UMAT子程序的案例
線彈性UMAT子程序驗證-對初級子程序用戶很有幫助 ￥3
umat子程序可以定義材料的本構(gòu)關(guān)系，abaqus 中自帶的材料模型通常為成熟的材料模型，當(dāng)新的材料行為被發(fā)現(xiàn)時，通常沒法應(yīng)用自帶的材料的模型，這時就需要用到umat子程序。現(xiàn)在以大家熟知的線彈性模型為例，以熟悉uamt子程序的語法結(jié)構(gòu)，并驗證其正確性。
問題描述：
一個0.2×0.2×1m的由Q235鋼組成的細桿，其一端固定，另一端面上施加100Mpa的拉力，我們都知道Q235鋼的屈服極限為235Mpa，因此這是個線彈性問題，可以用abaqus自帶的材料模型解決，但是為了熟悉umat子程序語法結(jié)構(gòu)并驗證其正確性，這里用umat子程序自定義材料本構(gòu)關(guān)系，對上述簡單拉伸問題進行模擬。
具體步驟如下：
1、建立部件
2、輸入材料參數(shù)
3、創(chuàng)建實例
4、定義分析步
5、定義邊界條件
6、劃分網(wǎng)格
7、提交作業(yè)
8、顯示結(jié)果
我們在細桿上任選一點，z方向的應(yīng)力為100Mpa，計算結(jié)果正確，從而驗證了umat子程序的正確性。
接下來，通過視頻的形式給大家詳細介紹本算例的umat子程序
展開 線性強化彈塑性umat子程序系列-子程序詳解 ￥5
對于這個問題，通過簡單計算可以發(fā)現(xiàn)試樣已經(jīng)發(fā)生塑性變形，通過自編的Umat子程序計算最后試樣應(yīng)力為509.1MPa。我們知道這個問題是有理論解的，下面我們來求理論解。
先求硬化系數(shù)H
總的應(yīng)變?yōu)?聯(lián)立后，試樣的軸向應(yīng)力為
基于Umat子程序的計算結(jié)果與理論值完全一致。
接下來請大家觀看Umat子程序逐句編寫視頻，也歡迎大家下載本次的Abaqus模型文件和*.for文件。
或者觀看視頻教程。
http://www.yqgqt.org.cn/college/video/c14014
基于Abaqus的umat子程序實現(xiàn)三維編織復(fù)合材料疲勞性能預(yù)測
（1）在Umat子程序中選用三維Hashin準(zhǔn)則進行纖維束的失效判斷，選用最大應(yīng)力準(zhǔn)則進行基體的失效判斷。
（2）確定單胞疲勞失效載荷、極限應(yīng)變，和纖維束平均應(yīng)力。
（3）計算時，計算模型需對界面和邊界條件采用一定的假設(shè)與簡化。提取加載點位移載荷數(shù)據(jù)。
圖1 單胞靜力分析載荷位移曲線
2.2 疲勞分析
（1）建立纖維束、界面、基體的疲勞折減模型，定義損傷因子，這些均在Umat子程序中加以實現(xiàn)。同樣的，在建立疲勞折減模型時也需引入一些的假設(shè)以簡化模型。在本此研究模型中，首先利用以下文獻中提供的單向復(fù)合材料疲勞壽命預(yù)測公式確定不同應(yīng)力水平下纖維束的疲勞壽命：
（2）依據(jù)計算所得的疲勞壽命，在Umat子程序中實現(xiàn)不同應(yīng)力水平下纖維束、界面和基體的剛度折減、強度折減。
（3）最后，利用子程序實現(xiàn)載荷步的循環(huán)加載，同時在每次循環(huán)中對單胞的疲勞狀態(tài)及失效狀態(tài)進行判斷，直至單胞達到最大疲勞壽命狀態(tài)。
2.3 結(jié)果展示
（1）單胞應(yīng)力水平與疲勞壽命的關(guān)系
本文分別依據(jù)以上方法，進行了不同應(yīng)力水平的單胞疲勞壽命預(yù)測，單胞應(yīng)力水平與疲勞壽命的關(guān)系如圖2所示。
圖2 單胞疲勞壽命預(yù)測結(jié)果
（2）損失演化
其中，纖維束的損傷演化如圖3所示。
展開 UMAT子程序（晶體塑性力學(xué)）經(jīng)典案例-單胞模型（inp+UMAT文件+子程序對應(yīng)的本構(gòu)模型文件） ￥10
這份資料是從事固體力學(xué)研究幾年經(jīng)典推薦教程，助你在有限元仿真理論部分有更深入的理解和認識，同時對有限元材料本構(gòu)模型的UMAT子程序的編寫、材料參數(shù)的設(shè)置、ABAQUS的前處理有更加深入的感悟。
附件中文件：inp單胞模型，UMAT晶體塑性經(jīng)典子程序，子程序對應(yīng)的經(jīng)典說明，固體力學(xué)國際經(jīng)典教程（Computational Methods for Plasticity），ABAQUS工程項目前處理經(jīng)典教程。
看好再購買，售出不退，謹慎入手
一種可用于形狀記憶合金（SMA）的UMAT子程序 ￥29.99
多尺度模型
子程序考慮微觀晶粒與宏觀單元的相互作用：
1) 宏觀單元由多個晶粒組成，每個晶粒的相變行為獨立計算；
2) 單元整體響應(yīng)為各晶粒響應(yīng)的加權(quán)平均，可模擬晶粒取向?qū)暧^行為的影響（如案例中 Element 1 及其不同晶粒的應(yīng)變差異）。
3、 案例介紹和結(jié)果對比
1. 案例介紹
為驗證 UMAT 子程序的有效性，構(gòu)建 NiTi 合金單向拉伸模型，參數(shù)如下：
幾何尺寸：矩形試件，長寬高均為1mm；
加載條件：位移控制加載，位移范圍0-0.05mm
材料參數(shù)：楊氏模量 E=40GPa，泊松比 ν=0.33，初始屈服應(yīng)力 σ0=353MPa，相變臨界應(yīng)力 σ_f=381MPa（正向）、σ_s=141MPa（反向）。
2.結(jié)果對比
實驗數(shù)據(jù)與模型預(yù)測結(jié)果如圖 1 所示（曲線趨勢與文獻 [鄒京辰等，2025] 一致）：
試件的力 - 位移曲線和應(yīng)力-應(yīng)變曲線均呈現(xiàn)典型超彈性特征：加載階段因奧氏體→馬氏體相變出現(xiàn)應(yīng)力平臺，卸載階段因反向相變應(yīng)力驟降；
曲線趨勢與文獻結(jié)果基本吻合，驗證了子程序對相變力學(xué)行為的精準(zhǔn)捕捉。
這個 UMAT 展示了如何在標(biāo)準(zhǔn)塑性框架內(nèi)嵌入相變效應(yīng)，為模擬如形狀記憶合金 (SMA)、相變誘發(fā)塑性 (TRIP) 鋼等智能材料或先進金屬提供了基礎(chǔ)。理解和應(yīng)用此代碼需要對彈塑性力學(xué)理論、ABAQUS UMAT 接口和特定材料的相變機制有深入的了解。
4、 代碼解釋以及案例文件（inp，umat子程序)
展開 非線性強化彈塑性umat子程序教程
其硬化曲線的方程為
自定義材料模型為
加載方式為
計算結(jié)果為
對于這個問題，通過簡單計算可以發(fā)現(xiàn)試樣已經(jīng)發(fā)生塑性變形，通過自編的Umat子程序計算最后試樣應(yīng)力為300MPa。我們知道這個問題是有理論解的，下面我們來求理論解。
假設(shè)z向應(yīng)力為，總的應(yīng)變?yōu)?聯(lián)立后，得方程
解方程，并取大于200的解為試樣的軸向應(yīng)力
基于Umat子程序的計算結(jié)果與理論值完全一致。
感興趣的朋友，請觀看非線性強化彈塑性umat子程序視頻教程，也歡迎大家下載本次的Abaqus模型文件和*.for文件。
視頻教程鏈接：
http://www.yqgqt.org.cn/college/video/c14070
展開 Johnson-cook 本構(gòu)模型 的umat子程序 ￥299
其中塑性模型相關(guān)的材料參數(shù)擬合可以參考網(wǎng)址：
abaqus中Johnson-Cook本構(gòu)模型理解 - 知乎 (zhihu.com)
損傷相關(guān)參數(shù)擬合可以參考網(wǎng)址：
abaqus中Johnson-Cook損傷失效模型 - 知乎 (zhihu.com)
umat子程序驗證
相同材料參數(shù)下，使用編寫的umat子程序與顯式求解器內(nèi)置的模型進行對比，模型材料參數(shù)取自文獻
模型尺寸為10*10*50（mm），沿著Z方向進行20%的拉伸。
材料參數(shù)分別為
楊氏模量：209Gpa
泊松比：0.29
塑性產(chǎn)熱分數(shù)：0.9
材料參數(shù)A：380
材料參數(shù)B：1200
材料參數(shù)C：0.012
材料參數(shù)m：2.5
材料參數(shù)n：0.55
材料參數(shù)epsilon_0：1.0
轉(zhuǎn)變溫度Tr：293K
熔融溫度Tm：2150K
材料參數(shù)d1-d5：0.1，0.7，1.5，0.005，0.84
模擬的結(jié)果
變形后顯式與隱式結(jié)果的應(yīng)力分布：
變形后顯式與隱式結(jié)果的等效塑性應(yīng)變分布：
可以看到顯式與編寫的隱式umat子程序具有良好的一致性，其中隱式計算需要大概600個增量步，而顯式則需要700萬的增量步，計算效率更高，但在高速沖擊，切屑問題更建議使用顯式程序進行計算
展開 Abaqus內(nèi)置LaRC05失效準(zhǔn)則子程序（附cohesive單元umat子程序開發(fā)教程）
Abaqus從2017版本開始，內(nèi)置了UVARM和DMGINI兩種子程序供用戶調(diào)用。UVARM子程序需要以“ABQ_LARC05_DMGCRT”作為材料名的前綴，可以用來評估是否滿足LaRC05強度指標(biāo)；DMGINI子程序需要以“ABQ_LARC05_DMGINI”作為材料名前綴，可以結(jié)合XFEM來分析裂紋萌生和擴展。
在abaqus中可以通過修改關(guān)鍵字來進行模型設(shè)置
UVARM設(shè)置
DMGINI設(shè)置
輸出變量的含義
計算結(jié)果
此外，如果有小伙伴希望學(xué)習(xí)cohesive單元的umat子程序開發(fā)相關(guān)技術(shù)，可以關(guān)注下圖的教程：
點擊鏈接進入了解詳情：Abaqus Cohesive單元介紹及其本構(gòu)的Umat子程序開發(fā)教程
最后，有ABAQUS子程序開發(fā)相關(guān)需求歡迎通過微信公眾號聯(lián)系我們。
微信公眾號：320科技工作室。
展開 下加載面修正劍橋模型及其在Abaqus中umat子程序的實現(xiàn)
這 個基本特征給程序的編程工作帶來了很大的方便，因為當(dāng)前的應(yīng)力點都處于下 加載屈服面上，因此不需要判斷此時的應(yīng)力狀態(tài)是否到達屈服面。
二。Abaqus的umat子程序實現(xiàn)
子程序編寫流程如下所示：
三。模型驗證
下面給出一個簡單的應(yīng)用 UMAT 子程序的算例，將模型劃分為多個單元， 分別采用位移的加載方式和力的加載方式，并分別與試驗數(shù)據(jù)和 ABAQUS 自帶 的修正劍橋模型計算結(jié)果進行對比，以此說明 UMAT 子程序可以應(yīng)用于多個單 元的模型，對于位移和力的加載方式均適用，并證明 UMAT 子程序適用于超固 結(jié)土體的力學(xué)特性分析和地基的固結(jié)沉降分析。
從而證明了編寫的子程序的有效性和實用性。
最后，大家有關(guān)于Abaqus二次開發(fā)的任何需求可以添加管理員微信號：CAE320，同時也歡迎大家關(guān)注“320科技工作室”的微信公眾號，掃一掃二維碼即可關(guān)注~~
展開 基于JC模型的umat子程序實例-Hopkinson壓桿的隱式求解器數(shù)值仿真實現(xiàn) ￥8
眾所周知，Abaqus中自帶的JC模型只能應(yīng)用于顯示求解器（Abaqus/Explicit），但是隱式求解器(Abaqus/Standard)通常具有更高的精度，而且隨著研究的進步JC模型也出現(xiàn)了新的修正形式，因此我們有必要基于Umat實現(xiàn)傳統(tǒng)的JC模型，以便有需要的朋友在此基礎(chǔ)上進一步創(chuàng)新，避免重復(fù)“造車”！本文基于JC模型的Umat子程序實現(xiàn)hopkinson壓桿的數(shù)值模擬，本文的一些參數(shù)設(shè)置和程序借鑒了盧劍鋒老師的碩士論文《沖擊載荷作用下材料和結(jié)構(gòu)力學(xué)行為有限元模擬》，且沒有考慮溫度軟化。
裝配圖如下
JC模型的方程形式如下
試樣的材料參數(shù)如下
入力桿和出力桿使用線彈性材料，彈性模量和泊松比分別為200GPa和0.3。
在入力桿的一端施加應(yīng)力脈沖
入力桿某點和出力桿某點上的應(yīng)力波時程曲線如下
建議大家按照前面提到的論文自己獨立把程序走一遍，也歡迎大家下載本次的*.cae文件和*.for文件，后期打算基于JC模型的umat子程序做一個侵徹的案例，感謝持續(xù)關(guān)注！
展開 Abaqus Cohesive單元的疲勞UMAT子程序
Abaqus內(nèi)置的內(nèi)聚力模型并未考慮循環(huán)載荷下的疲勞損傷，因此需要通過編寫UMAT子程序來實現(xiàn)循環(huán)載荷下的疲勞內(nèi)聚力模型。
Roe提出了一種不可逆的內(nèi)聚力模型來模擬界面的疲勞裂紋擴展行為。單調(diào)載荷下牽引力和位移之間符合指數(shù)關(guān)系
單調(diào)載荷下的損傷增量
同時，Roe提出了一種循環(huán)載荷下的損傷演化方程，將疲勞損傷與牽引力和累積位移聯(lián)系起來。
于是總損傷可以表示為
在Roe研究的基礎(chǔ)上，Emilio通過UEL編寫了考慮疲勞損傷的內(nèi)聚力模型，但是其只考慮了法向的疲勞損傷，同時UEL在ABAQUS中的前后處理都不太方便，因此本文在前人研究的基礎(chǔ)上(增加了兩個假設(shè)，1，卸載過程經(jīng)過原點；2，卸載時不產(chǎn)生疲勞損傷)編寫了考慮內(nèi)聚力模型法向和切向疲勞損傷的UMAT子程序。
通過子程序計算得到的結(jié)果如圖3-圖7所示
圖3 單向加載曲線
圖4 法向循環(huán)載荷下的響應(yīng)
圖5 切向循環(huán)載荷下的響應(yīng)
圖6 循環(huán)彎曲載荷下的損傷演化和裂紋擴展
圖7 循環(huán)載荷下的DCB試件裂紋擴展
展開 
Abaqus Cohesive單元的疲勞UMAT子程序
Abaqus內(nèi)置的內(nèi)聚力模型并未考慮循環(huán)載荷下的疲勞損傷，因此需要通過編寫UMAT子程序來實現(xiàn)循環(huán)載荷下的疲勞內(nèi)聚力模型。
Roe提出了一種不可逆的內(nèi)聚力模型來模擬界面的疲勞裂紋擴展行為。單調(diào)載荷下牽引力和位移之間符合指數(shù)關(guān)系
單調(diào)載荷下的損傷增量
同時，Roe提出了一種循環(huán)載荷下的損傷演化方程，將疲勞損傷與牽引力和累積位移聯(lián)系起來。
于是總損傷可以表示為
在Roe研究的基礎(chǔ)上，Emilio通過UEL編寫了考慮疲勞損傷的內(nèi)聚力模型，但是其只考慮了法向的疲勞損傷，同時UEL在ABAQUS中的前后處理都不太方便，因此本文在前人研究的基礎(chǔ)上(增加了兩個假設(shè)，1，卸載過程經(jīng)過原點；2，卸載時不產(chǎn)生疲勞損傷)編寫了考慮內(nèi)聚力模型法向和切向疲勞損傷的UMAT子程序。
通過子程序計算得到的結(jié)果如圖3-圖7所示
圖3 單向加載曲線
圖4 法向循環(huán)載荷下的響應(yīng)
圖5 切向循環(huán)載荷下的響應(yīng)
圖6 循環(huán)彎曲載荷下的損傷演化和裂紋擴展
圖7 循環(huán)載荷下的DCB試件裂紋擴展
最后，如果您有相關(guān)需求，歡迎通過微信公眾號聯(lián)系我們。
展開 Abaqus使用umat子程序的沙漏問題
在使用umat子程序時，采用縮減積分單元后，沙漏控制剛度是通過材料屬性中的彈性性質(zhì)定義的，這些剛度基于材料初始剪切模量的值。但是在使用umat時，Abaqus對程序輸入文件進行預(yù)處理時得不到剪切模量的數(shù)值，所以這時候必須通過hourglass stiffness定義具有沙漏模式的單元的沙漏控制剛度。
在Abaqus/CAE 的單元模塊即可定義沙漏控制模式，如圖所示。
abaqus子程序Umat基礎(chǔ)知識及實例解釋
UMAT子程序具有強大的功能，使用UMAT子程序：
(1)可以定義材料的本構(gòu)關(guān)系，使用ABAQUS材料庫中沒有包含的材料進行計算，擴充程序功能。
(2)幾乎可以用于力學(xué)行為分析的任何分析過程，幾乎可以把用戶材料屬性賦予ABAQUS中的任何單元；
(3)必須在UMAT中提供材料本構(gòu)模型的雅可比（Jacobian）矩陣，即應(yīng)力增量對應(yīng)變增量的變化率。
(4)可以和用戶子程序“USDFLD”聯(lián)合使用，通過“USDFLD”重新定義單元每一物質(zhì)點上傳遞到UMAT中場變量的數(shù)值。
今天跟大家分享一篇關(guān)于用戶子程序的基礎(chǔ)知識總結(jié)以及實例講解，我覺得是我看過眾多資料里面講解的比較好的，三年前收集的，實在不知道作者是誰，只能先謝謝了。希望對你們有幫助，因為umat應(yīng)該說是abaqus跟其他有限元軟件的最大不同之處了，是每個人學(xué)abaqus到最后無法避開的必經(jīng)之路，希望大家學(xué)好。其他也整理了好多學(xué)習(xí)Umat的資料，有需要的各位找我自取。
下面是pdf正文：
1
2
3
4
5
6
7
8
9
10
11
也算是自己之前學(xué)習(xí)的一點總結(jié)，可惜自己寫不了這么好。加油吧。內(nèi)容的pdf我也已經(jīng)上傳到了技術(shù)鄰。
其他的一些資料如下，好久之前整理的，有需要的請私聊自取哈。
—End—
發(fā)布于 00:18
展開 abaqus子程序Umat基礎(chǔ)知識及實例解釋
UMAT子程序具有強大的功能，使用UMAT子程序：
(1)可以定義材料的本構(gòu)關(guān)系，使用ABAQUS材料庫中沒有包含的材料進行計算，擴充程序功能。
(2)幾乎可以用于力學(xué)行為分析的任何分析過程，幾乎可以把用戶材料屬性賦予ABAQUS中的任何單元；
(3)必須在UMAT中提供材料本構(gòu)模型的雅可比（Jacobian）矩陣，即應(yīng)力增量對應(yīng)變增量的變化率。
(4)可以和用戶子程序“USDFLD”聯(lián)合使用，通過“USDFLD”重新定義單元每一物質(zhì)點上傳遞到UMAT中場變量的數(shù)值。
今天跟大家分享一篇關(guān)于用戶子程序的基礎(chǔ)知識總結(jié)以及實例講解，我覺得是我看過眾多資料里面講解的比較好的，三年前收集的，實在不知道作者是誰，只能先謝謝了。希望對你們有幫助，因為umat應(yīng)該說是abaqus跟其他有限元軟件的最大不同之處了，是每個人學(xué)abaqus到最后無法避開的必經(jīng)之路，希望大家學(xué)好。其他也整理了好多學(xué)習(xí)Umat的資料，有需要的各位找我自取。下面是pdf正文：
也算是自己之前學(xué)習(xí)的一點總結(jié)，可惜自己寫不了這么好。加油吧。pdf也已經(jīng)上傳到了技術(shù)鄰，大家可以下載。
其他的一些資料如下，好久之前整理的，懶得一個個區(qū)整理了，有需要的請私聊自取哈。
—End—
展開 