為什麼一座橋會被「對的風」吹垮？

為什麼一座橋會被「對的風」吹垮？

1940 年 11 月 7 日，美國華盛頓州的塔科馬海峽大橋（Tacoma Narrows Bridge）在僅約 64 公里/時的中等強度風中劇烈扭轉、最終斷裂墜海。風速並不算大，橋的靜力強度也綽綽有餘——它不是被「吹斷」的，而是被振動搖垮的。當外界擾動的節奏恰好踩中了結構自身的「脈搏」，能量便一波一波累積，振幅愈來愈大，直到材料無法承受。

這正是振動學（Vibration）的核心問題：任何具有質量與彈性的系統都有自己偏好的振動頻率，一旦外力以這個頻率持續推送，就可能發生共振（resonance），讓微小的力造成巨大的位移。 從手機馬達的嗡嗡聲、汽車過坑的彈跳、洗衣機脫水時的搖晃，到精密工具機的加工顫振（chatter）與航太結構的顫振（flutter），振動無所不在。理解它、預測它、控制它，是機械工程師的基本功。

從單自由度系統說起：質量、彈簧、阻尼

振動分析最基本的模型，是單自由度（Single Degree of Freedom, SDOF）彈簧—質量—阻尼系統：一個質量 $m$ 接在剛度為 $k$ 的彈簧上，並透過阻尼係數 $c$ 的黏滯阻尼器與地面相連。對位移 $x(t)$ 套用牛頓第二定律，得到振動學最重要的一條微分方程：

$$ m\ddot{x} + c\dot{x} + kx = F(t) $$

這裡 $\ddot{x}$ 是加速度、$\dot{x}$ 是速度。左邊三項分別代表慣性力（與加速度成正比）、阻尼力（與速度成正比、方向相反）、彈性恢復力（與位移成正比、把質量拉回平衡點）。右邊 $F(t)$ 是外加激振力。

若我們先把外力拿掉（$F(t)=0$），研究系統「被撥一下之後自己怎麼動」，就是自由振動（free vibration）。

無阻尼自由振動：固有頻率的誕生

把阻尼也暫時忽略（$c=0$），方程簡化為：

$$ m\ddot{x} + kx = 0 \quad\Longrightarrow\quad \ddot{x} + \frac{k}{m}x = 0 $$

這是標準的簡諧運動方程，其解為正弦振盪：

$$ x(t) = A\cos(\omega_n t) + B\sin(\omega_n t) $$

其中關鍵量是固有角頻率（natural angular frequency）：

$$ \omega_n = \sqrt{\frac{k}{m}} \quad (\text{rad/s}) $$

對應的固有頻率與週期為：

$$ f_n = \frac{\omega_n}{2\pi} = \frac{1}{2\pi}\sqrt{\frac{k}{m}} \quad (\text{Hz}), \qquad T_n = \frac{1}{f_n} = 2\pi\sqrt{\frac{m}{k}} $$

這個結果與你在普通物理學過的單擺、簡諧運動一脈相承。它告訴我們一件直覺卻深刻的事：系統愈硬（$k$ 大）振得愈快，質量愈重（$m$ 大）振得愈慢。 固有頻率只取決於系統本身的物理參數，與你一開始撥多大力無關——這正是「共振」之所以危險的根源：它是系統的內稟特性。

阻尼自由振動：能量如何被耗散

真實系統都有阻尼，振動會逐漸衰減。把阻尼放回去並定義兩個無因次量會讓分析變得優雅。阻尼比（damping ratio）：

$$ \zeta = \frac{c}{c_{cr}}, \qquad c_{cr} = 2\sqrt{km} = 2m\omega_n $$

其中 $c_{cr}$ 是臨界阻尼（critical damping）。原方程可改寫為標準型：

$$ \ddot{x} + 2\zeta\omega_n\dot{x} + \omega_n^2 x = 0 $$

依 $\zeta$ 的大小，系統行為分成三類：

• 欠阻尼（underdamped，$\zeta < 1$）：振盪式衰減，最常見。位移為 $$ x(t) = X_0 e^{-\zeta\omega_n t}\cos(\omega_d t - \phi) $$ 其中阻尼振盪頻率 $\omega_d = \omega_n\sqrt{1-\zeta^2}$ 略低於 $\omega_n$。包絡線 $e^{-\zeta\omega_n t}$ 控制衰減快慢。
• 臨界阻尼（critically damped，$\zeta = 1$）：最快回到平衡且不過衝，是汽車避震器與門弓器追求的理想狀態。
• 過阻尼（overdamped，$\zeta > 1$）：緩慢爬回平衡、不振盪，像泡在糖漿裡。

工程上常用對數衰減率（logarithmic decrement）從實驗波形反推阻尼比。量測相鄰兩個波峰振幅 $x_1, x_2$：

