对于承受很大变形的材料,例如泡沫。由于在计算过程中,某个单元可能变得非常扭曲以至于单元的体积出现负值,负体积可能发生在材料未达到失效标准前。LS-DYNA中计算得到负体积(negativevolume)会导致计算终止,除非在*control_timestep卡里面设置ERODE选项为1,而且在*control_termination里设置DTMIN项为任何非零的值,在这种情况下,出现负体积的单元会被删掉从而使计算继续进行(大多数情况)。有时在计算前即使设置了ERODE和DTMIN,还有可能出现负体积导致计算错误终止。
解决单元负体积常用方法如下:
-
简单的将材料应力-应变曲线在大应变时硬化。这种方法会非常有效;
-
有时候修改初始网格来适应特定的变形场将阻止负体积的形成。此外,负体积通常只对非常严重的变形情况是个问题,而且特别是仅发生在像泡沫这样的软的材料上面;
-
减小时间步缩放系数(timestepscale factor)。缺省的0.9可能不足以防止数值不稳定;
-
避免用全积分的体单元(单元类型2和3),这两种单元类型在大变形和扭曲的仿真中往往不是很稳定。全积分单元在大变形情况下的鲁棒性不如单点积分单元,因为单元的一个积分点可能出现负的Jacobian,而整个单元还维持正的体积。在计算中用全积分单元因计算出现负的Jacobian而终止会比单元积分单元来得快;
-
用缺省的单元方程(单点积分体单元)和类型4或者5的沙漏(hourglass)控制(将会刚化响应)。对泡沫材料首先的沙漏方程是:如果低速冲击type6,系数1.0;高速冲击type 2或者3;
-
对泡沫用四面体(tetrahedral)单