求教初始化渐变地应力问题出在哪儿

forsincerity

这是我编的一段初始化地应力的程序代码
;------------------------------加载及边界条件------------------------------------------
;设置边界条件
fix x range x -10.1 -9.9
fix x range x 9.9 10.1
fix z range z -10.1 -9.9
fix y range y -0.1 0.1
fix y range y 2.9 3.1
;设置初始应力状态(这基本不用解算就平衡了)
set grav 0 0 -9.8
ini density 2500
ini szz=-5.1891e6 grad 0,0,24500 range z -10 10
ini sxx=-2.71705e6 grad 0,0,12250 range z -10 10
ini syy =-2.71705e6 grad 0,0,12250 range z -10 10
apply szz=-4.9441e6 range z 9.9 10.1
apply szz =-5.4341e6 range z -10.1 -9.9
apply sxx=-2.71705e6 grad 0,0,12250 range z -10 10
apply sxx=-2.71705e6 grad 0,0,12250 range z -10 10
apply syy=-2.71705e6 grad 0,0,12250 range z -10 10
apply syy=-2.71705e6 grad 0,0,12250 range z -10 10
我做了一个10*3*10的模型，想让模型考虑重力的渐变，使模型上表面的地应力为4.9441MPa
下表面为5.4341MPa，水平应力为垂直应力的一半，但是在模型建好后，水平应力和考虑的不一样，这模型应力初始化问题出在哪儿了呢？
请大家帮忙！

songli269

你的Z向取的是20m，还是10m啊？你再仔细看看。

forsincerity

Z取得是20米，正负各10米，巷道建在之间了

叶飒2005

俺是新手！