new;
title;;
矿房
gen zone brick p0 0 0 0 p1 60 0 0 p2 0 20 0 p3 0 0 5 p4 60 20 0 p5 0 20 5 p6 60 0 35 p7 60 20 35 size 30 10 25 group 1;;页岩
gen zone brick p0 0 0 5 p1 60 0 35 p2 0 20 5 p3 0 0 8 p4 60 20 35 p5 0 20 8 p6 60 0 38 p7 60 20 38 size 30 10 25 group 2;;砂岩
gen zone brick p0 0 0 8 p1 60 0 38 p2 0 20 8 p3 0 0 13 p4 60 20 38 p5 0 20 13 p6 60 0 43 p7 60 20 43 size 30 10 25 group 3;;矿体
gen zone brick p0 0 0 13 p1 60 0 43 p2 0 20 13 p3 0 0 50 p4 60 20 43 p5 0 20 50 p6 60 0 50 p7 60 20 50 size 30 10 25 group 4;;白云岩
model mohr ;;设置模型本构关系
prop bulk 14.35e9 shear 3.4e9 coh 14.09e6 fric 42.82 tens 2.67e6 range group 1
prop bulk 10.97e9 shear 7.22e9 coh 29.78e6 fric 42.56 tens 3.02e6 range group 2
prop bulk 19.47e9 shear 11.68e9 coh 36.67e6 fric 41.94 tens 4.46e6 range group 3
prop bulk 21.61e9 shear 11.74e9 coh 37.49e6 fric 32.27 tens 5.61e6 range group 4
;;设置边界条件
fix x range x -.1 .1
fix x range x 59.9 60.1
fix y range y -.1 .1
fix y range y 19.9 20.1
fix z range z -.1 .1
;;设置初试应力状态
ini dens 2670 range group 1
ini dens 2680 range group 2
ini dens 3220 range group 3
ini dens 2810 range group 4
set grav 0 0 -10
app nstress -10802000 ran x 0 60 y 0 20 z 50
hist unbal
solve
