LS-DYNA中的接触界面模拟.docVIP

LS-DYNA中的接触界面模拟 1?引言 接触－碰撞问题属于最困难的非线性问题之一，因为在接触－碰撞问题中的响应是不平滑的。当发生碰撞时，垂直于接触界面的速度是瞬时不连续的。对于Coulcomb摩擦模型，当出现粘性滑移行为时，沿界面的切向速度也是不连续的。接触－碰撞问题的这些特点给离散方程的时间积分带来明显的困难。因此，方法和算法的适当选择对于数值分析的成功是至关重要的。虽然通用商业程LS-DYNA提供了大量的接触类型，可以对绝大多数接触界面进行合理的模拟，但用户在具体的工程问题中，面临接触类型的选择及棘手的接触参数控制等问题。基于以上，本文对LS-DYNA中的接触-碰撞算法作了简要的阐述，对接触类型作了详尽的总结归纳，并对接触界面的模拟提出了一些建议。 2?基本概念 基本概念：“slave”、“master”、“segment”。 在绝大多数的接触类型中，检查slave?nodes是否与master?segment产生相互作用(穿透或滑动，在Tied?Contacts?中slave限定在主面上滑动)。因此从节点的连接方式（或从面的网格单元形式）一般并不太重要。非对称接触算法中主、从定义的一般原则：粗网格表面定义为主面，细网格表面为从面；主、从面相关材料刚度相差悬殊，材料刚度大的一面为主面。平直或凹面为主面，凸面为从面。有一点值得注意的是，如有刚体包含在接触界面中，刚体的网格也必须适当，不可过粗。 3?接触算法 在LS-DYNA中有三种不同的算法处理碰撞、滑动接触界面，即：?动态约束法（kinematic?constraint?method）?罚函数法（penalty?method）?分布参数法（distributed?paramete?method） 3.1?Kinematic?Constraint?Method 采用碰撞和释放条件的节点约束法由Hughes?等于1976年提出，同年被Hallquit?首先应用在?DYNA2D中，后来扩展应用到?DYNA3D中。 其基本原理是：在每一时间步Δt修正构形之前，搜索所有未与主面(master?surface)接触的从节点(slave?node)，看是否在此Δt内穿透了主面。如是，则缩小Δt，使那些穿透主面的从节点都不贯穿主面，而使其正好到达主面。在计算下一Δt之前，对所有已经与主面接触的从节点都施加约束条件，以保持从节点与主面接触而不贯穿。此外还应检查那些和主面接触的从节点所属单元是否受到拉应力作用。如受到拉应力，则施加释放条件，使从节点脱离主面。 这种算法存在的主要问题是：如果主面网格划分比从面细，某些主节点（master?node）可以豪无约束地穿过从面(slave?surface)（这是由于约束只施加于从节点上），形成所谓的“纽结”(Kink)现象。当接触界面上的压力很大时，无论单元采用单点还是多点积分，这种现象都很容易发生。当然，好的网格划分可能会减弱这种现象。但是对于很多问题，初始构形上好的网格划分在迭代多次后可能会变得很糟糕，如爆炸气体在结构中的膨胀。由于节点约束算法较为复杂，目前在LS-DYNA程序中仅用于固连与固连一断开类型的接触界面（统称固连界面），主要用来将结构网格的不协调两部分联结起来。 3.2?Distributed?Parameter?Method 分配参数法也是发展较早的一种接触界面算法，Wilkins在1964年将该算法成功地应用到HEMP程序中，Burton等在1982年将其应用于TENSOR分析程序中。与节点约束法相比，这种算法具有较好的网格稳定性，因此被DYNA采用。目前，在LS-DYNA程序中用来处理接触一滑动界面的问题。该方法的基本原理是：将每一个正在接触的从单元(slave?element)的一半质量分配到被接触的主面面积上，同时根据每个正在接触的从单元的内应力确定作用在接受质量分配的主面面积上的分布压力。在完成质量和压力的分配后，修正主面的加速度。然后对从节点的加速度和速度施加约束，以保证从节点在主面上滑动，不允许从节点穿透主表面，从而避免了反弹现象。这种算法主要用来处理接触界面具有相对滑移而不可分开的问题。因此，在结构计算中，该算法并没有太多的用处。它最典型的应用是处理爆炸等问题，炸药爆炸产生的气体与被接触的结构之间只有相对滑动而没有分离。 3.3?Penalty?Method 该算法于1981年有Huag等人，1982年8月开始用于?DYNA2D中[2]。现在，罚函数法已发展为一种非常用的接触界面算法，在数值计算中被广泛应用。罚函数法的基本原理是：在每一个时间步首先检查各从节点是否穿透主面，如没有穿透不作任何处理。如果穿透，则在该从节点与被穿透主面间引入一个较大的界面接触力，其大小与穿透深度、主面的刚度成正比。这在物理