隧洞开挖

千里飞雪

求教各位大神，在进行隧洞开挖的静力计算时，模拟应力释放，得到的结果完全不同，而且又不提示错误，下面是命令流，求教各位大神！谢谢！


;位移归零
ini xdis 0 ydis 0 zdis 0
ini xvel 0 yvel 0 zvel 0
ini state=0
set large

model null range group 5
model null range group 6
model null range group 7
model null range group 8
model null range group 9
step 1
def relax
relax_coef=0.3
relax_m=1.0-relax_coef
p_gp=gp_head
loop while p_gp#null
if inrange ('chuchen',p_gp)=1 then
if inrange ('fangshui',p_gp)=1 then
if inrange ('erchen',p_gp)=1 then
if inrange ('weiyan',p_gp)=1 then
xf=-gp_xfunbal(p_gp)*relax_coef
yf=-gp_yfunbal(p_gp)*relax_coef
zf=-gp_zfunbal(p_gp)*relax_coef
pid=gp_id(p_gp)
command print pid
apply xforce xf range id pid
apply yforce yf range id pid
apply zforce zf range id pid
end command
end if
end if
p_gp=gp_next(p_gp)
endloop
end relax
save 开挖后结果保存.sav

千里飞雪

自己顶一下！