plot block group
gen zon radcyl p0 9 0 0 p1 15.5 0 0 p2 9 100 0 p3 9 0 6.5 size 5 20 10 4 dim 6 6 6 6 rat 1 1 1 1.2 group outrockys
gen zon cshell p0 9 0 0 p1 15 0 0 p2 9 100 0 p3 9 0 6 size 1 20 10 4 dim 6 6 6 6 rat 1 1 1 1 fill group inrockys
gen zon radcyl p0 9 0 0 p1 9 0 -4.1 p2 9 100 0 p3 15.5 0 0 size 5 20 10 4 dim 3.6 6 3.6 6 rat 1 1 1 1.2 group outrockyx
gen zon cshell p0 9 0 0 p1 9 0 -3.6 p2 9 100 0 p3 15 0 0 size 1 20 10 4 dim 3.6 6 3.6 6 rat 1 1 1 1 fill group inrockyx
gen zon reflect dip 90 dd 270 orig 9 0 0
gen zon brick p0 0 0 -4.1 p1 2.5 0 -4.1 p2 0 100 -4.1 p3 0 0 0 size 3 20 5 group outrock1
gen zon brick p0 15.5 0 -4.1 p1 100 0 -4.1 p2 15.5 100 -4.1 p3 15.5 0 0 size 20 20 5 group outrock1
gen zon brick p0 0 0 0 p1 2.5 0 0 p2 0 100 0 p3 0 0 6.5 size 3 20 5 group outrock2
gen zon brick p0 15.5 0 0 p1 100 0 0 p2 15.5 100 0 p3 15.5 0 6.5 size 20 20 5 group outrock2
gen zon brick p0 0 0 6.5 p1 2.5 0 6.5 p2 0 100 6.5 p3 0 0 60 size 3 20 10 group outrock3
gen zon brick p0 15.5 0 6.5 p1 100 0 6.5 p2 15.5 100 6.5 p3 15.5 0 60 size 20 20 10 group outrock3
gen zon brick p0 2.5 0 6.5 p1 15.5 0 6.5 p2 2.5 100 6.5 p3 2.5 0 60 size 10 20 10 group outrock3
gen zon brick p0 0 0 -60 p1 2.5 0 -60 p2 0 100 -60 p3 0 0 -4.1 size 3 20 8 group outrock4
gen zon brick p0 15.5 0 -60 p1 100 0 -60 p2 15.5 100 -60 p3 15.5 0 -4.1 size 20 20 8 group outrock4
gen zon brick p0 2.5 0 -60 p1 15.5 0 -60 p2 2.5 100 -60 p3 2.5 0 -4.1 size 10 20 8 group outrock4
gen zon reflect dip 90 dd 270 orig 0 0 0
plot block group
;模拟自重应力
set gravity 0 0 -10
fix z range z -60.01,-59.99
fix x range x -100.01,-99.99
fix x range x 99.99,100.01
fix y range y -0.01,0.01
fix y range y 99.99,100.01
model mohr
ini density 2400
prop bulk=6.0e9,shear=10.0e9,fric=45,coh=1.2e6,tension=2e6
set mech ratio=1e-7
save Gsol.sav
结果在定义model mohr模型的时候，就出现Warning zone******** has illegal geometry 的提示，请高手指教。谢谢
