阻尼、受迫振動與共振：從運動方程到品質因子

從自由振盪到受迫共振：一個微分方程的故事

把單擺拉開放手，它不會永遠擺下去；推鞦韆若推得「正是時候」，幅度卻能越來越大直到危險。這兩個日常現象的背後，是同一條二階線性常微分方程在不同邊界條件下的兩種面貌。本篇直接從運動方程出發，建立阻尼、受迫與共振的完整數學圖像。

考慮一個質量 $m$、彈簧常數 $k$ 的振子，額外受到正比於速度的黏滯阻尼力 $-b\dot{x}$ 與外加週期驅動力 $F(t)=F_0\cos\omega t$。牛頓第二定律給出

$$m\ddot{x} + b\dot{x} + kx = F_0\cos\omega t.$$

定義自然頻率 $\omega_0=\sqrt{k/m}$ 與阻尼比 $\zeta = b/(2\sqrt{mk})$（或阻尼率 $\gamma=b/2m=\zeta\omega_0$），可改寫為標準形式

$$\ddot{x} + 2\zeta\omega_0\,\dot{x} + \omega_0^2\,x = \frac{F_0}{m}\cos\omega t.$$

整個阻尼受迫振動的物理，都濃縮在這條方程的齊次解與特解之中。

阻尼自由振盪：三種命運

先關掉驅動力（$F_0=0$），求齊次解。代入試探解 $x\propto e^{\lambda t}$ 得特徵方程

$$\lambda^2 + 2\zeta\omega_0\lambda + \omega_0^2 = 0,\qquad \lambda_{\pm} = -\zeta\omega_0 \pm \omega_0\sqrt{\zeta^2-1}.$$

判別式的符號決定系統命運，這正是阻尼比 $\zeta$ 的物理意義所在：

• 欠阻尼（$\zeta<1$）：$\lambda_{\pm}=-\gamma\pm i\omega_d$，其中阻尼頻率 $\omega_d=\omega_0\sqrt{1-\zeta^2}$。解為

$$x(t)=A e^{-\gamma t}\cos(\omega_d t+\phi),$$

即在指數衰減的包絡內振盪，振幅每經一個週期衰減一個固定比例。

• 臨界阻尼（$\zeta=1$）：兩根重合於 $\lambda=-\omega_0$，解為 $x(t)=(A+Bt)e^{-\omega_0 t}$，是「不振盪卻最快回到平衡」的臨界情形，這也是車輛避震器與儀表指針的設計目標。

• 過阻尼（$\zeta>1$）：兩相異負實根，系統緩慢爬回平衡，比臨界阻尼還慢。

衡量欠阻尼系統「能量保存得多好」的無因次量是品質因子

$$Q = \frac{\omega_0}{2\gamma} = \frac{1}{2\zeta}.$$

$Q$ 物理上等於系統儲存能量與每弧度耗散能量之比的 $2\pi$ 倍除以 $2\pi$，亦即每振盪一弧度損失能量的倒數量度；$Q$ 越大代表阻尼越輕、振盪維持得越久。

受迫振動的穩態解

加回驅動力後，通解 = 齊次解（暫態，隨 $e^{-\gamma t}$ 消失）+ 特解（穩態）。我們關心長時間後的穩態。用複數法令 $\tilde{F}=F_0 e^{i\omega t}$、試探 $\tilde{x}=\tilde{X}e^{i\omega t}$，代入方程：

$$(-\omega^2 + 2i\zeta\omega_0\omega + \omega_0^2)\tilde{X} = \frac{F_0}{m}.$$

解出複振幅並取模，得到穩態振幅與相位落後角

$$X(\omega) = \frac{F_0/m}{\sqrt{(\omega_0^2-\omega^2)^2 + (2\zeta\omega_0\omega)^2}},\qquad \tan\delta = \frac{2\zeta\omega_0\omega}{\omega_0^2-\omega^2}.$$

這兩條曲線就是頻譜學、結構工程與電路分析的共同語言。相位 $\delta$ 從低頻的 $0$（位移與驅動同相）平滑過渡到高頻的 $\pi$（反相），而在 $\omega=\omega_0$ 時恰為 $\pi/2$。

共振：峰值在哪裡

「共振」常被籠統說成「驅動頻率等於自然頻率」，但嚴格而言，振幅 $X(\omega)$ 的峰值位置須令分母對 $\omega$ 微分為零。設 $g(\omega)=(\omega_0^2-\omega^2)^2+(2\zeta\omega_0\omega)^2$，

$$\frac{dg}{d\omega}=2(\omega_0^2-\omega^2)(-2\omega)+8\zeta^2\omega_0^2\omega=0.$$

解得振幅共振頻率

$$\omega_r = \omega_0\sqrt{1-2\zeta^2},$$

它略低於 $\omega_0$，且僅在 $\zeta<1/\sqrt{2}$ 時存在實數峰（否則振幅隨頻率單調遞減）。峰值振幅為

$$X_{\max}=\frac{F_0/m}{2\zeta\omega_0^2\sqrt{1-\zeta^2}}\approx \frac{F_0/k}{2\zeta}=Q\cdot\frac{F_0}{k}\quad(\zeta\ll1).$$

