CN110991049B - 一种基于过充锂离子电池的热失控仿真方法 - Google Patents

一种基于过充锂离子电池的热失控仿真方法 Download PDF

Info

Publication number
CN110991049B
CN110991049B CN201911230805.4A CN201911230805A CN110991049B CN 110991049 B CN110991049 B CN 110991049B CN 201911230805 A CN201911230805 A CN 201911230805A CN 110991049 B CN110991049 B CN 110991049B
Authority
CN
China
Prior art keywords
lithium battery
thermal runaway
heat generation
model
reaction
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201911230805.4A
Other languages
English (en)
Other versions
CN110991049A (zh
Inventor
胡广地
黄文才
郭峰
赵利民
赛影辉
胡坚耀
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN201911230805.4A priority Critical patent/CN110991049B/zh
Publication of CN110991049A publication Critical patent/CN110991049A/zh
Application granted granted Critical
Publication of CN110991049B publication Critical patent/CN110991049B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/10—Energy storage using batteries

Landscapes

  • Secondary Cells (AREA)
  • Battery Electrode And Active Subsutance (AREA)

Abstract

本发明公开了一种基于过充锂离子电池的热失控仿真方法，基于电池的电化学反应机理和产热特性原理，通过对锂离子电池进行热电特性试验和建模仿真计算，研究充电条件下电池温度变化特性，并根据该特性进而建立热滥用模型，从而分析出不同参数下的电池单体性质，实现过充热失控工况下电池内部化学反应变化的探究，解决了现有基于过充锂离子电池的热失控仿真方法缺乏对电化学反应机理和产热特性的全面分析，且未考虑不同电池单体的特异性的问题。

Description

一种基于过充锂离子电池的热失控仿真方法
技术领域
本发明涉及电池领域，具体涉及一种基于过充锂离子电池的热失控仿真方法。
背景技术
随着电子技术的发展，越来越多的电子产品采用锂电池作为可持续储能设备，同时由于不正当操作特别是过充造成的锂电池安全事故频繁发生。这种事故的内在机理是在锂电池充电过程中会伴随着一部分热量的产生，当充电达到额定电压和电量之后继续充电会造成热量的累积，累积的热量过大便会造成电池温度急剧升高引发副反应，严重情况下即会引起燃烧、爆炸等事故。过充现象在生活中非常普遍，无论是手机、电脑等小型电子设备，还是纯电动汽车、电动自行车等大型设备都存在过充问题，因此对锂电池过充现象而导致热失控的研究成为学术界的重要方向之一。
然而现有的锂离子电池过充的热失控仿真模型存在不少缺陷，尤其是缺乏对电化学反应机理和产热特性的全面分析，且未考虑不同电池单体的特异性。
发明内容
针对现有技术中的上述不足，本发明提供的一种基于过充锂离子电池的热失控仿真方法，基于电池的电化学反应机理和产热特性原理，通过对锂离子电池进行热电特性试验和建模仿真计算，研究充电条件下电池温度变化特性，并根据该特性进而建立热滥用模型，从而分析出不同参数下的电池单体性质，实现过充热失控工况下电池内部化学反应变化的探究，解决了现有基于过充锂离子电池的热失控仿真方法缺乏对电化学反应机理和产热特性的全面分析，且未考虑不同电池单体的特异性的问题。使用本发明提供的基于过充锂离子电池的热失控仿真方法，相比于现有基于过充锂离子电池的热失控仿真方法，能更有效地分析和预测电池的热行为。
为了达到上述发明目的，本发明采用的技术方案为：一种基于过充锂离子电池的热失控仿真方法，包括以下步骤：
S1：对锂电池样本进行充放电试验，得到锂电池样本充放电过程中的时域电压数据U(t)、时域电流数据I(t)和时域表面温度数据T(t)；
S2：对锂电池样本进行多重脉冲试验，得到锂电池样本的时域过电位变化数据OCVn(t)和时域剩余电量变化数据SOCn(t)，其中n为1～N的正整数，代表第n次试验，N为试验总次数；
S3：根据锂电池样本内部的各种材料具有各向异性导热系数的特征，按锂电池样本物质材料进行分层，建立锂电池分层结构模型；
S4：根据电池电化学反应热物理方程式，在三维空间中使用有限差分法对锂电池分层结构模型进行分格式建模，建立锂电池电化学产热模型；
S5：使用有限元法求解锂电池电化学产热模型，得到锂电池仿真时域电压和仿真时域温度数据；
S6：用锂电池样本充放电过程中的时域电压数据U(t)、时域电流数据I(t)、时域表面温度数据T(t)、时域过电位变化数据OCVn(t)和时域剩余电量变化数据SOCn(t)与锂电池仿真时域电压和仿真时域温度数据进行比对，验证锂电池电化学产热模型，判断温度最大误差Terror是否小于预设的温度误差阈值Tth，若是，则跳转至步骤S7，若否，则跳转至步骤S3；
S7、判断电压最大误差Uerror是否小于预设的电压误差阈值Uth，若是，则验证通过，并跳转至步骤S8，若否，则验证不通过，跳转至步骤S3；
S8：根据锂电池热失控反应方程组，在锂电池电化学产热模型基础上建立锂电池过充热失控模型进行过充模拟，得到不同高温阶段，锂电池逐步损毁时的表面温度分布特性以及电压变化曲线；
S9：对锂电池过充模拟得到的数据进行分析，得到锂电池热失控特征，从而预测真实环境中锂离子电池的热行为。
进一步地：步骤S3中建立锂电池分层结构模型时，对模型内部空间各处的电流密度进行相等近似处理。
进一步地：步骤S3中建立锂电池分层结构模型时，对各层材料模型进行介质均匀、密度一致的近似处理。
进一步地：步骤S4中在三维空间中使用有限差分法对锂电池分层结构模型进行分格式建模，建立锂电池电化学产热模型的步骤为：
S41：录入锂电池分层结构模型，并根据锂电池样本的材料属性设置参数；
S42：根据锂电池样本的电极和隔膜中的容积组分，录入电池的阳极、分隔膜和阴极中电解液、石墨、粘合剂、膜聚合物和镍钴锰酸锂的成分占比；
S43：对锂电池分层结构模型进行网格划分，形成相应的求解域；
S44：对求解域进行有限元建模，实现锂电池电化学产热模型的建立。
进一步地：步骤S8中的锂电池热失控反应方程组为：
Q＝Qrea+Qact+Qohm+Qp-Qdis(1)
Figure BDA0002303483470000031
Figure BDA0002303483470000032
Figure BDA0002303483470000033
Qp＝Qsei+Qneg+Qpos+Qele (5)
Qdis＝hA(T-T外)(6)
其中，Q为热失控过程中的产热速率，Qrea为可逆热，Qact为不可逆热，Qohm为欧姆热，Qp为副反应热，Qdis为对流热，σ为固相电导率，eff有效系数，
Figure BDA0002303483470000041
为固相的电位，
Figure BDA0002303483470000042
为液相的电位，
Figure BDA0002303483470000043
拉普拉斯算子，c2为液相电解质溶液中的浓度，k为电极反应速率常数，jLi为电极表面局部电流密度，Uj电路局部电流产生的电压，Qsei为固液膜的分解反应的产热速率，Qneg为负极材料分解反应的产热速率，Qpos为正极材料分解反应的产热速率，Qele为电解液分解的产热速率，h为对流传热系数，A为对流传热面积，T电池内环境温度，T外电池外环境温度。
进一步地：固液膜的分解反应的产热速率Qsei的计算公式为：
Qsei＝HseiWseiRsei (7)
Figure BDA0002303483470000044
Figure BDA0002303483470000045
其中，csei为亚稳态SEI的无量纲量；Rsei为膜结构中R型频率因子，Asei为膜结构A型频率因子；Ea,sei为膜结构反应的活化能，R为通用气体常数，Hsei是SEI膜的分解生热率，Wsei单位膜结构体积碳含量。
进一步地：负极材料分解反应的产热速率Qneg的计算公式为：
Qneg＝HnegWnegRneg(10)
Figure BDA0002303483470000046
Figure BDA0002303483470000047
其中，cneg为碳层中锂离子含量的无量纲量，Rneg为负极材料R型的频率因子，Aneg为负极材料A型的频率因子，Wneg为负极材料单位体积碳含量，Hneg为负极材料反应生热率，Ea,neg为负极材料的反应的活化能；R为通用气体常数。
进一步地：正极材料分解反应的产热速率Qpos的计算公式为：
Qpos＝HposWposRpos (13)
Figure BDA0002303483470000048
Figure BDA0002303483470000051
其中，Hpos为正极材料的反应生热率，Wpos为单位体积正极活性材料含量，Rpos为正极材料R型频率因子，Apos为正极材料A型频率因子，Ea,pos为正极材料反应的活化能，α为正极材料的转化率，R为通用气体常数。
进一步地：电解液分解的产热速率Qele的计算公式为：
Qele＝HeleWeleRele (16)
Figure BDA0002303483470000052
Figure BDA0002303483470000053
其中，Hele为电解液分解频率因子，Wele单位体积电解质含量，Rele为电解液R型频率因子，Aele为电解液A型频率因子，cele为电解质浓度的无量纲量，Ea,ele为电解液反应活化能，R为通用气体常数。
进一步地：步骤S8根据锂电池热失控反应方程组，在锂电池电化学产热模型基础上建立锂电池过充热失控模型进行过充模拟的步骤为：
S81：将锂电池热失控反应方程组设定为锂电池过充热失控模型的有限差分基本方程；
S82：对锂电池进行初始加热，将温度升高至80℃，使得锂电池的固液膜开始分解，记录此时模型仿真得到的锂电池表面温度分布特性以及电压变化曲线；
S83：对锂电池进一步加热至120℃，使得锂电池的负极材料开始反应分解，记录此时模型仿真得到的锂电池表面温度分布特性以及电压变化曲线；
S84：对锂电池进一步加热至140℃，使得锂电池的正极材料开始反应分解，记录此时模型仿真得到的锂电池表面温度分布特性以及电压变化曲线；
S85：对锂电池进一步加热至200℃，使得锂电池的电解液开始反应分解，记录此时模型仿真得到的锂电池表面温度分布特性以及电压变化曲线；
S86：对锂电池进一步加热直至锂电池的电解液完全分解，记录此时模型仿真得到的锂电池表面温度分布特性以及电压变化曲线。
本发明的有益效果为：首先根据锂电池样本内部的各种材料具有各向异性导热系数的特征，按锂电池样本物质材料进行了分层设置，而不是简单地将电池视为一个块状整体，充分考虑了不同材料的物理、化学特性；然后根据电池电化学反应热物理方程式，在三维空间中使用有限差分法对锂电池分层结构模型进行分格式建模，建立锂电池电化学产热模型，并用真实环境中的实验数据对该模型进行了验证，保障了模型的可靠性；然后在锂电池热失控反应方程组的设定中不仅考虑了电化学热耦合的基本方程，还引入了副反应方程，并依此对锂电池的五阶段热损毁的情况进行了严密的仿真和数据记录，数据全面而又可靠。本方法的模型也可灵活设置不同锂离子电池单体的参数，实现特异锂电池的过充热失控工况下电池内部化学反应变化的探究。
附图说明
图1为一种基于过充锂离子电池的热失控仿真方法的流程图；
图2为一种基于过充锂离子电池的热失控仿真方法的时域过电位变化数据OCV1(t)和时域剩余电量变化数据SOC1(t)的关系变化曲线；
图3为一种基于过充锂离子电池的热失控仿真方法的时域电压数据U(t)与锂电池仿真时域电压的关系变化曲线；
图4为一种基于过充锂离子电池的热失控仿真方法的时域表面温度数据T(t)与仿真时域温度数据的关系变化曲线；
图5为一种基于过充锂离子电池的热失控仿真方法的为环境温度为77℃时锂电池表面温度分布特性；
图6为一种基于过充锂离子电池的热失控仿真方法的为环境温度为175℃时锂电池表面温度分布特性。
具体实施方式
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
如图1所示，在本发明的一个实施例中，一种基于过充锂离子电池的热失控仿真方法，包括以下步骤：
S1：对锂电池样本进行充放电试验，得到锂电池样本充放电过程中的时域电压数据U(t)、时域电流数据I(t)和时域表面温度数据T(t)；
S2：对锂电池样本进行多重脉冲试验，得到锂电池样本的时域过电位变化数据OCVn(t)和时域剩余电量变化数据SOCn(t)，其中n为1～N的正整数，代表第n次试验，N为试验总次数；图2为第1次试验，时域过电位变化数据OCV1(t)和时域剩余电量变化数据SOC1(t)的关系曲线；
S3：根据锂电池样本内部的各种材料具有各向异性导热系数的特征，按锂电池样本物质材料进行分层，建立锂电池分层结构模型；
步骤S3中建立锂电池分层结构模型时，对模型内部空间各处的电流密度进行相等近似处理。
步骤S3中建立锂电池分层结构模型时，对各层材料模型进行介质均匀、密度一致的近似处理。
S4：根据电池电化学反应热物理方程式，在三维空间中使用有限差分法对锂电池分层结构模型进行分格式建模，建立锂电池电化学产热模型；
步骤S4中在三维空间中使用有限差分法对锂电池分层结构模型进行分格式建模，建立锂电池电化学产热模型的步骤为：
S41：在COMSOL Multiphysics软件中录入锂电池分层结构模型，并根据锂电池样本的材料属性设置参数，所述参数包括：极反应速率、阴极反应速率、扭曲布鲁格曼系数、初始电解质盐浓度、最大负极容量、初始固态锂离子浓度、电解质相体积分数、固相体积分数、电池材料厚度、粒子半径、电导率、扩散系数、电池材料密度、电池材料比热容和电池材料导热系数，如下表1和表2所示。
表1
Figure BDA0002303483470000081
表2
Figure BDA0002303483470000082
S42：根据锂电池样本的电极和隔膜中的容积组分，在COMSOL Multiphysics软件中录入电池的阳极、分隔膜和阴极中电解液、石墨、粘合剂、膜聚合物和镍钴锰酸锂的成分占比，如表3所示。
表3
参数 阳极 分隔膜 阴极
电解液 0.444 0.37 0.4
石墨 0.384 / /
粘合剂 0.172 / 0.17
膜聚合物 / 0.63 /
镍钴锰酸锂 / / 0.43
S43：利用COMSOL Multiphysics软件对锂电池分层结构模型进行网格划分，形成相应的求解域；
S44：利用COMSOL Multiphysics软件自带的锂电池场与温度场电化学反应热物理方程式，对求解域进行有限元建模，实现锂电池电化学产热模型的建立。
S5：使用有限元法求解锂电池电化学产热模型，得到锂电池仿真时域电压和仿真时域温度数据；
S6：用锂电池样本充放电过程中的时域电压数据U(t)、时域电流数据I(t)、时域表面温度数据T(t)、时域过电位变化数据OCVn(t)和时域剩余电量变化数据SOCn(t)与锂电池仿真时域电压和仿真时域温度数据进行比对，如图3～4所示，验证锂电池电化学产热模型，判断温度最大误差Terror是否小于预设的温度误差阈值Tth，若是，则跳转至步骤S7，若否，则跳转至步骤S3；
S7、判断电压最大误差Uerror是否小于预设的电压误差阈值Uth，若是，则验证通过，并跳转至步骤S8，若否，则验证不通过，跳转至步骤S3；
S8：根据锂电池热失控反应方程组，在锂电池电化学产热模型基础上建立锂电池过充热失控模型进行过充模拟，得到不同高温阶段，锂电池逐步损毁时的表面温度分布特性以及电压变化曲线；
步骤S8中的锂电池热失控反应方程组为：
Q＝Qrea+Qact+Qohm+Qp-Qdis(1)
Figure BDA0002303483470000101
Figure BDA0002303483470000102
Figure BDA0002303483470000103
Qp＝Qsei+Qneg+Qpos+Qele (5)
Qdis＝hA(T-T外)(6)
其中，Q为热失控过程中的产热速率，Qrea为可逆热，Qact为不可逆热，Qohm为欧姆热，Qp为副反应热，Qdis为对流热，σ为固相电导率，eff有效系数，
Figure BDA0002303483470000104
为固相的电位，
Figure BDA0002303483470000105
为液相的电位，
Figure BDA0002303483470000108
拉普拉斯算子，c2为液相电解质溶液中的浓度，k为电极反应速率常数，jLi为电极表面局部电流密度，Uj电路局部电流产生的电压，Qsei为固液膜的分解反应的产热速率，Qneg为负极材料分解反应的产热速率，Qpos为正极材料分解反应的产热速率，Qele为电解液分解的产热速率，h为对流传热系数，A为对流传热面积，T电池内环境温度，T外电池外环境温度。
固液膜的分解反应的产热速率Qsei的计算公式为：
Qsei＝HseiWseiRsei (7)
Figure BDA0002303483470000106
Figure BDA0002303483470000107
其中，csei为亚稳态SEI的无量纲量；Rsei为膜结构中R型频率因子，Asei为膜结构A型频率因子；Ea,sei为膜结构反应的活化能，R为通用气体常数，Hsei是SEI膜的分解生热率，Wsei单位膜结构体积碳含量。
负极材料分解反应的产热速率Qneg的计算公式为：
Qneg＝HnegWnegRneg (10)
Figure BDA0002303483470000111
Figure BDA0002303483470000112
其中，cneg为碳层中锂离子含量的无量纲量，Rneg为负极材料R型的频率因子，Aneg为负极材料A型的频率因子，Wneg为负极材料单位体积碳含量，Hneg为负极材料反应生热率，Ea,neg为负极材料的反应的活化能；R为通用气体常数。
正极材料分解反应的产热速率Qpos的计算公式为：
Qpos＝HposWposRpos (13)
Figure BDA0002303483470000113
Figure BDA0002303483470000114
其中，Hpos为正极材料的反应生热率，Wpos为单位体积正极活性材料含量，Rpos为正极材料R型频率因子，Apos为正极材料A型频率因子，Ea,pos为正极材料反应的活化能，α为正极材料的转化率，R为通用气体常数。
电解液分解的产热速率Qele的计算公式为：
Qele＝HeleWeleRele (16)
Figure BDA0002303483470000115
Figure BDA0002303483470000116
其中，Hele为电解液分解频率因子，Wele单位体积电解质含量，Rele为电解液R型频率因子，Aele为电解液A型频率因子，cele为电解质浓度的无量纲量，Ea,ele为电解液反应活化能，R为通用气体常数。
步骤S8根据锂电池热失控反应方程组，在锂电池电化学产热模型基础上建立锂电池过充热失控模型进行过充模拟的步骤为：
S81：通过COMSOL Multiphysics软件的UDF自定义域微分形式方程工具将锂电池热失控反应方程组设定为锂电池过充热失控模型的有限差分基本方程；
S82：在COMSOL Multiphysics软件中对锂电池进行初始加热，将温度升高至80℃，使得锂电池的固液膜开始分解，记录此时模型仿真得到的锂电池表面温度分布特性以及电压变化曲线，图5为温度为77℃时，锂电池表面温度分布特性图；
S83：对锂电池进一步加热至120℃，使得锂电池的负极材料开始反应分解，记录此时模型仿真得到的锂电池表面温度分布特性以及电压变化曲线；
S84：对锂电池进一步加热至140℃，使得锂电池的正极材料开始反应分解，记录此时模型仿真得到的锂电池表面温度分布特性以及电压变化曲线；
S85：对锂电池进一步加热至200℃，使得锂电池的电解液开始反应分解，记录此时模型仿真得到的锂电池表面温度分布特性以及电压变化曲线；图6为温度为175℃时，锂电池表面温度分布特性图；
S86：对锂电池进一步加热，直至锂电池的电解液完全分解，记录此时模型仿真得到的锂电池表面温度分布特性以及电压变化曲线。
S9：对锂电池过充模拟得到的数据进行分析，得到锂电池热失控特征，从而预测真实环境中锂离子电池的热行为。
本发明的有益效果为：首先根据锂电池样本内部的各种材料具有各向异性导热系数的特征，按锂电池样本物质材料进行了分层设置，而不是简单地将电池视为一个块状整体，充分考虑了不同材料的物理、化学特性；然后根据电池电化学反应热物理方程式，在三维空间中使用有限差分法对锂电池分层结构模型进行分格式建模，建立锂电池电化学产热模型，并用真实环境中的实验数据对该模型进行了验证，保障了模型的可靠性；然后在锂电池热失控反应方程组的设定中不仅考虑了电化学热耦合的基本方程，还引入了副反应方程，并依此对锂电池的五阶段热损毁的情况进行了严密的仿真和数据记录，数据全面而又可靠。本方法的模型也可灵活设置不同锂离子电池单体的参数，实现特异锂电池的过充热失控工况下电池内部化学反应变化的探究。

Claims (10)

1.一种基于过充锂离子电池的热失控仿真方法，其特征在于，包括以下步骤：
S1：对锂电池样本进行充放电试验，得到锂电池样本充放电过程中的时域电压数据U(t)、时域电流数据I(t)和时域表面温度数据T(t)；
S2：对锂电池样本进行多重脉冲试验，得到锂电池样本的时域过电位变化数据OCVn(t)和时域剩余电量变化数据SOCn(t)，其中n为1～N的正整数，代表第n次试验，N为试验总次数；
S3：根据锂电池样本内部的各种材料具有各向异性导热系数的特征，按锂电池样本物质材料进行分层，建立锂电池分层结构模型；
S4：根据电池电化学反应热物理方程式，在三维空间中使用有限差分法对锂电池分层结构模型进行分格式建模，建立锂电池电化学产热模型；
S5：使用有限元法求解锂电池电化学产热模型，得到锂电池仿真时域电压和仿真时域温度数据；
S6：用锂电池样本充放电过程中的时域电压数据U(t)、时域电流数据I(t)、时域表面温度数据T(t)、时域过电位变化数据OCVn(t)和时域剩余电量变化数据SOCn(t)与锂电池仿真时域电压和仿真时域温度数据进行比对，验证锂电池电化学产热模型，判断温度最大误差Terror是否小于预设的温度误差阈值Tth，若是，则跳转至步骤S7，若否，则跳转至步骤S3；
S7、判断电压最大误差Uerror是否小于预设的电压误差阈值Uth，若是，则验证通过，并跳转至步骤S8，若否，则验证不通过，跳转至步骤S3；
S8：根据锂电池热失控反应方程组，在锂电池电化学产热模型基础上建立锂电池过充热失控模型进行过充模拟，得到不同高温阶段锂电池逐步损毁时的表面温度分布特性以及电压变化曲线；
S9：对锂电池过充模拟得到的数据进行分析，得到锂电池热失控特征，从而预测真实环境中锂离子电池的热行为。
2.根据权利要求1所述的基于过充锂离子电池的热失控仿真方法，其特征在于，所述步骤S3中建立锂电池分层结构模型时，对模型内部空间各处的电流密度进行相等近似处理。
3.根据权利要求1所述的基于过充锂离子电池的热失控仿真方法，其特征在于，所述步骤S3中建立锂电池分层结构模型时，对各层材料模型进行介质均匀、密度一致的近似处理。
4.根据权利要求1所述的基于过充锂离子电池的热失控仿真方法，其特征在于，所述步骤S4中在三维空间中使用有限差分法对锂电池分层结构模型进行分格式建模，建立锂电池电化学产热模型的步骤为：
S41：录入锂电池分层结构模型，并根据锂电池样本的材料属性设置参数；
S42：根据锂电池样本的电极和隔膜中的容积组分，录入电池的阳极、分隔膜和阴极中电解液、石墨、粘合剂、膜聚合物和镍钴锰酸锂的成分占比；
S43：对锂电池分层结构模型进行网格划分，形成相应的求解域；
S44：对求解域进行有限元建模，实现锂电池电化学产热模型的建立。
5.根据权利要求1所述的基于过充锂离子电池的热失控仿真方法，其特征在于，所述步骤S8中的锂电池热失控反应方程组为：
Q＝Qrea+Qact+Qohm+Qp-Qdis (1)
Figure FDA0002303483460000021
Figure FDA0002303483460000022
Figure FDA0002303483460000023
Qp＝Qsei+Qneg+Qpos+Qele (5)
Qdis＝hA(T-T外) (6)
其中，Q为热失控过程中的产热速率，Qrea为可逆热，Qact为不可逆热，Qohm为欧姆热，Qp为副反应热，Qdis为对流热，σ为固相电导率，eff有效系数，
Figure FDA0002303483460000031
为固相的电位，
Figure FDA0002303483460000032
为液相的电位，
Figure FDA0002303483460000037
拉普拉斯算子，c2为液相电解质溶液中的浓度，k为电极反应速率常数，jLi为电极表面局部电流密度，Uj电路局部电流产生的电压，Qsei为固液膜的分解反应的产热速率，Qneg为负极材料分解反应的产热速率，Qpos为正极材料分解反应的产热速率，Qele为电解液分解的产热速率，h为对流传热系数，A为对流传热面积，T电池内环境温度，T处电池外环境温度。
6.根据权利要求5所述的基于过充锂离子电池的热失控仿真方法，其特征在于，所述固液膜的分解反应的产热速率Qsei的计算公式为：
Qsei＝HseiWseiRsei (7)
Figure FDA0002303483460000033
Figure FDA0002303483460000034
其中，csei为亚稳态SEI的无量纲量；Rsei为膜结构中R型频率因子，Asei为膜结构A型频率因子；Ea，sei为膜结构反应的活化能，R为通用气体常数，Hsei是SEI膜的分解生热率，Wsei单位膜结构体积碳含量。
7.根据权利要求5所述的基于过充锂离子电池的热失控仿真方法，其特征在于，所述负极材料分解反应的产热速率Qneg的计算公式为：
Qneg＝HnegWnegRneg (10)
Figure FDA0002303483460000035
Figure FDA0002303483460000036
其中，cneg为碳层中锂离子含量的无量纲量，Rneg为负极材料R型的频率因子，Aneg为负极材料A型的频率因子，Wneg为负极材料单位体积碳含量，Hneg为负极材料反应生热率，Ea，neg为负极材料的反应的活化能；R为通用气体常数。
8.根据权利要求5所述的基于过充锂离子电池的热失控仿真方法，其特征在于，所述正极材料分解反应的产热速率Qpos的计算公式为：
Qpos＝HposWposRpos (13)
Figure FDA0002303483460000041
Figure FDA0002303483460000042
其中，Hpos为正极材料的反应生热率，Wpos为单位体积正极活性材料含量，Rpos为正极材料R型频率因子，Apos为正极材料A型频率因子，Ea，pos为正极材料反应的活化能，α为正极材料的转化率，R为通用气体常数。
9.根据权利要求5所述的基于过充锂离子电池的热失控仿真方法，其特征在于，所述电解液分解的产热速率Qele的计算公式为：
Qele＝HeleWeleRele (16)
Figure FDA0002303483460000043
Figure FDA0002303483460000044
其中，Hele为电解液分解频率因子，Wele单位体积电解质含量，Rele为电解液R型频率因子，Aele为电解液A型频率因子，cele为电解质浓度的无量纲量，Ea，ele为电解液反应活化能，R为通用气体常数。
10.根据权利要求1所述的基于过充锂离子电池的热失控仿真方法，其特征在于，所述步骤S8根据锂电池热失控反应方程组，在锂电池电化学产热模型基础上建立锂电池过充热失控模型进行过充模拟的步骤为：
S81：将锂电池热失控反应方程组设定为锂电池过充热失控模型的有限差分基本方程；
S82：对锂电池进行初始加热，将温度升高至80℃，使得锂电池的固液膜开始分解，记录此时模型仿真得到的锂电池表面温度分布特性以及电压变化曲线；
S83：对锂电池进一步加热至120℃，使得锂电池的负极材料开始反应分解，记录此时模型仿真得到的锂电池表面温度分布特性以及电压变化曲线；
S84：对锂电池进一步加热至140℃，使得锂电池的正极材料开始反应分解，记录此时模型仿真得到的锂电池表面温度分布特性以及电压变化曲线；
S85：对锂电池进一步加热至200℃，使得锂电池的电解液开始反应分解，记录此时模型仿真得到的锂电池表面温度分布特性以及电压变化曲线；
S86：对锂电池进一步加热直至锂电池的电解液完全分解，记录此时模型仿真得到的锂电池表面温度分布特性以及电压变化曲线。
CN201911230805.4A 2019-12-05 2019-12-05 一种基于过充锂离子电池的热失控仿真方法 Active CN110991049B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911230805.4A CN110991049B (zh) 2019-12-05 2019-12-05 一种基于过充锂离子电池的热失控仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911230805.4A CN110991049B (zh) 2019-12-05 2019-12-05 一种基于过充锂离子电池的热失控仿真方法

Publications (2)

Publication Number Publication Date
CN110991049A CN110991049A (zh) 2020-04-10
CN110991049B true CN110991049B (zh) 2022-07-12

Family

ID=70090196

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911230805.4A Active CN110991049B (zh) 2019-12-05 2019-12-05 一种基于过充锂离子电池的热失控仿真方法

Country Status (1)

Country Link
CN (1) CN110991049B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111562499B (zh) * 2020-06-09 2022-11-11 温州长江汽车电子有限公司 一种新能源汽车锂动力电池热管理模拟方法
CN112380679A (zh) * 2020-11-02 2021-02-19 中国第一汽车股份有限公司 电池热失控的模拟方法、装置、设备及存储介质
CN113076518B (zh) * 2021-04-07 2022-04-15 湖南城市学院 一种锂离子软包电池单体内部温度压力估计方法
CN114035053A (zh) * 2021-11-11 2022-02-11 合肥工业大学 一种锂电池热失控预警方法及其系统
CN114279891B (zh) * 2021-12-29 2024-08-13 中国科学技术大学 一种锂离子电池热失控过程中产气速率原位测量方法
CN115062474B (zh) * 2022-06-17 2023-01-06 上海玫克生储能科技有限公司 一种锂电池的工况仿真方法、系统、终端设备及存储介质
CN115863810A (zh) * 2022-12-14 2023-03-28 国联汽车动力电池研究院有限责任公司 一种定量获取蓄电池过充过放过程中产热的方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103345561A (zh) * 2013-07-16 2013-10-09 上海航天电源技术有限责任公司 一种具有多层片状结构的锂离子电池热仿真分析方法
CN107145628A (zh) * 2017-03-31 2017-09-08 中南大学 基于电化学‑热耦合模型的预测锂电池循环寿命的方法
CN107390136A (zh) * 2017-08-15 2017-11-24 北京航空航天大学 一种老化锂离子电池热失控建模方法
CN109581064A (zh) * 2018-12-28 2019-04-05 蜂巢能源科技有限公司 测量电池低温直流电阻的方法、装置、介质
CN110457742A (zh) * 2019-06-26 2019-11-15 北京航空航天大学 一种锂离子电池电化学-热-机械耦合模型及建模方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10003111B2 (en) * 2013-05-13 2018-06-19 The Regents Of The University Of California Thermal interface composite material and method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103345561A (zh) * 2013-07-16 2013-10-09 上海航天电源技术有限责任公司 一种具有多层片状结构的锂离子电池热仿真分析方法
CN107145628A (zh) * 2017-03-31 2017-09-08 中南大学 基于电化学‑热耦合模型的预测锂电池循环寿命的方法
CN107390136A (zh) * 2017-08-15 2017-11-24 北京航空航天大学 一种老化锂离子电池热失控建模方法
CN109581064A (zh) * 2018-12-28 2019-04-05 蜂巢能源科技有限公司 测量电池低温直流电阻的方法、装置、介质
CN110457742A (zh) * 2019-06-26 2019-11-15 北京航空航天大学 一种锂离子电池电化学-热-机械耦合模型及建模方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Liang et al..A novel 3D-layered electrochemical-thermal coupled model strategy for the nail-penetration process simulation.《Journal of Power Sources》.2017,第342卷 *
黄文才 等.锂离子电池的高温热失控模拟.《电池》.2019,第49卷(第3期), *
黄文才 等.锂离子电池高温热模拟及热行为.《电池》.2018,第48卷(第6期), *

Also Published As

Publication number Publication date
CN110991049A (zh) 2020-04-10

Similar Documents

Publication Publication Date Title
CN110991049B (zh) 一种基于过充锂离子电池的热失控仿真方法
CN110968968A (zh) 一种基于针刺内部短路锂离子电池的热失控仿真方法
Xie et al. Challenges and opportunities toward fast-charging of lithium-ion batteries
CN107145628B (zh) 基于电化学-热耦合模型的预测锂电池循环寿命的方法
Leng et al. Hierarchical degradation processes in lithium-ion batteries during ageing
CN113705022B (zh) 一种锂离子电池低温充电性能的仿真方法及系统
CN106154172A (zh) 锂离子动力电池内短路程度的定量估算方法
CN108037463A (zh) 一种锂离子电池寿命预测方法
Li et al. Swelling force in lithium-ion power batteries
CN109991301A (zh) 一种基于电化学-力耦合模型的锂离子电池放电过程中电极颗粒间应力的确定方法
Sebastian et al. Adaptive fast charging methodology for commercial Li‐ion batteries based on the internal resistance spectrum
Gao et al. Influence of coupling of overcharge state and short-term cycle on the mechanical integrity behavior of 18650 Li-ion batteries subject to lateral compression
Nam et al. The characterization of dynamic behavior of Li-ion battery packs for enhanced design and states identification
CN109459463A (zh) 一种锂离子电池正极材料热存储稳定性的快捷评价方法
Balagopal et al. Effect of calendar ageing on SEI growth and its impact on electrical circuit model parameters in Lithium ion batteries
CN114996932A (zh) 一种锂离子电池过充热失控模型的建模方法
Zhang et al. Experimental investigation on the essential cause of the degrading performances for an overcharging ternary battery
Chen et al. Benign-to-malignant transition in external short circuiting of lithium-ion batteries
Leng et al. The effect of temperature on the electrochemistry in Lithium-ion batteries
Qiu et al. Understanding thermal and mechanical effects on lithium plating in lithium-ion batteries
Zhao et al. Cathode Chemistries and Electrode Parameters Affecting the Fast Charging Performance of Li-Ion Batteries
Yuan et al. Quantification of electrochemical-mechanical coupling in lithium-ion batteries
Soltani et al. Degradation behaviour analysis and end-of-life prediction of lithium titanate oxide batteries
CN110580389A (zh) 一种电池极耳的优化方法
Roy et al. Battery Health Management-A Perspective of Design, Optimization, Manufacturing, Fault Detection, and Recycling

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant