Q5u3~Q'e
new =vDpm,
config dyn F!2VTPm9z
set dyn off JeAyT48!M
;=============================================== : qV|rih_Q
;建模型 uC6 e2py<[
;=============================================== AgB$
w4
def setup ;定义函数setup =$ubSfx
yc=1 ;y方向尺寸 ju4wU; Nu
czz1=-13.16 +vPCr&40
czz11=1.9 ;四分之一车站z方向尺寸 " Lh XR
czz12=3 ;四分之一车站z方向尺寸 ^K
9jJS9K
czz13=6.8 ;四分之一车站z方向尺寸 Ye^ xV,U@
czz14=1.46 ;四分之一车站z方向尺寸 m`b:#z
ttx1=10.05 ;四分之一车站外围土体x方向尺寸 NH1ak(zHW
ttx2=-9.55 ;四分之一车站外围土体-x方向尺寸 2uMSeSx$
ht1=-3.1 A2I qn5
ht2=3.1 . TNJuuO
tt2=68.9 ;78.95 jO`L:D/C
tt1=-68.9 E~8J<g E
xh=33.74 ；车站以下土体高度 KqNsCT+j
g Eq6[G
end nE : Wl
setup .
ywVGBvJ
;------------------------------------------------------------------------------------------------- e'=#G$S?g
; 一、生成右侧土体 >R.!Qze\G
;------------------------------------------------------------------------------------------------- cCIs~*D
def gen_Sando1 akwS;| SZ
command u
N_< G
gen zone brick p0 0,0,ht1 p1 add ttx1,0,0 p2 add 0,yc,0 p3 add 0,0,ht2 & 0 4oMgH>Vd
p4 add ttx1,yc,0 p5 add 0,yc,ht2 p6 add ttx1,0,ht2 p7 add ttx1,yc,ht2 & p vRa
size 10 1 3 ;ratio 1 1 1 1.2 Ppx 4#j
.tK]-f2
end_command :z$+leNH \
end L_!} R
gen_Sando1 qVd s
2
group sand_o1 range x 0 10.05 y 0 1 z -3.1 0 _cJ\A0h^
;------------------------------------------------------------------------------------------------- t3!~= U
def gen_Sando1 ("=24R=a
command ^HHT>K-m
gen zone brick p0 ttx1,0,ht1 p1 add tt2,0,0 p2 add 0,yc,0 p3 add 0,0,ht2 & D@k#'KU
p4 add tt2,yc,0 p5 add 0,yc,ht2 p6 add tt2,0,ht2 p7 add tt2,yc,ht2 & s##XC^; p[
size 50 1 3 ;ratio 1.1 1 1 4 X0ku]
,{Z!T5 |
end_command /EL3 T t
end c{jTCkzq
gen_Sando1 e C&!yY2g
group sand_o1 range x 10.05 78.95 y 0 1 z -3.1 0 Owh:(EJ"d
;------------------------------------------------------------------------------------------------- lW]&a"1$
def gen_Sando2 T3-/+4$0v
command )=0@4
gen zone brick p0 ttx1,0,-5 p1 add tt2,0,0 p2 add 0,yc,0 p3 add 0,0,czz11 & qf%p#+: B3
p4 add tt2,yc,0 p5 add 0,yc,czz11 p6 add tt2,0,czz11 p7 add tt2,yc,czz11 & 5L\Im^
size 50 1 2 ;ratio 1.1 1 1 U {HBmSR
!ed0
end_command B5]nP .R
end B< ,AI7
gen_Sando2 jfZ)
group sand_o2 range x 10.05 78.95 y 0 1 z -5 -3.1 ;11.26 13.16 q c6d,z/
^5-SL?E
sT91>'&
;------------------------------------------------------------------------------------------------- <In+V
;------------------------------------------------------------------------------------------------- FXid=&T@0D
def gen_Sando3 quFNPdP
command A .<M*[{q
gen zone brick p0 ttx1,0,-8 p1 add tt2,0,0 p2 add 0,yc,0 p3 add 0,0,czz12 & 5"Y:^_8
p4 add tt2,yc,0 p5 add 0,yc,czz12 p6 add tt2,0,czz12 p7 add tt2,yc,czz12 & 0'R}'
size 50 1 3 ;ratio 1.1 1 1 YRj"]=
5N
P_M!h~
end_command ) =|8%IrB
end @%6"xnb `
gen_Sando3 vGp`P
group sand_o3 range x 10.05 78.95 y 0 1 z -8 -5 ;8.26 11.26 nB%[\LtZ?
;------------------------------------------------------------------------------------------------- $ u,`bX
def gen_Sando4 Kq:vTz&<
command _147d5
gen zone brick p0 ttx1,0,-14.8 p1 add tt2,0,0 p2 add 0,yc,0 p3 add 0,0,czz13 & Wk]E6yz6
p4 add tt2,yc,0 p5 add 0,yc,czz13 p6 add tt2,0,czz13 p7 add tt2,yc,czz13 & fc%C!^7
size 50 1 6 ;ratio 1.1 1 1 Bo/i =/7%
wd
Di5-A4
end_command Ia>~ph#]{`
end W1O Y}2kj
gen_Sando4 +qyx3c+
group sand_o4 range x 10.05 78.95 y 0 1 z -14.8 -8 ;1.46 8.26 ^]$rh.7&
Y,X0x-
;------------------------------------------------------------------------------------------------- eak+8URo
def gen_Sando5 P)UpUMt;k
command 'Y>@t6E4
gen zone brick p0 ttx1,0,-16.26 p1 add tt2,0,0 p2 add 0,yc,0 p3 add 0,0,czz14 & qkq^oHI
p4 add tt2,yc,0 p5 add 0,yc,czz14 p6 add tt2,0,czz14 p7 add tt2,yc,czz14 & /qXP \ a
size 50 1 2 ;ratio 1.1 1 1 ,h.hgyt
L>,x G.oG
end_command :tLMh08h
end v0+$d\mP4<
gen_Sando5 p[D,.0SuC
group sand_o5 range x 10.05 78.95 y 0 1 z -16.26 -14.8 ;0 1.46 QH6_nZY
;------------------------------------------------------------------------------------------------- 8{p#Nl?U1
;------------------------------------------------------------------------------------------------- 5K$d4 KT
BU%gXr4Ra
def gen_Sando6 . Kk'N
command 7T=:dv
gen zone brick p0 0,0,-50 p1 add ttx1,0,0 p2 add 0,yc,0 p3 add 0,0,xh & *G M.2``e
p4 add ttx1,yc,0 p5 add 0,yc,xh p6 add ttx1,0,xh p7 add ttx1,yc,xh &
C0j`H(
size 10 1 20 ;ratio 1.2 1 1 M V~3~h8
n*N`].r#{=
end_command CSMx]jbb
end t` zPx#])
gen_Sando6 =|q@ Q`DB
group sand_o6 range x 0 10.05 y 0 1 z -50 -16.26 qApf\o3[0
;------------------------------------------------------------------------------------------------- us^J!
s7
def gen_Sando6 4% 2MY\
command 'k9 1;T[
gen zone brick p0 ttx1,0,-50 p1 add tt2,0,0 p2 add 0,yc,0 p3 add 0,0,xh & (EOYJHZB!
p4 add tt2,yc,0 p5 add 0,yc,xh p6 add tt2,0,xh p7 add tt2,yc,xh & 0u ,nSvch
size 40 1 20 ;ratio 1.1 1 EBplr ,
2sT\+C&H
end_command BE," lX
end H`8}w{ ft&
gen_Sando6 w'P !<JaZ
group sand_o6 range x 10.05 78.95 y 0 1 z -50 -16.26 ;-33.74 0 ?) VBkA5j
;------------------------------------------------------------------------------------------------- nnMRp7LQ-
;------------------------------------------------------------------------------------------------- `Hs I)RmX
; 二、生成左侧土体 &fgfCZz'
;------------------------------------------------------------------------------------------------- :-1
i1d
;------------------------------------------------------------------------------------------------- FeuqqZ\=&
def gen_Sand11 n7'X.=o7
command Bfr$&?j#
gen zone brick p0 0,0,ht1 p1 add 0,yc,0 p2 add ttx2,0,0 p3 add 0,0,ht2 & oF(|NS^
p4 add ttx2,yc,0 p5 add ttx2,0,ht2 p6 add 0,yc,ht2 p7 add ttx2,yc,ht2 & E% \Ohs7
size 1 10 3 ;ratio 1 1 1 1.2 rFYw6&;vOi
a & tl@y1
end_command -<aN$O
end $U}GX'1LZ
gen_Sand11 EY[J;H_b
group sand_o7 range x -9.55 0 y 0 1 z -3.1 0 7bx!A+, t
;------------------------------------------------------------------------------------------------- X h"8uJD
def gen_Sand2 z4*`K4W
command lag%} ^
gen zone brick p0 ttx2,0,ht1 p1 add 0,yc,0 p2 add tt1,0,0 p3 add 0,0,ht2 & N(L?F):fT
p4 add tt1,yc,0 p5 add tt1,0,ht2 p6 add 0,yc,ht2 p7 add tt1,yc,ht2 & X##hSGQM
size 1 50 3 ;ratio 1 1.1 1 UWq[K&vQZ
$r mfE
end_command L2^M#G@t
end DsF<P@O 6
gen_Sand2 WhN~R[LE_
group sand_o7 range x -78.45 -9.55 y 0 1 z -3.1 0 m}?(c)ST
;------------------------------------------------------------------------------------------------- yX|0 R
H
def gen_Sand31 :LBRyBV
command 5<L+T
gen zone brick p0 ttx2,0,-5 p1 add 0,yc,0 p2 add tt1,0,0 p3 add 0,0,czz11 & y4 ~;H{!
p4 add tt1,yc,0 p5 add tt1,0,czz11 p6 add 0,yc,czz11 p7 add tt1,yc,czz11 & S(h+,+289
size 1 50 2 ;ratio 1 1.1 1 j43-YdCJ
D$G:#z*
end_command = $awUy
end hha^ :,
gen_Sand31 84u %_4/
group sand_o8 range x -78.45 -9.55 y 0 1 z -5 -3.1 ;11.26 13.16 cSs??i
D"q
;------------------------------------------------------------------------------------------------- K
na
def gen_Sand32 ld/\`s[i
command :I^I=A%Pe(
gen zone brick p0 ttx2,0,-8 p1 add 0,yc,0 p2 add tt1,0,0 p3 add 0,0,czz12 & k\wI^D
p4 add tt1,yc,0 p5 add tt1,0,czz12 p6 add 0,yc,czz12 p7 add tt1,yc,czz12 &
oN7JNMT
size 1 50 3 ;ratio 1 1.1 1 P|4qbm4%O,
#&ZwQw
end_command qx<h rC0Z&
end %,[p[`NRYR
gen_Sand32 !_My]>S
group sand_o9 range x -78.45 -9.55 y 0 1 z -8 -5 ;8.26 11.26 zeG_H} [2&
;------------------------------------------------------------------------------------------------- V50FX }i
def gen_Sand33 !`?*zf
command kA% "-$3
gen zone brick p0 ttx2,0,-14.8 p1 add 0,yc,0 p2 add tt1,0,0 p3 add 0,0,czz13 & "] 0sR
p4 add tt1,yc,0 p5 add tt1,0,czz13 p6 add 0,yc,czz13 p7 add tt1,yc,czz13 & $M 1 /74
size 1 50 6 ;ratio 1 1.1 1 *FrlzIAo m
XhEd9>#
end_command Y1Qg|U o
end h#!u" 'JW
gen_Sand33 O+Q t8,
group sand_o10 range x -78.45 -9.55 y 0 1 z -14.8 -8 ;1.46 8.26 LU4 k/
Auz.we s
;------------------------------------------------------------------------------------------------- aL%amL6CX
def gen_Sand34 ?Wz
rv&E2
command j[< }l&
gen zone brick p0 ttx2,0,-16.26 p1 add 0,yc,0 p2 add tt1,0,0 p3 add 0,0,czz14 & ( :2,Rr1"
p4 add tt1,yc,0 p5 add tt1,0,czz14 p6 add 0,yc,czz14 p7 add tt1,yc,czz14 & N]6M4j!
size 1 50 2 ;ratio 1 1.1 1 1rmK#ld"=Z
,/>hWAx
end_command yC]X&1,:z
end {@8TGHKv
gen_Sand34 ">^O{X\
group sand_o11 range x -78.45 -9.55 y 0 1 z -16.26 -14.8 ;0 1.46 4#^?- 6
;------------------------------------------------------------------------------------------------- lYq4f|5H}m
;------------------------------------------------------------------------------------------------- Z U^dLN- N
< _~>YJ
def gen_Sand4 1EvK\
command ) ]x/3J@
gen zone brick p0 0,0,-50 p1 add 0,yc,0 p2 add ttx2,0,0 p3 add 0,0,xh & * vMNv
p4 add ttx2,yc,0 p5 add ttx2,0,xh p6 add 0,yc,xh p7 add ttx2,yc,xh & 3A(sT}
size 1 10 20 42wa9UL<Ka
Y }pCBw
end_command s]z-d! G
end HX77XT y
gen_Sand4 6*n<emP
group sand_o12 range x -9.55 0 y 0 1 z -16.26 -50 tuF
hPqe {
;------------------------------------------------------------------------------------------------- -0/5 !
def gen_Sand5 DZI:zsf;5Q
command toCxY+"nbU
gen zone brick p0 ttx2,0,-50 p1 add 0,yc,0 p2 add tt1,0,0 p3 add 0,0,xh & =6sL }$
p4 add tt1,yc,0 p5 add tt1,0,xh p6 add 0,yc,xh p7 add tt1,yc,xh & ,>rr |O
size 1 50 20 ;ratio 1 1.1 1 c{dge/2yb
I*Q^$YnM
end_command XJG "Zr 9
end " +6:vhP5
gen_Sand5 ;/Z9M"!u[
group sand_o12 range x -78.45 -9.55 y 0 1 z -16.26 -50 ;-33.74 0 (79y!&9p
;------------------------------------------------------------------------------------------------- 1k
*gbXb
call Flacmodel_haitang.dat 1dy>a=W
;------------------------------------------------------------------------------------------------- ?Q ]{P]
plot add axes 6:v$g
;-------------------------------------------------------------------------------------- !5;A.f
;二、生成接触面的命令 -!k$ Z
;-------------------------------------------------------------------------------------- (A\p5@ht
interface 1 face range x -9.55 10.05 y 0 1 z -16.27 -16.25 K2u $1OKv
interface 2 face range x -9.55 10.05 y 0 1 z -3.11 -3.09 ;QA`2$Ow
interface 3 face range x -9.56 -9.54 y 0 1 z -5 -3.1 eKq`t.*Ft
interface 4 face range x -9.56 -9.54 y 0 1 z -8 -5 F( w
interface 5 face range x -9.56 -9.54 y 0 1 z -14.8 -8 lbCTc,xT
interface 6 face range x -9.56 -9.54 y 0 1 z -16.26 -14.8 ?x|8"*N
vIi#M 0@N
int 1 prop kn 1.570e9 ks 1.570e9 fric 10 coh 10000 ;prop属性 kn法向刚度 ks剪切刚度 fric为摩擦角 coh内聚力 JToc("V
int 2 prop kn 1.570e9 ks 1.570e9 fric 15 coh 55000 v~2$ 9x!9
int 3 prop kn 1.570e9 ks 1.570e9 fric 26 coh 8000 .9< i
int 4 prop kn 1.570e9 ks 1.570e9 fric 15 coh 55000 +x3T^G
int 5 prop kn 1.570e9 ks 1.570e9 fric 14 coh 24000 KfO$bmwmx
int 6 prop kn 1.570e9 ks 1.570e9 fric 17 coh 56000 %$)[ qa3
ini y add -2.0 rang group 1 * P#okwp
;=============================================== 9bL`0 L
;初始应力计算程序 %)Pn<! L
;=============================================== Aqwjs
3
;定义模型及性质 =X0"!y"
model mohr rang group sand_o1 i9q n_/<c
prop bulk 1.875e7 shear 0.967e7 coh 10000 fric 10 rang group sand_o1 us2X:X)
model mohr rang group sand_o2 CO` %eL ~
prop bulk 1.875e7 shear 0.967e7 coh 55000 fric 15 rang group sand_o2 2&f] v`|M|
model mohr rang group sand_o3 G8y:f%I!b
prop bulk 2.068e7 shear 0.793e7 coh 8000 fric 26 rang group sand_o3 t_-1sWeA!
model mohr rang group sand_o4 L@"1d.k_
prop bulk 1.875e7 shear 0.967e7 coh 55000 fric 15 rang group sand_o4 Yy$GfjJtL]
model mohr rang group sand_o5 TfD]`v`]
prop bulk 1.396e7 shear 0.758e7 coh 24000 fric 14 rang group sand_o5 LG0z|x(
model mohr rang group sand_o6 /$
-^k[%
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o6 m@r+M"!R
model mohr rang group sand_o7 x2|YrkGv
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o7 '8Q:}{
model mohr rang group sand_o8 f
xWW "B*A
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o8 vgc #IEx@
model mohr rang group sand_o9 Kd?TIeF E
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o9 dU7+rc2,CU
model mohr rang group sand_o10 ],lrT0_cT
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o10 ;^q@w
model mohr rang group sand_o11 fgs){ Ng`
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o11 TV~ <1vj
model mohr rang group sand_o12 1e 8J-Nkj
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o12 s<i& q {r
^ d i[J^
mo elas rang group 1 _%M5
T
;prop bulk 20.8e9 shear 15.58e9 rang group 1 ;正确 d+1q[,-
prop bulk 22.8e9 shear 17.08e9 rang group 1 ;只把柱子等效 X*'tJN $
ini dens 2500 rang group 1 om`x"x&6
ini dens 1900 rang group sand_o1 ;18000 I. [2 -~yf
ini dens 1990 rang group sand_o2 \"]vSx>
ini dens 1910 rang group sand_o3 +bG^SH2ke
ini dens 1990 rang group sand_o4 YceX)
ini dens 1970 rang group sand_o5 tSr.0'CE
ini dens 2010 rang group sand_o6 6;02_C]\o
ini dens 2010 rang group sand_o7 l ( EDe
ini dens 2010 rang group sand_o8 "k)}qI{
ini dens 2010 rang group sand_o9 ~ nQv
yM!$
ini dens 2010 rang group sand_o10 gEVN;G'B<=
ini dens 2010 rang group sand_o11 }tvLe3O
ini dens 2010 rang group sand_o12 )jM'
x&Vg
;边界条件 1@&i
ju5
fix x y z rang z -50.1 -49.9 @a08*"lbp
fix x rang x 78.85 79.05 PdUlwT? 8C
fix x rang x -78.55 -78.35 [jPUAr}
fix y rang y -.1 .1 0Q81$% @<
fix y rang y 0.9 1.1 dM%#DN8 l
;应力条件 4z^ ?3@:K
ini szz 0 grad 0 0 20e3 rang z -50 0 ;加权平均重度为1993 VJK?"m X
ini sxx 0 grad 0 0 10e3 rang z -50 0 K3uNR w
ini syy 0 grad 0 0 10e3 rang z -50 0 .?L&k|wX-
;--------------------------------------------------------------------------------- Uxla,CCp-
hist id 1 unbal cs]N%M^s
plot hist 1 ~uF%*
; ]SqLF!S (=
set grav 0 0 -10 'y(;: Kc
solve ratio 1.0e-4 9 ~bl
save station1.sav
