new
;-------------------------
;露天採礦場邊坡曲率穩定性影響
;-----------------------------
;創建採礦場
gen zone radcyl &
P0 0 0 25 P1 add 80 0 0 p2 add 0 0 -25 p3 add 0 80 0 &
dim 24.5 24.5 12 12 ratio 1 1 1 1.1
gen zone radcyl &
p0 0 0 0 p1 add 80 0 0 p2 add 0 0 -15 p3 add 0 80 0 &
dim 12 12 12 12 ratio 1 1 1 1.1 fill size 5 7 10 10
gen zone brick &
p0 0 -40 -15 p1 add 12 0 0 p2 add 0 40 0 p3 add 0 0 15 &
size 5 12 7
gen zone brick &
p0 12 -40 -15 p1 add 68 0 0 p2 add 0 40 0 p3 0 0 15 &
size 10 12 7 ratio 1.1 1 1
gen zone brick &
p0 12 -40 0 p1 add 68 0 0 p2 add 0 40 0 p3 add 12.5 0 25 &
p4 add 68 40 0 p5 add 12.5 40 25 p6 add 68 0 25 p7 &
add 68 40 25 size 10 12 10 rat 1.1 1 1
;設置 Mohr Coulomb模型及材量參數
mod mo
pro bulk 2e8 she 1e8 fric 45 coh 1e5 ten 1e5
;邊界條件
;重力加速度
set gravity 0 0 -9.81
fix x range x -.1 .1
fix x range x 79.9 80.1
fix y range y -40.1 -39.9
fix y range y 79.9 80.1
fix x range z -15.1 -14.9
fix y range z -15.1 -14.9
fix z range z -15.1 -14.9
;初始化重應力
ini syy -6.25e5 grad 0 0 2.5e4
ini sxx -6.25e5 grad 0 0 2.5e4
ini szz -6.25e5 grad 0 0 2.5e4
;監測
hist unbal
hist gp xdisp 24.5 -40 25
hist gp zdisp 24.5 -40 25
hist gp xdisp 24.5 0 25
hist gp ydisp 24.5 0 25
hist gp zdisp 24.5 0 25
hist gp ydisp 0 24.5 25
hist gp zdisp 0 24.5 25
;自動安全因子解算
plo con zdisp
solve
plo blo state
ret