CN111638462A - 一种soc-ocv分段拟合方法 - Google Patents

一种soc-ocv分段拟合方法 Download PDF

Info

Publication number
CN111638462A
CN111638462A CN202010289519.1A CN202010289519A CN111638462A CN 111638462 A CN111638462 A CN 111638462A CN 202010289519 A CN202010289519 A CN 202010289519A CN 111638462 A CN111638462 A CN 111638462A
Authority
CN
China
Prior art keywords
soc
battery
temperature
ocv
fitting
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
CN202010289519.1A
Other languages
English (en)
Other versions
CN111638462B (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN202010289519.1A priority Critical patent/CN111638462B/zh
Publication of CN111638462A publication Critical patent/CN111638462A/zh
Application granted granted Critical
Publication of CN111638462B publication Critical patent/CN111638462B/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/385Arrangements for measuring battery or accumulator variables
    • G01R31/387Determining ampere-hour charge capacity or SoC
    • G01R31/388Determining ampere-hour charge capacity or SoC involving voltage measurements
    • 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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Secondary Cells (AREA)
  • Tests Of Electric Status Of Batteries (AREA)

Abstract

本发明公开了一种SOC‑OCV分段拟合方法,首先,选定某动力电池然后该动力电池进行恒温脉冲充放电实验,记录实验数据,其中包括静置前的电池的正极温度Ta、电池体中间表面的温度Tb和电池的负极温度Tc;接着,基于常规的OCV‑SOC模型确定分段点,再进行分段拟合;然后,进行变环境温度高倍率、长时间连续充放电,拟合出电芯温度和Ta、Tb、Tc的关系模型,再进行考虑电芯温度的SOC‑OCV分段曲线拟合。本发明给出了两种分段曲线模型可供于不同的工况进行SOC估计,能够提高精度。

Description

一种SOC-OCV分段拟合方法
技术领域
本发明涉及汽车锂离子动力电池管理技术领域,尤其涉及一种SOC-OCV分段拟合方法。
背景技术
随着全球燃油汽车保有量迅速增长,能源以及燃油汽车带来的环境污染成为一个很大的问题,这对各个国家的能源结构和自然环境提出了很大的挑战。越来越多的国家,尤其是汽车工业强国,纷纷都加大清洁能源汽车的投资和研发力度,并给予相应的政策导向和扶持。因此,电动汽车在近几年成为解决此问题的一大选择并得到飞速发展。目前,对于电动汽车而言,锂离子电池技术是其中最重要的技术之一。动力电池的强时变非线性,以及诸多因,使得动力电池管理系统(BMS)对动力电池荷电状态(SOC)很难进行精确的估计。SOC-OCV关系曲线的拟合作为SOC估计中非常重要的一环,它所表示的电池荷电量(SOC)与电池开路电压(OCV)关系的准确度显得尤其重要。因此获得一种适用于多种情况的SOC-OCV关系曲线是精确估计电池SOC的基础。
目前电池SOC-OCV的测定方法,主要是SOC从0到100%之间每间隔10%的SOC测定一次电池的OCV,然后代入公式进行简单的拟合,并不能保证拟合的关系曲线的精度,而且考虑的影响因素不足。例如中国实用新型专利申请号为201811540158.2,名称为“一种SOC-OCV关系估算方法”中提出一种具有不温度修正的SOC-OCV曲线拟合方法;中国实用新型专利申请号为201811540158.2,名称为“一种锂离子电池SOC-OCV曲线的标定方法”中提出一种缩短了三元锂离子电池测试周期的SOC-OCV拟合方法。他们都未能提出一种能够考虑电池实际工作时各个部位的温度差异并且能够有效精确拟合SOC-OCV关系曲线的方法。本发明提出一种考虑电池各部位稳固差异、环境温度并SOC-OCV曲线分段拟合方法。
发明内容
本发明所要解决的技术问题是针对背景技术中所涉及到的缺陷,提供一种SOC-OCV分段拟合方法。
本发明为解决上述技术问题采用以下技术方案:
一种SOC-OCV曲线分段拟合方法,包含以下步骤:
步骤1),将电池温度与环境温度进行比较,获得其差值;
步骤2),当电池温度与环境温度的差值的绝对值小于预设的差值阈值时:
步骤2.1),采集分段拟合的试验数据,令M=5%:
步骤2.1.1),采用0.5C恒流放电的方法将电池放电至截至电压;
步骤2.1.2),按照标称充电电流给电池充电,充入当前温度下电池最大可用充电容量M的电量,记录此时电池的正极温度Ta、电池体中间表面的温度Tb和电池的负极温度Tc
步骤2.1.3),静置1小时以消除极化现象,然后测量此时的端电压作为电池在当前SOC状态下的充电开路电压OCV;
步骤2.1.4),令M=M+5%;
步骤2.1.5),重复步骤2.1.1)至步骤2.1.4),直至M大于100%;
步骤2.2),根据步骤2.1)中采集的试验数据,根据以下公式拟合SOC-OCV曲线:
UOC=k0+k1SOC+k2SOC2+k3SOC3+k4/SOC+k51nSOC+k61n(1-SOC)
式中,UOC为拟合后的电池开路电压,SOC为实验采集的电池荷电量,k1、k2、k3、k4、k5、k6均为需要拟合的多项式系数;
步骤2.3),根据以下公式计算步骤2.2)中拟合出的SOC-OCV曲线斜率的变化率:
Figure BDA0002449841970000021
其中,OCV为采集的电池开路电压,SOC为采集的电池荷电量,α为曲线斜率变化率的导数;
步骤2.4),根据α确定分段点A、B:
步骤2.4.1),当α<d时,取曲线上此时的点为A0;当α>d时,取曲线上此时的点为B0,d为预先设定的斜率阈值;
步骤2.4.2),选择采集的实验数据的样本点中与A0和B0最靠近的两个点分别作为分段点A,B,令A点的SOC值定为S1,B点的SOC值定为S2
步骤2.5),根据A、B点,进行SOC-OCV曲线的分段拟合:
步骤2.5.1),对于采集的实验数据,将其中SOC值小于等于S1的数据划分为第一部分,SOC值大于S1且小于S2的数据作为第二部分,SOC值大于等于S2的数据作为第三部分;
步骤2.5.2),将第一至第三部分的数据分别代入三个新的多项式进行拟合;
第一部分:UOC1=k10+k11SOC+k12SOC2+k13SOC3+k14SOC4
第二部分:UOC2=k20+k21SOC+k22SOC2+k23SOC3
第三部分:UOC3=k30+k31SOC+k32SOC2+k33SOC3+k34SOC4
其中,UOC1、UOC2、UOC3分别为第一部分、第二部分、第三部分数据拟合后的电池开路电压,kij为需要拟合的系数,i=1,2,3;j=1,2,3,4;
步骤2.5.3),将第一至第三部分数据拟合出的曲线整合到一起,得到SOC-OCV分段拟合曲线
Figure BDA0002449841970000031
Figure BDA0002449841970000032
步骤3),当电池温度与环境温度的差值的绝对值大于等于预设的差值阈值时:
步骤3.1),采集不同环境温度下的电池状态数据;
步骤3.1.1),令T0为预设的常温温度阈值、大于等于0℃且小于等于20℃,令N=1;
步骤3.1.2),将电池至于温度为T0的环境中,先以0.5C倍率进行电池放电至截止电压,电池静止一小时后,进行电池1.5C充电0.5小时后停止,不进行电池静置以保持电池各部分的温度差,分别测量电池的电芯温度Tbat、电池的正极温度Ta、电池体中间表面的温度Tb和电池的负极温度Tc
步骤3.1.3),令T0=T0+5,N=N+1;
步骤3.1.4),重复步骤3.1.2)至步骤3.1.3),直至N大于6;
步骤3.2),得出实时的电池的电芯温度
Figure BDA0002449841970000034
的拟合函数:
将步骤3.1)中采集到的数据,代入下列多项式中拟合,求出多项式系数β1、β2、β3、β4、β5、β6
Tbat=β1Ta2Ta 23Tb4Tb 25Tc6Tc 2
式中,β1、β2、β3、β4、β5、β6为需要拟合的系数;
则实时的电池的电芯温度
Figure BDA0002449841970000035
的拟合函数为;
Figure BDA0002449841970000033
步骤3.3),采集分段拟合的试验数据,令M=5%:
步骤3.3.1),采用0.5C恒流放电的方法将电池放电至截至电压;
步骤3.3.2),按照标称充电电流给电池充电,充入当前温度下电池最大可用充电容量M的电量,记录此时电池的正极温度Ta、电池体中间表面的温度Tb和电池的负极温度Tc
步骤3.3.3),静置1小时以消除极化现象,然后测量此时的端电压作为电池在当前SOC状态下的充电开路电压OCV;
步骤3.3.4),令M=M+5%;
步骤3.3.5),重复步骤3.3.1)至步骤3.3.4),直至M大于100%;
步骤3.4),根据步骤3.3)得到的数据,进行考虑温度时的SOC-OCV曲线的分段拟合;
步骤3.4.1),将步骤3.3)中采集的数据代入下列公式中,拟合出常用的SOC-OCV的曲线:
U'OC=k'0+k1'SOC+k'2SOC2+k'3SOC3+k'4/SOC+k5'1nSOC+k'61n(1-SOC)
式中,U'OC为拟合后的电池开路电压,k'0、k1'、k'2、k3'、k'4、k5'、k6'均为需要拟合的多项式系数;
步骤3.4.2),根据以下公式计算步骤3.4)中拟合出的SOC-OCV曲线斜率的变化率:
Figure BDA0002449841970000041
式中,α'为曲线斜率变化率的导数;
步骤3.4.3),根据α'确定分段点A'、B':
步骤3.4.3.1),当α'<d时,取曲线上此时的点为A'0;当α'>d时,取曲线上此时的点为B'0,d为预先设定的斜率阈值;
步骤3.4.3.2),选择采集的实验数据的样本点中与A'0和B'0最靠近的两个点分别作为分段点A',B',令A'点的SOC值定为S1',B'点的SOC值定为S'2
步骤3.4.4),根据A'、B'点,进行SOC-OCV曲线的分段拟合:
步骤3.4.4.1),对于采集的实验数据,将其中SOC值小于等于S1'的数据划分为第一部分,SOC值大于S1'且小于S'2的数据作为第二部分,SOC值大于等于S'2的数据作为第三部分;
步骤3.4.4.2),将三部分的数据分别代入考虑电芯温度的分段多项式进行拟合:
第一部分:
Figure BDA0002449841970000042
第二部分:
Figure BDA0002449841970000043
第三部分:
Figure BDA0002449841970000044
Figure BDA0002449841970000045
Figure BDA0002449841970000046
其中,
Figure BDA0002449841970000047
为考虑电芯温度的分段开路电压,SOC为电池荷电量,ko10、ko20、ko30、aij1、aij2为需要拟合的系数;
步骤3.4.5),将各部分拟合出的曲线整合到一起,得到SOC-OCV分段拟合曲线:
Figure BDA0002449841970000051
当电池各部分温度与环境的温度关系满足下列关系:|Tb-T0|≥d1℃,|Ta-Tb|≥d2℃,|Ta-Tc|≥d3℃,|Tb-Tc|≥d4℃中任意一项时启用步骤3)中的SOC-OCV关系曲线。d1、d2、d3、d4为1左右的值,可根据具体情况调整。
其余电池温度与环境温度差异不大的情况可使用第一种未考虑电芯温度步骤2)中的分段SOC-OCV曲线。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
给出了两种分段SOC-OCV关系模型及选择方法,供于不同情况选择来提高SOC估计的精确度。
两种曲线相比常规关系曲线都更能准确的反应电池SOC与OCV的关系。温差不大的情况下,第一种未考虑电芯温度曲线能够适当减小计算量。温度差异变大时可以采用第二种考虑电芯温度的曲线,考虑的电芯温度更能准确的反应出电池真实的工作温度差异,避免了了传热效率,电池壳体材料等因素对于测量的电池温度的可能。
附图说明
图1为本发明专利的整体流程框图;
图2为A0、B0、A、B点的确定及分段拟合后的关系曲线;
图3为测量温度Tbat,Ta,Tb,Tc在电池中的位置。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
本发明可以以许多不同的形式实现,而不应当认为限于这里所述的实施例。相反,提供这些实施例以便使本公开透彻且完整,并且将向本领域技术人员充分表达本发明的范围。在附图中,为了清楚起见放大了组件。
如图1所示,本发明公开了一种SOC-OCV曲线分段拟合方法,包含以下步骤:
步骤1),将电池温度与环境温度进行比较,获得其差值;
步骤2),当电池温度与环境温度的差值的绝对值小于预设的差值阈值时:
步骤2.1),采集分段拟合的试验数据,令M=5%:
步骤2.1.1),采用0.5C恒流放电的方法将电池放电至截至电压;
步骤2.1.2),按照标称充电电流给电池充电,充入当前温度下电池最大可用充电容量M的电量,记录此时电池的正极温度Ta、电池体中间表面的温度Tb和电池的负极温度Tc
步骤2.1.3),静置1小时以消除极化现象,然后测量此时的端电压作为电池在当前SOC状态下的充电开路电压OCV;
步骤2.1.4),令M=M+5%;
步骤2.1.5),重复步骤2.1.1)至步骤2.1.4),直至M大于100%;
步骤2.2),根据步骤2.1)中采集的试验数据,根据以下公式拟合SOC-OCV曲线:
UOC=k0+k1SOC+k2SOC2+k3SOC3+k4/SOC+k51nSOC+k61n(1-SOC)
式中,UOC为拟合后的电池开路电压,SOC为实验采集的电池荷电量,k1、k2、k3、k4、k5、k6均为需要拟合的多项式系数;
步骤2.3),根据以下公式计算步骤2.2)中拟合出的SOC-OCV曲线斜率的变化率:
Figure BDA0002449841970000061
其中,OCV为采集的电池开路电压,SOC为采集的电池荷电量,α为曲线斜率变化率的导数;
步骤2.4),根据α确定分段点A、B,如图2所示:
步骤2.4.1),当α<d时,取曲线上此时的点为A0;当α>d时,取曲线上此时的点为B0,d为预先设定的斜率阈值;
步骤2.4.2),选择采集的实验数据的样本点中与A0和B0最靠近的两个点分别作为分段点A,B,令A点的SOC值定为S1,B点的SOC值定为S2
步骤2.5),根据A、B点,进行SOC-OCV曲线的分段拟合:
步骤2.5.1),对于采集的实验数据,将其中SOC值小于等于S1的数据划分为第一部分,SOC值大于S1且小于S2的数据作为第二部分,SOC值大于等于S2的数据作为第三部分;
步骤2.5.2),将第一至第三部分的数据分别代入三个新的多项式进行拟合;
第一部分:UOC1=k10+k11SOC+k12SOC2+k13SOC3+k14SOC4
第二部分:UOC2=k20+k21SOC+k22SOC2+k23SOC3
第三部分:UOC3=k30+k31SOC+k32SOC2+k33SOC3+k34SOC4
其中,UOC1、UOC2、UOC3分别为第一部分、第二部分、第三部分数据拟合后的电池开路电压,kij为需要拟合的系数,i=1,2,3;j=1,2,3,4;
步骤2.5.3),将第一至第三部分数据拟合出的曲线整合到一起,得到SOC-OCV分段拟合曲线
Figure BDA0002449841970000071
Figure BDA0002449841970000072
步骤3),当电池温度与环境温度的差值的绝对值大于等于预设的差值阈值时:
步骤3.1),采集不同环境温度下的电池状态数据;
步骤3.1.1),令T0为预设的常温温度阈值、大于等于0℃且小于等于20℃,令N=1;
步骤3.1.2),将电池至于温度为T0的环境中,先以0.5C倍率进行电池放电至截止电压,电池静止一小时后,进行电池1.5C充电0.5小时后停止,不进行电池静置以保持电池各部分的温度差,分别测量电池的电芯温度Tbat、电池的正极温度Ta、电池体中间表面的温度Tb和电池的负极温度Tc,如图3所示;
步骤3.1.3),令T0=T0+5,N=N+1;
步骤3.1.4),重复步骤3.1.2)至步骤3.1.3),直至N大于6;
步骤3.2),得出实时的电池的电芯温度
Figure BDA0002449841970000075
的拟合函数:
将步骤3.1)中采集到的数据,代入下列多项式中拟合,求出多项式系数β1、β2、β3、β4、β5、β6
Tbat=β1Ta2Ta 23Tb4Tb 25Tc6Tc 2
式中,β1、β2、β3、β4、β5、β6为需要拟合的系数;
则实时的电池的电芯温度
Figure BDA0002449841970000073
的拟合函数为;
Figure BDA0002449841970000074
步骤3.3),采集分段拟合的试验数据,令M=5%:
步骤3.3.1),采用0.5C恒流放电的方法将电池放电至截至电压;
步骤3.3.2),按照标称充电电流给电池充电,充入当前温度下电池最大可用充电容量M的电量,记录此时电池的正极温度Ta、电池体中间表面的温度Tb和电池的负极温度Tc
步骤3.3.3),静置1小时以消除极化现象,然后测量此时的端电压作为电池在当前SOC状态下的充电开路电压OCV;
步骤3.3.4),令M=M+5%;
步骤3.3.5),重复步骤3.3.1)至步骤3.3.4),直至M大于100%;
步骤3.4),根据步骤3.3)得到的数据,进行考虑温度时的SOC-OCV曲线的分段拟合;
步骤3.4.1),将步骤3.3)中采集的数据代入下列公式中,拟合出常用的SOC-OCV的曲线:
U'OC=k'0+k1'SOC+k'2SOC2+k'3SOC3+k'4/SOC+k5'1nSOC+k'61n(1-SOC)
式中,U'OC为拟合后的电池开路电压,k'0、k1'、k'2、k3'、k'4、k5'、k6'均为需要拟合的多项式系数;
步骤3.4.2),根据以下公式计算步骤3.4)中拟合出的SOC-OCV曲线斜率的变化率:
Figure BDA0002449841970000081
式中,α'为曲线斜率变化率的导数;
步骤3.4.3),根据α'确定分段点A'、B':
步骤3.4.3.1),当α'<d时,取曲线上此时的点为A'0;当α'>d时,取曲线上此时的点为B'0,d为预先设定的斜率阈值;
步骤3.4.3.2),选择采集的实验数据的样本点中与A'0和B'0最靠近的两个点分别作为分段点A',B',令A'点的SOC值定为S1',B'点的SOC值定为S'2
步骤3.4.4),根据A'、B'点,进行SOC-OCV曲线的分段拟合:
步骤3.4.4.1),对于采集的实验数据,将其中SOC值小于等于S1'的数据划分为第一部分,SOC值大于S1'且小于S'2的数据作为第二部分,SOC值大于等于S'2的数据作为第三部分;
步骤3.4.4.2),将三部分的数据分别代入考虑电芯温度的分段多项式进行拟合:
第一部分:
Figure BDA0002449841970000082
第二部分:
Figure BDA0002449841970000083
第三部分:
Figure BDA0002449841970000084
Figure BDA0002449841970000085
Figure BDA0002449841970000091
其中,
Figure BDA0002449841970000092
为考虑电芯温度的分段开路电压,SOC为电池荷电量,ko10、ko20、ko30、aij1、aij2为需要拟合的系数;
步骤3.4.5),将各部分拟合出的曲线整合到一起,得到考虑温度的SOC-OCV分段拟合曲线:
Figure BDA0002449841970000093
当电池各部分温度与环境的温度关系满足下列关系:|Tb-T0|≥d1℃,|Ta-Tb|≥d2℃,|Ta-Tc|≥d3℃,|Tb-Tc|≥d4℃中任意一项时启用步骤3)中的SOC-OCV关系曲线。d1、d2、d3、d4为1左右的值,可根据具体情况调整。
其余电池温度与环境温度差异不大的情况可使用步骤2)中的分段SOC-OCV曲线。
本技术领域技术人员可以理解的是,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (1)

