LS-DYNA中的接触界面模拟

《LS-DYNA中的接触界面模拟》由会员分享，可在线阅读，更多相关《LS-DYNA中的接触界面模拟（9页珍藏版）》请在金锄头文库上搜索。

1、LS-DYNA中的接触界面模拟1引言接触一碰撞问题属于最困难的非线性问题之一，因为在接触一碰撞问题中的响应是不平 滑的。当发生碰撞时，垂直于接触界面的速度是瞬时丕连续的。对于Coulcomb摩擦模型， 当出现粘性滑移行为时，沿界面的切向速度也是不连续的。接触一碰撞问题的这些特点给离 散方程的时间积分带来明显的困难。因此，方法和算法的适当选择对于数值分析的成功是至 关重要的。虽然通用商业程序LS-DYNA提供了大量的接触类型，可以对绝大多数接触界面进行合 理的模拟,但用户在具体的工程问题中，面临接触类型的选择及棘手的接触参数控制等问题。基于以上，本文对LS-DYNA中的接触-碰撞算法作了简要的阐

2、述，对接触类型作了详 尽的总结归纳，并对接触界面的模拟提出了一些建议。2基本概念基本概念：“slave”、“master”、“segment”。在绝大多数的接触类型中，检查slave nodes是否与master segment产生相互作用(穿透 或滑动，在Tied Contacts中slave限定在主面上滑动)。因此从节点的连接方式(或从面的 网格单元形式)一般并不太重要。非对称接触算法中主、从定义的一般原则：粗网格表面定义为主面，细网格表面为从面；主、从面相关材料刚度相差悬殊，材料刚度大的一面为主面。 平直或凹面为主面，凸面为从面。有一点值得注意的是，如有刚体包含在接触界面中，刚体的网格也

3、必须适当，不可过粗。3接触算法在LS-DYNA中有三种不同的算法处理碰撞、滑动接触界面，即：动态约束法 (kinematic constraint method)罚函数法(penalty method)分布参数法(distributed paramete method)3.1 动态约束法 Kinematic Constraint Method采用碰撞和释放条件的节点约束法由Hughes等于1976年提出，同年被Hallquit首先 应用在DYNA2ID中，后来扩展应用到DYNA3D中。其基本原理是：在每一时间步 t修正构形之前，搜索所有未与主面(master surface)接触 的从节点(s

4、lave node)，看是否在此A t内穿透了主面。如是，则缩小 t,使那些穿透主面的 从节点都不贯穿主面，而使其正好到达主面。在计算下一 t之前，对所有已经与主面接触 的从节点都施加约束条件，以保持从节点与主面接触而不贯穿。此外还应检查那些和主面接 触的从节点所属单元是否受到拉应力作用。如受到拉应力，则施加释放条件，使从节点脱离 主面。这种算法存在的主要问题是：如果主面网格划分比从面细，某些主节点（master node） 可以豪无约束地穿过从面（slave surface）（这是由于约束只施加于从节点上），形成所谓的纽 结”（Kink）现象。当接触界面上的压力很大时，无论单元采用单点还是多

5、点积分，这种现象 都很容易发生。当然，好的网格划分可能会减弱这种现象。但是对于很多问题，初始构形上 好的网格划分在迭代多次后可能会变得很糟糕，如爆炸气体在结构中的膨胀。由于节点约束算法较为复杂，目前在LS-DYNA程序中仅用于固连与固连一断开类型的 接触界面（统称固连界面，主要用来将结构网格的不协调两部分联结起来。3.2 分布参数法 Distributed Parameter Method分配参数法也是发展较早的一种接触界面算法,Wilkins在1964年将该算法成功地应用 到HEMP程序中，Burton等在1982年将其应用于TENSOR分析程序中。与节点约束法相 比，这种算法具有较好的网格

6、稳定性，因此被DYNA采用。目前，在LS-DYNA程序中用 来处理接触一滑动界面的问题。该方法的基本原理是：将每一个正在接触的从单元（slave element）的一半质量分配到被 接触的主面面积上，同时根据每个正在接触的从单元的内应力确定作用在接受质量分配的主 面面积上的分布压力。在完成质量和压力的分配后，修正主面的加速度。然后对从节点的加 速度和速度施加约束，以保证从节点在主面上滑动，不允许从节点穿透主表面，从而避免了 反弹现象。这种算法主要用来处理接触界面具有相对滑移而不可分开的问题。因此，在结构计算中， 该算法并没有太多的用处。它最典型的应用是处理爆炸等问题，炸药爆炸产生的气体与被接

7、触的结构之间只有相对滑动而没有分离。3.3 罚函数法 Penalty Method该算法于1981年有Huag等人，1982年8月开始用于DYNA2D中。现在，罚函数法 已发展为一种非常实用的接触界面算法，在数值计算中被广泛应用。罚函数法的基本原理是：在每一个时间步首先检查各从节点是否穿透主面，如没有穿透 不作任何处理。如果穿透，则在该从节点与被穿透主面间引入一个较大的界面接触力，其大 小与穿透深度、主面的刚度成正比。这在物理上相当于在两者之间放置一法向弹簧，以限制 从节点对主面的穿透。接触力称为罚函数值。“对称罚函数法”则是同时对每个主节点也作 类似上述处理。对称罚函数法由于具有对称性、动量

