FLAC3D隧道耦合實例

1、z=Om)。分了 7層土，給定初始條件的代碼如下:config fluid A/ N( J8 c l5 m; & z P(建模過程省). M7 P8 p6 S3 q8 c4model fl_isotropicprop perm 1e-12 poro 0.48 range z 9.1 12.37$6p-h%|prop perm 1e-12 poro 0.55 range z -0.9 9.1prop perm 1e-12 poro 0.5 range z -7.4 -0.9prop perm 1e-12 poro 0.5 range z -11.4 -7.&u7S4*|6n;S;ccv|Upro

2、p perm 1e-12 poro 0.45 range z -18.4 -11.4 prop perm 1e-12 poro 0.45 range z -30.4 -18.4 ini fdens 1 e3:Y/b,u2 Y1 Y% i cini fmod 8.5e7/R9 b7 a1 d* w H) Y) aini sat 0 range z 12.6 15.1ini sat 1 range z -30.4 12.62Qe,q:|!nLini pp 0 grad 0 0 -1e4 range z -30.4 12.6?Q7k?6M7Q/h hfix pp range x -.1 .1fix

3、pp range x 39.9 40.1fix pp range y -.1 .1!b&R$ d* J ! C) X- Efix pp range y 119.9 1 20 . 14R+w6m8 B8 O- S, m+ B7 C8 m7 afix pp range z 12.5 12.7$M-V/&N# s: vfix pp range z -30.5 -30.3;material mecha nic parametersmodel mohr#g0V% u6 W, b. % f& def derive+J:6 d- I) F9 us_mod1=E_mod1/(2.0*(1.0+p_ratio1

