情况是这样的：在千米埋深处开挖一个巷道，垂直应力与水平应力已经给出，由三种不同的岩层：砂页岩：10m、泥岩:12m、破碎岩:8m（由上到下）(巷道在泥岩中) ，各岩层相关参数也已经知道
但是目前遇到了一个问题，在建模之后，开挖之前，显示塑性区的时候出现了下面的情况，改了好久，不知道问题出在了哪里请指教，感谢（最后，祝大家天天开心，万事如意呀）
n
plot block group
gen zone radcylinder dim 2 2 2 2 p0 0 0 0 p1 10 0 0 p2 0 30.0 0 p3 0 0 5 &
ratio 1.0 1.0 1.0 1.1 size 10 30 20 20 fill group uptunnel
group upexcav range cyl end1 0 0 0 end2 0 30.0 0 rad 2 group uptunnel
group uptunnel range group upexcav not
gen zone radtunnel dim 1.4 2. 1.4 2 p0 0 0 0 p1 0 0 -7 p2 0 30.0 0 &
p3 10 0 0 ratio 1.0 1.0 1.0 1.1 size 7 30 10 20 fill
group lowtunnel range group uptunnel not group upexcav not
group lowexcav range x 0 2 z -1.4 0 group lowtunnel
gen zone brick p0 0 0 -15 p1 10 0 -15 p2 0 30 -15 p3 0 0 -7 size 10 30 10 group 3
gen zone brick p0 0 0 5 p1 10 0 5 p2 0 30 5 p3 0 0 15 size 10 30 10 group 1
gen zone reflect norm -1,0,0 ori 0,0,0
save jianmo.sav
fix z range z -15
fix x range x -10
fix x range x 10
fix y range y 0
fix y range y 30
apply szz -27e6
ini syy -23.85e6
ini sxx -23.85e6
model mohr
set gravity 0 0 -10
prop bulk 0.83e9 shear 0.29e9 coh 1.0e5 fric 30 tens 1.0e4 range z -15 -7
prop bulk 2.66e9 shear 1.23e9 coh 1.2e6 fric 35 tens 1.0e6 range z -7 5
prop bulk 2.38e9 shear 1.63e9 coh 3.0e6 fric 40 tens 2.0e6 range z 5 15
initial density 2.0e3 range z -15 -7
initial density 2.5e3 range z -7 5
initial density 2.65e3 range z 5 15
set mech ratio=1e-4
plot block state
hist unb
step 1000
plot block state
save Gravsol.sav

塑性区太多，说明土体力学性质太差。