CN112285569A - 一种基于动态阈值模型的电动汽车故障诊断方法 - Google Patents

一种基于动态阈值模型的电动汽车故障诊断方法 Download PDF

Info

Publication number
CN112285569A
CN112285569A CN202011181863.5A CN202011181863A CN112285569A CN 112285569 A CN112285569 A CN 112285569A CN 202011181863 A CN202011181863 A CN 202011181863A CN 112285569 A CN112285569 A CN 112285569A
Authority
CN
China
Prior art keywords
battery
soc
parameter
model
fault diagnosis
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
CN202011181863.5A
Other languages
English (en)
Other versions
CN112285569B (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 Institute of Spacecraft System Engineering
Harbin Institute of Technology Weihai
Original Assignee
Beijing Institute of Spacecraft System Engineering
Harbin Institute of Technology Weihai
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 Institute of Spacecraft System Engineering, Harbin Institute of Technology Weihai filed Critical Beijing Institute of Spacecraft System Engineering
Priority to CN202011181863.5A priority Critical patent/CN112285569B/zh
Publication of CN112285569A publication Critical patent/CN112285569A/zh
Application granted granted Critical
Publication of CN112285569B publication Critical patent/CN112285569B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/367Software therefor, e.g. for battery testing using modelling or look-up tables
    • 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)
  • Secondary Cells (AREA)

Abstract

本发明提供一种基于动态阈值模型的电动汽车故障诊断方法,该方法用于电动汽车中电池系统故障诊断,在阈值模型建立和参数辨识算法两方面进行了改进,在不同温度下进行电路基础特性测试实验,得到等效电路模型参数;建立OCV‑SOC‑Q三维响应面模型;采用带遗忘因子的递推最小二乘法进行模型参数辨识,建立关于R0和τ的动态阈值模型。在实际故障诊断过程当中,利用双扩展卡尔曼滤波算法辨识参数和状态,得到电池R0和τ、容量及SOC;采用温度插值的方法确定参数参考值;确定参数阈值;生成残差;通过对比残差与阈值来判断电池是否发生故障。该方法不仅故障诊断率高,还能避免检测不及时、误警和漏警问题。

Description

