abaqus應力梯度的案例
LMS Virtual.Lab Durability_方法介紹11—高級分析：應力梯度修正
今天帶來LMS Virtual.Lab Durability高級分析：應力梯度修正
內容主要包括應力梯度，基于應力疊加的修正以及應用實例。
11LMS Virtual.Lab Durability高級分析—應力梯度修正.pdf
百度網盤鏈接http://pan.baidu.com/s/1pJuOgv5
（該目錄下“11LMS Virtual.Lab Durability高級分析—應力梯度修正.pdf“）
LMS Virtual.Lab Durability交流群，群號：捌叁捌伍叁柒捌零 歡迎各位入群討論交流。
展開 不銹鋼表面Fe-Al梯度涂層的ANSY殘余應力仿真分析
可見，在涂層邊緣處存在應力集中，尤其在涂層明銳界面特別顯著。這主要是由于在涂層邊緣處冷卻速度快，形成了較大的溫度梯度造成的。等效應力在第二層涂層處較大，主要是由于第二層涂層材料的熱膨脹系數較大，同時彈性模量較大造成的。徑向應力在梯度中間層產生拉應力，在基體內部產生壓應力，容易使涂層剝落。軸向應力沿厚度方向為均為拉應力，在基體與涂層的邊緣處存在應力突變。
圖4
圖5
圖6
圖（4）為對稱軸等效應力沿軸向的分布圖。可見在涂層與基體的界面處存在應力突變，這主要是由于涂層與基體熱膨脹系數不一致引起的，同時在第三層軸向應力最大，這與圖（1）的結果相吻合。
圖（5）—（6）分別為基體與第一層涂層界面處等效應力、徑向應力的分布圖。等效應力顯示在邊緣處應力突變，存在應力集中。
在求解分析之后，你可以用origin來繪制文中的曲線圖，以使曲線美觀。同時也可以研究涂層的變化順序，涂層的厚度，基體的尺寸等對結果的影響進而優化你的實驗。
展開 ABAQUS梯度晶體FGM二維模型
本案例介紹在Abaqus CAE內建立呈現不同梯度分布模式的二維Voronoi晶粒結構模型。
模型輪廓草圖預先在AutoCAD內建立，在“0”圖層上建立正方形，在“hole”圖層建立內部的孔，這里的孔采用的是正多邊形，以確保能以多邊形的邊長生成對應的梯度晶粒。圖形建立完成后，采用CAD二維圖形Voronoi劃分 V2.0插件進行梯度晶粒的生成，晶粒直徑參數設置為最大的晶粒尺寸，晶粒類型選取梯度適應，邊界模式勾選自動尺寸。
在Abaqus內建立對應尺寸的二維部件，部件內部的孔可以建立為圓形。將CAD內生成的梯度晶粒以dxf草圖的形式導入Abaqus，并用其對建立的部件進行分區。
分區完成后也可采用Random Material Partition插件對不同區域隨機設置材料及比例。
沿直線分布的FGM梯度晶體模型只需在CAD草圖建立時將邊界線用多段線分段繪制即可，每段的尺寸與對應位置的晶粒尺寸一致。
可對模型劃分網格，并進行后續的梯度晶粒結構仿真模擬分析。
展開 Abaqus三維梯度泰森多邊形插件：Voronoi FGM 3D（Mesh）- AbyssFish ￥298
<h1><strong>插件介紹</strong></h1><p>Voronoi FGM 3D (Mesh) V1.0 - AbyssFish 插件可在Abaqus軟件內生成梯度分布的三維泰森多邊形長方體模型。插件可用于梯度功能材料(Functionally Gradient Materials)、梯度納米金屬材料、梯度金屬結構等梯度晶體模型的建立。模型基于背景網格實現，通過單元集的劃分，將不同的晶格指定不同的材料類型。</p><div contenteditable="false" width="100%">
<figure class="figure-image" data-img="https://img.jishulink.com/202403/attachment/eca5c533e92442f7bf70cd803eb60e88.png" style="text-align: center">
<img src="https://img.jishulink.com/202403/attachment/eca5c533e92442f7bf70cd803eb60e88.png" data-mobile-src="https://img.jishulink.com/202403/attachment/eca5c533e92442f7bf70cd803eb60e88.png?image_process=/format,webp/resize,w_400" data-pc-src="https://img.jishulink.com/202403/attachment/eca5c533e92442f7bf70cd803eb60e88.png?
展開 
ABAQUS功能梯度材料FGM模型
功能梯度材料（FGM）作為一種新型復合材料，通過材料內部成分或微觀結構的梯度變化，優化特定性能適應復雜環境，被廣泛應用于高溫防護、結構優化、生物醫學、光電設備等領域。本案例介紹在ABAQUS內建立功能梯度材料模型。
首先采用CAD 功能梯度材料2D插件建立大小呈現梯度分布的AutoCAD模型。
將圓形部分通過拉伸生成三維模型。
同樣將矩形拉伸為三維，并與生成的三維圓柱進行差集操作。并將兩個模型分別導出為iges格式文件。
在ABAQUS內將兩個模型以部件的形式導入。
裝配到一起。
插件也可生成中心梯度分布的模型。
或兩種不同的材料在擴散狀態下產生的梯度分布模型。
展開 ABAQUS UEL-梯度損傷模型應用于4節點平面單元 ￥600
針對該問題，梯度損傷（Gradient-damage）模型的概念被提了出來。
本文詳細介紹了如何將梯度損傷模型應用于4節點平面單元，并在有限元模型中進行模擬。
ABAQUS提供了UEL（user defined element）給使用者進行開發。筆者利用UEL開發4節點平面單元，其邊界條件如下圖所示。其中，節點1的X、Y方向被限制住，節點2的Y方向被限制，節點4的X方向被限制，節點3、4的Y方向有豎向位移0.1mm。單元為100*100mm的二維正方形。
每個節點除了X和Y方向的位移，還帶有非局部應變（nonlocal strain）。
單個單元模型，
多個單元模型，
具體內容可參見知乎文章：
ABAQUS UEL-梯度損傷模型應用于4節點平面單元 - 知乎 (zhihu.com)
相應的input文件和uel代碼付費可見，
展開 ABAQUS三維功能梯度多孔結構材料FGM軸壓模擬
功能梯度多孔材料（FGM）通過梯度調控孔隙率，實現力學性能的連續分布，其彈性模量、強度等呈均勻變化。通過建立梯度多孔結構有限元模型，解析梯度參數對應力場及失效機制的影響，突破傳統試驗限制，優化設計。該研究對航空熱防護及生物醫用仿生植入體等功能化結構具有重要價值。本案例介紹在ABAQUS內建立三維梯度功能材料多孔結構模型，并對梯度結構模型進行軸心受壓力學仿真模擬。
三維梯度孔隙結構模型采用CAD球體功能梯度材料3D插件建立，模型建立完成后將梯度孔基體部分導出為iges格式。
將梯度多孔結構模型以部件的形式導入到ABAQUS內。
對模型設置材料屬性，這里采用EasyCDP插件快速生成C20混凝土塑性損傷材料模型并指派給部件。
設置軸心受壓載荷工況，將模型一端固定，另一端指定位移。
對模型劃分網格。
創建并提交作業，查看模擬結果。
展開 ABAQUS熱應力分析 附ABAQUS中初始地應力的施加下載
軋輥與Cu層的熱傳導系數
下載地址：ABAQUS中初始地應力的施加
【abaqus】個人筆記—應力奇異&應力平均&應力集中
【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腳本代碼（代碼中作了必要的簡單注釋）。
展開 ABAQUS中mises應力云圖顯示的最大值還不到屈服應力值為啥還有PEEQ值
ABAQUS中mises應力云圖顯示的最大值還不到屈服應力值為啥還有PEEQ值，PEEQ云圖有變形值
