4. 接触与冲击界面

冲击分析中，接触界面往往比材料本构更容易导致计算失败。显式动力学里，接触不是“加一个约束”那么简单：它决定了冲击波能否正确传递、靶板是否被穿透、能量是否守恒，以及计算能否在合理时间内完成。本节围绕 Abaqus/Explicit 中五类核心问题展开：接触算法选择、摩擦模型、侵彻接触、接触稳定化与穿透诊断。

一、通用接触与面面接触的适用边界

Abaqus/Explicit 提供两种主流接触定义方式：通用接触（General Contact）和面面接触（Surface-to-Surface Contact）。通用接触基于罚函数或运动学约束，自动识别模型中所有外表面之间的接触，适合复杂装配、自接触、碎片飞散等难以预先指定接触对的场景。面面接触则需要用户显式指定主面（Master Surface）和从面（Slave Surface），计算效率更高，适合接触区域明确、接触对数量有限的侵彻问题。

选择原则可以概括为：接触对不确定、存在自接触或大量碎片时用通用接触；接触区域清晰、需要精细控制主从关系时用面面接触。在侵彻分析中，弹体与靶板之间通常采用面面接触，而靶板破碎后的碎片之间则依赖通用接触或自接触来避免相互穿透。

二、摩擦模型：从库仑摩擦到罚函数实现

Abaqus/Explicit 默认摩擦系数为零，必须通过 *Friction 显式定义。最常用的是库仑摩擦模型：剪应力 τ 与接触压力 p 成正比，τ = μ·p，其中 μ 为摩擦系数。当剪应力达到临界值后，接触面开始塑性滑移；在此之前，罚函数允许微小的弹性滑移，以保持数值稳定。

冲击问题中，摩擦系数通常取 0.05～0.3。弹体与金属靶板之间可取 μ=0.1～0.15；混凝土侵彻中可取 μ=0.2～0.3。若滑移速率极高，可定义摩擦系数随滑移速率衰减，模拟高速冲击下的界面软化。关键字写法如下：

*Surface Interaction, name=IntProp-1
*Friction
0.15,
*Surface Behavior, pressure-overclosure=HARD

其中 pressure-overclosure=HARD 表示硬接触，即接触压力在闭合瞬间建立，不允许穿透。若模型出现接触抖动或收敛困难，可改用 SOFTENED 或指数型压力-过盈关系，但会引入一定穿透量。

三、侵彻接触：主从面选择与侵蚀准则

侵彻分析的核心是弹体与靶板的接触。主从面选择遵循三条经验：刚度大、网格粗的一侧为主面；材料软、网格细的一侧为从面；接触区域应尽量落在主面上。典型设置是弹体为主面、靶板为从面。若靶板发生单元删除（侵蚀），从面节点消失后接触会失效，此时需启用侵蚀接触（Eroding Contact）或通用接触中的节点-表面接触，使剩余单元继续参与接触计算。

单元删除准则通常基于等效塑性应变或应力三轴度。例如靶板材料设置 *Damage Initiation, criterion=DUCTILE，当等效塑性应变达到 0.8 时删除单元。删除后，新暴露的表面自动成为接触面，这要求接触定义覆盖内部表面，通用接触天然支持这一行为。

四、接触稳定化与穿透诊断

显式分析中，接触力突变会引入高频振荡，导致能量异常或计算终止。Abaqus 提供接触稳定化（Contact Stabilization），通过施加阻尼力抑制接触抖动。关键字为 *Contact Stabilization，可指定阻尼系数和衰减时间。阻尼系数过大会人为耗散能量，建议从 1e-4～1e-3 量级试算，并监控 ALLKE 与 ALLVD 的比例。

穿透诊断主要观察三个输出：STATUS（接触状态）、CSTRESS（接触应力）、CFORCE（接触力）。若 STATUS 在 0 和 1 之间频繁跳变，说明接触不稳定；若 CSTRESS 出现尖峰后骤降，可能发生了穿透。此外，可在历史输出中请求 COPEN（接触间隙）和 CPRESS（接触压力），当 COPEN 为负值且绝对值持续增大时，表明从面节点已侵入主面。

五、工程配置示例与判断依据

以弹体冲击金属靶板为例，给出可参考的接触配置：

判断接触是否正常，可检查能量平衡：ALLKE + ALLIE + ALLVD + ALLFD ≈ ETOTAL，其中 ALLFD 为摩擦耗散能。若 ALLFD 异常增大，说明摩擦系数过高或滑移过大；若 ETOTAL 持续漂移，说明接触力做功未正确计入。建议在冲击初期（0～50 μs）重点观察接触力曲线，若出现非物理尖峰，应降低接触稳定化阻尼或细化接触区域网格。