基于三维电化学热耦合析锂模型的锂离子电池参数设计

doi:10.19799/j.cnki.2095-4239.2022.0225

基于三维电化学热耦合析锂模型的锂离子电池参数设计

马勇^,¹, 李晓涵¹, 孙磊¹, 郭东亮¹, 杨景刚¹, 刘建军¹, 肖鹏¹, 钱广俊^,²

1.国网江苏省电力有限公司电力科学研究院，江苏南京 211103

2.清华大学汽车安全与节能国家重点实验室，北京 100084

Parameter design of lithium-ion batteries based on a three-dimensional electrochemical thermal coupling lithium precipitation model

MA Yong^,¹, LI Xiaohan¹, SUN Lei¹, GUO Dongliang¹, YANG Jinggang¹, LIU Jianjun¹, XIAO Peng¹, QIAN Guangjun^,²

1.State Grid Jiangsu Electric Power Co. , Ltd. Research Institute, Nanjing 211103, Jiangsu, China

2.State Key Laboratory of Automotive Safety and Energy, Tsinghua University, Beijing 100084, China

通讯作者: 钱广俊，博士，研究方向为锂离子电池组均衡算法和电化学阻抗谱模型，E-mail：qguangjun@163.com。

收稿日期: 2022-04-25 修回日期: 2022-05-14

基金资助:

国网江苏省电力有限公司“电动汽车动力电池状态检测与安全评估技术研究”项目. J2021042

Received: 2022-04-25 Revised: 2022-05-14

作者简介 About authors

马勇（1986—），男，硕士，研究方向为电力系统输变电设备状态评估技术，E-mail：ma.y@foxmail.com； E-mail：ma.y@foxmail.com。

摘要

锂离子电池负极析锂可能会诱发热失控，进而导致安全事故。而通过优化电池设计参数能够有效减少析锂副反应的发生，因此本工作提出一种基于三维电化学热耦合析锂模型的锂离子电池参数设计优化方法。首先，将模型参数进行分类，分别采用实验、精确测量、文献查找和参数辨识等方法获取相应的参数。同时加入可逆锂重嵌入机制和产热模型，建立三维电化学热耦合析锂模型。模型建立完成后，对模型精度进行验证，验证结果表明模型可以较好地模拟电池在常温和低温下端电压的变化，并且能够定量描述在低温大倍率充电期间电池内部的析锂程度、温度分布等非均一现象。最后，通过分析电极尺寸和极耳位置，研究电池设计参数对非均一析锂的影响。仿真结果表明：电极长度增加会导致电极区域温度差异和电流密度的不一致性增大，综合影响下使电池析锂时间略有提前，但对电池总体析锂程度影响较小；电池极耳位置处于长度方向的轴线对侧时能够有效缓解负极析锂，相对析锂程度降低了16.7%。

关键词： 锂离子电池 ; 三维电化学热耦合 ; 重嵌入机制 ; 析锂模型

Abstract

Lithium deposition on the negative electrode of lithium-ion batteries may induce thermal runaway that can lead to safety accidents. The occurrence of lithium precipitation side reactions can be reduced effectively by optimizing the battery design parameters. Therefore, this study proposes a parameter design optimization method for lithium-ion batteries based on a three-dimensional electrochemical, thermal coupled lithium precipitation model. First, the model parameters were classified, and the corresponding parameters were obtained using experiments, exact measurements, literature searches, and parameter identification, respectively. The reversible lithium re-embedding mechanism and heat production model were added to establish a three-dimensional electrochemical, thermal coupling lithium precipitation model. After constructing the model, the accuracy of the model was verified. The verification results showed that the model could simulate the changes in the terminal voltage of the battery at room temperature and low temperature and could quantitatively describe the nonhomogeneous phenomena, such as the degree of lithium precipitation and temperature distribution inside the battery during low temperature large rate charging. Finally, the effects of the battery design parameters on nonhomogeneous lithium precipitation were investigated by analyzing the electrode size and lug position. The simulation results showed that increasing the electrode length would increase the temperature difference in the electrode area and the inconsistency of the current density. The combined effect would advance the lithium precipitation time of the battery slightly, but the effect on the overall degree of lithium precipitation of the battery was relatively small. When the position of the battery tab was on the opposite side of the axis in the longitudinal direction, it could effectively alleviate lithium deposition on the negative electrode, and the relative lithium precipitation degree was reduced by 16.7%.

Keywords： lithium-ion battery ; three-dimensional electrochemical thermal coupling ; re-embedding mechanism ; lithium precipitation model

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

本文引用格式

马勇, 李晓涵, 孙磊, 郭东亮, 杨景刚, 刘建军, 肖鹏, 钱广俊. 基于三维电化学热耦合析锂模型的锂离子电池参数设计[J]. 储能科学与技术, 2022, 11(8): 2600-2611

MA Yong. Parameter design of lithium-ion batteries based on a three-dimensional electrochemical thermal coupling lithium precipitation model[J]. Energy Storage Science and Technology, 2022, 11(8): 2600-2611

