現在比較疑惑的地方在于GTN本身是個損傷模型，但它會導致屈服函數的改變，如果材料的應力應變關系用相關流動應力模型表示，如Hansel-Spittel高溫本構模型，這兩個模型會不會沖突？個人感覺應該不會，因為在abaqus設置多孔材料（porous metal plasticity）那里（就是GTN模型）設置完參數后也需要提供塑性應變和應力。那是否在二次開發（VUMAT和umat）中除了GTN模型，還需要提供相關的流動應力模型？按照這個思路我編寫了一個含有GTN和Hansel-Spittel高溫本構模型的vumat子程序，由于umat需要提供剛度矩陣看了相關論文還是不知道怎么整，編到最后放棄了。子程序大致思路按照論文《GTN模型的算法研究、程序開發及試驗驗證》編寫。最后做出來也不知道對不對，想請各位大神來教教我提供點意見，這啥都不懂一接手就高難度還沒教程太難了......

GTN VUMAT.rar

百度網盤鏈接：

鏈接：https://pan.baidu.com/s/1_AK23HgpzhkJusf7h_Escg 

提取碼：1234

一起學習，共同進步。。。