abaqus應力畸變的案例
解決橡膠大變形網格畸變，ABAQUS有絕招（上）
一直以來，這樣一個問題始終在困擾著CAE工程師及有限元仿真族們：在橡膠材料的大變形問題中，往往由于過于劇烈的材料變形，導致開始時劃分的網格畸變嚴重，不再適用于新的幾何形狀，局部應力集中，最終導致計算結果不收斂，分析無法繼續進行。如原本球形的橡膠材料最終變為立方體，又或者截面形貌發生劇烈變形等等，隨著變形的加劇，網格的不再適應對于大變形問題是一直是困擾工程師們的頭疼問題。
不需要的復雜拉格朗日-歐拉自適應網格技術，也不需要效果欠佳的自適應ALE劃分，今天，我們就通過一個實例來看看，如何運用ABAQUS的*map solution功能，聯合強大的前處理軟件ANSA，來快速、輕松、高效解決大變形中網格畸變的問題。
如圖所示，在ABAQUS/CAE中建立幾何模型并劃分網格定義接觸和邊界條件，在橡膠試樣頂部施加-Y方向50mm的位移，提交計算后可以看到，因為大變形導致的網格畸變使得計算無法完成。
可以看到隨著橡膠試樣被擠入左右的耳腔里，網格發生嚴重畸變，長寬比急劇增大，導致計算到0.89秒時，無法收斂而終止。我們可以將此分析分為以下兩步來完成：首先，對橡膠試樣下壓30mm（此時網格剛剛開始發生畸變，注意：不要用計算不收斂時的最后一步做為重構網格的初始模型，因為不收斂之前網格已經發生的劇烈的變形，局部應力集中，計算已經不再準確），隨后通過下壓30mm的odb結果文件導出最后一步中的part模型，在前處理軟件ANSA中進行網格重構。
重構后的網格較為規整，更加利于計算的收斂。接著將重構后的網格導入ABAQUS/CAE中進行前處理的重新定義，并寫入INP文件。
展開 解決橡膠大變形網格畸變，ABAQUS有絕招（下）
在《解決橡膠大變形網格畸變，ABAQUS有絕招（上）》中介紹到，運用ABAQUS中的map solution功能和強大的前處理軟件ANSA以及ABAQUS/CAE聯合操作，可以實現大變形中由于網格畸變而導致的計算不收斂問題，本文將進一步介紹，map solution的實現方法和橡膠大變形實例的操作。
在上篇文章中講到，試樣被下壓30mm后進行重構網格，之后導入ABAQUS/CAE中進行前處理的重新設置，最后將map solution命令寫入INP文件中，如下圖所示。
提取試樣頂部中間節點的載荷曲線可以看到，由于map solution的將節點應力映射到新的網格上，兩個模型的載荷曲線基本達到首尾相連的狀態，少量的誤差也是在可接受范圍之內的，這樣就通過映射網格節點數據的方法，對畸變的網格進行重新劃分，使得模型能夠成功的完成整個下壓50mm的大變形過程。
變形后應力云圖如下圖所示，可以看到，由于重構了網格，橡膠的填充效果比較理想，并沒有出現大的畸變和因此導致的不收斂情況發生。
可見，對于大變形中的網格畸變問題，我們只要運用ABAQUS的map solution功能，并結合強大的前處理軟件ANSA以及ABAQUS/CAE聯合作用，便可以在網格發生畸變之前進行重構，并映射前一步結果，再通過重啟動分析使計算繼續的進行下去，使得網格畸變不再是影響大變形不收斂的限制性因素。
來源：有限元在線的博客，版權歸作者所有。
展開 【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中mises應力云圖顯示的最大值還不到屈服應力值為啥還有PEEQ值
ABAQUS中mises應力云圖顯示的最大值還不到屈服應力值為啥還有PEEQ值，PEEQ云圖有變形值
ABAQUS中求解某部分單元的平均應力或平均應變 ￥10
1、參考模型：單向纖維的RVE模型；
2、腳本功能：針對指定的單元集合，在后處理中求解平均應力和平均應變。
3、應用的公式：一階均勻化計算方法。對于 RVE 模型的平均真應力和平均真應變，可通過對 RVE 內每一個單元的真應力 （真應變）取均值獲得。使用一階均勻化計算方法輸出的應力和應變適用于各種邊界條件，但需要對每個單元進行應力（應變）的輸出和計算。
abaqus模擬圓孔結構中應力集中分析 ￥19.89
<p>1.問題描述 </p><p><br></p><figure style="text-align: center;" class="ql-align-center">
<figure class="figure-image" contenteditable="false" data-img="https://img.jishulink.com/202502/attachment/f337440f18204bb49994cddb8926c825.png" style="display: inline-block;">
<img src="https://img.jishulink.com/202502/attachment/f337440f18204bb49994cddb8926c825.png" data-mobile-src="https://img.jishulink.com/202502/attachment/f337440f18204bb49994cddb8926c825.png?image_process=/format,webp" data-pc-src="https://img.jishulink.com/202502/attachment/f337440f18204bb49994cddb8926c825.png?image_process=/format,webp" data-initial-src="https://img.jishulink.com/202502/attachment/f337440f18204bb49994cddb8926c825.png">
</figure>
</figure><figure style="text-align: center;" class="ql-align-center">
展開 ABAQUS 噴丸殘余應力分析 ￥10
本案例適合哪些人學習：
1、學習型仿真工程師
2、理工科院校學生
3、對有限元分析感興趣的工程師
你會得到什么：
1、掌握噴丸三維模型的繪制
2、掌握顯示動力學分析相關的材料參數設置
3、理解顯示動力學分析步的建立
4、學習噴丸強化分析的相互關系的設置
5、了解顯示動力學網格的劃分
6、學習結果后處理的查看與對比
案例介紹：
所使用軟件為ABAQUS2018.
本案例完整得提供了分析相關所有的分析文件。
?
通過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腳本代碼（代碼中作了必要的簡單注釋）。
展開 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的工字梁靜應力分析
對于工字梁鋼架結構，通常都需要對鋼梁的應力及撓度進行分析，本案例通過Abaqus有限元分析軟件對此材料力學的問題進行了分析，相比利用材料力學知識的理論計算，不僅結果計算準確簡潔還能將結果可視化呈現，這也是本文采用Abaqus仿真分析的一大優勢。
2模型建立
本案例采用的工字型鋼梁的幾何尺寸如圖1所示，其中鋼梁的材料特性如表1所示。
圖1工字型鋼梁的結構尺寸
表1工字型鋼梁的材料特性
彈性模量E（N/m2）
泊松比μ
屈服強度fy（N/m2）
2.1E11
0.3
3.45E8
3.有限元建模分析
3.1單位制
這里不論是否使用abaqus軟件，只要是使用有限元仿真程序，均需要首先明確建模分析使用的單位制，保證單位制的統一。單位制的統一根據牛頓定律（F=ma）進行相互推導確認，這里不再贅述。
3.2ABAQUS建模
啟動ABAQUS/CAE，進入分析界面，依照有限元分析一般流程步驟：前處理→求解→后處理。我們通過工具欄的module模塊選擇part模塊進行工字梁幾何建模，按照圖1所示的幾何尺寸建立的模型如圖2所示。接著在property功能模塊定義工字梁的材料本構和截面屬性，并將截面屬性賦予相應的區域上。這里鋼梁的材料本構選用理想的線性彈塑性本構，即將切變模量（塑性應變）設置為0完成材料本構設置。應當注意要設置截面的屬性并給部件賦予上截面屬性。這是不同于ANSYS軟件分析流程的地方，ABAQUS不能把材料屬性直接賦予模型，而是先定義模型的截面，將材料屬性定義在截面上。之后定義裝配件并設置分析步。
展開 ABAQUS應力場
求大神告知abaqus中考慮計算出的應力場，是通過預定義場中的Initial-Mechanical-Stress，還是Step 1-Other-Field？
