abaqus輸出應(yīng)力應(yīng)變的案例
abaqus拉伸后處理（應(yīng)力-應(yīng)變，位移-力的輸出）
abaqus拉伸后處理（應(yīng)力-應(yīng)變，位移-力的輸出）
lashen.zip
py源代碼|平均應(yīng)力應(yīng)變位移輸出至Excel
/usr/bin/python
#-*-coding:utf-8-*-
import csv
import output_main
（3）以提取模型的分析時(shí)間、S11應(yīng)力和E11應(yīng)變數(shù)據(jù)為例，輸入下列代碼
time=[] #儲(chǔ)存時(shí)間數(shù)據(jù)
time=output_main.output_time(time)
stress=[] #儲(chǔ)存應(yīng)力數(shù)據(jù),其中'CONCRETE-1'部件名，'SET-1'單元集合名
stress=output_main.output_S11('CONCRETE-1','SET-1',stress)
strain=[] #儲(chǔ)存應(yīng)力數(shù)據(jù),其中'CONCRETE'部件名，'SET-1'單元集合名
stress=output_main.output_E11('CONCRETE-1','SET-1',strain)
（4）輸出數(shù)據(jù)至Excel
New=open('Output.csv','wb') #新建一個(gè)Excel文件儲(chǔ)存數(shù)據(jù)
New.write('time,stress,strain\n') #輸入Excel表頭分別為time，stress，strain
for i in range(len(stress)): #循環(huán)輸出數(shù)據(jù)
New.write('%s,%s,%s\n'%(time[i],stress[i],strain[i]))
New.close()
（5）保存新建txt文件，將.txt后綴修改為.py，在Abaqus中以腳本形式運(yùn)行該文件
output_main.output_S11中的S11用于控制輸出結(jié)果
結(jié)果控制參數(shù)如下
S11—x方向應(yīng)力；S22—y方向應(yīng)力；S33—z方向應(yīng)力；mises—mises應(yīng)力；Smax
展開 Python提取場(chǎng)輸出結(jié)果計(jì)算平均應(yīng)力應(yīng)變
###########################################################################
寫在前面，最近回看了一些自己之前編寫的和這個(gè)py代碼，發(fā)現(xiàn)自己寫的可能和公式表達(dá)的有偏差，在我的測(cè)試過程中，提取的是單元積分點(diǎn)處的應(yīng)力應(yīng)變值，由于只有一個(gè)積分點(diǎn)（這個(gè)積分點(diǎn)在質(zhì)心，與centroid提取得到的結(jié)果一樣），因此自然的把這個(gè)值當(dāng)成了單元的平均應(yīng)力或平均應(yīng)變，然后進(jìn)行計(jì)算，但是我現(xiàn)在覺得這并不是獲取單元平均應(yīng)力的方式，也就是代碼并沒有實(shí)現(xiàn)所謂的提取平均應(yīng)力應(yīng)變的功能，希望有大神可以指點(diǎn)迷津。
如果是有多個(gè)積分點(diǎn)的話，是不是應(yīng)該對(duì)每個(gè)積分點(diǎn)權(quán)重進(jìn)行積分，加權(quán)平均這樣得到單元的平均應(yīng)力，然后乘單元體積，將所有單元的值求和再除模型的總體積，就得到整個(gè)RVE模型的平均應(yīng)力。
###########################################################################
最近寫了一個(gè)簡(jiǎn)單的python讀取abaqus結(jié)果中的場(chǎng)輸出數(shù)據(jù)，想通過均勻化計(jì)算方法來計(jì)算所定義集合的平均應(yīng)力應(yīng)變曲線，之前是手動(dòng)提取了各個(gè)數(shù)據(jù)導(dǎo)出，然后用excel、matlab處理，但是很慢，而且很費(fèi)勁，于是就想著用Python來處理結(jié)果。
有需要的同學(xué)可以下載附件文件，打開abaqus，file→run script，選擇腳本文件即可運(yùn)行。
average.zip
均勻化計(jì)算方法：
參考文獻(xiàn)：馬思鳴. 精沖用碳鋼微觀組織對(duì)宏觀力學(xué)性能及精沖性能影響研究[D]. 上海：上海交通大學(xué)，2017.
我在網(wǎng)上找了挺久挺多的關(guān)于Python提取場(chǎng)輸出結(jié)果的，
主要對(duì)以下幾篇帖子進(jìn)行了參考：
http://forum.simwe.com/forum.php?
展開 基于HyperWorks膜單元輸出應(yīng)力應(yīng)變仿真分析與測(cè)試對(duì)標(biāo)的研究與應(yīng)用
計(jì)算機(jī)仿真技術(shù)越來越受到企業(yè)的認(rèn)可和重視，合理地建模是仿真分析成敗的關(guān)鍵，而驗(yàn)證建模的合理性，試驗(yàn)對(duì)標(biāo)是行之有效的方法，文中通過對(duì)焊接結(jié)構(gòu)、鑄造結(jié)構(gòu)、大型裝配體等零部件的應(yīng)力測(cè)試對(duì)標(biāo)為研究對(duì)象，以Altair公司的HypeWorks軟件為平臺(tái)，來尋找可行的對(duì)標(biāo)方法。
常亮_基于HyperWorks膜單元輸出應(yīng)力應(yīng)變仿真分析與測(cè)試對(duì)標(biāo)的研究與應(yīng)用.pdf
ABAQUS中求解某部分單元的平均應(yīng)力或平均應(yīng)變 ￥10
1、參考模型：?jiǎn)蜗蚶w維的RVE模型；
2、腳本功能：針對(duì)指定的單元集合，在后處理中求解平均應(yīng)力和平均應(yīng)變。
3、應(yīng)用的公式：一階均勻化計(jì)算方法。對(duì)于 RVE 模型的平均真應(yīng)力和平均真應(yīng)變，可通過對(duì) RVE 內(nèi)每一個(gè)單元的真應(yīng)力 （真應(yīng)變）取均值獲得。使用一階均勻化計(jì)算方法輸出的應(yīng)力和應(yīng)變適用于各種邊界條件，但需要對(duì)每個(gè)單元進(jìn)行應(yīng)力（應(yīng)變）的輸出和計(jì)算。
Abaqus平均應(yīng)力和應(yīng)變提取 ￥80
利用python腳本對(duì)ODB文件中單元集里所有積分點(diǎn)的應(yīng)力及應(yīng)變進(jìn)行自動(dòng)提取并計(jì)算平均值
能夠得到每一幀的應(yīng)力和應(yīng)變平均值，并保存到CSV文件中
所得到的應(yīng)力包括S11,S22,S33,S12,S13,S23以及Mises七個(gè)應(yīng)力平均值，以及E11,E22,E33,E12,E13,E23六個(gè)應(yīng)變平均值
ABAQUS提取單元平均應(yīng)力/應(yīng)變 ￥10
利用python讀取odb文件（可一次讀取多個(gè)odb）生成csv（excel）文件。提供源文件，注釋詳細(xì)，可根據(jù)需要進(jìn)行修改。
ABAQUS中對(duì)應(yīng)力、應(yīng)變的部分理解
對(duì)應(yīng)力的部分理解
對(duì)應(yīng)變的部分理解
轉(zhuǎn)自公眾號(hào)——ABAQUS大世界
旨在分享，若侵即刪.
Abaqus插件——平均應(yīng)力應(yīng)變提取 ￥60
通過該插件可實(shí)現(xiàn):
1）提取所有幀的任意單元集合的平均應(yīng)力（事先定義單元集合，如圖中的SET-1）
2）提取所有幀的任意區(qū)域的x、y、z方向的平均應(yīng)變（事先定義節(jié)點(diǎn)集合，如圖中的SET-2）
3）將以上數(shù)據(jù)保存至excel文件（excel文件名為odb文件名稱+_Stress_Strain.csv）
*************************注意事項(xiàng)******************************
1、插件使用過程中，如有任何問題請(qǐng)發(fā)郵件至shenz1hao@126.com
2、插件僅做學(xué)習(xí)交流使用，尊重原創(chuàng)者，切勿以營(yíng)利目的傳播
*****************************************************************
********************插件安裝及使用*******************************
1、電腦路徑下輸入 %homepath%\abaqus_plugins并回車
2、將Stress-Strain文件夾解壓至當(dāng)前目錄下
3、打開abaqus，菜單欄中點(diǎn)擊plug-ins，里面找出Stress-Strain
4、輸入相應(yīng)參數(shù)（hx、hy、hz表示x、y、z方向模型長(zhǎng)度，當(dāng)以上三參數(shù)取1時(shí)輸出的為該方向位移）
*****************************************************************
展開 讀取ABAQUS結(jié)果文件中的單元應(yīng)變能并輸出至excel文件的處理腳本 ￥1
腳本內(nèi)容如題。
ABAQUS學(xué)習(xí)筆記—對(duì)應(yīng)力應(yīng)變的部分理解
之前關(guān)于后處理的一些文章，由于一些原因全部刪除，故今天開始重新開始分享一些關(guān)于ABAQUS的一些知識(shí)，希望能夠?qū)Υ蠹矣兴鶐椭蚕Ｍ蠹夷軌蚶^續(xù)支持筆者。
那么今天，我們?cè)賹?duì)‘’ABAQUS中應(yīng)力應(yīng)變的部分理解的‘’內(nèi)容進(jìn)行講解。
在ABAQUS中，一般是把X軸當(dāng)成1軸，Y軸當(dāng)成2軸，Z軸當(dāng)成3軸，那么：
S11就是X軸向的應(yīng)力，正值為拉應(yīng)力，負(fù)值為壓應(yīng)力；
S22就是Y軸向的應(yīng)力，正值為拉應(yīng)力，負(fù)值為壓應(yīng)力；
S33就是Z軸向的應(yīng)力，正值為拉應(yīng)力，負(fù)值為壓應(yīng)力；
S12就是在YZ平面上，沿Y向的剪力；
S13就是在YZ平面上，沿Z向的剪力；
S23就是在XZ平面上，沿Z向的剪力；
由于剪力的對(duì)稱性:S11=S21,S13=S31,S23=S32；
由以上可知，S11，S22，S33為主應(yīng)力；S12，S13，S23為切應(yīng)力；
主應(yīng)力分別以σ1，σ2，σ3表示，按數(shù)值排序?yàn)椋害?≥σ2≥σ3。在ABAQUS中分別對(duì)應(yīng)為：Max.principal；Mid.principal；min.principal。這三個(gè)量在任何坐標(biāo)下都是不變量。
我們可利于最大應(yīng)力判斷一些情況：比如最大主應(yīng)力（拉應(yīng)力）大于混凝土的抗拉強(qiáng)度，則認(rèn)為混凝土開裂；通過顯示最大主應(yīng)力的法線方向，則可大致表示出裂縫的發(fā)展影響。
應(yīng)變中一些符號(hào)的含義
E——總應(yīng)變
EP——主應(yīng)變
EE——彈性形變
PE——塑性應(yīng)變分量
Eij——應(yīng)變分量
Ie——非彈性應(yīng)變分量
PEEQ——等效塑性應(yīng)變。若該值大于0，則認(rèn)為已經(jīng)屈服
注：在ABAQUS后處理中，盡量不要看Mises，其表示平均應(yīng)力，更適合金屬材料；對(duì)于鋼筋混凝土結(jié)構(gòu)，我們要看其單軸拉伸方向上的應(yīng)力和對(duì)應(yīng)的應(yīng)變。
展開 
一文搞清ABAQUS中真實(shí)應(yīng)力和真實(shí)應(yīng)變
應(yīng)變計(jì)算公式為：
②仿真中的真實(shí)應(yīng)力與真實(shí)應(yīng)變
在使用ABAQUS仿真時(shí)，如果我們的材料屬于塑性材料范疇，分析時(shí)涉及較大變形，在分析時(shí)必須將其應(yīng)力和應(yīng)變定義成真實(shí)應(yīng)力和真實(shí)應(yīng)變，我們就需要將實(shí)驗(yàn)數(shù)據(jù)中得到的名義數(shù)據(jù)轉(zhuǎn)化成真實(shí)數(shù)據(jù)。
總結(jié)來說，如果不使用真實(shí)應(yīng)力與真實(shí)應(yīng)變，我們就不能分析出我們想要的較為準(zhǔn)確結(jié)果，這里考慮的是材料的非線性問題。
ABAQUS變量解讀：教你讀懂應(yīng)力/應(yīng)變/損傷
結(jié)語
ABAQUS中的變量雖多，但歸根結(jié)底都是為描述材料行為和判斷結(jié)構(gòu)狀態(tài)服務(wù)的。理解每個(gè)變量的物理意義，結(jié)合實(shí)際分析目標(biāo)選擇合適的輸出，才能讓仿真結(jié)果真正“說話”。
如果你對(duì)某個(gè)變量還有疑問，或者想了解更深入的應(yīng)用場(chǎng)景，歡迎在評(píng)論區(qū)留言！
ABAQUS批量提交Job與Python讀取ODB結(jié)果應(yīng)力應(yīng)變數(shù)據(jù)
# coding: utf-8#微信公眾號(hào)：ABAQUS二次開發(fā)#作者：阿信老師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'] #需要輸出的場(chǎng)變量sub_f_S1 = f_S1.getSubset(region=region) S_Value = sub_f_S1.values#讀取所有單元的應(yīng)力值for i in range(len(S_Value)): Stress =S_Value[i].data print Stressodb.close()
若有緣，江湖再會(huì)。
阿信
2022.3.17
文章來源：ABAQUS二次開發(fā)
展開 abaqus中平面應(yīng)力應(yīng)變厚度對(duì)切削力的影響 ￥5
在鋁合金的二位正交切削仿真中，不同的平面應(yīng)力應(yīng)變厚度的對(duì)切削力的影響結(jié)果
以上為設(shè)定值為1的情況
