1. 冲击仿真概述与Abaqus/Explicit定位
冲击问题在工程中极为常见：汽车保险杠撞击刚性墙、手机跌落、弹体侵彻靶板、包装件运输跌落、叶片被飞鸟撞击。这类问题的共同特征是载荷在极短时间内释放、结构响应强烈非线性、并伴随接触状态突变。用有限元方法求解时，最核心的困难不是网格划分，而是时间积分方式的选择——这正是隐式与显式方法的分水岭。
显式动力学的基本思想
显式方法求解动力学方程 Mü + Cu̇ + Fint = Fext 时，把质量矩阵 M 处理为集中（对角）质量矩阵。对角矩阵求逆只需逐项取倒数，因此加速度可以直接写出：
ü(t) = M⁻¹ · [ F_ext(t) − F_int(t) − C·u̇(t) ]随后用中心差分法推进速度与位移：
u̇(t+Δt/2) = u̇(t−Δt/2) + (Δt_{t+Δt} + Δt_t)/2 · ü(t)
u(t+Δt) = u(t) + Δt_{t+Δt} · u̇(t+Δt/2)整个过程不需要组装整体刚度矩阵、不需要迭代求解线性方程组，每个时间增量只做向量运算，计算量与单元数近似成正比。代价是条件稳定性：时间增量必须小于临界值 Δtcr ≈ Le/cd，其中 Le 是单元特征长度，cd 是材料膨胀波速。以钢为例，弹性模量 210 GPa、密度 7800 kg/m³，波速约 5190 m/s；若单元尺寸 2 mm，则 Δtcr ≈ 3.9×10⁻⁷ s。冲击过程通常持续 1~10 ms，意味着需要数千到数万个增量步——这正是显式方法"步数多、每步便宜"的典型画像。
冲击问题的四个特征
- 短时强载：载荷持续时间与结构自振周期同量级，惯性效应不可忽略，准静态假设失效。
- 几何与材料非线性：大变形、大转动普遍存在；金属在高应变率下屈服强度上升，需用 Johnson-Cook 等率相关本构。
- 接触状态突变：撞击瞬间接触对从"未接触"跳到"接触"，法向力剧烈变化，隐式方法容易反复迭代甚至不收敛。
- 能量转换剧烈：动能迅速转化为内能、塑性耗散能与摩擦耗散能，能量平衡是判断结果可信度的关键指标。
隐式与显式方法对比
隐式方法（Abaqus/Standard）在每个增量步内求解 KΔu = R，需要迭代并判断收敛，时间步长可取得较大，适合长期准静态、低频响应问题。显式方法（Abaqus/Explicit）无需迭代，但步长受稳定性限制。二者并非优劣之分，而是适用场景不同：
| 对比维度 | Abaqus/Standard（隐式） | Abaqus/Explicit（显式） |
|---|---|---|
| 时间积分 | Newmark 隐式，需迭代 | 中心差分，显式推进 |
| 时间步长 | 由精度控制，可较大 | 受 Δtcr 限制，通常很小 |
| 每步成本 | 高（求解方程组） | 低（向量运算） |
| 接触处理 | 迭代收敛困难 | 自然处理状态突变 |
| 典型场景 | 静力、低频、长期响应 | 冲击、碰撞、侵彻、爆炸 |
Abaqus/Explicit 模块界面与求解流程
在 Abaqus/CAE 中，显式分析通过 Step 模块创建 Dynamic, Explicit 分析步，并设置时间长度（Time period）。与隐式分析不同，显式分析步没有增量步长输入项，软件自动根据网格与材料计算稳定增量；用户可指定质量缩放（Mass scaling）或时间缩放来调节计算成本。载荷、边界条件、接触均在对应模块中定义，接触推荐使用通用接触（General contact）以自动识别所有外表面。
求解流程可概括为：前处理建模 → 提交作业（Job）→ 显式积分推进 → 输出场变量与历史变量 → 后处理判读。其中能量历史输出（ALLKE、ALLIE、ALLPD、ALLSE、ALLAE）是显式分析特有的诊断工具：若沙漏能 ALLAE 超过内能的 5%~10%，说明网格或单元控制需要调整；若动能远大于内能，则可能加载过快或存在未约束刚体运动。
一个可复现的教学示例
为把上述概念落到具体操作，下面给出一个可执行的示例设置（参数为教学假设值，非实测结果）：
- 对象：100 mm × 100 mm × 5 mm 铝合金平板，密度 2700 kg/m³，弹性模量 70 GPa，泊松比 0.33，采用 Johnson-Cook 塑性参数。
- 网格：C3D8R 六面体单元，单元尺寸 1 mm，厚度方向 5 层，共约 50000 单元。
- 冲击条件：直径 20 mm 刚性球以 50 m/s 初速度垂直撞击板中心，板四边固支。
- 分析步：Dynamic, Explicit，Time period = 0.5 ms。
- 接触：General contact，法向硬接触，切向摩擦系数 0.2。
- 输出：场输出每 0.02 ms 一帧；历史输出记录撞击点位移、接触力与全部能量项。
提交后可在 Monitor 中观察稳定时间增量。若单元尺寸 1 mm、铝中波速约 5090 m/s，则 Δtcr ≈ 1.96×10⁻⁷ s，0.5 ms 约需 2550 个增量步。若计算时间过长，可对非关键区域施加质量缩放，例如将最小稳定增量提升到 5×10⁻⁷ s，步数降至约 1000 步，但需检查动能是否被人为放大。
理解显式方法"以步数换单步成本"的本质，以及冲击问题对接触突变与率相关本构的敏感性，是后续建模、调试与结果判读的基础。Abaqus/Explicit 正是为这类瞬态、强非线性、接触主导的问题而设计的求解器，其模块界面与求解流程也围绕"快速推进 + 能量诊断"这一核心逻辑展开。