4. 接触建模前处理实战
接触非线性分析的收敛困难，往往不是求解器设置的问题，而是在前处理阶段就已经埋下隐患。Abaqus/Standard 使用 Newton-Raphson 迭代求解接触问题，每一次迭代都要根据当前接触状态判断从面节点是否穿透主面、是否建立法向压力与切向摩擦。如果几何存在微小缝隙、面法向混乱、网格疏密悬殊，求解器就会在"开—闭—开"之间反复振荡，表现为增量步被反复削减、迭代次数超限甚至不收敛。本节把前处理拆成几何清理、部件装配、接触面提取、网格划分四个环节，逐一说明它们如何影响接触收敛，并给出可操作的判断依据。
几何清理：消除接触判断的"噪声源"
导入的 CAD 几何常带有微小面片、重复边、短边和尖角。这些特征在接触搜索中会被当作真实几何处理，导致主面被切分成大量碎片，接触压力分布出现虚假峰值。清理的目标是让接触区域的面尽量"干净、连续、可预测"。
具体操作上，先用 Tools → Query → Geometry Diagnostics 检查无效面、短边和微小面。对于长度小于接触容差量级的边（例如小于 0.01 mm 而模型特征尺寸为 100 mm），应通过 Virtual Topology 合并或忽略。虚拟拓扑不改变几何形状，只改变网格划分和接触搜索时对边的识别方式，因此不会引入几何误差。对于圆角与倒角，若其半径远小于接触区尺度且不承担传力作用，可考虑在几何层面抑制，避免网格在圆角处过度加密。
需要特别注意的是面法向。Abaqus 中壳单元和刚体面的法向决定主从关系与穿透方向。若导入几何的面法向不一致，接触对可能被错误建立。可在 Mesh → Verify → Normals 中统一法向，确保所有接触面的外法向指向实体外部。
部件装配：定位精度决定初始接触状态
装配阶段的核心是让接触面在初始时刻处于"刚好接触"或"留有合理间隙"的状态。若初始穿透过大，Abaqus 会在第一个增量步强行消除穿透，产生巨大的接触力，导致迭代发散；若初始间隙过大，则接触在加载初期不建立，结构可能先发生刚体位移。
推荐做法是：在装配时让从面与主面之间保留一个小于网格特征尺寸的微小间隙（例如 0.001～0.01 mm），并在接触定义中启用 Adjust only to remove overclosure，让求解器在分析开始时自动消除过盈。对于过盈配合，则应使用 Interference fit 选项，分步建立过盈量，而不是在几何上直接重叠。
装配定位还应检查坐标系一致性。若部件通过旋转、平移装配，需确认接触面的局部坐标系方向与预期一致，否则摩擦方向会出错。对于轴对称模型，所有部件必须共享同一对称轴。
接触面提取：主从选择与法向核对
接触面提取不是简单地"点两个面"，而是要遵循主从面的选择规则。一般原则是：主面选刚度大、网格粗、面积大的一侧；从面选刚度小、网格细、面积小的一侧。若两侧刚度相近，则选网格较粗的一侧作主面，以减少从面节点搜索主面单元的计算量。
在 Abaqus/CAE 中，通过 Interaction → Create → Surface-to-surface contact 建立接触对。选择面时建议使用 Surface 工具而非直接点选单元面，这样在网格重划后接触面定义不会丢失。对于复杂曲面，可使用 Mesh → Create Mesh Part 生成独立的接触面网格，便于检查面法向。
法向核对是容易被忽略的一步。在 View → Assembly Display Options → Surface 中打开面法向显示，确认主从面的法向相对（即主面外法向指向从面，从面外法向指向主面）。若法向同向，接触对会失效或产生错误穿透。
网格划分：接触区的局部控制策略
接触收敛对网格的敏感度远高于线性分析。核心原则是：接触区网格足够细以捕捉压力梯度，但主从面网格尺寸不能悬殊。经验上，主从面单元特征尺寸之比控制在 2:1 以内，超过 3:1 时从面节点可能同时与多个主面单元建立接触，导致接触力分配振荡。
具体操作：在接触区使用 Seed → Local Seeds 设置局部尺寸，例如全局尺寸 2 mm，接触区加密至 0.5 mm。对于从面，可进一步使用 Mesh → Element Type 选择 C3D8I 或 C3D10M 等适合接触的单元类型。C3D8I 为非协调模式单元，抗剪切自锁能力较好；C3D10M 为修正二次四面体单元，适合复杂几何但计算量较大。
网格质量同样关键。接触区的单元雅可比应大于 0.6，翘曲角小于 15°。若单元过度扭曲，接触搜索会因面片法向计算误差而失败。可在 Mesh → Verify 中检查 Analysis Checks，确保没有严重警告。
把四个环节串起来看：几何清理消除搜索噪声，装配定位给出合理初始状态，接触面提取保证主从关系正确，网格划分控制搜索精度与计算量。四者中任何一环出问题，都会在求解器的迭代日志中以"增量步削减、迭代次数超限、接触状态振荡"的形式暴露出来。前处理阶段多花十分钟检查，往往能省去求解阶段数小时的调试。