致lakewater等关于动力计算中的一些问题 ，请指点

司令大哥

如下：请指教和指点：
     1、本人在对一隧道在地震作用的动力计算中出现了如下的错误提示：，究竟什么原因呢？因为前面的静力计算、地震波的调整（基线校正和滤波）、输入地震荷载、设置边界条件和选择阻尼类型都已完成，当计算到solve age=10 时就出现了上述的错误提示。
     为了解决那个问题我首先以为是age设置有误，所以就改变了age，从10→20→5→2→1→0.5→0.015等等都还是提示以上的错误。
     本人百思不得其解，是不是网格尺寸不合要求？还是地震波的选取有误？盼指点，谢谢！
    2、本人还有一个问题，就是地震波的滤波，在陈博士的那本书上（第16章 阪神地震液化大变形分析）有这么一段话：由于地震波每次往返作用的周期大约为0.2s~1.0s，地震作用的频率约为1~5Hz，所以就过滤掉了地震波中频率大于5HZ的成分。但本人选取的地震波（注：汶川波）在20HZ处振幅仍然较强，那是否应选取20HZ作为截断频率呢？盼指点，谢谢！
    3、本人对动力计算中最后一步solve age=的取值有点疑问，即取值的大小如何确定》
    盼指点，谢谢！

xuelei

期待大虾的回复！

cloud328

1、这个可能是你动力计算过程中认为设置了迭代的时间步，动力计算最好不要去改时间步；
2、滤波不要人云亦云，本质就是把一部分高频能量滤掉。你可以首先对你的地震波进行FFT变化，因地震波能量跟振幅的平方成正比，故你可以把FFT变换后的曲线进行平方，然后积分，单位化，就可以看出0-多少赫兹的能量占据总能量的多少，再决定截断频率；
3、solve后面的时间指的是地震波的持续时间。比如EI波好像是54S，你可以选择算54S，也可以选择只选其中的0-35s。如果你选择0-35s就要把这段波拿出来重新滤波和基线校正，否则还是会出现0点位移漂移的现象。

lewisgogo

好东西 下来看看！谢谢楼主

司令大哥

感谢二楼版主的精彩解答
谢谢
谢谢
太谢谢了
我想再请问三楼cloud328
一个 问题，动力计算时，计算一秒（1s），动力响应就非常大，Z方向的位移达到5cm，加速度最大响应值达42m/s2，真是太不可思议了，我采用的是局部阻尼，换了另外一种地震波，响应值还是很大哈！！！
是不是和网格尺寸划分有关？我以前比较过，网格划分的大小对结果是有一定的影响，但一般不会超过2倍左右；
还是和阻尼的选取有关？
请cloud328 帮忙解答一下
谢谢！！！

司令大哥

好呢，马上把命令流贴出来
谢谢二楼和cloud328 ！！！
以下为动力计算的命令流
rest 4开挖支护.sav
config dynamic
set dyn on
ini xvel 0 yvel 0 zvel 0
ini xdis 0 ydis 0 zdis 0
ini state 0;将所有单元的“塑性指示器”归0，以表示没有塑性破坏
free x y z range z -39.9 -40.1
ini  density  2700
ini  density  2600 ran gro 围岩2
ini  density  2400 ran gro 围岩3
ini density 2500 range group 衬砌
set dyn multi on
table 1 read 0.15g加速度.dat
app xacc 1.0 hist table 1 range  z -39.9 -40.1
apply ff

set dyn time=0
set dyn damp local 0.1571
set hist_rep 100
hist reset
hist unbal
hist dytime

call常规动力监测.dat

pl set rot 20 0 30
pl con zdis out on magf 20 ran gr 衬砌 any gr 围岩1 any gr 围岩2 any gr 围岩3 any
pl add hist 1
set movie avi step 200 file donghua_wenchuanbo.avi
movie start ;输出同步动画

solve age 1.0
save 5动力计算_age1.0.sav
solve age 2.0
save 5动力计算_age2.0.sav
solve age 3.0
save 5动力计算_age3.0.sav
solve age 4.0
save 5动力计算_age4.0.sav
solve age 5.0
save 5动力计算_age5.0.sav
solve age 6.0
save 5动力计算_age6.0.sav
solve age 7.0
save 5动力计算_age7.0.sav
solve age 8.0
save 5动力计算_age8.0.sav
solve age 9.0
save 5动力计算_age9.0.sav
solve age 10.0
save 5动力计算_age10.0.sav
solve age 11.0
save 5动力计算_age11.0.sav

todayiswhat

6# 司令大哥

司令大哥，你好厉害啊！我现在也要做地震模拟，感觉真难啊！

张如林

有次我计算静力问题，不小心设置了config dy，就出现了那个警告。
去掉config dy就没事了

嘻嘻

司令大哥，你好，我现在也正在学习关于地震动力分析的模拟，不知大哥可以给我发几个你做过的地震分析的例子吗？谢谢大哥啦！！万分感激！！！

chen575982

为什么在动力计算时，还要设置一下围岩密度？有些不解

www

动力计算式FLAC3D最难学的一部分，也是最大特色之一