abaqus 應力比的案例
Abaqus熱應力分析、熱誘導振動分析簡單實例
1、一個金屬懸臂梁，一端固支，初始溫度20℃，溫度突變到120℃時由于膨脹及邊界約束而產生熱應力，進而引起振動，這種振動就是熱誘導振動。
2、熱誘導振動分析的成功應用不多見，在哈勃太空望遠鏡曾因熱誘導振動問題而發生故障。現在對航天器的分析中，熱誘導振動屬于難點和重點。國內曾有人對衛星天線做過準靜態熱誘導振動分析，也有人對空間站太陽能電池陣的桅桿做過基于模態的熱誘導振動分析（可能類似Abaqus中的線性攝動分析）。
3、熱應力分析與熱誘導振動分析進行耦合分析，還有難度，問題是多方面的。下面僅就準靜態非耦合的熱誘導振動分析為例，介紹由熱應力引起的振動。
4、懸臂梁材料屬性：
Conductity: 300W/(mK)
Density: 3000kg/m3
Elastic: E=3e10Pa, ν=0.3
Expansion: 3e-5 K-1
Specific Heat: 300J/(kgK)
5、分析結果
6、詳細步驟
見附件。
Abaqus熱應力分析、熱誘導振動分析簡單實例-kxh.part4.rar
Abaqus熱應力分析、熱誘導振動分析簡單實例-kxh.part1.rar
Abaqus熱應力分析、熱誘導振動分析簡單實例-kxh.part2.rar
Abaqus熱應力分析、熱誘導振動分析簡單實例-kxh.part3.rar
展開 Abaqus熱應力分析、熱誘導振動分析簡單實例
1、一個金屬懸臂梁，一端固支，初始溫度20℃，溫度突變到120℃時由于膨脹及邊界約束而產生熱應力，進而引起振動，這種振動就是熱誘導振動。
2、熱誘導振動分析的成功應用不多見，在哈勃太空望遠鏡曾因熱誘導振動問題而發生故障。現在對航天器的分析中，熱誘導振動屬于難點和重點。國內曾有人對衛星天線做過準靜態熱誘導振動分析，也有人對空間站太陽能電池陣的桅桿做過基于模態的熱誘導振動分析（可能類似Abaqus中的線性攝動分析）。
3、熱應力分析與熱誘導振動分析進行耦合分析，還有難度，問題是多方面的。下面僅就準靜態非耦合的熱誘導振動分析為例，介紹由熱應力引起的振動。
4、懸臂梁材料屬性：
Conductity: 300W/(mK)
Density: 3000kg/m3
Elastic: E=3e10Pa, ν=0.3
Expansion: 3e-5 K-1
Specific Heat: 300J/(kgK)
詳細步驟
見附件。
Abaqus熱應力分析、熱誘導振動分析簡單實例-kxh.part4.rar
Abaqus熱應力分析、熱誘導振動分析簡單實例-kxh.part2.rar
Abaqus熱應力分析、熱誘導振動分析簡單實例-kxh.part3.rar
Abaqus熱應力分析、熱誘導振動分析簡單實例-kxh.part1.rar
展開 Abaqus應力線性化-ASME Sec VIII Div 2_壓力容器分析設計
關于壓力容器分析設計的討論大多是基于ANSYS的應力線性化，而這方面Abaqus的公開資料不多，其實Abaqus早期版本就提供了在CAE界面下進行應力線性化的操作，為方便初學者使用Abaqus進行壓力容器分析設計，這篇文章介紹一下Abaqus應力線性化。
01. 壓力容器分析設計規范
目前最成熟、使用最多的壓力容器規范是由美國機械工程師協會（ASME）的鍋爐及壓力容器委員會（BPVC）制定的，我國的壓力容器相關規范有GB150、JB4732、JB4734等。
壓力容器的分析設計有別于傳統設計，主要是指通過有限元計算來校核壓力容器的設計方法，在ASME的壓力容器規范中是ASME Sec VIII Div 2的部分，相當于我國的JB4732。
分析設計的重要環節是應力線性化，為什么要進行應力線性化呢？其實主要是因為壓力容器的不同類型的故障（失效）模式是由不同類型的應力引起的，所以ASME的研究人員將它們進行了應力分類。
壓力容器的各種失效模式
ASME壓力容器規范的應力分類
如上圖所示，這些應力的類別大致分為三類：一次應力、二次應力和峰值應力，它們分別對應不同的故障模式。
一次應力與總塑性變形（gross plastic deformation）有關；
二次應力（在一次應力的基礎上）與增量塑性坍塌（incremental plastic collapse）有關；
峰值應力（在一次與二次應力的基礎上）與疲勞失效（fatigue failure）有關。
展開 有限元理論基礎及Abaqus內部實現方式研究系列15： 殼的剪切應力
我們關注CAE中的結構有限元，所以主要選擇了商用結構有限元軟件中文檔相對較完備的Abaqus來研究內部實現方式，同時對某些問題也會涉及其它的Nastran/Ansys等商軟。為了理解方便有很多問題在數學上其實并不嚴謹，同時由于水平有限可能有許多的理論錯誤，歡迎交流討論，也期待有更多的合作機會。
iSolver介紹視頻：
http://www.yqgqt.org.cn/college/video/c12884
==第15篇：殼的剪切應力 ==
自編有限元應力的校核除了Mises等合力外，也應該校核各個應力分量。材料力學中六個應力分量如下：
其中Tau11，Tau22，Tau33為正應力，Tau12,13,23為三個剪切應力，對殼來說，Tau33=0，Tau12為面內剪應力，Tau13,23即為本文所說的橫向剪切應力。
最近在做iSolver殼的應力分量和Abaqus比對時，發現Abaqus的橫向剪切應力和預想的不一致。iSolver按照常用的殼的理論得到的剪切應力是個與厚度無關的常量，但Abaqus的橫向剪切應力分量TSHR13，TSHR23，在各個截面方向積分點section point不一樣。
花了點時間細致的研究了一下，猜測Abaqus中剪切應力TSHR13、23是真實應力，但有限元理論和iSolver中計算的是板殼近似理論中平均剪切應力。本章將介紹殼單元中實際的和板殼近似理論中的剪切應力，也猜測了Abaqus的內部實現流程，最后通過一個算例來驗算Abaqus中的真實的剪切應力，并通過iSolver來計算板殼理論的平均剪切應力。
1.1 殼的真實的剪切應力
剪應力是材料由于抗拒面之間的滑動而產生的沿表面方向的應力。殼的中間層存在剪切應力，這個可以通過下面簡單的例子驗證。
展開 
ABAQUS-真實應力和名義應力轉化
ABAQUS-真實應力和名義應力轉化.doc
ABAQUS中mises應力云圖顯示的最大值還不到屈服應力值為啥還有PEEQ值
ABAQUS中mises應力云圖顯示的最大值還不到屈服應力值為啥還有PEEQ值，PEEQ云圖有變形值
ABAQUS學習筆記—對應力應變的部分理解
之前關于后處理的一些文章，由于一些原因全部刪除，故今天開始重新開始分享一些關于ABAQUS的一些知識，希望能夠對大家有所幫助，也希望大家能夠繼續支持筆者。
那么今天，我們再對‘’ABAQUS中應力應變的部分理解的‘’內容進行講解。
在ABAQUS中，一般是把X軸當成1軸，Y軸當成2軸，Z軸當成3軸，那么：
S11就是X軸向的應力，正值為拉應力，負值為壓應力；
S22就是Y軸向的應力，正值為拉應力，負值為壓應力；
S33就是Z軸向的應力，正值為拉應力，負值為壓應力；
S12就是在YZ平面上，沿Y向的剪力；
S13就是在YZ平面上，沿Z向的剪力；
S23就是在XZ平面上，沿Z向的剪力；
由于剪力的對稱性:S11=S21,S13=S31,S23=S32；
由以上可知，S11，S22，S33為主應力；S12，S13，S23為切應力；
主應力分別以σ1，σ2，σ3表示，按數值排序為：σ1≥σ2≥σ3。在ABAQUS中分別對應為：Max.principal；Mid.principal；min.principal。這三個量在任何坐標下都是不變量。
我們可利于最大應力判斷一些情況：比如最大主應力（拉應力）大于混凝土的抗拉強度，則認為混凝土開裂；通過顯示最大主應力的法線方向，則可大致表示出裂縫的發展影響。
應變中一些符號的含義
E——總應變
EP——主應變
EE——彈性形變
PE——塑性應變分量
Eij——應變分量
Ie——非彈性應變分量
PEEQ——等效塑性應變。若該值大于0，則認為已經屈服
注：在ABAQUS后處理中，盡量不要看Mises，其表示平均應力，更適合金屬材料；對于鋼筋混凝土結構，我們要看其單軸拉伸方向上的應力和對應的應變。
展開 ABAQUS地應力平衡
對于涉及開挖、填充的動態巖土工程問題，初始應力場是正確模擬其施工過程的先決條件。初始應力的加載必須滿足地應力平衡，而地應力平衡其實就是為了使數值模擬獲得一個存在初始應力，而無初始應變的狀態。當考慮自重場是產生地應力場的主要因素時，顯然重力是外力，初始應力場是內力，將提取出的內力施加于數值模擬后再施加重力，是內力和外力平衡，從而獲得較精確的沒有受到人為因素干擾情況下的數值模型的初始狀態。
在 ABAQUS 有限元軟件中，提供了 5 種定應力平衡方法，這里我就比較常用的，使用比較普遍的一種進行講解。用的例子還是上節ANSYS的例子。
接下來，小編就一個簡單的算例就行地應力平衡。在一個均質的地基上有一根立柱。
柱子的尺寸為1m×2m×10m，密度2400kg/m3，彈性模量2.8E8Pa（這里我為了方便觀察規律，降低這個彈性模量），泊松比0.167；均質地基的尺寸為30m×30m×15m，密度為2000kg/m3，彈性模量3E7Pa，泊松比0.3。模型見圖1.
1、將模型導入ABAQUS中。
2、創建靜力荷載步
3、殺死柱子單元
4、施加約束（四周施加法向連桿約束，底部施加全約束）
5、施加自重荷載
6、提交計算
7、查看后處理結果
豎向位移
豎向應力
8、應力場輸出
9、在工作目中找到stress_out.rpt文件
10、經過格式的轉化（文件的后面不要有空行）
11、在之前的inp文件中Step-1前面加上*Initial conditions,Type=stress,Input=stress_in.rpt這句命令，記得一定要把stress_in.rpt放在工作目錄里面。
展開 有限元理論基礎及Abaqus內部實現方式研究系列14： 殼的應力方向
我們關注CAE中的結構有限元，所以主要選擇了商用結構有限元軟件中文檔相對較完備的Abaqus來研究內部實現方式，同時對某些問題也會涉及其它的Nastran/Ansys等商軟。為了理解方便有很多問題在數學上其實并不嚴謹，同時由于水平有限可能有許多的理論錯誤，歡迎交流討論，也期待有更多的合作機會。
iSolver介紹視頻：
http://www.yqgqt.org.cn/college/video/c12884
==第14篇：殼的應力方向 ==
有限元中，物理量用的最多的是標量、矢量和二階張量。其中位移、坐標等都是矢量，而應變、應力等都是二階張量。矢量很容易理解，體的應力等二階張量直接就采用了全局坐標系的也不會有方向理解問題，但梁殼的應力結果很容易搞錯，后處理結果中的S11、S12等的方向有時會覺得和預想的不一致但又不明所以。同時，這個方向也是單元材料的方向，在自編程序時，如果一開始坐標系的定義就弄錯了，那么將直接導致和材料相關的剛度矩陣的錯誤，所以弄清應力的方向定義對自編程序和理解有限元結果都相當重要。
本章將簡單介紹一下數學上張量和Abaqus中殼的應力方向，并說明Abaqus這么選取的意義，最后通過自編程序iSolver來驗證殼的應力方向的正確性。
具體的驗證詳見下方視頻（帶配音）：
https://www.yqgqt.org.cn/college/video/c12884 20.14 理論系列文章14：殼的應力方向
1.1 數學上的張量方向
矢量的方向是一定的，但它的分量都是基于某個坐標系定義的，坐標系不同，那么分量結果也會不同。
矢量可以表示為：
顯然，分量和坐標系的選取有關。譬如我們一般的直角全局坐標系如下，那么分量就是普通的x、y、z三個分量值。
展開 ABAQUS熱應力分析 附ABAQUS中初始地應力的施加下載
軋輥與Cu層的熱傳導系數
下載地址：ABAQUS中初始地應力的施加
Abaqus中平面應力單元高斯積分點的順序
可以輸出umat接口中的變量coords進行查看
write(*,"(A,I4)") "npt = ", npt
write(*,"(A,3ES16.8)") "coords = ", coords
結果為：
npt = 1
coords = -5.77350269E-01 -5.77350269E-01 1.00000000E-02
npt = 2
coords = 5.77350269E-01 -5.77350269E-01 1.00000000E-02
npt = 3
coords = -5.77350269E-01 5.77350269E-01 1.00000000E-02
npt = 4
coords = 5.77350269E-01 5.77350269E-01 1.00000000E-02
因此Abaqus中平面應力單元高斯積分點的順序為：
展開 
希望對大家有幫助
觀看須知.txt
Abaqus熱應力分析、熱誘導振動分析簡單實例-kxh.part1.rar
Abaqus熱應力分析、熱誘導振動分析簡單實例-kxh.part1.rar
觀看須知.txt
Abaqus熱應力分析、熱誘導振動分析簡單實例-kxh.part1.rar
ANSYS與ABAQUS比較之實例8---帶孔平板的熱應力分析1
本博文是關于ANSYS與ABAQUS比較之系列博文，本例子使用ABAQUS做熱應力分析，后面會使用ANSYS對同一個問題做熱應力分析。
【問題描述】
一個帶孔平板結構如下圖
該平板上邊沿固定，左右兩邊是滾動支座支撐。該板的初始溫度是25度，現在要求當溫度升高到150度時，板中的應力分布。
已知：材料的彈性模量是2e9pa, 泊松比是0.3,熱膨脹系數是1.35e-5/度。
【問題分析】
1.
分析類型。這是一個平面應力問題，應力的產生是因為溫度的變化導致產生了熱應變，而該熱應變又被約束限制導致熱應力的產生。
2.
非線性考慮。只有一個物體，不存在接觸非線性；沒有材料非線性；沒有幾何非線性。總之，這就是一個最簡單的線彈性分析。
3. 幾何建模。由于該結構左右對稱，只取一半研究。
4.
邊界條件。除了常規的位移邊界條件以外，對該板施加預定義溫度場25度，而在第一個分析步修改該溫度場的溫度為150度。
【求解過程】
1. 創建部件
只取一半建模，它是一個二維的可變形部件。
2. 定義材料屬性
只需要定義彈性模量，泊松比及線膨脹系數。
3. 定義截面屬性
創建均質的實體截面，并將上述材料屬性賦予給它，然后將該截面屬性賦予給前面的部件。
4. 裝配部件
唯一的部件，導入到裝配即可。
5.設置分析步
兩個分析步。
新創建的分析步是最一般的靜力學通用分析步。
6.定義載荷和邊界條件
首先定義位移邊界條件，在初始分析步中，固定上邊，左右兩邊施加X方向的位移限制。
使用預定義場確定溫度。
對整塊板施加25度的初始溫度。
展開 基于ABAQUS的初始地應力平衡方法研究_代汝林.pdf
基于ABAQUS的初始地應力平衡方法研究_代汝林.pdf
abaqus拉伸后處理（應力-應變，位移-力的輸出）
abaqus拉伸后處理（應力-應變，位移-力的輸出）
lashen.zip
