CN106872904B - 基于离散滑模观测器的锂电池内部温度估计方法 - Google Patents
基于离散滑模观测器的锂电池内部温度估计方法 Download PDFInfo
- Publication number
- CN106872904B CN106872904B CN201710098924.3A CN201710098924A CN106872904B CN 106872904 B CN106872904 B CN 106872904B CN 201710098924 A CN201710098924 A CN 201710098924A CN 106872904 B CN106872904 B CN 106872904B
- Authority
- CN
- China
- Prior art keywords
- lithium battery
- temperature
- moment
- discrete
- mode observer
- 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
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 title claims abstract description 196
- 229910052744 lithium Inorganic materials 0.000 title claims abstract description 196
- 238000000034 method Methods 0.000 title claims abstract description 55
- 239000011159 matrix material Substances 0.000 claims description 34
- 230000014509 gene expression Effects 0.000 claims description 10
- 230000020169 heat generation Effects 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 7
- 230000005611 electricity Effects 0.000 claims 1
- 238000004519 manufacturing process Methods 0.000 abstract description 3
- 238000001914 filtration Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 230000032683 aging Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 239000011258 core-shell material Substances 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001566 impedance spectroscopy Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000003446 memory effect Effects 0.000 description 1
- 239000000178 monomer Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 229910001220 stainless steel Inorganic materials 0.000 description 1
- 239000010935 stainless steel Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/367—Software therefor, e.g. for battery testing using modelling or look-up tables
Landscapes
- Secondary Cells (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Battery Electrode And Active Subsutance (AREA)
Abstract
本发明公开了一种基于离散滑模观测器的锂电池内部温度估计方法。它包括以下步骤:首先,建立用于内部温度估计的锂电池离散状态空间模型;其次,实时采集工况下锂电池的工作电压、工作电流、表面温度和外界环境温度;然后,在线辨识锂电池离散状态空间模型参数;最后,构建离散滑模观测器,实现对锂电池内部温度的实时估计。本发明方法不仅具有良好的温度估计效果,同时能够严格保证观测器的收敛性,且对锂电池热模型建模误差、模型参数的摄动和外部扰动表现出较强的鲁棒性。
Description
技术领域
本发明涉及车用锂电池技术领域,具体涉及一种基于离散滑模观测器的锂电池内部温度估计方法。
背景技术
锂电池具有能量密度高、无记忆效应、环境污染低、循环寿命长等优点,在电动汽车上得到越来越多的应用。锂电池作为电动汽车的关键部件,其温度过高或过低对锂电池的性能、寿命和安全性有重要的影响,进而影响整车的动力性、经济性和安全性。
电池温度是电池管理系统进行热管理的一个重要参数,其数值对于估计电池剩余电量、制定热管理策略均具有重要的意义。
目前,传统的电池温度直接测量的方法主要有以下两种:一种方法是通过在电池表面安装温度传感器进行测量,得到电池表面温度,但是这种方法存在滞后性,测量得到的表面温度无法真实反映电池的内部状态,这是因为电池的热量是自内而外散发的,这个导热过程需要一段时间,因此电池内外存在温差;另一种是通过在电池内部内置温度传感器来测量电池内部温度,但这一方法目前只能在电池生产过程中进行,或在实验室环境下使用,对实际工况下的电池并不适用,这是因为在每个电池内部内置传感器不仅会增加成本,而且会降低电池的性能和安全性。
除了上述采用温度传感器的直接测量方法之外,还可以通过软测量的方式来估计出锂电池的内部温度。目前,用于估计锂电池内部温度方法有卡尔曼滤波算法、电池内阻-温度法、温度-阻抗谱法等。
中国发明专利(CN104462847A)于2014年12月23日公开的《一种电池内部温度实时预测方法》,它将电池分为内核和外壳,分别建立温度预测模型;通过实验的方法获取电池内部和表面材质比热容、电池内核-外壳以及外壳-外界的热阻参数、电池开路电压曲线、开路电压随温度变化曲线等信息;利用卡尔曼滤波方法对电池内部温度进行实时跟踪和修正,将实时电池表面和环境温度输入预测模型,实时预测电池内部温度。但是该方法存在不足:
1)卡尔曼滤波算法要求过程噪声和观测噪声是白噪声且噪声的均值、方差等统计特性已知,这在实际应用中是很难满足的,不仅因为噪声统计特性很难获得,而且白噪声也仅在理想条件下存在;
2)卡尔曼滤波算法对锂电池温度预测模型精度要求较高,当模型中的参数因温度、老化、寿命等因素的影响而发生变化时,会造成较大的温度预测误差,甚至发散。
中国发明专利(CN104865534A)于2015年08月26日公开了《一种单体电池内部温度估计方法》,该方法包括离线部分和在线部分,其中,离线部分包括:a1)获得电池单体在不同温度情况下的内阻特性;a2)根据所述内阻特性获取电池内部温度与电池内阻的标准关系;在线部分包括:b1)在线检测当前工作电池的电流、端电压及表面温度;b2)根据所述电流和端电压在线估计当前工作电池的电池内阻;b3)根据步骤b2)获得的电池内阻、步骤b1)获得的表面温度及步骤a2)获得的标准关系估计当前工作电池的内部温度。其存在的不足为:
仅考虑了电池内阻与电池内部温度之间的关系,采用离线方式获取电池内部温度与电池内阻的标准关系,而锂电池工作过程中,除了温度,电池内阻这一参数还受其他因素的影响而发生变化,比如循环次数、充放电倍率、荷电状态、老化程度等,因此会存在较大的估计误差。
发明内容
为了克服上述现有技术的不足,本发明提供了一种基于离散滑模观测器的锂电池内部温度估计方法。该方法不仅能够仅根据噪声的上界值和采集的锂电池电流、电压、表面温度,采用在线的方式对锂电池热模型的模型参数进行辨识,实现对锂电池内部温度的估计,而且该方法对锂电池热模型的建模误差和对因温度、老化、循环次数、充放电倍率等因素引起的模型参数变化等表现出较强的鲁棒性,同时能严格保证温度估计算法的收敛性,避免了估计发散现象的发生。
本发明的目的是这样实现的。本发明提供了一种基于离散滑模观测器的锂电池内部温度估计方法,主要步骤如下:
步骤1,根据锂电池内部与外壳之间的热平衡以及锂电池外壳与外界环境之间的热平衡,建立用于锂电池内部温度估计的离散状态空间模型:
x(k+1)=Adx(k)+Bdu(k)+w(k) (1)
y(k)=Cdx(k)+v(k) (2)
其中,
x(k)为锂电池在k时刻的状态变量,x(k)=[Tia(k) Tsa(k)]T,Tia(k)为k时刻锂电池内部温度Tin(k)与外界环境温度Tamb(k)之间的相对温度值,Tia(k)=Tin(k)-Tamb(k),Tsa(k)为k时刻锂电池表面温度Ts(k)与外界环境温度Tamb(k)之间的相对温度值,Tsa(k)=Ts(k)-Tamb(k);
x(k+1)为锂电池在(k+1)时刻的状态变量;
u(k)为锂电池在k时刻的输入,u(k)=Q(k),Q(k)为k时刻锂电池的生热速率;
y(k)为锂电池在k时刻的输出;
v(k)为锂电池在k时刻的观测噪声;
Ad为离散状态空间模型的系数矩阵,Ri为锂电池内部与锂电池外壳之间的等效热阻,Ro为锂电池外壳与外界环境之间的等效热阻,Cc为锂电池内部的等效比热容,Cs为锂电池外壳的比热容,Cs为常数,Δt为采样时间;
Bd为离散状态空间模型的输入矩阵,Cc为锂电池内部的等效比热容,Δt为采样时间;
Cd为离散状态空间模型的观测矩阵,Cd=[0 1],w(k)为锂电池在k时刻的过程噪声;
步骤2,分别通过电压传感器和电流传感器实时采集工况下k时刻锂电池的工作电压U(k)和k时刻锂电池的工作电流I(k),通过温度传感器实时采集工况下k时刻锂电池表面温度Ts(k)和外界环境温度Tamb(k);
步骤3,根据步骤2中采集到的k时刻锂电池工作电压U(k)、工作电流I(k)、表面温度Ts(k)和外界环境温度Tamb(k),在线辨识步骤1中锂电池离散状态空间模型中的参数,包括锂电池内部与锂电池外壳之间的等效热阻Ri、锂电池外壳与外界环境之间的等效热阻Ro和锂电池内部的等效比热容Cc;
步骤4,基于步骤1中建立的锂电池离散状态空间模型构建离散滑模观测器,并结合步骤2中采集到的k时刻锂电池工作电压U(k)、工作电流I(k)、表面温度Ts(k)和外界环境温度Tamb(k)以及步骤3中所辨识的锂电池离散状态空间模型参数Ri、Ro和Cc,实时获取Tia(k)的估计值从而实时估计锂电池内部温度,锂电池内部温度Tin(k)的实时估计值为
步骤4.1,构建离散滑模观测器的方程式如下:
其中,为x(k)的估计值,为x(k+1)的估计值,为y(k)的估计值,L为离散滑模观测器的增益矩阵,M为饱和函数增益矩阵,为饱和函数, 为边界层;
步骤4.2,将步骤2中采集到的k时刻锂电池工作电压U(k)、工作电流I(k)作为信号输入,利用步骤4.1中所构建的离散滑模观测器得到状态变量Tia(k)的实时估计值
步骤4.3,根据步骤4.2得到的状态变量Tia(k)的实时估计值和步骤3中采集到的k时刻外界环境温度Tamb(k),同时结合表达式实时估计锂电池内部温度,得到k时刻锂电池内部温度Tin(k)的实时估计值
优选地,步骤1中的k时刻锂电池的生热速率Q(k)为:
其中,I(k)为k时刻锂电池的工作电流,Vs为锂电池的体积,Eoc(k)为k时刻锂电池的电动势,U(k)为k时刻锂电池的工作电压,Ts(k)为k时刻锂电池的表面温度,α为锂电池电动势Eoc(k)的温度影响系数。
优选地,步骤1中的锂电池在k时刻的观测噪声v(k)有界,即满足以下表达式:
|v(k)|≤v0 (6)
其中,v0为v(k)上界参数,且v0为正实数。
优选地,步骤2中所述的电压传感器、电流传感器和温度传感器分别为霍尔电压传感器、霍尔电流传感器和热电偶。
优选地,步骤3中的在线辨识锂电池离散状态空间模型参数的方法为最小二乘法。
优选地,步骤1中的过程噪声w(k)和观测噪声v(k),以及步骤4中的离散滑模观测器的增益矩阵L满足以下表达式:
|w(k)-Lv(k)|≤Δ0 (7)
其中,Δ0为(w(k)-Lv(k))的上界参数,且Δ0为2×1的矩阵,且Δ0的每个元素都是正实数。
优选地,步骤4中的离散滑模观测器的增益矩阵L、饱和函数增益矩阵M、和边界层在进行选取时,满足以下表达式:
M≥Δ0 (9)
其中,Δ0为(w(k)-Lv(k))的上界参数,且Δ0为2×1的矩阵,且Δ0的每个元素都是正实数;v0是v(k)的上界参数,且v0为正实数。
本发明的有益的效果为:
1、对用于锂电池内部温度估计的模型精度要求不高,采用在线辨识方法对锂电池离散状态空间模型参数进行辨识,当温度、老化、循环次数、充放电倍率等因素引起模型参数发生变化时,仍能准确的估计锂电池内部温度,表现出较强的鲁棒性。
2、避免了卡尔曼滤波算法中因锂电池离散状态空间模型参数发生变化而引起的离散滑模观测器发散问题,能严格保证算法的收敛性。
3、避免了卡尔曼滤波算法中噪声统计特性的获取,仅需要噪声的上界值。
附图说明
图1是本发明锂电池内部温度估计方法流程示意图。
图2是锂电池内部-外壳-外界环境的热平衡示意图。
图3是锂电池的等效热模型示意图。
具体实施方式
下面将结合附图及实施例对本发明作进一步详细说明,但本发明的实施例不限于此。
图1是本发明锂电池内部温度估计方法流程示意图,由该图可见,本发明提供的一种基于离散滑模观测器的锂电池内部温度估计方法,包括以下步骤:
步骤1,根据锂电池内部-外壳的热平衡以及锂电池外壳-外界环境的热平衡,建立用于锂电池内部温度估计的锂电池离散状态空间模型。
锂电池内部温度估计的离散状态空间模型的表达式为:
x(k+1)=Adx(k)+Bdu(k)+w(k)
y(k)=Cdx(k)+v(k)
其中,
x(k)为锂电池在k时刻的状态变量,x(k)=[Tia(k) Tsa(k)]T,Tia(k)为k时刻锂电池内部温度Tin(k)与外界环境温度Tamb(k)之间的相对温度值,Tia(k)=Tin(k)-Tamb(k),Tsa(k)为k时刻锂电池表面温度Ts(k)与外界环境温度Tamb(k)之间的相对温度值,Tsa(k)=Ts(k)-Tamb(k);
x(k+1)为锂电池在(k+1)时刻的状态变量;
u(k)为锂电池在k时刻的输入,u(k)=Q(k),Q(k)为k时刻锂电池的生热速率;
y(k)为锂电池在k时刻的输出;
v(k)为锂电池在k时刻的观测噪声;
Ad为离散状态空间模型的系数矩阵,Ri为锂电池内部与锂电池外壳之间的等效热阻,Ro为锂电池外壳与外界环境之间的等效热阻,Cc为锂电池内部的等效比热容,Cs为锂电池外壳的比热容,Cs为常数,Δt为采样时间;
Bd为离散状态空间模型的输入矩阵,Cc为锂电池内部的等效比热容,Δt为采样时间;
Cd为离散状态空间模型的观测矩阵,Cd=[0 1],w(k)为锂电池在k时刻的过程噪声。
具体的,离散状态空间模型的建立分为以下步骤。
步骤1.1,图2为锂电池内部与外壳之间的热平衡以及锂电池外壳与外界环境之间的热平衡示意图,图3是锂电池的等效热模型,首先根据图2和图3分别建立锂电池热平衡微分方程:
其中,Tin为锂电池内部温度,Ts为锂电池表面温度,Tamb为外界环境温度,Ri为锂电池内部与锂电池外壳之间的等效热阻,Ro为锂电池外壳与外界环境之间的等效热阻,Cc为锂电池内部的等效比热容,Cs为锂电池外壳的比热容,在本实施例中,电池外壳材料为不锈钢,取值Cs=502.48J/(kg·K),Q为锂电池的生热速率。
在本实施例中,锂电池的生热速率Q的表达式为:
其中,I为锂电池的工作电流,放电时I取正值,充电时I取负值,Eoc为锂电池的电动势,U为锂电池的工作电压,Ts为锂电池的表面温度,Vs为锂电池的体积,在本实施例中,取值Vs=1.65×10-5m3,α为锂电池电动势Eoc的温度影响系数,在本实施例中,取值α=-5×10-4V/K。
步骤1.2,定义Tia为锂电池内部温度Tin与外界环境温度Tamb之间的相对温度值,即Tia=Tin-Tamb,Tsa为锂电池表面温度与外界环境温度之间的相对温度值,即Tsa=Ts-Tamb,则步骤1.1中的热平衡微分方程(1)(2)可以变化为:
步骤1.3,选取x=[Tia Tsa]T为状态变量,生热速率Q为输入量,Tsa为输出量,根据步骤1.2中的式(4)和式(5),建立用于锂电池内部温度估计的连续状态空间模型:
y=Cx (7)
其中,为状态变量x的微分,A为连续状态空间模型的系数矩阵,B为连续状态空间模型的输入矩阵,C为连续状态空间模型的观测矩阵,C=[0 1],u为锂电池的输入,u=Q,y为锂电池的输出。
步骤1.4,对步骤1.3中的式(6)和式(7)进行离散化处理,并增加过程噪声w(k)以补偿系统建模误差、参数摄动、外部扰动,增加观测噪声v(k)以补偿系统输出测量噪声,建立用于锂电池内部温度估计的离散状态空间模型:
x(k+1)=Adx(k)+Bdu(k)+w(k) (8)
y(k)=Cdx(k)+v(k) (9)
其中,x(k)为锂电池在k时刻的状态变量,x(k)=[Tia(k) Tsa(k)]T,x(k+1)为锂电池在(k+1)时刻的状态变量,Ad为离散化后的系数矩阵,I为2阶单位矩阵,Bd为离散化后的输入矩阵,Cd为离散化后的观测矩阵,Cd=C=[0 1],Δt为采样时间,在本实施例中,取值Δt=0.001s,观测噪声v(k)有界,即v(k)满足|v(k)|≤v0,v0为上界参数,在本实施例中,取值v0=0.1。
步骤2,分别使用霍尔电压传感器和霍尔电流传感器实时采集工况下k时刻锂电池的工作电压U(k)和k时刻锂电池的工作电流I(k),使用热电偶温度传感器实时采集工况下k时刻锂电池表面温度Ts(k)和外界环境温度Tamb(k)。
步骤3,根据步骤2中采集到的锂电池工作电压U(k)、工作电流I(k)、锂电池表面温度Ts(k)和外界环境温度Tamb(k),在线辨识步骤1中的锂电池离散状态空间模型参数,包括锂电池内部与锂电池外壳之间的等效热阻Ri、锂电池外壳与外界环境之间的等效热阻Ro和锂电池内部的等效比热容Cc。
在线辨识锂电池离散状态空间模型参数的方法包括最小二乘法、卡尔曼滤波法、模糊控制法或神经网络法。
在本实施例中,采用最小二乘法在线辨识锂电池离散状态空间模型中的参数Ri、Ro和Cc,过程如下。
步骤3.1,根据步骤1.2中的式(4)和式(5),先联立消去Tia,再进行离散化处理,可以得到如下的数学表达式::
Tsa(k+2)=aTsa(k+1)+bTsa(k)+cQ(k) (10)
其中,Tsa(k)、Tsa(k+1)和Tsa(k+2)分别为锂电池在k、k+1和k+2时刻的温度值,Q(k)为锂电池在k时刻的生热速率,
令θ=[a b c]T,z(k)=Tsa(k+2),采用最小二乘法对θ=[a b c]T进行辨识的步骤为:
(1)初始化
当k=0和1时,
(2)循环递推
当k≥2时,
根据递推得到的系统参数的最小二乘估计的值、Cs的值以及采样时间Δt的值,结合a、b和c的表达式,计算Ri、Ro和Cc的值。
步骤4,基于步骤1中建立的锂电池离散状态空间模型构建离散滑模观测器,并结合步骤2中采集到的k时刻锂电池工作电压U(k)、工作电流I(k)、表面温度Ts(k)和外界环境温度Tamb(k)以及步骤3中所辨识的锂电池离散状态空间模型参数Ri、Ro和Cc,实时获取Tia(k)的估计值从而实时估计锂电池内部温度,锂电池内部温度Tin(k)的实时估计值为 其过程如下:
步骤4.1,基于步骤1中建立的锂电池离散状态空间模型构建离散滑模观测器:
其中,为x(k)的估计值,为x(k+1)的估计值,为y(k)的估计值,L为离散滑模观测器的增益矩阵,M为饱和函数增益矩阵,为饱和函数, 为边界层。
步骤4.2,将步骤2中采集到的k时刻锂电池工作电压U(k)、工作电流I(k)作为信号输入,利用步骤4.1中所构建的离散滑模观测器得到状态变量Tia(k)的实时估计值
步骤4.3,根据步骤4.2得到的状态变量Tia(k)的实时估计值和步骤3中采集到的k时刻外界环境温度Tamb(k),同时结合表达式实时估计锂电池内部温度,得到k时刻锂电池内部温度Tin(k)的实时估计值
在本实施例中,离散滑模观测器的增益矩阵L、饱和函数增益矩阵M、和边界层按照以下公式进行选取:
|w(k)-Lv(k)|≤Δ0 (13)
M≥Δ0 (15)
其中,Δ0为(w(k)-Lv(k))的上界参数,在本实施例中,取值v0为v(k)的上界参数,在本实施例中,取值v0=0.1。
Claims (7)
1.一种基于离散滑模观测器的锂电池内部温度估计方法,包括对工况下锂电池表面温度和外界环境温度的采集,其特征在于,主要步骤如下:
步骤1,根据锂电池内部与外壳之间的热平衡以及锂电池外壳与外界环境之间的热平衡,建立用于锂电池内部温度估计的离散状态空间模型:
x(k+1)=Adx(k)+Bdu(k)+w(k) (1)
y(k)=Cdx(k)+v(k) (2)
其中,
x(k)为锂电池在k时刻的状态变量,x(k)=[Tia(k) Tsa(k)]T,Tia(k)为k时刻锂电池内部温度Tin(k)与外界环境温度Tamb(k)之间的相对温度值,Tia(k)=Tin(k)-Tamb(k),Tsa(k)为k时刻锂电池表面温度Ts(k)与外界环境温度Tamb(k)之间的相对温度值,Tsa(k)=Ts(k)-Tamb(k);
x(k+1)为锂电池在(k+1)时刻的状态变量;
u(k)为锂电池在k时刻的输入,u(k)=Q(k),Q(k)为k时刻锂电池的生热速率;
y(k)为锂电池在k时刻的输出;
v(k)为锂电池在k时刻的观测噪声;
Ad为离散状态空间模型的系数矩阵,Ri为锂电池内部与锂电池外壳之间的等效热阻,Ro为锂电池外壳与外界环境之间的等效热阻,Cc为锂电池内部的等效比热容,Cs为锂电池外壳的比热容,Cs为常数,△t为采样时间;
Bd为离散状态空间模型的输入矩阵,Cc为锂电池内部的等效比热容,△t为采样时间;
Cd为离散状态空间模型的观测矩阵,Cd=[0 1],w(k)为锂电池在k时刻的过程噪声;
步骤2,分别通过电压传感器和电流传感器实时采集工况下k时刻锂电池的工作电压U(k)和k时刻锂电池的工作电流I(k),通过温度传感器实时采集工况下k时刻锂电池表面温度Ts(k)和外界环境温度Tamb(k);
步骤3,根据步骤2中采集到的k时刻锂电池工作电压U(k)、工作电流I(k)、表面温度Ts(k)和外界环境温度Tamb(k),在线辨识步骤1中锂电池离散状态空间模型中的参数,包括锂电池内部与锂电池外壳之间的等效热阻Ri、锂电池外壳与外界环境之间的等效热阻Ro和锂电池内部的等效比热容Cc;
步骤4,基于步骤1中建立的锂电池离散状态空间模型构建离散滑模观测器,并结合步骤2中采集到的k时刻锂电池工作电压U(k)、工作电流I(k)、表面温度Ts(k)和外界环境温度Tamb(k)以及步骤3中所辨识的锂电池离散状态空间模型参数Ri、Ro和Cc,实时获取Tia(k)的估计值从而实时估计锂电池内部温度,锂电池内部温度Tin(k)的实时估计值为
步骤4.1,构建离散滑模观测器的方程式如下:
其中,为x(k)的估计值,为x(k+1)的估计值,为y(k)的估计值,L为离散滑模观测器的增益矩阵,M为饱和函数增益矩阵,为饱和函数, 为边界层;
步骤4.2,将步骤2中采集到的k时刻锂电池工作电压U(k)、工作电流I(k)作为信号输入,利用步骤4.1中所构建的离散滑模观测器得到状态变量Tia(k)的实时估计值
步骤4.3,根据步骤4.2得到的状态变量Tia(k)的实时估计值和步骤3中采集到的k时刻外界环境温度Tamb(k),同时结合表达式实时估计锂电池内部温度,得到k时刻锂电池内部温度Tin(k)的实时估计值
2.根据权利要求1所述的一种基于离散滑模观测器的锂电池内部温度估计方法,其特征在于,步骤1中的k时刻锂电池的生热速率Q(k)为:
其中,I(k)为k时刻锂电池的工作电流,Vs为锂电池的体积,Eoc(k)为k时刻锂电池的电动势,U(k)为k时刻锂电池的工作电压,Ts(k)为k时刻锂电池的表面温度,α为锂电池电动势Eoc(k)的温度影响系数。
3.根据权利要求1所述的一种基于离散滑模观测器的锂电池内部温度估计方法,其特征在于,步骤1中的锂电池在k时刻的观测噪声v(k)有界,即满足以下表达式:
|v(k)|≤v0 (6)
其中,v0为v(k)上界参数,且v0为正实数。
4.根据权利要求1所述的一种基于离散滑模观测器的锂电池内部温度估计方法,其特征在于,步骤2中所述的电压传感器、电流传感器和温度传感器分别为霍尔电压传感器、霍尔电流传感器和热电偶。
5.根据权利要求1所述的一种基于离散滑模观测器的锂电池内部温度估计方法,其特征在于,步骤3中的在线辨识锂电池离散状态空间模型参数的方法为最小二乘法。
6.根据权利要求1所述的一种基于离散滑模观测器的锂电池内部温度估计方法,其特征在于,步骤1中的过程噪声w(k)和观测噪声v(k),以及步骤4中的离散滑模观测器的增益矩阵L满足以下表达式:
|w(k)-Lv(k)|≤△0 (7)
其中,△0为(w(k)-Lv(k))的上界参数,且△0为2×1的矩阵,且△0的每个元素都是正实数。
7.根据权利要求1所述的一种基于离散滑模观测器的锂电池内部温度估计方法,其特征在于,步骤4中的离散滑模观测器的增益矩阵L、饱和函数增益矩阵M、和边界层在进行选取时,满足以下表达式:
M≥△0 (9)
其中,△0为(w(k)-Lv(k))的上界参数,且△0为2×1的矩阵,且△0的每个元素都是正实数;v0是v(k)的上界参数,且v0为正实数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710098924.3A CN106872904B (zh) | 2017-02-23 | 2017-02-23 | 基于离散滑模观测器的锂电池内部温度估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710098924.3A CN106872904B (zh) | 2017-02-23 | 2017-02-23 | 基于离散滑模观测器的锂电池内部温度估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106872904A CN106872904A (zh) | 2017-06-20 |
CN106872904B true CN106872904B (zh) | 2019-02-26 |
Family
ID=59169548
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710098924.3A Active CN106872904B (zh) | 2017-02-23 | 2017-02-23 | 基于离散滑模观测器的锂电池内部温度估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106872904B (zh) |
Families Citing this family (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108241128B (zh) * | 2018-01-09 | 2019-10-01 | 西南交通大学 | 一种质子交换膜燃料电池系统状态估计方法 |
CN109188285B (zh) * | 2018-07-24 | 2020-12-29 | 合肥工业大学 | 一种锂电池内外热阻估计方法 |
US11740607B2 (en) * | 2018-08-30 | 2023-08-29 | Abb Schweiz Ag | Method and system for monitoring condition of electric drives |
DE102018214984A1 (de) * | 2018-09-04 | 2020-03-05 | Robert Bosch Gmbh | Verfahren zur Ermittlung einer Umgebungstemperatur einer ersten elektrischen Energiespeichereinheit im Verbund mit zweiten elektrischen Energiespeichereinheiten sowie entsprechende Vorrichtung, Computerprogramm und maschinenlesbares Speichermedium |
CN109596241B (zh) * | 2018-11-28 | 2020-07-14 | 同济大学 | 电机转子在线温度融合估计方法 |
CN109921111B (zh) * | 2019-03-14 | 2020-08-04 | 上海大学 | 一种锂离子电池内部温度估测方法及系统 |
CN110118617A (zh) * | 2019-05-30 | 2019-08-13 | 上海元城汽车技术有限公司 | 电池模组的内部温度确定方法、装置及智能终端 |
CN112578298B (zh) | 2019-09-29 | 2022-03-15 | 比亚迪股份有限公司 | 电池温度估算方法、装置、电子设备及存储介质 |
CN110826206B (zh) * | 2019-10-28 | 2024-03-01 | 上海交通大学 | 电池内部三维温度无损软测量方法及系统 |
CN110816366B (zh) * | 2019-10-31 | 2022-12-09 | 上海交通大学 | 适用于单体电池内部的温度估算方法、系统、介质及设备 |
CN112964991B (zh) * | 2019-11-28 | 2022-07-15 | 比亚迪股份有限公司 | 电池内部温度信息处理方法、计算机设备和存储介质 |
CN111143974B (zh) * | 2019-12-06 | 2022-08-12 | 重庆大学 | 一种面向控制的锂电池热模型建立方法 |
CN111595485B (zh) * | 2020-05-07 | 2022-03-04 | 广东工业大学 | 一种基于降阶模型的锂离子电池在线温度分布观测器设计方法 |
CN112526349B (zh) * | 2020-11-30 | 2021-11-12 | 北京航空航天大学 | 基于电池荷电状态与温度联合估计的整车能量管理方法 |
CN112926186B (zh) * | 2021-01-26 | 2022-05-13 | 同济大学 | 一种超级电容器单体内部温度在线估计方法 |
CN113484770B (zh) * | 2021-06-10 | 2022-04-01 | 广东恒翼能科技有限公司 | 基于充放电数据在线测算电池内部核心温度的方法及系统 |
CN113611947B (zh) * | 2021-08-02 | 2022-10-18 | 湖南工学院 | 一种电动汽车电池温度控制装置 |
CN113791364A (zh) * | 2021-08-31 | 2021-12-14 | 上海交通大学 | 模型融合信号驱动的锂离子电池内部短路诊断方法与系统 |
CN113608127B (zh) * | 2021-10-08 | 2021-12-21 | 湖南大学 | 一种圆柱形锂离子电池热异常检测与定位方法 |
CN113659246B (zh) * | 2021-10-20 | 2022-01-25 | 中国气象科学研究院 | 适应于极地超低温环境的电池系统及其控温方法 |
CN113824114B (zh) * | 2021-11-22 | 2022-03-18 | 广东电网有限责任公司惠州供电局 | 基于滑模观测的配电网状态估计方法、装置、系统及介质 |
CN114899523B (zh) * | 2022-05-18 | 2023-05-02 | 浙江大学 | 一种锂离子电池单体热失控核心温度估计方法 |
CN115113052A (zh) * | 2022-06-27 | 2022-09-27 | 重庆大学 | 在线检测锂离子电池温度场及内热源场的方法 |
CN116387653B (zh) * | 2023-06-05 | 2023-08-22 | 钛玛科(北京)工业科技有限公司 | 一种用于锂电池的低损耗管理方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5491420A (en) * | 1993-03-01 | 1996-02-13 | Duracell Inc. | Battery tester with stacked thermochromic elements |
WO2014095279A1 (de) * | 2012-12-21 | 2014-06-26 | Continental Automotive Gmbh | Verfahren zur ermittlung einer temperatur einer zelle einer batterie, ermittlungsvorrichtung und batterie |
CN104865534A (zh) * | 2015-04-29 | 2015-08-26 | 同济大学 | 一种单体电池内部温度估计方法 |
CN105264709A (zh) * | 2013-06-14 | 2016-01-20 | Hrl实验室有限责任公司 | 用于感测电化学装置的内部温度的方法和设备 |
CN106324523A (zh) * | 2016-09-26 | 2017-01-11 | 合肥工业大学 | 基于离散变结构观测器的锂电池soc估计方法 |
-
2017
- 2017-02-23 CN CN201710098924.3A patent/CN106872904B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5491420A (en) * | 1993-03-01 | 1996-02-13 | Duracell Inc. | Battery tester with stacked thermochromic elements |
WO2014095279A1 (de) * | 2012-12-21 | 2014-06-26 | Continental Automotive Gmbh | Verfahren zur ermittlung einer temperatur einer zelle einer batterie, ermittlungsvorrichtung und batterie |
CN105264709A (zh) * | 2013-06-14 | 2016-01-20 | Hrl实验室有限责任公司 | 用于感测电化学装置的内部温度的方法和设备 |
CN104865534A (zh) * | 2015-04-29 | 2015-08-26 | 同济大学 | 一种单体电池内部温度估计方法 |
CN106324523A (zh) * | 2016-09-26 | 2017-01-11 | 合肥工业大学 | 基于离散变结构观测器的锂电池soc估计方法 |
Non-Patent Citations (1)
Title |
---|
"基于离散滑模观测器的锂电池荷电状态估计";孙冬等;《中国电机工程学报》;20150105;第35卷(第1期);第185-191页 |
Also Published As
Publication number | Publication date |
---|---|
CN106872904A (zh) | 2017-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106872904B (zh) | 基于离散滑模观测器的锂电池内部温度估计方法 | |
CN108919137B (zh) | 一种考虑不同电池状态的电池老化状态估计方法 | |
CN105319515B (zh) | 锂离子电池荷电状态和健康状态联合估算方法 | |
CN111581904B (zh) | 考虑循环次数影响的锂电池soc及soh协同估算方法 | |
CN106716158B (zh) | 电池荷电状态估算方法和装置 | |
CN103003709B (zh) | 电池状态估计装置以及电池状态估计方法 | |
CN111308379B (zh) | 一种基于局部恒压充电数据的电池健康状态估计方法 | |
CN104931784B (zh) | 基于频率的电池模型参数估计 | |
Hu et al. | Online model identification of lithium-ion battery for electric vehicles | |
CN105304952A (zh) | 用于车辆控制的降阶电化学电池模型 | |
CN111693877B (zh) | 一种锂离子电池的soc在线估测方法和系统 | |
CN110596606B (zh) | 一种锂电池剩余电量估计方法、系统及装置 | |
CN104849675A (zh) | 锂离子电池电化学和热耦合模型的获取方法 | |
CN105320033A (zh) | 用于车辆控制的温度相关的电化学电池模型 | |
CN105866700B (zh) | 一种锂离子电池快速筛选的方法 | |
CN103983920B (zh) | 一种建立电动车辆的动力电池的模型的方法 | |
JP2011521402A (ja) | 電気化学電池のスマート管理システム | |
CN114740386B (zh) | 一种基于健康状态的锂离子电池荷电状态估计方法 | |
CN110221221A (zh) | 锂离子电池荷电状态和健康状态联合估计方法 | |
CN105717460A (zh) | 一种基于非线性观测器的动力电池soc估计方法和系统 | |
CN109143097A (zh) | 一种计及温度和循环次数的锂离子电池soc估计方法 | |
CN113504481B (zh) | 一种电动汽车动力电池多状态联合估计系统及方法 | |
CN114729967A (zh) | 用于估计储能器的状态的方法 | |
CN105291873A (zh) | 在用于车辆控制的电池模型中的金属离子浓度的插值 | |
CN109828215A (zh) | 一种提升电池单体soc估算精度的方法和系统 |
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 |