大型磷酸铁锂电池高温热失控模拟研究

doi:10.19799/j.cnki.2095-4239.2020.0249

大型磷酸铁锂电池高温热失控模拟研究

梅文昕^,¹, 段强领¹, 王青山²^,³, 李妍²^,³, 李欣⁴, 朱金大⁴, 王青松^,¹

1.中国科学技术大学，安徽合肥 230026

2.国网江苏省电力有限公司经济技术研究院，江苏南京 210008

3.国网江苏省电力有限公司，江苏南京 210024

4.国网电力科学研究院有限公司，江苏南京 211000

Thermal runaway simulation of large-scale lithium iron phosphate battery at elevated temperatures

MEI Wenxin^,¹, DUAN Qiangling¹, WANG Qingshan²^,³, LI Yan²^,³, LI Xin⁴, ZHU Jinda⁴, WANG Qingsong^,¹

1.University of Science and Technology of China, Hefei 230026, Anhui, China

2.Economic and Technological Research Institute of State Grid Jiangsu Electric Power Co. Ltd. , Nanjing 210008, Jiangsu, China

3.State Grid Jiangsu Electric Power Co. Ltd. , Nanjing 210024, Jiangsu, China

4.State Grid Electric Power Research Institute Co. Ltd, Nanjing 211000, Jiangsu, China

通讯作者: 王青松，教授，主要从事火灾安全领域相关研究，E-mail：pinew@ustc.edu.cn。

收稿日期: 2020-07-15 修回日期: 2020-09-10 网络出版日期: 2021-01-08

基金资助:

国家电网公司总部科技项目. 5400-201940486A-0-0-00

Received: 2020-07-15 Revised: 2020-09-10 Online: 2021-01-08

作者简介 About authors

梅文昕（1995—），女，博士研究生，主要从事锂离子电池电化学-热-应力耦合模型的模拟研究，E-mail:heart@mail.ustc.edu.cn； E-mail：heart@mail.ustc.edu.cn。

摘要

高温是触发锂离子电池热失控的最直接原因，因此研究锂离子电池在高温加热中的热失控特征及其内在机制至关重要。本文选取109 A·h大型磷酸铁锂电池为研究对象，在COMSOL Multiphysics中建立了6种不同温度下(140 ℃、145 ℃、150 ℃、155 ℃、160 ℃、165 ℃)的烘箱热失控模型，模拟分析了电池在高温加热条件下的热失控特征和温度分布。研究结果表明在140 ℃和145 ℃时电池未发生热失控，其他工况下电池均发生热失控，且环境温度越高，电池发生热失控的时间越早，温升速率加快。此外，通过对热失控各副反应分解浓度分析得知，未发生热失控情况下只发生了SEI膜和负极的分解反应，而正极与电解液的反应是造成电池热失控的主要诱因。最后通过对比发生热失控和未发生热失控情况下电池的温度分布，发现未发生热失控条件下电池温度分布均匀，而发生热失控时电池温度均匀性变差。同时发现高环境温度下电池的热失控更为剧烈，温度分布极不均匀，且在热失控前后温度分布变化较快，预计电池材料发生的不可逆分解反应是导致电池损坏的主要原因。

关键词： 锂离子电池安全 ; 磷酸铁锂 ; 高温热失控模型 ; 副反应 ; 温度分布

Abstract

Elevated temperature is the most direct trigger of thermal runaway in lithium-ion batteries, so it is crucial to study the thermal runaway characteristics and mechanism of lithium-ion batteries at elevated temperatures. This paper presents the study of 109 A·h large-scale lithium iron phosphate power batteries, and an oven thermal runaway model at six different temperatures (140 ℃, 145 ℃, 150 ℃, 155 ℃, 160 ℃, 165 ℃) is presented via COMSOL Multiphysics software to simulate the thermal runaway characteristics and temperature distribution of the battery under high temperatures. The studies showed that the battery does not trigger thermal runaway at temperatures of 140 ℃ and 145 ℃, but does trigger thermal runaway at higher temperatures. The higher the oven temperature, the earlier the thermal runaway occurs, and the rate of temperature rise is also accelerated. Through the analysis of the decomposition concentration of each side reaction in the thermal runaway, it is observed that only the decomposition of the solid electrolyte interphase layer and anode occurred in the non-thermal runaway cases. The reaction between the cathode and the electrolyte is the main cause of thermal runaway. Finally, by comparing the thermal runaway cases with the non-thermal runaway cases, it was found that the temperature distribution of the battery is uniform in the case of a non-thermal runaway, while the temperature uniformity is poor in the case of a thermal runaway. At higher temperatures, the thermal runaway of the battery is more severe, where the temperature distribution is extremely uneven and changes rapidly before and after thermal runaway. In such thermal runaways, it is predicted that the electrode materials have undergone irreversible decomposition leading to battery damage.

