abaqus應力提取的案例
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平均應力和應變提取 ￥80
利用python腳本對ODB文件中單元集里所有積分點的應力及應變進行自動提取并計算平均值
能夠得到每一幀的應力和應變平均值，并保存到CSV文件中
所得到的應力包括S11,S22,S33,S12,S13,S23以及Mises七個應力平均值，以及E11,E22,E33,E12,E13,E23六個應變平均值
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提取單元平均應力/應變 ￥10
利用python讀取odb文件（可一次讀取多個odb）生成csv（excel）文件。提供源文件，注釋詳細，可根據需要進行修改。
在abaqus中采用python提取結果-平均應力 ￥5
Python腳本操作ODB文件進行計算結果數據的提取。
因對結果分析需求，需提取某單元集的每一分析步（包含每一子步）的平均應力，目前網上雖有很多代碼是關于單元集的平均應力的提取，但并未有針對每一子步都需要提取結果的代碼。故針對此需要編寫了python代碼。
Matlab“稍作修改”Abaqus-odb結果 | 提取三維裝配體Mises應力
數據及函數文件獲取方式：
Matlab“稍作修改”Abaqus-odb結果 | 提取三維裝配體Mises應力
https://mp.weixin.qq.com/s/l45KSgNxkJ-oCCrDMrUgHA
怎么在abaqus提取xfem下，動態三維應力強度因子？
求解答
ABAQUS怎么提取扣件支反力和結構層壓應力
ABAQUS怎么提取扣件支反力？軌道結構各結構層壓應力是提取S22還是最大主拉應力
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.
展開 應力曲線沿路徑提取 ￥8
通過鉆柱長度和角度繪制鉆柱簡化
序號
角度（°）
長度（m）
1
98.22
9.49
2
98.54
9.47
3
99.58
9.47
4
100.3
9.40
5
100.33
3.00
根據鉆柱簡化模型基本參數建立模型
對鉆柱添加邊界條件和載荷約束如圖2所示
載荷及邊界條件
鉆桿頭部
中間鉆桿
鉆桿尾部
約束
X軸自由度
Y軸自由度
Z軸自由度
Y軸自由度
X軸自由度
Y軸自由度
Z軸自由度
繞X軸轉動自由度
載荷
鉆壓5t
扭矩1000Nm
重力9.8m/s2（整個模型）
無
如果有需要文檔的同學，可以給我留言，備注信息。
切削殘余應力輸出提取
請問有人可以做切削過后的殘余應力提取嗎？就是殘余應力沿深度方向的分布。
ANSYS如何提取某一節點的應力時程 ￥100
在對結構進行時程分析后，我們經常提取的是全時程最大von Mises stress。
那么如何提取某一個節點的von Mises stress呢？
首先明確ANSYS的節點附加在單元上，可以通過選擇單元上節點的方法提取節點應力。
1 確定節點所在單元，顯示節點編號。
例單元號8560，節點號8678。
2 進入TimeHist Postpro, 定義變量。
3變量顯示。
付費內容為相關命令流。
殘余應力加載下的載荷深度曲線的提取！
由于殘余拉應力，殘余壓應力的存在，式樣會膨脹或者收縮，造成載荷曲線沒有從(0，0)點開始！
Python提取場輸出結果計算平均應力應變
###########################################################################
寫在前面，最近回看了一些自己之前編寫的和這個py代碼，發現自己寫的可能和公式表達的有偏差，在我的測試過程中，提取的是單元積分點處的應力應變值，由于只有一個積分點（這個積分點在質心，與centroid提取得到的結果一樣），因此自然的把這個值當成了單元的平均應力或平均應變，然后進行計算，但是我現在覺得這并不是獲取單元平均應力的方式，也就是代碼并沒有實現所謂的提取平均應力應變的功能，希望有大神可以指點迷津。
如果是有多個積分點的話，是不是應該對每個積分點權重進行積分，加權平均這樣得到單元的平均應力，然后乘單元體積，將所有單元的值求和再除模型的總體積，就得到整個RVE模型的平均應力。
###########################################################################
最近寫了一個簡單的python讀取abaqus結果中的場輸出數據，想通過均勻化計算方法來計算所定義集合的平均應力應變曲線，之前是手動提取了各個數據導出，然后用excel、matlab處理，但是很慢，而且很費勁，于是就想著用Python來處理結果。
有需要的同學可以下載附件文件，打開abaqus，file→run script，選擇腳本文件即可運行。
average.zip
均勻化計算方法：
參考文獻：馬思鳴. 精沖用碳鋼微觀組織對宏觀力學性能及精沖性能影響研究[D]. 上海：上海交通大學，2017.
我在網上找了挺久挺多的關于Python提取場輸出結果的，
主要對以下幾篇帖子進行了參考：
http://forum.simwe.com/forum.php?
展開 Ansys Workbench諧響應掃頻結果后處理，提取Von Mises掃頻曲線和應力幅值 ￥10
其作用是，對于多組件的產品掃頻仿真時，我們僅僅提取我們所關心的零部件的應力值和頻率。所以，我們要人為的選定一個selection集合。
在諧響應分析中插入后處理命令，并且確保其中cmsel命令后面的集合名稱，與named selections中的命名“body1”保持一致。
完成后處理Command命令修改后即可提交計算。計算完成后在結果文件夾中有txt文檔和Von Mises應力掃頻曲線。
