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=0.6e9 p_ratio1=0.27 E_mod2=0.8e9 p_ratio2=0.26 derive
prop bulk b_mod1 shear s_mod1 cohe 1.8e6 tens 0.8e6 fric 30 range z 4.5 20
prop bulk b_mod2 shear s_mod2 cohe 2.8e6 tens 1.0e6 fric 35 range z -40 4.5 ini dens=2300 set grav 0 0 -10
; boundary and initial conditions apply szz -1.4e6 range z 19.9 20.1 fix z range z -40.1 -39.1 fix x range x -45.1 -44.9 fix x range x 44.9 45.1 fix y range y 49.9 50.1 hist unbal
hist gp xdis 6.0,0,0 hist gp zdis 0,0,5
hist gp xdis 6.0,50,0 hist gp zdis 0,50,5 plot hist 3 solve
save tun_nature.sav
;对后面计算而言，模型建立时岩体在开挖前认为位移已经终了，因此需要对位移进行“清零”，而应力可以保留。假设隧道先开挖上断面，中间不设支护，直到进尺50m，那么此时位移和应力的分布情况可用如下命令流
ini xdis=0 ydis=0 zdis=0 plot re
model null range group 原岩 model null range group 二次衬砌 model null range group 初期支护 set large hist unbal plot hist 3 solve
save tun_ext1.sav
;如果在开挖后适时对隧道进行锚喷混凝土初期支护，该命令流如下： restore tun_nature.sav 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 plot hist 3
mo el range group 初期支护
prop bulk 1.33e9 shear 0.8e9 cohe 2.2e6 tens 1.2e6 fric 32 range z 4.5 7.0 x -7.0 7.0
prop bulk 10.9e9 shear 8.9e9 dens 2500 range group 初期支护 set large solve
save tun_ext2.sav
FLAC3D基础知识介绍解析
