;给土体赋弹性本构属性
model elastic range group 1
prop bulk 27.436e6 she 17.258e6 dens 2009
model elastic range group 2
prop bulk 28.205e6 she 17.742e6 dens 2115
model elastic range group 3
prop bulk 2.703e10 she 2.655e10 dens 2550
;固定边界
fix x range x -0.1 0.1
fix x range x 619.9 620.1
fix x y z range y -0.1 0.1
fix z range z -0.1 0.1
fix z range z -300.1 -299.9
;施加重力
set grav 0 -10 0
hist unb
hist gp xvel 200 80 -150
hist gp yvel 200 80 -150；无图形显示
hist gp zvel 200 80 -150；无图形显示
hist gp zdisp 200 80 -150；无图形显示
plot hist 1
plot hist 2
plot hist 3
plot hist 4
plot hist 5
solve
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
save 4-24chushi pingheng.sav
restore 4-24chushi pingheng.sav
；给土体赋摩尔库伦本构参数
model mohr range group 1
prop bulk 27.436e6 she 17.258e6 friction 33.2 coh 60e3
model mohr range group 2
prop bulk 28.205e6 she 17.742e6 friction 33.1 coh 76e3
model mohr range group 3
prop bulk 2.703e10 she 2.655e10 friction 50 coh 55e3
