flac3d动力分析的一些问题，请指教？

xjiang01

如下，请指教．
１．如何确切获得地震加速度－时间曲线？或何处有详细的地震加速度－时间数据可供参考引用？
２．动力分析时，对网格剖分有什么特殊要求？
３．若施加ff（free-field）边界，如何在后处理时将ff边界去除，即只显示实际模型？
４．假定地震持续时间为２０秒，则如何获得地震过程中，如第１５秒时刻点的力学响应图形或数值？
５．局部阻尼，瑞利阻尼的几个参数是如何具体确定的？即
    set dyn damp local 0.125              ;0.125 是如何确定的？
    set dyn damp rayleigh 0.05 22.8    ;0.05 22.8这两个数值是如何确定的？
６．施加ff边界或quiet边界时，有无必要先解除之前静力分析时设定的约束？
　　即：free x y z ran z -0.1 0.1  
　　　　apply ff　或　apply nquiet squiet dquiet range z -0.1 0.1
　　语句free x y z ran z -0.1 0.1 是否一定需要书写？apply　ff或apply nquiet squiet dquiet range z -0.1 0.1是否意味着会将之前静力分析时设定的约束自动解除，而施加ff或quiet边界条件？
７．进行动力分析前，需要如下语句：
　　ini xdis 0 ydis 0 zdis 0
　　ini xvel 0 yvel 0 zvel 0
　　ini state 0
　　请问语句ini　state　0是否必须书写？
８．flac3d是如何考虑动本构的？除了液化分析的finn模型外．
９．如果模型中增添结构性单元，然后再进行地震动力分析，请问flac3d是如何考虑土-结构动力相互作用的？

以上问题，盼指点，谢谢！

[ 本帖最后由 benjackxu 于 2006-12-4 16:46 编辑 ]

lakewater

xjiang01 wrote:
谢谢lakewater的热情指点．
尚有几个问题需要讨论：
２．不赞同斑竹的观点，动力分析时网格有特殊要求，见＂利用FLAC3D 分析某边坡地震稳定性＂（岩石力学与工程学报）；
４．这样做，那意味着动力计算分别进行了两次，计算耗时将急剧增加；
５．能否更详细描述这些阻尼参数的选取；
６．土－结构动力相互作用究竟具体如何操作？理论背景？
谢谢．

2. 这个要求是根据波动理论中来的，手册中讲得很清楚，很多有限元程序也会有这个要求，其实就是最小尺寸的要求；

3. 不对。sol age 15后再sol age 20，其dytime是从15s开始的，不存在什么“两次”的问题；

4. 阻尼参数还是先看手册的好，我也没有什么好的意见；

6. 操作：就是结构单元+ config dyn，理论背景其实还是flac的计算原理

starsmoon

楼主看来思考得很深刻啊，对于动力问题我也比较迷茫啊，向你请教。

xjiang01

starsmoon wrote:
楼主看来思考得很深刻啊，对于动力问题我也比较迷茫啊，向你请教。

呵呵，互相学习，共同进步．

lakewater

我来试着做答吧。
1. 地震反应分析一般有甲方提供地震波型，或者根据已有的波型（比如唐山波）按照分析要求，对最大加速度、历时、卓越周期进行调整；或者按照给出的关键参数进行人工合成地震波，这方面已有相关的研究介绍；

2. 没什么特殊要求；

3. ff边界是没有分组的，所以后处理可以这样
    plo con sz ran group 1 any group 2 any ...
    把所有已划分的组全部填上，就可以不显示ff边界了
这是笨方法！

4. 先sol age 15，然后save，或者进行后处理
   然后再sol age 20

5. 这个现在主要还是参考手册，瑞利阻尼的参数比较容易确定，因为很多有限元程序也是采用这种阻尼参数的。局部阻尼是参考了静力计算的阻尼公式；

6. 这个应该是必须的，后面的约束解除也是必须的；

7. 这个ini sate 0是将所有单元的“塑性指示器”归0，以表示没有塑性破坏；

8. Finn主要考虑液化，若不考虑液化则可以按照静力问题选择本构模型；

9. 手册中明确说明可以将结构单元进行动力分析，以进行土与结构的相互作用研究。

主要问题说清楚了，但我也有很多疑惑，希望大家共同交流！

xjiang01

谢谢lakewater的热情指点．
尚有几个问题需要讨论：
１．能否提供有＂确切数据＂的地震波（有可供引用的文献），而不是地震加速度－时间曲线；
２．不赞同斑竹的观点，动力分析时网格有特殊要求，见＂利用FLAC3D 分析某边坡地震稳定性＂（岩石力学与工程学报）；
３．此法甚好用，谢谢．
４．这样做，那意味着动力计算分别进行了两次，计算耗时将急剧增加；
５．能否更详细描述这些阻尼参数的选取；
６．土－结构动力相互作用究竟具体如何操作？理论背景？

