flac对几种阻尼的比较.doc

1 对几种阻尼的比较在考虑阻尼的影响时, 我做了一个很简单的例子,就是 10m × 10m × 10m 的立方块,每个方向划分为 10个单元,在底面输入正弦波,和以前给你看的一样。文件如下,我分了几种情况: config dynamic def wave wave=ampl*sin(*pi*freq*dytime) end set freq=4 ampl= wave gen zone brick p0000p11000p20100p30010 size 101010 me prop bulk 52e6 shear 11e6 ini dens 1760 ;固定底边界 fix z range z - set dyna off set grav 00 -10 hist unbal hist gp zvel 555 hist gp zdis 5510 plot bcon szz plot add hist 1 red cy500 save ;---------------- 动力计算--------------- ;动力计算前,应该先清零 ini xdisp=0 ydisp=0 zdisp=0 ini xvel=0 yvel=0 zvel=0 set dyn on ;----------- 侧面边界条件---------------------- fix x range x - fix x range fix y range y - fix y range ;------ 模型底面动力输入------------ app =1 his wave range z - ;--------- 设置阻尼条件----------- ;set dyn damp rayleigh ;set dyn damp local ;ini damp hyst hardin ;---选取的几个比较典型的节点进行记录----- hist reset 2 hist unbal hist dytime plot con sxx plot add hist 1 red ;-------x direction acceleration------- his gp 550;底部 his gp 555;中部 his gp 5510;顶部 solve age 1 save 1 、无阻尼时:波形很混乱 2 、瑞利阻尼: 已算出模型一阶自振频率为 3Hz ,所以有: set dyn damp rayleigh ,后面的 3Hz 为模型的一阶自振频率,但前面的参数ξ min该如何取值? 3 这个的结果还可以接受和解释。我监测了沿竖向方向的各个点( z= 0、 1、 2、 3、 4、 5、 6、 7、 8、 9、 10)的响应,响应在 0~ 7之间呈衰减趋势,但在 7~ 10 之间放大了,原因为土体自由面引起的波的反射叠加使响应增加的缘故。 3 、局部阻尼在瑞利阻尼中ξ min= 的计算结果比较能接受,因此由手册中提供的公式有: α L= pi× D= ,即 set