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殼單元質量矩陣研究。
展開 計算層合板的ABD陣/柔度陣/應力應變---abaqus插件 ￥29.9
Graphite/epoxy材料參數：
E1=181Gpa E2=10.3Gpa G12=7.17Gpa Nu12=0.28
例題中的剛度陣為：
插件輸出的剛度陣：
例題中的中面應變：
插件輸出的中面應變：
全局應變：
全局應力：
局部應變：
局部應力：
插件輸出的鋪層的全局和局部的應力應變：
對比上面的結果可以發現插件的計算結果和例題的計算結果是相吻合的。
由工程應力、應變計算真實應力、應變
設t時刻的作用于試樣標距段加載方向的載荷為F，工程應力為，工程應變為，真實應力為，真實應變為。規定試樣受拉伸長時，載荷、應力及應變取正值，反之，受壓縮短時，載荷、應力及應變取負值，則
式1
式2
式中表示0時刻至t時刻試樣長度的增量。
式3
式4
真實應力 的求解變換中利用了材料變形過程中體積不變的假設，即。
繪制應力-應變曲線時，往往不管拉伸或者壓縮，都將應力和應變繪制成正值。這樣，拉伸時，按照公式<1>至<4>計算出的應力和應變均為正值，不需要進行變換；壓縮時，按照上述方法計算出的應力和應度均為負值，需要進行變換。變換方式為：對上述公式中所有的應力和應變乘-1。按照這種規則，壓縮時應力、應變用<5>至公式<8>進行汁算。其中公式<5>中F取負值。
式5
式6
式7
式8
?
?
?
展開 由工程應力、應變計算真實應力、應變
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?
展開 Python提取場輸出結果計算平均應力應變
###########################################################################
寫在前面，最近回看了一些自己之前編寫的和這個py代碼，發現自己寫的可能和公式表達的有偏差，在我的測試過程中，提取的是單元積分點處的應力應變值，由于只有一個積分點（這個積分點在質心，與centroid提取得到的結果一樣），因此自然的把這個值當成了單元的平均應力或平均應變，然后進行計算，但是我現在覺得這并不是獲取單元平均應力的方式，也就是代碼并沒有實現所謂的提取平均應力應變的功能，希望有大神可以指點迷津。
如果是有多個積分點的話，是不是應該對每個積分點權重進行積分，加權平均這樣得到單元的平均應力，然后乘單元體積，將所有單元的值求和再除模型的總體積，就得到整個RVE模型的平均應力。
###########################################################################
最近寫了一個簡單的python讀取abaqus結果中的場輸出數據，想通過均勻化計算方法來計算所定義集合的平均應力應變曲線，之前是手動提取了各個數據導出，然后用excel、matlab處理，但是很慢，而且很費勁，于是就想著用Python來處理結果。
有需要的同學可以下載附件文件，打開abaqus，file→run script，選擇腳本文件即可運行。
average.zip
均勻化計算方法：
參考文獻：馬思鳴. 精沖用碳鋼微觀組織對宏觀力學性能及精沖性能影響研究[D]. 上海：上海交通大學，2017.
我在網上找了挺久挺多的關于Python提取場輸出結果的，
主要對以下幾篇帖子進行了參考：
http://forum.simwe.com/forum.php?
展開 ABAQUS中求解某部分單元的平均應力或平均應變 ￥10
1、參考模型：單向纖維的RVE模型；
2、腳本功能：針對指定的單元集合，在后處理中求解平均應力和平均應變。
3、應用的公式：一階均勻化計算方法。對于 RVE 模型的平均真應力和平均真應變，可通過對 RVE 內每一個單元的真應力 （真應變）取均值獲得。使用一階均勻化計算方法輸出的應力和應變適用于各種邊界條件，但需要對每個單元進行應力（應變）的輸出和計算。
平面應變單元CPE4R齒輪傳動接觸應力計算 ￥49.9
厚齒輪的應力符合平面應變狀態，可以采用平面應變單元CPE4R來進行快速接觸應力計算。
在sketch模塊建立非對稱結構齒輪的草圖，然后建立part，并在assembly模塊進行裝配。
非對稱齒輪草圖
齒輪裝配體
通過適當的結構設計，非對稱齒輪可以在定速轉動的情況下獲得按某規律的變化轉速，在工程上經常會用到。
非對稱齒輪傳動分析結果
非對稱齒輪應力云圖
非對稱齒輪齒合區域局部應力云圖
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科技工作室
展開 Abaqus通過USDFLD子程序進行泥巖的應變軟化模擬
本文基于泥巖的三軸壓縮試驗曲線，建立考慮應變軟化特性的泥巖彈塑性本構模型，使用Abaqus及其子程序對泥巖的三軸壓縮試驗進行了數值模擬。泥巖在受壓過程中主要經過了5個階段，即壓密階段、彈性變 形、應變硬化、應變軟化、殘余階段。泥巖應變軟化模型如下所示。
式中，ξ為強度參數，ξp為峰值強度參數，ξr為殘余階段強度參數，η為應變軟化參數，η*為殘余階段的應變軟化參數初始值。對于三軸壓縮試驗，η用塑性剪切應變來表示
塑性屈服準則采用Mohr-Coulomb準則，則粘聚力和內摩擦角的參數演化可以用下式表示
在巖石的塑性變形過程中會產生比較明顯的剪脹現象，而用來描述這一現象的較常用的力學參數就是剪脹角 Ψ，
上述模型可以通過USDFLD子程序進行實現，流程圖如下
有限元模型如下圖所示
計算得到的應力云圖及不同圍壓下的載荷位移響應如下所示
參考文獻：張力偉,賈善坡,鄒江濤,舒婧曦.泥巖的峰后軟化力學模型.中國科技論文,2016,11(21):2456-2461
有關于abaqus子程序開發的相關問題可以聯系扣扣1653004885或者關注cae320公眾號
展開 
Abaqus平均應力和應變提取 ￥80
利用python腳本對ODB文件中單元集里所有積分點的應力及應變進行自動提取并計算平均值
能夠得到每一幀的應力和應變平均值，并保存到CSV文件中
所得到的應力包括S11,S22,S33,S12,S13,S23以及Mises七個應力平均值，以及E11,E22,E33,E12,E13,E23六個應變平均值
ABAQUS 小應變分析(例3) 條形基礎或海洋淺基礎下壓模擬(Tresca 本構) ￥67
ABAQUS 小應變分析(例3) 條形基礎或海洋淺基礎下壓模擬(Tresca 本構)
條形基礎承載力是工程廣泛關注的問題，例如陸地條形基礎和海洋淺基礎。該模擬地基為飽和不排水的粘土，采用Tresca本構，粘土強度su = 15 kPa。條形基礎處理成剛體。最終數模結果顯示，條形基礎的無量綱承載力Nc0 = F/Asu 近似于 pi + 2 = 5.14, 與傳統理論解極好的契合。
建模過程及結果：
荷載及位移邊界條件
網格劃分
局部網格劃分
條形基礎的力位移曲線(已達到極限承載力)
地基的土體應力分布
地基的土體破壞模式
MATLAB/FORTRAN | 鍵基近場動力學（BBPD）動態松弛法實現準靜態單軸壓縮模擬（含預制裂隙），反力計算應力應變曲線 ￥119
鍵基 PD 理論基礎：嚴格遵循 BBPD 理論，涵蓋近場半徑（Horizon）確定、微模量計算及斷裂準則。
單軸壓縮工況：預設標準的單軸壓縮邊界條件，模擬材料在受壓狀態下的損傷演化。
應力應變曲線計算：通過反力計算試樣的應力應變曲線。
MATLAB/Fortran 編寫：代碼結構清晰，算法邏輯直觀，無須配置復雜的第三方環境，適合學習與二次開發。
損傷演化可視化：程序包含后處理模塊，可生成裂紋擴展路徑、損傷場分布圖。
參數可調：材料參數、幾何尺寸、離散間距及迭代終止條件均可靈活修改。
ABAQUS提取單元平均應力/應變 ￥10
利用python讀取odb文件（可一次讀取多個odb）生成csv（excel）文件。提供源文件，注釋詳細，可根據需要進行修改。
