CN110412483A - 一种基于互联观测器的车用锂电池荷电状态估计方法 - Google Patents
一种基于互联观测器的车用锂电池荷电状态估计方法 Download PDFInfo
- Publication number
- CN110412483A CN110412483A CN201910370131.1A CN201910370131A CN110412483A CN 110412483 A CN110412483 A CN 110412483A CN 201910370131 A CN201910370131 A CN 201910370131A CN 110412483 A CN110412483 A CN 110412483A
- Authority
- CN
- China
- Prior art keywords
- lithium battery
- state
- observer
- charge
- equation
- 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.)
- Pending
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60L—PROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
- B60L58/00—Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles
- B60L58/10—Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles for monitoring or controlling batteries
- B60L58/12—Methods or circuit arrangements for monitoring or controlling batteries or fuel cells, specially adapted for electric vehicles for monitoring or controlling batteries responding to state of charge [SoC]
-
- 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
-
- 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/385—Arrangements for measuring battery or accumulator variables
- G01R31/387—Determining ampere-hour charge capacity or SoC
-
- 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/396—Acquisition or processing of data for testing or for monitoring individual cells or groups of cells within a battery
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/60—Other road transportation technologies with climate change mitigation effect
- Y02T10/70—Energy storage systems for electromobility, e.g. batteries
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Power Engineering (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Secondary Cells (AREA)
Abstract
一种基于互联观测器的车用锂电池荷电状态估计方法,步骤包括:步骤1、构建锂电池的动力学方程,建立一种考虑液相动力学的锂电池扩展单粒子模型;步骤2、提出一种可同时实现对电池正、负电极浓度分布估计的锂电池双向互联观测器;步骤3、结合步骤1中的锂电池扩展单粒子模型和步骤2中所构建的互联观测器,提出并验证方法;本发明所提出的基于互联观测器的车用锂电池荷电状态估计方法能够同时对电池正、负电极的荷电状态进行估计,具有工程应用价值;解决了现有基于电化学模型的车用锂电池荷电状态估计仅能对一个电极荷电状态进行估计,同时锂电池出现老化时无法保证正、负电极锂离子摩尔数仍守恒,从而出现锂电池荷电状态估计误差较大的问题。
Description
技术领域
本发明涉及新能源汽车电池管理技术领域,具体涉及一种基于互联观测器的车用锂电池荷电状态估计方法。
背景技术
近年来,锂电池由于高能量密度、高输出功率及长寿命等优点在新能源汽车中具有广泛应用。作为构成锂离子动力电池组的重要单元,锂电池的精确建模和荷电状态估计对于开发高效的电池管理系统 (batterymanagement system,BMS)具有十分重要的意义。
现有技术中,基于等效电路模型(equivalent-circuitmodel,ECM) 的锂电池建模和内部参数识别方法以及不同的荷电状态估计方法比较常用。然而,ECM模型利用电阻、电容等元器件模拟锂电池输出电压,对于前期电池的实验具有很强的依赖性,且模型参数也不能完全对应锂电池内部实际的物理电化学状态变量。因此,为了更准确的表征锂电池内部的电化学行为与物理特性,需要构建一种基于电化学模型来对锂电池的荷电状态进行估计。
尽管已有基于电化学模型的锂电池荷电状态估计方法,但均假定锂电池正、负电极的锂离子摩尔数是恒定的,这使得两个电极的锂浓度状态存在一定的代数关系,从而使得该类电化学模型仅能对单个电极荷电状态是可观测的。然而当锂电池出现老化时,锂电池正、负极锂离子浓度的守恒就难以保证,这将导致锂电池的荷电状态估计出现较大误差。
综上所述,提出一种能够同时了解锂电池正、负电极的荷电状态估计方法是目前新能源汽车电池管理技术领域中研究的关键性问题。
发明内容
为克服上述现有技术的不足,本发明的目的在于提供一种基于互联观测器的车用锂电池荷电状态估计方法,基于多孔电极理论和浓度理论提出一种考虑液相动力学行为的锂电池扩展单粒子模型(ESPM 是enhanced single particle model的英文简称),在该模型基础上进一步简化,提出该方法可实现同时对锂电池的正、负电极荷电状态进行估计,有效提高了锂电池管理系统荷电状态估计的准确性;本发明解决了现有基于电化学模型的车用锂电池荷电状态估计仅能对一个电极荷电状态进行估计,同时锂电池出现老化时无法保证正、负电极锂离子摩尔数仍守恒,出现锂电池荷电状态估计误差较大的问题。
为实现上述目的,本发明采用的技术方案是:
一种基于互联观测器的车用锂电池荷电状态估计方法,其特征在于,包括以下步骤:
步骤1,构建锂电池的动力学方程,建立一种考虑液相动力学的锂电池扩展单粒子模型ESPM;
步骤2,提出一种可同时实现对正负电极浓度分布的锂电池双向互联观测器;
步骤3,结合步骤1中的锂电池扩展单粒子模型ESPM和步骤2 中所构建的互联观测器,提出并验证一种基于互联观测器的车用锂电池的荷电状态估计方法。
步骤1具体做法是:
锂电池是一种恒流等温电化学模型,构建一种考虑液相动力学行为的锂电池ESPM模型,其中主要包括的锂电池电化学动力学方程有:固相扩散方程、液相扩散方程、液相电荷守恒方程、固相电荷守恒方程、Bulter-Volmer动力学方程;应用有限差分法数值求解动力学方程,构建锂电池输出电压及荷电状态表达式分别为:
上式中,V(t)为电池端电压,为锂电池正电极固相电势,为锂电池负极固相电势,I(t)为电流输出电流密度,Rc为锂电池集电极接触电阻。SOC(t)为锂电池的电荷状态,Rs为固相粒子半径,为固相最大锂离子浓度。
步骤2具体做法是:
构建一种基于H∞理论框架的锂电池正、负极互联观测器,其结构不同于只能对单电极状态进行估计的观测器,它能够同时估计正、负极的浓度分布状态并使其相互提供反馈;为了降低锂电池H∞互联观测器的设计复杂度并满足其可观测性条件,给出相对应的假设条件。
步骤3具体做法是:
结合步骤1的扩展单粒子模型ESPM和步骤2中基于H∞互联观测器Obsv-2,用本发明提出的基于互联观测器的车用锂电池荷电状态方法对锂电池进行荷电状态估计,其结果在不同工况下与单电极观测器Obsv-1所计算的荷电状态估计进行对比分析,验证有效性。
所述的锂电池输出电压及荷电状态表达式的计算过程如下:
根据Fick第二定理,锂电池正、负电极的锂离子浓度为:
其边界控制条件包括:
其中,正、负极的反应电流密度表达式为:
在笛卡尔坐标系x轴上,ce(x,t)随着锂离子的流量密度梯度而变化,其动力学方程为:
液相浓度扩散方程的边界控制条件和浓度扩散连续条件如式(6) 和式(7):
式中,液相有效扩散系数
在x轴上只需要考虑电解液内锂离子传输引起的电势变化,则根据修正的欧姆定律,电解液的液相电势表示为:
式中,
沿x轴在锂电池厚度范围内对上式(8)积分两次可得液相电势差ΔΦe为
根据Bulter-Volmer方程以及式(4)的假设,可知η±(x,t)的计算式为:
取正负极电荷传输系数αa=αc=α=0.5,则交换电流密度定义为
又因为:
综合式(1),(9),(10)和(12),锂电池ESPM模型的终端电压V(t) 计算式为
为了便于标记和后续观测器设计,定义非线性函数(h(·)),则可得如下函数式:
因此,式(13)可简化为:
与温度变化相关的锂电池电化学参数包括和k±(T),这4个参数可以采用Arrhenius定律来标定温度对其的影响,表示为
式中,Tref=23℃(298K)是参考温度,和为活化能参数,此外与锂电池温度T和初始液相浓度有关的液相扩散系数De(T),以及依赖于温度的电解液离子电导率ke(T)的经验计算式为:
为了能够应用有限差分法数值求解上述模型,将偏微分方程(2) 中的固相浓度扩散方程在球形离子内部离散化为N+1个节点,同样地, 沿x轴将方程(5)离散化为M+1个节点,
定义状态向量x=[x1 x2 x3]T,其中,x1和x2分别是正、负电极的第N个节点的浓度状态变量,表示为 x3是电解液第M个节点的浓度状态变量, x3=[ce,1,ce,1,…,ce,M]T,y=V(t),
需要指出的是,和分别是正、负电极表面固相浓度,即和经过离散化之后,可得2N+M个常微分方程(ordinary differential equation,ODE),其状态空间方程的形式如下:
其中:u=I(t),状态变量系数矩阵列向量矩阵分别是包含电化学模型参数的状态变量和输入变量矩阵,注意,x1和x2分别包含正、负电极所有离散化节点的浓度值,若能求解方程(18)获得其数值解,则锂电池的荷电状态 (SOC)计算式为:
本发明所用的相关变量及参数所代表的含义如表1所示:
表1锂电池参数及相关符号
与现有技术相比,本发明的有益效果在于:
1)、与现有技术中基于正、负极锂离子守恒原理的单电极观测器荷电状态估计方法不同,基于H∞鲁棒控制理论提出的一种基于互联观测器的车用锂电池荷电状态估计方法,能够实现同时对锂电池正、负电极的荷电状态进行估计。
2)、即使本发明中提出的互联观测器结构中出现不准确的初始值,但本发明中提出的锂电池荷电状态估计方法仍能够使锂电池双电极的荷电状态估计值保持收敛。
3)、相比于HPPC和UDDS工况下的试验数据,对比分析了Obsv-1 以及Obsv-2模型计算出的锂电池在两种工况下的荷电状态响应曲线,结果发现:在HPPC工况下应用两种模型计算出的锂电池荷电状态估计值的最大相对误差分别为4.7%和3.4%;在UDDS工况下,应用两种模型计算出的锂电池荷电状态估计值最大相对误差分别为4.4%和 3.2%。
因此本发明所提出的一种基于互联观测器的车用锂电池荷电状态估计方法具有很好的适用性和工程价值。
附图说明
图1为本发明锂电池扩展单粒子模型(ESPM)示意图。
图2(a)为现有技术中锂电池的单电极浓度观测器。
图2(b)本发明的锂电池的互联观测器。
图3为HPPC工况下不同模型计算的锂电池(a)SOC估计值; (b)SOC误差曲线。
图4为UDDS工况下不同模型计算的锂电池(a)SOC估计值; (b)SOC误差曲线。
图5为UDDS工况下不同模型计算的锂电池SOC估计曲线。
图6为UDDS工况下不同模型计算的锂电池SOC误差曲线。
具体实施方式
为了便于理解本发明,下文将结合较佳的实施例对本发明作更全面、细致地描述,但本发明的保护范围并不限于以下具体的实施例。
一种基于互联观测器的车用锂电池荷电状态估计方法,其特征在于,包括以下步骤:
步骤1,构建锂电池的动力学方程,建立一种考虑液相动力学的锂电池扩展单粒子模型ESPM;
步骤2,提出一种可同时实现对正负电极浓度分布的锂电池双向互联观测器;
步骤3,结合步骤1中的锂电池扩展单粒子模型ESPM和步骤2 中所构建的互联观测器,提出并验证一种基于互联观测器的车用锂电池的荷电状态估计方法。
步骤1具体做法是:
锂电池是一种恒流等温电化学模型,构建一种考虑液相动力学行为的锂电池ESPM模型,其中主要包括的锂电池电化学动力学方程有:固相扩散方程、液相扩散方程、液相电荷守恒方程、固相电荷守恒方程、Bulter-Volmer动力学方程(这些方程是现有技术);应用有限差分法数值求解动力学方程,构建锂电池输出电压及荷电状态表达式分别为:
步骤2具体做法是:
提出一种基于H∞理论框架的锂电池正、负极互联观测器,其结构不同于只能对单电极状态进行估计的观测器,它但能同时估计正、负极的浓度分布状态并使其相互提供反馈;为了降低锂电池H∞互联观测器的设计复杂度并满足其可观测性条件,给出相对应的假设条件。
步骤3具体做法是:
结合步骤1的扩展单粒子模型ESPM和步骤2中基于H∞互联观测器Obsv-2,用本发明提出的基于互联观测器的车用锂电池荷电状态方法对锂电池进行荷电状态估计,其结果在不同工况下与单电极观测器Obsv-1所计算的荷电状态估计进行对比分析,验证有效性。
实施例1
本发明以索尼VTC4 18650锂电池为研究对象,其标称容量和电压分别为2Ah和3.6V,最大输出电压为4.2V,最小截止电压为2.5V。由此,本发明所提出的一种基于互联观测器的车用锂电池荷电状态估计方法,具体实施方式如下所述:
步骤1,由于电化学机理模型可以描述电池内部的电化学行为与物理特性,使其可以更准确地描述电池的运行状态,基于多孔电极理论和浓度理论,构建一种考虑液相动力学行为的锂电池扩展单粒子模型(ESPM),如图1所示,考虑到正、负极集流体产生的欧姆电势差,则锂电池终端电压V(t)的计算式为:
该模型考虑了电解液相动力学行为,耦合了温度和液相浓度变化对锂离子电池关键参数的影响,主要包括的锂电池电化学动力学方程有:固相扩散方程、液相扩散方程、液相电荷守恒方程、固相电荷守恒方程、Bulter-Volmer动力学方程,具体过程如下:
根据Fick第二定理,锂电池正、负电极的锂离子浓度为
其边界控制条件包括:
其中,正、负极的反应电流密度表达式为
在图1所示的笛卡尔坐标系x轴上,ce(x,t)随着锂离子的流量密度梯度而变化,其动力学方程为
液相浓度扩散方程的边界控制条件和浓度扩散连续条件如式(6) 和式(7):
式中,液相有效扩散系数
在x轴上只需要考虑电解液内锂离子传输引起的电势变化,则根据修正的欧姆定律,电解液的液相电势表示为
式中,
沿x轴在锂电池厚度范围内对上式(8)积分两次可得液相电势差ΔΦe为
根据Bulter-Volmer方程以及式(4)的假设,可知η±(x,t)的计算式为:
取正负极电荷传输系数αa=αc=α=0.5,则交换电流密度定义为
又因为:
综合式(1),(9),(10)和(12),锂电池ESPM模型的终端电压V(t) 计算式为
为了便于标记和后续观测器设计,定义非线性函数(h(·)),则可得如下函数式:
因此,式(13)可简化为
本发明中,与温度变化相关的锂电池电化学参数包括和k±(T),这4个参数可以采用Arrhenius定律来标定温度对其的影响,表示为
式中,Tref=23℃(298K)是参考温度,和为活化能参数。此外与锂电池温度T和初始液相浓度有关的液相扩散系数De(T),以及依赖于温度的电解液离子电导率ke(T)的经验计算式为:
为了能够应用有限差分法数值求解上述模型,将偏微分方程(2) 中的固相浓度扩散方程在球形离子内部离散化为N+1个节点,同样地, 沿x轴将方程(5)离散化为M+1个节点。
定义状态向量x=[x1 x2 x3]T,其中,x1和x2分别是正、负电极的第N个节点的浓度状态变量,表示为 x3是电解液第M个节点的浓度状态变量, x3=[ce,1,ce,1,…,ce,M]T,y=V(t)。
需要指出的是,和分别是正、负电极表面固相浓度,即和经过离散化之后,可得2N+M个常微分方程(ordinary differential equation,ODE),其状态空间方程的形式如下:
其中:u=I(t),状态变量系数矩阵列向量矩阵分别是包含电化学模型参数的状态变量和输入变量矩阵。注意,x1和x2分别包含正、负电极所有离散化节点的浓度值,若能求解方程(18)获得其数值解,则锂电池的荷电状态 (SOC)计算式为:
本发明所用的相关变量及参数所代表的含义如表1所示。
表1锂电池参数及相关符号
步骤2,基于H∞鲁棒控制理论提出一种可同时实现对正负电极浓度分布的锂电池双向互联观测器,具体过程如下:
从现有技术基于电化学模型的荷电状态估计方法中可知,如果单个电极的浓度状态可根据其开环模型的计算得到的话,则另一个电极(正或负)的状态是可观测的,该观测器结构如图2(a)所示;但是在该观测器结构中,不准确的初始值会导致开环模型的响应误差非常大,那么模型计算出的锂电池输出电压就是不准确的,进而导致锂电池的单电极浓度状态估计偏离于真实值,
为了解决该问题,本发明提出一种基于H∞理论框架的锂电池正、负极互联观测器,其结构类似现有技术中存在的单个电极观测器,但用此观测器能同时对锂电池正、负极的浓度状态进行估计并使其相互提供反馈,即本发明所提出的互联观测器结构如图2(b)所示;
首先,为降低锂电池H∞互联观测器的设计复杂度并满足其可观测性条件,假设锂电池电解液中的浓度和电势梯度是均匀分布的,即分别忽略偏微分方程(5)和(8),则可以从锂电池ESPM模型推导出单粒子模型(single-particle-model,SPM),SPM是一种简化的适合于观测器设计的锂电池电化学模型。若考虑系统输入高斯白噪声w(t),则锂电池SPM模型的状态方程可表示为:
式中:
对于方程(20)所描述的锂电池SPM模型,建立如下的互联观测器方程:
式中,L1和是H∞状态观测器的增益,角标ol代表开环模型的状态变量。
接着,定义系统状态误差为和锂电池正、负极的输出误差为同样地,
另外假设锂电池正、负极表面浓度的误差,e1,N和e2,N与状态误差向量e1和e2满足的关系式为e1,N=Ce1,e2,N=Ce2,其中C=[0,0,...,0,1]为满足输出电压的分布矩阵,且
由此根据方程(20),(21)和(22),可得正、负极状态观测器的误差动力学方程分别为
和
从而将系统(20)的状态误差问题归结为:求得输出误差反馈增益矩阵L1和L2,以保证在特定的衰减系数γ>0下,误差状态方程(23)和 (24)在w(t)下是渐进稳定的,即且在零初始条件下对所有满足:
本发明将H∞互联状态观测器的设计转化为寻找合适的对称正定矩阵P1>0,P2>0和适当维数矩阵Y1和Y2,进而获得状态观测器增益矩阵 L1和L2,使得下面的线性矩阵不等式(Linear Matrix Inequality,LMI)成立:
为了将式(25)转化为一个线性矩阵不等式,定义Y1=P1L1(即),Y2=P2L2(即),于是,不等式(26)和(27)可进一步转化为:
本发明利用MATLAB中的LMI工具箱求解式(28)和(29)可得 P1,Y1和P2,Y2,进而确定控制增益L1和L2;
步骤3,结合步骤1中的扩展单粒子模型和步骤2中基于H∞互联观测器(Obsv-2),提出一种基于互联观测器的车用锂电池荷电状态估计方法,从而对车用锂电池的荷电状态进行估计,其结果在不同工况下与现有技术中的单电极观测器(Obsv-1)所计算的荷电状态估计进行对比分析,验证本发明所提出的一种基于互联观测器的车用锂电池荷电状态估计方法的有效性;
本发明以索尼VTC4 18650锂电池为研究对象,其标称容量和电压分别为2Ah和3.6V;根据上述步骤2中提出的基于H∞互联观测器,对锂电池的双电极的浓度状态分布进行估计,然后结合公式(19)中锂电池的电极荷电状态表达式,提出本发明中的一种基于互联观测的车用锂电池的荷电状态估计方法;
为了验证本发明所提出一种基于互联观测器的车用锂电池荷电状态估计方法的有效性和准确性,在HPPC和UDDS工况下,对基于H∞互联观测器(Obsv-2)的锂电池荷电状态估计结果与现有技术中的单电极观测器(Obsv-1)所计算的荷电状态估计进行对比分析;其中的SOCref为应用安时积分法得到的SOC估计值,并将其作为SOC 参考值;另外,SOC误差计算式分别为:SOCerror-1=SOCref–SOCObsv-1, SOCerror-2=SOCref–SOCObsv-2;结果发现:当锂电池在室温23℃时,在 HPPC工况下应用两种模型计算出的锂电池荷电状态估计值的最大相对误差分别为4.7%和3.4%,具体如图3-4所示;在UDDS工况下,应用两种模型计算出的锂电池荷电状态估计值最大相对误差分别为 4.4%和3.2%,具体如图5-6所示;
本发明对现有技术中的单电极观测器Obsv-1、本发明中提出的互联观测器Obsv-2所计算的和荷电状态进行了仿真和对比分析,证明了本发明所提出的一种基于互联观测器的车用锂电池荷电状态估计方法的有效性,进而提高了锂电池管理系统中荷电状态估计的准确性,具有一定的工程应用价值。
Claims (5)
1.一种基于互联观测器的车用锂电池荷电状态估计方法,其特征在于,包括以下步骤:
步骤1,构建锂电池的动力学方程,建立一种考虑液相动力学的锂电池扩展单粒子模型ESPM;
步骤2,提出一种可同时实现对正负电极浓度分布的锂电池双向互联观测器;
步骤3,结合步骤1中的锂电池扩展单粒子模型ESPM和步骤2中所构建的互联观测器,提出并验证一种基于互联观测器的车用锂电池的荷电状态估计方法。
2.根据权利要求1所述的一种基于互联观测器的车用锂电池荷电状态估计方法,其特征在于,步骤1具体做法是:
锂电池是一种恒流等温电化学模型,构建一种考虑液相动力学行为的锂电池ESPM模型,其中主要包括的锂电池电化学动力学方程有:固相扩散方程、液相扩散方程、液相电荷守恒方程、固相电荷守恒方程、Bulter-Volmer动力学方程;应用有限差分法数值求解动力学方程,构建出锂电池输出电压及荷电状态表达式分别为:
上式中,V(t)为电池端电压,为锂电池正电极固相电势,为锂电池负极固相电势,I(t)为电流输出电流密度,Rc为锂电池集电极接触电阻。SOC(t)为锂电池的电荷状态,Rs为固相粒子半径,为固相最大锂离子浓度。
3.根据权利要求1所述的一种基于互联观测器的车用锂电池荷电状态估计方法,其特征在于,步骤2具体做法是:
构建一种基于H∞理论框架的锂电池正、负极互联观测器,其结构不同于只能对单电极状态进行估计的观测器,它能同时估计正、负极的浓度分布状态并使其相互提供反馈;为了降低锂电池H∞互联观测器的设计复杂度并满足其可观测性条件,给出相对应的假设条件。
4.根据权利要求1所述的一种基于互联观测器的车用锂电池荷电状态估计方法,其特征在于,步骤3具体做法是:
结合步骤1的扩展单粒子模型ESPM和步骤2中基于H∞互联观测器Obsv-2,用本发明提出的基于互联观测器的车用锂电池荷电状态方法对锂电池进行荷电状态估计,其结果在不同工况下与单电极观测器Obsv-1所计算的荷电状态估计进行对比分析,验证有效性。
5.根据权利要求2所述的一种基于互联观测器的车用锂电池荷电状态估计方法,其特征在于,所述的锂电池输出电压及荷电状态表达式的计算过程如下:
根据Fick第二定理,锂电池正、负电极的锂离子浓度为:
其边界控制条件包括:
其中,正、负极的反应电流密度表达式为:
在笛卡尔坐标系x轴上,ce(x,t)随着锂离子的流量密度梯度而变化,其动力学方程为:
液相浓度扩散方程的边界控制条件和浓度扩散连续条件如式(6)和式(7):
式中,液相有效扩散系数
在x轴上只需要考虑电解液内锂离子传输引起的电势变化,则根据修正的欧姆定律,电解液的液相电势表示为:
式中,
沿x轴在锂电池厚度范围内对上式(8)积分两次可得液相电势差ΔΦe为
根据Bulter-Volmer方程以及式(4)的假设,可知η±(x,t)的计算式为:
取正负极电荷传输系数αa=αc=α=0.5,则交换电流密度定义为
又因为:
综合式(1),(9),(10)和(12),锂电池ESPM模型的终端电压V(t)计算式为
为了便于标记和后续观测器设计,定义非线性函数(h(·)),则可得如下函数式:
因此,式(13)可简化为:
与温度变化相关的锂电池电化学参数包括和k±(T),这4个参数可以采用Arrhenius定律来标定温度对其的影响,表示为
式中,Tref=23℃(298K)是参考温度,和为活化能参数,此外与锂电池温度T和初始液相浓度有关的液相扩散系数De(T),以及依赖于温度的电解液离子电导率ke(T)的经验计算式为:
为了能够应用有限差分法数值求解上述模型,将偏微分方程(2)中的固相浓度扩散方程在球形离子内部离散化为N+1个节点,同样地,沿x轴将方程(5)离散化为M+1个节点,
定义状态向量x=[x1 x2 x3]T,其中,x1和x2分别是正、负电极的第N个节点的浓度状态变量,表示为 x3是电解液第M个节点的浓度状态变量,x3=[ce,1,ce,1,…,ce,M]T,y=V(t),
需要指出的是,和分别是正、负电极表面固相浓度,即和经过离散化之后,可得2N+M个常微分方程(ordinary differential equation,ODE),其状态空间方程的形式如下:
其中:u=I(t),状态变量系数矩阵A11,列向量矩阵B1,分别是包含电化学模型参数的状态变量和输入变量矩阵,注意,x1和x2分别包含正、负电极所有离散化节点的浓度值,若能求解方程(18)获得其数值解,则锂电池的荷电状态(SOC)计算式为:
本发明所用的相关变量及参数所代表的含义如表1所示:
表1锂电池参数及相关符号
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910370131.1A CN110412483A (zh) | 2019-05-06 | 2019-05-06 | 一种基于互联观测器的车用锂电池荷电状态估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910370131.1A CN110412483A (zh) | 2019-05-06 | 2019-05-06 | 一种基于互联观测器的车用锂电池荷电状态估计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110412483A true CN110412483A (zh) | 2019-11-05 |
Family
ID=68358176
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910370131.1A Pending CN110412483A (zh) | 2019-05-06 | 2019-05-06 | 一种基于互联观测器的车用锂电池荷电状态估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110412483A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110888056A (zh) * | 2019-11-20 | 2020-03-17 | 上海交通大学 | 适用于车载动力锂离子电池的在线soc观测器搭建方法及系统 |
CN111753402A (zh) * | 2020-05-29 | 2020-10-09 | 天津大学 | 燃料电池汽车辅助动力源锂电池的建模方法 |
CN111929581A (zh) * | 2020-06-05 | 2020-11-13 | 西安理工大学 | 一种动力锂电池内外部温度预测方法 |
CN114089191A (zh) * | 2021-11-17 | 2022-02-25 | 浙大城市学院 | 一种复合锂离子电池健康状况估计方法 |
CN114757026A (zh) * | 2022-04-08 | 2022-07-15 | 昆明理工大学 | 一种全工况多尺度动力锂电池电化学耦合建模方法 |
CN115236516A (zh) * | 2022-06-17 | 2022-10-25 | 上海玫克生储能科技有限公司 | 基于电化学模型的锂电池的预警方法 |
WO2023030024A1 (zh) * | 2021-09-01 | 2023-03-09 | 上海屹锂新能源科技有限公司 | 基于电化学模型的固态锂电池状态估计方法及系统 |
CN115935635A (zh) * | 2022-11-29 | 2023-04-07 | 上海玫克生储能科技有限公司 | 基于电化学模型的锂电池路端电压计算方法、装置及介质 |
CN116660752A (zh) * | 2023-05-29 | 2023-08-29 | 广州汽车集团股份有限公司 | 电池界面反应电压修正方法、装置、可读介质及电子设备 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107066722A (zh) * | 2017-04-06 | 2017-08-18 | 北京理工大学 | 一种基于电化学模型的动力电池系统荷电状态和健康状态的联合估计方法 |
-
2019
- 2019-05-06 CN CN201910370131.1A patent/CN110412483A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107066722A (zh) * | 2017-04-06 | 2017-08-18 | 北京理工大学 | 一种基于电化学模型的动力电池系统荷电状态和健康状态的联合估计方法 |
Non-Patent Citations (2)
Title |
---|
庞辉 等: "一种基于简化电化学模型的锂电池互联状态观测器", 《物理学报》 * |
庞辉: "基于电化学模型的锂离子电池多尺度建模及其简化方法", 《物理学报》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110888056A (zh) * | 2019-11-20 | 2020-03-17 | 上海交通大学 | 适用于车载动力锂离子电池的在线soc观测器搭建方法及系统 |
CN110888056B (zh) * | 2019-11-20 | 2022-05-06 | 上海交通大学 | 适用于车载动力锂离子电池的在线soc观测器搭建方法及系统 |
CN111753402A (zh) * | 2020-05-29 | 2020-10-09 | 天津大学 | 燃料电池汽车辅助动力源锂电池的建模方法 |
CN111753402B (zh) * | 2020-05-29 | 2022-03-29 | 天津大学 | 燃料电池汽车辅助动力源锂电池的建模方法 |
CN111929581B (zh) * | 2020-06-05 | 2022-10-21 | 西安理工大学 | 一种动力锂电池内外部温度预测方法 |
CN111929581A (zh) * | 2020-06-05 | 2020-11-13 | 西安理工大学 | 一种动力锂电池内外部温度预测方法 |
WO2023030024A1 (zh) * | 2021-09-01 | 2023-03-09 | 上海屹锂新能源科技有限公司 | 基于电化学模型的固态锂电池状态估计方法及系统 |
CN114089191A (zh) * | 2021-11-17 | 2022-02-25 | 浙大城市学院 | 一种复合锂离子电池健康状况估计方法 |
CN114757026A (zh) * | 2022-04-08 | 2022-07-15 | 昆明理工大学 | 一种全工况多尺度动力锂电池电化学耦合建模方法 |
CN114757026B (zh) * | 2022-04-08 | 2024-06-07 | 昆明理工大学 | 一种全工况多尺度动力锂电池电化学耦合建模方法 |
CN115236516A (zh) * | 2022-06-17 | 2022-10-25 | 上海玫克生储能科技有限公司 | 基于电化学模型的锂电池的预警方法 |
CN115935635A (zh) * | 2022-11-29 | 2023-04-07 | 上海玫克生储能科技有限公司 | 基于电化学模型的锂电池路端电压计算方法、装置及介质 |
CN115935635B (zh) * | 2022-11-29 | 2024-02-27 | 上海玫克生储能科技有限公司 | 基于电化学模型的锂电池路端电压计算方法、装置及介质 |
CN116660752A (zh) * | 2023-05-29 | 2023-08-29 | 广州汽车集团股份有限公司 | 电池界面反应电压修正方法、装置、可读介质及电子设备 |
CN116660752B (zh) * | 2023-05-29 | 2024-03-19 | 广州汽车集团股份有限公司 | 电池界面反应电压修正方法、装置、可读介质及电子设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110412483A (zh) | 一种基于互联观测器的车用锂电池荷电状态估计方法 | |
Zhang et al. | Lithium-ion battery pack state of charge and state of energy estimation algorithms using a hardware-in-the-loop validation | |
Sun et al. | Combined internal resistance and state-of-charge estimation of lithium-ion battery based on extended state observer | |
Tang et al. | Model migration based battery power capability evaluation considering uncertainties of temperature and aging | |
Tian et al. | Fractional-order model-based incremental capacity analysis for degradation state recognition of lithium-ion batteries | |
Ahmed et al. | Reduced-order electrochemical model parameters identification and soc estimation for healthy and aged li-ion batteries part i: Parameterization model development for healthy batteries | |
Kim | Nonlinear state of charge estimator for hybrid electric vehicle battery | |
Kim et al. | An on-board model-based condition monitoring for lithium-ion batteries | |
Tanim et al. | State of charge estimation of a lithium ion cell based on a temperature dependent and electrolyte enhanced single particle model | |
Allam et al. | An interconnected observer for concurrent estimation of bulk and surface concentration in the cathode and anode of a lithium-ion battery | |
Zhang et al. | Robust and adaptive estimation of state of charge for lithium-ion batteries | |
Corno et al. | Electrochemical model-based state of charge estimation for Li-ion cells | |
Ahmed et al. | Reduced-order electrochemical model parameters identification and state of charge estimation for healthy and aged Li-ion batteries—Part II: Aged battery model and state of charge estimation | |
Mohamed et al. | Electrical circuit model of a vanadium redox flow battery using extended Kalman filter | |
US9581988B2 (en) | Method and system for battery state of charge estimation | |
Moura | Estimation and control of battery electrochemistry models: A tutorial | |
Park et al. | Hybrid electrochemical modeling with recurrent neural networks for li-ion batteries | |
Li et al. | A novel state estimation approach based on adaptive unscented Kalman filter for electric vehicles | |
CN105891724A (zh) | 基于扩展单粒子模型的锂离子电池荷电状态在线估计方法 | |
Xu et al. | State of charge estimation for liquid metal battery based on an improved sliding mode observer | |
Shi et al. | State of charge estimation by joint approach with model-based and data-driven algorithm for lithium-ion battery | |
CN111812530A (zh) | 一种锂电池荷电状态估测方法 | |
Blondel et al. | Nonlinear circle-criterion observer design for an electrochemical battery model | |
Mu et al. | Modeling, evaluation, and state estimation for batteries | |
Lin et al. | Lithium-ion battery state of charge/state of health estimation using SMO for EVs |
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 |