Abaqus子程序開發的案例
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科技工作室。
展開 abaqus子程序開發：fortran-c-python混合編程開發
abaqus子程序開發語言不局限fortran和c了，還能支持python！
這是一套自研的開發框架，源于客戶需要獲取分析過程的中間數據作為python機器學習程序的輸入這樣的需求。這種需求乍一聽確實非主流、非常規。
后來仔細做了一些了解，這種需求有其合理性。第一，python在數據分析、機器學習方面有豐富的資源；第二，python程序嵌入abaqus求解器一起運行，那開發就會很方便高效；第三，借助python的pdb可以隨心所欲的設置斷點，調試程序。在此之前，python一般作為獨立程序，依賴數據文件的順序傳遞，實現與abaqus的協同。直白的講，abaqus算完了，形成數據文件再導入python。
于是火力全開研究解決方案，結果可謂“大快人心”：實現了abaqus的子程序對python的支持，而且還實現了fortran、c、python三種語言的混合編程，真正讓不同的編程語言在它最擅長的領域發揮作用。
以上介紹了abaqus子程序fortran-c-python混合編程，如有這方面的需求，歡迎私信聯系開展合作。
展開 abaqus 材料子程序開發
小妹最近想學習abaqus材料子程序開發，請問一下高手，需要哪些軟件方面的支持呢，是不是安裝visual studio和fotran就可以了，我的abaqus是6.7的，另外兩個軟件需要什么版本的呢。多謝多謝
ABAQUS子程序二次開發的軟件配置
由于工程問題的復雜度，Abaqus通用有限元軟件在求解某些專業問題上有所欠缺。針對此問題，Abaqus提供了二次開發的功能。 Abaqus的二次開發可以拓展Abaqus的應用范圍，提高對問題求解的復雜度和精確度。
Abaqus提供了兩種二次開發的接口：一種是子程序接口(User Subroutine)，一種是腳本接口(Abaqus Scripting Interface)。腳本語言是在Python語言的基礎上進行的定制開發，它擴充了Python的對象模型和數據類型，使Abaqus腳本接口的功能更加強大。子程序接口開發是在Fortran語言的基礎上進行的定制開發，可以根據需要開發新的材料模型，單元模型以及與其他應用軟件的數值交換。本文主要講述利用Fortran語言對Abaqus進行子程序開發的軟件配置。
Abaqus子程序的二次開發需要配置IVF和VS。做好相關設置后可通過本文提供的案例進行驗證以證明軟件的配置是否正確。
展開 
Abaqus通過USDFLD子程序進行泥巖的應變軟化模擬
本文基于泥巖的三軸壓縮試驗曲線，建立考慮應變軟化特性的泥巖彈塑性本構模型，使用Abaqus及其子程序對泥巖的三軸壓縮試驗進行了數值模擬。泥巖在受壓過程中主要經過了5個階段，即壓密階段、彈性變 形、應變硬化、應變軟化、殘余階段。泥巖應變軟化模型如下所示。
式中，ξ為強度參數，ξp為峰值強度參數，ξr為殘余階段強度參數，η為應變軟化參數，η*為殘余階段的應變軟化參數初始值。對于三軸壓縮試驗，η用塑性剪切應變來表示
塑性屈服準則采用Mohr-Coulomb準則，則粘聚力和內摩擦角的參數演化可以用下式表示
在巖石的塑性變形過程中會產生比較明顯的剪脹現象，而用來描述這一現象的較常用的力學參數就是剪脹角 Ψ，
上述模型可以通過USDFLD子程序進行實現，流程圖如下
有限元模型如下圖所示
計算得到的應力云圖及不同圍壓下的載荷位移響應如下所示
參考文獻：張力偉,賈善坡,鄒江濤,舒婧曦.泥巖的峰后軟化力學模型.中國科技論文,2016,11(21):2456-2461
有關于abaqus子程序開發的相關問題可以通過公眾號聯系我們.
公眾號: 320科技工作室
展開 abaqus二次開發子程序disp簡單案例 ￥5
abaqus二次開發子程序disp簡單案例，能夠快速了解disp子程序的組織架構。利用disp子程序使abaqus物體表面隨時間移動，即位移為時間的函數。
關鍵操作：
1、給物體面施加位移邊界條件——位移/轉角——選擇面——分布：用戶定義
2、新建作業——通用——添加disp.for文件——運行
Abaqus通過USDFLD子程序進行泥巖的應變軟化模擬
本文基于泥巖的三軸壓縮試驗曲線，建立考慮應變軟化特性的泥巖彈塑性本構模型，使用Abaqus及其子程序對泥巖的三軸壓縮試驗進行了數值模擬。泥巖在受壓過程中主要經過了5個階段，即壓密階段、彈性變 形、應變硬化、應變軟化、殘余階段。泥巖應變軟化模型如下所示。
式中，ξ為強度參數，ξp為峰值強度參數，ξr為殘余階段強度參數，η為應變軟化參數，η*為殘余階段的應變軟化參數初始值。對于三軸壓縮試驗，η用塑性剪切應變來表示
塑性屈服準則采用Mohr-Coulomb準則，則粘聚力和內摩擦角的參數演化可以用下式表示
在巖石的塑性變形過程中會產生比較明顯的剪脹現象，而用來描述這一現象的較常用的力學參數就是剪脹角 Ψ，
上述模型可以通過USDFLD子程序進行實現，流程圖如下
有限元模型如下圖所示
計算得到的應力云圖及不同圍壓下的載荷位移響應如下所示
參考文獻：張力偉,賈善坡,鄒江濤,舒婧曦.泥巖的峰后軟化力學模型.中國科技論文,2016,11(21):2456-2461
有關于abaqus子程序開發的相關問題可以聯系扣扣1653004885或者關注cae320公眾號
展開 【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用戶子程序安裝、開發、調試
引言
Abaqus為用戶提供了大量的用戶子程序，用戶可以通過子程序接口開發Abaqus原本不具備的單元、材料、載荷等。此外，利用子程序能幫助我們去理解、推測Abaqus的計算過程，這對國產自主軟件對標Abaqus功能開發具有一定的參考意義。
Abaqus 子程序開發日志01-umat之isotropic
對有限元進行了一些理論上的學習，現在正式嘗試對abaqus用戶子程序進行嘗試書寫，此貼只作為自己學習的記錄以及分享。
-------------------------------------------------------------------------------------------------------------------------------
umat子程序只用于通用靜力分析步使用來定義用戶自定義材料，這里首先對較為簡單的各向同性材料本構模型進行構建。
先說明一下umat的工作原理，再靜力分析時，每增加一個（極小的）微分應變增量，通過代入微分應力關于微分應變的雅可比矩陣B（j）方程，得到應力增量，再與原始應力累加得到新的應力（也可以說更新應力）。至此，此umat的工作到此結束。
對于各向同性的剛度矩陣方程為：
再來就是對矩陣中未知量進行求解，對其而言只需兩個獨立常量就可以表示。對于各向同性材料而言，G=E/2(1+NU),所以我們只要知道G、E、NU中的任意兩個量就可以得到另外一個量，習慣上我們常用E、NU來表示各向同性材料的材料屬性，所以我們這里就需要將矩陣中的未知量用E與NU表示出來就可以了。
求解的過程不再說明，此為材料力學的基礎，有需要可以查找相關資料。
求解完成后將其改寫成微分的形式得到雅可比矩陣后用程序進行表達出來，并使用abaqus提供的接口進行嵌套寫入用戶子程序，至此，全部工作完成。
先對其進行檢驗，abaqus中可以直接定義各向同性材料的材料屬性，在初期子程序的學習中檢驗umat的正確性提供了很好的對照。
這里就隨便建立一個長寬高為5；5；20的立方體對其進行施加一個單位為5的拉伸位移（較長方向拉伸）。
展開 Abaqus 修正劍橋模型的vumat子程序開發
考慮到莫 爾–庫侖準則有此特征，故將莫爾–庫侖準則與劍橋模型相結合，得
聯合式(1)，(4)，(5)，有改進的屈服函數：
根據上述理論，可以通過以下流程圖完成VUMAT子程序編寫
通過單胞模型計算得到得結果如下圖所示。
[1]袁克闊,陳衛忠,于洪丹,譚賢君,趙武勝,李香玲.考慮黏聚特性和拉壓不等效應的修正劍橋模型及數值實現[J].巖石力學與工程學報,2012,31(08):1574-1579.
Abaqus 子程序開發日志05-umat之fully anisotropic
此umat的材料參數設置參考abaqus，參數自上而下分別為：
使用方法和之前一致。
若選用縮減積分單元時要開啟沙漏控制。
完全各向異性umat：
fullyanisotropic-std.obj
共享庫：
standardU.dll
abaqus二次開發UNAT子程序實例 ￥12.4
二次開發的源程序及子程序可在“收費內容”中的附件中下載，編譯的話用txt文本就可以，編完了改一下后綴，將.txt改成.f然后直接調用子程序就可以了。需要注意的是你的電腦安裝了ABAQUS、VS和IVF
Abaqus 子程序開發日志07-vumat之anisotropic
正交各向異性（工程常數）：
EngineeringConstants-xplD.obj
EngineeringConstants-xpl.obj
共享庫：
explicitU.dll
explicitU-D.dll
參數輸入依次為：
E1,E2,E3,Nu23,Nu13,Nu12,G23,G13,G12
正交各向異性（剛度系數）：
OrthogonallyAnisotropic-xpl.obj
OrthogonallyAnisotropic-xplD.obj
共享庫：
explicitU.dll
explicitU-D.dll
?參數輸入依次為：
D1111,D1122,D2222,D1133,D2233,D3333,D1212,D1313,D2323
（按照σ1，σ2，σ3，┏12，┏13，┏23的形式輸入）
橫觀各向同性：
TransverselyIsotropic-xpl.obj
TransverselyIsotropic-xplD.obj
共享庫：
explicitU.dll
explicitU-D.dll
?參數輸入依次為：
Et(軸/橫向彈性模量）、Ep（面內彈性模量）、Nutp（面外泊松比，也就是說假設1為軸向，那么Nu12/13的值為Nutp）、Nup（面內泊松比）、Gt（軸/橫向剪切模量）
（這里的軸/橫向對abaqus來說是局部方向的1方向）
完全各向異性：
fullyanisotropic-xpl.obj
fullyanisotropic-xplD.obj
展開 abaqus剪切修正GTN模型的VUMAT子程序開發
這里對相應的算法進行簡要說明：
NH-GTN模型
屈服函數：
其中等效孔洞體積分數定義為：
孔洞體積分數包含新孔隙形核，原有空隙生長以及剪切相關的等效體積分數增加：
形核，生長，剪切相關體積分數的演化遵循：
其中：
剪切效應的修正，考慮應力狀態的影響
參數的物理含義如下
通過將文獻中的數值算法編程實現在VUMAT子程序中，可以用來實現對延性金屬材料在不同應力狀態下的損傷演化進行合理的數值預測，應用于金屬成型領域（沖壓，軋制，擠壓等）
預測修正后的模型應該在簡單拉伸情況下于abaqus自帶的GTN模型保持相同的損傷和其他狀態變量的分布，并在剪切情況中損傷發展顯著高于abaqus自帶的模型（自帶的模型忽略了剪切效應）。（為了進行對比使用于自帶的本構相同的硬化方式，模擬中使用了相同的質量縮放，但質量縮放容易產生數值振蕩，模擬的拉伸曲線存在波動。）
初步模擬結果：
拉伸情況（abaqus-VUMAT）
應力分情況
孔洞體積分數
剪切模型（abaqus-VUMAT）
不同變形時刻的應力分布
T=0.1s
局部放大圖
T=0.5s
局部放大圖
T=0.6s
局部放大圖
可以看到模型在拉伸預測中與原始模型保持一致，而在剪切修正后損傷發展顯著快于原始模型，利用作者提出的方法可以應用于復雜應力狀態下金屬材料的損傷分析，相關參數部分參考文獻，其中Kw=3.T1=0.2，T2=0.7.模擬結果符合文獻所提出方法的基本趨勢。
最后，如果有相關需要歡迎通過公眾號“320科技工作室”聯系我們。
展開 