振动力学:无阻尼单自由度系统
单自由度振动系统是最简单的一类振动系统,仅用一个坐标就可描述。从力学角度分析,一个实际的振动系统可由三个元件组成:惯性质量、弹性、阻尼,,它们分别描述系统的惯性、弹性、能耗机制。惯性元件是运动的实体,弹性元件提供振动回复力,阻尼元件在振动过程中消耗或吸收外界能量。
1. 单自由度系统振动方程
任何具有惯性和弹性的系统均可产生振动,单自由度系统是对这种系统的高度抽象概括。例如升降机吊篮、列车的一节车厢、高楼的一层、弹性体中的某点在某方向的振动,均可用单自由度系统振动模型描述,如图1所示。
惯性质量、弹性元件、阻尼元件的三个参数分别为:质量 m m m、刚度 k k k、粘性阻尼系数 c c c。特别地,采用线性阻尼可使运动方程得到简化,假设阻尼力与阻尼元件两端的速度成正比,即阻尼力 = c u ˙ ( t ) c\dot{u}(t) cu˙(t)。
图1 单自由度系统振动模型
以质量元件在平衡位置的坐标为坐标原点 O O O,对质量元件进行受力分析,根据牛顿第二定律可得:
m u ¨ ( t ) = − k [ u ( t ) + δ s ] − c u ˙ ( t ) + m g + f ( t ) ( 1.1 ) m \ddot{u}(t) = - k \left[ u(t) + \delta_s \right] - c\dot{u}(t) +mg + f(t) \qquad(1.1) mu¨(t)=−k[u(t)+δs]−cu˙(t)+mg+f(t)(1.1)
另一方面,根据静力平衡:
m g = k δ s mg = k\delta_s mg=kδs
于是式(1.1)进一步写为:
m u ¨ ( t ) + c u ˙ ( t ) + k u ( t ) = f ( t ) ( 1.3 ) m \ddot{u}(t) + c\dot{u}(t) +k u(t) = f(t) \qquad(1.3) mu¨(t)+cu˙(t)+ku(t)=f(t)(1.3)
式中 f ( t ) f(t) f(t)和 δ s \delta_s δs分别表示外力(方向设为与重力方向相同)和处于初始状态时的弹性体伸长量。
这就是单自由度系统振动方程的一般形式。它是二阶常系数线性非齐次微分方程。
若不受外力作用,则 f ( t ) = 0 f(t)=0 f(t)=0,于是单自由度系统自由振动方程为:
m u ¨ ( t ) + c u ˙ ( t ) + k u ( t ) = 0 ( 1.4 ) m \ddot{u}(t) + c\dot{u}(t) +k u(t) = 0 \qquad(1.4) mu¨(t)+cu˙(t)+ku(t)=0(1.4)
进一步,若 c = 0 c = 0 c=0,无阻尼自由振动方程为:
m u ¨ ( t ) + k u ( t ) = 0 ( 1.5 ) m \ddot{u}(t) +k u(t) = 0 \qquad(1.5) mu¨(t)+ku(t)=0(1.5)
2. 无阻尼单自由度振动:特征解和通解
式(1.5)是单自由度振动最简单的方程,根据微分方程知识,其解有如下形式:
u ( t ) = u ˉ e s t u(t) = \bar{u} \mathrm{e}^{st} u(t)=uˉest
式中, u ˉ \bar{u} uˉ和 s s s为常量。将上式代入式(1.5),得到:
( m s 2 + k ) u ˉ = 0 (ms^2 +k)\bar{u} =0 (ms2+k)uˉ=0
若系统存在非零解(即振动位移不为零),则满足条件:
m s 2 + k = 0 ( 2.4 ) ms^2 + k =0\qquad(2.4) ms2+k=0(2.4)
式(2.4)以 s s s为变量,称为特征方程,它的解 s 1 , 2 = ± j ω n s_{1,2} = \pm j\omega_{\rm n} s1,2=±jωn称为特征根。系统的固有圆频率,简称为固有频率,定义为:
ω n = k m ( 2.5 ) \omega_{\rm n} = \sqrt{\frac{k}{m}} \qquad (2.5) ωn=mk(2.5)
式中,各变量的单位为 k = [ N / m ] , m = [ k g ] , ω n = [ r a d / s ] k=\rm{[N/m]}, m=\rm{[kg]}, \omega_{\rm n} =\rm{[rad/s]} k=[N/m],m=[kg],ωn=[rad/s]。
方程(1.5)有两个特征解,即为 u ˉ 1 e s 1 t , u ˉ 2 e s 2 t \bar{u}_1 \mathrm{e}^{s_1 t}, \bar{u}_2 \mathrm{e}^{s_2 t} uˉ1es1t,uˉ2es2t,根据线性系统叠加原理,通解为两个特征解的线性叠加:
u ( t ) = a 1 cos ( ω n t ) + a 2 sin ( ω n t ) ( 2.6 ) u(t) = a_1 \cos(\omega_{\rm n} t) + a_2 \sin(\omega_{\rm n} t) \qquad (2.6) u(t)