當地震來臨，結構其實在「跳舞」：從靜力分析走向動力分析

當地震來臨，結構其實在「跳舞」：從靜力分析走向動力分析

讀完入門篇，你已經知道一座橋為什麼算得出來：平衡方程式、靜定與靜不定、贅餘度。但那整套思維有一個隱藏的前提——載重是「慢慢」加上去的。樓板的恆載、人群的活載，都是準靜態（quasi-static）的，結構在受力後安靜地停在一個變形位置，我們解一組代數方程式就結束了。

可是地震不是這樣。當斷層錯動，地面在零點幾秒內來回劇烈搖晃，建築物的每一層樓都被甩來甩去。此時讓結構受力的，不是某個固定的外力，而是它自己的質量乘上加速度所產生的慣性力（inertia force）。結構不再靜止，而是在「跳舞」——它會以自己特定的頻率擺盪，會與地震波發生共振，會在某些樓層放大、在某些樓層抵銷。

要算出這支舞，靜力分析的代數方程式 $[K]\{D\} = \{P\}$ 已經不夠了。我們必須走進結構動力學（structural dynamics），把時間 $t$ 與質量 $m$ 加進方程式。這篇進階文章，我們就從「靜力如何長出動力」談起，一路走到台灣耐震設計真正在用的反應譜分析（response spectrum analysis）。

從一條彈簧開始：單自由度系統

最簡單的動力模型是單自由度系統（single degree of freedom, SDOF）。想像一根質量集中在頂端的柱子，或一座只有屋頂有重量的單層廠房。把它的質量記為 $m$、側向勁度記為 $k$、阻尼係數記為 $c$，當地面以加速度 $\ddot{u}_g(t)$ 搖晃時，結構相對於地面的位移 $u(t)$ 滿足這條運動方程式（equation of motion）：

$$m\ddot{u} + c\dot{u} + ku = -m\ddot{u}_g(t)$$

這條式子值得逐項細看。左邊三項分別是慣性力（$m\ddot{u}$）、阻尼力（$c\dot{u}$）、彈性恢復力（$ku$）；右邊 $-m\ddot{u}_g$ 是地震透過質量灌進來的等效外力——質量越大，地震對你的「拳頭」就越重，這也是為什麼耐震設計常常要設法減輕結構自重。

如果暫時不考慮阻尼與外力（自由振動），方程式簡化為 $m\ddot{u} + ku = 0$，它的解是一個正弦波，振盪的自然圓頻率（natural circular frequency）為：

$$\omega_n = \sqrt{\frac{k}{m}}, \qquad T_n = 2\pi\sqrt{\frac{m}{k}}$$

其中 $T_n$ 是自然週期（natural period），代表結構自由擺盪一個來回所需的秒數。這是整個耐震分析中最重要的單一數字。週期短（剛硬）的結構擺得快，週期長（柔軟）的高樓擺得慢。台灣常用一個粗估公式 $T_n \approx 0.1N$（$N$ 為樓層數），所以一棟 10 層樓的建築週期大約 1 秒。

阻尼：能量去哪了

真實結構不會永遠擺盪下去，因為有阻尼（damping）把能量逐漸耗散掉。我們通常用阻尼比（damping ratio）$\zeta$ 來描述：

$$\zeta = \frac{c}{2\sqrt{km}} = \frac{c}{2m\omega_n}$$

對一般的鋼筋混凝土（RC）建築，工程上常假設 $\zeta \approx 0.05$，也就是 5% 臨界阻尼。這個看似不起眼的小數字，其實隱含了結構在反覆載重下從微裂縫摩擦、材料遲滯（hysteresis）到接頭滑動的全部能量耗散機制。

阻尼的重要性在於：地震是一種強迫振動（forced vibration）。當地震波某個頻率成分恰好接近結構的自然頻率時會發生共振（resonance），位移被大幅放大。理論上若完全無阻尼，共振時的放大倍率會趨於無限大；正是 $\zeta$ 把這個放大倍率壓到有限值。共振峰的放大倍率約為 $\frac{1}{2\zeta}$，當 $\zeta = 0.05$ 時約為 10 倍——這就是為什麼軟弱地盤上的長週期地震波，會對週期相近的高樓造成毀滅性破壞（1985 年墨西哥城地震、1999 年台灣 921 地震中部分中高樓的災情，都與這種頻率匹配有關）。

多自由度系統與模態：結構的「舞步分解」

真實建築當然不只一個自由度。一棟 $n$ 層樓的建築，若把每層樓板的水平位移當成一個自由度，就成為多自由度系統（multi degree of freedom, MDOF）。此時運動方程式升級為矩陣形式：

$$[M]\{\ddot{u}\} + [C]\{\dot{u}\} + [K]\{u\} = -[M]\{\iota\}\ddot{u}_g(t)$$

