纠结。。。谁能帮看下，一个FLAC3D沉降题 [复制链接]

我刚学FLAC3D，在做一个沉降的题。可是后面输出来的displacement里有些点x方向也有位移，实在找不出来哪里出错了。麻烦大家帮我看一下
基础做均布荷载算宽5m，网格51m，上面可压缩性土E=4MPa；下面岩石，E=100GPa。都是水平的。可压缩层2.5m，下面我定的10m。密度没给，我上网查的平均值
我是这么写的：
new
generate zone brick size 51 1 10 &
dimension 51 1 10 group 1
generate zone brick size 51 1 3 &
p0 0 0 10 &
dimension 51 1 2.5
model elastic
property bulk 4e6 range group 1
property bulk 1e11 range group 1 not
set grav 0 0 -9.8
initial dens 1800 range group 1
initial dens 2650 range group not
fix x ran x 0
fix x ran x 51
fix y ran y 0
fix y ran y 1
fix z ran z 0
apply szz -3e4 range z=12.4 12.6 x=23 28 y=0 1
hist id=1 gp zdis 23 0.5 12.5
hist id=2 gp zdis 25 0.5 12.5
hist id=3 gp zdis 28 0.5 12.5
solve
save 0.5.sav
不知道哪里有问题哇

楼主把你计算出来的图贴出来大家可以帮着看看，，从命令流上看你像没有什么问题

终于有人理我了。。。。。谢谢
模型图

第一个监控点的位移图

log 文件
Flac3D>print gp dis range id 1any id 2 any id 3 any
Flac3D>print gp dis range id 1 any id 2 any id 3 any
Gridpoint Displacement ...
   id        X-Dis         Y-Dis         Z-Dis
------  ------------- ------------- -------------
      1 (  0.0000e+000,  0.0000e+000,  0.0000e+000)
      2 (  5.0651e-007,  0.0000e+000,  0.0000e+000)
      3 (  0.0000e+000,  0.0000e+000,  0.0000e+000)
Flac3D>set log off
这个文件里z方向没有位移，反而有个x位移，不能理解

initial dens 2650 range group not 应该为initial dens 2650 range group 1 not
其他的没看出问题来，是不是应该把边界条件再看看