線性只是理想，非線性才是現實？當 P≠Ku

真實世界中幾乎不存在真正的線性力學行為。非線性分析（Nonlinear Analysis）是指結構的力-位移關係依賴於當前狀態（目前的位移、力、應力應變關係）的分析類型。

CAE 工程師最常遇到的挑戰，就是「非線性」分析，這也是在有限元素分析進階課程才深究的課題。如果你沒有做過非線性分析，可能難以想像，不知道也完全無法體會非線性分析的難處。它可能讓你忙個三天三夜，依舊得不到理想的結果？到底「非線性」是什麼，「非線性分析」又在解什麼呢？

嚴謹來說，線性行為只是非線性行為裡的一個特例。 在有限元素分析的理論基礎中，「線性問題其實是非線性問題的一個子集合」。換句話說，真實世界裡的力學行為本質上就是非線性的，只有在變形夠小、材料還未達降伏、邊界條件不會中途改變，這些條件都同時成立時，我們才能把它簡化成線性問題來計算。

這聽起來有點反直覺，畢竟大多數人學結構分析，都是先學線性，再學非線性。好像線性才是「正常情況」，非線性是「特殊狀況」。然而，我們從數學本質來看恰好相反，非線性分析是在求解一個「勁度矩陣會隨著受力或變形改變」的結構系統，而線性分析只是這個系統在特定假設下的簡化版本。 非線性分析中的勁度矩陣 K 值是會改變的，它會隨著材料狀態、幾何形狀、邊界條件的變化而即時改變，這也是為什麼非線性分析無法像線性分析一樣「一次計算就得到結果」，必須拆成一步一步的增量（increment）來逐步逼近答案。

這篇文章是「非線性分析」系列的入門，目標是先釐清核心觀念：為什麼我們說真實世界的力學行為，線性只是特例？非線性的『來源』又在哪裡？之後，我們再來深入討論材料、幾何、邊界條件三種非線性來源，這篇則是先了解整體概況。

線性分析的基本假設：P = Ku

在結構分析（無論是手算或有限元素法）中，最基礎的力-位移關係是虎克定律（Hooke’s Law）的線性矩陣方程式：

其中：

• P 是外加的力（或廣義力向量）
• u 是結構產生的位移（或廣義位移向量）
• K 是結構的勁度矩陣（stiffness matrix），代表結構「抵抗變形的能力」，和結構的幾何與材料有關。

在這個關係式是一個極為理想的假設，K 是一個固定不變的常數矩陣，完全不受 u 或 P 的大小影響。這代表什麼？代表你把外力加倍，位移就會剛好加倍；你把兩個載重案例的結果疊加起來，答案就是兩者的和（疊加原理)。變形量 u 與外力 P 成完美的線性正比關係。這是線性分析之所以計算快、結果好預測的根本原因。

然而，真實世界中幾乎不存在真正的線性力學行為。

當你拉伸一根金屬棒，原子間的鍵結力會隨距離改變；當你壓彎一張薄紙，幾何形狀的改變會瞬間讓結構失去抵抗力；當兩個零件碰撞在一起，傳遞力量的接觸面積隨時在變動。

「線性分析」是我們在變形極小、材料未降伏、邊界不變動的前提下，對物理世界所做的局部近似（Local Approximation）。這在很多「小變形、材料不降伏」的工程情境下是合理的簡化。一旦結構受力超出這個很小的理想範圍，剛度矩陣 K 就會開始動態變化，力學方程式必須是真實的非線性型態。

換個角度想：與其說「線性是常態，非線性是特例」，不如說線性分析是非線性分析在特定假設條件下的一個子集合。當這些假設（小變形、材料不降伏、邊界不變）同時成立時，非線性的方程式剛好退化成線性的形式；一旦有任何一個假設被打破，結構的真實行為就會顯露出它本來的非線性面貌。

當 K 不再是常數：P = K(P,u)·u

但真實世界的結構行為，往往不會這麼「簡單」。當結構產生的變形夠大、材料開始降伏、或是接觸邊界的狀態中途改變，勁度矩陣 K 就會跟著當下的受力或變形狀態改變，方程式因此變成：

注意這裡，K 從常數變成 K(P,u) —— K 隨著力量和位移變化，這是「非線性分析」與「線性分析」最核心的分水嶺。

也就是說，K 不再是一個固定的「抵抗變形能力」，而是會隨著結構目前的狀態（受力多大、變形多少）動態調整的量。 這也是為什麼非線性分析的結果無法用疊加原理計算——把兩個獨立的非線性結果加起來，並不等於同時施加兩個載重的結果，因為 K 在兩種情況下根本是不同的值。

關鍵定義：非線性分析（Nonlinear Analysis）是指結構的力-位移關係依賴於當前狀態（目前的位移、力、應力應變關係）的分析類型；線性分析則是這個關係與當前狀態無關的特例。

