CN107045105A - 一种锂离子电池组可用能量计算方法 - Google Patents

一种锂离子电池组可用能量计算方法 Download PDF

Info

Publication number
CN107045105A
CN107045105A CN201611076522.5A CN201611076522A CN107045105A CN 107045105 A CN107045105 A CN 107045105A CN 201611076522 A CN201611076522 A CN 201611076522A CN 107045105 A CN107045105 A CN 107045105A
Authority
CN
China
Prior art keywords
lithium
ion battery
capacity
batteries piles
monomer
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.)
Granted
Application number
CN201611076522.5A
Other languages
English (en)
Other versions
CN107045105B (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.)
Beijing Jiaotong University
CRRC Changchun Railway Vehicles Co Ltd
Original Assignee
Beijing Jiaotong University
CRRC Changchun Railway Vehicles Co Ltd
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 Beijing Jiaotong University, CRRC Changchun Railway Vehicles Co Ltd filed Critical Beijing Jiaotong University
Priority to CN201611076522.5A priority Critical patent/CN107045105B/zh
Publication of CN107045105A publication Critical patent/CN107045105A/zh
Application granted granted Critical
Publication of CN107045105B publication Critical patent/CN107045105B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/392Determining battery ageing or deterioration, e.g. state of health

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Battery Mounting, Suspending (AREA)
  • Secondary Cells (AREA)
  • Charge And Discharge Circuits For Batteries Or The Like (AREA)

Abstract

本发明涉及一种锂离子电池组可用能量计算方法。该方法包括如下步骤:S1、检测待测锂离子电池组(均分为n个电池组)中各锂离子电池单体的容量和中值电压,S2、求出任意一个电池组Yj的容量的分布函数;S3、根据分布函数求出Yj的期望值EYj和方差Var(Yj);S4、计算n个电池组的总容量Qsum的临界值Tn;S5、最后求得n个电池组的总能量En即为待测锂离子电池组可用能量。该方法利用对于同一批次的新锂离子电池或锂离子电池使用工况相同时,锂离子电池参数(容量)往往服从正态分布,因此,通过锂离子电池单体参数(容量)分布计算锂离子电池单体成组后的锂离子电池组可用能量,有利于锂离子电池制造厂商或锂离子电池使用者评价不同电压等级下的锂离子电池成组后的表现。

Description

