flac3d计算实例时从固定边界流下来啦！

我刚学flac3d不到一个星期，试着按照实例建了一个模型，可是计算的时候竟然发现模型像水一样从底层固定边界流下来了…… ;实体模型的建立
gen zone brick size 220 1 48 p0 0 0 0 p1 220 0 0 p2 0 1 0 p3 0 0 48.4
;指定材料性质
model mohr
;底板 5米砂岩
prop dens 2610 bulk 0.92e9 shear 0.86e9 cohesion 5.54e6 fric 40.9 tens 2.6e6 range x 0 220 y 0 1 z 0 5
;8号煤 3.8米
prop dens 1310 bulk 0.28e9 shear 0.12e9 cohesion 0.63e6 fric 34.6 tens 0.5e6 range x 0 220 y 0 1 z 5 8.8
;6号底板
;石灰岩 3.8米
prop dens 2710 bulk 3.5e9 shear 5.3e9 cohesion 6.4e6 fric 35 tens 4.0e6 range x 0 220 y 0 1 z 8.8 12.6
;钙质页岩 3.2米
prop dens 2720 bulk 2.5e9 shear 2.1e9 cohesion 6.61e6 fric 41.5 tens 2.1e6 range x 0 220 y 0 1 z 12.6 15.8
;石灰岩 7.5米
prop dens 2710 bulk 3.5e9 shear 5.3e9 cohesion 6.4e6 fric 35 tens 4.0e6 range x 0 220 y 0 1 z 15.8 23.3
;钙质页岩 3.6米
prop dens 2720 bulk 2.5e9 shear 2.1e9 cohesion 6.61e6 fric 41.5 tens 2.1e6 range x 0 220 y 0 1 z 23.3 26.9
;砂岩 3米
prop dens 2610 bulk 0.92e9 shear 0.86e9 cohesion 5.54e6 fric 40.9 tens 2.6e6 range x 0 220 y 0 1 z 26.9 29.9
;炭质页岩 0.6米
prop dens 2520 bulk 0.69e9 shear 0.68e9 cohesion 2.77e6 fric 39.8 tens 1.1e6 range x 0 220 y 0 1 z 29.9 30.5
;7号煤 0.8米
prop dens 1340 bulk 0.19e9 shear 0.08e9 cohesion 0.67e6 fric 30.8 tens 0.4e6 range x 0 220 y 0 1 z 30.5 31.3
;钙质页岩 1.4米
prop dens 2720 bulk 2.5e9 shear 2.1e9 cohesion 6.61e6 fric 41.5 tens 2.1e6 range x 0 220 y 0 1 z 31.3 32.7
;石灰岩 2.5米
prop dens 2710 bulk 3.5e9 shear 5.3e9 cohesion 6.4e6 fric 35 tens 4.0e6 range x 0 220 y 0 1 z 32.7 35.2
;页岩 1.5米
prop dens 2160 bulk 0.43e9 shear 0.64e9 cohesion 2.12e6 fric 33.2 tens 0.9e6 range x 0 220 y 0 1 z 35.2 36.7
;6号煤 1.7米
prop dens 1340 bulk 0.19e9 shear 0.08e9 cohesion 0.67e6 fric 30.8 tens 0.4e6 range x 0 220 y 0 1 z 36.7 38.4
;6号顶板
;页岩 6.0 米
prop dens 2160 bulk 0.43e9 shear 0.64e9 cohesion 2.12e6 fric 33.2 tens 0.9e6 range x 0 220 y 0 1 z 38.4 44.4
;石灰岩 2.1米
prop dens 2710 bulk 3.5e9 shear 5.3e9 cohesion 6.4e6 fric 35 tens 4.0e6 range x 0 220 y 0 1 z 44.4 46.5
;页岩 1.9米
prop dens 2160 bulk 0.43e9 shear 0.64e9 cohesion 2.12e6 fric 33.2 tens 0.9e6 range x 0 220 y 0 1 z 46.5 48.4
;设置初始边界条件
fix x range x -0.1 0.1
fix x range x 219.9 220.1
fix y range y -0.1 0.1
fix y range y 0.9 1.1
fix z range z -0.1 0.1
;设置应力条件
set grav 0 0 -9.81
apply szz -3.4e6 range z 48.3 48.5 x 0 220 y 0 1
set large
model null range x 20 200 y 0 1 z 6.8 8.8
hist unbal
solve
然后到20000+步之后就开始漏了……:L