非線性的三大來源

既然 K 會變，接下來自然的問題是：是什麼原因讓 K 改變？ 工程上把非線性的來源歸納成三大類：

非線性來源	K 為什麼改變	常見情境
材料非線性（Material）	沒有任何材料能無限度地維持線性彈性。當應力超過降伏點（Yield Point），材料進入塑性流動，切線模數（Tangent Modulus）大幅下降，剛度矩陣 K 隨之衰減。	金屬沖壓成型、橡膠大變形、高溫下的材料蠕變（Creep）與黏彈性（Viscoelasticity）時間累積效應。
幾何非線性（Geometric）	當結構產生大位移（Large Displacement）、大轉動（Large Rotation）或大應變（Large Strain）時，改變後的形狀會反過來影響剛度。	吉他弦拉緊後剛度增加（應力剛化）、薄板受拉時厚度減薄（Shell Thickness Update）、薄壁柱體受壓時突然發生的挫屈（Buckling）。
邊界條件非線性（Boundary Conditions）	現實中沒有永遠固定的約束。當零件互相接觸、分離或滑移時，傳遞力量的自由度會瞬間改變，導致剛度矩陣產生不連續的跳躍。	齒輪咬合、衝擊碰撞、負載方向隨結構變形而旋轉的隨動壓力（Follower Forces）。

這三種非線性來源可能個別或同時發生在一個分析裡（例如金屬沖壓成形，材料降伏、大變形、模具接觸三者同時出現）。

為什麼非線性問題不能「一次解完」

線性分析因為 K 是常數，所以只要解一次矩陣方程式 u = K⁻¹P 就能得到答案。但非線性分析中，當我們給予特定的力量 P 較計算位移 u，這時候需要知道 K(P,u) ，而 K(P,u) 本身又依賴答案 u，這形成一個雞生蛋、蛋生雞的問題：你需要先知道 u 才能算出正確的 K，但你需要正確的 K 才能算出 u。

正因為如此，非線性分析必須：

1. 把總載重拆成很多小步（load increment），一步一步慢慢施加載重，而不是一次全部加上去
2. 在每一小步裡，用目前已知最好的 K 估計值去猜一個 u，再檢查猜得準不準，不準就修正，直到誤差小到可以接受為止

這個「猜測、檢查、修正」的過程，就是非線性分析求解方法的核心邏輯（最常見的方法是 Newton-Raphson Method）。不過這套方法本身有很多細節，我們在這裡先不聊這個。這篇只需要先建立一個概念：因為 K 會變，所以非線性分析注定是一個「逐步逼近」的過程，而不是一次到位的計算。

如何判斷需使用「非線性分析」？

回到模擬分析實務上，到底什麼時候用線性分析就足夠，什麼時候又必須是非線性分析呢？

既然真實世界都是非線性的，那就所有分析都以非線性來進行，是不是就可以了呢？千萬別這樣做。畢竟我們是工程師，該做什麼要了解原因，有邏輯的分析，然後也要懂得Work Smart。模擬分析的目的是要解決工程問題，不是樣樣都要求要非常精準或者增加自己的困擾。在不同的工程研發製造階段做模擬，有不同的目的，大家要斟酌自己決定。在建立模型前，必須很清楚自己的分析哪些特質，我們這裡就模擬分析的常見狀況，給大家一個方向。

可以透過以下 Checklist 快速評估：

• 位移量檢查：結構的最大位移是否已經超過整體特徵尺寸（或板厚）的 10%？（若是， 開啟幾何非線性）
• 應力水準檢查：結構內部的最大 von Mises 應力是否接近或超過材料的降伏強度？（若是，開啟材料非線性）
• 接觸狀態檢查：模型中是否存在多個零件？受力過程中零件之間是否會產生接觸、滑移或分離？（若是，設定邊界條件非線性）
• 時間與溫度效應：結構是否在高溫環境下長時間承受固定載荷？（ 若是，開啟 Creep / 黏彈性模型）
• 特殊材料檢查：模型是否包含橡膠、矽膠或密封件等大應變材料？（ 若是，材料非線性（Hyperelastic 超彈性模型）。
• 幾何失穩檢查：結構是否為薄壁/長細比大，且承受壓應力或彎矩？（ 若是，開啟幾何非線性與 Buckling 挫屈分析）。

結論

「線性」可以是工程師快速估算所採用的分析方法；「非線性」才是真實世界的力學行為。唯有真正瞭解非線性的來源，CAE 工程師在處理和設定非線性分析時，才能游刃有餘。

常見問題

想學習結構非線性分析嗎？歡迎參考 Marc 或 MSC Nastran 高階結構分析軟體。