abaqus應力形式的案例
彈性力學基本方程之應力應變的矩陣形式
應力列陣（應力向量）
應力分量正負號規定
如果某一個面的外法線方向與坐標軸的正方向一致，這個面上的應力分量就以沿坐標軸正方向為正，與坐標軸反方向為負；相反如果某一個面的外法向方向與坐標軸的負方向一致，這個面上的應力分量就以沿坐標軸負方向為正，與坐標軸同向為負。應力分量及正方向如圖1
ABAQUS中輸出Cohesive單元的斷裂形式
Cohesive單元可以模擬三種基本斷裂形式：I型張開裂縫、II型滑移裂縫、III型撕開裂縫。其中后兩種破壞形式的驅動力為剪切力，所以可以稱為剪切破壞。
ABAQUS中其實已經提供了識別這兩種（拉伸和剪切）破壞形式的場輸出：MMIXDME、MMIXDMI。下面從定義、使用、結果演示等幾個方面進行介紹。
定義
幫助文檔中，二者定義如下：
MMIXDME，為損傷演化過程中混合斷裂模式的比例，定義為1-m1，位于單元積分點位置，單元未破壞時，其數值默認為-1；
MMIXDMI，為初始損傷時混合斷裂模式的比例，同樣定義為1-m1，位于單元積分點位置，單元未破壞時，其數值默認為-1；
后面看下m1的定義：
m1為Gn（一型張開斷裂能）與GT（三型斷裂能之和）的比值：
當m1為1，表示完全的拉伸破壞時，對應的MMIXDME和MMIXDMI數值為0；
當m1為0，表示完全的剪切破壞時，對應的MMIXDME和MMIXDMI數值為1；
通過這兩個參量，我們就可以判斷cohesive單元以哪種形式破壞為主。
使用方法
目前不支持GUI界面定義，Step模塊的場輸出無法找到MMIXDMI和MMIXDME。
展開 有限元理論基礎及Abaqus內部實現方式研究系列45：約束關系（1）-統一形式
Nastran的MPC的靈活度要遠遠超過Abaqus，Nastran的主節點可以選擇123自由度，也可以對每個從節點設置不同的自由度，還能主節點和從節點互相包含，Abaqus更多的是只負責80%的常用應用場景，復雜功能讓你編子程序，但事實上一線仿真工程師又有多少人愿意編子程序呢？這種做法導致雖然Abaqus無論從用戶體驗、非線性還是商業化都比Nastran好很多，但很多線性的工程復雜問題還是沒法替代Nastran。
（2） 另一類是Contact、Tie等的面之間的約束關系。在這方面Abaqus要明顯強于Nastran了。
我們將用統一的公式來求解這兩類關系，同時也從軟件實現層面說明一下針對這兩類情況的各自差異。分幾篇文章來介紹約束關系，本篇是約束關系（1）-統一形式，既然接觸僅是約束關系的一種，那么MPC、Tie、接觸等的求解過程也是很類似的，這里將介紹一下這些約束關系如何表達為統一形式。
3 統一形式的約束關系
在沒有約束關系時，如下圖情況，物體在體外力和面外力作用下變化。
有限元方程按照虛功原理求解，在物理上可解釋能量守恒原理，即在某一個時刻點，假定在外力作用下有個虛擬的位移，那么外力在虛擬位移下做的虛功=內部應變能的變化相同。
虛功原理中的每項都表示各自區域在虛位移下的能量變換
（1） fv是每單位體積內的力，外力fv和位移相乘表示單位體積內的虛功，所以對體積積分
（2） fs是每單位面積上的力，外力fv和位移相乘表示單位面積上的虛功，所以對面積積分，推論就是最后一項應力和應變是單位體積內的內能，所以對體積積分。
當存在約束關系時，在能量中加入約束關系相關的一項：
顯然的單位也是該約束關系所在區域在虛位移下的能量單位。
展開 怎樣將Abaqus中的材料設置導出成文本形式？
怎樣將Abaqus中的材料設置導出成文本形式？
在ABAQUS中如何采用DISP或者VDISP子程序模擬地基中地下水位的升降（以正弦波形式）？ ￥200
在ABAQUS中如何采用DISP或者VDISP子程序模擬地基中地下水位的升降（以正弦波形式）？
【abaqus】個人筆記—應力奇異&應力平均&應力集中
【abaqus】個人筆記—應力奇異&應力平均&應力集中
ABAQUS熱應力分析 附ABAQUS中初始地應力的施加下載
軋輥與Cu層的熱傳導系數
下載地址：ABAQUS中初始地應力的施加
Abaqus帶螺紋螺栓接觸應力分析淺析 Abaqus帶螺紋螺栓接觸應力分析淺析
目前的常規做法通常有兩種：1.簡化，用RBE2和beam梁來代替螺栓，這樣不能反映連接螺栓真實應力，圖1為某結構連接螺栓簡化的beam梁應力云圖，沒有接觸應力：
.直接做出來螺栓螺紋采用接觸分析，雖然得出的結果很精確，但這樣前處理工作量大（螺栓和螺紋用六面體網格建模）、計算量大（接觸收斂困難），如圖為某結構帶螺紋螺栓和連接件模型（圖2）和計算得出的結果（圖3）：
圖3 計算結果
那么，有什么好辦法可以不用簡化帶螺紋螺栓，不用直接做出帶螺紋螺栓，又能得到足夠精確的結果？
運用大型通用非線性有限元分析軟件Abaqus，只需要在接觸定義中設置跟實際螺紋形狀有關聯的參數，如牙角、螺距、螺栓小徑等，就可以模擬真實的連接螺栓接觸狀況。既可以得到足夠精確的分析結果，又節省了時間專注進行其他的分析設置。如圖4，為連接螺栓接觸來定義帶螺紋螺栓：
圖4 連接螺栓接觸來定義帶螺紋螺栓
圖5為某結構直徑10MM的帶螺紋的連接螺栓接觸壓力分布云圖：
圖5 某結構直徑10mm帶螺紋的連接螺栓接觸壓力分部云圖
展開 abaqus里的非線性薄層單元，零厚度cohesive單元，goodman接觸單元等的基本形式是什么？如何構建與應用？
在使用Abaqus，Comsol等軟件進行薄層區域的力學分析過程中，例如在研究水壓致裂、裂縫擴展，接觸粘結滑移的這類薄層力學性質時，我們經常需要采用應力-相對位移(σ-u)關系，而不是傳統本構描述的應力-應變(σ-ε)關系來描述，例如Abaqus里面的Cohesive單元，Goodman單元，以及Comsol里的彈性薄層（在后面我把這類單元統稱為增量非線性力學薄層）。這類單元厚度非常小甚至為0，薄層兩側的節點（單元）用一組力（應力）與相對位移的關系方程聯系起來，例如給出一個形式最為簡單的典型應力-位移方程
此方程描述了1,2,3方向（通常是法向和兩個切向）上相對位移與應力的關系，應力與相對位移呈線性關系，類似于“線性彈簧”。但是對于土-結構接觸、裂縫的張開閉合這類問題，線性方程已經不足以準確描述這些物理量之間的關系，這時就需要引入增量非線性方程來構建薄層單元。
引入增量非線性薄層的概念之前，首先介紹一下全量非線性薄層以理解非線性的概念，首先給出以下公式
這是一個全量非線性薄層，其非線性的表現可以用下面幾個例子體現，
對比①和②項，可以發現僅存在3方向上的位移變化的情況下，1，2方向上的力也會發生改變，體現了彈簧三個方向力學性質的非獨立性，對比①和③項，可以發現力的大小并不和位移大小成正比，也就是非線性特征。
所以對于增量非線性方程，就是把應力-位移關系方程寫成應力增量-位移增量的關系方程，例如
寫成微分形式的好處是，可以體現出應力路徑對位移結果的影響，也就是類似于“塑性”特征（所以所有的彈塑性本構也都是增量方程）。但是對于此類微分方程的求解，必須給定一個力的初始值。
展開 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 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模擬圓孔結構中應力集中分析 ￥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">
展開 