CN112485674B - 一种正向锂离子电池内短路热失控建模方法 - Google Patents

一种正向锂离子电池内短路热失控建模方法 Download PDF

Info

Publication number
CN112485674B
CN112485674B CN202011311126.2A CN202011311126A CN112485674B CN 112485674 B CN112485674 B CN 112485674B CN 202011311126 A CN202011311126 A CN 202011311126A CN 112485674 B CN112485674 B CN 112485674B
Authority
CN
China
Prior art keywords
short circuit
battery
lithium ion
model
internal short
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
CN202011311126.2A
Other languages
English (en)
Other versions
CN112485674A (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.)
Tsinghua University
Original Assignee
Tsinghua 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 Tsinghua University filed Critical Tsinghua University
Priority to CN202011311126.2A priority Critical patent/CN112485674B/zh
Publication of CN112485674A publication Critical patent/CN112485674A/zh
Application granted granted Critical
Publication of CN112485674B publication Critical patent/CN112485674B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/367Software therefor, e.g. for battery testing using modelling or look-up tables
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K13/00Thermometers specially adapted for specific purposes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/385Arrangements for measuring battery or accumulator variables
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/389Measuring internal impedance, internal conductance or related variables
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/50Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
    • G01R31/52Testing for short-circuits, leakage current or ground faults
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Secondary Cells (AREA)
  • Battery Electrode And Active Subsutance (AREA)

Abstract

本发明提出一种正向锂离子电池内短路热失控建模方法,属于电池技术领域技术领域。该方法首先获取锂离子电池的质量和关键几何参数,建立该电池的几何模型;然后再利用四个相同的锂离子电池制造四种内短路损伤,进而得到每种内短路损伤的内阻;在几何模型的基础上建立锂离子电池的三维热模型,并在电模型基础上考虑内短路过程建立内短路电模型;三维热模型和内短路电模型组成该电池的内短路热失控模型。本发明具有可以高效、准确、安全建立锂离子电池内短路模型的优点,用于探究内短路导致热失控的边界条件,指导电池与电池包设计,有效缩短开发时长。

Description

