FLAC隧道应力释放

FLAC 隧道应力释放(3 页)Good is good, but better carries it.精益求精，善益求善。flac3d 中应力释放的实现，在 flac3d中，很多时候在隧道模拟计算中，在隧道纵向取 1m 长计算，实际上是种平面应变计算，因此需要考虑岩体开挖后的应力释放问题，经本人多次试算，总结大家提出的各种方法，认为，下面方法可以实现应力释放的控制：range name chuzhi group 5range name weiyan group 3% {6 d4 w; J8;说明 其中 group 5 group 3 为需要释放应力的节点所共有的 2 个 group;开挖ini state 0ini xdis=0 ydis=0 zdis=0set large; y1 x-model null range group 5model null range group 6model null range group 7;应力释放step 1) @' o. def relax relax_coef=0.0 ;应力释放系数，取值 0－1之间+ relax_m=1.0-relax_coef* k8p_gp=gp_head loop while p_gp # null8 H5 if inrange('chuzhi',p_gp) = 1 then if inrange('weiyan',p_gp) = 1 then ;选择需要加反向节点力的节点 xf=-gp_xfunbal(p_gp)*relax_m yf=-gp_yfunbal(p_gp)*relax_m ~; {2 f5 k* p9 zf=-gp_zfunbal(p_gp)*relax_m ;计算所需加的反向节点力 pid=gp_id(p_gp) command apply xforce xf rang id pid apply yforce yf rang id pid apply zforce zf rang id pid ;加反向节点力 endcommand# f, B! s9 }1 Endif endif p_gp=gp_next(p_gp)2 I Endloopendrelax大家只需要根据自己的模型调整一下 group的名称，正确找到需要应力释放的节点，就可以很容易实现应力释放，并且可以如 ansys 里面一样做到控制应力释放率，并且不影响计算速度。比如需要释放 30％，只需要设置relax_coef=0.3 就可以了，非常容易操作。* 》~%》|:…）（}:《￥?$&~。/\%.{!.—）￥[【