模拟结果很奇怪

yebochuan

我用FLAC3D模拟煤层开挖，煤层厚度从2m增加到了12m，可是采厚覆岩破坏高度均不变，不知道哪里出了问题了，请教各位老师！
代码如下：
n
gen zone brick   ratio 1,1,1   size 50,36,10  & p0 0,0,0         p1 500,0,0   p2 0,360,0      p3 0,0,20   gr 铝土质泥岩
gen zone brick   ratio 1,1,1   size 50,36,1  & p0 0,0,20       p1 500,0,20  p2 0,360,20     p3 0,0,22   gr 4煤
gen zone brick   ratio 1,1,1   size 50,36,10  & p0 0,0,22        p1 500,0,22  p2 0,360,22     p3 0,0,42   gr 砂质泥岩
gen zone brick   ratio 1,1,1   size 50,36,16  & p0 0,0,42        p1 500,0,42  p2 0,360,42     p3 0,0,74   gr 细粒砂岩
gen zone brick   ratio 1,1,1   size 50,36,10  & p0 0,0,74        p1 500,0,74  p2 0,360,74     p3 0,0,94   gr 砂质泥岩
gen zone brick   ratio 1,1,1   size 50,36,4  & p0 0,0,94        p1 500,0,94  p2 0,360,94     p3 0,0,102   gr 粉砂岩
gen zone brick   ratio 1,1,1   size 50,36,10  & p0 0,0,102        p1 500,0,102  p2 0,360,102     p3 0,0,122  gr 砂质泥岩
gen zone brick   ratio 1,1,1   size 50,36,20  & p0 0,0,122        p1 500,0,122  p2 0,360,122     p3 0,0,162   gr  细粒砂岩
gen zone brick   ratio 1,1,1   size 50,36,24  & p0 0,0,162        p1 500,0,162  p2 0,360,162     p3 0,0,210   gr  砂质泥岩
gen zone brick   ratio 1,1,1   size 50,36,16  & p0 0,0,210       p1 500,0,210  p2 0,360,210    p3 0,0,242   gr  粗砾岩
gen zone brick   ratio 1,1,1   size 50,36,18  & p0 0,0,242       p1 500,0,242  p2 0,360,242    p3 0,0,278   gr  细石英砂岩

pl bl g
sav 2mnormmodel.sav



def para
;铝土质泥岩
p1=0.09
bulk1=3435e6
shear1=3876e6
coh1=3.9e6
fric1=37.5
tens1=1.9e6
dens1=2.55e3
;4煤
p2=0.29
bulk2=2381e6
shear2=1163e6
coh2=1.3e6
fric2=32.9
tens2=0.6e6
dens2=1.35e3
;砂质泥岩
p3=0.22
bulk3=4167e6
shear3=2869e6
coh3=3.5e6
fric3=37
tens3=1.8e6
dens3=2.42e3
;细粒砂岩
p4=0.13
bulk4=9302e6
shear4=9137e6
coh4=4.2e6
fric4=28
tens4=2.1e6
dens4=2.64e3
;粉砂岩
p5=0.14
bulk5=3372e6
shear5=3816e6
coh5=4.5e6
fric5=20.4
tens5=1.7e6
dens5=2.53e3
;粗砾岩
p6=0.19
bulk6=2135e6
shear6=1668e6
coh6=1e6
fric6=34
tens6=0.6e6
dens6=2.63e3
;细石英砂岩
p7=0.17
bulk7=7222e6
shear7=6111e6
coh7=5e6
fric7=35
tens7=1.3e6
dens7=2.41e3

end
para

mo m
;摩尔库伦破坏准则

set grav 0 0 -10
fix x range x -0.1 0.1
fix x range x 499.9 500.1
fix y range y -0.1 0.1
fix y range y 359.9 360.1
fix x y z range z -0.1 0.1
app nstress -1.24e6 ran x 0 500 y 0 360 z 278


prop dens=dens1 bulk=bulk1 sh=shear1 coh=coh1 fric=fric1 ten=tens1 range gr 铝土质泥岩
prop dens=dens2 bulk=bulk2 sh=shear2 coh=coh2 fric=fric2 ten=tens2 range gr 4煤
prop dens=dens3 bulk=bulk3 sh=shear3 coh=coh3 fric=fric3 ten=tens3 range gr 砂质泥岩
prop dens=dens4 bulk=bulk4 sh=shear4 coh=coh4 fric=fric4 ten=tens4 range gr 细粒砂岩
prop dens=dens5 bulk=bulk5 sh=shear5 coh=coh5 fric=fric5 ten=tens5 range gr 粉砂岩
prop dens=dens6 bulk=bulk6 sh=shear6 coh=coh6 fric=fric6 ten=tens6 range gr 粗砾岩
prop dens=dens7 bulk=bulk7 sh=shear7 coh=coh7 fric=fric7 ten=tens7 range gr 细石英砂岩

his unb
pl hi 1
so
sav 2mnormnaure.sav


ini zdis=0 ydis=0 xdis=0
ini sta 0;对四面体的塑性标志设置为0 没塑性破坏
set large ;大变形

def ji
loop n (0,29)
jilu=n+1
jilu1=100+n*10
jilu2=jilu1+10
name='dig' + string(jilu) + '.sav'
command
model null range gr 4煤 x jilu1 jilu2 y 80 280 z 20 22
step 2000

sav name
endcommand
end_loop
end
ji


求助！！