在当前的快充场景中，电池有可能会出现负极析锂现象^[1-2]。这种现象会导致电池管理系统对电池单体监测异常^[3]，从而造成电池过热和热稳定性降低等问题，甚至可能引发热失控^[4-5]。因此，电池的充电安全问题亟须解决^[6-7]。

抑制电池充电时发生析锂可以从电池设计方面入手^[8-9]。Sturm等^[10]基于多维结构模型，以负极过电位和温度评估了极片和极耳的设计对负极析锂的影响，并通过设定过电位阈值将高倍率充电模式下的析锂风险降到最低。Zhao等^[11]建立了三维模型，通过多极耳的设计有效降低了长电极的集流体带来的电子传递损失，提高了能量密度。Samba等^[12]对软包电池进行三维建模，从电流密度分布、电势、温度分布等角度分析了极耳位置对电池带来的影响。其研究结果表明，对称结构的设计使得电池性能更加均匀，电流密度梯度更小，能够使活性材料最大程度得到利用。然而，常规的电化学热耦合模型没有考虑负极发生析锂后的可逆锂重嵌入机制，无法很好地描述电池在充电过程中实际的析锂状况。因此，添加可逆锂重嵌入机制的多维结构电化学热耦合析锂模型将能够更好地分析电池内部非均一析锂问题^[13]。

本工作提出一种基于三维电化学热耦合析锂模型的锂离子电池参数设计优化方法。首先，分析各种参数的获取方式，并加入可逆锂重嵌入机制，建立三维电化学热耦合析锂模型。之后，对模型精度进行验证。最后，通过分析电极尺寸和极耳位置，研究结构设计对非均一析锂的影响。

1 三维电化学热耦合析锂模型

1.1　模型建立

本工作所使用的模型是在Newman和Doyle等^[14-15]建立的伪二维(pseudo two dimensional，P2D)机理模型基础上进一步建立完成的。

电池在充放电过程中会产热，造成电池温度上升，从而影响电池内部粒子输运特性及反应动力学参数。因此，需要考虑电池产热对其内部反应状态的影响。根据能量守恒定律，电池内部能量满足式(1)，即：

ρ C_{p} \frac{\partial T}{\partial t} = λ_{x} \frac{\partial^{2} T}{\partial x^{2}} + λ_{y} \frac{\partial^{2} T}{\partial y^{2}} + λ_{z} \frac{\partial^{2} T}{\partial z^{2}} + q

(1)

其中，ρ为电池密度；C_p 为电池比热容；λ为不同方向的电池导热系数；q为总生热率。边界条件如式(2)所示：

- λ (\frac{\partial T}{\partial n}) = h (T_{a m b} - T_{\infty})

(2)

其中，n为电池表面法线矢量方向；h为对流换热系数；T_amb为电池所处环境温度；T_∞为电池温度。

本工作参数获取方式主要为拆解测量、实验标定、参数辨识、文献参考等，模型具体参数和分类如表1所示。模型中设计参数部分，如极片厚度L、集流体厚度L_cc、正负极厚度L_p/n、隔膜厚度L_sep、活性物质总反应面积A_total等参数可通过工具测量获得。而正负极均衡电势曲线U_{ref, p/n}则通过制作纽扣电池，以极小电流的充放电实验来获取。固相颗粒最大嵌锂浓度c_{s, max}，初始时刻化学计量比x_0/100、y_0/100等一般通过参数辨识进行获取。正负极固相颗粒半径r_{s, n/p}、电解质初始盐浓度cl₀、孔隙率ε等由电池厂商测定后提供。模型中的热物性参数，如比热容C_p 、对流换热系数h等通过实验进行测定。液相参数中的液相扩散系数D_e、液相离子电导率与温度T及电解质浓度c有关，可以直接引用文献中的参考值^[16]。正负极材料的熵热系数通过对COMSOL案例库中已录入的数据进行适当调整即可得到有效的dU/dT曲线。由于固相扩散系数D_s、反应速率常数k与温度T密切相关且符合Arrhenius公式，可以通过不同温度下的脉冲实验标定得到相应数值和活化能E_a。

表1 电化学热耦合模型参数

Table 1 Electrochemical thermal coupling model parameters

