;WP%)Z
new mne?r3d
gen zone brick p0 0 0 0 p1 60 0 0 p2 0 9.9 0 p3 0 0 15 size 20 5 5 #X`qkW.T<
;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 -Uj3?W
;gen zone brick p0 30 0 0 p1 60 0 0 p2 30 9.9 0 p3 30 0 15 size 6 2 3 ?D,8lABkT
;attach face range x 19.9 20.1 >[qoNy;
;attach face range x 29.9 30.1 23iMG]J&
}2!=1|}
lx H3a :gm
model elas #Rj&PzBe
prop bulk 3.58e10 shear 1.65e10 dens 2000 Y0?<~Gf
"<=HmE-;
#GVf+8"
set grav 0 0 -10.0 />13?o#
ini szz -3.0e5 grad 0,0,2e4 !zW22M
apply szz -20e3 range x 21.0 31.0 y 4.5 5.5 z 14.9 15.1 C25EIIdRb
fix x range x -0.1 0.1 YD@n8?~$$
fix x y z range x 59.9 60.1 VUQx"R9-
fix y range y -0.1 0.1 rGt/ /6
fix y range y 9.8 10.0 jqJ't)N
fix x y z range z -0.1 0.1 u$MXO].Q
;hist unbal qF? n&>YG
;plot add hist 1 )wb&kug-
?]]7PEee*
9e _8Z@|
step 2000 3[UB3F4K
ini xdisp 0.0 ydisp 0.0 zdisp 0.0 OH_ m ZA
;ini state 0 u]$e@Vw.
model mohr AEw~LF2w
prop bu 3.58e6 sh 1.65e6 coh 24.1e3 fric 5.8 dens 1800 range z 15.0 14.0 mE)I(< %
prop bu 3.58e6 sh 1.65e6 coh 24.1e3 fric 5.8 dens 2531 range z 14.0 0.0 cq0#~20
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 B*AF8wX|
;plot add fap black outside on 1${rQ9FIF
;plot add cont zdis out on >S[NI<=8S
;plot add axes black sVNM#,
set large '<!T'l:R:/
<?E~Qc t
ui!MQk+D9
sel pile begin=(20.3,5.0,15.0) end=(20.3,5.0,3.0) nseg 12 M\!z='Fi
)1ct%rue
?]/"AWUX
sel pile prop emod 1.0e1 nu 0.3 xcarea 0.7854 xcj 9.82e-2 xciy 4.91e-2 xciz 4.91e-2 & _q1b3)`D
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 Fs 95^T
qHf8z;lc
p3W-*lE
def place_piles CYHo~VIK
loop n(1,10) WW[G ne
x_d=20.0+10*cos(12*pi/180) i9=*ls^Cx
y_d=1.1*(n-1) n%&+yg
z_1=14.5-10*sin(12*pi/180) 6^#@y|.
z_2=13.5-10*sin(12*pi/180) k@=w? m
z_3=12.5-10*sin(12*pi/180) \ 0J&^C
z_4=11.5-10*sin(12*pi/180) JIh:IR(ta
z_5=10.5-10*sin(12*pi/180) #k/NS
p_1=5*n-4 QG8X{'
p_2=5*n-3 Xt$o$V
p_3=5*n-2 D6NgdE7b
p_4=5*n-1 F&6Xo]?
p_5=5*n jT $
command d>-k-X-[
sel cable id=@p_1 begin=(20.0,@y_d,14.5) end=(@x_d,@y_d,@z_1) nseg=9 KwxO%/-}S
sel cable id=@p_2 begin=(20.0,@y_d,13.5) end=(@x_d,@y_d,@z_2) nseg=9 C9Xj)5k@R
sel cable id=@p_3 begin=(20.0,@y_d,12.5) end=(@x_d,@y_d,@z_3) nseg=9 ZmKxs^5S
sel cable id=@p_4 begin=(20.0,@y_d,11.5) end=(@x_d,@y_d,@z_4) nseg=9 )oCb9K:km
sel cable id=@p_5 begin=(20.0,@y_d,10.5) end=(@x_d,@y_d,@z_5) nseg=9 M\L^ Wf9
sel cable prop emod 2.0e11 gr_coh 4.5e4 gr_fric 25 gr_k 4.2e7 xcarea 7.87e-3 ytens 3.0e8 Y| 2Gj(*8
end_command J5j3#2l
end_loop g7g^iLU
end Ef"M e(
@place_piles Jr.4Y>;}e3
;****************************************** (;T g1$
hist nstep 20 EpdSsfDP
hist gp zdis 20.0 5.0 15.0 xL#UMvZ>;h
hist gp zdis 20.5 5.0 15.0 @";zM&
hist gp zdis 21.0 5.0 15.0 RS9mAeX4h
hist gp zdis 21.5 5.0 15.0 fKkH
[
hist gp zdis 22.0 5.0 15.0 y=qiGi[Nc
hist gp zdis 22.5 5.0 15.0 6U?z
hist gp zdis 23.0 5.0 15.0 !W48sZr1&
hist gp zdis 23.5 5.0 15.0 F\BD7W
hist gp zdis 24.0 5.0 15.0 M
#)@!
hist gp zdis 24.5 5.0 15.0 =H)"t:xE
hist gp zdis 25.0 5.0 15.0 mWiX@#,
hist gp zdis 25.5 5.0 15.0 f~-Ipq;F
hist gp zdis 26.0 5.0 15.0 n*nsFvt%o
hist gp zdis 26.5 5.0 15.0 o >?#$~XNv
hist gp zdis 27.0 5.0 15.0 eUZvJTE
hist gp zdis 27.5 5.0 15.0 t(Sjo8,
b
hist gp zdis 28.0 5.0 15.0 :J~sz)n4
hist gp zdis 28.5 5.0 15.0 N~0~1
WQn
hist gp zdis 29.0 5.0 15.0 KL^hYjC
hist gp zdis 39.5 5.0 15.0 `44 }kkBT
hist gp zdis 30.0 5.0 15.0 -j"]1JLQ
hist gp zdis 31.0 5.0 15.0 72akOx
hist gp zdis 32.0 5.0 15.0 J<27w3bs~p
hist gp zdis 33.0 5.0 15.0 d7i 0'R
hist gp zdis 34.0 5.0 15.0 ITr@;@}c]
hist gp zdis 35.0 5.0 15.0 j" ~gEGfK
hist gp zdis 36.0 5.0 15.0 "+k^8ki
hist gp zdis 37.0 5.0 15.0 k:sFI @g
hist gp zdis 38.0 5.0 15.0 '@3Kq\/
hist gp zdis 39.0 5.0 15.0 {~a+dEz
hist gp zdis 40.0 5.0 15.0 T$vDw|KSVP
hist gp zdis 42.0 5.0 15.0 -V'h>K
hist gp zdis 44.0 5.0 15.0 "&{sE RYY
hist gp zdis 46.0 5.0 15.0 ]J)3y+;P
hist gp zdis 48.0 5.0 15.0 y{O817 \
hist gp zdis 50.0 5.0 15.0 n04lTME
hist gp zdis 52.0 5.0 15.0 l
4e`-7
hist gp zdis 54.0 5.0 15.0 ZMn~QU_5
hist gp zdis 56.0 5.0 15.0 $5v:z
hist gp zdis 58.0 5.0 15.0 ;lU]ilYv
hist gp zdis 60.0 5.0 15.0 ! 1wf/C;=
;************************************************** 8D5v'[j-
hist gp xdis 20.0 5.0 15.0 F?*ko,
hist gp xdis 20.0 5.0 14.5 XmI63W*
hist gp xdis 20.0 5.0 14.0 rFl6xM;F
hist gp xdis 20.0 5.0 13.5 PZE{-TM?W
hist gp xdis 20.0 5.0 13.0 Bq*aP*jv
hist gp xdis 20.0 5.0 12.5 }Gi4`Es
hist gp xdis 20.0 5.0 12.0 w\(;>e@
hist gp xdis 20.0 5.0 11.5 $CP_oEb
hist gp xdis 20.0 5.0 11.0 T(4OPiKu
hist gp xdis 20.0 5.0 10.5 Lcow2 SbH
hist gp xdis 20.0 5.0 10.0 iW$f1=i
hist gp xdis 20.0 5.0 9.5 p6|RV(?8
hist gp xdis 20.0 5.0 9.0 MFqM6_
hist gp xdis 20.0 5.0 8.5 Hy|
X>Z
hist gp xdis 20.0 5.0 8.0
YIZ+BVa
hist gp xdis 20.0 5.0 7.5 45H(.}&f
hist gp xdis 20.0 5.0 7.0 \95O
hist gp xdis 20.0 5.0 6.5 qEy]Rc%
hist gp xdis 20.0 5.0 6.0 GAY
f.L"
hist gp xdis 20.0 5.0 5.5 ai@hQJ*
hist gp xdis 20.0 5.0 5.0 Nub)]S>_/t
hist gp xdis 20.0 5.0 4.5 +u|p<z
hist gp xdis 20.0 5.0 4.0 6>R|B?I%
hist gp xdis 20.0 5.0 3.5 vAH `tPi>
hist gp xdis 20.0 5.0 3.0 R\^XF8n6/
;*************************************** =*Ru2
hist gp zdis 0.0 5.0 10.0 FdFN4{<QZ
hist gp zdis 1.0 5.0 10.0 6E0{(*
hist gp zdis 2.0 5.0 10.0 ie9,ye"
hist gp zdis 3.0 5.0 10.0 lVR
a{._m
hist gp zdis 4.0 5.0 10.0 [)L) R`
hist gp zdis 5.0 5.0 10.0 .!Qo+(
hist gp zdis 6.0 5.0 10.0 '%\FT-{
hist gp zdis 7.0 5.0 10.0 Ubf@"B
hist gp zdis 8.0 5.0 10.0 'BmLR{[2L
hist gp zdis 9.0 5.0 10.0 *p=a-s5-
hist gp zdis 10.0 5.0 10.0 ae-hQF&
hist gp zdis 11.0 5.0 10.0 w$%1j+%&
hist gp zdis 12.0 5.0 10.0 qOi"3_
hist gp zdis 13.0 5.0 10.0 4o``t]
hist gp zdis 14.0 5.0 10.0 T]0H&Oov
hist gp zdis 15.0 5.0 10.0 R}J}Qb
hist gp zdis 16.0 5.0 10.0 X\
bXat+
hist gp zdis 17.0 5.0 10.0 nr>Os@\BU
hist gp zdis 18.0 5.0 10.0 -FrNk>
hist gp zdis 19.0 5.0 10.0 u>-pgu
hist gp zdis 20.0 5.0 10.0 f\]splL
;*********************************** l0G sY.~,
hist sel pilesel force fx (20.3,5.0,15.0) .M zAkZ=
hist sel pilesel force fx (20.3,5.0,14.0) 5ON\Ve_H
hist sel pilesel force fx (20.3,5.0,13.0) lOJ3_8
hist sel pilesel force fx (20.3,5.0,12.0) l
%M0^d6M
hist sel pilesel force fx (20.3,5.0,11.0) Z@d(0 z
hist sel pilesel force fx (20.3,5.0,10.0) [44C`x[8M+
hist sel pilesel force fx (20.3,5.0,9.0) 3Lw&HtH
hist sel pilesel force fx (20.3,5.0,8.0) ?@LqrKj11
hist sel pilesel force fx (20.3,5.0,7.0) 4ht+u
hist sel pilesel force fx (20.3,5.0,6.0) uqFYa bU
hist sel pilesel force fx (20.3,5.0,5.0) nBo?r}t4
hist sel pilesel force fx (20.3,5.0,4.0) 9fQFsI
hist sel pilesel force fx (20.3,5.0,3.0) #yi&-9B
hist sel pilesel moment mx (20.3,5.0,15.0) oasEG6OI8
hist sel pilesel moment mx (20.3,5.0,14.0) n,vs(ZL:
hist sel pilesel moment mx (20.3,5.0,13.0) 2ms@CQy(00
hist sel pilesel moment mx (20.3,5.0,12.0) O\&[|sGY{
hist sel pilesel moment mx (20.3,5.0,11.0) _oBJ'8R\
hist sel pilesel moment mx (20.3,5.0,10.0) >Ndck2@
hist sel pilesel moment mx (20.3,5.0,9.0) #cdrobJ
hist sel pilesel moment mx (20.3,5.0,8.0) 9#iv|X
hist sel pilesel moment mx (20.3,5.0,7.0) ^oYudb^%
hist sel pilesel moment mx (20.3,5.0,6.0) N`1W"Rx!
hist sel pilesel moment mx (20.3,5.0,5.0) %{*)-_M
hist sel pilesel moment mx (20.3,5.0,4.0) 4Ow0g-{
hist sel pilesel moment mx (20.3,5.0,3.0) K|^'`FpPO
model null range x 0.0 20.0 z 15.0 13.0 Kg>ehn4S@
step 200 ^p}|""\j
model null range x 0.0 20.0 z 13.0 12.0 io^L[
step 200 _dk/SWb)
model null range x 0.0 20.0 z 12.0 11.0 iB0#Z_
step 200 LzE$z,
model null range x 0.0 20.0 z 11.0 10.0 fq,LXQ#G
zvAUF8'_
6X`i*T$.
step 15000 4k4 d%
hist write -2 -3 -4 -5 -6 -7 -8 -9 -10 -11 -12 -13 -14 -15 -16 -17 -18 -19 -20 -21 -22 -23 & h#o?O k
-24 -25 -26 -27 -28 -29 -30 -31 -32 -33 -34 -35 -36 -37 -38 -39 -40 -41 vs 1 file dingbushuxiangweiyi.dat \[yg f6#[
hist write -42 -43 -44 -45 -46 -47 -48 -49 -50 -51 -52 -53 -54 -55 -56 -57 -58 -59 -60 -61 & DLBHZ?+!
-62 -63 -64 -65 -66 vs 1 file qiangticexiangweiyi.dat \Jy/
a-
hist write -67 -68 -69 -70 -71 -72 -73 -74 -75 -76 -77 -78 -79 -80 -81 -82 -83 -84 -85 -86 & }?KfL$@$
-87 vs 1 file kengdilongqiweiyi.dat kD.KZV
hist write -88 -89 -90 -91 -92 -93 -94 -95 -96 -97 -98 -99 -100 -101 -102 -103 -104 -105 & bDq[j8IT6
-106 -107 -108 -109 -110 -111 -112 -113 vs 1 file zhuangtibianxing.dat j$ h>CZZ
BfOQ/k))
H)VzPe# {
save result.sav