做了一个深埋隧道的模型，120米埋深，高地应力，上边界至地表，侧压力系数1.5，solve后地表应力不为0，貌似不对啊，求改正！
set log on
gen zon radcyl p0 0 0 0 p1 15 0 0 p2 0 2 0 p3 0 0 13.5 size 6 2 16 6 dim 10.5 9 10.5 9 rat 1 1 1 1 group 围岩
gen zon cshell p0 0 0 0 p1 10.5 0 0 p2 0 2 0 p3 0 0 9 size 4 2 16 4 dim 8.5 7 8.5 7 rat 1 1 1 1 group 预支护
gen zon cshell p0 0 0 0 p1 8.5 0 0 p2 0 2 0 p3 0 0 7 size 2 2 16 2 dim 8.2 6.7 8.2 6.7 rat 1 1 1 1 group 初期支护
gen zon cshell p0 0 0 0 p1 8.2 0 0 p2 0 2 0 p3 0 0 6.7 size 2 2 16 10 dim 7.6 6.1 7.6 6.1 rat 1 1 1 1 group 二次衬砌 fill group 原岩
gen zon radcyl p0 0 0 0 p1 0 0 -11.5 p2 0 2 0 p3 15 0 0 size 6 2 16 6 dim 7 10.5 7 10.5 rat 1 1 1 1 group 围岩2
gen zon cshell p0 0 0 0 p1 0 0 -7 p2 0 2 0 p3 10.5 0 0 size 4 2 16 4 dim 5 8.5 5 8.5 rat 1 1 1 1 group 仰拱预支护
gen zon cshell p0 0 0 0 p1 0 0 -5 p2 0 2 0 p3 8.5 0 0 size 2 2 16 4 dim 4.7 8.2 4.7 8.2 rat 1 1 1 1 group 仰拱初期支护
gen zon cshell p0 0 0 0 p1 0 0 -4.7 p2 0 2 0 p3 8.2 0 0 size 2 2 16 10 dim 4.1 7.6 4.1 7.6 rat 1 1 1 1 group 仰拱二次衬砌 fill group 仰拱原岩
gen zone reflect normal -1 0 0
gen zone radtun p0 0 0 0 p1 60 0 0 p2 0 2 0 p3 0 0 60 size 8 2 8 20 dim 15 13.5 15 13.5 rat 1 1 1 1.05 group 围岩3
gen zone reflect normal -1 0 0 range group 围岩3
gen zone radtun p0 0 0 0 p1 -60 0 0 p2 0 2 0 p3 0 0 -60 size 8 2 8 20 dim 15 11.5 15 11.5 rat 1 1 1 1.05 group 围岩4
gen zone reflect normal -1 0 0 range group 围岩4
gen zon brick p0 0 0 60 p1 60 0 60 p2 0 2 60 p3 0 0 120 size 8 2 20 rat 1 1 1 group 围岩5
gen zone reflect normal -1 0 0 range group 围岩5
save tun_model.sav
plot show
;
;设定材料参数
;mohr-coulomb model
mo m
def derive
s_mod1=E_mod1/(2.0*(1.0+p_ratio1))
b_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio1))
s_mod2=E_mod2/(2.0*(1.0+p_ratio2))
b_mod2=E_mod2/(3.0*(1.0-2.0*p_ratio2))
end
set E_mod1=3.6e9 p_ratio1=0.3
set E_mod2=19.3e9 p_ratio2=0.38
derive
pro bulk b_mod1 shear s_mod1 fric 30 cohe 1e10 tens 1e10 range z -60 120 x -60 60;
in density 2700
set grav 0 0 -10
;边界条件
fix x range x -60.1 -59.9
fix x range x 59.9 60.1
fix z range z -59.9 -60.1
fix y range y -0.1 0.1
fix y range y 1.9 2.1
;初始应力
;k0=1.5
ini szz -3240000 grad 0,0,27000
ini sxx -4860000 grad 0,0,40500
ini syy -4860000 grad 0,0,40500
apply sxx -4860000 grad 0,0,40500 range x -60.1 -59.9
apply sxx -4860000 grad 0,0,40500 range x -0.1 0.1
apply szz -4.86e6 range z -60.1 -59.9;底部
apply syy -4860000 grad 0,0,40500 range y 1.9 2.1
apply syy -4860000 grad 0,0,40500 range y -0.1 0.1
history unbal
set mech ratio 1.0e-5
plot hist 1
solve
mo m
pro bulk b_mod1 shear s_mod1 fric 30 cohe 0.224e6 tens 0.5e6 density 2700 range z -60 120 x -60 60
ini state 0
ini xdis=0 ydis=0 zdis=0
ini xvel 0 yvel 0 zvel 0
solve
save qinling.sav
plot surf yellow fap red