参数	参数意义	单位	负极	隔膜	正极
L	极片厚度	$μ m$	100^①	15^①	51^①
L_cc	集流体厚度	$μ m$	15.3^①	—	9^①
i_1C	1 C倍率下平均电流密度	A/m²	26.67^①
c_s,max	固相颗粒最大嵌锂浓度	mol/m³	25255^②	—	110995^②
x_0/100	负极初始时刻化学计量比	1	0.0073/0.9501^②
y_0/100	正极初始时刻化学计量比	1	0.9329/0.2904^②
A_total	活性物质总反应面积	m²	5.49^①
r_s	固相颗粒半径	$μ m$	12^③	—	10^③
κ_e	液相离子电导率	S/m	$\begin{array}{l} κ_{e} = 1 \times 10^{- 4} c (- 10.5 + 0.074 T - 3.96 \times 10^{- 5} T^{2} + 6.68 \times 10^{- 4} c - \\ 1.78 \times 10^{- 5} T + 2.8 \times 10^{- 8} T^{2} + {4.94 \times 10^{- 7} c^{2} - 8.86 \times 10^{- 10} c^{2} T)}^{2} \end{array}$
σ_s	固相电导率	S/m	100^④	—	10^④
ε_e	液相体积分数	1	0.357^③	0.4^③	0.444^③
σ_s^eff	固相有效电导率	S/m	ε_s^1.5σ_s	—	ε_s^1.5σ_s
cl₀	电解质初始盐浓度	mol/m³	1000^③
E_a	活化能	J/mol	21350^⑤	—	20500^⑤
k_ref	参考温度下的反应速率常数	m/s	3e^-11④,⑤	—	4e^-11④,⑤
D_e	液相扩散系数	m²/s		$D_{e} = 1 \times 10^{- 8.83 - \frac{54}{(T - 5 \times 10^{- 3} c - 299)} - 2.2 \times 10^{- 4} c}$
D_s,ref	参考温度下的固相扩散系数	m²/s	3.9e^-14④,⑤	—	4.5e^-13④,⑤
T_ref	参考温度	K		298.15
α_a、α_c	传递系数	1	0.5	—	0.5
C_p	比热容	J/(kg·K)		1150^①
h	对流换热系数	W/(m²·K)		22^①
F	法拉第常数	C/mol		96485
R	理想气体常数	J/(mol·K)		8.314

注：①为实际测量值，②为参数辨识数据，③为电池厂商提供数据，④为文献参考值^[17-18]，⑤为标定数据。

新窗口打开| 下载CSV

之后，对某软包电池进行三维结构建模，电池标称容量为24 Ah，工作电压区间为2.5～4.3 V，工作温度区间为-20～55 ℃。图1左侧为三维模型几何结构，模型考虑了x、y方向的极片尺寸，并引入了集流体和极耳的几何结构。三维模型中电化学部分采用P2D模型理论，模型中正负极由固相颗粒和电解液组成固液两相，隔膜区域充满电解液。其中，正负极集流体箔材属性为Al和Cu，仅起到电子载体作用。

图1

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

图1 锂离子电池三维模型结构示意图

Fig. 1 Schematic diagram of the three-dimensional model structure of the lithium-ion battery

同时考虑三维模型中集流体长度方向上的电势分布，满足式(3)

\frac{\partial}{\partial x} ϕ = i_{a p p} ρ / A_{c c}

(3)

其中，ρ为集流体的电阻率；A_cc为集流体沿长度方向横截面积。三维传热模型中，电流密度在空间上存在分布差异，造成局部产热功率q不同。相应区域的热物性参数如表2所示。

表2 三维几何模型热物性参数

Table 2 Thermal physical parameters of the three-dimensional geometric model

参数	单位	正极		隔膜	负极
参数	单位	NCM	集流体	隔膜	石墨	集流体
比热容C_p	J/(kg·K)	1100	900	1978	881	285
导热系数λ	W/(m·K)	1.58	238	0.344	1.04	400
密度ρ	kg/m³	2628	2700	1108	2270	8960
对流换热系h	W/(m²·K)	大面：1.1 极耳尾部：10

新窗口打开| 下载CSV

三维模型中局部活性区域的反应速率k、扩散系数D_s、D_e和电导率κ等参数受局部温度影响，利用Arrhenius公式进行控制。

进一步，在该三维电化学热耦合模型中对多孔电极负极添加析锂副反应方程建立析锂模型。在只含有锂析出副反应的模型基础上考虑了可逆锂重嵌入机制，搭建更加符合实际的析锂模型。除了正常嵌锂反应动力学理论外，析锂反应速率也满足Bulter-Volmer方程，如式(4)所示：

j_{p l} = j_{p l, 0} [e x p (\frac{α_{a, p l}}{R T} η_{L i}) - e x p (\frac{α_{c, p l}}{R T} η_{L i})]

(4)

其中，j_{pl, 0}为析锂交换电流密度；α_{a, pl}和α_{c, pl}为电化学传递系数，分别取值为0.3、0.7；η_Li为析锂反应过电势，当η_Li<0时，析锂副反应触发。析锂反应交换电流密度j_{pl, 0}如式(5)所示：

j_{p l, 0} = F k_{p l} {(\frac{c_{e}}{c_{e, r e f}})}^{α_{a, p l}}

(5)

其中，k_pl为锂沉积反应速率常数。析锂过电势η_Li满足式(6)：

η_{L i} = ϕ_{s} - ϕ_{e} - U_{e, L i} - Δ ϕ_{f i l m, L i}

(6)

其中，U_{e, Li}为析锂反应平衡电势，一般取值0 V； $Δ ϕ_{f i l m, L i}$ 为固体电解质界面(solid electrolyte interface，SEI)膜和析锂产物内层与外层之间的电势差V。锂析出后，一旦负极η_Li＞0时，可逆锂会重新溶解，其动力学反应如式(7)所示：

