什么是电池热失控？能否通过仿真进行预测？

当锂离子电池内部温度超过约130°C时，SEI层开始分解；超过200°C时，正极与电解液之间的放热反应将发生失控。通过基于有限元法（FEM）的电化学-热耦合分析，可以预测热失控的传播时间，并用于隔热屏障的设计。

电池热失控的发热模型通常使用哪些方程？

主要采用阿伦尼乌斯型反应速率方程 k=A·exp(-Ea/RT) 以及基于发热量 Q=I²R+I·T·(dU/dT) 的能量守恒方程作为控制方程。其中，SEI分解、负极-电解液反应、正极分解和电解液分解的四阶段反应模型是标准建模方法。

如何对电芯之间的热传播进行建模？

以单个电芯的热失控作为触发条件，通过包含电芯间热传导、辐射和对流的三维非稳态有限元分析，计算整个模组的温度传播过程。通过改变隔热屏障材料的厚度和材质进行参数化分析，可实现防止热蔓延的设计。

能否通过仿真替代UL 2580或GB 38031的测试？

虽然无法完全替代实际测试，但虚拟测试可用于测试条件的前期筛选、隔热材料设计的优化以及最坏工况的识别。尽管法规认证仍需实物测试，但仿真可显著减少所需测试的次数。

电池热失控模拟

分类: 連成解析 > 電気化学連成 | 更新 2026-04-12

理论与物理

🧑‍🎓

电池的热失控就是会爆炸的那个吗？能用仿真预测吗？

🎓

简单来说，锂离子电池内部温度超过约130°C时，SEI（固体电解质界面）层开始分解。然后暴露出来的负极锂会与电解液直接反应放热，导致温度进一步升高。超过200°C时，正极的氧气释放和电解液的放热反应会失控，有时在几秒钟内温度就能超过800°C。

🧑‍🎓

诶，几秒钟就800°C！？如果装在车里的电池发生这种情况，那可就危险了…

🎓

没错。一个电芯的热失控如果通过热传播影响到相邻电芯，就会导致整个模组→整个电池包连锁性地发生火灾。这就是“延烧（thermal propagation）”，也是电动汽车火灾的主要原因。因此，通过CAE预测传播时间，并优化隔热屏障的厚度和冷却设计就变得至关重要。

热失控的机制

🧑‍🎓

能再详细讲讲温度上升的“阶段”吗？

🎓

温度上升是通过4个阶段的连锁反应进行的。在仿真模型中，通常用Arrhenius公式分别描述这4个阶段。

阶段	温度范围	反应	发热量参考
Stage 1	80〜120°C	SEI膜分解	~250 J/g
Stage 2	120〜200°C	负极-电解液反应（嵌入Li暴露）	~1700 J/g
Stage 3	180〜250°C	正极分解（NMC, LFP差异很大）	~800 J/g
Stage 4	200°C〜	电解液分解（有机溶剂汽化・燃烧）	~400 J/g

🧑‍🎓

也就是说，SEI层的分解是导火索，然后反应会接连连锁发生。听说LFP比NMC更安全，是因为Stage 3的正极分解比较温和吗？

🎓

没错。LiFePO₄（LFP）的Fe-O键很强，氧气释放温度较高。因此正极引起的发热较少，不易发生热失控。但并非“绝对不会失控”，在外部短路或过充电等情况下，仍有可能被迫进入Stage 2及以后。仿真中会通过改变不同正极材料的Arrhenius常数来建模。

控制方程

🧑‍🎓

热失控仿真的控制方程是怎样的？

🎓

基本是在能量守恒方程（非稳态热传导方程）的基础上，加上化学反应发热项和焦耳发热项：

$$ \rho c_p \frac{\partial T}{\partial t} = \nabla \cdot (k \nabla T) + Q_{\text{rxn}} + Q_{\text{ohm}} $$

其中 $\rho$ 是密度，$c_p$ 是比热容，$k$ 是热导率。右边第一项是热传导，$Q_{\text{rxn}}$ 是化学反应发热，$Q_{\text{ohm}}$ 是焦耳发热。

Arrhenius反应速率模型

🧑‍🎓

