跳转至

线性振动

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

简谐振动

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

运动方程推导

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

d2xdt2+ω2x=0d^2xdt^2 + omega^2 x = 0

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

F=−kxF = -kx

根据牛顿第二定律:

F=md2xdt2F = md^2xdt^2

代入回复力表达式:

md2xdt2=−kxmd^2xdt^2 = -kx

整理得:

d2xdt2+kmx=0d^2xdt^2 + kmx = 0

令 ω=kmomega = k/m,即得简谐振动的运动方程.

解的形式

该方程的通解为:

x(t)=Acos⁡(ωt+ϕ)x(t) = A cos(omega t + phi)

其中:

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

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

解:
易得运动方程为:

d2xdt2+kmx=0d^2xdt^2 + kmx = 0

角频率:

ω=kmomega = km

因此,运动方程为:

x(t)=Acos⁡(kmt+ϕ)x(t) = A cosleft(kmt + phiright)

其中 AA 和 ϕphi 由初始条件决定.

例 2

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

解:

M=−mglsin⁡θ=Iβ=ml2d2θdt2M = -mglsintheta = Ibeta = ml^2 d^2dt^2

对于小角度摆动,sin⁡θ≈θsintheta ~ theta,单摆的运动方程为:

d2θdt2+glθ=0d^2dt^2 + gltheta = 0

角频率:

ω=glomega = gl

因此,运动方程为:

θ(t)=Acos⁡(glt+ϕ)theta(t) = A cosleft(glt + phiright)

其中 AA 和 ϕphi 由初始条件决定.

小角度近似

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

简谐振动的判断

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

周期与频率

振动的周期 TT 和频率 ff 分别为:

T=2πω,f=1T=ω2π.T = 2omega, quad f = 1T = 2pi.

相位

振动的相位 θtheta 定义为:

θ=ωt+ϕ.theta = omega t + phi.

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

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

Δθ=θ1−θ2=ϕ1−ϕ2.Delta theta = theta_1 - theta_2 = phi_1 - phi_2.

若 Δθ=0Delta theta = 0,则两振动 同相,若 Δθ=πDelta theta = pi,则两振动 反相. 若 0<Δθ<π0 < Delta theta < pi,则称 x1x_1 比 x2x_2 超前,若 π<Δθ<2πpi < Delta theta < 2pi,则称 x1x_1 比 x2x_2 落后.

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

x(t)=Acos⁡(ωt+ϕ)x(t) = A cos(omega t + phi)v(t)=−Aωcos⁡(ωt+ϕ+π2)=Aωsin⁡(ωt+ϕ)v(t) = -Aomega cos(omega t + phi + 2) = Aomega sin(omega t + phi)a(t)=−Aω2cos⁡(ωt+ϕ+π)=−Aω2cos⁡(ωt+ϕ)a(t) = -Aomega^2 cos(omega t + phi + pi) = -Aomega^2 cos(omega t + phi)

由此可见,速度 vv 比位移 xx 超前 π2pi/2,加速度 aa 与位移 xx 反相.你可以想想为什么会有这样的相位关系.

能量分析

动能与势能

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

  • 动能:
Ek=12mv2=12mω2A2sin2⁡(ωt+ϕ)E_k = 12mv^2 = 12momega^2A^2sin^2(omega t + phi)

在一个周期内的平均值:

Ek―=1T∫0TEkdt=14kA2E_k = 1Tint_0^T E_k dt = 14kA^2
  • 势能:
Ep=12kx2=12kA2cos2⁡(ωt+ϕ)E_p = 12kx^2 = 12kA^2cos^2(omega t + phi)

在一个周期内的平均值:

Ep―=1T∫0TEpdt=14kA2E_p = 1Tint_0^T E_p dt = 14kA^2

总能量守恒

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

E=Ek+Ep=12kA2=12mω2A2E = E_k + E_p = 12kA^2 = 12momega^2A^2

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

简谐振动的几何表示

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

简正模

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

用复数表示简谐振动

简谐振动

x(t)=Acos⁡(ωt+φ)x(t) = A cos(omega t + varphi)

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

x~(t)=Aei(ωt+φ) x(t)=Ae^i(omega t+varphi)

其中,AA 为振幅,ωomega 为角频率,φvarphi 为初相位.

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

x~(t)=A~eiωtx(t)=Ae^iomega t

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

x(t)=Re(x~(t))=Re(A~eiωt)=Acos⁡(ωt+φ)x(t)=Re(x(t))=Re(Ae^iomega t)=Acos(omega t+varphi)

用 x~(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)cases v(t)=dx(t)dt=iomega Ae^iomega t=iomega x(t) a(t)=dv(t)dt=d^2x(t)dt^2=-omega^2 Ae^iomega t=-omega^2 x(t) cases
例题

已知一个线性三原子分子 A2BA_2B 的纵向振动模型:质量为 mm 的两个 AA 原子位于两端,质量为 MM 的 BB 原子位于中间;相邻原子间由劲度系数均为 kk 的弹簧连接(AA–BB 与 BB–AA).求简正模与对应角频率.