j_{s t} = j_{s t, 0} [e x p (\frac{α_{a, s t}}{R T} η_{L i}) - e x p (\frac{α_{c, s t}}{R T} η_{L i})]

(7)

其中，j_{st, 0}为锂溶解交换电流密度。

SEI膜阻和析锂产物造成电势差满足Ohm定律，即：

Δ ϕ_{f i l m, p l} = j_{p l} R_{f i l m}

(8)

其中，在析锂期间，R_film随析锂程度的增加而增厚，增厚部分主要因析出的锂金属覆盖在固相颗粒上。另外，析出的金属锂与电解液发生反应，产物主要为Li₂CO₃、LiF等，造成SEI膜增厚。因而，膜阻R_film可表示为初始SEI膜阻R₀和新增SEI膜阻△R_film之和，即：

R_{f i l m} = R_{0} + Δ R_{f i l m}

(9)

其中，膜阻R_film表达式可由电导率σ和δ之间关系计算，即：

R_{f i l m} = \frac{δ_{0} + Δ δ}{σ}

(10)

其中，δ、σ分别为SEI膜厚度和SEI膜电导率。

新增SEI膜厚度△δ变化率与析锂反应电流密度之间关系满足下式，即：

\frac{d δ}{d t} = - \frac{M j_{p l}}{F ρ}

(11)

其中，M为金属锂摩尔质量，ρ为锂金属密度。

表3为文中析锂模型中金属锂沉积及可逆锂重嵌入反应动力学重要参数。部分参数来源于文献，没有引用的参数在实际范围内进行评估，以符合实验数据。

表3 析锂与可逆锂重嵌入动力学反应相关参数

Table 3 Parameters related to the kinetic reactions of analytical and reversible lithium re-embedding

参数	物理意义	单位	负极
α_{a, pl}、α_c,_pl	析锂/锂溶解传递系数	1	0.3、0.7^[19]
k_pl/st	锂析出/溶解反应速率常数	m/s	6.5×10^-6
E_{a, pl/st}	析锂与锂溶解反应活化能	J/mol	50000
z₁、z₂、z₃	可逆锂/不可逆锂/SEI占比	1	0.57、0.375、0.05
σ_SEI	SEI膜电导率	S/m	5×10^-6[19]
δ_SEI	SEI膜t厚度	m	式(10)
ρ_SEI	SEI 膜密度	kg/m³	1690^[19]
ρ_Li	锂金属密度	kg/m³	534^[19]
M_SEI	SEI膜摩尔质量	kg/mol	0.162^[19]
M_Li	锂金属摩尔质量	kg/mol	6.94×10^-3[19]

新窗口打开| 下载CSV

在COMSOL Multiphysics 5.4软件平台中基于上述控制方程对锂离子电池建立三维电化学热耦合析锂模型。

1.2　模型验证

模型建立完成后，通过常温25 ℃下不同倍率恒流充电仿真与实验结果对比，验证模型准确性。其中，25 ℃下选用1/3 C、1/2 C、1 C、5/2 C倍率，从0% SOC状态下充电至截止电压4.2 V。仿真结果如图2(a)所示，图中为25 ℃下不同倍率充电电压曲线(带标记为实验数据，实线为仿真数据)。对比结果表明，常温下模型仿真数据和实验数据较吻合，且端电压的RMSE小于10 mV，表明该三维电化学热耦合模型能够较好地模拟电池充电过程中的端电压变化，模型精度较高。

图2

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

图2 25 ℃不同倍率充电时端电压变化和电池采样位置

Fig. 2 Variation of terminal voltage and battery sampling position during charging at different multipliers at 25 ℃

进一步，利用三维模型模拟电极局部区域温度差异，选取电极区域5个测温点，采样点位置如图2(b)所示。图3(a)和3(b)分别为常温1 C和5/2 C恒流充电过程中电极区域5个点温度变化情况。采样点温度表现为T2>T1>T3>T5>T4，其中1、2位置靠近极耳，该处电流密度较大，产热量大，因此温度要高于其他位置。图3(a)中1 C倍率下模型与实际采样点温度最大误差不超过0.5 ℃，图3(b)中5/2 C倍率的模型与实际采样点温度最大误差小于1.5 ℃，造成误差的主要原因来自两方面：一方面是热物性参数测试的精确性；另一方面由于实验时受到温箱内部气体的对流散热影响，导致实验温度值相对仿真值更低。但模型曲线整体上的趋势与实际一致，不会影响后续的分析结果。总体而言，1 C倍率下各区域温度差异不大，靠近负极极耳处的位置2与底部位置4之间最大温差仅为0.4 K，原因为该款电池为软包叠片型，外部有足够的散热面积用于散热。在大倍率5/2 C下，位置2与位置4之间最大温差也仅为2.3 K。

图3

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

图3 不同倍率充电时电极区域温度验证

Fig. 3 Temperature verification of electrode area at different multiplier charging

