ABAQUS UMAT的案例
初識ABAQUS UMAT二次開發（二）
圖5 從abaqus計算結果中導出的豎桿應力-應變曲線
對比圖5和圖2可知，模型中豎桿的應變-應變曲線與輸入的本構關系完全一致，可見UMAT程序正確。
（本示例UMAT源碼、CAE模型均可在公眾號QQ群文件中免費獲得）
4
結語
Abaqus UMAT用戶子程序接口功能非常強大，編寫一個漂亮的UMAT子程序需要一定的力學基礎（尤其是有限元、連續介質力學）。本文旨在讓大家對UMAT有一個初步的了解，以及使用UMAT需要用到哪些基本知識，好在今后的學習中有側重點。后續的推文中，小編還會進一步結合我們有限元理論模塊的進度，進一步講解Abaqus UMAT的使用。
展開 ABAQUS UMAT for cohesive element compression ￥150
本文和我的另一篇文章ABAQUS UMAT for cohesive element tension behaviour 比較相像，所以不重復介紹abaqus建模過程，主要介紹理論和子程序。
本文利用 cohesive element 模擬兩個粘接的物體受剪破壞，利用的本構模型是compression cut-off模型（和tension cut-off模型原理一樣，只是cut-off的面在compression段），有限元實現算法是帶有返回功能的前歐拉法（Crisfield在1991出版的Non-linear finite element analysis of solid and structure中的第6章有提到），軟件是abaqus，子程序為UMAT。
這種有限元模型主要應用在在膠合破壞的預測，比如磚結構的水泥砂漿、纖維復合材料加固結構中復合材料的剝落、多層玻璃的脫膠以及夾層板的滑移等等。可以說，與膠有關的結構都可以試著利用cohesive單元模擬，所以開發cohesive單元的本構模型是很多領域的熱點，因為相比起固體單元，cohesive單元的發展相對較晚。
理論：請先看我的另一篇文章ABAQUS UMAT for cohesive element tension behaviour，然后在結合下文。
compression cut-off返回算法：
abaqus模擬的結果：
最大主應力：
最大塑性應變：
完整視頻模擬結果詳見：https://zhuanlan.zhihu.com/p/113157257
展開 初識Abaqus UMAT 二次開發
使用UMAT必須在子程序中定義材料本構模型的雅克比矩陣，即應力增量對應變增量的變化率。
在Abaqus軟件中UMAT需要使用Fortran語言（一般都采用Fortran77語法、格式）進行編寫。在WIN10操作系統下需要配置好Abaqus、intel Fortran、Visual Studio的連接，在B站、技術鄰上有很多免費的教程，在此不再累述。
【7月20-21日 上海】Abaqus UMAT用戶子程序二次開發技術培訓
理論和實際操作結合，一次性將abaqus UMAT用戶子程序開發課程講解透徹！
至于python的自動化和GUI部分內容的操作練習，如有此方面需求，將有老師帶領學習。
背景介紹：
Abaqus是當今最流行的通用結構有限元軟件之一，以其強大的非線性和良好的用戶體驗著稱。但在實際工程應用中，設計對象專業針對性強，通用的工具軟件往往解決不了專業問題，需要通過二次開發對商業軟件進行定制和優化，集成高精度的計算方法和專有技術，以更好的解決專業問題，甚至形成自主品牌的CAE軟件。
ABAQUS不僅提供標準的有限元分析程序，而且具有良好的開放性，可利用它提供的用戶子程序接口生成非標準的分析程序來滿足用戶的需要，在實際工程中得到廣泛應用。
Abaqus UMAT是ABAQUS提供給用戶定義自材料屬性的程序接口。通過UMAT，用戶可以定義材料的本構關系，建立Abaqus材料庫沒有的材料本構模型，滿足實際工程分析的需要。
本次培訓主要包括：
（1） Abaqus二次開發整體介紹；
（2） UMAT接口原理和使用方法；
（3） 材料本構模型和基礎理論
（4） 基于Fortran和Matlab的UMAT編寫示例和練習。
通過本次培訓，學員可以：
（1） 學員在老師指導下完成用戶子程序開發環境配置；
（2） 建立UMAT編寫相關的理論基礎；
（3） 可自行從零開始構建UMAT算例；
（4） 學員實際項目問題現場交流；
（5） 和專家及其它行業用戶深入溝通，并獲得長期的技術支持。
培訓方式：
培訓主要采用：
（1） PPT介紹；
（2） 老師操作演示；
（3） 學員上機練習。
三者結合的方式直觀的學習各個知識點。
講師介紹：
技術鄰ID：SnowWave02 團隊，3個講師同時授課。
展開 
ABAQUS UMAT-混凝土受拉狀態下塑性損傷模型的簡單實現 ￥600
本文利用ABAQUS UMAT子程序，簡單實現了混凝土受拉狀態下的破壞。本構模型的實現算法摘抄自DeBorst的書籍《Nonlinear Finite Element Analysis of Solids and Structures》，基本如下：
為了簡化模型，筆者將書中損傷部分做了簡化，不再采用損傷屈服面進行判定。損傷影子w的計算直接由塑性等效應變確定。
在ABAQUS中建立100*100*100的立方體塊，試件的底部固定，頂部反復加載-卸載，通過UMAT得到的模擬結果如下：
Abaqus中UMAT二次開發詳細教程由入門到精通
一．Abaqus中UMAT的使用背景
ABAQUS中雖然提供了很多材料本構模型，但是有時候我們的分析的實際情況是比較復雜的，比如當我們要計算一種特定的橡膠材料的時候，已有的橡膠本構比如Mooney等都有著這樣、那樣的不足，很難滿足我們的要求。又比如我們要計算土壤應力的時候，土壤本身受壓不受拉，有著很強的非線性因素，這時候就需要我們建立足夠精確的土壤本構模型。如果你正在處理的問題，ABAQUS的常規本構不能滿足你的要求，那么UMAT將是你的首選。
ABAQUS UMAT UEL子程序咨詢解疑 ￥400
ABAQUS UMAT UEL子程序咨詢解疑
有限元理論基礎及Abaqus內部實現方式研究系列9： 編寫線性UMAT Step By Step
介紹六面體單元的基礎理論和Abaqus中C3D8R六面體單元的剛度矩陣的修正方式，采用這些修正方式可以得到和Abaqus六面體單元完全一致的剛度矩陣。
https://www.yqgqt.org.cn/content/post/430177
第八篇：UMAT用戶子程序開發步驟。介紹基于Fortran和Matlab兩種方式的Abaqus的UMAT的開發步驟，對比發現開發步驟基本相同，同時采用Matlab更加高效和靈活。
https://www.yqgqt.org.cn/content/post/432848
展開 ABAQUS中UMAT中的經典循環塑性本構模型及相應子程序代碼 ￥20
ABAQUS中UMAT中的循環塑性模型，包含非線性各向同性強化彈塑性、線性各向同性強化彈塑性、線性隨動強化彈塑性模型，包含CAE文件、UMAT文件等。
ABAQUS二次開發-UMAT入門
ABAQUS中UMAT使用背景
ABAQUS中雖然提供了很多的單元，很多的材料本構模型，但是有時候我們的分析的實際情況是比較復雜的，比如當我們要計算一種特定的橡膠材料的時候，已有的橡膠本構比如Mooney等都有著這樣、那樣的不足，很難滿足我們的要求。又比如我們要計算土壤應力的時候，土壤本身受壓不受拉，有著很強的非線性因素，這時候就需要我們建立足夠精確的土壤本構模型。如果你正在處理的問題，ABAQUS的常規本構不能滿足你的要求，那么UMAT將是你的首選。
(2)
UMAT參數說明
UMAT是一個較為復雜的程序，涉及到的參數較多，首先我們要明白UMAT完成了什么功能，簡單的說，就一句話，開發者需要定義材料積分點的Jacobian矩陣。由開發者完成應力-應變之間的關系曲線。下面我們對一些參數做一些介紹，其中kstep為分析步參數、kinc為增量步參數、time為主程序傳入進來的分析時間參數noel為單元編號、npt為單元積分點編號。stran為應變。dstran為應變增量。ddsdde為最關鍵Jacobian矩陣。
下面我們以最簡單的等溫線彈性本構為例與ABAQUS自帶的本構進行對比計算。我們要建立的模型如下圖所示：
在分析中我們要保證兩個本構的邊界條件一致、載荷一致、網格一致。我們所寫的本構這里給出主要內容,如下所示：
使用命令行模式提交計算，兩個本構計算結果對比如下所示：
由上兩圖可知，兩個本構計算的結果最大應力都為6.552Mpa。由此我們可知在ABAQUS中其自帶的等溫彈性本構一定也是如此。
展開 Abaqus使用umat子程序的沙漏問題
在使用umat子程序時，采用縮減積分單元后，沙漏控制剛度是通過材料屬性中的彈性性質定義的，這些剛度基于材料初始剪切模量的值。但是在使用umat時，Abaqus對程序輸入文件進行預處理時得不到剪切模量的數值，所以這時候必須通過hourglass stiffness定義具有沙漏模式的單元的沙漏控制剛度。
在Abaqus/CAE 的單元模塊即可定義沙漏控制模式，如圖所示。
Abaqus子程序代碼分享
Engineering Fracture Mechanics, 185: 210-226 (2017)
17、UELcode.zip
18、ABAQUS UMAT SUBROUTINE FOR CMSG PLASTICITY
ABAQUS UMAT子程序用于CMSG可塑性
ABAQUS User Material (UMAT) subroutine with the constitutive formulation of the conventional mechanism-based strain gradient (CMSG) plasticity theory.
ABAQUS用戶材料（UMAT）子例程，具有常規基于機理的應變梯度（CMSG）可塑性理論的本構公式。
Paper: E. Martínez-Pa?eda and C. Betegón. Modeling damage and fracture within strain-gradient plasticity. International Journal of Solids and Structures, 59: 208-215 (2015)
18、CodeCMSG.zip
展開 利用Abaqus的UMAT子程序仿真木材蠕變現象
下面我們將以木材蠕變為例，介紹下在ABAQUS UMAT中如何實現蠕變仿真。
2. 本構理論
文獻[1]給出了木材蠕變過程中本構：
3. 算例
3.1 模型
考慮懸臂梁模型，如下圖。
3.2 邊界條件
根據蠕變的定義，模型必須現有一個穩定的載荷，因此可以分成兩個step。第一個step，完成力加載，第二個step保持載荷，實現蠕變變形的生長。
需要指出的是，蠕變通常需要在較長的時間尺度上才能有明顯的效果。比如我們要觀察100天的變形情況，那么這個時候，ABAQUS設定的總時間還是1，在UMAT里面要乘以相應的系數，給出物理時間，才能有效的實現蠕變效果。
3.3 結果
最終得到懸臂梁端部位移如下圖所示。從圖中可以看出，在最初的幾天，蠕變變形較大，隨著時間的推移，變形增加的幅度放緩，符合蠕變的特點。
參考文獻
[1] 《濕度變化和荷載作用下膠合木曲梁的工作性能研究》
[2] 《旋切板膠合木的蠕變及其對結構穩定性的影響》
最后，有abaqus二次開發相關需求歡迎大家通過我們的微信公眾號聯系我們。
微信公眾號：320科技工作室
展開 下加載面修正劍橋模型及其在Abaqus中umat子程序的實現
Abaqus的umat子程序實現
子程序編寫流程如下所示：
三。模型驗證
下面給出一個簡單的應用 UMAT 子程序的算例，將模型劃分為多個單元， 分別采用位移的加載方式和力的加載方式，并分別與試驗數據和 ABAQUS 自帶 的修正劍橋模型計算結果進行對比，以此說明 UMAT 子程序可以應用于多個單 元的模型，對于位移和力的加載方式均適用，并證明 UMAT 子程序適用于超固 結土體的力學特性分析和地基的固結沉降分析。
從而證明了編寫的子程序的有效性和實用性。
最后，大家有關于Abaqus二次開發的任何需求可以添加管理員微信號：CAE320，同時也歡迎大家關注“320科技工作室”的微信公眾號，掃一掃二維碼即可關注~~
展開 Abaqus——UMAT學習
本來想制作一個入門級abaqus子程序UMAT的程序，線彈性材料比較簡單。講起來方便，明白。確實是這樣。那我就先弄一個這個來說明一下。 我這里就簡單說一下2D情況，你把剛度矩陣DDSDDE直接賦值，然后求一下應力STRESS，值得注意的是 應力的順序：sigma（i）=D（i，j）*E(j) ,遍歷i，j=1,3。我們知道小變形abaqus給出的是Cauchy應力，所以是對稱的。它的排布是先主應力再剪應力。 你可能需要 一個fortran入門的教材，我有一個，但是由于版權，我不上傳。需要的我私發一下（只發前十份）。供大家學習。
Mooney_Rivlin Material Model in abaqus 是我應該寫的，應該是這周末吧。到此
展開 