flac对几种阻尼的比较

《flac对几种阻尼的比较》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、对几种阻尼的比较在考虑阻尼的影响时，我做了一个很简单的例子，就是10m×10m×10m的立方块，每个方向划分为10个单元，在底面输入正弦波，和以前给你看的一样。文件如下，我分了几种情况：configdynamicdefwavewave=ampl*sin(2.0*pi*freq*dytime)endsetfreq=4ampl=0.1386wavegenzonebrickp0000p11000p20100p30010size101010mepropbulk52e6shear11e6inidens1760;固定底边界fixzrangez-0.010.01s

2、etdynaoffsetgrav00-10histunbalhistgpzvel555histgpzdis5510plotbconszzplotaddhist1redcy500savestatic.sav;----------------动力计算---------------;动力计算前，应该先清零inixdisp=0ydisp=0zdisp=0inixvel=0yvel=0zvel=0setdynon;-----------侧面边界条件----------------------fixxrangex-0.010.01fixxrangex9.9910

3、.01fixyrangey-0.010.01fixyrangey9.9910.01;------模型底面动力输入－－－－－－－－－－－－appxacc=1hiswaverangez-0.010.01;---------设置阻尼条件-----------;setdyndamprayleigh0.13;setdyndamplocal0.2;inidamphysthardin0.234;－－－选取的几个比较典型的节点进行记录－－－－－histreset5histunbalhistdytimeplotconsxxplotaddhist1red;-------

4、xdirectionacceleration-------hisgpxacc550；底部hisgpxacc555；中部hisgpxacc5510；顶部solveage1savedynamic.sav1、无阻尼时：波形很混乱2、瑞利阻尼：已算出模型一阶自振频率为3Hz，所以有：setdyndamprayleigh0.13，后面的3Hz为模型的一阶自振频率，但前面的参数ξmin该如何取值？5这个的结果还可以接受和解释。我监测了沿竖向方向的各个点（z＝0、1、2、3、4、5、6、7、8、9、10）的响应，响应在0～7之间呈衰减趋势，但在7～10之间放大了，

5、原因为土体自由面引起的波的反射叠加使响应增加的缘故。3、局部阻尼在瑞利阻尼中ξmin＝0.08的计算结果比较能接受，因此由手册中提供的公式有：αL＝pi×D=0.25，即setdyndamplocal0.25，响应如无阻尼一样，很杂乱。我也试算过别的阻尼值，响应都很杂乱。这里的问题同样是参数值如何选取？54、滞后阻尼响应同样很杂乱。采用不同的模型（default、S3、S4与Hardin）结果差别不大。从上可以看出：瑞利阻尼的优点在于保持波形，而局部阻尼和滞后阻尼则不能；但瑞利阻尼的最明显不足在于他显著降低了计算时步。以上例子中无阻尼和局部阻尼以及滞

6、后阻尼时的时步均为1.142e－3，而瑞利阻尼随ξmin的增加而降低，ξmin＝0.08时的时步为5.423e－4问题：1、各阻尼中参数如何取值？2、出瑞利阻尼外其余的集中情况波形杂乱是什么原因？应该还有其他的原因？如波在介质中的传播的问题（关于这点我也正在查看相关资料）或我在动力分析阶段仍采用的是弹性模型？还请指点！5附上频率计算：configdynamicgenzonebrickp0000p11000p20100p30010size101010mepropbulk52e6shear11e6inidens1760;固定底边界fixzrangez-0

7、.010.01setdynaonsetgrav00-10histunbalhistdytimehistgpzvel5510histgpzdis5510plotbconszzplotaddhist1redsolveage1savefreq.sav;3Hz5