初学者多多指教，我的思路是做一个大的山脉网格模型，然后在下面构建采空区模型，然后把采空区模型设为空，从而得到竖直和水平位移等值线。却不知道要注意些什么，请大神指点
现在用下面代码来实验 （里面开挖部分语法不懂，求助）
new
;研究区矩形=====================
generate zone brick size 44 44 12 p0 = (47.0919,15.80734,0)&
p1 = (2677.6602,15.80734,0) p2 = (47.0919,2646.37557,0)&
p3 = (47.0919,15.80734,120)
;read D:\flacbox\mesh.sav
model elastic
prop shear 3e8 bulk 5e8
def mountain
gp = gp_head
loop while gp #null
zz = sprt(gp_xpos(gp)^2 + gp_ypos(gp)^2)
;傅里叶条件求和
dz = 0.06 * sin(0.2*zz + 100.0)
;伪随机表面
dz = dz + 0.06 * sin(0.22^zz - 20.3)
dz = dz - 0.04 * sin(0.33*zz + 33.3);topoglgy.
gp_zpos(gp) = 0.5*gp_zpos(gp)*(1.0+dz)
gp = gp_next(gp)
end_loop
end
mountain
;边界条件=======================
fix x range x 47.09 47.10
fix x range x 2677.65 2677.67
fix y range y 15.807 15.808
fix y range y 2646.375 2646.376
fix z range z -0.1 0.1
set grav 0,0,-10
ini density = 2000
ini szz = -2.0e6 (grad 0,0,2.0e4)
ini sxx = -4.0e6 (grad 0,0,4.0e4) syy = -4.0e6 (grad 0,0,4.0e4)
solve
;挖三块===============================
generate zone brick p0 1392.6947, 496.181, 0 &
p1 1783.725, 192.14335, 0 p2 1619.8115, 1403.15808, 0 &
p3 1392.6947, 496.181, 20 p4 2186.8456, 1317.11615, 0 &
p5 1619.8115, 1403.15808, 20 p6 1783.725, 192.14335, 20 &
p7 2186.8456, 1317.11615, 20 size 8,8,5 group 1
; 2 3 7 1===
gen zone brick p0 933.721, 1113.24163, 0 &
p1 1392.6947, 496.181, 0 p2 1136.79, 1492.21374, 0 &
p3 933.721, 1113.24163, 20 p4 1514.8379, 983.95258, 0 &
p5 1136.79, 1492.21374, 20 p6 1392.6947, 496.181, 20 &
p7 1514.8379, 983.95258, 20 size 8,8,5 group 2
;1 7 6 0====
gen zone brick p0 1136.79, 1492.21374, 0 &
p1 1514.8379, 983.95258, 0 p2 1129.8124, 1663.21506, 0 &
p3 1136.79, 1492.21374, 20 p4 1619.8115, 1403.15808, 0 &
p5 1129.8124, 1663.21506, 20 p6 1514.8379, 983.95258, 20 &
p7 1619.8115, 1403.15808, 20 size 8,8,5 group 3
model null range group 1
model null range group 2
model null range group 3
save d:\mod ;!!!!!!!!!!!!!!!!!!!!此步得到的模型看不出来已经挖掉group的区域 为什么呢????
set large
initial xdisplacement = 0 ydisplacement = 0 zdisplacement = 0
step = 1000 ;计算步数
save d:\1000.sav
;输出位移等值线图======================
plot creat t
;plot add contour disp ;添加 当前视图等值线 位移量
plot add contour zdisplacement
plot show
