abaqus計算應變能的案例
有限元理論基礎及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中的偽應變能（artificial strain energy:ALLAE）
經常在參考書上看到關于沙漏控制的描述：當偽應變能（ALLAE）超過內能（ALLIE）的10%以后，沙漏效應過大，仿真無效。
于是在6.13幫助文檔查了一下關于ALLAE的描述，具體目錄為Getting Started with Abaqus: Interactive Edition中的10.5.2節，見附圖1。
可見ALLAE是abaqus用來控制沙漏效應的一種手段，當該手段與真實內能ALLIE的比值超過2%后，仿真失真，應該細化或調整網格以減小沙漏效應。
關于能量的進一步描述參見Getting Started with Abaqus: Interactive Edition中的9.6節Energy balance，相關文檔已上傳。
Getting Started with Abaqus_ Interactive Edition (6.13).pdf
附圖 1：
展開 計算層合板的ABD陣/柔度陣/應力應變---abaqus插件 ￥29.9
Graphite/epoxy材料參數：
E1=181Gpa E2=10.3Gpa G12=7.17Gpa Nu12=0.28
例題中的剛度陣為：
插件輸出的剛度陣：
例題中的中面應變：
插件輸出的中面應變：
全局應變：
全局應力：
局部應變：
局部應力：
插件輸出的鋪層的全局和局部的應力應變：
對比上面的結果可以發現插件的計算結果和例題的計算結果是相吻合的。
讀取ABAQUS結果文件中的單元應變能并輸出至excel文件的處理腳本 ￥1
腳本內容如題。
abaqus做DCB模擬時如何導出應變能釋放率曲線，可以不編程嗎？
寫論文需要R曲線（斷裂韌性-裂紋長度），但是不會編程
abaqus批處理提交計算3個完成后自動退出不能繼續進行的解決辦法
abaqus批處理提交計算3個完成后自動退出不能繼續進行的解決辦法.pdf
一個奇怪的問題（很多人沒有遇到過）：
Abaqus進行批處理時，批處理文件中有超過3個job時，第三個完成后自動退出，后續的不再進行了，必須設置多個批處理文件才能完成更多作業（每個批處理文件中不超過3個job），最初的設置如下（大家通常使用的）：
只能計算出A1/A2/A3， A4/A5不再進行計算。
解決辦法：
ABAQUS斷裂模擬收徒 ，保證快速學會各種ABAQUS斷裂模擬方法 1200/人（將享有各種插件以及程序，價值3000+、專門定制視頻、全程親自教學、各種模型調試及解答問題等等，傾囊相教）
展開 