下面是我的数值模型：
new
gen zo radc p0 (0,0,0) p1 (36.3,0,0) p2 (0,35,0) p3 (0,0,8)&
d (6.3,5.3,6.3,5.3) size 6 7 6 5 rat 1 1 1 1 group 围岩1
gen zo cs p0 (0,0,0) p1 6.3 0 0 p2 0 35 0 p3 0 0 5.3 &
d 6 5 6 5 size 3 7 6 6 rat 1 1 1 1 gro 初期支护 fill group 原岩
gen zo reflect nor (-1,0,0) origin (0,0,0)
gen zo b p0 10 0 8 p1 36.3 0 8 p2 10 35 8 p3 10 0 20 p4 36.3 35 8 &
p5 10 35 25 p6 36.3 0 30 p7 36.3 35 35 size 5 7 5 rat 1 1 1 gro 围岩2
gen zo w p0 (10,35,8) p1 (-36.3,35,8) p2 (10,0,8) p3 (10,35,25) p4 (-36.3 0 8) &
p5 10 0 20 size 18 7 6 rat 1 1 1 group 仰坡
gen zo b p0 (36.3,35,0) p1 (36.3,0,0) p2 (-36.3,35,0) p3 (36.3,35,-15) size 7 20 5 rat 1 1 1 &
group 围岩3
attach face ran z 7.9 8.1
attach face ran x -36.3 36.3
attach face ran z -0.1 0.1
save 模型.sav
model mohr
pro bulk 1e7 shear 3.34e6 coh 20e10 fr 15 ten 0
fix z x y ran z -15
fix x ran x -36.3
fix x ran x 36.3
fix y ran y 0
fix y ran y 35
ini dens 1600
set grav 0 0 -10
solve elas
save 初始应力场.sav
ini xdis 0 ydis 0 zdis 0
model null range gro 初期支护
model null range gro 原岩
step 100
save 开挖后.sav
model ela ran gro 初期支护
pro bulk 2.5e10 shear 1.875e10 ran gro 初期支护
hist unbal
plo hi 1
solve
save 最终.sav
