分层开挖，挖一层打一层土钉，然后计算到第五层的时候就开始出问题了，坑壁水平位移太大，好像土钉没加似得，帮我看看呗
new
gen zone brick p0 40 -30 -39 p1 100 -30 -39 p2 40 0 -39 p3 40 -30 -18.5 size 30 15 5 group b8
gen zone brick p0 40 -30 -18.5 p1 100 -30 -18.5 p2 40 0 -18.5 p3 40 -30 -14 size 30 15 9 group b7
gen zone brick p0 40 -30 -14 p1 100 -30 -14 p2 40 0 -14 p3 40 -30 -10.5 size 30 15 7 group b6
gen zone brick p0 40 -30 -10.5 p1 100 -30 -10.5 p2 40 0 -10.5 p3 40 -30 -9 size 30 15 3 group b5
gen zone brick p0 40 -30 -9 p1 100 -30 -9 p2 40 0 -9 p3 40 -30 -6 size 30 15 6 group b4
gen zone brick p0 40 -30 -6 p1 100 -30 -6 p2 40 0 -6 p3 40 -30 -3 size 30 15 6 group b3
gen zone brick p0 40 -30 -3 p1 100 -30 -3 p2 40 0 -3 p3 40 -30 -2 size 30 15 2 group b2
gen zone brick p0 40 -30 -2 p1 100 -30 -2 p2 40 0 -2 p3 40 -30 0 size 30 15 4 group b1
plot block group
fix x range x 39.9 40.1
fix x range x 99.9 100.1
fix x y z range z -38.9 -39.1
fix y range y -0.1 0.1
fix y range y -29.9 -30.1
model ela
pro bulk 6.7308e7 shear 4.2339e7 range group b1
pro bulk 5.8333e7 shear 3.3333e7 range group b2
pro bulk 5.2083e7 shear 2.9762e7 range group b3
pro bulk 2.5e7 shear 1.1539e7 range group b4
pro bulk 2.5e7 shear 1.1539e7 range group b5
pro bulk 2.25e7 shear 1.0385e7 range group b6
pro bulk 1.27e7 shear 4.7e6 range group b7
pro bulk 1.27e7 shear 4.7e6 range group b8
ini density 2000 range group b1
ini density 2000 range group b2
ini density 2000 range group b3
ini density 1990 range group b4
ini density 1990 range group b5
ini density 1970 range group b6
ini density 2000 range group b7
ini density 2000 range group b8
set gravity 0 0 -10
plot add contour szz
solve
save mod1.sav
model mohr
pro bulk 6.7308e7 shear 4.2339e7 fric 38 cohesion 0 range group b1
pro bulk 5.8333e7 shear 3.3333e7 fric 32 cohesion 0 range group b2
pro bulk 5.2083e7 shear 2.9762e7 fric 30 cohesion 0 range group b3
pro bulk 2.5e7 shear 1.1539e7 fric 23 cohesion 2.5e4 range group b4
pro bulk 2.5e7 shear 1.1539e7 fric 15 cohesion 3e4 range group b5
pro bulk 2.25e7 shear 1.0385e7 fric 25 cohesion 2.4e4 range group b6
pro bulk 1.27e7 shear 4.7e6 fric 18 cohesion 5e3 range group b7
pro bulk 1.27e7 shear 4.7e6 fric 12 cohesion 1e4 range group b8
plot block group
plot add contour szz
solve
；第一层开挖
new
;初始节点三方向位移0，四面体塑性标志0（无塑性破坏）
ini xdisp 0 ydisp 0 zdisp 0
ini state 0
model null range x 80 100 z 0 -1
def ins_cable1
y1=0
loop i(101,121)
x1=80
x2=68.724
z1=-1
z2=-5.1
command
sel cable id=i begin x1 y1 z1 end x2 y1 z2 nseg 12
sel cable id=i pro emod 2e11 xcarea 4.9e-4 gr_per 4.71e-1 yten 1.47e5 gr_k 7e6 gr_c 2e4 range cid 1 12
;sel cable id=i pretension 9.8e4 range cid 1 12
endcommand
y1=y1-1.5
endloop
end
ins_cable1
;sel shell id=1001 range x 79.9 80.1 y 0 -30 z 0 -1
;sel node fix y z yr range x 79.9 80.1
;sel shell id=1001 prop iso=(2.5e10,0.20) thick=0.1
plot block group
plot add contour szz
solve
下面都跟第一层差不多，只是z坐标改下，每次挖1.5米深，挖到10米，你们看看这是咋回事
