地下洞室开挖初始应力场计算

X.g

图为计算的初始应力场，中间应力有很大问题，不知哪位大神可以指导一下啊？
命令流如下：model mohr
ini dens 2800
prop bulk 16.67e9 shear 10e9 coh 15e6 fric 50 tension 5.5e6

ini dens 2300 range group 37 any
prop bulk 5e9 shear 2.31e9 coh 0.7e6 fric 35 range group 37 any

ini dens 2000 range group 5 any group 9 any group 20 any group 31 any group 33 any group 34 any
prop bulk 3.33e9 shear 0.714e9 coh 0.12e6 fric 29 range group 5 any group 9 any group 20 any &
group 31 any group 33 any group 34 any

apply sxx=-23.8e6 grad 0,0,3.4877e4 range x -360.5 -360.3 z 0 682.397
apply sxx=-23.8e6 grad 0,0,2.4895e4 range x 500.8 500.9 z 0 956
apply syy=-32.55e6 grad 0,0,3.1814e4 range y -266.1 -265.9 z 0 1023.139
apply syy=-32.55e6 grad 0,0,5.2905e4 range y 373.4 373.6 z 0 615.258
apply szz=-32.2e6 grad 0,0,3.93066e4 range z 0 819.2

set gravity 0 -9.8 0

hist id=1 unbal
step 10000