请教：深埋 隧道 初始应力场（S-B法），结果不理想。

fanyonghui618

请教：深埋 隧道 初始应力场（S-B法），埋深170米，S-B法模拟初始应力，结果不理想。

set log on

gen zon radcyl p0 0 0 0 p1 12 0 0 p2 0 2 0 p3 0 0 11 size 6 2 14 6 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 14 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 3 2 14 8 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 -9 p2 0 2 0 p3 12 0 0 size 6 2 14 6 dim 5 8.5 5 8.5 rat 1 1 1 1 group 围岩2
gen zon cshell p0 0 0 0 p1 0 0 -5 p2 0 2 0 p3 8.5 0 0 size 2 2 14 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 3 2 14 8 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 40 0 0 p2 0 2 0 p3 0 0 30 size 7 2 7 18 dim 12 11 12 11 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 -40 0 0 p2 0 2 0 p3 0 0 -30 size 7 2 7 18 dim 12 9 12 9 rat 1 1 1 1.05 group 围岩4
gen zone reflect normal -1 0 0 range group 围岩4

gen zon brick p0 40 0 0 p1 60 0 0 p2 40 2 0 p3 40 0 30 size 7 2 7 rat 1.05 1 1 group 围岩5
gen zone reflect normal 0 0 -1 range group 围岩5
gen zone reflect normal -1 0 0 range group 围岩5
gen zon brick p0 0 0 30 p1 40 0 30 p2 0 2 30 p3 0 0 170 size 7 2 35 rat 1 1 1.05 group 围岩6
gen zone reflect normal -1 0 0 range group 围岩6
gen zon brick p0 40 0 30 p1 60 0 30 p2 40 2 30 p3 40 0 170 size 7 2 35 rat 1.05 1 1.05 group 围岩7
gen zone reflect normal -1 0 0 range group 围岩7
gen zon brick p0 0 0 -30 p1 0 0 -50 p2 0 2 -30 p3 40 0 -30 size 8 2 7 rat 1.05 1 1 group 围岩8
gen zone reflect normal -1 0 0 range group 围岩8
gen zon brick p0 40 0 -30 p1 40 0 -50 p2 40 2 -30 p3 60 0 -30 size 8 2 7 rat 1.05 1 1.05 group 围岩9
gen zone reflect normal -1 0 0 range group 围岩9

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))
end

set E_mod1=11e9 p_ratio1=0.33
derive

pro bulk b_mod1 shear s_mod1 friction 25 coh 1e15 ten 1e15 range z -50 170 x -60 60
;prop bulk 5e8 shear 1e9 coh 1e6 tens 0.3e6 fric 30 density 2000 range z -120 170 x -60 60
ini density 2700
set grav 0 0 -10
;边界条件

;初始应力

apply sxx -1e9 grad 0 0 -450e5 range x -59.9 -60.1
apply sxx -1e9 grad 0 0 -450e5 range x 59.9 60.1
apply syy -1e9 grad 0 0 -450e5 range y -0.1 0.1
apply syy -1e9 grad 0 0 -450e5 range y 1.9 2.1
apply szz -1e8 grad 0 0 -450e4 range z -50 170


history unbal
step 9000
plot hist 1

mo m
pro bulk b_mod1 shear s_mod1 fric 25 cohe 2.25e5 tens 0.18e3 range z -50 170 x -60 60
ini density 2700

step 800
ini state 0
ini xdis=0 ydis=0 zdis=0
ini xvel=0 yvel=0 zvel=0
save qinling.sav

goodmangis

不理想是指什么？
apply 中grad 0 0 -450e4 为负值，会造成越往上越大，结果好像不对。另外，z的梯度递增从0开始的吧