abaqus讀取應(yīng)變的案例
ABAQUS批量提交Job與Python讀取ODB結(jié)果應(yīng)力應(yīng)變數(shù)據(jù)
問題2
如何從ODB結(jié)果文件中讀取計算結(jié)果數(shù)據(jù)
第二個問題是大家經(jīng)常要面對的，常規(guī)的數(shù)據(jù)提取采用ABAQUS自帶的后處理工具已經(jīng)能滿足需求了，但是如果我們需要進行一些寫出文本，數(shù)據(jù)篩選，數(shù)據(jù)的二次計算處理。。。，我們就需要用程序語言的幫助了。
針對今天的問題——
提取結(jié)果并寫出到文本
，如果用自帶的工具，可以導(dǎo)出abaqus.rpt文件，不過這種方式并不能滿足我們較為苛刻的要求，用python程序會的自由度會更大一些。
另外需要說明的是，如果我們僅僅只是看看某單元或者節(jié)點的各種計算結(jié)果，或者簡要的數(shù)據(jù)處理，臨時用一用，可以采用如下方式，不必費心寫復(fù)雜的程序，不值當(dāng)浪費時間，
青春易逝，頭發(fā)寶貴
。當(dāng)然了，通過阿信這里給出的基本程序改一改，能省幾根頭發(fā)。
限于時間關(guān)系，這里只列出基本的應(yīng)力提取程序，寫入文本比較簡單，隨意百度即可實現(xiàn)，不再給出。
展開 讀取ABAQUS結(jié)果文件中的單元應(yīng)變能并輸出至excel文件的處理腳本 ￥1
腳本內(nèi)容如題。
Abaqus二次開發(fā)讀取變形后的節(jié)點坐標(biāo)并輸出到txt文件中。讀取模態(tài)頻率到txt 文件中。 ￥10
有時候在abaqus中，我們需要知道某一個集合中的節(jié)點變形后的坐標(biāo)。以此為輸入來進行一些研究。這里我用一個自編函數(shù)，將節(jié)點集合變形后的坐標(biāo)寫入到一個txt文件中。格式化的寫入文件，方便用此坐標(biāo)來進行計算研究。
同樣有時候，需要將計算的模態(tài)頻率值提取出來。同樣用一個函數(shù)將模態(tài)頻率提取出來，放進txt文件中，方便后續(xù)研究。
讀取的節(jié)點結(jié)果如下圖所示：
Abaqus利用region讀取結(jié)果數(shù)據(jù)
Abaqus利用region讀取結(jié)果數(shù)據(jù)
長安CAE
Abaqus計算完成后，為輸出某部分模型的結(jié)果數(shù)據(jù)，需要用到region，并且配套使用的是getSubset函數(shù)。
示例一創(chuàng)建一個變量center，將一個叫做PUNCH的節(jié)點set賦值給這個變量，通過整個模型的位移場displacement得到該節(jié)點set的輸出數(shù)據(jù)，代碼如下：
center = odb.rootAssembly.instances['PART-1-1'].nodeSets['PUNCH']
centerDisplacement = displacement.getSubset(region=center)
centerValues = centerDisplacement.values
for v in centerValues:
print v.nodeLabel, v.data
其中的displacement是已經(jīng)創(chuàng)建好的位移場數(shù)據(jù)，centerDisplacement變量得到該節(jié)點集的所有數(shù)據(jù)，centerValues變量得到該節(jié)點集的數(shù)值數(shù)據(jù)，通過循環(huán)輸出節(jié)點號和當(dāng)前數(shù)據(jù)。
展開 
Abaqus利用region讀取結(jié)果數(shù)據(jù)
單元質(zhì)心點
轉(zhuǎn)自公眾號——ABAQUS大世界
旨在分享，若侵即刪.
Abaqus二次開發(fā)之Python腳本讀取場輸出位移
操作腳本如下：
結(jié)果如下：
Node label: 1000
Displacement in X direction: -7.891572e-34
Displacement in Y direction: -76.45553
Displacement: 76.4555282593
注：本文中涉及的模型文件和完整腳本文件可在Up的公眾號：易木木響叮當(dāng)，中回復(fù)“讀取場輸出數(shù)據(jù)”自動獲取。
ANSYS與Abaqus球坐標(biāo)系下的結(jié)果讀取
ANSYS與Abaqus球坐標(biāo)系下的結(jié)果讀取
1 概述
采用ANSYS和Abaqus軟件計算的結(jié)果通常默認(rèn)的結(jié)果是在總體笛卡爾坐標(biāo)系下產(chǎn)生的結(jié)果，這對于應(yīng)力或者應(yīng)變等分量的分析有時候不方便，比如對于一個圓筒體，比較關(guān)心其徑向應(yīng)力和環(huán)向應(yīng)力，而這個結(jié)果直接讀取使不可能的，需要一定的轉(zhuǎn)換。
這就是結(jié)果坐標(biāo)系轉(zhuǎn)換。
在軟件里，應(yīng)力分量表示為sx,xy,xz（ANSYS），s11,s22,s22（Abaqus），當(dāng)其轉(zhuǎn)換到柱坐標(biāo)或者球坐標(biāo)時，對應(yīng)的應(yīng)力分量就發(fā)生變化，sx和s11均表示徑向應(yīng)力。
2 ANSYS
建立一個球體模型，如圖1，加載求解，得到其總體坐標(biāo)系下的sx應(yīng)力分量。
圖1
在后處理器中，將結(jié)果坐標(biāo)系轉(zhuǎn)換為球坐標(biāo)系，采用的命令為：RSYS。查詢ANSYS幫助文檔，如圖2：
圖2 RSYS
0,1,2分別代表笛卡爾坐標(biāo)系，柱坐標(biāo)系，球坐標(biāo)系。
輸入命令：RSYS，2
顯式結(jié)果sx為圖3，此時的sx應(yīng)力分量為徑向應(yīng)力。
圖3
3 Abaqus
建立模型加載求解，得到s11應(yīng)力分量如圖4.
圖4
轉(zhuǎn)換結(jié)果坐標(biāo)系，Visualization模塊下選擇 Tools--Create Coordinate Aystem，按指定方法建立局部坐標(biāo)系，然后選擇Result-Option，選擇Transformation標(biāo)簽，User-specified，就可以看到新建立的坐標(biāo)系，選擇新建的坐標(biāo)系即可完成坐標(biāo)轉(zhuǎn)換。
如圖5，圖6
圖5
建立球坐標(biāo)系的時候根據(jù)Abaqus窗口下方的提示進行操作。
圖6
最終轉(zhuǎn)換為徑向應(yīng)力的顯式結(jié)果，如圖7
圖7
展開 Abaqus利用腳本讀取結(jié)果并創(chuàng)建新的輸出場
Abaqus利用腳本讀取結(jié)果并創(chuàng)建新的輸出場
長安CAE
Abaqus的后處理有時候采用Python腳本會比較方便，此次實例通過腳本語言讀取已經(jīng)存在的Odb結(jié)果文件，并且通過程序?qū)崿F(xiàn)不同載荷步下的結(jié)果場數(shù)據(jù)運算，以得到新的場數(shù)據(jù)，并顯示其云圖。
abaqus子程序VUSDFLD——考慮應(yīng)變率與應(yīng)變軟化效應(yīng)的軟土模型 ￥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提出的西澳模型，用于求解軟黏土體劇烈變形后的強度變化，可應(yīng)用于的大變形計算。</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中求解某部分單元的平均應(yīng)力或平均應(yīng)變 ￥10
1、參考模型：單向纖維的RVE模型；
2、腳本功能：針對指定的單元集合，在后處理中求解平均應(yīng)力和平均應(yīng)變。
3、應(yīng)用的公式：一階均勻化計算方法。對于 RVE 模型的平均真應(yīng)力和平均真應(yīng)變，可通過對 RVE 內(nèi)每一個單元的真應(yīng)力 （真應(yīng)變）取均值獲得。使用一階均勻化計算方法輸出的應(yīng)力和應(yīng)變適用于各種邊界條件，但需要對每個單元進行應(yīng)力（應(yīng)變）的輸出和計算。
有限元理論基礎(chǔ)及Abaqus內(nèi)部實現(xiàn)方式研究系列27： Abaqus內(nèi)部計算和顯示的應(yīng)變
（1）顯示應(yīng)變：Abaqus計算完畢后得到導(dǎo)入結(jié)果，在后處理中查看，應(yīng)變E11=8.528e-1,E22=-5.173e-1如下：
（2）計算應(yīng)變：Abaqus中采用UMAT子程序，利用我們的子程序調(diào)試插件DUS調(diào)試UMAT，在Visual Studio中查看dStran的值，發(fā)現(xiàn)在計算完應(yīng)變后，進入UMAT時，E11=8.528e-1,E22=-5.173e-1，調(diào)試如下：
可以發(fā)現(xiàn)殼單元Abaqus的計算應(yīng)變和顯示應(yīng)變一樣，猜測都是對數(shù)應(yīng)變。
1.5.3 iSolver的應(yīng)變
iSolver中采用自帶材料進行計算，材料參數(shù)和UMAT的輸入完全一致。
為了計算和Abaqus完全一致，iSolver也采用對數(shù)應(yīng)變計算方式，得到的應(yīng)變顯示如下，可發(fā)現(xiàn)和Abaqus完全一致。
==總結(jié)==
由上可以看到，在實際計算中，對體單元，Abaqus和iSolver都采用變形率積分方式來計算應(yīng)變，對殼單元，Abaqus和iSolver都采用對數(shù)應(yīng)變。一般理論書都認(rèn)為Abaqus是因為對數(shù)應(yīng)變計算復(fù)雜才采用別的應(yīng)變，但個人認(rèn)為應(yīng)該不是這個原因，因為Abaqus對體單元為了顯示對數(shù)應(yīng)變，依然重新計算了一遍，說明Abaqus體單元采用變形率是有其它原因的，具體什么原因我也沒研究清楚，歡迎探討。
如果有任何其它疑問或者項目合作意向，也歡迎聯(lián)系我們：
snowwave02 From www.yqgqt.org.cn
email: snowwave02@qq.com
以往的系列文章：
1.7.1 ========第一階段========
第一篇：S4殼單元剛度矩陣研究。
http://www.yqgqt.org.cn/content/post/338859
第二篇：S4殼單元質(zhì)量矩陣研究。
展開 
Abaqus平均應(yīng)力和應(yīng)變提取 ￥80
利用python腳本對ODB文件中單元集里所有積分點的應(yīng)力及應(yīng)變進行自動提取并計算平均值
能夠得到每一幀的應(yīng)力和應(yīng)變平均值，并保存到CSV文件中
所得到的應(yīng)力包括S11,S22,S33,S12,S13,S23以及Mises七個應(yīng)力平均值，以及E11,E22,E33,E12,E13,E23六個應(yīng)變平均值
ABAQUS提取單元平均應(yīng)力/應(yīng)變 ￥10
利用python讀取odb文件（可一次讀取多個odb）生成csv（excel）文件。提供源文件，注釋詳細，可根據(jù)需要進行修改。
ABAQUS中對應(yīng)力、應(yīng)變的部分理解
對應(yīng)力的部分理解
對應(yīng)變的部分理解
轉(zhuǎn)自公眾號——ABAQUS大世界
旨在分享，若侵即刪.
ABAQUS變量解讀：教你讀懂應(yīng)力/應(yīng)變/損傷
有限元后處理直接與數(shù)據(jù)圖片處理、論文撰寫相關(guān)，除了典型的應(yīng)力張量與應(yīng)變張量外，ABAQUS還提供了大量可供使用者讀取的其他應(yīng)力/應(yīng)變/損傷參數(shù)，這都有助于結(jié)果的分析。今天喵星人就教你讀懂其中的應(yīng)力、應(yīng)變及損傷的后處理細節(jié)。
一、應(yīng)力相關(guān)
根據(jù)用戶手冊及后處理分類，ABAQUS提供了三類典型的后處理變量：
1.不變量
不變量的定義是指張量在坐標(biāo)旋轉(zhuǎn)下保持不變的量。這些量反映了材料內(nèi)在的力學(xué)狀態(tài)，與觀察方向無關(guān)，因此常用于判斷材料的屈服、破壞或變形行為。
喵星人認(rèn)為以下幾種應(yīng)力相關(guān)不變量相對比較重要：
Mises：基于第四強度理論，用戶手冊定義如下：
Tresca：基于第三強度理論，用戶手冊定義如下：
Tresca equivalent stress, defined as the maximum difference between principal stresses.
Pressure：靜水壓力，注意正值為壓，負值為拉，用戶手冊定義如下：
編輯
跳轉(zhuǎn)
Third Invariant：第三應(yīng)力不變量，用戶手冊定義如下：
編輯
跳轉(zhuǎn)
編輯
跳轉(zhuǎn)
這里，喵星人給出更加簡潔的定義：
TRIAX：應(yīng)力三軸度，可用于評估斷裂行為，尤其適用于延性損傷模型，用戶手冊定義如下：
編輯
跳轉(zhuǎn)
2.主應(yīng)力
主應(yīng)力是指該點應(yīng)力狀態(tài)中僅存在正應(yīng)力、無剪應(yīng)力的特殊方向上的應(yīng)力值。
展開 