8、守恒准确，不需要碰撞和释放条件，因此很少引起 Hourglass （沙漏）效应，噪声小。对称罚函数法在每一个时间步对从节点和主节点循环处理一遍，算法相同。下面以从节 点n为例详细描述该算法的基本步骤：$1搜索所有从动点，确定从动点是否穿透主面。2如否，不作处理，搜索结束；如是，则在从节点与主面上的接触点间附加 一法向接触力F。n3处理摩擦力。4将接触力F和摩擦力投影到总体坐标，组集到总体载荷向量中。n4接触类型在具体介绍各种类型的接触前，先阐述几个基本的概念。在壳单元中，自动接触通过法向投影中面的1/2 “Contact Thickness”来确定接触面。这 就是“shell thicknes

9、s offsets”。接触厚度可以在接触的定义中明确指定。如果接触厚度没有 指定，则等于壳的厚度（在单面接触中，为壳厚度或单元边长的最小值）。相同的，在梁的 接触中，接触面从梁的基线偏置梁截面等效半径距离。因此，在有限元几何建模时，为考虑 壳厚、梁截面尺寸必须在壳、梁的part间有适当的间隙，否则会有初始穿透现象发生（即发 生不真实的接触现象）。虽然LS-DYNA可以通过移动穿透的从节点到主面上来消除初始穿 透，但是并不是所有的初始穿透都能检查出。DYNA中大多数的接触有一个“极限穿透深度”，如侵彻超过这个深度则从节点被释放， 接触力置为0。这主要用在自动接触中，防止过大接触力的产生而引起数值

10、不稳定性。然而 在有些情况下，因为这个阈值过早达到而使接触失效（常发生在非常薄的壳单元中）。此时 应采取的措施是放大接触厚度因子或设置接触厚度为大于壳厚度的一个值，或者改变接触刚 度的计算方法（如改为Soft=1）。LS-DYNA中的接触允许从节点与主段间压缩载荷的传递。如接触摩擦激活，也允许切 向载荷的传递。Coulomb摩擦列式用来处理从静到动摩擦的转换，这种转换要求一个衰减系 数、静摩擦系数大于动摩擦系数。关于接触搜索方法，这里仅给出几个简单的要点，详细描述见Theoretical Manual of LS-DYNAo DYNA 中有两种搜索方法：Incremental Search T

11、echnique与 Bucket SortoIncremental Search AlgorithmsGlobal Bucket Sort搜索方向仅在主段正方向从节点的穿透主面正、负方向检查穿透搜索步骤 对每一个从节点的：找出最接近的主节点；搜索接近的主段（不止一个）；搜索相邻的主段；局部利用Incremental Search确穿透检查；定最接近的主段；施加作用力。穿透检查；施加作用力。主面要求主面连续主面可以不连续特点简单、速度快非常有效，但耗时大所有的非自动LS-DYNA中的接触类型大体上可以分为五大类：One-Way Contact （单向接触）Two-Way Contact （双向接

12、触）Single Contact （单面接触）EntityTied Contac （固连接触）在以上接触类型中，前四种接触类型的接触算法均采用罚函数法。固一连接触有的采用 的罚函数法，有的采用动约束法，少部分采用分布参数法。4.1 单向接触 One-Way Treatment of ContactOne-Way、Two-Way是对接触搜索来讲的。One-way仅检查从节点是否穿透主面，而不 检查主节点。在Two-Way Contact中从节点与主节点是对称的，从节点与主节点都被检查是 否穿透相应的主面或从面。LS-DYNA中的Node_To_Surface接触类型都属于单向接触，另外还有特别注

13、明为单向 接触的_Surface_To_Surface接触类型： *Contact_Nodes_To_surface *Contact_Automatic_Nodes_To_Surface *Contact_Froming_Nodes_To_Surface (自动接触类型、主要用于金属拉压成形) *Contact_Constraint_Nodes_To_Surface (现已很少用) *Contact_Eroding_Nodes_To_Surface *Contact_One_Way_Surface_To_Surface Contact_One_Way_Automatic_Surface_To

14、_Surface由于在单向接触中，仅有从节点被检查是否穿透主面，而不考虑主节点，因此在使用时 必须注意，应保证在接触过程中主节点不会穿过从面。同样的原因，单向接触要比双向接触 运行速度快得多，因此仍被广泛应用。在以下情况中使用单向接触是合适的：主面是刚体 相对细的网格(从)与相对平滑、粗的网格(主)接触 beam_to_surface、Shell edge_to_surface 接触。beam node、Shell edge node 作从点。 在接触分析中，由于问题的复杂性，判断接触发生的方向有时是很困难的，因此分析中应尽量使用自动接触(不需要人工干预接触方向)。但当面的方向在整个分析过程中

15、都能确 定的情况下，下面的非自动接触类型是非常有效的：*Contact_Nodes_To_Surface(5)*Contact_One_Way_Surface_To_Surface(10) *Contact_Constraint_Nodes_To_Surface(18)*Contact_Eroding_Nodes_To_Surface(16)4.2 双向接触 Two-Way Treatment of Contact主、从面的定义与算法处理上是完全对称的。因此主面、从面可以随意定义。计算资源 大约是单向的2倍。LS-DYNA中绝大多数_Surface_To_Surface接触都是双向接触类型。双向接触除对主节点的搜索外，其它方面同单向接触是完全一样的。与前述接触类型5、 18、16相对应的双向接触为：*Contact_Surface_To_Surface(3)*Contact_Constraint_Surfaces_To_Surface(17) *Contact_Eroding_Surface_To_Surface(14)在 Crash Analysis 中，*Contact_Automatic_Surface_To_Surface(a3)推荐使用。在金属的 拉压成形分析中推荐使用 *Contact_