首先給大家介紹一下非線性行為的含義，載荷引起結(jié)構(gòu)剛度發(fā)生變化，其中典型情況包括：金屬材料的彈塑性，大變形問題（細(xì)長桿+長薄板），狀態(tài)變化（接觸+單元生死）。一般非線性分為：幾何非線性，材料非線性和狀態(tài)非線性。

ANSYS程序的方程求解器計(jì)算一系列的聯(lián)立線性方程來預(yù)測工程系統(tǒng)的響應(yīng)。然而，非線性結(jié)構(gòu)的行為不能直接用這樣一系列的線性方程表示。需要一系列的帶校正的線性近似來求解非線性問題。

一、純粹的增量式求解

一種近似的非線性求解是將載荷分成一系列的載荷增量?？梢栽趲讉€(gè)載荷步內(nèi)或者在一個(gè)載荷步的幾個(gè)子步內(nèi)施加載荷增量。在每一個(gè)增量的求解完成后，繼續(xù)進(jìn)行下一個(gè)載荷增量之前程序調(diào)整剛度矩陣以反映結(jié)構(gòu)剛度的非線性變化。

但是，純粹的增量近似不可避免地要隨著每一個(gè)載荷增量積累誤差，導(dǎo)種結(jié)果最終失去平衡。如圖1所示

二、全牛頓-拉普森迭代

ANSYS程序通過使用牛頓－拉普森平衡迭代克服了這種困難，它迫使在每一個(gè)載荷增量的末端解達(dá)到平衡收斂（在某個(gè)容限范圍內(nèi)）。

下圖圖2描述了在單自由度非線性分析中牛頓－拉普森平衡迭代的使用。在每次求解前，NR方法估算出殘差矢量，這個(gè)矢量是回復(fù)力（對應(yīng)于單元應(yīng)力的載荷）和所加載荷的差值，然后使用非平衡載荷進(jìn)行線性求解，且核查收斂性。如果不滿足收斂準(zhǔn)則，重新估算非平衡載荷，修改剛度矩陣，獲得新解。持續(xù)這種迭代過程直到問題收斂。

三、弧長法

對某些物理意義上不穩(wěn)定系統(tǒng)的非線性靜態(tài)分析,如果你僅僅使用NR方法，正切剛度矩陣可能變?yōu)榻抵榷剃嚕瑢?dǎo)致嚴(yán)重的收斂問題。這樣的情況包括獨(dú)立實(shí)體從固定表面分離的靜態(tài)接觸分析，結(jié)構(gòu)或者完全崩潰或者“突然變成”另一個(gè)穩(wěn)定形狀的非線性彎曲問題。對這樣的情況，你可以激活另外一種迭代方法：弧長方法，來幫助穩(wěn)定求解?；￠L方法導(dǎo)致NR平衡迭代沿一段弧收斂，從而即使當(dāng)正切剛度矩陣的傾斜為零或負(fù)值時(shí)，也往往阻止發(fā)散。這種迭代方法以圖形表示如圖3所示。

ANSYS程序還提供了一系列命令來增強(qiáng)問題的收斂性，如線性搜索、自動載荷步長及二分等，可被激活來加強(qiáng)問題的收斂性，如果不能得到收斂，那么程序或者繼續(xù)計(jì)算下一個(gè)載荷步或者終止（依據(jù)你的指示）。

線性搜索選項(xiàng)

ANSYS的自動求解控制，將根據(jù)需要關(guān)閉或打開線性搜索。對大多數(shù)接觸問題，線性搜索打開。這個(gè)收斂增強(qiáng)工具用程序計(jì)算出的比例因子(具有0和1之間的值)乘以計(jì)算出的位移增量。如圖4所示。

牛頓－拉普森選項(xiàng)

在存在非線性時(shí)，ANSYS的自動求解控制將應(yīng)用完全牛頓－拉普森選項(xiàng)。

完全牛頓－拉普森法：程序使用完全的牛頓－拉普森方法。在這種處理方法中，每進(jìn)行一次平衡迭代，就修改剛度矩陣一次。?

修正的牛頓－拉普森法：使用修正的牛頓－拉普森方法。

不對稱矩陣完全牛頓－拉普森方法：應(yīng)用完全牛頓－拉普森方法，剛度矩陣在每一次平衡迭代中都修正。