这个是模型的代码：
;煤层5.8m，变电所宽4.5m，直墙高2m，拱高1.5m，碹体厚0.5m，
generate zone radcylinder p0 0 0 0 p1 3.5 0 0 p2 0 3 0 p3 0 0 3 dimension 2.75 2 2.75 2 size 5 3 10 1 group 煤层
generate zone cshell p0 0 0 0 p1 2.75 0 0 p2 0 3 0 p3 0 0 2 dimension 2.25 1.5 2.25 1.5 size 1 3 10 4 group 碹体 fill group 变电所
generate zone brick p0 3.5 0 0 p1 add 20 0 0 p2 add 0 3 0 p3 add 0 0 3 size 20 3 5 group 煤层
generate zone brick p0 0 0 -2.0 p1 add 2.25 0 0 p2 add 0 3 0 p3 add 0 0 2.0 size 4 3 4 group 变电所
generate zone brick p0 2.25 0 -2.0 p1 add 0.5 0 0 p2 add 0 3 0 p3 add 0 0 2.0 size 1 3 4 group 碹体
generate zone brick p0 2.75 0 -2.0 p1 add 0.75 0 0 p2 add 0 3 0 p3 add 0 0 2.0 size 1 3 4 group 煤层
generate zone brick p0 3.5 0 -2.0 p1 add 20 0 0 p2 add 0 3 0 p3 add 0 0 2.0 size 20 3 4 group 煤层generate zone brick p0 0 0 -2.8 p1 add 2.25 0 0 p2 add 0 3 0 p3 add 0 0 0.8 size 4 3 1 group 煤层
generate zone brick p0 2.25 0 -2.8 p1 add 0.5 0 0 p2 add 0 3 0 p3 add 0 0 0.8 size 1 3 1 group 煤层
generate zone brick p0 2.75 0 -2.8 p1 add 0.75 0 0 p2 add 0 3 0 p3 add 0 0 0.8 size 1 3 1 group 煤层
generate zone brick p0 3.5 0 -2.8 p1 add 20 0 0 p2 add 0 3 0 p3 add 0 0 0.8 size 20 3 1 group 煤层
;直接顶 3.7m
generate zone brick p0 0 0 3 p1 add 3.5 0 0 p2 add 0 3 0 p3 add 0 0 3.7 size 5 3 5 group 直接顶
generate zone brick p0 3.5 0 3 p1 add 20 0 0 p2 add 0 3 0 p3 add 0 0 3.7 size 20 3 5 group 直接顶;基本顶 8.4m
generate zone brick p0 0 0 6.7 p1 add 3.5 0 0 p2 add 0 3 0 p3 add 0 0 8.4 size 5 3 10 group 基本顶
generate zone brick p0 3.5 0 6.7 p1 add 20 0 0 p2 add 0 3 0 p3 add 0 0 8.4 size 20 3 10 group 基本顶
;直接底 2.3m
generate zone brick p0 0 0 -5.1 p1 add 2.25 0 0 p2 add 0 3 0 p3 add 0 0 2.3 size 4 3 4 group 直接底
generate zone brick p0 2.25 0 -5.1 p1 add 0.5 0 0 p2 add 0 3 0 p3 add 0 0 2.3 size 1 3 4 group 直接底
generate zone brick p0 2.75 0 -5.1 p1 add 0.75 0 0 p2 add 0 3 0 p3 add 0 0 2.3 size 1 3 4 group 直接底
generate zone brick p0 3.5 0 -5.1 p1 add 20 0 0 p2 add 0 3 0 p3 add 0 0 2.3 size 20 3 4 group 直接底;基本底 1.1m
generate zone brick p0 0 0 -6.2 p1 add 2.25 0 0 p2 add 0 3 0 p3 add 0 0 1.1 size 4 3 2 group 基本底
generate zone brick p0 2.25 0 -6.2 p1 add 0.5 0 0 p2 add 0 3 0 p3 add 0 0 1.1 size 1 3 2 group 基本底
generate zone brick p0 2.75 0 -6.2 p1 add 0.75 0 0 p2 add 0 3 0 p3 add 0 0 1.1 size 1 3 2 group 基本底
generate zone brick p0 3.5 0 -6.2 p1 add 20 0 0 p2 add 0 3 0 p3 add 0 0 1.1 size 20 3 2 group 基本底
;下覆岩层 10.7m
generate zone brick p0 0 0 -16.9 p1 add 2.25 0 0 p2 add 0 3 0 p3 add 0 0 10.7 size 4 3 10 group 下覆岩层
generate zone brick p0 2.25 0 -16.9 p1 add 0.5 0 0 p2 add 0 3 0 p3 add 0 0 10.7 size 1 3 10 group 下覆岩层
generate zone brick p0 2.75 0 -16.9 p1 add 0.75 0 0 p2 add 0 3 0 p3 add 0 0 10.7 size 1 3 10 group 下覆岩层
generate zone brick p0 3.5 0 -16.9 p1 add 20 0 0 p2 add 0 3 0 p3 add 0 0 10.7 size 20 3 10 group 下覆岩层generate zone reflect origin 0 0 0 normal -1 0 0
save model.sav
model mohr
set gravity 0 0 -10
prop dens 2500 bu 3.1e9 sh 2.0e9 fric 28 coh 2.0e6 range group 基本顶
prop dens 2500 bu 3.2e9 sh 1.8e9 fric 30 coh 1.0e6 range group 直接顶
prop dens 2500 bu 1.3e9 sh 0.6e9 fric 28 coh 0.6e6 range group 煤层
prop dens 2500 bu 1.3e9 sh 0.6e9 fric 28 coh 0.6e6 range group 变电所
prop dens 2500 bu 1.3e9 sh 0.6e9 fric 28 coh 0.6e6 range group 碹体
prop dens 2500 bu 3.2e9 sh 1.8e9 fric 30 coh 1.0e6 range group 直接底
prop dens 2500 bu 3.1e9 sh 2.0e9 fric 28 coh 2.0e6 range group 基本底
prop dens 2500 bu 3.2e9 sh 1.8e9 fric 30 coh 1.0e6 range group 下覆岩层
fix x range x -23.4 -23.6
fix y range y -0.1 0.1
fix x range x 23.4 23.6
fix y range y 2.9 3.1
fix z range z -16.8 -17
apply szz -5e6 range z 15 15.2
ini szz -5.3775e6 grad 0 0 2.5e4
ini sxx -6.453e6 grad 0 0 3e4
ini syy -6.453e6 grad 0 0 3e4
;set large
hist unbal
pl cont szz ou on
solve
建好的模型图为
但计算初始地应力时一直不收敛，求大神们指点！！！不甚感激，本人菜鸟一个！
