计算隧道最佳支护时间，一般是等效于在围岩释放不同比例的荷载后施作支护，如下面这个网页所做。
http://www.doc88.com/p-0418017102072.html
step 1
def relax
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

所以问题是“在加反向节点力后，直接计算再施加支护，再计算？还是加方向节点力后，直接加支护计算？”

这个用不了 从inrange开始就错了 我的是5.0