电池热失控模拟
电池热失控的理论基础
电池热失控就是爆炸的那个吗？ 能通过模拟预测吗？
简单来说，锂离子电池内部温度超过约130°C时，SEI（固体电解质界面）层开始分解。露出来的负极锂直接与电解液反应放热，温度进一步升高。200°C以上时，正极放热反应和电解液放热反应失控，数秒内温度可超过800°C。
数秒就能到800°C! 如果装在车里就危险了…
正是。1个电池失控后热量传到相邻电池，模块整体→电池组整体火灾呈链式蔓延。这叫"热传播（thermal propagation）"，是EV火灾的主要原因。所以CAE用来预测传播时间，优化断热屏障厚度和冷却设计就很重要了。
热失控的机制
能详细说说温度上升的"阶段"吗？
温度上升是4阶段连锁反应进行的。模拟模型也标准地用Arrhenius式分别表述这4个阶段。
| 阶段 | 温度范围 | 反应 | 放热量目安 |
|---|---|---|---|
| 第1阶段 | 80～120°C | SEI膜分解 | ～250 J/g |
| 第2阶段 | 120～200°C | 负极-电解液反应（镶嵌Li露出） | ～1700 J/g |
| 第3阶段 | 180～250°C | 正极分解（NMC、LFP差异大） | ～800 J/g |
| 第4阶段 | 200°C～ | 电解液分解（有机溶媒气化·燃烧） | ～400 J/g |
SEI层分解是引火芯吧。LFP比NMC安全，是因为第3阶段正极分解反应温和？
完全正确。LiFePO₄（LFP）的Fe-O键很强，放出氧气温度高。所以正极引起的放热少，难以失控。但不是"绝对不会失控"，外部短路或过充时被迫进入第2阶段之后也会有问题。模拟中按正极材对Arrhenius常数分别建模。
支配方程
热失控模拟的支配方程是什么？
基本就是能量守恒式（非稳态热传导方程）加上化学反应放热项和焦耳发热项：
其中 $\rho$ 是密度，$c_p$ 是比热，$k$ 是热导率。右边第1项是热传导，$Q_{\text{rxn}}$ 是化学反应放热，$Q_{\text{ohm}}$ 是焦耳发热。
Arrhenius反应速率模型
$Q_{\text{rxn}}$ 内部是什么？ Arrhenius式是什么…
每个阶段 $i$ 的反应速率用Arrhenius式记述，乘以反应热 $H_i$ 后累加就是 $Q_{\text{rxn}}$：
各参数的意义如下：
- $A_i$：频率因子（前指数因子）[1/s]
- $E_{a,i}$：活性化能 [J/mol]
- $R$：气体常数 = 8.314 J/(mol·K)
- $c_i$：反应物的无量纲浓度（0～1，消耗后减小）
- $n_i$：反应阶数
- $H_i$：反应焓 [J/m³]
$\exp(-E_a/RT)$ 这部分是关键。温度 $T$ 上升时反应速率指数增长，放热增加，温度进一步上升……就是失控循环！
没错，这就是Semenov理论所说的"热爆炸"的本质。放热跟不上散热的瞬间，温度就会因正反馈急剧上升。
| 阶段 | $E_a$ [kJ/mol] | $A$ [1/s] | $H$ [J/g] |
|---|---|---|---|
| SEI分解 | 135 | 1.67×10¹⁵ | 257 |
| 负极-电解液 | 110～140 | 2.50×10¹³ | 1714 |
| 正极分解（NMC） | 170～200 | 1.75×10⁹ | 792 |
| 电解液分解 | 140～160 | 5.14×10²⁵ | 420 |
这些值一般从DSC（差示扫描量热计）或ARC（加速速率量热计）的实验数据拟合得出。
电池的放热模型
$Q_{\text{ohm}}$ 的焦耳发热是什么？ 正常工作时的发热和失控时的发热不同？
正常工作时的发热模型叫Bernardi式。这对热失控模拟的初值条件也很重要：
- 第1项 $I^2 R_{\text{int}}$：焦耳发热（电阻发热）。总是正的，大电流时越大。内阻 $R_{\text{int}}$ 是SOC和温度的函数。
- 第2项 $I \cdot T \cdot dU_{\text{OCV}}/dT$：熵发热（可逆热）。充放电时符号相反。$dU_{\text{OCV}}/dT$ 是OCV（开路电压）的温度系数。
例如快速充电（2C～3C）时，第1项焦耳发热占主导，电池温度能升到60～80°C。冷却不足会接近SEI分解的引发温度。
原来是把正常工作时的发热模型和失控时的Arrhenius反应模型结合在一起啊。
电池间热传播模型
1个电池失控后，热量如何传到相邻电池的？怎样建模？
电池间的传播模型考虑3种传热模式，全都要算：
- 热传导 $q_{\text{cond}}$：电池间直接接触或通过屏障材传导。最占主导。通过屏障材的热导率 $k$ 和厚度 $d$ 控制。
- 对流 $q_{\text{conv}}$：高温气体泄放（喷出）导致模块内强制对流。泄放气温度达600～1000°C。
- 辐射 $q_{\text{rad}}$：按 $\sigma T^4$ 的Stefan-Boltzmann辐射。800°C以上不能忽略。
实务中设置"触发电池"（钉刺或内部短路），看其旁边电池温度超过失控起始温度需要多长时间，称为"传播时间"。这个时间能否确保乘员逃生时间（通常5分钟以上）是设计基准。
屏障厚度稍微改变就会大幅影响传播时间啊。这种问题很适合模拟。
完全同意。某车企的例子，把陶瓷纤维屏障从5mm改到8mm，传播时间从约40秒延长到120秒以上。那80秒的差异关乎乘员安全。这种参数化分析正是FEM的长处。
热失控的"引火芯"——130°C会发生什么
理解锂离子电池热失控的关键是温度与反应的"连锁"。SEI膜在80～120°C开始分解，隔膜在130～150°C熔融导致内部短路，正极在180～200°C释放氧气——这些放热反应相继触发，温度数秒内超过800°C。但失控的起始温度受电池制造偏差和劣化状态影响很大，"所有电池都在同一温度失控"并不成立。这种个体差异的建模，是热失控模拟最难的课题之一。
电池热失控的数值计算方法
理论懂了。这个方程在计算机上怎样实际求解？
电化学-热耦合的非稳态问题，空间离散用FEM（有限元法），时间离散用隐式方案（后退Euler法或Crank-Nicolson法）是标准做法。Arrhenius项的强非线性在数值上很头疼，这是热失控模拟特有的难点。
FEM表述
FEM的表述怎样？ 和结构分析的热传导FEM类似吗？
用形函数 $N_i$ 近似温度场，变换为弱形式的过程和普通热传导FEM一样。用Galerkin法得到：
- $[C] = \int_\Omega \rho c_p [N]^T [N] \, d\Omega$：热容量矩阵
- $[K] = \int_\Omega [B]^T k [B] \, d\Omega$：热导率矩阵
- $\{Q\}$：放热向量（$Q_{\text{rxn}} + Q_{\text{ohm}}$ 的积分）
关键是右边 $\{Q\}$ 强烈依赖温度 $\{T\}$。因为有Arrhenius项 $\exp(-E_a/RT)$，放热量是温度的非线性函数。所以迭代法（Newton-Raphson法）必须的。
时间积分方案
失控开始温度变化急剧。时间步长怎样决定？
这是实务中最常烦扰的地方。失控前 $\Delta t = 1～10$ 秒就够了，但失控开始后温度梯度陡峭，非得 $\Delta t = 0.001～0.01$ 秒才跟得上。
| 方案 | 稳定性 | 精度 | 推荐场景 |
|---|---|---|---|
| 后退Euler（1阶隐式） | 无条件稳定 | $O(\Delta t)$ | 稳健但扩散大。初期筛选向导 |
| Crank-Nicolson（2阶隐式） | 无条件稳定 | $O(\Delta t^2)$ | 精度重视。失控前的缓和升温适合 |
| BDF2（2阶后退差分） | A-稳定 | $O(\Delta t^2)$ | Stiff系统强悍。失控过渡期有效 |
实务中**自适应时间步控制**必需。温度变化量 $|\Delta T|$ 超过阈值（比如5°C/步）时自动把 $\Delta t$ 减半，缓和阶段则增大。COMSOL Multiphysics等商业软件都有这个功能。
非线性耦合的处理
用Newton-Raphson法迭代，Arrhenius项是指数函数，收敛困难吗？
问得好。Arrhenius项对温度的敏感性极高。雅可比矩阵（切线刚性矩阵）中反应速率的温度导数：
失控区域 $\partial Q/\partial T$ 极大，雅可比条件数恶化。对策有：
- 荷载步细分：失控起始自动加入子步
- 线性搜索法：Newton增量加阻尼防止越界
- 反应物浓度钳制：强制 $c_i$ 在 [0, 1] 范围内防非物理值
网格策略
网格划分要注意什么？ 电池是卷绕结构，有几百层…
实务中不可能把几百层都单独划网。得用多级方法：
| 尺度 | 模型 | 典型单元数 | 网格类型 |
|---|---|---|---|
| 电极层（μm） | 1D电化学模型（P2D模型） | 100～1,000 | 1D线段单元 |
| 单体电池（mm） | 均质各向异性3D FEM | 10,000～50,000 | 六面体推荐 |
| 模块（cm） | 3D FEM（电池+屏障+壳体） | 100,000～500,000 | 四面体/六面体混合 |
| 整包（m） | 等效电路模型+简化3D | 500,000～2,000,000 | 粗六面体 |
电池间屏障的厚度方向最少要4层单元。厚度1～3mm的薄屏障温度梯度陡峭，这里省网格的话传播时间预测精度大幅下降。
电池热失控的实务应用
分析流程
实际上分析"电池会不会失控"用什么步骤？
标准流程5步：
- 电池级校准：用ARC试验数据拟合Arrhenius常数（$A_i, E_{a,i}, H_i$）
- 正常工作热模型验证：充放电循环中电池温度实测值与Bernardi式模型的一致性确认
- 触发事件建模：钉刺（局部短路电阻设定）、过充（SOC强制升高）、外部加热（加热器边界条件）
- 模块级传播分析：从触发电池到相邻电池温度的时间序列评估
- 参数化研究：屏障材厚度·材质、冷却流量、电池间间隙等设计变量变化，传播时间优化
第1步校准最困难吧。没实验数据就没法开始？
没错。ARC或DSC实验必需最低限度。但文献数据也积累了，NMC111/NMC622/NMC811/LFP等各正极材的代表Arrhenius常数已公开。初期可用文献值验证，再用自家电池数据微调是现实的做法。
材料参数的获取
Arrhenius常数以外还需什么材料参数？
| 参数 | 取得方法 | 注意事项 |
|---|---|---|
| 热导率（面方向/厚度方向） | 激光闪光法、热磁盘法 | 各向异性超大（40倍差异） |
| 比热 $c_p$ | DSC | 温度依存性要考虑（尤其电解液相变） |
| 内阻 $R_{\text{int}}(T, \text{SOC})$ | EIS（电化学阻抗谱） | 需要SOC和温度的2变量数据表 |
| OCV温度系数 $dU/dT$ | 低速充放电的OCV-T测量 | 正极材差异大 |
| 屏障材的热物性 | 厂家规格书+验证试验 | 高温（>300°C）劣化要考虑 |
虚拟试验（UL 2580 / GB 38031）
UL 2580和GB 38031是规格试验吧？ 模拟能代替吗？
规格认证必须实物试验，完全代替不了。但虚拟试验在3个场景很有效：
- 事前筛选：10种断热设计全部实物试？ 成本太高。模拟先筛到3个候选，再实物试验。
- 最坏情况特定："哪个位置电池失控传播最快？"FEM网格化搜索，实物不可能全试。
- 规格变更对应：新规格出现时，设计改动方向模拟即时评估。
| 规格 | 对象 | 主试验项目 | 模拟应用 |
|---|---|---|---|
| UL 2580 | EV用电池包 | 钉刺、过充、外火 | 传播时间预测、泄放气行为 |
| GB 38031 | 中国市场EV电池 | 热传播试验（5分钟非蔓延） | 屏障材参数化设计 |
| UN ECE R100 | 国际统一基准 | 火灾安全性 | 外火暴露模拟 |
| SAE J2464 | 试验手册 | 滥用试验方法标准化 | 试验条件的根据设定 |
GB 38031的"5分钟非蔓延"，就是1个电池失控5分钟内不蔓延到整包？
准确说是"触发电池热失控后到电池包外部出现火焰或爆炸为止，5分钟以上的余地"。乘员逃生时间的最低限保障。模拟中以这5分钟为目标值优化屏障设计。
断热屏障设计
断热屏障设计，厚一点更安全但重，薄点轻但危险…
正是矛盾。模拟找最优解的典型参数化设计例：
| 屏障材 | 热导率 [W/(m·K)] | 密度 [kg/m³] | 5mm厚传播时间 |
|---|---|---|---|
| 气凝胶片材 | 0.015 | 150 | 约180秒 |
| 陶瓷纤维 | 0.03 | 200 | 约120秒 |
| 云母片 | 0.4 | 2800 | 约50秒 |
| 无屏障（仅空气层） | 0.026（空气） | — | 约30秒 |
气凝胶断热最优但贵。陶瓷纤维性价比好，很多OEM选它。加主动冷却（冷却水路）还能进一步延长传播时间。
EV火灾的"阻止蔓延"——电池包设计者的斗争
热失控模拟实务的最大目标是"1个电池失控也别蔓延到整包"的设计。某国内车企，根据模拟结果把屏障材从陶瓷纤维改成气凝胶复合体，厚度从8mm减到5mm，电池包重量轻3kg，同时传播时间延长2倍。但最终设计确定前经过1000次以上的模拟和20次以上的实物钉刺试验。"不要模拟的话100次以上的实物试验会必要"——设计者的实感是。
电池热失控的软件对比
工具对比
电池热失控模拟具体用什么软件？
主流工具和特点对比如下：
| 工具 | 开发商 | 优点 | 电池热失控对应 |
|---|---|---|---|
| COMSOL Multiphysics | COMSOL AB | 多物理耦合简单。Arrhenius反应可GUI定义 | Battery Design Module（电化学+热+结构） |
| Ansys Fluent | Ansys Inc. | CFD耦合。泄放气行为分析 | MSMD（多尺度多区域）电池模型 |
| Simcenter STAR-CCM+ | Siemens | 多面体网格。大规模电池包分析 | Battery Simulation Module（ESC/NTGK模型） |
| GT-AutoLion | Gamma Technologies | 电化学详细模型（P2D） | 电池级高精度预测。系统模型连接 |
| LS-DYNA | LSTC / Ansys | 碰撞+热失控耦合。机械触发再现 | *BATTERY_ECHEM + 热solver耦合 |
LS-DYNA进来了！ 碰撞分析的软件？
车辆碰撞→电池壳体变形→电池压壊→内部短路→热失控，这一系列过程在LS-DYNA 1个模型追踪。机械触发到热失控无缝连接，其他工具难做到。车企碰撞安全部队常用它。
选型指南
结局哪个选好？ 预算也有限…
用途选就对：
- 电池级电化学+热详细分析 → COMSOL（直观GUI）或GT-AutoLion（P2D精度）
- 模块/电池包蔓延分析主目的 → Ansys Fluent或STAR-CCM+（大规模3D热+CFD）
- 碰撞→热失控一贯 → LS-DYNA一択
- 成本抑制 → OpenFOAM自作solver加Arrhenius项。但开发时间要有心理准备
多尺度分析的工具使分
热失控模拟工具选型常忽视的点是"能跨尺度计算吗"。电池内部电化学反应是微米尺度，模块热传导是厘米尺度，电池包冷却流道是毫米尺度——这些都用1个模型算计算成本极高。实务做法是"电化学用GT-AutoLion的P2D模型，蔓延用Fluent/STAR-CCM+的3D模型，系统全体用MATLAB/Simulink"这样工具分工，数据相互转递。这种多工具流程是通常做法。
电池热失控的前沿研究
机器学习代理模型
AI/机器学习预测热失控的研究也有吗？
有，大2方向：
- 代理模型：FEM数千案例回数据作成、（DNN）過程回帰（GPR）近似模型構築。参数优化的評価関数使、毎回FEM走100～1000倍高速化。
- PINN（物理）：保存式損失関数組込物理則満解学習。実験数据少場合物理的妥当外挿可。失控域的急激温度変化PINN捕捉的研究段階。
代理模型材最適设计面白！
数字孪生·BMS联动
車走時熱失控予測？
的考方。BMS（管理系统）取得電圧·電流·温度数据、事前構築模型入力、「何分熱失控至」推定。
例温度異常上昇場合、模型「条件120秒後SEI分解開始」予測、BMS自動的急速放電停止→冷却強化→乗員警告実行。行CAE安全制御直結。
、仿真设计運用！
电池热失控的故障排查
常见错误与对策
先生、熱失控仿真実際何故障遭遇？
最常是以下4：
| 症状 | 原因 | 対策 |
|---|---|---|
| 失控開始付近収束失敗 （"小"） | Arrhenius項的急激温度感度悪化 | 適応時間刻的下限 $10^{-5}$ s 程度設定。Newton反復最大回数50→200増加。探索法有効化 |
| 温度負値 | 反応物浓度 $c_i$ 流程放熱異常 | $c_i$ [0, 1] 追加。COMSOL: "Constraints" 下限設定 |
| 伝播時間実験大幅 | 材高温劣化（熱導率的温度依存性）無視 | 300°C以上的 $k(T)$ 入力。700°C超焼結熱導率急増 |
| 温度上昇遅 | 熱導率異方性設定（面方向厚方向逆） | 面方向 $k_{\parallel} \approx 30$ W/(m·K)、厚方向 $k_{\perp} \approx 0.7$ W/(m·K) 確認 |
異方性設定… 面方向厚方向40倍違的。
。銅·箔的集電体面方向熱通高値。厚方向·電解液的低熱導率支配的。的異方性等方性近似、失控開始時刻30%以上。
最後「覚」？
3。
- ：自社的ARC/DSC数据信頼熱失控仿真不可能。文献値初期验证用的。
- 材的温度依存性忘：室温熱物性模型組的最危険失敗。高温域材料特性劇変。
- 最悪案例探：「真中失控最危険」限。冷却结构近失控冷媒蒸発冷却喪失案例。分析漏。
！ 熱失控仿真全体像。COMSOL的Battery Module始！
！ COMSOL的Application Gallery "Lithium-Ion Battery Pack Thermal Runaway Propagation" 。的記事説明4段階Arrhenius模型的実装、初歩最適。
更多
错误