abaqus各向異性實體
關注創建者：王靖雯 創建時間：2023-02-27
abaqus各向異性實體的視頻教程
abaqus各向異性實體的實例教程
請問各位大佬，各向異性材料的塑性階段怎么設置參數呀
ABAQUS各向異性彈塑性PUCK損傷VUMAT ￥5999
該子程序為ABAQUS用戶自定義材料模型（VUMAT），用于模擬復合材料的非線性力學行為。其核心功能包含三部分：首先基于正交各向異性彈性本構更新應力，通過材料屬性計算剛度矩陣并響應應變增量；其次實現彈塑性修正，采用J2流動理論判斷屈服狀態，通過牛頓迭代求解塑性變形并更新應力；最后建立漸進損傷模型，分別針對纖維方向（拉伸/壓縮失效）和基體方向（通過180°平面搜索臨界斷裂面）定義損傷初始判據，結合斷裂能與特征長度控制損傷演化過程。該模型通過21個狀態變量跟蹤材料歷史響應，包括塑性應變、損傷變量及主應變等。
該子程序專用于預測纖維增強復合材料（如碳纖維、玻璃纖維層合板）的漸進失效行為，適用于航空航天結構（機翼蒙皮、整流罩）、新能源汽車電池包防護結構、風電葉片等領域的強度分析與失效預測。其優勢在于精確模擬從初始彈性響應、塑性變形到最終斷裂的全過程，尤其擅長處理沖擊載荷、復雜應力狀態下的損傷演化問題。
展開 由于砌體結構所采用的砌體材料具有明顯的正交各項異性，故先從正交各向異性彈性入手，根據彈性理論中的正交各向異性彈性理論，建立砌體的正交各向異性彈性本構模型，并將該彈性本構模型寫入Abaqus的材料子程序UMAT中，與Abaqus中自帶的正交各向異性彈性本構模型進行對比驗證，為后續砌體的正交各向異性彈塑性本構模型做好準備。
一、正交各向異性彈性基本理論
砌體的彈性各向異性主要是由其不同彈性特性的材料組分引起的（同樣研究復合材料時也可能會遇到相同問題）。當通過不同的方向測量砌體，會得到不同的砌體的彈性特性。屬于典型的正交各向異性材料，本文先從其平面正交各向異性彈性特性入手。
在正交各向異性材料的分析中，需要使用兩個坐標系統：材料坐標系統與整體坐標系統。以砌體為例，材料坐標是指由平行于砂漿接縫（1軸）和垂直于砂漿接縫（2軸）所形成的坐標系統。整體坐標系統指的是在結構體系下，平行于水平面（x軸）與垂直于水平面（y軸）所形成的坐標系統。材料坐標與整體坐標間的夾角為θ，二者的關系如下圖1所示：
圖1 正交各向異性材料的材料坐標(1-2)與整體坐標(x-y)示意圖
正交各項異性材料具有三個互相垂直坐標軸的材料彈性對稱性，將坐標軸x、y和z分別垂直于三個材料對稱，并要求繞這些軸轉動180°之后彈性性能不發生改變，由此XX中的常數具有一定的關系。
展開 一個基體里面只有一根熱導率各向異性的纖維，怎么在不和基體合并的情況下正常運算？或者說應該怎么將兩者綁定呢？
這是我第一次實踐各向異性超彈性本構子程序UANISOHYPER_INV，中間走了幾步彎路，好在最后問題都解決了。把這個過程記錄下來，為后人鑒。
1 超彈性本構
剛接觸超彈性本構的時候，很不適應。因為我之前研究的本構，都會給出非常明確的應力應變關系。比如最簡單的：應力=剛度矩陣×彈性應變。
超彈性本構一般不這么給，給的都是應變能和不變量之間的關系。比如這樣：
對于新的東西，我本能地用原有的知識體系去套。于是開始拼命的去檢索相關文獻，試圖找到超彈性本構應力和應變的關系到底是怎么表達的。
結果呢就是，撲街。大家都在各種秀張量，秀應變能，秀不變量。我一度認為這些人閑著沒事，凈搞形式。
但是也不是一無所獲，文獻中的蛛絲馬跡都指向了UANISOHYPER_INV子程序。
2 UANISOHYPER_INV子程序
UANISOHYPER_INV子程序是干啥的？它就是專門用來定義各向異性超彈性本構的。那么自然的你就會想，是不是還有專門定義各向同性超彈性本構的呢？當然，這個子程序叫UHYPER。
看懂了UANISOHYPER_INV子程序的設定，你就會恍然大悟，原來真的不需要定義應力應變關系，只要知道應變能和不變量的關系就行了啊。
UANISOHYPER_INV子程序的基本結構如下：
其中主要變量的介紹如下：
也就是說，在UANISOHYPER_INV子程序中定義出應變能、應變能對不變量的導數即可。
于是我按照幫助文檔的提示，一步步完成了子程序編寫。但是在測試的時候，問題接踵而來。
1 無法提交計算
做了一個簡單拉伸算例。但是提交計算時候，總是報錯：
報錯信息告訴我，可壓縮性材料不能用雜交單元。這個確實不能用，但是我好像沒有用啊。
展開 
abaqus各向異性實體的相關專題、標簽、搜索
abaqus各向異性實體的最新內容
該子程序為ABAQUS用戶自定義材料模型（VUMAT），用于模擬復合材料的非線性力學行為。其核心功能包含三部分：首先基于正交各向異性彈性本構更新應力，通過材料屬性計算剛度矩陣并響應應變增量；其次實現彈塑性修正，采用J2流動理論判斷屈服狀態，通過牛頓迭代求解塑性變形并更新應力；最后建立漸進損傷模型，分別針對纖維方向（拉伸/壓縮失效）和基體方向（通過180°平面搜索臨界斷裂面）定義損傷初始判據，結合斷裂能與特征長度控制損傷演化過程
1.ABAQUS三維hill48彈塑性模型VUmat子程序
2.彈性階段為正交各項異性材料
3.hill48和正交各項異性材料參數參考ABAQUS靜力模塊自帶的模型參數
4.發貨方式為百度網盤鏈接，包含子程序及上面跑的兩個模型相關文件，包含Cae，inp文件，odb文件等
5.ABAQUS版本為2024，低版本可以利用導入inp
寫在前文
材料的線彈性本構模型能夠很好的描述處于工作荷載水平下的材料性能情況，后續材料的塑性理論也需要在彈性本構模型的基礎上進行開展。由于砌體結構所采用的砌體材料具有明顯的正交各項異性，故先從正交各向異性彈性入手，根據彈性理論中的正交各向異性彈性理論，建立砌體的正交各向異性彈性本構模型，并將該彈性本構模型寫入Abaqus的材料子程序UMAT中，與Abaqus中自帶的正交各向異性彈性本構模型進行對比驗證
各向異性 復材層板 的 振型分析
建模介紹：
（1）復材層板，長30cm, 寬2cm, 厚0.5cm，共計10層的，每層厚度為0.05cm。
（2）鋪層角度為0度和90度交替，零度方向與長邊方向相同。
邊界條件：單邊固支
模擬結果：
能看到固有頻率、模態振型、應力、應變及位移云圖
復材層板鋪層角度的實現：
網格劃分細節
一個基體里面只有一根熱導率各向異性的纖維，怎么在不和基體合并的情況下正常運算？或者說應該怎么將兩者綁定呢？
請問各位大佬，各向異性材料的塑性階段怎么設置參數呀
這是我第一次實踐各向異性超彈性本構子程序UANISOHYPER_INV，中間走了幾步彎路，好在最后問題都解決了。把這個過程記錄下來，為后人鑒。
1 超彈性本構
剛接觸超彈性本構的時候，很不適應。因為我之前研究的本構，都會給出非常明確的應力應變關系。比如最簡單的：應力=剛度矩陣×彈性應變。
超彈性本構一般不這么給，給的都是應變能和不變量之間的關系。比如這樣：
對于新的東西
近在做有關木材的抗拔性能的分析，遇到了一系列的問題，比如各向異性材料參數填寫，材料方向分配，基于HILL屈服準則的材料參數等，經過一段時間的文獻調研，資料查詢，已經基本搞清楚，鑒于網上相關帖子比較少，所以特地寫一篇帖子，以供參考！
1. 3D模型建立。利用相關軟件（solidworks，creo或inventor），建立3D模型Part-1,Part-2，并導出STEP中性格式。
2. ABAQUS
