FLAC3D 隧道渗流

在开挖计算以前位移已清零，为何位移会是这样？？？？

config fluid
model fl_iso
prop poros 0.3 perm 4.4574e-10
ini fmod 2e3
ini fdensity 1e3
ini pp -100e4 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
ini szz -4400e3 gra 0 0 22e3
ini syy -2200e3 gra 0 0 11e3
ini sxx -2200e3 gra 0 0 11e3
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 block1.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
model null ran gro excavation any gro liner any
ini pp 0.0 ran gro excavation any gro liner any
fix pp ran gro excavation any gro liner any
plo con pp
solve
save seeapge001.sav