gen zone brick p0 0 0 0 p1 30 0 0 p2 0 3 0 p3 0 0 12 size 60 6 24 group 1 ;;底板（12m泥岩）
gen zone brick p0 0 0 12 p1 30 0 12 p2 0 3 12 p3 0 0 14 size 60 6 4 group 2 ;2m3-1煤
gen zone brick p0 0 0 14 p1 30 0 14 p2 0 3 14 p3 0 0 26 size 60 6 24 group 3 ;顶板1（12m中砂岩）
pl bl gr
model mohr ;;赋予模型莫尔-库伦本构模型准则
prop bulk 1.97e10 she 1.14e10 fric 20.9 coh 4e6 ten 4e6 ran group 1
prop bulk 3.1e6 she 2.8e6 fric 23 coh 0.45e6 ten 0.32e6 ran group 2
prop bulk 1.65e10 she 1.25e10 fric 30.5 coh 8e6 ten 25.7e6 ran group 3
;密度
ini dens 2400 ran group 1
ini dens 1290 ran group 2
ini dens 2600 ran group 3
; ----------------------- 施加边界条件 ----------------
attach face range z 11.9 12.1
attach face range z 13.9 14.1
set grav 0 0 -9.8
fix x range x -0.1 0.1
fix x range x 29.9 30.1
fix y range y -0.1 0.1
fix y range y 2.9 3.1
fix z range z -0.1 0.1
app nstress -2.72e9 ran x 0 30 y 0 3 z 26 ;;在计算模型的顶部(Z=60)、X=0~30、Y=0~2的范围内，设定均布载荷8.0MPa，以“-”表示压应力
app nstress 1.17e9 ran x 0 y 0 3 z 0 26 ;;在计算模型的左部(X=0)、z=0~42、Y=0~2的范围内，设定向右均布载荷，以“-”表示压应力
app nstress -1.17e8 ran x 30 y 0 3 z 0 26 ;;在计算模型的左部(X=0)、z=0~42、Y=0~2的范围内，设定向左均布载荷，以“-”表示压应力
hist unbal ;;设定模型的最大不平衡力
set large
solve
model null ran x 12.5 17.5 y 0 3 z 12 14.5
initial xdis = 0 ydis = 0 zdis = 0
set mechanical force 50 ;设置最大不平衡力为50N，一旦小于此值则终止计算
solve
plot show