一种基于动态阈值模型的电动汽车故障诊断方法
技术领域
本发明主要涉及新能源汽车系统相关技术领域,具体是一种基于动态阈值模型的电动汽车故障诊断方法。
背景技术
当前电动汽车技术发展十分迅速,正处于市场推广的上升阶段。电动车作为一种新兴车型,相比传统车型,车辆安全及零部件质量问题比较突出。而电池作为电动汽车的核心零部件,其故障状态及运行寿命直接影响汽车使用功能及安全,同时也是生产厂商及4S店十分关心的内容。
动力电池系统故障诊断是电池管理系统的核心功能之一,常用的电池系统故障诊断方法是基于模型的方法,该方法共分残差生成和残差评价两步。残差生成是通过构建故障模型,得到参数或状态的估计值并与相应参数或状态的参考值做差生成残差,再通过对比残差与故障阈值来判断故障是否发生。鉴于电池模型参数是其状态估计的基础,通过残差生成残差是故障诊断的典型方法,但目前研究仍存在如下问题:
1)、在残差生成过程中,尽管参数的估计值易确定,但对应的参考值通常是固定值,实际上,电池在一定温度下的残差估计值对应的参考值也应是对应温度下的值,忽略温度影响易造成生成的残差过大或过小,均影响故障诊断精度;
2)、残差评价过程中,故障阈值通常也是固定值,而电池残差受温度、老化和电池荷电状态(State of charge,SOC)影响,若忽略SOC、温度和老化对残差的影响而采用固定阈值进行残差评价,易造成故障误警率和漏警率过高问题。因此,如何获得准确的参数参考值以生成残差且随着电池工作状态而更新阈值进行故障检测是当然的难点。
发明内容
为解决目前技术的不足,本发明结合现有技术,从实际应用出发,提供一种基于动态阈值模型的电动汽车故障诊断方法。
本发明的技术方案如下:
一种基于动态阈值模型的电动汽车故障诊断方法,所述方法用于电动汽车中电池系统故障诊断,所述方法包括:构建动态阈值模型以及基于阈值模型进行电池系统故障诊断,其中,
构建动态阈值模型包括如下步骤:
步骤S11、对锂电池动力电池进行充放电老化循环实验、容量实验、开路电压OCV实验以及混合动力脉冲能力实验,分别用于获取对应循环下的容量Q、OCV-SOC关系式和电池等效电路模型参数;
步骤S12、建立电池全寿命周期下的OCV-SOC-Q三维响应面模型;
步骤S13、建立基于欧姆内阻参数Ro以及时间常数参数τ的一阶等效电路模型,并进行模型参数辨识,建立动态阈值模型;
步骤S14、在不同温度下对动力电池进行上述步骤,得到不同温度下的阈值模型;基于阈值模型进行电池系统故障诊断包括如下步骤:
步骤S21、记录电池充放电过程的端电压信号、电流信号和温度信号;
步骤S22、基于建立的一阶等效电路模型,辨识电池模型参数和状态,得到电池内阻参数Ro、时间常数参数τ、容量Q和SOC的估计值;
步骤S23、采用温度差值的方法,确定参数参考值;
步骤S24、根据步骤S22得到的参数估计值和步骤S23得到的参数参考值,得到故障诊断残差;
步骤S25、分别确定参数Ro及参数τ的阈值;
步骤S26、对比残差和阈值判断电池是否发生故障。
进一步,步骤S11中,在老化实验过程中,在第0个循环后及每50个老化循环后进行一次容量实验、开路电压OCV实验以及混合动力脉冲能力实验。
进一步,步骤S12中,根据第0个循环后的OCV实验建立OCV与SOC的如下关系式:
Figure BDA0002750406200000031
式中,Uoc表示开路电压OCV;α012...α6为拟合系数,可通过OCV实验拟合得到;z为电池SOC,其值在电池完全充满时为100%,完全放光时为0,其他时刻的值可通过安时积分法计算:
Figure BDA0002750406200000032
其中z0为SOC初始值,iL为电池电流;
对每50个老化循环后的OCV实验数据进行处理得到每50个循环对应容量下的Uoc(z),然后将系数α012...α6进一步解析成容量Q的二次函数以得到电池全寿命周期下的OCV-SOC-Q三维响应面模型Uoc(Q,z):
Figure BDA0002750406200000033
式中,上标T表示矩阵的转置;Ω为7×3的拟合系数矩阵。
进一步,步骤S13中,一阶RC等效电路模型的数学表达式如下:
Figure BDA0002750406200000034
式中,Up为电池极化电容两端的极化电压,Ut为端电压,Rp为极化内阻,Cp为极化电容,Ro为欧姆内阻,时间常数τ=Rp×Cp
对第0次循环后的混合动力脉冲能力特性实验数据采用带遗忘因子的递推最小二乘法进行模型参数辨识,得到第0次循环后电池SOC从100%下降至0全过程对应的参数R0,0和τ0的估计值,即得到SOC分别与Ro,0和τ0的映射关系;
然后再对第50个循环后的混合动力脉冲能力特性实验数据进行参数辨识,得到对应循环下电池放电过程的SOC与Ro,50和τ50的映射关系,将对应SOC点下的Ro,50和τ50的估计值分别与Ro,0和τ0作差得到参数误差ΔRo,50和Δτ50,然后建立参数误差和SOC的如下拟合式:
Figure BDA0002750406200000041
Figure BDA0002750406200000042
式中,β01...β6和γ01...γ6分别为ΔRo,50和Δτ50与SOC的拟合系数;
而后将β01...β6和γ01...γ6分别进一步解析成容量Q的二次函数,得到全寿命周期内的ΔR0(Q,z)和Δτ(Q,z):
Figure BDA0002750406200000043
Figure BDA0002750406200000044
式中,A1和A2均为7×3的拟合系数矩阵;
定义参数Ro对应的阈值模型J1(Q,z)为:
J1(Q,z)=1.1*ΔRo(Q,z) (8)
定义参数τ对应的阈值模型J2(Q,z)为:
J2(Q,z)=1.1*Δτ(Q,z) (9)。
进一步,步骤S14中,将不同批次的动力电池分别在0℃、10℃、20℃、30℃、40℃和50℃环境为目的下进行步骤S11至步骤S13,得到不同温度下SOC分别与Ro,0和τ0的映射关系及三个模型:Uoc(Q,z)、J1(Q,z)和J2(Q,z),用上标Ta表示不同温度下的三个模型,即
Figure BDA0002750406200000053
Ta取集合[-10℃ 0℃ 10℃ 20℃ 30℃ 40℃ 50℃]中的某一值。
进一步,步骤S22中,基于建立的一阶等效电路模型,采用双扩展卡尔曼滤波算法辨识电池模型参数和状态,得到电池Ro、τ、Q和SOC,具体过程如下:
选取极化电压Up和SOC为状态向量,x=[UP SOC]T;参数向量θ=[Ro Rp Cp Q]T;电流为输入向量,u=[iL];模型端电压为输出向量y=[Ut],状态噪声为wk,参数噪声为rk,观测噪声为vk,参数和状态联合估计所用离散状态空间形式为:
Figure BDA0002750406200000051
式中,下标k表示第k个采样时刻;wk、rk、和vk为互相独立的零均值高斯白噪声,且误差协方差矩阵分别为:∑w、∑r和∑v
定义雅克比系数矩阵:
Figure BDA0002750406200000052
式中,Δt为采样间隔,上标-表示先验估计;
完成上述状态向量、参数向量、输入、输出和状态空间方程后,按照双扩展卡尔曼滤波算法完成参数和状态的联合估计,得到当前温度下每个采样时刻的Ro,k、τk、Qk和SOCk
进一步,步骤S23中,在采样时刻k,基于电池当前的温度Tex、SOCk和容量估计值Qk确定对应的参数参考值Ro,0,k和τ0,k的具体确定过程如下:
根据电池当前温度Tex读数找到Ta集合中与其最接近的两个温度Ta1和Ta2,且Ta1>Ta2,根据Ta1和Ta2温度下的SOC与Ro,0的映射关系,结合估计出的SOCk值,确定出Ta1和Ta2温度下的Ro,0,k,Ta1和Ro,0,k,Ta2,则在k采样时刻,Ro,k对应的参考值Ro,0,k通过下式确定:
Figure BDA0002750406200000061
通过相同的方式可确定k采样时刻,τk对应的参考值τ0,k,r
Figure BDA0002750406200000062
式中,
Figure BDA0002750406200000063
Figure BDA0002750406200000064
分别为SOCk时Ta1和Ta2温度下τ0的值。
进一步,步骤S24中,根据参数估计值和参数参考值,得到如下两个故障诊断残差:
Figure BDA0002750406200000065
进一步,步骤S25中,首先根据电池温度Tex读数找到Ta集合中与其最接近的两个温度Ta1和Ta2,且Ta1>Ta2,结合估计出的SOCk值,确定出Ta1和Ta2温度下的阈值模型
Figure BDA0002750406200000066
中的阈值
Figure BDA0002750406200000067
Figure BDA0002750406200000068
Figure BDA0002750406200000069
中的阈值
Figure BDA00027504062000000610
Figure BDA00027504062000000611
最终利用得到k时刻的阈值如下:
Figure BDA00027504062000000612
进一步,步骤S26中,通过对比残差和阈值来判断电池是否发生故障:
Figure BDA00027504062000000613
时,电池系统出现故障,反之,系统无故障。
本发明的有益效果:
(1)本发明所提方法在残差生成过程,充分考虑温度和SOC影响,且考虑到老化、温度和SOC影响,提出了动态阈值模型代替传统单一恒定阈值,故障诊断率高;
(2)本发明提出采用欧姆内阻和时间常数两个参数分别生成残差,当其中一个残差超过对应阈值即可判定电池系统出现故障,避免了单一残差进行故障诊断带来的检测不及时、误警和漏警问题。
附图说明
附图1为本发明流程框图。
附图2为本发明一阶等效电路图。
具体实施方式
结合附图和具体实施例,对本发明作进一步说明。应理解,这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解,在阅读了本发明讲授的内容之后,本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样落于本申请所限定的范围。
本发明提供了一种基于动态阈值模型的电动汽车故障诊断方法,该方法用于电动汽车中电池系统故障诊断,本发明所提方法需先构建动态阈值模型,然后基于阈值模型进行电池系统故障诊断。
其中,动态阈值模型的构建过程如下:
步骤S11、确定所用锂离子动力电池的型号及技术参数,根据电池企业提供的手册或相关标准对编号为01的电池在-10℃的温度下展开老化循环实验,在老化实验过程中,需在第0个循环后(即新电池)及每50个老化循环后进行一次电池基础特性测试实验。电池特性测试实验包括容量测试、开路电压(Open circuit voltage,OCV)测试和混合动力脉冲能力特性(Hybrid pulse power characteristic,HPPC)测试,分别用于获取对应循环下的容量Q、OCV-SOC关系式和电池等效电路模型参数。
步骤S12、根据第0个循环后的OCV实验建立OCV与SOC的如下关系式:
Figure BDA0002750406200000071
式中,Uoc表示开路电压OCV;α012...α6为拟合系数,可通过OCV实验拟合得到;z为电池SOC,其值在电池完全充满时为100%,完全放光时为0,其他时刻的值可通过安时积分法计算:
Figure BDA0002750406200000081
其中z0为SOC初始值,iL为电池电流。
对每50个老化循环后的基础特性测试中的OCV实验数据进行处理得到每50个循环对应容量下的Uoc(z),然后将系数α012...α6进一步解析成容量Q的二次函数以得到电池全寿命周期(即每个容量点)下的OCV-SOC-Q三维响应面模型Uoc(Q,z):
Figure BDA0002750406200000082
式中,上标T表示矩阵的转置;Ω为7×3的拟合系数矩阵。
步骤S13、建立图2所示的一阶RC等效电路模型,模型数学表达式如下:
Figure BDA0002750406200000083
式中,Up为电池极化电容两端的极化电压,Ut为端电压,Rp为极化内阻,Cp为极化电容,Ro为欧姆内阻。时间常数τ=Rp×Cp
对第0次循环后的HPPC实验数据采用带遗忘因子的递推最小二乘法进行模型参数辨识,得到第0次循环后电池SOC从100%下降至0全过程对应的参数R0,0和τ0的估计值,即得到SOC分别与Ro,0和τ0的映射关系。然后再对第50个循环后的HPPC实验数据进行参数辨识得到对应循环下电池放电过程的SOC与Ro,50和τ50的映射关系。将对应SOC点下的Ro,50和τ50的估计值分别与Ro,0和τ0作差得到参数误差ΔRo,50和Δτ50,然后建立参数误差和SOC的如下拟合式:
Figure BDA0002750406200000084
Figure BDA0002750406200000085
式中,β01...β6和γ01...γ6分别为ΔRo,50和Δτ50与SOC的拟合系数。
因不同循环对应不同老化状态(即不同容量),可利用此方法获得其余每隔50个老化循环下的参数误差与SOC的拟合式,而后将β01...β6和γ01...γ6分别进一步解析成容量Q的二次函数,得到全寿命周期内的ΔR0(Q,z)和Δτ(Q,z):
Figure BDA0002750406200000091
Figure BDA0002750406200000092
式中,A1和A2均为7×3的拟合系数矩阵。
定义参数Ro对应的阈值模型J1(Q,z)为:
J1(Q,z)=1.1*ΔRo(Q,z) (8)
定义参数τ对应的阈值模型J2(Q,z)为:
J2(Q,z)=1.1*Δτ(Q,z) (9)
步骤S14、将编号为02~07的六只同批次动力电池分别在0℃、10℃、20℃、30℃、40℃和50℃环境为目的下进行上述步骤S11至步骤S13,得到不同温度下SOC分别与Ro,0和τ0的映射关系及三个模型:Uoc(Q,z)、J1(Q,z)和J2(Q,z)。为后续叙述方便,用上标Ta表示不同温度下的三个模型,即
Figure BDA0002750406200000093
Figure BDA0002750406200000094
Ta取集合[-10℃ 0℃ 10℃ 20℃ 30℃ 40℃ 50℃]中的某一值。
电池系统实际故障诊断过程步骤如下:
步骤S21、记录电池充放电过程的端电压信号Ut、电流信号iL和温度信号Tex
步骤S22、基于建立的一阶等效电路模型,采用应用较多的双扩展卡尔曼滤波算法辨识电池模型参数和状态,得到电池Ro、τ、Q和SOC,具体过程如下:
(i)、选取极化电压Up和SOC为状态向量,x=[UP SOC]T;参数向量θ=[Ro Rp Cp Q]T;电流为输入向量,u=[iL];模型端电压为输出向量y=[Ut],状态噪声为wk,参数噪声为rk,观测噪声为vk,参数和状态联合估计所用离散状态空间形式为:
Figure BDA0002750406200000101
式中,下标k表示第k个采样时刻;wk、rk、和vk为互相独立的零均值高斯白噪声,且误差协方差矩阵分别为:∑w、∑r和∑v
定义雅克比系数矩阵:
Figure BDA0002750406200000102
式中,Δt为采样间隔,上标-表示先验估计。
完成上述状态向量、参数向量、输入、输出和状态空间方程后,即可按照双扩展卡尔曼滤波算法完成参数和状态的联合估计,得到当前温度下每个采样时刻的Ro,k、τk、Qk和SOCk
步骤S23、在采样时刻k,基于电池当前的温度Tex、SOCk和容量估计值Qk确定对应的参数参考值Ro,0,k和τ0,k的具体确定过程如下:
根据电池当前温度Tex读数找到Ta集合中与其最接近的两个温度Ta1和Ta2,且Ta1>Ta2,根据Ta1和Ta2温度下的SOC与Ro,0的映射关系,结合步骤S22中估计出的SOCk值,确定出Ta1和Ta2温度下的Ro,0,k,Ta1和Ro,0,k,Ta2,则在k采样时刻,Ro,k对应的参考值Ro,0,k通过下式确定:
Figure BDA0002750406200000103
通过相同的方式可确定k采样时刻,τk对应的参考值τ0,k,r
Figure BDA0002750406200000104
式中,
Figure BDA0002750406200000105
Figure BDA0002750406200000106
分别为SOCk时Ta1和Ta2温度下τ0的值。
步骤S24、根据步骤S22得到的参数估计值和步骤S23得到的参数参考值,可得到如下两个故障诊断残差:
Figure BDA0002750406200000111
步骤S25、确定故障诊断阈值过程与参数参考值确定过程相似,首先根据电池温度Tex读数找到Ta集合中与其最接近的两个温度Ta1和Ta2,且Ta1>Ta2,结合步骤II中估计出的SOCk值,确定出Ta1和Ta2温度下的阈值模型
Figure BDA0002750406200000112
中的阈值
Figure BDA0002750406200000113
Figure BDA0002750406200000114
Figure BDA0002750406200000115
中的阈值
Figure BDA0002750406200000116
Figure BDA0002750406200000117
最终利用得到k时刻的阈值如下:
Figure BDA0002750406200000118
步骤S26、通过对比残差和阈值来判断电池是否发生故障:
Figure BDA0002750406200000119
时,电池系统出现故障,反之,系统无故障。本发明在阈值模型建立和参数辨识算法两方面进行了改进。在不同温度下进行电路基础特性测试实验,得到等效电路模型参数;建立OCV-SOC-Q三维响应面模型;采用带遗忘因子的递推最小二乘法进行模型参数辨识,建立关于R0和τ的动态阈值模型。在实际故障诊断过程当中,利用双扩展卡尔曼滤波算法辨识参数和状态,得到电池R0和τ、容量及SOC;采用温度插值的方法确定参数参考值;确定参数阈值;生成残差;通过对比残差与阈值来判断电池是否发生故障。该方法不仅故障诊断率高,还能避免检测不及时、误警和漏警问题。