一种锂离子电池组可用能量计算方法
技术领域
本发明涉及新能源车辆的电池管理领域,具体说是一种锂离子电池组可用能量计算方法。
背景技术
近年来,锂离子电池材料和电池制造技术飞速发展。国内外电池厂家不断地生产出大容量、高比能量和高安全性的电池,使得锂离子电池在电动汽车、储能和轨道交通等领域得到广泛的应用。
尽管锂离子电池制造工艺上的进步在一定程度上减小了电池出厂时的差异,但锂离子电池在大规模串联成组后,由于电池的木桶效应,锂离子电池的不一致性使得锂离子电池组的实际可用能量远小于锂离子电池单体的能量。同样,锂离子电池组在使用过程中,电池老化的不一致使得锂离子电池参数的不一致加剧,造成锂离子电池组的能量利用率迅速衰减。
发明内容
针对现有技术中存在的缺陷,本发明的目的在于提供一种锂离子电池组可用能量计算方法,该方法利用对于同一批次的新锂离子电池或锂离子电池使用工况相同时,锂离子电池参数(容量)往往服从正态分布,因此,通过锂离子电池单体参数(容量)分布计算锂离子电池单体成组后的锂离子电池组可用能量,有利于锂离子电池制造厂商或锂离子电池使用者评价不同电压等级下的锂离子电池成组后的表现。
为达到以上目的,本发明采取的技术方案是:
一种锂离子电池组可用能量计算方法,包括如下步骤:
S1、分别检测待测锂离子电池组中N个锂离子电池单体的容量和中值电压,其中,待测锂离子电池组由N个锂离子电池单体串联成组,且平均分n个电池组;
S2、求出任意一个电池组的容量Yj(j=1,2,……n)的分布函数FY(y),y代表锂离子电池单体的容量;
S3、根据Yj的分布函数FY(y)(具体为通过数值模拟的方法)
求出Yj的期望值EYj和方差Var(Yj);
S4、根据式(17)计算n个电池组的总容量Qsum的临界值Tn
式(17)中,φ为标准正态分布的分布函数;α为置信区间;
S5、将Tn代入式(10)中的Qsum,求得En即为待测锂离子电池组可用能量(即待测锂离子电池组最大可用能量的理论值):
En=Uav·Qsum (10);
式(10)中,Uav为待测锂离子电池组中所有锂离子电池单体的中值电压均值。
在上述技术方案基础上,步骤S2中,任意一个电池组的容量Yj(j=1,2,……n)的分布函数FY(y)表示为:
FY(y)=P(Y1≤y)=1-(1-F(y))N/n (11);
Yj的密度函数fY(y)为:
fY(y)=F'Y(y)=N(1-F(y))N/n-1f(y)/n (12);
式(11)和(12)中,
F(y)为锂离子电池单体容量的分布函数;
f(y)为锂离子电池单体容量的概率密度函数:
式(1)中,μ和σ2分别为锂离子电池单体的容量的均值和方差。
本发明的有益效果如下:
1.本发明利用统计学方法,通过获取某一批次锂离子电池的参数分布特性,可计算出在一定置信区间下,任意电压等级下锂离子电池组的可用能量。
2.本发明中所提供的锂离子电池组可用能量的计算方法,并非针对某一特定的已知所有锂离子电池单体参数的电池组进行可用能量计算,而是根据大量锂离子电池单体的参数分布特性,计算在锂离子电池单体随机成组情况下,锂离子电池组最大可用能量的理论值。
3.本发明不仅适用于新的锂离子电池,也适用于梯次利用锂离子电池再成组的锂离子电池组可用能量计算。
4.本发明适用于不同类型的锂离子电池,影响本发明计算结果的是中锂离子电池的容量分布参数,而与锂离子电池的类型无关。
附图说明
本发明有如下附图:
图1为锂离子电池单体的容量测试结果直方图;
图2为n个电池组的总能量En与电池组数n的关系图。
具体实施方式
以下结合附图对本发明作进一步详细说明。
本发明所述的一种锂离子电池组可用能量计算方法,可实现同一类型锂离子电池随机组合情况下的锂离子电池组可用能量的计算。
实施例1
一、首先,针对某一批次同一类型锂离子电池单体,对锂离子电池单体进行容量测试和中值电压测试,确定所有锂离子电池单体的容量和中值电压。
容量测试和中值电压测试的具体测试方法如下:
(1)以0.05C(C表示为电流大小的比率,若锂离子电池额定容量为60Ah,则1C=60A)倍率对锂离子电池单体进行恒流放电,当锂离子电池单体电压到达放电截止电压时,放电停止;
(2)静置2h;
(3)以0.05C倍率对锂离子电池单体进行恒流充电,当锂离子电池单体电压到达充电截止电压时,充电停止,此阶段锂离子电池单体充电容量记为Qc
(4)静置2h;
(5)以0.05C倍率对锂离子电池单体进行恒流放电,当锂离子电池单体电压到达放电截止电压时,放电停止;此阶段锂离子电池单体放电容量记为Qd
(6)静置2h;
(7)以0.05C倍率对锂离子电池单体恒流充电,当充电容量等于Qc/2时停止充电;
(8)静置2h;
(9)测量锂离子电池单体电压V,即为锂离子电池单体的中值电压;
(10)通常情况下,Qc和Qd非常接近,可以认为Qc=Qd,因此,锂离子电池单体的容量为Q=Qc=Qd,锂离子电池类型不同,锂离子电池的充电截止电压和放电截止电压也不同,进行容量测试时,充电截止电压和放电截止电压的参数通常由电池厂提供。
二、锂离子电池组可用能量计算方法的推导
若有N个锂离子电池单体,锂离子电池单体在上述容量测试中的容量为Qi(i=1,…,N),锂离子电池单体容量的均值和方差分别为μ、σ2。若锂离子电池单体容量服从正态分布,则锂离子电池单体容量的概率密度函数f(y)为:
锂离子电池单体容量的分布函数为F(y)。
当N个锂离子电池单体串联成组时,锂离子电池组的容量QN为:
N个锂离子电池单体串联成组时,锂离子电池组的能量EN为:
式(3)中,Umi表示为第i个锂离子电池单体的中值电压。
设将N个锂离子电池单体平均分为n个电池组,则每个电池组的锂离子电池单体数量为N/n个。第j个电池组的容量可表示为:
n个电池组的总能量En可表示为:
式(5)中,Uavj表示为第j个电池组中锂离子电池单体的中值电压平均值。
考虑如下情况,在由大量锂离子电池单体样本中随机选出特定数量的锂离子电池单体串联成组,若串联的锂离子电池单体数量越多,串联电池组中锂离子电池单体中值电压的平均值就会越接近于整个锂离子电池单体样本的中值电压平均值。由于所有的锂离子电池单体服从相同的分布,则每一个电池组当中的Uavj相同。设每一个电池组当中的Uavj均为Uav,因此,n个电池组的总能量可表示为:
设:
因此:
设:
则:
En=Uav·Qsum (10)
需要指出的是,Yj(j=1,…,n)是一组独立同分布的数,因此,Yj的分布函数FY(y)可表示为:
Yj的密度函数fY(y)为:
fY(y)=F'Y(y)=N(1-F(y))N/n-1f(y)/n (12)
设Yj的期望值和方差分别为EYj和Var(Yj),EYj和Var(Yj)数值可由Yj的分布函数FY(y)通过数值模拟的方法求出。
假设在置信区间为α条件下有一临界值Tn使得:
即:
可推导出:
根据中心极限定理可得:
因此,在置信区间为α条件下:
式(17)中,Tn是输出数据,其余参数为输入数据;将根据式(17)求得的Tn代入式(10)中的Qsum,求得En,即为锂离子电池组可用能量。
三、下面针对上述锂离子电池组可用能量计算方法,结合实际数据说明其有益效果。
1、锂离子电池组可用能量计算
按照步骤一的方法对同一批次的95块磷酸铁锂电池单体的容量和中值电压进行测试,容量结果直方图如图1所示,95块磷酸铁锂电池单体的中值电压的平均值Uav为3.3V。
测试过程中,磷酸铁锂电池单体的充电截止电压3.65V,放电截止电压为2.5V。图1的测试结果显示,95块磷酸铁锂电池单体容量均值为49.46Ah,方差为1.499Ah。
根据上述测试结果,假设上述磷酸铁锂电池单体的数量N=100,设将这N个串联的磷酸铁锂电池单体平均分为n个电池组,分别计算在不同的置信区间(α=0.99和α=0.95)下,n个电池组的总能量En,En与电池组数n的关系图如图2所示。
En的具体计算方法如下:
S1、分别检测待测锂离子电池组中N个锂离子电池单体的容量和中值电压,其中,待测锂离子电池组由N个锂离子电池单体串联成组,且平均分n个电池组;
S2、求出任意一个电池组的容量Yj(j=1,2,……n)的分布函数FY(y),y代表锂离子电池单体的容量;
任意一个电池组的容量Yj(j=1,2,……n)的分布函数FY(y)表示为:
FY(y)=P(Y1≤y)=1-(1-F(y))N/n (11);
Yj的密度函数fY(y)为:
fY(y)=F'Y(y)=N(1-F(y))N/n-1f(y)/n (12);
式(11)和(12)中,
F(y)为锂离子电池单体容量的分布函数;
f(y)为锂离子电池单体容量的概率密度函数:
式(1)中,μ和σ2分别为锂离子电池单体的容量的均值和方差。
S3、根据Yj的分布函数FY(y)(具体为通过数值模拟的方法)
求出Yj的期望值EYj和方差Var(Yj);
S4、根据式(17)计算n个电池组的总容量Qsum的临界值Tn
式(17)中,φ为标准正态分布的分布函数;α为置信区间;
S5、将Tn代入式(10)中的Qsum,求得En即为待测锂离子电池组可用能量(即待测锂离子电池组最大可用能量的理论值):
En=Uav·Qsum (10);
式(10)中,Uav为待测锂离子电池组中所有锂离子电池单体的中值电压均值。
2、锂离子电池组可用能量的验证
从以上95块磷酸铁锂电池单体中选取24块磷酸铁锂电池单体,对按照步骤二所述锂离子电池组可用能量的计算方法得到的结果进行验证。24块磷酸铁锂电池单体的容量均值为49.27Ah,容量标准差为1.82Ah。验证实验共分为3组进行,实验1、实验2、实验3(a)和实验3(b)依次进行:
实验1:24块磷酸铁锂电池单体串联成电池组;
实验2:24块磷酸铁锂电池单体随机均分成2组,组内串联成电池组;
实验3:(a)24块磷酸铁锂电池单体随机均分成3组,组内串联成电池组;
(b)24磷酸铁锂电池单体按照容量大小排序,均分为3组,组内串联成电池组。
进行电池组测试时,电池单体均处于满电状态,以0.05C倍率对电池组进行充放电实验,当任意一个电池单体达到充电截止电压或放电截止电压时,充电或放电停止,充电或放电过程之间,电池组静置2h。实验过程中,由电池组测试设备记录电池组放电过程的放电容量和放电能量。实验结果如表1所示。
同时按照上述步骤S1—S5的方法计算获得各组实验的En
以能量利用率作为评价指标,具体计算式为:
式(18)中,En代表实验获得的各电池组的放电能量之和或计算获得的可用能量,Etotal表示电池单体的放电能量之和,在步骤一中(1)-(9)测试电池单体容量的过程中,可同步得到电池单体的放电能量。
测试结果显示,在置信区间为0.95条件下,步骤S1—S5的方法计算获得的能量利用率误差小于1%。
表1电池组实验结果
本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。

