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

whutzll · 发表于 2011-9-27 15:26:22

正需要这方面的参考资料呢，多谢了

hustzhdm · 发表于 2011-9-28 11:04:05

顶一个，谢谢楼主分享！！！

gaochao20wo · 发表于 2011-10-25 22:32:38

非常感谢，明年我毕业了，我也会把我的论文的命令流上传，向师兄们学习看齐

gaochao20wo · 发表于 2011-10-26 22:29:54

请各位大哥指点一下啊，我的沉降位移很小啊，

gaochao20wo · 发表于 2011-10-28 11:11:37

plot block group

range  name chuzhi  group 3
range  name weiyan  group 19

def excavate
cut=0
loop while cut <=10.0
cut1 =cut+2
cut2=cut-2

;开挖
command
model null range group 3 z cut cut1          ;限定开挖长度
model null range group 8 z cut cut1
model null range group 14 z cut cut1
set hist_rep 1
set large
step 1

his  unbal
his  gp ydisp 17.464 6.1606 2    ;拱顶节点
his  gp xdisp 24.686 1.0216 2    ;上台阶右拱角节点
his  gp xdisp 10.243 1.0216 2    ;上台阶左拱角节点
end_command

if cut =10.0
command
set plot jpg
set plot quality 100
plot cont syy
plot hard file 7.jpg
plot cont sxx
plot hard file 8.jpg
plot cont ydis
plot hard file 9.jpg
end_command
end_if

if cut >=2 then
command
model elas range group 3 z cut2 cut
prop bulk=1.167e8 shear=0.875e8 range group 3 z cut2 cut
ini density 2500 range group 3 z cut2 cut
set hist_rep 1
step 1

def relax
  relax_coef=0.4
  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 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
relax

his  unbal
his  gp ydisp 17.464 6.1606 2    ;拱顶节点
his  gp xdisp 24.686 1.0216 2    ;上台阶右拱角节点
his  gp xdisp 10.243 1.0216 2    ;上台阶左拱角节点
end_command
end_if
cut =cut+2
end_loop

command
set plot jpg
set plot quality 100
plot cont syy
plot hard file 10.jpg
plot cont sxx
plot hard file 11.jpg
plot cont ydis
plot hard file 12.jpg
end_command

end
excavate
plot block group

gaochao20wo · 发表于 2011-10-28 11:14:46

请教一下各位大哥，上面的程序中应力释放对么？我是在初支后面加应力释放，开挖后没加应力释放，这样对么？是不是要在开挖后面加应力释放，初支后再加应力释放？

gongzhenyu · 发表于 2011-10-28 16:08:20

很好的程序，收藏了

stuy_1001823 · 发表于 2012-2-22 19:36:31

挺好的东西

:D:D

caopangziabc · 发表于 2012-2-29 22:01:04

step 1 能实现吗，我改了几个应力释放系数，但zdis结果一样啊

小珍鹏 · 发表于 2012-3-27 09:24:08

借鉴下，很不错。

:lol

gao9941 · 发表于 2012-4-26 13:16:38

兄台你太犀利了

haineiwushuang · 发表于 2012-11-26 14:07:33

厉害，谢谢

gongzhenyu · 发表于 2012-11-28 00:02:48

这个好，应力释放，一直都是纠结的一个问题。

vwxy295 · 发表于 2013-2-24 17:47:32

不错，顶一个

xiaojun · 发表于 2013-3-27 14:18:07

谢谢！

leonandg · 发表于 2013-4-15 15:43:39

学渣苦，学渣累，学渣学习了

过客1314 · 发表于 2014-10-13 11:51:13

论坛真是好东西呀！

yamcsu · 发表于 2014-10-20 11:49:17

qiandao 发表于 2008-6-6 14:27
姚兄，类似ANSYS，在施作二衬时，应该把加的这个反力撤除（remove）吧？不过撤除的时候我老出错，还麻烦指 ...

请问你现在明白了么？？我也不明白，麻烦指点一下。

京京宋007 · 发表于 2015-8-2 11:14:36

但是我做出来的结果怎么感觉不对呢

cumtdenglei · 发表于 2015-9-21 11:31:52

不太懂，再看看

账号		自动登录	找回密码
密码			注册

[前后处理] flac3d中应力释放的实现（总结加实现）

本帖子中包含更多资源

浏览过的版块