在mpm88中使用Monaghan状态方程

我发现在mpm88中使用Momaghan状态方程,也就是在更新应力的时候(J[p]-1)需要改为(J[p]**7-1)就算不下去了,请问这是什么原因