應力結果的案例
axial stress+ bending stress=long stress，axial stress（軸向應力）和normal Stress SNMAX（最大法向應力）的結果不相等 ￥20
繼上一篇關于梁單元結構建模，optistruct求解后，hyperview查看應力，只有axial stress和long stress,沒有Von mises stress 等應力結果的原因后，但篇中沒有關于對axial stress和long stress的解釋，其實在查看結果時發現，axial stress（軸向應力）和normal Stress SNMAX（最大法向應力）的結果不相等，但是在幫助文件了翻找了好久也沒有關于axial stress等應力的相關解釋，卡在這好幾天。
突然想到Nastran的關鍵字與Optistruct的一樣，能在Optistruct運行也可以在Nastran中運行求解。
展開 hyperview二次開發-批量輸出應力結果
一、hyperview二次開發-批量輸出應力結果（此腳本源代碼50元）
分析結果文件：NASTRAN的.op2文件
功能：可輸出多個comp在不同工況下的最大應力值及安全系數，保存在.csv的表格文件中，思路如下
二、需要定制功能或其他腳本（價格可議）
hypermesh前處理、hyperview云圖動畫、hypergrash曲線，自動輸出分析報告等均可滿足您的需求，
可聯系Q1700403934
https://www.yqgqt.org.cn/college/video/c13045
workbench怎么展示多方向的結果應力圖
問題：workbench在出圖的時候，不管是模型示意圖，載荷施加圖還是應力結果圖，有時候一張圖顯示不全，效果不是很直觀，就需要多個方向的去截圖顯示。
如圖所示，點擊紅圈部分，會有one，horizontal，vertical，four四種顯示模式，可以是默認的一個界面，橫向兩個界面，豎向兩個界面，四個界面。這樣可以多方向的展示你的結果，每個圖框可以調整方向，還可以單獨做個截面等。
單元應力求解方法簡述、結果應力種類選項值的含義以及計算舉例 ￥1
后處理求解結果應力種類選項值的含義與計算舉例：
unveraged
areraged
nodal difference
nodal fraction
elmemntal difference
elemental mean
elemental fraction
ABAQUS Fortran基于提供的代碼 P53，完善 8 節點單元程序，并增加對應力結果的處理。
計算節點應力，給出并實現至少一種應力處理方案，提供處理 前后的應力結果（可用表格和云圖表示），可與其它軟件對比； 3.提交總結報告（包括方法/方案描述、帶詳細注釋的代碼、程序框圖、算例描述、結果比較分析等）、可編譯源代碼、可執行文件、 數據文件、結果文件
program p53
!-----------------------------------------------------------------------------
! program 5.3 plane strain of an elastic solid using uniform
! 8-node quadrilateral elements numbered in the x direction
!-----------------------------------------------------------------------------
use new_library ; use geometry_lib ; implicit none
integer::nels,nxe,neq,nband,nn,nr,nip,nodof=2,nod=8,nst=3,ndof,loaded_nodes,&
i,k,iel,ndim=2
real::aa,bb,e,v,det ; character(len=15) :: element = 'quadrilateral'
!
展開 WORKBENCH 后處理中各種應力結果說明
在ANSYS WORKBENCH的結構分析后處理中，我們經常會關注應力。在選擇一個對象并查看某種應力后，會在其細節視圖中出現一個積分點結果的顯示選項，說明要觀察應力的哪種結果，如下圖。
這里面有七種查看應力的方式。那么這些方式分別是什么含義呢？
（1）unveraged---------沒有平均化的應力。此時在單元內部，基于積分點的應力外推算該單元幾個節點的應力。此時如果一個節點周圍毗鄰幾個單元，那么這幾個單元在同一點處會有不同的應力值。
（2）areraged--------節點的平均化應力。在對所有單元進行計算，得到其節點的應力后，此時對于共享節點，對該點的幾個應力進行平均，得到該點的應力。（默認）
（3）nodal difference---節點應力差的最大值。對于共享節點，還沒有進行應力平均時，它有幾個應力，對這幾個應力排序，得到最大值，最小值；用最大值減去最小值，得到的值稱為nodal difference.
(4) nodal fraction------對于一個共享節點，用（3）除以（2），得到一個比率，就是nodal fraction.
（5）elemental difference-----在一個單元內部操作。一個單元的節點應力的最大值，最小值，用最大值減去最小值，得到一個值，稱為elemental difference.
(6) elemental mean-----在一個單元內部操作。在節點應力平均后，對于單元的所有的節點應力，再一次平均，得到單元內部的elemental mean.
(7) elemental fraction------在一個單元內部，用（5）除以（6），得到elemental fraction.
展開 Python提取場輸出結果計算平均應力應變
###########################################################################
寫在前面，最近回看了一些自己之前編寫的和這個py代碼，發現自己寫的可能和公式表達的有偏差，在我的測試過程中，提取的是單元積分點處的應力應變值，由于只有一個積分點（這個積分點在質心，與centroid提取得到的結果一樣），因此自然的把這個值當成了單元的平均應力或平均應變，然后進行計算，但是我現在覺得這并不是獲取單元平均應力的方式，也就是代碼并沒有實現所謂的提取平均應力應變的功能，希望有大神可以指點迷津。
如果是有多個積分點的話，是不是應該對每個積分點權重進行積分，加權平均這樣得到單元的平均應力，然后乘單元體積，將所有單元的值求和再除模型的總體積，就得到整個RVE模型的平均應力。
###########################################################################
最近寫了一個簡單的python讀取abaqus結果中的場輸出數據，想通過均勻化計算方法來計算所定義集合的平均應力應變曲線，之前是手動提取了各個數據導出，然后用excel、matlab處理，但是很慢，而且很費勁，于是就想著用Python來處理結果。
有需要的同學可以下載附件文件，打開abaqus，file→run script，選擇腳本文件即可運行。
average.zip
均勻化計算方法：
參考文獻：馬思鳴. 精沖用碳鋼微觀組織對宏觀力學性能及精沖性能影響研究[D]. 上海：上海交通大學，2017.
我在網上找了挺久挺多的關于Python提取場輸出結果的，
主要對以下幾篇帖子進行了參考：
http://forum.simwe.com/forum.php?
展開 有限元分析中的一些問題--應力集中結果的可信性
如果結構中沒有應力集中，答案是肯定的。
如果結構中存在應力集中，則結果未必會收斂。
為了說明這一點，我們選取了一個平面應力問題。它是一個角支座，其圖形及尺寸如下。在角支座上鉆了兩個孔，現在我們固定左上邊的孔，而在右下方孔的第四象限半圓上施加壓力。并通過不斷的加密網格來考慮計算結果的可信性。
生成的有限元模型如下
固定左上邊的孔，并對右下方孔施加右下方向的壓力，當單元尺寸取5mm時候，應力云圖如下
可見，此時最大應力發生在拐角處，是34.383MPa.
單元尺寸全局細分到3mm，結果是
最大應力是44.44MPa.
單元尺寸全局細分到1mm，結果是
最大應力是74.004MPa.
單元尺寸全局細分到0.4mm，結果是
最大應力是112.873MPa.
可見，結果并沒有收斂的趨勢。
如果我們進一步細分網格，會發現數據無限增大，不會收斂。
實際上，理論證明，在該拐角處如果是直角，而沒有倒圓角的話，應力集中系數會趨向無窮大，所以在實踐設計中絕對禁止出現這種直角。
這也意味著，如果我們在有限元分析前進行模型簡化時，絕不可輕易將一些倒角隨便刪除，否則會出現奇怪的結果。
內容轉自宋博士的博客
展開 ANSYS Workbench中的7種應力結果如何理解？
在ANSYS Workbench的結構分析后處理中，我們經常會關注應力。在選擇一個對象并查看某種應力后，會在其細節視圖中出現一個積分點結果的顯示選項，說明要觀察應力的哪種結果，如下圖。
這里面有七種查看應力的方式。那么這些方式分別是什么含義呢？
由于應力是我們做結構力學分析時最為關注的對象，因此弄清楚ANSYS所給的應力究竟是什么意思也就變得非常重要。這七種應力的含義及相互關系如下圖。
從上圖中可以看到，在計算出積分點的應力以后，其它應力都是在其基礎上推算出來的。下面說明每一個選項的推算過程。
（1）unveraged---------沒有平均化的應力。此時在單元內部，基于積分點的應力根據形函數推算該單元幾個節點的應力。因為它是在積分點應力的基礎上做的第一次運算，所以相對準確。此時如果一個節點周圍毗鄰幾個單元，那么這幾個單元在同一點處會有不同的應力值。
（2）areraged--------節點的平均化應力。在對所有單元進行計算，得到其節點的應力后，此時對于共享節點，對該點的幾個應力進行平均，得到該點的應力。
（3）nodal difference------節點應力差的最大值。對于共享節點，還沒有進行應力平均時，它有幾個應力，對這幾個應力排序，得到最大值，最小值；用最大值減去最小值，得到的值稱為nodal difference.
（4） nodal fraction------對于一個共享節點，用（3）除以（2），得到一個比率，就是nodal fraction.
所以，（2）（3）（4）都是對于共享節點，在不同的單元間進行計算的。
（5）elmemntal difference-----在一個單元內部操作。
展開 四面體網格，六面體網格，低高階單元，對比研究
01 幾何模型
02 部分網格展示
04 用低階六面體單元進行仿真計算
某兩點的位移隨節點數的變化趨勢：
某應力梯度較小位置的應力隨節點數的變化趨勢：
某應力梯度較大位置的應力隨節點數的變化趨勢：
05 用高階六面體單元進行仿真計算
某兩點的位移隨節點數的變化趨勢：
某應力梯度較小位置的應力隨節點數的變化趨勢：
某應力梯度較大位置的應力隨節點數的變化趨勢：
06 六面體單元的相關結論
01 位移結果可靠，節點數和單元階數的影響較小；
02 應力梯度較小位置的應力結果可靠，節點數和單元階數的影響較小；
03 應力梯度較大位置的應力結果不可靠，節點數和單元階數的影響較大；
07 四面體單元仿真計算與相關結論
01 高階四面體單元的位移結果可靠，節點數的影響較小；
02 低階四面體單元的位移結果不可靠，建議不要使用；
03 高階單元在應力梯度較小位置的應力結果可靠，節點數的影響較小；
04 低階單元在應力梯度較小位置的應力結果不可靠，建議不要使用；
05 應力梯度較大位置的應力結果不可靠，節點數和單元階數的影響較大；
08 總結論
01 在結構有限元分析中，建議不要使用低階四面體單元；
02 對于位移結果來說，六面體單元，高階四面體單元的求解都是可靠的，并且節點數影響較小。
03 在應力梯度較小位置，六面體單元，高階四面體單元的求解都是可靠的，并且節點數影響較小。
展開 模態瞬態疲勞分析操作流程（nastran、Ncode、Femfat）
1.模態應力結果保存在op2結果文件中，結果中包含所有求解的模態階次的應力結果。
2.模態應力結果有幾個頻率，就需要設置幾個LC工況。
3.模態位移時間歷程結果保存在pch文件中，femfat無法直接識別，可以通過Ncode將其轉化為femfat可以識別的文件，如Ncode dac格式文件。
Ncode和femfat都可以進行模態瞬態疲勞分析。有限元求解設置是一致的，包括模態應力和模態位移時間歷程結果的生成。Ncode可以直接識別nastran的模態位移結果pch文件，而femfat無法直接識別pch格式的結果文件，因此需要通過中間轉換來完成，比如hypergraph、Ncode等工具或者自己寫個轉換工具。
模態瞬態疲勞分析和多通道靜態疲勞分析的設置基本是一致的。
（文章來源于CAE數值優化輕量化 ，作者方永利）
展開 
Abaqus梁單元基礎知識 附ABAQUS基礎入門與案例精通下載
以下是梁單元的命名規律：
由于空間梁單元除了拉壓、彎曲自由度外，還具有扭轉（翹曲）自由度，所以一般相同邊界載荷條件下，平面梁單元與空間梁單元計算結果會有一些差異，因此，在選擇梁單元時要根據實際情況選擇。
梁單元按節點數量分為兩類：2節點梁單元、3節點梁單元
具有不同積分點的梁單元分類如下：對于單個單元來說，積分點數量越多，單個單元具有更好的柔度，越適合模擬大彎曲變形的結構，如海底光纜。
本文以工字梁作為建模單元：在定義工字梁截面屬性時，I 的作用如下：定義單元橫截面軸在截面內與截面底部的距離。
I=0.2
I=0.6
以下為部分工字梁單元輸出結果：Abaqus梁單元計算結果具有豐富的計算結果（幾十種結果類型），能夠滿足科研、常規工程的計算需求。
單元變形結果（打開beam截面渲染）
單元Mises應力結果（打開beam截面渲染）
本例用工字梁單元建模，下面是工字梁的截面積分點分布情況：
baqus能夠基于beam 截面渲染，給出梁單元截面內的應力分布情況：
單元截面Mises應力結果（必須打開beam截面渲染）
單元截面主應力結果（必須打開beam截面渲染）
單元截面彎曲應力結果（必須打開beam截面渲染）
單元截面翹曲應力結果（必須打開beam截面渲染）
在考慮溫度的情況下，工字梁截面溫度（場變量）積分點的分布情況如下：
應用梁單元做熱分析計算時，應該可以參照截面應力輸出的方法輸出截面溫度分布。
下載地址：ABAQUS基礎入門與案例精通
展開 Ansys Workbench ACT插件，由窗口選中體單元，提取體積和表面積，計算幾何特征尺寸 ￥20
Ansys Workbench ACT插件，由窗口選中體單元，提取體積和表面積，計算幾何特征尺寸
問題：
在FKM關于結構疲勞評估計算方法中指出：零部件特征尺寸，影響疲勞結果評估。原因是材料的應力壽命曲線是由標準試樣進行試驗測試獲得的。當零部件的特征尺寸與測試樣件不一致時，需要考慮零部件的特征尺寸這一因素。（一般而言，當零部件的尺寸大于材料標準測試樣件時，零部件的表面或內部缺陷發生的概率會增加，從而導致零部件尺寸越大，疲勞壽命越低）
對與規則幾何形狀的零部件，有相應的經典公式提供特征尺寸的計算；例如圓形細長桿的特征尺寸是直徑；薄板零部件的特征尺寸是板厚等；但是實際工作中的零部件幾何形狀千差萬別，沒有統一的經典公式可以提供特征尺寸的計算；在FKM手冊中給出了一個通用公式，用于估計零部件疲勞危險區域的局部特征尺寸；
FKM關于循環載荷的疲勞評估中，提及可以使用循環載荷下的有限元應力結果進行疲勞損傷估計。此時，除了需要由應力結果估計危險疲勞區域，提取危險點的應力結果外，還需要給出危險疲勞區域的特征尺寸。在Ansys Workbench中，用戶可以方便的查看應力結果云圖，從而大體評估出危險疲勞區域。并且用戶可以通過選取高應力區域的單元體，再通過特征尺寸一般計算公式，來估計高應力區域的特征尺寸，進行進行合理的FKM疲勞評估。
但是，Ansys Workbench中，當用戶選中了某個/某些體單元后，在選擇信息欄中并不能直接給出單元體積和表面的有效信息輸出。并且通過查詢資料，即使在APDL經典界面中對與體單元也是僅僅只能輸出體積（沒有體單元表面的輸出）；并且對與FKM特征尺寸的一般計算公式中，關于表面積A，也并不是指每個體單元所有面的表面積的總和。
展開 有限元分析的一些基本考慮---位移解應變解和應力解
即使將一個公共節點的多個應力進行平均，以代表該節點的應力值，該平均過程稱為“平滑”。
總之，求解節點應力的步驟是：
（1）根據總體方程，得到節點的位移解。
（2）根據幾何方程，得到單元高斯點的應變解。
（3）根據物理方程，得到單元高斯點的應力解。
（4）在某一個單元內，基于形函數，將高斯點的應力外推到該單元的所有節點。
（5）對于某一個公共節點，將該節點關聯的所有單元所推出的該節點的應力解進行平均，最終得到該節點的應力解。
在ANSYS WORKBENCH的后處理中，如果我們加入了一個應力對象，我們可以看到其細節視圖中有下列選項---積分點結果選項，如下圖
那么這里面的7項是什么含義呢？
下面闡釋這七項的意義。
(1)unaveraged:顯示沒有進行平均的應力結果。
(2)averaged：顯示平均后的應力結果。
(3) nodal difference：對于公共節點，計算其相連各單元計算得到的非平均應力的差的最大值。
(4)nodal fraction：計算公共節點的nodal difference與節點平均值的比值。
(5)elemental difference：對于一個單元上的所有節點，計算其非平均結果的最大差值。
(6)elemental fraction：計算element difference與單元平均值的比值。
(7)elemental mean:根據平均化的應力結果來計算單元的平均值。
這樣，我們在瀏覽應力結果時，應根據需要來選擇我們需要查看的對象。
展開 【OptiStruct要領】掃頻/定頻疲勞以及隨機振動疲勞
需要注意的是在掃頻/定頻疲勞中，一個FATEVNT中僅支持單個FATLOAD，不支持應力疊加。
FATPARM
同樣加上SWEEP 關鍵字，NF/DF 指定參與掃頻的頻率個數/增量；如果當前取到的頻率點沒有對應的應力結果，則用前后的頻率結果進行插值得到其應力。如果定義了NF則忽略DF。STSUBID指定靜力工況，用于引入平均應力。
2）損傷/壽命計算
a、定頻疲勞的損傷計為損傷量和總時長的乘積。
's o 損傷量可以由該頻率下的應力結果和材料的 SN 曲線確定在定頻疲勞中 FATSEQ 直接指定T(總時長)
b、掃頻疲勞的損傷則是計算頻率段內各采樣頻率上的損傷（=單循環損傷*循環次數），然后再疊加。
從上面的定義可知，重點在于確定循環次數以及單次循環的損傷量。
單次循環的損傷量可以由該頻率下的應力結果和材料的SN曲線確定。在掃頻中，按照NF/DF對直接求解得到的頻響函數應力結果進行采樣得到指定頻率下的應力結果。應用該應力結果計算單次循環產生的損傷。如果存在平均應力，可考慮平均應力修正。
通過兩種不同掃頻方式計算相對應的循環次數。
SRUNIT= HZPS 以Hz/s 的速度進行掃頻
f是某個需要計算的頻率（NF/DF指定的）
SRUNIT=OCTPRM 以oct/m 的速度進行掃頻
得到這兩個參數后，就可以計算整個掃頻/定頻中的疲勞損傷/壽命。
展開 