大家好,我在模拟盾构循环开挖,考虑了应力释放,命令流如下:
;施加释放荷载
command
step 1
endcommand
relax_m=1.0-relax_coef
;管片轮廓界面释放
p_gp=gp_head
loop while p_gp # null
if inrange('cz1',p_gp) = 1 then
if inrange('fnt',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
;删除释放荷载
p_gp=gp_head
loop while p_gp # null
if inrange('cz1',p_gp) = 1 then
pidd=gp_id(p_gp)
command
apply remove gp rang id pidd
apply remove gp rang id pidd
apply remove gp rang id pidd
endcommand
endif
p_gp=gp_next(p_gp)
endloop