$Q_{\text{rxn}}$ 具体是什么？经常听到Arrhenius公式…

🎓

用Arrhenius公式描述每个阶段 $i$ 的反应速率，再乘以反应热 $H_i$ 后求和。这就是 $Q_{\text{rxn}}$：

$$ Q_{\text{rxn}} = \sum_{i} H_i \cdot A_i \exp\!\left(-\frac{E_{a,i}}{RT}\right) c_i^{n_i} $$

各参数的含义如下：

$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公式。它作为热失控仿真的初始条件也很重要：

$$ Q_{\text{cell}} = I^2 R_{\text{int}} + I \cdot T \frac{dU_{\text{OCV}}}{dT} $$

第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反应模型结合起来。

电芯间热传播模型

🧑‍🎓

一个电芯失控后，火势蔓延到相邻电芯的过程是如何建模的？

🎓

电芯间的传播模型需要考虑三种传热模式：

$$ q_{\text{prop}} = q_{\text{cond}} + q_{\text{conv}} + q_{\text{rad}} $$

热传导 $q_{\text{cond}}$：通过电芯间直接接触或屏障材料进行的传导。占主导地位。通过屏障材料的热导率 $k$ 和厚度 $d$ 来控制。
对流 $q_{\text{conv}}$：高温气体喷出（Vent）引起的模组内强制对流。喷出气体温度可达600〜1000°C。
辐射 $q_{\text{rad}}$：与 $\sigma T^4$ 成正比的斯蒂芬-玻尔兹曼辐射。超过800°C时不可忽略。

实际工作中，会设定一个“触发电芯”（如针刺或内部短路），然后评估其相邻电芯温度超过失控开始温度所需的时间，即“传播时间”。这个传播时间能否确保乘员逃生时间（通常5分钟以上）是设计基准。

🧑‍🎓

稍微改变一下屏障材料的厚度，传播时间似乎就会有很大变化。这确实是适合仿真的问题。

🎓

是的。某汽车制造商的案例显示，仅将陶瓷纤维系屏障从5mm改为8mm，传播时间就从约40秒延长到120秒以上。这80秒的差距关系到乘员的安全。这种参数化分析正是FEM的强项。

Coffee Break 闲谈

热失控的“导火索”——130°C时会发生什么

理解锂离子电池热失控的关键在于温度与反应的“连锁”。SEI膜开始分解的80〜120°C，隔膜熔融导致内部短路的130〜150°C，正极释放氧气的180〜200°C——这些放热反应会接连触发，温度有时在几秒钟内就能超过800°C。然而，发热的时机强烈依赖于电池的制造差异和劣化状态，因此“所有电池都在相同温度下开始失控”的情况并不存在。对这种个体差异的建模，是热失控仿真中最困难的课题之一。

各项的物理意义

蓄热项 $\rho c_p \partial T / \partial t$：表示电池材料（正极・负极・电解液・集流体）的热惯性。对于卷绕结构，由于存在各向异性，有时会在面方向和厚度方向设置不同的有效热容。
热传导项 $\nabla \cdot (k \nabla T)$：圆柱电芯的面方向（ $k_r \approx$ 30 W/(m·K)）和径向（$k_z \approx$ 0.7 W/(m·K)）热导率差异超过40倍，具有很强的各向异性。软包电芯的面方向和厚度方向也有很大差异。
化学反应项 $Q_{\text{rxn}}$：Arrhenius型的4阶段反应模型。温度上升→反应加速→温度进一步上升的正反馈是热失控的本质。
电发热项 $Q_{\text{ohm}}$：正常运行时用Bernardi公式（$I^2R + I T \, dU/dT$）描述。内部短路时，通过短路电阻的局部大电流焦耳发热占主导。

假设条件与适用范围

均质化假设：将卷绕内部的电极叠层结构视为均质各向异性材料。忽略电芯内部的局部反应不均匀性。
Arrhenius型反应速率：用简单的一级反应近似每个阶段。实际反应包含多个并行路径，但通过对DSC/ARC数据的拟合可以得到有效的模型。
反应物浓度假设：追踪无量纲浓度 $c_i$（0〜1）的消耗。有时也用未反应率表示。
内部短路模型：在再现针刺试验时，短路部位的电阻值和接触面积对结果影响很大。必须进行实验标定。
适用范围：气体喷出后的流体力学（喷出气体行为）需要另外用CFD求解，或作为简化的对流换热系数给出。

