跳转至

线性振动

线性振动是物理学中一种重要的运动形式,广泛存在于各种物理系统中。它通常指系统在平衡位置附近的小幅度振动,其运动方程可以用线性微分方程描述。

简单谐振动

如果某一物量在平衡点附近往复振动,且随时间的变化是余弦(或正弦)的函数形式,我们就把这种振动称为简谐振动。

运动方程推导

简单谐振动是最基本的线性振动形式,其运动方程为:

d2xdt2+ω2x=0

这一方程可以通过牛顿第二定律推导得到。例如,对于一个质量为 m 的物体,受到的回复力 F 满足胡克定律:

F=kx

根据牛顿第二定律:

F=md2xdt2

代入回复力表达式:

md2xdt2=kx

整理得:

d2xdt2+kmx=0

ω=km,即得简单谐振动的运动方程。

解的形式

该方程的通解为:

x(t)=Acos(ωt+ϕ)

其中: - A 为振幅,振子离开平衡位置的最大位移的绝对值,由初始条件决定。是求解上述微分方程时引入的常量; - ϕ 为初相位,初始时刻振子的相位,由初始条件决定。是求解上述微分方程时引入的常量; - ω 为角频率(圆频率),在简谐振动的几何表示中,表示旋转矢量的角速度; - x 为位移,表示振子相对于平衡位置的偏移量。

例1

设一个质量为 m 的物体连接在弹簧上,弹簧劲度系数为 k。求该系统的角频率 ω 及其运动方程。

解: 易得运动方程为:

d2xdt2+kmx=0

角频率:

ω=km

因此,运动方程为:

x(t)=Acos(kmt+ϕ)

其中 Aϕ 由初始条件决定。

例2

设一个质量为 m 的物体连接在摆绳上,摆长为 l。求该单摆系统的角频率 ω 及其运动方程。

解:

M=mglsinθ=Iβ=ml2d2θdt2

对于小角度摆动, $$ \sin\theta \approx \theta $$,单摆的运动方程为:

d2θdt2+glθ=0

角频率:

ω=gl

因此,运动方程为:

θ(t)=Acos(glt+ϕ)

其中 Aϕ 由初始条件决定。

小角度近似

在单摆的例子中,我们使用了小角度近似 $$ \sin\theta \approx \theta $$。这一近似仅在摆动角度较小时成立(通常小于约 10 度)。对于较大的摆动角度,运动方程将变得非线性,不能再用简单谐振动的形式描述。

简谐振动的判断

并非所有往复振动都是简谐振动。简谐振动要求回复力与位移成正比且方向相反,即满足胡克定律的形式。对于非线性系统或大幅度振动,运动方程可能不再是线性的,不能用简谐振动的形式描述。我们往往也可以通过得到相应的特征微分方程来判断振动是否为简谐振动。

周期与频率

振动的周期 T 和频率 f 分别为:

T=2πω,f=1T=ω2π.

相位

振动的相位 θ 定义为:

θ=ωt+ϕ.

相位表示振动在某一时刻所处的位置状态。规定 0<θ<2π

对于两个同频简谐振动,其相位差 Δθ 为:

Δθ=θ1θ2=ϕ1ϕ2.

Δθ=0,则两振动同相,若 Δθ=π,则两振动反相。 若 0<Δθ<π,则称 x1x2 超前,若 π<Δθ<2π,则则称 x1x2 落后

x, v, a 三者的相位关系

x(t)=Acos(ωt+ϕ)v(t)=Aωcos(ωt+ϕ+π2)=Aωsin(ωt+ϕ)a(t)=Aω2cos(ωt+ϕ+π)=Aω2cos(ωt+ϕ)

由此可见,速度 v 比位移 x 超前 π2,加速度 a 与位移 x 反相。你可以想想为什么会有这样的相位关系。

能量分析

动能与势能

简单谐振动的能量包括动能和势能,其表达式分别为:

  • 动能:
Ek=12mv2=12mω2A2sin2(ωt+ϕ)

在一个周期内的平均值:

Ek=1T0TEkdt=14kA2
  • 势能:
Ep=12kx2=12kA2cos2(ωt+ϕ)

在一个周期内的平均值:

Ep=1T0TEpdt=14kA2

总能量守恒

总能量为动能与势能之和:

E=Ek+Ep=12kA2=12mω2A2

总能量在振动过程中保持不变,动能和势能在振动过程中相互转化。

简谐振动的几何表示

简谐振动可以通过旋转矢量的投影来几何表示。设一个长度为 A 的旋转矢量以角速度 ω 绕原点逆时针旋转,其在水平轴上的投影即为简谐振动的位移 x(t)

简正模

在多自由度系统中,线性振动可以分解为若干个独立的简正模(Normal Modes)。每个简正模对应一个特定的频率和振动模式,系统的总振动可以看作这些简正模的叠加。

用复数表示简谐振动

