1. 接触调试工程总览
在LS-DYNA的显式动力学分析中，接触定义往往是"看起来最简单、调起来最麻烦"的部分。一个整车正碰模型可能有上百个接触对，一个电池包挤压模型可能同时存在壳-壳、壳-实体、自接触与绑定接触。接触参数写错一个数量级，计算结果可能从"合理变形"变成"穿透飞散"，而求解器未必报错。因此，接触调试不是建模的附属步骤，而是一项需要独立方法论的工程活动。
1.1 LS-DYNA接触算法家族
LS-DYNA的接触算法可以按"搜索方式"和"约束方式"两个维度理解。按搜索方式，主要分为罚函数法（Penalty）、拉格朗日乘子法（Lagrange Multiplier）和罚函数+拉格朗日混合法；按接触面离散方式，分为节点-表面（NTS）、表面-表面（STS）、单面（Single Surface）和自动通用接触（Automatic General）。工程中最常用的是 *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE 与 *CONTACT_AUTOMATIC_SINGLE_SURFACE。
罚函数法的核心思想是：当从节点穿透主面时，在接触界面施加一个与穿透量成正比的法向力 F = k·δ，其中 k 为接触刚度，δ 为穿透深度。刚度越大，穿透越小，但过大刚度会引起数值振荡甚至发散。LS-DYNA默认通过 SOFT、SLSFAC、PENMAX 等参数控制这一平衡。理解这一点，是后续所有调试动作的基础。
1.2 复杂接触问题的典型特征
所谓"复杂接触"，通常不是指接触对数量多，而是指以下一种或多种情况同时出现：
- 初始穿透与初始间隙并存：CAD装配公差导致部分区域网格重叠，部分区域又存在微小间隙。LS-DYNA默认会尝试消除初始穿透，但若穿透量超过单元厚度的30%，强行消除可能引入巨大初始应力。
- 接触刚度与材料刚度不匹配：软材料（如橡胶、泡沫）与硬材料（如钢）接触时，统一刚度会导致软材料被"弹飞"或硬材料穿透。
- 自接触与多体接触耦合：薄壁结构在屈曲后发生自接触，同时与外部刚体碰撞，接触状态在时间步内频繁切换。
- 摩擦与热-力耦合：高速滑动产生摩擦生热，温度改变材料屈服强度，进而改变接触压力分布。
这些特征在求解器中的表现往往是"不报错但结果不对"。例如，接触力曲线出现高频毛刺，或总能量中滑移能（sliding energy）持续增长超过内能的5%，都提示接触定义需要调整。
1.3 调试工作的核心目标与交付物
接触调试的核心目标可以概括为三条：穿透可控、能量可信、结果可复现。穿透可控指最大穿透量不超过单元特征尺寸的10%~20%（工程经验值，需按项目要求调整）；能量可信指沙漏能低于总能量的5%、滑移能无异常突增；结果可复现指同一K文件在不同计算平台或不同LS-DYNA版本上给出一致趋势。
为达成这些目标，调试工作应产出以下交付物：
| 交付物 | 内容要点 | 判断依据 |
|---|---|---|
| 接触参数表 | 接触类型、SOFT、SLSFAC、FS、FD、PENMAX、DEPTH | 与材料、网格尺寸匹配 |
| 调试日志 | 每次修改的参数、修改原因、观察到的现象 | 可追溯、可回退 |
| 能量曲线 | 内能、动能、沙漏能、滑移能随时间变化 | 沙漏能<5%，滑移能无突增 |
| 穿透云图 | 最大穿透量、穿透位置、穿透发生时刻 | 穿透量<单元尺寸10%~20% |
| 可复现K文件 | 包含所有接触卡、控制卡、输出卡 | 独立运行结果一致 |
*CONTROL_CONTACT 的 NSBCS、INTERM 等参数会影响接触搜索频率和穿透消除策略。修改这些参数前，务必记录原始值，并在小模型上验证效果。
本章后续内容将围绕上述目标展开：先建立接触算法的分类认知，再通过典型故障现象反推参数调整方向，最终形成一套可落地的调试流程。读者应带着自己的模型问题阅读，边学边改，而不是等全部学完再动手。