谢谢．

xjiang01

另动力分析时材料参数是否需要将静力分析时的材料参数进行适当调整？

starsmoon

先收藏起来，等学习动力部分再来学习。

xjiang01

lakewater wrote:
2. 这个要求是根据波动理论中来的，手册中讲得很清楚，很多有限元程序也会有这个要求，其实就是最小尺寸的要求；

3. 不对。sol age 15后再sol age 20，其dytime是从15s开始的，不存在什么“两次”的问题；

4. 阻尼参数还是先看手册的好，我也没有什么好的意见；

6. 操作：就是结构单元+ config dyn，理论背景其实还是flac的计算原理

关于土－结构动力相互作用，理论背景是flac的计算原理．但究竟是如何处理的？还请斑竹明示．谢谢．即是否结构性单元也设置阻尼等？
另关于网格，不是最小尺寸，应当是最大尺寸要求，估计斑竹笔误．
再次感谢斑竹对此问题的关注！

freddie_1999

我的看法,请大侠指点

2.对于中心差分法,因其是显式差分,对网格有比较严格的要求,称作为'有条件收敛',它的时间步的选择由最小的单元尺寸决定,临界时间步为delta(T)=L/C(其中,L-单元的最小边长,C-弹性波的速度).而对于隐式差分,比如NEWMARK法,则称作为'无条件收敛'.
   在FLAC中,单元的边长由弹性波波长决定
5.阻尼分为内阻尼和外阻尼,工程中一般用rayleigh做近似.FLAC中的rayleigh要求两个参数,分别为最小阻尼比和最小频率(不是圆频率),阻尼比一般由经验得到,岩土介质的阻尼比的范围为0.02--0.05,如果考虑了结构单元,结构单元的阻尼在0.02--0.10.
   至于最小频率的获得比较麻烦,最笨的方法是先对模型做弹性无阻尼的计算,记录某处的速度或位移,由曲线得到最小频率,然后在恢复正常计算,其实这样的过程就是有限元软件中的模态的计算过程.
   局部阻尼实质和动力阻尼是不同的,它的作用仅仅是为了保证静力计算的收敛.
7.　　ini xvel 0 yvel 0 zvel 0这句是必须的

其他的观点同LAKEWATER斑竹

riken

１．能否提供有＂确切数据＂的地震波（有可供引用的文献），而不是地震加速度－时间曲线；
去美国加州伯克力大学地震研究中心http://nisee.berkeley.edu/去看看，
另外princeton和bafullo两个大学也有相关的研究机构
http://www.princeton.edu/~dynaflow/vitae/Prevost.htm  
在美国的各类相关著名实验室都有相关数据或试用软件可供下载

本人也做动力问题，希望与各位高手交流学习 30726262

believewhy

原帖由 freddie_1999 于 2005-9-12 15:11 发表
我的看法,请大侠指点

2.对于中心差分法,因其是显式差分,对网格有比较严格的要求,称作为'有条件收敛',它的时间步的选择由最小的单元尺寸决定,临界时间步为delta(T)=L/C(其中,L-单元的最小边长,C-弹性波的速度).而 ...

请问如果动力计算的时候存在结构单元，如何设置结构单元的阻尼比，通过什么命令来设置阻尼参数？

zhuhaha

晕，动力计算这么长时间？都说FLAC计算非线性比有限元快，那么计算动力问题，它和LS-DYNA相比哪个更好呢？

lbt18872

没有十全十美的事物哦，大家在享受其带来便利的一方面的同时，也应该想到起不便利的一些地方，这些都很正常的哦，呵呵

时间的主人

看了上面的帖子,收益匪浅.我也有个问题:
1.我在做一个隧道动力开挖的例子,其中也设置了锚杆单元,那么是否也需要为锚杆单元这一部分单独设定:set dyn,是否也需要为锚杆部分单独设定瑞雷阻尼比?

freddie_1999

原帖由 时间的主人 于 2007-1-14 13:20 发表
看了上面的帖子,收益匪浅.我也有个问题:
1.我在做一个隧道动力开挖的例子,其中也设置了锚杆单元,那么是否也需要为锚杆单元这一部分单独设定:set dyn,是否也需要为锚杆部分单独设定瑞雷阻尼比?

需要的，但是对于频率的取值大小可能需要进一步讨论

zhuhaha

看一遍学习一遍，多看一遍，就多学习一遍！！！！！！！！！！！！！！！！！

flankness

我的动力timestep＝3.639e-9，怎么那么小啊，如果set dytime＝20 s，那不是要算很久吗？
怎么改变timestep？觉得我的网格单元不是很小啊，单元总数也才700多个。请帮帮忙？
我是用轴对称分析的

[ 本帖最后由 flankness 于 2007-5-10 11:44 编辑 ]

suizhirong

看得头都大了，不过确实学到不少东西。问个问题：滞后阻尼怎么取参数？