ansys在面上應力積分的案例
將應力拖到屈服面上程序代碼
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中平面應力單元高斯積分點的順序為：
展開 一個斷裂分析實例（J積分和應力強度因子）
（參看圖1.1）
（2） 載荷為均布拉力q=1200MPa時的J積分值。（參看圖1.1）
（3） 當載荷為轉速n=600r/min的KI ，并計算沿裂紋尖端不同路徑的積分值，與KI比較。
二、 求解
求解中統一采用國際單位制，長度m，壓力、應力與彈性模量Pa，密度Kg/m3，轉速rad/s。
對圓盤的1/4進行ANSYS建模，網格劃分如圖2.1。單元類型為6節點三角形單元plane2。裂紋附近單元邊長為0.0002m。載荷施加如圖2.2，扇形兩條半徑（裂紋處除外）上施加對稱位移邊界條件，弧上加均布拉力。裂紋處無位移約束。
計算的應力強度因子和J積分結果如表2.1。前面給出的J積分由于坐標系錯誤，做法不對，現在改正1200Mpa下的J積分結果。但是，結果差別不大。因為，即使在局部坐標系下，J積分中用到的 XG, YG, ZG的坐標還是在全球坐標系下的。坐標系的改變只對J積分的第一項有一點影響（對Y坐標積分這部分，Y坐標變了）。
新加兩個附件：改正后求解1200MPa下J積分的命令流，及兩張路徑定義的圖片命令流里面涉及到路徑定義的命令是GUI方式進行的。所以要分開看。
simwe_060715_001.rar
。
展開 單元積分點應力如何外插至節點上 | 數值實現篇
繼上次的推文：有限元計算過程中積分點應力如何外插至節點處？【公式推導篇】，本次分享單元積分點應力外插至節點處的數值實現過程。
數值實現
借助以上理論，我們可以基于matlab平臺編制以下代碼段：
% 將積分點應力外插至單元節點上，這里只列舉了Q4的情況
for i = 1:3
StressElem(e,:,i) = [1+0.5*sqrt(3) -0.5 1-0.5*sqrt(3) -0.5;
-0.5 1+0.5*sqrt(3) -0.5 1-0.5*sqrt(3);
1-0.5*sqrt(3) -0.5 1+0.5*sqrt(3) -0.5;
-0.5 1-0.5*sqrt(3) -0.5 1+0.5*sqrt(3)]*...
[stress(e,1,i);stress(e,2,i);stress(e,3,i);stress(e,4,i)];
end
對標Abaqus
模型材料參數為普通的線彈性材料，單元類型選擇CPS4，網格劃分及邊界條件設置如下：
在結果對標過程中，可以先對比自研程序與Abaqus的節點位移場：
Abaqus位移場結果
自研程序位移場結果
在位移場一致的前提下，我們再來對標應力結果。以常見的mises應力為例：
Abaqus位移應力場結果
自研程序應力場結果
結果是一致的，說明了程序的正確性。
展開 
分享一個斷裂分析實例（J積分和應力強度因子）
（參看圖1.1）
（2） 載荷為均布拉力q=1200MPa時的J積分值。（參看圖1.1）
（3） 當載荷為轉速n=600r/min的KI ，并計算沿裂紋尖端不同路徑的積分值，與KI比較。
二、 求解
求解中統一采用國際單位制，長度m，壓力、應力與彈性模量Pa，密度Kg/m3，轉速rad/s。
對圓盤的1/4進行ANSYS建模，網格劃分如圖2.1。單元類型為6節點三角形單元plane2。裂紋附近單元邊長為0.0002m。載荷施加如圖2.2，扇形兩條半徑（裂紋處除外）上施加對稱位移邊界條件，弧上加均布拉力。裂紋處無位移約束。
計算的應力強度因子和J積分結果如表2.1。前面給出的J積分由于坐標系錯誤，做法不對，現在改正1200Mpa下的J積分結果。但是，結果差別不大。因為，即使在局部坐標系下，J積分中用到的 XG, YG, ZG的坐標還是在全球坐標系下的。坐標系的改變只對J積分的第一項有一點影響（對Y坐標積分這部分，Y坐標變了）。
新加兩個附件：改正后求解1200MPa下J積分的命令流，及兩張路徑定義的圖片命令流里面涉及到路徑定義的命令是GUI方式進行的。所以要分開看。
展開 有限元計算過程中積分點應力如何外插至節點處？【公式推導篇】
注：由于技術鄰排版風格有限，故部分內容顯示不全，感興趣的小伙伴可點擊原文進行閱覽：
有限元計算過程中積分點應力如何外插至節點處？【公式推導篇】
https://mp.weixin.qq.com/s/47byQ3b3e5UpbUp7Krs2mQ
本次分享的是：有限元計算過程中，單元積分點應力如何外推至節點？
有關積分點與節點的概念可點擊跳轉閱讀歷史推文：有限元基本概念-【節點和積分點】，現科普一下Q4單元、Q8單元、Q9單元的形函數和高斯積分方案。
Q4單元
Q8/9單元
應力外插
核心理念：坐標系的轉換。
假設是母單元的自然坐標系，是由高斯積分點控制的坐標系（術語可能不專業），假設高斯積分方案為。坐標系轉換關系：
單元內任一點的應力，由4個高斯積分點應力進行插值時，可表示為
其中，是基于高斯積分點的形函數，第一個積分點的坐標在母單元坐標系下為（-1，-1），根據上述的坐標系轉換的方式，在高斯積分點的坐標系下，第一個單元節點在高斯積分點坐標系下坐標為，將此坐標值代入第一個形函數，得，相同的道理，可推導至四個節點在4個形函數下的外插矩陣：
對于Q8、Q9單元，依然可采用高斯積分方案（減縮積分）。
展開 有限元中單元積分點與節點應力相互轉換(CPE4為例)（ABAQUS）
（注意：變量是a,b,c,d,而不是x,y.所以方程組是線性的）
第一個積分點的應力和坐標：S1，(X1，Y1);
第二個積分點的應力和坐標：S2，(X2，Y2);
第三個積分點的應力和坐標：S3，(X3，Y3);
第四個積分點的應力和坐標：S4，(X4，Y4);
現在的問題是：應力分量S1，S2，S3，S4是已知的，我們需要知道真實的積分點的坐標信息嗎？
答案：不需要，只需要知道積分點在整個單元相對位置即可。即等參元中的坐標。（教材中有）
等參元的長和寬都為2.
而有限元中的積分是高斯積分，積分點的位置是固定的。由查表可知：
上表是一維的高斯積分點的坐標，后面的加權系數不用管（我們不求積分）。由一維可以猜出二維（兩個一維）。二維有4個積分點，所以我們對應一維選第二行的數據。
展開 ansys 如何添加圓柱面上小部分小角度的徑向均布載荷
ansys 如何添加圓柱面上小部分小角度的徑向均布載荷，也就是說在圓柱面上的一小段，比如說120mm的圓柱，在其中間60mm的一段上，60度的扇形面上添加均布的徑向載荷？
ansys J-積分
J—積分計算方法
J 積分_命令流.doc
J積分_GUI具體步驟.doc
J積分_基于ANSYS的J積分計算與分析.pdf
斷裂／裂縫／裂紋／應力強度因子和Ｊ積分計算實例（原創，如轉載，請注明出處）
材料：復合材料
分析類型：斷裂力學
技術難點：斷裂 應力強度因子 Ｊ積分
完成人：技術鄰ANSYS專家
網址：http://www.yqgqt.org.cn/content/other/402981
模擬過程：
斷裂力學計算應力強度因子和Ｊ積分
圖１ 裂縫模型
圖２ Ｋ１應力強度因子
圖３ Ｋ２應力強度因子
圖４ Ｋ３應力強度因子
圖５ Ｊ積分
圖６ 裂縫前沿集中應力
展開 ANSYS中的LAREA命令——在面上的兩個關鍵點之間生成最短的線
注：在面上的兩個關鍵點P1和P2之間生成一條最短的線，生成的線也位于面內。P1和P2關鍵點也可以在面的同一側，且到面的距離相等，這種情況下則生成一條與面相平行的線。
2.操作路徑
Main Menu> Preprocessor> Modeling> Create> Lines> Lines> Overlaid on Area
3.實例
輸入命令：
/PREP7
CYL4,0,0,1,,,,2
LAREA,6,7,2 !在平面上生成最短的線L11
LAREA,1,5,4 !在圓柱面上生成最短的線L12
則生成的圖線如圖1所示
圖1 生成的圖線
4.參考資料
ANSYS HELP 15.0
展開 
【JY】ANSYS Workbench在減隔震應用分析中的單元積分技術筆記
黏滯阻尼器的固流耦合分析：
對于ABAQUS的單元介紹已經做了詳盡，個人感覺固體力學上ABAQUS還是上手比較方便，而多場耦合、快速建模預估Workbench會方便一些，因人而異：
【JY】有限單元分析的常見問題及單元選擇
ANSYS Workbench就像一個科技界的“瑞士軍刀”，集合了各種強大的單元技術，為減隔震元件提供全面且準確的分析支持。近期對于ANSYS Workbench進行了學習，本文將對ANSYS Workbench 各類單元技術做一個筆記總結，便于為減隔震元件分析提供理論基礎。（畢竟Workbench大部分時候會自動匹配相應所需技術）
B-bar方法完全積分
Workbench中的B-bar方法是一種常用于處理低階單元完全積分的技術，也被稱為選擇性減積分策略。它是針對有限元分析（FEA）中的一種改進方法，旨在提高計算效率和準確性。
在傳統的有限元分析中，低階單元（如線性單元）在處理不可壓縮材料或近似不可壓縮材料時，常常遇到體積鎖定問題。體積鎖定是指在近似不可壓縮材料的有限元模擬中，由于體積應變被過度限制，導致計算結果偏離實際情況的現象。為了解決這個問題，B-bar方法被引入到ANSYS Workbench中。
B-bar方法的核心思想是在低階單元的完全積分過程中進行選擇性減積分。它通過將高斯積分點處的體積應變替換為單元的平均體積應變，實現了對應變的軟化處理，從而防止了體積鎖定的發生。這種選擇性減積分的策略可以在保證計算精度的同時，提高計算的收斂性和效率。
需要注意的是，B-bar方法并不能解決剪切鎖定問題，這是另一種常見的有限元分析問題。對于彎曲主導的問題，剪切鎖定可能導致結果的失真。因此，在處理這類問題時，用戶需要采用其他方法，如使用增強應變公式等。
展開 ANSYS后處理中的應力與屈服準則！
ansys后處理該看的那些應力
01
應力
材料發生形變時，內部產生了大小相等但方向相反的反作用力，抵抗外力把分布內力在一點的集度稱為應力 (Stress)，應力與微面積的乘積即微內力或物體由于外因（受力、濕度變化等）而變形時，在物體內各部分之間產生相互作用的內力，以抵抗這種外因的作用，并力圖使物體從變形后的位置回復到變形前的位置。我們分析后查看應力，目的就是在于確定該結構的承載能力是否足夠。那么承載能力是如何定義的呢？比如混凝土、鋼材，應該就是用萬能壓力機進行的單軸破壞試驗吧。也就是說，我們在ANSYS計算中得到的應力，總是要和單軸破壞試驗得到的結果進行比對的。所以，當有限元模型本身是一維或二維結構時，通過查看某一個方向，如plnsol,s,x 等，是有意義的。但三維實體結構中，應力分布要復雜得多，不能僅用單一方向上的應力來代表結構此處的確切應力值——就出現了強度理論學說。
材料力學中的四種強度理論
01
最大拉應力強度理論
該理論認為，材料破壞的主要因素是最大拉應力，無論何種狀態，只要最大拉應力達到材料的單向拉伸斷裂時的最大拉應力，則材料斷裂。其中，某點的最大拉應力數值，就是其第一主應力數值。
展開 ANSYS后處理中的應力與屈服準則
ansys后處理該看的那些應力
01
應力
材料發生形變時，內部產生了大小相等但方向相反的反作用力，抵抗外力把分布內力在一點的集度稱為應力 (Stress)，應力與微面積的乘積即微內力或物體由于外因（受力、濕度變化等）而變形時，在物體內各部分之間產生相互作用的內力，以抵抗這種外因的作用，并力圖使物體從變形后的位置回復到變形前的位置。我們分析后查看應力，目的就是在于確定該結構的承載能力是否足夠。那么承載能力是如何定義的呢？比如混凝土、鋼材，應該就是用萬能壓力機進行的單軸破壞試驗吧。也就是說，我們在ANSYS計算中得到的應力，總是要和單軸破壞試驗得到的結果進行比對的。所以，當有限元模型本身是一維或二維結構時，通過查看某一個方向，如plnsol,s,x 等，是有意義的。但三維實體結構中，應力分布要復雜得多，不能僅用單一方向上的應力來代表結構此處的確切應力值——就出現了強度理論學說。
材料力學中的四種強度理論
01
最大拉應力強度理論
該理論認為，材料破壞的主要因素是最大拉應力，無論何種狀態，只要最大拉應力達到材料的單向拉伸斷裂時的最大拉應力，則材料斷裂。
展開 應力集中問題與ANSYS驗證
在工程上，應力集中的程度用局部最大應力σmax與該截面上的名義應力σnom的比值來表示，即
Ktσ=σmax/σnom
Ktσ稱為理論應力集中系數。下面，我們將通過一個典型應力集中問題——帶孔平板，使用ANSYS軟件求出最大應力和應力分布圖，并與彈性理論計算的結果進行比較：
根據彈性力學知識，孔邊環向正應力的大小是無孔時的3倍，隨著遠離孔邊而極速趨近于q。
ANSYS求解：
Step1：在SCDM中創建平面模型。
由于我們使用平面應力模型計算，所以建模時必須要將橫截面建立在xy平面上。建立一個邊長為20mm×10mm的平面模型，中間孔的直徑為2mm。我們將模型分為四部分，方便在每部分的邊界上設置Path，從而繪制應力曲線。由于該模型同時關于X軸和Y軸對稱，我們也可以使用四分之一模型建模。此處筆者使用完整模型。建立完成以后，使用share命令共享拓撲，然后點擊菜單欄Workbench→ANSYS transfer→2020R1進入Workbench。
Step2：設置分析類型（2D）。
在Project Schematic中的空白處點擊右鍵，選擇Properties，打開Properties of Project Schematic。單擊項目中的A3（Geometry）欄，在Propertiesof Project Schematic A3: Geometry中將AnalysisType切換為2D。（若Analysis Type為3D，則導入平面幾何后軟件將使用殼單元計算。）
展開 