解除节点应力的方法

fc3535

在进行隧道开挖时宜考虑围岩的应力释放，西南交大姚长春博士提出可以通过在衬砌与围岩边界上的施加节点力达到这种效果。隧道衬砌拼装模拟之后应将施加的节点力消除。经本人实验，解除施加的节点力的方法有以下两种，且效果差不多。
1）将应力释放率设为1（姚的方法），命令流如下：
def relax
        relax_m=1.0-relax_coef
        p_gp=gp_head
        loop while p_gp#null
                if inrange('kaiwa',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 range id pid
                                        apply yforce yf range id pid
                                        apply zforce zf range id pid
                                endcommand
                        endif
                endif
                p_gp=gp_next(p_gp)
        endloop
end
set relax_coef=1.0

2）使用remove命令
def remove_nf
        p_gp=gp_head
        loop while p_gp#null
                if inrange('kaiwa',p_gp)=1 then
                        if inrange('weiyan',p_gp)=1 then                               
                                pid=gp_id(p_gp)
                                command
                                        apply remove gp range id pid                                       
                                endcommand
                        endif
                endif
                p_gp=gp_next(p_gp)
        endloop
end

fc3535

自己顶一个。交流方能成长，希望大家提出宝贵的想法！

enochwang

解除倒还容易，如果是多次分步应力释放，该怎么实现呢？尝试过多次，也不知对错