?)60JWOJ1
new MV~-']2u
gen zone cyl p0 0 0 0 p1 1 0 0 p2 0 2 0 p3 0 0 1 size 4 5 4 r=`]L-}V
gen zone reflect norm 1,0,0 u`Sg' ro
gen zone reflect norm 0,0,1 HZG<aY="
bSKe@4C
model ss oD0N<Ln}
prop bulk 1.19e6 shear 1.1e6 dens 1500 FTA[O.tiG
prop coh 2.72e10 fric 44 ten 2e5 ctab 1 ftab 2 #U=}Pv~wM
table 1 0 100e3 0.05 50e3 0.1 25e3 1 25e3 |.q K69
table 2 0 23 0.05 20 0.1 18 1 18 gmU0/z3&
f:).wi
Ld
fix x y z range y -.1 .1 Gp PlO]
fix x y z range y 1.9 2.1 v4YY6?4
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 {e3XmVAI
ini yvel 1e-7 range y -.1 .1 ,)CRozC\}K
>o#^)LN
hist gp ydisp 0,0,0 ; 4;_<CB
hist zone syy 0,0,0 ; ~kkwPs2V
hist zone syy 0,1,0 ; +pq/:h
step 1000 Z,? T`[4B
2f=7`1RCD
plot hist -2 vs 1 --32kuF&(
hist write 1 2 file 单压.txt