CN113884901A - 一种电池表面温度分布估计方法及系统 - Google Patents
一种电池表面温度分布估计方法及系统 Download PDFInfo
- Publication number
- CN113884901A CN113884901A CN202111183859.7A CN202111183859A CN113884901A CN 113884901 A CN113884901 A CN 113884901A CN 202111183859 A CN202111183859 A CN 202111183859A CN 113884901 A CN113884901 A CN 113884901A
- Authority
- CN
- China
- Prior art keywords
- battery
- heat
- temperature
- model
- resistance
- 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
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000009826 distribution Methods 0.000 title claims abstract description 52
- 238000012546 transfer Methods 0.000 claims abstract description 44
- 238000004364 calculation method Methods 0.000 claims abstract description 19
- 230000001052 transient effect Effects 0.000 claims abstract description 17
- 238000012360 testing method Methods 0.000 claims description 33
- 230000010287 polarization Effects 0.000 claims description 16
- 238000007599 discharging Methods 0.000 claims description 12
- 230000020169 heat generation Effects 0.000 claims description 11
- 238000007600 charging Methods 0.000 claims description 10
- 239000003990 capacitor Substances 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 9
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 230000002441 reversible effect Effects 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 230000007613 environmental effect Effects 0.000 claims description 3
- 230000008878 coupling Effects 0.000 abstract description 14
- 238000010168 coupling process Methods 0.000 abstract description 14
- 238000005859 coupling reaction Methods 0.000 abstract description 14
- 238000013461 design Methods 0.000 abstract description 6
- 238000012544 monitoring process Methods 0.000 abstract description 4
- 230000004043 responsiveness Effects 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 abstract description 2
- 238000012795 verification Methods 0.000 description 10
- HBBGRARXTFLTSG-UHFFFAOYSA-N Lithium ion Chemical compound [Li+] HBBGRARXTFLTSG-UHFFFAOYSA-N 0.000 description 4
- 125000004122 cyclic group Chemical group 0.000 description 4
- 229910001416 lithium ion Inorganic materials 0.000 description 4
- 238000010278 pulse charging Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical group [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000001808 coupling effect Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 229910052744 lithium Inorganic materials 0.000 description 2
- 230000026676 system process Effects 0.000 description 2
- 239000002253 acid Substances 0.000 description 1
- OJIJEKBXJYRIBZ-UHFFFAOYSA-N cadmium nickel Chemical compound [Ni].[Cd] OJIJEKBXJYRIBZ-UHFFFAOYSA-N 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 210000005069 ears Anatomy 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003487 electrochemical reaction Methods 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000000178 monomer Substances 0.000 description 1
Images
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/374—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC] with means for correcting the measurement for temperature or ageing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Secondary Cells (AREA)
Abstract
本发明公开了一种电池表面温度分布估计方法及系统，电池温度估计方法首次提出三热源电‑热耦合模型，包括三热源传热模型和电阻分布式等效电流模型，分别用以计算电池的热特性以及电特性，以实现特征温度点的瞬时温度估计；具体包括步骤：建立三热源传热模型以及电阻分布式等效电路模型；不同工况下测量并获取电池电压值、电流值以及特征点的温度值；分别对模型中相关参数进行辨识；在稳态和瞬态工况下验证模型的精度和响应性；进一步提高了电池管理系统的温度实时监测能力。本发明所述方法中仅选取单体电池中四个特征温度点进行分析，大幅降低计算复杂度，易于实际应用，对电池表面温度分布实时估计及热管理结构设计具有重要意义。
Description
技术领域
本发明属于电池管理技术领域，涉及一种电池表面温度分布估计方法及系统。
背景技术
锂离子电池作为储能装置被广泛地运用于电动汽车及移动机器人中，其性能受到温度影响较大。随着电池本体技术的发展，锂离子电池能量密度不断提高，热安全问题越来越值得重视。由于电池内正负两极集流体的汇流作用，锂离子电池的产热不均导致温度分布不均，在一些工况下，单体电池上的温度差异甚至超过10℃，极大影响了电池的安全性以及使用寿命。因此，建立一个能快速估算电池表面温度分布的热模型显得尤为重要。
在电池温度分布研究中，由于电池是一个复杂的多物理场耦合系统，不同因素之间相互关联，内部状态不断变化且难以测量，这给电池温度估计的实时性和准确性增加了难度。
当前电池热模型研究可主要分为两类，电-热耦合模型及电化学-热耦合模型，包含电化学特性子模型、电特性子模型以及传热子模型。电化学-热耦合模型主要是从电化学反应层面上描述电池工作特性，基于能量守恒进而得到电池的产热，结合传热模型得到温度分布规律。该模型有较高的精度，但由于其模型复杂度及参数辨识的难度较大而不适合用于实时的温度估计；电-热耦合模型通过对电池的宏观特性等进行分析(如：电压、电流、电阻、荷电状态(SOC)等)计算产热量，该模型计算速度快及准确性较好，但由于其一般通过总内阻计算热量，往往只适用于整体温度估计而无法考虑温度的分布规律。在传热子模型方面，以三维有限元模型以及一维集总参数热模型为代表，其优缺点较为明显，三维模型能较好反应出电池温度分布规律但计算量大，一维模型计算速度快但难以反应温度的空间分布规律。
现有方法中未能实现高精度、速度快的单体电池温度分布规律估计，而难以对电池在瞬态工况下的温度特性进行在线评估。因此，如何同时在小计算量的情况下，对于电池在实际使用中的瞬态温度分布进行精确预估，也是本领域技术人员亟待解决的重点问题。
发明内容
为了克服上述现有技术的缺点，本发明的目的在于提供一种电池表面温度分布估计方法及系统，旨在解决现有技术中不能实现单体电池温度分布的快速估计，计算复杂度高，不利于实际应用的缺陷性技术问题。
为了达到上述目的，本发明采用以下技术方案予以实现：
本发明公开了一种单体电池温度分布的快速估计方法，包括以下步骤：
1)在脉冲循环充放电工况下测量并获取电池特征点的温度值；在不同温度下进行混合脉冲功率特性HPPC试验并获取电压值和电流值；
2)根据步骤1)获得的温度值、电压值和温度值建立三热源传热模型以及电阻分布式等效电路模型；
3)根据最小二乘算法离线辨识步骤2)建立的三热源传热模型和电阻分布式等效电路模型中的未知参数；
4)根据步骤3)得到离线辨识的未知参数计算电池特征点的瞬时温度，实现电池表面温度的分布估计。
优选地，在步骤2)中，电阻分布式等效电路模型的参数包括总欧姆内阻R0、极化内阻R1以及极化电容C1；电阻分布式等效电路模型在一阶RC等效电路基础上将总欧姆内阻R0等效为不同分布式内阻Re1、Re2、Rp和Rn，将不同分布式内阻通过串联和并联结合对应不同热源。优选地，在步骤2)中，三热源传热模型包括三个独立热源Q1、Q2和Q3，三个独立热源对应的热容为Ct1、Ct2和Ct3，以及三个独立热源对应的传热热阻Rx1、Rx2、Ry1、Ry2和Rh。
优选地，在步骤3)中，包括如下步骤：
a)、根据正、负极耳附近温度T2和T3完成对两极耳处热源Q2和Q3及其关联热阻Rx1和Rx2的估计；
b)、根据混合脉冲功率特性HPPC试验数据估计电池总欧姆内阻R0、极化内阻R1和电容C1；
c)、根据估计的正、负极耳处热源Q2和Q3推算欧姆内阻分量Re1和Re2，计算得到中间热源Q1的值；
d)、根据电池本体中间温度T1以及下部温度T4估计相应部位的传热热阻Ry1、Ry2和Rh参数，其中Rh为电池与空气之间的传热热阻。
优选地，在步骤4)中，包括如下步骤：
a)、在电池温度的瞬态计算中，根据初始温度、电流和电池荷电状态SOC，计算电池总欧姆内阻R0以及各欧姆内阻分量Re1和Re2；
b)、根据生热公式计算三个热源Q1、Q2和Q3的瞬时生热量，再根据三热源传热模型计算四个特征点瞬时温度T1、T2、T3和T4；
其中，中间本体热源Q1包括欧姆热、可逆热和极化热，正、负极耳两处热源Q2和Q3以欧姆热形式体现；
c)、完成一个时间步长的温度计算后，更新电池荷电状态SOC，并转入更新测量时刻t＝t+1，重新执行步骤b)再次计算瞬时生热量以及瞬时温度，循环至预设时间结束或电池荷电状态SOC达到限值；
其中，电池荷电状态的计算方法如下：
I表示为电流；Qc表示为标称电荷容量。
优选地，电池表面温度分布估计过程中，在恒流放电工况和混合道路循环工况下验证三热源传热模型及电阻分布式等效电路模型在瞬态温度的精度和工况适应性。
优选地，在脉冲循环脉冲充放电过程中，充放电过程中产生的可逆热相互抵消，能够使电池热源保持稳定利于参数辨识。
优选地，通过拟合不同环境温度下的混合脉冲功率特性HPPC测试数据，获取总欧姆内阻随温度及电池荷电状态SOC的变化规律；
通过拟合熵热系数，得到电池产生的可逆热随电池荷电状态SOC的变化规律。
优选地，监控步骤b)中四个特征温度点的实时温度，能够通过插值估算得到电池单体的温度分布。
本发明还提出了一种电池表面温度分布估计方法的估计系统，包括：
实际参数获取模块，用于在脉冲循环充放电工况下测量并获取电池特征点的温度值，在不同温度下进行混合脉冲功率特性HPPC试验并获取电压值和电流值；
模型建立模块，用于建立三热源传热模型以及电阻分布式等效电路模型；
模型参数获取模块，用于获取三热源传热模型以及电阻分布式等效电路模型的未知参数，再根据未知参数计算电池特征点的瞬时温度。
与现有技术相比，本发明具有以下有益效果：
1)、本发明公开了一种电池表面温度分布估计方法，基于三热源电-热耦合模型估计单体电池温度分布，该方法根据实际应用需要，首次提出三热源传热模型和电阻分布式等效电路模型，对三热源传热模型进行创新设计，通过将单体电池的不均匀产热等效为三个独立热源，实现迅速、准确地单体电池特征点温度瞬时估计。本发明提供的上述方法能够实现高精度的瞬时温度分布估计，提高了电池热管理系统的实时温度估计能力。
2)、本发明简化算法设计，通过采集单体电池特征点在脉冲循环充放电工况和不同温度下进行混合脉冲功率特性HPPC试验，获得电池单体的温度值、电压值和电流值，在不同测试条件下分别辨识得到三热源传热模型的未知参数和电阻分布式等效电路模型的未知参数，避免电-热耦合作用带来的参数辨识难度，简化了温度分布计算量，大幅降低计算复杂度，易于实际应用。
进一步地，本发明的算法步骤设计合理，在对不同特征点的瞬时温度计算中，直接利用不同热源以及热阻进行计算，避免了现有技术中采用有限元计算划分大量网格获得电池温度分布，大大减少了计算量。
本发明还提出了一种电池表面温度分布估计方法的估计系统，采用模块化思想来实现初至拾取，使各个模块之间相互独立，便于实现对各个模块的统一管理。
附图说明
图1为本发明所提出的方法进行电池温度估计的流程图；
图2为电池三热源电-热耦合模型图；其中，(a)为三热源传热模型，(b)为电阻分布式电模型；
图3为本实施例中25℃下HPPC试验电压特性曲线；
图4为本实施例中2C倍率下循环脉冲充放电下电特性；
图5为本实施例中1C和3C倍率下参数辨识结果验证；其中，(a)为1C倍率实验结果验证，(b)为3C倍率实验结果验证；
图6为本实施例中恒流放电温度估计结果；其中，(a)为1C倍率温度验证，(b)为2C倍率温度验证；
图7为本实施例中混合道路循环工况下模型验证结果；其中，(a)为电压验证，(b)为温度验证。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
下面结合附图对本发明做进一步详细描述：
参见图1所示流程图，为本发明提供的一种电池表面温度分布估计方法，包括以下步骤：
1)、在脉冲循环充放电工况以及不同温度下进行混合脉冲功率特性HPPC试验下测量并获取单体电池特征点的温度值、电池电压和电流测量值，电池特征点的温度值是四个特征点表面温度；
2)、根据步骤1)获得的温度值、电压值和温度值建立建立三热源传热模型以及分布式一阶RC等效电路模型；
在步骤2)中，电阻分布式等效电路模型的参数包括总欧姆内阻R0、极化内阻R1以及极化电容C1；电阻分布式等效电路模型在一阶RC等效电路基础上将总欧姆内阻R0等效为不同分布式内阻Re1、Re2、Rp和Rn，将不同分布式内阻通过串联和并联结合对应不同热源；
在步骤2)中，三热源传热模型包括三个独立热源Q1、Q2和Q3，三个独立热源对应的热容为Ct1、Ct2和Ct3，以及三个独立热源对应的传热热阻Rx1、Rx2、Ry1、Ry2和Rh；正负两极耳处分别对应两个热源Q2和Q3，电池本体处对应一个热源Q1；
3)、根据最小二乘算法离线辨识步骤2)建立的电阻分布式等效电路模型和三热源传热模型中的未知参数包括以下步骤：
a)、根据正、负极耳附近温度T2和T3完成对两极耳处热源Q2和Q3及其关联热阻Rx1和Rx2的估计；
b)、根据混合脉冲功率特性HPPC试验数据估计电池总欧姆内阻R0、极化内阻R1和电容C1；
c)、根据估计的正、负极耳处热源Q2和Q3推算欧姆内阻分量Re1和Re2，计算得到中间热源Q1的值；
d)、根据电池本体中间温度T1以及下部温度T4估计相应部位的传热热阻Ry1、Ry2和Rh参数，其中Rh为电池与空气之间的传热热阻。
具体的：在上述步骤中，辨识电池总欧姆内阻R0：
其中ΔU0为瞬间电压变化，I为测试电流；
根据电池一阶RC等效电路模型推导得出电压满足以下方程：
根据最小二乘法辨识电池极化电压和极化内阻：
其中Uexp(t)为试验值，U(t)为拟合值，J表示为电压方差的最小值。辨识热源Q2,Q3以及热阻Rx1+Rh,Rx2+Rh：
其中Ct2,Ct3为对应热容Tamb表示为环境温度；Rx1,2表示为电池内部传热热阻；Rh表示为电池对环境的散热热阻；Q2,3表示为正负极处对应热源；Tc表示为电池初始温度；
根据总欧姆内阻R0、辨识热源Q2和Q3辨识电路模型中内阻分量Re1和Re2：
其中，RP表示为正极处内阻；Rn表示为负极处内阻；i表示为总电流；
辨识热源Q1以及热阻Ry1和Ry2：
辨识热阻Rh：
优选地，步骤3)中，参数辨识过程采用稳态工况测试数据，包括不同环境温度下的混合脉冲功率特性(HPPC)测试以及0.5SOC下的循环脉冲充放电测试：
1)HPPC测试：电池充满电情况下，以1C放电6分钟，搁置30分钟，再执行20秒的脉冲充放电一次，搁置10分钟；循环以上步骤直到SOC降至0.1，并根据电压曲线变化计算电流中内阻值以及电容值。在不同环境温度下重复以上步骤，完成不同温度下的等效电路内阻以及电容估计；
2)脉冲循环充放电测试：当SOC为0.5时，执行充放电各10秒，循环100次，；在1C、2C、3C倍率下重复以上步骤，根据电池各区域的不同温度完成对热阻的估计。
优选地，步骤3)中，参数辨识算法为最小二乘参数辨识算法。
优选地，步骤3)中，参数辨识主要分为两部分：根据正负极耳附近温度T2和T3完成对两极耳处热源Q2和Q3及其关联热阻Rx1和Rx2的估计；根据HPPC试验数据估计电池欧姆内阻R0、极化内阻R1和电容C1；根据估计的极耳处热源推算欧姆内阻分量Re1和Re2，计算得到中间热源Q1的值；根据电池本体中间温度T1以及下部温度T4估计相应部位的传热热阻Ry1、Ry2和Rh参数，其中Rh为电池与空气之间的传热热阻。
4)、根据步骤3)得到的未知参数计算电池特征点的瞬时温度，通过插值得到整体的温度分布，利用恒流放电工况以及混合道路循环工况数据验证三热源传热模型及电阻分布式等效电路模型精度以及工况适应性，实现电池表面温度的分布估计及验证。
包括以下步骤：
a)、在电池温度的瞬态计算中，根据初始温度、电流和电池荷电状态SOC，计算电池总欧姆内阻R0以及各欧姆内阻分量Re1和Re2；
b)、根据生热公式计算三个热源Q1、Q2和Q3的瞬时生热量，再根据三热源传热模型计算四个特征点瞬时温度T1、T2、T3和T4；
其中，中间本体热源Q1包括欧姆热、可逆热和极化热，正、负极耳两处热源Q2和Q3以欧姆热形式体现；
c)、完成一个时间步长的温度计算后，更新电池荷电状态SOC，并转入更新测量时刻t＝t+1，重新执行步骤b)再次计算瞬时生热量以及瞬时温度，循环至预设时间结束或电池荷电状态SOC达到限值；
其中，电池荷电状态的计算方法如下：
I表示为电流；Qc表示为标称电荷容量。
优选地，步骤4)中，混合道路循环工况由两种典型城市道路工况循环(UDDS、NEDC)组成，先进行2C倍率充电至截止电压，之后转入三个UDDS(Urban Dynamometer DrivingSchedule)工况测试，再执行2C倍率充电至截止电压，最后完成三个NEDC(New EuropeanDriving Cycle)工况测试。
优选地，步骤4)中，一般只需监控四个特征温度点的实时温度，即可通过插值估算得到电池单体的温度分布。
初始化t-1时刻的参数如电流I(t-1)、温度T(t-1)、电荷状态SOC(t-1)；
具体的，在以上步骤中：
i)、在电池瞬态温度计算之前，根据经验公式对不同温度及SOC下的电池内阻进行拟合。根据经验公式计算电池总内阻以及分阻值在当前温度以及SOC下的熵热系数，并利用多项式拟合计算：
ii)、根据拟合公式计算熵热系数：
其中b0,b1,b2,b3,b4为拟合系数；
计算瞬态热源Q2(t)和Q3(t)：
计算瞬态温度T2和T3：
计算瞬态热源Q1(t)：
计算瞬态温度T1和T4：
iii)预设时间及SOC判断：若满足t＜tn或0＜SOC＜1，其中tn为预设时间，则SOC更新：
时间更新：t＝t+1，并返回第ii)步，继续循环直到系统过程结束；否则，系统过程结束；
所述预设时间tn为根据实际运行工况提前预设的量；
所述电池为锂电池、镍镉电池、镍氢电池或铅酸电池等。
其中，电池特征点温度估计是基于三热源电热耦合模型实现的，如图2所示。其中电模型包括由开路电压源Voc及电池欧姆内阻R0、极化内阻R1、极化电容C1组成的Rint模型部分，以及分布式电阻Re1和Re2。测量电流记为I，测量端电压为V0；热模型包括三个独立热源Q1、Q2和Q3，对应的热容为Ct1、Ct2和Ct3，以及传热热阻Rx1、Rx2、Ry1、Ry2和Rh；
参见图2可知，三个热源Q1、Q2和Q3与电池表面温度T1、T2、T3和T4存在关系，温度估计表达式如下所示：
根据预设时间步长，将方程在时间域进行离散，进而进行迭代更新各点的温度，其离散化温度估计表达式如下所示：
实施例：
本实施例采用的电池为三元锂软包电池，标称容量为10Ah，上限电压为4.2V，下限电压为2.75V。为模拟单体电池在不同工况下的适应性，采取稳态工况以及瞬态工况以验证所述方法。
对电模型中电阻及电容进行参数辨识，分别在25℃、30℃、35℃、40℃下进行混合脉冲功率特性(HPPC,Hybrid Pulse Power Characteristic)测试，25℃下电池电压特性由图3所示；对热模型中热源及热阻进行参数辨识，以2C倍率脉冲循环充放电试验的温度作为作为数据进行辨识，每个循环充电放电各10s，共循环100次，由图4可见，之后对辨识出的数据代入1C和3C倍率下试验数据进行验证。不同倍率下循环脉冲充放电的估计结果参见图5。由图5可见，试验值与计算值吻合程度较好，验证了参数辨识结果的稳定性；1C倍率下的温度波动要大于3C倍率，当t>1000s时，整体温度趋于稳定，说明了在该工况下内阻的稳定性。
为验证在内阻随温度以及SOC瞬时变化的工况下对表面温度的估计，进行了不同倍率下恒流放电工况以及混合道路循环工况下的实验，在这两种工况下，估计时间小于1s，从而说明了该方法用于单体温度估计的快速性。
1C和2C恒流放电条件下各点温度估计结果参见图6。由图6可见，最高温度出现在靠近正极的测点1，而最低温度出现在电池底部的测点4，从试验值和计算值之间的对比来看，1C倍率下最大误差为0.9℃而2C倍率下温度最大误差为1.4℃，估计温度值可以实现良好的精度；同时可得1C倍率下最大温差为3.6℃，而在1C倍率下，最大温差则达到了6.8℃。
为验证瞬态情况下电池温度估计的响应性以及精度，在混合道路循环工况下进行测试，电压和温度估计结果参见图7。由图7可见，测点1和测点4温度试验值与计算值吻合性较好，过程中的最大偏差不超过1℃，均方根误差不超过0.3℃，且可以在1s内计算完成四个测点的温度，验证了长时间瞬态工况下温度估计的准确性以及响应性。
本发明从实际安全角度出发，研究了单体电池在不同工况下的温度分布，在短时间的情况下得到电池特征区域的温度，提高了目前温度估计的运算效率，利于扩展到实际的实时监控。
本发明提出一种单体电池三热源热模型以及电阻分布式等效电路模型，利用模型耦合对特征温度点进行估计，属于模型的创新设计；通过电模型与热模型之间的耦合作用，可实现在瞬时工况下单体电池温度分布的实时估计，应用到实际中也可以对电池单体热不均匀性进行定量诊断。此外，本申请通过对算法的创新设计，选取了脉冲循环充放电试验数据和HPPC试验数据，先分别对热模型以及电模型中的参数进行辨识，之后再结合起来即可实现对整个电-热耦合模型的参数辨识，而不需要因为耦合作用同时对两个模型中的所有参数进行相关的辨识，从而有效降低辨识难度和计算复杂度。此外，现有温度分布估计方法需要对多个节点进行计算，而本发明中则直接利用四个特征点的温度进行插值计算即可获得单体电池的温度分布。
本发明公开了一种电池表面温度分布估计方法，基于三热源电-热耦合模型估计单体电池表面温度分布的方法，该方法包括步骤：建立三热源传热模型以及电阻分布式电模型；不同工况下测量并获取电池特征点的温度值以及电压值和电流值；分别对三热源传热模型及电阻分布式等效电路模型中相关参数进行辨识；在稳态和瞬态工况下验证模型的精度和响应性。本发明根据实际应用需要，首次提出三热源电-热耦合模型，实现对单体电池温度分布的快速估计，过程中仅选取单体电池中四个温度特征点进行分析，大幅降低计算复杂度，易于实际应用，对电池组热管理设计及在线监测具有重要意义。
本发明还提出了一种电池表面温度分布估计方法的估计系统，包括：
实际参数获取模块，用于在脉冲循环充放电工况下测量并获取电池特征点的温度值，在不同温度下进行混合脉冲功率特性HPPC试验并获取电压值和电流值；
模型建立模块，用于建立三热源传热模型以及电阻分布式等效电路模型；
模型参数获取模块，用于获取三热源传热模型以及电阻分布式等效电路模型的未知参数，再根据未知参数计算电池特征点的瞬时温度。
以上内容仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明权利要求书的保护范围之内。
Claims (10)
1.一种电池表面温度分布估计方法，其特征在于，包括以下步骤：
1)在脉冲循环充放电工况下测量并获取电池特征点的温度值；在不同温度下进行混合脉冲功率特性HPPC试验并获取电压值和电流值；
2)根据步骤1)获得的温度值、电压值和温度值建立三热源传热模型以及电阻分布式等效电路模型；
3)根据最小二乘算法离线辨识步骤2)建立的三热源传热模型和电阻分布式等效电路模型中的未知参数；
4)根据步骤3)得到离线辨识的未知参数计算电池特征点的瞬时温度，实现电池表面温度的分布估计。
2.根据权利要求1所述的电池表面温度分布估计方法，其特征在于，在步骤2)中，电阻分布式等效电路模型的参数包括总欧姆内阻R0、极化内阻R1以及极化电容C1；电阻分布式等效电路模型在一阶RC等效电路基础上将总欧姆内阻R0等效为不同分布式内阻Re1、Re2、Rp和Rn，将不同分布式内阻通过串联和并联结合对应不同热源。
3.根据权利要求2所述的电池表面温度分布估计方法，其特征在于，在步骤2)中，三热源传热模型包括三个独立热源Q1、Q2和Q3，三个独立热源对应的热容为Ct1、Ct2和Ct3，以及三个独立热源对应的传热热阻Rx1、Rx2、Ry1、Ry2和Rh。
4.根据权利要求3所述的电池表面温度分布估计方法，其特征在于，在步骤3)中，包括如下步骤：
a)、根据正、负极耳附近温度T2和T3完成对两极耳处热源Q2和Q3及其关联热阻Rx1和Rx2的估计；
b)、根据混合脉冲功率特性HPPC试验数据估计电池总欧姆内阻R0、极化内阻R1和电容C1；
c)、根据估计的正、负极耳处热源Q2和Q3推算欧姆内阻分量Re1和Re2，计算得到中间热源Q1的值；
d)、根据电池本体中间温度T1以及下部温度T4估计相应部位的传热热阻Ry1、Ry2和Rh参数，其中Rh为电池与空气之间的传热热阻。
5.根据权利要求4所述的电池表面温度分布估计方法，其特征在于，在步骤4)中，包括如下步骤：
a)、在电池温度的瞬态计算中，根据初始温度、电流和电池荷电状态SOC，计算电池总欧姆内阻R0以及各欧姆内阻分量Re1和Re2；
b)、根据生热公式计算三个热源Q1、Q2和Q3的瞬时生热量，再根据三热源传热模型计算四个特征点瞬时温度T1、T2、T3和T4；
其中，中间本体热源Q1包括欧姆热、可逆热和极化热，正、负极耳两处热源Q2和Q3以欧姆热形式体现；
c)、完成一个时间步长的温度计算后，更新电池荷电状态SOC，并转入更新测量时刻t＝t+1，重新执行步骤b)再次计算瞬时生热量以及瞬时温度，循环至预设时间结束或电池荷电状态SOC达到限值；
其中，电池荷电状态的计算方法如下：
I表示为电流；Qc表示为标称电荷容量。
6.根据权利要求5所述的电池表面温度分布估计方法，其特征在于，电池表面温度分布估计过程中，在恒流放电工况和混合道路循环工况下验证三热源传热模型及电阻分布式等效电路模型在瞬态温度的精度和工况适应性。
7.根据权利要求5所述的电池表面温度分布估计方法，其特征在于，在脉冲循环充放电过程中，充放电过程中产生的可逆热相互抵消，能够使电池产热保持稳定。
8.根据权利要求5所述的电池表面温度分布估计方法，其特征在于，通过拟合不同环境温度下的混合脉冲功率特性HPPC测试数据，获取总欧姆内阻随温度及电池荷电状态SOC的变化规律；
通过拟合熵热系数，得到电池产生的可逆热随电池荷电状态SOC的变化规律。
9.根据权利要求5所述的电池表面温度分布估计方法，其特征在于，监控步骤b)中四个特征温度点的实时温度，能够通过插值估算得到电池单体的温度分布。
10.采用权利要求1～9任意一项所述的电池表面温度分布估计方法的估计系统，其特征在于，包括：
实际参数获取模块，用于在脉冲循环充放电工况下测量并获取电池特征点的温度值，在不同温度下进行混合脉冲功率特性HPPC试验并获取电压值和电流值；
模型建立模块，用于建立三热源传热模型以及电阻分布式等效电路模型；
模型参数获取模块，用于获取三热源传热模型以及电阻分布式等效电路模型的未知参数，再根据未知参数计算电池特征点的瞬时温度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111183859.7A CN113884901A (zh) | 2021-10-11 | 2021-10-11 | 一种电池表面温度分布估计方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111183859.7A CN113884901A (zh) | 2021-10-11 | 2021-10-11 | 一种电池表面温度分布估计方法及系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113884901A true CN113884901A (zh) | 2022-01-04 |
Family
ID=79006173
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111183859.7A Pending CN113884901A (zh) | 2021-10-11 | 2021-10-11 | 一种电池表面温度分布估计方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113884901A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116432384A (zh) * | 2023-02-14 | 2023-07-14 | 华中科技大学 | 一种二元电流引线电磁热耦合求解方法和系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110161423A (zh) * | 2019-06-26 | 2019-08-23 | 重庆大学 | 一种基于多维度耦合模型的动力锂电池状态联合估计方法 |
| CN110795851A (zh) * | 2019-10-30 | 2020-02-14 | 西安理工大学 | 一种考虑环境温度影响的锂离子电池建模方法 |
| CN111929581A (zh) * | 2020-06-05 | 2020-11-13 | 西安理工大学 | 一种动力锂电池内外部温度预测方法 |
-
2021
- 2021-10-11 CN CN202111183859.7A patent/CN113884901A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110161423A (zh) * | 2019-06-26 | 2019-08-23 | 重庆大学 | 一种基于多维度耦合模型的动力锂电池状态联合估计方法 |
| CN110795851A (zh) * | 2019-10-30 | 2020-02-14 | 西安理工大学 | 一种考虑环境温度影响的锂离子电池建模方法 |
| CN111929581A (zh) * | 2020-06-05 | 2020-11-13 | 西安理工大学 | 一种动力锂电池内外部温度预测方法 |
Non-Patent Citations (3)
| Title |
|---|
| LIANG GUO等: "Analysis on internal defects of electrical equipments in substation analysis on internal defects of electrical equipments in substation using heating simulation for infrared diagnose using heating simulation for infrared diagnose", 《2014 CHINA INTERNATIONAL CONFERENCE ON ELECTRICITY DISTRIBUTION (CICED)》 * |
| ZHECHEN GUO等: "A Three-Heat-Source Electro-Thermal Coupled Model for Fast Estimation of the Temperature Distribution of a Lithium-Ion Battery Cell", 《IEEE TRANSACTIONS ON TRANSPORTATION ELECTRIFICATION》 * |
| 姜水生等: "锂离子电池电-热耦合模型分析及其温度场仿真研究", 《热科学与技术》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116432384A (zh) * | 2023-02-14 | 2023-07-14 | 华中科技大学 | 一种二元电流引线电磁热耦合求解方法和系统 |
| CN116432384B (zh) * | 2023-02-14 | 2024-02-27 | 华中科技大学 | 一种二元电流引线电磁热耦合求解方法和系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106909716B (zh) | 计及容量损耗的磷酸铁锂电池建模及soc估计方法 | |
| Li et al. | Energy efficiency of lithium-ion battery used as energy storage devices in micro-grid | |
| CN105954679B (zh) | 一种锂电池荷电状态的在线估计方法 | |
| Kang et al. | Comparison of comprehensive properties of Ni-MH (nickel-metal hydride) and Li-ion (lithium-ion) batteries in terms of energy efficiency | |
| CN104849675A (zh) | 锂离子电池电化学和热耦合模型的获取方法 | |
| CN110795851A (zh) | 一种考虑环境温度影响的锂离子电池建模方法 | |
| CN111537885B (zh) | 一种串联电池组的多时间尺度短路电阻估计方法 | |
| CN111062137A (zh) | 一种锂离子电池性能预测模型及其构建方法和应用 | |
| CN108508370A (zh) | 一种基于温度校正的开路电压-安时积分soc估计方法 | |
| CN109752660B (zh) | 一种无电流传感器的电池荷电状态估计方法 | |
| CN111610450A (zh) | 一种锂离子电池充电产热来源的估算方法 | |
| CN110442901A (zh) | 一种锂离子电池电化学简化模型及其参数的获取方法 | |
| Seger et al. | Open data model parameterization of a second-life Li-ion battery | |
| CN112068000A (zh) | 一种考虑动力电池耐久性影响的峰值功率预测方法 | |
| CN111965544A (zh) | 基于电压及电流双约束的车用并联动力电池的最小包络线soc估计方法 | |
| Kim et al. | Hysteresis modeling for model-based condition monitoring of lithium-ion batteries | |
| CN113325327B (zh) | 基于内阻测试的动力电池瞬态产热率的测算方法 | |
| Bellache et al. | Aging characterization of lithium iron phosphate batteries considering temperature and direct current undulations as degrading factors | |
| CN113515860B (zh) | 一种锂离子电池组运行状态的数值模拟方法及系统 | |
| CN113884901A (zh) | 一种电池表面温度分布估计方法及系统 | |
| CN114089193A (zh) | 电池的温度和负极电位在线估计方法、装置和计算机设备 | |
| Deng et al. | An improved algorithm of SOC testing based on open-circuit voltage-ampere hour method | |
| CN110031772B (zh) | 一种锂离子电池等效内阻的实时估算方法 | |
| CN112684342A (zh) | 一种钠离子电池充/放电过程中产热率的评估方法 | |
| CN114184958B (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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20220104 |
|
| RJ01 | Rejection of invention patent application after publication |