这是我的第一步开挖程序：
new
title nail-wall
plot ske
gen zone brick size 60 2 40 p0 0 0 0 p1 30 0 0 p2 0 1 0 p3 0 0 20
model mohr
prop bulk 40.9e6 shear 18.9e6 coh 15e3 fric 25 dil=15
ini dens 1900
fix x z range z -0.1 0.1
fix x rang x -0.1 0.1
fix x range x 29.9 30.1
fix y
ini szz=-3.8e5 grad 0,0,1.9e4
ini sxx=-1.04e5 grad 0,0,0.57e4
ini syy=-1.04e5 grad 0,0,0.57e4
set grav 0,0,-10
hist unbal
def laycable
x1=6.1*cos((3.14*20)/180)+9.0
z1=19.5-sin((3.14*20)/180)
z2=17.7-sin((3.14*20)/180)
z3=15.9-sin((3.14*20)/180)
z4=13.1-sin((3.14*20)/180)
z5=11.3-sin((3.14*20)/180)
end
laycable
solve
save step0.sav
set large
ini stat 0
ini xdis 0 ydis 0 zdis 0
ini xve 0 yve 0 zve 0
;第一层开挖
mod null range x 0 9 z 18.2 20
sel cable id=10 begin=(9.0,0.5,19.5) end=(x1,0.5,z1) nseg=12
sel cable id=10 prop xcarea=8.17e-3 emod=38.25e9 ytens=3.4e8 gr_k=30e6 gr_coh=10e3 gr_per=0.32 gr_fric=20
sel shell id=10 range x 8.99 9.01 z 18.2 20
sel shell id=10 prop iso=(26e9,0.3) thick 0.1
solve
save step1.sav
计算的结果，土钉出现了压应力
