1. 冲击分析工程导论

冲击问题的共同特征是:载荷在极短时间内完成输入,结构响应以应力波的形式在材料中传播,惯性效应和应变率效应不可忽略,且往往伴随大变形、接触碰撞甚至材料失效。这类问题如果沿用静力学或隐式动力学思路,会遇到两个根本困难:一是时间步长受稳定性条件限制,隐式方法在强非线性接触下收敛困难;二是冲击过程本身持续时间常在毫秒到微秒量级,需要捕捉高频瞬态响应。Abaqus/Explicit 正是为这一类问题设计的显式动力学求解器。

显式动力学与隐式方法的定位差异

显式方法基于中心差分法,不需要组装和求解整体刚度矩阵,也不进行迭代收敛判断。它通过“稳定时间增量”控制每一步的推进:单元最小特征长度除以材料波速,得到临界时间步长,实际步长通常取该值的 0.9 倍左右。这意味着显式方法的计算成本与最小单元尺寸强相关,而不是与模型自由度总数强相关。隐式方法则相反,它需要求解大型线性方程组并迭代满足平衡,适合静态、准静态和低频振动问题。

因此,选择求解器的判断依据可以归纳为:载荷持续时间是否短到惯性效应主导、是否涉及复杂接触和材料退化、是否需要捕捉应力波传播。如果答案是肯定的,显式动力学通常是更稳妥的选择。反之,如果问题是长期蠕变、低频模态或纯静力平衡,隐式方法更高效。

Abaqus/Explicit 求解流程

一个完整的显式冲击分析通常包含以下环节:前处理建模、求解设置、提交计算、结果后处理与验证。前处理阶段需要定义几何、材料本构、截面属性、装配、分析步、接触、边界条件和载荷。求解设置中,最关键的是时间步长控制、质量缩放策略和输出频率。提交计算后,需要监控能量平衡、稳定时间增量和接触状态。后处理阶段则关注应力波传播、塑性应变分布、接触力时程和失效单元比例。

下面这张图展示了显式冲击分析从输入到结果判读的完整信号链与数据流关系,包括模型输入、求解器内部推进、关键监控量和结果输出。

Abaqus/Explicit 冲击分析:输入—求解—监控—输出链路 模型输入 • 几何与网格(最小单元尺寸) • 材料本构(弹性/塑性/损伤) • 接触定义(通用接触/面面) • 边界条件与初始速度 • 分析步时间与输出频率 • 质量缩放与沙漏控制 显式求解推进 中心差分:u(t+Δt) 由前两步显式 稳定步长:Δt ≤ L_min / c_d 无需整体刚度矩阵求逆 接触搜索与罚函数/运动学约束 单元内力计算与沙漏抑制 逐增量推进至分析步结束 求解监控 • 能量平衡(ALLKE/ALLIE) • 稳定时间增量变化 • 接触状态与穿透量 • 沙漏能与伪应变能 • 质量缩放比例 • 单元删除与失效比例 结果输出与判读 场输出 • 应力波传播云图 • 等效塑性应变 PEEQ • 损伤变量与失效单元 • 位移/速度/加速度场 历史输出 • 接触力—时间曲线 • 参考点位移/速度 • 能量时程曲线 • 单元应力/应变时程 判读依据 • 能量比是否合理 • 沙漏能 < 5% 内能 • 接触穿透是否可接受 • 网格收敛性验证 典型工程场景 • 弹体侵彻靶板 • 汽车碰撞吸能 • 跌落冲击与包装 • 爆炸冲击波响应 迭代验证与模型修正

典型工程场景与建模要点

冲击分析在工程中覆盖面很广。弹体侵彻靶板问题需要关注材料在高应变率下的本构响应、单元侵蚀准则和接触算法;汽车碰撞吸能问题需要处理薄壁结构的屈曲、焊点失效和整车多体接触;跌落冲击问题通常关注包装材料的缓冲性能和内部组件的加速度峰值;爆炸冲击波问题则涉及流固耦合和远场压力加载。这些场景的共同点是:载荷时间短、非线性强、接触状态复杂。

以弹体侵彻靶板为例,一个合理的教学示例可以设定为:钢制弹体以 500 m/s 初速度撞击 10 mm 厚铝合金靶板,弹体直径 20 mm,靶板尺寸 200 mm × 200 mm。材料模型采用 Johnson-Cook 本构描述应变率强化和温度软化,失效准则采用等效塑性应变阈值。接触采用通用接触,摩擦系数取 0.15。分析步时间设为 200 μs,输出间隔 2 μs。这些参数属于教学示例设定,实际工程中需要根据材料试验数据标定。

求解过程中需要重点观察三个量:一是稳定时间增量是否在合理范围,如果突然下降几个数量级,通常意味着网格畸变或接触穿透;二是能量平衡,动能与内能之和应接近初始动能,沙漏能应控制在总内能的 5% 以内;三是接触力时程,侵彻过程中接触力会出现明显峰值和振荡,峰值大小与弹靶材料阻抗匹配相关。

课程任务总览与学习路径

本课程围绕显式动力学冲击分析展开,从求解器定位、建模流程、材料本构、接触设置、网格策略、能量监控到结果判读,逐步深入。学习过程中需要建立三个核心意识:第一,时间步长和网格尺寸是显式分析的成本与精度杠杆;第二,能量监控是判断结果可信度的第一道防线;第三,接触定义和材料失效参数对冲击响应的影响往往超过网格密度本身。

显式分析不是“不需要验证”的分析。相反,由于它不进行迭代收敛判断,结果是否可信更依赖能量平衡、沙漏控制、网格收敛性和接触穿透量等辅助指标。建议在每次提交计算后,先看能量曲线和稳定时间增量曲线,再看应力云图和接触力时程。
质量缩放可以显著提高稳定时间增量,但会引入附加惯性效应。如果缩放后总质量增加超过 1%,或动能与内能之比出现异常,需要重新评估缩放策略。对于冲击问题,质量缩放应谨慎使用,优先通过调整网格尺寸来控制计算成本。