LS-DYNA的接触算法简介

（2）DistributedParameter Method     

分配参数法也是发展较早的一种接触界面算法,与节点约束法相比，这种算法具有较好的网格稳定性，因此被DYNA 采用。目前，在LS-DYNA程序中用来处理接触一滑动界面的问题。

该方法的基本原理是：将每一个正在接触的从单元(slave element)的一半质量分配到被接触的主面面积上，同时根据每个正在接触的从单元的内应力确定作用在接受质量分配的主面面积上的分布压力。在完成质量和压力的分配后，修正主面的加速度。然后对从节点的加速度和速度施加约束，以保证从节点在主面上滑动，不允许从节点穿透主表面，从而避免了反弹现象。

这种算法主要用来处理接触界面具有相对滑移而不可分开的问题。因此，在结构计算中，该算法并没有太多的用处。它最典型的应用是处理爆炸等问题，炸药爆炸产生的气体与被接触的结构之间只有相对滑动而没有分离。

（3）Penalty Method     

现在，罚函数法已发展为一种非常用的接触界面算法，在数值计算中被广泛应用。罚函数法的基本原理是：在每一个时间步首先检查各从节点是否穿透主面，如没有穿透不作任何处理。如果穿透，则在该从节点与被穿透主面间引入一个较大的界面接触力，其大小与穿透深度、主面的刚度成正比。这在物理上相当于在两者之间放置一法向弹簧，以限制从节点对主面的穿透。接触力称为罚函数值。“对称罚函数法”则是同时对每个主节点也作类似上述处理。

      对称罚函数法由于具有对称性、动量守恒准确，不需要碰撞和释放条件，因此很少引起Hourglass 效应，噪声小。

     对称罚函数法在每一个时间步对从节点和主节点循环处理一遍，算法相同。下面以从节点ns 为例详细描述该算法的基本步骤：