1. 汽车碰撞安全仿真概述

汽车碰撞安全仿真,本质上是把一次真实撞车过程搬到计算机里“重演”:用有限元方法把车身、乘员约束系统、障碍物离散成成千上万个单元,用显式时间积分求解碰撞瞬间的大变形、接触与失效,再通过假人模型或人体生物力学模型评估伤害风险。它之所以成为现代汽车开发的核心手段,是因为实车碰撞试验成本高、周期长、可重复性差,而仿真可以在概念设计阶段就快速比较几十种结构方案,把“撞坏再改”变成“算准再试”。

被动安全法规与 NCAP 评价体系

被动安全的目标是事故不可避免时,尽量降低乘员与行人的伤害。它由两类要求共同约束:一类是强制性法规,另一类是新车评价规程。法规是准入门槛,NCAP 是市场竞争力,两者工况、评价指标并不完全一致,仿真工程师必须同时满足。

常见法规工况包括正面 100% 重叠刚性壁障碰撞(如 GB 11551、ECE R94)、正面 40% 重叠可变形壁障偏置碰撞、侧面移动变形壁障碰撞(GB 20071、ECE R95)、以及行人保护(GB/T 24550)等。法规通常以假人特定部位的伤害指标为判据,例如头部 HIC36/HIC15、胸部压缩量、大腿轴向力、腹部力等,超过限值即不通过。

NCAP 则更严格、更全面。以中国 C-NCAP 为例,包含正面 100% 重叠刚性壁障、正面 50% 重叠移动渐进变形壁障(MPDB)、侧面柱碰撞、侧面移动变形壁障、鞭打试验以及行人保护、主动安全等板块,用星级和得分率呈现。Euro NCAP 与 C-NCAP 在壁障质量、速度、假人类型(如 THOR、WorldSID)上存在差异,同一车型往往需要分别建模、分别校核。下表给出典型工况的对比,便于建立整体印象。

工况典型速度主要评价对象关键伤害指标
正面 100% 刚性壁障50 km/h前排乘员HIC、胸部加速度、大腿力
正面 40% 偏置可变形壁障64 km/h前排乘员、车体侵入HIC、胸压、小腿力
侧面移动变形壁障50 km/h侧碰侧乘员肋骨变形、腹部力、骨盆力
侧面柱碰撞32 km/h侧碰侧乘员头部、胸部、骨盆
行人保护20–40 km/h行人头部、腿部HIC、小腿弯矩、膝部剪切

需要强调的是,法规与 NCAP 的限值、假人标定、壁障刚度都是公开但不断更新的,仿真中必须使用对应版本的假人模型与壁障模型,否则结果无法用于认证或评分。

LS-DYNA 在碰撞仿真中的角色

碰撞是典型的瞬态、强非线性、大变形问题:接触状态在毫秒内反复切换,材料进入塑性甚至断裂,结构出现屈曲与折叠。隐式方法在这种工况下收敛困难,而 LS-DYNA 以显式中心差分法为主,配合丰富的接触算法、材料本构和沙漏控制,成为行业事实标准。它的核心能力体现在:

  • 显式时间积分,稳定时间步由最小单元尺寸与材料波速决定,适合毫秒级瞬态求解;
  • 自动接触与罚函数接触,可处理车身自接触、焊点失效后的部件分离;
  • 材料库覆盖金属塑性、复合材料、泡沫、橡胶、玻璃及点焊、胶粘等连接模型;
  • 与假人模型(LS-DYNA 关键字格式)、气囊折叠与展开(CV/ALE 或 FPM)配合完成乘员约束系统分析。

在整车开发流程中,LS-DYNA 通常承担三类任务:整车碰撞耐撞性分析、乘员约束系统与假人响应分析、以及行人保护与低速碰撞。它与前后处理工具(如 ANSA、HyperMesh、LS-PrePost)以及优化工具(LS-OPT、optiSLang)共同构成仿真链。

典型工作流程

一个完整的碰撞仿真项目,可以概括为“几何—网格—连接—边界—求解—评价—迭代”的闭环。下面这张图展示了从 CAD 到设计决策的数据流与关键判断点,其中每一步的输出都是下一步的输入,而结果评价会反向驱动结构或约束系统的修改。

整车碰撞仿真典型工作流程与数据流 CAD 几何 白车身/内外饰/底盘 几何清理与网格 目标尺寸 5–8 mm 连接建模 焊点/胶粘/铰链 材料与属性 *MAT / 应变率 假人与壁障 THOR/WorldSID 边界与初始 速度/重力/接触 控制卡片 *CONTROL_TIMESTEP LS-DYNA 求解 显式中心差分 结果提取 d3plot / binout 伤害指标评价 HIC / 侵入量 / 力 满足法规 与 NCAP? 结构/约束优化 厚度·材料·吸能盒 否 迭代重算 设计冻结与试验验证 滑车/实车对标 是 关键输入:几何、网格尺寸、材料曲线、焊点失效参数、壁障与假人模型、初速度与接触定义。 关键输出:B 柱/门槛侵入量、加速度时程、假人 HIC 与胸部压缩量、能量曲线与沙漏能占比。

流程中有几个容易被忽视但决定成败的环节。第一是网格尺寸与时间步:显式求解的稳定时间步约为 Δt = L/c,其中 L 是最小单元特征长度,c 是材料声速。若最小单元只有 2 mm,钢中声速约 5000 m/s,则 Δt 约 0.4 μs;一个 100 ms 的碰撞过程需要约 25 万步,计算量巨大。因此工程上常把非关注区域网格放大,或使用质量缩放(*CONTROL_TIMESTEP 的 DT2MS)人为提高时间步,但必须监控质量增加比例,一般控制在 5% 以内。

第二是接触与沙漏。碰撞中大量部件相互挤压,若接触定义遗漏,会出现穿透导致结果失真;而单点积分单元容易产生沙漏模式,需通过 *HOURGLASS 控制并检查沙漏能,通常要求沙漏能不超过总内能的 5%。第三是材料应变率效应:碰撞应变率可达 10²–10³ s⁻¹,低碳钢屈服强度可提升 30%–50%,必须使用 *MAT_PIECEWISE_LINEAR_PLASTICITY 等带应变率曲线的模型,否则吸能量会被低估。

教学示例:某简化正面 100% 刚性壁障模型,整车质量 1500 kg,初速度 50 km/h(13.89 m/s),总动能约 144.7 kJ。若仿真结束时总内能约 130 kJ、沙漏能约 4 kJ、剩余动能约 8 kJ,则沙漏能占比约 3%,能量平衡合理,结果可用于初步评价。以上数值为教学假设,不代表任何实际车型结果。

结果评价阶段,除了看假人伤害指标是否低于限值,还要看车体变形模式是否符合预期:前纵梁应发生有序的压溃折叠,而不是整体弯曲失稳;乘员舱应保持完整,A 柱、B 柱、门槛的侵入量应控制在目标范围内。若发现某处过早失效或吸能不足,通常回到连接建模或材料卡片修改,再重新求解,形成迭代闭环。最终,仿真结论需要与滑车试验或实车碰撞试验对标,用试验数据修正模型,才能让仿真真正具备预测能力。

注意:法规与 NCAP 的工况参数、假人类型和限值会随版本更新,建模前务必确认所依据的法规版本与评价规程年份,避免用旧模型套新要求。