{keZ_ 2
new LD$5KaOW
config dyn ~P4C`Q1PT#
set dyn off K`4rUEf}V"
;=============================================== ]P4Wf V
d
;建模型 &OMlW _FHR
;=============================================== 0Y.z
def setup ;定义函数setup `zd,^.i5~
yc=1 ;y方向尺寸 54
lD+%E
czz1=-13.16 p~pD`'%
czz11=1.9 ;四分之一车站z方向尺寸 @KN+)q P
czz12=3 ;四分之一车站z方向尺寸 Z oW1Cc&p
czz13=6.8 ;四分之一车站z方向尺寸 H?$dnwR
czz14=1.46 ;四分之一车站z方向尺寸 wj)LOA0
ttx1=10.05 ;四分之一车站外围土体x方向尺寸 7#V7D6j1
ttx2=-9.55 ;四分之一车站外围土体-x方向尺寸 ;Rrh$Ag
ht1=-3.1 ke4E 1T-1n
ht2=3.1 !FEc:qH
tt2=68.9 ;78.95 s\*L5{kiSl
tt1=-68.9 =~D? K9o
xh=33.74 ；车站以下土体高度 i2F7O" f.
BZq#OA p
end oV|O`n
setup mG~_*8}e<
;------------------------------------------------------------------------------------------------- _jWs(OmJ
; 一、生成右侧土体 XlmX3RU
;------------------------------------------------------------------------------------------------- B}e/MlX3M
def gen_Sando1 <Xf6?nyZ(
command q6zVu (
gen zone brick p0 0,0,ht1 p1 add ttx1,0,0 p2 add 0,yc,0 p3 add 0,0,ht2 & W5*ldX Xk
p4 add ttx1,yc,0 p5 add 0,yc,ht2 p6 add ttx1,0,ht2 p7 add ttx1,yc,ht2 & TO wd+]B
size 10 1 3 ;ratio 1 1 1 1.2 UhxM85M;x
> I2rj2M#
end_command -IPo/?}
end ]\E"oZ
gen_Sando1 5N U{y+
group sand_o1 range x 0 10.05 y 0 1 z -3.1 0 ER0
Yl
;------------------------------------------------------------------------------------------------- vygzL U^
def gen_Sando1 5mH [|_
command PmR].Ohzi
gen zone brick p0 ttx1,0,ht1 p1 add tt2,0,0 p2 add 0,yc,0 p3 add 0,0,ht2 & r`H}f#.KR
p4 add tt2,yc,0 p5 add 0,yc,ht2 p6 add tt2,0,ht2 p7 add tt2,yc,ht2 & agIqca;
size 50 1 3 ;ratio 1.1 1 1 : 4ryi&Y
_&gi4)q
end_command Y(Q
0m|3P
end S{c/3k ~
gen_Sando1 /6x&%G:m#
group sand_o1 range x 10.05 78.95 y 0 1 z -3.1 0 ?
SFBUX(p
;------------------------------------------------------------------------------------------------- 144Y.
def gen_Sando2 DfXkLOGik
command 5`;SI36"
gen zone brick p0 ttx1,0,-5 p1 add tt2,0,0 p2 add 0,yc,0 p3 add 0,0,czz11 & /{T&l* '
p4 add tt2,yc,0 p5 add 0,yc,czz11 p6 add tt2,0,czz11 p7 add tt2,yc,czz11 & ; pnF%co9
size 50 1 2 ;ratio 1.1 1 1 bI):-2&s}
TI>5g(:3\
end_command _v!7
|&\
end X+X:nL.t
gen_Sando2 ZmA }i`
group sand_o2 range x 10.05 78.95 y 0 1 z -5 -3.1 ;11.26 13.16 rrD6x>
]so/AdT9hA
Xp<q`w 0I,
;------------------------------------------------------------------------------------------------- *x36;6~W;
;------------------------------------------------------------------------------------------------- ah>c)1DA*H
def gen_Sando3 ~7m` p3W@
command L/k35 x8
gen zone brick p0 ttx1,0,-8 p1 add tt2,0,0 p2 add 0,yc,0 p3 add 0,0,czz12 & fSo8O
p4 add tt2,yc,0 p5 add 0,yc,czz12 p6 add tt2,0,czz12 p7 add tt2,yc,czz12 & 4:/]Y=)x
size 50 1 3 ;ratio 1.1 1 1 MZgaQU g
Kb~nC6yJc
end_command H.5
6
end ,=/9Ld2w9
gen_Sando3 (tK_(gO
group sand_o3 range x 10.05 78.95 y 0 1 z -8 -5 ;8.26 11.26 a k(s@@k
;------------------------------------------------------------------------------------------------- w$]G$e
def gen_Sando4 P1#g{f
command r!r08y f
gen zone brick p0 ttx1,0,-14.8 p1 add tt2,0,0 p2 add 0,yc,0 p3 add 0,0,czz13 & l*Cul VX
p4 add tt2,yc,0 p5 add 0,yc,czz13 p6 add tt2,0,czz13 p7 add tt2,yc,czz13 & |J $A%2 7
size 50 1 6 ;ratio 1.1 1 1 { FJM c O=
b
74 !Zw
end_command 7HkO:/
end _u]Z+H"
gen_Sando4 UZdnsG7
group sand_o4 range x 10.05 78.95 y 0 1 z -14.8 -8 ;1.46 8.26 K\xz|Gq
lrrTeE*
;------------------------------------------------------------------------------------------------- ~OD6K`s3
def gen_Sando5
f<o|5r
command I|X` 9
gen zone brick p0 ttx1,0,-16.26 p1 add tt2,0,0 p2 add 0,yc,0 p3 add 0,0,czz14 & ,!=
sGUQ)
p4 add tt2,yc,0 p5 add 0,yc,czz14 p6 add tt2,0,czz14 p7 add tt2,yc,czz14 & K,4Ig!
size 50 1 2 ;ratio 1.1 1 1 fB2ILRc
ZCCCuB
end_command 2&hv6Y1
end R
iZ)FW
gen_Sando5 E|=x+M1sH
group sand_o5 range x 10.05 78.95 y 0 1 z -16.26 -14.8 ;0 1.46 ;qmnG3;Q
;------------------------------------------------------------------------------------------------- #2=l\y-#
;------------------------------------------------------------------------------------------------- ~WrpJjI[
Et6 j6gmif
def gen_Sando6 .Xlo-gHk
command :zW? O#aL-
gen zone brick p0 0,0,-50 p1 add ttx1,0,0 p2 add 0,yc,0 p3 add 0,0,xh & /xn|d#4
p4 add ttx1,yc,0 p5 add 0,yc,xh p6 add ttx1,0,xh p7 add ttx1,yc,xh & <bcf"0A
size 10 1 20 ;ratio 1.2 1 1 lMv6QL\>'
2/SUEnaLy_
end_command 2d ts}G
end mnTF40l
gen_Sando6 "&_+!TBg,
group sand_o6 range x 0 10.05 y 0 1 z -50 -16.26 }T_"Vg q
;------------------------------------------------------------------------------------------------- !t % 1G.
def gen_Sando6 6 (rm%c
command 8vx
ca]DcV
gen zone brick p0 ttx1,0,-50 p1 add tt2,0,0 p2 add 0,yc,0 p3 add 0,0,xh &
?K-4T
p4 add tt2,yc,0 p5 add 0,yc,xh p6 add tt2,0,xh p7 add tt2,yc,xh & r CJ$Pl9R
size 40 1 20 ;ratio 1.1 1 Zma;An6
Z~:)hwF
end_command xI,3( A.
end CvIuH=,
gen_Sando6 FM]clC;X?
group sand_o6 range x 10.05 78.95 y 0 1 z -50 -16.26 ;-33.74 0 d=1\= d/K
;------------------------------------------------------------------------------------------------- O>E}Lu;|
;------------------------------------------------------------------------------------------------- Wvb ~ j
; 二、生成左侧土体 RHsVG &<j
;------------------------------------------------------------------------------------------------- 2j
f!o
;------------------------------------------------------------------------------------------------- 4s{=/,f
def gen_Sand11 A &}]:4@{
command D;sG9H ky
gen zone brick p0 0,0,ht1 p1 add 0,yc,0 p2 add ttx2,0,0 p3 add 0,0,ht2 & k$=L&id
p4 add ttx2,yc,0 p5 add ttx2,0,ht2 p6 add 0,yc,ht2 p7 add ttx2,yc,ht2 & J;sQvPHV8
size 1 10 3 ;ratio 1 1 1 1.2 ^"Bhp:o2
=Ti!9 _~
end_command A 6z2KVk
end W@!qp
gen_Sand11 =8_b&4.:&
group sand_o7 range x -9.55 0 y 0 1 z -3.1 0 w,hm_aDq
;------------------------------------------------------------------------------------------------- ^*jw e^
def gen_Sand2 $H*8H`
command VjM/'V5
gen zone brick p0 ttx2,0,ht1 p1 add 0,yc,0 p2 add tt1,0,0 p3 add 0,0,ht2 & ]ty$/{hx'
p4 add tt1,yc,0 p5 add tt1,0,ht2 p6 add 0,yc,ht2 p7 add tt1,yc,ht2 & j32*9
size 1 50 3 ;ratio 1 1.1 1 nscnG5'{+
=2q#- ,t
end_command {?Slo5X|
end m{oe|UVcmr
gen_Sand2 (~Z &U
group sand_o7 range x -78.45 -9.55 y 0 1 z -3.1 0 [kJ;Uxncz~
;------------------------------------------------------------------------------------------------- p$"~v A .
def gen_Sand31 !S~)U{SSK
command "%K'~"S#Q,
gen zone brick p0 ttx2,0,-5 p1 add 0,yc,0 p2 add tt1,0,0 p3 add 0,0,czz11 & r?KRK?I
p4 add tt1,yc,0 p5 add tt1,0,czz11 p6 add 0,yc,czz11 p7 add tt1,yc,czz11 & Q``1^E'
size 1 50 2 ;ratio 1 1.1 1 E2|M#Y
"s@Hg1
end_command 5$
How!
end x~k3kj
gen_Sand31 x,: k /]
group sand_o8 range x -78.45 -9.55 y 0 1 z -5 -3.1 ;11.26 13.16 IF W(nB(
;------------------------------------------------------------------------------------------------- r'ilJ("
def gen_Sand32 f0eQq;D$K
command PU% Zay
gen zone brick p0 ttx2,0,-8 p1 add 0,yc,0 p2 add tt1,0,0 p3 add 0,0,czz12 & yR5XcPo KI
p4 add tt1,yc,0 p5 add tt1,0,czz12 p6 add 0,yc,czz12 p7 add tt1,yc,czz12 & :TVo2Zm[@
size 1 50 3 ;ratio 1 1.1 1 +:Xg7H*
B<~ AUf*y
end_command 7( 5 d$ W
end ;nl J D#
gen_Sand32 p$F` 9_bZ
group sand_o9 range x -78.45 -9.55 y 0 1 z -8 -5 ;8.26 11.26 ~ LI }
;------------------------------------------------------------------------------------------------- }S"gZ6
def gen_Sand33 SN|!FW.*:
command xlcCL?qQj
gen zone brick p0 ttx2,0,-14.8 p1 add 0,yc,0 p2 add tt1,0,0 p3 add 0,0,czz13 & /#5 rt &q
p4 add tt1,yc,0 p5 add tt1,0,czz13 p6 add 0,yc,czz13 p7 add tt1,yc,czz13 & 75>%!mhM
size 1 50 6 ;ratio 1 1.1 1 e O~p"d-|
~ M7X]
end_command iwIn3R,
end 0aq{ Y7sYU
gen_Sand33 Dh m ;K$T
group sand_o10 range x -78.45 -9.55 y 0 1 z -14.8 -8 ;1.46 8.26 3N]ushMO
/@Jg [na
;------------------------------------------------------------------------------------------------- RJ'za1@z;b
def gen_Sand34 zvN7aG
command A46dtFD{
gen zone brick p0 ttx2,0,-16.26 p1 add 0,yc,0 p2 add tt1,0,0 p3 add 0,0,czz14 & [cFD\"gJAr
p4 add tt1,yc,0 p5 add tt1,0,czz14 p6 add 0,yc,czz14 p7 add tt1,yc,czz14 & 2cu2S"r
size 1 50 2 ;ratio 1 1.1 1 TlO=dLR7d
LQqba4$
end_command !C4)P3k
end AW;xlY= g
gen_Sand34 gWro])3
group sand_o11 range x -78.45 -9.55 y 0 1 z -16.26 -14.8 ;0 1.46 3Ch 4 2<
;------------------------------------------------------------------------------------------------- t .&JPTK-H
;------------------------------------------------------------------------------------------------- ZT"vVX- )G
V (XU^}b#
def gen_Sand4 S^q^=q0F
command m
Urb
gen zone brick p0 0,0,-50 p1 add 0,yc,0 p2 add ttx2,0,0 p3 add 0,0,xh & n=f?Q= h\3
p4 add ttx2,yc,0 p5 add ttx2,0,xh p6 add 0,yc,xh p7 add ttx2,yc,xh & ?gu!P:lZS
size 1 10 20 )a!f")@uz
@WXRZEz
end_command 1A 93ol=
end /g0' +DP
gen_Sand4 %X}D(_
group sand_o12 range x -9.55 0 y 0 1 z -16.26 -50 zp4@T)
;------------------------------------------------------------------------------------------------- 5T"h7^}e
def gen_Sand5 8 bw, dBN
command ao7M([ff
gen zone brick p0 ttx2,0,-50 p1 add 0,yc,0 p2 add tt1,0,0 p3 add 0,0,xh & vWW Q/^
p4 add tt1,yc,0 p5 add tt1,0,xh p6 add 0,yc,xh p7 add tt1,yc,xh &
w@mCQ$
size 1 50 20 ;ratio 1 1.1 1 WCaMPz
cEXd#TlY~X
end_command C(sz/x?11
end {
zL4dJw
gen_Sand5 .$f0!`
t
group sand_o12 range x -78.45 -9.55 y 0 1 z -16.26 -50 ;-33.74 0 J;V#a=I
;------------------------------------------------------------------------------------------------- +pUYFDwFx
call Flacmodel_haitang.dat =C#z Px,
;------------------------------------------------------------------------------------------------- * W"Pv,:
plot add axes *?Oh%.HgF
;-------------------------------------------------------------------------------------- l"&iSq!3=
;二、生成接触面的命令 <%4pvn8d?&
;-------------------------------------------------------------------------------------- amQiH!}8R
interface 1 face range x -9.55 10.05 y 0 1 z -16.27 -16.25 8z- wdO\
interface 2 face range x -9.55 10.05 y 0 1 z -3.11 -3.09 1)e[F#|
interface 3 face range x -9.56 -9.54 y 0 1 z -5 -3.1 u3cg&l EgT
interface 4 face range x -9.56 -9.54 y 0 1 z -8 -5 -R$ Q`Xw
interface 5 face range x -9.56 -9.54 y 0 1 z -14.8 -8 ?!t O'}?
interface 6 face range x -9.56 -9.54 y 0 1 z -16.26 -14.8 .K_50 %s
i*xVD`x ~
int 1 prop kn 1.570e9 ks 1.570e9 fric 10 coh 10000 ;prop属性 kn法向刚度 ks剪切刚度 fric为摩擦角 coh内聚力 !! 6@r|.
int 2 prop kn 1.570e9 ks 1.570e9 fric 15 coh 55000 t[({KbIy
int 3 prop kn 1.570e9 ks 1.570e9 fric 26 coh 8000 K-(C5 "j_
int 4 prop kn 1.570e9 ks 1.570e9 fric 15 coh 55000 =%77~q-HL
int 5 prop kn 1.570e9 ks 1.570e9 fric 14 coh 24000 4X:mb}(
int 6 prop kn 1.570e9 ks 1.570e9 fric 17 coh 56000 ,S}wOjb@
ini y add -2.0 rang group 1 u w>y*OLU+
;=============================================== C:!&g~{cKi
;初始应力计算程序 FIJ]`
;=============================================== (h&=N a~
;定义模型及性质 wO %617Av
model mohr rang group sand_o1 D.F1^9Q
prop bulk 1.875e7 shear 0.967e7 coh 10000 fric 10 rang group sand_o1 8M4GforP
model mohr rang group sand_o2
8`<3rj
prop bulk 1.875e7 shear 0.967e7 coh 55000 fric 15 rang group sand_o2 [8b{Yba z
model mohr rang group sand_o3 ?U:c\TA,m
prop bulk 2.068e7 shear 0.793e7 coh 8000 fric 26 rang group sand_o3 *EU1`q*
model mohr rang group sand_o4 1"l48NL L|
prop bulk 1.875e7 shear 0.967e7 coh 55000 fric 15 rang group sand_o4 XU['lr&,W
model mohr rang group sand_o5 HT_nxe`E
prop bulk 1.396e7 shear 0.758e7 coh 24000 fric 14 rang group sand_o5 7&dF=/:X@
model mohr rang group sand_o6 +n YF9z2
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o6 >)`*:_{
model mohr rang group sand_o7 }Zp [f6^Q
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o7 DI :
model mohr rang group sand_o8 M VE:JNm
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o8 gE23C*!'&:
model mohr rang group sand_o9 C;Kq_/l
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o9 Bgs U:eKe
model mohr rang group sand_o10
" v' %M({
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o10 SphP@J<ONW
model mohr rang group sand_o11 yTWicW7i
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o11
4f213h
model mohr rang group sand_o12 ?R @u'4yK
prop bulk 2.178e7 shear 1.183e7 coh 56000 fric 17 rang group sand_o12 |i (@1 l
!'bZ|j%
mo elas rang group 1 >b*}Td~J
;prop bulk 20.8e9 shear 15.58e9 rang group 1 ;正确 ) Cm95,Y
prop bulk 22.8e9 shear 17.08e9 rang group 1 ;只把柱子等效 C61E=$
ini dens 2500 rang group 1 |kHzp^S
ini dens 1900 rang group sand_o1 ;18000 )uCa]IR
ini dens 1990 rang group sand_o2 @B svk9}
ini dens 1910 rang group sand_o3 SG)|4$"
ini dens 1990 rang group sand_o4 aqP"Y9l
ini dens 1970 rang group sand_o5 5n=~l [O
ini dens 2010 rang group sand_o6 vf!lhV-UG+
ini dens 2010 rang group sand_o7 YQ-V^e6
ini dens 2010 rang group sand_o8 V4<f4|IL
ini dens 2010 rang group sand_o9 M r~IVmtf
ini dens 2010 rang group sand_o10 wem hP8!gc
ini dens 2010 rang group sand_o11 |x ir93 |
ini dens 2010 rang group sand_o12 |v" &Y
;边界条件 _10I0Z0
fix x y z rang z -50.1 -49.9 _ dVA^m
fix x rang x 78.85 79.05 T \TKgO=)
fix x rang x -78.55 -78.35 aslb^
fix y rang y -.1 .1 P.Qz>c^-C
fix y rang y 0.9 1.1 5BM rn0
;应力条件 *kNXju
ini szz 0 grad 0 0 20e3 rang z -50 0 ;加权平均重度为1993 RPTIDA))
ini sxx 0 grad 0 0 10e3 rang z -50 0 E`q)vk
ini syy 0 grad 0 0 10e3 rang z -50 0 xf3/J{ n3
;--------------------------------------------------------------------------------- 6 Znt
hist id 1 unbal 1[(/{CClB
plot hist 1 [58qC:
; 7!g4 `@!5M
set grav 0 0 -10 k|rbh.Q
solve ratio 1.0e-4 Yx"~_xA/u
save station1.sav
