mm\Jf
new X{\>TOk
gen zone brick p0 0 0 0 p1 60 0 0 p2 0 9.9 0 p3 0 0 15 size 20 5 5 9+U%k(9
;gen zone brick p0 20 0 0 p1 60 0 0 p2 20 9.9 0 p3 20 0 15 size 20 3 5 ratio 1.1 048BQ
;gen zone brick p0 30 0 0 p1 60 0 0 p2 30 9.9 0 p3 30 0 15 size 6 2 3 GrjL9+|x
;attach face range x 19.9 20.1 _,_>B8
;attach face range x 29.9 30.1 :XcU @m
P!'Sx;C^f
p>h B &h
model elas ]R\k@a|G
prop bulk 3.58e10 shear 1.65e10 dens 2000 |[(4h
5c($3Pno=
Dg%zN i2GS
set grav 0 0 -10.0 VS@o_fUx)
ini szz -3.0e5 grad 0,0,2e4 g
#u1.|s&p
apply szz -20e3 range x 21.0 31.0 y 4.5 5.5 z 14.9 15.1 A;m)/@
fix x range x -0.1 0.1 )Xg,;^
fix x y z range x 59.9 60.1 tKrr5SRb
fix y range y -0.1 0.1 ,,S5 8\x
fix y range y 9.8 10.0 xekW-=#a7-
fix x y z range z -0.1 0.1 yZ)GP!cM4c
;hist unbal ,Y=r]
fk
;plot add hist 1 ?~uTbNR
p8hF`D~
YW9 [^
step 2000 _$PZID
ini xdisp 0.0 ydisp 0.0 zdisp 0.0 .9+"rK}u
;ini state 0 P;bl+a'gu
model mohr aAiSP+#
prop bu 3.58e6 sh 1.65e6 coh 24.1e3 fric 5.8 dens 1800 range z 15.0 14.0 ^%-$8sV
prop bu 3.58e6 sh 1.65e6 coh 24.1e3 fric 5.8 dens 2531 range z 14.0 0.0 VI0^Zq!6R
prop bu 8.6e8 sh 4.0e8 coh 5.06e5 fric 7.0 dens 1900 range x 20.0 20.8 z 3.0 15.0 G<OC99;8
;plot add fap black outside on bPEf2Z
G4
;plot add cont zdis out on )+RTA
y [k
;plot add axes black csz/[*
set large //]g78]=O
dX{|-;6vm
8<t6_* f
sel pile begin=(20.3,5.0,15.0) end=(20.3,5.0,3.0) nseg 12 )dgXS//Y
<KqZ.7XfB
\.2?951}
sel pile prop emod 1.0e1 nu 0.3 xcarea 0.7854 xcj 9.82e-2 xciy 4.91e-2 xciz 4.91e-2 & S-WD?BFC
per 3.14 cs_sk 1.0e1 cs_scoh 1.0e1 cs_sfric 0.0 cs_nk 1.0e1 cs_ncoh 1.0e1 cs_nfric 0.0 Or=
[2@Wg
c<