1.一种SOC-OCV分段拟合方法,其特征在于,包含以下步骤:
步骤1),将电池温度与环境温度进行比较,获得其差值;
步骤2),当电池温度与环境温度的差值的绝对值小于预设的差值阈值时:
步骤2.1),采集分段拟合的试验数据,令M=5%:
步骤2.1.1),采用0.5C恒流放电的方法将电池放电至截至电压;
步骤2.1.2),按照标称充电电流给电池充电,充入当前温度下电池最大可用充电容量M的电量,记录此时电池的正极温度Ta、电池体中间表面的温度Tb和电池的负极温度Tc
步骤2.1.3),静置1小时以消除极化现象,然后测量此时的端电压作为电池在当前SOC状态下的充电开路电压OCV;
步骤2.1.4),令M=M+5%;
步骤2.1.5),重复步骤2.1.1)至步骤2.1.4),直至M大于100%;
步骤2.2),根据步骤2.1)中采集的试验数据,根据以下公式拟合SOC-OCV曲线:
UOC=k0+k1SOC+k2SOC2+k3SOC3+k4/SOC+k51nSOC+k61n(1-SOC)
式中,UOC为拟合后的电池开路电压,SOC为实验采集的电池荷电量,k1、k2、k3、k4、k5、k6均为需要拟合的多项式系数;
步骤2.3),根据以下公式计算步骤2.2)中拟合出的SOC-OCV曲线斜率的变化率:
Figure FDA0002449841960000011
其中,OCV为采集的电池开路电压,SOC为采集的电池荷电量,α为曲线斜率变化率的导数;
步骤2.4),根据α确定分段点A、B:
步骤2.4.1),当α<d时,取曲线上此时的点为A0;当α>d时,取曲线上此时的点为B0,d为预先设定的斜率阈值;
步骤2.4.2),选择采集的实验数据的样本点中与A0和B0最靠近的两个点分别作为分段点A,B,令A点的SOC值定为S1,B点的SOC值定为S2
步骤2.5),根据A、B点,进行SOC-OCV曲线的分段拟合:
步骤2.5.1),对于采集的实验数据,将其中SOC值小于等于S1的数据划分为第一部分,SOC值大于S1且小于S2的数据作为第二部分,SOC值大于等于S2的数据作为第三部分;
步骤2.5.2),将第一至第三部分的数据分别代入三个新的多项式进行拟合;
第一部分:UOC1=k10+k11SOC+k12SOC2+k13SOC3+k14SOC4
第二部分:UOC2=k20+k21SOC+k22SOC2+k23SOC3
第三部分:UOC3=k30+k31SOC+k32SOC2+k33SOC3+k34SOC4
其中,UOC1、UOC2、UOC3分别为第一部分、第二部分、第三部分数据拟合后的电池开路电压,kij为需要拟合的系数,i=1,2,3;j=1,2,3,4;
步骤2.5.3),将第一至第三部分数据拟合出的曲线整合到一起,得到SOC-OCV分段拟合曲线
Figure FDA0002449841960000021
Figure FDA0002449841960000022
步骤3),当电池温度与环境温度的差值的绝对值大于等于预设的差值阈值时:
步骤3.1),采集不同环境温度下的电池状态数据;
步骤3.1.1),令T0为预设的常温温度阈值、大于等于0℃且小于等于20℃,令N=1;
步骤3.1.2),将电池至于温度为T0的环境中,先以0.5C倍率进行电池放电至截止电压,电池静止一小时后,进行电池1.5C充电0.5小时后停止,不进行电池静置以保持电池各部分的温度差,分别测量电池的电芯温度Tbat、电池的正极温度Ta、电池体中间表面的温度Tb和电池的负极温度Tc
步骤3.1.3),令T0=T0+5,N=N+1;
步骤3.1.4),重复步骤3.1.2)至步骤3.1.3),直至N大于6;
步骤3.2),得出实时的电池的电芯温度
Figure FDA0002449841960000023
的拟合函数:
将步骤3.1)中采集到的数据,代入下列多项式中拟合,求出多项式系数β1、β2、β3、β4、β5、β6
Tbat=β1Ta2Ta 23Tb4Tb 25Tc6Tc 2
式中,β1、β2、β3、β4、β5、β6为需要拟合的系数;
则实时的电池的电芯温度
Figure FDA0002449841960000024
的拟合函数为;
Figure FDA0002449841960000025
步骤3.3),采集分段拟合的试验数据,令M=5%:
步骤3.3.1),采用0.5C恒流放电的方法将电池放电至截至电压;
步骤3.3.2),按照标称充电电流给电池充电,充入当前温度下电池最大可用充电容量M的电量,记录此时电池的正极温度Ta、电池体中间表面的温度Tb和电池的负极温度Tc
步骤3.3.3),静置1小时以消除极化现象,然后测量此时的端电压作为电池在当前SOC状态下的充电开路电压OCV;
步骤3.3.4),令M=M+5%;
步骤3.3.5),重复步骤3.3.1)至步骤3.3.4),直至M大于100%;
步骤3.4),根据步骤3.3)得到的数据,进行考虑温度时的SOC-OCV曲线的分段拟合;
步骤3.4.1),将步骤3.3)中采集的数据代入下列公式中,拟合出常用的SOC-OCV的曲线:
U'OC=k'0+k′1SOC+k'2SOC2+k'3SOC3+k'4/SOC+k′51nSOC+k'61n(1-SOC)
式中,U'OC为拟合后的电池开路电压,k'0、k′1、k'2、k′3、k'4、k′5、k′6均为需要拟合的多项式系数;
步骤3.4.2),根据以下公式计算步骤3.4)中拟合出的SOC-OCV曲线斜率的变化率:
Figure FDA0002449841960000031
式中,α'为曲线斜率变化率的导数;
步骤3.4.3),根据α'确定分段点A'、B':
步骤3.4.3.1),当α'<d时,取曲线上此时的点为A'0;当α'>d时,取曲线上此时的点为B'0,d为预先设定的斜率阈值;
步骤3.4.3.2),选择采集的实验数据的样本点中与A'0和B'0最靠近的两个点分别作为分段点A',B',令A'点的SOC值定为S′1,B'点的SOC值定为S'2
步骤3.4.4),根据A'、B'点,进行SOC-OCV曲线的分段拟合:
步骤3.4.4.1),对于采集的实验数据,将其中SOC值小于等于S′1的数据划分为第一部分,SOC值大于S′1且小于S'2的数据作为第二部分,SOC值大于等于S'2的数据作为第三部分;
步骤3.4.4.2),将三部分的数据分别代入考虑电芯温度的分段多项式进行拟合:
第一部分:
Figure FDA0002449841960000032
第二部分:
Figure FDA0002449841960000033
第三部分:
Figure FDA0002449841960000041
Figure FDA0002449841960000042
Figure FDA0002449841960000045
其中,
Figure FDA0002449841960000043
为考虑电芯温度的分段开路电压,SOC为电池荷电量,ko10、ko20、ko30、aij1、aij2为需要拟合的系数;
步骤3.4.5),将各部分拟合出的曲线整合到一起,得到SOC-OCV分段拟合曲线:
Figure FDA0002449841960000044
CN202010289519.1A 2020-04-14 2020-04-14 一种soc-ocv分段拟合方法 Active CN111638462B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010289519.1A CN111638462B (zh) 2020-04-14 2020-04-14 一种soc-ocv分段拟合方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010289519.1A CN111638462B (zh) 2020-04-14 2020-04-14 一种soc-ocv分段拟合方法

