高手们，求助深部开挖位移分析问题

haozhenggong

我在模拟一个深部开挖，但是在开挖后得到的y向正好相反的，两帮应该是向采空区移动。而图显示的是向外侧移动，这到底是哪里出错了呢。这是平衡代码：
m m
pro bulk 18.7e9   shear 11.2e9     cohesion 8.7e6    fric 34  ten 3.0e6   range group 粉砂岩
pro bulk 21.0e9   shear 11.0e9     cohesion 8.2e6    fric 35  ten 2.6e6   range group 细砂岩
pro bulk 21.6e9   shear 11.2e9     cohesion 8.3e6    fric 34  ten 3.0e6   range group 中砂岩
pro bulk 3.4e9    shear 1.5e9      cohesion 1.5e6    fric 28  ten 1.6e6   range group 煤层
pro bulk 3.4e9    shear 1.5e9      cohesion 1.5e6    fric 28  ten 1.6e6   range group 煤层工作面
pro bulk 18.3e9   shear 8.4e9      cohesion 8.5e6    fric 35  ten 2.2e6   range group 泥岩
pro bulk 26.7e9   shear 8.9e9      cohesion 8.4e6    fric 38  ten 2.5e6   range group 砂质泥岩
pro bulk 4.5e9    shear 1.0e9      cohesion 0.9e6    fric 40 ten 0.2e6    range group 陷落柱
ini density 2400 range group 粉砂岩
ini density 2300 range group 细砂岩
ini density 2300 range group 中砂岩
ini density 1400 range group 煤层
ini density 1400 range group 煤层工作面
ini density 2200 range group 泥岩
ini density 2100 range group 砂质泥岩
ini density 1200 range group 陷落柱
;确定边界条件
===================================================
plot set rot 0 0 0
plot set center auto
plot con szz outline on
plot show
fix x y z               range z -0.1    0.1
fix x                   range x -0.1    0.1
fix x                   range x 1399.9  1400.1
fix y                   range y -0.1    0.1
fix y                   range y 599.9   600.1
set grav 0 0 -9.8
ini szz -1.5e7 grad 0,0,2.5e4
ini syy -0.75e7 grad 0,0,1.25e4
ini sxx -0.75e7 grad 0,0,1.25e4
apply szz -0.75e7 range z 299.9 300.1
set large
hist unbal    ;
solve
ini xdisp=0 ydisp=0 zdisp=0
ini xvel=0  yvel=0  zvel=0

漫步听海

重力加速度的方向应该是0 -9.8 0

haozhenggong

这个图是我吧平面的转了90度， 重力加的是对的 2# 漫步听海

yangxin546

你算自重场为什么开大变形？
你的开挖步骤呢？

tanhao365

不会是上下那两个面没加约束吧（Y向）？

haozhenggong

5# tanhao365
加约束了啊