new
;矿体
gen zone brick p0 0 0 0 p1 80 0 0 p2 0 30 0 p3 0 45 45 p4 80 30 0 p5 0 75 45 p6 80 45 45 p7 80 75 45 size 20 20 2 0 group KT
gen zone brick p0 80 0 0 p1 320 0 0 p2 80 30 0 p3 80 45 45 p4 320 30 0 p5 80 75 45 p6 320 45 45 p7 320 75 45 size 18 20 20 ratio 1.1 1 1 group KT
gen zone brick p0 -240 0 0 p1 0 0 0 p2 -240 30 0 p3 -240 45 45 p4 0 30 0 p5 -240 75 45 p6 0 45 45 p7 0 75 45 size 18 20 20 ratio 0.91 1 1 group KT
gen zone brick p0 0 45 45 p1 80 45 45 p2 0 75 45 p3 0 180 180 p4 80 75 45 p5 0 210 180 p6 80 180 180 p7 80 210 180 size 20 20 18 ratio 1 1 1.1 group KT
gen zone brick p0 0 -135 -135 p1 80 -135 -135 p2 0 -105 -135 p3 0 0 0 p4 80 -105 -135 p5 0 30 0 p6 80 0 0 p7 80 30 0 size 20 20 18 ratio 1 1 0.91 group KT
gen zone brick p0 -240 -135 -135 p1 0 -135 -135 p2 -240 -105 -135 p3 -240 0 0 size 18 20 18 ratio 0.91 1 0.91 group KT
gen zone brick p0 80 -135 -135 p1 320 -135 -135 p2 80 -105 -135 p3 80 0 0 size 18 20 18 ratio 1.1 1 0.91 group KT
gen zone brick p0 -240 45 45 p1 0 45 45 p2 -240 75 45 p3 -240 180 180 size 18 20 18 ratio 0.91 1 1.1 group KT
gen zone brick p0 80 45 45 p1 320 45 45 p2 80 75 45 p3 80 180 180 size 18 20 18 ratio 1.1 1 1.1 group KT
;下盘
gen zone brick p0 0 -105 -135 p1 80 -105 -135 p2 0 210 -135 p3 0 30 0 p4 80 210 -135 p5 0 210 0 p6 80 30 0 p7 80 210 0 size 20 18 18 ratio 1 1.1 0.91 &
group XP
gen zone brick p0 0 30 0 p1 80 30 0 p2 0 210 0 p3 0 75 45 p4 80 210 0 p5 0 210 45 p6 80 75 45 p7 80 210 45 size 20 18 20 ratio 1 1.1 1 group XP
gen zone w p0 0 210 45 p1 0 75 45 p2 80 210 45 p3 0 210 180 size 18 20 18 ratio 0.91 1 1.1 group XP
gen zone brick p0 -240 -105 -135 p1 0 -105 -135 p2 -240 210 -135 p3 -240 30 0 p4 0 210 -135 p5 -240 210 0 p6 0 30 0 p7 0 210 0 size 18 18 18 ratio 0.91 1.1 0.91 &
group XP
gen zone brick p0 80 -105 -135 p1 320 -105 -135 p2 80 210 -135 p3 80 30 0 p4 320 210 -135 p5 80 210 0 p6 320 30 0 p7 320 210 0 size 18 18 18 ratio 1.1 1.1 0.91 &
group XP
gen zone brick p0 -240 30 0 p1 0 30 0 p2 -240 210 0 p3 -240 75 45 p4 0 210 0 p5 -240 210 45 p6 0 75 45 p7 0 210 45 size 18 18 20 ratio 0.91 1.1 1 &
group XP
gen zone brick p0 80 30 0 p1 320 30 0 p2 80 210 0 p3 320 75 45 p4 320 210 0 p5 80 210 45 p6 320 75 45 p7 320 210 45 size 18 18 20 ratio 1.1 1.1 1 &
group XP
gen zone w p0 -240 210 45 p1 -240 75 45 p2 0 210 45 p3 -240 210 180 p4 0 75 45 p5 0 210 180 size 18 18 18 ratio 0.91 0.91 1.1 group XP
gen zone w p0 80 210 45 p1 80 75 45 p2 320 210 45 p3 80 210 180 p4 320 75 45 p5 320 210 180 size 18 18 18 ratio 0.91 1.1 1.1 group XP
;上盘
gen zone brick p0 0 -135 0 p1 80 -135 0 p2 0 0 0 p3 0 -135 45 p4 80 0 0 p5 0 45 45 p6 80 -135 45 p7 80 45 45 size 20 18 20 ratio 1 0.91 1 group SP
gen zone w p0 0 -135 0 p1 0 0 0 p2 80 -135 0 p3 0 -135 -135 size 18 20 18 ratio 0.91 1 1.1 group SP
gen zone brick p0 -240 -135 0 p1 0 -135 0 p2 -240 0 0 p3 -240 -135 45 p4 0 0 0 p5 -240 45 45 p6 0 -135 45 p7 0 45 45 size 18 18 20 ratio 0.91 0.91 1 group SP
gen zone brick p0 80 -135 0 p1 320 -135 0 p2 80 0 0 p3 80 -135 45 p4 320 0 0 p5 80 45 45 p6 320 -135 45 p7 320 45 45 size 18 18 20 ratio 1.1 0.91 1 group SP
gen zone w p0 -240 -135 0 p1 -240 0 0 p2 0 -135 0 p3 -240 -135 -135 size 18 18 18 ratio 0.91 0.91 1.1 group SP
gen zone w p0 80 -135 0 p1 80 0 0 p2 320 -135 0 p3 80 -135 -135 size 18 18 18 ratio 0.91 1.1 1.1 group SP
gen zone brick p0 0 -135 45 p1 80 -135 45 p2 0 45 45 p3 0 -135 180 p4 80 45 45 p5 0 180 180 p6 80 -135 180 p7 80 180 180 size 20 18 18 ratio 1 0.91 1.1 group SP
gen zone brick p0 80 -135 45 p1 320 -135 45 p2 80 45 45 p3 80 -135 180 p4 320 45 45 p5 80 180 180 p6 320 -135 180 p7 320 180 180 size 18 18 18 ratio 1.1 0.91 1.1 group SP
gen zone brick p0 -240 -135 45 p1 0 -135 45 p2 -240 45 45 p3 -240 -135 180 p4 0 45 45 p5 -240 180 180 p6 0 -135 180 p7 0 180 180 size 18 18 18 ratio 0.91 0.91 1.1 group SP
plot add surface yellow
plot add axes black
plot show
;
model mohr
prop bulk 2.24e11 shear 1.68e11 coh 5.72e6 fri 30.6 ten 3.18e6 range group SP
prop bulk 2.42e11 shear 1.89e11 coh 6.43e6 fri 32.6 ten 3.72e6 range group KT
prop bulk 3.29e11 shear 2.07e11 coh 10.7e6 fri 36.9 ten 4.31e6 range group XP
ini dens 2706 range group SP
ini dens 2710 range group KT
ini dens 2635 range group XP
;应力条件
apply sxx -16.04e6 grad 0 0 1.95e4 ran x -240.1 -239.9
apply sxx -16.04e6 grad 0 0 1.95e4 ran x 319.9 320.1
apply syy -30.16e6 grad 0 0 3.62e4 ran y -135.1 -134.9
apply syy -30.16e6 grad 0 0 3.62e4 ran y 209.9 210.1
apply szz -16.59e6 grad 0 0 1.86e4 ran z -135 180
set grav 0 0 -9.81
solve elas
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
我查看了好多次 都找不出原因
