金屬彈塑性各向異性的案例
ABAQUS各向異性彈塑性PUCK損傷VUMAT ￥5999
其核心功能包含三部分：首先基于正交各向異性彈性本構更新應力，通過材料屬性計算剛度矩陣并響應應變增量；其次實現彈塑性修正，采用J2流動理論判斷屈服狀態，通過牛頓迭代求解塑性變形并更新應力；最后建立漸進損傷模型，分別針對纖維方向（拉伸/壓縮失效）和基體方向（通過180°平面搜索臨界斷裂面）定義損傷初始判據，結合斷裂能與特征長度控制損傷演化過程。該模型通過21個狀態變量跟蹤材料歷史響應，包括塑性應變、損傷變量及主應變等。
該子程序專用于預測纖維增強復合材料（如碳纖維、玻璃纖維層合板）的漸進失效行為，適用于航空航天結構（機翼蒙皮、整流罩）、新能源汽車電池包防護結構、風電葉片等領域的強度分析與失效預測。其優勢在于精確模擬從初始彈性響應、塑性變形到最終斷裂的全過程，尤其擅長處理沖擊載荷、復雜應力狀態下的損傷演化問題。
展開 ABAQUS三維hill48彈塑性模型VUmat子程序（彈性為正交各向異性） ￥388
1.ABAQUS三維hill48彈塑性模型VUmat子程序
2.彈性階段為正交各項異性材料
3.hill48和正交各項異性材料參數參考ABAQUS靜力模塊自帶的模型參數
4.發貨方式為百度網盤鏈接，包含子程序及上面跑的兩個模型相關文件，包含Cae，inp文件，odb文件等
5.ABAQUS版本為2024，低版本可以利用導入inp文件的方式運行及修改
6.可以免費答疑三次，后續添加你自己的模型或者相關參數等輔導都可以優惠。
各向同性硬化彈塑性umat開發
1 說明
該本構完全從文檔《Writing User Subroutines with Abaqus》中摘抄而來，采用Fortran77格式編寫。
2 本構理論
3 與Abaqus自帶本構的對比
4 源代碼
iso_hardening_plasticity.f
各向同性硬化von Mises率無關彈塑性本構理論以及umat源代碼 ￥99
各向同性硬化von Mises率無關彈塑性本構理論以及umat源代碼
1 本構理論
1.1 率形式
對于各向同性線彈性材料，其本構方程為：
式中假設了應變張量可以分解為彈性應變和塑性應變兩部分：
因此塑性本構的關鍵在于計算塑性應變的演化。對于率無關彈塑性的本構理論，需要確定以下三個部分：
(1)：屈服條件
(2)：流動法則
(3)：硬化法則
在此采用的是 von Mises 屈服條件：
式中后繼屈服應力是等效塑性應變的函數：
流動法則為：
式中流動方向的表達式為：
硬化法則為：
1.2 Return-mapping算法
上述的本構方程均為率形式。在增量步中，給定增量應變：
首先假設該增量應變全為彈性應變，計算試驗狀態下的一些物理量：
試驗狀態下的應力
試驗狀態下的屈服函數值：
利用該試驗屈服函數值來判斷在該增量步下是否發生了塑性屈服。如果：
則說明試驗狀態即為真實狀態，即可進行更新：
反之則需要進行塑性更正，即需要計算塑性乘子的增量，利用以下非線性方程組進行計算：
可以將該非線性方程組簡化至一個非線性方程，過程如下，將該方程組中的第一式分解為球量和偏量兩部分：
因此可以計算應力為：
將上式中的第二式整理得到：
可以得到兩個張量的方向相同：
因此偏應力可以用試驗狀態的信息表示出來：
代入到最后一個一致性方程中可得：
即可利用牛頓迭代法對上述非線性方程進行求解，得到塑性乘子增量。
展開 
各向同性彈塑性本構的vumat源代碼：通過修改umat ￥99
更新狀態變量
alpha_n1 = alpha_n0
strain_pl_n1 = strain_pl_n0
Dp_1 = Dp_0
statev_n1(1) = alpha_n1
statev_n1(2:7) = strain_pl_n1
statev_n1(8) = Dp_1
endif
end subroutine plastic_iso_vumat
3 算例
3.1 單單元拉伸測試
對單個單元進行單軸拉伸，邊界條件如下：
von Mises應力對比結果如下（左圖為Abaqus材料庫計算，右圖為vumat子程序計算結果）：
等效塑性應變對比結果如下（左圖為Abaqus材料庫計算，右圖為vumat子程序計算結果）：
反力曲線對比如下：
塑性耗散曲線對比如下：
3.2 圓棒拉伸測試
對一圓棒骨料進行單軸拉伸，其邊界條件如下：
von Mises應力對比結果如下（左圖為Abaqus材料庫計算，右圖為vumat子程序計算結果）：
等效塑性應變對比結果如下（左圖為Abaqus材料庫計算，右圖為vumat子程序計算結果）：
反力曲線對比如下：
算例cae模型
abaqus_cae.zip
展開 