利用晶体塑性有限元模拟疲劳裂纹扩展

大家好，我在abaqus中利用umat与udmgini子程序，结合xfem来模拟裂纹的疲劳扩展，模型可正常计算无报错，但计算结果裂纹始终不会扩展，我有尝试调整分析步，边界条件，晶粒的晶体取向以及udmgini子程序的损伤起始准则，但结果始终不尽如人意。我不知该如何着手调整我的模型，希望能得到大家的帮助。