ls-dyna中JOHOSON_COOK关键字在撞击压缩生热过程中出现负体积和复杂声速怎样处理？ 200

kaolvwendu new.k

      目前我在做一个模拟火炸药撞击感度测试的仿真，其中主要内容是把火炸药薄片夹在上下击柱间并置于套筒中，对上击柱赋予3kg的质量，随后让它以不同的速度与炸药薄片进行碰撞。大致如图所示。

      建模后上下击柱留出了一点间隙，如图：

      目前的第一阶段，我不考虑炸药的化学反应，只考虑材料的压缩生热，因此我对中间的绿色药片部分使用的是JOHOSON_COOK关键字，搭配EOS_GRUNEISEN状态方程；上下击柱和套筒都使用一种钢材，也用JOHOSON_COOK关键字，搭配EOS_GRUNEISEN状态方程。

       但是在实际计算中，计算过程会出现复杂声速的报错，

 *** Error 40510 (SOL+510)

     complex sound speed in solid element # 97194 cycle 14118

 *** Error 40510 (SOL+510)

     complex sound speed in solid element # 97214 cycle 14118

   14118 t 6.6122E+01 dt 3.28E-06 write d3plot file            08/03/23 15:03:52

   14118 t 6.6122E+01 dt 3.28E-06 write d3plot file            08/03/23 15:03:52

结果文件导入prepost中后是这样的情况（其他部分都比较正常，我只贴出药片和下击柱部分）：

上面是药片，这是下击柱：

我也尝试过更换材料模型，将中间的绿色药片也换成跟摩擦板一样的钢材质，结果并没有什么变化。

而当我把中间的药片换成MAT_PIECEWISE_LINEAR_PLASTICITY材料模型后，则会出现负体积，但是看起来好了很多：

 *** Error 40509 (SOL+509)

     negative volume in solid element # 93174 cycle 13982

   13982 t 6.4856E+01 dt-1.01E-04 write d3plot file            08/03/23 15:39:30

   13982 t 6.4856E+01 dt-1.01E-04 write d3plot file            08/03/23 15:39:30

我想解决这个问题，但又不知从何着手，下面是我用的材料模型具体参数，我用的单位制是cm-g-us

上面这个是钢质摩擦板的参数。

这个是hmx炸药的参数。

这个是摩擦板的状态方程。

这个是药片的状态方程。

我现在完全不知道问题出在哪里，细化过网格后并没有改善，之前使用其他材料模型我记得有成过，但是没有保存（毕竟药片用其他材料意义不大），希望可以得到大佬帮助，指出问题所在。k文件附上，希望能得到回应，想做生热方面内容的朋友也可以下载查看。