1、
隧洞埋深取80
米，建立50*50*50模型，隧洞中心点坐标为（0 0 -25） 模型顶部施加上部岩体自重应力作为顶部边界条件，暂取地下水为充满整个模型，命令流如下，为什么施加衬砌之后计算难以收敛，且查看衬砌的smax最大主应力超过几十兆帕，远远超于极限抗拉强度，是哪里出错了？哪个兄弟帮忙看下啊，或者提供下可能错在哪里的思路 急用，万分感谢。调试好多遍模型了，还是得不到理想的结果，谢谢了，可能对流固耦合了解不够。config fluid
;----------
流体模型参数----------model fl_iso
prop perm 4.08e10-11 poro 0.3
set fl biot off
ini fdensity 1e3
ini sat 1.0
ini fmod 2e9 ftens -1e-3
ini pp 0 grad 0 0 -1e4 range z -50 0
fix pp 0 range z -0.1 0.1
;----------
力学参数----------def derive0
b_mod0=Emod/(3.0*(1.0-2.0*p_ratio))
s_mod0=Emod/(2.0*(1.0+p_ratio))
end
set Emod 8e9
set p_ratio 0.25
derive0
model mohr
pro den 2400 bulk b_mod0 shear s_mod0 coh 8e5 fric 40
;----------
位移边界条件----------fix x range x -24.9 -25.1
fix x range x 24.9 25.1
fix y range y -0.1 0.1
fix y range y 49.9 50.1
fix x y z range z -49.9 -50.1
;----------
初始条件----------apply nstress -1.5e6 range z -0.1 0.1
ini szz -1.5e6 grad 0 0 2.6e4 range z -50 0
ini sxx -1.5e6 grad 0 0 2.6e4 range z -50 0 ;
侧压力系数k取0.8ini syy -1.5e6 grad 0 0 2.6e4 range z -50 0
set grav 0 0 -10
set mech ratio 5e-5
Solve
ini xdis=0 ydis=0 zdis=0
ini xvel=0 yvel=0 zvel=0
;----------
隧洞开挖及施加反力模拟应力释放---------- model null range gr
隧洞 any gr 衬砌1 any gr 衬砌2 any ;----------
应力释放---------set fluid off mech on
fix pp
def sup_stress
lumda=0.7
o_press=-1.36e6
s_press=o_press*lumda
o_grad=2.4e4
s_grad=o_grad*lumda
end
sup_stress
apply nstress s_press grad 0 0 s_grad range x -3.6 3.6 z -28.6 -21.4
set mech ratio 1e-4
Solve
衬砌施作apply nstress 0 range x -3.6 3.6 z -28.6 -21.4
set fluid on mech on
free pp
model fl_null range gr
隧洞 any gr 衬砌1 any gr 衬砌2 any fix pp 0 range gr
隧洞 any gr 衬砌1 any gr 衬砌2 any model el range gr
衬砌1 any gr 衬砌2 any model fl_iso range gr
衬砌1 any gr 衬砌2 any prop bulk 4.642e8 shear 8.974e8 poro 0.05 perm 4.08e10-12 range gr
衬砌1 any gr 衬砌2 any set mech ratio 1e-4
solve