一种正向锂离子电池内短路热失控建模方法
技术领域
本发明属于电池技术领域技术领域,特别提出一种正向锂离子电池内短路热失控建模方法。
背景技术
锂离子电池内短路是重要的安全问题。内短路是电池内部正极和负极直接或间接接触后,发生的放电并伴随放热的现象,是锂离子电池热失控的共性环节,因此建立锂离子电池内短路导致热失控的模型十分关键。该模型可以用于研究内短路导致热失控的边界条件,进而指导电池的设计和电池包的设计。
目前,常用的锂离子电池内短路建模方法,主要基于等效电阻内短路实验,获取电池内短路的热特征与电特征,该方法由于外电路产热、短路面积过大等原因,偏离真实内短路程度较高,很难在内短路后触发电池的热失控事故。同时,该方法没有触发真实的内短路,缺少对于电池真实内短路损伤程度的测量。因此,所以建立的锂离子电池内短路模型,无法真实预测电池内短路后造成电池热失控的行为。
发明内容
本发明的目的是为克服已有技术的不足之处,提出一种正向锂离子电池内短路热失控建模方法。本发明具有可以高效、准确、安全建立锂离子电池内短路模型的优点,用于探究内短路导致热失控的边界条件,指导电池与电池包设计,有效缩短开发时长。
本发明提出一种正向锂离子电池内短路热失控建模方法,其特征在于,该方法包括以下步骤:
1)任意选定一款待建模的锂离子电池,测量该电池的质量m和关键几何参数;
基于关键几何参数建立该电池的三维几何模型,计算该电池电芯的体积V,计算该电池的密度ρ=m/V;
2)获取四个与步骤1)选定的电池型号相同且全新的锂离子电池,将该四个锂离子电池均放电至0V并拆解,然后将得到的四个锂离子电池卷芯静置、风干,至电解液完全挥发后,分别在每个卷芯内部制造一种内短路损伤,得到分别对应四种内短路损伤的四个卷芯;其中,所述四种内短路损伤分别为:正极-负极内短路损伤、正极-铜内短路损伤、铝-负极内短路损伤和铝-铜内短路损伤;
3)将步骤2)得到的每个卷芯夹紧,并使用直流内阻仪器连接该卷芯正负极并测量电阻,得到每个卷芯的内短路内阻,分别记为RISC1、RISC2、RISC3、RISC4,分别代表四种内短路损伤的内阻;
4)获取一个与步骤1)选定的电池型号相同且全新的锂离子电池,对该锂离子电池进行绝热热失控测试,记录测试过程中的该电池不同时刻的温度T(t)构成的温度时间曲线和不同时刻的电压V(t)构成的电压时间曲线;
对温度时间曲线进行分析,计算温度时间曲线的微分,获得dT/dt(t)构成的温升速率曲线,根据温升速率曲线,获得该电池热失控起始温度T1;根据温度时间曲线获得热失控最高温度T2;根据电压时间曲线获得电池隔膜失效温度Ts
进一步,计算该电池所包含的热失控化学能:
HTR=CP×m×(T2-T1)
其中,CP是电池的比热容;
5)获取一个与步骤1)选定的电池型号相同且全新的锂离子电池,对该电池进行常温充放电测试、混合动力脉冲能力特性测试,并建立该锂离子电池的电模型,计算该电池使用过程中的端电压Ut与电流I,
6)在步骤1)建立的几何模型基础上,建立锂离子电池的三维热模型,其中,三维热模型的电池温度T分布的控制方程为:
Figure BDA0002789820980000021
其中,λi为电池的导热系数,下标i={x,y,z}分别代表xyz三个导热方向,Q为电池的生热速率;
Q=Qrev+Qirr+QISC+QTR
其中,Qrev为电池可逆产热,Qirr为电池不可逆产热,QISC为内短路产热,QTR为热失控产热;
Figure BDA0002789820980000022
其中I为电池电流,Tavg为电池平均温度,
Figure BDA0002789820980000023
为熵系数;
Qirr=I2R0
其中,R0为电池的欧姆内阻;
Figure BDA0002789820980000024
其中,RISC为电池短路内阻,RISC=ξISCRISCj,RISCj由步骤3)获得,下标j={1,2,3,4}分别代表步骤3)的四种内短路损伤;ξISC为内短路系数,
Figure BDA0002789820980000031
Figure BDA0002789820980000032
其中,CTR为热失控速率,
Figure BDA0002789820980000033
CTR为热失控速率参数;
该模型中,边界条件为
Figure BDA0002789820980000034
即为绝热边界条件;
7)在步骤5)建立的电模型基础上将短路内阻串联连接至电模型中的端电压部分,建立内短路电模型,其中端电压Ut与电流I的关系为Ut=IRISC
8)锂离子电池内短路热失控模型建立完毕,该模型包括步骤6)建立的三维热模型和步骤7)建立的内短路电模型;
将设定的模型计算步长Δt,初始荷电状态SOC,电池初始温度T0,以及初始短路内阻RISC0输入锂离子电池内短路热失控模型,该模型输出不同时刻的端电压Ut(t)构成的端电压时间曲线及电池不同位置在不同时刻的温度T(t)构成的的温度时间曲线T(t)。
本发明的特点及有益效果:
本发明建立了考虑热失控的电-热耦合内短路模型,可以同时模拟电池内短路后的电行为与热行为,进一步判断内短路是否会造成电池的热失控,并模拟电池热失控后的行为,建立了电池内短路全过程模型,更全面仿真模拟电池的安全性能,用于研究内短路导致热失控的边界条件,进而指导电池的设计和电池包的设计。该方法明确了四种内短路损伤的制造及定量测试方法,并应用四种内短路损伤对电池的内短路安全性能进行判断。建模中,考虑了随着电池温度升高,隔膜融化导致短路内阻下降的现象,使得对于内短路后的电压行为预测更为准确,进一步提高了所获得的热失控边界的精度。该方法无需进行复杂的、危险的内短路模拟实验,即可建立较为准确的内短路模型,有效提高了建模的效率与安全性,进一步缩短电池设计、电池包设计与管理方法开发的周期。
附图说明
图1是本发明的一种正向锂离子电池内短路热失控建模方法的整体流程图;
图2是本实施例建立的软包电池几何模型,该电池由主体、极耳、铝塑膜外壳组成;
图3是本实施例中制作第一种内短路损伤电池的方法,即正极-负极内短路损伤;
图4是本实施例中制作第二种内短路损伤电池的方法,即正极-铜内短路损伤;
图5是本实施例中制作第三种内短路损伤电池的方法,即铝-负极内短路损伤;
图6是本实施例中制作第四种内短路损伤电池的方法,即铝-铜内短路损伤;
图7是本实施例中的电池绝热热失控测试结果,包括温度时间曲线、电压时间曲线、电压温度曲线及温升速度温度曲线;
图8是本实施例使用的一阶RC等效电路模型;
图9是本实施例使用的一阶RC等效电路模型中的参数;
图10是本实施例使用的一阶RC等效电路内短路模型;
图11是使用本发明建立的模型,预测的四种内短路损伤下电池的电压时间曲线与电池最高温度时间曲线;
图12是使用本发明建立的模型,研究电池内短路热失控边界条件的结果。
具体实施方式
本发明提出一种正向锂离子电池内短路热失控建模方法,下面结合附图及具体实施例详细说明如下:
本发明提出一种正向锂离子电池内短路热失控建模方法,整体流程如图1所示,包括以下步骤:
1)任意选定一款待建模的锂离子电池(型号无特殊要求),测量该电池的质量m和关键几何参数;
对于软包及方形电池,其关键几何参数包括电池的长度L、厚度B及高度H;对于圆柱形电池,其关键几何参数包括电池的高度H和直径D;
基于关键几何参数建立该电池的三维几何模型,并计算该电池电芯的体积V,计算表达式为:对于软包及方形电池V=L×B×H,对于圆柱形电池V=0.25π×D2×H,进而计算该电池的密度ρ=m/V。
2)获取四个与步骤1)选定的电池型号相同且全新的锂离子电池,将该四个锂离子电池均放电至0V并拆解,然后将得到的四个锂离子电池卷芯静置、风干,至电解液完全挥发后,分别在每个卷芯内部制造一种内短路损伤,得到分别对应四种内短路损伤的四个卷芯。
所述四种内短路损伤制造方法分别为:
将第一个锂离子电池卷芯拆解,在隔膜中心位置切割一个直径为2~4mm的圆形孔洞,保证该孔洞两侧材料分别是正极活性材料和负极活性材料,将卷芯重新组装并复原,得到了正极-负极内短路损伤卷芯。
将第二个锂离子电池卷芯拆解,在隔膜中心位置切割一个直径为2~4mm的圆形孔洞,保证该孔洞两侧材料分别是正极活性材料和负极活性材料,使用刀轻轻刮除或者使用少量N-甲基吡咯烷酮擦拭掉对应圆形孔洞位置的负极活性材料,使对应该位置的铜集流体露出,将卷芯重新组装并复原,得到了正极-铜内短路损伤卷芯。
将第三个锂离子电池卷芯拆解,在隔膜中心位置切割一个直径为2~4mm的圆形孔洞,保证该孔洞两侧材料分别是正极活性材料和负极活性材料,使用刀轻轻刮除或者使用少量N-甲基吡咯烷酮擦拭掉对应圆形孔洞位置的正极活性材料,使对应该位置的铝集流体露出,将卷芯重新组装并复原,得到了铝-负极内短路损伤卷芯。
将第四个锂离子电池卷芯拆解,在隔膜中心位置切割一个直径为2~4mm的圆形孔洞,保证该孔洞两侧材料分别是正极活性材料和负极活性材料,使用刀轻轻刮除或者使用少量N-甲基吡咯烷酮擦拭掉对应圆形孔洞位置的正极活性材料及负极活性材料,使对应该位置的铝集流体和铜集流体分别露出,将卷芯重新组装并复原,得到了铝-铜内短路损伤卷芯。
3)将步骤2)得到的每个卷芯夹紧,并使用直流内阻仪器连接该卷芯正负极并测量电阻,得到每个卷芯的内短路内阻,分别记为RISC1、RISC2、RISC3、RISC4,分别代表四种内短路损伤的内阻。
4)获取一个与步骤1)选定的电池型号相同且全新的锂离子电池,对该锂离子电池进行绝热热失控测试,记录测试过程中的该电池不同时刻的温度T(t)构成的温度时间曲线和不同时刻的电压V(t)构成的电压时间曲线。
对温度时间曲线进行分析,计算温度时间曲线的微分,获得dT/dt(t)构成的温升速率曲线,根据温升速率曲线,获得该电池热失控起始温度T1(温度上升速率dT/dt首次大于1℃/s的时刻对应的电池温度)。根据温度时间曲线获得热失控最高温度T2(温度曲线中的温度最大值)。根据电压时间曲线获得电池隔膜失效温度Ts(电池电压曲线发生突降的时刻对应的电池温度)。
进一步,计算该电池所包含的热失控化学能HTR=CP×m×(T2-T1),m是步骤1)测量的电池质量;CP是电池的比热容,可以通过测量得到,也可以采用文献值。
5)获取一个与步骤1)选定的电池型号相同且全新的锂离子电池,对该电池进行常温充放电测试、混合动力脉冲能力特性测试,并建立该锂离子电池的电模型,计算该电池使用过程中的端电压Ut与电流I,此处锂离子电池的电模型可以采用等效电路模型、电化学模型等。
本实施例中电模型采用一阶RC等效电路模型,表达式如下:
U1=IR1×[1-exp(-t/τ1)]
Ut=OCV(SOC)-IR0-U1
其中,R0为该锂离子电池的欧姆阻抗,R1为极化阻抗,τ1为时间常数,OCV为相应SOC下电池的开路电压,这些参数都可以通过混合动力脉冲能力特性测试、充放电测试辨识得到。
6)在步骤1)建立的几何模型基础上,建立锂离子电池的三维热模型,其中,三维热模型的电池温度T分布的控制方程为:
Figure BDA0002789820980000061
其中ρ为电池密度,由步骤1)计算得到,λi为电池的导热系数,下标i={x,y,z}代表不同的导热方向,可通过测试得到,也可以采用文献值,Q为电池的生热速率。
Q=Qrev+Qirr+QISC+QTR,由四部分组成,Qrev为电池可逆产热,Qirr为电池不可逆产热,QISC为内短路产热,QTR为热失控产热。其中,前三项的产热与电流有关,电流在电模型中与内短路阻值有关,因此前三项都与内短路损伤程度有关,不是确定的。第四项与内短路损伤无关。
Figure BDA0002789820980000062
其中I为电池电流,Tavg为电池平均温度,通过计算电池三维温度场平均值得到,
Figure BDA0002789820980000063
为熵系数,可通过实验测量,也可以采用文献值。
Qirr=I2R0,R0为电池的欧姆内阻。
Figure BDA0002789820980000064
RISC为电池短路内阻,RISC=ξISCRISCj,RISCj由步骤3)获得,下标j={1,2,3,4}分别代表步骤3)获得的四种内短路损伤,ξISC为内短路系数,反应电池当前隔膜融化情况,
Figure BDA0002789820980000065
其中Ts为步骤4)获得的电池隔膜失效温度。
Figure BDA0002789820980000066
其中V为步骤1)计算的电池体积,HTR为步骤4)计算的电池热失控化学能,CTR为热失控速率,
Figure BDA0002789820980000067
CTR为热失控速率参数,一般采用文献值5~12,T1为步骤4获得的该电池热失控起始温度。
该模型中,边界条件为
Figure BDA0002789820980000068
即为绝热边界条件。
7)在步骤5)建立的电模型基础上将短路内阻串联连接至电模型中的端电压部分,建立内短路电模型,其中端电压Ut与电流I的关系为Ut=IRISC,RISC为步骤6)所计算的电池短路内阻。
8)锂离子电池内短路热失控模型建立完毕,该模型包括步骤6)建立的三维热模型和步骤7)建立的内短路电模型;
锂离子电池内短路热失控模型的使用方法为:设定模型计算步长Δt,初始荷电状态SOC,电池初始温度T0,以及初始短路内阻RISC0。上述四个参数输入步骤7)建立的内短路电模型中,计算当前时刻的端电压Ut与电流I,并将模型计算步长Δt,初始荷电状态SOC,电池初始温度T0,以及初始短路内阻RISC0和内短路模型中计算得到的当前时刻的端电压Ut与电流I输入至步骤6)建立的三维热模型,计算当前时刻下的温度分布T及短路内阻RISC,然后将T和RISC输入至步骤7)建立的内短路模型,实现循环计算,直至内短路电模型计算得到的电池的端电压Ut小于5mV,计算结束,输出不同时刻的端电压Ut(t)构成的端电压时间曲线及电池不同位置在不同时刻的温度T(t)构成的的温度时间曲线T(t)。
本实施例中,设定SOC=100%,T0=25℃,RISC0=RISCj,下标j={1,2,3,4}代表步骤3)获得的四种内短路损伤,将上述参数输入锂离子电池内短路热失控模型,该模型输出四种内短路损伤下的电池的电压时间曲线Ut(t)与温度时间曲线T(t),用于描述内短路后的电池行为,其中每种内短路损伤下温度时间曲线的最大值为电池最高温度Tmax,将该最高温度作为对应内短路损伤的内短路安全性的评价指标。
下面结合一个实施例对本发明进一步详细说明如下:
本实施例提出一种正向锂离子电池内短路热失控建模方法,包括以下步骤:
1)选定一款1Ah软包锂离子电池,测量电池的质量m=25g和电池主体部分的关键几何参数,其关键几何参数包括电池的长度L=0.048m、厚度B=0.007m及高度H=0.058m,基于所测量参数建立该电池的三维几何模型,如图2所示,该模型中可以忽略电池极耳等飞主体部件,仅考虑卷芯所在的电池主体,并在电池中央建立直径为4mm的圆形内短路区域,并计算该电芯的体积为V=1.95×10-5m3,及密度S=1302kg/m3
2)获取四个与步骤1)选定的电池型号相同且全新的锂离子电池,将该四个锂离子电池均放电至0V并拆解,然后将得到的四个锂离子电池卷芯静置、风干,至电解液完全挥发后,分别在每个卷芯内部制造一种内短路损伤,共四种内短路损伤卷芯。
四种内短路损伤制造方法分别为:
将第一个锂离子电池卷芯拆解,在隔膜中心位置切割一个直径为4mm的圆形孔洞,保证该孔洞两侧材料分别是正极活性材料和负极活性材料,将卷芯重新组装并复原,得到了正极-负极内短路损伤卷芯,如图3所示。
将第二个锂离子电池卷芯拆解,在隔膜中心位置切割一个直径为4mm的圆形孔洞,保证该孔洞两侧材料分别是正极活性材料和负极活性材料,使用刀轻轻刮除或者使用少量N-甲基吡咯烷酮擦拭掉对应圆形孔洞位置的负极活性材料,使对应该位置的铜集流体露出,将卷芯重新组装并复原,得到了正极-铜内短路损伤卷芯,如图4所示。
将第三个锂离子电池卷芯拆解,在隔膜中心位置切割一个直径为4mm的圆形孔洞,保证该孔洞两侧材料分别是正极活性材料和负极活性材料,使用刀轻轻刮除或者使用少量N-甲基吡咯烷酮擦拭掉对应圆形孔洞位置的正极活性材料,使对应该位置的铝集流体露出,将卷芯重新组装并复原,得到了铝-负极内短路损伤卷芯,如图5所示。
将第四个锂离子电池卷芯拆解,在隔膜中心位置切割一个直径为4mm的圆形孔洞,保证该孔洞两侧材料分别是正极活性材料和负极活性材料,使用刀轻轻刮除或者使用少量N-甲基吡咯烷酮擦拭掉对应圆形孔洞位置的正极活性材料及负极活性材料,使对应该位置的铝集流体和铜集流体分别露出,将卷芯重新组装并复原,得到了铝-铜内短路损伤卷芯,如图6所示。
3)将步骤2)得到的每个卷芯夹紧,并使用直流内阻仪器连接该卷芯正负极,测量电阻,得到四个卷芯的内短路内阻,分别记为RISC1、RISC2、RISC3、RISC4,分别对应四种内短路损伤的内阻。
本实施例测量得到的四种内短路损伤的内阻值如下表所示:
表1本实施例的四种内短路损伤的内阻值表
符号 案例测试结果/Ω 一般范围/Ω
R<sub>ISC1</sub> 163.0 >100
R<sub>ISC2</sub> 49.40 30~100
R<sub>ISC3</sub> 0.015 <0.05
R<sub>ISC4</sub> 3.069 0.5~10
4)获取一个与步骤1)选定的电池型号相同且全新的锂离子电池,对该锂离子电池进行绝热热失控测试,记录测试过程中的不同时刻的温度T(t)构成的温度时间曲线和不同时刻的电压V(t)构成的电压时间曲线,如图7所示。
对温度时间曲线进行分析,计算温度时间曲线的微分,获得dT/dt构成的温升速率曲线,根据该曲线,获得该电池热失控起始温度T1为169.1℃(温度上升速率dT/dt首次大于1℃/s的时刻对应的电池温度)。根据温度曲线获得热失控最高温度T2为830.5℃(温度曲线中的温度最大值)。根据电压曲线获得电池隔膜失效温度Ts为162.9℃(电池电压曲线发生突降的时刻对应的电池温度)。
进一步,计算该电池所包含的热失控化学能HTR=18188.5J,CP是电池的比热容,采用文献值1100J/kg·K。
5)获取一个与步骤1)选定的电池型号相同且全新的锂离子电池,对该电池进行常温充放电测试、混合动力脉冲能力特性测试,并建立该锂离子电池的电模型,计算该电池使用过程中的端电压Ut与电流I,此处可以采用等效电路模型、电化学模型等。
本实施例中电模型可采用一阶RC等效电路模型,如图8所示。
U1=IR1×[1-exp(-t/τ1)]
Ut=OCV(SOC)-IR0-U1
其中,R0为该锂离子电池的欧姆阻抗,R1为极化阻抗,τ1为时间常数,OCV为相应SOC下电池的开路电压,这些参数都可以通过混合动力脉冲能力特性测试、充放电测试辨识得到,如图9所示,其中图9a)为通过测试数据及模型辨识得到的极化阻抗R1-SOC曲线,图9b)为通过测试数据及模型辨识得到的欧姆阻抗R0-SOC曲线,图9c)为通过测试数据及模型辨识得到的时间常数τ1-SOC曲线,图9d)为通过测试数据及模型辨识得到的开路电压OCV-SOC曲线。
6)在步骤1)建立的几何模型基础上,建立锂离子电池的三维热模型,电池温度T分布的控制方程为:
Figure BDA0002789820980000091
其中ρ为电池密度,由步骤1)计算得到,λi为电池的导热系数,下标i={x,y,z}代表不同的导热方向,可通过测试得到,也可以采用文献值,Q为电池的生热速率。
Q=Qrev+Qirr+QISC+QTR,由四部分组成,Qrev为电池可逆产热,Qirr为电池不可逆产热,QISC为内短路产热,QTR为热失控产热。
Figure BDA0002789820980000092
其中I为电池电流,Tavg为电池平均温度,通过计算三维温度场平均值得到,
Figure BDA0002789820980000093
为熵系数,可通过实验测量,也可以采用文献值。
Qirr=I2R0,R0为电池的欧姆内阻。
Figure BDA0002789820980000094
RISC为电池短路内阻,RISC=ξISCRISCj,RISCj由步骤3获得,下标j={1,2,3,4}分别代表步骤3)获得的四种内短路损伤,ξISC为内短路系数,反应电池当前隔膜融化情况,
Figure BDA0002789820980000095
其中Ts为步骤4获得的电池隔膜失效温度。
Figure BDA0002789820980000096
其中V为步骤1计算的电池体积,UTR为步骤4)计算的电池热失控化学能,CTR为热失控速率,
Figure BDA0002789820980000097
CTR为热失控速率参数,一般采用文献值5~12,T1为步骤4)获得的该电池热失控起始温度。
该模型中,边界条件为
Figure BDA0002789820980000098
即为绝热边界条件。
7)在步骤5)建立的电模型基础上将短路内阻串联连接至电模型中的端电压部分建立内短路电模型,如图10所示,其中端电压Ut与电流I的关系为Ut=IRISC,RISC为电池短路内阻为步骤6)所计算的短路内阻。
8)锂离子电池内短路热失控模型建立完毕,该模型包括步骤6)建立的三维热模型和步骤7)建立的内短路电模型;
该模型使用方法为:设定模型计算步长Δt,初始荷电状态SOC,电池初始温度T0,以及初始短路内阻RISC0。上述四个参数输入步骤7)建立的内短路电模型中,计算当前时刻的端电压Ut与电流I,并将模型计算步长Δt,初始荷电状态SOC,电池初始温度T0,以及初始短路内阻RISC0和内短路模型中计算得到的当前时刻的端电压Ut与电流I输入至步骤6)建立的三维热模型,计算当前时刻下的温度分布T及短路内阻RISC,然后将T和RISC输入至步骤7)建立的内短路模型,实现循环计算,直至内短路电模型计算得到的电池的端电压Ut小于5mV,计算结束,输出不同时刻的端电压Ut(t)构成的端电压时间曲线及电池不同位置在不同时刻的温度T(t)构成的的温度时间曲线T(t)。
本实施例中,设定Δt=10ms,SOC=100%,T0=25℃,RISC0=RISCj,下标j={1,2,3,4}代表步骤3)获得的四种内短路损伤,将上述参数输入锂离子电池内短路热失控模型,该模型输出四种内短路损伤下的电池的端电压时间曲线与温度时间曲线,如图11所示,其中图11a)为第一种内短路损伤下的电池端电压时间曲线与第二种内短路损伤下的电池端电压时间曲线,图11b)为第一种内短路损伤下的电池最高温度时间曲线与第二种内短路损伤下的电池最高温度时间曲线,第一种内短路损伤为正极-负极类型内短路,第二种内短路损伤为正极-铜类型内短路,其内短路阻值较大,电压下降缓慢,且没有明显的温度升高;图11c)为第三种内短路损伤下的电池端电压时间曲线与第四种内短路损伤下的电池端电压时间曲线,图11d)为第三种内短路损伤下的电池最高温度时间曲线与第四种内短路损伤下的电池最高温度时间曲线,第三种内短路损伤为铝-负极类型内短路,第四种内短路损伤为铝-铜类型内短路,铝-负极内短路发生后,电池释放出较多电能且在局部区域累积出发了电池的热失控,温度不可控上升至1300℃左右,铝-铜类型内短路后发生后电池电压瞬间失效,电池最高温度上升至150℃左右,但未触发电池热失控。
选取SOC={0,10,20,30,40,50,60,70,80,90,100}%,选取RISC0={0.05,0.1,0.15,0.2,…,1.95,2}Ω,对于每一个SOC取值,将其与每一个RISC0取值进行组合输入锂离子电池内短路热失控模型,即将11*40组SOC、RISC0取值组合输入锂离子电池内短路热失控模型,根据的输出结果,可计算得到该电池的内短路-热失控边界,如图12所示。

