ABAQUS屈服應力的案例
ABAQUS中mises應力云圖顯示的最大值還不到屈服應力值為啥還有PEEQ值
ABAQUS中mises應力云圖顯示的最大值還不到屈服應力值為啥還有PEEQ值，PEEQ云圖有變形值
針對平面應力問題的YLD2000-2D屈服準則及其在ABAQUS中UMAT子程序的實現
Barlat在2003年提出了專門針對平面應力問題的各向異性屈服準則，該屈服準則對于各向異性材料具有很高的精度，得到了廣泛的應用。
YLD2000-2D屈服面示意圖
Yld2000-2d屈服準則由下式給出
其中
矩陣X′和X″的元素分別由柯西應力的下列線性變換獲得
L′和L″的分量由下式求得
積分算法采用徑向返回算法，該方法是穩健和精確的。
當彈性試算超出屈服面時，則需要進行塑性修正
使其滿足
公式9可以通過牛頓法進行迭代求解。
計算的應力應變曲線如下圖所示
B, F. Barlat A , et al. "Plane stress yield function for aluminum alloy sheets—part 1: theory." International Journal of Plasticity 19. 9(2003):1297-1319.
王海波, 萬敏, 閻昱,等. 屈服準則在有限元軟件中實現的正確性驗證[J]. 固體力學學報, 2010, 031(002):173-180.
最后，有需要歡迎通過微信公眾號聯系我們。
微信公眾號：320科技工作室。
展開 ANSYS后處理中的應力與屈服準則！
此時，形式上將主應力的某一綜合值與材料單向拉伸軸向拉壓許用應力比較，這個綜合值就是等效應力——equivalent stress。
03
最大切應力理論
該理論認為，引起材料屈服的主要因素是最大切應力，不論何種狀態，只要最大切應力達到材料單向拉伸屈服時的最大切應力，則認為材料屈服。
04
畸變能理論
該理論認為，彈性體在外力作用下產生變形，荷載做功、彈性體變形儲能，稱之為應變能（分為畸變能和體積的改變能）。引起材料屈服的主要因素是畸變能密度，無論何種狀態，只要畸變能密度達到材料單向拉伸屈服時的畸變能密度，材料就屈服。
對于各向同性材料的屈服準則
01
屈雷斯加屈服準則
當材料中的最大剪應力達到某一臨界值時，材料發生屈服。該臨界值取決于材料在變形條件下的性質，而與應力狀態無關。
展開 ANSYS后處理中的應力與屈服準則
其中，某點的最大拉應力數值，就是其第一主應力數值。
02
最大拉應變理論
該理論認為，引起材料破壞的主要因素，是最大拉應變。無論何種狀態，只要最大拉應變達到材料拉伸斷裂時的最大應變值，則材料斷裂。此時，形式上將主應力的某一綜合值與材料單向拉伸軸向拉壓許用應力比較，這個綜合值就是等效應力——equivalent stress。
03
最大切應力理論
該理論認為，引起材料屈服的主要因素是最大切應力，不論何種狀態，只要最大切應力達到材料單向拉伸屈服時的最大切應力，則認為材料屈服。
04
畸變能理論
該理論認為，彈性體在外力作用下產生變形，荷載做功、彈性體變形儲能，稱之為應變能（分為畸變能和體積的改變能）。引起材料屈服的主要因素是畸變能密度，無論何種狀態，只要畸變能密度達到材料單向拉伸屈服時的畸變能密度，材料就屈服。
對于各向同性材料的屈服準則
01
屈雷斯加屈服準則
當材料中的最大剪應力達到某一臨界值時，材料發生屈服。
展開 
ANSYS后處理中的應力與屈服準則
其中，某點的最大拉應力數值，就是其第一主應力數值。
02
最大拉應變理論
該理論認為，引起材料破壞的主要因素，是最大拉應變。無論何種狀態，只要最大拉應變達到材料拉伸斷裂時的最大應變值，則材料斷裂。此時，形式上將主應力的某一綜合值與材料單向拉伸軸向拉壓許用應力比較，這個綜合值就是等效應力——equivalent stress。
03
最大切應力理論
該理論認為，引起材料屈服的主要因素是最大切應力，不論何種狀態，只要最大切應力達到材料單向拉伸屈服時的最大切應力，則認為材料屈服。
04
畸變能理論
該理論認為，彈性體在外力作用下產生變形，荷載做功、彈性體變形儲能，稱之為應變能（分為畸變能和體積的改變能）。引起材料屈服的主要因素是畸變能密度，無論何種狀態，只要畸變能密度達到材料單向拉伸屈服時的畸變能密度，材料就屈服。
對于各向同性材料的屈服準則
01
屈雷斯加屈服準則
當材料中的最大剪應力達到某一臨界值時，材料發生屈服。
展開 將應力拖到屈服面上程序代碼
Mises屈服的Chaboche運動硬化背應力詳細推導過程 ￥15
鑒于目前國內外很多書上關于Chaboche運動硬化背應力只是給出了一個表達式的具體形式而其具體推導過程并沒有給出，這造成很多小伙伴在科研的時候很疑惑，因此我把詳細的推導過程給了出來，希望能幫助到大家。
Abaqus 采用YLD2000-2D屈服準則的UMAT子程序
Barlat在2003年提出了專門針對平面應力問題的各向異性屈服準則，該屈服準則對于各向異性材料具有很高的精度，得到了廣泛的應用。
YLD2000-2D屈服面示意圖
Yld2000-2d屈服準則由下式給出
其中
矩陣X′和X″的元素分別由柯西應力的下列線性變換獲得
L′和L″的分量由下式求得
積分算法采用徑向返回算法，該方法是穩健和精確的。
當彈性試算超出屈服面時，則需要進行塑性修正
使其滿足
公式9可以通過牛頓法進行迭代求解。
計算的應力應變曲線如下圖所示
B, F. Barlat A , et al. "Plane stress yield function for aluminum alloy sheets—part 1: theory." International Journal of Plasticity 19. 9(2003):1297-1319.
王海波, 萬敏, 閻昱,等. 屈服準則在有限元軟件中實現的正確性驗證[J]. 固體力學學報, 2010, 031(002):173-180.
展開 【abaqus】個人筆記—應力奇異&應力平均&應力集中
【abaqus】個人筆記—應力奇異&應力平均&應力集中
ABAQUS熱應力分析 附ABAQUS中初始地應力的施加下載
軋輥與Cu層的熱傳導系數
下載地址：ABAQUS中初始地應力的施加
Abaqus帶螺紋螺栓接觸應力分析淺析 Abaqus帶螺紋螺栓接觸應力分析淺析
目前的常規做法通常有兩種：1.簡化，用RBE2和beam梁來代替螺栓，這樣不能反映連接螺栓真實應力，圖1為某結構連接螺栓簡化的beam梁應力云圖，沒有接觸應力：
.直接做出來螺栓螺紋采用接觸分析，雖然得出的結果很精確，但這樣前處理工作量大（螺栓和螺紋用六面體網格建模）、計算量大（接觸收斂困難），如圖為某結構帶螺紋螺栓和連接件模型（圖2）和計算得出的結果（圖3）：
圖3 計算結果
那么，有什么好辦法可以不用簡化帶螺紋螺栓，不用直接做出帶螺紋螺栓，又能得到足夠精確的結果？
運用大型通用非線性有限元分析軟件Abaqus，只需要在接觸定義中設置跟實際螺紋形狀有關聯的參數，如牙角、螺距、螺栓小徑等，就可以模擬真實的連接螺栓接觸狀況。既可以得到足夠精確的分析結果，又節省了時間專注進行其他的分析設置。如圖4，為連接螺栓接觸來定義帶螺紋螺栓：
圖4 連接螺栓接觸來定義帶螺紋螺栓
圖5為某結構直徑10MM的帶螺紋的連接螺栓接觸壓力分布云圖：
圖5 某結構直徑10mm帶螺紋的連接螺栓接觸壓力分部云圖
展開 
ABAQUS-真實應力和名義應力轉化
ABAQUS-真實應力和名義應力轉化.doc
針對某袋除塵器整體進行ABAQUS有限元分析，考慮九項載荷工況，分析設備靜應力、熱應力、變形及熱膨脹數值 ￥15
某袋除塵殼體結構選型如下：
箱體板厚5mm
箱體角柱：角鋼L90*56*8
箱體加強筋：角鋼L90*56*6
花板厚6mm
花板下加強筋：橫向為扁鋼80*6，縱向為扁鋼100*6
箱體中間支撐管：鋼管Φ60*5
圖1 袋除塵殼體結構示意圖
2、 建立模型
按照殼體結構示意圖建立幾何模型如圖2所示。
圖2 建立幾何模型
三、約束條件及載荷
立柱底部約束如圖3所示。
圖3 立柱底部邊界約束
載荷：
（1）自重（軟件考慮）；
（2） 頂部載荷：檢修載（按400kg/m2）；
（3） 花板處載荷：濾袋、濾籠、濾袋積灰（積灰厚度按5mm）共3.06t；
（4） 灰斗積灰重：滿灰9.6t；
（5） 保溫載荷：按25kg/m2；
（6） 負壓11000Pa或正壓8000Pa兩種工況分別施加；
（7） 煙道及檢修平臺載荷：上煙道（出氣端）900kg，下煙道（進氣端）
400kg，上中下三層檢修平臺檢修載荷均為400×2.85×3.25=3705kg。
注：此項載荷殼體和鋼支架各占一半。
（8） 灰斗卸灰口載荷（方向按照幾何模型坐標系）：FX=4700N，FY=3500N，FZ=-4700N，MX=3690N.m，MY=4800N.m，MZ=5540N.m。
（9） 頂部牛腿處檢修荷載：單個牛腿處載荷為1t，頂板為260×260，轉化為面壓添加，面壓為1×10×1000/260/260=0.148N/mm2。
下圖4所示為載荷添加圖示：
（a）負壓11000Pa （b）正壓8000Pa （c）花板處載荷
展開 ABAQUS中求解某部分單元的平均應力或平均應變 ￥10
1、參考模型：單向纖維的RVE模型；
2、腳本功能：針對指定的單元集合，在后處理中求解平均應力和平均應變。
3、應用的公式：一階均勻化計算方法。對于 RVE 模型的平均真應力和平均真應變，可通過對 RVE 內每一個單元的真應力 （真應變）取均值獲得。使用一階均勻化計算方法輸出的應力和應變適用于各種邊界條件，但需要對每個單元進行應力（應變）的輸出和計算。
通過Abaqus python腳本批量獲取節點的應力 ￥25
背景
有限單元法計算單元積分點的應力應變值，而對于節點的應力應變值是通過外插得到的，Abaqus中云圖顯示的就是經過插值和平均后的節點的值。通過工具欄的Query-Probe values可以查看單元或節點的應力應變等結果。
對于自動化的后處理場景，通常需要自動批量地獲取單元/節點的結果，通常都需要通過python腳本來實現。通過類似odb.steps['Step-1'].frames[-1].fieldOutputs['S']的場輸出可以比較方便地直接獲得單元的積分點應力，但沒有直接的API可以獲取節點的應力應變等結果。
如果需要獲取部件表面節點應力，可以通過創建路徑+XYData的方式實現，但想要獲得最大節點應力，則該方式不便實現。
2. 通過python腳本獲取節點應力結果
本文通過fieldOutput.getSub()函數獲取所有單元的節點結果，并對每一節點關聯的多個單元的節點值進行平均后得到節點的結果。以下以某個簡單的odb結果進行驗證。
（1）批量獲得節點的mises應力值
（2）批量獲得節點的X方向正應力值
（3）批量獲得節點的最大主應力值
（4）獲取節點的最大mises應力及編號
3. 獲取節點應變等結果
只需將腳本程序中的應力場改為應變成E等即可，此處不再演示。
以下為本文的python腳本代碼（代碼中作了必要的簡單注釋）。
展開 