1. 接触非线性问题导论

在Abaqus的求解体系中，接触非线性与几何非线性、材料非线性并列为三大非线性来源。很多初学者会把接触简单理解为“两个零件碰到一起”，但在有限元里，接触本质上是一类边界条件随求解过程不断变化的问题：接触面之间的间隙、法向压力、切向摩擦、接触区域的大小和位置，都只有在变形过程中才能确定。正因为这种“边界未知、状态可变”的特性，接触分析往往比同规模的线性静力分析更难收敛，也更依赖建模经验。

从力学角度看，接触问题包含三个必须同时处理的核心要素：法向不可穿透、切向摩擦传递以及接触状态的切换。法向不可穿透要求两个物体不能互相嵌入，Abaqus通过接触压力来阻止穿透；切向摩擦决定接触面之间能否传递剪力，常用库仑摩擦模型描述；接触状态切换则指同一对接触面在不同时刻可能处于“张开”“粘结”“滑移”三种状态之一。理解这三要素，是后续学习接触对定义、接触属性设置和收敛调试的基础。

接触非线性在工程中极为常见。例如：法兰螺栓连接中，螺栓预紧后法兰面之间的压力分布；齿轮啮合时齿面接触应力的传递；橡胶密封圈被压入沟槽后的回弹与密封压力；轮胎与地面的滚动接触；以及金属成形中模具与坯料的挤压。这些场景的共同点是：载荷通过接触面传递，接触区域随载荷变化，且往往伴随大变形或摩擦滑移。若用线性分析近似，会严重高估结构刚度或低估局部应力。

在Abaqus中，接触建模主要围绕“接触对”或“通用接触”展开。接触对需要指定主面与从面，并选择接触属性；通用接触则通过自动探测所有可能接触的面，适合复杂装配体。无论哪种方式，都需要定义法向行为和切向行为。法向常用“硬接触”，即接触压力为零时允许分离，压力大于零时不允许穿透；切向则常用“罚摩擦”或“拉格朗日摩擦”。这些设置直接决定求解器能否稳定迭代。

学习Abaqus接触非线性，建议按“先理解状态、再掌握定义、后调试收敛”的路径推进。第一步要能判断一个工程问题是否属于接触非线性，并识别接触面在哪里；第二步要会在Interaction模块中创建接触对或通用接触，设置法向硬接触和切向罚摩擦；第三步要能读懂MSG和DAT文件中的接触警告，例如“overclosure”“zero pivot”“too many attempts”，并据此调整接触刚度、摩擦系数或求解控制参数。软件界面方面，接触相关操作集中在Interaction模块，接触属性在Property模块中定义，而接触诊断信息则主要在Job模块的MSG文件中查看。

下面给出一个可直接参考的入门级接触设置示例，用于模拟一个橡胶块被刚性压头压缩的过程。该示例仅用于说明参数含义，实际数值需根据具体问题调整。

*Surface Interaction, name=IntProp-1
*Friction
0.3,
*Surface Behavior, pressure-overclosure=HARD
*Contact Pair, interaction=IntProp-1, type=SURFACE TO SURFACE
Surf-Rigid, Surf-Rubber

其中，*Friction下的0.3表示摩擦系数μ；pressure-overclosure=HARD表示法向硬接触，即接触压力为零时允许分离，压力大于零时不允许穿透；*Contact Pair指定主面为刚性压头表面，从面为橡胶块表面。提交计算后，应在MSG文件中搜索“CONTACT”关键字，观察接触状态变化和迭代次数。若出现大量“OPEN”与“STICK”交替，说明接触状态不稳定，可尝试减小初始时间增量或调整接触刚度。