FLAC3D流固耦合例子

;导入模型
;==================================
def _prop ;4级围岩
c_el=9e9 ;弹性模量
c_bo=0.3 ;泊松比
c_bu=c_el/(3*(1-2*c_bo)) ;排水体积模量
c_sh=c_el/(2*(1+c_bo)) ;剪切模量
c_po=0.3 ;孔隙率
c_pe=2.22e-10
c_coh=1.5e6
c_fric=40.0
c_tension=1e6
c_dens=2400 ;密度
end
_prop
;注浆
def _zprop ;注浆
c_zel=14e9 ;弹性模量
c_zbo=0.3 ;泊松比
c_zbu=c_el/(3*(1-2*c_bo)) ;排水体积模量
c_zsh=c_el/(2*(1+c_bo)) ;剪切模量
c_zpo=0.3 ;孔隙率
c_zpe=2.22e-11
c_zcoh=1.5e6
c_zfric=40.0
c_ztension=1e6
c_zdens=2400
end
_zprop
;初期衬砌(弹性材料)
def ini_chu
c_elchu=10e9 ;弹性模量
c_bochu=0.30 ;泊松比
c_buchu=c_el/(3*(1-2*c_bo)) ;排水体积模量
c_shchu=c_el/(2*(1+c_bo)) ;剪切模量
c_pochu=0.2 ;孔隙率
c_pechu=1e-12
end
ini_chu
;============================================
;设置材料模型
model elastic
pro bulk c_bu shear=c_sh
ini dens c_dens
;============================================
;注浆静力模型
model elastic range group 7
pro bulk c_zbu shear=c_zsh range group 7
ini density c_zdens r g 7
;============================================
;设置边界条件
fix x range x -20.1 -19.9
fix x range x 19.9 20.1
fix y range y -0.01 0.01
fix y range y 0.59 0.61
fix z range z -15.1 -14.9
;============================================
config fluid
;设置初始应力状态，静力
set gravity 0 0 -10
;设置初始应力状态,水位4m
ini pp 10e4 grad 0 0 -1e4
model fl_iso
prop perm c_pe poros c_po biot_c 1
set fl biot off
ini fmod 2e9 fdensity 1e3 ftens -1e10
;注浆水力模型
model fl_iso range group 7
prop perm c_zpe poros c_zpo range group 7
;流体平衡
set mech off
set fluid on
pl con pp ou on
solve age 1
;力学平衡
ini fmod 0
set mech on
set fluid off
hist unbal
pl add hist 1
solve
;求解
ini fmod 2e9
set mech on fluid on
solve
save 初始应力.sav