ABAQUS/Standard vs. 非線性 (基礎篇)
公眾號:CAE分析技術共享
眾所周知,ABAQUS是一款功能非常強大的通用有限元軟件,一提到ABAQUS大家首先想到的肯定是其強大的非線性求解功能。那么本篇文章我們將一起研究一下ABAQUS/Standard是怎么求解非線性問題的。在說非線性之前,我們先簡單的過一遍線性問題,以便更好的理解非線性,以及非線性的求解。
01. 什么是線性問題
線性問題的定義:在外部載荷的作用下,結構的響應是線性的。比如假設一個線性彈簧在10N的載荷下,伸長1cm。那么當施加20N載荷時,它將會伸長2cm。
ABAQUS/Standard在線性分析中,結構的剛度矩陣只計算一次。通過將新的載荷向量乘以剛度矩陣的逆矩陣(又稱柔度矩陣),可得到結構對其他載荷的相應。
02 什么是非線性
非線性結構問題是指結構的剛度隨其變形而改變的問題。所有的物理結構均是非線性的。線性只是一種近似,很顯然對許多有限元模型比如加工過程 ,沖壓,碰撞等分析,線性分析是遠遠不夠的。
由于結構的剛度會時刻改變,所以我們不能像線性問題一樣,還是用直接法計算的結構剛度矩陣。在非線性分析中,結構的剛度矩陣必須進行許多次迭代生成計算,因此,這使得非線性的分析相比線性分析昂貴的多。
非線性的來源:
1. 材料非線性,比如材料的應力-應變關系成非線性
2. 邊界非線性,比如結構在變形的過程過中生接觸,使得整個結構的邊界條件發生變化
3. 幾何非線性,大轉動擾動,突然的翻轉,屈曲現象等。
非線性問題求解的實質就是求解靜態平衡方程,【P - I = 0】。外載荷P必須與內應力I互相平衡。
ABAQUS/Standard應用迭代算法獲得非線性問題的解答,通過把總的外部載荷,分成若干個小的載荷,增量的施加給結構,逐步的獲取最終的解答。
ABOUT STEPS,INCREMENTS,ITERATIONS!
這三個概念是Standard求解器非常重要的概念,弄清楚三者的含義,以及他們之間的關系非常重要。下面我們通過一個實際的小例子來看一下,假設圓柱受到一個剛性板約10000N的載荷向下壓。

分析步:這里我們可以將這個問題分為2個STEP來計算。
1. STEP1施加一個很小的載荷比如10N,使得上板與下板之間的接觸關系穩定。
2. STEP2經過第一步,接觸已穩定建立,這時就可以將剩余的載荷在STEP2中施加到圓柱上。
增量步:我們假設在STEP2中,我們設置了1000個增量步,也就意味著,剩余的載荷將會分成1000份來逐步施加到圓柱上。
迭代步:迭代步是在一個增量步中尋找平衡解答的一次嘗試。在迭代結束后,如果模型不處于平衡狀態,那么就會采用新一輪的迭代步進行嘗試。一個增量步中需要多次的迭代才能找到可以接受的解答。
04 平衡迭代和收斂判斷
ABAQUS/Standard求解器對于一個小的載荷增量△P,結構的非線性如下圖所示:
平衡迭代:
已知△P和結構在U0處的切向剛度K0, 計算出修正位移Ca。Ca = △P / K0
利用Ca,將結構更新到Ua處。Ua = U0 + Ca
基于結構的更新構形Ua,并形成新的剛度Ka,隨之計算出新的內部作用力Ia。 Ia = I(Ua)
計算迭代殘差力Ra;Ra = P – Ia
收斂判斷:
對迭代殘差力Ra的比較,如果Ra小于平均力的0.5%,那么接受此解。
對修正位移Ca的比較,如果Ca小于總的增量位移,那么接受此解。
當倆個判斷標準均滿足時,Solver接受此解,并寫入結果文件,進行下一個增量步的計算。
如果上述條件不滿足時,將會進行下一次迭代,如下圖所示。此時solver采用Ua處的切向剛度Ka來計算Cb,并更新Ib,獲得Rb后對Rb以及Cb進行收斂判斷

為了更清楚了解平衡迭代的過程,將其步驟整理成流程表,如下所示:

05 Static, general分析步
本小節,我們將進入ABAQUS/CAE界面。觀察下Static,general分析步在做非線性的時候有哪些設置,以及這些設置的含義。
[Basic 選項卡]
Description:描述,可以寫一些對分析步的描述信息,方便他人觀察
Time period: 默認1,并非實際物理意義上的1秒。意思就是這個分析步的整體為1。
NLgeom: 幾何非線性開關,做非線性的時候打開
Automatic stabiization: 自動穩定控制,不推薦打開。這個功能將會引入人工能來幫助求解器穩定求解。
[Increment 選項卡]
Type: Automatic, 自動增量控制選項,推薦。
Maximum number of increments: 最大增量步數,推薦設大一點,這個選項對計算速度精度沒有任何影響,只是控制增量步數,如果你的計算中增量步數超過設定值,軟件就會拋出錯誤。
Increment Size:增量步大小,intial size確定初始增量步大小,推薦設置的小一點,幫助求解器更容易獲得收斂的解。Min Size,增量步最小值,當增量步小于這個值,軟件將自動中止分析。Max Size,增量步最大值,控制增量步的最大值。

[Other 選項卡]
Method:Solver求解方法,直接法和迭代法
Solution Technique: 求解技術,分為Full
Newton和Quasi-Newton, 兩者之間的最大差距就是,Full Newton在每個迭代步中結構的剛度矩陣都會重新生成,而Quasi-Newton默認下,每8個增量步才會重新生成剛度矩陣。
Extrapolation: 外推。默認Linear,非線性分析步的第一個增量步之后的每一個增量上,ABAQUS/Standard通過從前面的增量外推解來評估當前增量的解。默認情況下,使用100%的線性外推。如果Ti ≤ DE * Ti-1,則放棄外推,DE默認為0.1。
【PS】如喜歡本文請掃一下二維碼,關注一下公眾號謝謝。
工程師必備
- 項目客服
- 培訓客服
- 平臺客服
TOP




















