CN105548901B - 一种轨道交通用钛酸锂电池功率状态预测方法 - Google Patents

一种轨道交通用钛酸锂电池功率状态预测方法 Download PDF

Info

Publication number
CN105548901B
CN105548901B CN201610009188.5A CN201610009188A CN105548901B CN 105548901 B CN105548901 B CN 105548901B CN 201610009188 A CN201610009188 A CN 201610009188A CN 105548901 B CN105548901 B CN 105548901B
Authority
CN
China
Prior art keywords
battery
soc
point
lithium titanate
rail traffic
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
Application number
CN201610009188.5A
Other languages
English (en)
Other versions
CN105548901A (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 Beijiao New Energy Technology Co Ltd
Original Assignee
Beijing Beijiao New Energy Technology Co Ltd
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 Beijiao New Energy Technology Co Ltd filed Critical Beijing Beijiao New Energy Technology Co Ltd
Priority to CN201610009188.5A priority Critical patent/CN105548901B/zh
Publication of CN105548901A publication Critical patent/CN105548901A/zh
Application granted granted Critical
Publication of CN105548901B publication Critical patent/CN105548901B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/382Arrangements for monitoring battery or accumulator variables, e.g. SoC
    • G01R31/3842Arrangements for monitoring battery or accumulator variables, e.g. SoC combining voltage and current measurements

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Secondary Cells (AREA)

Abstract

本发明涉及一种轨道交通用钛酸锂电池功率状态预测方法,基于电池改进模型,通过OCV‑SOC测试数据获得基准倍率下不同SOC的参数值,求取每个SOC点下典型倍率的充放电极化电压,基于复合形法测算电池特性参数,并根据该特性参数集合计算电池最大允许充放电电流和电池功率状态。本发明,提出一种能够实时准确预测电池功率状态的方法,所预测的电池功率状态不仅精度高,同时可直接用于系统级仿真,可充分利用电池功率能力以提升整车的动力特性,进而为轨道交通用钛酸锂电池运行过程中的最优控制策略提供理论依据。

Description

一种轨道交通用钛酸锂电池功率状态预测方法
技术领域
本发明涉及动力电池功率状态预测领域,具体说是一种轨道交通用钛酸锂电池功率状态预测方法。尤指基于高精度电池改进模型的轨道交通用钛酸锂电池功率状态预测方法。
背景技术
为了应对日益严重的环境污染和能源危机,电动车辆的开发及普及迫在眉睫。动力电池作为电动车辆动力系统的核心部件,其性能对车辆运行的安全性和高效性至关重要。
近年来,锂电池凭借自身在比功率、比能量及循环寿命方面的优势,被广泛应用于新能源车辆和电网储能领域。相对于采用石墨作为负极的功率型磷酸铁锂电池和锰酸锂电池,采用钛酸锂材料作为负极的钛酸锂电池具有更出色的倍率充放电特性、温度特性和更为稳定的结构,更适合用于频繁启停的电力驱动系统,成为未来电动车辆动力系统,尤其是轨道交通车辆的最佳选择之一。
作为电力驱动系统的储能部分,电池不仅需要在运行中提供一定量的能量,还要根据环境要求输出足够的功率。而电池功率受到电流、电压、荷电状态SOC、温度及时间间隔方面的约束条件的限制。当然,电池功率也会随着电池的老化而改变。通常将t秒内电池能够在安全运行区域内放电的最大功率称为电池t秒内最大放电功率,最大充电功率的定义类似。
目前关于电池功率状态预测的方法主要分为两种:基于特性图表的方法和基于动态电池模型的方法。
基于特性图表的方法主要通过在电池的不同SOC下施加特定脉冲组激励,从而得到相应的电压响应并进行功率预测。电池厂通常采用这类方法,以日本电动车辆协会标准JEVS D713 20035《混合动力电动汽车用密闭型镍氢电池的输出密度及输入密度试验方法》(简称JEVS)、美国Freedom-CAR项目《功率辅助型混合动力汽车用动力电池测试手册》(简称HPPC)为主,其优点是求取过程直接,方法简单。当然,其缺点就是仅考虑了电池静态特性,未能将实际运行中动态激励施加后初始极化状态的影响考虑在内,因此在动态工况中预测精度较低。
基于动态电池模型的方法则充分考虑电池动态特性,通过电池模型的参数值来预测电池当前状态下的功率状态。这类方法间的区别主要在于所采用的电池模型的类型。
作为预测电池功率状态的基础,电池模型的准确度直接影响到最终的预测精度。建立电池模型需要兼顾复杂程度和准确度,在保证模型精度的基础上使模型简单化。
目前常用的电池模型包括:电化学模型和电学模型,其中:
电化学模型充分表征了电池的材料特性以及反应机理,可以为电池优化设计提供依据,然而电化学模型中包含多个分数偏微分方程及大量未知变量,增加了模型的复杂程度,而且无法和动力系统其它部分直接连接。虽然电化学模型可以从多个角度进行简化,但简化的同时模型精度降低,并且简化后的模型仍难以满足应用的要求。
电学模型通过电学元件的组合来刻画电池动态特性,拥有更简单的结构和更少的未知参数,包括一阶RC电路模型等,更适合于电池性能仿真。虽然很多文献中针对电池在不同倍率激励下的表现提出了各种改进的电学模型,但这些模型均无法兼顾复杂程度和准确度,且均在窄倍率范围内进行验证,并未涉及钛酸锂电池。
因此有必要在兼顾复杂程度和准确度的电池改进模型的基础上,建立一种动态电池功率状态预测方法以准确描述钛酸锂电池高倍率应用下的电学行为和不同时刻的功率状态。
发明内容
针对现有技术中存在的缺陷,本发明的目的在于提供一种轨道交通用钛酸锂电池功率状态预测方法,在建立一种兼顾复杂程度和准确度的电池改进模型的基础上,进而提出一种基于动态电池改进模型的功率预测方法。
为达到以上目的,本发明采取的技术方案是:
一种轨道交通用钛酸锂电池功率状态预测方法,其特征在于:基于电池改进模型,通过OCV-SOC测试数据获得基准倍率下不同SOC的参数值,求取每个SOC点下典型倍率的充放电极化电压,基于复合形法测算电池特性参数,并根据该特性参数集合计算电池最大允许充放电电流和电池功率状态。
在上述技术方案的基础上,基准倍率取为1C。
在上述技术方案的基础上,所述电池改进模型,在传统一阶RC电路模型结构的基础上加入倍率因子的影响,由受控电压源Uocv、串联电阻RΩ、RC网络和极化修正单元四部分构成;
所述受控电压源Uocv代表电池的开路电压OCV随电池SOC的变化;
所述串联电阻RΩ体现了电池在充放电过程中产生的能量损失;
所述RC网络为RpCp并联网络,反映了电池内部由电化学反应不平衡所引起的过电位,通常称为极化电压,记为Up
受控电压源Uocv,串联电阻RΩ,电容Cp与SOC存在非线性关系,具体关系为:
式(1)中,a1、a2、b1、b2、c1、c2为常数,在不同SOC的参数基础上通过线性内插法得到。
在上述技术方案的基础上,所述极化修正单元由Butler-Volmer方程简化形式构成,用于修正电池电流对极化电阻Rp的影响程度,具体关系为:
式(2)中,f1(SOC)、f2(SOC)表示不同SOC下Butler-Volmer方程简化形式的系数,而Rpc代表基准倍率下获得的极化电阻值,Io为通过电池的电流值。
在上述技术方案的基础上,所述OCV-SOC测试数据包括20个充电脉冲和20个放电脉冲,每个脉冲持续时间均为180秒;
为使电池达到平衡状态,在每个脉冲后均会有1小时的静置时间。
在上述技术方案的基础上,所述每个SOC点下典型倍率的充放电极化电压通过典型倍率下的充放电实验和基准倍率下的参数值获得,具体如式(11)所示:
Up(t)=Uocv(t)-Uo(t)-RΩIo(t) (11)
其中Uo代表电池端电压。
在上述技术方案的基础上,所述基于复合形法测算得到的电池特性参数,包括:Butler-Volmer方程简化形式系数的充电部分和放电部分;
Butler-Volmer方程简化形式系数的充电部分讨论的区间为SOC∈[10,75],当SOC∈[75,90]时,采用基准倍率下获得的参数值,
Butler-Volmer方程简化形式系数的放电部分讨论的区间为SOC∈[10,90]。
在上述技术方案的基础上,所述电池最大允许充放电电流和电池功率状态通过电池上限电压和电池下限电压确定。
在上述技术方案的基础上,电池最大允许充放电电流的表达式为:
其中,Uch_lim代表电池充电上限电压,Udch_lim代表电池充电上限电压,Δt表示时间间隔,t3,t4分别对应于电池功率预测过程的起始和结束时刻,Ich_max代表电池最大允许充电电流,,Idch_max代表电池最大允许放电电流;
电池功率状态的表达式为:
其中,Pch_lim和Pdch_lim分别代表电池功率状态的充电部分和放电部分。
本发明所述的轨道交通用钛酸锂电池功率状态预测方法,基于一种电池改进模型,提出一种能够实时准确预测电池功率状态的方法,所预测的电池功率状态不仅精度高,同时可直接用于系统级仿真,可充分利用电池功率能力以提升整车的动力特性,进而为轨道交通用钛酸锂电池运行过程中的最优控制策略提供理论依据。有益效果为:
1.所使用的电池改进模型将电化学Butler-Volmer方程简化形式与一阶RC电路模型相结合,不仅继承了理论描述的优越性,更可直接用于系统级仿真且精度高;
2.首次提出一种适用于钛酸锂电池的动态功率状态预测方法,能够准确刻画钛酸锂电池不同时刻下的功率状态,不仅为电池厂确定电池的限制条件提供了基础,更为钛酸锂电池的广泛应用中的管理策略提供了理论依据。
附图说明
本发明有如下附图:
图1本发明采用模型的结构示意图;
图2钛酸锂电池典型OCV-SOC曲线;
图3复合形法流程图
图4恒电流工况下提出模型仿真电压与实验数据的对比;
图5恒电流工况下仿真电压的误差分布;
图6FUDS(Federal Urban Dynamic Schedule)变电流工况下提出模型仿真电压与实验数据的对比;
图7FUDS变电流工况下仿真电压的误差分布;
图8功率状态预测预测结果;
图9本发明整体流程图。
具体实施方式
以下结合附图对本发明作进一步详细说明。
如图9所示,为本发明的整体流程图,本发明所述的轨道交通用钛酸锂电池功率状态预测方法,基于电池改进模型,通过OCV-SOC测试数据获得基准倍率(本发明中取为1C)下不同SOC的参数值,求取每个SOC点下典型倍率的充放电极化电压,基于复合形法测算电池特性参数,并根据该特性参数集合计算电池最大允许充放电电流和电池功率状态。
在上述技术方案的基础上,如图1所示,所述电池改进模型,在传统一阶RC电路模型结构的基础上加入倍率因子的影响,由受控电压源Uocv、串联电阻RΩ、RC网络和极化修正单元四部分构成;
所述受控电压源Uocv代表电池的开路电压OCV随电池SOC的变化;
所述串联电阻RΩ体现了电池在充放电过程中产生的能量损失;
所述RC网络为RpCp并联网络,反映了电池内部由电化学反应不平衡所引起的过电位,通常称为极化电压,记为Up
受控电压源Uocv,串联电阻RΩ,电容Cp与SOC存在非线性关系,具体关系为:
式(1)中,a1、a2、b1、b2、c1、c2为常数,在不同SOC的参数基础上通过线性内插法得到。
在上述技术方案的基础上,所述极化修正单元由Butler-Volmer方程简化形式构成,用于修正电池电流对极化电阻Rp的影响程度,具体关系为:
式(2)中,f1(SOC)、f2(SOC)表示不同SOC下Butler-Volmer方程简化形式的系数,而Rpc代表基准倍率下获得的极化电阻值,Io为通过电池的电流值。
所述Butler-Volmer方程简化形式由Butler-Volmer方程原始形式经过多次变形所得。Butler-Volmer方程原始形式从电池内部反应的角度解释了电池过电位与电流密度之间的关系,又称为稳态电化学极化方程式,是建立电化学模型的基础公式之一,其表达式为:
其中,R表示摩尔气体常量;T代表开尔文温度;F是法拉第常量;J0表示交换电流密度,α代表传递系数,J表示电流密度,η代表电池过电势,也称为电池极化电压。
然而,式(3)包含多个未知变量且无法从方程中直接得到由电流激励变化引起的过电位值,因此需要对该方程进行形式变换后才能使用。基于式(4)表示的电流与电流密度的关系,式(3)可以化简为式(5),其中S表示有效面积。
Io(t)=J(t)×S (4)
考虑到电化学反应发生过程中一般有α≈0.5,此时,式(5)可以被进一步化简,得到如式(6)所示的关系式。
式(6)中虽然仍包含部分无法外部获得的电化学变量,但可以直接得到由电池电流变化引起的过电位,同时,交换电流密度和有效面积两个变量与电池状态紧密相关,反映了极化过程的复杂程度。可知,其他变量一定时,电池电流的增加会使电池严重偏离平衡态,从而引起更高的过电位。而当电池电流减小并接近零时,过电位也趋于零,表明此时电池接近平衡态。
当电池放电电流突然由I1变为I2时,可认为Δt<ε,其中ε为一个有限小变量,则有式(7)所示不等关系,此时ΔSOC≈0。串联电阻作为电池电极、电解液等电阻特性的集中体现,其数值大小基本与电流无关。由式(6)可得电流变化对电池极化的影响,如式(8)所示,其中Upmax表示电池极化电压建立完全后的值。进而,可以得到Rp与Io的关系式,即本发明中Butler-Volmer方程简化形式。
ΔSOC≤max{I1Δt,I2Δt}≤max{I1ε,I2ε} (7)
在上述技术方案的基础上,所述OCV-SOC测试数据包括20个充电脉冲和20个放电脉冲,每个脉冲持续时间均为180秒;
为使电池达到平衡状态,在每个脉冲后均会有1小时的静置时间。
静置结束时刻的电压值即作为当前SOC下的受控电压源Uocv,图2为钛酸锂电池充放电两种状态下的电池开路电压,可以看出开路电压与SOC密切相关,且钛酸锂电池充电和放电状态下的开路电压具有较高的一致性。考虑到实际采样频率的限制,基于式(9)选取通过1s内的电压跳变获得串联电阻的数值。由于电路结构仅采用一个RC网络表示电池极化过程,因此电压响应采用式(10)进行最小二乘拟合,进而获得不同SOC下的Rpc和Cp的数值。需要注意的是,本发明采用两组参数来描述电池充电与放电过程的差异。
RΩ=ΔUo/ΔIo (9)
其中,ΔUo和ΔIo分别表示1秒内的电压变化和电流变化,而A表示需拟合的极化电压幅值,τ表示极化电压建立过程中的时间常数。
在上述技术方案的基础上,所述每个SOC点下典型倍率的充放电极化电压通过典型倍率下的充放电实验和基准倍率下的参数值获得,具体如式(11)所示:
Up(t)=Uocv(t)-Uo(t)-RΩIo(t) (11)
其中Uo代表电池端电压。
在上述技术方案的基础上,所述复合形法,是求解约束优化问题的一种重要的直接解法,其具体过程如图3所示:
1.选择复合形的顶点数K,给定反射比a,扩张比r,收缩比p和压缩比q,同时给定反射比的极限值δ,将sign清零,并在可行域内构造初始复合形;
2.计算K个顶点的目标函数值,找出最坏点XH,次坏点XG,最好点XL,其中点的好坏程度由其目标函数值的大小确定,最好点为目标函数值最小的点;
3.判断是否满足终止条件,若是,则跳出循环,否则继续下一步;
4.以式(12)计算除最坏点XH外其余各点的形心XC:
5.以式(13)计算反射点XR:
XR=XC+a·(XC-XH) (13)
并检查反射点XR是否在可行域,若是则转向步骤6,否则以a=0.5a代入后,再次以上式进行计算;
6.计算反射点XR的目标函数值,
若f(XR)<f(XH);以XR代替XH重构复合形后转向步骤7,
若f(XR)≥f(XH);转向步骤9,
7.以式(14)计算扩张点XE
XE=XR+r·(XR-XC) (14)
并检查扩张点XE是否在可行域,若是转向步骤8,否则以XR代替XH,并转向步骤2;
8.计算扩张点XE的目标函数值,
若f(XE)<f(XH);以XE代替XH重构复合形后转向步骤2;
若f(XE)≥f(XH);以XR代替XH,并转向步骤2;
9.判断反射比a是否大于δ;
若是,则令a=0.5a后转向步骤5;
若否且sign不为1,则表示之前未用次坏点代替过坏点,因此用次坏点代替过坏点,令sign=1,并转向步骤2;
若否且sign为1,则表示之前已经用次坏点代替过坏点,则转向步骤10;
10.以式(15)计算收缩点XK
XK=XH+p·(XC-XH) (15)
检查收缩点XK是否在可行域,若是转向步骤11,否则转向步骤12;
11.计算收缩点XK的目标函数值,
若f(XK)<f(XH);以XK代替XH重构复合形后转向步骤2;
若f(XK)≥f(XH);转向步骤12;
12.压缩所有顶点并转向步骤2,
Xj=XL-q·(XL-Xj),j=1,2,3,4 and j≠L (16)
由于复合形的形状不必保持规则的图形,对目标函数和约束函数无特殊要求,因此这种方法适应性强,在机械优化设计中应用广泛。复合形法适用于仅含不等式约束的问题计算简单,不必求导,占内存小。
在上述技术方案的基础上,所述基于复合形法测算得到的电池特性参数,包括:Butler-Volmer方程简化形式系数的充电部分和放电部分;
鉴于不同倍率充电过程中,电池从恒流充电转变为恒压充电时处于不同的SOC,且随着倍率的增加电池更早的进入恒压充电阶段,因此本发明的Butler-Volmer方程简化形式系数的充电部分仅讨论SOC∈[10,75]。而当SOC∈[75,90]时,由于恒压阶段电流激励的逐渐降低接近于零,认为电池参数基本不变,仍采用基准倍率下获得的参数值。而Butler-Volmer方程简化形式系数的放电部分讨论的区间为SOC∈[10,90]。
在上述技术方案的基础上,所述电池最大允许充放电电流和电池功率状态通过电池上限电压和电池下限电压确定。
在上述技术方案的基础上,电池最大允许充放电电流的表达式为:
其中,Uch_lim代表电池充电上限电压,Udch_lim代表电池充电上限电压,Δt表示时间间隔,t3,t4分别对应于电池功率预测过程的起始和结束时刻,Ich_max代表电池最大允许充电电流,,Idch_max代表电池最大允许放电电流;
电池功率状态的表达式为:
其中,Pch_lim和Pdch_lim分别代表电池功率状态的充电部分和放电部分。
如上所述技术方案,发明人分别对高倍率恒流工况和FUDS动态工况进行了模型准确性验证,从图4和图5中可以看出,本发明采用的模型的恒流工况仿真结果与实验数据具有高度的一致性,误差在240秒之内小于0.5%,在300秒内误差小于1.5%。从图6和图7可以看出,FUDS变电流工况下,本发明采用的模型仿真结果同样与实验数据具有高度的一致性,误差在2250秒内误差小于1.5%。
在验证模型准确性的基础上,本发明对动态运行条件下的功率型钛酸锂电池进行了功率估计的准确性的验证。从图8可以看出,基于改进模型的充放电功率估计结果均具备较高精度,大部分仿真误差均在10%以内。
由以上对比分析可知,本发明采用的钛酸锂电池模型具备较高精度,能够很好的模拟电池的动静态特性,本发明提出的电池峰值功率的估算方案不仅精度高,同时可直接用于系统级仿真,可充分利用电池的功率能力以提升整车的动力特性,进而为轨道交通用钛酸锂电池运行过程中的最优控制策略提供理论依据。
以上为本发明较佳的具体实现方案,除此之外还有其他实现方案,需要说明的是在没有脱离本发明构思的前提下任何显然意见的替换均在本发明保护范围之内。
本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。

Claims (9)

1.一种轨道交通用钛酸锂电池功率状态预测方法,其特征在于:基于电池改进模型,通过OCV-SOC测试数据获得基准倍率下不同SOC的参数值,求取每个SOC点下典型倍率的充放电极化电压,基于复合形法测算电池特性参数,并根据该特性参数集合计算电池最大允许充放电电流和电池功率状态;
所述复合形法的具体过程如下所示:
1.选择复合形的顶点数K,给定反射比a,扩张比r,收缩比p和压缩系数q,同时给定反射比的极限值δ,将sign清零,并在可行域内构造初始复合形;
2.计算K个顶点的目标函数值,找出最坏点XH,次坏点XG,最好点XL,其中点的好坏程度由其目标函数值的大小确定,最好点为目标函数值最小的点;
3.判断是否满足终止条件,若是,则跳出循环,否则继续下一步;
4.以式(12)计算除最坏点XH外其余各点的形心XC
5.以式(13)计算反射点XR
XR=XC+a*(XC-XH) (13)
并检查反射点XR是否在可行域,若是转向步骤6,否则以a=0.5a代入后,再次以上式进行计算;
6.计算反射点XR的目标函数值,
若f(XR)<f(XH);以XR代替XH重构复合形后转向步骤7,
若f(XR)≥f(XH);转向步骤9,
7.以式(14)计算扩张点XE
XE=XR+r*(XR–XC) (14)
并检查扩张点XE是否在可行域,若是转向步骤8,否则以XR代替XH,并转向步骤2;
8.计算扩张点XE的目标函数值,
若f(XE)<f(XH);以XE代替XH重构复合形后转向步骤2;
若f(XE)≥f(XH);以XR代替XH,并转向步骤2;
9.判断反射比a是否大于δ;
若是,则a=0.5a后转向步骤5;
若否且sign不为1,则表示之前未用次坏点代替过坏点,因此则用次坏点代替过坏点,令sign=1,并转向步骤2;
若否且sign为1,则表示之前已经用次坏点代替过坏点,则转向步骤10;
10.以式(15)计算收缩点XK
XK=XH+p*(XC-XH) (15)
检查收缩点XK是否在可行域,若是转向步骤11,否则转向步骤12;
11.计算收缩点XK的目标函数值,
若f(XK)<f(XH);以XK代替XH重构复合形后转向步骤2;
若f(XK)≥f(XH);转向步骤12;
12.压缩所有顶点并转向步骤2,
Xj=XL-q·(XL-Xj),j=1,2,3,4and j≠L (16);
其中Xj代表复合形的顶点,其中j=1…K;f(XR)、f(XH)、f(XE)、f(XK)代表算法中的目标函数。
2.如权利要求1所述的轨道交通用钛酸锂电池功率状态预测方法,其特征在于:基准倍率取为1C。
3.如权利要求1所述的轨道交通用钛酸锂电池功率状态预测方法,其特征在于:所述电池改进模型,在传统一阶RC电路模型结构的基础上加入倍率因子的影响,由受控电压源Uocv、串联电阻RΩ、RC网络和极化修正单元四部分构成;
所述受控电压源Uocv代表电池的开路电压OCV随电池SOC的变化;
所述串联电阻RΩ体现了电池在充放电过程中产生的能量损失;
所述RC网络为RpCp并联网络,反映了电池内部由电化学反应不平衡所引起的过电位,通常称为极化电压,记为Up
受控电压源Uocv,串联电阻RΩ,电容Cp与SOC存在非线性关系,具体关系为:
式(1)中,a1、a2、b1、b2、c1、c2为常数,在不同SOC的参数基础上通过线性内插法得到。
4.如权利要求3所述的轨道交通用钛酸锂电池功率状态预测方法,其特征在于:所述极化修正单元由Butler-Volmer方程简化形式构成,用于修正电池电流对极化电阻Rp的影响程度,具体关系为:
式(2)中,f1(SOC)、f2(SOC)表示不同SOC下Butler-Volmer方程简化形式的系数,而Rpc代表基准倍率下获得的极化电阻值,Io为通过电池的电流值。
5.如权利要求1所述的轨道交通用钛酸锂电池功率状态预测方法,其特征在于:所述OCV-SOC测试数据包括20个充电脉冲和20个放电脉冲,每个脉冲持续时间均为180秒;
为使电池达到平衡状态,在每个脉冲后均会有1小时的静置时间。
6.如权利要求1所述的轨道交通用钛酸锂电池功率状态预测方法,其特征在于:所述每个SOC点下典型倍率的充放电极化电压通过典型倍率下的充放电实验和基准倍率下的参数值获得,具体如式(11)所示:
Up(t)=Uocv(t)-Uo(t)-RΩIo(t) (11)
其中Uo代表电池端电压;Uocv为受控电压源,RΩ为串联电阻,Io为通过电池的电流值。
7.如权利要求1所述的轨道交通用钛酸锂电池功率状态预测方法,其特征在于:所述基于复合形法测算得到的电池特性参数,包括:Butler-Volmer方程简化形式系数的充电部分和放电部分;
Butler-Volmer方程简化形式系数的充电部分讨论的区间为SOC∈[10,75],当SOC∈[75,90]时,采用基准倍率下获得的参数值,
Butler-Volmer方程简化形式系数的放电部分讨论的区间为SOC∈[10,90]。
8.如权利要求1所述的轨道交通用钛酸锂电池功率状态预测方法,其特征在于:所述电池最大允许充放电电流和电池功率状态通过电池上限电压和电池下限电压确定。
9.如权利要求8所述的轨道交通用钛酸锂电池功率状态预测方法,其特征在于:电池最大允许充放电电流的表达式为:
其中,Uch_lim代表电池充电上限电压,Udch_lim代表电池放电下限电压,Δt表示时间间隔,t3,t4分别对应于电池功率预测过程的起始和结束时刻,Ich_max代表电池最大允许充电电流,Idch_max代表电池最大允许放电电流;RΩ为串联电阻,T代表开尔文温度,τ表示极化电压建立过程中的时间常数,f1(SOC)和f2(SOC)为Butler-Volmer方程简化形式的系数,Uocv为受控电压源,Up为RC网络两端的极化电压,Rpc代表基准倍率下获得的极化电阻值;
电池功率状态的表达式为:
其中,Pch_lim和Pdch_lim分别代表电池功率状态的充电部分和放电部分。
CN201610009188.5A 2016-01-07 2016-01-07 一种轨道交通用钛酸锂电池功率状态预测方法 Active CN105548901B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610009188.5A CN105548901B (zh) 2016-01-07 2016-01-07 一种轨道交通用钛酸锂电池功率状态预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610009188.5A CN105548901B (zh) 2016-01-07 2016-01-07 一种轨道交通用钛酸锂电池功率状态预测方法

Publications (2)

Publication Number Publication Date
CN105548901A CN105548901A (zh) 2016-05-04
CN105548901B true CN105548901B (zh) 2019-03-26

Family

ID=55828214

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610009188.5A Active CN105548901B (zh) 2016-01-07 2016-01-07 一种轨道交通用钛酸锂电池功率状态预测方法

Country Status (1)

Country Link
CN (1) CN105548901B (zh)

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106383324B (zh) * 2016-12-07 2019-07-09 上海动力储能电池系统工程技术有限公司 一种基于容量衰减机理分解分析的锂离子电池寿命预测方法
CN109991553B (zh) * 2017-12-29 2021-01-15 微宏动力系统(湖州)有限公司 一种电池soc的估算方法
CN108226804A (zh) * 2018-01-17 2018-06-29 杭州六创电动汽车科技有限公司 一种电动汽车锂电池sop估算方法
CN108845273B (zh) * 2018-08-30 2020-11-03 北京经纬恒润科技有限公司 一种动力电池功率状态估算功能测试方法和装置
CN109085509A (zh) * 2018-09-11 2018-12-25 山东鲁能智能技术有限公司 锂离子电池开路电压与soc关系的参数辨识方法及系统
CN109738809A (zh) * 2019-01-10 2019-05-10 上海动力储能电池系统工程技术有限公司 动力与储能电池功率特性的估计方法
CN109490790B (zh) * 2019-01-23 2021-05-07 哈尔滨理工大学 采用补偿脉冲法的锂动力电池功率特性测试方法及装置
CN109991547B (zh) * 2019-03-29 2021-06-18 深圳猛犸电动科技有限公司 锂离子电池包最大允许充放电电流估计方法及装置
CN110208700B (zh) * 2019-04-09 2020-07-10 清华大学 一种直流微网中储能系统虚拟电池开路电压的计算方法
CN110244228B (zh) * 2019-07-17 2022-04-15 东软睿驰汽车技术(沈阳)有限公司 一种电池控制方法及装置
CN111177992B (zh) * 2019-12-16 2024-03-29 中车工业研究院有限公司 基于电化学理论和等效电路模型的电池模型及其构建方法
CN111060823A (zh) * 2019-12-24 2020-04-24 南京航空航天大学 一种基于dp模型的低温环境下电池sop在线估计方法
CN111025172B (zh) * 2019-12-31 2022-03-01 国联汽车动力电池研究院有限责任公司 一种实现锂离子电池充放电最大允许功率快速测量的方法
CN111044909B (zh) * 2019-12-31 2022-02-08 北京理工大学 电池的状态预测方法及装置
CN111413630B (zh) * 2020-03-11 2022-06-14 合肥国轩高科动力能源有限公司 一种锂电池的脉冲放电功率修正方法
CN111426962B (zh) * 2020-03-23 2022-08-09 华为技术有限公司 电池参数的辨识方法和终端
CN111532143B (zh) * 2020-04-24 2024-05-31 北京北交新能科技有限公司 一种轨道交通用无网自走行辅助蓄电池储能系统
CN112684347B (zh) * 2020-12-15 2022-08-02 华北电力大学(保定) 一种多重约束条件下的锂电池的可用功率预测方法及系统
CN114167289B (zh) * 2021-12-01 2024-01-16 中国第一汽车股份有限公司 一种电池输出功率的确定方法、装置及车辆
CN116593905B (zh) * 2023-06-30 2023-09-08 广汽埃安新能源汽车股份有限公司 电池功率状态sop计算方法、装置、电子设备和存储介质

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8466684B2 (en) * 2006-06-16 2013-06-18 Chevron Technology Ventures Llc Determination of battery predictive power limits
JP4802945B2 (ja) * 2006-08-31 2011-10-26 トヨタ自動車株式会社 二次電池の制御システムおよびそれを搭載したハイブリッド車両
US8346495B2 (en) * 2010-04-22 2013-01-01 Battelle Energy Alliance, Llc Systems, methods and computer-readable media to model kinetic performance of rechargeable electrochemical devices
CN103077291B (zh) * 2013-01-25 2016-05-18 华北电力大学 可设置初始荷电状态的电池充放电过程数字仿真方法
CN104899439B (zh) * 2015-06-02 2019-05-03 吉林大学 一种锂离子电池机理建模方法

Also Published As

Publication number Publication date
CN105548901A (zh) 2016-05-04

Similar Documents

Publication Publication Date Title
CN105548901B (zh) 一种轨道交通用钛酸锂电池功率状态预测方法
CN106909716B (zh) 计及容量损耗的磷酸铁锂电池建模及soc估计方法
Zhang et al. Battery modelling methods for electric vehicles-A review
Świerczyński et al. Selection and performance-degradation modeling of limo $ _ {2} $/Li $ _ {4} $ Ti $ _ {5} $ O $ _ {12} $ and lifepo $ _ {4} $/C battery cells as suitable energy storage systems for grid integration with wind power plants: an example for the primary frequency regulation service
CN103529399B (zh) 一种基于铅酸电池改进型pngv模型的模拟方法
CN104965179B (zh) 一种锂离子蓄电池的温度组合电路模型及其参数识别方法
CN101226229B (zh) 用于在电化学系统中对扩散建模的方法和装置
CN109586373A (zh) 一种电池充电方法和装置
CN107329094A (zh) 动力电池健康状态估算方法及装置
CN110061531A (zh) 储能电池的均衡方法
CN103983833B (zh) 基于gm(1,1)灰色模型的电池开路电压预测方法
CN105207243A (zh) 一种用于风电场实时功率预报修正的电池能量管理方法
CN105116338A (zh) 一种基于soc补偿器的并联型电池系统建模方法
Tian et al. Lithium-ion battery charging optimization based on electrical, thermal and aging mechanism models
CN113156316A (zh) 盐水电池soc估算算法
CN105677979A (zh) 基于Butler-Volmer方程的钛酸锂电池改进模型
Boi et al. A non-linear electrical model for iron doped sodium metal halides batteries
Li et al. Evaluation and analysis of circuit model for lithium batteries
CN108427080A (zh) 混合动力船舶的动力电池组的荷电状态计算方法
US20230032837A1 (en) Method for determining the value of a parameter related to the state of health of an electrochemical cell in a battery, system for electronic management of a battery, and corresponding battery
Nájera Álvarez Study and Analysis of the Behavior of LFP and NMC Electric Vehicle Batteries concerning their Ageing and their Integration into the Power Grid
Li et al. High-accuracy parameters identification of non-linear electrical model for high-energy lithium-ion capacitor
Celen et al. Equivalent Circuit Modeling of All-Solid-State Battery by using DC-IR Data
Raman et al. Computationally efficient and accurate modeling of Li-ion battery
Nejad Adaptive techniques for estimation and online monitoring of battery energy storage devices

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant