abaqus顯示積分點(diǎn)
關(guān)注創(chuàng)建者：王靖雯 創(chuàng)建時(shí)間：2023-02-27
abaqus顯示積分點(diǎn)的實(shí)例教程
方法
在ABAQUS CAE的場輸出中選擇的坐標(biāo)點(diǎn)是節(jié)點(diǎn)的坐標(biāo)，而節(jié)點(diǎn)是從積分點(diǎn)插值出來的，單元積分點(diǎn)的信息相對真實(shí)。所以最好是獲取積分點(diǎn)的信息，其中積分點(diǎn)的坐標(biāo)無法在CAE中獲取，需要在關(guān)鍵字中添加。具體在每個(gè)分析步的單元輸出下面添加COORD，如果需要輸出節(jié)點(diǎn)的坐標(biāo)也可以在節(jié)點(diǎn)場輸出下面添加COORD（這和CAE中場輸出選擇節(jié)點(diǎn)坐標(biāo)的效果是一致的）。具體如下圖：
2.注意
在ODB結(jié)果中創(chuàng)建場輸出時(shí)會(huì)附帶著一份XYZ坐標(biāo)，這個(gè)應(yīng)該也可以當(dāng)做單元的坐標(biāo)，，但是我比較過這個(gè)附帶的坐標(biāo)和單元的COORD輸出的坐標(biāo)，有時(shí)候有點(diǎn)差別，可能是數(shù)據(jù)精度的問題。
展開 總結(jié)inp中添加關(guān)鍵字
輸出單元的積分點(diǎn)坐標(biāo)：*EL FILE
COORD
輸出節(jié)點(diǎn)坐標(biāo)：*NODE FILE
COORD
原貼出處：https://www.researchgate.net/post/How-to-find-integration-point-coordinates-in-Abaqus-CAE
這是帖子討論的，但是我的嘗試是兩個(gè)COORD生成的結(jié)果文件是一樣的，都是節(jié)點(diǎn)坐標(biāo)
可以輸出umat接口中的變量coords進(jìn)行查看
write(*,"(A,I4)") "npt = ", npt
write(*,"(A,3ES16.8)") "coords = ", coords
結(jié)果為：
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中平面應(yīng)力單元高斯積分點(diǎn)的順序?yàn)椋?/span>
展開 經(jīng)常有小伙伴問獲取積分點(diǎn)坐標(biāo)的方法，今天給大家介紹三種獲取積分點(diǎn)坐標(biāo)的方式，希望能給你們帶來幫助。
1 通過abaqus子程序獲取積分點(diǎn)坐標(biāo)
Abaqus一些子程序中可以直接獲取積分點(diǎn)坐標(biāo)，例如我們熟知的UMAT子程序中包含COORD參數(shù)，即為積分點(diǎn)坐標(biāo)。順帶一提的是，當(dāng)打開了幾何非線性時(shí)，該積分點(diǎn)是當(dāng)前構(gòu)形下的坐標(biāo)，如果未打開幾何非線性則為初始坐標(biāo)。
2通過history output輸出積分點(diǎn)坐標(biāo)
Abaqus可以直接在歷程變量history output中輸出積分點(diǎn)坐標(biāo)。直接在history output中勾選COORD選項(xiàng)，但是這里需要注意的是，Domain中的Set集合如果是node set，這里輸出來的是節(jié)點(diǎn)坐標(biāo)，當(dāng)這里是element set的時(shí)候，輸出來的才是積分點(diǎn)坐標(biāo)。
3通過等參單元映射函數(shù)計(jì)算
等參元中，為了方便計(jì)算，把整體坐標(biāo)映射到自然坐標(biāo)，然后在自然坐標(biāo)下進(jìn)行高斯積分。如果知道了自然坐標(biāo)下的高斯積分點(diǎn)，通過映射函數(shù)反算，便能得到整體坐標(biāo)下的高斯積分點(diǎn)坐標(biāo)。以四邊形等參單元為例，其以自然呢坐標(biāo)為變量的插值形函數(shù)如下
坐標(biāo)變換采取同樣的插值函數(shù)(叫做等參的原因)，整體坐標(biāo)和自然坐標(biāo)的關(guān)系式如下，如果知道自然坐標(biāo)下的高斯積分點(diǎn)，直接通過此公式計(jì)算其在整體坐標(biāo)下的坐標(biāo)。
展開 如題 為什么我查詢的時(shí)候坐標(biāo)都是0
abaqus顯示積分點(diǎn)的相關(guān)專題、標(biāo)簽、搜索
abaqus顯示積分點(diǎn)的最新內(nèi)容
Abaqus中平面應(yīng)力單元高斯積分點(diǎn)的順序11個(gè)月前
可以輸出umat接口中的變量coords進(jìn)行查看
write(*,"(A,I4)") "npt = ", npt
write(*,"(A,3ES16.8)") "coords = ", coords
結(jié)果為：
npt = 1
coords = -5.77350269E-01 -5.77350269E-01 1.00000000E-02
npt = 2
1. 方法
在ABAQUS CAE的場輸出中選擇的坐標(biāo)點(diǎn)是節(jié)點(diǎn)的坐標(biāo)，而節(jié)點(diǎn)是從積分點(diǎn)插值出來的，單元積分點(diǎn)的信息相對真實(shí)。所以最好是獲取積分點(diǎn)的信息，其中積分點(diǎn)的坐標(biāo)無法在CAE中獲取，需要在關(guān)鍵字中添加。具體在每個(gè)分析步的單元輸出下面添加COORD，如果需要輸出節(jié)點(diǎn)的坐標(biāo)也可以在節(jié)點(diǎn)場輸出下面添加COORD（這和CAE中場輸出選擇節(jié)點(diǎn)坐標(biāo)的效果是一致的）。具體如下圖：
按照正常的理解，毫無.疑問，abaqus 全積分一定是采用了2x2x2=8個(gè)積分點(diǎn)。
從后處理結(jié)果來看，似乎也是如此，每個(gè)單元存在8個(gè)積分點(diǎn)。
然而，如果自己動(dòng)手跑一遍程序，就會(huì)發(fā)現(xiàn)事實(shí)遠(yuǎn)非如此，采用全積分計(jì)算得到的結(jié)果與abaqus 存在差異，原因何在？
事實(shí)賞，abaqus C3D8 采用的選擇積分方式（selective intergation schema），即對于偏應(yīng)變，采用
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
如題 為什么我查詢的時(shí)候坐標(biāo)都是0
在ABAQUS中,當(dāng)需要獲取節(jié)點(diǎn)上的應(yīng)力時(shí),可以在后處理中建立路徑或者用查詢功能等獲取.
但是當(dāng)需要大量的節(jié)點(diǎn)上應(yīng)力數(shù)據(jù)時(shí),很多人會(huì)用Python編程進(jìn)行大批量的提取應(yīng)力.但是提取出來的應(yīng)力為單元積分點(diǎn)上的應(yīng)力.無法獲取節(jié)點(diǎn)上的應(yīng)力.同時(shí)在ABAQUS中的子程序中,也是對積分點(diǎn)上的數(shù)據(jù)進(jìn)行操作.
本文基于個(gè)人興趣同時(shí)想要更加了解有限元背后原理和公式的想法.近日進(jìn)行了一些初步的探索.希望大家批評指正
經(jīng)常有小伙伴問獲取積分點(diǎn)坐標(biāo)的方法，今天給大家介紹三種獲取積分點(diǎn)坐標(biāo)的方式，希望能給你們帶來幫助。
1 通過abaqus子程序獲取積分點(diǎn)坐標(biāo)
Abaqus一些子程序中可以直接獲取積分點(diǎn)坐標(biāo)，例如我們熟知的UMAT子程序中包含COORD參數(shù)，即為積分點(diǎn)坐標(biāo)。順帶一提的是，當(dāng)打開了幾何非線性時(shí)，該積分點(diǎn)是當(dāng)前構(gòu)形下的坐標(biāo)，如果未打開幾何非線性則為初始坐標(biāo)。
