flac3d中应力释放的实现（总结加实现）

chchyao

flac3d中应力释放的实现，在flac3d中，很多时候在隧道模拟计算中，在隧道纵向取1m长计算，实际上是种平面应变计算，因此需要考虑岩体开挖后的应力释放问题，经本人多次试算，总结大家提出的各种方法，认为，下面方法可以实现应力释放的控制：
range  name chuzhi  group 5
range  name weiyan  group 3
;说明 其中group 5  group 3 为需要释放应力的节点所共有的2个group

;开挖
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
;应力释放
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

大家只需要根据自己的模型调整一下group的名称，正确找到需要应力释放的节点，就可以很容易实现应力释放，并且可以如ansys里面一样做到控制应力释放率，并且不影响计算速度。比如需要释放30％，只需要设置relax_coef=0.3 就可以了，非常容易操作。

如果有问题欢迎加我qq讨论  158536860

lessse

借鉴ansys果然不错，兄台貌似很猛啊

一米阳光！

兄台不是貌是很猛，是本来就很猛。牛人一个。

benjackxu

原帖由 chchyao 于 2008-3-24 13:55 发表
flac3d中应力释放的实现，在flac3d中，很多时候在隧道模拟计算中，在隧道纵向取1m长计算，实际上是种平面应变计算，因此需要考虑岩体开挖后的应力释放问题，经本人多次试算，总结大家提出的各种方法，认为，下面方法 ...

希望论坛多出现姚兄这样有心、有创新的软件使用者。

shilw

来赞美一个！*^_^*

wangyi2002598

赞一个

xbwei33

姚兄厉害，感谢

hbw8456

呵呵，发帖子根写论文一样，同样的内容要我说应该写出一个精华帖来，以让更多的人参与讨论，广发言论~~~
至少给个例子，贴几张图，再和其他几种方法做个比较~~~
这个工作是有点烦琐，想发个好论文不也很辛苦嘛~~~~

baomuqiqi83song

很有参考价值，谢谢

looloozheng

牛人啊!赞一个

中南土木人

呵呵！多谢楼主了啊！

freeflywyx

收藏
楼主强

kang8107

姚兄厉害，感谢

syx2006

牛人啊   我要加油向师兄看齐 嘿嘿:lol

looloozheng

我是根据分布开挖过程中值实际监测值，来发推应力释放的结果，根据支护时间来采用setp 多少步来解决应力释放！不知道这样可行否？？仅个人之言

lhjhh1973

高，实在是高！

qiandao

姚兄，类似ANSYS，在施作二衬时，应该把加的这个反力撤除（remove）吧？不过撤除的时候我老出错，还麻烦指点一下。

chen023536

高手啊！！！

070542

不错，顶

longyutian

不错，顶一个