flac3d中应力释放的实现

flac3d中应力释放的实现，在flac3d中，很多时候在隧道模拟计算中，在隧道纵向取1m长计算，实际上是种平面应变计算，因此需要考虑岩体开挖后的应力释放问题，经本人多次试算，总结大家提出的各种方法，认为，下面方法可以实现应力释放的控制：rangenamechuzhigroup5#@#~3I#x7^:Nj,bX&I+Nrangenameweiyangroup3;说明其中group5group3为需要释放应力的节点所共有的2个group;y,y7B2d:P0Z!U+s(Q*Q;开挖4}4Cg0`1v/J0e6K:`dinistate0inixdis=0ydis=0zdis=01x-[*?*g1|5GM8w%z_0Asetlarge-a8XQ#_%X8CC5Z,Wmodelnullrangegroup5modelnullrangegroup6modelnullrangegroup7*Q0H7?(E+A&D(]7b;应力释放step1:v/n*{5l5b2K8aQ)z7Ndefrelax/C'~,V5|$y4g8mrelax_coef=0.0;应力释放系数，取值0－1之间5F#\!Q*c%G#f0B$@relax_m=1.0-relax_coefp_gp=gp_head/\;Y3E/]:Nloopwhilep_gp#null2I5g3c#I?I(E0{%x:oifinrange('chuzhi',p_gp)=1thenifinrange('weiyan',p_gp)=1then;选择需要加反向节点力的节点xf=-gp_xfunbal(p_gp)*relax_myf=-gp_yfunbal(p_gp)*relax_mzf=-gp_zfunbal(p_gp)*relax_m;计算所需加的反向节点力$N5G#M9`)M5U8S)|pid=gp_id(p_gp)commandapplyxforcexfrangidpidapplyyforceyfrangidpidapplyzforcezfrangidpid;加反向节点力!n;A0D/CM0h0?,Y+`+B;e/dendcommandendifn(y%r8y,^)Bendif%\7_4A9P!J.n2|p_gp=gp_next(p_gp)0V8oX$c7a4@5jv&s3`7\endloop+A;U+Z9V0^end;k9~,^+I!S0X3y.Srelax-c#z)['g:m2Y大家只需要根据自己的模型调整一下group的名称，正确找到需要应力释放的节点，就可以很容易实现应力释放，并且可以如ansys里面一样做到控制应力释放率，并且不影响计算速度。比如需要释放30％，只需要设置relax_coef=0.3就可以了，非常容易操作。*z&l!]0J0{8c4q+_$D9u