CN106096281B - 一种考虑单体不一致的锂离子动力电池组在线参数估计方法 - Google Patents
一种考虑单体不一致的锂离子动力电池组在线参数估计方法 Download PDFInfo
- Publication number
- CN106096281B CN106096281B CN201610417242.XA CN201610417242A CN106096281B CN 106096281 B CN106096281 B CN 106096281B CN 201610417242 A CN201610417242 A CN 201610417242A CN 106096281 B CN106096281 B CN 106096281B
- Authority
- CN
- China
- Prior art keywords
- monomer
- moment
- weight vector
- battery pack
- parameter
- 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
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Secondary Cells (AREA)
Abstract
本发明涉及一种考虑单体不一致的锂离子动力电池组在线参数估计方法,用以获取锂离子动力电池组中各个单体的参数,包括欧姆内阻、极化电阻和时间常数,包括以下步骤:1)获取锂离子动力电池组的总体平均参数;2)分别构造欧姆内阻权向量A、极化电阻权向量B和时间常数权向量C并根据总体平均参数分别获取锂离子动力电池组中各个单体的参数。与现有技术相比,本发明具有考虑单体不一致性、计算简单等优点。
Description
技术领域
本发明涉及动力电池组领域,尤其是涉及一种考虑单体不一致的锂离子动力电池组在线参数估计方法。
背景技术
目前的电池管理***(Battery Management System,BMS)对电池内部参数进行估计时,往往将整个电池组看成均匀的。实际上,电池在制作和装配过程中,由于工艺和材质的不均匀等,使得同一类型、规格、型号的电池在出厂时的参数不可能完全一致。而在电池的使用过程中,每个电池单体的充放电过程、自放电程度、通风条件和温度环境等的差异在一定程度上也会增加电池参数的不一致性电池组内部各单体之间存在不一致。对于成组的电池来说,由于其内阻不一致性,在充电过程中,相同的充电电流下,内阻大的电池升压高,容易导致过充;而在放电过程中,相同的放电电流下,内阻大的电池压降高,容易导致过放。因此,以平均参数代替各个单体参数进行管理是不合理的。在电池的使用过程中,随着循环次数的增加,电池的不一致性会加剧,电池组很有可能会因为某一节单体的损坏而提前报废,甚至引起安全隐患,导致***使用和维护成本增加。这一问题已经成为制约电动汽车产业化发展和应用的瓶颈问题之一。
由于单体不一致的问题,我们必须要对每一节单体电池进行参数估计,对于BMS来说,电池信息显示必须实时。因此,如果我们直接将传统参数辨识方法在每节单体电池上实施,将导致计算量过大,从而就对处理器提出较高的要求,势必会增加BMS的成本。因此,本发明所提的方法是对现有估计算法进行改进和优化,使之能在不增加计算量或者只增加少量计算量的前提下估计出考虑一致性的参数值。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种考虑单体不一致性、计算简单的考虑单体不一致的锂离子动力电池组在线参数估计方法。
本发明的目的可以通过以下技术方案来实现:
一种考虑单体不一致的锂离子动力电池组在线参数估计方法,用以获取锂离子动力电池组中各个单体的参数,包括欧姆内阻、极化电阻和时间常数,包括以下步骤:
1)获取锂离子动力电池组的总体平均参数;
2)分别构造欧姆内阻权向量A、极化电阻权向量B和时间常数权向量C并根据总体平均参数分别获取锂离子动力电池组中各个单体的参数。
所述的步骤1)中通过最小二乘法或自适应滤波法获取锂离子动力电池组的总体平均参数。
当通过最小二乘法获取锂离子动力电池组的总体平均参数时包括以下步骤:
11)根据多个电池单体串联的电池组结构建立电池平均模型,该模型为:
其中,UL为端电压,Uoc为开路电压,UThm为平均极化电压,IL为电流,R0m为平均欧姆内阻,RThm为平均极化电阻,CThm为平均极化电容;
12)构建电池平均模型的差分方程,并采用递推公式,获取动力电池组的总体平均参数,包括平均欧姆内阻、平均极化电阻和平均时间常数,
所述的差分方程为:
Ud,k+1=aUd,k+bIL,k+1+cIL,k
递推公式为:
zk+1=Ud,k+1
hk+1=[Ud,k+1,IL,k+1,IL,k]
动力电池组的总体平均参数为:
其中,Ud,k、Ud,k+1分别为电池组k时刻和k+1时刻的平均动态电压向量,IL,k、IL,k+1分别为电池组k时刻和k+1时刻的电流,τm为平均时间常数且τm=RThmCThm,Δt为采样周期,Pk、Pk+1分别为k和k+1时刻的协方差矩阵,Lk+1为k+1时刻的增益矩阵。
所述的步骤2)中,欧姆内阻权向量A、极化电阻权向量B和时间常数权向量C的表达式为:
其中,Udn为电池组中第n个单体的平均动态电压,UBn为电池组中第n个单体基于权向量B的动态电压,UCn为电池组中第n个单体基于权向量C的动态电压。
所述的步骤2)中,
欧姆内阻权向量A的递推公式为:
Up,k=Udmean,kAk-1
Ak=Ak-1+gAUdmean,k(Ud,k-Up,k)
极化电阻权向量B的递推公式为:
UBp,k=UBmean,kBk-1
Bk=Bk-1+gBUBmean,k(UB,k-UBp,k)
时间常数权向量C的递推公式为:
UCp,k=UCmean,kCk-1
Ck=Ck-1+gCUCmean,k(UC,k-UCp,k)
其中,Up,k为电池组k时刻的预测动态电压,Udmean,k为k时刻所有电池组单体动态电压的算数平均值,Ud,k为电池组k时刻的平均动态电压向量,gA为欧姆增益,UBp,k为电池组k时刻的UB的预测值,UBmean,k为k时刻电池组各单体对应UB的算数平均值,gB为极化增益,UCp,k为电池组k时刻的UC的预测值,UCmean,k为k时刻电池组各单体对应UB的算数平均值,gC为时间增益。
与现有技术相比,本发明具有以下优点:
一、考虑单体不一致性:本发明分别根据每节单体的端电压,荷电状态等参数估计出每节电池的参数,而不是把整包电池看成均匀状态来求出平均参数。
二、计算简单:本发明从原理出发,推导出考虑单体不一致性的权向量,计算时只需将平均参数与权向量相乘即可得到每节单体的参数,与目前的方法相比,大大减小了计算量。
附图说明
图1为本发明算法流程图。
图2为本发明实施例涉及的电池模型图。
图3为本发明实施例的两种电池模型示意图,其中,图(3a)为本发明实施例的电池串联模型示意图,图(3b)为本发明实施例的电池平均模型示意图。
图4为本发明实施例的电池参数曲线,其中,图(4a)为本发明实施例中两节电池动态电压曲线,图(4b)为本发明实施例中两节电池欧姆内阻曲线。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例:
为了使本发明的目的、技术方案及创新点更加清晰,以下结合附图及实施例对本发明作进一步的阐释。应当理解,此处所描述的具体实施例仅用于解释本发明,并不用于限定本发明。
如图1所示,图为本发明的算法流程图。
图2所示为本发明实施例涉及的电池模型图,其中Uoc为开路电压,UL为端电压,IL为电流R0为欧姆内阻,RTh和CTh分别为极化电阻和极化电容,ITh为极化电容的电流,UTh为极化电压,模型的数学关系为:
为简化起见,假设一个简单的电池模组,其中有两节电池串联,其阻抗模型如图3a所示。首先,可以构建一个上述两节电池的平均阻抗模型,如图3b所示,R0m为平均欧姆内阻,RThm为平均极化电阻,CThm平均极化电容。通过一般的参数估计方法,我们可以求出电池组的平均参数,此处以最小二乘法为例。电池组的差分方程如下式:
Ud,k+1=aUd,k+bIL,k+1+cIL,k
上式中,Ud,k为k时刻平均端电压和平均开路电压之差,即电池组平均动态电压,IL,k为k时刻的电流,a、b、c为:
上式中,R0m为平均欧姆内阻,RThm为平均极化电阻,τm为平均时间常数,e为自然对数,△t为采样周期。最小二乘法的递推公式如下:
上式中,为k时刻的参数估计值,zk+1=Ud,k+1,hk+1=[Ud,k+1,IL,k+1,IL,k]。设定好初值之后可以递推出a、b、c的值,然后根据下式可以求出平均参数:
求出平均参数之后,下面求单体参数,如图3a所示,有:
R01:R02=U01:U02
上式中U01和U02分别为两节电池欧姆内阻R01和R02上的电压。但是在实际计算中U01和U02的求解比较困难,而实际情况中发现欧姆内阻之比与动态电压之比比较接近。本发明的实施例,如图4a和4b所示,在动态电压曲线上随机选取了两个时刻,两节电池动态电压的比值为1.06和1.07,两节电池欧姆内阻的比值为1.05。
因此,我们可以推断:
R01:R02≈Ud1:Ud2
根据上述比例关系,我们可以构造一个权向量A:
式中,
根据平均欧姆内阻和向量A,我们可以求出两节电池的欧姆内阻:
[R01,R02]=R0mA
这样就可以求出两节电池的欧姆内阻,然后根据动态电压Ud和电流IL可以求出极化电压UTh:
根据模型的数学关系可以推出:
由于采样周期△t一般远远小于时间常数τ,因此有:
从而得到:
上式整理得:
令:
UB=UTh,k+1/UTh,k
这样我们就可以构造一个权向量B:
式中:
根据平均时间常数τm和权向量B,我们可以求出:
这样就求出了两节电池的时间常数τ1和τ2,由电池模型的数学关系得:
因此有两节电池的极化电阻之比:
令:
这样我们就可以构造一个权向量C:
根据RThm和权向量C,我们可以求出:
[RTh1,RTh2]=RThmC
然后,根据两节电池的时间常数和极化电阻的值可以求出两节电池的极化电容:
以上是两节电池的的推导,同理可以推出多节电池的权向量:
如果直接利用上述权向量进行求解,由于电池电压的波动比较大,在参数估计过程中会造成参数比较大的波动。首先,我们必须对权向量进行归一化处理,权向量实质上反应了所有电池参数的比率,当参数全部相同时,A=B=C=[1,1,1,……,1],权向量中的元素大于1说明对应的参数大于平均参数,小于1说明对应的参数小于平均参数。
其次,为了减缓波动,我们可以采用递推的方式得到权向量。下面以基于动态电压Ud求权向量Ak为例,说明权向量的求解方法。k为当前时刻,首先利用上一时刻的权向量预测当前的时刻的动态电压值:
Up,k=Udmean,kAk-1
上式中Udmean,k为当前时刻的平均动态电压值,然后利用预测的动态电压值和真实动态电压值之差来修正当前时刻的权向量:
Ak=Ak-1+gAUdmean,k(Ud,k-Up,k)
上式中,gA为增益,增益越大,修正量越大,算法越快,也可能出现过修调的情况,因此精度越低。
同理可以求出权向量Bk和Ck的递推公式:
上两式中,UBmean,k和UCmean,k分别为k时刻UB和UC的平均值,gB和gC分别为权向量B和C的增益。
Claims (2)
1.一种考虑单体不一致的锂离子动力电池组在线参数估计方法,用以获取锂离子动力电池组中各个单体的参数,包括欧姆内阻、极化电阻和时间常数,其特征在于,包括以下步骤:
1)获取锂离子动力电池组的总体平均参数,通过最小二乘法或自适应滤波法获取锂离子动力电池组的总体平均参数,当通过最小二乘法获取锂离子动力电池组的总体平均参数时包括以下步骤:
11)根据多个电池单体串联的电池组结构建立电池平均模型,该模型为:
其中,UL为端电压,Uoc为开路电压,UThm为平均极化电压,IL为电流,R0m为平均欧姆内阻,RThm为平均极化电阻,CThm为平均极化电容;
12)构建电池平均模型的差分方程,并采用递推公式,获取动力电池组的总体平均参数,包括平均欧姆内阻、平均极化电阻和平均时间常数,
所述的差分方程为:
Ud,k+1=aUd,k+bIL,k+1+cIL,k
递推公式为:
zk+1=Ud,k+1
hk+1=[Ud,k+1,IL,k+1,IL,k]
动力电池组的总体平均参数为:
其中,Ud,k、Ud,k+1分别为电池组k时刻和k+1时刻的平均动态电压向量,IL,k、IL,k+1分别为电池组k时刻和k+1时刻的电流,τm为平均时间常数且τm=RThmCThm,Δt为采样周期,Pk、Pk+1分别为k和k+1时刻的协方差矩阵,Lk+1为k+1时刻的增益矩阵;
2)分别构造欧姆内阻权向量A、极化电阻权向量B和时间常数权向量C并根据总体平均参数分别获取锂离子动力电池组中各个单体的参数,所述的步骤2)中,欧姆内阻权向量A、极化电阻权向量B和时间常数权向量C的表达式为:
其中,Udn为电池组中第n个单体的平均动态电压,UBn为电池组中第n个单体基于权向量B的动态电压,UCn为电池组中第n个单体基于权向量C的动态电压。
2.根据权利要求1所述的一种考虑单体不一致的锂离子动力电池组在线参数估计方法,其特征在于,所述的步骤2)中,
欧姆内阻权向量A的递推公式为:
Up,k=Udmean,kAk-1
Ak=Ak-1+gAUdmean,k(Ud,k-Up,k)
极化电阻权向量B的递推公式为:
UBp,k=UBmean,kBk-1
Bk=Bk-1+gBUBmean,k(UB,k-UBp,k)
时间常数权向量C的递推公式为:
UCp,k=UCmean,kCk-1
Ck=Ck-1+gCUCmean,k(UC,k-UCp,k)
其中,Up,k为电池组k时刻的预测动态电压,Udmean,k为k时刻所有电池组单体动态电压的算数平均值,Ud,k为电池组k时刻的平均动态电压向量,gA为欧姆增益,UBp,k为电池组k时刻的UB的预测值,UBmean,k为k时刻电池组各单体对应UB的算数平均值,gB为极化增益,UCp,k为电池组k时刻的UC的预测值,UCmean,k为k时刻电池组各单体对应UB的算数平均值,gC为时间增益。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610417242.XA CN106096281B (zh) | 2016-06-15 | 2016-06-15 | 一种考虑单体不一致的锂离子动力电池组在线参数估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610417242.XA CN106096281B (zh) | 2016-06-15 | 2016-06-15 | 一种考虑单体不一致的锂离子动力电池组在线参数估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106096281A CN106096281A (zh) | 2016-11-09 |
CN106096281B true CN106096281B (zh) | 2018-10-26 |
Family
ID=57846877
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610417242.XA Active CN106096281B (zh) | 2016-06-15 | 2016-06-15 | 一种考虑单体不一致的锂离子动力电池组在线参数估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106096281B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109031145B (zh) * | 2018-08-10 | 2020-12-11 | 山东大学 | 一种考虑不一致性的串并联电池组模型及实现方法 |
CN109212429B (zh) * | 2018-08-14 | 2021-02-09 | 国网冀北电力有限公司秦皇岛供电公司 | 多参数加权判定蓄电池性能的方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103793605A (zh) * | 2014-01-26 | 2014-05-14 | 浙江省计量科学研究院 | 基于粒子群算法的磷酸铁锂动力电池等效电路模型参数估计的方法 |
-
2016
- 2016-06-15 CN CN201610417242.XA patent/CN106096281B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103793605A (zh) * | 2014-01-26 | 2014-05-14 | 浙江省计量科学研究院 | 基于粒子群算法的磷酸铁锂动力电池等效电路模型参数估计的方法 |
Non-Patent Citations (4)
Title |
---|
Model-based dynamic multi-parameter method for peak power estimation of lithium–ion batteries;FengChun Sun et al;《Applied Energy》;20120328;378-386 * |
电动汽车锂离子电池等效电路模型的参数辨识研究;李思;《中国优秀硕士学位论文全文数据库(电子期刊)工程科技Ⅱ辑》;20150715;第3.1-3.2节 * |
考虑电池组不一致性的储能***建模及仿真;范刘洋等;《电力***自动化》;20160210;第40卷(第3期);110-115 * |
锂离子动力电池参数辨识与状态估计;魏学哲等;《同济大学学报(自然科学版)》;20080229;第36卷(第2期);231-235 * |
Also Published As
Publication number | Publication date |
---|---|
CN106096281A (zh) | 2016-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
von Lüders et al. | Modeling of lithium plating and lithium stripping in lithium-ion batteries | |
US10027134B2 (en) | Active equalization method and system of lithium iron phosphate battery pack | |
CN106291378B (zh) | 一种电动汽车动力电池soh的测算方法 | |
Kang et al. | Comparison of comprehensive properties of Ni-MH (nickel-metal hydride) and Li-ion (lithium-ion) batteries in terms of energy efficiency | |
WO2022105104A1 (zh) | 一种基于多新息递推贝叶斯算法的电池模型参数辨识方法 | |
CN110940920B (zh) | 锂电池在预定soc下的不析锂的最大充电电流的获取方法 | |
EP3786653B1 (en) | Battery state estimation method | |
KR102169774B1 (ko) | 리튬-이온 배터리 시스템의 개개의 전극들의 용량 및 총 용량을 추정하기 위한 방법 및 시스템 | |
EP3475712B1 (en) | Method and device for estimating a voltage of a battery | |
EP3141919B1 (en) | Apparatus and method for estimating open circuit voltage | |
CN204269787U (zh) | 一种锂离子电池低温性能一致性的检测*** | |
KR20200049596A (ko) | 리튬 이온 배터리를 급속 충전하기 위한 시스템 및 방법 | |
US10935607B2 (en) | Control device for secondary battery, control method of secondary battery, battery system, and motor-driven vehicle | |
US20130314042A1 (en) | Method for Ascertaining the Open Circuit Voltage of a Battery, Battery with a Module for Ascertaining the Open Circuit Voltage and a Motor Vehicle Having a Corresponding Battery | |
CN113779794B (zh) | 计及微观约束的锂离子电池sop估计方法及*** | |
CN105518927B (zh) | 用于确定电池组的温度的方法 | |
CN110854959B (zh) | 电池***被动均衡方法 | |
EP3988954B1 (en) | Method for detecting internal short-circuited cell | |
EP4152022B1 (en) | Method for determining full-charge capacity of battery pack, method for determining state of health of battery pack, system, and apparatus | |
CN111987377B (zh) | 电池的快充方法、电池管理***以及电池的快充装置 | |
CN106096281B (zh) | 一种考虑单体不一致的锂离子动力电池组在线参数估计方法 | |
CN105811028A (zh) | 一种锂离子电池***的soc状态估计方法 | |
Tran et al. | State of charge estimation of lithium ion batteries using an extended single particle model and sigma-point Kalman filter | |
Gao et al. | An optimal charging strategy of lithium-ion batteries based on polarization and temperature rise | |
Liu et al. | Improving battery safety for electric vehicles through the optimization of battery design parameters |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |