Abaqusvumat实现变形和损伤断裂

1）各向异性塑性本构（Hill48）和损伤演化方程（扩展的GTN模型）已经给出，根据方程，写清楚本构实现的过程，积分方法和增量方程，给一个word; 2) 给一个实现计算的VUMAT子程序，包含详细注释； 3）给一个一个单元运行的单向拉伸和剪切的应力应变曲线、孔洞参数演化曲线；4）关于VUMAT公式推导过程和子程序实现方法的答疑；5）本构方程是基于参考文献确定的，个别地方不确定，需要结合参考文献跟你讨论并最终确定；6）9月30日前给出这些结果。