FLAC3D流固耦合急急急！！！

渗流场出错！！！！
命令流如下：
plot set rotation 20 0 30
plo blo gro
config fluid
model fl_iso
prop poros 0.3 perm 4.4574e-10
ini fmod 2e9
ini fdensity 1e3
ini pp 49.5e4 gra 0 0 -1e4
fix pp ran z 49.4 49.6
model elas
prop bulk 1e10 shear 1e10
ini dens=2200
fix x range x -77.6 -77.4
fix x range x 42.4 42.6
fix y range y -.1 .1
fix y range y 29.9 30.1
fix x y z range z -49.6 -49.4
set grav 0 0 -9.81
;网格三维显示
plot set rotation 30 0 30
plot cont zdisp outl on
hist unbal
set mech ratio 1e-4
solve
save block.sav
ini xdis 0.0 ydis 0.0 zdis 0.0
def derive
s_mod=E_mod/(2.0*(1.0+p_ratio))
b_mod=E_mod/(3.0*(1.0-2.0*p_ratio))
end
set E_mod=5e9 p_ratio=0.3
derive
model mohr
prop bulk b_mod shear s_mod cohe 0.5e6 fric 35 tens 0.5e6
model null ran gro outsiderock not y 0 15
ini pp 0.0 ran gro outsiderock not y 0 15
fix pp ran gro outsiderock not y 0 15

plo con pp
solve
save seepage001.sav