1、非线性有限元第 10章 接触 -碰撞计算固体力学第 10章 接触 -碰撞1 引言2 接触界面方程 3 摩擦模型 4 弱形式 5 有限元离散 6 显式方法 1 引言在制造加工过程的仿真中包括接触和碰撞。例如,在薄金属板的成型中,模具和工件之间接触面的模拟;汽车碰撞和挤压成型的模拟,都需要有接触滑移界面。接触问题的失效形式是相互作用的接触面发生破坏。对于接触 -碰撞,将展示 Lagrangian网格的控制方程和有限元程序。在接触中,物体的控制方程与前面介绍的方程是一致的,但在接触界面上,需要增加动力学和运动学的条件。关键条件是 不可侵彻性条件 :即两个物体不能互相侵入的条件,不可侵彻性不能表示为一
2、个简单的方程,所以,发展了几种简单的方法。其中两个基于最近点映射的形式:1)显式动态问题适用的率形式;2)主要适用于隐式方法的平衡解答。此外,经典的 Coulomb摩擦模型和界面本构模型。1 引言给出控制方程的弱形式,处理接触界面约束的 4种方法:1 Lagrange乘子法,2 罚方法,3 增广的 Lagrangian法,4 摄动的 Lagrangian法。由 Lagrange乘子法,在接触问题的离散化中,在接触界面上乘子必须是近似的。 乘子必须满足法向面力是压力的约束 。在罚方法中,面力不等式源于 Heaviside分步函数,该函数被嵌入在罚力之中。1 引言接触 -碰撞问题是属于最困难的非线
3、性问题之一,因为在接触 -碰撞问题中的响应是不平滑的。当发生瞬时接触时,垂直于接触界面的速度是瞬时不连续的。对于 Coulomb摩擦模型,当出现粘性滑移行为时,沿着界面的切向速度是不连续的。接触 -碰撞问题的这些特性给离散方程的时间积分带来了明显的困难,削弱了 Newton算法的功能。因此,选择适当的方法和算法是至关重要的,并且在获得强健的求解程序中,规则化的技术是非常有用的。2 接触界面方程模拟接触 -碰撞问题的标记 接触界面包含两个物体表面的交界。 接触界面包括两个物体处于接触的两个物理表面,它们是重合的,在数值计算中,两个表面一般不重合,分为主控和从属表面。接触界面是时间的函数,确定它是
4、接触 -碰撞问题解答的重要部分。在通用有限元软件中,接触算法能够处理多个物体的相互作用,然而多个物体的接触包含成对物体的相互作用。因此,从考虑两个物体的问题入手。2 接触界面方程在主控接触表面的每一点建立局部坐标系统,可以构造相切于主控物体表面的单位矢量:物体 A的法线给出为 在接触界面上 有即两个物体的法线方向相反。以局部分量的形式表示速度场在三维问题中希腊字母下角标的取值范围为 2;当问题是二维时,接触表面成为一条线,取值为 1。 2 接触界面方程物体由标准场方程控制:质量、动量和能量的守恒,应变度量,以及本构方程。接触增加了条件:1 在界面上,两个物体 不可相互侵入 和 面力满足动量守恒
5、 ;2 横跨接触界面的 法向面力不能为拉力 。按照要求分类:1 对于 位移和速度 的要求作为 运动学 条件,2 对于 面力 的要求作为 动力学 条件。 不可侵彻性条件2 接触界面方程一对物体的不可侵彻性条件可以表示为交集为零 例如,如果物体在旋转中,对于P 点接触 Q 点是可能的,而一个不同的相对运动可能导致 P 点与S 点接触。结论是,除了以一般的形式,找不到其它的方程表示P 点没有侵入物体 A 的事实。 两个物体不允许重叠,这可以视为一个协调条件。对于大位移问题,不可侵彻性条件是高度非线性的,并且一般不能以位移的形式表示为一个代数方程或者微分方程。其困难源于在一个任意运动中,不可能预先估计到两个物体的哪些点将发生接触。不可侵彻性条件运动学2 接触界面方程由于以位移的形式表示交集为零的公式是不可能的,所以,在接触过程的每一阶段中以率形式或者增量形式表示不可侵彻性方程是很方便的。其率形式应用到物体 A和 B上发生接触的部分,即是位于接触表面上的那些点 两个物体的相互侵彻速率 利用 点乘 得到上两式 并且利用法线是正交于与平面相切的单位矢量的事实。