基于matlab的阻尼振动.pdf

基于matlab 的阻尼振动 孟利平 （船院 2009011521） 摘要： 分析单自由度阻尼系统的阻尼系统对其固有振动模态的影响，用传统方 法通过matlab实现其振动规律曲线，并改变其初始条件，实现一组曲线的脉 冲过度函数，改变传统运算算法，使运算简便明了。 关键字： 阻尼振动，复数运算，简便 引言： 通过传统方法，用基本公式分析阻尼振动，并根据传统算法编程，用 matlab实现其振动规律曲线，并绘制曲线。由于传统计算公式要考虑复数运 算，其要避开复数运算，所以程序显得繁琐，改变其中算法公式，用matlab 直接实现其算法，并绘制相应曲线。改变其初始条件，实现一组曲线的脉冲 过度函数，并用matlab绘制其过渡曲线。 基本知识： 1 不论是弹簧振子还是单摆，由于外界的摩擦和介质阻力总是存在，， 在振动过程中要不断克服外界阻力做功，消耗能量，振幅就会逐渐减小， 经过一段时间，振动就会完全停下来。这种振幅越来越小的振动叫做阻 尼振动。 振幅随时间减小的振动称为阻尼振动．因为振幅与振动的能量有关， 阻尼振动也就是能量不断减少的振动．阻尼振动是非简谐运动．阻尼振 动系统属于耗散系统。 阻尼振动 振动系统因受阻力而作振幅减小的运动。 一 . 阻尼振动的动力学方程 现象：振幅随时间减小 原因：阻尼 动力学分析： 假设：振动速度较小时，摩擦力正比于质点的速率。即： 假设：振动速度较小时，摩擦力正比于质点的速率。即： 阻尼力 F Cv 阻力系数 C r 对物块应用牛顿第二定律： 2 kxCv ma 2 d x dx m  C kx 0 dt2 dt 化解为二阶线性常系数齐次方程，即阻尼振动的动力学方程： 2 d x dx 2 2  x 0 dt2 dt 0 k 其中： 固有角频率 0 m 阻尼系数  C 2m 解得： t x Ae cos(t) 其中： 振幅 Aet 角频率  2 2 2π 2 2    T 2π   0 0  Aet   振幅 为随时间的推移，呈指数递减， 越大，振动衰减越快； 越小，振幅衰减越慢。随时间的推移，质点坐标单调地趋于零。质点运动是 非周期的，甚至不是往复的。将质点移开平衡位置后释放，质点便慢慢