Niordson和Hutchinson(2003)的应变梯度塑性框架在Abaqus UEL子程序中的实现问题

我在UEL中实现了Niordson和Hutchinson(2003)的速率无关应变梯度塑性框架。我能够验证两个刚性表面之间无限宽薄层剪切的结果(Niordson和Hutchinson, 2003)。然而，对于更复杂的有限宽板剪切问题，UEL子程序不能收敛。我怀疑收敛问题是由内部弹塑性边界条件和屈服标准引起的。在这方面我没有什么疑问。

1.Niordson和Hutchinson(2003)提到，每当在前一个增量的高斯点处获得负塑性增量时，高斯点被弹性卸载。在这种情况下应力更新方法是什么?会是dQ=0 d =0 d = L*d(因为d ^p =0)同样，在弹性区域q=\sigma_e。因此，在弹性卸载(负塑性应变增量)期间，还应该赋值Q=\sigma_e。

2.在这样一个步骤中，如果在之前的增量中发现了负的塑性增量，那么对于该高斯点，刚度矩阵K_p也应该设置为一个非常大的数字，例如K_p = 10^ 8e。

3.文中还提到，当Q达到Qy时，塑性屈服再次开始。这是否意味着一旦对高斯点进行弹性卸载，就可以根据上一项(2)的条件求解系统=\s Q达到Qy。在这段时间里Q=\sigma_e一直都是吗?同样，Qy是否对应于开始弹性卸载时的增量Q值?

如果有人能说明这个问题，那就太有帮助了。

多谢。