Claims (10)

1.一种基于动态阈值模型的电动汽车故障诊断方法,其特征在于,所述方法用于电动汽车中电池系统故障诊断,所述方法包括:构建动态阈值模型以及基于阈值模型进行电池系统故障诊断,其中,
构建动态阈值模型包括如下步骤:
步骤S11、对锂电池动力电池进行充放电老化循环实验、容量实验、开路电压OCV实验以及混合动力脉冲能力实验,分别用于获取对应循环下的容量Q、OCV-SOC关系式和电池等效电路模型参数;
步骤S12、建立电池全寿命周期下的OCV-SOC-Q三维响应面模型;
步骤S13、建立基于欧姆内阻参数Ro以及时间常数参数τ的一阶等效电路模型,并进行模型参数辨识,建立动态阈值模型;
步骤S14、在不同温度下对动力电池进行上述步骤,得到不同温度下的阈值模型;基于阈值模型进行电池系统故障诊断包括如下步骤:
步骤S21、记录电池充放电过程的端电压信号、电流信号和温度信号;
步骤S22、基于建立的一阶等效电路模型,辨识电池模型参数和状态,得到电池内阻参数Ro、时间常数参数τ、容量Q和SOC的估计值;
步骤S23、采用温度差值的方法,确定参数参考值;
步骤S24、根据步骤S22得到的参数估计值和步骤S23得到的参数参考值,得到故障诊断残差;
步骤S25、分别确定参数Ro及参数τ的阈值;
步骤S26、对比残差和阈值判断电池是否发生故障。
2.根据权利要求1所述的基于动态阈值模型的电动汽车故障诊断方法,其特征在于,步骤S11中,在老化实验过程中,在第0个循环后及每50个老化循环后进行一次容量实验、开路电压OCV实验以及混合动力脉冲能力实验。
3.根据权利要求2所述的基于动态阈值模型的电动汽车故障诊断方法,其特征在于,步骤S12中,根据第0个循环后的OCV实验建立OCV与SOC的如下关系式:
Figure FDA0002750406190000021
式中,Uoc表示开路电压OCV;α012...α6为拟合系数,可通过OCV实验拟合得到;z为电池SOC,其值在电池完全充满时为100%,完全放光时为0,其他时刻的值可通过安时积分法计算:
Figure FDA0002750406190000022
其中z0为SOC初始值,iL为电池电流;
对每50个老化循环后的OCV实验数据进行处理得到每50个循环对应容量下的Uoc(z),然后将系数α012...α6进一步解析成容量Q的二次函数以得到电池全寿命周期下的OCV-SOC-Q三维响应面模型Uoc(Q,z):
Figure FDA0002750406190000023
式中,上标T表示矩阵的转置;Ω为7×3的拟合系数矩阵。
4.根据权利要求3所述的基于动态阈值模型的电动汽车故障诊断方法,其特征在于,步骤S13中,一阶RC等效电路模型,模型数学表达式如下:
Figure FDA0002750406190000024
式中,Up为电池极化电容两端的极化电压,Ut为端电压,Rp为极化内阻,Cp为极化电容,Ro为欧姆内阻,时间常数τ=Rp×Cp
对第0次循环后的混合动力脉冲能力特性实验数据采用带遗忘因子的递推最小二乘法进行模型参数辨识,得到第0次循环后电池SOC从100%下降至0全过程对应的参数R0,0和τ0的估计值,即得到SOC分别与Ro,0和τ0的映射关系;
然后再对第50个循环后的混合动力脉冲能力特性实验数据进行参数辨识得到对应循环下电池放电过程的SOC与Ro,50和τ50的映射关系,将对应SOC点下的Ro,50和τ50的估计值分别与Ro,0和τ0作差得到参数误差ΔRo,50和Δτ50,然后建立参数误差和SOC的如下拟合式:
Figure FDA0002750406190000031
Figure FDA0002750406190000032
式中,β01...β6和γ01...γ6分别为ΔRo,50和Δτ50与SOC的拟合系数;
而后将β01...β6和γ01...γ6分别进一步解析成容量Q的二次函数,得到全寿命周期内的ΔR0(Q,z)和Δτ(Q,z):
Figure FDA0002750406190000033
Figure FDA0002750406190000034
式中,A1和A2均为7×3的拟合系数矩阵;
定义参数Ro对应的阈值模型J1(Q,z)为:
J1(Q,z)=1.1*ΔRo(Q,z) (8)
定义参数τ对应的阈值模型J2(Q,z)为:
J2(Q,z)=1.1*Δτ(Q,z) (9)。
5.根据权利要求4所述的基于动态阈值模型的电动汽车故障诊断方法,其特征在于,步骤S14中,将不同批次的动力电池分别在0℃、10℃、20℃、30℃、40℃和50℃环境为目的下进行步骤S11至步骤S13,得到不同温度下SOC分别与Ro,0和τ0的映射关系及三个模型:Uoc(Q,z)、J1(Q,z)和J2(Q,z),用上标Ta表示不同温度下的三个模型,即
Figure FDA0002750406190000035
Ta取集合[-10℃ 0℃ 10℃ 20℃ 30℃ 40℃ 50℃]中的某一值。
6.根据权利要求5所述的基于动态阈值模型的电动汽车故障诊断方法,其特征在于,步骤S22中,基于建立的一阶等效电路模型,采用双扩展卡尔曼滤波算法辨识电池模型参数和状态,得到电池Ro、τ、Q和SOC,具体过程如下:
选取极化电压Up和SOC为状态向量,x=[UP SOC]T;参数向量θ=[Ro Rp Cp Q]T;电流为输入向量,u=[iL];模型端电压为输出向量y=[Ut],状态噪声为wk,参数噪声为rk,观测噪声为vk,参数和状态联合估计所用离散状态空间形式为:
Figure FDA0002750406190000041
式中,下标k表示第k个采样时刻;wk、rk、和vk为互相独立的零均值高斯白噪声,且误差协方差矩阵分别为:∑w、∑r和∑v
定义雅克比系数矩阵:
Figure FDA0002750406190000042
式中,Δt为采样间隔,上标-表示先验估计;
完成上述状态向量、参数向量、输入、输出和状态空间方程后,按照双扩展卡尔曼滤波算法完成参数和状态的联合估计,得到当前温度下每个采样时刻的Ro,k、τk、Qk和SOCk
7.根据权利要求6所述的基于动态阈值模型的电动汽车故障诊断方法,其特征在于,步骤S23中,在采样时刻k,基于电池当前的温度Tex、SOCk和容量估计值Qk确定对应的参数参考值Ro,0,k和τ0,k的具体确定过程如下:
根据电池当前温度Tex读数找到Ta集合中与其最接近的两个温度Ta1和Ta2,且Ta1>Ta2,根据Ta1和Ta2温度下的SOC与Ro,0的映射关系,结合估计出的SOCk值,确定出Ta1和Ta2温度下的Ro,0,k,Ta1和Ro,0,k,Ta2,则在k采样时刻,Ro,k对应的参考值Ro,0,k通过下式确定:
Figure FDA0002750406190000051
通过相同的方式可确定k采样时刻,τk对应的参考值τ0kr
Figure FDA0002750406190000052
式中,
Figure FDA0002750406190000053
Figure FDA0002750406190000054
分别为SOCk时Ta1和Ta2温度下τ0的值。
8.根据权利要求7所述的基于动态阈值模型的电动汽车故障诊断方法,其特征在于,步骤S24中,根据参数估计值和参数参考值,得到如下两个故障诊断残差:
Figure FDA0002750406190000055
9.根据权利要求8所述的基于动态阈值模型的电动汽车故障诊断方法,其特征在于,步骤S25中,首先根据电池温度Tex读数找到Ta集合中与其最接近的两个温度Ta1和Ta2,且Ta1>Ta2,结合估计出的SOCk值,确定出Ta1和Ta2温度下的阈值模型
Figure FDA0002750406190000056
中的阈值
Figure FDA0002750406190000057
Figure FDA0002750406190000058
Figure FDA0002750406190000059
中的阈值
Figure FDA00027504061900000510
Figure FDA00027504061900000511
最终利用得到k时刻的阈值如下:
Figure FDA00027504061900000512
10.根据权利要求9所述的基于动态阈值模型的电动汽车故障诊断方法,其特征在于,步骤S26中,通过对比残差和阈值来判断电池是否发生故障:当
Figure FDA00027504061900000513
时,电池系统出现故障,反之,系统无故障。
CN202011181863.5A 2020-10-29 2020-10-29 一种基于动态阈值模型的电动汽车故障诊断方法 Active CN112285569B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011181863.5A CN112285569B (zh) 2020-10-29 2020-10-29 一种基于动态阈值模型的电动汽车故障诊断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011181863.5A CN112285569B (zh) 2020-10-29 2020-10-29 一种基于动态阈值模型的电动汽车故障诊断方法