析锂现象主要发生于大倍率快充和低温工况下，模型中输运性能参数、反应动力学参数均对温度敏感，并符合Arrhenius公式。在低温0 ℃下仿真验证该模型对端电压及温度的响应能力。验证结果如图4(a)和4(b)所示，模型在0 ℃下不同倍率充电电压及位置3处的温度均能和实验数据较好地吻合。电压在低SOC区间误差较大，主要原因为电池动力学参数对温度更加敏感，但整体上0 ℃仿真电压曲线基本能描述电压的变化情况。实际测量中，热电偶采集温度为大面中心位置温度，接近模型仿真中测温点3处温度，因此，选取此处进行温度验证。结果显示，在1/2 C小倍率下温度较为接近，1 C及以上倍率存在一定误差，但温度最大误差小于1 ℃。低温验证结果表明，模型能较为准确地描述低温充电下电池端电压及温度变化，模型精度较高。

图4

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

图4 0 ℃不同倍率充电时模型验证

Fig. 4 Model validation at different multiplier charging at 0 ℃

同时，电池析锂动力学反应和局部温度密切相关。排除产热过于严重的极耳区域，可以用电极区域 $Δ T$ 变化来描述充电温度差异的变化规律， $Δ T$ 为局部温度与电极区域平均温度之差，即：

Δ T = T - T_{a v e r a g e}

(12)

如图5所示，在充电初期10 s时电极区域最大温差 $Δ T$ 仅为0.05 K，发生在极耳附近。随着充电的进行，在300 s时温度 $Δ T$ 增长到0.18 K，充电截止时刻仍达到0.15 K。低温下不同倍率仿真验证表明，三维电化学热耦合模型可以有效模拟电池在低温下电极局部温度间的差异。

图5

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

图5 0 ℃下1 C倍率充电期间电池温度分布差异变化情况

Fig. 5 Variation of battery temperature distribution during 1 C multiplier charging at 0 ℃

进一步，为分析电池局部析锂情况，图6显示了充电期间负极-隔膜处的固液相电势差分布图。达到析锂触发时间t_onset=205 s时，靠近极耳区域开始出现φ_s-φ_e<0 V，造成靠近极耳区域首先触发析锂副反应，其余区域暂未出现。

图6

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

图6 负极-隔膜处 φ_s-φ_e 分布

Fig. 6 φ_s-φ_e distribution at the negative electrode and diaphragm

在三维析锂模型理论中，金属锂析出后副反应产物包括3部分，为不可逆锂、可逆锂和用于SEI膜增厚的锂。图7为隔膜-负极边界处副反应产物浓度面内分布情况。在达到t_onset=205 s之前，析锂反应暂未触发，副反应产物浓度为0。t=300 s时，界面处副反应产物浓度已经呈现分布非均一现象，靠近极耳处浓度较大。随着时间的推移，界面处副反应产物浓度增加，截止时刻达到了1.09×10³ mol/m³。

图7

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

图7 隔膜与负极边界处副反应产物浓度变化情况

Fig. 7 Variation of by-product concentration at the boundary between the diaphragm and the negative electrode

对极片厚度z轴方向分析，如图8所示。负极厚度方向析锂程度存在差异，在充电截止时刻，负极区域析锂副产物浓度最大的位置在靠近隔膜处，产物浓度达到1.09×10³ mol/m³，而靠近集流体处析锂最轻微，产物浓度仅为220 mol/m³。

图8

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

图8 负极区域副反应产物浓度

Fig. 8 Concentration of side reaction products in the negative region

最后，对负极整个区域副反应产物浓度进行积分处理。图9为副反应产物中各组分对全电池容量衰减率的影响。3种产物分别造成0.2%、0.69%、3.03%的全电池容量衰减，充电结束因析锂整体上造成电池发生了3.92%的容量衰减率。

图9

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

图9 0 ℃下1 C倍率充电副反应产物造成容量变化情况

Fig. 9 Capacity change due to 1 C multiplier charging side reaction products at 0 ℃

基于上述三维电化学热耦合析锂模型可以定量描述在低温大倍率充电期间电池内部的析锂程度、温度分布等非均一现象。充电期间，极耳区域电流密度较大，该区域负极电位最先达到析锂电位0 V，过早析锂。另外，极耳区域过大的电流密度也造成该区域温度高于其他区域。上述两现象综合影响电池局部析锂的发生，总体上靠近极耳区域先发生析锂，且极耳区域最终析锂程度最严重。

1.3　电池参数设计优化

1.3.1　电极尺寸对析锂影响

为探究电极尺寸对析锂的影响，对1 C倍率充电过程中负极-隔膜处固液电势差进行分析。由图10中不同电极长度的电池局部开始析锂时间、图11中析锂时负极-隔膜处固液电势差分布情况可知，长度为20 cm时，在t_onset=205 s时靠近极耳区域固液相电势差开始小于0。底部最大值仅为2.28 mV，整个界面处电势差分布差异可以忽略。电极长度增加到40 cm时，靠近极耳处t_onset提前到193 s。分布差异也明显增大。电极长度增加到60 cm时t_onset提前到175 s开始析锂，80 cm时t_onset提前153 s，界面处最大差异达到33.4 mV。由此说明，电极长度的增加会使得电池内部析锂触发时间提前。

图10

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

图10 不同电极长度的电池局部开始析锂时间

Fig. 10 Time of partial onset of lithium precipitation for cells with different electrode lengths

图11

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

图11 不同长度电池局部开始析锂时隔膜处 φ_s-φ_e 分布

Fig. 11 φ_s-φ_e distribution at the diaphragm when lithium precipitation begins partially in cells of different lengths

进一步，对局部析锂开始触发时刻负极-隔膜界面电流密度分布分析。图12为电流密度i_z分布情况，由图12(a)可知，长度为20 cm电池达到t_onset=205 s时，i_z_{, max}为0.7 A/m²。随着电极长度增加，长度为40 cm电池开始析锂时，界面i_z_{, max}也随之增大到0.74 A/m²。长度为60 cm和80 cm，界面处i_z_{, max}分别为0.81 A/m²和0.9 A/m²，这也造成了尺寸增加后，电池内部提前达到析锂触发时刻t_onset。随着长度的增加，电极区域温度差异增大和电流密度不一致性程度增大。

图12

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

图12 不同长度电池局部开始析锂时负极-隔膜处电流密度 i_z

Fig. 12 Current i_z density at the anode-diaphragm when lithium precipitation starts par tially in cells of different lengths

图13为不同长度电极析锂量情况，即容量衰减率。在局部温度和电流密度多因素综合影响下，长度越长，析锂开始时间越早。

图13

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

图13 不同长度电极整体析锂程度

Fig. 13 Different lengths of electrodes in the overall lithium precipitation range

1.3.2　极耳位置对电池析锂影响

极耳位置的设计将影响靠近极耳区域的电流密度和温度的分布，对局部析锂的发生存在一定影响。本工作选取了长度均为40 cm的具有3种不同极耳位置的电池进行分析，以探究极耳位置对各项电性能的影响，尤其是局部析锂。极耳位置设计如图14所示，分为I、II、III三种类型，并对其在0 ℃下1 C倍率恒流充电。不同极耳类型的电池局部开始析锂时间如图15所示。

图14

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

图14 3种不同极耳位置电池结构

Fig. 14 3 different pole lug position cell structure

图15

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

图15 不同极耳类型的电池局部开始析锂时间

Fig. 15 Time of partial onset of lithium precipitation for cells with different lug types

图16为3种极耳类型电池在局部开始出现析锂时负极-隔膜界面的固液相电势差 $Δ φ$ 分布情况。I类电池在达到析锂触发时刻t_onset=193 s时，靠近极耳区域开始出现 $Δ φ$ <0，该区域首先析锂，但底部的 $Δ φ$ 为8.67 mV， $Δ φ$ 不一致性较大。然而，II类电池在t_onset为201 s时靠近极耳区域先析锂，此时中部区域 $Δ φ_{m a x}$ 接近4.02 mV。III类电池t_onset为221 s时极耳区域局部析锂，中部区域 $Δ φ_{m a x}$ 为3.07 mV， $Δ φ$ 不一致性远低于I类，接近II类。

图16

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

图16 不同极耳位置电池局部开始析锂时隔膜处 φ_s-φ_e 分布

Fig. 16 φ_s-φ_e distribution at the diaphragm when lithium precipitation starts at different lug positions

进一步对局部析锂开始触发时负极-隔膜界面电流密度i_z 分布进行分析。由电流密度分布图可知，3种极耳设计下的电池i_z_,max数值大小差异不大，如图17所示。其中，I类电池中界面处 $Δ i_{z}$ 为0.13 A/m²，另外两种电池较小，II类电池 $Δ i_{z}$ 为0.08 A/m²，III类电池则更小，仅为0.05 A/m²。说明III类电池在电流密度分布上更加均匀，也使得电池内部各区域产热均匀，温度差异更小。

图17

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

图17 不同类型极耳位置电池局部开始析锂时负极-隔膜边界处电流密度 i_z

Fig. 17 Current i_z density at the anode-diaphragm boundary when lithium precipitation begins locally in cells with different types of lug positions

极耳位置不同导致电极区域温度和电流密度不一致性程度有所差异，在多种因素的耦合下电池整体析锂程度也有所不同。如图18所示，I类电池和II类电池在温度、电流密度多种因素综合影响下析锂量几乎相同，析锂量达到3.77%。然而，III类电池在温度分布、电流密度分布不一致性要优于I和II类电池，充电截止时刻析锂量降低到3.14%，相对析锂程度降低16.7%。因此，III类电池将极耳置于长度方向轴线对侧将有利于减缓析锂程度。

图18

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

图18 不同极耳位置设计电池整体析锂程度

Fig. 18 Different lug positions to design the overall lithium precipitation degree of the battery

2 结论

本工作以P2D模型为基础，通过耦合产热模型，并在模型中加入了考虑可逆锂重嵌入机制的析锂副反应，建立了三维电化学热耦合析锂模型。通过多种方式获取模型参数后，对模型进行验证，验证结果表明了该模型在常温和低温下条件下可以较好地模拟电池端电压的变化，并且能定量描述在低温大倍率充电期间电池内部的析锂程度、温度分布等非均一现象。此外，通过仿真对电极尺寸、极耳位置等因素的分析，研究了结构设计对非均一析锂的影响。

