汽车碰撞模拟仿真：64 km/h 正面刚性壁碰撞减速度与能量吸收曲线

项目过程

本项目来自一款在研紧凑型 SUV 的整车安全开发阶段，委托方需要确认 100% 正面刚性壁碰撞工况下的前端吸能结构表现。整车整备质量 m=1500 kg，初始车速 v=64 km/h（17.78 m/s），对应初始动能 E_k,0=½mv²≈237 kJ，与图中参考线一致。计算目的有三点：①提取 B 柱下方减速度时程并与乘员伤害限值对照；②核算前端结构吸能完整性；③输出能量时程曲线，确认无异常沙漏能与接触能。

本项目采用 LS-DYNA 显式求解器。整车网格约 80 万单元，前端吸能盒、防撞梁、纵梁外板用 Belytschko-Tsay 壳单元（默认 5 积分厚度点），白车身内部用 1 积分单元并开启沙漏能控制（<5%）。材料模型前端用 54 号 *MAT_PIECEWISE_LINEAR_PLASTICITY（双线性随动强化），白车身用 24 号 *MAT_PIECEWISE_LINEAR_PLASTICITY_TITLE，含应变率项。整车 1 g 重力初速度通过 *INITIAL_VELOCITY_GENERATION 加在 64 km/h 方向，刚性壁 *RIGIDWALL_PLANAR 距车头 5 mm，时间步 0.5 μs，整车仿真 0~120 ms。

案例分析

减速度时程曲线与速度衰减曲线完全对应：t=0~20 ms 内前端压溃开始，减速度上升至 40 g（≈390 m/s²）峰值，恰好出现在 t≈25 ms 处；与此同时速度从 64 km/h 近似线性下降到 32 km/h（约 17.78 m/s → 8.9 m/s）。20 ms 之后减速度维持约 25 g 水平，到 60 ms 速度归零，120 ms 时整车相对刚性壁完全停止。40 g 峰值乘以等效压溃距离 ≈5.4 m（64 km/h 的动能 237 kJ 除以平均 44 kN 力）即得整车大致压溃行程，与 LS-DYNA 结果一致。

能量时程曲线方面，动能从初始 237 kJ 在前 40 ms 内单调下降，约 60 ms 时归零；内能（被结构吸收部分）从 0 单调上升，约 60 ms 后接近 237 kJ 水平并保持稳定，全程与动能曲线镜像。值得注意的是内能曲线在 t=20 ms 附近有微小”凹陷”——这是前防撞梁与吸能盒之间焊接接头断裂导致的局部能量释放，但波动幅度 <3%，不影响整车吸能完整性评估。

总结教训

本项目最容易忽视的细节是接触定义中的 *CONTACT_AUTOMATIC_SINGLE_SURFACE 与 *CONTACT_AUTOMATIC_NODES_TO_SURFACE 混用。客户最初建模时整车内骨架用单面接触，导致 A 柱内板与仪表板骨架出现穿透告警，最终改用 *CONTACT_AUTOMATIC_GENERAL 加厚度偏移 0.4 mm 才解决。

另一坑是惯性释放（*INERTIA_RELIEVE）与重力边界冲突。本项目用 v=64 km/h 初始速度，重力也用 1 g 全方向，加速度时程才能正确反映结构反应；如果误关重力，会出现减速度峰值偏小 ~5%。

项目经验

经验一：40 g 峰值对应的乘员伤害指标 HIC_15 在简化头部减速度估算下约为 720，低于 Euro NCAP 头部限值 1000 与 C-NCAP 五星限值 700~800，结构方案通过。若峰值 >45 g，需重点加强前纵梁吸能盒的诱导褶皱位置。

经验二：前端比吸能 SEA=E_absorbed/m_structure 是评价轻量化水平的核心指标。本项目 SEA≈32 kJ/kg，超过 Euro NCAP 五星门槛（25 kJ/kg），说明铝合金前端结构应用得当。

经验三：整车能量曲线的”平台后凸起”是判断沙漏能异常的金标准。本项目沙漏能在 60 ms 后仍 <2.5%，无异常；若该值超过 5%，必须加密 1 积分单元或更换为全积分单元。

有限元仿真