gruneisen状态方程的奇怪现象

木风10086

做一个碰撞分析，材料用的JC模型和gruneisen状态方程，每次计算都出错，后来发现初始的pressure就已经非常大，所以一开始计算结构就不断膨胀，很快失效。可是检查很久也没发现状态方程和JC模型哪里出问题了，用的都是国际单位制。后来建立了个简单模型，两个方块碰，同样的材料参数就没有这个问题。不知道有没有哪位有过类似的经历，给指点一二。
*MAT_JOHNSON_COOK
$#     mid        ro         g         e        pr       dtf        vp    rateop
         3 7800.00008.0700E+102.1000E+11  0.300000     0.000     0.000     0.000
$#       a         b         n         c         m        tm        tr      epso
7.9000E+8 4.9700E+8  0.261000  0.014500  1.030000 1800.0000 294.00000  1.000000
$#      cp        pc     spall        it        d1        d2        d3        d4
470.00000-9.000E+11  2.000000     0.000  0.800000     0.000     0.000     0.000
$#      d5      c2/p      erod     efmin
     0.000     0.000         0 1.0000E-6
*EOS_GRUNEISEN
$#   eosid         c        s1        s2        s3     gamao         a        e0
         1 4569.0000  1.490000     0.000     0.000  2.170000  0.460000     0.000
$#      v0
     0.000

木风10086

初始的pressure不应该默认为0吗？