應力提取的案例
切削殘余應力輸出提取
請問有人可以做切削過后的殘余應力提取嗎？就是殘余應力沿深度方向的分布。
在abaqus中采用python提取結果-平均應力 ￥5
Python腳本操作ODB文件進行計算結果數據的提取。
因對結果分析需求，需提取某單元集的每一分析步（包含每一子步）的平均應力，目前網上雖有很多代碼是關于單元集的平均應力的提取，但并未有針對每一子步都需要提取結果的代碼。故針對此需要編寫了python代碼。
ABAQUS怎么提取扣件支反力和結構層壓應力
ABAQUS怎么提取扣件支反力？軌道結構各結構層壓應力是提取S22還是最大主拉應力
ANSYS如何提取某一節點的應力時程 ￥100
在對結構進行時程分析后，我們經常提取的是全時程最大von Mises stress。
那么如何提取某一個節點的von Mises stress呢？
首先明確ANSYS的節點附加在單元上，可以通過選擇單元上節點的方法提取節點應力。
1 確定節點所在單元，顯示節點編號。
例單元號8560，節點號8678。
2 進入TimeHist Postpro, 定義變量。
3變量顯示。
付費內容為相關命令流。
Ansys Workbench諧響應掃頻結果后處理，提取Von Mises掃頻曲線和應力幅值 ￥10
其作用是，對于多組件的產品掃頻仿真時，我們僅僅提取我們所關心的零部件的應力值和頻率。所以，我們要人為的選定一個selection集合。
在諧響應分析中插入后處理命令，并且確保其中cmsel命令后面的集合名稱，與named selections中的命名“body1”保持一致。
完成后處理Command命令修改后即可提交計算。計算完成后在結果文件夾中有txt文檔和Von Mises應力掃頻曲線。
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 nCode DesigenLife的創建與計算原則淺述
焊縫單元厚度推薦為：
亦或定義為：
三、角焊縫計算要點和應力提取
在ANSYS nCode DesignLife中通過高級編輯“Advanced edit”進行焊縫求解中的“EntityDataType”、“WeldResultLocation”、“ WeldEndElements”的配置。
求解引擎屬性中“EntityDataType”控制采用的數據類型，如圖2所示。
① Stress：直接采用應力作為疲勞評估計算。
② ForceMoment：采用節點力和力矩提取疲勞計算的應力。
③ Displacement：采用位移量提取進行疲勞計算的應力。
WeldResultLocation=NodeOnElement/MidElementEdge考慮基于單元節點或者單元邊長的中間點進行評估，如圖3所示。
WeldEndElements用于在分析中考慮焊線端部單元。
圖2
圖3
限于篇幅以下僅針對基于“EntityDataType=ForceMoment”進行節點力和力矩進行應力計算要點的說明，其他方法可以參閱相關技術文檔。
“EntityDataType=ForceMoment”進行節點力和力矩提取要點：
① 垂直于焊接邊的正應力被提取（Weld Top面），平行的正應力和剪應力不進行計算。
② 對于焊趾和焊根單元的應力提取，是毗鄰焊縫單元的邊的中間點的均值應力，方向垂直于該邊。
③ 焊喉單元應力提取的是焊喉單元兩條邊的均值。
④ 參與計算的焊趾單元是與焊喉單元共享節點的且位于焊線前部的那些殼體單元。
展開 應力曲線沿路徑提取 ￥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（整個模型）
無
如果有需要文檔的同學，可以給我留言，備注信息。
ABAQUS實用子程序SPRINC提取主應力
調用的時候直接在UMAT子程序中輸入：CALL SPRINC(S,PS,LSTR,NDI,NSHR)
需要提供給該實用子程序的變量包括：
實用子程序返回UMAT的變量
實例：
CALL SPRINC(S,PS,LSTR,NDI,NSHR) 可以直接寫成
CALL SPRINC(STRESS,PS,1,3,3)或CALL SPRINC(STRESS,PS,1, NDI,NSHR)
這行代碼的意思是將UMAT中的應力分量數組STRESS提供給實用子程序SPRINC計算；1表示LSTR=1，表明提供給實用子程序SPRINC的為應力分量數組；后面的兩個“3”分別表示正應力分量個數為3和切應力分量個數為3，也可以不用改變直接使用NDI和NSHR，因為在三維模型的UMAT中NDI和NSHR已經被賦予了值為3。PS數組儲存返回的計算得到的主應力，這個需要自己在umat程序中定義數組PS.
下面是關于返回主應力順序的驗證：
CALL SPRINC(STRESS,PS,1,3,3)
STATEV(4)=PS(1)
STATEV(5)=PS(2)
STATEV(6)=PS(3)
需要注意以下點：
1.PS數組一般需要在程序中自定義數組變量；
2.S數組一般也需要自定義數組變量，且需要為S數組賦予應力分量或應變分量的值；
3.如果SPRINC或SPRIND返回的主應力為零，檢查輸入的變量或數組是否未定義或者未賦值，其次檢查輸出數組是否未定義。
展開 殘余應力加載下的載荷深度曲線的提取！
由于殘余拉應力，殘余壓應力的存在，式樣會膨脹或者收縮，造成載荷曲線沒有從(0，0)點開始！
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）文件。提供源文件，注釋詳細，可根據需要進行修改。
Matlab“稍作修改”Abaqus-odb結果 | 提取三維裝配體Mises應力
數據及函數文件獲取方式：
Matlab“稍作修改”Abaqus-odb結果 | 提取三維裝配體Mises應力
https://mp.weixin.qq.com/s/l45KSgNxkJ-oCCrDMrUgHA
怎么在abaqus提取xfem下，動態三維應力強度因子？
求解答
