以下是命令流，开挖之后隧道拱顶出现向上的位移，不知是哪里出错了，请各位大神指正，谢谢!!!!
gen zon radcyl p0 9 0 0 p1 18 0 0 p2 9 1 0 p3 9 0 9 &
size 4 1 8 4 dim 3 3 3 3 rat 1 1 1 1.2 group outr
gen zon cshell p0 9 0 0 p1 12 0 0 p2 9 1 0 p3 9 0 3 &
size 1 1 8 4 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1 group liny fill group inr
gen zon reflect dip 90 dd 90 orig 9 0 0
gen zon reflect dip 0 dd 0 orig 9 0 0
gen point id 1 -6 0 -3
gen point id 2 -6 1 -3
gen zon radcyl p0 -9 0 -3 p1 0 0 -3 p2 -9 1 -3 p3 -9 0 6 &
size 4 24 8 4 dim 3 3 3 3 rat 1 1 1 1.2 group outr
gen zon cshell p0 -9 0 -3 p1 point 1 p2 -9 1 -3 p3 -9 0 0 &
size 1 1 8 4 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1 group liny fill group inr
gen zon radcyl p0 -9 0 -3 p1 -9 0 -12 p2 -9 1 -3 p3 0 0 -3 &
size 4 1 8 4 dim 3 3 3 3 rat 1 1 1 1.2 group outr
gen zon cshell p0 -9 0 -3 p1 -9 0 -6 p2 -9 1 -3 p3 -6 0 -3 & ;这里不对
size 1 1 8 4 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1 group liny fill group inr
gen zon radcyl p0 -9 0 -3 p1 -9 0 6 p2 -9 1 -3 p3 -18 0 -3 &
size 4 1 8 4 dim 3 3 3 3 rat 1 1 1 1.2 group outr
gen zon cshell p0 -9 0 -3 p1 -9 0 0 p2 -9 1 -3 p3 -12 0 -3 &
size 1 1 8 4 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1 group liny fill group inr
gen zon radcyl p0 -9 0 -3 p1 -18 0 -3 p2 -9 1 -3 p3 -9 0 -12 &
size 4 24 8 4 dim 3 3 3 3 rat 1 1 1 1.2 group outr
gen zon cshell p0 -9 0 -3 p1 -12 0 -3 p2 -9 1 -3 p3 -9 0 -6 &
size 1 24 8 4 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1 group liny fill group inr
gen zon brick p0 18 0 -9 p1 36 0 -9 p2 18 1 -9 p3 18 0 9 size 10 1 8 group outr1
gen zon brick p0 18 0 9 p1 36 0 9 p2 18 1 9 p3 18 0 13 size 10 1 2 group outr2
gen zon brick p0 0 0 9 p1 18 0 9 p2 0 1 9 p3 0 0 13 size 8 12 2 group outr2
gen zon brick p0 18 0 -24 p1 36 0 -24 p2 18 1 -24 p3 18 0 -9 size 10 1 10 group outr3
gen zon brick p0 0 0 -24 p1 18 0 -24 p2 0 1 -24 p3 0 0 -9 size 10 1 10 group outr3
gen zon brick p0 -36 0 6 p1 0 0 6 p2 -36 1 6 p3 -36 0 13 size 18 12 3 group outr2
gen zon brick p0 -36 0 -12 p1 -18 0 -12 p2 -36 1 -12 p3 -36 0 6 size 10 1 8 group outr1
gen zon brick p0 -36 0 -24 p1 0 0 -24 p2 -36 1 -24 p3 -36 0 -12 size 18 1 8 group outr3
attach face
plo blo gro
set gravity 0 0 -10
fix z range z -24.01 -23.99
fix x range x -36.01 -35.99
fix x range x 35.99 36.01
fix y range y -0.01 0.01
fix y range y 0.99 1.01
model mohr
ini density 1800
prop bulk=3e6 shear=1e6 fric=20 coh=2.5e4
solve
plot cont szz
plot cont zdisp
save chushi.sav
ini xdisp=0 ydisp=0 zdisp=0
ini xvel=0 yvel=0 zvel=0
model null range group inr y 0 1
model elas range group liny y 0 1
ini density 2500 range group liny y 0 1
prop bulk=17e9 shear=6.8e9 range group liny y 0 1
solve
plot cont szz
plot cont zdisp
save step1.sav
