无反射边界条件和负体积问题

本人在做磨料射流切割岩石的仿真过程中，对岩石四周施加了无反射边界条件。

数值模型计算过程中，总是提示岩石单元出现负体积，造成数值模型无法计算。

如果去除无反射边界条件且正常计算。后处理中查看产生负体积的岩石单元并未出现大变形。

尝试过提高岩石单元的硬化程度，修改时间步长，加密或放大网格尺寸，控制沙漏等操作均为解决上述问题。

请问各位经验丰富的网友指点一下，不胜感激。