Publications (2)

Publication Number Publication Date
CN112285569A true CN112285569A (zh) 2021-01-29
CN112285569B CN112285569B (zh) 2022-02-01

Family

ID=74353488

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011181863.5A Active CN112285569B (zh) 2020-10-29 2020-10-29 一种基于动态阈值模型的电动汽车故障诊断方法

Country Status (1)

Country Link
CN (1) CN112285569B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113687251A (zh) * 2021-08-23 2021-11-23 重庆邮电大学 一种基于双模型的锂离子电池组电压异常故障诊断方法
CN115201687A (zh) * 2022-07-13 2022-10-18 西南交通大学 基于在线宽频带阻抗的电池模型参数和SoC联合估计方法
CN116295934A (zh) * 2023-05-18 2023-06-23 中国华能集团清洁能源技术研究院有限公司 一种基于动态阈值的电池温度异常检测方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104297691A (zh) * 2014-09-29 2015-01-21 上海吉能电源系统有限公司 电池组健康状态诊断系统和方法
CN107831443A (zh) * 2017-10-20 2018-03-23 开沃新能源汽车集团有限公司 基于相关系数的电池系统短路故障诊断方法
CN108957349A (zh) * 2018-08-17 2018-12-07 北京航空航天大学 一种锂离子电池故障检测方法和系统
CN109143096A (zh) * 2018-09-04 2019-01-04 石家庄科林电气股份有限公司 电动自行车电池参数检测装置及检测方法
CN111257753A (zh) * 2020-03-10 2020-06-09 合肥工业大学 一种电池系统故障诊断方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104297691A (zh) * 2014-09-29 2015-01-21 上海吉能电源系统有限公司 电池组健康状态诊断系统和方法
CN107831443A (zh) * 2017-10-20 2018-03-23 开沃新能源汽车集团有限公司 基于相关系数的电池系统短路故障诊断方法
CN108957349A (zh) * 2018-08-17 2018-12-07 北京航空航天大学 一种锂离子电池故障检测方法和系统
CN109143096A (zh) * 2018-09-04 2019-01-04 石家庄科林电气股份有限公司 电动自行车电池参数检测装置及检测方法
CN111257753A (zh) * 2020-03-10 2020-06-09 合肥工业大学 一种电池系统故障诊断方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
RUI XIONG: "A Sensor Fault Diagnosis Method for a Lithium-IonBattery Pack in Electric Vehicles", 《IEEE TRANSACTIONS ON POWER ELECTRONICS》 *
刘征宇: "基于幅值平方相干谱的电池故障诊断方法", 《中 国 电 机 工 程 学 报 》 *
陈如尹: "基于改进阈值比较法的电动汽车电池故障在线诊断方法", 《电气技术》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113687251A (zh) * 2021-08-23 2021-11-23 重庆邮电大学 一种基于双模型的锂离子电池组电压异常故障诊断方法
CN115201687A (zh) * 2022-07-13 2022-10-18 西南交通大学 基于在线宽频带阻抗的电池模型参数和SoC联合估计方法
CN115201687B (zh) * 2022-07-13 2023-08-29 西南交通大学 基于在线宽频带阻抗的电池模型参数和SoC联合估计方法
CN116295934A (zh) * 2023-05-18 2023-06-23 中国华能集团清洁能源技术研究院有限公司 一种基于动态阈值的电池温度异常检测方法及系统
CN116295934B (zh) * 2023-05-18 2023-08-04 中国华能集团清洁能源技术研究院有限公司 一种基于动态阈值的电池温度异常检测方法及系统

