车辆碰撞仿真详解
理论与物理
车辆碰撞仿真
老师，汽车的碰撞安全没有FEM就无法设计，对吧？
没错。在现代汽车开发中，在进行实车碰撞试验之前，会使用FEM进行数百至数千次仿真。碰撞安全的设计是由FEM主导的。
碰撞的分类
| 碰撞类型 | 标准 | 速度 | 特征 |
|---|---|---|---|
| 正面碰撞（全宽） | FMVSS 208, Euro NCAP | 56 km/h | 全宽刚性壁碰撞 |
| 正面碰撞（偏置） | Euro NCAP, IIHS | 64 km/h | 40%偏置ODB |
| 侧面碰撞 | FMVSS 214, Euro NCAP | 50 km/h | 可变形壁障侧面碰撞 |
| 后面碰撞 | FMVSS 301 | 80 km/h | 防止燃油泄漏 |
| 柱状侧面碰撞 | Euro NCAP | 32 km/h | 电线杆等细长障碍物 |
| 行人保护 | Euro NCAP | — | 引擎盖上的头部冲击 |
竟然有这么多种碰撞模式啊。
一个车型要仿真20到50个碰撞工况。每个工况使用数百万单元的整车模型进行显式求解，计算50到200毫秒。计算资源是巨大的。
FEM模型的规模
典型的整车碰撞模型：
| 项目 | 数值 |
|---|---|
| 单元数量 | 300万～1000万 |
| 节点数量 | 100万～500万 |
| 材料模型数量 | 50～200 |
| 接触定义数量 | 数百 |
| 计算时间 | 4～24小时（100～200 CPU） |
| 结果文件大小 | 10～100 GB |
1000万个单元！规模真惊人。
包含BIW（白车身）、覆盖件、底盘、动力总成、内饰、座椅、假人、安全气囊……所有部件。网格生成可能需要数周，计算设置也可能需要数天。
碰撞安全的设计思想
能量吸收是碰撞安全的基本概念：
1. 前部碰撞吸能区 — 通过受控屈曲吸收能量
2. 乘员舱 — 不变形的高刚性笼式结构
3. 约束系统 — 通过安全带、安全气囊使乘员减速
“应该溃缩的部分”和“不能溃缩的部分”是设计的核心呢。
用FEM来仿真这种“受控屈曲”。通过FEM优化碰撞盒的加强筋形状、板厚、材料，以达到目标的能量吸收量和减速度脉冲。
总结
要点：
- 使用FEM仿真20～50个碰撞工况 — 在实车试验之前
- 300万～1000万单元的整车模型 — LS-DYNA显式求解
- 能量吸收与乘员舱变形限制 — 碰撞安全的设计思想
- 碰撞盒的受控屈曲 — 通过FEM优化
- 符合Euro NCAP、FMVSS等标准 — 多种场景
碰撞安全工程学由Hugh DeHaven创立
被誉为现代碰撞安全工程学之父的Hugh DeHaven在1942年提出了“可变碰撞区”的概念。汽车在撞击障碍物时，通过使发动机舱发生有意的变形来吸收能量，从而保护客舱，这一思想是现今所有汽车上实现的溃缩区设计的原型。福特于1956年首次在量产车上采用DeHaven理论的带衬垫仪表板，也源于同样的思路。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。您有过急刹车时身体被向前甩出的经历吗？那种“被带走的感觉”正是惯性力。物体越重，越难启动，一旦启动也越难停止。地震时建筑物摇晃，也是因为地面突然移动，而建筑物的质量“被落下”。静力分析中此项设为零，那是基于“缓慢施力，加速度可忽略”的假设。对于冲击载荷或振动问题，此项绝对不能省略。
- 刚度项（弹性恢复力）：$Ku$ 或 $\nabla \cdot \sigma$。拉弹簧时会感觉到“想要恢复原状的力”吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——铁棒和橡皮筋，用相同的力拉，哪个伸得更长？当然是橡皮筋。这种“不易伸长性”就是杨氏模量 $E$，它决定了刚度。常见的误解是：“刚度高=强度高”。刚度是“不易变形的程度”，强度是“不易破坏的程度”，这是两个不同的概念。
- 外力项（载荷项）：体积力 $f_b$（重力等）和表面力 $f_s$（压力、接触力等）。可以这样想——桥上卡车的重量是“作用在整个内部上的力”（体积力），轮胎压路面的力是“只作用在表面上的力”（表面力）。风压、水压、螺栓预紧力……这些都是外力。这里容易犯的错误是：弄错载荷方向。本想施加“拉力”却成了“压力”——听起来像笑话，但在3D空间中坐标系发生旋转时，确实会发生。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。试着弹一下吉他的弦。声音会一直持续吗？不，会逐渐变小。这是因为振动能量通过空气阻力和弦的内部摩擦转化成了热能。汽车的减震器也是同样的原理——特意吸收振动能量来改善乘坐舒适性。如果阻尼为零会怎样？建筑物在地震后会一直摇晃下去。实际上不会这样，所以设置适当的阻尼很重要。
假设条件与适用范围
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 输入为mm时，载荷·弹性模量也需统一为MPa/N系 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意工程应变与对数应变的区别（大变形时） |
| 弹性模量 $E$ | Pa | 钢：约210 GPa，铝：约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm系中为tonne/mm³（钢约为 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm系用N，m系也用N统一 |
数值解法与实现
碰撞仿真的FEM
请告诉我碰撞仿真的技术细节。
单元类型
- BIW（白车身） — 壳单元（主要为Quad4, HEX8R）
- 覆盖件 — 壳单元
- 保险杠、纵梁 — 壳单元+实体单元
- 假人 — 壳单元+实体单元+1D单元（关节）
- 安全气囊 — 壳单元+气体模型（ALE/CPM）
材料模型
- 钢板 — MAT24（弹塑性）+ 应变率相关性（Cowper-Symonds）
- 铝 — MAT24 或 MAT125
- 树脂 — MAT24 或 MAT89
- CFRP — MAT54/58（渐进损伤）
- 橡胶 — MAT77（Ogden超弹性）
- 泡沫 — MAT57/63（可压缩泡沫）
应变率相关性很重要呢。
碰撞时的应变率为 $10 \sim 1000$ /s。钢的屈服强度会因应变率而上升20～50%。忽略此效应会低估能量吸收。Cowper-Symonds法则：
接触
碰撞模型需要数百个接触定义。LS-DYNA的*CONTACT_AUTOMATIC_GENERAL（全局自动接触）是标准。使用罚函数法防止穿透。
总结
1000万单元模型2小时即可求解的时代
现代汽车整车碰撞模型的规模达到单元数700万～1200万，材料定义5000以上，接触对200以上。截至2024年，使用256核（例如：AMD EPYC 9354·128核×2节点）运行LS-DYNA MPP，100ms的全宽正面碰撞分析大约2～4小时即可完成。丰田和大众等公司通过在一夜批处理中同时执行多个试验模式的“夜间运行”体制，大幅缩短了开发周期时间。
线性单元（1次单元）
节点间线性插值。计算成本低，但应力精度低。注意剪切锁定（可通过减缩积分或B-bar法缓解）。
2次单元（带中间节点）
可以表现曲线变形。应力精度大幅提高，但自由度约增加2～3倍。推荐：应力评估很重要的情况。
完全积分 vs 减缩积分
完全积分：有过约束（锁定）风险。减缩积分：沙漏模式（零能模式）风险。根据情况选择。
自适应网格
基于误差指标（ZZ估计量等）的自动细化。有效提高应力集中区域的精度。有h法（单元细分）和p法（增加阶次）。
牛顿·拉弗森法
非线性分析的标准方法。每次迭代更新切线刚度矩阵。在收敛半径内具有二次收敛性，但计算成本高。
修正牛顿·拉弗森法
切线刚度矩阵使用初始值或每隔几次迭代更新。每次迭代成本低，但收敛速度为线性。
收敛判定标准
力残差范数: $||R|| / ||F_{ext}|| < \epsilon$（通常 $\epsilon = 10^{-3}$～$10^{-6}$）。位移增量范数: $||\Delta u|| / ||u|| < \epsilon$。能量范数: $\Delta u \cdot R < \epsilon$
载荷增量法
不一次性施加全部载荷，而是分小步增加。弧长法（Riks法）可以超越载荷-位移关系的极值点进行追踪。
直接法 vs 迭代法的比喻
直接法是“用笔算精确解联立方程”的方法——可靠但大规模问题耗时过长。迭代法是“反复猜测逼近正确答案”的方法——最初是粗略的答案，但每次迭代精度都会提高。就像查字典时，从第一页开始按顺序找（直接法）不如先估计大概位置翻开，再前后调整（迭代法）更高效，原理相同。
网格阶次与精度的关系
1次单元是“用直尺近似曲线”——用直线折线表现，精度有限。2次单元是“柔性曲线”——可以表现曲线变化，即使网格密度相同，精度也显著提高。但是，每个单元的计算成本增加，因此需要根据总体的成本效益来判断。
实践指南
碰撞仿真实务
请告诉我碰撞仿真的工作流程。
工作流程
1. 接收CAD数据 — 整合各部件的CAD
2. 网格生成 — 使用HyperMesh/ANSA生成壳网格（5～10 mm）
3. 材料定义 — 根据材料试验数据设置MAT24等
4. 连接部位的建模 — 点焊（*CONSTRAINED_SPOTWELD）、粘接、螺栓
5. 假人布置 — WorldSID/THOR等认证假人模型
6. 约束系统 — 安全带（*ELEMENT_SEATBELT）、安全气
相关主题
なった
詳しく
報告