其中 $[M]$ 是質量矩陣、$[K]$ 是入門篇提過的整體勁度矩陣、$[C]$ 是阻尼矩陣，$\{\iota\}$ 是影響向量（influence vector，標示哪些自由度感受到地面運動）。

這個 $n \times n$ 的耦合方程組看起來很可怕，但結構動力學有一個極優雅的武器：模態分析（modal analysis）。先解自由振動下的特徵值問題（eigenvalue problem）：

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

解出來會得到 $n$ 組特徵值 $\omega_i^2$（對應 $n$ 個自然頻率）與特徵向量 $\{\phi_i\}$（對應 $n$ 個振態 / 模態形狀，mode shape）。每一個振態就是結構的一種「基本舞步」：第一振態（first mode）通常是整棟樓像懸臂梁一樣同方向擺動的形狀，週期最長；高階振態則是樓層彼此反向、形狀更複雜、週期更短。

最關鍵的數學性質是模態正交性（modal orthogonality）：不同振態對質量與勁度矩陣都正交，

$$\{\phi_i\}^T [M] \{\phi_j\} = 0, \quad \{\phi_i\}^T [K] \{\phi_j\} = 0 \quad (i \neq j)$$

這個性質讓我們能把原本耦合的 $n$ 維方程組，解耦（decouple）成 $n 個彼此獨立的 SDOF 方程式。換句話說，一棟 20 層樓的複雜振動，可以被拆解成 20 個簡單彈簧各自擺盪後再疊加。這就是模態疊加法（mode superposition）的威力，也是所有商用軟體（ETABS、SAP2000、MIDAS）做動力分析的底層原理。

實務上還有個好消息：低階的少數幾個振態往往就涵蓋了大部分的有效質量參與（effective modal mass participation）。耐震規範通常要求納入的振態，其累積有效質量需達總質量的 90% 以上，因此一棟規則的建築可能只算前 3 到 6 個振態就足夠精確。

反應譜：把無數地震濃縮成一條曲線

工程師不可能對每棟樓都拿幾十條真實地震歷時去做逐步積分（雖然 time-history analysis 確實存在且用於重要結構）。我們需要一個更精煉的工具——反應譜（response spectrum）。

反應譜的構想極為巧妙：拿一條真實地震加速度紀錄，丟給一系列「自然週期不同、但阻尼比都是 5%」的 SDOF 系統，記錄每個系統在這場地震中經歷的最大反應（最大位移、最大加速度），然後把「最大反應 vs. 自然週期」畫成一條曲線。這條曲線就告訴你：對一個週期為 $T$ 的結構，這場地震最多會逼它承受多大的加速度。

設計用的反應譜不是單一場地震，而是把該地區許多可能地震、考慮地盤條件後包絡平滑化的結果。台灣《建築物耐震設計規範及解說》就規定了設計水平譜加速度係數 $S_{aD}$ 隨週期 $T$ 變化的形狀：短週期段為一段平台（加速度敏感區），長週期段則隨 $T$ 增大而下降（位移敏感區）。地震力的計算大致是：

$$V = \frac{S_{aD} \, I}{1.4 \alpha_y} \, W$$

其中 $V$ 是設計地震力（基底剪力 base shear）、$W$ 是建築總重、$I$ 是用途係數、$\alpha_y$ 是與結構系統韌性相關的折減係數。注意這裡再次出現入門篇強調的觀念：贅餘度與韌性高的靜不定系統，可以用較大的 $\alpha_y$ 把設計地震力折減下來——因為我們允許它在大地震下進入非線性、靠塑鉸耗能，而不要求它全程保持彈性。

看一個例子

讓我們用一個具體數字感受「週期決定地震力」這件事。

問題：某棟單層 RC 廠房，屋頂集中質量 $m = 1.0 \times 10^5\ \text{kg}$（即重量約 $981\ \text{kN}$），四根柱子合計提供側向勁度 $k = 4.0 \times 10^6\ \text{N/m}$。阻尼比 $\zeta = 0.05$。試求其自然週期，並估算它在設計反應譜下承受的基底剪力。

第一步：求自然頻率與週期。

$$\omega_n = \sqrt{\frac{k}{m}} = \sqrt{\frac{4.0 \times 10^6}{1.0 \times 10^5}} = \sqrt{40} \approx 6.32\ \text{rad/s}$$

$$T_n = \frac{2\pi}{\omega_n} = \frac{2\pi}{6.32} \approx 0.99\ \text{s}$$

這棟廠房的週期約 1 秒，落在一般設計反應譜由平台轉入下降段的附近。

第二步：從反應譜讀取譜加速度。

假設根據場址與地盤類別查得，在 $T \approx 1.0\ \text{s}$ 時設計譜加速度為 $S_a = 0.40g$（這裡採簡化值；實際需按規範微分區計算）。

