FLAC3D隧道耦合實(shí)例.doc

1、作一個(gè)盾構(gòu)隧道開(kāi)挖，地下水位設(shè)為 z=12.6m，（隧道中心線設(shè)為z=0m）。分了 7 層土，給定初始條件的代碼如下： 4 s m- j! m3 p$ Q8 qconfig fluid. A/ N( J8 c l5 m; & z P（建模過(guò)程省） . M7 P8 p6 S3 q8 c4 Rmodel fl_isotropicprop perm 1e-12 poro 0.48 range z 9.1 12.63 7 $ p- h% |* u8 B- r+ D: U: |$ |4 O5 Wprop perm 1e-12 poro 0.55 range z -0.9 9.1|3 e: S& l% C

2、* V& d6 p4 tprop perm 1e-12 poro 0.5 range z -7.4 -0.9- f x9 ? G J c$ l* _5 O/ : Bprop perm 1e-12 poro 0.5 range z -11.4 -7.4& u7 S * |6 n; S; c c v | U 4 X: / m& j2 S, K- v1 |4 xprop perm 1e-12 poro 0.45 range z -18.4 -11.4 K. E c: O, H( R4 pprop perm 1e-12 poro 0.45 range z -30.4 -18.4 2 7 j ?% R:

3、 # Cini fdens 1e3: Y/ b, u2 Y1 Y% i cini fmod 8.5e7/ R9 b7 a1 d* w H Y a8 P. a4 U/ m+ s+ Wini sat 0 range z 12.6 15.1ini sat 1 range z -30.4 12.62 Q e, q: |! n Lini pp 0 grad 0 0 -1e4 range z -30.4 12.6 ? Q7 k ?6 M7 Q/ h h3 s2 tfix pp range x -.1 .1fix pp range x 39.9 40.1 . t4 H, Q* f; K. D% z$ Jfi

4、x pp range y -.1 .1! b& R$ d* J ! C X- Efix pp range y 119.9 120.14 R+ w6 m8 B8 O- S, m+ B7 C8 m7 a fix pp range z 12.5 12.7$ M- V/ & N# s: v fix pp range z -30.5 -30.3;material mechanic parameters& P* % w. % 3 q( : s y % Q8 c & % ! model mohr# g0 V% u6 W, b. k% f& 0 X, w0 n* def derive+ J: 6 d- I F

5、9 us_mod1=E_mod1/(2.0*(1.0+p_ratio1b_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio11 V9 X* k D. b5 y; I% c. fs_mod2=E_mod2/(2.0*(1.0+p_ratio2b_mod2=E_mod2/(3.0*(1.0-2.0*p_ratio2! & ; B V( 4 ?& as_mod3=E_mod3/(2.0*(1.0+p_ratio3( M/ o2 c# s4 j, G7 N qb_mod3=E_mod3/(3.0*(1.0-2.0*p_ratio3 H M3 f* p1 Gs_mod4=E_mod4/

6、(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_ratio4s_mod5=E_mod5/(2.0*(1.0+p_ratio5# q4 x c% n4 W4 K# E5 v# b_mod5=E_mod5/(3.0*(1.0-2.0*p_ratio58 L2 M; i u$ ss_mod6=E_mod6/(2.0*(1.0+p_ratio6b_mod6=E_mod6/(3.0*(1.0-2.0*p_ratio6- W9 l i3 q0 n0 a; f2 _s_mod7=E_mod7/(2.0*(1.0+

7、p_ratio7 5 T/ R7 |, e bb_mod7=E_mod7/(3.0*(1.0-2.0*p_ratio7# M$ e$ f9 + 4 F: F2 w# b; - 6 y8 b end( S/ y9 ! Q: Cset E_mod1=2.40e6 p_ratio1=0.25 E_mod2=5.0e6 p_ratio2=0.30 E_mod3=11.0e6 p_ratio3=0.25 &E_mod4=8.5e6 p_ratio4=0.25 E_mod5=11.5e6 p_ratio5=0.27 E_mod6=9.25e6 p_ratio6=0.22 &, X z i* y3 N$ B

8、8 I% aE_mod7=12.40e6 p_ratio7=0.20 m u, s0 d, B* m# d3 U; Cderive9 i. N; k, Q# A; G* vprop 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.6B D B 7 A0 d- # s% d8 Qprop bulk b_mod1 shear s_mod1

9、 cohe 13.2e3 fric 10.0 ten 74.861e3 range z 0 9.1 F% |6 b7 ?, lprop 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 8 B+ c7 w0 e8 g. 4 T F S+ Gprop bulk b_mod4 shear s_mod4 co

10、he 22.0e3 fric 20.0 ten 60.445e3 range z -11.4 -7.4& Q o. 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* I0 M+ n/ D5 u8 i C4 3 L O! p;boundary conditions! y Q#

11、 Z9 % C7 Q$ f# y( 8 G$ Y7 s8 v# Ofix x range x -0.1 0.1 h( 6 + k% V0 bfix x range x 39.9 40.1fix y range y -0.1 0.10 N, A, u5 A( V& A/ ?. c( v J fix y range y 119.9 120.1fix x range z -30.5 -30.33 y ! ( v3 D6 |fix y range z -30.5 -30.3 # P8 F4 m8 X $ Nfix z range z -30.5 -30.3 m! L3 g% 7 v$ Y$ V3 _i

12、nterface 1 prop kn 3e9 ks 1e9 fric 20 coh 3e5 interface 2 prop kn 3e9 ks 1e9 fric 20 coh 3e5;stress conditions; w( a3 b- Z( P: x. 1 O 7 , Kset grav 0 0 -10 W6 j1 h* - $ y, E7 s% U Z$ Kini dens 1.87e3 range z 12.6 15.1ini dens 1.87e3 range z 9.1 12.6( K2 F; l2 e% v( X r |, ! l- ( p o8 x! 8 C2 D5 3 P/

13、 xini dens 1.76e3 range z 0 9.1 $ B9 g3 M1 A ! E ini dens 1.76e3 range z -0.9 0ini dens 1.84e3 range z -7.4 -0.9! e& g9 O; U- g7 h; B8 j4 c 0 1 T R& c9 k J& Y* Kini dens 2.0e3 range z -11.4 -7.42 u% M |; n2 K0 8 7 Y ini dens 1.89e3 range z -18.4 -11.4ini dens 1.93e3 range z -30.4 -18.4 U& j. o7 |! r

14、/ c% ! C - t8 Y1 y$ l7 S; q$ Q. L0 3 cini 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 0 X+ M9 h/ m( Zini szz -27.164e4 grad 0 0 1.84e4 range z -7.4 -0.

15、9 ini szz -25.98e4 grad 0 0 2.0e4 range z -11.4 -7.4ini szz -27.234e4 grad 0 0 1.89e4 range z -18.4 -11.44 o5 Z! U- J- H$ m ini szz -26.498e4 grad 0 0 1.93e4 range z -30.4 -18.4! R! |8 h/ 4 Z/ u0 Qini sxx -28.237e4 grad 0 0 1.87e4 range z 12.6 15.1ini sxx -28.237e4 grad 0 0 1.87e4 range z 9.1 12.6!

16、p% c. D d3 p! G7 C0 z9 Lini sxx -27.236e4 grad 0 0 1.76e4 range z 0 9.1ini sxx -27.236e4 grad 0 0 1.76e4 range z -0.9 0, i q8 l c. h$ n. e; N ini sxx -27.164e4 grad 0 0 1.84e4 range z -7.4 -0.9! Q; r: M: J iini sxx -25.98e4 grad 0 0 2.0e4 range z -11.4 -7.4ini sxx -27.234e4 grad 0 0 1.89e4 range z -

17、18.4 -11.4 ini sxx -26.498e4 grad 0 0 1.93e4 range z -30.4 -18.4- I8 ?4 i1 b0 8 t& dini syy -14.905e4 grad 0 0 9.871e3 range z 12.6 15.17 V. K 2 h7 / s9 W 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.1ini syy -16.834e4 grad 0 0 1.239e4 range z -0

18、.9 0ini syy -16.783e4 grad 0 0 1.296e4 range z -7.4 -0.94 h. q! E- $ $ D K6qini syy -19.117e4 grad 0 0 9.806e3 range z -11.4 -7.4 ; U- R U5 J; B! n v2 Dini syy -21.550e4 grad 0 0 7.672e3 range z -18.4 -11.43 _, z6 _- H- 5 d9 f3 ini syy -21.252e4 grad 0 0 7.834e3 range z -30.4 -18.47Z-_&:Y,M（取控制點(diǎn)省） 4 Q/ a h K |9 j2 S! tsolvesave iniconditions.sav初始平衡后的 PP如下圖。請(qǐng)問(wèn)：z=0 以上都出現(xiàn)了負(fù)孔隙水壓，還有z=12.6m 以上無(wú)水，為什么出現(xiàn)孔隙水壓了呢？給定土層初始地 應(yīng)力 時(shí)， szz和 syy 中的水壓力部分是 FLAC 自動(dòng)加上的嗎？ ; p2 M% i K3 H i: k側(cè)壓力系數(shù)如何求？ # Z* G7 F+ r4 ?: J* m( x-