不同电极尺寸的分析结果表明，电极长度增加会致使电极区域温度差异增大和电流密度不一致性程度增大，综合影响下使电池析锂时间略有提前，但对电池总体析锂程度影响较小。不同极耳位置的分析结果表明，当极耳位置处于电极长度方向的轴线对侧处时，电池温度和电流密度分布更加均匀，可以延缓负极析锂起始时间，能够有效缓解负极析锂，相对析锂程度降低了16.7%。

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

INUI Y, KOBAYASHI Y, WATANABE Y, et al. Simulation of temperature distribution in cylindrical and prismatic lithium ion secondary batteries[J]. Energy Conversion and Management, 2007, 48(7): 2103-2109.

[本文引用: 1]

[2]

GRANDJEAN T, BARAI A, HOSSEINZADEH E, et al. Large format lithium ion pouch cell full thermal characterisation for improved electric vehicle thermal management[J]. Journal of Power Sources, 2017, 359: 215-225.

[本文引用: 1]

[3]

SHARMA G, SOOD V K, ALAM M S, et al. Comparison of common DC and AC bus architectures for EV fast charging stations and impact on power quality[J]. eTransportation, 2020, 5: doi: 10.1016/j.etran.2020.100066.

[本文引用: 1]

[4]

HALES A, PROSSER R, BRAVO DIAZ L, et al. The cell cooling coefficient as a design tool to optimise thermal management of lithium-ion cells in battery packs[J]. eTransportation, 2020, 6: doi: 10.1016/j.etran.2020.100089.

[本文引用: 1]

[5]

FLECKENSTEIN M, BOHLEN O, ROSCHER M A, et al. Current density and state of charge inhomogeneities in Li-ion battery cells with LiFePO₄ as cathode material due to temperature gradients[J]. Journal of Power Sources, 2011, 196(10): 4769-4778.

[本文引用: 1]

[6]

GUO C L, ZHU K J, CHEN C C, et al. Characteristics and effect laws of the large-scale electric Vehicle's charging load[J]. eTransportation, 2020, 3: doi: 10.1016/j.etran.2020.100049.

[本文引用: 1]

[7]

孙占宇. 基于电化学模型的车用锂离子电池安全快速充电算法[J]. 汽车安全与节能学报, 2017, 8(1): 97-101.

[本文引用: 1]

SUN Z Y. Safe fast charging algorithm of lithium ion battery based on an electrochemical model[J]. Journal of Automotive Safety and Energy, 2017, 8(1): 97-101.

[本文引用: 1]

[8]

DIXON J, BELL K. Electric vehicles: Battery capacity, charger power, access to charging and the impacts on distribution networks[J]. eTransportation, 2020, 4: doi: 10.1016/j.etran.2020.100059.

[本文引用: 1]

[9]

彭敏, 申文静, 罗兆东. 层叠式锂离子电池二维热模型研究[J]. 电源技术, 2018, 42(9): 1312-1315.

[本文引用: 1]

PENG M, SHEN W J, LUO Z D. Two-dimensional thermal modeling for laminated lithium ion battery[J]. Chinese Journal of Power Sources, 2018, 42(9): 1312-1315.

[本文引用: 1]

[10]

STURM J, RHEINFELD A, ZILBERMAN I, et al. Modeling and simulation of inhomogeneities in a 18650 nickel-rich, silicon-graphite lithium-ion cell during fast charging[J]. Journal of Power Sources, 2019, 412: 204-223.

[本文引用: 1]

[11]

ZHAO W, LUO G, WANG C Y. Effect of tab design on large-format Li-ion cell performance[J]. Journal of Power Sources, 2014, 257: 70-79.

[本文引用: 1]

[12]

SAMBA A, OMAR N, GUALOUS H, et al. Impact of tab location on large format lithium-ion pouch cell based on fully coupled tree-dimensional electrochemical-thermal modeling[J]. Electrochimica Acta, 2014, 147: 319-329.

[本文引用: 1]

[13]

FEAR C, PARMANANDA M, KABRA V, et al. Mechanistic underpinnings of thermal gradient induced inhomogeneity in lithium plating[J]. Energy Storage Materials, 2021, 35: 500-511.

[本文引用: 1]

[14]

DOYLE M, FULLER T F, NEWMAN J. Modeling of galvanostatic charge and discharge of the lithium/polymer/insertion cell[J]. Journal of the Electrochemical Society, 1993, 140(6): 1526-1533.

[本文引用: 1]

[15]

DOYLE M, NEWMAN J, GOZDZ A S, et al. Comparison of modeling predictions with experimental data from plastic lithium ion cells[J]. Journal of the Electrochemical Society, 1996, 143(6): 1890-1903.

[本文引用: 1]

[16]

