cohesive vumat 子程序的案例
cohesive界面單元VUMAT子程序 ￥20
1 基于cohesive單元的裂縫擴展模擬
1.1 cohesive單元簡介
cohesive單元（粘性單元）定義： 通過預置裂紋邊（面）的方式來模擬二維（三維）裂紋，即在預判的裂紋區域加入一層厚度為0的cohesive單元。
cohesive單元簡界：如圖1所示，cohesive單元由頂面，中面和底面組成。cohesive模型通過損傷起始準則以及演化方法來判斷其損傷情況。損傷發生在cohesive單元上，當滿足損傷判據時cohesive單元開始進入損傷，根據預先定義的損傷演化模型，當單元損傷累積到完全失效時cohesive單元的中層;一分為二，使物體產生幾何上的不連續，從而形成開裂。
圖1 cohesive單元示意圖
1.2 斷裂力學原理
cohesive單元的損傷時由牽引-分離(traction-separation)定理描述的，在宏觀上牽引力與位移的關系如圖2所示。
展開 Abaqus復合材料仿真指南：UMA/VUMAT子程序、Cohesive單元入門（資源分享，錯過后悔
六.案例/程序分享
1.Abaqus復合材料層壓結構快速建模插件（附多個版本執行程序）
同樣是復合材料力學老師做的一個復合材料建模插件，超好用！
2.復合材料沖擊動力學合集及VUMAT子程序
包含計算合集中所有低速沖擊模型，VUMAT子程序，子程序使用方法以及子程序注釋，老師的主頁還有很多視頻演示，非常良心。
3.UMAT/VUMAT從入門到進階子程序集合————復合材料系列
帖主分享的自己整理的子程序合集，包含基于不同失效準則和不同損傷模型的UMAT/VUMAT源代碼。
4.基于Python腳本提取復合材料應力應變均勻化有效彈性模量
基于均勻化理論的確定復合材料結構均勻應力應變場，進一步求解單向加載情況下的均勻模量的腳本文件。
5.ABAQUS復合材料建模插件CM TOOL2.0（GUI界面）
插件名為CM TOOL，本插件件基于Python結合ABAQUS二次開發，實現復合材料幾何模型的建快速建立。
6.ABAQUS復合材料UMA/VUMAT子程序豪華大禮包
包括七種失效準則、五種損傷演化模式，30個UMAT/VUMAT子程序
七.文檔分享
1.UMAT復合材料損傷子程序詳解.pdf
各種常用Abaqus UMAT子程序，以及相應變量參數的介紹和解釋。
2.Abaqus復合材料分析要點詳解-kxh.pdf
文檔介紹了體單元、殼單元的復合材料層壓板靜力分析，多種定義鋪層方向、堆疊方向的方法，復合材料截錐殼鋪層方向的定義方法等。
展開 Abaqus內置LaRC05失效準則子程序（附cohesive單元umat子程序開發教程）
Abaqus從2017版本開始，內置了UVARM和DMGINI兩種子程序供用戶調用。UVARM子程序需要以“ABQ_LARC05_DMGCRT”作為材料名的前綴，可以用來評估是否滿足LaRC05強度指標；DMGINI子程序需要以“ABQ_LARC05_DMGINI”作為材料名前綴，可以結合XFEM來分析裂紋萌生和擴展。
在abaqus中可以通過修改關鍵字來進行模型設置
UVARM設置
DMGINI設置
輸出變量的含義
計算結果
此外，如果有小伙伴希望學習cohesive單元的umat子程序開發相關技術，可以關注下圖的教程：
點擊鏈接進入了解詳情：Abaqus Cohesive單元介紹及其本構的Umat子程序開發教程
最后，有ABAQUS子程序開發相關需求歡迎通過微信公眾號聯系我們。
微信公眾號：320科技工作室。
展開 Hill VUMat 子程序（ ￥20
subroutine vumat(
C Read only -
1 nblock, ndir, nshr, nstatev, nfieldv, nprops, lanneal,
2 stepTime, totalTime, dt, cmname, coordMp, charLength,
3 props, density, strainInc, relSpinInc,
4 tempOld, stretchOld, defgradOld, fieldOld,
3 stressOld, stateOld, enerInternOld, enerInelasOld,
6 tempNew, stretchNew, defgradNew, fieldNew,
C Write only -
5 stressNew, stateNew, enerInternNew, enerInelasNew )
C
include 'vaba_param.inc'
C
C All arrays dimensioned by (*) are not used in this algorithm
dimension props(nprops), density(nblock),
1 coordMp(nblock,*),
2 charLength(*), strainInc(nblock,ndir+nshr),
3 relSpinInc(*), tempOld(*),
4 stretchOld(*), defgradOld(*),
5 fieldOld(nblock,nfieldv
展開 
VUMAT子程序使用方法
VUMAT子程序使用方法
在ABAQUS屬性模塊（property）中
1. props ( * )
* 數值代表在在User Material 中填寫的對應數值，例如E1 = props(1) 對應的是對一個數值130000，即后續在子程序中出現E1代表數值130000。
2. 子程序中stateNew(k,*) 與Depvar相對應
②代表變量stateNew(k,*) 的總數，填寫方法為：檢查.for文件中stateNew(k,*)中*的最大值，填寫要求是大于這個最大值均可。例如.for文件中stateNew(k,*)中*的最大值為9，即填寫比9大的數值都是可以的，例如9,10,11 …..均可。③對應單元刪除狀態變量，即找出單元刪除判定條件，滿足對應的損傷變量大于等于1，一般取纖維拉伸損傷變量（纖維方向應力S11大于0下的損傷）
if (statenew(k,*) .ge. one ) then
stateNew(k,**) = zero
end if
此時**就是單元刪除變量，即填寫③中。
②和③的區別，②的狀態變量初始值（除③對應的數值）均為0；③的狀態變量默認值是1，及初始值stateNew(k,11)=1
3.在Job模塊中
選擇General ＞ User subroutine file 文件夾中對應的子程序 .for 文件，確保目錄文件名均為英文名稱，然后提交任務即可。
展開 求購VUMAT子程序
求得基于應變的Hashin,Chang-Chang三維損傷失效判據的VUMAT子程序，有償，私聊，謝謝
VUMAT子程序簡介
VUMAT子程序簡介
有償求vumat子程序
有償求vumat子程序，復合材料漸進損傷子程序，帶能量演化(線性剛度退化)
ABAQUS用戶子程序VUMAT
Abaqus用戶子程序VUMAT.pdf
vpsc-vumat子程序的開發與應用文章推薦
結合作者的理論（尤其是分段線性化和應力驅動的求解思路）我們可以把獨立的vpsc子程序編寫進abaqus里面，為了避免復雜的雅可比推導，以及適用各種復雜的變形工況，推薦使用abaqus的顯式求解器，即vumat程序
以下展示一個使用vpsc-鎂合金本構模型，模擬包含1個單元，單元包好100個晶粒在RD方向壓縮20%的模擬效果（原始模型參數取自vpsc官方案例,為了減少計算時間使用高應變率進行計算，應力高于vpsc模擬）：
等效塑性應變：
第一個晶粒的累計剪切滑移：
發生孿晶次數;
變形后的形狀演化：
Abaqus修正GTN模型的VUMAT子程序
本文在GTN模型中引入剪應力的影響，編寫了相關的VUMAT子程序。
GTN模型的屈服函數可以用下式表示
其中q1，q2是模型參數，取q1=1.5，q2=1,σ0為等效應力，p為靜水應力，q為Mises等效應力；f為空洞的體積分數。
p和q可以通過徑向返回算法得到
應變控制的孔洞形核系數
GTN模型可以通過以下4個方程進行描述
Nahshon and Hutchinson考慮了剪應力對模型的影響
于是孔隙體積分數的演化可以通過下式描述
仿真計算得到的結果如下圖所示
有問題私信或者關注cae320公眾號
ABAQUS的vumat子程序
粘超彈性子程序有大神代做嘛？有需要的請聯系QQ615329992，
開源Johnson-Cook損傷vumat子程序
其次是計算效率與數值穩定性極佳，它的數學形式簡潔高效，非常適合顯式動力學子程序（如 VUMAT）進行大規模并行計算，不易發生數值發散。最后是完美閉環了“力-熱-損傷”的耦合，它不僅能算應力，還能同步算出溫度升高以及材料的受損程度，在模擬金屬穿透、飛濺、切屑形成等斷裂失效行為時，具有無與倫比的仿真精度和視覺逼真度。
這里分享一個經典的vumat子程序，方便大家學習Johnson-Cook的相關理論模型：
原始鏈接：https://github.com/mauroarcidiacono/Abaqus-VUMAT-Johnson-Cook/tree/main
代碼由Arcidiacono, Mauro F. and Rahimi, Salaheddin等人開發
! ########################################################################! User subroutine to model the Johnson-Cook plasticity, damage and the ! Taylor-Quinney conversion of mechanical work into heat during plastic! deformation.!! Abaqus version: Abaqus 2022! Intel Fortran Compiler 2021.11! Visual Studio 2019!! Author: Mauro Francisco Arcidiacono! ########################################################################!
展開 Abaqus Vumat子程序計算復合材料損傷 ￥10
用于復合材料三維實體單元，3D Hashin損傷準則
Abaqus 子程序開發日志07-vumat之anisotropic
（沒測試，自己寫了個umat轉vumat的接口，直接替換的。各向同性的是全寫的，非替換）