可見輕阻尼下，穩態振幅可放大到靜態位移 $F_0/k$ 的約 $Q$ 倍——這正是塔科馬窄橋與酒杯碎裂背後的放大機制。值得注意的是，能量（功率）共振峰值恰在 $\omega=\omega_0$，與振幅共振頻率 $\omega_r$ 並不完全重合；這是共振討論中常被混淆的細節。

驅動力對系統做功的平均功率為

$$\langle P\rangle = \tfrac{1}{2}F_0\,\omega X\sin\delta = \frac{F_0^2\,\gamma\,\omega^2/m}{(\omega_0^2-\omega^2)^2+(2\gamma\omega)^2},$$

其在 $\omega=\omega_0$ 達極大，半高全寬 $\Delta\omega\approx 2\gamma$，於是

$$Q\approx \frac{\omega_0}{\Delta\omega},$$

把品質因子與共振峰的「尖銳程度」直接聯繫起來。

定量小範例

設 $m=0.50\,\text{kg}$、$k=200\,\text{N/m}$、阻尼係數 $b=2.0\,\text{kg/s}$，驅動力幅 $F_0=5.0\,\text{N}$。

1. 自然頻率：$\omega_0=\sqrt{k/m}=\sqrt{200/0.50}=20\,\text{rad/s}$。
2. 阻尼比：$\zeta=b/(2\sqrt{mk})=2.0/(2\sqrt{0.50\times200})=2.0/(2\times10)=0.10$。屬欠阻尼。
3. 品質因子：$Q=1/(2\zeta)=5.0$。
4. 振幅共振頻率：$\omega_r=\omega_0\sqrt{1-2\zeta^2}=20\sqrt{1-0.02}\approx19.80\,\text{rad/s}$（與 $\omega_0$ 僅差約 1%）。
5. 靜態位移：$F_0/k=5.0/200=0.025\,\text{m}$。
6. 峰值振幅：$X_{\max}=\dfrac{F_0/k}{2\zeta\sqrt{1-\zeta^2}}=\dfrac{0.025}{2\times0.10\times\sqrt{0.99}}\approx0.126\,\text{m}$。

驗證放大倍率約為 $0.126/0.025\approx5.0$，與 $Q=5.0$ 一致，再次印證輕阻尼下峰值放大倍率 $\approx Q$。

深入探討（研究所視角）

阻尼受迫振子可在分析力學框架下重新理解。無阻尼受迫部分對應拉格朗日量 $L=\tfrac12 m\dot x^2-\tfrac12 kx^2+F(t)x$，由 Euler–Lagrange 方程 $\frac{d}{dt}\frac{\partial L}{\partial \dot x}-\frac{\partial L}{\partial x}=0$ 即得無阻尼運動方程。然而黏滯阻尼是非保守耗散力，無法直接寫入標準 $L$；通行做法是引入 Rayleigh 耗散函數 $\mathcal{F}=\tfrac12 b\dot x^2$，把 Euler–Lagrange 方程推廣為 $\frac{d}{dt}\frac{\partial L}{\partial \dot x}-\frac{\partial L}{\partial x}+\frac{\partial \mathcal{F}}{\partial \dot x}=0$，從而還原 $m\ddot x+b\dot x+kx=F(t)$。另一條進路是 Caldirola–Kanai 哈密頓量 $H=\frac{p^2}{2m}e^{-bt/m}+\tfrac12 k x^2 e^{bt/m}$，以顯含時間的方式把耗散編入正則結構，這在開放量子系統與耗散量子化的討論中具有重要地位。

頻域觀點則把穩態響應寫成轉移函數與輸入的乘積。定義 $\chi(\omega)=1/[m(\omega_0^2-\omega^2)+ib\omega]$，即線性響應理論中的廣義感受率（susceptibility），則 $\tilde X=\chi(\omega)\tilde F$。$\chi$ 的實部與虛部並非獨立，受 Kramers–Kronig 關係約束——這是因果性（響應不能早於激勵）在數學上的必然推論，從電介質極化率到光學折射率的色散關係皆由此而來。$\chi$ 的虛部正比於前述吸收功率，是漲落–耗散定理（fluctuation–dissipation theorem）的雛形：系統對外力的耗散性響應，與其平衡態的自發熱漲落由同一個 $\mathrm{Im}\,\chi$ 描述。

把單一振子推廣到連續介質，即得受迫阻尼波動方程，共振演化為駐波與本徵模態；耦合大量振子並送入量子化，則每個簡正模成為一個量子諧振子，受迫項對應相干態的位移算符，阻尼對應與熱庫耦合的 Lindblad 主方程。當回復力不再線性（如 Duffing 振子 $\ddot x+2\gamma\dot x+\omega_0^2 x+\beta x^3=F_0\cos\omega t$），共振峰會「折彎」產生遲滯與跳躍、出現次諧波與混沌，這是非線性動力學與當代 MEMS／光力學共振器研究的核心課題。從一條二階方程出發，竟能一路延伸到場論、開放量子系統與混沌，這正是阻尼受迫振動作為物理「最小模型」的迷人之處。