gen zone radbrick size 3 5 3 7 p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 ratio 1 1 1 1.2
plo sur yellow
plo add axes red
gen zone reflect dip 90 dd 90 origin 0 0 0
gen zone reflect dip 0 dd 0 origin 0 0 0
attach face range y -0.1 0.1
attach face range x -0.1 0.1
;依次联接y=0 x=0的交界面
Model mohr
pro bulk 8e9 shear 7e9 fric 25 coh 7.5e6 ten 2.6e6
set grav 0 10 0
ini dens 2000
fix x range x -10.1 -9.9
fix x range x 9.9 10.1
fix y range y 4.9 5.1
fix z range z 9.9 10.1
fix z range z -10.1 -9.9
hist unbal
hist gp zdisp -2 0 -6
;step 1000
;initial xdis 0.0 ydis 0.0 zdis 0.0
;set grav 0 -10 0
step 1000
;plo hist 1
plo con dis
plot con sxy ou on magf 20
plo add axes red
range name trench1 x=-5 3 y=0 5 z=-3 5
model null range trench1
;开挖
