摩擦力、空氣阻力與終端速度：耗散力的數學機制

從「等速墜落」談起：耗散力如何重塑運動

跳傘者離機後並不會無限加速；幾秒之內，速度便趨於一個穩定值。這個看似違反「自由落體」直覺的現象，正是耗散力（dissipative force）的標誌：當阻力隨速度增長到足以抵銷重力時，加速度歸零，運動進入終端速度（terminal velocity）。本文不停留在「摩擦會減速」的定性描述，而是直接從受力方程出發，推導終端速度、刻劃趨近的動力學，並釐清乾摩擦與流體阻力在數學結構上的根本差異。

對已具普物基礎的學習者而言，關鍵不在於「有沒有阻力」，而在於阻力作為速度的函數 $f(v)$ 具有何種形式——線性、平方、抑或與速度大小無關——這直接決定了運動方程是否可解、解的形式，以及系統是否存在吸引子（attractor）。

乾摩擦：庫倫定律的非解析本質

乾摩擦（Coulomb friction）的經驗律可寫為：靜摩擦上限 $f_s \le \mu_s N$，動摩擦 $f_k = \mu_k N$，其中 $N$ 為正向力、$\mu$ 為摩擦係數。值得強調的是，動摩擦的大小與滑動速度幾乎無關，但其方向恆與相對速度反向：

$$\vec{f}_k = -\mu_k N\,\frac{\vec{v}}{|\vec{v}|}$$

這個 $\vec{v}/|\vec{v}|$ 的因子使方程在 $\vec{v}=0$ 處不可微，這是乾摩擦數學上的「惡性」來源：它不是光滑的力場，而是一個帶有符號函數 $\operatorname{sgn}(v)$ 的非線性項。考慮一維滑塊在水平面減速：

$$m\frac{dv}{dt} = -\mu_k m g\,\operatorname{sgn}(v)$$

對 $v>0$，$dv/dt = -\mu_k g$ 為常數，故 $v(t)=v_0-\mu_k g\,t$，物體在有限時間 $t^\* = v_0/(\mu_k g)$ 內精確停止——這與流體阻力的指數漸近截然不同。乾摩擦使系統在有限時間到達靜止並「黏住」，這種有限時間奇異性（finite-time singularity）正是 stick-slip 振盪與摩擦音（如琴弓摩擦）的根源。

流體阻力：兩種極限與雷諾數

物體在流體中運動的阻力由雷諾數 $\mathrm{Re}=\rho v L/\eta$ 決定其機制（$\rho$ 為流體密度、$L$ 特徵尺度、$\eta$ 動黏度）。

低 $\mathrm{Re}$（黏滯主導，$\mathrm{Re}\lesssim 1$）：阻力線性正比於速度，史托克斯定律（Stokes' law）對球體給出 $f = 6\pi\eta r v$，即 $f = bv$。適用於霧滴、微生物、密立根油滴實驗。

高 $\mathrm{Re}$（慣性主導，$\mathrm{Re}\gtrsim 10^3$）：阻力來自前方流體動量轉移，正比於速度平方：

$$f = \tfrac{1}{2}C_d\,\rho A\,v^2 \equiv c v^2$$

其中 $C_d$ 為阻力係數、$A$ 為迎風截面積。跳傘者、雨滴、汽車皆屬此域。注意這兩種律並非矛盾，而是同一物理在不同 $\mathrm{Re}$ 的漸近極限。

終端速度與趨近動力學的推導

考慮鉛直下墜、平方阻力的情形，向下為正：

$$m\frac{dv}{dt} = mg - cv^2$$

終端速度由穩態 $dv/dt=0$ 立得：

$$v_T = \sqrt{\frac{mg}{c}} = \sqrt{\frac{2mg}{C_d\rho A}}$$

此式揭示為何重物落得快：$v_T \propto \sqrt{m/A}$，即正比於「面積荷重」的平方根。要解出完整時間演化，將方程改寫並分離變數。令 $k=c/m$，則 $\dot v = g - kv^2 = g(1-v^2/v_T^2)$。積分

$$\int \frac{dv}{1-(v/v_T)^2} = \int g\,dt$$

利用 $\int \frac{dx}{1-x^2}=\operatorname{artanh}(x)$，得

$$v(t) = v_T\tanh\!\left(\frac{g\,t}{v_T}\right)$$

雙曲正切的性質完美對應物理：$t\to 0$ 時 $\tanh x\approx x$，回到自由落體 $v\approx gt$；$t\to\infty$ 時 $\tanh\to 1$，速度漸近 $v_T$ 但永不真正達到——與乾摩擦的有限時間停止形成鮮明對比。特徵時間尺度為 $\tau = v_T/g$。

對線性阻力 $m\dot v = mg-bv$，解則為 $v(t)=v_T(1-e^{-t/\tau})$，$v_T=mg/b$，$\tau=m/b$，呈純指數弛豫。兩種阻力律都產生「吸引子」式的終端速度，但趨近的函數形式（$\tanh$ vs. 指數）不同，這在實驗上可用以辨別主導機制。