第三步：估算基底剪力。

最簡化的估法是把譜加速度乘上質量（單振態、近似全部質量參與）：

$$V \approx m \cdot S_a = 1.0 \times 10^5 \times 0.40 \times 9.81 \approx 3.92 \times 10^5\ \text{N} = 392\ \text{kN}$$

也就是說，這場設計地震對廠房施加的等效水平力約 $392\ \text{kN}$，相當於建築自重（$981\ \text{kN}$）的 40%。

對比思考：如果我們把柱子加粗，讓勁度變成 4 倍（$k = 1.6 \times 10^7\ \text{N/m}$），週期會縮短一半到約 $0.49\ \text{s}$，落入反應譜的短週期平台段，此時 $S_a$ 可能反而上升到 $0.55g$ 左右——把結構做得更硬，地震力可能變更大！ 這個反直覺的結果正是靜力思維永遠看不到的：地震力不是固定外力，而是隨結構自身週期變動的「對手」。耐震設計因此往往不是一味追求剛硬，而是在勁度、週期、韌性之間做精巧的權衡。

重點回顧

1. 靜力長出動力的關鍵是兩個新角色：質量與時間。 運動方程式 $m\ddot{u} + c\dot{u} + ku = -m\ddot{u}_g$ 把慣性力與阻尼力加入入門篇的彈性恢復力。
2. 自然週期 $T_n = 2\pi\sqrt{m/k}$ 是耐震分析最核心的數字，決定結構與地震波會不會共振、以及在反應譜上讀到多大的地震力。
3. 阻尼比 $\zeta$（RC 建築常取 0.05）把共振放大倍率壓到有限，約為 $1/(2\zeta)$；頻率匹配的長週期地震波對高樓特別危險。
4. 模態分析靠特徵值問題與模態正交性，把複雜的 MDOF 振動解耦成多個獨立 SDOF，是所有結構分析軟體的動力學底層；低階少數振態通常即涵蓋 90% 以上有效質量。
5. 反應譜把無數地震濃縮成「最大反應 vs. 週期」的一條曲線；把結構做硬反而可能讓地震力變大，耐震設計是勁度、週期與韌性的權衡，而非單純求剛。

深入探討（研究所視角）

進入研究所，動力分析會從「假設線彈性、用反應譜」往三個更艱難也更真實的方向延伸。

非線性歷時分析與遲滯模型。 反應譜分析骨子裡仍是線彈性疊加（再乘上折減係數近似非線性）。但對隔減震建築、不規則結構或重要設施，規範允許甚至要求非線性歷時分析（nonlinear response history analysis）：選取多組符合場址設計譜的真實或人造地震歷時，對結構模型逐時步積分（常用 Newmark-$\beta$ 法）。此時每個塑鉸都要賦予遲滯模型（hysteresis model，如 Bilinear、Takeda、Clough 模型）來描述反覆載重下的勁度退化與能量耗散，而能量耗散的面積（遲滯迴圈所圍的面積）正是評估耐震性能的核心指標。模態正交性在非線性階段不再成立，疊加法失效，只能直接積分耦合方程組。

阻尼矩陣的建構難題。 $[M]$ 與 $[K]$ 都能由物理直接組裝，但 $[C]$ 其實沒有清晰的物理對應。實務上常用瑞利阻尼（Rayleigh damping）假設 $[C] = a_0[M] + a_1[K]$，再選兩個關鍵振態指定其阻尼比來反推 $a_0, a_1$。然而這個假設會讓其他振態的阻尼比偏離設計值，高階振態甚至可能被過度阻尼或激發數值假象，是計算結構動力學中持續被研究的議題。

主動／被動控制與性能設計。 現代耐震已從「讓結構硬扛」轉向「主動管理結構的動力特性」。被動控制如調諧質量阻尼器（tuned mass damper, TMD，台北 101 的著名巨型風阻尼球即屬此類）刻意調整一個附加質量的週期去抵銷主結構振動；隔震（base isolation）則用柔軟的隔震層大幅拉長週期、避開地震主要能量頻段。這些都直接操作本文的 $\omega_n$ 與 $\zeta$。而性能基礎耐震設計（performance-based seismic design, PBSD）更把分析目標從「不倒塌」細分為立即使用、生命安全、避免崩塌等多個性能水準，要求工程師對結構在不同地震回歸期下的反應做機率性評估。

值得延伸的思考：當你下次走進一棟高樓、感覺到它在強風中極輕微地晃動時，不妨想想——那不是結構出了問題，而是它正以自己的自然週期溫和地呼吸。真正危險的從來不是「會晃」，而是「晃的節奏剛好和地震對上」。理解這支看不見的舞，正是從靜力分析跨入動力分析、從「算得出來」進階到「撐得過去」的分水嶺。