abaqus輸出應變的案例
abaqus拉伸后處理（應力-應變，位移-力的輸出）
abaqus拉伸后處理（應力-應變，位移-力的輸出）
lashen.zip
讀取ABAQUS結果文件中的單元應變能并輸出至excel文件的處理腳本 ￥1
腳本內容如題。
【Ls-dyna】Hypermesh&Ls-dyna聯合仿真時如何設置輸出單元應變和查看單元應變？
通常，使用Ls-dyna進行瞬態分析時，默認的輸出控制中是不包含單元的應變值的。但是，如果關心單元上的應變，或者需要對單元應變有更全面的了解，那么就需要在建模和分析中進行輸出控制，使計算結果中包含單元的應變。那么，問題來了，如何在Hypermesh中如何設置可以輸出單元的應變？如何在LS-POST中顯示單元的應變云圖？ 接下來，依次解答。
如何在Hypermesh中如何設置可以輸出單元的應變？
這里需要首先說明的是，本人使用Ls-dyna以來，一直使用的前處理不是Ls-prepost前處理，而是前處理軟件Hypermesh，后處理則使用Ls-prepost。所以，這里只說明Hypermesh&Ls-dyna聯合仿真時如何在Hypermesh中設置控制應變輸出。
前處理軟件工作界面
后處理軟件工作界面
首先啟動Hypermesh，點擊“user Profiles”，選擇Ls-dyna，其他默認，進入Ls-dyna分析模塊。
然后，在軟件的面板區域選擇“analysis”，并點擊“control card”。
點擊面板區域的Next，一直到出現“database-extent-binary”，點擊該按鈕
在關鍵字*Database-Extent-Binary的定義中，將第一行第四個參數【STRFLAG】的值設置為1，表示在二進制結果文件d3plot中輸出單元的應變。
在Hypermesh中按照上述方式就可以在d3plot中輸出單元的應變，由于d3plot是二進制文件，所以只能借助后處理LS-Prepost查看應變結果。
那么，如何在LS-POST查看應變，顯示應變云圖呢？
打開LS-PrePost，點擊【Fcomp】，選擇【strain】，就可以觀察計算后的應變云圖。
展開 輸出方法：abaqus最值隨時間變化
最近做項目，需要輸出abaqus中最大應變值隨時間變化的曲線，但是翻了很多帖子都沒找到相關的，所以后來自己寫了個Python代碼，用以提取最大值。需要用到這個功能的同學把代碼里的注釋刪掉直接運行就行。
viewportName = session.currentViewportName
j=getInput('how much frame','0') #要提取多少幀的最值，默認為0
#這一步也可以用在step上，我的需求始終是step-1，所以沒寫那部分代碼，和幀數的代碼同理
a=int(j)+1
i=1
while i<a:
cobject = session.viewports[viewportName].odbDisplay.setFrame(step=0, frame=i)
#有興趣的朋友可以print一下cobject1，看看odbDisplay里面包含了哪些函數
i+=1
cobject = session.viewports[viewportName].odbDisplay.contourOptions
#同上，想學習的朋友可以print一下cobject,
print("step=1", "frame=", i-1)
print(cobject.autoMaxValue)#輸出最大值
print(cobject.autoMinValue)#輸出最小值
值得注意的是，代碼里step和frame的值都比實際值要小1，比如我要提取本來是step-1的數值，但是代碼里要寫成step=0。當然我的代碼并不完整，這樣輸出的max/min只顯示在abaqus底部信息欄，還需要在Excel中處理數據。
展開 
py源代碼|平均應力應變位移輸出至Excel
—最大主應力；Smin—最小主應力；E11—x方向應變；E22—y方向應變；E33—z方向應變；Emax—最大主應變；U11—x方向位移；U22—y方向位移；U33—z方向位移；RF1—x方向支反力；RF2—y方向支反力；RF3—z方向支反力
output_constants.py、output_main.py獲取方式如下
公眾號(有限元與力學)任意兩篇文章轉發朋友圈12h后截圖加客服(QQ:416140933)索取
Python提取場輸出結果計算平均應力應變
###########################################################################
寫在前面，最近回看了一些自己之前編寫的和這個py代碼，發現自己寫的可能和公式表達的有偏差，在我的測試過程中，提取的是單元積分點處的應力應變值，由于只有一個積分點（這個積分點在質心，與centroid提取得到的結果一樣），因此自然的把這個值當成了單元的平均應力或平均應變，然后進行計算，但是我現在覺得這并不是獲取單元平均應力的方式，也就是代碼并沒有實現所謂的提取平均應力應變的功能，希望有大神可以指點迷津。
如果是有多個積分點的話，是不是應該對每個積分點權重進行積分，加權平均這樣得到單元的平均應力，然后乘單元體積，將所有單元的值求和再除模型的總體積，就得到整個RVE模型的平均應力。
###########################################################################
最近寫了一個簡單的python讀取abaqus結果中的場輸出數據，想通過均勻化計算方法來計算所定義集合的平均應力應變曲線，之前是手動提取了各個數據導出，然后用excel、matlab處理，但是很慢，而且很費勁，于是就想著用Python來處理結果。
有需要的同學可以下載附件文件，打開abaqus，file→run script，選擇腳本文件即可運行。
average.zip
均勻化計算方法：
參考文獻：馬思鳴. 精沖用碳鋼微觀組織對宏觀力學性能及精沖性能影響研究[D]. 上海：上海交通大學，2017.
我在網上找了挺久挺多的關于Python提取場輸出結果的，
主要對以下幾篇帖子進行了參考：
http://forum.simwe.com/forum.php?
展開 基于HyperWorks膜單元輸出應力應變仿真分析與測試對標的研究與應用
常亮_基于HyperWorks膜單元輸出應力應變仿真分析與測試對標的研究與應用.pdf
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?
展開 ABAQUS中求解某部分單元的平均應力或平均應變 ￥10
1、參考模型：單向纖維的RVE模型；
2、腳本功能：針對指定的單元集合，在后處理中求解平均應力和平均應變。
3、應用的公式：一階均勻化計算方法。對于 RVE 模型的平均真應力和平均真應變，可通過對 RVE 內每一個單元的真應力 （真應變）取均值獲得。使用一階均勻化計算方法輸出的應力和應變適用于各種邊界條件，但需要對每個單元進行應力（應變）的輸出和計算。
abaqus二次開發：后處理批量提取場輸出和歷程輸出結果（源碼帶注釋） ￥198
abaqus二次開發：后處理批量提取場輸出和歷程輸出結果（購買后需要加V私聊獲取使用教程）
購買后，請加V：wzd_1021_
更多插件內容，請私聊
同時提供定制hypermesh/hyperview/abaqus前后處理插件開發。
Abaqus-高清圖片輸出 ￥4.99
<p>在使用 Abaqus 的過程中，輸出高清圖片是一個常見的需求。</p><p>高清圖片的本質是 DPI（<span style="background-color: rgba(18, 18, 18, 0);">Dots Per Inch，圖像每英寸長度內的像素點數）達到一定數值。</span>網上的很多資料，僅僅是調整輸出圖片的格式，或者調整 Abaqus 內的顯示設置，并沒有從根本上改善圖片的 DPI。</p><p>本帖提供了操作簡便的處理思路，可以精準的控制圖片的清晰度及尺寸，適用于 原始文件（cae文件）的幾何模型、裝配模型以及結果文件（odb文件）的應力云圖、變形圖等圖片的高清輸出。</p><p>高清圖片導出--圖文流程，見付費部分。</p>
展開 
有限元理論基礎及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】輸出任意參量平均值的插件
image_process=/format,webp/resize,w_579" data-initial-src="https://img.jishulink.com/202405/attachment/ff59c2bd1e5247c6965f3dcd3a0a898c.png">
</figure>
</div><p class="ql-align-justify">即模型中每一個單元的應力（應變）對單元體積積分后，除以模型整個體積。上述應力應變曲線也證實，采用該方法能夠得到較為真實的數據。</p><p class="ql-align-justify">如果進提取某個方向的平均應力-應變，也許采用作用力反作用關系更快些；然而，有時候需要不同參量，如von mises應力等的平均值，此時需要通過python進行批量處理。</p><p class="ql-align-justify">****基于以上需求，做了各abaqus后處理插件，可以提取任意參量的平均值。
展開 【ABAQUS模態動力學】Composite&abaqus 預應力模態分析&輸出單元剛度矩陣
從上面這個理解出發，ABAQUS預應力模態只要在frequency分析步之前進行General,Static分析步,打開NLGeom選項(分析過程中剛度矩陣會不斷變化)。
提取單元剛度矩陣：
【ABAQUS 二次開發筆記】輸出單元剛度矩陣 - hayden_william - 博客園
以上均為我的一點理解，不一定完全正確，本文僅作為個人學習記錄之用，其他概不負責。
ABAQUS批量提交Job與Python讀取ODB結果應力應變數據
# coding: utf-8#微信公眾號：ABAQUS二次開發#作者：阿信老師CAE#email：axin_cae@163.com#2022.03.17
from odbAccess import *from visualization import *
frame_num = 20 #定義需要輸出的幀odbName = 'X:/***/***.odb' #odb路徑和名稱elemSetName = '***' #單元setodb = openOdb(path=odbName)#打開odbregion = odb.rootAssembly.elementSets[elemSetName]frameRps = odb.steps['Step-1'].framesf_S1 = frameRps[frame_num].fieldOutputs['S'] #需要輸出的場變量sub_f_S1 = f_S1.getSubset(region=region) S_Value = sub_f_S1.values#讀取所有單元的應力值for i in range(len(S_Value)): Stress =S_Value[i].data print Stressodb.close()
若有緣，江湖再會。
阿信
2022.3.17
文章來源：ABAQUS二次開發
展開 