UMAT/VUMAT子程序開發(fā)的案例
UMAT/VUMAT從入門到進階子程序集合————復合材料系列 ￥66.59
UMAT/VUMAT從入門到進階子程序集合——Progressive damage analysis of Laminated Composites（PDALAC）層合板漸進損傷分析
==前情須知：以下內(nèi)容為外網(wǎng)搬運，屬于難得的高質(zhì)量子程序集合，由于完全開源，內(nèi)容不涉及侵權，不涉及抄襲，內(nèi)容收費并非版權收費==
適用人群：
1、關注復合材料漸進損傷分析的研究者
2、ABAQUS二次開發(fā)子程序編寫復合材料本構需求者
3、具備一定的自主學習能力和FORTRAN程序閱讀能力者
4、具備一定的英文閱讀能力的人（單詞不認識ctrl+c\ctrl+v查一下都不會的別碰）
5、不適合不具備獨立思考能力的人
6、歡迎各位大佬交流，如果有質(zhì)量比較好的代碼互相分享，可以直接加我QQ互幫互助，計算力學相關來者不拒
由于技術鄰發(fā)帖內(nèi)容有限，有需求的可以下載PDF深入了解本貼內(nèi)容。
UMATVUMAT從入門到進階子程序集合——Progressive damage analysis of Laminated Composites（PDALAC）層合板漸進損傷分析系列.pdf
以下為付費內(nèi)容部分展示
新手區(qū)：
進階區(qū)：
展開 UMAT/VUMAT從入門到進階子程序集合————復合材料系列 ￥66.59
UMAT/VUMAT從入門到進階子程序集合——Progressive damage analysis of Laminated Composites（PDALAC）層合板漸進損傷分析
==前情須知：以下內(nèi)容為外網(wǎng)搬運，屬于難得的高質(zhì)量子程序集合，由于完全開源，內(nèi)容不涉及侵權，不涉及抄襲，內(nèi)容收費并非版權收費==
適用人群：
1、關注復合材料漸進損傷分析的研究者
2、ABAQUS二次開發(fā)子程序編寫復合材料本構需求者
3、具備一定的自主學習能力和FORTRAN程序閱讀能力者
4、具備一定的英文閱讀能力的人（單詞不認識ctrl+c\ctrl+v查一下都不會的別碰）
5、不適合不具備獨立思考能力的人
6、歡迎各位大佬交流，如果有質(zhì)量比較好的代碼互相分享，可以直接加我QQ互幫互助，計算力學相關來者不拒
由于技術鄰發(fā)帖內(nèi)容有限，有需求的可以下載PDF深入了解本貼內(nèi)容。
UMATVUMAT從入門到進階子程序集合——Progressive damage analysis of Laminated Composites（PDALAC）層合板漸進損傷分析系列.pdf
以下為付費內(nèi)容部分展示
新手區(qū)：
進階區(qū)：
展開 vpsc-vumat子程序的開發(fā)與應用文章推薦
結(jié)合作者的理論（尤其是分段線性化和應力驅(qū)動的求解思路）我們可以把獨立的vpsc子程序編寫進abaqus里面，為了避免復雜的雅可比推導，以及適用各種復雜的變形工況，推薦使用abaqus的顯式求解器，即vumat程序
以下展示一個使用vpsc-鎂合金本構模型，模擬包含1個單元，單元包好100個晶粒在RD方向壓縮20%的模擬效果（原始模型參數(shù)取自vpsc官方案例,為了減少計算時間使用高應變率進行計算，應力高于vpsc模擬）：
等效塑性應變：
第一個晶粒的累計剪切滑移：
發(fā)生孿晶次數(shù);
變形后的形狀演化：
復合材料螺栓連接視頻模型及UMAT/VUMAT子程序 ￥200
1.復合材料螺栓連接視頻教程；
2.復合材料螺栓連接UMAT子程序；
3.復合材料螺栓連接VUMAT子程序；
4.復合材料螺栓連接接觸視頻案例
Abaqus內(nèi)置LaRC05失效準則子程序（附cohesive單元umat子程序開發(fā)教程）
Abaqus從2017版本開始，內(nèi)置了UVARM和DMGINI兩種子程序供用戶調(diào)用。UVARM子程序需要以“ABQ_LARC05_DMGCRT”作為材料名的前綴，可以用來評估是否滿足LaRC05強度指標；DMGINI子程序需要以“ABQ_LARC05_DMGINI”作為材料名前綴，可以結(jié)合XFEM來分析裂紋萌生和擴展。
在abaqus中可以通過修改關鍵字來進行模型設置
UVARM設置
DMGINI設置
輸出變量的含義
計算結(jié)果
此外，如果有小伙伴希望學習cohesive單元的umat子程序開發(fā)相關技術，可以關注下圖的教程：
點擊鏈接進入了解詳情：Abaqus Cohesive單元介紹及其本構的Umat子程序開發(fā)教程
最后，有ABAQUS子程序開發(fā)相關需求歡迎通過微信公眾號聯(lián)系我們。
微信公眾號：320科技工作室。
展開 Abaqus 子程序開發(fā)日志07-vumat之a(chǎn)nisotropic
（沒測試，自己寫了個umat轉(zhuǎn)vumat的接口，直接替換的。各向同性的是全寫的，非替換）
基于Abaqus的修正劍橋模型的vumat子程序開發(fā)
考慮到莫 爾–庫侖準則有此特征，故將莫爾–庫侖準則與劍橋模型相結(jié)合，得
聯(lián)合式(1)，(4)，(5)，有改進的屈服函數(shù)：
根據(jù)上述理論，可以通過以下流程圖完成VUMAT子程序編寫
通過單胞模型計算得到得結(jié)果如下圖所示。
[1]袁克闊,陳衛(wèi)忠,于洪丹,譚賢君,趙武勝,李香玲.考慮黏聚特性和拉壓不等效應的修正劍橋模型及數(shù)值實現(xiàn)[J].巖石力學與工程學報,2012,31(08):1574-1579.
最后，有需要歡迎通過微信公眾號聯(lián)系我們。
微信公眾號：320科技工作室。
Abaqus 修正劍橋模型的vumat子程序開發(fā)
考慮到莫 爾–庫侖準則有此特征，故將莫爾–庫侖準則與劍橋模型相結(jié)合，得
聯(lián)合式(1)，(4)，(5)，有改進的屈服函數(shù)：
根據(jù)上述理論，可以通過以下流程圖完成VUMAT子程序編寫
通過單胞模型計算得到得結(jié)果如下圖所示。
[1]袁克闊,陳衛(wèi)忠,于洪丹,譚賢君,趙武勝,李香玲.考慮黏聚特性和拉壓不等效應的修正劍橋模型及數(shù)值實現(xiàn)[J].巖石力學與工程學報,2012,31(08):1574-1579.
abaqus剪切修正GTN模型的VUMAT子程序開發(fā)
這里對相應的算法進行簡要說明：
NH-GTN模型
屈服函數(shù)：
其中等效孔洞體積分數(shù)定義為：
孔洞體積分數(shù)包含新孔隙形核，原有空隙生長以及剪切相關的等效體積分數(shù)增加：
形核，生長，剪切相關體積分數(shù)的演化遵循：
其中：
剪切效應的修正，考慮應力狀態(tài)的影響
參數(shù)的物理含義如下
通過將文獻中的數(shù)值算法編程實現(xiàn)在VUMAT子程序中，可以用來實現(xiàn)對延性金屬材料在不同應力狀態(tài)下的損傷演化進行合理的數(shù)值預測，應用于金屬成型領域（沖壓，軋制，擠壓等）
預測修正后的模型應該在簡單拉伸情況下于abaqus自帶的GTN模型保持相同的損傷和其他狀態(tài)變量的分布，并在剪切情況中損傷發(fā)展顯著高于abaqus自帶的模型（自帶的模型忽略了剪切效應）。（為了進行對比使用于自帶的本構相同的硬化方式，模擬中使用了相同的質(zhì)量縮放，但質(zhì)量縮放容易產(chǎn)生數(shù)值振蕩，模擬的拉伸曲線存在波動。）
初步模擬結(jié)果：
拉伸情況（abaqus-VUMAT）
應力分情況
孔洞體積分數(shù)
剪切模型（abaqus-VUMAT）
不同變形時刻的應力分布
T=0.1s
局部放大圖
T=0.5s
局部放大圖
T=0.6s
局部放大圖
可以看到模型在拉伸預測中與原始模型保持一致，而在剪切修正后損傷發(fā)展顯著快于原始模型，利用作者提出的方法可以應用于復雜應力狀態(tài)下金屬材料的損傷分析，相關參數(shù)部分參考文獻，其中Kw=3.T1=0.2，T2=0.7.模擬結(jié)果符合文獻所提出方法的基本趨勢。
最后，如果有相關需要歡迎通過公眾號“320科技工作室”聯(lián)系我們。
展開 【12月7日-9日 西安】Abaqus復合材料分析技術培訓班-三天
基礎班第二天主要講解復合材料失效理論、面內(nèi)漸進損傷分析、層間/界面漸進損傷分析、三維VUMAT子程序開發(fā)等內(nèi)容；通過提高班課程的學習，學員能夠掌握Abaqus初始損傷分析、漸進損傷分析、基于Cohesive內(nèi)聚力模型的層間/界面漸進損傷分析、基于VCCT虛擬裂紋閉合技術的層間/界面裂紋擴展分析以及VUMAT子程序開發(fā)與應用等。
高級班安排在第三天，主要講解多種復合材料子程序的開發(fā)、Puck失效理論子程序實現(xiàn)、材料非線性表征、率相關特定表征、子程序調(diào)試等方面的內(nèi)容。通過該課程學習，能夠掌握復合材料UVARM子程序、USDFLD子程序、VUSDFLD子程序、UMAT子程序、VUMAT子程序的開發(fā)，能夠基于VUMAT子程序進一步針對Puck失效理論、材料非線性、率相關性進行深入開發(fā)。該課程幾乎全部涉及程序開發(fā)，適合已掌握Abaqus基礎復合材料建模，在失效理論、程序開發(fā)、論文發(fā)表方面有更深入需求的研究生和科研人員。
學員可根據(jù)自身的條件和需求選擇適合自己的課程進行學習。
【基礎班課程大綱】：
【高級班課程大綱】：
【培訓講義】：
該培訓采用獨家內(nèi)部講義，該講義分為20個部分，共計256頁，每一部分形成獨立章節(jié)，從基本建模操作、理論講解、公式推導到程序編寫全面覆蓋。
展開 Abaqus 子程序開發(fā)日志01-umat之isotropic
分別用自帶的各向同性材料屬性與umat子程序進行計算。
自帶的材料屬性的設定以及模型的建立、邊界條件的設定不再說明，只對umat的使用進行說明。
此次直接使用編譯好的二進制文件，在材料屬性的輸入過程，在通用中找到用戶自定義材料，輸入設定的參數(shù)，第一個為楊氏模量，第二個為泊松比。如圖所示：
點擊ok就好，其余的和平時操作一致，只需在最后提交job時提供用戶子程序。
然后提交計算就好。結(jié)果如圖所示：
到此結(jié)束。cae文件如下：
umattest.cae
若選用縮減積分單元時要開啟沙漏控制。
各向同性umat：
isotropic-std.obj
共享庫：
standardU.dll
展開 
LS-DYNA用戶子程序（UMAT）開發(fā)案例
最近在International Journal of Impact Engineering（中科院一區(qū)）上看到一篇文章，題目是”Plasticity based material model for concrete subjected to dynamic loadings（https://doi.org/10.1016/j.ijimpeng.2020.103581）“，作者開發(fā)了一個適用于描述動態(tài)荷載作用下混凝土力學響應的塑性模型，并在LS-DYNA中通過用戶自定義子程序（UMAT）二次開發(fā)的方式實現(xiàn)了模型。論文中展示的模擬效果圖如下：
模擬鋼筋混凝土板的爆炸響應
我和這個作者要到了他開發(fā)的UMAT源代碼，在這里給大家共享一下，對于想學習LS-DYNA子程序開發(fā)的絕對是不可多得的資料（免費共享，不負責答疑哦）
（希望大家下載了關注我一下就行，有空可以看看我發(fā)的各種教學視頻和資料貼，985高校在讀博士，長期穩(wěn)定更新各種內(nèi)容）
UMAT.f
展開 技術鄰Abaqus復合材料高級培訓定制（長期有效）
基礎班第二天主要講解復合材料失效理論、面內(nèi)漸進損傷分析、層間/界面漸進損傷分析、三維VUMAT子程序開發(fā)等內(nèi)容；通過提高班課程的學習，學員能夠掌握Abaqus初始損傷分析、漸進損傷分析、基于Cohesive內(nèi)聚力模型的層間/界面漸進損傷分析、基于VCCT虛擬裂紋閉合技術的層間/界面裂紋擴展分析以及VUMAT子程序開發(fā)與應用等。
基礎班課程大部分都基于Abaqus/CAE進行操作，一小部分需要開發(fā)程序，適合尚不熟悉Abaqus復合材料分析的在校學生及工程單位CAE分析人員。
高級班主要講解多種復合材料子程序的開發(fā)、Puck失效理論子程序實現(xiàn)、材料非線性表征、率相關特定表征等方面的內(nèi)容。通過該課程學習，能夠掌握復合材料UVARM子程序、USDFLD子程序、VUSDFLD子程序、UMAT子程序、VUMAT子程序的開發(fā)，能夠基于VUMAT子程序進一步針對Puck失效理論、材料非線性、率相關性進行深入開發(fā)。該課程幾乎全部涉及程序開發(fā)，適合已掌握Abaqus基礎復合材料建模，在失效理論、程序開發(fā)、論文發(fā)表方面有更深入需求的研究生和科研人員。
【培訓要求】：
培訓期間，學員應自帶筆記本電腦，電腦中預先安裝好Abaqus以及與Abaqus版本對應的Visual Studio和Intel Visual Fortran，并完成Abaqus的子程序驗證。
展開 Abaqus 子程序開發(fā)日志05-umat之fully anisotropic
此umat的材料參數(shù)設置參考abaqus，參數(shù)自上而下分別為：
使用方法和之前一致。
若選用縮減積分單元時要開啟沙漏控制。
完全各向異性umat：
fullyanisotropic-std.obj
共享庫：
standardU.dll
Abaqus復合材料Larc05失效準則Umat子程序開發(fā)
之前在https://www.yqgqt.org.cn/content/post/1279990 這個帖子里介紹了Abaqus內(nèi)置的Larc05子程序調(diào)用方法，因為內(nèi)置的Larc05子程序是通過Udmgini子程序結(jié)合Xfem來實現(xiàn)的，因此還是存在一定的局限性的。本帖考慮了復合材料的剪切非線性行為，在材料本構中引入了Hahn-Tsai的剪切非線性模型
式中，β為剪切非線性因子。
通過編寫Umat子程序，對復合材料的損傷行為進行了模擬，計算結(jié)果如下。
90度鋪層損傷和載荷位移曲線
0度鋪層損傷和載荷位移曲線
45度鋪層損傷和載荷位移曲線
由上圖可以看出，引入剪切非線性后，45度鋪層試驗件的拉伸曲線表現(xiàn)出明顯的非線性行為
有關于子程序二次開發(fā)或者復材仿真的問題可以聯(lián)系QQ1653004885或者關注CAE320公眾號
展開 