4id3P{aU
new M'\pkzx
gen zone brick p0 0 0 0 p1 60 0 0 p2 0 9.9 0 p3 0 0 15 size 20 5 5 D3 %l4.h
;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 [urH a
;gen zone brick p0 30 0 0 p1 60 0 0 p2 30 9.9 0 p3 30 0 15 size 6 2 3 ')C|`(hs
;attach face range x 19.9 20.1 iY0,WT}&n
;attach face range x 29.9 30.1 cwu$TP A>
`]K,'i{ R
L3B8IDq
model elas ;c>>$lr
prop bulk 3.58e10 shear 1.65e10 dens 2000 0SJ{@ *
} c{Fa&
7'_nc!ME
set grav 0 0 -10.0 =a?a@+
ini szz -3.0e5 grad 0,0,2e4 LIr(mB"Y0
apply szz -20e3 range x 21.0 31.0 y 4.5 5.5 z 14.9 15.1 ':,>eL#+uV
fix x range x -0.1 0.1 R]CZw;zS_
fix x y z range x 59.9 60.1 X,>(Y8
fix y range y -0.1 0.1 HR[ Q
?rg
fix y range y 9.8 10.0 U:qF/%w
fix x y z range z -0.1 0.1 'Z\{D*=V8
;hist unbal ]@*tfz\YaH
;plot add hist 1 "w*@R 8v
xe(MHNrj
&}zRH}s;
step 2000 U+4HG
ini xdisp 0.0 ydisp 0.0 zdisp 0.0 b\1+kB/8
;ini state 0 LkaG8#m1R
model mohr 'oC$6l'rQ
prop bu 3.58e6 sh 1.65e6 coh 24.1e3 fric 5.8 dens 1800 range z 15.0 14.0 d4ic9u*D
prop bu 3.58e6 sh 1.65e6 coh 24.1e3 fric 5.8 dens 2531 range z 14.0 0.0 )*!1bgXQ
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 M} O[`Fx{W
;plot add fap black outside on k?^%hO>[
;plot add cont zdis out on s,84* 6u
;plot add axes black ,q8(]n 4
set large %z1{Kus
(-bRj#
z8b
_ _%Br
sel pile begin=(20.3,5.0,15.0) end=(20.3,5.0,3.0) nseg 12 M|U';2hZN:
+``>,O6
%v]7BV^%6
sel pile prop emod 1.0e1 nu 0.3 xcarea 0.7854 xcj 9.82e-2 xciy 4.91e-2 xciz 4.91e-2 & q z=yMIy=
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 Ru)(dvk}S
b![t6-f^z
e@[9C(5E"
def place_piles hgj CXl
loop n(1,10) >RM
0=bO
x_d=20.0+10*cos(12*pi/180) HKpD 2M
y_d=1.1*(n-1) G-2EQ .
z_1=14.5-10*sin(12*pi/180) +jcdf }
z_2=13.5-10*sin(12*pi/180) DZJ eup?Z
z_3=12.5-10*sin(12*pi/180) 4w @v# H@
z_4=11.5-10*sin(12*pi/180) h=h oV5d@
z_5=10.5-10*sin(12*pi/180) GSQfg
p_1=5*n-4 DeA @0HOxh
p_2=5*n-3 7. %f01/i
p_3=5*n-2 wOjv[@d
p_4=5*n-1 - dl}_
p_5=5*n BpH|/7
command
LlU' _}>
sel cable id=@p_1 begin=(20.0,@y_d,14.5) end=(@x_d,@y_d,@z_1) nseg=9 BApa^j\?
sel cable id=@p_2 begin=(20.0,@y_d,13.5) end=(@x_d,@y_d,@z_2) nseg=9 ,MRvuw0P
sel cable id=@p_3 begin=(20.0,@y_d,12.5) end=(@x_d,@y_d,@z_3) nseg=9 9^oo-,Su_
sel cable id=@p_4 begin=(20.0,@y_d,11.5) end=(@x_d,@y_d,@z_4) nseg=9 jUY F.K&
sel cable id=@p_5 begin=(20.0,@y_d,10.5) end=(@x_d,@y_d,@z_5) nseg=9 ] yo_wGiwY
sel cable prop emod 2.0e11 gr_coh 4.5e4 gr_fric 25 gr_k 4.2e7 xcarea 7.87e-3 ytens 3.0e8 f b /qoZ
end_command
+q7qK*
end_loop l x7Kw%
end w-lrnjs
@place_piles ;JYoW{2
;****************************************** yP
x\ltG3
hist nstep 20 ]+AAT=B<!
hist gp zdis 20.0 5.0 15.0 {&uT3*V 1
hist gp zdis 20.5 5.0 15.0 9 >%+bA(
hist gp zdis 21.0 5.0 15.0 ^ b@! dS
hist gp zdis 21.5 5.0 15.0 }gCG&7C
hist gp zdis 22.0 5.0 15.0 -=}b;Kf -
hist gp zdis 22.5 5.0 15.0 r WJ*e Y
hist gp zdis 23.0 5.0 15.0 B gf=\7;5
hist gp zdis 23.5 5.0 15.0 TNx _Rc}
hist gp zdis 24.0 5.0 15.0 i9KTX%s5^
hist gp zdis 24.5 5.0 15.0 g+.0c=G(
hist gp zdis 25.0 5.0 15.0 <p09oZ {6
hist gp zdis 25.5 5.0 15.0 [ qiOd!
hist gp zdis 26.0 5.0 15.0 U>oW~Z
hist gp zdis 26.5 5.0 15.0 Im6U_JsNZh
hist gp zdis 27.0 5.0 15.0 Q2q| *EL
hist gp zdis 27.5 5.0 15.0 D {>\- ]\
hist gp zdis 28.0 5.0 15.0
t7&Dwmck9
hist gp zdis 28.5 5.0 15.0 sqT^t!
hist gp zdis 29.0 5.0 15.0 97!H`|u <
hist gp zdis 39.5 5.0 15.0 : aH%bk
hist gp zdis 30.0 5.0 15.0 \OP9_J(*
hist gp zdis 31.0 5.0 15.0 A$Wx#r7)
hist gp zdis 32.0 5.0 15.0 6>s=Ci ZB
hist gp zdis 33.0 5.0 15.0 !Pw$4 8cg
hist gp zdis 34.0 5.0 15.0 XYts8}y5
hist gp zdis 35.0 5.0 15.0 "i&fp:E0
hist gp zdis 36.0 5.0 15.0 {f-XyF1`
hist gp zdis 37.0 5.0 15.0 `6P?G|'
hist gp zdis 38.0 5.0 15.0 {b2 aL7
hist gp zdis 39.0 5.0 15.0 z<t>h zl 7
hist gp zdis 40.0 5.0 15.0 > <X $#
hist gp zdis 42.0 5.0 15.0 Zlrh C= 0
hist gp zdis 44.0 5.0 15.0 {(%~i3 7
hist gp zdis 46.0 5.0 15.0 0:. 6rp
hist gp zdis 48.0 5.0 15.0 ":V%(c
hist gp zdis 50.0 5.0 15.0 #J\s%60p t
hist gp zdis 52.0 5.0 15.0 dKb ^x^
hist gp zdis 54.0 5.0 15.0 ~zMDY F"&
hist gp zdis 56.0 5.0 15.0 *(icR
hist gp zdis 58.0 5.0 15.0 b )Tl*
hist gp zdis 60.0 5.0 15.0 Vx kEe z'|
;************************************************** X CzXS .
hist gp xdis 20.0 5.0 15.0 `&H04x"Y$>
hist gp xdis 20.0 5.0 14.5 @O'I)(To
hist gp xdis 20.0 5.0 14.0 b TiBmS
hist gp xdis 20.0 5.0 13.5 ZE qE$:
hist gp xdis 20.0 5.0 13.0 9Yn)t#G'`F
hist gp xdis 20.0 5.0 12.5 O)`Gzx*ShU
hist gp xdis 20.0 5.0 12.0 s0 47"Q
hist gp xdis 20.0 5.0 11.5 g** 5z'7
hist gp xdis 20.0 5.0 11.0 '<W<B!HP5Z
hist gp xdis 20.0 5.0 10.5 e9lOk)`t
hist gp xdis 20.0 5.0 10.0 z#RuwB+
hist gp xdis 20.0 5.0 9.5 J]dW1boT@
hist gp xdis 20.0 5.0 9.0 ^@K
WYAAW5
hist gp xdis 20.0 5.0 8.5 rEmwKZ F'
hist gp xdis 20.0 5.0 8.0 Si]X
rub
hist gp xdis 20.0 5.0 7.5 -n-Z/5~ X
hist gp xdis 20.0 5.0 7.0 (V!0'9c
hist gp xdis 20.0 5.0 6.5 p= vV4 C:
hist gp xdis 20.0 5.0 6.0 'aZAS Pn[
hist gp xdis 20.0 5.0 5.5 _\UIc;3Gl
hist gp xdis 20.0 5.0 5.0 2
^oGwx @
hist gp xdis 20.0 5.0 4.5 Wa<-AZnh
hist gp xdis 20.0 5.0 4.0 xdqK.Z%
hist gp xdis 20.0 5.0 3.5 fQO
""qh
hist gp xdis 20.0 5.0 3.0 U:\p$ hL9
;*************************************** c`ftd>]
hist gp zdis 0.0 5.0 10.0 Sj@15 W
hist gp zdis 1.0 5.0 10.0 ** n y!
hist gp zdis 2.0 5.0 10.0 94y9W#
hist gp zdis 3.0 5.0 10.0 V,m3-=q
hist gp zdis 4.0 5.0 10.0 xvB8YW"
hist gp zdis 5.0 5.0 10.0 { l@WCR
hist gp zdis 6.0 5.0 10.0 n_}aZB3;U
hist gp zdis 7.0 5.0 10.0 T=>vh*J
hist gp zdis 8.0 5.0 10.0 qW 3x {L$c
hist gp zdis 9.0 5.0 10.0 m2Q$+p@
hist gp zdis 10.0 5.0 10.0 ~XKZXGw
hist gp zdis 11.0 5.0 10.0 EWO / u.z
hist gp zdis 12.0 5.0 10.0 SwX@I6huM
hist gp zdis 13.0 5.0 10.0 p jd o|
hist gp zdis 14.0 5.0 10.0 oBC]UL;8xJ
hist gp zdis 15.0 5.0 10.0 > e/ r2U
hist gp zdis 16.0 5.0 10.0 z>p ]/Sa
hist gp zdis 17.0 5.0 10.0 {pC\\}
hist gp zdis 18.0 5.0 10.0 g8'~e{= (
hist gp zdis 19.0 5.0 10.0 `6}Yqh))
hist gp zdis 20.0 5.0 10.0 [1U{ci&=p
;*********************************** "O``7HA}
hist sel pilesel force fx (20.3,5.0,15.0) y]
y9'5_
hist sel pilesel force fx (20.3,5.0,14.0) % 0zS
hist sel pilesel force fx (20.3,5.0,13.0) S}b~_}
hist sel pilesel force fx (20.3,5.0,12.0) 6uq UiRs()
hist sel pilesel force fx (20.3,5.0,11.0) 9$wAm89
hist sel pilesel force fx (20.3,5.0,10.0) <S&]$?`{Wi
hist sel pilesel force fx (20.3,5.0,9.0) 5e8xKL
hist sel pilesel force fx (20.3,5.0,8.0) ki#bPgT
hist sel pilesel force fx (20.3,5.0,7.0) WGPD8.
hist sel pilesel force fx (20.3,5.0,6.0) J)KnE2dw5
hist sel pilesel force fx (20.3,5.0,5.0) F<gMUDB
hist sel pilesel force fx (20.3,5.0,4.0) #"<?_fao~
hist sel pilesel force fx (20.3,5.0,3.0) J
3B`Krh
hist sel pilesel moment mx (20.3,5.0,15.0) hJLT!33:
hist sel pilesel moment mx (20.3,5.0,14.0) {!r#f(?uT
hist sel pilesel moment mx (20.3,5.0,13.0) _ ~[M+IO
hist sel pilesel moment mx (20.3,5.0,12.0) ;?`@"YG)
hist sel pilesel moment mx (20.3,5.0,11.0) iu|v9+
hist sel pilesel moment mx (20.3,5.0,10.0) "[.ne)/MC
hist sel pilesel moment mx (20.3,5.0,9.0) F 3s?&T)[G
hist sel pilesel moment mx (20.3,5.0,8.0) DN<M ?u]
hist sel pilesel moment mx (20.3,5.0,7.0) ?<6@^X"
hist sel pilesel moment mx (20.3,5.0,6.0) AOAO8%|I
hist sel pilesel moment mx (20.3,5.0,5.0) \OY}GR Kt
hist sel pilesel moment mx (20.3,5.0,4.0) /?U!y?t&@
hist sel pilesel moment mx (20.3,5.0,3.0) %I0}4$
model null range x 0.0 20.0 z 15.0 13.0 v^TkDf(Oz
step 200 e[8UH =`|
model null range x 0.0 20.0 z 13.0 12.0 <]'|$8&jY
step 200
QAUykS8
model null range x 0.0 20.0 z 12.0 11.0 o} {-j
step 200 t#~XL CE
model null range x 0.0 20.0 z 11.0 10.0 *6P 'q4 )
e=L*&X
\%4|t,en
step 15000 hkF^? AJ
hist write -2 -3 -4 -5 -6 -7 -8 -9 -10 -11 -12 -13 -14 -15 -16 -17 -18 -19 -20 -21 -22 -23 & D J_DonO]
-24 -25 -26 -27 -28 -29 -30 -31 -32 -33 -34 -35 -36 -37 -38 -39 -40 -41 vs 1 file dingbushuxiangweiyi.dat M $uf:+F
hist write -42 -43 -44 -45 -46 -47 -48 -49 -50 -51 -52 -53 -54 -55 -56 -57 -58 -59 -60 -61 & sG1BNb_
-62 -63 -64 -65 -66 vs 1 file qiangticexiangweiyi.dat ST%
T =_q
hist write -67 -68 -69 -70 -71 -72 -73 -74 -75 -76 -77 -78 -79 -80 -81 -82 -83 -84 -85 -86 &
mV;3I LO
-87 vs 1 file kengdilongqiweiyi.dat N|<bVq%
hist write -88 -89 -90 -91 -92 -93 -94 -95 -96 -97 -98 -99 -100 -101 -102 -103 -104 -105 & [<S^c[47U
-106 -107 -108 -109 -110 -111 -112 -113 vs 1 file zhuangtibianxing.dat A2 BRbwr>
-N4z-ozhC
GXYj+ qJ
save result.sav
