GTN模型在abaqus的二次开发

  现在比较疑惑的地方在于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

一起学习，共同进步。。。