$$ \delta = \ln\frac{x_1}{x_2} = \frac{2\pi\zeta}{\sqrt{1-\zeta^2}} $$

小阻尼時 $\delta \approx 2\pi\zeta$，這是實驗室裡識別系統阻尼最方便的手法之一。

強迫振動與共振：振動工程的主戰場

當外力是簡諧形式 $F(t) = F_0\cos(\omega t)$，我們關心穩態響應（steady-state response）的振幅。定義頻率比 $r = \omega/\omega_n$（激振頻率除以固有頻率），穩態振幅為：

$$ X = \frac{F_0/k}{\sqrt{(1-r^2)^2 + (2\zeta r)^2}} $$

注意 $F_0/k$ 就是把力靜態加上去時的位移（靜變形 $X_{st}$）。兩者相除得到無因次的放大係數（Magnification Factor, MF）：

$$ \text{MF} = \frac{X}{X_{st}} = \frac{1}{\sqrt{(1-r^2)^2 + (2\zeta r)^2}} $$

這條曲線是整個振動學的靈魂：

• 當 $r \to 0$（激振很慢），MF $\to 1$，系統幾乎像靜力一樣跟著外力走。
• 當 $r \to \infty$（激振很快），MF $\to 0$，質量因慣性「跟不上」而幾乎不動。
• 當 $r \approx 1$（激振頻率接近固有頻率），分母第一項 $(1-r^2)^2$ 趨近於零，MF 急遽放大，這就是共振。阻尼愈小，共振峰愈尖、放大愈劇烈。在無阻尼極限下振幅理論上趨於無窮，這正是塔科馬大橋的物理本質。

精確的峰值出現在 $r_{peak} = \sqrt{1-2\zeta^2}$，峰值放大為 $\text{MF}_{max} = \dfrac{1}{2\zeta\sqrt{1-\zeta^2}}$。低阻尼時近似 $\text{MF}_{max}\approx \dfrac{1}{2\zeta}$——例如 $\zeta = 0.02$（2% 阻尼）時，共振可把外力放大 25 倍！這也說明為什麼阻尼在工程上如此珍貴。

激振力與位移之間還有相位差 $\phi = \arctan\dfrac{2\zeta r}{1-r^2}$：低頻時位移與力同相，共振時恰好落後 $90°$，高頻時落後接近 $180°$（質量與力反向）。這個相位翻轉是辨識共振是否發生的重要指標。

振動控制：三條主要路線

理解了共振，工程師的任務就是避開它或抑制它。三種策略：

1. 調離共振區（避頻）：透過改變 $m$ 或 $k$ 移動固有頻率，使工作頻率遠離 $r=1$。一般要求 $r > \sqrt{2}$（此時 MF < 1，傳遞力小於施力，即進入隔振有效區），這是隔振設計的黃金門檻。

2. 隔振（vibration isolation）：在振源與基礎之間放置軟彈簧與阻尼（如機器底座的避震墊）。傳遞率（transmissibility） 描述傳到基礎的力佔施力的比例：

$$ T_R = \frac{\sqrt{1+(2\zeta r)^2}}{\sqrt{(1-r^2)^2 + (2\zeta r)^2}} $$

只有當 $r > \sqrt{2}$ 時 $T_R < 1$ 才真正達到隔振效果。有趣的是，在隔振區增加阻尼反而會略微提高傳遞率——這是設計者必須權衡的取捨：阻尼幫你度過啟動時掃過共振的瞬間，卻在高頻區拖後腿。

3. 吸振（dynamic vibration absorber, DVA）：在主系統上加裝一個調諧質量—彈簧子系統（tuned mass damper, TMD），讓子系統在特定頻率「替主系統振動」、把能量吸走。台北 101 大樓頂端那顆 660 噸的金球，就是用來抵抗強風與地震引發的搖晃的調諧質量阻尼器。

看一個例子

一台 80 kg 的工業風扇安裝在四個並聯彈簧上，每個彈簧剛度 $k_1 = 1.0\times10^5\ \text{N/m}$，系統阻尼比 $\zeta = 0.05$。風扇轉速為 1200 rpm，因葉片不平衡產生簡諧激振，靜變形響應 $X_{st} = 0.5\ \text{mm}$。試評估其振動狀態。

步驟一：等效剛度與固有頻率。 四個彈簧並聯： $$ k = 4k_1 = 4.0\times10^5\ \text{N/m} $$ $$ \omega_n = \sqrt{\frac{k}{m}} = \sqrt{\frac{4.0\times10^5}{80}} = \sqrt{5000} \approx 70.7\ \text{rad/s} $$ $$ f_n = \frac{\omega_n}{2\pi} \approx 11.25\ \text{Hz} $$

