耦合损伤的弹塑性本构umat子程序代码计算损伤累积和疲劳寿命是否正确？

最近在复现一篇论文，涉及到圆柱滚子轴承的损伤累积和疲劳寿命计算，施加Pmax=2500MPa，滚子：R=9，但是算出来的寿命和论文中的小10倍，论文中是1.2e6，我算出来3.7e5。不知道哪里有问题

图1.损伤累积曲线（横坐标寿命，纵坐标损伤值）

图2.mises应力图（上为滚子：弹性，下为滚道：弹塑性用umat子程序算损伤）

图3.寿命云图（不知道为什么寿命只有一个值，不同单元的寿命是相同的，这里是不是也有问题？）

图4.接触压力CPRESS云图

图5.CPRESS的曲线图是这样，横坐标是x的坐标

参考文献是上述链接