abaqus真實應變的案例
一文搞清ABAQUS中真實應力和真實應變
應變計算公式為：
②仿真中的真實應力與真實應變
在使用ABAQUS仿真時，如果我們的材料屬于塑性材料范疇，分析時涉及較大變形，在分析時必須將其應力和應變定義成真實應力和真實應變，我們就需要將實驗數據中得到的名義數據轉化成真實數據。
總結來說，如果不使用真實應力與真實應變，我們就不能分析出我們想要的較為準確結果，這里考慮的是材料的非線性問題。
工程應力應變和真實應力應變 附常用材料應力應變數據下載
通常處理方法是：實驗采集的數據轉換成工程應力應變數據①，再通過上述公式轉換成真實的應力應變曲線②，通過真實應變減去彈性應變，得到最終的塑性應變。
實驗數據處理方法：將計算好的工程應變應力分別輸入EXCEL表格中，插入計算公式：Ln（1+A2）即可計算出真實應變，代入公式：B2*(1+A2)并下拉即可得到真實應力，假定第三行為最大彈性應變，真實應變減去彈性應變得到有效塑性應變。
有效塑性應變真實應力曲線即是我們處理好的可以導入有限元軟件的材料模型數據。
下載地址：常用材料應力應變數據
展開 名義應力應變與真實應力應變
名義應力應變與真實應力應變
在進行結構或者構件分析時，材料屬性往往是最為重要的。我們在材料試驗測試時，一般測出的試驗曲線是名義應力—應變曲線，即所謂的工程應力和工程應變之間的關系。還有一張應力應變，我們稱為真實應力應變，那么工程應力和真實應力有什么區別？
首先請看這張圖:
這里面的Stress 和 Strain 就是指的工程應力和工程應變。只是由于在進行應變計算時，并未考慮測試構件的長度伸長或者截面縮小，這相當于沒有考慮非線性的影響。
但其實我們可以看到，在斷口處A（這個面積才代表真正的受應力面）是非常小的，因而材料的真實強度時上升了的（是指單位體積或者單位面積上的，不是結構上的）。
真實應力的定義
考慮到上述情況，真實應力被定義了出來：
在有限變形中，只有Δl→0時，拉伸與壓縮的應變才是相同的，即：
以及：
其中：l為當前長度；l0為初始長度；ε為真實應變或對數應變。
與真實應變相對應的是真實應力，定義為：
其中：F是施加在材料上的力；A是當前面積。
如果給出真實應力和真實應變的曲線，那么在拉伸和壓縮下，承受有限變形的金屬有相同的應力—應變關系。
應力應變的轉換
在一些有限元軟件中，必須輸入真實應力—應變關系，MARC和ABAQUS都是這樣的有限元軟件，尤其是在定義塑性數據時。
展開 名義應力應變與真實應力應變
我們在材料試驗測試時，一般測出的試驗曲線是名義應力—應變曲線，即所謂的工程應力和工程應變之間的關系。還有一張應力應變，我們稱為真實應力應變，那么工程應力和真實應力有什么區別？
首先請看這張圖:
這里面的Stress 和 Strain 就是指的工程應力和工程應變。只是由于在進行應變計算時，并未考慮測試構件的長度伸長或者截面縮小，這相當于沒有考慮非線性的影響。
但其實我們可以看到，在斷口處A（這個面積才代表真正的受應力面）是非常小的，因而材料的真實強度時上升了的（是指單位體積或者單位面積上的，不是結構上的）。
真實應力的定義
考慮到上述情況，真實應力被定義了出來：
在有限變形中，只有Δl→0時，拉伸與壓縮的應變才是相同的，即：
以及：
其中：l為當前長度；l0為初始長度；ε為真實應變或對數應變。
與真實應變相對應的是真實應力，定義為：
其中：F是施加在材料上的力；A是當前面積。
如果給出真實應力和真實應變的曲線，那么在拉伸和壓縮下，承受有限變形的金屬有相同的應力—應變關系。
應力應變的轉換
在一些有限元軟件中，必須輸入真實應力—應變關系，MARC和ABAQUS都是這樣的有限元軟件，尤其是在定義塑性數據時。這時需要對試驗給出的材料數據進行轉換。
展開 
由工程應力、應變計算真實應力、應變
設t時刻的作用于試樣標距段加載方向的載荷為F，工程應力為，工程應變為，真實應力為，真實應變為。規定試樣受拉伸長時，載荷、應力及應變取正值，反之，受壓縮短時，載荷、應力及應變取負值，則
式1
式2
式中表示0時刻至t時刻試樣長度的增量。
式3
式4
真實應力 的求解變換中利用了材料變形過程中體積不變的假設，即。
繪制應力-應變曲線時，往往不管拉伸或者壓縮，都將應力和應變繪制成正值。這樣，拉伸時，按照公式<1>至<4>計算出的應力和應變均為正值，不需要進行變換；壓縮時，按照上述方法計算出的應力和應度均為負值，需要進行變換。變換方式為：對上述公式中所有的應力和應變乘-1。按照這種規則，壓縮時應力、應變用<5>至公式<8>進行汁算。其中公式<5>中F取負值。
式5
式6
式7
式8
?
?
?
展開 由工程應力、應變計算真實應力、應變
材料的名義應力、應變與真實應力、應變轉換公式的推導
名義應力、應變也叫工程(Engineering)應力、應變。
CAE軟件需要采用基于變形后的應力、應變，即真實的應力、應變。
其中
為試件當前截面面積，
為試件當前長度。
LS-DYNA 材料的真實應力應變曲線是怎么回事？ ￥3.99
LS-DYNA 材料的真實應力應變曲線是怎么回事？
關于《LS-DYNA 材料的真實應力應變曲線是怎么回事？》的補充-Python后處理腳本 ￥1.99
《LS-DYNA 材料的真實應力應變曲線是怎么回事？》貼對于怎么計算仿真得出的應力應變曲線說的很含糊，且將文本文件中的數據粘貼至excel也較為繁瑣，在給學生講課的時候，普遍反映不感興趣。為簡化此問題，特地制作一個python腳本，用于自動讀取spcforc文件中的數據，并輸出名義應力應變曲線，增加學習的趣味性。代碼如下。腳本在附件中，有興趣的朋友可以嘗試下載。
記得點贊，謝謝。
熱處理消除Q235鋼焊接殘余應力的研究 附Q235鋼真實應力應變曲線研究下載
下載地址：Q235鋼真實應力應變曲線研究
abaqus子程序VUSDFLD——考慮應變率與應變軟化效應的軟土模型 ￥25
<p><strong>【注意】本貼子只包含子程序文件</strong></p><p>基于<a href="https://www.yqgqt.org.cn/qa/6302" rel="noopener noreferrer" target="_blank">abaqus子程序</a>VUSDFLD編寫的由Einav與Randolph提出的西澳模型，用于求解軟黏土體劇烈變形后的強度變化，可應用于的大變形計算。</p><div contenteditable="false" width="100%"><img src="https://img.jishulink.com/upload/202306/f69f50d42a81489ea1cb5e7a03da5c14.png" title="8$U(VZ82]O{OEMQB}[P(ZMB.png" alt="8$U(VZ82]O{OEMQB}[P(ZMB.png" style="max-width:760px;" data-mobile-src="https://img.jishulink.com/upload/202306/f69f50d42a81489ea1cb5e7a03da5c14.png?image_process=/format,webp" data-pc-src="https://img.jishulink.com/upload/202306/f69f50d42a81489ea1cb5e7a03da5c14.png?
展開 
ABAQUS中求解某部分單元的平均應力或平均應變 ￥10
1、參考模型：單向纖維的RVE模型；
2、腳本功能：針對指定的單元集合，在后處理中求解平均應力和平均應變。
3、應用的公式：一階均勻化計算方法。對于 RVE 模型的平均真應力和平均真應變，可通過對 RVE 內每一個單元的真應力 （真應變）取均值獲得。使用一階均勻化計算方法輸出的應力和應變適用于各種邊界條件，但需要對每個單元進行應力（應變）的輸出和計算。
有限元理論基礎及Abaqus內部實現方式研究系列27： Abaqus內部計算和顯示的應變
（1）顯示應變：Abaqus計算完畢后得到導入結果，在后處理中查看，應變E11=8.528e-1,E22=-5.173e-1如下：
（2）計算應變：Abaqus中采用UMAT子程序，利用我們的子程序調試插件DUS調試UMAT，在Visual Studio中查看dStran的值，發現在計算完應變后，進入UMAT時，E11=8.528e-1,E22=-5.173e-1，調試如下：
可以發現殼單元Abaqus的計算應變和顯示應變一樣，猜測都是對數應變。
1.5.3 iSolver的應變
iSolver中采用自帶材料進行計算，材料參數和UMAT的輸入完全一致。
為了計算和Abaqus完全一致，iSolver也采用對數應變計算方式，得到的應變顯示如下，可發現和Abaqus完全一致。
==總結==
由上可以看到，在實際計算中，對體單元，Abaqus和iSolver都采用變形率積分方式來計算應變，對殼單元，Abaqus和iSolver都采用對數應變。一般理論書都認為Abaqus是因為對數應變計算復雜才采用別的應變，但個人認為應該不是這個原因，因為Abaqus對體單元為了顯示對數應變，依然重新計算了一遍，說明Abaqus體單元采用變形率是有其它原因的，具體什么原因我也沒研究清楚，歡迎探討。
如果有任何其它疑問或者項目合作意向，也歡迎聯系我們：
snowwave02 From www.yqgqt.org.cn
email: snowwave02@qq.com
以往的系列文章：
1.7.1 ========第一階段========
第一篇：S4殼單元剛度矩陣研究。
http://www.yqgqt.org.cn/content/post/338859
第二篇：S4殼單元質量矩陣研究。
展開 ABAQUS批量提交Job與Python讀取ODB結果應力應變數據
批量提交的核心是需要等待當前計算任務的結束，上圖中若干個job的提交代碼如下：
# coding: utf-8
#微信公眾號：ABAQUS二次開發
#作者：阿信老師CAE
#email：axin_cae@163.com
#2022.03.17
from abaqus import *
from abaqusConstants import *
jobList = mdb.jobs.keys() #獲取所有計算任務的jobName
for jobName in jobList:
mdb.jobs[jobName].submit(consistencyChecking=OFF) #提交計算
mdb.jobs[jobName].waitForCompletion() #等待計算完成
print jobName , "is completed"
不過在實際的操作過程中，我們可能會需要避開一些job不提交，或者job窗口事先并沒有建立job，等等，總之實際問題永遠比任何教程都復雜，不過只需要靈活面對就可以了，處理起來并不難。
展開 Abaqus通過USDFLD子程序進行泥巖的應變軟化模擬
本文基于泥巖的三軸壓縮試驗曲線，建立考慮應變軟化特性的泥巖彈塑性本構模型，使用Abaqus及其子程序對泥巖的三軸壓縮試驗進行了數值模擬。泥巖在受壓過程中主要經過了5個階段，即壓密階段、彈性變 形、應變硬化、應變軟化、殘余階段。泥巖應變軟化模型如下所示。
式中，ξ為強度參數，ξp為峰值強度參數，ξr為殘余階段強度參數，η為應變軟化參數，η*為殘余階段的應變軟化參數初始值。對于三軸壓縮試驗，η用塑性剪切應變來表示
塑性屈服準則采用Mohr-Coulomb準則，則粘聚力和內摩擦角的參數演化可以用下式表示
在巖石的塑性變形過程中會產生比較明顯的剪脹現象，而用來描述這一現象的較常用的力學參數就是剪脹角 Ψ，
上述模型可以通過USDFLD子程序進行實現，流程圖如下
有限元模型如下圖所示
計算得到的應力云圖及不同圍壓下的載荷位移響應如下所示
參考文獻：張力偉,賈善坡,鄒江濤,舒婧曦.泥巖的峰后軟化力學模型.中國科技論文,2016,11(21):2456-2461
有關于abaqus子程序開發的相關問題可以通過公眾號聯系我們.
公眾號: 320科技工作室
展開 