TANG S Q, WANG Z X, GUO H J, et al. Systematic parameter acquisition method for electrochemical model of 4.35 V LiCoO₂ batteries[J]. Solid State Ionics, 2019, 343: doi: 10.1016/j.ssi.2019.115083.

[本文引用: 1]

[17]

YE Y H, SHI Y X, CAI N S, et al. Electro-thermal modeling and experimental validation for lithium ion battery[J]. Journal of Power Sources, 2012, 199: 227-238.

[本文引用: 1]

[18]

李斌, 常国峰, 林春景, 等. 车用动力锂电池产热机理研究现状[J]. 电源技术, 2014, 38(2): 378-381.

[本文引用: 1]

LI B, CHANG G F, LIN C J, et al. Research on heat generate mechanism of Li-ion batteries for electric vehicles[J]. Chinese Journal of Power Sources, 2014, 38(2): 378-381.

[本文引用: 1]

[19]

REN D S, SMITH K, GUO D X, et al. Investigation of lithium plating-stripping process in Li-ion batteries at low temperature using an electrochemical model[J]. Journal of the Electrochemical Society, 2018, 165(10): doi: 10.1149/2.0661810jes.

[本文引用: 6]

... 在当前的快充场景中，电池有可能会出现负极析锂现象^[1-2].这种现象会导致电池管理系统对电池单体监测异常^[3]，从而造成电池过热和热稳定性降低等问题，甚至可能引发热失控^[4-5].因此，电池的充电安全问题亟须解决^[6-7]. ...

... 抑制电池充电时发生析锂可以从电池设计方面入手^[8-9].Sturm等^[10]基于多维结构模型，以负极过电位和温度评估了极片和极耳的设计对负极析锂的影响，并通过设定过电位阈值将高倍率充电模式下的析锂风险降到最低.Zhao等^[11]建立了三维模型，通过多极耳的设计有效降低了长电极的集流体带来的电子传递损失，提高了能量密度.Samba等^[12]对软包电池进行三维建模，从电流密度分布、电势、温度分布等角度分析了极耳位置对电池带来的影响.其研究结果表明，对称结构的设计使得电池性能更加均匀，电流密度梯度更小，能够使活性材料最大程度得到利用.然而，常规的电化学热耦合模型没有考虑负极发生析锂后的可逆锂重嵌入机制，无法很好地描述电池在充电过程中实际的析锂状况.因此，添加可逆锂重嵌入机制的多维结构电化学热耦合析锂模型将能够更好地分析电池内部非均一析锂问题^[13]. ...

... 本工作所使用的模型是在Newman和Doyle等^[14-15]建立的伪二维(pseudo two dimensional，P2D)机理模型基础上进一步建立完成的. ...

... 本工作参数获取方式主要为拆解测量、实验标定、参数辨识、文献参考等，模型具体参数和分类如表1所示.模型中设计参数部分，如极片厚度L、集流体厚度L_cc、正负极厚度L_p/n、隔膜厚度L_sep、活性物质总反应面积A_total等参数可通过工具测量获得.而正负极均衡电势曲线U_{ref, p/n}则通过制作纽扣电池，以极小电流的充放电实验来获取.固相颗粒最大嵌锂浓度c_{s, max}，初始时刻化学计量比x_0/100、y_0/100等一般通过参数辨识进行获取.正负极固相颗粒半径r_{s, n/p}、电解质初始盐浓度cl₀、孔隙率ε等由电池厂商测定后提供.模型中的热物性参数，如比热容C_p 、对流换热系数h等通过实验进行测定.液相参数中的液相扩散系数D_e、液相离子电导率与温度T及电解质浓度c有关，可以直接引用文献中的参考值^[16].正负极材料的熵热系数通过对COMSOL案例库中已录入的数据进行适当调整即可得到有效的dU/dT曲线.由于固相扩散系数D_s、反应速率常数k与温度T密切相关且符合Arrhenius公式，可以通过不同温度下的脉冲实验标定得到相应数值和活化能E_a. ...

... 注：①为实际测量值，②为参数辨识数据，③为电池厂商提供数据，④为文献参考值^[17-18]，⑤为标定数据. ...

... Parameters related to the kinetic reactions of analytical and reversible lithium re-embedding

Table 3

参数	物理意义	单位	负极
α_{a, pl}、α_c,_pl	析锂/锂溶解传递系数	1	0.3、0.7^[19]
k_pl/st	锂析出/溶解反应速率常数	m/s	6.5×10^-6
E_{a, pl/st}	析锂与锂溶解反应活化能	J/mol	50000
z₁、z₂、z₃	可逆锂/不可逆锂/SEI占比	1	0.57、0.375、0.05
σ_SEI	SEI膜电导率	S/m	5×10^-6[19]
δ_SEI	SEI膜t厚度	m	式(10)
ρ_SEI	SEI 膜密度	kg/m³	1690^[19]
ρ_Li	锂金属密度	kg/m³	534^[19]
M_SEI	SEI膜摩尔质量	kg/mol	0.162^[19]
M_Li	锂金属摩尔质量	kg/mol	6.94×10^-3[19]