abaqus提取應力場的案例
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平均應力和應變提取 ￥80
利用python腳本對ODB文件中單元集里所有積分點的應力及應變進行自動提取并計算平均值
能夠得到每一幀的應力和應變平均值，并保存到CSV文件中
所得到的應力包括S11,S22,S33,S12,S13,S23以及Mises七個應力平均值，以及E11,E22,E33,E12,E13,E23六個應變平均值
python提取abaqus場輸出結果
最近寫了一個簡單的python讀取abaqus結果中的場輸出數據，就在這里簡單的講一下整個流程。后續采用abaqus讀取場輸出結果應用也挺廣泛。整個過程的流程圖可以用下圖來概括，因為不是計算機專業，考慮沒那么細致，很多python函數也是即查即用，當然程序后續根據要求也可以不斷更改。
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插件——平均應力應變提取 ￥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時輸出的為該方向位移）
*****************************************************************
展開 ABAQUS應力場
求大神告知abaqus中考慮計算出的應力場，是通過預定義場中的Initial-Mechanical-Stress，還是Step 1-Other-Field？
在abaqus中采用python提取結果-平均應力 ￥5
Python腳本操作ODB文件進行計算結果數據的提取。
因對結果分析需求，需提取某單元集的每一分析步（包含每一子步）的平均應力，目前網上雖有很多代碼是關于單元集的平均應力的提取，但并未有針對每一子步都需要提取結果的代碼。故針對此需要編寫了python代碼。
Abaqus_Python 后處理結點數據提取、場變量云圖輸出 ￥1.99
每個case需要輸出幾個路徑上的部分節點值，而且需要針對不同的case、不同的幀、不同的視角截圖、保存云圖，如果人工出圖和提取節點數據，工作量會非常大，費時會遠遠超出計算所用時間；并且，同事還要求我根據case名稱、組件名稱（模型有多個組件組成，且網格是在part上畫的，所以同一個節點編號可能對應多個節點）等保存成不同的文件。
結點溫度數據沒有保存，因此沒有輸出。下圖是自動輸出的png格式云圖。
比較倉促，代碼也比價粗糙。回頭做一個詳細的介紹。
主要代碼見下。
ABAQUS提取單元平均應力/應變 ￥10
利用python讀取odb文件（可一次讀取多個odb）生成csv（excel）文件。提供源文件，注釋詳細，可根據需要進行修改。
怎么在abaqus提取xfem下，動態三維應力強度因子？
求解答
ABAQUS怎么提取扣件支反力和結構層壓應力
ABAQUS怎么提取扣件支反力？軌道結構各結構層壓應力是提取S22還是最大主拉應力
Matlab“稍作修改”Abaqus-odb結果 | 提取三維裝配體Mises應力
數據及函數文件獲取方式：
Matlab“稍作修改”Abaqus-odb結果 | 提取三維裝配體Mises應力
https://mp.weixin.qq.com/s/l45KSgNxkJ-oCCrDMrUgHA
Abaqus CAE Python后處理提取每一幀最大等效應力
使用Python語言對Abaqus CAE后處理結果進行分析，并提取一個分析步中每一幀的最大等效應力，其中Python代碼如下：
from abaqus import *
from abaqusConstants import *
from odbAccess import *
import visualization
myFile=open('DATA.txt','w')
print('********************************\n')
myFile.write('********************************\n')
myOdb=openOdb(path='viewer_tutorial.odb')
myStepValue=myOdb.steps.values()
for step in myStepValue:
print('The current step is: %s.\n'%step.name)
myFile.write('The current step is: %s.\n'%step.name)
frameID=0
for frame in step.frames:
print('The current frame is: %d.\n'%frameID)
myFile.write('The current frame is: %d.
展開 ABAQUS殘余應力預定義場設置案例
六、小結
本文通過一個簡單的輪轂鑄造后所得殘余應力，并且在Abaqus中進行預定義場設置的方法和流程，僅供參考，實際制造過程的模型，網格劃分以及材料，裝配關系也更為復雜。
abaqus二次開發：后處理批量提取場輸出和歷程輸出結果（源碼帶注釋） ￥198
abaqus二次開發：后處理批量提取場輸出和歷程輸出結果（購買后需要加V私聊獲取使用教程）
購買后，請加V：wzd_1021_
更多插件內容，請私聊
同時提供定制hypermesh/hyperview/abaqus前后處理插件開發。
