FLAC隧道应力释放.doc

《FLAC隧道应力释放.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、flac3d中应力释放的实现，在flac3d中，很多时候在隧道模拟计算中，在隧道纵向取1m长计算，实际上是种平面应变计算，因此需要考虑岩体开挖后的应力释放问题，经本人多次试算，总结大家提出的各种方法，认为，下面方法可以实现应力释放的控制：:t+]:^;h-C;Mrange  namechuzhi  group5.w9I$S#V#R5E-Frange  nameweiyan  group3;说明其中group5  group3为需要释放应力的节点所共有的2个group*l;3m2o.

2、;p#I  O/M0_(N7j2Y!f9i+;开挖

3、,P8t+a5N  i0e1v/Z;iinistate07f1t"p(X8F;n%yinixdis=0ydis=0zdis=00q3T0Z-E!W!A5i7f8vsetlargemodelnullrangegroup5modelnullrangegroup6modelnullrangegroup7;M'^%]-z,U6y+C;应力释放*i2@6N.w2t  @5kstep1defrelax1q.`(Z7r*C  relax_coef=0.0  ;应力释放系数，取值0－1之间  relax_m=1.0-relax_coefp_gp=gp_

4、head+j+n8W)^$n7E.a.r  loopwhilep_gp#null    ifinrange('chuzhi',p_gp)=1then ifinrange('weiyan',p_gp)=1then  ;选择需要加反向节点力的节点8K*x%{9?8D#Q0_      xf=-gp_xfunbal(p_gp)*relax_m1P&!B.q:s&k&u      yf=-gp_yfunbal(p_gp)*relax_m      zf=-gp_zfunbal(p_gp)*relax_m;计算所需加的反向节点力5P%x+w  

5、F*g.v2~3k      pid=gp_id(p_gp))F$D6K&R4O1x1j      command7T  ?+a#x#g6N      applyxforcexfrangidpid+c0b2*b,P+y      applyyforceyfrangidpid;H.{0y5G0^8D$t9V7W  x/t:p      applyzforcezfrangidpid  ;加反向节点力2}+Q&d*s'D'S,_;z      endcommand  Endif,Q,J"f&j%l  endif,h)y8w;D2u7z+X)

6、[    p_gp=gp_next(p_gp)  Endloop)L!J5b&z5@7Qend+i8b9@(b5P+Frelax-h:I+s$Y9w&V/O${)P%A&  B2q,Z&^#r;M大家只需要根据自己的模型调整一下group的名称，正确找到需要应力释放的节点，就可以很容易实现应力释放，并且可以如ansys里面一样做到控制应力释放率，并且不影响计算速度。比如需要释放30％，只需要设置relax_coef=0.3就可以了，非常容易操作。$@7l1j(Q1@9`