坝基的位移向上！

xfjt

加了孔隙水压力后，坝基的位移向上而且位移量级很大，不知道原因请大侠们指点！命令流如下：
config fl
model mohr range group 1
prop bulk 1.8e7 shear 1e7 coh 1e3 fric 30 range group 1
mode mohr range group 2
prop bulk 1.2e7 shear 9e6 coh 6.6e3 fric 26 range group 2
mode mohr range group 3
prop bulk 1.4e7 shear 1.1e7 coh 8.3e3 fric 28 range group 3
mode mohr range group 4
prop bulk 6e9 shear 8e8 coh 1.2e4 fric 32 range group 4
ini density 1570 range group 1
ini density 1530 range group 2
ini density 1650 range group 3
ini density 2100 range group 4
mode fl_iso
prop perm 1.6e-10 poros 0.31 range group 1
prop perm 3.7e-10 poros 0.34 range group 2
prop perm 3.7e-10 poros 0.34 range group 3
prop perm 3e-12  poros 0.17 range group 4
ini fmod 2e8
set gravity 0 -10 0
ini fdens 1000 ftens=0
ini sat 1 range y 0 13
ini sat 0 range y 13 16.87
ini pp 1.3e5 grad 0 -10000 0 range y 0 13
;ini pp 0   range y 13 16.87
;fix pp  range x -0.01 0.01
;fix pp  range x 159.99 160.01
fix pp  range z -50.75 -50.73
fix pp  range z 110.73 110.75
;fix pp  range y 12.9 13.1
;fix pp 1.3e5 range y -0.01 0.01
fix x range x -0.01 0.01
fix x range x 159.99 160.01
fix z range z -50.75 -50.73
fix z range z 110.73 110.75
fix y range y -0.01 0.01
set gravity 0 -10 0
set fluid on mech off
;set mech ratio 1e-4
step 1000
save gongkuang-1-1.sav
ini fmod 0
set fl off mech on
solve

xfjt

如果不加 set fl on mech off 直接set fl off mech on 计算，结果是这样的

xfjt

差别怎么这么大呀？是那里的问题？

fangzhenkeji

肯定是存在超孔压