ansys提取積分點應力的案例
Abaqus中平面應力單元高斯積分點的順序
可以輸出umat接口中的變量coords進行查看
write(*,"(A,I4)") "npt = ", npt
write(*,"(A,3ES16.8)") "coords = ", coords
結果為：
npt = 1
coords = -5.77350269E-01 -5.77350269E-01 1.00000000E-02
npt = 2
coords = 5.77350269E-01 -5.77350269E-01 1.00000000E-02
npt = 3
coords = -5.77350269E-01 5.77350269E-01 1.00000000E-02
npt = 4
coords = 5.77350269E-01 5.77350269E-01 1.00000000E-02
因此Abaqus中平面應力單元高斯積分點的順序為：
展開 單元積分點應力如何外插至節點上 | 數值實現篇
繼上次的推文：有限元計算過程中積分點應力如何外插至節點處？【公式推導篇】，本次分享單元積分點應力外插至節點處的數值實現過程。
數值實現
借助以上理論，我們可以基于matlab平臺編制以下代碼段：
% 將積分點應力外插至單元節點上，這里只列舉了Q4的情況
for i = 1:3
StressElem(e,:,i) = [1+0.5*sqrt(3) -0.5 1-0.5*sqrt(3) -0.5;
-0.5 1+0.5*sqrt(3) -0.5 1-0.5*sqrt(3);
1-0.5*sqrt(3) -0.5 1+0.5*sqrt(3) -0.5;
-0.5 1-0.5*sqrt(3) -0.5 1+0.5*sqrt(3)]*...
[stress(e,1,i);stress(e,2,i);stress(e,3,i);stress(e,4,i)];
end
對標Abaqus
模型材料參數為普通的線彈性材料，單元類型選擇CPS4，網格劃分及邊界條件設置如下：
在結果對標過程中，可以先對比自研程序與Abaqus的節點位移場：
Abaqus位移場結果
自研程序位移場結果
在位移場一致的前提下，我們再來對標應力結果。以常見的mises應力為例：
Abaqus位移應力場結果
自研程序應力場結果
結果是一致的，說明了程序的正確性。
展開 有限元中單元積分點與節點應力相互轉換(CPE4為例)（ABAQUS）
在ABAQUS中,當需要獲取節點上的應力時,可以在后處理中建立路徑或者用查詢功能等獲取.
但是當需要大量的節點上應力數據時,很多人會用Python編程進行大批量的提取應力.但是提取出來的應力為單元積分點上的應力.無法獲取節點上的應力.同時在ABAQUS中的子程序中,也是對積分點上的數據進行操作.
本文基于個人興趣同時想要更加了解有限元背后原理和公式的想法.近日進行了一些初步的探索.希望大家批評指正. 本文基本不涉及原理公式，只在轉換積分點和節點的應力時列出公式。盡可能簡介易懂。
一: 單元類型及節點數目與位移,應變,應力階次的關系
本節內容基于有限元教材及一些網上資料.
(1)有限元求解的思路是:
一: 建立單元節點力與節點位移關系式.
二: 將彈性體上的外載荷等效移置到節點上.
三: 在節點上建立力的平衡方程,求得節點位移.
四: 通過彈性力學基本方程,可求得單元的應力和應變.
(2) 四節點矩形單元
以四節點矩形單元為例，在此只表達有限元教材中的結論，具體公式可參考有限元教材。
（3）ABAQUS中的CPE4單元
CPE4: A 4-node bilinear plane strain quadrilateral.
該單元有四個節點，同時有四個積分點。
對于每個應力分量（注意：在此只看一個應力分量），單元內任一點（x,y）的應力表達式為：
stress=a*x*y+b*x+c*y+d （1）
該表達式有四個未知量：a,b,c,d。
若知道四個積分點的應力分量。將每個積分點帶入上式，則會形成包含四個方程的線性方程組。
展開 有限元計算過程中積分點應力如何外插至節點處？【公式推導篇】
注：由于技術鄰排版風格有限，故部分內容顯示不全，感興趣的小伙伴可點擊原文進行閱覽：
有限元計算過程中積分點應力如何外插至節點處？【公式推導篇】
https://mp.weixin.qq.com/s/47byQ3b3e5UpbUp7Krs2mQ
本次分享的是：有限元計算過程中，單元積分點應力如何外推至節點？
有關積分點與節點的概念可點擊跳轉閱讀歷史推文：有限元基本概念-【節點和積分點】，現科普一下Q4單元、Q8單元、Q9單元的形函數和高斯積分方案。
Q4單元
Q8/9單元
應力外插
核心理念：坐標系的轉換。
假設是母單元的自然坐標系，是由高斯積分點控制的坐標系（術語可能不專業），假設高斯積分方案為。坐標系轉換關系：
單元內任一點的應力，由4個高斯積分點應力進行插值時，可表示為
其中，是基于高斯積分點的形函數，第一個積分點的坐標在母單元坐標系下為（-1，-1），根據上述的坐標系轉換的方式，在高斯積分點的坐標系下，第一個單元節點在高斯積分點坐標系下坐標為，將此坐標值代入第一個形函數，得，相同的道理，可推導至四個節點在4個形函數下的外插矩陣：
對于Q8、Q9單元，依然可采用高斯積分方案（減縮積分）。
展開 
ANSYS如何提取某一節點的應力時程 ￥100
在對結構進行時程分析后，我們經常提取的是全時程最大von Mises stress。
那么如何提取某一個節點的von Mises stress呢？
首先明確ANSYS的節點附加在單元上，可以通過選擇單元上節點的方法提取節點應力。
1 確定節點所在單元，顯示節點編號。
例單元號8560，節點號8678。
2 進入TimeHist Postpro, 定義變量。
3變量顯示。
付費內容為相關命令流。
ansys workbench rst 文件應力、位移和坐標結果提取
采用python語言提取rst 文件結果提取
ANSYS瞬態分析全時程結構響應最大值的提取方法（變形、應力、應變、能量） ￥100
</p><p>同樣的方法，可以提取全時程最大的位移、應力、應變、能量等結果。</p><p><img src="https://img.jishulink.com/upload/202302/623025b5c0d646b9973cd2adc6c6037f.png" alt="1.png"></p><p>收費內容為相關命令流。</p>
Ansys Workbench諧響應掃頻結果后處理，提取Von Mises掃頻曲線和應力幅值 ￥10
問題：
Ansys workbench進行諧響應仿真計算的后處理結果中，提供了單一頻率下的Von Mises應力查看功能和應力頻響曲線功能，但是應力頻響曲線的應力列表中沒有Von Mises應力查看項。因為Von Mises應力太常用，所以這就給我們在整個掃頻范圍內，定位Von Mises應力的最大頻率和應力值帶來一定的困難。如下所示。
需求：
希望后處理結果中可以在應力響應曲線中，有一項Von Mises應力選項。實現每個掃頻點的最大Von Mises應力和掃頻頻率的曲線圖顯示，從而一眼就看出產品在整個掃頻范圍內，哪個頻率下結構的等效應力最大。而后再通過應力云圖查看這個頻率下的Von Mises應力。
解決方法：
利用APDL命令實現。簡要流程為：首先，讀取每一個掃頻點的最大Von Mises應力值。記下應力值、頻率值和最大節點號。再統計記錄的所有掃頻點的Von Mises應力值，提取整個掃頻過程中最大應力值及其頻率。并將結果寫出到txt文件。進一步提取這個最大Von Mises應力點對應的整個掃頻范圍內的Von Mises應力曲線。
這個樣就可以在txt文檔中直接看到所有掃頻點下，結構的等效應力幅值；以及全頻段中最大Von Mises應力所在節點的等效應力掃頻曲線圖。
效果展示如下：
在結果文件夾中，會生成一個txt結果文件和一張Von Mises應力曲線圖。如此我們可以直觀注意到，在當前掃頻范圍內，結構在78.95Hz時應力最大約為17.552Mpa。
結果后處理問題示例：
Ansys workbench進可以查看某個頻率下的 Von Mises應力幅值
Ansys workbench進掃頻應力響應曲線中，應力選項卻沒有Von Mises應力選型，只能按三個方向來分別查看。
展開 Ansys 查看高斯點上的應力
許多時候我們需要在ANSYS中查看高斯點上的應或者和應變，然而我們看到的節點上的應力或者應變通常是由高斯點上的應力或者應變外插而來，這時候我們就需要用到ERESX這個命令了。
ERESX命令使用格式：ERESX，Key(GUI: Main>solution > Load Step Opts > Output Ctrls > Integration Pt或Main Menu > Preprocessor > Loads > Load
Step Opts > Output
Ctrls > Integration Pt)
Key為外插法控制鍵，有DEFA,YES和NO三個選項，分別對應著三種情況：
DEFA(默認設置)：除了具有塑性、蠕變或膨脹等非線性特性的單元意外，將積分點的結果進行外插擴展到所有單元的節點上。
YES: 將積分點的結果進行外插擴展到所有單元的節點上，僅將線性結果數據通過外插法擴展到這些具有塑性、蠕變或膨脹非線性特性的單元上。
NO: 將積分點上的結果復制(不是外插)到所有單元的節點上。
顯然，當我們不確定ANSYS是如何外推的，想直接查看高斯點上的應力、應變或其它結果的時候，我們就可以直接使用ERESX,no這個命令來查看了。
注意：對于非線性的數據ANSYS總是采用復制的方式擴展到節點上，而不是外推法,當 然，你也可以用ERESX,yes來采用外推法；這個命令同樣可以在prep7中使用；
轉載來源于
http://blog.sina.com.cn/s/blog_934e096a0102wkyb.html
展開 