4、)b_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio1)s_mod2=E_mod2/(2.0*(1.0+p_ratio2)b_mod2=E_mod2/(3.0*(1.0-2.0*p_ratio2)s_mod3二E_mod3/(2.0*(1.0+p_ratio3)b_mod3二E_mod3/(3.0*(1.0-2.0*p_ratio3)s_mod4=E_mod4/(2.0*(1.0+p_ratio4)/0 3 y/ H c: d; b3 h( qb_mod4=E_mod4/(3.0*(1.0-2.0*p_ratio4)s_mod5=E_mod5/(2.0*(1.0+p_rat

5、io5)b_mod5=E_mod5/(3.0*(1.0-2.0*p_ratio5)s_mod6=E_mod6/(2.0*(1.0+p_ratio6)b_mod6=E_mod6/(3.0*(1.0-2.0*p_ratio6)s_mod7=E_mod7/(2.0*(1.0+p_ratio7)5 T/ R7 |, e) bb_mod7=E_mod7/(3.0*(1.0-2.0*p_ratio7)end(S/yE_mod2=5.0e6 p_ratio2=0.30set E_mod1=2.40e6 p_ratio1=0.25E_mod3=11.0e6 p_ratio3=0.25 &E_mod4=8.5e

6、6 p_ratio4=0.25E_mod5=11.5e6 p_ratio5=0.27E_mod6=9.25e6 p_ratio6=0.22 &,X)zi* yE_mod7=12.40e6 p_ratio7=0.2deriveprop bulk b_mod5 shear s_mod5 cohe 18.0e3 fric 18.0 ten 55.398e3 range z 12.6 15.1prop bulk b_mod3 shear s_mod3 cohe 6.3e3 fric 21.0 ten 16.412e3 range z 9.1 12.6prop bulk b_mod1 shear s_m

7、od1 cohe 13.2e3 fric 10.0 ten 74.861e3 range z 0 9.1)F%|6b7?,prop bulk b_mod1 shear s_mod1 cohe 13.2e3 fric 10.0 ten 74.861e3 range z -0.9 0prop bulk b_mod2 shear s_mod2 cohe 15.3e3 fric 10.0 ten 86.771e3 range z -7.4 -0.9. Q* I5 S$ H2 x prop bulk b_mod4 shear s_mod4 cohe 22.0e3 fric 20.0 ten 60.445

8、e3 range z -11.4 -7.4&Qo.d&W!n;q7 7 Z4 x5 prop bulk b_mod6 shear s_mod6 cohe 3.0e3 fric 25.0 ten 6.434e3 range z -18.4 -11.4prop bulk b_mod7 shear s_mod7 cohe 3.0e3 fric 25.0 ten 6.434e3 range z -30.4 -18.4;boundary conditionsQ# Z9 % C7 Q $ ffix x range x -0.1 0.1fix x range x 39.9 40.1fix y range y

9、 -0.1 0.10N, A, u5 A( V& A/ ?. c( v Jfix y range y 119.9 120.1fix x range z -30.5 -30.3fix y range z -30.5 -30.3 # P8 F4 m8 X $ Nfix z range z -30.5 -30.3m!L3g% 7 v$ Y$ V3 _interface 1 prop kn 3e9 ks 1e9 fric 20 coh 3e5interface 2 prop kn 3e9 ks 1e9 fric 20 coh 3e5;stress conditions;w(a3b-Z(P:x. 1 O

10、 7 , Kset grav 0 0 -10)W6j1h*- ini dens 1.87e3 range z 12.6 15.1ini dens 1.87e3 range z 9.1 12.(K62F;l2e%v(X)r)|ini dens 1.76e3 range z 0 9.1$B9g3M1A ! Eini dens 1.76e3 range z -0.9 0ini dens 1.84e3 range z -7.4 -0.! 9e&g9O;U-g7h;B8j4cini dens 2.0e3 range z -11.4 -7.24u%M|;n2K08 7 Yini dens 1.89e3 r

11、ange z -18.4 -11.4ini de ns 1.93e3 range z -30.4 -18.4ini szz -28.237e4 grad 0 0 1.87e4 range z 12.6 15.1 ini szz -28.237e4 grad 0 0 1.87e4 range z 9.1 12.6 ini szz -27.236e4 grad 0 0 1.76e4 range z 0 9.1ini szz -27.236e4 grad 0 0 1.76e4 range z -0.9 X+ M9 h/ m0(Zini szz -27.164e4 grad 0 0 1.84e4 ra

12、nge z -7.4 -0.9 ini szz -25.98e4 grad 0 0 2.0e4 range z -11.4 -7.4 ini szz -27.234e4 grad 0 0 1.89e4 range z -18.4 -114o5Z.!U-J-H$m4 ini szz -26.498e4 grad 0 0 1.93e4 range z -30.4 -18.4ini sxx -28.237e4 grad 0 0 1.87e4 range z 12.6 15.1 ini sxx -28.237e4 grad 0 0 1.87e4 range z 9.1 12.!p%c.D 6d3p!G

13、7C0 z9 L ini sxx -27.236e4 grad 0 0 1.76e4 range z 0 9.1 ini sxx -27.236e4 grad 0 0 1.76e4 range z -0.9 , i q8 l c.0h$n.e;N ini sxx -27.164e4 grad 0 0 1.84e4 range z -7.4 -0.!Q;r:M:9Ji ini sxx -25.98e4 grad 0 0 2.0e4 range z -11.4 -7.4 ini sxx -27.234e4 grad 0 0 1.89e4 range z -18.4 -11.4 ini sxx -2

14、6.498e4 grad 0 0 1.93e4 range z -30.4 -18.4ini syy -14.905e4 grad 0 0 9.871e3 range z 12.6 157V.K2h7/1s9W ini syy -13.597e4 grad 0 0 8.833e3 range z 9.1 12.6 ini syy -16.834e4 grad 0 0 1.239e4 range z 0 9.1 ini syy -16.834e4 grad 0 0 1.239e4 range z -0.9 0 ini syy -16.783e4 grad 0 0 1.296e4 range z

15、-7.4 -0.4h.q!E9-$DK6q ini syy -19.117e4 grad 0 0 9.806e3 range z -11.4 -7.4;U-R)U5J;B!nv2D ini syy -21.550e4 grad 0 0 7.672e3 range z -18.4 -11.3_,z6_-H- 45d9f3 ini syy -21.252e4 grad 0 0 7.834e3 range z -30.4 -18.4(取控制點省) 4 Q/ ah K) |9 j2 S! tsolvesave iniconditions.sav初始平衡后的PP如下圖。請問：z=0以上都出現(xiàn)了負孔隙水壓

16、，還有 z=12.6m以上無水，為什么出 現(xiàn)孔隙水壓了呢？給定土層初始地應力時，szz和syy中的水壓力部分是FLAC自動力口 上的嗎？側壓力系數(shù)如何求？fmod應該取多少？我開始取 e9級，結果計算時間長的嚇人，頂部 位移還出現(xiàn)正值。請高手幫忙看看。不勝感激！FLAC3D 3.00Steo B711 Model Perspective1041:45 Sat Dec 09 23Q6Center:X： 3.D28e+-D0lY： 4,E23e+001Z: -2J54e+0013 st: 1.856e*002Rotation:X: 2D.0Q0Y： 0.0002: 30 000Mag.:1Ang.; 22.500Contour of Pore ProcureMaifac = 0 000e+000J 2600e+005tO-1.0000e+005 -rOOOOe+OOSto-S.dOOQe+DM-5.0D00e+D04to 0.0000e+000C.DOOOe+OOOto 5.00008+(X)4 5.Q0Q0e+0(Mtc 1.0000e+Q051.00006+005 to 1.50006 tO 051.5000e+005b： 2.i)00Gef005Interval = SOe+WSketdiM a Tfac = 0.0006+0302.aOOOe+OC5t. 2.5000005