flac对几种阻尼的比较.docVIP

对几种阻尼的比较 在考虑阻尼的影响时，我做了一个很简单的例子，就是10m×10m×10m的立方块，每个方向划分为10个单元，在底面输入正弦波，和以前给你看的一样。文件如下，我分了几种情况： config dynamic def wave wave=ampl*sin(2.0*pi*freq*dytime) end set freq=4 ampl=0.1386 wave gen zone brick p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 size 10 10 10 m e prop bulk 52e6 shear 11e6 ini dens 1760 ;固定底边界 fix z range z -0.01 0.01 set dyna off set grav 0 0 -10 hist unbal hist gp zvel 5 5 5 hist gp zdis 5 5 10 plot bcon szz plot add hist 1 red cy 500 save static.sav ;----------------动力计算--------------- ;动力计算前，应该先清零 ini xdisp=0 ydisp=0 zdisp=0 ini xvel=0 yvel=0 zvel=0 set dyn on ;-----------侧面边界条件---------------------- fix x range x -0.01 0.01 fix x range x 9.99 10.01 fix y range y -0.01 0.01 fix 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 reset hist unbal hist dytime plot con sxx plot 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 1 save dynamic.sav 1、无阻尼时：波形很混乱 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 dynamic gen zone brick p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 size 10 1