newconfig fluidtitle基于流固耦合作用下的双龙富水隧道稳定性研宄set fluid offset log onset logfile yang 1 .loggenzonradcylpOOOOpl 9.00 0p2 0 50 0 p3 0 0 8 size4 20 64 dim6 5 6 5 rat 1 1 I 1 group围岩gen zon cshell pOOOOpl 6.0 0 0 p2 0 50 0 p3 0 0 5.0 size 4 2064 dim 5.6 4.6 5.6 4.6 rat 1 1 1 1 group初期支护gen zon cshell pO 0 0 0 pi 5.6 0 0 p2 0 50 0 p3 0 0 4.6 size 4 2064 dim 5.0 4.0 5.0 4.0 rat 1 1 1 1 group 二次衬砲fill group原岩gen zon radcyl pOOOOpl 0 0 ?8.0 p2 0 50 0 p3 9.0 0 0 size 4 20 64 dim 3 6 3 6 rat 1 1 1 1 group围岩2gen zon cshell pOOO Opl 0 0 -3.0 p2 0 50 0 p3 6.0 0 0 size 4 2064 dim 2.6 5.6 2.6 5.6 rat 1 1 1 1 group仰拱初期支护gen zon cshell pO 0 0 Opl 0 0 -2.6 p2 0 50 0 p3 5.6 0 0 size 4 2064 dim2 5 2 5 rat 1 1 1 1 group仰拱二次衬砲fill group仰拱原岩gen zone reflect normal -10 0gen zone radtun pO 0 0 0 pi 45 0 0 p2 0 50 0 p3 0 0 20 size 3 20 3 12 dim 9 8 9 8 rat 1 1 1 1.1 group围岩3gen zon reflect dip 0 ori 0 0 0 range x09y0 50z8 20gen zon reflect dip 0 ori 0 0 0 range x 9 45 y 0 50 z 0 20gen zon reflect dip 90 dd 270 ori 0 0 0 range x09y0 50z8 20gen zon reflect dip 90 dd 270 ori 0 0 0 range x 0 9 y 0 50 z -8 -20gen zon reflect dip 90 dd 270 ori 0 0 0 range x 9 45 y 0 50 z -20 20save shuitun一model.savmodel fl_isoprop perm 1.23e-9 poro 0.45 range z 4.5 20prop perm 4.70e-10 poro 0.4 range z -20 4.5set fl biot offini fdensity le3ini sat 1.0ini food 2.0e9 ftens -le-3;假设围岩岩体符合mohr-coulomb本构模型,给围岩陚参数命令流如下,;mohr-coulomb modelmodel mohrdef derives_modl=E一modl/(2.0*(l .0+p—ratio 1))b_modl=E_modl/(3.0*(1.0-2.0*p_ratiol))s_mod2=E—mod2/(2.0*(l .0+p_ratio2))bjmod2~E_mod2/(3.0*( 1.0-2.0*p_ratio2))endset E_modl=0.6e9 p_ratio 1=0.27 E_mod2=0.8e9 p__ratio2=0.26deriveprop bulk b_jnodl shear s_modl cohe 1.8e6 tens 0.8e6 firic 30 range z 4.5 20 prop bulk b_jnod2 shear s_mod2 cohe 2.8e6 tens 1.0e6 fric 35 range z -20 4.5ini dens=2300set grav 0 0 -10;boundary and initial conditionsapply szz -1.4e6 range z 19.9 20.1ini pp 0 grad 0 0 -1.0e4fix pp 0 range z 19.9 20.1fix z range z -20.1 -19.9fix x range x -45.1 -44.9fix x range x 44.9 45.1fix y range y 49.9 50.1hist unbalhist gp xdis 6.0,0,0hist gp zdis 0,0,5hist gp xdis 6.0,50,0hist gp zdis 0,50,5hist gp pp 6 0 0hist gp pp 0 0 5hist gp disp id=82hist gp disp id=162hist gp disp id=4243hist gp disp id=l 1621plot hist 3solvesave shuitun_nature.savrestore shuitun_nature.savset fluid onini xdis=0 ydis=0 zdis=0set fluid onset mech onplot remodel null range group 原岩any group 二次衬砲any group 初期支护any model fl_null range group 原岩any group 二次衬砌any group 初期支护any fix pp 0 range group 原岩any group 二次衬砌any group 初期支护anyhist unbalsolvesave shuitun_extl .savnewtitle基于流固耦合作用下的双龙富水隧道稳定性研究set log onset logfile yang.loggenzonradcylpO000pi 9.000p20500p3 008 size4 2064dim6565rat 1111 group围岩gen zon cshell pO 0 0 0 pi 6.0 0 0 p2 0 50 0 p3 0 0 5.0 size 4 2064 dim 5.6 4.6 5.6 4.6 rat 1 1 1 1 group初期支护gen zon cshell pO 0 0 0 pi 5.6 0 0 p2 0 50 0 p3 0 0 4.6 size 4 2064 dim 5.0 4.0 5.0 4.0 rat 1 1 1 1 group 二次衬御fill group 原岩gen zon radcyl pOOOOpl 0 0 -8.0 p2 0 50 0 p3 9.0 0 0 size 4 2064 dim 3 6 3 6 rat 1 1 1 1 group围岩2gen zon cshell pOOOOpl 0 0 -3.0 p2 0 50 0 p3 6.0 0 0 size 4 20 6 4 dim 2.6 5.6 2.6 5.6 rat 1 1 1 1 group仰拱初期支护gen zon cshell pOOOOpl 0 0 -2.6 p2 0 50 0 p3 5.6 0 0 size 4 20 6 4 dim 2 5 2 5 rat 1 1 1 1 group仰拱二次衬砲fill group仰拱原岩gen zone reflect normal ?10 0gen zone radtun pO 0 0 0 pi 45 0 0 p2 0 50 0 p3 0 0 20 size 3 20 3 12 dim 9898 rat 1 1 1 1.1 group围岩3gen zon reflect dip 0 on 0 0 0 range x09y0 50z8 20gen zon reflect dip 0 ori 0 0 0 range x 9 45 y 0 50 z 0 20gen zon reflect dip 90 dd 270 ori 0 0 0 range x09y0 50z8 20gen zon reflect dip 90 dd 270 ori 0 0 0 range x 0 9 y 0 50 z -8 -20gen zon reflect dip 90 dd 270 ori 0 0 0 range x 9 45 y 0 50 z -20 20save tun一model.sav;假设围岩岩体符合mohr?coulomb本构模型,给围岩赋参数命令流如下,;mohr-coulomb modelmodel mohrdef derives一modl=E一modl/(2.0*(1.0+p一ratiol))b一mod 1 =E_modl/(3.0*( 1.0-2.0*p一ratio 1))s_mod2=E_mod2/(2.0*(l .0+p_ratio2))b_mod2=E_mod2/(3.0*(1.0-2.0*p_ratio2))endset E_modl^0.6e9 p_ratiol=0.27 E_mod2=0.8e9 p_ratio2=0.26deriveprop bulk b_modl shear sjtnodl cohe 1.8e6 tens 0.8e6 fric 30 range z 4.5 20prop bulk b_mod2 shear s_mod2 cohe 2.8e6 tens 1.0e6 fric 35 range z -20 4.5ini dens=2300set grav 0 0 -10;boundary and initial conditionsfix z range z -20.1 -19.9fix x range x -45,1 -44.9fix x range x 44.9 45.1fix y range y 49.9 50.1hist unbalhist gp disp id=82 hist gp disp id=162 hist gp disp id=4243 hist gp disp id=l 1621 plot hist 3solvesave tunjiature.sav。