《最新FLAC3D隧道耦合实例.doc》由会员分享，可在线阅读，更多相关《最新FLAC3D隧道耦合实例.doc（20页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、Four short words sum up what has lifted most successful individuals above the crowd: a little bit more.-author-dateFLAC3D隧道耦合实例FLAC3D隧道耦合实例作一个盾构隧道开挖，地下水位设为z=12.6m，（隧道中心线设为z=0m）。分了7层土，给定初始条件的代码如下：4 sm- j! m3 p$ Q8 q$ V. J: ?. Y5 3 z, o* c* gconfig fluid. A/ N( J8 c l5 m; & z P（建模过程省）. M7 P8 p6 S3 q8
2、c4 Rmodel fl_isotropic# x2 X% M3 n |- ?% % w6 C1 j) Q0 J7 r$ p; J2 u. n8 iprop 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 O5 Q 5 : v0 H) |3 e: S& l% C* V& d6 p4 tprop perm 1e-12 poro 0.5 range z -7.4 -0.9- f x9 ?) GJ c$ l*
3、_5 O/ : B( D1 L; L+ M- q* : F* Fprop perm 1e-12 poro 0.5 range z -11.4 -7.4& u7 S* |6 n; S; c cv |U4 X: / m& j2 S, K- v1 |4 xprop perm 1e-12 poro 0.45 range z -18.4 -11.4: _ f* P: K8 c) O ) K. E c: O, H( R4 pprop perm 1e-12 poro 0.45 range z -30.4 -18.47 S; # VS; k- 6 k; c& M1 _7 pm2 7 j ?% R: # Cin
4、i fdens 1e3: Y/ b, u2 Y1 Y% i c! S; 5 $ A- P$ fini fmod 8.5e7/ R9 b7 a1 d* w H) Y) a8 P. a4 U/ m+ s+ Wini sat 0 range z 12.6 15.1) _( C, I/ I9 i7 y R5 Zini sat 1 range z -30.4 12.62 Q e, q: |! n L7 GQ6 K* Q( H pini pp 0 grad 0 0 -1e4 range z -30.4 12.6 ? Q7 k?6 M7 Q/ h h3 s2 t& hx( D6 v/ Q! Y3 T( D3
5、 U9 k: ?3 r; ; F3 ufix pp range x -.1 .1) 9 Z3 T0 ) H/ w. u% 0 G) fix pp range x 39.9 40.1 . t4 H, Q* f; K. D% z$ Jfix 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 afix pp range z 12.5 12.7$ M- V/ & N# s: vfix pp range z -30.5 -30.3 . D; R; _2 E6 o1
6、 h1 M7 8 o n! L8 V7 t& ;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) F9 us_mod1=E_mod1/(2.0*(1.0+p_ratio1)7 2 U+ R$ v7 W* q2 Qb_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio1)1 V9 X* k D. b5 y; I% c. fs_mod2=E_mod2/(2.0*(1.0+p_
7、ratio2)3 A) bR* s X7 w6 Z7 Vb_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/(2.0*(1.0+p_ratio4)/ 0 3 y/ Hc: d; b3 h( qb_mod4=E_mod4/(3.0*(1.0-2.0*p_ratio4) T/ a: h0 r) O* k
8、& u$ s_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_ratio5)8 L2 M; i u$ ss_mod6=E_mod6/(2.0*(1.0+p_ratio6)% F; bm8 a/ ?# em$ b_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) 5 T/ R7 |, e) bb_mod7=E_mod7/(3.0*(1.0-2
9、.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.30 E_mod3=11.0e6 p_ratio3=0.25 &# l. X- q/ M3 j5 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) zi* y3 N$ B8 I% a E_mod7=12.40e6 p_ratio7=
10、0.2: L. q1 . r* |& 0 m u, s0 d, B* m# d3 U; Cderive % V8 b+ uT5 l, y- P8 t9 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.1; m* g) j/ P8 4 I5 gprop bulk b_mod3 shear s_mod3 cohe 6.3e3 fric 21.0 ten 16.412e3 range z 9.1 12.60 l u P& q8 6 B) D) B 7
11、 A0 d- # s% d8 Qprop bulk b_mod1 shear s_mod1 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 0 q. t. x1 L7 n2 Oh! vprop bulk b_mod2 shear s_mod2 cohe 15.3e3 fric 10.0 ten 86.771e3 range z -7.4 -0.9 . Q* I5
12、S$ H2 x8 B+ c7 w0 e8 g. 4 TF S+ Gprop bulk b_mod4 shear s_mod4 cohe 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.4$ B 6 ST, r; v/ n( qprop bulk b_mod7 shear s_mod7 cohe 3.0e3 fric 25.0 ten
13、6.434e3 range z -30.4 -18.47 w% d, X e iC3 m; m* | * I0 M+ n/ D5 u8 i C4 3 L) O! p;boundary conditions! y Q# 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.17 m; v. D+ J$ Pfix y range y -0.1 0.10 N, A, u5 A( V& A/ ?. c( v Jfix y range y 119.9 120.11 U5 y%
14、S J+ 1 Z; u* h7 sfix 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 _- c: q$ s2 F: H+ fY( CXinterface 1 prop kn 3e9 ks 1e9 fric 20 coh 3e5f4 Y z, z3 Z4 yinterface 2 prop kn 3e9 ks 1e9 fric 20 coh 3e5, ?7 r& J, T! u; W S$ Bd
15、! M4 C: X5 V;stress conditions; w( a3 b- Z( P: x. 1 O 7 , K9 o# u. l G1 E/ bset grav 0 0 -10) W6 j1 h* - $ y, E7 s% U Z$ K! S( . F2 dl8 Zini dens 1.87e3 range z 12.6 15.1+ P7 A0 w( C6 u+ x. ini 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/ xini dens 1.76e3 range z
16、 0 9.1 $ B9 g3 M1 A! Eini dens 1.76e3 range z -0.9 03 4 R i) ; G8 V* N3 Z sini 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.4# C+ v- i) V! e oini dens 1.93e3 range z -30.4 -18.4)
17、U& j. o7 |! r/ ) 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% y! e0 V/ x/ w9 N9 s J, ?/ eini szz -28.237e4 grad 0 0 1.87e4 range z 9.1 12.67 b% k$ 9 B. n1 M2 i B8 Qini szz -27.236e4 grad 0 0 1.76e4 range z 0 9.16 X; w, M# I* # T. V( b! s: Eini szz -27.236e4
18、 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.9% e3 6 a2 b2 t# I2 ini szz -25.98e4 grad 0 0 2.0e4 range z -11.4 -7.45 t, a. V- 2 O$ e/ e0 E+ d) D! ini szz -27.234e4 grad 0 0 1.89e4 range z -18.4 -11.44 o5 Z! U- J- H$ mini szz -26.498e4 grad 0 0 1.93e4
19、 range z -30.4 -18.41 N+ E2 I- O5 C9 K8 ?- C! R! |8 h/ 4 Z/ u0 Qini sxx -28.237e4 grad 0 0 1.87e4 range z 12.6 15.16 h( c* q1 3 y. 8 H! Tini sxx -28.237e4 grad 0 0 1.87e4 range z 9.1 12.6! p% c. Dd3 p! G7 C0 z9 Lini sxx -27.236e4 grad 0 0 1.76e4 range z 0 9.1, u5 H6 o3 H+ ?ini sxx -27.236e4 grad 0 0
20、 1.76e4 range z -0.9 0, i q8 l c. h$ n. e; Nini sxx -27.164e4 grad 0 0 1.84e4 range z -7.4 -0.9! Q; r: M: Jiini sxx -25.98e4 grad 0 0 2.0e4 range z -11.4 -7.4. G2 g x* k- C) h# 8 Q! Dini sxx -27.234e4 grad 0 0 1.89e4 range z -18.4 -11.41 9 4 f+ + s9 kini sxx -26.498e4 grad 0 0 1.93e4 range z -30.4 -
21、18.45 z X/ Z6 Q! u ( - e- 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 Wini syy -13.597e4 grad 0 0 8.833e3 range z 9.1 12.65 ?# ? K m# w ) Lini syy -16.834e4 grad 0 0 1.239e4 range z 0 9.15 a9 i% & s( b1 Z& v z% z* t( Cini syy -16.834e4 grad 0 0 1.239e4 ra
22、nge z -0.9 0& b- ?! M9 x. ini 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 range z -11.4 -7.4 ; U- R) U5 J; B! nv2 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
23、-18.4- q; ) _4 X1 k1 T5 S7 Z- _& : Y, M（取控制点省）4 Q/ a h K) |9 j2 S! tsolve! _1 _ d9 j9 r3 T/ N! psave iniconditions.sav) b) d, X4 ! y4 X4 v& ln* V, _% K1 V! |& , 0 |初始平衡后的PP如下图。请问：8 ; S3 , q/ z=0以上都出现了负孔隙水压，还有z=12.6m以上无水，为什么出现孔隙水压了呢？( R0 Y% ( w/ O# C9 t0 p给定土层初始地应力时，szz和syy中的水压力部分是FLAC自动加上的吗？; p2 M%
24、i K3 H) i: k侧压力系数如何求？# Z* G7 F+ r4 ?: J* m( x- nfmod应该取多少？我开始取e9级，结果计算时间长的吓人，顶部位移还出现正值。+ e: % + X5 k1 Y1 C9 d( g请高手帮忙看看。不胜感激！你施加的孔压边界没问题，初始孔压加错了3 o$ R# ! J/ X8 # p应为ini pp 1.26e5 grad 0 0 -1e4 range z 12.6 -30.4, ( F- 5 P7 1 K! r- B地层中的孔压分布是难以确定的! Y1 zi / Z: O既然是人为施加的孔压，那边在初始平衡时就不要让它变了U& J& & M O7 j( S, j% V用INI FMOD 0% K9 O! L2 ! x8 c1 A) x即可“ini pp 0 grad 0 0 -1e4 range z -30.4 12.6”该命令流有问题，在Z=12.6时，其孔隙水压力并不为0，而是等于-12.6e4-