Claims (2)

1.一种锂离子电池组可用能量计算方法,其特征在于,包括如下步骤:
S1、分别检测待测锂离子电池组中N个锂离子电池单体的容量和中值电压,其中,待测锂离子电池组由N个锂离子电池单体串联成组,且平均分n个电池组;
S2、求出任意一个电池组的容量Yj(j=1,2,……n)的分布函数FY(y),y代表锂离子电池单体的容量;
S3、根据Yj的分布函数FY(y)求出Yj的期望值EYj和方差Var(Yj);
S4、根据式(17)计算n个电池组的总容量Qsum的临界值Tn
式(17)中,φ为标准正态分布的分布函数;α为置信区间;
S5、将Tn代入式(10)中的Qsum,求得En即为待测锂离子电池组可用能量:
En=Uav·Qsum (10);
式(10)中,Uav为待测锂离子电池组中所有锂离子电池单体的中值电压均值。
2.如权利要求1所述的锂离子电池组可用能量计算方法,其特征在于:步骤S2中,任意一个电池组的容量Yj(j=1,2,……n)的分布函数FY(y)表示为:
FY(y)=P(Y1≤y)=1-(1-F(y))N/n (11);
Yj的密度函数fY(y)为:
fY(y)=F'Y(y)=N(1-F(y))N/n-1f(y)/n (12);
式(11)和(12)中,
F(y)为锂离子电池单体容量的分布函数;
f(y)为锂离子电池单体容量的概率密度函数:
式(1)中,μ和σ2分别为锂离子电池单体的容量的均值和方差。
CN201611076522.5A 2016-11-29 2016-11-29 一种锂离子电池组可用能量计算方法 Active CN107045105B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611076522.5A CN107045105B (zh) 2016-11-29 2016-11-29 一种锂离子电池组可用能量计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611076522.5A CN107045105B (zh) 2016-11-29 2016-11-29 一种锂离子电池组可用能量计算方法

