abaqus+franc3d裂纹扩展分析中应力强度因子求解问题？

大佬们好，本人正在用abaqus+franc3d软件做CT试样的裂纹扩展分析，现在在应力强度因子的求解方面遇到一些问题。在abaqus中进行前处理后（边界条件设置为上端UR3自由，下端U1，UR3自由，载荷施加在X轴垂直向下的拉力），导入到franc3d中插入一个3mm左右初始裂纹进行应力强度因子的求解，但是得到的结果比解析解偏大（尝试过图四的方法），如图2中所示，解析解的数值大概在840MPA√mm,而franc3d的结果在935MPA√mm左右，相差较大。求各位大佬们解答一下。