《FLAC隧道应力释放》由会员分享,可在线阅读,更多相关《FLAC隧道应力释放(2页珍藏版)》请在金锄头文库上搜索。
1、精品文档,仅供学习与交流,如有侵权请联系网站删除flac3d中应力释放的实现,在flac3d中,很多时候在隧道模拟计算中,在隧道纵向取1m长计算,实际上是种平面应变计算,因此需要考虑岩体开挖后的应力释放问题,经本人多次试算,总结大家提出的各种方法,认为,下面方法可以实现应力释放的控制:: t+ : ; h- C; Mrangename chuzhigroup 5. w9 I$ S# V# R5 E- Frangename weiyangroup 3% 6 d4 w; J8;说明 其中group 5group 3 为需要释放应力的节点所共有的2个group* l; 3 m2 o. |; p# I
2、O/ M0 _( N7 j2 Y! f9 i+ ;开挖, P8 t+ a5 Ni0 e1 v/ Z; iini state 07 f1 t p( X8 F; n% yini xdis=0 ydis=0 zdis=00 q3 T0 Z- E! W! A5 i7 f8 vset large; y1 x-model null range group 5model null range group 6model null range group 7; M % - z, U6 y+ C;应力释放* i2 6 N. w2 t5 kstep 1) o. def relax1 q. ( Z7 r* Crela
3、x_coef=0.0 ;应力释放系数,取值01之间+relax_m=1.0-relax_coef* k8p_gp=gp_head+ j+ n8 W) $ n7 E. a. r loop while p_gp # null8 H5 if inrange(chuzhi,p_gp) = 1 then if inrange(weiyan,p_gp) = 1 then;选择需要加反向节点力的节点8 K* x% 9 ?8 D# Q0 _ xf=-gp_xfunbal(p_gp)*relax_m1 P& ! B. q: s& k& u yf=-gp_yfunbal(p_gp)*relax_m; 2 f5 k
4、* p9 zf=-gp_zfunbal(p_gp)*relax_m ;计算所需加的反向节点力5 P% x+ wF* g. v2 3 k pid=gp_id(p_gp) F$ D6 K& R4 O1 x1 j command7 T?+ a# x# g6 N apply xforce xf rang id pid+ c0 b2 * b, P+ y apply yforce yf rang id pid; H. 0 y5 G0 8 D$ t9 V7 Wx/ t: p apply zforce zf rang id pid ;加反向节点力2 + Q& d* s D S, _; z endcommand
5、# f, B! s9 1 Endif, Q, J f& j% l endif, h) y8 w; D2 u7 z+ X) p_gp=gp_next(p_gp)2 I Endloop) L! J5 b& z5 7 Qend + i8 b9 ( b5 P+ Frelax- h: I+ s$ Y9 w& V/ O$ ) P% A& B2 q, Z& # r; M大家只需要根据自己的模型调整一下group的名称,正确找到需要应力释放的节点,就可以很容易实现应力释放,并且可以如ansys里面一样做到控制应力释放率,并且不影响计算速度。比如需要释放30,只需要设置relax_coef=0.3 就可以了,非常容易操作。$ 7 l1 j( Q1 9 * 【精品文档】第 页