ABAQUS/Standard vs. 非線性 (基礎篇)

  公眾號:CAE分析技術共享

  眾所周知,ABAQUS是一款功能非常強大的通用有限元軟件,一提到ABAQUS大家首先想到的肯定是其強大的非線性求解功能。那么本篇文章我們將一起研究一下ABAQUS/Standard是怎么求解非線性問題的。在說非線性之前,我們先簡單的過一遍線性問題,以便更好的理解非線性,以及非線性的求解。

01. 什么是線性問題

    線性問題的定義:在外部載荷的作用下,結構的響應是線性的。比如假設一個線性彈簧在10N的載荷下,伸長1cm。那么當施加20N載荷時,它將會伸長2cm。

    ABAQUS/Standard在線性分析中,結構的剛度矩陣只計算一次。通過將新的載荷向量乘以剛度矩陣的逆矩陣(又稱柔度矩陣),可得到結構對其他載荷的相應。

ABAQUS/Standard vs. 非線性 (基礎篇)的圖1

     

   

02 什么是非線性

    非線性結構問題是指結構的剛度隨其變形而改變的問題。所有的物理結構均是非線性的。線性只是一種近似,很顯然對許多有限元模型比如加工過程 ,沖壓,碰撞等分析,線性分析是遠遠不夠的。

    由于結構的剛度會時刻改變,所以我們不能像線性問題一樣,還是用直接法計算的結構剛度矩陣。在非線性分析中,結構的剛度矩陣必須進行許多次迭代生成計算,因此,這使得非線性的分析相比線性分析昂貴的多。

非線性的來源:

    1.    材料非線性,比如材料的應力-應變關系成非線性

    2.    邊界非線性,比如結構在變形的過程過中生接觸,使得整個結構的邊界條件發生變化

    3.    幾何非線性,大轉動擾動,突然的翻轉,屈曲現象等。

ABAQUS/Standard vs. 非線性 (基礎篇)的圖2



   
03 非線性問題的求解    

    非線性問題求解的實質就是求解靜態平衡方程,【P - I = 0】。外載荷P必須與內應力I互相平衡。

ABAQUS/Standard vs. 非線性 (基礎篇)的圖3

    

    ABAQUS/Standard應用迭代算法獲得非線性問題的解答,通過把總的外部載荷,分成若干個小的載荷,增量的施加給結構,逐步的獲取最終的解答。

ABOUT STEPS,INCREMENTS,ITERATIONS!

    這三個概念是Standard求解器非常重要的概念,弄清楚三者的含義,以及他們之間的關系非常重要。下面我們通過一個實際的小例子來看一下,假設圓柱受到一個剛性板約10000N的載荷向下壓。

ABAQUS/Standard vs. 非線性 (基礎篇)的圖4

ABAQUS/Standard vs. 非線性 (基礎篇)的圖5

    分析步:這里我們可以將這個問題分為2個STEP來計算。

            1.    STEP1施加一個很小的載荷比如10N,使得上板與下板之間的接觸關系穩定。

            2.    STEP2經過第一步,接觸已穩定建立,這時就可以將剩余的載荷在STEP2中施加到圓柱上。

    增量步:我們假設在STEP2中,我們設置了1000個增量步,也就意味著,剩余的載荷將會分成1000份來逐步施加到圓柱上。

    迭代步:迭代步是在一個增量步中尋找平衡解答的一次嘗試。在迭代結束后,如果模型不處于平衡狀態,那么就會采用新一輪的迭代步進行嘗試。一個增量步中需要多次的迭代才能找到可以接受的解答。



04 平衡迭代和收斂判斷

ABAQUS/Standard求解器對于一個小的載荷增量△P,結構的非線性如下圖所示:

ABAQUS/Standard vs. 非線性 (基礎篇)的圖6

平衡迭代:

    • 已知△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進行收斂判斷

ABAQUS/Standard vs. 非線性 (基礎篇)的圖7

ABAQUS/Standard vs. 非線性 (基礎篇)的圖8

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

ABAQUS/Standard vs. 非線性 (基礎篇)的圖9

ABAQUS/Standard vs. 非線性 (基礎篇)的圖10


 

05 Static, general分析步

    本小節,我們將進入ABAQUS/CAE界面。觀察下Static,general分析步在做非線性的時候有哪些設置,以及這些設置的含義。

ABAQUS/Standard vs. 非線性 (基礎篇)的圖11

[Basic 選項卡]

    • Description:描述,可以寫一些對分析步的描述信息,方便他人觀察

    • Time period: 默認1,并非實際物理意義上的1秒。意思就是這個分析步的整體為1。

    • NLgeom: 幾何非線性開關,做非線性的時候打開

    • Automatic stabiization: 自動穩定控制,不推薦打開。這個功能將會引入人工能來幫助求解器穩定求解。

      ABAQUS/Standard vs. 非線性 (基礎篇)的圖12

ABAQUS/Standard vs. 非線性 (基礎篇)的圖13[Increment 選項卡]

    • Type: Automatic, 自動增量控制選項,推薦。

    • Maximum number of increments: 最大增量步數,推薦設大一點,這個選項對計算速度精度沒有任何影響,只是控制增量步數,如果你的計算中增量步數超過設定值,軟件就會拋出錯誤。

    • Increment Size:增量步大小,intial size確定初始增量步大小,推薦設置的小一點,幫助求解器更容易獲得收斂的解。Min Size,增量步最小值,當增量步小于這個值,軟件將自動中止分析。Max Size,增量步最大值,控制增量步的最大值。

ABAQUS/Standard vs. 非線性 (基礎篇)的圖14ABAQUS/Standard vs. 非線性 (基礎篇)的圖15

[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】如喜歡本文請掃一下二維碼,關注一下公眾號謝謝。

ABAQUS/Standard vs. 非線性 (基礎篇)的圖16

登錄后免費查看全文
立即登錄
App下載
技術鄰APP
工程師必備
  • 項目客服
  • 培訓客服
  • 平臺客服

TOP

14
4
15