###########################################################################

寫在前面，最近回看了一些自己之前編寫的和這個py代碼，發現自己寫的可能和公式表達的有偏差，在我的測試過程中，提取的是單元積分點處的應力應變值，由于只有一個積分點（這個積分點在質心，與centroid提取得到的結果一樣），因此自然的把這個值當成了單元的平均應力或平均應變，然后進行計算，但是我現在覺得這并不是獲取單元平均應力的方式，也就是代碼并沒有實現所謂的提取平均應力應變的功能，希望有大神可以指點迷津。

如果是有多個積分點的話，是不是應該對每個積分點權重進行積分，加權平均這樣得到單元的平均應力，然后乘單元體積，將所有單元的值求和再除模型的總體積，就得到整個RVE模型的平均應力。

###########################################################################

最近寫了一個簡單的python讀取abaqus結果中的場輸出數據，想通過均勻化計算方法來計算所定義集合的平均應力應變曲線，之前是手動提取了各個數據導出，然后用excel、matlab處理，但是很慢，而且很費勁，于是就想著用Python來處理結果。

有需要的同學可以下載附件文件，打開abaqus，file→run script，選擇腳本文件即可運行。

average.zip

均勻化計算方法：

參考文獻：馬思鳴. 精沖用碳鋼微觀組織對宏觀力學性能及精沖性能影響研究[D]. 上海：上海交通大學，2017.

我在網上找了挺久挺多的關于Python提取場輸出結果的，
主要對以下幾篇帖子進行了參考：
http://forum.simwe.com/forum.php?mod=viewthread&tid=1283175&highlight=%E5%B9%B3%E5%9D%87%E5%BA%94%E5%8A%9B%E5%BA%94%E5%8F%98
http://forum.simwe.com/forum.php?mod=viewthread&tid=1075473&highlight=Python%E5%90%8E%E5%A4%84%E7%90%86

http://www.yqgqt.org.cn/content/post/296436

http://www.yqgqt.org.cn/content/post/1189653

http://www.yqgqt.org.cn/content/post/542728

但是我找到的案例中輸出的都是某一幀所有單元的或者某一單元所有幀的結果，
而我要想求的是所有單元每一幀的結果，以第二篇帖子中的簡單模型為案例，集合Set-1定義為整個模型，編寫程序計算輸出了其應力應變數據，結果如下：