2D中侧压力系数为2的地应力的构造问题

liuzhiqiang06

埋深200m，密度2000，侧压力系数为2，在模型顶部施加4MPa的压力模拟上部200m土体重量，网格10m*10m，应该说最下面为4.2MPa的力，可是syy下部大于4.2MPa，sxx图在最上部应该为8MPa，下部为8.4MPa，但是结果却有差别，不知道为什么，贴上命令流，请高手指点！
config
grid 10,10
model elastic
group 'User:rock' notnull
model mohr notnull group 'User:rock'
prop density=2000.0 bulk=1.33333E9 shear=8E8 cohesion=200000.0 friction=30.0 dilation=0.0 tension=0.0 notnull group 'User:rock'
apply syy=-4e6 j=10
ini syy -4.2e6 var 0,2000
ini sxx -8.4e6 var 0,4000
ini szz -8.4e6 var 0,4000
fix  x i 1
fix  x i 11
fix  y j 1
set gravity=10.0
solve

liuzhiqiang06

自己顶一下  望高手指点