有限元分析丨接触（二）
接触设置特别重要，这篇文章没任何原创，只能学习过程中记录。 目录2 接触算法2.1 法向接触算法2.1.1 罚函数算法和增广拉格朗日算法（惩罚方法）2.1.2 法向拉格朗日算法2.1.3 MPC (多点约束)算法2.2 切向接触算法2.3 渗透量 3 接触设置3.1 接触行为3.1.1不对称3.1.2 对称3.2 裁剪设置2 接触算法为了更全面的认识接触，接下里从接触算法角度认识接触。前文提到接触是数值接触而非物理上认知中的几何接触。 接触算法我自己用的乱七八糟...努力试着争取写明白。 首先看一下Workbench帮助文档关于接触算法内容，可按照法向接触和切向接触进行分类说明。2.1 法向接触算法法向接触算法包括：① 罚函数算法（Pure Penalty）② 增广拉格朗日算法（Augmented Lagrange）③ 法向拉格朗日算法（Normal Lagrange）④ MPC算法。注意：其中①、②是基于惩罚方法，存在穿透，惩罚函数是一种前置接触的数学方法，而不是限制它。 2.1.1 罚函数算法和增广拉格朗日算法（惩罚方法）对于存在接触关系的几何模型是不能互相穿透的，但罚函数算法和增广拉格朗日接触算法需要接触面和目标面之间存在小穿透量Xp，接触刚度为Kn。则：纯罚法和增广拉格朗日法的主要区别在于增广拉格朗日法增加了接触力的计算。由于附加了λ项，增广拉格朗日法对接触刚度Kn的大小不敏感。注意：① 对于“软的”（如橡胶）和硬材料接触存在变形时的，适合选择罚函数和法向拉格朗日算法；② 接触界面弹簧刚度越高，穿透量就越小；③ 当Kn无限大时，则穿透量Xp接近0，这对罚函数是无法实现的；④ 但如果Xp足够小或可忽略，则是可以办到的，这时就认为求解的结果是精确的。2.1.2 法向拉格朗日算法使用法向拉格朗日中接触压力被作为额外的自由度参与求解，以强制满足接触界面变形协调性，不涉及接触刚度和穿透问题。注意：这种算法缺点是存在震荡问题，收敛困难，计算耗时长（所有接触算法中耗时最长的）。而使用罚函数算法相对来说更容易收敛，因为接触状态的改变不再存在阶跃。2.1.3 MPC (多点约束)算法对于绑定和不分离接触方式，可以使用MPC接触方式，MPC内部添加了约束方程来“连接”接触面之间的位移。这种方法不是基于惩罚或拉格朗日算法。是一种直接、有效的连接接触区域表面的方法。基于MPC算法的绑定接触支持大变形效应。注意：使用罚函数和增广拉格朗日算法使用积分点进行接触探测，而拉格朗日算法和MPC算法使用节点进行探测，节点探测数量少；可以通过局部网格细化可以使节点探测获和积分点接触探测相同效果。2.2 切向接触算法对于无摩擦、粗糙和摩擦这3种接触类型，接触发生在切向方向上。同法向不 穿透条件类似，如果在切向上两个实体是粘接在一起（不存在相对滑动），则在切向方向一般采用罚函数来防止相对滑动，切向引入刚度Kt的弹簧，切向力和滑动距离之间满足：罚函数算法增广拉格朗日算法法向拉格朗日算法MPC算法收敛性很好如果穿透量较大，迭代次数较多如果存在抖振，迭代次数较多。很好接触刚度影响大影响小无穿透不可控制一定程度控制≈00接触类型所有接触类型只适用于绑定与不分离求解器Iterative或Direct求解器只能用Direct求解器Iterative或Direct求解器对称与非对称对称或非对称非对称接触检测高斯积分点接触节点2.3 渗透量 Penetration 罚函数算法和增广拉格朗日算法是允许渗透的，对于渗透进行一些简单说明。渗透量如何获得？以及如何使用Contact Tool使用后续内容以案例形式进行说明。（给自己挖坑了）3 接触设置物体接触之前，接触面和目标面是互相分离的，至少有一个物体处于无约束的状态，在受到外界载荷作用时，必然出现刚体运动。对于线性接触类型，可以通过设置Pinball Region，保证接触面和目标面之间的初始间隙忽略;对于非线性接触，接触面和目标面之间的初始间隙可以通过以下方法来消除。3.1 接触行为接触行为包括：① 系统自动控制（Program Controlled）② 不对称（Asymmetric）③ 对称（Symmetric）④ 自动非对称（Auto Asymmetric）。3.1.1不对称（Asymmetric）即一侧为接触面，一侧为目标面，系统默认设置；3.1.2 对称（Symmetric）即两侧互为接触面以及目标面；注意：① 只有罚函数算法和增广拉格朗日算法支持对称行为；法向拉格朗日算法和MPC算法支持非对称行为；② 当接触面和目标面不易区分情况下，推荐使用对称行为； ③ 网格质量差，十分粗糙时，推荐使用对称行为；④ 使用非对称行为时，接触结果仅使用于Contact Side当使用对称行为时，实际接触结果是两个接触面Contact Side的平均值；⑤ 对称行为计算成本大，代价大。3.2 裁剪设置Trim contact用于自动减少接触对数量，提高求解速度。可供选择包括：On，一般情况选择为On；Off：当存在手工创建的接触区域时，设置为Off。注意：有大挠度滑动必须设置为Off，否则在原始接触边缘外出现穿透现象。如果设置为On，一定要定义合理的Trim Tolerance，以保证足够的接触区域。Trim Tolerance为裁剪容差，但Trim选项设置为program controlled，仅对自动探测接触可用。当Trim选项设置为On时，对自动探测接触和手工定义接触均可用。来源：认真的假装VS假装的认真
