自接触（Self-contact）是 Abaqus 中同一零件不同区域发生相互接触的仿真技术，广泛用于橡胶大变形折叠、薄壁壳体屈曲、金属冲压、织物压缩、软管挤压等工况。自接触相比普通零件间接触，更容易出现初始穿透、接触搜索失败、迭代震荡、不收敛、虚假接触、应力异常等问题。本文从自接触基础定义、法向 / 切向核心参数、接触控制、分析步设置、网格与几何预处理、Standard 与 Explicit 差异化调参思路，系统讲解自接触仿真的参数调整策略，给出工程调试流程与常见坑点，帮助 CAE 工程师高效完成自接触模型调试。
一、自接触基础选型：通用接触 vs 接触对自接触
Abaqus 提供两种自接触定义方式，建模之初就要选定，两种方式参数调整逻辑差异很大。
1. 通用接触 General Contact（推荐用于自接触）
直接勾选包含自接触，自动识别所有外表面，支持复杂折叠、多区域相互接触；适合大变形、多位置随机接触场景。可单独对自接触区域指派接触属性、接触算法、排除不需要的接触面，减少无效搜索，提升稳定性。
注意：通用接触默认会自动排除刚体自身自接触，刚体不会与自身发生接触。
2. 接触对 Contact Pair 自接触
将同一个面同时作为主面与从面；自接触强制采用平衡主从算法（balanced master-slave），适合接触区域已知、接触位置固定的模型。缺点：需要手动创建表面，复杂变形下容易漏接触，不适合接触位置不可预知的大折叠问题。
二、接触属性核心参数调整
2.1 法向接触压力 - 过盈关系（Normal Behavior）
硬接触 Hard Contact（默认）
优点：物理最准确，无额外穿透；缺点：大变形自接触极易迭代震荡、收敛困难。 调参策略：
l 标准隐式 Standard 遇到不收敛，优先改为软化接触 Softened contact（分段线性压力过盈曲线），允许微小穿透，平滑接触刚度突变；
l 或者切换罚函数 Penalty，设置合理接触刚度，刚度不宜过大（容易震荡）、不宜过小（穿透超标）。罚刚度一般取相邻单元刚度 0.1~1 倍。
增强拉格朗日 Augmented Lagrange：兼顾硬接触精度与罚函数稳定性，是隐式自接触工程首选方案。
接触粘性阻尼 Contact Damping
自接触发生反复拍打、接触分离交替（如橡胶折叠回弹），极易出现数值振荡。开启接触粘性阻尼，阻尼系数从较小值开始试算，推荐 0.01~0.1，用于抑制接触瞬间的节点震荡，提升收敛。 ⚠️ 阻尼不要过大，否则会引入虚假能量，改变仿真结果。
2.2 切向摩擦行为（Tangential Behavior）
摩擦是自接触不收敛的高频诱因，摩擦参数按下面顺序逐步调试：
1. 先无摩擦试算：模型收敛性测试阶段，先设置摩擦系数 = 0，排除摩擦带来的切向非线性，确认模型可以跑完；
2. 收敛稳定后，再逐步增加摩擦系数；
3. 弹性滑移 Elastic slip：摩擦较大的自接触模型，适当放大弹性滑移值，缓解切向刚度突变；
4. 摩擦模型选择：默认库伦摩擦，若出现粘滑交替震荡，可尝试指数衰减摩擦模型。
三、接触控制 Contact Controls（自接触收敛最重要的调参项）
3.1 自动稳定 Automatic Stabilization（Standard 隐式核心参数）
大变形自接触，接触建立前存在刚体位移、接触状态突变，打开自动稳定，引入局部粘性稳定阻尼。
l 稳定系数 Stabilization factor 推荐初始取 0.001~0.01；
l 调试原则：能收敛就尽量取小值，完成计算后关闭稳定，重新验算，确认稳定引入的附加阻尼不显著影响结果；稳定系数过大，结果失真。
3.2 接触搜索与阈值参数
1. 初始过盈调整 Initial Overclosure 导入几何存在微小面重叠，Abaqus 默认会无应变调整从节点消除小穿透；过大初始穿透不能依靠接触求解消除，必须在 CAD / 网格阶段修复。
2. 邻近搜索距离 Proximity Distance：限制接触搜索范围，减少无效面搜索，提升计算效率，同时避免远距离虚假接触；
3. 反向角 Maximum reverse angle：自接触时，两个面法向夹角超过该阈值则不判定为接触，用于排除同一连续面自身误识别（比如曲面连续相邻单元，误判自接触），常用 30°~60°。
3.3 接触离散与滑移公式
l 有限滑移 Finite sliding：自接触必须开启有限滑移，小滑移 Small sliding 只适合微小相对位移，折叠大变形模型直接报错或漏接触；
l 离散方式：面到面 Surface-to-surface 优先；节点到面 Node-to-surface 适合粗糙网格，但精度偏低；
l 通用接触中，必要时开启 Edge-to-edge 边边接触，处理壳体 / 薄膜边缘相互折叠接触。
四、分析步与求解器参数配套调整
4.1 增量步设置（Standard 隐式）
自接触属于强非线性，增量步不能过大：
1. 减小初始增量步 Initial increment size，建议 0.01~0.05；
2. 设置最小增量步下限，防止增量步无限切分；
3. 增加每个增量步最大迭代次数，默认 16，可上调至 25~30；
4. 分步加载策略：拆分多个分析步，第一步轻载荷、缓慢建立接触，第二步再加满载荷，避免一次性大变形导致接触突变。
4.2 Explicit 显式自接触调参要点
显式多用于冲击、动态折叠自接触：
1. 接触厚度：壳单元自接触务必定义接触厚度；
2. 接触搜索域：调整接触初始化搜索上下距离，保证折叠后能搜到接触面；
3. 粘性阻尼用于抑制接触冲击振荡；
4. 慎用质量缩放，质量缩放会改变惯性，影响动力学结果。
区分：Standard 适合准静态缓慢折叠；Explicit 适合动态冲击、高速挤压自接触；准静态大变形也可以用 Explicit 低速加载。
五、几何、网格预处理（很多自接触问题根源不在参数，而在网格）
1. 网格质量优先：自接触区域单元不能畸变，避免狭长单元、高度扭曲单元；接触区域网格适度加密，保证折叠接触时，至少 2~3 个单元覆盖接触区域；
2. 法线检查：自接触面单元法线方向混乱，极易造成接触识别失效，提前检查并统一单元法向；
3. 排除无效自接触面：连续光滑的同一外表面，本身不会自接触，在通用接触中做接触排除，减少大量无效接触搜索，降低计算成本与虚假接触；
4. 几何清理：CAD 模型的微小碎面、缝隙、重叠，在导入 Abaqus 前修复，微小几何缺陷会引发大量虚假接触。
Abaqus 自接触仿真调参遵循先几何网格，再接触算法，最后摩擦与稳定的原则。优先采用通用接触 + 有限滑移，隐式推荐增强拉格朗日法；调试阶段先无摩擦、小载荷、开启自动稳定，收敛成功后再逐步还原真实摩擦与边界条件。 自接触参数不是固定值，参数调整本质是平衡求解稳定性与物理精度；稳定阻尼、罚刚度这类数值参数不能一味加大，完成计算后必须评估数值引入的附加效应，避免得到数值上收敛但物理失真的仿真结果。
