疲劳裂纹的模拟方法？

各位大神好，本人最近在研究疲劳裂纹扩展方面的课题，现已用abaqus中的扩展有限元建立了含初始裂纹的结构，并利用direct cyclic分析步进行低周循环模拟。但模拟出的结果并不理想，裂纹只扩展了1-2个单元，就算再增加循环次数，也不会继续扩展。而且direct cyclic分析步适用于低周循环，我模拟列车工况应为高周疲劳，这也是一个问题。所以想请问各位大神有其他的模拟方法吗，最好是能与abaqus联立，还是说只能二次开发了？