Claims (4)

1.一种正向锂离子电池内短路热失控建模方法,其特征在于,该方法包括以下步骤:
1)任意选定一款待建模的锂离子电池,测量该电池的质量m和关键几何参数;
基于关键几何参数建立该电池的三维几何模型,计算该电池电芯的体积V,计算该电池的密度ρ=m/V;
2)获取四个与步骤1)选定的电池型号相同且全新的锂离子电池,将步骤2)获取的四个锂离子电池均放电至0V并拆解,然后将得到的四个锂离子电池卷芯静置、风干,至电解液完全挥发后,分别在每个卷芯内部制造一种内短路损伤,得到分别对应四种内短路损伤的四个卷芯;其中,所述四种内短路损伤分别为:正极-负极内短路损伤、正极-铜内短路损伤、铝-负极内短路损伤和铝-铜内短路损伤;
3)将步骤2)得到的每个卷芯夹紧,并使用直流内阻仪器连接每个卷芯正负极并测量电阻,得到每个卷芯的内短路内阻,分别记为RISC1、RISC2、RISC3、RISC4,分别代表四种内短路损伤的内阻;
4)获取一个与步骤1)选定的电池型号相同且全新的锂离子电池,对步骤4)获取的锂离子电池进行绝热热失控测试,记录测试过程中的步骤4)获取的锂离子电池不同时刻的温度T(t)构成的温度时间曲线和不同时刻的电压V(t)构成的电压时间曲线;
对温度时间曲线进行分析,计算温度时间曲线的微分,获得dT/dt(t)构成的温升速率曲线,根据温升速率曲线,获得步骤4)获取的锂离子电池热失控起始温度T1;根据温度时间曲线获得热失控最高温度T2;根据电压时间曲线获得电池隔膜失效温度Ts
进一步,计算步骤4)获取的锂离子电池所包含的热失控化学能:
HTR=CP×m×(T2-T1)
其中,CP是电池的比热容;
5)获取一个与步骤1)选定的电池型号相同且全新的锂离子电池,对步骤5)获取的锂离子电池进行常温充放电测试、混合动力脉冲能力特性测试,并建立步骤5)获取的锂离子电池的电模型,计算步骤5)获取的锂离子电池使用过程中的端电压Ut与电流I,
6)在步骤1)建立的三维几何模型基础上,建立锂离子电池的三维热模型,其中,三维热模型的电池温度T分布的控制方程为:
Figure FDA0003304974170000011
其中,λi为电池的导热系数,下标i={x,y,z}分别代表xyz三个导热方向,Q为电池的生热速率;
Q=Qrev+Qirr+QISC+QTR
其中,Qrev为电池可逆产热,Qirr为电池不可逆产热,QISC为内短路产热,QTR为热失控产热;
Figure FDA0003304974170000021
其中I为电池电流,Tavg为电池平均温度,
Figure FDA0003304974170000022
为熵系数;
Qirr=I2R0
其中,R0为电池的欧姆内阻;
Figure FDA0003304974170000023
其中,RISC为电池短路内阻,RISC=ξISCRISCj,RISCj由步骤3)获得,下标j={1,2,3,4}分别代表步骤3)的四种内短路损伤;ξISC为内短路系数,
Figure FDA0003304974170000024
Figure FDA0003304974170000025
其中,cTR为热失控速率,
Figure FDA0003304974170000026
CTR为热失控速率参数;
该模型中,边界条件为
Figure FDA0003304974170000027
即为绝热边界条件;
7)在步骤5)建立的电模型基础上将内短路内阻串联连接至电模型中的端电压部分,建立内短路电模型,其中端电压Ut与电流I的关系为Ut=IRISC
8)锂离子电池内短路热失控模型建立完毕,该模型包括步骤6)建立的三维热模型和步骤7)建立的内短路电模型;
将设定的锂离子电池内短路热失控模型计算步长Δt,初始荷电状态SOC,电池初始温度T0,以及初始短路内阻RISC0输入锂离子电池内短路热失控模型,锂离子电池内短路热失控模型输出不同时刻的端电压Ut(t)构成的端电压时间曲线及电池不同位置在不同时刻的温度T(t)构成的温度时间曲线。
2.如权利要求1所述的方法,其特征在于,所述步骤1)中关键几何参数具体如下:
对于软包及方形电池,关键几何参数包括电池的长度L、厚度B及高度H;对于圆柱形电池,关键几何参数包括电池的高度H和直径D。
3.如权利要求1所述的方法,其特征在于,所述步骤2)中四种内短路损伤制造方法分别为:
将第一个锂离子电池卷芯拆解,在隔膜中心位置切割一个直径为2~4mm的圆形孔洞,使得该孔洞两侧材料分别是正极活性材料和负极活性材料,将卷芯重新组装并复原,得到正极-负极内短路损伤卷芯;
将第二个锂离子电池卷芯拆解,在隔膜中心位置切割一个直径为2~4mm的圆形孔洞,使得该孔洞两侧材料分别是正极活性材料和负极活性材料,使用刀轻轻刮除或者使用少量N-甲基吡咯烷酮擦拭掉对应圆形孔洞位置的负极活性材料,使对应该位置的铜集流体露出,将卷芯重新组装并复原,得到正极-铜内短路损伤卷芯;
将第三个锂离子电池卷芯拆解,在隔膜中心位置切割一个直径为2~4mm的圆形孔洞,使得该孔洞两侧材料分别是正极活性材料和负极活性材料,使用刀轻轻刮除或者使用少量N-甲基吡咯烷酮擦拭掉对应圆形孔洞位置的正极活性材料,使对应该位置的铝集流体露出,将卷芯重新组装并复原,得到铝-负极内短路损伤卷芯;
将第四个锂离子电池卷芯拆解,在隔膜中心位置切割一个直径为2~4mm的圆形孔洞,使得该孔洞两侧材料分别是正极活性材料和负极活性材料,使用刀轻轻刮除或者使用少量N-甲基吡咯烷酮擦拭掉对应圆形孔洞位置的正极活性材料及负极活性材料,使对应该位置的铝集流体和铜集流体分别露出,将卷芯重新组装并复原,得到铝-铜内短路损伤卷芯。
4.如权利要求1所述的方法,其特征在于,所述步骤5)中电模型采用一阶RC等效电路模型,表达式如下:
U1=IR1×[1-exp(-t/τ1)]
Ut=OCV(SOC)-IR0-U1
其中,R1为极化阻抗,τ1为时间常数,OCV为相应SOC下电池的开路电压。
CN202011311126.2A 2020-11-20 2020-11-20 一种正向锂离子电池内短路热失控建模方法 Active CN112485674B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011311126.2A CN112485674B (zh) 2020-11-20 2020-11-20 一种正向锂离子电池内短路热失控建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011311126.2A CN112485674B (zh) 2020-11-20 2020-11-20 一种正向锂离子电池内短路热失控建模方法

