contact的soft如何选择

LS-DYNA 程序中， 罚函数法 用于 接触力 的计算。 在罚函数法中，一个弹簧被放在两个物体之间，其间的接触力由下式给出

k = 接触界面刚度、 δ= 物体间的穿透量

理想情况下，接触过程中接触面间不应该有穿透

k=无穷大

 需要尽量避免初始单元穿透——可使用LS-PrePost提前检查

两种方法来计算弹簧刚度值：

           1,  根据接触面的材料特性和尺寸自动计算接触刚度。也就是SOFT=0

           2,  基于接触节点的质量和积分时间步长，它与接触材料常数无关， SOFT=1、2

故需要经常检验 d3hsp文件以确保求解时间步长 小于 接触稳定所需时间步长

图中基于材料属性和单元尺寸的最小时间步长为 0.56918e-6 【也就是timestep】小于 稳定接触所需的(0.185e-5)，因此模型是 OK的。

1， 自动计算接触刚度。

       接触刚度由接触面单元的大小和材料属性决定。因为这种方法依赖材料常数和单元大小，所以当接触面材料的刚度参数在同一个数量级上时非常有效。在材料不相似的情况下，接触可能会失效，因为LS-DYNA会粗略地把SLAVE和MASTER中最小的刚度认为是接触刚度，这样做也许刚度太小。这种情况经常会碰到，比如疏松的薄膜材料和铁材料间的接触。因此，在碰撞分析中，我们不建议使用SOFT=0 这个选项，除非你有丰富的经验确保不会发生问题。

2， 基于接触节点的质量和积分时间步长，它与接触材料常数无关

       很适合处理不同材料间的接触方式，这个刚度值等于节点的质量除以时间步长的平方。可以通过比例缩放因子来确保计算的稳定。通常，在金属和金属接触的情况下，SOFT=0或SOFT=1相当。当疏松薄膜材料和金属相接触时SOFT=1选项经常提供大于一到两个数量级的接触刚度。这个选项也用在材料不同的接触分析中。

       SOFT=2选项，其接触刚度的算法和SOFT=1一样，SOFT=2涉及基于面的接触算法，它起源于Belyischko和他的同事开发的Pinball接触。这种接触算法的处理方式是segment-to-segment接触而不是普通的node-to-segment接触。当两个四节点单元接触时，接触力加载到八个节点上以阻止单元穿透。这种处理能更加真实地离散接触力，有时候在处理复杂的接触时十分有效。最近SOFT=2选项还提供了MPP计算端口。SOFT=2选项中还不能处理梁单元。此外，SOFT=2只适合surface-to-surface接触和单面接触，而不适合nodes-to-surface接触。另外，如果希望定义segment-edge-to-segment接触，CardA中可选参数EDGE应该小心使用