微分方程式 · 阻尼

阻尼

阻尼帶走振盪的能量。調整 ζ,比較振盪式的回復與不跨越平衡位置的回復。時間以無阻尼週期 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

相關概念

在相圖頁面的跡–行列式圖上,阻尼振子對應水平線 Δ = ω₀²。

延伸閱讀: 維基百科: Damping (英文); OpenStax, University Physics Volume 1, 15.5 Damped Oscillations (英文).