Publications (2)

Publication Number Publication Date
CN112485674A CN112485674A (zh) 2021-03-12
CN112485674B true CN112485674B (zh) 2021-12-10

Family

ID=74932427

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011311126.2A Active CN112485674B (zh) 2020-11-20 2020-11-20 一种正向锂离子电池内短路热失控建模方法

Country Status (1)

Country Link
CN (1) CN112485674B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113484760B (zh) * 2021-07-12 2022-06-24 杭州华塑科技股份有限公司 一种电池热失控识别方法、装置、设备及存储介质
CN114091404B (zh) * 2022-01-20 2022-05-17 天津力神电池股份有限公司 一种锂离子电池的电-热模型
CN114509678B (zh) * 2022-02-10 2024-06-18 北京理工大学重庆创新中心 一种基于大数据的动力电池热失控预测方法
CN115544839A (zh) * 2022-08-05 2022-12-30 天目湖先进储能技术研究院有限公司 锂电池电化学-热-力-短路-热失控耦合模型构建方法
CN116562107B (zh) * 2023-07-06 2023-09-05 中汽研(天津)汽车工程研究院有限公司 动力电池系统热扩散仿真方法、装置及介质

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104035048A (zh) * 2014-06-20 2014-09-10 上海出入境检验检疫局工业品与原材料检测技术中心 一种锂离子电池过充安全性能的热电检测方法及其装置
CN104346524A (zh) * 2014-09-16 2015-02-11 清华大学 一种锂离子电池热失控的建模方法
CN105510843A (zh) * 2016-01-08 2016-04-20 北京北交新能科技有限公司 一种锂离子电池的低温频域电-热模型的构建方法
CN107192914A (zh) * 2017-04-18 2017-09-22 宁德时代新能源科技股份有限公司 锂离子动力电池内短路检测方法
CN108646186A (zh) * 2018-04-19 2018-10-12 浙江大学城市学院 一种基于电化学特征的单体锂离子电池三维热模型
CN109521367A (zh) * 2018-10-31 2019-03-26 华南理工大学 一种锂离子动力电池过充/过放检测系统
CN110376530A (zh) * 2019-08-08 2019-10-25 清华大学 电池内短路检测装置及方法
CN110968968A (zh) * 2019-12-05 2020-04-07 西南交通大学 一种基于针刺内部短路锂离子电池的热失控仿真方法
CN111199106A (zh) * 2020-01-06 2020-05-26 上海空间电源研究所 一种电池绝热热失控过程参数获取方法
CN111597747A (zh) * 2020-05-18 2020-08-28 华南理工大学 一种多点触发三元锂动力电池模组热失控仿真及预测方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9478829B2 (en) * 2013-05-16 2016-10-25 Ec Power, Llc Rechargeable battery with multiple resistance levels
US10330739B2 (en) * 2016-08-24 2019-06-25 The Boeing Company Detecting internal short circuits in batteries
CN110083871B (zh) * 2019-03-27 2021-03-30 北京航空航天大学 一种基于锂离子电池热失控预测模型的热失控模拟方法和装置

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104035048A (zh) * 2014-06-20 2014-09-10 上海出入境检验检疫局工业品与原材料检测技术中心 一种锂离子电池过充安全性能的热电检测方法及其装置
CN104346524A (zh) * 2014-09-16 2015-02-11 清华大学 一种锂离子电池热失控的建模方法
CN105510843A (zh) * 2016-01-08 2016-04-20 北京北交新能科技有限公司 一种锂离子电池的低温频域电-热模型的构建方法
CN107192914A (zh) * 2017-04-18 2017-09-22 宁德时代新能源科技股份有限公司 锂离子动力电池内短路检测方法
CN108646186A (zh) * 2018-04-19 2018-10-12 浙江大学城市学院 一种基于电化学特征的单体锂离子电池三维热模型
CN109521367A (zh) * 2018-10-31 2019-03-26 华南理工大学 一种锂离子动力电池过充/过放检测系统
CN110376530A (zh) * 2019-08-08 2019-10-25 清华大学 电池内短路检测装置及方法
CN110968968A (zh) * 2019-12-05 2020-04-07 西南交通大学 一种基于针刺内部短路锂离子电池的热失控仿真方法
CN111199106A (zh) * 2020-01-06 2020-05-26 上海空间电源研究所 一种电池绝热热失控过程参数获取方法
CN111597747A (zh) * 2020-05-18 2020-08-28 华南理工大学 一种多点触发三元锂动力电池模组热失控仿真及预测方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
A study on effect of lithium ion battery design variables upon features of thermal-runaway using mathematical model and simulation;Chan Ho Lee 等;《Journal of Power Sources》;20151231;第29卷;第498-511页 *
Analysis of internal short-circuit in a lithium ion cell;Shriram Santhanagopalan 等;《J.Power Sources.》;20091231;第194卷;第550-557页 *
三元锂离子动力电池针刺热失控实验与建模;张明轩 等;《汽车工程》;20151231;第37卷(第7期);第743-750页 *
循环老化对锂离子电池在绝热条件下的产热及热失控影响;毛亚 等;《储能科学与技术》;20181130;第7卷(第6期);第1120-1127页 *
锂离子电池热失控传播研究进展;杨国明 等;《船电技术》;20151231;第35卷(第9期);第48-51页 *

