abaqus單元積分點的案例
ABAQUS輸出單元積分點坐標
方法
在ABAQUS CAE的場輸出中選擇的坐標點是節點的坐標，而節點是從積分點插值出來的，單元積分點的信息相對真實。所以最好是獲取積分點的信息，其中積分點的坐標無法在CAE中獲取，需要在關鍵字中添加。具體在每個分析步的單元輸出下面添加COORD，如果需要輸出節點的坐標也可以在節點場輸出下面添加COORD（這和CAE中場輸出選擇節點坐標的效果是一致的）。具體如下圖：
2.注意
在ODB結果中創建場輸出時會附帶著一份XYZ坐標，這個應該也可以當做單元的坐標，，但是我比較過這個附帶的坐標和單元的COORD輸出的坐標，有時候有點差別，可能是數據精度的問題。
展開 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中平面應力單元高斯積分點的順序為：
展開 abaqus怎么獲取單元積分點上的坐標呢
如題 為什么我查詢的時候坐標都是0
abaqus C3D8 單元 計算中采用了多少個積分點？
按照正常的理解，毫無.疑問，abaqus 全積分一定是采用了2x2x2=8個積分點。
從后處理結果來看，似乎也是如此，每個單元存在8個積分點。
然而，如果自己動手跑一遍程序，就會發現事實遠非如此，采用全積分計算得到的結果與abaqus 存在差異，原因何在？
事實賞，abaqus C3D8 采用的選擇積分方式（selective intergation schema），即對于偏應變，采用8個積分，對于球應變，采用中心點積分。這樣計算得到的結果才能與abaqus 完全對標，亦可從abaqus 幫助文檔得到答案。
展開 
有限元中單元積分點與節點應力相互轉換(CPE4為例)（ABAQUS）
在ABAQUS中,當需要獲取節點上的應力時,可以在后處理中建立路徑或者用查詢功能等獲取.
但是當需要大量的節點上應力數據時,很多人會用Python編程進行大批量的提取應力.但是提取出來的應力為單元積分點上的應力.無法獲取節點上的應力.同時在ABAQUS中的子程序中,也是對積分點上的數據進行操作.
本文基于個人興趣同時想要更加了解有限元背后原理和公式的想法.近日進行了一些初步的探索.希望大家批評指正. 本文基本不涉及原理公式，只在轉換積分點和節點的應力時列出公式。盡可能簡介易懂。
一: 單元類型及節點數目與位移,應變,應力階次的關系
本節內容基于有限元教材及一些網上資料.
(1)有限元求解的思路是:
一: 建立單元節點力與節點位移關系式.
二: 將彈性體上的外載荷等效移置到節點上.
三: 在節點上建立力的平衡方程,求得節點位移.
四: 通過彈性力學基本方程,可求得單元的應力和應變.
(2) 四節點矩形單元
以四節點矩形單元為例，在此只表達有限元教材中的結論，具體公式可參考有限元教材。
（3）ABAQUS中的CPE4單元
CPE4: A 4-node bilinear plane strain quadrilateral.
該單元有四個節點，同時有四個積分點。
對于每個應力分量（注意：在此只看一個應力分量），單元內任一點（x,y）的應力表達式為：
stress=a*x*y+b*x+c*y+d （1）
該表達式有四個未知量：a,b,c,d。
若知道四個積分點的應力分量。將每個積分點帶入上式，則會形成包含四個方程的線性方程組。
展開 關于有限單元法中節（結）點與積分點的幾點釋疑
節點和積分點是有限單元法（FEM）的兩個基本概念，初涉有限元計算的同志往往在這點上產生混淆，假設導師面試的時候，問單元應力是什么，若回答不慎，將貽笑大方，得不償失。本文試圖以簡略易懂的說法來闡述節點和積分點的區別。
1.節點位移是有限元法的基本未知量。節點構筑了問題域的幾何離散化形狀，節點是形函數的零點，通常形函數是以節點為依據進行假設的。形函數決定了單元內部各點運動的位移模式（常用帕斯卡三角形來選擇單元位移模式），這樣就形成了數學上所說的插值。
有限元法的原理就是將問題域分割成N多小單元，在每個單元內采用簡單的函數來近似表達單元的真實位移，將各單元再連接起來，就可以近似描述整個問題域的運動。因此，有限元法從根本上就是精確的，而不是準確的。
2.積分點是單元進行數值積分的已知量。有限元法中一般采用高斯積分，但是積分方法不限于高斯積分，如果有人用了Irons積分或者Hammer積分，請不要驚訝。在形成單元剛度矩陣和進行節點應力磨平的時候，需要高斯積分。
以等參單元為例，其剛度矩陣
，這個就需要數值積分來快速計算，高斯點坐標及權系數如表4.2[王勖成]所示。
老師授課時一般對常應力單元進行推導，而常應力單元只有一個積分點，被積函數是常數，因此體現不出高斯積分來。很多老師對高斯積分在單元剛度矩陣的應用不予細述，導致部分同學對單元積分點認識不足。
3.單元應力指的是高斯積分點的應力，而非節點上的應力。有了位移模式，再通過虛功原理得到單元剛度矩陣，然后聚合總剛，求解平衡方程，就會把基本未知量——節點位移求出來了。通過節點位移得到單元應變結果，利用物理方程求得單元應力結果。
在等參元中，單元中n+1階(n=p-m)高斯積分點上的應變或應力近似解比其它部位具有較高的精度，因此我們稱(n+1)階高斯積分點是等參元中的最佳應力點。
展開 單元積分點應力如何外插至節點上 | 數值實現篇
繼上次的推文：有限元計算過程中積分點應力如何外插至節點處？【公式推導篇】，本次分享單元積分點應力外插至節點處的數值實現過程。
數值實現
借助以上理論，我們可以基于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位移應力場結果
自研程序應力場結果
結果是一致的，說明了程序的正確性。
展開 ABAQUS 輸出節點坐標和積分點坐標
總結inp中添加關鍵字
輸出單元的積分點坐標：*EL FILE
COORD
輸出節點坐標：*NODE FILE
COORD
原貼出處：https://www.researchgate.net/post/How-to-find-integration-point-coordinates-in-Abaqus-CAE
這是帖子討論的，但是我的嘗試是兩個COORD生成的結果文件是一樣的，都是節點坐標
展開 Abaqus中獲取積分點坐標的三種方法
經常有小伙伴問獲取積分點坐標的方法，今天給大家介紹三種獲取積分點坐標的方式，希望能給你們帶來幫助。
1 通過abaqus子程序獲取積分點坐標
Abaqus一些子程序中可以直接獲取積分點坐標，例如我們熟知的UMAT子程序中包含COORD參數，即為積分點坐標。順帶一提的是，當打開了幾何非線性時，該積分點是當前構形下的坐標，如果未打開幾何非線性則為初始坐標。
2通過history output輸出積分點坐標
Abaqus可以直接在歷程變量history output中輸出積分點坐標。直接在history output中勾選COORD選項，但是這里需要注意的是，Domain中的Set集合如果是node set，這里輸出來的是節點坐標，當這里是element set的時候，輸出來的才是積分點坐標。
3通過等參單元映射函數計算
等參元中，為了方便計算，把整體坐標映射到自然坐標，然后在自然坐標下進行高斯積分。如果知道了自然坐標下的高斯積分點，通過映射函數反算，便能得到整體坐標下的高斯積分點坐標。以四邊形等參單元為例，其以自然呢坐標為變量的插值形函數如下
坐標變換采取同樣的插值函數(叫做等參的原因)，整體坐標和自然坐標的關系式如下，如果知道自然坐標下的高斯積分點，直接通過此公式計算其在整體坐標下的坐標。
展開 ABAQUS讀懂用戶手冊系列—修煉Cohesive內功：內聚力單元/接觸基礎知識點
<p>Cohesive作為ABAQUS中常用的粘結技術，無論在模擬粘結界面（例如新舊混凝土疊合面、復合材料粘結界面）或是全局粘結單元（例如模擬細觀混凝土開裂）具有較廣泛的應用。今天喵星人從官方的用戶手冊中選取了幾個Cohesive基礎而又關鍵的知識點，幫助大家修煉Cohesive內功。</p><p class="ql-align-center"><strong>01</strong>內聚力單元/接觸區別</p><p>對于內聚力單元/接觸的相同之處，用戶手冊指出：</p><p>The formulae and laws that govern cohesive constitutive behavior are very similar for cohesive contact and cohesive elements. The similarities extend to the linear elastic traction-separation model, damage initiation criteria, and damage evolution laws.</p><p><strong>喵星人翻譯：</strong></p><p>控制內聚本構行為的公式和法則在內聚力接觸與內聚力單元中極為相似。包括適用于線性彈性牽引-分離模型、損傷起始準則以及損傷演化法則。</p><p>可見，內聚力接觸與內聚力單元的本構模型基本相同，在ABAQUS中的操作差異在于內聚力單元需要在材料的屬性中輸入，內聚力接觸則是在相互作用的接觸中輸入。
展開 一維單元模擬混凝土構件開裂的解決方案
圖1 鋼筋混凝土結構裂縫
本案例使用ABAQUS對一根鋼筋混凝土受拉構件進行裂縫估算分析，根據鋼筋應力狀態計算等效裂縫寬度，并采用Python腳本在ODB結果文件中創建裂縫場變量，實現在ABAQUS中使用一維單元進行快速分析，并在后處理模塊 顯示等效裂縫寬度的目標。
案例涉及的相關技術：
①ABAQUS梁單元Rebar積分點插入；
②利用Python腳本提取ABAQUS場變量數據；
③利用Python腳本創建ABAQUS場變量數據。
計算報告編寫采用操作引導式，希望能為讀者使用ABAUQS場變量創建提供有益參考。操作分析要點為：
①ABAQUSABAQUS梁單元Rebar積分點插入；
②ABAQUS場變量輸出Python腳本getSubset()函數應用；
③ABAQUS場變量編輯Python腳本addData()()函數應用。
二、計算任務
1.模型裝配及接觸連接
計算模型取自《混凝土結構：混凝土結構設計原理》（第六版）習題8-3。
計算模型為鋼筋混凝土屋架下弦按軸心受拉構件，見圖2。模型只包含1個part。截面寬200mm，截面高160mm。因為案例模型較為簡單，混凝土梁采用B21單元模擬，鋼筋通過在與混凝土單元共節點建立鋼筋箱型截面單元實現?；炷敛牧蠟镃40，鋼筋為HRB400，為實現材料單軸滯回損傷模擬，混凝土受壓骨架線上升段采用Hognested曲線，下降段為直線；鋼筋均采用《混凝土結構設計規范》（GB 50010-2010）附錄C中鋼筋本構關系，材料本構示意圖見圖3，材料參數見表1。
展開 
Abaqus疑難雜癥——INP文件自定義輸出
以下內容轉載本人的公眾號：易木木響叮當
本期內容延續上一期內容，繼續詳解INP文件的功能，內容涉及分析結果的輸出以及單元剛度矩陣、質量矩陣的打印。
將分析結果輸出到DAT文件
可在INP文件中加入以下內容：
*STEP,NAME=step-1
......
*NODE PRINT,NSET=node-11
U1
*EL PRINT,ELSET=element-11
S11
*END STEP
注：位移對應節點、應力對應單元積分點，在打印輸出時應注意區別。軟件后處理時可以看到節點上的應力結果是因為ABAQUS對單元積分點上的應力做了差值平均。以下列出常用關鍵詞的用法，大家可根據自己需要學習。
展開 