各位大侠, 我想做个钻孔分级扩孔开挖的数值模拟，自己写了一小段程序，先求初始应力平衡，但是结果明显有问题，由于是初学者，不知道问题出在哪里，希望给位大侠有空给查看一下，希望不吝赐教,谢谢！
程序如下：
new
gen zon radcyl p0 0 ,0,0 p1 8,0,0 p2 0,2,0 p3 0,0,8 dim 1 1 1 1 size 1 15 10 15 ratio 1 1 1 1.15 1 group 4;
gen zon cs p0 0,0,0 p1 1,0,0 p2 0,2,0 p3 0,0,1 dim 0.7 0.7 0.7 0.7 size 2,5,10,2 group 3 ；第三级开挖区域
gen zon cs p0 0,0,0 p1 0.7,0,0 p2 0,2,0 p3 0,0,0.7 dim 0.4 0.4 0.4 0.4 size 2,5,10,2 group 2 ；第二级开挖区域
gen zon c p0 0,0,0 p1 0.4,0,0 p2 0,2,0 p3 0,0,0.4 size 2,5,10 group 1;第一级开挖区域
gen zone reflect orig 0 0 0 norm 0 0 1
pl bl gr
model mohr
prop bulk 19e6 shear 1e7 coh 10e3 fric 18 ten 0
ini dens 1790
fix x range x -0.1 0.1;
fix x range x 7.9 8.1;
fix y range y -0.1 0.1;
fix y range y 1.9 2.1;
fix z range z -8.1 -7.9;
set grav 0 0 -10
history n=50;
hist unbal;
pl hist 1
pl add con szz ou on
solve elas
save yingli.sav
理论上z方向应力应该是成层分布的，但是我计算的结果明显有误，请大侠们指教
