本期內(nèi)容將詳細(xì)講解Abaqus二次開發(fā)之非線性強(qiáng)化Umat代碼，通過這次的學(xué)習(xí)，只要是應(yīng)用Mises屈服法則，率無關(guān)的任何非線性強(qiáng)化均可適用。

模型介紹：

彈性模量20GPa,泊松比0.3，延續(xù)上一節(jié)Umat的模型，硬化曲線方程：

材料設(shè)置：

編寫流程：

計算彈性雅克比矩陣、更新彈性試應(yīng)力、讀取狀態(tài)變量（彈性應(yīng)變和塑性應(yīng)變以及等效塑性應(yīng)變）、編寫計算強(qiáng)化系數(shù)和等效Mises應(yīng)力的子程序（方便隨時調(diào)用）、屈服后計算等效塑性增量、計算當(dāng)前應(yīng)力、彈性應(yīng)變、塑性應(yīng)變、調(diào)用Mises子程序、計算彈塑性雅克比矩陣、更新狀態(tài)變量。

本次代碼中用到了13個狀態(tài)變量，1~6個為彈性應(yīng)變，7~12個為塑性應(yīng)變，13為等效塑性應(yīng)變。

代碼講解：

部分代碼如下（篇幅有限）：

Props(1)代表?xiàng)钍夏Ａ浚?/p>

Props(2)代表泊松比，

Props(3)代表初始屈服應(yīng)力，

Props(4)~Props(6)非線性強(qiáng)化相關(guān)系數(shù)。

38~69行：計算彈性雅克比矩陣，更新彈性試應(yīng)力；

73~77行：讀取狀態(tài)變量；

83行：調(diào)用計算非線性強(qiáng)化系數(shù)子程序；

88~89行：調(diào)用等效Mises應(yīng)力子程序；

97~103行：

107~119行：（其中112~119行的循環(huán)結(jié)構(gòu)要看懂）

123~133行：

138~145行：調(diào)用等效Mises應(yīng)力，更新彈塑性雅克比矩陣；

150~154行：將彈性應(yīng)變、塑性應(yīng)變、等效塑性應(yīng)變儲存在狀態(tài)變量中。

160~185行：編寫計算非線性強(qiáng)化系數(shù)子程序；

187~217行：編寫等效Mises應(yīng)力子程序。

本次的代碼原本要為大家驗(yàn)證一下的，打開Abaqus的時候才發(fā)現(xiàn)子程序關(guān)聯(lián)失常，就把代碼先寫出來吧，有興趣大家可以自行驗(yàn)證，然后一起交流。

注：篇幅有限，故在本文中僅展示部分代碼，源代碼請在我的公眾號：“易木木響叮當(dāng)”內(nèi)回復(fù)：“非線性強(qiáng)化”，自動獲取。