FLAC3D动力计算中的问题

magicdonkey

图1是一个城门洞型的隧洞开挖完毕的破坏区图2是是施加完自由场边界后的模型示意图
图3是在模型底部输入加速度时程后的计算5s所得的破坏区分布
图4是中间一片模型的破坏区分布

请问为何用FLAC3D计算后，得到的破坏区是如此一个规律呢？

yinyunteng

free x y z ran z -133.83 -1就可以了
free y ran y -10.1 -9.9   
free y ran y 19.9 20.1   
free x ran x -105 -104
free x ran x  104  105就不需要了吧？

magicdonkey

顶起!希望得到解答

magicdonkey

不是我恶意灌水，只是希望能引起大家的关注，所以再次顶起！
谢谢指教！

深夜来访电话

开始学习动力计算，好没有入门

aizhongluo

希望楼主贴出源代码，那样才能分析原因

ZZYWLS

这种情况可能是边界条件没有设置正确。。。。引起大面积破坏

猪猪加油

看不到，但是也支持一下

yinyunteng

参数设置对了没有？最好把程序贴出来，大家看看。

magicdonkey

附上动力计算步骤代码，请大家斧正！
new

rest Exca-Range1st.sav  ;开挖完毕的SAV文件

set dyn on
free x y z ran z -133.83 -1
free y ran y -10.1 -9.9   
free y ran y 19.9 20.1   
free x ran x -105 -104
free x ran x  104  105

table 1 read EL-CENTRO-filtered.dat
apply xacc 1.0 hist table 1 range z -133.83 -133.63
apply yacc 0.85 hist table 1 range z -133.83 -133.63  
apply zacc 0.65 hist table 1 range z -133.83 -133.63

apply ff

set dyn damp local 0.157
hist dytime    ;1
hist gp disp   0     -10.0 -133.73   ;2
hist gp disp   0     -10.0   28.97   ;3   
hist gp disp  -12.02 -10.0    0      ;4
hist gp disp   0     -10.0  -19.25   ;5
hist zone smax 0     -10.0   28.97   ;6
hist zone smin 0     -10.0   28.97   ;7
hist zone smax -12.02 -10.0  0       ;8
hist zone smin -12.02 -10.0 0        ;9
hist zone smax 0 -10.0 -19.25        ;10
hist zone smin 0 -10.0 -19.25        ;11

set dyn time 0
set dyn multi on
solve age 5   ;帖子中的截图是计算到2.2时刻的塑性区截图。

liugang1986

不懂
帮顶

magicdonkey

11# yinyunteng 后面的四行是因为这段代买有时候我也用于静态边界设置，所以在使用自由场边界时就没有删掉，因为施加APPLY FF后，四周的边界条件自动转换为了自由状态。所以这四行在这里应该没有作用。

yzhanguo

不是很懂
没有研究过动力分析
看来是落伍了

xnjdlhy

FALC 3D能动力计算啊！