解: 设三原子沿直线位移为 x1,x2,x3x_1,x_2,x_3(分别对应 A,B,AA,B,A),取平衡位置为零.
弹簧形变为 x1−x2x_1-x_2 与 x3−x2x_3-x_2,于是运动方程为

{mx¨1=−k(x1−x2),Mx¨2=−k(x2−x1)−k(x2−x3)=−k(2x2−x1−x3),mx¨3=−k(x3−x2).cases mddot x_1=-k(x_1-x_2), Mddot x_2=-k(x_2-x_1)-k(x_2-x_3)=-k(2x_2-x_1-x_3), mddot x_3=-k(x_3-x_2). cases

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

{mω2a1=k(a1−a2),Mω2a2=k(2a2−a1−a3),mω2a3=k(a3−a2).cases momega^2 a_1=k(a_1-a_2), Momega^2 a_2=k(2a_2-a_1-a_3), momega^2 a_3=k(a_3-a_2). cases

利用对称性分两类:

  1. 整体平移模:a1=a2=a3a_1=a_2=a_3,此时弹簧不伸长,ω0=0omega_0=0.

  2. 反对称伸缩模:a1=−a3a_1=-a_3,由对称性得 a2=0a_2=0.
    代入第一式得 mω2a1=ka1momega^2 a_1=k a_1,因此

ω−=km,(a1,a2,a3)∝(1,0,−1).omega_- = km,quad (a_1,a_2,a_3)propto(1,0,-1).
  1. 对称伸缩模:a1=a3=aa_1=a_3=a,a2=ba_2=b.
    由第一式 mω2a=k(a−b)momega^2 a=k(a-b),第二式 Mω2b=2k(b−a)Momega^2 b=2k(b-a).
    联立消去 bb 可得
ω+2=km+2kM,(a1,a2,a3)∝(1,−2mM,1).omega_+^2=km+2kM,quad (a_1,a_2,a_3)proptoleft(1,-2mM,1right).

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

例题:受迫振动中「功率共振」的频率

对线性受迫振动

x¨+2γx˙+ω02x=F0mcos⁡(Ωt),ddot x+2gammadot x+omega_0^2x=F_0mcos(Omega t),

证明稳态下平均耗散功率 ⟨P⟩langle Prangle 取最大值时的驱动频率为 Ω=ω0Omega=omega_0.

解: 稳态响应 x=Acos⁡(Ωt−δ)x=Acos(Omega t-delta).
阻尼耗散功率 Pd=bx˙2P_d=bdot x^2,周期平均

⟨P⟩=b⟨x˙2⟩=b⋅12A2Ω2.langle Prangle=blangledot x^2rangle=b*frac12A^2Omega^2.

用 b=2mγb=2mgamma 与

A2=(F0/m)2(ω02−Ω2)2+(2γΩ)2A^2=(F_0/m)^2(omega_0^2-Omega^2)^2+(2gammaOmega)^2

得

⟨P⟩∝Ω2(ω02−Ω2)2+(2γΩ)2.langle Pranglepropto Omega^2(omega_0^2-Omega^2)^2+(2gammaOmega)^2.

对 ΩOmega 求极值可得最大值发生在 Ω=ω0Omega=omega_0.

阻尼振动

运动方程

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

d2xdt2+2βdxdt+ω2x=0d^2xdt^2 + 2beta dxdt + omega^2 x = 0

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

解的形式

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

  1. 欠阻尼振动(β<ωbeta < omega):
x(t)=Ae−βtcos⁡(ωdt+ϕ)x(t) = A e^-beta t cos(omega_d t + phi)

其中 ωd=ω2−β2omega_d = omega^2 - beta^2 为阻尼振动的角频率.

图例

欠阻尼振动

  1. 临界阻尼(β=ωbeta = omega):
x(t)=(C1+C2t)e−βtx(t) = (C_1 + C_2 t)e^-beta t
图例

临界阻尼振动

  1. 过阻尼(β>ωbeta > omega):
x(t)=C1er1t+C2er2tx(t) = C_1 e^r_1 t + C_2 e^r_2 t

其中 r1r_1 和 r2r_2 为两个负实根.

图例

过阻尼振动

能量衰减

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

E(t)=E0e−2βtE(t) = E_0 e^-2beta t

受迫振动

运动方程

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

d2xdt2+2βdxdt+ω2x=F0cos⁡(Ωt)d^2xdt^2 + 2beta dxdt + omega^2 x = F_0 cos(Omega t)

其中:

  • F0F_0 为外力的振幅;
  • ΩOmega 为外力的角频率.

稳态解

稳态解为:

x(t)=Acos⁡(Ωt−δ)x(t) = A cos(Omega t - delta)

其中:

  • A=F0(ω2−Ω2)2+(2βΩ)2A = F_0/(omega^2 - Omega^2)^2 + (2betaOmega)^2 为稳态振幅;
  • δ=arctan⁡(2βΩω2−Ω2)delta = arctanleft(2betaOmega/omega^2 - Omega^2right) 为相位差.

共振现象

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

共振条件下的振幅为:

Ares=F02βωA_res = F_02betaomega

应用

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

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


评论