new
gen zon radcyl p0 0 0 0 p1 30 0 0 p2 0 50 0 p3 0 0 30 &
size 15 25 15 8 dim 7.4 8 7.4 8 rat 1 1 1 1.2 group 围岩
gen zon cshell p0 0 0 0 p1 7.4 0 0 p2 0 50 0 p3 0 0 8 &
size 15 25 15 8 dim 7.1 7.7 7.1 7.7 rat 1 1 1 1 group 初期支护
gen zon cshell p0 0 0 0 p1 7.1 0 0 p2 0 50 0 p3 0 0 7.7 &
size 15 25 15 8 dim 6.5 7.1 6.5 7.1 rat 1 1 1 1 group 二次衬砌 fill group 原岩
gen zon radcyl p0 0 0 0 p1 0 0 -30 p2 0 50 0 p3 30 0 0 &
size 15 25 15 8 dim 4 7.4 4 7.4 rat 1 1 1 1.2 group 围岩2
gen zon cshell p0 0 0 0 p1 0 0 -4.0 p2 0 50 0 p3 7.4 0 0 &
size 15 25 15 8 dim 3.7 7.1 3.7 7.1 rat 1 1 1 1 group 仰拱初期支护
gen zon cshell p0 0 0 0 p1 0 0 -3.7 p2 0 50 0 p3 7.1 0 0 &
size 15 25 15 8 dim 3.1 6.5 3.1 6.5 rat 1 1 1 1 group 仰拱二次衬砌 fill group 仰拱原岩
gen zone reflect normal -1 0 0
save tunnel-model1.sav
plot block group
model mohr ;莫尔～库仑模型
ini density 1960 ;围岩的密度
prop bulk=2.5e9 shear=3.78e8 fric=26 coh=3.04e5 tension=8.15e6
set gravity 0 0 -10 ;设置重力加速度为z方向-10
apply szz -3.53e6 range z 29.9 30.1
fix z range z -30.1,-29.9 ;设置底边界
fix x range x 30.1,29.9 ;设置左边界
fix x range x -29.9, -30.1 ;设置右边界
fix y range y -0.1 0.1 ;设置前边界
fix y range y 50.1,49.9 ;设置后边界
hist unbal
step 4000 ;求解计算
save gsol1.sav ;计算结果保存在Gsol.sav文件中
plot block group
ini xdis=0 ydis=0 zdis=0
plot re
model null range group 原岩
model null range group 二次衬砌
model null range group 初期支护
step 100
hist unbal
mo el range group 初期支护
prop bulk 1.43e10 shear 1.3e10 fric 70 coh 3e6 dens 2200 tension 5e6 range group 初期支护(这条命令输入后就会显示warning-Property not available to any zone in range，求高手帮助万分感激)
