flac对几种阻尼比较.docVIP

对几种阻尼的比较在考虑阻尼的影响时，我做了一个很简单的例子，就是10m×10m×10m的立方块，每个方向划分为10个单元，在底面输入正弦波，和以前给你看的一样。文件如下，我分了几种情况：config dynamicdef wavewave=ampl*sin(2.0*pi*freq*dytime)endset freq=4 ampl=0.1386wavegen zone brick p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 size 10 10 10 m eprop bulk 52e6 shear 11e6ini dens 1760;固定底边界fix z range z -0.01 0.01set dyna offset grav 0 0 -10hist unbalhist gp zvel 5 5 5hist gp zdis 5 5 10plot bcon szzplot add hist 1 redcy 500save static.sav;----------------动力计算---------------;动力计算前，应该先清零ini xdisp=0 ydisp=0 zdisp=0ini xvel=0 yvel=0 zvel=0set dyn on;-----------侧面边界条件----------------------fix x range x -0.01 0.01fix x range x 9.99 10.01fix y range y -0.01 0.01fix y range y 9.99 10.01;------模型底面动力输入－－－－－－－－－－－－app xacc=1 his wave range z -0.01 0.01;---------设置阻尼条件-----------;set dyn damp rayleigh 0.1 3;set dyn damp local 0.2;ini damp hyst hardin 0.234;－－－选取的几个比较典型的节点进行记录－－－－－hist resethist unbal hist dytimeplot con sxxplot add hist 1 red;-------x direction acceleration-------his gp xacc 5 5 0 ；底部 his gp xacc 5 5 5 ；中部 his gp xacc 5 5 10 ；顶部solve age 1save dynamic.sav1、无阻尼时：波形很混乱2、瑞利阻尼：已算出模型一阶自振频率为3Hz，所以有：set dyn damp rayleigh 0.1 3，后面的3Hz为模型的一阶自振频率，但前面的参数ξmin该如何取值？这个的结果还可以接受和解释。我监测了沿竖向方向的各个点（z＝0、1、2、3、4、5、6、7、8、9、10）的响应，响应在0～7之间呈衰减趋势，但在7～10之间放大了，原因为土体自由面引起的波的反射叠加使响应增加的缘故。3、局部阻尼在瑞利阻尼中ξmin＝0.08的计算结果比较能接受，因此由手册中提供的公式有：αL＝pi×D=0.25，即set dyn damp local 0.25，响应如无阻尼一样，很杂乱。我也试算过别的阻尼值，响应都很杂乱。这里的问题同样是参数值如何选取？4、滞后阻尼响应同样很杂乱。采用不同的模型（default、S3、S4与Hardin）结果差别不大。从上可以看出：瑞利阻尼的优点在于保持波形，而局部阻尼和滞后阻尼则不能；但瑞利阻尼的最明显不足在于他显著降低了计算时步。以上例子中无阻尼和局部阻尼以及滞后阻尼时的时步均为1.142e－3，而瑞利阻尼随ξmin的增加而降低，ξmin＝0.08时的时步为5.423e－4问题：1、各阻尼中参数如何取值？2、出瑞利阻尼外其余的集中情况波形杂乱是什么原因？应该还有其他的原因？如波在介质中的传播的问题（关于这点我也正在查看相关资料）或我在动力分析阶段仍采用的是弹性模型？还请指点！附上频率计算：config dynamicgen zone brick p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 size 10 10 10 m eprop bulk 52e6 shear 11e6ini dens 1760;固定底边界fix z range z -0.01 0.01set dyna onset grav 0 0 -10hist unbalhist dytimehist gp zvel 5 5 10hist gp zdis 5 5 10plot bcon szz