Also Published As

Publication number Publication date
CN112285569B (zh) 2022-02-01

Similar Documents

Publication Publication Date Title
CN112285569B (zh) 一种基于动态阈值模型的电动汽车故障诊断方法
CN111610456B (zh) 一种区分电池微短路和小容量故障的诊断方法
CN111198327B (zh) 一种单体电池内短路故障自检测方法
CN110261779B (zh) 一种三元锂电池荷电状态与健康状态在线协同估计方法
CN107368619B (zh) 一种扩展卡尔曼滤波soc估算方法
CN111965547B (zh) 一种基于参数辨识法的电池系统传感器故障诊断方法
CN107748336B (zh) 锂离子电池的荷电状态在线估计方法及系统
CN108732508B (zh) 一种锂离子电池容量的实时估计方法
CN112285570B (zh) 一种基于衰减记忆滤波器的电动汽车故障诊断方法
CN111929602B (zh) 一种基于容量估计的单体电池漏电或微短路定量诊断方法
CN112147514B (zh) 基于rls的锂电池全工况自适应等效电路模型
CN112462282A (zh) 基于机理模型的用于确定电池组实时荷电状态的方法
CN111044924B (zh) 一种全工况电池剩余容量确定方法及系统
CN113484771A (zh) 一种锂离子电池宽温度全寿命soc及容量估计的方法
CN115877232A (zh) 一种基于卡尔曼滤波的锂离子电池内部温度估计方法
CN115792638A (zh) 基于电池模型参数辨识的soc-内短路联合估计方法
CN112009252B (zh) 一种动力电池系统故障诊断及容错控制方法
CN114720881A (zh) 一种基于改进初值带遗忘因子递推最小二乘法的锂电池参数辨识方法
CN112213644B (zh) 电池荷电状态估算方法及电池管理系统
CN116400248B (zh) 电池内短路故障的确定方法及装置
CN109298340B (zh) 一种基于可变时间尺度的电池容量在线估计方法
CN110927582B (zh) 一种基于多种sigma点的锂电池SOC估计方法
CN111965548A (zh) 一种基于状态估计法的电池系统传感器故障诊断方法
CN117007969A (zh) 可重构电池系统中单个模组soc确定方法、系统及设备
CN114152826B (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