求助，关于深部煤层开采模拟的问题 [复制链接]

本人小硕一枚，目前在做千米深井煤层开采模拟研究，在模拟过程中遇到了很多问题，不知道该如何解决，请求各位大神指导一二，不胜感激！
问题一：我建模型的时候，煤层开挖2m，顶板下沉了近6m，地表下沉1.5左右。地表下沉达到我的要求了，但是顶板下沉过大~顶板调到下沉1.85m左右后达到我要求了，地表下沉最大又只有几十毫米了，我已经将各个岩层的弹性模量调到几百MPa了，难道还要继续将岩层弹性模量调小吗？
问题二：在做煤层条带开采塑性区分析的时候，1000m情况下塑性区宽度11m，参数不变采深800m塑性区变成3m，按理说应该7~9m才符合实际情况~需要调整什么参数？

命令流让我们给你指教指教

多多交流
，有意者私下交流

我也是。。。。同求第一个问题，。。。开挖2米 下沉很大。。。而且算不完锕。。。请大神指教  
new
gen zone brick p0 (0,0,0)   p1 (300,0,0)   p2 (0,100,0)   p3 (0,0,10)  size 30 10 10  group 1 ;粉砂质砂岩
gen zone brick p0 (0,0,10)  p1 (300,0,10)  p2 (0,100,10)  p3 (0,0,25)  size 30 10 15  group 2 ;泥岩
gen zone brick p0 (0,0,25)  p1 (300,0,25)  p2 (0,100,25)  p3 (0,0,30)  size 30 10 10  group 3 ;粉砂质砂岩
gen zone brick p0 (0,0,30)  p1 (300,0,30)  p2 (0,100,30)  p3 (0,0,32)  size 30 10 10  group 4 ;C5煤
gen zone brick p0 (0,0,32)  p1 (300,0,32)  p2 (0,100,32)  p3 (0,0,35)  size 30 10 10  group 5 ;泥岩
gen zone brick p0 (0,0,35)  p1 (300,0,35)  p2 (0,100,35)  p3 (0,0,43)  size 30 10 20  group 6 ;细粒砂岩
gen zone brick p0 (0,0,43)  p1 (300,0,43)  p2 (0,100,43)  p3 (0,0,54)  size 30 10 11  group 7 ;粉砂质砂岩
gen zone brick p0 (0,0,54)  p1 (300,0,54)  p2 (0,100,54)  p3 (0,0,67)  size 30 10 13  group 8 ;泥岩
gen zone brick p0 (0,0,67)  p1 (300,0,67)  p2 (0,100,67)  p3 (0,0,77)  size 30 10 20  group 9 ;砂质泥岩
gen zone brick p0 (0,0,77)  p1 (300,0,77)  p2 (0,100,77)  p3 (0,0,87)  size 30 10 20  group 10;中砂岩
gen zone brick p0 (0,0,87)  p1 (300,0,87)  p2 (0,100,87)  p3 (0,0,107) size 30 10 20  group 11;粉砂质砂岩
gen zone brick p0 (0,0,107) p1 (300,0,107) p2 (0,100,107) p3 (0,0,122) size 30 10 30  group 12;灰岩岩

model mohr;摩尔-库仑模型的材料参数
prop bulk 1.1e10  shear 7.3e9   coh 4e6   fric 32 tens 1.8e6  ran group 1 ;粉砂质砂岩
prop bulk 1e10    shear 6e9     coh 3e6   fric 36 tens 2e6    ran group 2 ;泥岩
prop bulk 1.1e10  shear 7.3e9   coh 4e6   fric 32 tens 1.8e6  ran group 3 ;粉砂质砂岩
prop bulk 2.8e9   shear 1.2e9   coh 1.6e6 fric 38 tens 1e6    ran group 4 ;C5煤
prop bulk 1e10    shear 6e9     coh 3e6   fric 36 tens 2e6    ran group 5 ;泥岩
prop bulk 1.34e10 shear 1.05e10 coh 7e6   fric 38 tens 3e6    ran group 6 ;细粒砂岩
prop bulk 1.1e10  shear 7.3e9   coh 4e6   fric 32 tens 1.8e6  ran group 7 ;粉砂质砂岩
prop bulk 1e10    shear 6e9     coh 3e6   fric 36 tens 2e6    ran group 8 ;泥岩
prop bulk 5.33e10 shear 3.2e10  coh 1.6e6 fric 32 tens 2e6    ran group 9 ;砂质泥岩
prop bulk 7.28e10 shear 5.46e10 coh 2e6   fric 36 tens 3e6    ran group 10 ;中砂岩
prop bulk 1.1e10  shear 7.3e9   coh 4e6   fric 32 tens 1.8e6  ran group 11 ;粉砂质砂岩
prop bulk 4.8e10  shear 9.6e9   coh 5e6   fric 35 tens 6e6    ran group 12 ;灰岩

ini dens 2560 ran group 1;粉砂质砂岩
ini dens 2400 ran group 2;泥岩
ini dens 2560 ran group 3;粉砂质砂岩
ini dens 1500 ran group 4;C5煤
ini dens 2400 ran group 5;泥岩
ini dens 2450 ran group 6;细粒砂岩
ini dens 2560 ran group 7;粉砂质砂岩
ini dens 2400 ran group 8;泥岩
ini dens 2410 ran group 9;砂质泥岩
ini dens 2700 ran group 10;中砂岩
ini dens 2560 ran group 11;粉砂质砂岩
ini dens 2750 ran group 12;灰岩

fix x ran x -0.1 0.1
fix x ran x 299.9 300.1
fix y ran y -0.1 0.1
fix y ran y 99.9 100.1
fix z ran z -0.1 0.1

set gravity 0 0 -9.81
app nstress -8e6 ran x 0 300 y 0 100 z 122

h unb

solve
save xgxiaofalu.sav

开挖150m的
rest xgxiaofalu.sav
set large
ini xdis = 0 ydis = 0 zdis = 0
model null ran x 50 200 y 0 100 z 30 32

hist gp zdisp 30 40 35
hist gp zdisp 50 40 35
hist gp zdisp 100 40 35
hist gp zdisp 150 40 35
hist gp zdisp 200 40 35

solve
save kw150.sav