_?M34&.X
new V$VqYy9 *
gen zone cyl p0 0 0 0 p1 1 0 0 p2 0 2 0 p3 0 0 1 size 4 5 4 ?>cx;"xF
gen zone reflect norm 1,0,0 LdwWB
`L
gen zone reflect norm 0,0,1 ri1D*CS
z6]dF"N
model ss q.U` mtS
prop bulk 1.19e6 shear 1.1e6 dens 1500 s]50Y-C
prop coh 2.72e10 fric 44 ten 2e5 ctab 1 ftab 2 ~m8".Z"
table 1 0 100e3 0.05 50e3 0.1 25e3 1 25e3 rCGXHbj%
table 2 0 23 0.05 20 0.1 18 1 18 +w[vYKSZm
L3GA]TIf
fix x y z range y -.1 .1 Ci4`,
fix x y z range y 1.9 2.1 VdjS\VYe,
fix x y z range cylinder end1 0,0,0 end2 0,2,0 radius 1.1 cyl end1 0,0,0 end2 0,2,0 radius 0.9 not Yrs7F.Y"
ini yvel 1e-7 range y -.1 .1 NQz*P.q
bGOOC?[UX
hist gp ydisp 0,0,0 ; /W1!mih
hist zone syy 0,0,0 ; <qT[
hist zone syy 0,1,0 ; ?1*Ka
step 1000 jzI70+E
de1cl<
plot hist -2 vs 1 Y#S<:,/sb?
hist write 1 2 file 单压.txt