Ba6xkEd
new Pz 'Hqvd
gen zone brick p0 0 0 0 p1 60 0 0 p2 0 9.9 0 p3 0 0 15 size 20 5 5 j66@E\dN
;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 v4$,Vt:7
;gen zone brick p0 30 0 0 p1 60 0 0 p2 30 9.9 0 p3 30 0 15 size 6 2 3 ,3j7Y5v
;attach face range x 19.9 20.1 H/$q]i*#K
;attach face range x 29.9 30.1 hs4r5[
*"ShE=\p
*C BCQp[$
model elas <Va>5R_d<
prop bulk 3.58e10 shear 1.65e10 dens 2000 bhmjH(.t
(
~>Q2DS
.kIf1-(<U
set grav 0 0 -10.0 C#Jj;Gd
ini szz -3.0e5 grad 0,0,2e4 S\A[Z&k0
apply szz -20e3 range x 21.0 31.0 y 4.5 5.5 z 14.9 15.1 %vXQ Sz
fix x range x -0.1 0.1 hd~rC*I
fix x y z range x 59.9 60.1 DFonK{
fix y range y -0.1 0.1 Oq5k4
fix y range y 9.8 10.0 UwU]l17~
fix x y z range z -0.1 0.1 3fPv71NVtt
;hist unbal UL%ihWq
;plot add hist 1 A=K1T]o
(R-Q9F+;
MF`'r#@:wa
step 2000 3[0:,^a
ini xdisp 0.0 ydisp 0.0 zdisp 0.0 yKJ^hv"#
;ini state 0 Ei-OuDM;)
model mohr E`|qFG<
prop bu 3.58e6 sh 1.65e6 coh 24.1e3 fric 5.8 dens 1800 range z 15.0 14.0 (XJQ$n
prop bu 3.58e6 sh 1.65e6 coh 24.1e3 fric 5.8 dens 2531 range z 14.0 0.0 r.^&%D
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 7 SZR#L
;plot add fap black outside on ${wE5^ky
;plot add cont zdis out on :+Kesa:E
;plot add axes black MeX1y]<It
set large n&]w* (,
BpT&vbY
m!_ghD{5h
sel pile begin=(20.3,5.0,15.0) end=(20.3,5.0,3.0) nseg 12 67x^{u7
jq)|Uq'6
jH1~Ve+q9
sel pile prop emod 1.0e1 nu 0.3 xcarea 0.7854 xcj 9.82e-2 xciy 4.91e-2 xciz 4.91e-2 & bed+Ur&
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 vknFtpx
9YABr>
?
YC'~8\x3z
def place_piles $b} +5
loop n(1,10) @Hh"Y1B
x_d=20.0+10*cos(12*pi/180) *$VurqLn
y_d=1.1*(n-1) In&vh9Lw
z_1=14.5-10*sin(12*pi/180) M6iKl
z_2=13.5-10*sin(12*pi/180) %a `dOEO
z_3=12.5-10*sin(12*pi/180) 7Cf(y'w^
z_4=11.5-10*sin(12*pi/180) sXxO{aeev
z_5=10.5-10*sin(12*pi/180) bSLj-vp
p_1=5*n-4 GHY>DrXO1u
p_2=5*n-3 R{_IrYk
p_3=5*n-2 U4gJ![>5j
p_4=5*n-1 mQd?Tyvn
p_5=5*n }3 }=tN5
command =HHg:"
sel cable id=@p_1 begin=(20.0,@y_d,14.5) end=(@x_d,@y_d,@z_1) nseg=9 ([~`{,sv
sel cable id=@p_2 begin=(20.0,@y_d,13.5) end=(@x_d,@y_d,@z_2) nseg=9 _=5ZB_I
sel cable id=@p_3 begin=(20.0,@y_d,12.5) end=(@x_d,@y_d,@z_3) nseg=9 Mk[`HEO
sel cable id=@p_4 begin=(20.0,@y_d,11.5) end=(@x_d,@y_d,@z_4) nseg=9 Kdm5O@tq
sel cable id=@p_5 begin=(20.0,@y_d,10.5) end=(@x_d,@y_d,@z_5) nseg=9 YqgW8EM
sel cable prop emod 2.0e11 gr_coh 4.5e4 gr_fric 25 gr_k 4.2e7 xcarea 7.87e-3 ytens 3.0e8 SO/]d70HG
end_command k6BgY|0gC
end_loop pZxL?N!
end j&.BbcE45
@place_piles 5ov%(QI
;****************************************** 1tW:(~=a;
hist nstep 20 :(Bi{cw
hist gp zdis 20.0 5.0 15.0 Fev3CV$
hist gp zdis 20.5 5.0 15.0 Z&,}Fgl!F
hist gp zdis 21.0 5.0 15.0 qS
al~
hist gp zdis 21.5 5.0 15.0 =6&D4~R
hist gp zdis 22.0 5.0 15.0 )v~]lk,o
hist gp zdis 22.5 5.0 15.0 [2V/v
hist gp zdis 23.0 5.0 15.0 cmI#R1\
hist gp zdis 23.5 5.0 15.0 I.!/R`
hist gp zdis 24.0 5.0 15.0 ub5hX{uT
hist gp zdis 24.5 5.0 15.0 b`zf&Mn
hist gp zdis 25.0 5.0 15.0 U@nwSfp:G
hist gp zdis 25.5 5.0 15.0 ]6 wi
hist gp zdis 26.0 5.0 15.0 A]$+
`uS\
hist gp zdis 26.5 5.0 15.0 :8rCCop
Uv
hist gp zdis 27.0 5.0 15.0 k#xpY!'7
hist gp zdis 27.5 5.0 15.0 `'WLGQG
hist gp zdis 28.0 5.0 15.0 *\", qMp
hist gp zdis 28.5 5.0 15.0 Kf#!IY][
hist gp zdis 29.0 5.0 15.0 3G^Ed)JvE
hist gp zdis 39.5 5.0 15.0 5eA]7$ic
hist gp zdis 30.0 5.0 15.0 *.g?y6d
hist gp zdis 31.0 5.0 15.0 99K+7G\{
hist gp zdis 32.0 5.0 15.0 \<**SSN
hist gp zdis 33.0 5.0 15.0 N&=2 /
hist gp zdis 34.0 5.0 15.0 <J-Z;r(gQN
hist gp zdis 35.0 5.0 15.0 Sj?sw]3
hist gp zdis 36.0 5.0 15.0 _+iz?|U
hist gp zdis 37.0 5.0 15.0 [[Z>(d$8
hist gp zdis 38.0 5.0 15.0 7`HUwu
hist gp zdis 39.0 5.0 15.0 sfE8b/Z8
hist gp zdis 40.0 5.0 15.0 /&7Yi_]r
hist gp zdis 42.0 5.0 15.0 HU9y{H
hist gp zdis 44.0 5.0 15.0 lvi:I+VgA
hist gp zdis 46.0 5.0 15.0 (_ah~VnO
hist gp zdis 48.0 5.0 15.0 JB@VP{
hist gp zdis 50.0 5.0 15.0 avu,o
hist gp zdis 52.0 5.0 15.0 UI C? S
hist gp zdis 54.0 5.0 15.0 ;!?K.,N:N
hist gp zdis 56.0 5.0 15.0 vq1&8=
hist gp zdis 58.0 5.0 15.0 o"[bIXf-h
hist gp zdis 60.0 5.0 15.0 ,np`:fBMy
;************************************************** $:!T/*p*
hist gp xdis 20.0 5.0 15.0 ;0}2@Q2@ZK
hist gp xdis 20.0 5.0 14.5 4sjr\9IDC
hist gp xdis 20.0 5.0 14.0 mC92J@m/L!
hist gp xdis 20.0 5.0 13.5 i5Q<~;Z+
hist gp xdis 20.0 5.0 13.0 L0tKIpk
hist gp xdis 20.0 5.0 12.5 VhgEG(Ud
hist gp xdis 20.0 5.0 12.0 NCt sx /C
hist gp xdis 20.0 5.0 11.5 i&)C,
hist gp xdis 20.0 5.0 11.0 6a?p?I K^
hist gp xdis 20.0 5.0 10.5 2]=I'U<E!
hist gp xdis 20.0 5.0 10.0 o[hP&9>q
hist gp xdis 20.0 5.0 9.5 D|9fHMg%
hist gp xdis 20.0 5.0 9.0 Xkom@F~]
hist gp xdis 20.0 5.0 8.5 vWs c{9
hist gp xdis 20.0 5.0 8.0 ton`ji\^
hist gp xdis 20.0 5.0 7.5 Q%Q?q)x
hist gp xdis 20.0 5.0 7.0 =fK'Ep[
hist gp xdis 20.0 5.0 6.5 <t% A)L%
hist gp xdis 20.0 5.0 6.0 om?CFl
hist gp xdis 20.0 5.0 5.5 VY@hhr1s~
hist gp xdis 20.0 5.0 5.0 u V7Hsg9l
hist gp xdis 20.0 5.0 4.5 T0%TeFY
hist gp xdis 20.0 5.0 4.0 tYZGf xj
hist gp xdis 20.0 5.0 3.5 J|S^K kC
hist gp xdis 20.0 5.0 3.0 <|3v@
;*************************************** _9|@nUD
hist gp zdis 0.0 5.0 10.0 C}_:K)5q
hist gp zdis 1.0 5.0 10.0 ]xEE7H]\h
hist gp zdis 2.0 5.0 10.0 Y{RB\}f(
hist gp zdis 3.0 5.0 10.0 yuEOQ\!(u
hist gp zdis 4.0 5.0 10.0 ~79Qg{+]N
hist gp zdis 5.0 5.0 10.0 p]Zabky
hist gp zdis 6.0 5.0 10.0 Tj5@OcA$
hist gp zdis 7.0 5.0 10.0 J5_Y\@
hist gp zdis 8.0 5.0 10.0 vfJk?
(
hist gp zdis 9.0 5.0 10.0 WG}CPkj
hist gp zdis 10.0 5.0 10.0 4uAafQ`@H
hist gp zdis 11.0 5.0 10.0 y51D-vj
hist gp zdis 12.0 5.0 10.0
"B3:m-'
hist gp zdis 13.0 5.0 10.0 E^a`IA
hist gp zdis 14.0 5.0 10.0 \/'n[3x
hist gp zdis 15.0 5.0 10.0 IQe[ CcM
hist gp zdis 16.0 5.0 10.0 5C1Rub)
hist gp zdis 17.0 5.0 10.0 CL :M>(
hist gp zdis 18.0 5.0 10.0 K"j=_%{
hist gp zdis 19.0 5.0 10.0 Ag0_^
hist gp zdis 20.0 5.0 10.0 5KE%@,k k
;*********************************** 8p{
hist sel pilesel force fx (20.3,5.0,15.0) Ml?)Sc"\7
hist sel pilesel force fx (20.3,5.0,14.0) U^BXCu1km
hist sel pilesel force fx (20.3,5.0,13.0) x=Jn&4q
hist sel pilesel force fx (20.3,5.0,12.0) 2_n*u^X:_
hist sel pilesel force fx (20.3,5.0,11.0) 6xh#;+e}
hist sel pilesel force fx (20.3,5.0,10.0) 3Lki7QW`
hist sel pilesel force fx (20.3,5.0,9.0) _PUm
Pom.
hist sel pilesel force fx (20.3,5.0,8.0) ID#qKFFW
hist sel pilesel force fx (20.3,5.0,7.0) Gj`Y2X2r
hist sel pilesel force fx (20.3,5.0,6.0) &xroms"S=
hist sel pilesel force fx (20.3,5.0,5.0) ?3,tG z)
hist sel pilesel force fx (20.3,5.0,4.0) j%jd@z ]@
hist sel pilesel force fx (20.3,5.0,3.0) OB^?cA>
hist sel pilesel moment mx (20.3,5.0,15.0) 2Q|*xd4B^
hist sel pilesel moment mx (20.3,5.0,14.0) 5dw@g4N %^
hist sel pilesel moment mx (20.3,5.0,13.0) UMQW#$~C{g
hist sel pilesel moment mx (20.3,5.0,12.0) FNCLGAiZ
hist sel pilesel moment mx (20.3,5.0,11.0) 3}{5
X'
hist sel pilesel moment mx (20.3,5.0,10.0) UQ])QTrZFi
hist sel pilesel moment mx (20.3,5.0,9.0) x*8f3^ wE
hist sel pilesel moment mx (20.3,5.0,8.0) h\*rv5\M
hist sel pilesel moment mx (20.3,5.0,7.0) O)%kl
hist sel pilesel moment mx (20.3,5.0,6.0) zN/~a)
hist sel pilesel moment mx (20.3,5.0,5.0) [.xk
hist sel pilesel moment mx (20.3,5.0,4.0) (!5}" fj
hist sel pilesel moment mx (20.3,5.0,3.0) cjC6\.+l3
model null range x 0.0 20.0 z 15.0 13.0 DN':-PK
step 200 l~TIFmHkh%
model null range x 0.0 20.0 z 13.0 12.0 u%T$XG
step 200 |!5T+H{Sj
model null range x 0.0 20.0 z 12.0 11.0 \Rqh|T<D
step 200 9w;J7jgOT!
model null range x 0.0 20.0 z 11.0 10.0 =^y{@[p`(
vr;Br-8
Z !25xqNCd
step 15000 w })Pedg
hist write -2 -3 -4 -5 -6 -7 -8 -9 -10 -11 -12 -13 -14 -15 -16 -17 -18 -19 -20 -21 -22 -23 & J?quYlS
-24 -25 -26 -27 -28 -29 -30 -31 -32 -33 -34 -35 -36 -37 -38 -39 -40 -41 vs 1 file dingbushuxiangweiyi.dat umZ
g}|C_
hist write -42 -43 -44 -45 -46 -47 -48 -49 -50 -51 -52 -53 -54 -55 -56 -57 -58 -59 -60 -61 & cN}A rv
-62 -63 -64 -65 -66 vs 1 file qiangticexiangweiyi.dat *jw$d8q2
hist write -67 -68 -69 -70 -71 -72 -73 -74 -75 -76 -77 -78 -79 -80 -81 -82 -83 -84 -85 -86 &
%%cSvPcz
-87 vs 1 file kengdilongqiweiyi.dat $1zeY6O
hist write -88 -89 -90 -91 -92 -93 -94 -95 -96 -97 -98 -99 -100 -101 -102 -103 -104 -105 & Cmx2/N
-106 -107 -108 -109 -110 -111 -112 -113 vs 1 file zhuangtibianxing.dat )-KE4/G
F%Umau*1
m_02"'
save result.sav