config fluid
impgrid flacmodel_haitang.flac3d
set fluid off
;初始应力
pl bl g
pause
model fl_iso
model elas
prop bulk 5.1e12 shear 2.35e12
pro por 0.5 perm 1.34e-11 range group 1 any group 2 any
;初始化材料密度
ini dens 2100 range group 1
ini dens 1860 range group 2
ini dens 2460 range group 3
ini sat 0.15 range group 1 any group 2 any
ini sat 0 range group 3
;设置边界条件
fix x y z range y -50.1 -49.9
fix z range z -0.1 0.1
fix x range x -0.1 0.1
fix z range z 149.9 150.1
fix x range x 199.9 200.1
set gravity 0 -10.0 0
history unb
pl con syy
solve
save chushiyingli.sav
pause
n
set log on
;初始孔压
rest chushiyingli.sav
;位移和速度归零
ini xdisp 0
ini ydisp 0
ini zdisp 0
ini xvel 0
ini yvel 0
ini zvel 0
config fluid
set fluid on mech off
model fl_iso
prop poros 0.5 perm 1.34e-11 range group 1 any group 2 any
model fl_null range group 3
ini fden 1000
ini fmod 2e3 ftens -1e5
ini sat 0 range group 3
apply pp 50e4 grad 0 10e3 0 range x 0 200 y -50 0 z 0
pl con pp ou on
pl add flow scale 0.08 red
set grav 0 -10 0
pause
step 400000
save chushikongya.sav
算完后group 3的饱和度怎么会变成1了呢,不是已经设为0了吗