CN110208704A - 一种基于电压滞后效应的锂电池建模方法和系统 - Google Patents
一种基于电压滞后效应的锂电池建模方法和系统 Download PDFInfo
- Publication number
- CN110208704A CN110208704A CN201910355957.0A CN201910355957A CN110208704A CN 110208704 A CN110208704 A CN 110208704A CN 201910355957 A CN201910355957 A CN 201910355957A CN 110208704 A CN110208704 A CN 110208704A
- Authority
- CN
- China
- Prior art keywords
- lithium battery
- voltage
- model
- constant
- integrated model
- 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
Links
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
-
- 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
- G01R31/388—Determining ampere-hour charge capacity or SoC involving voltage measurements
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Secondary Cells (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
Abstract
本发明提供了一种基于电压滞后效应的锂电池建模方法和系统,该方法包括:步骤A:对锂电池进行某温度下的恒流恒压充电和不同温度、倍率下的脉冲充放电测试实验,获取锂电池充放电端电压与锂电池荷电状态之间的函数关系,进而分析确定开路电压中的滞后电压响应;步骤B:建立锂电池综合模型,所述锂电池综合模型融合了具有由锂电池荷电状态获取的锂电池开路电压的二阶非线性的等效电路模型以及描述动态滞后电压的滞后模型;步骤C:分别对锂电池综合模型中的动态参数进行辨识;步骤D:将所获得的动态参数进行仿真验证。该建模方法为锂电池模型的精度与实时性提供保障,能够为基于模型的SOC估计方法提供更精确的信息,进一步提高了BMS的性能。
Description
技术领域
本发明涉及电动汽车动力电池管理技术领域,具体涉及一种基于电压滞后效应的锂电池建模方法和系统。
背景技术
近百年来随着道路上汽车数量的暴增,给自然环境、工作生活带来了越来越多的危害:能源危机、雾霾频发、交通堵塞等,汽车反而成为了社会进一步发展的障碍。在能源和环境的双重危机下,结合国情与汽车产业发展现状,我国不断推广电动汽车(包括纯电动汽车EV、插电式混合动力汽车PHEV)和燃料汽车FCEV等环保类型的汽车,致力于实现零排放、零污染和能源转型。
锂离子电池由于其能量密度高、成本低、循环寿命长已经广泛应用在电动汽车上作为能量存储单元,为了保证电池组的安全可靠,需要动力电池管理系统(BMS)监控电池运行过程中的实时状态以控制充放电过程。同时,一个安全可靠的BMS还可以通过电池均衡,提高续航里程,延长电池使用寿命,否则,将会对电池造成不可逆的损害。通常,BMS采用模型来描述动力电池在使用过程中的变化特征,通过所建立的数学模型可以表示电池内部的化学反应原理和外部电压电流之间的关系,获取电池荷电状态信息(SOC),得到剩余续航里程,避免过充过放,因此精确的锂离子电池模型可以为基于模型的SOC估计方法提供更精确的信息,对于BMS的性能至关重要。SOC的不准确会降低用户满意度,是企业亟待解决的技术难题。
由于电化学电池非线性特征明显,目前广泛使用的有3种电池模型:电化学模型、等效电路模型(ECM)和黑箱模型(数据驱动模型和神经网络模型)。其中ECM由于未知参数少、计算简单应用较为广泛。模型中包含有:被控电压源表示电池开路电压(OCV),欧姆电阻表示直流内阻和RC并联电路表示电池内部极化反应。其中OCV-SOC关系对SOC估计过程中初值的获取至关重要,然而,由于电池OCV-SOC曲线平坦,存在滞后电压响应,直接影响电池模型的准确性与实时性,即使是较小的OCV误差也会增大SOC估计误差。
由于OCV-SOC关系通常是根据先验的充放电测试确定的,因此必须保证实验的精确同时进行参数校准,基于容量增量的分析方法已被用来测量OCV和OCV滞后,但是当初始SOC值不正确,会抵消充电和放电的OCV-SOC曲线,得到的滞后电压部分通常存在较大误差。因此,准确估计OCV之后电压对SOC估计至关重要,而现有的OCV测试都没有考虑时间对于滞后参数间耦合的作用,因此本发明除了建立了电池的综合模型外,也提出了新的滞后测试方法,用来获取精确的滞后电压响应。
发明内容
针对现有的锂离子电池模型准确性低和实时性差导致增大SOC估计误差从而影响影响锂电池使用寿命以及存在安全隐患等问题,本发明提供了一种基于电压滞后效应的锂电池建模方法,通过获取锂电池充放电端电压与SOC的函数关系来确定开路电压的滞后电压响应,为锂电池模型的精度与实时性提供保障,能够为基于模型的SOC估计方法提供更精确的信息,进一步提高了BMS的性能。本发明还涉及一种基于电压滞后效应的锂电池建模系统。
本发明技术方案如下:
一种基于电压滞后效应的锂电池建模方法,其特征在于,所述方法包括下述步骤:
步骤A:对锂电池进行某温度下的恒流恒压充电以获取锂电池容量信息,再进行不同温度、倍率下的脉冲充放电测试实验获取锂电池充放电端电压与锂电池荷电状态之间的函数关系,进而分析确定开路电压中的滞后电压响应;
步骤B:建立锂电池综合模型,所述锂电池综合模型融合了具有由锂电池荷电状态获取的锂电池开路电压的二阶非线性的等效电路模型以及描述动态滞后电压的滞后模型;
步骤C:分别对锂电池综合模型中的动态参数进行辨识;
步骤D:将所获得的动态参数进行仿真验证。
进一步地,所述步骤A中的恒流恒压充电和脉冲充放电测试实验包括下述步骤:
步骤A1:在某温度下静置3-5小时后,首先采用恒流—恒压的方式将锂电池充满电,继而恒流放电,然后静置3-5小时,重复3次确定锂电池的实际容量;
步骤A2:静置3-5小时;
步骤A3:然后在不同温度、倍率下采用与充电电流相同的电流值进行脉冲放电,每放锂电池容量的特定百分比后静置2小时,至截止电压;
步骤A4:最后再采用与步骤A1相同的充电电流进行脉冲充电,每充锂电池容量的特定百分比后静置2小时,至截止电流。
进一步地,所述步骤B中是根据Play算子的包络函数的输入输出关系将具有不同阈值和加权值的play算子进行线性加权叠加进而获取的滞后模型。
进一步地,所述步骤B中的等效电路模型具有开路电压、滞后电压相应、电化学极化电压以及浓差极化电压。
进一步地,所述步骤C中的分别对锂电池综合模型中的动态参数进行辨识包括:对建立的锂电池综合模型进行仿真,利用参数估计算法辨识锂电池综合模型的动态参数,对锂电池综合模型的不同参数利用最小二乘法进行自动优化,并通过信赖域反射算法约束优化,直至仿真输出与步骤A的实验结果匹配。
进一步地,所述步骤D中的将所获得的动态参数进行仿真验证包括:对锂电池进行恒流放电,记录锂电池的端电压数据,然后以同样的电流作为锂电池综合模型的输入,设置仿真时间与实验时间一致,将仿真步长设定为固定步长,并与实验数据采样时间间隔一致,获得锂电池综合模型的电压估计值结果。
一种基于电压滞后效应的锂电池建模系统,其特征在于,包括依次连接的实验模块、建模模块、参数辨识模块和参数验证模块,
所述实验模块:对锂电池进行某温度下的恒流恒压充电和不同温度、倍率下的脉冲充放电测试实验,获取锂电池充放电端电压与锂电池荷电状态之间的函数关系,进而分析确定开路电压中的滞后电压响应;
所述建模模块:建立锂电池综合模型,所述锂电池综合模型融合了具有由锂电池荷电状态获取的锂电池开路电压的二阶非线性的等效电路模型以及描述动态滞后电压的滞后模型;
所述参数辨识模块:分别对锂电池综合模型中的动态参数进行辨识;
所述参数验证模块:将所获得的动态参数进行仿真验证。
进一步地,所述建模模块是根据Play算子的包络函数的输入输出关系将具有不同阈值和加权值的play算子进行线性加权叠加进而获取的滞后模型,所述等效电路模型具有开路电压、滞后电压相应、电化学极化电压以及浓差极化电压。
进一步地,所述参数辨识模块是对建立的锂电池综合模型进行仿真,利用参数估计算法辨识锂电池综合模型的动态参数,对锂电池综合模型的不同参数利用最小二乘法进行自动优化,并通过信赖域反射算法约束优化,直至仿真输出与实验模块的实验结果匹配。
进一步地,所述参数验证模块是对锂电池进行恒流放电,记录锂电池的端电压数据,然后以同样的电流作为锂电池综合模型的输入,设置仿真时间与实验时间一致,将仿真步长设定为固定步长,并与实验数据采样时间间隔一致,获得锂电池综合模型的电压估计值结果。
本发明技术效果如下:
本发明提供了一种基于电压滞后效应的锂电池建模方法,该方法考虑了动态滞后电压特征,首先设计合理有效的锂电池充放电测试实验,获得电池在不同环境温度、不同放电倍率下的电池电压、电流等数据,获取锂电池充放电端电压与锂电池荷电状态的函数关系,进而分析确定开路电压中的滞后电压响应,避免了实验方法对于数据以及模型精度的影响,考虑了滞后时间增加时参数之间的耦合作用,可以获得更精确的电池OCV及其滞后响应;再建立等效电路模型和滞后模型并进一步融合为锂电池综合模型,等效电路模型为具有由锂电池荷电状态获取的锂电池开路电压的二阶非线性模型,即为简化的二阶GNL模型,滞后模型即为PI模型其用来描述动态滞后电压,与现有的ECM电池建模方法相比,本发明能够减小由于锂电池开路电压的不对称滞后所带来的模型误差,相比电化学模型模型参数少,计算简单、精度高,从理论上能够有效提高基于模型的SOC估计精度,并保证初始电压误差和电流测量值误差下的鲁棒性。该方法设计了电池充放电测试实验,获得电池在不同环境温度、不同放电倍率下的电池电压、电流等数据,然后针对锂电池电压滞后特征和模型参数固定、适用范围有限的问题,利用模型化的思想对电池的开路电压、电流、内阻等参数进行非线性建模,获得基于数据的电池电压滞后特征。该方法主要包括分析不同应力下的电池特性,包括电池不同参数之间的变化规律,建立动力锂电池电压滞后综合模型。本发明建立的模型能够准确描述锂电池的充电与放电特性,模型拟合精度高,解决了现有的锂离子电池模型由于准确性低和实时性差导致增大SOC估计误差从而影响影响锂电池使用寿命以及存在安全隐患等问题,为锂电池模型的精度与实时性提供保障,能够为基于模型的SOC估计方法提供更精确的信息,有助于设计动力电池管理系统和开发高级功能算法,进一步提高了BMS的性能。
本发明还涉及一种基于电压滞后效应的锂电池建模系统,该系统与上述的基于电压滞后效应的锂电池建模方法相对应,可理解为是实现上述方法的系统,该系统包括依次连接的实验模块、建模模块、参数辨识模块和参数验证模块,各模块协同工作,实现电池充放电测试实验、滞后模型及等效电路模型进而融合锂电池综合模型、参数辨识和参数修正。通过获取锂电池充放电端电压与SOC的函数关系来确定开路电压的滞后电压响应,为锂电池模型的精度与实时性提供保障,能够为基于模型的SOC估计方法提供更精确的信息。
附图说明
图1为本发明基于电压滞后效应的锂电池建模方法的流程图。
图2为充放电电压与SOC的函数关系图。
图3为锂电池的等效电路模型。
图4为PI模型中广义加权play算子的输入输出关系。
图5为锂电池综合模型中的动态参数进行辨识方法流程图。
图6为锂电池综合精度验证端电压对比曲线。
图中各标号列示如下:
1—充电电压测量值;2—放电电压测量值;OCVcha—充电工况下的SOC-OCV曲线;OCVdis—放电工况下的SOC-OCV曲线;OCVave—充放电的平均电压。
具体实施方式
下面结合附图对本发明的技术方案做详细的说明。
本发明涉及一种基于电压滞后效应的锂电池建模方法,也就是说,本发明的目的在于,提供一种锂电池(比如磷酸铁锂电池)滞后电压模型的建模方法,用于确定开路电压的滞后电压响应,为锂电池模型的精度与实时性提供保障。该方法的流程图如图1所示,包括下述步骤:
步骤A:也称为脉冲充放电测试实验步骤,是对锂电池进行不同温度、倍率下的恒流恒压充电和脉冲充放电测试实验,也就是说,是对锂电池进行某温度下的恒流恒压充电以获取锂电池容量信息,再进行不同温度、倍率下的脉冲充放电测试实验获取锂电池充放电端电压与锂电池荷电状态之间的函数关系,进而分析确定开路电压中的滞后电压响应;如图1所示,在恒流I、倍率C-rate下对动力电池(即锂电池)进行脉冲充放电测试,获取到的锂电池充放电端电压与锂电池荷电状态的函数关系可以理解为是放入至电池充放电测试系统中,再将电荷状态SOC、电流I、电压U、温度T等输入建模步骤的模型;
步骤B:也称为建模步骤,是建立锂电池综合模型,所述锂电池综合模型融合了具有由锂电池荷电状态获取的锂电池开路电压的二阶非线性的等效电路模型以及描述动态滞后电压的滞后模型;滞后模型也可称为是滞后PI(Prandtl–Ishlinskii)模型;
步骤C:也称为参数辨识步骤,分别对锂电池综合模型中的动态参数进行辨识;
步骤D:也称为参数验证,将所获得的动态参数进行仿真验证。
其中,步骤A所述的恒流恒压充电和脉冲充放电测试实验方面,现有的电池测试标准不包括OCV测量的测试方法和OCV滞后的识别。目前获取滞后使用的测试方法主要有:低电流充电/放电,增量电流充电/放电,其滞后的评估方法大多存在错误,根本原因包括:静置时间短、电池容量随测试方法不同、充放电端电压发生偏移等等。实际上,电池的开路电压与放电容量存在一一对应的关系,在滞后测试方法中,增加步数(或减小步长)后由于极化反应减弱,电池总的放电容量就会增加。本发明经过实验验证后建立了一个准确的电池滞后特征测试方法,完整的测试过程包括:
步骤A1:在某温度下静置3-5小时后,首先采用恒流—恒压的方式将锂电池充满电,继而恒流放电,然后静置3-5小时,重复3次确定锂电池的实际容量;具体地,比如在T1温度下静置4h后,首先采用恒流-恒压的方式将电池充满电,充电电流为1C,截止电流为1/20C;
步骤A2:静置3-5小时;优选静置4h;
步骤A3:然后在不同温度、倍率下采用与充电电流相同的电流值进行脉冲放电,每放锂电池容量的特定百分比后静置2小时,至截止电压;具体地,在不同温度、倍率下,可采用1C电流放电,每放电池容量的5%静置2h,至截止电压;
步骤A4:最后再采用与步骤A1相同的充电电流进行脉冲充电,每充锂电池容量的特定百分比后静置2小时,至截止电流。具体地,在不同温度、倍率下,可采用1C电流充电,每充电池容量的5%静置2h,至截止电流。
实验方法对于OCV特征的获取和滞后的判断极其重要,本发明中的测试方法同时考虑了实验时间和分辨率,能够作为研究应用的重要参考。
如图2所示,在对锂电池进行完整的充放电测试之后,锂电池充放电端电压与SOC的函数关系。OCVcha为充电工况下的SOC-OCV曲线,OCVdis为放电工况下的SOC-OCV曲线,OCVave为充放电的平均电压。实验发现:在任何SOC值处,充电曲线的电压值都要高于放电曲线。即充电电压测量值1高于放电电压测量值2。这是由于电池内阻的存在,在充电过程中充电电流经过内阻导致电池的端电压大于电池的电动势,相反,在放电过程中内阻上的电压又使得端电压的值小于电池的电动势,因此,取它们的平均值OCVave作为等效电路模型中关于开路电压的SOC-OCV函数曲线不包含电池内阻上的电压影响,从而获得精确的SOC-OCV函数关系。
步骤B的建模步骤,是建立锂电池综合模型,其包括等效电路模型和滞后模型。其中,等效电路模型(Equivalent Circuit Model,ECM)如图3所示,等效电路模型是采用电容、电阻、电压源等理想元件组成电路网络来描述电池的伏安特性,各参数的物理含义明确,可以很好的模拟电池动态行为的基本特点。这种模型属于半机理的灰箱模型,易于计算、实时性较高,便于分析、通用性强。
这是一个动态的线性时变(LTV)系统,用来模拟电池的瞬态行为,既能较好的描述电池动态特性,观测电池参数变化,又可以在线准确估测SOC。在图3所示的等效电路模型中,左图里的Rsd是锂电池静置的自放电电阻,电容CQ表示电池内部容量;在前面的电池性能测试分析中已经获得了电池开路滞后电压VOCV,主要表现为充放电过程中的SOC-OCV曲线不一致,因而在该等效电路模型右图表示的充放电状态中采用Ucd表示电池的滞后电压,两个串联的RC并联电路分别表示电池的电化学极化内阻Re和电化学极化电容Ce,浓差极化电阻Rd和浓差极化电容Cd,电化学极化是活性极化现象,浓差极化是由于电池内部的离子扩散;R0表示电池的欧姆内阻,主要来源为电池电解液的内阻;Ui为电池端电压。
等效电路模型由二阶GNL模型简化而来,首选获得锂电池的开路电压为:
OCVk=aSOCk 8+bSOCk 7+cSOCk 6+dSOCk 5+eSOCk 4+fSOCk 3+gSOCk 2+hSOCk+mlogSOCk
其中,OCVk为k时刻电池的开路电压,SOCk为k时刻电池的荷电状态,a、b、c、d、e、f、g、h、m为系数。
根据该二阶等效电路模型列写KVL方程如下:
UOCV=Uo+Ucd+Ui(t)+Ue(t)+Ud(t) (1)
Uo=RoI(t) (2)
经拉氏变化后的电池等效电路方程为:
其中,Ui为电池端电压,Ucd为滞后电压响应,Re为电化学极化内阻,为电化学极化电容,Rd为浓差极化电阻,Cd浓差极化电容,R0表示电池的欧姆内阻。滞后电压Ucd可以用一个非线性滞后函数描述,滞后函数是一个有关于SOC变化率和当前电压的函数,与放电/充电状态有关。
本发明的滞后模型是根据Play算子的包络函数的输入输出关系将具有不同阈值和加权值的play算子进行线性加权叠加等过程获取的。PI模型相对于其他滞后模型的主要优点是表达式简单、可以通过解析的方法获得其逆,并降低了建模复杂度,从而提高了实时应用中的效率。PI模型将具有不同阈值和加权值的初等线性play算子进行线性加权叠加,其play算子是一种具有对称性和与速率无关特性的基本滞后算子。
然而,由于线性play算子是对称并且凸的,因此在描述非线性饱和的滞后回路时,通常有三种方式:1)用非线性算子代替线性算子;2)使用非对称无内存叠加运算符级联Play算子;3)添加非线性无记忆函数/非线性输入函数。本发明采用第二种方式利用不同的非线性包络函数逼近非对称滞后循环,可以表征不对称且饱和的滞后回线。
如图4所示,是PI模型中广义加权play算子的输入输出关系,输入x的范围为[ζr,ζl],随着x的增大输出ω沿着曲线γr增大,反之,则沿着曲线γl减小。两个曲线γr和γl就是play算子的包络函数,连续递增且γl≤γr。当阈值γ≥0,Play算子的输出为:
其中gmr(x,ωm)=max{m0(x-r),min[m0(x+r),ωm]},对于任何分段均匀输入x(t)∈C[0,T],它在每个子空间[ti,ti+1]中都是单调的i=0,1,…,N-1;m0∈R+为引入的斜坡向量来调整滞后形状;r是play算子的阈值:
采用广义play算子的加权叠加来描述非线性滞后电压的线性输入函数。本发明利用不同的非线性包络函数逼近非对称滞后循环,可以表征不对称和饱和滞后回线;将有限个play算子的加权叠加后得到PI模型:
其中p(r)是可积密度函数,H为非递减Lipschitz连续函数,满足:
q是常数且为正,实际应用中,将式(8)中的模型离散化得到有限个play算子的离散化表达式:
其中n为play算子的个数,第i个算子的阈值ri=αi,那么常数α、ρ和τ由充放电实验数据获得。该模型使用不同的非线性包络函数以获得非对称滞后函数。
步骤C中的分别对锂电池综合模型中的动态参数进行辨识,优选为:首先在Matlab/Simscape中对建立的锂电池综合模型进行仿真,利用Parameter Estimation参数估计算法估计锂电池综合模型的动态参数,对锂电池综合模型的不同参数利用最小二乘法进行自动优化,并通过Trust-Region-Reflective(信赖域反射算法)约束优化,直至仿真输出与步骤A的实验结果完全匹配,最后验证模型精度。在参数辨识过程中,将实验中的电流、电压数据作为该仿真模型的输入,仿真后获得模型的仿真电压输出,参数的辨识过程就是将模型仿真输出电压匹配到测量电压的过程。电池的等效电路模型是为了更好地模拟电池在充放电过程中的动态变化,因此ECM的输出必须能够动态匹配电池实际输出,这样该电路模型才能够完整地替代实际电池进行分析控制。其中的参数优化方法为非线性最小二乘法,采用的算法为Trust-Region-Reflective(信赖域反射算法),属于约束优化中的求解算法。如图5所示的参数辨识流程为:首先输入脉冲电流,为仿真模型参数设置初值与有效区间,然后将仿真输出结果与采集到的锂电池输出电压进行对比验证直至误差小于限定阈值则认为完全匹配。Em(SOC)即为VOCV(SOC),R/C即为Re、Rd或Ce、Cd等参数,这种参数估计方法将参数估计作为一个优化问题,能够同时估计和验证模型中的多个参数,并制定参数的边界。
步骤D中的将所获得的动态参数进行仿真验证具体为:对锂电池进行恒流放电,记录锂电池的端电压数据,然后以同样的电流作为锂电池综合模型的输入,设置仿真时间与实验时间一致,将仿真步长设定为固定步长,并与实验数据采样时间间隔一致,获得锂电池综合模型的电压估计值结果。如图6所示,在使用了脉冲放电进行了参数估计之后,如何评价参数估计结果的优劣,就需要对模型参数进行实验验证。本发明采用恒流放电实验进行电池模型参数的验证,具体的做法是先对电池进行1/3C(0.56A)恒流放电,记录电池的端电压数据,然后以同样的电流作为模型的输入,设置仿真时间与实验时间一致,将仿真步长设定为固定步长,并与实验数据采样时间间隔一致,获得模型的电压估计值结果。图6所示为电池端电压的测量值与仿真输出对比曲线,在整个放电过程中,锂电池综合模型输出能够紧紧跟踪电池实际的端电压,验证了模型参数的准确性,也说明了动态参数能够确保模型的精确度。
本发明还涉及一种基于电压滞后效应的锂电池建模系统,该系统与上述的基于电压滞后效应的锂电池建模方法相对应,可理解为是实现上述方法的系统,该系统包括依次连接的实验模块、建模模块、参数辨识模块和参数验证模块。其中,实验模块:对锂电池进行某温度下的恒流恒压充电和不同温度、倍率下的脉冲充放电测试实验,获取锂电池充放电端电压与锂电池荷电状态之间的函数关系,进而分析确定开路电压中的滞后电压响应;建模模块:建立锂电池综合模型,所述锂电池综合模型融合了具有由锂电池荷电状态获取的锂电池开路电压的二阶非线性的等效电路模型以及描述动态滞后电压的滞后模型;参数辨识模块:分别对锂电池综合模型中的动态参数进行辨识;参数验证模块:将所获得的动态参数进行仿真验证。
优选地,建模模块是根据Play算子的包络函数的输入输出关系将具有不同阈值和加权值的play算子进行线性加权叠加进而获取的滞后模型,可参考上述公式(6)—(10);等效电路模型可参考图3所示,具有开路电压、滞后电压相应、电化学极化电压以及浓差极化电压,模型计算可参考上述公式(1)—(5)。
参数辨识模块是对建立的锂电池综合模型进行仿真,利用参数估计算法辨识锂电池综合模型的动态参数,对锂电池综合模型的不同参数利用最小二乘法进行自动优化,并通过信赖域反射算法约束优化,直至仿真输出与实验模块的实验结果匹配。参数验证模块是对锂电池进行恒流放电,记录锂电池的端电压数据,然后以同样的电流作为锂电池综合模型的输入,设置仿真时间与实验时间一致,将仿真步长设定为固定步长,并与实验数据采样时间间隔一致,获得锂电池综合模型的电压估计值结果。
本发明的锂电池建模系统通过各模块协同工作,实现电池充放电测试实验、动力电池、滞后电压模型、等效电路模型、参数辨识和参数修正。实验模块实现电池充放电测试实验,获得电池在不同环境温度、不同放电倍率下的电池电压、电流等数据,然后针对锂电池电压滞后特征和模型参数固定、适用范围有限的问题,由建模模块利用模型化的思想对电池的开路电压、电流、内阻等参数进行非线性建模,获得基于数据的电池电压滞后特征。该方法主要包括分析不同应力下的电池特性,包括电池不同参数之间的变化规律,建立动力锂电池电压滞后综合模型。再由参数辨识模块在Matlab/Simulink中搭建电池Simscape模型并使用Parameter Estimation完成电池模型中参数的动态辨识。本发明建立的模型能够准确描述电池的充电与放电特性,模型拟合精度高,有助于设计动力电池管理系统和开发高级功能算法,提高了BMS的性能。
应当指出,以上所述具体实施方式可以使本领域的技术人员更全面地理解本发明创造,但不以任何方式限制本发明创造。因此,尽管本说明书参照附图和实施例对本发明创造已进行了详细的说明,但是,本领域技术人员应当理解,仍然可以对本发明创造进行修改或者等同替换,总之,一切不脱离本发明创造的精神和范围的技术方案及其改进,其均应涵盖在本发明创造专利的保护范围当中。
Claims (10)
1.一种基于电压滞后效应的锂电池建模方法,其特征在于,所述方法包括下述步骤:
步骤A:对锂电池进行某温度下的恒流恒压充电以获取锂电池容量信息,再进行不同温度、倍率下的脉冲充放电测试实验获取锂电池充放电端电压与锂电池荷电状态之间的函数关系,进而分析确定开路电压中的滞后电压响应;
步骤B:建立锂电池综合模型,所述锂电池综合模型融合了具有由锂电池荷电状态获取的锂电池开路电压的二阶非线性的等效电路模型以及描述动态滞后电压的滞后模型;
步骤C:分别对锂电池综合模型中的动态参数进行辨识;
步骤D:将所获得的动态参数进行仿真验证。
2.根据权利要求1所述的锂电池建模方法,其特征在于,所述步骤A中的恒流恒压充电和脉冲充放电测试实验包括下述步骤:
步骤A1:在某温度下静置3-5小时后,首先采用恒流—恒压的方式将锂电池充满电,继而恒流放电,然后静置3-5小时,重复3次确定锂电池的实际容量;
步骤A2:静置3-5小时;
步骤A3:然后在不同温度、倍率下采用与充电电流相同的电流值进行脉冲放电,每放锂电池容量的特定百分比后静置2小时,至截止电压;
步骤A4:最后再采用与步骤A1相同的充电电流进行脉冲充电,每充锂电池容量的特定百分比后静置2小时,至截止电流。
3.根据权利要求1或2所述的锂电池建模方法,其特征在于,所述步骤B中是根据Play算子的包络函数的输入输出关系将具有不同阈值和加权值的play算子进行线性加权叠加进而获取的滞后模型。
4.根据权利要求3所述的锂电池建模方法,其特征在于,所述步骤B中的等效电路模型具有开路电压、滞后电压相应、电化学极化电压以及浓差极化电压。
5.根据权利要求3所述的锂电池建模方法,其特征在于,所述步骤C中的分别对锂电池综合模型中的动态参数进行辨识包括:对建立的锂电池综合模型进行仿真,利用参数估计算法辨识锂电池综合模型的动态参数,对锂电池综合模型的不同参数利用最小二乘法进行自动优化,并通过信赖域反射算法约束优化,直至仿真输出与步骤A的实验结果匹配。
6.根据权利要求5所述的锂电池建模方法,其特征在于,所述步骤D中的将所获得的动态参数进行仿真验证包括:对锂电池进行恒流放电,记录锂电池的端电压数据,然后以同样的电流作为锂电池综合模型的输入,设置仿真时间与实验时间一致,将仿真步长设定为固定步长,并与实验数据采样时间间隔一致,获得锂电池综合模型的电压估计值结果。
7.一种基于电压滞后效应的锂电池建模系统,其特征在于,包括依次连接的实验模块、建模模块、参数辨识模块和参数验证模块,
所述实验模块:对锂电池进行某温度下的恒流恒压充电和不同温度、倍率下的脉冲充放电测试实验,获取锂电池充放电端电压与锂电池荷电状态之间的函数关系,进而分析确定开路电压中的滞后电压响应;
所述建模模块:建立锂电池综合模型,所述锂电池综合模型融合了具有由锂电池荷电状态获取的锂电池开路电压的二阶非线性的等效电路模型以及描述动态滞后电压的滞后模型;
所述参数辨识模块:分别对锂电池综合模型中的动态参数进行辨识;
所述参数验证模块:将所获得的动态参数进行仿真验证。
8.根据权利要求7所述的锂电池建模系统,其特征在于,所述建模模块是根据Play算子的包络函数的输入输出关系将具有不同阈值和加权值的play算子进行线性加权叠加进而获取的滞后模型,所述等效电路模型具有开路电压、滞后电压相应、电化学极化电压以及浓差极化电压。
9.根据权利要求7或8所述的锂电池建模系统,其特征在于,所述参数辨识模块是对建立的锂电池综合模型进行仿真,利用参数估计算法辨识锂电池综合模型的动态参数,对锂电池综合模型的不同参数利用最小二乘法进行自动优化,并通过信赖域反射算法约束优化,直至仿真输出与实验模块的实验结果匹配。
10.根据权利要求9所述的锂电池建模系统,其特征在于,所述参数验证模块是对锂电池进行恒流放电,记录锂电池的端电压数据,然后以同样的电流作为锂电池综合模型的输入,设置仿真时间与实验时间一致,将仿真步长设定为固定步长,并与实验数据采样时间间隔一致,获得锂电池综合模型的电压估计值结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910355957.0A CN110208704B (zh) | 2019-04-29 | 2019-04-29 | 一种基于电压滞后效应的锂电池建模方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910355957.0A CN110208704B (zh) | 2019-04-29 | 2019-04-29 | 一种基于电压滞后效应的锂电池建模方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110208704A true CN110208704A (zh) | 2019-09-06 |
CN110208704B CN110208704B (zh) | 2021-08-06 |
Family
ID=67786712
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910355957.0A Active CN110208704B (zh) | 2019-04-29 | 2019-04-29 | 一种基于电压滞后效应的锂电池建模方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110208704B (zh) |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110795851A (zh) * | 2019-10-30 | 2020-02-14 | 西安理工大学 | 一种考虑环境温度影响的锂离子电池建模方法 |
CN110826023A (zh) * | 2019-11-13 | 2020-02-21 | 国网电力科学研究院有限公司 | 一种电池系统运行数据分辨率提升方法 |
CN111104759A (zh) * | 2019-12-09 | 2020-05-05 | 国联汽车动力电池研究院有限责任公司 | 锂电池焊接形状的设计仿真方法 |
CN111367262A (zh) * | 2020-04-17 | 2020-07-03 | 广东戈兰玛汽车系统有限公司 | 新能源汽车台架动力系统的bms仿真方法及系统 |
CN111506977A (zh) * | 2019-12-11 | 2020-08-07 | 安徽贵博新能科技有限公司 | 一种动力电池建模方法 |
CN111537887A (zh) * | 2020-04-27 | 2020-08-14 | 南京航空航天大学 | 考虑迟滞特性的混合动力系统电池开路电压模型优化方法 |
CN111914503A (zh) * | 2020-08-04 | 2020-11-10 | 重庆大学 | 一种锂离子电池功率输入电热耦合模型建立方法 |
CN112180278A (zh) * | 2020-09-28 | 2021-01-05 | 天津大学 | 考虑电压迟滞特性的电动汽车动力电池性能无损检测方法 |
CN112379289A (zh) * | 2020-10-31 | 2021-02-19 | 浙江锋锂新能源科技有限公司 | 一种锂离子电池最大电流的测试方法 |
CN112379280A (zh) * | 2020-11-10 | 2021-02-19 | 南京理工大学 | 基于恒压恒流充电曲线的电池模型参数与ocv-soc关系确定方法 |
CN112462264A (zh) * | 2019-12-19 | 2021-03-09 | 淮南师范学院 | 一种基于asr-ukf的矿用锂电池soc估算方法 |
CN112731160A (zh) * | 2020-12-25 | 2021-04-30 | 东莞新能安科技有限公司 | 电池滞回模型训练方法、估算电池soc的方法和装置 |
CN112881921A (zh) * | 2021-01-08 | 2021-06-01 | 恒大新能源汽车投资控股集团有限公司 | 电池等效电路模型参数辨识方法、装置、设备及存储介质 |
CN113011007A (zh) * | 2021-02-26 | 2021-06-22 | 山东大学 | 一种锂离子动力电池热模型参数快速辨识方法及系统 |
CN113111579A (zh) * | 2021-04-02 | 2021-07-13 | 华北电力大学(保定) | 一种自适应天牛须优化神经网络的锂电池等效电路模型参数辨识方法 |
CN113125971A (zh) * | 2021-03-18 | 2021-07-16 | 深圳逸驰云动科技有限公司 | 一种确定电池开路电压的方法、装置及充电桩 |
CN113341319A (zh) * | 2021-04-19 | 2021-09-03 | 北京交通大学 | 基于参数插值获取任意温度和倍率下放电曲线的方法 |
CN113687251A (zh) * | 2021-08-23 | 2021-11-23 | 重庆邮电大学 | 一种基于双模型的锂离子电池组电压异常故障诊断方法 |
CN113985286A (zh) * | 2021-10-14 | 2022-01-28 | 合肥国轩高科动力能源有限公司 | 一种锂离子电池不同温度下soc-ocv测试方法 |
CN114236394A (zh) * | 2021-12-15 | 2022-03-25 | 天津市捷威动力工业有限公司 | 一种电芯最大电流的测试方法和应用 |
CN115101138A (zh) * | 2022-06-29 | 2022-09-23 | 上海玫克生储能科技有限公司 | 基于参数辨识的锂电池设计优化方法及装置、存储介质 |
WO2022242653A1 (zh) * | 2021-05-17 | 2022-11-24 | 蜂巢能源科技股份有限公司 | 确定电芯充放电性能的方法、装置、存储介质及电子设备 |
CN115494406B (zh) * | 2022-10-28 | 2023-03-24 | 北京航空航天大学 | 电池参数辨识方法、装置及电子设备 |
CN116593895A (zh) * | 2023-06-16 | 2023-08-15 | 中国科学技术大学 | 一种基于应变的锂离子电池组电流检测方法及系统 |
CN116680938A (zh) * | 2023-08-04 | 2023-09-01 | 通达电磁能股份有限公司 | 一种锂电池电热耦合模型的建模方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1243933A1 (en) * | 2001-03-23 | 2002-09-25 | General Motors Corporation | State of charge prediction method and apparatus for a battery |
CN101371155A (zh) * | 2005-12-13 | 2009-02-18 | 科巴西斯有限责任公司 | 充电电压滞后估算器的电池状态 |
CN105842627A (zh) * | 2016-02-01 | 2016-08-10 | 北京理工大学 | 基于数据模型融合的动力电池容量和荷电状态的估计方法 |
CN108736083A (zh) * | 2017-04-25 | 2018-11-02 | 通用汽车环球科技运作有限责任公司 | 用于评估含有表现出电压滞后的材料的电池的方法和设备 |
-
2019
- 2019-04-29 CN CN201910355957.0A patent/CN110208704B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1243933A1 (en) * | 2001-03-23 | 2002-09-25 | General Motors Corporation | State of charge prediction method and apparatus for a battery |
CN101371155A (zh) * | 2005-12-13 | 2009-02-18 | 科巴西斯有限责任公司 | 充电电压滞后估算器的电池状态 |
CN105842627A (zh) * | 2016-02-01 | 2016-08-10 | 北京理工大学 | 基于数据模型融合的动力电池容量和荷电状态的估计方法 |
CN108736083A (zh) * | 2017-04-25 | 2018-11-02 | 通用汽车环球科技运作有限责任公司 | 用于评估含有表现出电压滞后的材料的电池的方法和设备 |
Non-Patent Citations (1)
Title |
---|
张华辉 等: "动力锂离子电池荷电态估计的建模与仿真", 《哈尔滨工程大学学报》 * |
Cited By (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110795851B (zh) * | 2019-10-30 | 2024-05-14 | 西安理工大学 | 一种考虑环境温度影响的锂离子电池建模方法 |
CN110795851A (zh) * | 2019-10-30 | 2020-02-14 | 西安理工大学 | 一种考虑环境温度影响的锂离子电池建模方法 |
CN110826023A (zh) * | 2019-11-13 | 2020-02-21 | 国网电力科学研究院有限公司 | 一种电池系统运行数据分辨率提升方法 |
CN110826023B (zh) * | 2019-11-13 | 2023-08-25 | 国网电力科学研究院有限公司 | 一种电池系统运行数据分辨率提升方法 |
CN111104759A (zh) * | 2019-12-09 | 2020-05-05 | 国联汽车动力电池研究院有限责任公司 | 锂电池焊接形状的设计仿真方法 |
CN111104759B (zh) * | 2019-12-09 | 2023-09-29 | 国联汽车动力电池研究院有限责任公司 | 锂电池焊接形状的设计仿真方法 |
CN111506977A (zh) * | 2019-12-11 | 2020-08-07 | 安徽贵博新能科技有限公司 | 一种动力电池建模方法 |
CN112462264A (zh) * | 2019-12-19 | 2021-03-09 | 淮南师范学院 | 一种基于asr-ukf的矿用锂电池soc估算方法 |
CN111367262A (zh) * | 2020-04-17 | 2020-07-03 | 广东戈兰玛汽车系统有限公司 | 新能源汽车台架动力系统的bms仿真方法及系统 |
CN111537887A (zh) * | 2020-04-27 | 2020-08-14 | 南京航空航天大学 | 考虑迟滞特性的混合动力系统电池开路电压模型优化方法 |
CN111537887B (zh) * | 2020-04-27 | 2021-10-01 | 南京航空航天大学 | 考虑迟滞特性的混合动力系统电池开路电压模型优化方法 |
CN111914503B (zh) * | 2020-08-04 | 2023-10-20 | 重庆大学 | 一种锂离子电池功率输入电热耦合模型建立方法 |
CN111914503A (zh) * | 2020-08-04 | 2020-11-10 | 重庆大学 | 一种锂离子电池功率输入电热耦合模型建立方法 |
CN112180278B (zh) * | 2020-09-28 | 2022-10-18 | 天津大学 | 考虑电压迟滞特性的电动汽车动力电池性能无损检测方法 |
CN112180278A (zh) * | 2020-09-28 | 2021-01-05 | 天津大学 | 考虑电压迟滞特性的电动汽车动力电池性能无损检测方法 |
CN112379289A (zh) * | 2020-10-31 | 2021-02-19 | 浙江锋锂新能源科技有限公司 | 一种锂离子电池最大电流的测试方法 |
CN112379289B (zh) * | 2020-10-31 | 2024-03-29 | 浙江锋锂新能源科技有限公司 | 一种锂离子电池最大电流的测试方法 |
CN112379280A (zh) * | 2020-11-10 | 2021-02-19 | 南京理工大学 | 基于恒压恒流充电曲线的电池模型参数与ocv-soc关系确定方法 |
CN112379280B (zh) * | 2020-11-10 | 2022-12-27 | 南京理工大学 | 基于恒压恒流充电曲线的电池模型参数与ocv-soc关系确定方法 |
CN112731160A (zh) * | 2020-12-25 | 2021-04-30 | 东莞新能安科技有限公司 | 电池滞回模型训练方法、估算电池soc的方法和装置 |
CN112881921A (zh) * | 2021-01-08 | 2021-06-01 | 恒大新能源汽车投资控股集团有限公司 | 电池等效电路模型参数辨识方法、装置、设备及存储介质 |
CN113011007A (zh) * | 2021-02-26 | 2021-06-22 | 山东大学 | 一种锂离子动力电池热模型参数快速辨识方法及系统 |
CN113125971A (zh) * | 2021-03-18 | 2021-07-16 | 深圳逸驰云动科技有限公司 | 一种确定电池开路电压的方法、装置及充电桩 |
CN113111579A (zh) * | 2021-04-02 | 2021-07-13 | 华北电力大学(保定) | 一种自适应天牛须优化神经网络的锂电池等效电路模型参数辨识方法 |
CN113111579B (zh) * | 2021-04-02 | 2022-09-09 | 华北电力大学(保定) | 一种自适应天牛须优化神经网络的锂电池等效电路模型参数辨识方法 |
CN113341319B (zh) * | 2021-04-19 | 2023-02-21 | 北京交通大学 | 基于参数插值获取任意温度和倍率下放电曲线的方法 |
CN113341319A (zh) * | 2021-04-19 | 2021-09-03 | 北京交通大学 | 基于参数插值获取任意温度和倍率下放电曲线的方法 |
WO2022242653A1 (zh) * | 2021-05-17 | 2022-11-24 | 蜂巢能源科技股份有限公司 | 确定电芯充放电性能的方法、装置、存储介质及电子设备 |
CN113687251A (zh) * | 2021-08-23 | 2021-11-23 | 重庆邮电大学 | 一种基于双模型的锂离子电池组电压异常故障诊断方法 |
CN113687251B (zh) * | 2021-08-23 | 2024-07-30 | 重庆邮电大学 | 一种基于双模型的锂离子电池组电压异常故障诊断方法 |
CN113985286B (zh) * | 2021-10-14 | 2024-03-08 | 合肥国轩高科动力能源有限公司 | 一种锂离子电池不同温度下soc-ocv测试方法 |
CN113985286A (zh) * | 2021-10-14 | 2022-01-28 | 合肥国轩高科动力能源有限公司 | 一种锂离子电池不同温度下soc-ocv测试方法 |
CN114236394A (zh) * | 2021-12-15 | 2022-03-25 | 天津市捷威动力工业有限公司 | 一种电芯最大电流的测试方法和应用 |
CN115101138A (zh) * | 2022-06-29 | 2022-09-23 | 上海玫克生储能科技有限公司 | 基于参数辨识的锂电池设计优化方法及装置、存储介质 |
CN115494406B (zh) * | 2022-10-28 | 2023-03-24 | 北京航空航天大学 | 电池参数辨识方法、装置及电子设备 |
CN116593895A (zh) * | 2023-06-16 | 2023-08-15 | 中国科学技术大学 | 一种基于应变的锂离子电池组电流检测方法及系统 |
CN116593895B (zh) * | 2023-06-16 | 2024-02-23 | 中国科学技术大学 | 一种基于应变的锂离子电池组电流检测方法及系统 |
CN116680938A (zh) * | 2023-08-04 | 2023-09-01 | 通达电磁能股份有限公司 | 一种锂电池电热耦合模型的建模方法及系统 |
CN116680938B (zh) * | 2023-08-04 | 2023-10-20 | 通达电磁能股份有限公司 | 一种锂电池电热耦合模型的建模方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN110208704B (zh) | 2021-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110208704A (zh) | 一种基于电压滞后效应的锂电池建模方法和系统 | |
Cao et al. | Multi-timescale parametric electrical battery model for use in dynamic electric vehicle simulations | |
CN103926538B (zh) | 基于aic准则的变阶数rc等效电路模型及实现方法 | |
CN106569136B (zh) | 一种电池健康状态在线估计方法及系统 | |
CN105807226B (zh) | 基于等效电路模型的锂离子电池soc预测方法及装置 | |
Zheng et al. | A novel capacity estimation method based on charging curve sections for lithium-ion batteries in electric vehicles | |
CN105206888B (zh) | 一种锂离子电池内部温度监测方法 | |
CN102590751B (zh) | 动力电池包一致性评估方法和装置 | |
CN103323775B (zh) | 用于电池模块的平衡监控及测试系统 | |
CN104392080B (zh) | 一种锂电池分数阶变阶等效电路模型及其辨识方法 | |
CN105548896A (zh) | 基于n-2rc模型的动力电池soc在线闭环估计方法 | |
CN102680795B (zh) | 一种二次电池内阻的实时在线估计方法 | |
Wang et al. | An improved coulomb counting method based on dual open‐circuit voltage and real‐time evaluation of battery dischargeable capacity considering temperature and battery aging | |
CN103744026A (zh) | 基于自适应无迹卡尔曼滤波的蓄电池荷电状态估算方法 | |
CN104965179A (zh) | 一种锂离子蓄电池的温度组合电路模型及其参数识别方法 | |
CN104714177A (zh) | 锂电池管理系统的检测方法及系统 | |
CN102565710A (zh) | 用于估计蓄电池健康状态的方法和装置 | |
Feng et al. | A graphical model for evaluating the status of series‐connected lithium‐ion battery pack | |
CN103197251A (zh) | 一种动力锂电池二阶rc等效模型的辨识方法 | |
CN104849672A (zh) | 基于等效电路模型的锂电池动态阻抗参数识别方法 | |
CN103744028A (zh) | 一种基于ukf的蓄电池荷电状态soc估算方法 | |
CN103632018A (zh) | 一种基于Simscape平台的电池建模方法 | |
CN106872899A (zh) | 一种基于降维观测器的动力电池soc估计方法 | |
Li et al. | Discharge capacity estimation for Li-ion batteries based on particle filter under multi-operating conditions | |
CN103713264A (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 |