`Up3p24
new 8*lVO2
gen zone brick p0 0 0 0 p1 60 0 0 p2 0 9.9 0 p3 0 0 15 size 20 5 5 $2\OBc=
;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 Q
xKC5`1
;gen zone brick p0 30 0 0 p1 60 0 0 p2 30 9.9 0 p3 30 0 15 size 6 2 3 mhVSZhx|
;attach face range x 19.9 20.1 N5o jXX!l%
;attach face range x 29.9 30.1 \|Us/_h
GsxrqIaD
eAPGy-
model elas (rvK@
prop bulk 3.58e10 shear 1.65e10 dens 2000 E QMn'>
p'!cGJL
8I`>tY
set grav 0 0 -10.0 4s%vx]E
ini szz -3.0e5 grad 0,0,2e4 L%Me
wU0TZ
apply szz -20e3 range x 21.0 31.0 y 4.5 5.5 z 14.9 15.1 g&X$)V4C
fix x range x -0.1 0.1 4y:yFTp
fix x y z range x 59.9 60.1 |&3[YZY
fix y range y -0.1 0.1 `cCsJm$V"
fix y range y 9.8 10.0 &h')snp:#
fix x y z range z -0.1 0.1 YG$Y4h"
@"
;hist unbal E]i3E[T
;plot add hist 1 +%oXPG?
'
FF@I^O
,d$V-~2,
step 2000 HW)> `
ini xdisp 0.0 ydisp 0.0 zdisp 0.0 QT&2&#Z
;ini state 0 G
D$o|l]\
model mohr KHI-m9(
prop bu 3.58e6 sh 1.65e6 coh 24.1e3 fric 5.8 dens 1800 range z 15.0 14.0 ]*D=^kA0[
prop bu 3.58e6 sh 1.65e6 coh 24.1e3 fric 5.8 dens 2531 range z 14.0 0.0 %;~Vc{Xxt/
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 }f/ 1
;plot add fap black outside on 0p)#!$
;plot add cont zdis out on 1=}+NK!
;plot add axes black ( hp 52Vse
set large :H@Q`g u
F8S% \i
Il `35~a
sel pile begin=(20.3,5.0,15.0) end=(20.3,5.0,3.0) nseg 12 pxDkf|*
krFuEaO
!g'kWE[
sel pile prop emod 1.0e1 nu 0.3 xcarea 0.7854 xcj 9.82e-2 xciy 4.91e-2 xciz 4.91e-2 & i]Njn k
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 il|e5TD^
j KGfm9|zj
iz.J._&
def place_piles 'pa[z5{k+
loop n(1,10) j]Kpwf<NS
x_d=20.0+10*cos(12*pi/180) jb|al[p\
y_d=1.1*(n-1)
LhKaqR{
z_1=14.5-10*sin(12*pi/180) 5<^$9('
z_2=13.5-10*sin(12*pi/180) dG2k4 O
z_3=12.5-10*sin(12*pi/180) *eK\W00
z_4=11.5-10*sin(12*pi/180) &e0BL z
z_5=10.5-10*sin(12*pi/180) db -h=L|
p_1=5*n-4 >g+yw1nC
p_2=5*n-3 fh,kbn==r?
p_3=5*n-2 _)XQb1]
p_4=5*n-1 oR}cE
Sr
p_5=5*n B:J([@\'
command L7VD ZCV
sel cable id=@p_1 begin=(20.0,@y_d,14.5) end=(@x_d,@y_d,@z_1) nseg=9 n@[_lNa4GD
sel cable id=@p_2 begin=(20.0,@y_d,13.5) end=(@x_d,@y_d,@z_2) nseg=9 LDc?/
Z1
sel cable id=@p_3 begin=(20.0,@y_d,12.5) end=(@x_d,@y_d,@z_3) nseg=9 D<U^FT
sel cable id=@p_4 begin=(20.0,@y_d,11.5) end=(@x_d,@y_d,@z_4) nseg=9 M#Kke9%2
sel cable id=@p_5 begin=(20.0,@y_d,10.5) end=(@x_d,@y_d,@z_5) nseg=9 nvY%{Zf$}
sel cable prop emod 2.0e11 gr_coh 4.5e4 gr_fric 25 gr_k 4.2e7 xcarea 7.87e-3 ytens 3.0e8 _UI*W&*
end_command G9 v'a&
end_loop 1cq"H/N
end RUco3fZ
@place_piles
!?)ky `S3
;****************************************** ,a
2(h
hist nstep 20 1 etl:gcEC
hist gp zdis 20.0 5.0 15.0 ![1+=F!
hist gp zdis 20.5 5.0 15.0 kD;1+lNz
hist gp zdis 21.0 5.0 15.0 M~+}ss
hist gp zdis 21.5 5.0 15.0 =>3wI'I
hist gp zdis 22.0 5.0 15.0 IyK^` y
hist gp zdis 22.5 5.0 15.0 \no6]xN;
hist gp zdis 23.0 5.0 15.0 QiKci%=SX
hist gp zdis 23.5 5.0 15.0 BA(erf>
hist gp zdis 24.0 5.0 15.0 Z uFk}R"x
hist gp zdis 24.5 5.0 15.0 H?oBax:
hist gp zdis 25.0 5.0 15.0 Z mVw5G
q
hist gp zdis 25.5 5.0 15.0 C-h9_<AwJQ
hist gp zdis 26.0 5.0 15.0 v'RpsCov
hist gp zdis 26.5 5.0 15.0 .bYZkO:oy
hist gp zdis 27.0 5.0 15.0 7XzhKA6
hist gp zdis 27.5 5.0 15.0 27J!oin$
hist gp zdis 28.0 5.0 15.0 C fs2tN
hist gp zdis 28.5 5.0 15.0 Y=83r]%
hist gp zdis 29.0 5.0 15.0 yM}~]aQ y
hist gp zdis 39.5 5.0 15.0 aQ.QkMZ
hist gp zdis 30.0 5.0 15.0 G.")Bg
hist gp zdis 31.0 5.0 15.0 Z=ZTSl
hist gp zdis 32.0 5.0 15.0 Lh@0|k
hist gp zdis 33.0 5.0 15.0 YN]xI
hist gp zdis 34.0 5.0 15.0 I-Q@v`
hist gp zdis 35.0 5.0 15.0 H2kib4^i
hist gp zdis 36.0 5.0 15.0 `D-P}hDm!
hist gp zdis 37.0 5.0 15.0 0Wkk$0h9
hist gp zdis 38.0 5.0 15.0 ]y0bgKTK
hist gp zdis 39.0 5.0 15.0 *qa.hqas
hist gp zdis 40.0 5.0 15.0 Sy@)Q[A
hist gp zdis 42.0 5.0 15.0 Nwwn #+
hist gp zdis 44.0 5.0 15.0 [x%[N)U3
hist gp zdis 46.0 5.0 15.0 =y^`yv 3
hist gp zdis 48.0 5.0 15.0 2]of4
hist gp zdis 50.0 5.0 15.0 \+M6R<Qw
hist gp zdis 52.0 5.0 15.0 *c<6 Er>s
hist gp zdis 54.0 5.0 15.0 zsM2R"[X
hist gp zdis 56.0 5.0 15.0 |E?PQ?P
hist gp zdis 58.0 5.0 15.0 dvAG}<
hist gp zdis 60.0 5.0 15.0 \+)aYP2Hu
;************************************************** yopC
<k
hist gp xdis 20.0 5.0 15.0 5bB\i79$
hist gp xdis 20.0 5.0 14.5 I=[Ir8};
hist gp xdis 20.0 5.0 14.0 9| g]M:{
hist gp xdis 20.0 5.0 13.5 ObPXVqG"?
hist gp xdis 20.0 5.0 13.0 o;TS69|D
hist gp xdis 20.0 5.0 12.5 `.z"Q%uz
hist gp xdis 20.0 5.0 12.0 Y1Q240
hist gp xdis 20.0 5.0 11.5 LwlO)|E
hist gp xdis 20.0 5.0 11.0 oyT`AYa
hist gp xdis 20.0 5.0 10.5 x{O) n
hist gp xdis 20.0 5.0 10.0 !F:ANoaS
hist gp xdis 20.0 5.0 9.5 Rtu"#XcBw+
hist gp xdis 20.0 5.0 9.0 @
D,]v:
hist gp xdis 20.0 5.0 8.5 v!`:{)2C
hist gp xdis 20.0 5.0 8.0 &HQ_e$1
hist gp xdis 20.0 5.0 7.5 T>hrKn.!D:
hist gp xdis 20.0 5.0 7.0 a,|Hn
hist gp xdis 20.0 5.0 6.5 G! Y
l0Zr
hist gp xdis 20.0 5.0 6.0 R$ra=sL`
hist gp xdis 20.0 5.0 5.5 S,Z~-j
hist gp xdis 20.0 5.0 5.0 29AE B
hist gp xdis 20.0 5.0 4.5 nR]*RIp5
hist gp xdis 20.0 5.0 4.0 38 ]}+Bb
hist gp xdis 20.0 5.0 3.5 ;kI)j
?
hist gp xdis 20.0 5.0 3.0 AQ&vq$
;*************************************** 8NkyT_\
hist gp zdis 0.0 5.0 10.0 b>-h4{B[
hist gp zdis 1.0 5.0 10.0 v M'!WVs
hist gp zdis 2.0 5.0 10.0 7qSnP30}
hist gp zdis 3.0 5.0 10.0 a&Qr7tTY"
hist gp zdis 4.0 5.0 10.0 Z9G4in8
hist gp zdis 5.0 5.0 10.0 Kf-XL),3l
hist gp zdis 6.0 5.0 10.0 %4wEAi$I
hist gp zdis 7.0 5.0 10.0 \G#_z|'dN
hist gp zdis 8.0 5.0 10.0 ?y-@c]
hist gp zdis 9.0 5.0 10.0 GrUpATIx
hist gp zdis 10.0 5.0 10.0 -5oYGLS$y3
hist gp zdis 11.0 5.0 10.0 v2IcDz`}7
hist gp zdis 12.0 5.0 10.0 (mr*Thy`@
hist gp zdis 13.0 5.0 10.0 l 88n*O
hist gp zdis 14.0 5.0 10.0 QQ`tSYgex
hist gp zdis 15.0 5.0 10.0 9u)p9)^-.v
hist gp zdis 16.0 5.0 10.0 6FYL},.R
hist gp zdis 17.0 5.0 10.0 IL`LIJ:O
hist gp zdis 18.0 5.0 10.0 S{RRlR6Z
hist gp zdis 19.0 5.0 10.0 !1ZItJ74#
hist gp zdis 20.0 5.0 10.0 / Xq|SO
;*********************************** w&@zJ [
hist sel pilesel force fx (20.3,5.0,15.0) mGDy3R90
hist sel pilesel force fx (20.3,5.0,14.0) Egz6rRCvg
hist sel pilesel force fx (20.3,5.0,13.0) R4X9g\KpAt
hist sel pilesel force fx (20.3,5.0,12.0) (*M0'5
hist sel pilesel force fx (20.3,5.0,11.0) x<`^4|<
hist sel pilesel force fx (20.3,5.0,10.0) Y<IuwS
hist sel pilesel force fx (20.3,5.0,9.0) rd>>=~vx=/
hist sel pilesel force fx (20.3,5.0,8.0) \2!.
hist sel pilesel force fx (20.3,5.0,7.0) `YUeVz>q?
hist sel pilesel force fx (20.3,5.0,6.0) ]x).C[^
hist sel pilesel force fx (20.3,5.0,5.0) %"cOX
hist sel pilesel force fx (20.3,5.0,4.0) =),O ;M
hist sel pilesel force fx (20.3,5.0,3.0) %K+hG=3O
hist sel pilesel moment mx (20.3,5.0,15.0) JMk2OK{0
hist sel pilesel moment mx (20.3,5.0,14.0) NT2XG&$W>
hist sel pilesel moment mx (20.3,5.0,13.0) 4X",:B}
hist sel pilesel moment mx (20.3,5.0,12.0) :_<&LO]Q
hist sel pilesel moment mx (20.3,5.0,11.0) q44vI
hist sel pilesel moment mx (20.3,5.0,10.0) /0cm7[a ?
hist sel pilesel moment mx (20.3,5.0,9.0) S|xwYaoy%
hist sel pilesel moment mx (20.3,5.0,8.0) 9B![l=Gh
hist sel pilesel moment mx (20.3,5.0,7.0)
lI%RdA[
hist sel pilesel moment mx (20.3,5.0,6.0) ,na}' A@a`
hist sel pilesel moment mx (20.3,5.0,5.0) Rp;"]Q&b
hist sel pilesel moment mx (20.3,5.0,4.0) |)xWQ KzA
hist sel pilesel moment mx (20.3,5.0,3.0) $}IG+,L
model null range x 0.0 20.0 z 15.0 13.0 ju5o).!bg
step 200 +1I7K|M
model null range x 0.0 20.0 z 13.0 12.0 {Mj- $G"
step 200 8bQXC+bK
model null range x 0.0 20.0 z 12.0 11.0 3-E-\5I
step 200 \crh`~?>
model null range x 0.0 20.0 z 11.0 10.0 #eN2{G=4+
G)^/#d#&
!4`:(G59
step 15000 !4!qHJISa
hist write -2 -3 -4 -5 -6 -7 -8 -9 -10 -11 -12 -13 -14 -15 -16 -17 -18 -19 -20 -21 -22 -23 & !Pz#czo
-24 -25 -26 -27 -28 -29 -30 -31 -32 -33 -34 -35 -36 -37 -38 -39 -40 -41 vs 1 file dingbushuxiangweiyi.dat iITMBS`}
hist write -42 -43 -44 -45 -46 -47 -48 -49 -50 -51 -52 -53 -54 -55 -56 -57 -58 -59 -60 -61 & *P5\T4!+d
-62 -63 -64 -65 -66 vs 1 file qiangticexiangweiyi.dat m]*a;a'}#
hist write -67 -68 -69 -70 -71 -72 -73 -74 -75 -76 -77 -78 -79 -80 -81 -82 -83 -84 -85 -86 & +8W5amk.P|
-87 vs 1 file kengdilongqiweiyi.dat :Tv>)N
hist write -88 -89 -90 -91 -92 -93 -94 -95 -96 -97 -98 -99 -100 -101 -102 -103 -104 -105 & H)h^|A/vO
-106 -107 -108 -109 -110 -111 -112 -113 vs 1 file zhuangtibianxing.dat lxCAZa\
VxS3lR=
,}khu
save result.sav