步驟二：激振頻率與頻率比。 1200 rpm 換算： $$ f = \frac{1200}{60} = 20\ \text{Hz}, \qquad \omega = 2\pi f \approx 125.7\ \text{rad/s} $$ $$ r = \frac{\omega}{\omega_n} = \frac{125.7}{70.7} \approx 1.78 $$

步驟三：放大係數與穩態振幅。 $$ \text{MF} = \frac{1}{\sqrt{(1-1.78^2)^2 + (2\times0.05\times1.78)^2}} = \frac{1}{\sqrt{(-2.17)^2 + (0.178)^2}} \approx \frac{1}{2.18} \approx 0.46 $$ $$ X = \text{MF}\times X_{st} \approx 0.46 \times 0.5\ \text{mm} \approx 0.23\ \text{mm} $$

步驟四：判讀。 因為 $r \approx 1.78 > \sqrt{2} \approx 1.41$，系統運轉在隔振有效區，振幅被壓到靜變形之下，設計合理。但要警惕：風扇啟動或停機時轉速會由 0 緩慢掃過 $f_n \approx 11.25\ \text{Hz}$，那一瞬間 $r=1$、MF 可達 $\approx 1/(2\zeta)=10$，振幅暴增至約 5 mm。因此啟停過程必須夠快、且阻尼不可太低，避免在共振點停留太久。

重點回顧

• 任何「質量＋彈性」系統都有固有頻率 $\omega_n = \sqrt{k/m}$，它只取決於系統本身，與激振大小無關；當外力頻率接近 $\omega_n$ 就會共振。
• 阻尼比 $\zeta = c/(2\sqrt{km})$ 決定自由振動是欠阻尼（振盪衰減）、臨界阻尼（最快不過衝）還是過阻尼（緩慢爬回）。
• 強迫振動的放大係數 $\text{MF}=1/\sqrt{(1-r^2)^2+(2\zeta r)^2}$，共振峰約為 $1/(2\zeta)$，低阻尼時放大極為驚人。
• 隔振的黃金門檻是 $r>\sqrt{2}$，此時傳遞率才小於 1；高頻隔振區增加阻尼反而不利，需權衡。
• 振動控制三策略：避頻（改 $m,k$）、隔振（軟彈簧＋阻尼）、吸振（調諧質量阻尼器，如台北 101 的金球）。

深入探討（研究所視角）

單自由度模型雖然優雅，真實結構卻是多自由度（Multi-DOF）甚至連續體系統。對 $n$ 個自由度，運動方程寫成矩陣形式：

$$ [M]\{\ddot{x}\} + [C]\{\dot{x}\} + [K]\{x\} = \{F(t)\} $$

其中 $[M]、[C]、[K]$ 分別是質量、阻尼、剛度矩陣。無阻尼自由振動的解導向特徵值問題：

$$ \left([K] - \omega^2[M]\right)\{\phi\} = \{0\} $$

求解可得 $n$ 個固有頻率 $\omega_i$（特徵值）與對應的模態振型（mode shapes） $\{\phi_i\}$（特徵向量）。每個模態都是一個獨立的「振動花樣」——第一模態通常是整體最柔軟的擺動，高階模態則出現節點（node）與更複雜的變形。

關鍵技巧是模態分析（modal analysis）：利用振型的正交性 $\{\phi_i\}^T[M]\{\phi_j\}=0\ (i\neq j)$，把耦合的矩陣方程解耦成 $n$ 個獨立的 SDOF 方程，每個都用我們前面學的方法求解，再疊加回去。這讓上千自由度的有限元素（FEM）模型得以高效分析，也是實驗模態分析（EMA，用敲擊鎚與加速度計識別結構動態特性）的理論基礎。

對連續體（如樑、板、軸），系統有無限多個自由度，振動由偏微分方程描述。以等截面歐拉—伯努利樑（Euler–Bernoulli beam）的橫向自由振動為例：

$$ EI\frac{\partial^4 w}{\partial x^4} + \rho A\frac{\partial^2 w}{\partial t^2} = 0 $$

其中 $EI$ 為抗彎剛度、$\rho A$ 為單位長度質量。分離變數後可解出與邊界條件相關的固有頻率序列。例如兩端簡支樑：

$$ \omega_i = (i\pi)^2\sqrt{\frac{EI}{\rho A L^4}}, \quad i = 1,2,3,\dots $$

最後值得一提的是自激振動（self-excited vibration）——塔科馬大橋的真正元兇並非單純共振，而是氣流與結構耦合產生的顫振（flutter）與渦旋脫落（vortex shedding）：結構的運動本身會改變氣動力，形成正回饋，把能量持續注入系統，即使外界沒有週期性激振也會發散。這類流固耦合（fluid-structure interaction）問題進入非線性動力學與穩定性理論的範疇，是航太、風工程與旋轉機械（如工具機顫振、轉子動力學）的研究前沿，也提醒我們：線性振動理論是強大的起點，但絕非世界的全部。