UMAT 收敛问题请教

lungo999

最近在debug Yoshida-Uemori材料模型，UMAT采用的是Implicit 积分方法，采用最近点投影算法，但对一个单元的计算时发现，如果以0.002（即2000个微应变）个应变增量计算时结果收敛，如果增大应变增量步的话，在弹塑性转换时容易出现伪加载和伪卸载（即应力有抖动）。用Write语句写出结果发现，在欧拉迭代过程中,等效塑性应变增量delta(Delta P)出现负值。
请问，大家有没有遇到过相似的情况，个人认为DDSDDE的推导应该是正确的，有没有哪位可以推荐下解决方法，主要是想保证大应变增量下模型的收敛。

msyexplore

在欧拉迭代过程中,等效塑性应变增量delta(Delta P)出现负值,说明你的推导公式可能有问题，或者说你的本构太复杂，而你的初值选的不好，这样你在newton迭代的时候出现负值要退出循环，给一个新的初值

lungo999

回复 2# msyexplore

非常感谢，我觉得是本构复杂的原因，因为最近点投影方法在大应变增量下并不能无条件收敛，不知道楼上所说的新的初值该如何给定，因为abaqus提供的应变增量，而其他的变量都是上一步收敛的变量，如何给定新的初值呢？

msyexplore

回复 3# lungo999


    比如说，先假设你材料不发生硬化，也就是把屈服应力"冰住"，基于这个，推倒其他变量的增量，总的来说，就是你有很多变量，先假设几个不变，推倒其他变量的初值。
或者你是使用newton+line search方法改进你的迭代法

轩落子

学习了！