解决LS-DYNA中负体积方法

对于承受很大变形的材料，例如泡沫。由于在计算过程中，某个单元可能变得非常扭曲以至于单元的体积出现负值，负体积可能发生在材料未达到失效标准前。LS-DYNA中计算得到负体积(negativevolume)会导致计算终止，除非在*control_timestep卡里面设置ERODE选项为1，而且在*control_termination里设置DTMIN项为任何非零的值，在这种情况下，出现负体积的单元会被删掉从而使计算继续进行（大多数情况）。有时在计算前即使设置了ERODE和DTMIN，还有可能出现负体积导致计算错误终止。

解决单元负体积常用方法如下：

1. 简单的将材料应力－应变曲线在大应变时硬化。这种方法会非常有效；

2. 有时候修改初始网格来适应特定的变形场将阻止负体积的形成。此外，负体积通常只对非常严重的变形情况是个问题，而且特别是仅发生在像泡沫这样的软的材料上面；

3. 减小时间步缩放系数(timestepscale factor)。缺省的0.9可能不足以防止数值不稳定；

4. 避免用全积分的体单元(单元类型2和3)，这两种单元类型在大变形和扭曲的仿真中往往不是很稳定。全积分单元在大变形情况下的鲁棒性不如单点积分单元，因为单元的一个积分点可能出现负的Jacobian，而整个单元还维持正的体积。在计算中用全积分单元因计算出现负的Jacobian而终止会比单元积分单元来得快；

5. 用缺省的单元方程(单点积分体单元)和类型4或者5的沙漏(hourglass)控制(将会刚化响应)。对泡沫材料首先的沙漏方程是：如果低速冲击type6，系数1.0；高速冲击type 2或者3；

6. 对泡沫用四面体(tetrahedral)单