Publications (2)

Publication Number Publication Date
CN111638462A true CN111638462A (zh) 2020-09-08
CN111638462B CN111638462B (zh) 2021-05-11

Family

ID=72329665

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010289519.1A Active CN111638462B (zh) 2020-04-14 2020-04-14 一种soc-ocv分段拟合方法

Country Status (1)

Country Link
CN (1) CN111638462B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112305434A (zh) * 2020-09-15 2021-02-02 深圳新恒业电池科技有限公司 电池容量的测量方法、装置、终端及存储介质
CN113466696A (zh) * 2021-05-17 2021-10-01 北京交通大学 一种基于电压曲线变换的电池组单体状态估计方法
CN113740754A (zh) * 2021-09-06 2021-12-03 北京西清能源科技有限公司 一种检测电池组不一致性的方法与系统
CN117686918A (zh) * 2024-01-31 2024-03-12 深圳市卓芯微科技有限公司 电池soc预测方法、装置、电池管理设备及存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100121591A1 (en) * 2008-11-13 2010-05-13 Lockheed Martin Corporation Method and apparatus that detects state of charge (soc) of a battery
CN107741569A (zh) * 2017-11-16 2018-02-27 温州大学 一种基于分段扩展卡尔曼滤波的锂电池荷电状态的估算方法
CN108169682A (zh) * 2017-12-14 2018-06-15 江苏大学 一种基于气液动力学模型的锂离子电池soc精确估算方法
CN108829911A (zh) * 2018-04-16 2018-11-16 西南科技大学 一种开路电压与soc函数关系优化方法
CN109239608A (zh) * 2018-08-10 2019-01-18 安徽力高新能源技术有限公司 一种实时修正锂电池soc-ocv曲线的方法
CN110618389A (zh) * 2019-09-25 2019-12-27 宝能汽车有限公司 电池soc-ocv曲线的测试方法及测试电池soc-ocv曲线的设备

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100121591A1 (en) * 2008-11-13 2010-05-13 Lockheed Martin Corporation Method and apparatus that detects state of charge (soc) of a battery
CN107741569A (zh) * 2017-11-16 2018-02-27 温州大学 一种基于分段扩展卡尔曼滤波的锂电池荷电状态的估算方法
CN108169682A (zh) * 2017-12-14 2018-06-15 江苏大学 一种基于气液动力学模型的锂离子电池soc精确估算方法
CN108829911A (zh) * 2018-04-16 2018-11-16 西南科技大学 一种开路电压与soc函数关系优化方法
CN109239608A (zh) * 2018-08-10 2019-01-18 安徽力高新能源技术有限公司 一种实时修正锂电池soc-ocv曲线的方法
CN110618389A (zh) * 2019-09-25 2019-12-27 宝能汽车有限公司 电池soc-ocv曲线的测试方法及测试电池soc-ocv曲线的设备

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张佳倩等: "基于模型参数辨识和扩展卡尔曼滤波的锂电池荷电状态估计", 《工业控制计算机》 *
黄玉莎等: "基于CARMA模型的动力电池荷电状态估计", 《控制与信息技术》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112305434A (zh) * 2020-09-15 2021-02-02 深圳新恒业电池科技有限公司 电池容量的测量方法、装置、终端及存储介质
CN112305434B (zh) * 2020-09-15 2022-05-06 创普斯(深圳)新能源科技有限公司 电池容量的测量方法、装置、终端及存储介质
CN113466696A (zh) * 2021-05-17 2021-10-01 北京交通大学 一种基于电压曲线变换的电池组单体状态估计方法
CN113740754A (zh) * 2021-09-06 2021-12-03 北京西清能源科技有限公司 一种检测电池组不一致性的方法与系统
CN113740754B (zh) * 2021-09-06 2023-10-13 北京西清能源科技有限公司 一种检测电池组不一致性的方法与系统
CN117686918A (zh) * 2024-01-31 2024-03-12 深圳市卓芯微科技有限公司 电池soc预测方法、装置、电池管理设备及存储介质
CN117686918B (zh) * 2024-01-31 2024-05-07 深圳市卓芯微科技有限公司 电池soc预测方法、装置、电池管理设备及存储介质

