flac3d中应力释放实现方法与总结

flac3d中应力释放的实现，在flac3d中，很多时候在隧道模拟计算中，在隧道纵向取1m长计算，实际上是种 平面应变计算，因此需要考虑岩体开挖后的应力释放问题，经本人多次试算，总结大家提出的各种方法，认为，下面方法可以实现应力释放的控制：
range name chuzhi group 5
range name weiyan group 3
;说明 其中group 5 group 3 为需要释放应力的节点所共有的2个group
;开挖
ini state 0
ini xdis=0 ydis=0 zdis=0
set large
model null range group 5
model null range group 6
model null range group 7
;应力释放
step 1
def relax. N4 k) p9 O# X, l s" u
relax_coef=0.0 ;应力释放系数，取值0－1之间
relax_m=1.0-relax_coef
p_gp=gp_head
loop while p_gp # null
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
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
endif
endif
p_gp=gp_next(p_gp
endloop
end
relax
大家只需要根据自己的模型调整一下group的名称，正确找到需要应力释放的节点，就可以很容易实现应力释放，并且可以如ansys里面一样做到控制应力释放率，并且不影响计算速度。比如需要释放30％，只需要设置relax_coef=0.3 就可以了，非常容易操作。