SPH仿真中近似方法由经典的改为重整化的之后为什么计算出现问题了？

目前正在做SPH方法的精密切削仿真，由于切屑形状不好（很直，几乎不弯曲），于是想把SPH的计算方法从默认的经典算法改为Renormalization的计算方法，但是改动过后出现了如图中的问题：

出现这个问题之后计算时间无限长，无法计算了：

找了很久也没有找到解决方法，所以想来此求助。附上了K文件，为了便于排错，所以这里没加热分析，只有结构分析。

test.k