602 岩 土 力 学 2011年

行迭代计算，计算结果如图1～5所示。

表1 Duncan-Chang本构模型试验参数

Table1 Duncan-Chang model parameters

c / kPa

φ / (°)

Rf K n Kb m Kur

1 186

338.8 36.83 0.81 741.3 0.24 199.5 0.26

表2 NHRI本构模型试验参数

Table 2 NHRI model parameters

c / kPa

φ

/ (°) / (°)

Rf K n Cb nd rd Kur864

0 56.25 12 0.7 540 0.3 0.01 0.6 0.66

图3 NHRI体应变与轴向应变模拟与试验结果 Fig.3 NHRI volumetric strain and axial strain

model and test results

表3 砂土本构模型试验参数 Table 3 Sand model parameters

弹性参数

G

临界状态参数

状态相关参数

µ

M eΓ λ ξ m D0 n h1

h2

2.55

22 0.05 1.75 0.512 0.025 0.7 1.28 0.61 1.1 2.60

图4 SAND偏应力与轴向应变模拟与试验结果

Fig.4 SAND deviatoric stress and axial strain

model and test results

图1 Duncan-Chang偏应力与轴向应变模拟与试验结果 Fig.1 Duncan-Chang deviatoric stress and axial strain

model and test result

图5 SAND体应变与轴向应变模拟与试验结果 Fig.5 SAND volumetric strain and axial strain

model and test results

从模拟过程、计算结果与试验结果比较看，本文二次开发子程序运行稳定，结果正确、可靠。图

1、2、4显示，拟合σa-εa曲线时，这3类模型模拟结果均较好，共同点是在剪切初始阶段模拟值偏

图2 NHRI偏应力与轴向应变模拟与试验结果 Fig.2 NHRI deviatoric stress and axial strain

model and test results

大而后期阶段偏小。由于Duncan-Chang模型为弹性非线性模型，在偏应力与轴向应变平面上，仅近似反映土体的非线性，不能反映土体的剪胀性。图