本構模型開發的案例
flac3d的interface界面單元非線性本構模型開發代碼 ￥12
<p>基于fish語言的flac3d的interface界面單元非線性本構模型開發實例</p>
#PLAXIS#用戶自定義本構模型（UDSM）編譯經驗
很多有限元軟件都提供二次開發接口，供用戶把自己的本構模型編入有限元軟件，PLAXIS也不例外。
不過這是個VIP功功能……目前國內正版用戶本來就不多、VIP用戶更少，其中有二次開發需求的少之又少了。
如果想學習二次開發的話，首先用官方提供的兩個模型的代碼示例玩一玩，將其編譯成軟件可以調用的dll文件。編譯成功后，再考慮自己寫代碼編譯。
總結了一點編譯的經驗，如果有同學用得上，就再好不過了。
PLAXIS 自定義本構模型開發經驗.pdf
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書寫格式問題就可以通過一次次改錯而解決掉了。
展開 基于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的混凝土損傷本構模型與LSDYNA的JHC本構模型分析與研究
圖1混凝土材料本構參數設置
分析：在損傷系數的定義中，應特別注意以下幾點，
1.ABAQUS的混凝土損傷本構模型采用的是非關聯的流動法則，其中系數Dilation Angle，即膨脹角控制著塑性勢函數開口的大小。膨脹角越小，材料越容易破壞，那么相應的結構計算機構就偏向安全，但膨脹角越小就越不容易收斂。因此，膨脹角的取值應當適中，本案例中混凝土本構參數中的膨脹角取值一般在30~35之間，取30。
2.Eccentricity（塑性勢偏移量）決定了塑性勢函數趨近其漸近線的速率。該參數的引入主要是為了保證塑性勢函數的連續、光滑及塑性勢函數在頂點處的可導性。本案例取值0.1。
3.Viscosity Parameter（黏度系數）是為了使材料模型在軟化階段更容易收斂，仍然保持0.1。
3.2基于ANSYS/LSDYNA的混凝土JHC損傷本構模型
對于混凝土材料的本構模型眾多學者進行了深入分析研究以期望獲得一個更加準確描述混凝土材料在壓縮拉伸等力學變化過程中的斷裂行為。除去上述本構損傷模型以外，還有一種專門用來描述混凝土材料的本構模型JHC本構模型。然而，Abaqus自帶的材料模型中并沒有JHC本構，其提供了內置的子程序以供調用。為方便分析進行，本文借助LSDYAN平臺對該本構模型各參數含義進行分析以了解此種本構模型的優勢之處，LSDYNA中對該JHC本構參數的定義界面如圖2所示。JHC本構模型是LSDYNA軟件材料庫中常用于模擬脆性材料的方程之一，尤其是方程中對材料的逐漸累積損傷的計算使得其能夠準確模擬脆性材料的大變形、高應變率效應問題。JHC本構包括應力應變模型、損傷失效模型、靜水壓力模型以及多項式狀態方程[1-2]。
展開 結構分析的常用本構關系（本構模型）
微信 leslie_wj
01 本構關系的定義
本構關系屬于材料的屬性，其實就是材料的應力應變關系，也就是內力和變形的關系。
02 線彈性本構：彈性，并且應力和應變線性相關。
03 非線性彈性本構：彈性，應力和應變非線性相關。
04 理想彈塑性本構：彈性階段應力和應變線性相關，塑性階段應力保持不變。
05 線性強化彈塑性本構：彈性階段應力和應變線性相關，塑性階段應力和應變線性相關。
06 剛塑性本構：彈性可忽略。
一個有意思的材料本構模型設計方案，拉伸變形采用von Mises屈服，壓縮側 cap屈服本構模型設計。
分享這個代碼的主要原因：一方面，它很適合做玻璃、非晶材料、壓痕問題中的壓力敏感塑性分析；另一方面，它也是學習 cap 模型、致密化硬化和隱式本構積分的一個很好的范例。論文結果表明，這一模型能夠較好復現實驗載荷—位移曲線以及壓痕致密化分布，不過需要明確指出的是，當前模型暫時還沒有考慮剪切硬化，因此更適合用于理解“壓痕致密化”這一核心機制，而不是直接覆蓋所有復雜失效問題。作為一份用于科研復現和二次開發的代碼，我覺得它很有參考價值。
約束混凝土cdp塑性損傷本構，mander混凝土本構模型 ￥10
約束混凝土本構，mander混凝土本構，自己做的箍筋約束方柱和圓柱本構模型，表格只要輸入相關參數，自動生成ABAQUS塑性損傷本構關系。
巖土非線性本構UMAT開發 ￥666
巖土非線性本構UMAT開發
YLD2004本構模型 ￥199
其主要應用場景包括以下幾個方面：
車身設計：YLD2004模型可用于描述汽車車身零部件的各向異性塑性行為，進行強度、剛度和疲勞壽命等方面的分析和設計。
金屬成形加工：YLD2004模型可用于描述金屬材料的塑性變形行為，進行成形加工過程的模擬和優化。
建筑結構設計：YLD2004模型可用于描述鋼結構的各向異性塑性行為，進行結構的強度、穩定性和疲勞壽命等方面的分析和設計。
航空航天領域：YLD2004模型可用于描述航空航天結構材料的各向異性塑性行為，進行結構的強度、穩定性和疲勞壽命等方面的分析和設計。
其他領域：YLD2004模型還可用于其他領域的金屬材料塑性分析和設計，例如機械制造、電子設備等領域。
YLD2004模型模型包含18個材料參數用于確定屈服面，以及其他參數去頂硬化和加卸載等：
這些參數需要通過試驗或數值模擬獲得。其中，確定屈服面參數需要進行單軸拉伸、單軸壓縮、剪切等試驗；確定硬化規則參數需要進行多次加載和卸載試驗以測定材料的塑性硬化行為；確定加載面參數需要進行不同方向的應力應變試驗。
而更精確的模擬往往以更高的數值計算為代價，通過原始模型結合線搜索可以實現更快的數值收斂
以YLD2004為本構模型進行單向拉伸加載模擬，結果取下：
展開 ABAQUS umat 非線性等向硬化本構模型（Voce 硬化模型） ￥129
<p class="ql-align-justify">本資源包含一份 PDF 文檔和可直接編譯運行的 Fortran UMAT 代碼，具體內容為：</p><p class="ql-align-justify">非線性等向硬化本構模型（Voce硬化模型） + 隱式積分 + 徑向返回</p><p class="ql-align-justify">完整公式推導 + Fortran 源碼直接編譯</p><p class="ql-align-justify">完整的算法一致切線模量推導與實現</p><p class="ql-align-justify">PDF 包含規范化的本構方程、隱式積分、徑向返回與一致切線模量推導，可供初學者學習。配套 UMAT 代碼可直接在 ABAQUS 編譯運行，采用全隱式積分搭配一致切線模量，收斂速度極快、計算精度極高，適合初學者快速入門。</p><p class="ql-align-justify">下圖展示了部分PDF內容，及umat計算結果與abaqus內置模型對比，可以發現umat收斂速度極快，與abaqus內置模型幾乎一致。
展開 
方鋼管混凝土本構模型 ￥15
方鋼管混凝土本構模型，方鋼管約束混凝土本構模型，mander混凝土本構模型，自己做的方鋼管混凝土本構模型，表格只要輸入相關參數，自動生成混凝土塑性損傷本構關系，塑性損傷本構模型。B站有鋼管混凝土軸壓驗證操作詳細視頻：https://www.bilibili.com/video/BV19R4y147gb?spm_id_from=333.337.search-card.all.click
Abaqus橡膠本構模型選擇
圖
6
試驗數據本構模型識別
圖
7
選擇可能的本構模型
其中，圖7中Test setup項可以默認；后面一個是可能相關的本構模型，可以根據數據大體判斷勾選。然后點擊OK開始根據數據進行本構模型識別。
4、在計算完成之后，會出現兩種結果，如圖8~圖10所示。
圖
8
根據數據擬合的曲線
圖
9
不同本構模型的識別結果
1
圖
10
不同本構模型的識別結果
2
由圖8可知試驗數據與不同本構模型的曲線相似度，圖9和圖10可直接判斷哪個本構模型更合適，如圖9的unstable可能不如圖10的stable本構模型合適。
然后再回到圖4中，在strain energy potential中下拉選擇比較合適的本構模型即可。
三、其他說明
需要補充幾點說明：
1、當材料行為是不可壓縮（泊松比=0.5）或非常接近于不可壓縮（泊松比＞0.475）時，則不能用常規單元來模擬（平面應力情況除外），因為此時單元中的壓應力是不確定的；
2、如圖11所示，考慮均勻靜水壓力作用下的一個單元，材料若不可壓縮，則其體積在均勻壓力作用下并不改變，單元內部的變形是非確定量，壓應力無法由單元內部積分點處的應變得到，或者無法從節點位移得到節點力；
圖
11
承受靜水壓力下的單元
3、對于具有不可壓縮材料性質的任何單元，一個純位移的數學公式是不確定的。Abaqus中采用雜交單元來處理，雜交單元包含一個可直接確定單元壓應力的附加自由度。
展開 ABAQUS umat 非線性混合硬化本構模型（Chaboche 硬化模型 ） ￥239
<p>本資源包含一份 PDF 文檔和可直接編譯運行的 Fortran UMAT 代碼，具體內容為：</p><p>Chaboche硬化本構模型 + 隱式積分 + 徑向返回</p><p>完整公式推導 + Fortran 源碼直接編譯</p><p>任意個數背應力分量 + 解析一致切線模量</p><p>PDF 包含規范化的本構方程、隱式積分、徑向返回與一致切線模量推導，可供初學者學習。配套 UMAT 代碼可直接在 ABAQUS 編譯運行，采用全隱式積分搭配一致切線模量，收斂速度極快、計算精度極高，適合初學者快速入門。</p><p>下圖展示了部分PDF內容，及umat計算結果與abaqus內置模型對比，可以發現umat收斂速度極快，與abaqus內置模型幾乎一致。
展開 橡膠的本構模型和參數
請問用abaqus分析橡膠的常用本構模型有哪些？材料參數大概是個什么范圍？