Also Published As

Publication number Publication date
CN112485674A (zh) 2021-03-12

Similar Documents

Publication Publication Date Title
CN112485674B (zh) 一种正向锂离子电池内短路热失控建模方法
Abaza et al. Experimental study of internal and external short circuits of commercial automotive pouch lithium-ion cells
Liu et al. Safety issues caused by internal short circuits in lithium-ion batteries
Abada et al. Safety focused modeling of lithium-ion batteries: A review
Jaguemont et al. Characterization and modeling of a hybrid-electric-vehicle lithium-ion battery pack at low temperatures
Lu et al. A method of cell-to-cell variation evaluation for battery packs in electric vehicles with charging cloud data
CN111062137B (zh) 一种锂离子电池性能预测模型及其构建方法和应用
Wang et al. Influence of connecting plate resistance upon LiFePO4 battery performance
Xie et al. An improved resistance-based thermal model for prismatic lithium-ion battery charging
CN112470325A (zh) 检测电池中的内部短路(isc)的方法和系统及实现这种isc检测方法的电池电芯
Jin et al. Comparison of Li-ion battery degradation models for system design and control algorithm development
Ma et al. Electro-thermal modeling of a lithium-ion battery system
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
CN108363016B (zh) 基于人工神经网络的电池微短路定量诊断方法
Yazdanpour et al. A circuit-based approach for electro-thermal modeling of lithium-ion batteries
Lazreg et al. Soc estimation of li-ion battery pack for light electric vehicles using enhanced coulomb counting algorithm
Zhang et al. Non-isothermal Ragone plots of Li-ion cells from datasheet and galvanostatic discharge tests
Liao et al. A dynamic equivalent circuit model of LiFePO 4 cathode material for lithium ion batteries on hybrid electric vehicles
CN112198437A (zh) 模拟锂离子电池内短路热失控的刺片、样品电芯及方法
Mo et al. Working temperature effects on mechanical integrity of cylindrical lithium-ion batteries
Watrin et al. Multiphysical lithium-based battery pack modeling for simulation purposes
Wang et al. Development of a one-dimensional thermal-electrochemical model of lithium ion battery
Keyser et al. Thermal characterization of advanced lithium-ion polymer cells
Saqli et al. Battery Pack Thermal Modeling, Simulation and electric model Identification
Ohshima et al. Thermal behavior of small lithium‐ion secondary battery during rapid charge and discharge cycles

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