abaqus 減縮積分
關(guān)注創(chuàng)建者：王靖雯 創(chuàng)建時(shí)間：2023-02-27
abaqus 減縮積分的視頻教程
選擇積分與 Abaqus 梁?jiǎn)卧獌?nèi)核
首先討論選擇積分與降階積分的理論爭(zhēng)議（Bathe vs 王勖成），分析零能模式與剪切鎖死的數(shù)學(xué)機(jī)制。隨后通過(guò) Timoshenko 梁經(jīng)典例題，對(duì)比經(jīng)典梁、精確積分與縮減積分三種結(jié)果，解釋 25% 誤差來(lái)源。接著深入 Abaqus 梁?jiǎn)卧碚摚榻B中心線描述、變形梯度分解、四元數(shù)大轉(zhuǎn)動(dòng)更新及虛功方程。最后說(shuō)明普通梁、開(kāi)口薄壁梁與混合梁?jiǎn)卧倪x型邏輯，并引入張量分析基礎(chǔ)。
￥65 52分鐘 1播放
查看
ABAQUS材料斷裂與失效系列 之 圍道積分運(yùn)算
本專題分一下幾部分進(jìn)行講解： 1.1 靜態(tài)裂縫的兩種模擬方法（達(dá)索09教程） 1.2 圍道積分運(yùn)算（達(dá)索09教程） 1.3 圍道積分相關(guān)的版本更新 1.4 示例1：無(wú)限空間內(nèi)圓形裂縫的圍道積分運(yùn)算的五種方法五種方法： 1.4.1 軸對(duì)稱尖銳裂縫模型 1.4.2 軸對(duì)稱鈍角裂縫模型 1.4.3 三維尖銳裂縫模型 1.4.4 擴(kuò)展有限元（XFEM）方法 1.4.5 子模型方法
￥50 3小時(shí)15分鐘 3034播放
查看
ABAQUS中利用圍線積分法輸出裂紋尖端積分值
基于ABAQUS的官方幫助文檔，詳細(xì)講解了ABAQUS中裂紋體的創(chuàng)建方法+裂紋尖端奇異性單元的構(gòu)件方法（多種單元類型計(jì)算精度對(duì)比）+裂尖積分值的輸出設(shè)置
￥6 34分鐘 167播放
查看
abaqus 減縮積分的實(shí)例教程
公式排版、代碼排版效果不佳，所以上傳的圖片，見(jiàn)諒
2階8節(jié)點(diǎn)減縮積分平面應(yīng)變單元子程序UELMAT源代碼及計(jì)算算例
如果使用了減縮積分的線性單元，即使不是在純彎曲加載模式下，其得到的應(yīng)力應(yīng)變值相比理論預(yù)示值應(yīng)該要?。ㄎ彝茰y(cè)的^_^，沒(méi)空詳細(xì)證實(shí)），所以用這樣的數(shù)據(jù)構(gòu)造的切線剛度矩陣相比其他單元構(gòu)造的切線剛度矩陣要小，這也許就是通常所說(shuō)的出現(xiàn)沙漏問(wèn)題的單元“太軟”的緣故。
結(jié)語(yǔ)：本文算不得什么，只是從公式上加深了商業(yè)軟件使用者對(duì)沙漏這一現(xiàn)象的了解，稍微知其所以然罷了。如果要進(jìn)一步探究如何防止沙漏問(wèn)題，要構(gòu)造怎樣的位移模式，需要更多功夫，可見(jiàn)如下博文：
FEMer，公眾號(hào)：易木木響叮當(dāng)
減縮積分單元、沙漏控制與自定義單元：與Abaqus C3D8R單元的精度對(duì)比之旅
參考資料：
《有限元分析及應(yīng)用》曾攀，清華大學(xué)出版社，2004.
《基于ABAQUS的有限元分析和應(yīng)用》莊茁等，清華大學(xué)出版社2008.
展開(kāi) 之前算過(guò)一個(gè)關(guān)于裂紋擴(kuò)展的問(wèn)題，當(dāng)時(shí)創(chuàng)建裂紋選擇的是contour intergral，后來(lái)又有人咨詢我裂紋尖端J積分的計(jì)算問(wèn)題。我才恍然大悟，其實(shí)圍道積分方法還是適用于計(jì)算裂紋尖端在某時(shí)刻的J積分，至于動(dòng)態(tài)擴(kuò)展問(wèn)題，還是交給XFEM吧（雖然也不太好）。
計(jì)算了幾種情況下的裂紋尖端J積分，包括直裂紋、斜裂紋以及裂紋尖端傾斜等三種情況。
部分試件的應(yīng)力分布及J積分結(jié)果如圖所示：
方法
在ABAQUS CAE的場(chǎng)輸出中選擇的坐標(biāo)點(diǎn)是節(jié)點(diǎn)的坐標(biāo)，而節(jié)點(diǎn)是從積分點(diǎn)插值出來(lái)的，單元積分點(diǎn)的信息相對(duì)真實(shí)。所以最好是獲取積分點(diǎn)的信息，其中積分點(diǎn)的坐標(biāo)無(wú)法在CAE中獲取，需要在關(guān)鍵字中添加。具體在每個(gè)分析步的單元輸出下面添加COORD，如果需要輸出節(jié)點(diǎn)的坐標(biāo)也可以在節(jié)點(diǎn)場(chǎng)輸出下面添加COORD（這和CAE中場(chǎng)輸出選擇節(jié)點(diǎn)坐標(biāo)的效果是一致的）。具體如下圖：
2.注意
在ODB結(jié)果中創(chuàng)建場(chǎng)輸出時(shí)會(huì)附帶著一份XYZ坐標(biāo)，這個(gè)應(yīng)該也可以當(dāng)做單元的坐標(biāo)，，但是我比較過(guò)這個(gè)附帶的坐標(biāo)和單元的COORD輸出的坐標(biāo)，有時(shí)候有點(diǎn)差別，可能是數(shù)據(jù)精度的問(wèn)題。
展開(kāi) 
abaqus 減縮積分的相關(guān)專題、標(biāo)簽、搜索
abaqus 減縮積分的最新內(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
個(gè)人筆記，錯(cuò)誤難免，懇請(qǐng)指出，共同進(jìn)步。
參考資料見(jiàn)文后，文中的引用以“作者+頁(yè)碼”、“作者名年份+頁(yè)碼”等方式呈現(xiàn)。
引言：
莊茁P65對(duì)沙漏現(xiàn)象的描述如下圖：
本文試圖基于純彎曲加載下線性減縮積分的應(yīng)變公式，對(duì)沙漏現(xiàn)象的產(chǎn)生機(jī)理進(jìn)行淺淺的理論闡述。
我們?cè)谇耙黄┪闹泻?jiǎn)述了有限元中的數(shù)值積分機(jī)理：
數(shù)峰青，公眾號(hào)：數(shù)峰青
有限元筆記#1：什么是剪切自鎖
1. 方法
在ABAQUS CAE的場(chǎng)輸出中選擇的坐標(biāo)點(diǎn)是節(jié)點(diǎn)的坐標(biāo)，而節(jié)點(diǎn)是從積分點(diǎn)插值出來(lái)的，單元積分點(diǎn)的信息相對(duì)真實(shí)。所以最好是獲取積分點(diǎn)的信息，其中積分點(diǎn)的坐標(biāo)無(wú)法在CAE中獲取，需要在關(guān)鍵字中添加。具體在每個(gè)分析步的單元輸出下面添加COORD，如果需要輸出節(jié)點(diǎn)的坐標(biāo)也可以在節(jié)點(diǎn)場(chǎng)輸出下面添加COORD（這和CAE中場(chǎng)輸出選擇節(jié)點(diǎn)坐標(biāo)的效果是一致的）。具體如下圖：
公式排版、代碼排版效果不佳，所以上傳的圖片，見(jiàn)諒
按照正常的理解，毫無(wú).疑問(wèn)，abaqus 全積分一定是采用了2x2x2=8個(gè)積分點(diǎn)。
從后處理結(jié)果來(lái)看，似乎也是如此，每個(gè)單元存在8個(gè)積分點(diǎn)。
然而，如果自己動(dòng)手跑一遍程序，就會(huì)發(fā)現(xiàn)事實(shí)遠(yuǎn)非如此，采用全積分計(jì)算得到的結(jié)果與abaqus 存在差異，原因何在？
事實(shí)賞，abaqus C3D8 采用的選擇積分方式（selective intergation schema），即對(duì)于偏應(yīng)變，采用
2階8節(jié)點(diǎn)減縮積分平面應(yīng)變單元子程序UELMAT源代碼及計(jì)算算例
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)力.但是提取出來(lái)的應(yīng)力為單元積分點(diǎn)上的應(yīng)力.無(wú)法獲取節(jié)點(diǎn)上的應(yīng)力.同時(shí)在ABAQUS中的子程序中,也是對(duì)積分點(diǎn)上的數(shù)據(jù)進(jìn)行操作.
本文基于個(gè)人興趣同時(shí)想要更加了解有限元背后原理和公式的想法.近日進(jìn)行了一些初步的探索.希望大家批評(píng)指正
可在ABAQUS中對(duì)減縮積分單元引入少量的人工“沙漏剛度”以限制沙漏模式的擴(kuò)展。當(dāng)模型中有更多的單元時(shí)，這種剛度在限制沙漏模式方面是更有效的，這意味著只要采用合理的細(xì)網(wǎng)格，線性減縮積分單元會(huì)給出可接受的結(jié)果。對(duì)許多應(yīng)用而言，采用細(xì)網(wǎng)格的線性減縮積分單元所產(chǎn)生的誤差是在一個(gè)可接受的范圍內(nèi)的。這個(gè)結(jié)果說(shuō)明當(dāng)用這類單元來(lái)模擬承受彎曲載荷的結(jié)構(gòu)時(shí)，在厚度方向上至少應(yīng)采用四個(gè)單元。
