Abaqus VUMAT对正交各向异性线弹性材料没有给出相同的结果

亲爱的iM万博manbetx平台echanica社区:

嗨！我试图用简单的例子验证Abaqus VUMAT，所以我写了一个正交各向异性材料的VUMAT，并将结果与Abaqus内置材料模型(弹性->工程常数)进行比较。材料常数在1.png和2.png中给出。所有其他设置从部分到网格是相同的。内置材料模型的应力分布为3.png, VUMAT的应力分布为4.png。这种差异是不可忽视的。VUMAT代码也以.txt文件的形式附上。

我在代码中遗漏了什么吗?我会很感激你的帮助。谢谢!如果问题解决了，我会为此写一篇教程，并将其发布在iMechanica和Github上。万博manbetx平台

注:Fortran代码是基于Rutger's repo的VUMAT:https://github.com/rutger-kok/composite_cdm_tan