Also Published As

Publication number Publication date
CN111638462B (zh) 2021-05-11

Similar Documents

Publication Publication Date Title
CN111638462B (zh) 一种soc-ocv分段拟合方法
CN107991623B (zh) 一种考虑温度和老化程度的电池安时积分soc估计方法
CN110501652B (zh) 一种退役锂电池可用容量快速评估方法及评估装置
CN111239611B (zh) 一种基于单体电池容量校准packsoc的计算方法
CN105467328A (zh) 一种锂离子电池荷电状态估计方法
CN111766530B (zh) 锂离子蓄电池单体寿命检测方法
CN108445422B (zh) 基于极化电压恢复特性的电池荷电状态估算方法
CN112485673B (zh) 一种基于动态多安全约束的电池充放电峰值功率预测方法
CN111366864B (zh) 一种基于固定压升区间的电池soh在线估计方法
CN110442901A (zh) 一种锂离子电池电化学简化模型及其参数的获取方法
CN112327180A (zh) 一种锂离子电池自放电评价及其电压均衡调节方法
CN114865117B (zh) 锂离子电池电极嵌锂量检测方法、装置及电池管理系统
CN104681851A (zh) 一种汽车用锂离子动力电池配组方法
CN108872867B (zh) 一种电池自放电测试方法
CN112986828A (zh) 预估电池健康度的方法
CN112819196B (zh) 一种等效库伦效率的测定方法及动力电池循环寿命的预测方法
CN109143096A (zh) 电动自行车电池参数检测装置及检测方法
CN116754971A (zh) 电池析锂检测的方法、设备及存储介质
CN115331743A (zh) 基于实验类比法的高倍率工况电化学模型建模方法
CN114879053A (zh) 一种储能磷酸铁锂电池寿命预测方法
JP7290126B2 (ja) 全固体電池の製造方法
CN111948554B (zh) 一种降低锂离子电池力学退化的方法
CN111832169B (zh) 一种电芯寿命模型自动校正方法
CN114609525A (zh) 基于分数阶容积卡尔曼滤波的动力电池soc估计方法
CN114035075A (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