FLAC3D隧道耦合實例

1、作一個盾構隧道開挖，地下水位設為z=12.6m,（隧道中心線設為z=0m）。分了 7層土，給定初始條件的代碼如下：4 s m-j! m3 p$ Q8 qconfig fluid. A/ N（ J8 c 15 m;& z P（過程省）.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* V& d6 p

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

3、ens 1e3: Y/ b, u2 Y1 Y% i cini fmod 8.5e7/ R9 b7 al 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 A. t4 H, Q* f; K. D% z$ Jfix pp rang

4、e 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 afix pp range z 12.5 12.7$ M- V/ & N# s: vfix pp range z -30.5 -30.3;material mechanic parameters& P* a% w. % 3 q(: s y% Q8 c & % A! model mohr# g0 V% u6 W, b. k% f& 0 X, w0 n* def derive+ J: 6 d- I F9 us_mod1=E

5、_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* pl Gs_mod4=E_mod4/(2.0*(1.0+

6、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+p_ratio7

7、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 bend( S/ y9 ! Q: Cset E_mod1=2.40e6 p_ratio1=0.25 E_mod2=5.0e6 p_ratio2=0.30E_mod3=11.0e6 p_ratio3=0.25 & 一一E_mod4=8.5e6 p_ratio4=0.25 E_mod5=11.5e6 p_ratio5=0.27E_mod6=9.25e6 p_ratio6=0.22 &, x-z i* y3 N$ B8 T% aE_m

8、od7=12.40e6 p_ratio7=0.20 m u, s0 d, B* m# d3 U; C derive9 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 z9.1 12.6 一一B D B 7 A0 d- # s% d8 Qprop bulk b_mod1 shear s_mod1 cohe

9、13.2e3 fric 10.0 ten 74.861e3 rangez 0 9.1 F% |6 b7 ?, l 一prop bulk b_mod1 shear s_mod1 cohe 13.2e3 fric 10.0 ten 74.861e3 rangez -0.9 0 一一prop bulk b_mod2 shear s_mod2 cohe 15.3e3 fric 10.0 ten 86.771e3 rangez -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 cohe 2

10、2.0e3 fric 20.0 ten 60.445e3 rangez -11.4 -7.4&Q o. d& W! nq7 7 Z4 x5 prop bulk b mod6 shear s mod6 cohe 3.0e3 fric 25.0 ten 6.434e3 range z -18.4 -11.4 一一prop 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# Z9

11、% 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 Jfix 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 _inter

12、face 1 prop kn 3e9 ks 1e9 fric 20 coh 3e5interface 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; 12 e% v( X r |,! l-( p o8 x! 8 C2 D5 3 P/ xini d

13、ens 1.76e3 range z 0 9.1 $ B9 g3 M1 A ! Eini 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 c0 1 T R& c9 k J& Y* Kini dens 2.0e3 range z -11.4 -7.42 u% M |; n2 K0 8 7 Yini dens 1.89e3 range z -18.4 -11.4ini dens 1.93e3 range z -30.4 -18.4 U& j. o7 |! r/ c% ! C-t

14、8 Y1 y$ 17 S; q$ Q. L0 八3 cini szz -28.237e4 grad 0 0 1.87e4 range z 12.6 15.1ini szz -28.237e4 grad 0 0 1.87e4 range z 9.1 12.6ini 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.9ini szz -25.

15、98e4 grad 0 0 2.0e4 range z -11.4 -7.4ini szz -27.234e4 grad0 01.89e4range z -18.4 -11.44 o5Z! U-J- H$mini szz -26.498e4 grad0 01.93e4range z -30.4 -18.4! R! |8 h/ 4 Z/ u0 Qini sxx -28.237e4 grad0 01.87e4range z 12.6 15.1ini sxx -28.237e4 grad0 01.87e4range z 9.1 12.6! p% c.D d3p! G7C0z9 Lini sxx -2

16、7.236e4 grad0 01.76e4range z 0 9.1ini sxx -27.236e4 grad0 01.76e4range z -0.9 0, i q8 l c.h$ n.e; Nini sxx -25.98e4 grad 0 0 2.0e4 range z -11.4 -7.4ini sxx -27.234e4 grad 0 0 1.89e4 range z -18.4 -11.4ini 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

17、 9.871e3 range z 12.6 15.17 V. K 2 h7 / s9 Wini syy -13.597e4 grad 0 0 8.833e3 range z 9.1 12.6ini syy -16.834e4 grad 0 0 1.239e4 range z 0 9.1ini syy -16.834e4 grad 0 0 1.239e4 range z -0.9 0ini syy -16.783e4 grad 0 0 1.296e4 range z -7.4 -0.94 h. q! E- $ $ D K6 qini syy -19.117e4 grad 0 0 9.806e3

18、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.47 Z- _& : Y, M（取控制點省）4 Q/ a h K |9 j2 S! tsolvesave iniconditions.sav初始平衡后的PP如下圖。請問：z=0以上都出現(xiàn)了負孔隙水壓，還有 z=12.6m以上無水，為什么出現(xiàn)孔隙水壓了呢？給定土層初始地.時，szz和syy中的水壓力部分是FLAC自動加上的 嗎? ; p2 M% i K3 Hi: k側壓力系數如何求？ # Z* G7 F+ r4 ?: J* m( x- nfmod應該取多少？我開始取e9級，計算.長的嚇人，頂部-還出現(xiàn)正 值。請看看。不勝感激!( Wum 國 Pine -I 皿0M TXCMKfAI 即2賓：”！:二3+、工里知心匕：皿m_ IMOKEKt 5_ SC