abaqus應力應變提取的案例
Abaqus平均應力和應變提取 ￥80
利用python腳本對ODB文件中單元集里所有積分點的應力及應變進行自動提取并計算平均值
能夠得到每一幀的應力和應變平均值，并保存到CSV文件中
所得到的應力包括S11,S22,S33,S12,S13,S23以及Mises七個應力平均值，以及E11,E22,E33,E12,E13,E23六個應變平均值
ABAQUS提取單元平均應力/應變 ￥10
利用python讀取odb文件（可一次讀取多個odb）生成csv（excel）文件。提供源文件，注釋詳細，可根據需要進行修改。
Abaqus插件——平均應力應變提取 ￥60
通過該插件可實現:
1）提取所有幀的任意單元集合的平均應力（事先定義單元集合，如圖中的SET-1）
2）提取所有幀的任意區域的x、y、z方向的平均應變（事先定義節點集合，如圖中的SET-2）
3）將以上數據保存至excel文件（excel文件名為odb文件名稱+_Stress_Strain.csv）
*************************注意事項******************************
1、插件使用過程中，如有任何問題請發郵件至shenz1hao@126.com
2、插件僅做學習交流使用，尊重原創者，切勿以營利目的傳播
*****************************************************************
********************插件安裝及使用*******************************
1、電腦路徑下輸入 %homepath%\abaqus_plugins并回車
2、將Stress-Strain文件夾解壓至當前目錄下
3、打開abaqus，菜單欄中點擊plug-ins，里面找出Stress-Strain
4、輸入相應參數（hx、hy、hz表示x、y、z方向模型長度，當以上三參數取1時輸出的為該方向位移）
*****************************************************************
展開 Python提取場輸出結果計算平均應力應變
###########################################################################
寫在前面，最近回看了一些自己之前編寫的和這個py代碼，發現自己寫的可能和公式表達的有偏差，在我的測試過程中，提取的是單元積分點處的應力應變值，由于只有一個積分點（這個積分點在質心，與centroid提取得到的結果一樣），因此自然的把這個值當成了單元的平均應力或平均應變，然后進行計算，但是我現在覺得這并不是獲取單元平均應力的方式，也就是代碼并沒有實現所謂的提取平均應力應變的功能，希望有大神可以指點迷津。
如果是有多個積分點的話，是不是應該對每個積分點權重進行積分，加權平均這樣得到單元的平均應力，然后乘單元體積，將所有單元的值求和再除模型的總體積，就得到整個RVE模型的平均應力。
###########################################################################
最近寫了一個簡單的python讀取abaqus結果中的場輸出數據，想通過均勻化計算方法來計算所定義集合的平均應力應變曲線，之前是手動提取了各個數據導出，然后用excel、matlab處理，但是很慢，而且很費勁，于是就想著用Python來處理結果。
有需要的同學可以下載附件文件，打開abaqus，file→run script，選擇腳本文件即可運行。
average.zip
均勻化計算方法：
參考文獻：馬思鳴. 精沖用碳鋼微觀組織對宏觀力學性能及精沖性能影響研究[D]. 上海：上海交通大學，2017.
我在網上找了挺久挺多的關于Python提取場輸出結果的，
主要對以下幾篇帖子進行了參考：
http://forum.simwe.com/forum.php?
展開 
ANSYS瞬態分析全時程結構響應最大值的提取方法（變形、應力、應變、能量） ￥100
</p><p>同樣的方法，可以提取全時程最大的位移、應力、應變、能量等結果。</p><p><img src="https://img.jishulink.com/upload/202302/623025b5c0d646b9973cd2adc6c6037f.png" alt="1.png"></p><p>收費內容為相關命令流。</p>
展開 基于Python腳本提取復合材料應力應變均勻化有效彈性模量 ￥75
通過均勻化方法使用Python腳本從細觀尺度得到材料宏觀的等效模量
https://www.bilibili.com/video/BV1sb4y1273j
具體腳本使用方法及理論介紹詳見視頻教程
ABAQUS中求解某部分單元的平均應力或平均應變 ￥10
1、參考模型：單向纖維的RVE模型；
2、腳本功能：針對指定的單元集合，在后處理中求解平均應力和平均應變。
3、應用的公式：一階均勻化計算方法。對于 RVE 模型的平均真應力和平均真應變，可通過對 RVE 內每一個單元的真應力 （真應變）取均值獲得。使用一階均勻化計算方法輸出的應力和應變適用于各種邊界條件，但需要對每個單元進行應力（應變）的輸出和計算。
ABAQUS實用子程序SPRINC提取主應力
實用子程序SPRINC
在ABAQUS中使用UMAT子程序時有時會使用到最大主應力進行計算。通過查閱幫助文檔，ABAQUS實用子程序SPRINC可以在UMAT中計算最大主應力和最大主應變，SPRIND可以計算最大主應力和最大主應變的方向。
下面是ABAQUS幫助文檔關于實用子程序SPRINC的介紹：
SPRINC (calculate principal values)
Interface
CALL SPRINC(S,PS,LSTR,NDI,NSHR)
Variables to be provided to the utility routine
S
Stress or strain tensor.
LSTR
An identifier. LSTR=1 indicates that S contains stresses; LSTR=2 indicates that S contains strains.
NDI
Number of direct components.
NSHR
Number of shear components.
Variables returned from the utility routine
PS(I), I=1,2,3
The three principal values.
展開 ABAQUS中對應力、應變的部分理解
對應力的部分理解
對應變的部分理解
轉自公眾號——ABAQUS大世界
旨在分享，若侵即刪.
ABAQUS學習筆記—對應力應變的部分理解
之前關于后處理的一些文章，由于一些原因全部刪除，故今天開始重新開始分享一些關于ABAQUS的一些知識，希望能夠對大家有所幫助，也希望大家能夠繼續支持筆者。
那么今天，我們再對‘’ABAQUS中應力應變的部分理解的‘’內容進行講解。
在ABAQUS中，一般是把X軸當成1軸，Y軸當成2軸，Z軸當成3軸，那么：
S11就是X軸向的應力，正值為拉應力，負值為壓應力；
S22就是Y軸向的應力，正值為拉應力，負值為壓應力；
S33就是Z軸向的應力，正值為拉應力，負值為壓應力；
S12就是在YZ平面上，沿Y向的剪力；
S13就是在YZ平面上，沿Z向的剪力；
S23就是在XZ平面上，沿Z向的剪力；
由于剪力的對稱性:S11=S21,S13=S31,S23=S32；
由以上可知，S11，S22，S33為主應力；S12，S13，S23為切應力；
主應力分別以σ1，σ2，σ3表示，按數值排序為：σ1≥σ2≥σ3。在ABAQUS中分別對應為：Max.principal；Mid.principal；min.principal。這三個量在任何坐標下都是不變量。
我們可利于最大應力判斷一些情況：比如最大主應力（拉應力）大于混凝土的抗拉強度，則認為混凝土開裂；通過顯示最大主應力的法線方向，則可大致表示出裂縫的發展影響。
應變中一些符號的含義
E——總應變
EP——主應變
EE——彈性形變
PE——塑性應變分量
Eij——應變分量
Ie——非彈性應變分量
PEEQ——等效塑性應變。若該值大于0，則認為已經屈服
注：在ABAQUS后處理中，盡量不要看Mises，其表示平均應力，更適合金屬材料；對于鋼筋混凝土結構，我們要看其單軸拉伸方向上的應力和對應的應變。
展開 在abaqus中采用python提取結果-平均應力 ￥5
Python腳本操作ODB文件進行計算結果數據的提取。
因對結果分析需求，需提取某單元集的每一分析步（包含每一子步）的平均應力，目前網上雖有很多代碼是關于單元集的平均應力的提取，但并未有針對每一子步都需要提取結果的代碼。故針對此需要編寫了python代碼。
ABAQUS變量解讀：教你讀懂應力/應變/損傷
其中第一、第二、第三主應力分別是此狀態下最大、中間、最小的應力值。
ABAQUS提供了以下幾種主應力：
Max/Mid/Min Principal：第一、二、三主應力，分別對應最大、中間、最小主應力。在判斷以脆性材料為主的第一強度理論時有奇效。例如在砌體結構主壓破壞識別中，如下圖所示：
Max/Min In-Plane Principal：平面問題最大/最小主應力。
Max Principal(abs)：絕對值最大主應力。
既然講到強度理論，喵星人就帶大家回憶一下材料力學四大強度理論：
3.應力張量
如果你關心具體方向上的應力分量，可以選擇S11、S22、S12等張量分量。具體方向與單元類型相關，細節可參考我之前的文章《喵星人教你看懂不同類型單元的應力方向》。
二、應變相關
根據用戶手冊及后處理分類，ABAQUS提供了三類典型的后處理變量：
1.不變量
喵星人認為以下幾種應變相關不變量相對比較重要：
PEEQ：等效塑性應變，描述材料塑性變形的絕對值累積，用戶手冊定義如下：
典型的等效塑性應變區域如下圖：
PEMAG：塑性應變幅值，用戶手冊定義如下：
可見，相比PEEQ，PEMAG不考慮累積效應，兩種僅在單調加載下相同。
AC YIELD：標識“正在屈服”的區域，輸出0或1。適合用單色云圖顯示塑性區，如下圖所示。
2.
展開 ABAQUS批量提交Job與Python讀取ODB結果應力應變數據
問題2
如何從ODB結果文件中讀取計算結果數據
第二個問題是大家經常要面對的，常規的數據提取采用ABAQUS自帶的后處理工具已經能滿足需求了，但是如果我們需要進行一些寫出文本，數據篩選，數據的二次計算處理。。。，我們就需要用程序語言的幫助了。
針對今天的問題——
提取結果并寫出到文本
，如果用自帶的工具，可以導出abaqus.rpt文件，不過這種方式并不能滿足我們較為苛刻的要求，用python程序會的自由度會更大一些。
另外需要說明的是，如果我們僅僅只是看看某單元或者節點的各種計算結果，或者簡要的數據處理，臨時用一用，可以采用如下方式，不必費心寫復雜的程序，不值當浪費時間，
青春易逝，頭發寶貴
。當然了，通過阿信這里給出的基本程序改一改，能省幾根頭發。
限于時間關系，這里只列出基本的應力提取程序，寫入文本比較簡單，隨意百度即可實現，不再給出。
展開 一文搞清ABAQUS中真實應力和真實應變
應變計算公式為：
②仿真中的真實應力與真實應變
在使用ABAQUS仿真時，如果我們的材料屬于塑性材料范疇，分析時涉及較大變形，在分析時必須將其應力和應變定義成真實應力和真實應變，我們就需要將實驗數據中得到的名義數據轉化成真實數據。
總結來說，如果不使用真實應力與真實應變，我們就不能分析出我們想要的較為準確結果，這里考慮的是材料的非線性問題。
Matlab“稍作修改”Abaqus-odb結果 | 提取三維裝配體Mises應力
數據及函數文件獲取方式：
Matlab“稍作修改”Abaqus-odb結果 | 提取三維裝配體Mises應力
https://mp.weixin.qq.com/s/l45KSgNxkJ-oCCrDMrUgHA
