微分方程 · 阻尼

阻尼

阻尼使振动损失能量。调整阻尼比ζ,比较物块来回振荡与不越过平衡位置的返回过程。时间以无阻尼周期T₀为单位。

01

阻尼与本页画面

阻尼是让振动衰减下去的能量损耗。带阻尼器的弹簧振子满足mx′′+cx′+kx=0m x'' + c x' + k x = 0,写成x′′+2ζω0x′+ω02x=0x'' + 2\zeta\omega_0 x' + \omega_0^2 x = 0,其中固有角频率ω0=k/m\omega_0 = \sqrt{k/m},阻尼比ζ=c/(2mk)\zeta = c/(2\sqrt{mk})。ζ小于1是欠阻尼(来回振荡),等于1是临界阻尼,大于1是过阻尼(慢慢爬回)。竖直挂着的物块受的重力只把平衡位置下移mg/k,x从这个位置量起。本页让物块从x = 1静止放开,问:阻尼越大,停得越快吗?

特征根决定运动;临界重根还产生t·e^(λt)项
λ2+2ζω0λ+ω02=0,λ=ω0(−ζ±ζ2−1),σ=−max⁡Re⁡λ\lambda^2 + 2\zeta\omega_0\lambda + \omega_0^2 = 0, \qquad \lambda = \omega_0\left(-\zeta \pm \sqrt{\zeta^2 - 1}\right), \qquad \sigma = -\max \operatorname{Re}\lambda

画面右边是特征根平面。两根之积是ω02\omega_0^2,所以欠阻尼时它们在半径ω₀的圆上:实部−ζω0-\zeta\omega_0决定衰减,虚部ωd=ω01−ζ2\omega_d = \omega_0\sqrt{1-\zeta^2}决定振动的角频率,所以阻尼也让每次来回变慢。ζ = 1时两根在−ω₀相遇;过了临界,一个根往左跑,另一个退回0附近,ζ较大时约为−ω0/(2ζ)=−k/c-\omega_0/(2\zeta) = -k/c。运动中最慢的那部分按e−σte^{-\sigma t}衰减,这个渐近衰减率σ就是最右边的根到虚轴的距离,恰好在ζ = 1时最大。x(t)图上的蓝色虚线画出了它:欠阻尼时是包络±e−ζω0t/1−ζ2\pm e^{-\zeta\omega_0 t}/\sqrt{1-\zeta^2}(曲线在两个极值之间与它相切),过阻尼时是慢根那一项。

“停下”要有一个标准。本页用T1%T_{1\%}:从这一刻起|x|始终不超过放开时的1%,这样的最早时刻。x(t)下方的窄图把这条带放大。极值在tn=nπ/ωdt_n = n\pi/\omega_d,∣x(tn)∣=e−ζω0tn|x(t_n)| = e^{-\zeta\omega_0 t_n},由此精确判断哪些峰伸出带外;最后一次回到带内的时刻再在闭式解上用二分法求出,所以这个读数是数值的。ζ稍小一点、多一个峰伸出带外时,稳定时间就跳一次,所以T1%T_{1\%}随ζ的曲线是锯齿形。它的最小值不在临界阻尼,而在ζ∗=ln⁡100/π2+ln⁡2100≈0.826\zeta^* = \ln 100/\sqrt{\pi^2 + \ln^2 100} \approx 0.826:这时唯一一次越过平衡位置恰好碰到−1%,用时0.667 T₀,临界阻尼要1.057 T₀。临界阻尼只在“不越过x = 0”的运动里最快,而且只对从静止放开而言;带宽不同,ζ*也不同(2%带约为0.780)。

02

值得注意的

  • 两种“最快”不是一回事:渐近衰减率σ恰好在ζ = 1最大;稳定时间在稍小一点的ζ*最短。字面上的0.826不是ζ*:它越过1.0015%,稳定时间跳到0.896 T₀。本页把ζ*存成精确值。
  • 阻尼比已经很大时再加倍,爬回来的时间也约加倍:慢根的时间常数约为c/k。
  • 三种情形的x(t)都是精确解;键入的阻尼比按写下的数精确判断情形,所以0.99999999999999999是欠阻尼。液缸里的液体只是示意。
03

操作

  • 推动ζ推子(在0.1、1、10处有卡位),或键入分数、小数,或按ζ = 1、ζ*。根可以沿根轨迹拖动;到−ω₀会被吸住,也只有在那里才能换到另一支。拖动红色的时间游标,物块跟着移动。
  • 演示模式打开全屏视图:← →按推子的一步改变ζ,H隐藏衰减率和稳定时间。
04

相关概念

在相图一页的迹–行列式地图上,阻尼振子是一条横线Δ = ω₀²。

延伸阅读:维基百科:阻尼;OpenStax:University Physics Volume 1, 15.5 Damped Oscillations(英文)。