一种基于气液动力学电池模型在线估算开路电压的方法及
系统
技术领域
[0001] 本
发明涉及动力电池开路电压在线估算领域，特别是涉及一种基于气液动力学电池模型在线估算开路电压的方法及系统。
背景技术
[0002] 随着国家经济的飞速发展和社会不断进步，大众物质生活
水平得到不断提升，与此同时作为人们出行的主要交通工具——
汽车快速发展，全球的石油存储量日益减少，因
此世界各国都在大力鼓励发展新
能源汽车，尤其是电动汽车。动力电池作为混合动力和电
动汽车的核屯、动力部件，直接决定着新能源汽车的整车性能，不仅制约其由科研到产业
化、规模化发展，而且决定着新能源汽车的市场价格。
电池管理系统(BMS)是控
制动力电池的核心，其主要功能是依靠对电池状态变量例如端电压，
电流，
温度的采集，实现对电池充放电，电池状态估算，电池
热管理等多方面的协调控制。一个高效的电池管理系统能够保证
电动车充放电安全，提高车辆的行驶里程，延长电池包的使用寿命，降低电池全生命周期使用成本等。电池管理系统中的一个关键变量就是动力电池的开路电压，开路电压估算的准
确性将直接决定电池状态估计的
精度。
[0003] 现有的开路电压估算模型主要有等效
电路模型和电化学模型。等效电路模型以最简单的一阶RC电路模型为例，虽然具有简单直观的特点，但是其估算精度较低，而高阶RC电路模型尽管估算的精确度有一定的提高，但是相应的计算复杂度以及计算时间也显著提
高，很难在车辆上实现对开路电压的在线估计。电化学模型主要从微观的
角度上反映动力
电池内部的化学机理，同时运用相当复杂的偏微分方程，描述微观信息(离子的浓度，分布)与电池宏观信息(电压，电流，
电阻)的关系，尽管相当精确，但是需要经过复杂的计算，很难在电动车辆上实现电路电压的在线估计。
[0004] 综上所述，为保证高效实时的对动力电池的电量状态(SOC)、
健康状态(SOH)以及功率状态(SOP)进行估算的同时作出相应策略以满足驾驶人员的需求，设计一种实时，准确的电池开路电压估算模型对于提高电池管理系统的性能意义重大。
发明内容
[0005] 本发明的主要目的是提供一种基于气液动力学电池模型在线估算开路电压的方法及系统，将气液动力学电池模型动态估算特性强与微步长安时积分法精度高相结合，实
现对电池开路电压在线高精度估算，在估算电池开路电压之前无需对模型参数和初值进行
辨识或识别，该估算方法将极大地提高模型对电池的适应性。
[0006] 为了实现上述发明目的，本发明提出一种技术解决方案为：一种基于气液动力学电池模型在线估算开路电压的方法，包括以下步骤：
[0007] 步骤S1：通过气液动力学开路电压电池模型估算电池第一开路电压；
[0008] 步骤S2：依据第一开路电压查“开路电压-容量”表得到第一容量；
[0009] 步骤S3：以第一容量为初值，利用微步长安时积分法得到第二容量；
[0010] 步骤S4：依据第二容量查“容量-开路电压”表得到第二开路电压；
[0011] 步骤S5：通过气液动力学端电压电池模型估算电池端电压；
[0012] 步骤S6：计算气液动力学开路电压电池模型的
扩展卡尔曼滤波
算法的系数矩阵和协方差矩阵并更新气液动力学开路电压电池模型的扩展卡尔曼滤波算法的系数矩阵和协
方差矩阵；
[0013] 循环上述步骤S1至S6的过程，完成气液动力学开路电压电池模型参数在线辨识和电池开路电压在线估算。
[0014] 上述方案中，所述步骤S1中气液动力学开路电压电池模型为：
[0015] P2＝P0(i)-k3v(i)-k4|v(i)|v(i)
[0016]
[0017]
[0018]
[0019] 其中，P0为端电压；k1，k2，k3，k4为模型参数，记Para＝[k1，k2，k3，k4]为模型参数矩阵，v为电流，充电时v>0,放电时v<0，P1为初始开路电压，P3为第一开路电压，i为计数序号。
[0020] 检测动力电池端电压和电流，模型参数矩阵Para＝[k1，k2，k3，k4]的初始值为随机赋予，输入到气液动力学开路电压电池模型中，估算第一开路电压P3；所述步骤S2中，依据第一开路电压P3查“开路电压-容量”表得到第一容量First_Ah。
[0021] 上述方案中，所述步骤S3中以第一容量First_Ah为初值，利用微步长安时积分法通过以下公式得到第二容量Next_Ah：
[0022] Next_Ah＝First_Ah+v(i+1)×Δt/3600 公式十一
[0023] 公式十一中，Δt为电流
采样时间间隔，Δt＝t(i+1)-t(i)。所述步骤S4依据第二容量Next_Ah查“容量-开路电压”表得到第二开路电压Next_OCV。
[0024] 上述方案中，所述步骤S5中气液动力学端电压电池模型为：
[0025]
[0026] 公式十三中，P0_g为估算的端电压，设充电时v>0，放电时v<0。
[0027] 上述方案中，所述步骤S6中计算气液动力学开路电压电池模型的扩展卡尔曼滤波算法的系数矩阵和协方差矩阵通过下面公式计算：
[0028]
[0029] 公式十四中，C为协方差矩阵；Pn为中间传递矩阵，H为P0_g关于参数k1,k2,k3,k4的偏导数矩阵；R为增量因子；I为单位矩阵；K为系数矩阵。
[0030] 上述方案中，所述步骤S6中更新气液动力学开路电压电池模型的参数和扩展卡尔曼滤波算法的协方差矩阵具体为：
[0031] 模型参数矩阵更新为：Para＝Parα+K×(P0(i+1)-P0g) 公式十五
[0032] 协方差矩阵更新为：C＝Pn。
[0033] 一种实现所述基于气液动力学电池模型在线估算开路电压的方法的系统，包括信号采集模
块、开路电压估算模块和显示模块；
[0034] 所述
信号采集模块用于采集电池的电流和电压；
[0035] 所述信号采集模块与开路电压估算模块连接并将采集的电流和电压信号传送到开路电压估算模块，所述开路电压估算模块根据气液动力学开路电压电池模型计算出开路
电压值；
[0036] 所述开路电压估算模块与显示模块连接，将电池电流、电压和开路电压值发送给显示模块显示。
[0037] 上述方案中，所述信号采集模块包括电流
传感器和电压传感器。
[0038] 上述方案中，所述开路电压估算模块包括
单片机。
[0039] 与
现有技术相比，有益效果为：本发明无需对模型参数和初值进行离线辨识或识别，极大地提高模型对电池的适应性；求解系数矩阵K时，(H×P×H′+R)-1为求一个实数的逆矩阵，计算过程非常简单，保证在线估算实时性；将气液动力学电池模型动态估算特性强与微步长安时积分法精度高相结合，实现对电池开路电压在线高精度估算。
附图说明
[0041] 图2为气液动力学电池模型原理图；
[0042] 图3为开路电压-容量表的绘图；
[0043] 图4为容量-开路电压表的绘图；
[0045] 图6为本发明估算电池的端电压和开路电压效果图；
[0046] 图7为本发明系统框图。
具体实施方式
[0047] 下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附
图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
[0048] 在本发明的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“轴向”、“径向”、“竖直”、“水平”、“内”、“外”等指示的方位或
位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为
对本发明的限制。此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中，“多个”的含义是两个或两个以上，除非另有明确具体的限定。
[0049] 在本发明中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。
[0050] 图1所述为本发明所述基于气液动力学电池模型在线估算开路电压的方法的一种较佳实施方式，所述基于气液动力学电池模型在线估算开路电压的方法包括如下步骤：
[0051] 步骤S1，通过气液动力学开路电压电池模型估算电池第一开路电压包括：
[0052] 推导所得的气液动力学开路电压电池模型：
[0053] 将气液动力学模型的充放气过程对应起电池的充放电过程，如图2，所述气液动力学电池模型为一个筒状密闭容器1，在容器的顶部安装有管道2和
阀门3，所述容器内装有Vw体积的液体4，剩下的容积V是压强为P，物质的量为n，
密度为ρ的气体5；其中，气体流动的阻力系数为μ，管道外部管口6的气体压强为P0，气体7溶解于液体中的物质的量为nj；
[0054] 设在t1时刻，气液动力学电池模型处于稳态，此时气体压强为P1、气体物质的量为n1、气体溶解于液体中的物质的量为nj1；打开所述容器的阀门，时间长度为Δt，模型中气体向外放出，流速为v，气体流动的阻力系数为μ，外部管口压强为P0，在t2＝t1+Δt时刻关闭阀门，此时气体压强为P2、气体物质的量为n2；经过数十小时到t3时，气液动力学电池模型再次达到稳态，此时容器内气体压强为P3，气体物质的量为nj3；
[0055] 气液动力学电池模型参数与电池参数对应关系：外部管口压强P0对应电池端电压，气体流速v对应电池电流，气体压强P1对应电池初始开路电压，气体压强P3对应电池第一开路电压；P2为中间变量。
[0056] 在模型推导之前需要作出合理假设：其一，因为在物理方程中所有参数具有实际物理意义，因此所有参数非负；其二，因为Δt时间很短，因此在打开阀门时从液体中析出的气体忽略不计；
[0057] 推导气液动力学电池模型的关键是建立P3与P1、v和P0之间的显函数关系；
[0058] 在t1时刻，可以写出气体间隙填充
溶解度方程：
[0059]
[0061] 有效间隙度，
[0062] bm：气体分子的范德华体积，
[0063] R：热力学常数；
[0064] 在t2时刻，可以写出理想气体
状态方程为：
[0065] P2V＝n2RT 公式二
[0066] 气体流动的伯努利方程：
[0067]
[0068] 在t3时刻，可以写出理想气体状态方程为：
[0069] P3V＝n3RT 公式四
[0070] 气体间隙填充溶解度方程：
[0071]
[0072] 在阀门关闭瞬间，容器内气体物质的量为n2；关闭阀门之后，从液体中析出的气体物质的量为nj1-nj3；因此在t3时物质的量关系为n3＝n2+(nj1-nj3)；在t2到t3时间段，容器内改变的压强为 也就是 将t2和t3时刻气体间隙填充溶解度方程代入得 将该式恒等变换得
[0073]
[0074] 为了简化方程，设 公式六化简为：
[0075]
[0076] 因为k2+P3＞0和k2+P1＞0，因此公式七可
整理为：
[0077]
[0078] 设a＝1， 因为ac＜0，且b2-4ac＞0，因此上式有且仅有唯一正实根，为：
[0079]
[0080] 以上推导过程以打开阀门放气过程为例，其实，打开阀门从外界
泵入气体进入容器内的推到结果与上述结果一致，此处不再赘述；但是，气体泵入容器内的伯努利方程为：
[0081]
[0082] 因此，可以将充气过程与放气过程的伯努利方程统一为公式十：
[0083]
[0084] 设充电时v>0,放电时v<0；
[0085] 对推导的气液动力学电池模型进行简化，设
[0086] 因此，推导所得的气液动力学开路电压电池模型：
[0087] P2＝P0(i)-k3v(i)-k4|v(i)|v(i)
[0088]
[0089]
[0090]
[0091] 此处，P0为端电压；记Para＝[k1,k2,k3,k4]为模型参数，v为电流，充电时v>0,放电时v<0，P1为初始开路电压，P3为第一开路电压,i为计数序号(i＝1、2、3…)；初始Para＝[k1,k2,k3,k4]＝[0.005,0.005,0.005,0.005],P1＝4.0；以上初始值为随机赋予，不局限于本
专利列出的具体值；
[0092] 检测动力电池端电压和电流，输入到气液动力学开路电压电池模型中，估算第一开路电压P3；
[0093] 步骤S2，如图3依据第一开路电压P3查“开路电压-容量”表得到第一容量First_Ah；
[0094] 步骤S3，以第一容量为初值，利用微步长安时积分法得到第二容量Next_Ah：
[0095] Next_Ah＝First_Ah+v(i+1)×Δt/3600 公式十一
[0096] 公式十一中，Δt为电流采样时间间隔，Δt＝t(i+1)-t(i)。
[0097] 步骤S4，如图4，依据第二容量Next_Ah查“容量-开路电压”表得到第二开路电压Next_OCV；
[0098] 步骤S5，通过气液动力学端电压电池模型估算电池端电压：
[0099] 推导所得的气液动力学端电压电池模型为：
[0100] 根据公式七和公式十恒等变换得到公式十二，
[0101]
[0102] 将步骤S4得到的Next_OCV赋值给P3，将步骤S2得到的P3赋值给P1，将v(i+1)赋值给v(i)得到公式十三：
[0103]
[0104] 公式十三中，P0_g为估算的端电压；设充电时v>0，放电时v<0；
[0105] 步骤S6，所述计算气液动力学开路电压电池模型的扩展卡尔曼滤波算法的系数矩阵和协方差矩阵包括：
[0106]
[0107] 此处，C为协方差矩阵；Pn为中间传递矩阵，H为P0_g关于参数k1,k2,k3,k4的偏导数矩阵，H′为H矩阵的转置；R为增量因子；I为单位矩阵；K为系数矩阵；
[0108] 初始协方差矩阵C＝[1,0,0,0；0,1,0,0；0,0,1,0；0,0,0,1],R＝0.001，以上初始值为随机赋予，不局限于本专利列出的具体值；
[0109]
[0110] 所述更新气液动力学开路电压电池模型的参数和扩展卡尔曼滤波算法的协方差矩阵具体算法为：
[0111] Para＝Para+K×(P0(i+1)-P0_g) 公式十五
[0112] C＝Pn。
[0113] 如图5，循环上述步骤S1到6过程，完成模型参数在线辨识和电池开路电压在线估算。
[0114] 如图6是本发明估算电池的端电压和开路电压效果图，图中圆形实线为实验测试的端电压，叉形实线为估算的端电压，方形实线为估算的开路电压，圆形实线与叉形实线几乎是重合的，表明采用本发明方案能够现实电池端电压高精度估算，因为在步骤一中估算
的第一开路电压P3和步骤五中估算的端电压P0_g采用同一组参数，且开路电压估算模型与
端电压估算模型均由公式七和公式十推导得出，因此本发明方法对电池的充电和放电端电
压(方形实线所示)均具有非常高的估算精度。
[0115] 如图7所示一种实现所述基于气液动力学电池模型在线估算开路电压的方法的系统，包括信号采集模块、开路电压估算模块和显示模块；所述信号采集模块用于采集电池的电流和电压；所述信号采集模块与开路电压估算模块连接并将采集的电流和电压信号传送
到开路电压估算模块，所述开路电压估算模块根据开路电压估算方程计算出开路电压值；
所述开路电压估算模块与显示模块连接，将电池电流、电压和开路电压值发送给显示模块
显示。
[0116] 所述信号采集模块包括电流传感器和电压传感器。
[0117] 所述开路电压估算模块包括单片机，所述单片机优选为STM32。将基于气液动力学电池模型在线估算开路电压的方法在
硬件上实现，可以在Keil uVision5开发平台上运用C
语言编写的代码在STM32单片机上实现的。
[0118] 所述开路电压估算模块具体为：
[0119] 首先加载STM32单片机库函数文件，运用库函数配置STM32单片机寄存器，编写时钟函数、
定时器函数、延迟函数、存储函数、数据校验函数、OCV估算函数和主函数等；
[0120] ①将电流传感器连接到信号采集卡上，采集卡可以直接采集
单体电池电压，优选的，单体电池电压范围在0—5V以内；
[0121] ②采集卡与STM单片机串口相连，通讯方式选择RS-232，将电池的电流、电压信号传给单片机；
[0122] ③STM32单片机主函数读取电池的电流、电压信号，调用OCV估算函数算出当前输入下的开路电压值；将电池电流、电压和算出的开路电压值写入内存卡中，并把电池电流、电压和算出的开路电压值发送给上位机的显示模块显示；
[0123] ④如此循环第①—③步，完成
电池组实时开路电压估算。
[0124] 所述上位机是基于Microsoft Visual Studio平台开发的，用于显示电池组端电压、开路电压、所有
串联单体电池的开路电压和拟合的电池最低开路电压；
[0125] 所述单片机包括：2n位单片机，n＝1,2,3...,以及各种ARM
内核的运算单元；
[0126] 运用的信号通讯协议包括：RS-485、CAN、TCP、modbus、MPI、串口通信等。
[0127] 本发明将气液动力学电池模型动态估算特性强与微步长安时积分法精度高相结合，实现对电池开路电压在线高精度估算，在估算电池开路电压之前无需对模型参数和初
值进行辨识或识别，该估算方法将极大地提高模型对电池的适应性。
[0128] 应当理解，虽然本
说明书是按照各个实施例描述的，但并非每个实施例仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说
明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以
理解的其他实施方式。
[0129] 上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施例的具体说明，它们并非用以限制本发明的保护范围，凡未脱离本发明技艺精神所作的等效实施例或变更
均应包含在本发明的保护范围之内。
