abaqus 厚度積分的案例
ABAQUS直裂紋、斜裂紋圍道積分計算裂紋尖端J積分
之前算過一個關于裂紋擴展的問題，當時創建裂紋選擇的是contour intergral，后來又有人咨詢我裂紋尖端J積分的計算問題。我才恍然大悟，其實圍道積分方法還是適用于計算裂紋尖端在某時刻的J積分，至于動態擴展問題，還是交給XFEM吧（雖然也不太好）。
計算了幾種情況下的裂紋尖端J積分，包括直裂紋、斜裂紋以及裂紋尖端傾斜等三種情況。
部分試件的應力分布及J積分結果如圖所示：
HyperWorks-Abaqus二次開發之自動抽中面+自動分組+識別厚度+賦予厚度屬性
[圖片]
abaqus顯式接觸厚度問題
我做了一個落錘沖擊梁的模型，落錘用的離散剛體，加工成殼，用動力隱式可以計算，但用動力顯式就報錯。
動力隱式容易不收斂，所以現在不知道在顯式中怎么改。
For two-dimensional models, if a non-unity thickness is specified for two-dimensional solid elements and these elements are involved in an interaction such as contact, the same thickness should be specified for the out-of-plane thickness of the corresponding surface under *surface interaction.
展開 ABAQUS輸出單元積分點坐標
方法
在ABAQUS CAE的場輸出中選擇的坐標點是節點的坐標，而節點是從積分點插值出來的，單元積分點的信息相對真實。所以最好是獲取積分點的信息，其中積分點的坐標無法在CAE中獲取，需要在關鍵字中添加。具體在每個分析步的單元輸出下面添加COORD，如果需要輸出節點的坐標也可以在節點場輸出下面添加COORD（這和CAE中場輸出選擇節點坐標的效果是一致的）。具體如下圖：
2.注意
在ODB結果中創建場輸出時會附帶著一份XYZ坐標，這個應該也可以當做單元的坐標，，但是我比較過這個附帶的坐標和單元的COORD輸出的坐標，有時候有點差別，可能是數據精度的問題。
展開 
ABAQUS 輸出節點坐標和積分點坐標
總結inp中添加關鍵字
輸出單元的積分點坐標：*EL FILE
COORD
輸出節點坐標：*NODE FILE
COORD
原貼出處：https://www.researchgate.net/post/How-to-find-integration-point-coordinates-in-Abaqus-CAE
這是帖子討論的，但是我的嘗試是兩個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中獲取積分點坐標的三種方法
經常有小伙伴問獲取積分點坐標的方法，今天給大家介紹三種獲取積分點坐標的方式，希望能給你們帶來幫助。
1 通過abaqus子程序獲取積分點坐標
Abaqus一些子程序中可以直接獲取積分點坐標，例如我們熟知的UMAT子程序中包含COORD參數，即為積分點坐標。順帶一提的是，當打開了幾何非線性時，該積分點是當前構形下的坐標，如果未打開幾何非線性則為初始坐標。
2通過history output輸出積分點坐標
Abaqus可以直接在歷程變量history output中輸出積分點坐標。直接在history output中勾選COORD選項，但是這里需要注意的是，Domain中的Set集合如果是node set，這里輸出來的是節點坐標，當這里是element set的時候，輸出來的才是積分點坐標。
3通過等參單元映射函數計算
等參元中，為了方便計算，把整體坐標映射到自然坐標，然后在自然坐標下進行高斯積分。如果知道了自然坐標下的高斯積分點，通過映射函數反算，便能得到整體坐標下的高斯積分點坐標。以四邊形等參單元為例，其以自然呢坐標為變量的插值形函數如下
坐標變換采取同樣的插值函數(叫做等參的原因)，整體坐標和自然坐標的關系式如下，如果知道自然坐標下的高斯積分點，直接通過此公式計算其在整體坐標下的坐標。
展開 批量0厚度cohesive單元的Abaqus插件（極速版） ￥300
新功能要求Abaqus2018及以上，其他功能要求Abaqus6.14及以上。</span>
</div>
</div><div contenteditable="false" width="100%">
<div data-offset-key="9a73f-0-0" class="public-DraftStyleDefault-block public-DraftStyleDefault-ltr" style="max-width: 100%;">
<span data-offset-key="9a73f-0-0">全局插入有厚度單元目前似乎無應用場景，不打算增加這個功能。新生成的有厚度單元從Matrix中劃分出來，非Matrix區域的節點和單元保持不變 (如果remesh 編號可能會發生變化）。可以選擇對分割出界面單元后的matrix區域重新劃分網格（remesh），remesh功能要求Matrix區域的網格必須全是三角形網格。
展開 Abaqus創建零厚度cohesive單元
有些情況下需要建立零厚度的cohesive單元以進行裂紋擴展的計算模擬，下面以第一種方法進行零厚度的cohesive單元的建模。
新建一個二維模型，如下圖所示：
假設在模型的中間有一層零厚度的cohesive單元層，做一個切割處理，如下圖所示：
切割出一個很小的縫隙之后進行網格劃分，如下所示：
完成這一步之后需要通過菜單欄Mesh-creat mesh part將單元變成孤立網格。
在此基礎上即可創建零厚度的單元。點擊菜單欄Mesh下面的Edit，如下所示：
然后框選出中間縫隙的上面一層節點，之后再選擇下面一層的一個單元面，將縫隙的上面的所有節點投影到縫隙的下層。
投影完的效果如下，中間的縫隙不見了。
打開節點號的顯示，如下所示，可以發現在中間的一層位置處相同位置有兩個節點，也即是該層單元為零厚度單元
再之后，通過Assign Element Type可修改中間這一層重合的單元，修改為cohesive單元。
長安CAE的博客
http://blog.sina.com.cn/zuoerninan
展開 abaqus模擬研究不同強度的再生磚混凝土和鋁管厚度軸壓性能的差異 ￥9.9
1.1項目概況
該課題研究不同強度的再生磚混凝土和鋁管厚度軸壓性能的差異。
1.2項目要求
以上述參數進行有限元分析，并提取其荷載-縱向應變關系與試驗數據進行比。
1.3單位制
在CAE項目計算以及報告中使用的基本單位系統如表格 01所示。
表格 11單位系統
序號
單位類型
單位名稱
符號
1
長度
毫米
mm
2
時間
秒
s
3
質量
噸
T
4
溫度
攝氏度
℃
1.4軟件介紹
本項目在實施過程中主要采用的CAE軟件包括： Abaqus.
(1) 有限元求解軟件采用ABAQUS2020
ABAQUS是一套功能強大的工程模擬的有限元軟件，其解決問題的范圍從相對簡單的線性分析到許多復雜的非線性問題。ABAQUS包括一個豐富的、可模擬任意幾何形狀的單元庫。
展開 批量0厚度cohesive單元的Abaqus插件（終極版） ￥300
（3）提交運行，生成新的Part：Cohesive，生成后的新part和單元集合如下：
新part中的cohesive單元如下：
插件試用及購買
百度網盤下載插件，解壓后放在Abaqus插件目錄中，即可使用，具體Abaqus插件的安裝方法可以百度，非常簡單。插件未激活前，只能用于單元數量不超過2000的part， 激活后的插件，單元數量和使用時間沒有限制。如果使用過程中遇到buge，歡迎反饋問題，提交過buge的用戶購買有優惠。
適用Abaqus6.14及以上版本。
使用過程遇到問題，聯系QQ：1913905197，或本貼留言。
知乎下載鏈接
https://zhuanlan.zhihu.com/p/401803056
展開 
Abaqus降魔篇之圍道積分小探
1 分析目的
平面應變的一塊板，在扭矩的作用下，分析裂紋端的應力強度因子以及斷裂能
模擬采用CPE4R（使用CPE8R往往會得到更精確的結果）
2 模型建立
然后進行切分，主要是切出預制裂紋面以及圍道外邊
在這里注意切那個圓的時候一定要將半徑點選到裂紋上，以便劃分網格
這里的圓半徑為0.5mm
3 建立材料以及相關屬性
4 建立模型，預制裂紋
選擇之前的切分線
選擇裂紋端點
這里的0.25是根據劃分網格的種子點確定的（4個點）
而collapsed的選項是要求對重疊單元的裂紋控制只對單一點有效
5 設置邊界條件
a首先是施加扭矩端的耦合
b然后是扭矩施加
另一端同理可得
c 限制位移
6 計算吧！
然后你可以看到應力強度因子以及斷裂能
或者你比較喜歡看趨勢
在或者你對歷史輸出里面的5感到困惑
那就看下圖
展開 abaqus中平面應力應變厚度對切削力的影響 ￥5
在鋁合金的二位正交切削仿真中，不同的平面應力應變厚度的對切削力的影響結果
以上為設定值為1的情況
abaqus怎么獲取單元積分點上的坐標呢
如題 為什么我查詢的時候坐標都是0
abaqus C3D8 單元 計算中采用了多少個積分點？
按照正常的理解，毫無.疑問，abaqus 全積分一定是采用了2x2x2=8個積分點。
從后處理結果來看，似乎也是如此，每個單元存在8個積分點。
然而，如果自己動手跑一遍程序，就會發現事實遠非如此，采用全積分計算得到的結果與abaqus 存在差異，原因何在？
事實賞，abaqus C3D8 采用的選擇積分方式（selective intergation schema），即對于偏應變，采用8個積分，對于球應變，采用中心點積分。這樣計算得到的結果才能與abaqus 完全對標，亦可從abaqus 幫助文檔得到答案。
展開 