abaqus網格積分
關注創建者：王靖雯 創建時間：2023-02-27
abaqus網格積分的視頻教程
ABAQUS材料斷裂與失效系列 之 圍道積分運算
本專題分一下幾部分進行講解： 1.1 靜態裂縫的兩種模擬方法（達索09教程） 1.2 圍道積分運算（達索09教程） 1.3 圍道積分相關的版本更新 1.4 示例1：無限空間內圓形裂縫的圍道積分運算的五種方法五種方法： 1.4.1 軸對稱尖銳裂縫模型 1.4.2 軸對稱鈍角裂縫模型 1.4.3 三維尖銳裂縫模型 1.4.4 擴展有限元（XFEM）方法 1.4.5 子模型方法
￥50 3小時15分鐘 3034播放
查看
選擇積分與 Abaqus 梁單元內核
首先討論選擇積分與降階積分的理論爭議（Bathe vs 王勖成），分析零能模式與剪切鎖死的數學機制。隨后通過 Timoshenko 梁經典例題，對比經典梁、精確積分與縮減積分三種結果，解釋 25% 誤差來源。接著深入 Abaqus 梁單元理論，介紹中心線描述、變形梯度分解、四元數大轉動更新及虛功方程。最后說明普通梁、開口薄壁梁與混合梁單元的選型邏輯，并引入張量分析基礎。
￥65 52分鐘 1播放
查看
ABAQUS中利用圍線積分法輸出裂紋尖端積分值
基于ABAQUS的官方幫助文檔，詳細講解了ABAQUS中裂紋體的創建方法+裂紋尖端奇異性單元的構件方法（多種單元類型計算精度對比）+裂尖積分值的輸出設置
￥6 34分鐘 167播放
查看
abaqus網格積分的實例教程
之前算過一個關于裂紋擴展的問題，當時創建裂紋選擇的是contour intergral，后來又有人咨詢我裂紋尖端J積分的計算問題。我才恍然大悟，其實圍道積分方法還是適用于計算裂紋尖端在某時刻的J積分，至于動態擴展問題，還是交給XFEM吧（雖然也不太好）。
計算了幾種情況下的裂紋尖端J積分，包括直裂紋、斜裂紋以及裂紋尖端傾斜等三種情況。
部分試件的應力分布及J積分結果如圖所示：
參考文獻：《Physically based crystal plasticity FEM including geometrically necessary dislocations: Numerical implementation and applications in micro-forming》
GND 演化方程依賴依賴于剪切應變率的梯度或者塑性變形梯度的旋度，而標準FEM/VUMAT 只告訴你每個積分點本身的 γ˙a、Fp，不會直接給梯度。以往廣泛應用的數值方案通常是：先把 積分點的數據外推到節點，再用線性形函數求梯度，然而這類方案只能用特定單元（如 C3D8），對自適應網格、復雜接觸不友好。
該文章提出的一個mesh-free的方案，該方案的主要優勢是不改單元、不加 DOF，只在材料子程序內部，用鄰近積分點的數據做一次局部重構，就算出梯度，該策略對某個積分點 x，附近有一團“鄰居積分點” xI，作者把它們當成 mesh-free 的“節點”，對每個場變量 u(x)（可以是 γ˙a，Fp 的分量）做 MLS 擬合，如下圖所示：
權函數使用立方樣條，有緊支撐，距離越近權越大：
在實現上作者提到，立方支撐三個方向尺寸約為5個單元尺寸，最多取最鄰近60個（3D）或者30個（2D）積分點，作者指出：當鄰域尺寸比網格尺寸還小的時候，這個非局域模型就自然退化為傳統的局域模型。也就是說，鄰域尺寸本身就扮演了“材料內在長度標度”的角色。
為了提高計算效率秒作者使用了一個“時間滯后 + 公共塊”的策略對GND進行更新。
展開 方法
在ABAQUS CAE的場輸出中選擇的坐標點是節點的坐標，而節點是從積分點插值出來的，單元積分點的信息相對真實。所以最好是獲取積分點的信息，其中積分點的坐標無法在CAE中獲取，需要在關鍵字中添加。具體在每個分析步的單元輸出下面添加COORD，如果需要輸出節點的坐標也可以在節點場輸出下面添加COORD（這和CAE中場輸出選擇節點坐標的效果是一致的）。具體如下圖：
2.注意
在ODB結果中創建場輸出時會附帶著一份XYZ坐標，這個應該也可以當做單元的坐標，，但是我比較過這個附帶的坐標和單元的COORD輸出的坐標，有時候有點差別，可能是數據精度的問題。
展開 總結inp中添加關鍵字
輸出單元的積分點坐標：*EL FILE
COORD
輸出節點坐標：*NODE FILE
COORD
原貼出處：https://www.researchgate.net/post/How-to-find-integration-point-coordinates-in-Abaqus-CAE
這是帖子討論的，但是我的嘗試是兩個COORD生成的結果文件是一樣的，都是節點坐標
可以輸出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中平面應力單元高斯積分點的順序為：
展開 
abaqus網格積分的相關專題、標簽、搜索
abaqus網格積分的最新內容
<ul><li>允許<a href="https://zhida.zhihu.com/search?content_id=224701271&content_type=Article&match_order=1&q=Abaqus&zhida_source=entity" rel="noopener noreferrer" target="_blank" style="
參考文獻：《Physically based crystal plasticity FEM including geometrically necessary dislocations: Numerical implementation and applications in micro-forming》
GND 演化方程依賴依賴于剪切應變率的梯度或者塑性變形梯度的旋度，而標準FEM/VUMAT
[圖片]
可以輸出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
插件介紹
Random Agg ITZ Pore 3D (Mesh) V1.0 - AbyssFish 插件可在Abaqus內參數化建立包含水泥漿基體、粗細骨料、界面過渡區（ITZ）、孔隙在內的多相材料混凝土細觀背景網格模型。
模型說明
插件采用材料映射單元的方式，將不同相材料賦值到網格單元，實現三維混凝土細觀有限元模型。
1. 方法
在ABAQUS CAE的場輸出中選擇的坐標點是節點的坐標，而節點是從積分點插值出來的，單元積分點的信息相對真實。所以最好是獲取積分點的信息，其中積分點的坐標無法在CAE中獲取，需要在關鍵字中添加。具體在每個分析步的單元輸出下面添加COORD，如果需要輸出節點的坐標也可以在節點場輸出下面添加COORD（這和CAE中場輸出選擇節點坐標的效果是一致的）。具體如下圖：
<p>在基于實際混凝土斷面圖像進行混凝土細觀有限元模型重建研究方面，主要可采用兩種方式實現：一是根據圖像數據建立實體模型；另一種是采用材料映射單元的方式將不同組分建立背景網格。</p><p><span style="color: rgb(25, 27, 31);">本文將基于以上兩種方式，通過混凝土切片圖片建立二維混凝土細觀有限元模型，并對模型進行軸壓模擬分析。</span></p><div contenteditable
<p class="ql-align-justify">內容記錄帖子，不包含課程內容：請勿購買！</p><p class="ql-align-justify">關于SHPB數值模擬的研究已較為深入，模擬優勢主要在于可通過修正參數使模擬結果與實際一致，以此為基礎對材料的動態破壞過程及更為復雜的工況進行模擬研究，主要研究對象主要分為混凝土、巖石、金屬、陶瓷等材料，并通過<a href="https://
來源：
虛擬Abaqus仿真現實世界
編輯：心印玅經
大部分有限元工程師更愿意花費更多的時間劃分六面體網格，可見六面體網格在分析時是有優勢的，本文分享支架導入的方式對獲取六面體網格的影響，其他較復雜模型可能也同樣適用，如果你學會了，又剛好適合你的模型，那將為你省去很多的時間。
關于該方法，是我在最近仿真冠脈支架時發現的，我使用了不同的3種外觀的支架都是可以滿足使用的
