flac3d中應(yīng)力釋放的實(shí)現(xiàn)，在flac3d中，很多時(shí)候在隧道模擬計(jì)算中，在隧道縱向取1m長(zhǎng)計(jì)算，實(shí)際上是種 平面應(yīng)變計(jì)算，因此需要考慮巖體開(kāi)挖后的應(yīng)力釋放問(wèn)題，經(jīng)本人多次試算，總結(jié)大家提出的各種方法，認(rèn)為，下面方法可以實(shí)現(xiàn)應(yīng)力釋放的控制：
range name chuzhi group 5
range name weiyan group 3
;說(shuō)明 其中g(shù)roup 5 group 3 為需要釋放應(yīng)力的節(jié)點(diǎn)所共有的2個(gè)group
;開(kāi)挖
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
;應(yīng)力釋放
step 1
def relax. N4 k) p9 O# X, l s" u
relax_coef=0.0 ;應(yīng)力釋放系數(shù)，取值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 ;選擇需要加反向節(jié)點(diǎn)力的節(jié)點(diǎn)
xf=-gp_xfunbal(p_gp)*relax_m
yf=-gp_yfunbal(p_gp)*relax_m
zf=-gp_zfunbal(p_gp)*relax_m ;計(jì)算所需加的反向節(jié)點(diǎn)力
pid=gp_id(p_gp)
command
apply xforce xf rang id pid
apply yforce yf rang id pid
apply zforce zf rang id pid ;加反向節(jié)點(diǎn)力
endcommand
endif
endif
p_gp=gp_next(p_gp
endloop
end
relax
大家只需要根據(jù)自己的模型調(diào)整一下group的名稱，正確找到需要應(yīng)力釋放的節(jié)點(diǎn)，就可以很容易實(shí)現(xiàn)應(yīng)力釋放，并且可以如ansys里面一樣做到控制應(yīng)力釋放率，并且不影響計(jì)算速度。比如需要釋放30％，只需要設(shè)置relax_coef=0.3 就可以了，非常容易操作。