子程序，cohesive單元，沖擊的案例
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科技工作室。
展開 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所示。
展開 cohesive界面單元UMAT子程序 ￥20
1 基于cohesive單元的裂縫擴展模擬
1.1 cohesive單元簡介
cohesive單元（粘性單元）定義： 通過預置裂紋邊（面）的方式來模擬二維（三維）裂紋，即在預判的裂紋區域加入一層厚度為0的cohesive單元。
cohesive單元簡界：如圖1所示，cohesive單元由頂面，中面和底面組成。cohesive模型通過損傷起始準則以及演化方法來判斷其損傷情況。損傷發生在cohesive單元上，當滿足損傷判據時cohesive單元開始進入損傷，根據預先定義的損傷演化模型，當單元損傷累積到完全失效時cohesive單元的中層;一分為二，使物體產生幾何上的不連續，從而形成開裂。
圖1 cohesive單元示意圖
1.2 斷裂力學原理
cohesive單元的損傷時由牽引-分離(traction-separation)定理描述的，在宏觀上牽引力與位移的關系如圖2所示。
圖2 牽引力-牽引位移示意圖
在損傷前，應力與應變滿足如下關系:ABAQUSlass="p1">
公式中，上標0單元分層損傷起始位移，f代表最終開裂位移，max代表加載過程中最大位移，并且滿足損傷不可逆性。
損傷起始準則用于判斷材料是否開始損傷。常見的損傷判據如下表。
當材料開始滿足損傷起始準則之后，表現出整體剛度退化的現象。并且損傷演化定義了剛度退化的規律和裂縫的擴展，同時損傷演化分為位移和能量兩種模式：位移即為圖2的橫軸，能量為圖中曲線圍成的面積。
展開 Cohesive單元的疲勞載荷子程序
有Cohesive單元的疲勞子程序應該怎么根據自己的模型修改子程序
Abaqus Cohesive單元的疲勞UMAT子程序
Abaqus內置的內聚力模型并未考慮循環載荷下的疲勞損傷，因此需要通過編寫UMAT子程序來實現循環載荷下的疲勞內聚力模型。
Roe提出了一種不可逆的內聚力模型來模擬界面的疲勞裂紋擴展行為。單調載荷下牽引力和位移之間符合指數關系
單調載荷下的損傷增量
同時，Roe提出了一種循環載荷下的損傷演化方程，將疲勞損傷與牽引力和累積位移聯系起來。
于是總損傷可以表示為
在Roe研究的基礎上，Emilio通過UEL編寫了考慮疲勞損傷的內聚力模型，但是其只考慮了法向的疲勞損傷，同時UEL在ABAQUS中的前后處理都不太方便，因此本文在前人研究的基礎上(增加了兩個假設，1，卸載過程經過原點；2，卸載時不產生疲勞損傷)編寫了考慮內聚力模型法向和切向疲勞損傷的UMAT子程序。
通過子程序計算得到的結果如圖3-圖7所示
圖3 單向加載曲線
圖4 法向循環載荷下的響應
圖5 切向循環載荷下的響應
圖6 循環彎曲載荷下的損傷演化和裂紋擴展
圖7 循環載荷下的DCB試件裂紋擴展
展開 Abaqus Cohesive單元的疲勞UMAT子程序
Abaqus內置的內聚力模型并未考慮循環載荷下的疲勞損傷，因此需要通過編寫UMAT子程序來實現循環載荷下的疲勞內聚力模型。
Roe提出了一種不可逆的內聚力模型來模擬界面的疲勞裂紋擴展行為。單調載荷下牽引力和位移之間符合指數關系
單調載荷下的損傷增量
同時，Roe提出了一種循環載荷下的損傷演化方程，將疲勞損傷與牽引力和累積位移聯系起來。
于是總損傷可以表示為
在Roe研究的基礎上，Emilio通過UEL編寫了考慮疲勞損傷的內聚力模型，但是其只考慮了法向的疲勞損傷，同時UEL在ABAQUS中的前后處理都不太方便，因此本文在前人研究的基礎上(增加了兩個假設，1，卸載過程經過原點；2，卸載時不產生疲勞損傷)編寫了考慮內聚力模型法向和切向疲勞損傷的UMAT子程序。
通過子程序計算得到的結果如圖3-圖7所示
圖3 單向加載曲線
圖4 法向循環載荷下的響應
圖5 切向循環載荷下的響應
圖6 循環彎曲載荷下的損傷演化和裂紋擴展
圖7 循環載荷下的DCB試件裂紋擴展
最后，如果您有相關需求，歡迎通過微信公眾號聯系我們。
展開 Abaqus纖維復合材料蜂窩板落錘沖擊仿真模型
內插0厚度cohesive單元以模擬分層
模擬過程采用puck子程序，有錄制整個建模操作視頻，可贈送復合材料層合板快速建模插件及蜂窩建模插件!
cae ￥20
image_process=/format,webp" data-initial-src="https://img.jishulink.com/202505/attachment/49623d80bdb74936898c3463aebb8345.png" data-extentions-extra-ocr-id="e6cb4a74c55e38de39a7e4f229d3e914"></figure>
</figure><div contenteditable="false" width="100%">
Abaqus纖維復合材料蜂窩板落錘沖擊仿真模型!
</div><div contenteditable="false" width="100%">
內插0厚度cohesive單元以模擬分層
</div><div contenteditable="false" width="100%">
模擬過程采用puck子程序，有錄制整個建模操作視頻，可贈送復合材料層合板快速建模插件及蜂窩建模插件!
</div><div contenteditable="false" width="100%">
cae,inp文件及ODB文件，操作視頻(注意：并未含puck子程序，僅作學習參考)
</div><p><br></p>
展開 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
文檔介紹了體單元、殼單元的復合材料層壓板靜力分析，多種定義鋪層方向、堆疊方向的方法，復合材料截錐殼鋪層方向的定義方法等。
展開 批量插入嵌入0厚度cohesive粘結單元模擬沖擊載荷下裂紋擴展特征
<p>本帖給出了通過批量插入0厚度的cohesive單元，模擬在沖擊條件下，小球開裂和基體開裂兩種情況下的效果圖，對整個模型都插入cohesive單元或者對局部插入cohesive單元或者其他特定條件下的插入都可以使用程序插入，在小球沖擊情況下，裂紋的擴展過程在下面給出來了，由于參數是隨便假設的，所以效果可能不是很理想，對于不同行業，可以通過調整實際的參數可以很好地模擬裂紋的擴展，</p><p>批量插入cohesive單元：</p><div contenteditable="false" width="100%"><img src="https://img.jishulink.com/upload/201806/1528620192337_六面體單元間嵌入cohesive單元.jpg" title="六面體單元間嵌入cohesive單元.jpg" alt="六面體單元間嵌入cohesive單元.jpg" style="max-width: 760px; width: 484px; height: 274px;" width="484" height="274" data-mobile-src="https://img.jishulink.com/upload/201806/1528620192337_六面體單元間嵌入cohesive單元.jpg?image_process=/format,webp/quality,q_40/resize,w_400" data-pc-src="https://img.jishulink.com/upload/201806/1528620192337_六面體單元間嵌入cohesive單元.jpg?
展開 批量嵌入cohesive單元源程序
批量嵌入cohesive單元程序（2D/3D）
截止目前位置，在論壇和群里也逛了好長時間了，發現現在很多很多的朋友都需求批量嵌入cohesive單元程序（不論是出于好奇，還是科研需要，或者其他…），個人在這上面也是下足了功夫，為了回饋壇友的幫助，同時也為了幫助現在急需程序的科研人員，現在超低價出售批量嵌入coh單元程序（不要又問：既然是回饋和幫助，為什么還要收費？這個問題不想回答，你覺得為什么呢，你可以不買的，不喜勿噴！）。
1. 2D實體單元間批量嵌入cohesive單元（嵌入區域：全部區域，部分區域，界面嵌入等，適用單元：四邊形單元或三角形單元）
效果圖展示：
2. 3D實體單元間批量嵌入cohesive單元（嵌入區域：全部區域，部分區域，集合嵌入等，適用單元：六面體單元或四面體單元）
效果圖展示：
程序價格在199-599之間，價格在不同時間段可能有輕微浮動，需要的同學請速度，可以郵箱我1057593923@qq.com，告知你的需求（也可以加qq，由于比較忙，很多時候看不到），看到后會第一時間回復大家！
Tip：
1. 售出的程序都是源程序；
2. 即便極個別不是源程序，但也沒有注冊要求，可以一程序多機器使用；
3. 程序購買后，只能用于科研學習使用，不得用于其他商業用途；
4. 個性化定制，程序功能算是比較全的了，請購買前描述你的問題，看看是否適合，避免不必要的浪費；
5. 程序的總體價格在199-599之間（需求功能不同，價格不同），這個價格相當人性了吧，畢竟都是源程序啊；
6. 程序編寫語言MATLAB或python，不懂沒關系，會用就行了嘛；
展開 cohesive單元破壞數提取程序 ￥200
cohesive單元破壞數提取程序
基于cohesive單元模型的裂紋長度提取源程序 ￥300
本程序可以提取每個時間步下裂紋的長度。
復合材料沖擊動力學合集及VUMAT子程序 ￥58
包含計算合集中所有碳纖維層合板低速沖擊cae模型，VUMAT子程序，子程序使用方法以及子程序注釋
三維hashin應變子程序模擬效果
基于Vumat子程序的復合材料壓縮沖擊破壞預測
利用這些理論，通過在子程序中構建本構方程，可實現在ABAQUS等有限元軟件中得到與實驗結果較為接近的仿真結果。同樣的，耐壓縮沖擊性能也是實際工程應用中經常需要考慮的一個性能，但由于復合材料壓縮的破壞形式較為復雜，且對于壁板結構，壓縮過程中又涉及到了穩定性問題，所以對其進行有效的預測具有較大的難度。
本期，我們就講講如何在ABAQUS中基于Vumat子程序來實現對復合材料在壓縮沖擊載荷作用下，有效地預測其破壞模式。
1. 方法概述—“兩步走”
對于本文案例類似的結構，其在壓縮沖擊載荷作用下的破壞模式預測大致可分為兩大步，即穩定性分析和破壞模式預測。
1. 首先，我們需要建立一個屈曲分析步，引入擾動，在Model-edit keyword中添加相關關鍵字，并輸出節點位移文件。
2. 然后，復制模型，修改分析步類型及邊界條件，并在Vumat子程序中構建材料本構，引入失效準則，利用ABAQUS顯式運算的方法，不停迭代產生損傷后的剛度矩陣，直至結構發生破壞。同樣的，在此步驟中也需要對關鍵字進行編輯。
由于壓縮沖擊的破壞模式多種多樣，若需模擬某些特定的損傷模式，可在子程序中進行單元刪除的定義。
2. 案例
本文以碳纖維增強復合材料薄壁方管為例，利用Vumat子程序預測其在沿軸長方向的壓縮沖擊載荷作用下的破壞模式。
2.1 問題描述
試驗件由[0°/90°]的碳纖維布采用RTM工藝固化成型。實驗時，試驗件放置在試驗機底座上，上壓頭以1m/s的速度對試驗件進行壓縮沖擊，其有限元模型如圖1所示。
展開 Abaqus 通過VUMAT子程序模擬復合材料的低速沖擊行為
單元的應力應變關系
[1]. 徐穎, 溫衛東與崔海坡, 復合材料層合板低速沖擊逐漸累積損傷預測方法. 材料科學與工程學報, 2006(01): 第77-81頁.
[2]. 王躍全, 童明波與朱書華, 三維復合材料層合板漸進損傷非線性分析模型. 復合材料學報, 2009(05): 第159-166頁.
最后，有相關需求，歡迎通過微信公眾號聯系我們。
微信公眾號：320科技工作室。