Publications (2)

Publication Number Publication Date
CN107045105A true CN107045105A (zh) 2017-08-15
CN107045105B CN107045105B (zh) 2019-10-18

Family

ID=59543507

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611076522.5A Active CN107045105B (zh) 2016-11-29 2016-11-29 一种锂离子电池组可用能量计算方法

Country Status (1)

Country Link
CN (1) CN107045105B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111142028A (zh) * 2020-02-20 2020-05-12 广东天波信息技术股份有限公司 一种锂电池电量的测量方法及装置
CN112462286A (zh) * 2020-10-13 2021-03-09 中车长春轨道客车股份有限公司 一种基于能量对锂电池健康状态的估计方法
CN116774044A (zh) * 2023-05-23 2023-09-19 山东科技大学 一种退役锂电池组的健康状态评价方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001045611A (ja) * 1999-08-02 2001-02-16 Mazda Motor Corp ハイブリッド自動車のパワートレンの故障判断装置
CN101131414A (zh) * 2006-08-25 2008-02-27 比亚迪股份有限公司 一种评价球镍中值电压性能的方法
CN101930056A (zh) * 2009-06-24 2010-12-29 力博特公司 一种电池备电时间预测方法
CN103499792A (zh) * 2013-07-18 2014-01-08 浙江工业大学 电动汽车动力电池集群可用容量的预测方法
CN105044605A (zh) * 2015-06-26 2015-11-11 江苏北新能源科技有限公司 一种获取电池中值电压的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001045611A (ja) * 1999-08-02 2001-02-16 Mazda Motor Corp ハイブリッド自動車のパワートレンの故障判断装置
CN101131414A (zh) * 2006-08-25 2008-02-27 比亚迪股份有限公司 一种评价球镍中值电压性能的方法
CN101930056A (zh) * 2009-06-24 2010-12-29 力博特公司 一种电池备电时间预测方法
CN103499792A (zh) * 2013-07-18 2014-01-08 浙江工业大学 电动汽车动力电池集群可用容量的预测方法
CN105044605A (zh) * 2015-06-26 2015-11-11 江苏北新能源科技有限公司 一种获取电池中值电压的方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111142028A (zh) * 2020-02-20 2020-05-12 广东天波信息技术股份有限公司 一种锂电池电量的测量方法及装置
CN111142028B (zh) * 2020-02-20 2022-02-08 广东天波信息技术股份有限公司 一种锂电池电量的测量方法及装置
CN112462286A (zh) * 2020-10-13 2021-03-09 中车长春轨道客车股份有限公司 一种基于能量对锂电池健康状态的估计方法
CN116774044A (zh) * 2023-05-23 2023-09-19 山东科技大学 一种退役锂电池组的健康状态评价方法
CN116774044B (zh) * 2023-05-23 2024-02-20 山东科技大学 一种退役锂电池组的健康状态评价方法