量纲分析与单位制

变量	SI单位	典型值・注意事项
活化能 $E_a$	J/mol	SEI: 135 kJ/mol、正极: 170〜200 kJ/mol
频率因子 $A$	1/s	10⁹〜10²⁵量级。不同阶段差异很大
反应焓 $H$	J/kg 或 J/m³	注意体积基准与质量基准的转换
内阻 $R_{\text{int}}$	Ω	SOC・温度依赖。约1〜100 mΩ
电芯间热导率	W/(m·K)	屏障材料: 0.03〜0.5、铝壳体: 150〜200

数值解法与实现

🧑‍🎓

理论我大概明白了。这些方程在计算机上具体是怎么求解的？

🎓

这是一个电化学-热耦合的非稳态问题，通常空间离散用FEM（有限元法），时间离散用隐式格式（后向欧拉法或Crank-Nicolson法）。Arrhenius项具有很强的非线性，这是热失控仿真特有的数值难点。

FEM公式化

🧑‍🎓

FEM的公式化是怎样的？热传导的FEM和结构分析类似吗？

🎓

用形函数 $N_i$ 近似温度场，并转换到弱形式，这一步和通常的热传导FEM是一样的。应用伽辽金法得到：

$$ [C]\left\{\dot{T}\right\} + [K]\{T\} = \{Q(\{T\})\} $$

$[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)$，发热量是温度的非线性函数。因此必须使用迭代法（牛顿-拉夫逊法）。

时间积分格式

🧑‍🎓

失控开始时温度变化很剧烈吧？时间步长是怎么决定的？

🎓

这是实际工作中最常遇到的难点。热失控前 $\Delta t = 1\text{〜}10$ 秒就足够了，但失控开始后温度梯度变得陡峭，需要将 $\Delta t$ 细化到 0.001\text{〜}0.01 秒才能跟上变化。

格式	稳定性	精度	推荐场景
后向欧拉（1阶隐式）	无条件稳定	$O(\Delta t)$	鲁棒但扩散较大。适合初期筛选
Crank-Nicolson（2阶隐式）	无条件稳定	$O(\Delta t^2)$	重视精度。适合失控前缓慢升温阶段
BDF2（2阶后向差分）	A-稳定	$O(\Delta t^2)$	对刚性系统强。在失控过渡期有效

实际工作中，自适应时间步长控制是必须的。当温度变化量 $|\Delta T|$ 超过阈值（例如5°C/步）时，自动将 $\Delta t$ 减半，在平缓阶段则增大。COMSOL Multiphysics等软件的标准功能就支持这个。

非线性耦合的处理

🧑‍🎓

您刚才说用牛顿-拉夫逊法收敛，但Arrhenius项是指数函数，会不会很难收敛？

🎓

问得很到位。Arrhenius项对温度的敏感性极高。因为雅可比矩阵（切线刚度矩阵）中包含反应速率对温度的微分：

$$ \frac{\partial Q_{\text{rxn}}}{\partial T} = \sum_i H_i A_i c_i^{n_i} \frac{E_{a,i}}{RT^2} \exp\!\left(-\frac{E_{a,i}}{RT}\right) $$

在失控区域，$\partial Q/\partial T$ 变得非常大，导致雅可比矩阵的条件数恶化。对策包括：

载荷步细分：检测到失控开始时自动增加子步
线搜索法：对牛顿更新量施加阻尼以防止超调
反应物浓度钳制：将 $c_i$ 强制约束在 [0, 1] 范围内，防止出现非物理值

网格策略

🧑‍🎓

划分网格时需要注意什么？电池卷绕结构有几百层吧？

🎓

在实际工作中，对卷绕结构的所有层进行单独网格划分在计算成本上是不可能的。因此采用以下多阶段方法：

尺度	模型	典型单元数	网格类型
电极层（μm）	1D电化学模型（P2D模型）	100〜1,000	1D线单元
单电芯（mm）	均质各向异性体3D FEM	10,000〜50,000	六面体