流固耦合ansys命令流

new

config fluid

title

基于流固耦合作用下的双龙富水隧道稳定性研宄

set fluid off

set log on

set logfile yang 1 .log

genzonradcylpOOOOpl 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围岩2

gen 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 0

gen 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围岩3

gen zon reflect dip 0 ori 0 0 0 range x09y0 50z8 20

gen zon reflect dip 0 ori 0 0 0 range x 9 45 y 0 50 z 0 20

gen zon reflect dip 90 dd 270 ori 0 0 0 range x09y0 50z8 20

gen zon reflect dip 90 dd 270 ori 0 0 0 range x 0 9 y 0 50 z -8 -20

gen zon reflect dip 90 dd 270 ori 0 0 0 range x 9 45 y 0 50 z -20 20

save shuitun一model.sav

model fl_iso

prop perm 1.23e-9 poro 0.45 range z 4.5 20

prop perm 4.70e-10 poro 0.4 range z -20 4.5

set fl biot off

ini fdensity le3

ini sat 1.0

ini food 2.0e9 ftens -le-3

;假设围岩岩体符合mohr-coulomb本构模型,给围岩陚参数命令流如下,

;mohr-coulomb model

model mohr

def derive

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

end

set E_modl=0.6e9 p_ratio 1=0.27 E_mod2=0.8e9 p__ratio2=0.26

derive

prop 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.5

ini dens=2300

set grav 0 0 -10

;boundary and initial conditions

apply szz -1.4e6 range z 19.9 20.1

ini pp 0 grad 0 0 -1.0e4

fix pp 0 range z 19.9 20.1

fix z range z -20.1 -19.9

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

hist gp pp 6 0 0

hist gp pp 0 0 5

hist gp disp id=82

hist gp disp id=162

hist gp disp id=4243

hist gp disp id=l 1621

plot hist 3

solve

save shuitun_nature.sav

restore shuitun_nature.sav

set fluid on

ini xdis=0 ydis=0 zdis=0

set fluid on

set mech on

plot re

model 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 初期支护any

hist unbal

solve

save shuitun_extl .sav

new

title