Also Published As

Publication number Publication date
CN107045105B (zh) 2019-10-18

Similar Documents

Publication Publication Date Title
CN106291372B (zh) 一种新的锂离子动力电池剩余寿命预测方法
CN103545567B (zh) 一种快速分选锂离子电池的方法
CN103792495B (zh) 基于德尔菲法和灰色关联理论的电池性能评价方法
CN107607881A (zh) 一种锂离子动力电池自放电一致性的评价方法
CN108089133A (zh) 储能系统电池组一致性检测方法及检测装置
CN106824831A (zh) 一种提高锂离子电池一致性的动力型电池的制造方法
CN101504443B (zh) 一种锂离子电池放电容量的预测方法
CN106154165A (zh) 一种大容量电池储能系统性能的评估方法和评估系统
CN103412264B (zh) 蓄电池组内单体电池一致性的评价方法
CN108037462A (zh) 蓄电池健康状况量化方法及系统
CN102879743B (zh) 锂离子电池剩余电量值实时分析计算法
CN107607874A (zh) 快速充/放电锂离子电池的三点式筛选方法
CN107843846A (zh) 一种锂离子电池健康状态估计方法
CN110031771A (zh) 一种描述电池一致性的方法
CN103956513B (zh) 一种大容量锂动力电池的配组方法
CN105652106A (zh) 一种储能系统综合性能测试方法及装置
CN103579700A (zh) 一种锂离子电池分选配组方法
CN107064806A (zh) 一种锂电池容量衰减率计算方法
CN107377422A (zh) 一种单体电池的分选方法
CN111458649A (zh) 一种电池模组健康度快速检测方法
CN104502844A (zh) 一种基于交流阻抗的动力锂电池劣化程度诊断方法
CN111036575A (zh) 一种基于温度变化分析的锂离子电池分选方法
CN103487758B (zh) 一种锂离子电池配组方法
CN107045105B (zh) 一种锂离子电池组可用能量计算方法
CN107132481A (zh) 一种识别电池组中单体一致的方法及系统

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