命令流: dr^pzM!N
o4o&}
new .x$+R%5U
def _parm 4pV.R5:
`B) ~
;-----------------几何参数---------------------------------------------------- {5-4^|!
X3(:)zUL
p_d1=0.15 ;-----------桩径------------------------------- 8wwqV{O7
tslope=7.5 ;-----------坡顶------------------------------- /u!I2DF
dslope=-7.5 ;-----------坡脚------------------------------- '[I_Iu#,
ngrad01=1.5 ;-----------坡率------------------------------- [0yKd?e
s1=3.0*p_d1 ;-----------桩排距S1=3D------------------------- .a:"B\B`
s2=3.0*p_d1 ;-----------桩间距S2=3D---------------------- qWKpnofa
s=2.0*s2 ;-----------模型厚度S=2S2---------------------- u@}((V
sl1=s-0.1 Wd~aSz9
sl2=s+0.1 dn$1OhN8M
yY4*/w7*j4
;-----------------土体材料参数---------------------------------------------------- "^ydoRZ
)^
R]3!v
s_Emod=1.0e7 $Eg|Qc-1
s_nu=0.3 'Z+w\0}@
s_fri=20.0 3qiE#+dC
s_coh=10e3 ~;aSX1
s_tens=0 ;{ XKZ}
soil_dens=2000 T2Z;)e$m_
soil_bulk=s_Emod/3/(1-2*s_nu) O;ZU{VY
soil_shear=s_Emod/2/(1+s_nu) ,G
e7
9(
Tc,Bv7:
7S7gU\qOj
;-----------------桩位---------------------------------------------------- m'KY;C
jiYYDGs77
;-----------------1# 桩---------------------- bRJYw6oA<
x_beg01=-s1/2 _2q4Aaza
y_beg01=0.5*s2 "BK'<j^q
z_beg01=10.0+5.0*(7.5+x_beg01)/7.5 \2U^y4K.
YM
0f_G=
x_end01=-s1/2 Pr{? A]dQ
y_end01=0.5*s2 Ljq/f&
c
z_end01=0 g[@Kd
C#emmg!a\
;-----------------2# 桩---------------------- -. L)-%wIV
x_beg02=-s1/2 _Hl[Fit<j1
y_beg02=1.5*s2 {0q;:7Bt
z_beg02=10.0+5.0*(7.5+x_beg01)/7.5 a~Sf~ka
R5ZnkPEA
x_end02=-s1/2 sVZb[|zSri
y_end02=1.5*s2 R/waWz\D
z_end02=0 hc5M)0d
jx!)N>
;-----------------3# 桩---------------------- ?1?^>M
x_beg03=s1/2 Ra'0 ^4t
y_beg03=0.5*s2 A)2vjM9}K
z_beg03=10.0+10.0*(7.5+x_beg03)/15 +)^F9LPl
+[UFf3(ON
x_end03=s1/2 VFT
G3,kI
y_end03=0.5*s2 bEV
9l
z_end03=0 [VT&
<P3r+ 1|R
;-----------------4# 桩---------------------- >>rW-&
x_beg04=s1/2 miCt)Qd
y_beg04=1.5*s2 WiH%URFB
z_beg04=10.0+10.0*(7.5+x_beg03)/15 VjtI1I
02Ur'|
x_end04=s1/2 R@5jEf
y_end04=1.5*s2 ::g"dRS<v
z_end04=0 ~cz]Rhq
)pjd*+V
_J6|ju\
end {q~Bss{z
_parm 8X!UtHml
F:%= u
=
;-----------------建立模型---------------------------------------------------- $yIcut7
U
u(ysN4`
gen zone brick p0 -17.5 0 0 p1 -7.5 0 0 p2 -17.5 s 0 p3 -17.5 0 10 size 20 3 20 ratio 1 1 1 group 21 $JypVA(CX
gen zone brick p0 7.5 0 0 p1 17.5 0 0 p2 7.5 s 0 p3 7.5 0 20 size 20 3 20 ratio 1 1 1 group 22 W8,XSUl
;DL|%-%;$r
gen zone brick p0 dslope 0 0 p1 tslope 0 0 p2 dslope s 0 p3 dslope 0 10 p4 tslope s 0 p5 dslope s 10 p6 tslope 0 20.0 p7 tslope s 20.0 size 20 3 20 ratio 1 1 1 group 12 )-Zpr1kD
}Jk.c~P)
group slope range group 12 any group 21 any group 22 any P!,\V\TY]
attach face >-r\]/^
7I/Sfmqy"O
model mohr ,ZJI]Q=!
prop shear soil_shear bulk soil_bulk fric s_fri coh s_coh tens s_tens CM>/b3nOW
ini dens soil_dens 4))u*c/,
^TyusfOz
def pileproperty qBpv[m
d_pile=p_d1 y%|nE((
A_pile=pi*(d_pile)^2/4 ;截面积 !aeL*`;
Iy_pile=pi*(d_pile)^4/64 ;Y轴惯性矩 n]}+ :
Iz_pile=pi*(d_pile)^4/64 ;Z轴惯性矩 V,<,;d fR
Ij_pile=Iy_pile+Iz_pile ;惯性积 0Tp?ED_
C_pile=pi*d_pile ;周长 /SY40;k:
emod_pile=2.5e10 ;微型桩弹性模量 $N5VoK
nu_pile=0.2 ;微型桩泊松比 R$K.;
end pileproperty M/5/Tp
pileproperty K)]7e?:Wu
%/oOM\}++
sel pile id=1 begin=(x_beg01,y_beg01,z_beg01) end=(x_end01,y_end01,z_end01) nseg=15 +o?.<[>!GR
sel pile id=1 begin=(x_beg02,y_beg02,z_beg02) end=(x_end02,y_end02,z_end02) nseg=15 82~UI'f \
sel pile id=1 begin=(x_beg03,y_beg03,z_beg03) end=(x_end03,y_end03,z_end03) nseg=15 8d-; ;V
sel pile id=1 begin=(x_beg04,y_beg04,z_beg04) end=(x_end04,y_end04,z_end04) nseg=15 0'TqW9P
&"._%S58V
sel pile id=1 prop emod=emod_pile nu=nu_pile & .i"v([eQ
xcarea=A_pile per=C_pile & Z9i,#/
xcj=Ij_pile xciy=Iy_pile xciz=Iz_pile & eE>3=1d]w
cs_sk=1.3e11 cs_scoh=1.5e10 cs_sfric=20 & ?#qA>:2,
cs_nk=1.3e9 cs_ncoh=1.0e4 cs_nfric=20 cs_ngap=off z<+".sD'
K-)*S\<}
;-------------------------边界条件-------------------------------- .hg<\-:_
"}\2zub9
fix x range x -17.6 -17.4 *wwLhweQ5W
fix x range x 17.4 17.6 ]"dZE2!
fix y range y -0.1 0.1 H #E
fix y range y sl1 sl2 5u:+hB
fix x y z range z -0.1 0.1 pCb@4nb
`V w9j,G
;----------------------求解安全系数FOS----------------------- k1L GT&
1i_~ZzX8
set grav 0,0,-10 @!Il!+^3
en=Z[ZIPO
solve fos file 0.5l.sav