思路：1.用surf2flac的代码，转换出地表，solve一下
2.在山地模型底下粘上一层标准的20米厚的矩形六面体网格
3.把那20米厚的网格中间挖出一块（煤矿采空区）
4.设定计算歩数，计算并显示Z轴位移等值线
模型：
mesh.rar (98 K) 下载次数:49 代码：
代码
Q.txt (2 K) 下载次数:38 new
restore d:\Box1\mesh group 1 ;读取转换后模型
generate zone brick size 16 16 10 & ;生成“底”
p0 = (10,10,-20) &
p1 = (1610,10,-20) &
p2 = (10,1610,-20) &
p3 = (10,10,0) group 2
attach face range z = -0.01,0.01 ;融合
;模型及参数============================
model mohr
property bulk = 2.08e8 shear = 9.6e7 friction = 25
property cohesion = 1e10 tension = 1e10
;边界条件=======================
fix x range x 9.99 10.01
fix x range x 1609.99 1610.01
fix y range y 9.99 10.01
fix y range y 1609.99 1610.01
fix z range z -19.99 20.01
set grav 0,0,-10
ini density = 1860
solve
save d:\Box1\z1.sav
;开挖===========================
model null range x=410,710 y=1010,1410 z=-20,0
model null range x=710,1310 y=410,1310 z=-20,0
save z2.sav
initial xdisplacement = 0 ydisplacement = 0 zdisplacement = 0
property cohesion = 1e5 tension = 1e4
;定义剩余煤层
model mohr range z = -20,0
property bulk = 3.846e8 shear = 7.258e8 friction = 32
property cohesion = 5e5 tension = 2e5
ini density = 2550
set large
step = 500 ;计算步数
save z3.sav
;显示开挖位移变化=====================
plot create DispCont
plot add contour zdisp; 添加 等值线 z位移量
plot add axes
plot show
得到效果图：
问题大大的啊！！！！！》》》
怎么山顶的位移量最大，而挖出的采空区上的地表却没有碗型沉陷？是不是我模型参数的问题啊！求大神们帮忙解答！~
