_^SNI ~
new \ 9 ap$
config dyn i~K~Czmok+
set dyn off 1
#EmZ{*
;=============================================== =q6yb@
;建模型 {`J7>K
;=============================================== bz}T}n j
def setup ;定义函数setup T \0e8"iZ
yc=1 ;y方向尺寸 ](T*f'LN
czz1=-13.16 q= 96Ci _a
czz11=1.9 ;四分之一车站z方向尺寸 A`OU} 'v?L
czz12=3 ;四分之一车站z方向尺寸 >TL0hBaaR
czz13=6.8 ;四分之一车站z方向尺寸 B'~.>, fg
czz14=1.46 ;四分之一车站z方向尺寸 [$] JvF
ttx1=10.05 ;四分之一车站外围土体x方向尺寸 /j S
ttx2=-9.55 ;四分之一车站外围土体-x方向尺寸 c&'T By
ht1=-3.1 . 5ingB3%
ht2=3.1 qP zgGbmD9
tt2=68.9 ;78.95 0Q>Yoa
11
tt1=-68.9 eRv3ZHH
xh=33.74 ；车站以下土体高度 (_-z m)F7
w l9E
end ,Hj=]e2?
setup 3!* J;Y
;------------------------------------------------------------------------------------------------- rJ<v1Yb
; 一、生成右侧土体 <Pf W
;------------------------------------------------------------------------------------------------- }1(F~6RH
def gen_Sando1 x>Q#B v y
command {},G xrQm
gen zone brick p0 0,0,ht1 p1 add ttx1,0,0 p2 add 0,yc,0 p3 add 0,0,ht2 & 1}V_:~7
p4 add ttx1,yc,0 p5 add 0,yc,ht2 p6 add ttx1,0,ht2 p7 add ttx1,yc,ht2 &
MNJ$/l)h
size 10 1 3 ;ratio 1 1 1 1.2 ?[VS0IBS
iCw~4KG
end_command 8G{} r
end @-u/('vpB
gen_Sando1 gy* N)iv%
group sand_o1 range x 0 10.05 y 0 1 z -3.1 0 \/lS!+~'']
;------------------------------------------------------------------------------------------------- [>6:xGSe9X
def gen_Sando1 8[zP2L!-
command nk6xavQji
gen zone brick p0 ttx1,0,ht1 p1 add tt2,0,0 p2 add 0,yc,0 p3 add 0,0,ht2 & DmD*,[rD
p4 add tt2,yc,0 p5 add 0,yc,ht2 p6 add tt2,0,ht2 p7 add tt2,yc,ht2 & $Cf_RFH0
size 50 1 3 ;ratio 1.1 1 1 >e Jk)qM
O{%y `|m
end_command 4}?Yp e-
end %JiA ,
gen_Sando1 <4( rY9
group sand_o1 range x 10.05 78.95 y 0 1 z -3.1 0 bh_i*DJ]
;------------------------------------------------------------------------------------------------- =zI
eZ 7
def gen_Sando2 5N '
QG<jE
command zXMIDrq
gen zone brick p0 ttx1,0,-5 p1 add tt2,0,0 p2 add 0,yc,0 p3 add 0,0,czz11 & !&19%C4
p4 add tt2,yc,0 p5 add 0,yc,czz11 p6 add tt2,0,czz11 p7 add tt2,yc,czz11 & YFsEuaV
size 50 1 2 ;ratio 1.1 1 1 % ]
4$vya+mAk5
end_command #PtV=Ee1
end 6AzH'H F
gen_Sando2 rEZa%)XJ
group sand_o2 range x 10.05 78.95 y 0 1 z -5 -3.1 ;11.26 13.16 l} h<2
RCgs3JIE+2
#msk'MVt
;------------------------------------------------------------------------------------------------- ^V> sNR
;------------------------------------------------------------------------------------------------- HH[ ?LKd<
def gen_Sando3 z!t &zkAK
command b/ur!2yr
gen zone brick p0 ttx1,0,-8 p1 add tt2,0,0 p2 add 0,yc,0 p3 add 0,0,czz12 & (wv DiW5
p4 add tt2,yc,0 p5 add 0,yc,czz12 p6 add tt2,0,czz12 p7 add tt2,yc,czz12 & e9N 1xB
size 50 1 3 ;ratio 1.1 1 1 5H? `a7q N
"1`c^
end_command &.Yh _
end | \Qr
c f
gen_Sando3 !y?g$e`
group sand_o3 range x 10.05 78.95 y 0 1 z -8 -5 ;8.26 11.26 R+, tn,<<
;------------------------------------------------------------------------------------------------- wCc:HfmjJ
def gen_Sando4 o),i2
command ,uuQj]Dac+
gen zone brick p0 ttx1,0,-14.8 p1 add tt2,0,0 p2 add 0,yc,0 p3 add 0,0,czz13 & >#c]rk:
p4 add tt2,yc,0 p5 add 0,yc,czz13 p6 add tt2,0,czz13 p7 add tt2,yc,czz13 & ,?i#NN5p
size 50 1 6 ;ratio 1.1 1 1 ^=Up U B
hC5iv J
end_command {f )" ,#
end DA[-(
s
gen_Sando4 hwko l W
group sand_o4 range x 10.05 78.95 y 0 1 z -14.8 -8 ;1.46 8.26 sa*]q~ a
3P'.)=}
;------------------------------------------------------------------------------------------------- CWt,cwFW
def gen_Sando5 bx Eb2D
command t~X wF(";
gen zone brick p0 ttx1,0,-16.26 p1 add tt2,0,0 p2 add 0,yc,0 p3 add 0,0,czz14 & )SUT+x(DU
p4 add tt2,yc,0 p5 add 0,yc,czz14 p6 add tt2,0,czz14 p7 add tt2,yc,czz14 & ^vo^W:
size 50 1 2 ;ratio 1.1 1 1 0.|tKetHq
v@"xEf1n[
end_command z26zl[.
end Y.jg
}oV
gen_Sando5 Wc03Sv&FZ
group sand_o5 range x 10.05 78.95 y 0 1 z -16.26 -14.8 ;0 1.46 $eRxCX?b2
;------------------------------------------------------------------------------------------------- * F~"4g
;------------------------------------------------------------------------------------------------- 3vmLftZE}
% E~4 Ur
def gen_Sando6 K'n^,
t
command
0qZ{:}`3
gen zone brick p0 0,0,-50 p1 add ttx1,0,0 p2 add 0,yc,0 p3 add 0,0,xh & 7P:0XML}
p4 add ttx1,yc,0 p5 add 0,yc,xh p6 add ttx1,0,xh p7 add ttx1,yc,xh & b*r1Jn"h
size 10 1 20 ;ratio 1.2 1 1 pk(<],0]X
Q Q1+u Y
end_command 56&s'
end _W tSZmW?
gen_Sando6 p;Bd zV>
group sand_o6 range x 0 10.05 y 0 1 z -50 -16.26 BI,K?D&W-
;------------------------------------------------------------------------------------------------- kB"S h_:m
def gen_Sando6 :3{@LOil^
command =@ V4V} ?
gen zone brick p0 ttx1,0,-50 p1 add tt2,0,0 p2 add 0,yc,0 p3 add 0,0,xh & B>#zrCD
p4 add tt2,yc,0 p5 add 0,yc,xh p6 add tt2,0,xh p7 add tt2,yc,xh & 8uS1HE \%
size 40 1 20 ;ratio 1.1 1 #j iQa"
VLu_SXlo*
end_command M)Tv (7
end Lg_y1Mu7o
gen_Sando6 @+B
.<@V
group sand_o6 range x 10.05 78.95 y 0 1 z -50 -16.26 ;-33.74 0 `H ^Nc\P#
;------------------------------------------------------------------------------------------------- r/:s2 oQ
;------------------------------------------------------------------------------------------------- J'^BxN&
; 二、生成左侧土体 {HvR24#
;------------------------------------------------------------------------------------------------- E9%xSMS8@
;-------------------------------------------------------------------------------------------------
QH9(l
def gen_Sand11 f x%'7/+
command ,N <;!6e
gen zone brick p0 0,0,ht1 p1 add 0,yc,0 p2 add ttx2,0,0 p3 add 0,0,ht2 & <w.V !"!
p4 add ttx2,yc,0 p5 add ttx2,0,ht2 p6 add 0,yc,ht2 p7 add ttx2,yc,ht2 & ZEUd?"gaR
size 1 10 3 ;ratio 1 1 1 1.2 t<}'/
)
VH+^G)^) W
end_command RAXqRP,iw
end mcS/-DaN?
gen_Sand11 u|BD%5+J
group sand_o7 range x -9.55 0 y 0 1 z -3.1 0 m!N_TOl-^
;------------------------------------------------------------------------------------------------- f1hi\p0q
def gen_Sand2 R b\=\
command tG{e(
gen zone brick p0 ttx2,0,ht1 p1 add 0,yc,0 p2 add tt1,0,0 p3 add 0,0,ht2 & w0^( jMQe^
p4 add tt1,yc,0 p5 add tt1,0,ht2 p6 add 0,yc,ht2 p7 add tt1,yc,ht2 & u%VO'}Gz
size 1 50 3 ;ratio 1 1.1 1 (3c,;koRR
burEo.=
end_command 1Qhx$If~
end 'rp(k\ pY
gen_Sand2 ^X"G~#v=q
group sand_o7 range x -78.45 -9.55 y 0 1 z -3.1 0 0ydAdgD
;------------------------------------------------------------------------------------------------- zu^?9k
def gen_Sand31 JS}W 4 N
command ZCbxL.fFz
gen zone brick p0 ttx2,0,-5 p1 add 0,yc,0 p2 add tt1,0,0 p3 add 0,0,czz11 & E' JVf%)
p4 add tt1,yc,0 p5 add tt1,0,czz11 p6 add 0,yc,czz11 p7 add tt1,yc,czz11 & 7+#^:;19`
size 1 50 2 ;ratio 1 1.1 1 IP&En8W+
<s]K~ Vo
end_command N;Hrc6nin^
end 4 h :Oo
gen_Sand31 H<X4R
group sand_o8 range x -78.45 -9.55 y 0 1 z -5 -3.1 ;11.26 13.16 hLYSYMUb
;------------------------------------------------------------------------------------------------- 5}e-\:J >B
def gen_Sand32 S;i^ucAF
command XrFyN(p
gen zone brick p0 ttx2,0,-8 p1 add 0,yc,0 p2 add tt1,0,0 p3 add 0,0,czz12 & V:Q fI
p4 add tt1,yc,0 p5 add tt1,0,czz12 p6 add 0,yc,czz12 p7 add tt1,yc,czz12 & [)Xu60? Q
size 1 50 3 ;ratio 1 1.1 1 8[(c'rl|)|
xm/v :hl=
end_command 2s?j5 Sd
end K]N ^6ome
gen_Sand32 hSp[BsF`,
group sand_o9 range x -78.45 -9.55 y 0 1 z -8 -5 ;8.26 11.26 Dn<2.!ZKQ
;------------------------------------------------------------------------------------------------- " *kWM
def gen_Sand33 NAx( Qi3
command
2Z7smDJ
gen zone brick p0 ttx2,0,-14.8 p1 add 0,yc,0 p2 add tt1,0,0 p3 add 0,0,czz13 & u?Iop/b
p4 add tt1,yc,0 p5 add tt1,0,czz13 p6 add 0,yc,czz13 p7 add tt1,yc,czz13 & 08yTTt76t
size 1 50 6 ;ratio 1 1.1 1 @hQ+pG@s
@UkcvhH
end_command _+z@Qn?#6h
end V <:kS
gen_Sand33 2=(=Wjk.
group sand_o10 range x -78.45 -9.55 y 0 1 z -14.8 -8 ;1.46 8.26 ehO F@IA_
}I#;~|v~<
;------------------------------------------------------------------------------------------------- H P*x?|4
def gen_Sand34 0*B_$E06
command [-s0'z
gen zone brick p0 ttx2,0,-16.26 p1 add 0,yc,0 p2 add tt1,0,0 p3 add 0,0,czz14 & e`<= & w
p4 add tt1,yc,0 p5 add tt1,0,czz14 p6 add 0,yc,czz14 p7 add tt1,yc,czz14 & >pG]#Z g
size 1 50 2 ;ratio 1 1.1 1 wf6ZzG:
:*1|ERGoay
end_command PrDvRWM
end Y\dK- M{$
gen_Sand34 F! c%& Z
group sand_o11 range x -78.45 -9.55 y 0 1 z -16.26 -14.8 ;0 1.46 i.0}qS?
;------------------------------------------------------------------------------------------------- kx]f`b
;------------------------------------------------------------------------------------------------- oopTo51,a
Fm*n>^P@Y
def gen_Sand4 XH1so1h
command xfos>|0N
gen zone brick p0 0,0,-50 p1 add 0,yc,0 p2 add ttx2,0,0 p3 add 0,0,xh & d3 fE[/oU
p4 add ttx2,yc,0 p5 add ttx2,0,xh p6 add 0,yc,xh p7 add ttx2,yc,xh & 67/hhO
size 1 10 20 75Jh(hd(
GB^Ch YOb
end_command wgkh} b
end !@ai=p
gen_Sand4 31Zl"-<#-
group sand_o12 range x -9.55 0 y 0 1 z -16.26 -50 '`/1?,=
;------------------------------------------------------------------------------------------------- /D 2v 1
def gen_Sand5 76zi)f1f
command .;/@k%>
gen zone brick p0 ttx2,0,-50 p1 add 0,yc,0 p2 add tt1,0,0 p3 add 0,0,xh & yY`<t
p4 add tt1,yc,0 p5 add tt1,0,xh p6 add 0,yc,xh p7 add tt1,yc,xh & hh
<=D .u
size 1 50 20 ;ratio 1 1.1 1 7! <cU
I&Yu=v/_
end_command 2ai \("?
end 1i2jYDB"
gen_Sand5 'ZT E"KT
group sand_o12 range x -78.45 -9.55 y 0 1 z -16.26 -50 ;-33.74 0 r\ ` R$
;------------------------------------------------------------------------------------------------- G)5% f\&
call Flacmodel_haitang.dat Ax=Rb
B"
;------------------------------------------------------------------------------------------------- ?Xh=rx_
plot add axes gPQ2i])"Q
;-------------------------------------------------------------------------------------- ~?[%uGI0h
;二、生成接触面的命令 tA}O'x
;-------------------------------------------------------------------------------------- WH/r$.&
interface 1 face range x -9.55 10.05 y 0 1 z -16.27 -16.25 .CQ
IN] iD
interface 2 face range x -9.55 10.05 y 0 1 z -3.11 -3.09 C Ij3D"
interface 3 face range x -9.56 -9.54 y 0 1 z -5 -3.1 v(h
interface 4 face range x -9.56 -9.54 y 0 1 z -8 -5 f o4j^,`
interface 5 face range x -9.56 -9.54 y 0 1 z -14.8 -8 2[qO;js
interface 6 face range x -9.56 -9.54 y 0 1 z -16.26 -14.8 n CGLuZn
BU<A+Pe>
int 1 prop kn 1.570e9 ks 1.570e9 fric 10 coh 10000 ;prop属性 kn法向刚度 ks剪切刚度 fric为摩擦角 coh内聚力 uD Q
d48>
int 2 prop kn 1.570e9 ks 1.570e9 fric 15 coh 55000 5s;H F |2x
int 3 prop kn 1.570e9 ks 1.570e9 fric 26 coh 8000 ;w}ZI<ou
int 4 prop kn 1.570e9 ks 1.570e9 fric 15 coh 55000 B=u@u([.
int 5 prop kn 1.570e9 ks 1.570e9 fric 14 coh 24000 |DwI%%0(F
int 6 prop kn 1.570e9 ks 1.570e9 fric 17 coh 56000 !OPa
`kSh
ini y add -2.0 rang group 1 PZeVjL?E
;=============================================== o-GlBXI;
;初始应力计算程序 hgfCM
;=============================================== 5~aSkg,MD
;定义模型及性质 `|
L+a~~
model mohr rang group sand_o1 % ]r@vjeyd
prop bulk 1.875e7 shear 0.967e7 coh 10000 fric 10 rang group sand_o1
/DQoM@X
model mohr rang group sand_o2 z"=#<C
prop bulk 1.875e7 shear 0.967e7 coh 55000 fric 15 rang group sand_o2 "Z& {
model mohr rang group sand_o3 hi`\3B
prop bulk 2.068e7 shear 0.793e7 coh 8000 fric 26 rang group sand_o3 -P(q<T2MV'
model mohr rang group sand_o4 T%
Kj >-
prop bulk 1.875e7 shear 0.967e7 coh 55000 fric 15 rang group sand_o4 ! Hdg
$,
model mohr rang group sand_o5 HGh`O\f8
prop bulk 1.396e7 shear 0.758e7 coh 24000 fric 14 rang group sand_o5 j:O=9
model mohr rang group sand_o6 Z +(V'e;
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o6 O292 JA
model mohr rang group sand_o7 q 1Rk'k4+
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o7 $*9h\W-)`Q
model mohr rang group sand_o8 lRZt)) 3
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o8 P7 H -Dw
model mohr rang group sand_o9 .FXq4who
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o9 )+u|qT3%
model mohr rang group sand_o10 ]Tx8ImD#)A
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o10 HZkC3$
model mohr rang group sand_o11 =5[}&W
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o11 )l\BZndf
model mohr rang group sand_o12 l`[*b_
Xt
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o12 ~{=+dQ
5m=3{lBi
mo elas rang group 1 UyEyk$6SU
;prop bulk 20.8e9 shear 15.58e9 rang group 1 ;正确 n;xtUw6 \
prop bulk 22.8e9 shear 17.08e9 rang group 1 ;只把柱子等效 TStu)6%`
ini dens 2500 rang group 1 }f; Zx)!
ini dens 1900 rang group sand_o1 ;18000 iNT 1lk
ini dens 1990 rang group sand_o2 r*p<7
ini dens 1910 rang group sand_o3 zX{ [Z
ini dens 1990 rang group sand_o4 .B6$U>>NS^
ini dens 1970 rang group sand_o5 }ytc oIuLf
ini dens 2010 rang group sand_o6 YaFQy0t%/5
ini dens 2010 rang group sand_o7 ,CA,7Mu:
ini dens 2010 rang group sand_o8 %K\_gR}V
ini dens 2010 rang group sand_o9 PHyS^J`
ini dens 2010 rang group sand_o10 X%h1r`h&
ini dens 2010 rang group sand_o11 &b i Bm
ini dens 2010 rang group sand_o12 r[Qk-}@vp
;边界条件 9V0iV5?( P
fix x y z rang z -50.1 -49.9 U*(m'Ea
fix x rang x 78.85 79.05 gk >A
fix x rang x -78.55 -78.35 kC!7 <%(
fix y rang y -.1 .1 /IxMRi=
fix y rang y 0.9 1.1 AVZ -g/<
;应力条件 15)=>=1mR.
ini szz 0 grad 0 0 20e3 rang z -50 0 ;加权平均重度为1993 CD +,&id
ini sxx 0 grad 0 0 10e3 rang z -50 0 V2u^sy
ini syy 0 grad 0 0 10e3 rang z -50 0 s~6?p%
2]
;--------------------------------------------------------------------------------- eg3zp gZ
hist id 1 unbal QDzFl1\P
plot hist 1 Y 'Yo c
; =SRp
set grav 0 0 -10 S "!nM]2L
solve ratio 1.0e-4 l=Jbu c
save station1.sav
