弹塑性模型计算？

abaqus子程序怎么将uel和umat结合计算弹塑性模型？原先是将umat计算的材料矩阵和更新后的应力应变传输到uel中进行计算，但是对测试的单元计算出是结果不是很正确。

有没有好的参考文献呢？