規則說變就變：邊界條件非線性

邊界條件非線性是：結構的邊界條件本身——接觸狀態、支承勁度、載重方向，隨著分析過程改變，讓勁度矩陣 K 變為非線性。

我們在前面幾篇文章，了解到材料非線性透過 D 矩陣連續變化讓 K 改變，幾何非線性則是透過形狀連續變化（K1、K2）讓 K 改變；邊界非線性走的又是另一條路，而且這條路有一個前兩者都沒有的特徵：邊界條件的改變，往往不是平滑連續的，而是突然的跳變（例如：接觸從分離瞬間變成閉合、彈簧勁度隨變形量不連續切換）。這篇要談的邊界非線性有三種類型：接觸、非線性支承、非線性載重。

接觸：間隙的開與合

接觸問題在實際的物理系統裡很常見，金屬成形製程裡工件與模具的界面、齒輪轉動間的介面、汽車碰撞模擬，都是典型的接觸問題。

接觸問題有兩個關鍵特徵：間隙的開合，以及摩擦。間隙描述的是兩個物體之間接觸（間隙閉合）或分離（間隙打開）的狀態，這個狀態取決於物體的位移；摩擦則影響物體接觸之後彼此之間的作用關係，取決於接觸力的大小，以及接觸面的摩擦係數（庫倫摩擦）。因為間隙和摩擦的行為必須逐步累加才能追蹤，牽涉間隙與摩擦的分析必須以增量（increment）方式進行，每個增量步裡也經常需要疊代，才能讓間隙-摩擦的行為穩定下來。

這正是邊界非線性「不連續」特性最直接的表現：兩個物體的間隙狀態，只有「開」和「合」兩種，沒有中間過渡的模糊地帶，這一步還沒接觸、下一步突然就接觸了，勁度矩陣在這個瞬間會出現不連續的跳變，跟材料非線性裡 D 矩陣隨應力平滑變化的方式不同。

非線性支承：彈簧與基礎

除了接觸，結構的支承條件本身也可能是非線性的。常見的處理方式有兩種：非線性彈簧，以及彈性基礎。

線性和非線性的彈簧都可以在模型裡指定；在非線性問題中，彈簧的勁度、以及彈性基礎的等效彈簧勁度，可以透過使用者副程式（user subroutine）來修改 （這裡提到的是 Marc 軟體內的設定）。非線性彈簧選項裡，彈簧的增量力是：

其中 K 是彈簧勁度，Δu₂ 是彈簧第二端自由度的位移增量，Δu₁ 是彈簧第一端自由度的位移增量。彈簧的線性與非線性資料透過模型定義選項輸入；針對非線性彈簧，也可以依照先前已經發生的變形量，用使用者副程式指定 K 的數值。動態分析裡，彈簧選項也可以用來定義阻尼器；在熱傳分析或電學分析裡，彈簧選項則可以用來定義熱連結或電連結。

在無摩擦的非線性彈性基礎選項裡，結構的元素可以指定成支承在一個彈性基礎上，這個基礎以每單位面積的增量力支撐結構：

其中 K 是基礎的等效彈簧勁度（每單位表面面積），Δuₙ 是表面上某一點沿著 ΔPₙ 方向的增量位移。要輸入非線性基礎資料，需要用對應的模型定義選項；要依照基礎先前已經發生的變形量，指定非線性等效彈簧勁度 K 的數值，同樣可以用使用者副程式。

這兩條式子有一個共通點：K 本身是位移（或已經累積的變形量）的函數，這跟材料非線性裡 D 矩陣依賴應力狀態的邏輯很類似，只是這裡改變的不是材料內部的應力應變關係，而是結構跟外部支承之間的連接方式。

非線性載重：保守與非保守

第三種邊界非線性，發生在載重系統本身是非保守的情況下，跟隨力（follower force）和摩擦滑動（friction/sliding）效應都屬於這一類。

結構受力時，會不會形成非線性，關鍵在於這個力是不是隨著結構的位移而改變。這類力可以是保守的，例如離心力場，雖然離心力的方向和大小會隨結構旋轉、變形而改變，但它仍然是一個保守力場；也可以是非保守的，例如作用在懸臂梁上的跟隨力，這種力的方向會隨著結構的變形持續調整，跟隨結構末端的姿態改變。而且跟隨力也可能在局部是非保守的，但整合到整個結構上卻是一個保守的載重系統，受壓力的圓筒就是這樣的例子。

這裡有個容易忽略的重點：離心力雖然會隨位移改變、方向也會跟著旋轉，卻仍然屬於保守力；反而是懸臂梁上的跟隨力，看起來只是簡單地「方向跟著端點轉」，卻是非保守的。這代表「隨載重方向會不會變」不是判斷保守與否的依據，兩者的差異在於這個力做的功，經過一個完整的變形循環之後，能不能回到原本的能量狀態。

為什麼邊界非線性特別難處理

邊界非線性的本質，讓它們在數值上特別不容易處理：

不連續性是這類非線性天生的特性，這使得用增量線性近似的方式求解變得困難。力學裡最嚴重的一些非線性，正是由邊界條件的非線性所帶來的。也因此，充分意識到這些潛在的問題區域、對背後的原理有良好的理解，非常重要——這樣的認知與理解，能讓你有能力驗證數值計算出來的結果，並在初次嘗試失敗時採取其他的解決方式。

這段話直接點出邊界非線性和前兩種非線性來源（材料、幾何）最大的差異：材料非線性和幾何非線性，K 的變化雖然也很劇烈，但基本上是連續的函數關係；邊界非線性的 K，卻可能在某個瞬間直接跳變（間隙從開到合、彈簧勁度隨位移超過某個門檻突然切換），這種不連續性，正是它特別難處理的根源。

和模擬分析的關係?

接觸問題有專門的元素和選項處理。 Marc 提供兩種模擬接觸問題的方式：特殊的間隙-摩擦元素，以及 CONTACT 選項——這裡不深入這兩種方式各自的計算機制，那屬於求解方法的範疇，留給之後的系列討論。

非線性支承需要透過特定選項輸入。 線性與非線性彈簧資料，透過 SPRINGS 模型定義選項輸入；如果彈簧勁度要依照先前的變形量動態調整，需要用到 USPRNG 使用者副程式。非線性基礎資料則透過 FOUNDATION 模型定義選項輸入，同樣可以搭配 USPRNG 副程式指定隨變形量變化的等效彈簧勁度。

剛性連結的啟用與停用，可以用調適性線性約束處理。 在分析過程中，因為剛性連結被啟用或停用而產生的非線性，可以透過調適性線性約束（adaptive linear constraints）來建模——這代表邊界非線性不只是接觸、彈簧、載重這三種情況，結構內部的約束關係本身，也可能在分析過程中動態改變。

跟隨力需要搭配特定參數。 這部分在上一篇幾何非線性已經談過 FOLLOW FOR 參數的使用方式，這裡要補充的是：跟隨力屬於邊界非線性的原因，不是因為它牽涉到幾何形狀的改變，而是因為它本身是一種「非保守的載重系統」——這是邊界非線性的分類邏輯，跟幾何非線性著重的「K1、K2隨形狀改變」是不同的切入角度，兩篇文章雖然都提到跟隨力，但關注的重點不同。

常見問題

想學習結構非線性分析嗎？歡迎報名 CAE 訓練課程。