Keywords： lithium ion battery safety ; lithium iron phosphate ; oven thermal runaway model ; side reaction ; temperature distribution

PDF (4213KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

梅文昕, 段强领, 王青山, 李妍, 李欣, 朱金大, 王青松. 大型磷酸铁锂电池高温热失控模拟研究[J]. 储能科学与技术, 2021, 10(1): 202-209

MEI Wenxin. Thermal runaway simulation of large-scale lithium iron phosphate battery at elevated temperatures[J]. Energy Storage Science and Technology, 2021, 10(1): 202-209

随着锂离子电池尺寸、能量密度的不断增大，其安全性也越来越成为消费者关注的重点^[1]。锂离子电池在受到电滥用、机械滥用以及热滥用的情况下会因内部热量的积累造成高温，引发热失控^[2]，进而造成火灾爆炸事故。可以看出高温是触发锂离子电池热失控的最直接原因，因此科学认识高温作用下诱发的锂离子电池热失控特征至关重要。目前对锂离子电池高温热失控的研究局限于三元电池，如杜光超等^[3]分析了圆柱形高镍三元锂离子电池在高温下的热失控参数，黄文才等^[4]针对三元电池建立了高温加热条件下的三维分层物理热失控模型，研究了不同传热系数下电池的热失控行为。相比于三元电池，目前关于磷酸铁锂电池在高温下的热失控研究还比较缺乏。特别是，对于大型磷酸铁锂电池的高温热失控行为特征、内在机制的认识更是严重不足。鉴于此，本文选取大型磷酸铁锂电池(109 A·h)为研究对象，对其建立不同加热温度下的烘箱热失控模型，通过对比发生热失控和未发生热失控两种条件，得出磷酸铁锂电池高温热失控行为。

1 三维烘箱热失控模型的建立

锂离子电池的烘箱热失控模型即模拟电池在烘箱中恒温加热的情况下的热失控发展规律。锂离子电池热失控各副反应阶段见图1。在90～120 ℃时^[5]，负极表面形成的SEI膜发生分解，当SEI膜分解完成后，负极失去保护，嵌入在负极的锂继而与电解液发生反应。随着温度的持续升高，隔膜开始熔化，且正极与电解液反应，且伴随着电解液的自身分解，这时会发生隔膜破裂、正负极短路，而导致热失控的发生。之后负极与黏结剂发生反应，引发电解液着火，造成电池的燃烧。在烘箱热失控模拟中，电池热失控的最高温度较绝热热失控模拟低，也不会发生燃烧现象，因此只考虑前4种副反应的发生。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 锂离子电池热失控各副反应历程

Fig.1 The process of thermal runaway side reactions of lithium ion battery

1.1　模型描述

本文选取的电池为109 A·h的磷酸铁锂/石墨电池，其外形及几何尺寸如图2(a)和(b)所示。本文通过COMSOL Multiphysics中的PDE模块和固体传热模块的实时耦合，建立了耦合副反应产热的大型磷酸铁锂/石墨动力电池的三维烘箱热失控模型。模型中对电池划分网格的示意图见图2(c)。将109 A·h的电池电芯部分看作一个整体，基于集总模型假设，电池内部的能量守恒方程见式(1)^[6]

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 109 A·h磷酸铁锂电池的(a) 外形；(b) 几何及尺寸；(c) 网格

Fig.2 The (a) configuration; (b) geometry and dimension, as well as mesh of the 109 A·h lithium iron phosphate battery

ρ C_{p} \frac{\partial T}{\partial t} = λ \nabla^{2} T + q_{s i d e} / V_{c e l l} - q_{d i s s} / V_{b a t t}

(1)

式中， $ρ$ 和 $C_{p}$ 为电池的密度和比热容； $T$ 为温度； $t$ 为时间； $λ$ 为电池各个方向的导热系数； $q_{s i d e}$ 为电池热失控副反应产热； $q_{d i s s}$ 为散热项； $V_{c e l l}$ 为电芯的体积(不包括极耳和集流体)； $V_{b a t t}$ 为电池的体积。在热失控模型中不考虑电芯的可逆热和不可逆热。散热考虑热对流 $q_{c o n v}$ 和热辐射 $q_{r a d}$ ，见式(2)

q_{d i s s} = q_{c o n v} + q_{r a d} = A_{b a t t} h (T - T_{a m b}) + A_{b a t t} ε σ_{B} (T^{4} - T_{a m b}^{4})

(2)

式中， $A_{b a t t}$ 为电池表面与空气的换热面积； $h$ 和 $ε$ 分别为对流换热系数和辐射率； $T$ 和 $T_{a m b}$ 分别为电池表面温度和环境温度； $σ_{B}$ 为Stefan-Boltzmann常数。

1.2　热失控副反应

热失控副反应产热 $q_{s i d e}$ 包含四部分，即SEI膜分解放热 $q_{S E I}$ ，负极与电解液反应放热 $q_{a n}$ ，正极与电解液反应放热 $q_{c a}$ ，电解液分解热 $q_{e l e}$ ，满足式(3)

q_{s i d e} = q_{S E I} + q_{a n} + q_{c a} + q_{e l e}

(3)

1.2.1　SEI膜分解

SEI膜分解放热量可通过式(4)计算^[7]，其反应物浓度与分解速率的关系通过式(5)和(6)得出

q_{S E I} = H_{S E I} m_{S E I} R_{S E I}

(4)

R_{S E I} = A_{S E I} e x p (- \frac{E_{a, S E I}}{R T}) c_{S E I}

(5)

\frac{d c_{S E I}}{d t} = - R_{S E I}

(6)

式中， $H_{S E I}$ 为SEI膜的化学反应生成焓； $m_{S E I}$ 为SEI膜质量； $R_{S E I}$ 为SEI膜分解反应速率； $c_{S E I}$ 为SEI的归一化浓度； $A_{S E I}$ 和 $E_{a, S E I}$ 分别为指前因子和活化能。

1.2.2　负极与电解液反应

负极中的嵌入锂与电解液分解产生的热量可由式(7)计算^[7]，而负极与电解液分解反应速率除了与负极嵌入锂浓度 $c_{a n}$ 有关外，还会因发生SEI膜的重生反应而与SEI膜厚度 $δ_{S E I}$ 有关，见式(8)～(10)。

q_{a n} = H_{a n} m_{a n} R_{a n}

(7)

R_{a n} = A_{a n} e x p (- \frac{E_{a, a n}}{R T}) c_{a n} e x p (- \frac{δ_{S E I}}{δ_{S E I, r e f}})

(8)

\frac{d δ_{S E I}}{d t} = R_{a n}

(9)

\frac{d c_{a n}}{d t} = - R_{a n}

(10)

式中， $δ_{S E I, r e f}$ 为参考值，其值为1。

1.2.3　正极与电解液反应

正极与电解液反应放热量见式(11)，反应速率与正极转化率 $c_{c a}$ 之间的关系可由式(12)和(13)确定^[7]

q_{c a} = H_{c a} m_{c a} R_{c a}

(11)

R_{c a} = A_{c a} e x p (- \frac{E_{a, c a}}{R T}) c_{c a} (1 - c_{c a})

(12)

\frac{d c_{c a}}{d t} = R_{c a}

(13)

1.2.4　电解液分解反应

电解液自身分解放热量见式(14)，反应速率与电解液归一化浓度 $c_{e l e}$ 之间的关系由式(15)和(16)确定

q_{e l e} = H_{e l e} m_{e l e} R_{e l e}

(14)

R_{e l e} = A_{e l e} e x p (- \frac{E_{a, e l e}}{R T}) c_{e l e}

(15)

\frac{d c_{e l e}}{d t} = - R_{e l e}

(16)

1.3　热失控模型参数

电池几何尺寸为实验测量，各副反应动力学相关参数以及材料热物性参数均取自文献，并在同一数量级下调整得出，具体模型参数详见表1。

表1 热失控模型参数

Table 1 The parameters of thermal runaway model

参数	符号(单位)	负极	SEI	正极	电解液	来源
活化能 /J·mol^-1	E_a	1.3508E5	1.3508E5	1.396E5	2.74E5	[7]
指前因子/s^-1	A	2.5E13	1.667E15	6.667E13	5.14E25	[7]
化学反应生成焓/J·g^-1	H	1714	257	314	155	[7]
各组分质量/g	m	102	102	240	112	a
归一化浓度初始值	c₀	0.75	0.15	0.04	1	[8]
密度/kg·m^-3	$ρ$	2600				[9]
比热容 /J·(kg·K)^-1	$C_{p}$	1100				[9]
导热系数 /W·(m·K)^-1	$λ$	$λ_{x} = λ_{z} = 21$ , $λ_{y} = 1.1$				[9]
对流换热系数 /W·(m^-2·K^-1)	$h$	8.7				[10]
辐射率	$ε$	0.8				[10]

注：a表示没有具体的参考文献，根据文献取值在同一范围内调整得到。

新窗口打开| 下载CSV

1.4　热失控模型的验证

在155 ℃的烘箱模拟中，本研究模拟得到的电池热失控过程平均温度结果与文献^[10]中对磷酸铁锂电池的热失控模拟结果对比见图3。从图3中可以看出电池平均温度曲线吻合良好，热失控触发时间和热失控最高温度基本一致。不同之处在于热失控前的温升速率：文献中的电池初始温升速率较快，而本文模拟的电池初始温升速率较慢。这是由于文献中的磷酸铁锂电池为55 A·h，电池尺寸为62 mm×100 mm×50 mm(宽×高×厚)，其尺寸小于本文模拟的109 A·h电池，小尺寸的电池具有较大的比表面积，Kim等^[7]的研究表明，小尺寸电池由于其较大的体积比热交换面积，前期温升大于大尺寸电池。而在热失控后期由于大尺寸电池具有较大的散热面积，因此温度下降较小尺寸电池快。综上，认为本文模型具有较高的精确度，可保证文中模拟研究结果的可信性。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 本文温度模拟结果与文献数据^[10]对比

Fig.3 The temperature comparison between simulation in this paper and in the reference^[10]

2 不同烘箱温度下磷酸铁锂电池热失控机制

模拟磷酸铁锂电池在140、145、150、155、160和165 ℃六个烘箱温度下的热失控机制，包括烘箱加热温度对电池温度分布、各副反应浓度的影响。

2.1　不同烘箱温度下磷酸铁锂电池热失控过程

2.1.1　未发生热失控过程

图4为六个烘箱温度下磷酸铁锂电池最高温度随时间的变化曲线。在140 ℃时随着加热时间的延长，电池温度逐渐升高，并在后期趋于平稳；而145 ℃工况下电池温度在5160 s左右出现一个微小的峰值(约160 ℃)。二者的温度变化规律可由图5中SEI膜和负极分解情况来说明：在两个工况下，SEI膜和负极均未反应完全，且并未激发正极与电解液的反应，SEI和负极的分解反应产生的热量不足以触发电池热失控；而与140 ℃相比，145 ℃时SEI和负极分解速率较快，且反应物浓度消耗较大，因此产生了一个温度峰值。综合上述分析，可以看出在烘箱温度为140 ℃和145 ℃时，电池并未发生热失控。而与黄文才等^[4]对NCM电池的高温热失控研究对比可知，NCM电池在145 ℃加热时已发生热失控，且热失控最高温度达到近200 ℃，也表明磷酸铁锂电池的热稳定性高于三元电池。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 不同烘箱温度下电池平均温度曲线

Fig.4 The battery temperature profiles at different oven temperature

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 140 ℃和145 ℃下SEI(a)和负极(b)归一化浓度的变化

Fig.5 The normalized concentration of SEI (a) and anode (b) at 140 ℃ and 145 ℃

2.1.2　发生热失控过程

从图4看出，其他4个工况(150 ℃、155 ℃、160 ℃、165 ℃)下电池均发生热失控，图中标注为该4个工况下电池的最高温度。表2给出了不同烘箱温度下电池的最高温度及对应时间。随着烘箱温度的提高，电池发生热失控的时间提前，温升速率加快。但是从图4和表2中还可以发现，当烘箱温度从150 ℃增加到155 ℃时，电池最高温度增加了21 ℃，而当烘箱温度从160 ℃增加到165 ℃时，电池最高温度只增加了7 ℃，且热失控开始时间也提前较少。这是由于电池体系中各反应产物浓度有限，在反应物浓度消耗殆尽时，电池温度便不会继续升高，因此随着烘箱温度进一步的提升，电池热失控最高温度不会无限升高。与Kim等对LCO体系电池的高温热失控行为对比可知，LCO电池在150 ℃下热失控最高温度达到近350 ℃，而本文对磷酸铁锂电池的热失控模拟结果表明最高温度仅有208 ℃，进一步说明了磷酸铁锂电池的高安全性和热稳定性。

表2 不同烘箱温度下电池最高温度及对应时间

Table 2 The maximum temperature and its corresponding time at different oven temperature

烘箱温度/℃	140	145	150	155	160	165
电池最高温度/℃	150	162	208	229	240	247
电池最高温度对应时间/s	6095	5175	4635	3610	3030	2655

新窗口打开| 下载CSV

为了进一步说明电池热失控进程，图6代表性地给出了在150 ℃和160 ℃两个烘箱温度下，磷酸铁锂电池的平均温度，以及各副反应的归一化浓度。两条点划线分别代表两种工况下电池最高温度对应的时间。从图6中可以看出两个温度下的电池均发生热失控，我们可以从图中得出各温度下的热失控反应历程，以160 ℃烘箱温度为例进行说明：在大约110 ℃(1000 s)，SEI膜开始分解，但此时分解速度较慢；之后在约115 ℃(1140 s)开始，负极嵌入锂在失去SEI膜的保护后也开始分解，两个反应相互促进，造成了在150 ℃(2300 s)左右SEI膜和负极的快速分解，这也正是热失控起始温度点。SEI膜首先分解完成，而后负极与电解液反应结束；此时二者的快速消耗也加剧了正极与电解液的反应，见图6(e)。对比图6(a)和(e)可以发现电池热失控最高温度出现的时间点恰好是正极快速分解的时间点，且正极反应迅速，在290 s内全部反应完毕，可见正极与电解液的反应是造成热失控的主要导火索；而后紧接着电解液的自身分解也进一步加剧了热失控的发生，但是电解液归一化浓度在反应结束只减少到0.8，说明电解液分解并不完全。对比150 ℃和160 ℃两个烘箱温度下磷酸铁锂电池的热失控反应进程，我们发现在较高的烘箱温度下，各副反应分解提前，且反应速率较快。并且对比图6(f)中电解液分解反应历程可以发现在150 ℃烘箱温度下电解液只反应了很少的一部分，这是由于电池最高温度仅有204 ℃，刚刚达到电解液分解温度(200 ℃)，因此电解液一开始分解就由于较低的温度而停止。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 150 ℃和160 ℃烘箱温度下模拟结果：(a) 电池平均温度；(b) SEI归一化浓度， $c_{S E I}$ ；(c) SEI厚度， $t_{S E I}$ ；(d) 负极归一化浓度， $c_{a n}$ ；(e) 正极转化率， $c_{c a}$ ；(f) 电解液归一化浓度， $c_{e l e}$

Fig.6 Simulation results at 150 ℃ and 160 ℃ oven temperature: (a) the average temperature of the battery; (b) normalized concentration of SEI, $c_{S E I}$ ; (c) normalized thickness of SEI, $t_{S E I}$ ; (d) normalized concentration of anode, $c_{a n}$ ; (e) degree of conservation of cathode, $c_{c a}$ ; (f) normalized concentration of electrolyte, $c_{e l e}$

2.2　不同烘箱温度下磷酸铁锂电池温度分布

2.2.1　未发生热失控过程

为了明确不同烘箱温度下磷酸铁锂电池的温度分布，图7给出了未发生热失控情况下(140 ℃和145 ℃烘箱温度)电池在不同时刻的三维温度分布图。分别选取了加热1000、3000、5000 s，电池达到最高温度以及加热结束时的温度分布为代表进行分析。可以看出未发生热失控时两个烘箱温度条件下电池的温度分布一致且全程呈均匀分布，即在加热前期，电池温度低于烘箱环境温度，使得电池表面温度高于内部温度；而在加热至5000 s时，电池温度超过烘箱环境温度，电池内部热量的积累使得内部温度高于表面温度，直至加热结束一直维持同样的温度分布，这是由于电池的温度依然未降至烘箱温度之下。由此可见，在未发生热失控的情况下，电池温度分布一直比较均匀。

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7 烘箱温度140 ℃和145 ℃，电池未发生热失控情况下不同时刻的温度分布

Fig.7 The battery temperature distribution at different characteristic time at 140 ℃ and 145 ℃ with no thermal runaway cases

2.2.2　发生热失控过程

图8为发生热失控的情况下150 ℃和160 ℃两个工况下磷酸铁锂电池的温度分布图，选择了6个时间点，分别为热失控前20 s，前10 s，热失控时，热失控后10 s，后20 s，以及加热结束(7200 s)。可以看出在150 ℃烘箱加热条件下，电池在热失控前后40 s内的温度分布基本一致，且温度变化较小。电池内部各副反应加速分解导致的热量积累，使得电池内部温度升高而电池表面温度低于内部温度，且温度分布不均匀。直到加热结束，电池温度降到152 ℃，基本与环境温度一致，电池温度又恢复均匀分布。160 ℃工况下电池的温度分布变化更快且不均匀程度增大：在热失控前20 s，热量已聚集在电池内部，内部温度开始升高；在热失控前10 s，副反应的继续进行使得电池温度持续升高，且较高的环境温度使得电池内部的热量无法散失并且表现出在电池表面的温度升高；直至发生热失控，电池温度达到峰值，且温度分布极度不均匀，意味着电池内已发生不可逆的损失；而在热失控过后的20 s内电池的温度虽略有降低，但是内部依然维持高温。与150 ℃工况不同的是，在加热结束后160 ℃工况下电池表面温度高于内部温度，这是由于此时环境温度(160 ℃)高于电池温度(159 ℃)，因此电池表面温度高于内部温度，但是依然呈现出均匀的温度分布。由此可以看出，高环境温度下电池的热失控更为剧烈，温度分布极不均匀，且在热失控前后温度分布变化较快，表明电池材料已发生不可逆的分解导致电池的损坏。且通过与图7对比可以看出，在已发生热失控的情况下，电池的温度分布极不均匀，也表明热失控时电池内部已经损坏。

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8 烘箱温度150 ℃和160 ℃，电池发生热失控情况下不同时刻的温度分布

Fig.8 The battery temperature distribution at different characteristic time at 150 ℃ and 160 ℃ with thermal runaway cases

3 结论

本文针对109 A·h大型磷酸铁锂电池建立了三维烘箱热失控模型，并与文献数据对比验证了模型的精确性。进一步模拟研究了在不同烘箱温度下电池的热失控行为特征及其内在机制。通过对比发生热失控和未发生热失控情况下各副反应分解浓度和电池的温度分布可以得出：环境温度在140 ℃和145 ℃时，电池未发生热失控，仅发生了SEI膜和负极的分解反应，且温度分布比较均匀；在环境温度为150 ℃及以上工况下，电池均发生热失控，这一过程中正极与电解液的反应是热失控的主要致因。此外，高环境温度下电池的热失控表现更为剧烈，温度分布也极不均匀，且在热失控前后温度分布变化较快，预计电池材料发生的不可逆分解反应是导致电池损坏的主要原因。本文建立的大型磷酸铁锂电池热失控模型，可以直观的揭示磷酸铁锂电池高温加热条件的热失控微观机理，对后续其他滥用条件触发的热失控研究具有指导意义。

符号说明

$A$	指前因子，s^-1
$A b a t t$	电池表面与空气的换热面积，m²
$C p$	比热容，J/(kg·K)
$c$	归一化浓度
$E a$	活化能，J/mol
$H$	化学反应生成焓，J/g
$h$	对流换热系数，W/(m²·K)
$m$	各反应物质量，g
$q$	反应物分解放热量，W
$q c o n v$	对流散热量，W
$q d i s s$	总散热量，W
$q r a d$	辐射散热量，W
$q s i d e$	热失控副反应总产热，W
$R$	化学反应速率，s^-1
$T$	电池温度，℃
$T a m b$	环境温度，℃
$t$	时间，s
$V b a t t$	电池总体积，m³
$V c e l l$	电芯总体积(电池不包括极耳的体积)，m³
$δ S E I$	SEI膜厚度，m
$ε$	辐射率
$λ x, z$	平行于极片方向的导热系数，W/(m·K)
$λ y$	垂直于极片方向的导热系数，W/(m·K)
$ρ$	电池密度，kg/m³
$σ B$	Stefan-Boltzmann常数，W/(m²·K⁴)
下标
an	负极
ca	正极
ele	电解液
SEI	SEI膜
0	初始值

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

羡学磊, 董海斌, 张少禹, 等. 三元锂离子动力电池热失控及火灾特性研究[J]. 储能科学与技术, 2020, 9(1): 239-248.

[本文引用: 1]

XIAN Xuelei, DONG Haibin, ZHANG Shaoyu, et al. Thermal runaway and fire characteristics of NCM lithium-ion power battery[J]. Energy Storage Science and Technology, 2020, 9(1): 239-248.

[本文引用: 1]

[2]

WANG Q, MAO B, STOLIAROV S I, et al. A review of lithium ion battery failure mechanisms and fire prevention strategies[J]. Progress in Energy and Combustion Science, 2019, 73: 95-131.

[本文引用: 1]

[3]

杜光超, 郑莉莉, 张志超, 等. 圆柱形高镍三元锂离子电池高温热失控实验研究[J]. 储能科学与技术, 2020, 9(1): 249-256.

[本文引用: 1]

DU Guangchao, ZHENG Lili, ZHANG Zhichao, et al. Experimental study on high temperature thermal runaway of cylindrical high nickel ternary lithium-ion batteries[J]. Energy Storage Science and Technology, 2020, 9(1): 249-256.

[本文引用: 1]

[4]

黄文才, 胡广地, 邓宇翔, 等. 锂离子电池的高温热失控模拟[J]. 电池, 2019, 49(3): 204-207.

[本文引用: 2]

HUANG Wencai, HU Guangdi, DENG Yuxiang, et al. Analysis on high temperature thermal runaway simulation of Li-ion battery[J]. Battery Bimonthly, 2019, 49(3): 204-207.

[本文引用: 2]

[5]

RICHARD M, DAHN J. Accelerating rate calorimetry study on the thermal stability of lithium intercalated graphite in electrolyte. I. Experimental[J]. Journal of the Electrochemical Society, 1999, 146 (6): 2068-2077.

[本文引用: 1]

[6]

MEI W, CHEN H, SUN J, et al. The effect of electrode design parameters on battery performance and optimization of electrode thickness based on the electrochemical-thermal coupling model[J]. Sustainable Energy & Fuels. 2019, 3: 148-165.

[本文引用: 1]

[7]

KIM G H, PESARAN A, SPOTNITZ R. A three-dimensional thermal abuse model for lithium-ion cells[J]. Journal of Power Sources, 2007, 170 (2): 476-489.

[本文引用: 4]

[8]

HATCHARD T D, MACNEIL D D, BASU A, et al. Thermal model of cylindrical and prismatic lithium-ion cells[J]. Journal of the Electrochemical Society, 2001, 148 (7): A755-A761.

[9]

冯旭宁. 车用锂离子动力电池热失控诱发与扩展机理、建模与防控[D]. 北京: 清华大学, 2016.

FENG Xuning. Thermal runaway initiation and propagation of lithium-ion traction battery for electric vehicle: Test, modeling and prevention[D]. Beijing: Tsinghua University, 2016.

[10]

GUO G F, LONG B, CHENG B, et al. Three-dimensional thermal finite element modeling of lithium-ion battery in thermal abuse application[J]. Journal of Power Sources. 2010, 195 (8): 2393-2398.

[本文引用: 3]

... 随着锂离子电池尺寸、能量密度的不断增大，其安全性也越来越成为消费者关注的重点^[1].锂离子电池在受到电滥用、机械滥用以及热滥用的情况下会因内部热量的积累造成高温，引发热失控^[2]，进而造成火灾爆炸事故.可以看出高温是触发锂离子电池热失控的最直接原因，因此科学认识高温作用下诱发的锂离子电池热失控特征至关重要.目前对锂离子电池高温热失控的研究局限于三元电池，如杜光超等^[3]分析了圆柱形高镍三元锂离子电池在高温下的热失控参数，黄文才等^[4]针对三元电池建立了高温加热条件下的三维分层物理热失控模型，研究了不同传热系数下电池的热失控行为.相比于三元电池，目前关于磷酸铁锂电池在高温下的热失控研究还比较缺乏.特别是，对于大型磷酸铁锂电池的高温热失控行为特征、内在机制的认识更是严重不足.鉴于此，本文选取大型磷酸铁锂电池(109 A·h)为研究对象，对其建立不同加热温度下的烘箱热失控模型，通过对比发生热失控和未发生热失控两种条件，得出磷酸铁锂电池高温热失控行为. ...

... 图4为六个烘箱温度下磷酸铁锂电池最高温度随时间的变化曲线.在140 ℃时随着加热时间的延长，电池温度逐渐升高，并在后期趋于平稳；而145 ℃工况下电池温度在5160 s左右出现一个微小的峰值(约160 ℃).二者的温度变化规律可由图5中SEI膜和负极分解情况来说明：在两个工况下，SEI膜和负极均未反应完全，且并未激发正极与电解液的反应，SEI和负极的分解反应产生的热量不足以触发电池热失控；而与140 ℃相比，145 ℃时SEI和负极分解速率较快，且反应物浓度消耗较大，因此产生了一个温度峰值.综合上述分析，可以看出在烘箱温度为140 ℃和145 ℃时，电池并未发生热失控.而与黄文才等^[4]对NCM电池的高温热失控研究对比可知，NCM电池在145 ℃加热时已发生热失控，且热失控最高温度达到近200 ℃，也表明磷酸铁锂电池的热稳定性高于三元电池. ...

... 锂离子电池的烘箱热失控模型即模拟电池在烘箱中恒温加热的情况下的热失控发展规律.锂离子电池热失控各副反应阶段见图1.在90～120 ℃时^[5]，负极表面形成的SEI膜发生分解，当SEI膜分解完成后，负极失去保护，嵌入在负极的锂继而与电解液发生反应.随着温度的持续升高，隔膜开始熔化，且正极与电解液反应，且伴随着电解液的自身分解，这时会发生隔膜破裂、正负极短路，而导致热失控的发生.之后负极与黏结剂发生反应，引发电解液着火，造成电池的燃烧.在烘箱热失控模拟中，电池热失控的最高温度较绝热热失控模拟低，也不会发生燃烧现象，因此只考虑前4种副反应的发生. ...

... 本文选取的电池为109 A·h的磷酸铁锂/石墨电池，其外形及几何尺寸如图2(a)和(b)所示.本文通过COMSOL Multiphysics中的PDE模块和固体传热模块的实时耦合，建立了耦合副反应产热的大型磷酸铁锂/石墨动力电池的三维烘箱热失控模型.模型中对电池划分网格的示意图见图2(c).将109 A·h的电池电芯部分看作一个整体，基于集总模型假设，电池内部的能量守恒方程见式(1)^[6] ...

... SEI膜分解放热量可通过式(4)计算^[7]，其反应物浓度与分解速率的关系通过式(5)和(6)得出 ...

... 负极中的嵌入锂与电解液分解产生的热量可由式(7)计算^[7]，而负极与电解液分解反应速率除了与负极嵌入锂浓度

c_{a n}

有关外，还会因发生SEI膜的重生反应而与SEI膜厚度

δ_{S E I}

有关，见式(8)～(10). ...

... 正极与电解液反应放热量见式(11)，反应速率与正极转化率

c_{c a}

之间的关系可由式(12)和(13)确定^[7] ...

... 在155 ℃的烘箱模拟中，本研究模拟得到的电池热失控过程平均温度结果与文献^[10]中对磷酸铁锂电池的热失控模拟结果对比见图3.从图3中可以看出电池平均温度曲线吻合良好，热失控触发时间和热失控最高温度基本一致.不同之处在于热失控前的温升速率：文献中的电池初始温升速率较快，而本文模拟的电池初始温升速率较慢.这是由于文献中的磷酸铁锂电池为55 A·h，电池尺寸为62 mm×100 mm×50 mm(宽×高×厚)，其尺寸小于本文模拟的109 A·h电池，小尺寸的电池具有较大的比表面积，Kim等^[7]的研究表明，小尺寸电池由于其较大的体积比热交换面积，前期温升大于大尺寸电池.而在热失控后期由于大尺寸电池具有较大的散热面积，因此温度下降较小尺寸电池快.综上，认为本文模型具有较高的精确度，可保证文中模拟研究结果的可信性. ...

... [10]对比The temperature comparison between simulation in this paper and in the reference<sup>[<xref ref-type="bibr" rid="R10">10</xref>]</sup>

Fig.3

2 不同烘箱温度下磷酸铁锂电池热失控机制

模拟磷酸铁锂电池在140、145、150、155、160和165 ℃六个烘箱温度下的热失控机制，包括烘箱加热温度对电池温度分布、各副反应浓度的影响. ...

... [10]Fig.32 不同烘箱温度下磷酸铁锂电池热失控机制

模拟磷酸铁锂电池在140、145、150、155、160和165 ℃六个烘箱温度下的热失控机制，包括烘箱加热温度对电池温度分布、各副反应浓度的影响. ...

〈

〉