第十一章接触..doc

第十一章 接触 许多工程问题涉及两个或多个部件间的接触。在这类问题中，当两个物体彼此接触时，物体间存在沿接触面法向的相互作用力。如果接触面间存在摩擦，沿接触面的切线方向也会产剪力以抵抗物体间切向运动（滑动）。接触模拟通常的目标是确定接触面积及计算所产生的接触压力。 在有限元中，接触条件是一类特殊的不连续的约束，它允许力从模型的一部分传递到另一部分。因为只有当两个表面接触时才用到接触条件，所以这种约束是不连续的。当两个接触的面分开时，就不再存在约束作用了。因此，分析方法必须能够判断什么时候两个表面是接触的并且采用相应的接触约束。同样，分析方法也必须能判断什么时候两个表面分开并解除接触约束。 在ABAQUS接触分析过程中，必须在模型的各个部件上创建可能接触的面。一对彼此可能接触的面，称为接触对，必须被标识。最后各接触面服从的本构模型必须定义。这些接触面间的相互作用的定义包括诸如摩擦等行为。 11.1 接触面间的相互作用 接触面之间的相互作用包含两个部分：一部分是接触面的法向作用，另一部分是接触面的切向作用。切向部分包括接触面间的相对运动（滑动）和可能的摩擦剪应力。 11.1.1 接触面法向性质 两个面之间分开的距离称为间隙。当两个面之间的间隙变为零，接触约束就起作用了。在接触问题的公式中，对接触面之间相互传递接触压力的大小未作任何限制。当接触压力变为零或负值时，接触面分离，约束就被撤出。这个行为称为“硬”接触。图11.1中的接触压力－间隙关系中描述了这种行为。 图11.1 硬接触的接触压力与间隙的关系 当接触条件从“开”（正的间隙）到“闭”（间隙为零）时，接触压力发生剧烈的变化，有时使得接触计算很难完成。在后面的章节中将讨论克服接触计算困难的若干技术。另外一些信息可参见ABAQUS/Standard 用户手册22.1.1节。和ABAQUS/Standard讲义的接触部分。 11.1.2 表面的滑动 除了要确定在某一点是否发生接触外，分析中还必须计算两个表面间的相互滑动关系。这个问题可能是一个很复杂的计算；因此，在分析时ABAQUS对小滑动量和有限滑动量做了区分。接触面间小滑移量问题的计算量较小。什么是“小滑移”通常很难定义，但一般所遵循的原则是，当一点与一表面接触时，只要这点滑动量不超过一个典型的单元尺度的很小部分，就可以近似地认为是“小滑动”。 11.1.3 摩擦 当表面接触时，就像传递法向力一样，接触面间要传递切向力。所以分析时需要考虑阻止面之间相对滑动趋势的摩擦力。库仑摩擦是常用的描述接触面的相互作用的摩擦模型。这个模型用摩擦系数??来描述两个表面间的摩擦行为。乘积?p给出了接触面之间摩擦剪应力的极限值，这里p是两接触面之间的接触压力。直到接触面之间的剪应力达到摩擦剪应力的极限?p时，接触面间才发生相对滑动。大多数表面的?通常小于单位1。图11-2中的实线描述了库仑摩擦模型的行为：当它们粘结在一起即剪应力小于?p时，表面间的相对运动(滑移)量为零。 图11.2 摩擦特性 在分析过程中，在粘结和滑移两种状态间的不连续性，可能导致收敛问题。只有在摩擦力对模型的响应有显著的影响时才应该在分析中考虑摩擦。如果在有摩擦的接触分析中出现收敛问题，首先必须尝试改进的方法之一就是重新进行没有摩擦的分析。 模拟真实的摩擦行为可能是非常困难的；因此在默认情况下， ABAQUS使用一个允许“弹性滑动”的罚摩擦公式，见图11-2中的虚线。“弹性滑动”是指表面粘结在一起时所发生的小量的相对运动。ABAQUS会自动选择罚刚度（虚线的斜率），从而这个允许的“弹性滑动”的滑动值只有单元特征长度非常小的部分那么大。罚摩擦公式适用于对大多数问题，其中包括大部分金属成型问题。在那些必须包括理想的粘结－滑动摩擦行为的问题中，可以使用“Lagrange”摩擦公式。使用“Lagrange”摩擦公式需要花费更多的计算机资源，其原因是在使用“Lagrange”擦公式时ABAQUS需对每个摩擦接触的表面节点额外增加变量。另外，解的收敛会很慢，通常也需要更多的迭代。在这本指南中不讨论这种摩擦公式。 通常刚开始滑动与滑动中的摩擦系数是不同的。前者称为静摩擦系数，后者称为动摩擦系数。在 ABAQUS/Standard中用指数衰减规律来模拟静和动摩擦系数的变化。在这本指南中不讨论这个摩擦公式。 在模型中考虑了摩擦，就会在求解的方程组中增添了不对称项。如果?值0.2，不对称项的值及其影响非常小，一般而言，采用正规的、对称求解器法求解的效果还是很好（接触面的曲率很大除外）。在摩擦系数较大时，会自动要调用非对称求解器求解，因为它将改进收敛速度。非对称求解器所需