微分方程 · 共振

共振

共振是振子在接近固有频率的周期力驱动下出现的大幅响应。振幅一直增长,直到外力每个周期送进的能量等于阻尼耗散的能量。

共振与本页画面

共振是振子受到周期力推动、而周期力的频率接近它自身固有频率时出现的大幅振动。本页的物块挂在弹簧上,阻尼器让它减慢,外力F0cos⁡ωtF_0\cos\omega t推着它。页面逐个周期记下振幅,标出每一次推动什么时候在送进能量,再把结果和各个驱动频率下的稳态响应对照。

受迫阻尼振子
x′′+2ζω0x′+ω02x=F0mcos⁡ωtx'' + 2\zeta\omega_0 x' + \omega_0^2 x = \frac{F_0}{m}\cos\omega t

驱动频率写成r=ω/ω0r = \omega/\omega_0,位移以静态伸长xst=F0/kx_{\mathrm{st}} = F_0/k为单位,时间以固有周期T0T_0为单位。有阻尼时,从静止出发的运动是以驱动频率进行的稳态振动,加上一段逐渐衰减的自由运动。黄色的点标出每个驱动周期里∣x∣|x|的最大值。自由运动衰减时,它们趋向稳态振幅,途中也可能起伏。

记录下方的条是功率P=FvP = F v。力和物块运动方向相同时,推一下就送进能量。有阻尼、r=1r = 1的稳态里,物块正好比力晚四分之一周期,于是力和速度同步,每一次推动都在送进能量。rr偏离1的稳态里,每个周期都有一段力和运动反向,条上画成斜线。离r=1r = 1越远,这一段越接近半个周期,正负功几乎抵消。

振幅为什么不再长?设想一个振幅为AA、比力晚四分之一周期的正弦振动,r=1r = 1的稳态就是这样。每个周期外力送进πF0A\pi F_0 A,阻尼系数为cc的阻尼器耗散πcωA2\pi c\omega A^2。送进随AA增长,耗散随A2A^2增长,所以两者只在一个振幅上持平:A=xst/(2ζ)A = x_{\mathrm{st}}/(2\zeta)。振幅增长的过程中,能量账逐个周期列出实际的送进、耗散和储能的变化。

稳态振幅
Axst=1(1−r2)2+(2ζr)2\frac{A}{x_{\mathrm{st}}} = \frac{1}{\sqrt{(1 - r^2)^2 + (2\zeta r)^2}}

阻尼越小,峰越高越窄。ζ>0\zeta > 0时,平均输入功率不低于峰值一半的那段驱动频率,宽度正好是2ζ2\zeta,就是图上标出的半功率宽度。位移的峰低于r=1r = 1,在r=1−2ζ2r = \sqrt{1 - 2\zeta^2}处,前提是0<2ζ2<10 < 2\zeta^2 < 1。

阻尼很小时,自由运动按自己的频率ωd=ω01−ζ2\omega_d = \omega_0\sqrt{1 - \zeta^2}振荡。它和驱动频率接近但不相等时,两者时而同步、时而错开,记录就一涨一落形成拍。阻尼耗尽自由运动后,拍也就消失了。没有阻尼时,自由运动永远不会消失,运动也就不会稳定下来。r≠1r \ne 1时,曲线给出一个周期运动的振幅,但运动并不趋向它;r=1r = 1时,振幅会无限增长。

值得注意

  • 红箭头画的是力F/kF/k,用位移的刻度,所以最长时正好是xstx_{\mathrm{st}}。在共振的大振幅旁边,它显得很小。
  • 曲线和蓝点表示稳态。黄点是这次运动逐个周期测得的值,运动还在调整时,它们不一定落在曲线上。

操作

  • 改变rr或ζ\zeta不会让运动重来。物块从当前状态接着走,有阻尼时会趋向新的稳态。“从静止重来”才清空记录。
  • 在响应图或相位图上任意位置按下,驱动频率就设到那里。靠近r=1r = 1时会停在正好1。

相关概念