熱力耦合本構開發的案例
Abaqus形狀記憶合金熱力耦合本構模型umat子程序分享
This manual provides a comprehensive review of user subroutine SMA_UM (Shape Memory Alloy — Unified Model). SMA_UM is a FORTRAN coded numerical implementation of an SMA thermomechanical constitutive model. The current version of SMA_UM implements the unified constitutive model presented by Lagoudas et al. (1996), which unifies and generalizes to three dimensions the constitutive models presented by Tanaka (1986), Boyd and Lagoudas (1996) and Liang and Rogers (1992). The model is numerically implemented using return mapping algorithms. Detailed description of the implementation is provided in the research paper by Qidwai and Lagoudas (2000). The current implementation of the subroutine follows the specifications for user-material subroutines by ABAQUS and is intended primarily for use
展開 巖土非線性本構UMAT開發 ￥666
巖土非線性本構UMAT開發
晶體塑性耦合連續損傷本構框架
經典文章推薦
《Necking behavior of AA 6022-T4 based on the crystal
plasticity and damage models
是最經典的耦合晶體塑性理論和連續損傷的文章之一，損傷力學有兩種主要方法。第一種是Gurson提出的基于微觀力學的損傷模型。在基于微觀力學的方法中，損傷演化通過孔隙成核、生長和聚結來描述。對空穴成核和生長進行了建模，必須使用實驗數據確定相關系數。另一種方法是連續損傷力學（CDM）。在CDM框架中，使用應力、壓力、溫度和應力三軸性確定斷裂應變。在這些研究之后，提出了許多改進的模型，以包括洛德角和各向異性損傷的影響，
作者在研究中使用的損傷模型基于連續損傷力學（CDM）。然而，通過結合CPFEM可以預測孔隙的萌生、生長和聚結行為。此外，材料因損傷而弱化用于描述頸縮后承載能力的突然下降，通過顯式時間積分方案進行了分析，這為通過CPFEM預測頸縮行為提供了可能性。然而，沒有預測頸縮形狀和載荷位移曲線。為了準確預測頸縮和載荷位移曲線，使用隱式時間積分方案進行了分析，可以獲得更合理的載荷位移曲線。此外，還進行了實驗，并與分析結果進行了比較。最后，新提出了四種不同的帶系數校準的損傷模型，并提出了一種最能描述頸縮行為的模型。
作者使用的四類連續損傷模型理論如下
（1）最大塑性應變損傷模型：該模型將損傷定義為當主塑性應變大于某一臨界值時開始和累積的損傷。此模型寫為：
ε1f.ini是損傷萌生塑性應變值，ε1f.ini是最大塑性應變值，D是損傷因子，M是損傷指數（通常取值大于1.0有利于流動應力平滑過渡）
（2）等效塑性應變損傷模型：該模型將損傷定義為當等效塑性應變大于某一臨界值時開始和累積的損傷。
展開 lsdyna材料本構模型二次開發經驗分享（umat41）
二、lsdyna二次開發前期儲備
2.1 軟件安裝（最基本最簡單）
（1）必須擁有對應版本的ls-dynalib文件包。需要根據不同的系統平臺（32位/64位、SMP/MPP）、版本（971R5.1.1/R711/R8.0.0）得到對應的lib包。
（2）裝IFC之前需要先裝MVS
（3）有lstc授權文件，即ls‐dyna求解器可求解。
下面是版本對應關系：
2.2 力學知識儲備（最難）
因為做lsdyna本構模型二次開發和直接使用內置本構進行計算難易程度差距很大，采用內置本構進行計算不用過多了解本構底層邏輯。而二次開發自己的本構需要對整套本構的內在邏輯有很好的理解，包括基本的應力應變關系、應力偏量、靜水應力、應變率、應力不變量、應力偏量不變量等等。如果有做二次開發的打算，建議提早學習一下以上知識點，可以查閱相關書籍，個人建議如果想速成可以在B站上去學習，有一位女老師講的線上網課非常受用。
2.3 Fortran語言基礎（相對較容易）
二次開發對于編程的要求是很低的，只需要掌握最基本的即可。比如用到最多的條件語句里的比較：gt為大于、le為小于等于、ge為大于等于。其他的也都和上邊這些最基本的一樣，在具備以上兩點之后可以在lsdyna手冊里看一下具體代碼，先讀一遍，主要是學習套路和編程語言。
如果以上三點你都基本具備了，那么就可以自己嘗試去根據推導的本構去編一下umat代碼了。
如果umat代碼寫完了，那么恭喜你，可以進入下一道難關了：編譯生成求解器。如果你是初學者，那么寫完之后的代碼肯定會錯誤百出，不過不要慌，這很好解決。因為在編譯的時候如果你的語法有錯誤，他都會提示的，在vs里打開行號就可以清晰地看到具體哪一行出錯了，fortran語法以及umat書寫格式問題就可以通過一次次改錯而解決掉了。
展開 
木材三維彈塑性本構子程序開發
問題介紹
木材的本構模型是采用連續體單元建模模擬木材彈塑性響應的基礎，然而木材復雜的力學性質常常為其本構模型的建立帶來困難。木材力學性質的復雜性主要表現在：
不同方向的強度值和剛度值各不相同；
2. 同一方向的抗拉強度和抗壓強度之間存在差異；
3. 不同形式荷載作用下材料的響應不同，壓力作用下材料的表現以延性為主， 而拉力和剪力作用下材料的破壞呈脆性。
木材在復雜應力狀態下的彈塑性本構模型。以經典彈塑性力學為框架，該本構模型建立在如下四個基本假設的基礎之上：
木材在彈性階段是理想的橫觀各向同性材料；
2. 材料的屈服符合簡化的 Hashin 屈服準則；
3. 材料在受拉和受剪屈服之前是理想線彈性的，屈服之后進入塑性流動階段；
4. 材料受壓初始屈服之前是理想線彈性的，屈服之后進入應變硬化階段，隨 著屈服面的轉移到達最終屈服面后進入完全的塑性流動。
二。子程序編寫流程
本工作室在三維hashing模型的基礎上，利用Abaqus軟件平臺，開發了完整的木材的彈塑性本構umat子程序，包含木材完整的彈性、塑性、強化以及軟化階段。編寫子程序的流程如下：
三。結果驗證
通過如下圖的木材模型進行驗證：
該模型在受壓、受剪及受拉的工況下，應力應變曲線如下所示：
該子程序還有以下特征：
能計算靜力非線性
2. 收斂性好
3. 能計算復雜應力狀態
附件為本子程序參考的文獻，供大家學習探討~
2. 木材的力學性質試驗研究及數值模擬方法.pdf
最后，大家有關于Abaqus二次開發的相關需求可以添加管理員扣扣：3045552826，微信：CAE320，同時也歡迎大家關注“320科技工作室”的微信公眾號，掃一掃二維碼即可關注~~
展開 flac3d的interface界面單元非線性本構模型開發代碼 ￥12
<p>基于fish語言的flac3d的interface界面單元非線性本構模型開發實例</p>
基于VB的ANSYS二次開發之Duncan-Chang本構模型算法介紹
所謂的本構模型是指材料的應力應變關系，比如遵循虎克定律的線彈性應力應變關系就是一種常用的本構模型。ANSYS為用戶提供了許多本構模型，但在某些特殊領域，現有的本構模型卻很少，完全不能滿足分析要求。為了解決這個問題，ANSYS為用戶提供了usermat等UPFs用戶子程序，這些用戶子程序擁有強大的二次開發功能，可以實現各種復雜本構模型的開發。但是，對于一些簡單的本構模型，用戶也可以利用APDL語言進行開發，比如Duncan-Chang本構模型。
Duncan-Chang本構模型介紹
對于應力應變關系復雜的材料而言，變形主要是由顆粒間位置的變化引起的，不同應力水平下相同的應力增量引起的應變增量并不相同，從而表現出復雜的非線性特征。為了反映材料特性，人們提出了許多本構模型，以期更好的反映材料的應力應變規律。
Duncan-Chang本構模型主要優點是可以利用常規三軸固結排水剪試驗確定模型參數，因此在工程實踐中得到了廣泛應用。Duncan-Chang本構模型是非線性彈性模型，彈性矩陣中的彈性系數不再是常量，而是隨應力狀態而改變。由于不考慮塑性變形，所以一般只適用于載荷不大的情況（即不太接近破壞的情況）。Duncan-Chang模型有E-V和E-B模型兩類。
當然，該模型也存在一些缺陷，如無法反應不同應力路徑的影響、加卸載判斷不明確等，不可避免造成了計算分析誤差，長期以來許多學者試圖來對其進行修正。有限元軟件通常采用傳統塑性理論的應力符號，以拉應力為正，下面對Duncan-Chang模型采用有限元軟件的應力形式進行說明。
Duncan-Chang本構模型算法
該模型是Duncan和Chang根據Konder關于巖土材料的三軸試驗的偏應力與軸向應變近似呈雙曲線的假定而提出的。
展開 ABAQUS增材制造AlSi10Mg修正JC本構二次開發（論文復現）
[圖片]
Abaqus考慮拉壓不對稱的樹脂彈塑性損傷本構vumat子程序開發
復合材料在航空航天領域的應用非常廣泛。在研究復合材料失效機理的過程中，學者們提出了許多宏觀和細觀尺度上的失效準則。其中應用比較廣泛的宏觀失效準則包括hashin準則、puck準則等。復合材料在細觀尺度上的失效行為通常通過代表體積單元（RVE）模型來研究。
RVE代表體積單元
RVE模型由纖維和樹脂構成，一般假設纖維是橫觀各向同性線彈性材料，樹脂則為彈塑性材料。本文通過在屈服準則中引入拉壓非對稱參量，研究了樹脂的拉壓不對稱彈塑性損傷行為。
由于樹脂的屈服行為與靜水壓力相關，這里采用下式所示的拋物面屈服準則。
式中J2為偏應力的第二不變量，I1為應力第一不變量，σt和σc為拉壓屈服應力
采用非關聯塑性流動準則，如下所示。
式中σvm為mises等效應力，P為靜水壓力，α為材料參數
損傷萌生準則如下所示
式中J2和I1為無損應力下的不變量。
為了降低模型的網格依賴性，損傷演化采用特征長度相關的指數模型
式中，rm為損傷內變量，am為特征長度相關的材料參數。
Melro的文章中給出了通過Simpson積分和弦截法計算Am的方法，實際計算發現通過該方法計算的am效果不是太理想，因此本文未對am進行迭代，直接采用其初值進行仿真計算，如下所示。
計算流程如下
計算流程圖
根據上文的彈塑性損傷模型編寫了vumat子程序，并通過單胞模型進行了驗證，計算結果如下圖所示。
abaqus單胞模型
拉伸載荷下的應力應變曲線
壓縮載荷下的應力應變曲線
展開 建筑黑科技來襲 之Abaqus材料本構庫開發、驗證及工程實踐
下面以韓林海本構模型為例，分別用方/圓試驗對其進行驗證。
鋼管混凝土本構-試驗介紹
韓林海本構方鋼管滯回-數值與試驗對比
韓林海本構方鋼管滯回-滯回動畫
韓林海本構圓鋼管滯回-數值與試驗對比
韓林海本構方鋼管滯回-滯回動畫
方/圓鋼管約束混凝土本構數值計算得到的滯回曲線和試驗滯回曲線相差較小，表明所開發的本構具有良好的計算精度，可用于工程實踐。
3
同濟剪力墻二維損傷本構
Abaqus軟件內嵌了塑性損傷本構，但是該本構是塑型驅動損傷演化滯回規則，容易高估滯回環的耗能，導致計算所得滯回環過于飽滿。
研發室團隊與同濟大學合作，開發了用于剪力墻的二維損傷本構。
展開 ABAQUS基于VUHARD的JC本構二次開發：以二維切削仿真為例（1） ￥600
[圖片]
ABAQUS基于VUHARD子程序的JC本構二次開發：以二維切削仿真為例(3)
[圖片]
ABAQUS基于VUHARD子程序的JC本構二次開發：以二維切削仿真為例（2）
[圖片]
