在mpm88中使用Monaghan状态方程

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

Hi @Otis, 抱歉回复晚了。请问你说的算不下去是什么意思?

当那个幂取7次方的时候,计算就会卡住,然后流体信息会出现nan

hello~ 可以贴一下完整的代码不~