定量小範例：雨滴的終端速度

求一顆半徑 $r=1.0\,\text{mm}$ 雨滴的終端速度。取水密度 $\rho_w=1000\,\text{kg/m}^3$、空氣密度 $\rho=1.2\,\text{kg/m}^3$、球體阻力係數 $C_d\approx 0.47$、$g=9.8\,\text{m/s}^2$。

步驟一：質量 $m=\rho_w\cdot\tfrac{4}{3}\pi r^3 = 1000\times\tfrac{4}{3}\pi(10^{-3})^3 \approx 4.19\times10^{-6}\,\text{kg}$。

步驟二：迎風面積 $A=\pi r^2 = \pi(10^{-3})^2 \approx 3.14\times10^{-6}\,\text{m}^2$。

步驟三：代入

$$v_T=\sqrt{\frac{2mg}{C_d\rho A}}=\sqrt{\frac{2\times4.19\times10^{-6}\times9.8}{0.47\times1.2\times3.14\times10^{-6}}}$$

分子 $=8.21\times10^{-5}$，分母 $=1.77\times10^{-6}$，比值 $\approx 46.4$，故 $v_T\approx 6.8\,\text{m/s}$。

驗證 $\mathrm{Re}$：$\mathrm{Re}=\rho v_T (2r)/\eta = 1.2\times6.8\times2\times10^{-3}/(1.8\times10^{-5})\approx 900$，落在慣性主導域，確認採用平方律自洽——若誤用史托克斯線性律將得到嚴重偏高的速度。此即理論的自我一致性檢驗。

深入探討（研究所視角）

耗散力使系統脫離保守力學的核心：拉格朗日量 $L=T-V$ 與哈密頓量 $H=T+V$ 的標準表述要求力可由位能梯度導出，而摩擦與阻力無法寫成 $-\nabla V$。處理之道有二。其一是 Rayleigh 耗散函數（Rayleigh dissipation function）$\mathcal{F}=\tfrac{1}{2}\sum_{ij}c_{ij}\dot q_i\dot q_j$，修正後的尤拉–拉格朗日方程為

$$\frac{d}{dt}\frac{\partial L}{\partial \dot q_i}-\frac{\partial L}{\partial q_i}=-\frac{\partial \mathcal{F}}{\partial \dot q_i}$$

對線性阻力此式精確；$2\mathcal{F}$ 恰為系統耗散功率 $P=\vec{f}\cdot\vec{v}$。平方阻力則需推廣的非二次型耗散勢，數學上不再保證唯一性。其二是 Caldeira–Leggett 模型：將耗散視為系統與一個由無窮多諧振子構成之熱浴（bath）的耦合，對整體寫下保守哈密頓量，再積掉熱浴自由度。這在凝態與量子耗散理論中是處理摩擦與布朗運動的標準框架，並直接連結到 Langevin 方程 $m\dot v = -\gamma v + \xi(t)$ 與漲落–耗散定理（fluctuation–dissipation theorem），後者規定阻尼係數 $\gamma$ 與隨機力 $\xi$ 的關聯 $\langle\xi(t)\xi(t')\rangle=2\gamma k_B T\,\delta(t-t')$，將宏觀耗散與微觀熱漲落統一。

在連續介質層面，平方阻力的物理根源是 Navier–Stokes 方程在高雷諾數下的非線性對流項 $(\vec{u}\cdot\nabla)\vec{u}$，它導致邊界層分離、尾流渦旋與壓差阻力（form drag）。阻力係數 $C_d$ 隨 $\mathrm{Re}$ 的「阻力危機」（drag crisis）——球體在 $\mathrm{Re}\sim 3\times10^5$ 附近 $C_d$ 驟降——正源於邊界層由層流轉湍流、分離點後移，這是高爾夫球凹痕設計的物理基礎。

微觀摩擦學（nanotribology）則揭示乾摩擦的庫倫律並非基本定律，而是大量原子尺度黏滑事件的統計平均。Prandtl–Tomlinson 模型以單一原子在週期勢中受彈簧牽引，重現了與速度近乎無關的動摩擦與能量耗散的滯後迴圈，並預測在極低速度下出現對數速度依賴。從這個視角看，宏觀的 $\mu_k$ 是無數微觀能壘躍遷的湧現量，而真實接觸面積（由 Greenwood–Williamson 粗糙接觸理論刻劃）才是 $f_k\propto N$ 線性關係的微觀來源——亞莫頓定律（Amontons' law）的「正比於正向力、與表觀面積無關」由此得到統計力學的解釋。耗散力雖被普物課程當作經驗修正項，其完整理論卻橫跨流體力學、統計物理與量子多體，是連結宏觀現象與微觀機制的絕佳範例。