線性強化彈塑性umat子程序系列-增量迭代法

捕獲.png

一些非線性問題可歸納為一個如下的數學表達式

捕獲.png

對于力學問題,我們可以把P看作外載荷向量,q看作位移向量,Q(q)是關于q的非線性表達式。對于這樣的非線性問題,一般的有限元程序都是通過增量迭代法求解。增量迭代法的核心思想是,將最終的狀態看成是一個加載過程,將載荷分成多個增量,逐級加載,然后在每個增量步內多次迭代,收斂后進行下一個增量步。

1、增量法

將{P}荷載分成為m個荷載增量(相等或不等)

捕獲.png

,即總荷載為

捕獲.png

每次施加一個荷載增量,在第i步加載后,荷載為

捕獲.png

每一荷載增量產生一個位移增量

捕獲.png

線性強化彈塑性umat子程序系列-增量迭代法的圖7和應力增量

捕獲.png

在第i步加載后,位移、應力分別為

捕獲.png

第m步加載后,得到最終位移、應力。     

增量法的關鍵在于:已知前一個增量步的相關信息,如何由荷載增量

捕獲.png

線性強化彈塑性umat子程序系列-增量迭代法的圖11計算位移增量

捕獲.png

線性強化彈塑性umat子程序系列-增量迭代法的圖13和應力增量

捕獲.png

線性強化彈塑性umat子程序系列-增量迭代法的圖15,進而求出位移

捕獲.png

線性強化彈塑性umat子程序系列-增量迭代法的圖17和應力

捕獲.png

線性強化彈塑性umat子程序系列-增量迭代法的圖19的問題,這個問題通常應用牛頓-拉普森迭代法求解,接下來介紹這種方法。

         

捕獲.png

2、修正的牛頓-拉普森迭代法(mN-R)

在載荷

捕獲.png

線性強化彈塑性umat子程序系列-增量迭代法的圖22時,位移為

捕獲.png

線性強化彈塑性umat子程序系列-增量迭代法的圖24,下一個增量為

捕獲.png

線性強化彈塑性umat子程序系列-增量迭代法的圖26,下一個增量步結束后載荷為

捕獲.png

線性強化彈塑性umat子程序系列-增量迭代法的圖28,在已知以上條件后用mN-R方法計算下一個增量結束后的位移。

捕獲.png

計算切線剛度

捕獲.png

線性強化彈塑性umat子程序系列-增量迭代法的圖31

計算不平衡力

捕獲.png

根據非平衡力計算位移修正量

捕獲.png

位移修正量為

捕獲.png

一次迭代后修正的位移為

捕獲.png

判斷是否收斂,若不收斂繼續迭代,直到近似收斂于真實解

線性強化彈塑性umat子程序系列-增量迭代法的圖36

增量迭代法3.png

附件有個小算例,希望能幫助大家理解增量迭代法

同時也歡迎觀看本次的視頻教程

http://www.yqgqt.org.cn/college/video/c14014

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

TOP

5
1
2