隧道模型，请各位老师看看是什么问题

aei001

做了一个隧道的模型
在计算初始应力的时候就问题百出了
希望各位老师能看看
谢谢

模型用ansys建的
用mesh200单元画了平面网格
之后拉伸成体
用的solid45单元

之后导到flac里面的
计算初始应力命令流如下

mo mohr
pro bulk 1.36e10 shear 0.703e10 fric 45 coh 1.3e6 ten 1.1e6 range z 0 40
fix y range y -0.1 0.1
fix x range x -40.1 -39.9
fix x range x 39.9 40.1
fix z range z -0.1 0.1
fix z range z 39.9 40.1
ini dens 2500
set grav 0 -10 0
ini sxx -9.7e5 grad 0 9.7e3 0
ini szz -9.7e5 grad 0 9.7e3 0
ini syy -25.e5 grad 0 2.5e4 0
apply syy -5e5 range y 79.9 80.1
hist unbal
plot hist 1
solve
save 531-ini-1645-step1.sav
ret

加了y方向的ini力的话，云图就像1里面那样

wylong1979

ini sxx -9.7e5 grad 0 9.7e3 0
!W d2h8?&K8Bini szz -9.7e5 grad 0 9.7e3 0|Simwe.com|仿真|设计|有限元|虚拟仪器&Q$@,~2e"W.T"k
ini syy -25.e5 grad 0 2.5e4 0|Simwe.com|仿真|设计|有限元|虚拟仪器5f7n)q.v,`
apply syy -5e5 range y 79.9 80.1仿真分析,


应力梯度是怎么得来的？y方向的应该是2.5e5