FLAC3D隧道耦合实例.doc

《FLAC3D隧道耦合实例.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、作一个盾构隧道开挖，地下水位设为z=12.6m，（隧道中心线设为z=0m）。分了7层土，给定初始条件的代码如下：4s  m-j!m3p$Q8q$V.J:?.Y5@3z,o*c*gconfigfluid.A/N(J8c'l5m;~&z'P（建模过程省）.M7P8p6S3q8c4Rmodelfl_isotropic#x2X%M3n"

2、-?%^%w6C1j)Q0J7r$p;J2u.n8i  propperm1e-12poro0.48rangez9.112.63~7~$p-h%

3、*u8B-r+D:U:

4、$

5、4O5W  propperm1e-12poro0.55ran

6、gez-0.99.1"O5Q'@5]:v0H  )

7、3e:S&l%C*V&d6p4t  propperm1e-12poro0.5rangez-7.4-0.9-f"x9?)^"G  J"c$l*_5O/`:B  (D1L;L+M-q*`:F*F  propperm1e-12poro0.5rangez-11.4-7.4&u7S  ~*

8、6n;S;c"c  v"

9、  U4X:@/m&j2S,K-v1

10、4x  propperm1e-12poro0.45rangez-18.4-11.4:_'f*P:K8c)O)K.E"c:O,H(R4p  propperm1e-12p

11、oro0.45rangez-30.4-18.47S;~#V  S;k-6k;c&M1_7p  m2@7j"?%R:{#C  inifdens1e3:Y/b,u2Y1Y%i'c!S;]5@$A-P$f  inifmod8.5e7/R9b7a1d*w"H)Y)a8P.a4U/m+s+W  inisat0rangez12.615.1)_(C,I/I9i7y'R5Z  inisat1rangez-30.412.62Q"e,q:

12、!n'L7{  G  Q6K*Q(H"p  inipp0grad00-1e4rangez-30.412.6'?'Q7k  ?6M7Q/h'

13、h3s2t  &h  x(D6v/Q!Y3T(D3@  'U9k:?3r;{;F3u  fixpprangex-.1.1)]9Z3T0[)H/w.u%`0G)@  fixpprangex39.940.1"^.t4H,Q*f;K.D%z$J  fixpprangey-.1.1!b&R$d*J  @!C)X-E  fixpprangey119.9120.14R+w6m8B8O-S,m+{"B7C8m7a  fixpprangez12.512.7$M-V/~&N#s:v  fixpprangez-30.5-30.3.D;R;_2E6o  1h1M7@8o"n!L 

14、 ]8V7t&}  ;materialmechanicparameters&P*^%w.`%`3q([:s"y%Q8c'~&[%^!@  modelmohr#g0V%u6W,b.k%f&}0]'X,w0n*`  defderive+J:@6d-I)F9u  s_mod1=E_mod1/(2.0*(1.0+p_ratio1))7~2U+R$v7W*q2Q  b_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio1))1V9X*k"D.b5y;I%c.f  s_mod2=E_mod2/(2.0*(1.0+p_ratio2))3A)b  R*s"X

15、7w6Z7V  b_mod2=E_mod2/(3.0*(1.0-2.0*p_ratio2))!&[;B)V(`4?&a  s_mod3=E_mod3/(2.0*(1.0+p_ratio3))(M/o2c#^"s4j,G7N'q  b_mod3=E_mod3/(3.0*(1.0-2.0*p_ratio3))"H)M3f*p1G  s_mod4=E_mod4/(2.0*(1.0+p_ratio4))/}0]3y/H  c:d;b3h(q  b_mod4=E_mod4/(3.0*(1.0-2.0*p_ratio4))"T/a:h0r)O*k&u$~  s_mod

16、5=E_mod5/(2.0*(1.0+p_ratio5))#q4x'c%n4W4K#E5v#  b_mod5=E_mod5/(3.0*(1.0-2.0*p_ratio5))8L2^"M;i"u$s  s_mod6=E_mod6/(2.0*(1.0+p_ratio6))%F;b  m8a/?#]  e  m$]  b_mod6=E_mod6/(3.0*(1.0-2.0*p_ratio6))-W9l"i3q0n0a;f2_  s_mod7=E_mod7/(2.0*(1.0+p_ratio7))"{5T/R7

17、,e)b  b_mod7=E_mod7/(3.0*(

18、1.0-2.0*p_ratio7))#M$e$f9[