简谐振动

x(t)=Acos(ωt+φ)

也可以用复数的实部和虚部表示:

x~(t)=Aei(ωt+φ) 

其中,A 为振幅,ω 为角频率,φ 为初相位。

上式的右端又可以写为 (Aeiφ)eiωt=A~eiωt,其中 Aeiφ=A~ 是一个复常数,表示振动的初始状态,称为复振幅

x~(t)=A~eiωt

如果我们取上式的实部,就得到了简谐振动的标准形式:

x(t)=Re(x~(t))=Re(A~eiωt)=Acos(ωt+φ)

x~(t) 表示简谐振动,则速度和加速度为

{v~(t)=dx~(t)dt=iωA~eiωt=iωx~(t)a~(t)=dv~(t)dt=d2x~(t)dt2=ω2A~eiωt=ω2x~(t)
例题

已知一个线性三原子分子 A2B 的纵向振动模型:质量为 m 的两个 A 原子位于两端,质量为 MB 原子位于中间;相邻原子间由劲度系数均为 k 的弹簧连接(ABBA)。求简正模与对应角频率。

解:设三原子沿直线位移为 x1,x2,x3(分别对应 A,B,A),取平衡位置为零。 弹簧形变为 x1x2x3x2,于是运动方程为

{mx¨1=k(x1x2),Mx¨2=k(x2x1)k(x2x3)=k(2x2x1x3),mx¨3=k(x3x2).

取简正模形式 xj(t)=ajeiωt,得到代数方程

{mω2a1=k(a1a2),Mω2a2=k(2a2a1a3),mω2a3=k(a3a2).

利用对称性分两类:

1) 整体平移模a1=a2=a3,此时弹簧不伸长,ω0=0

2) 反对称伸缩模a1=a3,由对称性得 a2=0。 代入第一式得 mω2a1=ka1,因此

ω=km,(a1,a2,a3)(1,0,1).

3) 对称伸缩模a1=a3=aa2=b。 由第一式 mω2a=k(ab),第二式 Mω2b=2k(ba)。 联立消去 b 可得

ω+2=km+2kM,(a1,a2,a3)(1,2mM,1).

因此该分子有一个零频平移模与两个非零振动模 ω,ω+

例题:受迫振动中“功率共振”的频率

对线性受迫振动 $$ \ddot x+2\gamma\dot x+\omega_0^2x=\frac{F_0}{m}\cos(\Omega t), $$ 证明稳态下平均耗散功率 P 取最大值时的驱动频率为 Ω=ω0

解:稳态响应 x=Acos(Ωtδ)。 阻尼耗散功率 Pd=bx˙2,周期平均

P=bx˙2=b12A2Ω2.

b=2mγ 与 $$ A2=\frac{(F_0/m)2}{(\omega_02-\Omega2)2+(2\gamma\Omega)2} $$

得 $$ \langle P\rangle\propto \frac{\Omega2}{(\omega_02-\Omega2)2+(2\gamma\Omega)^2}. $$

Ω 求极值可得最大值发生在 Ω=ω0

阻尼振动

运动方程

当系统受到阻尼力作用时,振动会逐渐衰减,其运动方程为:

d2xdt2+2βdxdt+ω2x=0

其中 β 为阻尼系数,描述阻尼的强弱。

解的形式

根据阻尼的大小,解的形式不同:

  1. 欠阻尼振动β<ω):
x(t)=Aeβtcos(ωdt+ϕ)

其中 ωd=ω2β2 为阻尼振动的角频率。

图例

欠阻尼振动

  1. 临界阻尼β=ω):
x(t)=(C1+C2t)eβt
图例

临界阻尼振动

  1. 过阻尼β>ω):
x(t)=C1er1t+C2er2t

其中 r1r2 为两个负实根。

图例

过阻尼振动

能量衰减

欠阻尼振动中,系统的总能量随时间指数衰减:

E(t)=E0e2βt

强迫振动

运动方程

当系统受到周期性外力作用时,其运动方程为:

d2xdt2+2βdxdt+ω2x=F0cos(Ωt)

其中: - F0 为外力的振幅; - Ω 为外力的角频率。

稳态解

稳态解为:

x(t)=Acos(Ωtδ)

其中: - A=F0(ω2Ω2)2+(2βΩ)2 为稳态振幅; - δ=arctan(2βΩω2Ω2) 为相位差。

共振现象

当驱动力频率接近系统的固有频率时,即 Ωω,系统的振幅达到最大值,称为共振

共振条件下的振幅为:

Ares=F02βω

应用

线性振动在工程、物理和日常生活中有广泛的应用,例如:

  • 钟摆的运动:简单谐振动的经典例子;
  • 弹簧-质量系统:描述机械振动;
  • 电路中的交流电振荡:电感-电容回路的振荡;
  • 建筑物的抗震设计:利用阻尼减小振幅。


评论