gen zon radcyl p0 0 0 0 p1 8 0 0 p2 0 4 0 p3 0 0 8 size 4 8 16 16 dim 1.4 1.4 1.4 rat 1 1 1 1.1 group outsiderock
gen zon cshel p0 0 0 0 p1 1.4 0 0 p2 0 4 0 p3 0 0 1.4 size 2 8 16 4 dim 1.2 1.2 1.2 rat 1 1 1 1 &
group concretline fill group insiderock
gen zon brick p0 0 0 8 p1 8 0 8 p2 0 4 8 p3 0 0 24 size 8 8 16 group outsiderock1
gen zon brick p0 8 0 0 p1 24 0 0 p2 8 4 0 p3 8 0 24 size 16 8 24 group outsiderock2
gen zon reflect dip 90 dd 90 orig 0 0 0
gen zon reflect dip 0 dd 0 orig 0 0 0
plot block group
plot add axes red
m e
plot con szz ou on
plot add axes red
fix z range z -24.1 -23.9
fix z range z 24.1 23.9
fix x range x 23.9 24.1
fix x range x -24.1 -23.9
fix y ran y 0.1 -0.1
fix y ran y 3.9 4.1
ini szz -2.3e6 gradient 0 0 0.023e6 range z 24 -24
ini sxx -3.47e6 gradient 0 0 0.0297e6 range z 24 -24
ini syy -2.02e6 gradient 0 0 0.0303e6 range z 24 -24
ini density 2300
prop bulk=1.39e11 shear=0.57e11
set gravity 0 0 -10
set mech ratio=1e-6
solve
initial xdisp=0 ydisp=0 zdisp=0
ini xvel 0 yvel 0 zvel 0
m m
ini density 2300
prop bulk=1.39e9 shear=0.57e9 fric=34 coh=0.15e6 tension=1.5e6
mod null range group concretline any group insiderock any
;;;;;;;;;;;;;;;;;;;;混凝土衬砌
mod e range group concretline any
ini density 2200 range group concretline any
prop bulk=14e9 shear=8.4e9 range group concretline any
plo blo state
set mech ratio=1e-6
solve
save zhihusol.sav
