命令流如下，先有e模型，为什么求解得：zero stiffness in grid point-1
new
plo sh
plo add axes y
plo add sur red
;模型1建立,模型X=600m，Y=150m，Z=200m
;生成矿体1（Ore1）网格，矿体倾角50°，矿体宽30m即采场宽
;为便于后面模拟采场宽度，矿体分10m(R1)-5m(R2)-5m(L2)-10m(L1)每层建模
;生成矿体（Ore-R1-10）
gen zo brick P0 150 0 0 P1 160 0 0 P2 150 100 0 P3 234 0 100 P4 160 100 0 P5 234 100 100 P6 244 0 100 P7 244 100 100 size 2 20 20 group Ore-R1-10
;生成矿体（Ore-R2-5）
gen zo brick P0 160 0 0 P1 165 0 0 P2 160 100 0 P3 244 0 100 P4 165 100 0 P5 244 100 100 P6 249 0 100 P7 249 100 100 size 1 20 20 group Ore-R2-5
;生成矿体（Ore-L2-5）
gen zo brick P0 165 0 0 P1 170 0 0 P2 165 100 0 P3 249 0 100 P4 170 100 0 P5 249 100 100 P6 254 0 100 P7 254 100 100 size 1 20 20 group Ore-L2-5
;生成矿体（Ore-L1-10）
gen zo brick P0 170 0 0 P1 180 0 0 P2 170 100 0 P3 254 0 100 P4 180 100 0 P5 254 100 100 P6 264 0 100 P7 264 100 100 size 2 20 20 group Ore-L1-5
;生成左半围岩1（LeftRock-2）网格
gen zo brick P0 0 0 0 P1 150 0 0 P2 0 100 0 P3 0 0 100 size 15 20 20 group LeftRock-1
;生成左半围岩2（LeftRock-1）网格
gen zo wedge P0 150 100 100 P1 234 100 100 P2 150 0 100 P3 150 100 0 P4 234 0 100 P5 150 0 0 size 6 20 20 group LeftRock-2
;生成右半围岩（RightRock-2）网格
gen zo wedge P0 180 0 0 P1 264 0 0 P2 180 100 0 P3 264 0 100 P4 264 100 0 P5 264 100 100 size 6 20 20 group RightRock-2
;生成右半围岩（RightRock-1）网格
gen zo brick P0 264 0 0 P1 414 0 0 P2 264 100 0 P3 264 0 100 size 15 20 20 group RightRock-1
;生成顶板围岩（TopRock1）网格
gen zo brick P0 0 0 100 P1 150 0 100 P2 0 100 100 P3 0 0 200 size 15 20 10 group TopRock1
;生成顶板围岩（TopRock2）网格
gen zo brick P0 150 0 100 P1 234 0 100 P2 150 100 100 P3 150 0 200 size 6 20 10 group TopRock2
;生成顶板围岩（TopRock3）网格
gen zo brick P0 234 0 100 P1 264 0 100 P2 234 100 100 P3 234 0 200 size 6 20 10 group TopRock3
;生成顶板围岩（TopRock4）网格
gen zo brick P0 264 0 100 P1 414 0 100 P2 264 100 100 P3 264 0 200 size 15 20 10 group TopRock4
;把在group中部分染成红色、黑色、蓝色和橙色
plot add block group red yellow blue orange
;连接毗邻的网格
attach face range z -0.1 200.1
;边界条件
;定义模型材料，材料为砂岩，体积模量为27GPa 剪切模量7.0GPa
mo e
;model mohr
;prop bulk 2.7e10 shear 7.0e9 fric 35
;约束边界
fix z range z -0.1 0.1
fix x range x -0.1 0.1
fix x range x 413.9 414.1
fix y range y -0.1 0.1
fix y range y 99.9 100.1
;为进行开挖，进行模型初始化密度2700kg/m3
initial density=2700
;加载重力加速度为10
set gravity 0 0 -10
;施加初始应力
;initial sxx -1.215e7 grad 0,0,2.7e4
;initial syy -1.215e7 grad 0,0,2.7e4
;initial szz -1.215e7 grad 0,0,2.7e4
;施加上部压力2700*250*10=6.75e6
;apply szz -6.75e6 range z 199.9 200.1
;绘制任意剖面位移等值线图-XX面（命令详解P14）
plot create TrenchXX
plot set plane dip=90 dd=0 origin=300 75 100
plot add boundary behind
plot add contour disp
plot add axes black
pl sh
pl grid red vel blue
save model-1.sav
set large
step 100000
save model-1.sav
