FLAC隧道应力释放

flac3d中应力释放的实现，在flac3d中，很多时候在隧道模拟计算中，在隧道纵向取1m长计算，实际上是种平面应变计算，因此需要考虑岩体开挖后的应力释放问题，经本人多次试算，总结大家提出的各种方法，认为，下面方法可以实现应力释放的控制：rangenamechuzhigroup5rangenameweiyangroup3%{6d4w;J8;说明其中group5group3为需要释放应力的节点所共有的2个group;开挖inistate0inixdis=0ydis=0zdis=0setlarge;y1x-modelnullrangegroup5modelnullrangegroup6modelnullrangegroup7;应力释放step1)@'o.defrelaxrelax_coef=0.0;应力释放系数，取值0－1之间+relax_m=1.0-relax_coef*k8p_gp=gp_headloopwhilep_gp#null8H5ifinrange('chuzhi',p_gp)=1thenifinrange('weiyan',p_gp)=1then;选择需要加反向节点力的节点xf=-gp_xfunbal(p_gp)*relax_myf=-gp_yfunbal(p_gp)*relax_m~;{2f5k*p9zf=-gp_zfunbal(p_gp)*relax_m;计算所需加的反向节点力pid=gp_id(p_gp)commandapplyxforcexfrangidpidapplyyforceyfrangidpidapplyzforcezfrangidpid;加反向节点力endcommand#f,B!s9}1Endifendifp_gp=gp_next(p_gp)2IEndloopendrelax大家只需要根据自己的模型调整一下group的名称，正确找到需要应力释放的节点，就可以很容易实现应力释放，并且可以如ansys里面一样做到控制应力释放率，并且不影响计算速度。比如需要释放30％，只需要设置relax_coef=0.3就可以了，非常容易操作。*