abaqus網格積分點
關注創建者：王靖雯 創建時間：2023-02-27
abaqus網格積分點的視頻教程
ABAQUS-ECC和FRP網格加固混凝土梁四點彎曲CBM論文復現
ABAQUS-ECC和FRP網格加固混凝土梁四點彎曲CBM論文復現。兩個混凝土梁受彎試件荷載-位移曲線如何吻合。保姆及教程，1.3小時，從建模開始錄制。
￥240 1小時12分鐘 82播放
查看
abaqus腳本插件104-基于空間點云坐標重構六面體網格節點（2025-10-01）
abaqus腳本插件104-基于空間點云坐標重構六面體網格節點（2025-10-01）
免費 12分鐘 5播放
查看
abaqus腳本插件104-基于空間點云坐標重構六面體網格節點（2025-09-24）
abaqus腳本插件104-基于空間點云坐標重構六面體網格節點（2025-09-24）
免費 5分鐘 3播放
查看
abaqus網格積分點的實例教程
參考文獻：《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中平面應力單元高斯積分點的順序為：
展開 經常有小伙伴問獲取積分點坐標的方法，今天給大家介紹三種獲取積分點坐標的方式，希望能給你們帶來幫助。
1 通過abaqus子程序獲取積分點坐標
Abaqus一些子程序中可以直接獲取積分點坐標，例如我們熟知的UMAT子程序中包含COORD參數，即為積分點坐標。順帶一提的是，當打開了幾何非線性時，該積分點是當前構形下的坐標，如果未打開幾何非線性則為初始坐標。
2通過history output輸出積分點坐標
Abaqus可以直接在歷程變量history output中輸出積分點坐標。直接在history output中勾選COORD選項，但是這里需要注意的是，Domain中的Set集合如果是node set，這里輸出來的是節點坐標，當這里是element set的時候，輸出來的才是積分點坐標。
3通過等參單元映射函數計算
等參元中，為了方便計算，把整體坐標映射到自然坐標，然后在自然坐標下進行高斯積分。如果知道了自然坐標下的高斯積分點，通過映射函數反算，便能得到整體坐標下的高斯積分點坐標。以四邊形等參單元為例，其以自然呢坐標為變量的插值形函數如下
坐標變換采取同樣的插值函數(叫做等參的原因)，整體坐標和自然坐標的關系式如下，如果知道自然坐標下的高斯積分點，直接通過此公式計算其在整體坐標下的坐標。
展開 
abaqus網格積分點的相關專題、標簽、搜索
abaqus網格積分點的最新內容
參考文獻：《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
1. 方法
在ABAQUS CAE的場輸出中選擇的坐標點是節點的坐標，而節點是從積分點插值出來的，單元積分點的信息相對真實。所以最好是獲取積分點的信息，其中積分點的坐標無法在CAE中獲取，需要在關鍵字中添加。具體在每個分析步的單元輸出下面添加COORD，如果需要輸出節點的坐標也可以在節點場輸出下面添加COORD（這和CAE中場輸出選擇節點坐標的效果是一致的）。具體如下圖：
按照正常的理解，毫無.疑問，abaqus 全積分一定是采用了2x2x2=8個積分點。
從后處理結果來看，似乎也是如此，每個單元存在8個積分點。
然而，如果自己動手跑一遍程序，就會發現事實遠非如此，采用全積分計算得到的結果與abaqus 存在差異，原因何在？
事實賞，abaqus C3D8 采用的選擇積分方式（selective intergation schema），即對于偏應變，采用
In CAE, you can request coordinate outputs by going to the Field Output Request->Edit->Volume/Thickness/Coordinates->COORD, Current nodal coordinates. This way, you will be able to obtain nodal coordinates
如題 為什么我查詢的時候坐標都是0
在ABAQUS中,當需要獲取節點上的應力時,可以在后處理中建立路徑或者用查詢功能等獲取.
但是當需要大量的節點上應力數據時,很多人會用Python編程進行大批量的提取應力.但是提取出來的應力為單元積分點上的應力.無法獲取節點上的應力.同時在ABAQUS中的子程序中,也是對積分點上的數據進行操作.
本文基于個人興趣同時想要更加了解有限元背后原理和公式的想法.近日進行了一些初步的探索.希望大家批評指正
經常有小伙伴問獲取積分點坐標的方法，今天給大家介紹三種獲取積分點坐標的方式，希望能給你們帶來幫助。
1 通過abaqus子程序獲取積分點坐標
Abaqus一些子程序中可以直接獲取積分點坐標，例如我們熟知的UMAT子程序中包含COORD參數，即為積分點坐標。順帶一提的是，當打開了幾何非線性時，該積分點是當前構形下的坐標，如果未打開幾何非線性則為初始坐標。
