abaqus混合單元計算的案例
ABAQUS UMAT后歐拉法實現cohesive單元的拉-剪混合本構 ￥1500
本文介紹了如何用后歐拉算法模擬2D cohesive單元的復雜破壞，采用的本構模型是hyperbolic曲線。這種曲線可以將受拉破壞和受剪破壞耦合在一起，很好地反應了膠合接觸面地特性。
主要參考文獻為：Caballero, A., Willam, K.J. and Carol, I., 2008. Consistent tangent formulation for 3D interface modeling of cracking/fracture in quasi-brittle materials. Computer Methods in Applied Mechanics and Engineering, 197(33-40), pp.2804-2822.
Willam, K.J 和 Carol, I. 長期致力于模擬巖石、混凝土和磚結構地開裂破壞。
本文的很多數學推導是通過軟件Wolfram Mathematica獲得的。
本文所用的算法是Caballero（2008）中的算法1：monolithic iteration strategy without substepping。
具體建模過程和建模結果見知乎文章：https://zhuanlan.zhihu.com/p/113538156
umat子程序和代碼對應的
詳細解答見附件。
展開 Abaqus仿真計算中的單元選擇
目前第一、二期直播已結束（聯系文末客服看回放），第三期直播<Abaqus仿真計算中的單元選擇>，已經開啟報名，歡迎參加~
”
對于有限元分析的網格模型，不僅需要較高的網格質量，還需要擁有合適的單元類型。ABAQUS為用戶提供了豐富的單元庫，幾乎可以模擬實際工程中任意幾何形狀的有限元模型，在對一個問題進行分析時，可以根據情況選擇使用。
如何才能選取出適合于分析的單元類型呢？
本次分享首先介紹ABAQUS中對于單元的分類，每種單元特定的使用范圍，各種單元類型的節點數目、單元形狀、插值函數階次以及單元構造的方式。
展開 巧用單元提高Abaqus計算效率：帶扭曲的軸對稱單元-橡膠阻尼器內摩擦生熱分析 ￥49.99
Abaqus有非常豐富的單元庫，其中就有軸對稱單元，比如CAX4（I/R/H/T），當一個回轉結構具有某種載荷對稱性時，可以用它將三維模型縮減為軸對稱模型來分析，能減少大量的內存和分析時間，而同樣的模型規模，3D實體單元要更耗費計算資源。
那么，回轉結構受到側向彎曲或軸向扭轉的載荷時，有沒有類似的單元可以用呢？
橡膠阻尼器的內摩擦生熱分析-節點溫度云圖
比如，假設上圖中的阻尼器不再是長方體，而是回轉體，且發生軸向扭曲變形，那么能不能用軸對稱單元來建模呢？
答案是可以的，在Abaqus的軸對稱單元系里還有一種可考慮Twist的單元，即帶字母G標識的那種類型，能夠在分析時充分考慮回轉體的整體扭轉變形。
首先，我們可以在part模塊使用Axisymmetric建立環形塊狀阻尼器的回轉截面；然后在mesh模塊劃分好四邊形網格；最后，定義單元類型為CGAX4T，即帶扭曲的4節點軸對稱位移-溫度耦合單元。
這里的橡膠阻尼器材料本構采用的是超彈性模型，應變能描述形式為Neo Hooke，再結合時域黏彈性Prony參數與非彈性變形能耗散比，來計算阻尼器周期性扭轉過程中的材料內摩擦生熱。
阻尼器上、下兩個端面的節點分別使用位于回轉軸上的兩個參考點來耦合，固定下端面參考點，并在上端面參考點施加軸向的周期性扭角位移。
阻尼器的回轉結構與網格-單元
雖然建模時只考慮了回轉截面，但是帶扭曲的軸對稱單元可以將回轉體發生扭轉時的整體結構響應考慮在內，這是因為這種單元多了一個扭轉自由度5，拿本例中的位移-溫度耦合單元CGAX4T來說，該單元的節點具有1、2、5和11四個自由度。
展開 abaqus單元平均應力的計算
abaqus計算結果可得到節點應力，那位大牛知道如何得到單元（或二維計算模型）的平均應力？
abaqus C3D8 單元 計算中采用了多少個積分點？
按照正常的理解，毫無.疑問，abaqus 全積分一定是采用了2x2x2=8個積分點。
從后處理結果來看，似乎也是如此，每個單元存在8個積分點。
然而，如果自己動手跑一遍程序，就會發現事實遠非如此，采用全積分計算得到的結果與abaqus 存在差異，原因何在？
事實賞，abaqus C3D8 采用的選擇積分方式（selective intergation schema），即對于偏應變，采用8個積分，對于球應變，采用中心點積分。這樣計算得到的結果才能與abaqus 完全對標，亦可從abaqus 幫助文檔得到答案。
展開 ABAQUS 二次開發（UEL+UMAT C3D8 單元動力隱式計算）
概述：采用UEL接口二次開發實現八節點單元，考慮BBAR修正，避免體積自鎖，對標ABAQUS自帶的C3D8單元，計算的剛度矩陣、質量矩陣和阻尼矩陣均與ABAQUS保持一致。并且采用UMAT子程序進行應力和應變數據的可視化，計算的應力應變數據同樣與ABAQUS保持一致，可視化效果同ABAQUS。以方塊的受動力簡諧荷載為例，采用上述程序，應用動力隱式計算分析步，最終計算的位移、應變等時程曲線均與ABAQUS保持一致。
（）模型信息
模型尺寸為10x10x10，彈性模量1e10，密度2000，泊松比0.25，荷載和邊界條件示意圖為：
一面的所有節點均固定。另一面所有的節點施加簡諧荷載，簡諧荷載曲線為（詳細的參數見附件）：
計算的網格圖為：
網格尺寸為1，共計10x10x=1000個單元。總的計算時長為1s，這只固定增量步長為0.01s，所以總增量步數為100。
（）計算結果
以上面網格圖中中間角點為例，提取加載向位移時程曲線如下圖：
0.89s時刻x向應力云圖為：
注：左一為ABAQUS計算結果，中間為umat實現線彈性各向同性本構計算結果，右一為UEL+UMAT計算結果，后面的云圖也是這樣排布。
0.89s時刻y向應力云圖為：
寫文字好麻煩，不寫了！放個動圖算了，有興趣的移步附件。
展開 第一篇梁單元的軸力圖 （理論計算、ABAQUS仿真、ANSYS仿真方法） ￥10
第一篇梁單元的軸力圖
（理論計算、ABAQUS仿真、ANSYS仿真方法）
篇幅內容僅針對自我學習總結展示，并希望給軟件初學者帶來一定啟發。
結構有限元仿真中有兩種一維單元：桁架與梁
桁架單元：僅承受軸力作用；如二力桿。由于只在軸向承受拉/壓載荷，所以只需要定義截面面積；應力和變形均與截面形狀無關。ABAQUS 6.14-4中對應單元為truss T2D2；ANSYS 18.0中對應單元為link180。
梁單元：可承受軸向拉/壓載荷，具有承受扭轉和彎曲的能力。由于可承受扭轉、彎曲等組合變形，梁單元需要定義截面形狀。ABAQUS與ANSYS對應均為beam單元。
孫訓芳先生的《材料力學》例題2-1：一等直桿及其受力情況如下圖，試作桿的軸力圖。
由于桁架單元僅能承受拉/壓載荷；而梁單元可承受拉、壓、彎曲、扭轉的組合變形，梁單元可承受的載荷類型更為復雜，故此篇通篇采用梁單元作為分析。
展開 ANSYS與ABAQUS關于梁單元后處理的計算與理論值比較（糾錯）- CAE夢想很偉大
ANSYS與ABAQUS關于梁單元后處理的計算與理論值比較（推薦）- CAE夢想很偉大
本文原創，若是轉載，請注明出處和筆名CAE-夢想很偉大。
感謝abaqus襄陽對于本文中錯誤Mises應力的問題的糾正。
本文目的
本文以工程項目中出現的評估問題為原型，以懸臂梁為例，對abaqus的mises應力在評估梁單元的如何獲得正確性進行說明。以理論計算為主，聯合ansys 和ansys workbench的計算結果，縱向評估正確的abaqus查看梁單元的正確用法beam-stress。
雖然本文可能小題大做，但是對于新手和一般不了解beam-mises的工程師，都希望引起足夠的重視。若是有任何異議，請大家留言，也歡迎大家留言討論。
具體內容如下
以10×10mm矩形截面，長度100mm的矩形管為例進行說明。
載荷：軸向載荷為10000N，彎矩為100N.m。通過理論計算
理論計算結果
軸向正應力為 ，
彎曲最大應力為
疊加組合應力
最大組合應力100+60=160
最小組合應力100-60=40
下面對比有限元計算結果與理論值比對，如表格所示
可以知道ANSYS、WB、ABAQUS顯示結果均與理論值一致。但是需要注意的是，ABAQUS需要修改截面顯示設置，需要考慮TOP和BOTTOM同時顯示數據，才能獲得正確的MISES結果。
ABAQUS的Mises不同截面激活設置顯示形式的比較如圖4所示。
展開 