一种锂离子电池机理建模方法
一种锂离子电池机理建模方法
【技术领域】
[0001] 本发明属于电动汽车锂离子动力电池技术领域，更具体地说，本发明涉及一种锂 离子电池机理建模方法。
【背景技术】
[0002] 动力电池的性能对电动汽车的整车动力性、续驶里程和安全性影响很大。电池的 荷电状态和健康状态估计的准确性、寿命的预测精度、安全稳定的运行都与动力电池的建 模精度密切相关，因此，建立精确的电池模型对电池的研宄具有重要意义。
[0003] 目前，常用的电池模型有经验模型、多物理场耦合模型和电化学模型。经验模型不 考虑物理、化学反应原理，依据所采集的实验数据，拟合出锂离子电池的特性曲线，建立电 池的经验模型，该模型结构简单，但预测能力差；多物理场耦合模型同时考虑了锂离子电池 的内部电场和温度场反应和分布过程，基于电池内部电化学反应过程的生热原理和传热学 原理，分析锂离子电池内部的热特性，该模型主要用于模拟电池在正常工作状态下的内部 温度变化情况，只适用于小功率情况，使用范围有限；电化学模型是基于多孔电极原理、浓 溶液原理和欧姆定律等理论建立的，由偏微分方程、常微分方程、代数方程组成，该模型精 度高，可以较为准确的描述电池内部电化学反应机理和外部反应特性，但模型参数相互耦 合，且与电池结构、尺寸、所用材料等因素有关，模型求解过程复杂，计算量大，给实车上应 用带来困难。
[0004] 在锂离子电池模型中，由于正极和负极活性区域的表面积、正极和负极固相锂离 子扩散系数与正极和负极反应速率常数无法直接测量，且会随着锂离子电池的使用而变 化，所以需要对这些参数进行辨识。目前常用的最小二乘法和极大似然参数辨识方法，要求 目标函数连续可导，极容易使辨识结果陷入局部最小的情况，从而难以获得全局最优的辨 识结果；基于神经网络的参数辨识方法可以高精度逼近非线性函数，但是神经网络结构的 确定、样本数据的选择和神经网络训练算法等问题尚未有效解决，影响了该方法的应用；群 体智能优化算法从生物生命演化过程中得到启示，模拟生物进化行为，是一种多智能体并 行的算法，鲁棒性好。
【发明内容】
[0005] 本发明所要解决的技术问题是克服锂离子电池电化学模型结构复杂、参数难以辨 识，经验模型精度低的缺点，本发明提出了一种基于菌群觅食优化算法的锂离子电池单粒 子模型建模方法。
[0006] 为解决上述技术问题，本发明是采用如下技术方案实现的，结合【附图说明】如下：
[0007] 一种锂离子电池机理建模方法，包括以下步骤：
[0008] 步骤1 :建立锂离子电池单粒子模型；
[0009] 所述锂离子电池单粒子模型描述锂离子电池内的机理反应过程和物理特性，锂离 子电池单粒子模型包括固相锂离子扩散方程、巴特勒-伏尔摩Butler-Volmer动力学方程 和端电压表达式；
[0010] 步骤2 :采用三参数抛物线方法简化锂离子电池单粒子模型中的固相锂离子扩散 方程；
[0011] 所述锂离子电池单粒子模型中的固相锂离子扩散方程为二阶偏微分方程，采用三 参数抛物线方法将其简化为常微分方程组和代数方程；
[0012] 步骤3 :采用菌群觅食优化算法辨识锂离子电池单粒子模型中的未知参数；
[0013] 所述锂离子电池单粒子模型中的未知参数为：正极活性区域的表面积Sp，负极活 性区域的表面积Sn，正极固相锂离子扩散系数Ds，p，负极固相锂离子扩散系数Ds，n，正极反应 速率常数kp，负极反应速率常数kn;
[0014] 步骤4 :拟合锂离子单粒子模型的正极开路电压表达式。
[0015] 技术方案中所述固相锂离子扩散方程为：
[0017] 式⑴的边界条件为：
[0020] 其中，i=p，n分别代表锂离子电池的正极和负极；
[0021] y为固相锂离子浓度；
[0022] t为时间；
[0023] Ds，i为固相锂离子扩散系数；
[0024] rf(〇，Ri)为电极内球状粒子的径向距离；
[0025] 氏为电极内球状粒子的半径；
[0026] 上为固相粒子表面的锂尚子孔壁通量。
[0027] 所述固相粒子表面的锂离子孔壁通量i的定义为：
[0029] 其中，I为锂离子电池的充放电电流；
[0030] F为法拉第常数；
[0031] Si*电极内活性区域的表面积；
[0032] 所述巴特勒-伏尔摩Butler-Volmer动力学方程为：
[0034] 其中，匕为电极反应速率常数；
[0035] cs，imax为固相最大锂离子浓度；
[0036] 为固相粒子表面锂离子浓度；
[0037] c为液相锂离子浓度；
[0038] R为普适气体常数；
[0039]T为温度；
[0040] 1为电极过电压；
[0041] 求解式（5)，得电极过电压表达式为：
[0044] 由于锂离子电池单粒子模型忽略了与液相扩散相关的反应过程，液相电压为零， 电极过电压与电极固相电压、电极开路电压之间的关系为：
[0046] 其中，（i>si为电极固相电压；
[0047] 仏为电极开路电压，其表达式是关于电极固相粒子表面锂离子浓度的非线性函 数。
[0048] 锂离子电池正极固相电压与负极固相电压之间的差值为锂离子电池端电压，为了 计算方便，引入荷电状态变量9 i:
[0050] 所述锂离子电池端电压表达式为：
[0051] v= <i>s；p-<i)s；n= (up( 0p)-un( 0 n)) + (np-nn) (9)
[0052] 式（l)-(4)、（6)、（8)-(9)为以电流为输入，端电压为输出的锂离子电池单粒子模 型数学表达式。
[0053] 技术方案中所述锂离子电池单粒子模型固相锂离子扩散方程为二阶偏微分方程， 采用三参数抛物线方法将其简化为常微分方程组和代数方程；
[0054] 采用三参数抛物线方法表示固相锂离子浓度为：
[0056] 其中，a(t)、b(t)、c(t)是待求解的系数；
[0057]式（10)满足式⑵中r= 0处的边界条件，将式（10)代入式⑴和式（3)，得：
[0060] 根据式（10)，变量固相锂离子平均浓度1(0、固相粒子表面锂离子浓度和 固相锂离子体积平均浓度通量_)可以表示为：
[0064] 由式（13)_(15)，求得a(t)、b(t)和c(t)的表达式为：
[0068] 则固相粒子锂离子平均浓度满足：
[0070] 固相锂离子体积平均浓度通量满足：
[0072] 固相粒子表面锂离子浓度为：
[0074] 式（19)-(21)为锂离子电池单粒子模型简化的固相锂离子扩散方程。
[0075] 技术方案中锂离子电池单粒子模型中的未知参数向量为0 = (Sp，Sn，Ds，p,Ds，n,kp，kn) ; 0为锂离子电池单粒子模型中的未知参数向量；
[0076] 所述采用菌群觅食优化算法辨识锂离子电池单粒子模型中的未知参数过程中，目 标函数为：
[0078] 其中，yi为锂离子电池充放电实验第i个采样点的端电压；
[0079] f(Xi，0 )为锂离子电池单粒子模型第i个采样点的端电压；
[0080] XiS充放电实验第i个采样点的输入电流；
[0081] 0 = (Sp，Sn，Ds，p,Ds，n,kp，kn)为锂离子电池单粒子模型中待辨识的参数向量；
[0082] 将目标函数作为菌群觅食优化算法的适应度函数，基于实测的锂离子电池充放电 电流和端电压数据，得到使式（22)中目标函数最小的参数向量0即为参数辨识结果。
[0083] 技术方案中所述菌群觅食优化算法模拟了细菌觅食的趋药性、聚集、复制和消散 四种智能行为；
[0084] 所述趋药性行为包括翻转行为和前进行为；翻转行为是细菌沿着任意一个新的方 向运动，前进行为是细菌沿着与上一步相同的方向运动；
[0085] -次翻转行为优化后参数向量0值为：
[0086] 0i(j+l，k，1) = 0i(j，k，1)+C⑴伞⑴(23)
[0087] 其中，0W，k，1)为第j(0 <j彡Nc)步趋药性，第k(0 <k彡Nre)步复制，第1 (0 < 1 <NJ步消散行为对应的第i个参数向量0 ;
[0088] 0i(j+l，k, 1)为一次翻转行为优化后第i个参数向量0值；
[0089] N。为趋药性行为总步数；
[0090] 为复制行为总步数；
[0091] Ned为消散行为总步数；
[0092] 巾（i)为单位长度的 随机方向向量；
[0093] C(i)为参数向量优化步长；
[0094] 所述聚集行为是菌群觅食过程中，细菌个体之间通过释放引诱剂和排斥剂来完成 菌群的聚集行为；
[0095] 细菌间聚集行为的数学表达式为：
[0097] 其中，JJ0 )为每两个细菌间的适应度函数；
[0098] S为参数变量0的个数；
[0099] dattraet为引诱剂扩散深度；
[0100] ?attrac;t为引诱剂的扩散率的大小；
[0101 ]D为待辨识参数的维数；
[0102] 9 为参数向量9在第m维空间的分量；
[0103] 化为第i个参数向量9在m维空间的分量；
[0104] h,epellant为排斥剂扩散高度；
[0105] ?Mpellant为排斥剂扩散率的大小；
[0106] 所述复制行为是一个周期的趋药性行为后，根据细菌健康函数值的大小进行复制 与淘汰，健康函数值较小的细菌开始复制行为，生成与父代完全相同的子代，没有复制的细 菌将被淘汰，复制细菌的个数与淘汰细菌的个数相等，以此来维持菌群的个数不变；
[0107] 健康函数为趋药性行为中参数向量0的适应度函数值的和：
[0109] 其中，表示第i个参数向量0的健康函数，，越大表示参数向量0优化效 果越差；
[0110] J(i，j，k，l)表示对应第i个参数向量0的适应度函数值；
[0111] 所述消散行为是指细菌被外力杀死或者被驱散到新的区域中的行为；消散行为破 坏了细菌的趋药性过程，但是，这使得细菌可能寻找到食物更加丰富的区域。
[0112] 技术方案中所述菌群觅食优化算法辨识参数的过程为：
[0113] 步骤1 :在参数优化范围内，随机给出参数向量0 = (Sp,Sn,Ds，p,Ds，n,kp,kn)的初 始值；
[0114] 步骤2 :通过细菌趋药性行为优化参数向量0，使目标函数减小：
[0115] 首先，根据式（22)计算目标函数值。考虑菌群聚集行为，根据式（24)对菌群适应 度函数进行修正，加快参数辨识速度；
[0116] 然后，参数向量0随机选择一个方向进行翻转行为优化，根据式（23)计算优化后 的参数向量0，并计算目标函数值J(0)。如果J(0)减小，则在该方向上继续前进行为优 化，直到J(0)不再减小，或参数向量0在该方向上达到了最大的前进行为优化次数，否则 参数向量0将随机选择另外一个方向进行翻转行为优化；
[0117] 步骤3 :参数向量0优化后，通过细菌复制行为对参数向量0进行复制与淘汰：
[0118] 首先，根据式（25)计算参数向量0的健康函数值，并对各参数向量0的健康函 数值按从大到小的原则进行排序；
[0119] 然后，对健康函数值低的参数向量0进行复制，将健康函数值高的参数向量0淘 汰，淘汰的参数向量0的个数为*5； =|，以保证参数向量0的总数不变；
[0120] 步骤4:参数向量0复制与淘汰后，通过细菌的消散行为将参数向量0按消散概 率ped重新随机分布到寻优区间，使算法逃逸出局部极值，从而求出全局最优点；
[0121] 步骤5 :参数向量0重新分布后，计算参数向量0对应的目标函数值，使目标函 数最小的参数向量0即为参数辨识结果。
[0122] 技术方案中所述锂离子电池单粒子模型的正极开路电压表达式的拟合步骤为：
[0123] (1)采用0. 1C电流对锂离子电池进行恒流放电实验，测得锂离子电池端电压
[0124] (2)在矩阵实验室Matlab中搭建锂离子电池单粒子模型，将正、负极开路电压的 经验公式代入到所搭建的锂离子电池单粒子模型中；采用0. 1C电流对锂离子电池单粒子 模型进行恒流放电仿真实验，得到锂离子电池单粒子模型端电压V;
[0125] (3)当锂离子电池进行小电流充放电实验时，端电压与开路电压差近似相等，而且 负极开路电压较小，则锂离子电池单粒子模型正极开路电压通过式（28)计算：
[0127] 其中，Up( 0p)是锂离子电池单粒子模型正极开路电压；
[0128] 是正极开路电压经验值；
[0129] 是通过恒流放电实验测得的锂离子电池端电压；
[0130] V是锂离子电池单粒子模型端电压；
[0131] 将拟合得到锂离子电池单粒子模型正极开路电压表达式代入到锂离子电池单粒 子模型中，得到锂离子电池单粒子模型端电压；
[0132] (4)将（3)中得到的锂离子电池单粒子模型端电压与（1)中实验测得的锂离子电 池端电压比较，若差值小于±〇.〇5V，则（3)中拟合得到的锂离子电池单粒子模型正极开路 电压表达式为所求，否则，用该表达式代替正极开路电压经验公式，返回步骤（2)，重复以上 步骤，直至端电压差小于±〇.〇5V。
[0133] 技术方案中建立锂离子电池单粒子模型的条件为：
[0134] 1)假设锂离子电池电极由多个具有相同大小和动力学特性的球状粒子组成，并且 电流通过电极时在所有活性粒子内均匀分布；
[0135] 2)假设在固相颗粒内或颗粒之间电压降为零；
[0136] 3)假设在整个锂离子电池内部液相锂离子浓度恒定，并且在时间和空间上均匀分 布；
[0137]4)忽略液相电压对电池端电压的影响；
[0138] 5)忽略锂离子电池充放电过程中产生的热量。
[0139] 与现有技术相比本发明的有益效果是：
[0140] (1)锂离子电池单粒子模型忽略液相电压对端电压的影响，采用三参数抛物线方 法将锂离子电池单粒子模型中表示固相扩散方程的偏微分方程组简化为常微分方程组和 代数方程，简化了锂离子电池单粒子模型的结构。
[0141] (2)采用菌群觅食优化算法辨识锂离子电池单粒子模型中的未知参数，辨识速度 快，得到了全局最优解。
[0142] (3)锂离子电池单粒子模型能够较准确地模拟锂离子电池的放电行为，锂离子电 池单粒子模型端电压数据与锂离子电池端电压数据可以很好地吻合，最大误差在±〇.〇7V 左右。此外，锂离子电池单粒子模型可以为锂离子电池状态估计，寿命预测，特性分析提供 理论支持。
【附图说明】
[0143] 下面结合附图对本发明作进一步的说明：
[0144] 图1为锂离子电池单粒子模型建模流程图；
[0145] 图2为锂离子电池单粒子模型示意图；
[0146] 图3为锂离子电池单粒子模型仿真模型结构框图；
[0147] 图4为菌群觅食优化算法辨识锂离子电池单粒子模型参数流程图；
[0148] 图5为适应度函数随趋药性行为步数变化曲线；
[0149] 图6为磷酸铁锂电池正极开路电压与正极荷电状态变量之间的关系曲线；
[0150] 图7为磷酸铁锂电池负极开路电压与负极荷电状态变量之间的关系曲线；
[0151] 图8为1. 6Ah磷酸铁锂电池端电压曲线与单粒子模型端电压曲线；
[0152] 图9为1. 6Ah磷酸铁锂电池单粒子模型端电压误差曲线。
【具体实施方式】
[0153] 下面结合附图对本发明作详细的描述：
[0154] 锂离子电池单粒子模型建模流程图参阅图1所示，本发明基于菌群觅食优化算法 的锂离子电池单粒子模型建模步骤如下：
[0155] 1、根据锂离子电池工作原理，建立锂离子电池单粒子模型。
[0156] 锂离子电池主要由正极、负极和隔膜组成，正极活性物质为富锂金属氧化物与相 应金氧化物的混合物质，负极活性物质是石墨或近似石墨结构的碳，隔膜是一种经特殊成 型的高分子薄膜，有微孔结构，锂离子可以自由通过，而电子不能通过。锂离子电池充电 时，正极上有锂离子生成，生成的锂离子通过电池隔膜运动到电池负极，而作为负极的碳呈 现层状结构，它有很多微孔，运动到负极的锂离子就嵌入到碳层的微孔中，嵌入的锂离子越 多，充电量越高。当对锂离子电池进行放电时，嵌在负极碳层中的锂离子析出回到正极，回 到正极的锂离子越多，放电量越高。
[0157] 锂离子电池单粒子模型是指利用一个球状粒子特性来代表整个电极特性而建立 的一种简化的锂离子电池机理模型，锂离子电池单粒子模型示意图如图2所示。建立锂离 子电池单粒子模型的条件为：
[0158] (1)假设锂离子电池电极由多个具有相同大小和动力学特性的球状粒子组成，并 且电流通过电极时在所有活性粒子内均匀分布；
[0159] (2)假设在固相颗粒内或颗粒之间电压降为零；
[0160] (3)假设在整个锂离子电池内部液相锂离子浓度恒定，并且在时间和空间上均匀 分布；
[0161] (4)忽略液相电压对电池端电压的影响；
[0162] (5)忽略锂离子电池 充放电过程中产生的热量。
[0163] 依据上述假设条件，根据锂离子电池内的机理反应过程和物理特性，锂离子电池 单粒子模型包括固相扩散方程，巴特勒-伏尔摩（Butler-Volmer)动力学方程和端电压表 达式。
[0164] 锂离子在电极固相中的扩散为非稳态扩散，根据菲克第二定律，在球形电极上，固 相锂离子扩散方程为：
[0166] 式⑴的边界条件为：
[0169] 其中，i=p，n分别代表锂离子电池的正极和负极；为固相锂离子浓度；t为时 间；Ds;i为固相锂离子扩散系数；r# (〇，Ri)为电极内球状粒子的径向距离况为电极内球 状粒子的半径；1为固相粒子表面的锂离子孔壁通量。
[0170] 固相粒子表面的锂离子孔壁通量1的定义为：
[0172] 其中，I为锂离子电池的充放电电流；F为法拉第常数A为电极内活性区域的表 面积。巴特勒-伏尔摩（Butler-Volmer)动力学方程描述固相活性粒子表面与电解液溶液 临界面处的电化学反应过程，该动力学方程为：
[0174] 其中，&为电极反应速率常数；Cs，imax为固相最大锂离子浓度；为固相粒子表面 锂离子浓度；C为液相锂离子浓度；R为普适气体常数；T为温度；n电极过电压。
[0175]求解式（5)，得电极过电压表达式为：
[0178] 由于锂离子电池单粒子模型忽略了与液相扩散相关的反应过程，液相电压为零， 电极过电压与电极固相电压、电极开路电压之间的关系为：
[0180] 其中，为电极固相电压；为电极开路电压，是关于电极固相粒子表面锂 离子浓度的非线性函数。
[0181] 锂离子电池正极固相电压与负极固相电压之间的差值为锂离子电池端电压，为了 计算方便，引入荷电状态变量9 i:
[0183] 锂离子电池端电压表达式为：
[0184] v= <i>s；p-<i)s；n= (up( 0p)-un( 0 n)) + (np-nn) (9)
[0185] 式（l)-(4)、（6)、（8)-(9)即为以电流为输入，端电压为输出的锂离子电池单粒子 模型数学表达式，锂离子电池单粒子模型仿真模型结构框图参阅图3所示。
[0186] 2、采用三参数抛物线方法简化锂离子电池单粒子模型中的固相扩散方程。
[0187] 锂离子电池单粒子模型固相扩散方程为二阶偏微分方程，采用三参数抛物线方法 将其简化为常微分方程组和代数方程。采用三参数抛物线方法表示固相锂离子浓度为：
[0189] 其中，a(t)、b(t)、c(t)是待求解的系数。
[0190] 显然，式（10)满足式⑵中r= 0处的边界条件，将式（10)代入式⑴和式（3)， 得：
[0193] 根据式（10)，变量固相锂离子平均浓度己(0、固相粒子表面锂离子浓度和固 相锂离子体积平均浓度通量砂)可以表示为：
[0197] 由式（13)-(15)，可求得3(〇、13(〇、(^)的表达式为：
[0201] 结合式（1)-⑶求解，得：
[0205] 式（19)-(21)即为锂离子电池单粒子模型简化的固相扩散方程。
[0206] 3、采用菌群觅食优化算法辨识锂离子电池单粒子模型中的未知参数。
[0207] 锂离子电池单粒子模型中的正极和负极活性区域的表面积、正极和负极固相锂 离子扩散系数、正极和负极反应速率常数无法直接测量，并且会随着锂离子电池的使用 而变化，所以需要对这些参数进行辨识。本发明采用菌群觅食优化算法辨识参数0 = (Sp，Sn，Ds，p,Ds，n,kp，kn) 〇
[0208] 细菌在觅食时希望在最短的时间内获得最多的营养。为了达到这个目标，细菌 要决定觅食策略，通过感应自身周围的化学物质浓度，寻找食物所在的区域，找到后决 定进入该区域还是继续寻找，进入后什么时候离开该区域。通过模拟细菌的觅食策略， Passino在文献Biomimicryofbacterialforagingfordistributedoptimizationand control(ControlSystems,IEEE, 2002, 22(3) : 52-67.)中提出 了菌群觅食优化算法，算法 流程图参阅图4所示。
[0209] 在菌群觅食优化算法辨识锂离子电池单粒子模型未知参数过程中，目标函数为：
[0211] 其中，yi为锂离子电池充放电实验第i个采样点的端电压；f(Xi，0)为锂离子电 池单粒子模型第i个采样点的端电压；^为充放电实验第i个采样点的输入电流；0 = (Sp，Sn，Ds，p,Ds，n,kp，kn)为锂离子电池单粒子模型中待辨识的参数向量。
[0212] 将目标函数作为菌群觅食优化算法的适应度函数，基于实测的锂离子电池端电压 和充放电电流数据，采用菌群觅食优化算法，得到使式（22)中目标函数最小的参数向量0 即为参数辨识结果。菌群觅食优化算法模拟了细菌觅食的趋药性、聚集、复制和消散四种智 能行为。
[0213] ⑴趋药性行为
[0214] 由于细菌经常生存在化学引诱剂环境中，它对环境的应激反应称为趋药性行为。 细菌的基本趋药性行为包括：翻转行为和前进行为。翻转行为是细菌沿着任意一个新的方 向运动，而前进行为是细菌沿着与上一步相同的方向运动。
[0215] -次翻转行为优化后参数向量0值为：
[0216] 0i(j+l，k，1) = 0i(j，k，1)+C⑴伞⑴(23)
[0217] 其中，0i(j，k，l)为第j(0<j<N。）步趋药性，第k(0<k<Nj步复制，第1(0 < 1 <NJ步消散行为对应的第i个参数向量0 ; 0 1 (j+1,k, 1)为一次翻转行为优化后第 i个参数向量9 ;N。为趋药性行为总步数；NM为复制行为总步数；Nm为消散行为总步数； (i)为单位长度的随机方向向量；C(i)为参数向量优化步长。
[0218] ⑵聚集行为
[0219] 在菌群觅食过程中，细菌个体之间通过释放引诱剂和排斥剂来完成菌群的聚集行 为。弓丨诱剂使细菌聚集在一起，甚至出现"抱团"现象；排斥剂使每个细菌都有一定的位置， 使其能在该位置上获取营养，以维持生存。
[0220] 细菌间聚集行为的数学表达式为：
[0222] 其中，上。（0 )为每两个细菌间的适应度函数；S为参数变量0的个数；dattMet为 引诱剂扩散深度，即释放引诱剂的量；《attM;t为引诱剂的扩散率的大小；D为寻优空间的维 数，即待辨识参数的维数；0m为参数向量0在第m维空间的分量；％为第i个参数向量0 在m维空间的分量；hMIrellant为排斥剂扩散高度，即释放排斥剂的量；《 MIrellant为排斥剂扩散 率的大小。
[0223] ⑶复制行为
[0224] 一个周期的趋药性行为后，根据细菌健康函数值的大小进行复制与淘汰，健康函 数值较小的细菌开始复制行为，生成与父代完全相同的子代，没有复制的细菌将被淘汰，复 制细菌的个数与淘汰细菌的个数相等，以此来维持细菌的个数不变。
[0225] 健康函数为趋药性行为中参数向量0的适应度函数值的和：
[0227] 其中，/Lm表示第i个参数向量9的健康函数，越大表示参数向量9优化效 果越差；J(i,j,k, 1)表示对应第i个参数向量0的适应度函数值。
[0228] ⑷消散行为
[0229] 消散行为是指细菌被外力杀死或者被驱散到新的区域中的行为。消散行为破坏了 细菌的趋药性过程，但是，这使得细菌可能寻找到食物更加丰富的区域。
[0230] 菌群觅食优化算法辨识参数的过程可以归纳为：
[0231] 步骤1在参数优化范围内，随机给出参数向量0 = (Sp,Sn,Ds，p,Ds，n,kp,kn)的初始 值；
[0232] 步骤2通过细菌趋药性行为优化参数向量0，使目标函数减小：
[0233] 首先，根据式（22)计算目标函数值。考虑菌群聚集行为，根据式（24)对菌群适应 度函数进行修正，加快参数辨识速度；
[0234] 然后，参数向量0随机选择一个方向进行翻转行为优化，根据式（23)计算优化后 的参数向量0，并计算目标函数值J(0)。如果J(0)减小，则在该方向上继续前进行为优 化，直到J(0)不再减小，或参数向量0在该方向上达到了最大的前进行为优化次数，否则 参数向量0将随机选择另外一个方向进行翻转行为优化；
[0235] 步骤3参数向量0优化后，通过细菌复制行为对参数向量0进行复制与淘汰：
[0236] 首先，根据式（25)计算参数向量0健康函数值，并对各参数向量0的健康函数 值按从大到小的原则进行排序；
[0237]然后，对健康函数值低的参数向量0进行复制，将健康函数值高的参数向量0淘 汰，淘汰的参数向量Q个数为，以保证参数向量Q的总数不变；
[0238] 步骤4参数向量0复制与淘汰后，通过细菌的消散行为将参数向量0按消散概 率ped重新随机分布到寻优区间，使算法逃逸出局部极值，从而求出全局最优点；
[0239] 步骤5参数向量0重新分布后，计算参数向量0对应的目标函数值，使目标函数 最小的参数向量0即为参数辨识结果。
[0240] 4、拟合锂离子电池单粒子模型的正极开路电压表达式。
[0241] 锂离子电池负极的活性物质为石墨或近似石墨结构的碳，结构稳定，其开路电压 与负极荷电状态变量之间的关系可以用经验公式表示。锂离子电池正极的活性物质为富锂 金属氧化物与相应金属氧化物的混合物质，并且因为锂离子电池的正极开路电压决定了电 池的开路电压差，所以需要拟合锂离子电池单粒子模型中的正极开路电压表达式。正极开 路电压Up(0p)是正极荷电状态变量0p的非线性函数，本发明通过锂离子电池小电流恒流 放电实验，基于实验测得的锂离子电池端电压与锂离子电池单粒子模型端电压之间的差， 结合正极开路电压经验公式，拟合锂离子电池单粒子模型正极开路电压表达式。
[0242] 当锂离子电池进行小电流充放电实验时，端电压与开路电压差近似相等，而且负 极开路电压较小，则锂离子电池单粒子模型正极开路电压可以通过式（26)计算：
[0244] 其中，Up(0p)是锂离子电池单粒子模型正极开路电压；是正极开路电压经 验值；la是通过恒流放电实验测得的锂离子电池端电压；V是锂离子电池单粒子模型端 电压。
[0245] 以磷酸铁锂电池为例，根据文献Modelingofacommercialgraphite/ LiFeP04cell(JournaloftheElectrochemicalSociety,2011， 158(5):A562_A571)，正、 负极开路电压经验公式为：
[0248] 其中，W吹）为负极开路电压经验值。
[0249] 锂离子电池单粒子模型正极开路电压表达式具体拟合过程为：
[0250] (1)采用0. 1C电流对锂离子电池进行恒流放电实验，测得锂离子电池端电压。
[0251] (2)在矩阵实验室（Matlab)中搭建锂离子电池单粒子模型，将正、负极开路电压 的经验公式代入到所搭建的锂离子电池单粒子模型中。采用〇. 1C电流对锂离子电池单粒 子模型进行恒流放电仿真实验，得到锂离子电池单粒子模型端电压V。
[0252] (3)根据式（26)，拟合得到锂离子电池单粒子模型正极开路电压表达式，并将其 代入到锂离子电池单粒子模型中，得到锂离子电池单粒子模型端电压。
[0253] (4)将（3)中得到的锂离子电池单粒子模型端电压与⑴中实验测得的锂离子电 池端电压比较，若差值小于±〇.〇5V，则（3)中拟合得到的锂离子电池单粒子模型正极开路 电压表达式即为所求，否则，用该表达式代替正极开路电压经验公式，返回步骤（2)，重复以 上步骤，直至端电压差小于±〇.〇5V。
[0254] 下面以一个具体实例解释说明本发明。
[0255] 本实例以正极为LiFeP04，负极为1^(：6的1.6Ah锂离子电池为研宄对象，具体实施 过程如下：
[0256] 1、根据锂离子电池内部电化学反应机理，得到锂离子电池单粒子模型的数学表达 式为式（1)-(4)、（6)、（8)-(9)。
[0257] 2、采用三参数抛物线方法化简锂离子电池单粒子模型，得到化简后的锂离子电池 单粒子模型数学表达式为式（6)、（8) _(9)、（19)-(21)。
[0258] 3、采用菌群觅食优化算法辨识锂离子电池单粒子模型中的未知参数0 = (Sp，Sn，Ds，p,Ds，n,kp，kn) 〇
[0259] 首先对磷酸铁锂电池进行0. 1C恒流放电实验得到电池端电压数据。
[0260] 然后确定锂离子电池单粒子模型中待辨识参数0 = (Sp，Sn，Ds,p，Ds,n，kp，kn)的范 围：SpG[0? 1，l]、SnG[0? 1，l]、Ds，pG[0? 5X1(T14,lXl(T13]、Ds，nG[0? 5X1(T16, 1XKT15]、 kpG[0? 1X10 -11，1X1(T10]、knG[0? 1X10 -12, 1XKT11] 〇
[0261] 其次参数初始化：设辨识参数初始值为61,;,(〇)，待辨识参数向量0的维数D= 6 ;根 据文献"菌群优化算法的研宄"（哈尔滨工业大学，2009.)取参数向量0的个数S= 20， 趋药性行为总步数凡=60,趋药性行为中前进行为优化最大步数NS= 4,复制行为总步数 Nra= 8,消散行为总步数Ned= 2,复制参数向量0组数
消散概率口《1= 〇. 25,引诱 剂扩散深度dattrart= 0.05,引诱剂的扩散率的大小《 attrart= 0. 1，排斥剂扩散高度h_llant = 0.05,排斥剂的扩散率的大小 〇repellant4〇
[0262] 最后根据图4菌群觅食优化算法辨识锂离子电池单粒子模型参数流程，结合实验 数据进行参数辨识，辨识结果为：
[0263] Sp= 0? 4945m2，Sn= 0? 4937m2，Ds，p= 8. 13Xl(T14m2 .s'Ds，n= 6. 64Xl(T16m2 .s' kp= 4. 25Xl〇-nm2.5 ?m〇r0.5 ?s'kn= 2. 09X10_12m2.5 ?m〇r0.5 ?s-1。
[0264] 参数向量适应度函数随趋药性行为步数变化曲线如图5所示，适应度函数值随着 趋药性行为步数的增加不断减小最后趋于稳定。表1为磷酸铁锂电池单粒子模型参数值。
[0265] 表 1
[0266]
[0267] 4、拟合锂离子电池单粒子模型的正极开路电压表达式。
[0268] 经过多次循环迭代和仿真实验，得到磷酸铁锂电池单粒子模型的正极开路电压表 达式为：
[0270] 磷酸铁锂电池单粒子模型正极开路电压与正极荷电状态之间的关系曲线参阅图6 所示，磷酸铁锂电池单粒子模型负极开路电压与负极荷电状态之间的关系参阅图7所示。
[0271] 用上面得到的磷酸铁锂电池单粒子模型模拟电池0. 1C恒流放电行为，端电压曲 线参阅图8所示，误差曲线参阅图9所示。参阅图8，在放电过程中，当t< 10000s时，磷酸 铁锂电池单粒子模型可以准确地模拟电池放电的压降行为，端电压误差在±0.01V以内。 在放电截止时，磷酸铁锂电池正极锂离子浓度增大，根据式（8)，正极荷电状态随之增大，根 据式（29)，正极荷电状态对锂离子电池正极开路电压的影响为指数形式，磷酸铁锂电池正 极开路电压会产生急速下降的趋势。因为正极开路电压决定电池端电压，所以放电截止时 端电压产生陡降行为，磷酸铁锂电池单粒子模型端电压误差增大，但整体上，磷酸铁锂电池 单粒子模型的放电端电压可以与磷酸铁锂电池实际放电端电压的采样点较好地吻合。
【主权项】
1. 一种锂离子电池机理建模方法，其特征在于，包括以下步骤： 步骤1 :建立锂离子电池单粒子模型； 所述锂离子电池单粒子模型描述锂离子电池内的机理反应过程和物理特性，锂离子电 池单粒子模型包括固相锂离子扩散方程、巴特勒-伏尔摩Butler-Volmer动力学方程和端 电压表达式； 步骤2:采用三参数抛物线方法简化锂离子电池单粒子模型中的固相锂离子扩散方 程； 所述锂离子电池单粒子模型中的固相锂离子扩散方程为二阶偏微分方程，采用三参数 抛物线方法将其简化为常微分方程组和代数方程； 步骤3 :采用菌群觅食优化算法辨识锂离子电池单粒子模型中的未知参数； 所述锂离子电池单粒子模型中的未知参数为：正极活性区域的表面积Sp，负极活性区 域的表面积Sn，正极固相锂离子扩散系数Ds，p，负极固相锂离子扩散系数D s，n，正极反应速率 常数kp，负极反应速率常数kn; 步骤4 :拟合锂离子单粒子模型的正极开路电压表达式。2. 根据权利要求1所述的一种锂离子电池机理建模方法，其特征在于： 所述固相锂离子扩散方程为：(1) 式（1)的边界条件为：(2) (3) 其中，i = ρ，η分别代表锂离子电池的正极和负极； Cy为固相锂离子浓度； t为时间； Dsii为固相锂离子扩散系数； rie (OJi)为电极内球状粒子的径向距离； Ri为电极内球状粒子的半径； 1为固相粒子表面的锂离子孔壁通量。 所述固相粒子表面的锂离子孔壁通量Ji的定义为：(4) 其中，I为锂离子电池的充放电电流； F为法拉第常数； S#电极内活性区域的表面积； 所述巴特勒-伏尔摩Butler-Volmer动力学方程为：其中，1^为电极反应速率常数；cEUmax为固相最大锂离子浓度； 固相粒子表面锂离子浓度； c为液相锂离子浓度； R为普适气体常数； T为温度； n i为电极过电压； 求解式（5)，得电极过电压表达式为：由于锂离 子电池单粒子模型忽略了与液相扩散相关的反应过程，液相电压为零，电极 过电压与电极固相电压、电极开路电压之间的关系为： >ι=Φ,J-UiKn (7) 其中，为电极固相电压； 仏为电极开路电压，其表达式是关于电极固相粒子表面锂离子浓度的非线性函数； 锂离子电池正极固相电压与负极固相电压之间的差值为锂离子电池端电压，为了计算 方便，引入荷电状态变量9i:(8) 所述锂离子电池端电压表达式为： ν= Φ3,Ρ_Φ3,η= (υρ( θΡ)-υη( θη)) + (ηρ-ηη) (9) 式（1)-(4)、（6)、（8)-(9)为以电流为输入，端电压为输出的锂离子电池单粒子模型数 学表达式。3.根据权利要求1所述的一种锂离子电池机理建模方法，其特征在于： 采用三参数抛物线方法表示固相锂离子浓度为：(10) 其中，a (t)、b (t)、c (t)是待求解的系数； 式（10)满足式⑵中r = 0处的边界条件，将式（10)代入式⑴和式（3)，得：根据式（10)，变量固相锂离子平均浓度己(0、固相粒子表面锂离子浓度和固相 锂离子体积平均浓度通量获t)可以表示为：由式（13)-(15)，求得a(t)、b(t)和c(t)的表达式为：则固相粒子锂离子平均浓度满足：(19) 固相锂离子体积平均浓度通量满足：(20) 固相粒子表面锂离子浓度为：(21) 式（19)-(21)为锂离子电池单粒子模型简化的固相锂离子扩散方程。4. 根据权利要求1所述的一种锂离子电池机理建模方法，其特征在于： 锂离子电池单粒子模型中的未知参数向量为Θ = (Sp，Sn，Ds，p，Ds， n，kp，kn) ; Θ为锂离子 电池单粒子模型中的未知参数向量； 所述采用菌群觅食优化算法辨识锂离子电池单粒子模型中的未知参数过程中，目标函 数为：(22) 其中，Yi为锂离子电池充放电实验第i个采样点的端电压； f(Xi，Θ )为锂离子电池单粒子模型第i个采样点的端电压； Xi为充放电实验第i个采样点的输入电流； 将目标函数作为菌群觅食优化算法的适应度函数，基于实测的锂离子电池充放电电流 和端电压数据，得到使式（22)中目标函数最小的参数向量Θ即为参数辨识结果。5. 根据权利要求4所述的一种锂离子电池机理建模方法，其特征在于： 所述菌群觅食优化算法模拟了细菌觅食的趋药性、聚集、复制和消散四种智能行为； 所述趋药性行为包括翻转行为和前进行为；翻转行为是细菌沿着任意一个新的方向运 动，前进行为是细菌沿着与上一步相同的方向运动； 一次翻转行为优化后参数向量Θ值为： Θ 1 (j+1，k，1) = Θ i (j，k，I) +C ⑴ Φ ⑴(23) 其中，9\」，1^，1)为第」_(0<_]_彡凡）步趋药性，第1^0<1^彡〇步复制，第1(0 < I < NJ步消散行为对应的第i个参数向量θ ; Θ 1 (j+1, k, 1)为一次翻转行为优化后第i个参数向量Θ值； N。为趋药性行为总步数； Nm为复制行为总步数； Ned为消散行为总步数； Φ (i)为单位长度的随机方向向量； C(i)为参数向量优化步长； 所述聚集行为是菌群觅食过程中，细菌个体之间通过释放引诱剂和排斥剂来完成菌群 的聚集行为； 细菌间聚集行为的数学表达式为：其中，U Θ )为每两个细菌间的适应度函数； S为参数变量Θ的个数； dattrac；t为引诱剂扩散深度； Wattaart为引诱剂的扩散率的大小； D为待辨识参数的维数； 为参数向量Θ在第m维空间的分量； <为第i个参数向量Θ在m维空间的分量； Iwllant为排斥剂扩散高度； ω ^llant为排斥剂扩散率的大小； 所述复制行为是一个周期的趋药性行为后，根据细菌健康函数值的大小进行复制与淘 汰，健康函数值较小的细菌开始复制行为，生成与父代完全相同的子代，没有复制的细菌将 被淘汰，复制细菌的个数与淘汰细菌的个数相等，以此来维持菌群的个数不变； 健康函数为趋药性行为中参数向量Θ的适应度函数值的和：(25) 其中，7LflM表示第i个参数向量Θ的健康函数，越大表示参数向量Θ优化效果越 差； J(i，j，k, 1)表示对应第i个参数向量Θ的适应度函数值； 所述消散行为是指细菌被外力杀死或者被驱散到新的区域中的行为；消散行为破坏了 细菌的趋药性过程，但是，这使得细菌可能寻找到食物更加丰富的区域。6.根据权利要求4所述的一种锂离子电池机理建模方法，其特征在于： 所述菌群觅食优化算法辨识参数的过程为： 步骤1 :在参数优化范围内，随机给出参数向量Θ = (sp, Sn, Ds，p, Ds，n, kp, kn)的初始值； 步骤2 :通过细菌趋药性行为优化参数向量θ，使目标函数减小： 首先，根据式（22)计算目标函数值。考虑菌群聚集行为，根据式（24)对菌群适应度函 数进行修正，加快参数辨识速度； 然后，参数向量Θ随机选择一个方向进行翻转行为优化，根据式（23)计算优化后的参 数向量Θ，并计算目标函数值J(0)。如果J(0)减小，则在该方向上继续前进行为优化， 直到J(9)不再减小，或参数向量Θ在该方向上达到了最大的前进行为优化次数，否则参 数向量Θ将随机选择另外一个方向进行翻转行为优化； 步骤3:参数向量Θ优化后，通过细菌复制行为对参数向量Θ进行复制与淘汰： 首先，根据式（25)计算参数向量Θ的健康函数值，并对各参数向量Θ的健康函数值 按从大到小的原则进行排序； 然后，对健康函数值低的参数向量Θ进行复制，将健康函数值高的参数向量Θ淘汰， 淘汰的参数向量Θ的个数为，以保证参数向量Θ的总数不变； 步骤4:参数向量Θ复制与淘汰后，通过细菌的消散行为将参数向量Θ按消散概率Ped重新随机分布到寻优区间，使算法逃逸出局部极值，从而求出全局最优点； 步骤5:参数向量Θ重新分布后，计算参数向量Θ对应的目标函数值，使目标函数最 小的参数向量Θ即为参数辨识结果。7. 根据权利要求1所述的一种锂离子电池机理建模方法，其特征在于： 所述锂离子电池单粒子模型的正极开路电压表达式的拟合步骤为： (1) 采用0.1 C电流对锂离子电池进行恒流放电实验，测得锂离子电池端电压Pcc,; (2) 在矩阵实验室Matlab中搭建锂离子电池单粒子模型，将正、负极开路电压的经验 公式代入到所搭建的锂离子电池单粒子模型中；采用0.1 C电流对锂离子电池单粒子模型 进行恒流放电仿真实验，得到锂离子电池单粒子模型端电压V ; (3) 当锂离子电池进行小电流充放电实验时，端电压与开路电压差近似相等，而且负极 开路电压较小，则锂离子电池单粒子模型正极开路电压通过式（28)计算： Ull(Oll)=U^dp)+ (Vocr-V) (28) 其中，Up( Θ p)是锂离子电池单粒子模型正极开路电压； WA)是正极开路电压经验值； Kct是通过恒流放电实验测得的锂离子电池端电压； V是锂离子电池单粒子模型端电压； 将拟合得到锂离子电池单粒子模型正极开路电压表达式代入到锂离子电池单粒子模 型中，得到锂离子电池单粒子模型端电压； (4) 将（3)中得到的锂离子电池单粒子模型端电压与（1)中实验测得的锂离子电池 端电压比较，若差值小于±〇.〇5V，则（3)中拟合得到的锂离子电池单粒子模型正极开路电 压表达式为所求，否则，用该表达式代替正极开路电压经验公式，返回步骤（2)，重复以上步 骤，直至端电压差小于±〇.〇5V。8. 根据权利要求1所述的一种锂离子电池机理建模方法，其特征在于： 建立锂离子电池单粒子模型的条件为： 1)假设锂离子电池电极由多个具有相同大小和动力学特性的球状粒子组成，并且电流 通过电极时在所有活性粒子内均匀分布； 2) 假设在固相颗粒内或颗粒之间电压降为零； 3) 假设在整个锂离子电池内部液相锂离子浓度恒定，并且在时间和空间上均匀分布； 4) 忽略液相电压对电池端电压的影响； 5) 忽略锂离子电池充放电过程中产生的热量。
【专利摘要】本发明属于电动汽车锂离子动力电池技术领域，涉及一种锂离子电池机理建模方法；克服了锂离子电池电化学模型结构复杂、参数难以辨识，经验模型精度低的缺点；包括以下步骤：1)建立锂离子电池单粒子模型；2)采用三参数抛物线方法简化锂离子电池单粒子模型中的固相扩散方程；3)采用菌群觅食优化算法辨识锂离子电池单粒子模型中的未知参数；4)拟合锂离子单粒子模型的正极开路电压表达式；本发明采用三参数抛物线方法，简化了锂离子电池单粒子模型的结构；采用菌群觅食优化算法辨识锂离子电池单粒子模型中的未知参数，辨识速度快，得到了全局最优解；本发明为锂离子电池状态估计，寿命预测，特性分析提供理论支持。
【IPC分类】H01M10/0525, G06F19/00
【公开号】CN104899439
【申请号】CN201510296543
【发明人】马彦, 殷明月, 李炳思, 周秀文, 应振华, 茹敬佩, 何鹏才, 孙延帅, 王留, 陈虹
【申请人】吉林大学
【公开日】2015年9月9日
【申请日】2015年6月2日
【技术领域】
[0001] 本发明属于电动汽车锂离子动力电池技术领域，更具体地说，本发明涉及一种锂 离子电池机理建模方法。
【背景技术】
[0002] 动力电池的性能对电动汽车的整车动力性、续驶里程和安全性影响很大。电池的 荷电状态和健康状态估计的准确性、寿命的预测精度、安全稳定的运行都与动力电池的建 模精度密切相关，因此，建立精确的电池模型对电池的研宄具有重要意义。
[0003] 目前，常用的电池模型有经验模型、多物理场耦合模型和电化学模型。经验模型不 考虑物理、化学反应原理，依据所采集的实验数据，拟合出锂离子电池的特性曲线，建立电 池的经验模型，该模型结构简单，但预测能力差；多物理场耦合模型同时考虑了锂离子电池 的内部电场和温度场反应和分布过程，基于电池内部电化学反应过程的生热原理和传热学 原理，分析锂离子电池内部的热特性，该模型主要用于模拟电池在正常工作状态下的内部 温度变化情况，只适用于小功率情况，使用范围有限；电化学模型是基于多孔电极原理、浓 溶液原理和欧姆定律等理论建立的，由偏微分方程、常微分方程、代数方程组成，该模型精 度高，可以较为准确的描述电池内部电化学反应机理和外部反应特性，但模型参数相互耦 合，且与电池结构、尺寸、所用材料等因素有关，模型求解过程复杂，计算量大，给实车上应 用带来困难。
[0004] 在锂离子电池模型中，由于正极和负极活性区域的表面积、正极和负极固相锂离 子扩散系数与正极和负极反应速率常数无法直接测量，且会随着锂离子电池的使用而变 化，所以需要对这些参数进行辨识。目前常用的最小二乘法和极大似然参数辨识方法，要求 目标函数连续可导，极容易使辨识结果陷入局部最小的情况，从而难以获得全局最优的辨 识结果；基于神经网络的参数辨识方法可以高精度逼近非线性函数，但是神经网络结构的 确定、样本数据的选择和神经网络训练算法等问题尚未有效解决，影响了该方法的应用；群 体智能优化算法从生物生命演化过程中得到启示，模拟生物进化行为，是一种多智能体并 行的算法，鲁棒性好。
【发明内容】
[0005] 本发明所要解决的技术问题是克服锂离子电池电化学模型结构复杂、参数难以辨 识，经验模型精度低的缺点，本发明提出了一种基于菌群觅食优化算法的锂离子电池单粒 子模型建模方法。
[0006] 为解决上述技术问题，本发明是采用如下技术方案实现的，结合【附图说明】如下：
[0007] 一种锂离子电池机理建模方法，包括以下步骤：
[0008] 步骤1 :建立锂离子电池单粒子模型；
[0009] 所述锂离子电池单粒子模型描述锂离子电池内的机理反应过程和物理特性，锂离 子电池单粒子模型包括固相锂离子扩散方程、巴特勒-伏尔摩Butler-Volmer动力学方程 和端电压表达式；
[0010] 步骤2 :采用三参数抛物线方法简化锂离子电池单粒子模型中的固相锂离子扩散 方程；
[0011] 所述锂离子电池单粒子模型中的固相锂离子扩散方程为二阶偏微分方程，采用三 参数抛物线方法将其简化为常微分方程组和代数方程；
[0012] 步骤3 :采用菌群觅食优化算法辨识锂离子电池单粒子模型中的未知参数；
[0013] 所述锂离子电池单粒子模型中的未知参数为：正极活性区域的表面积Sp，负极活 性区域的表面积Sn，正极固相锂离子扩散系数Ds，p，负极固相锂离子扩散系数Ds，n，正极反应 速率常数kp，负极反应速率常数kn;
[0014] 步骤4 :拟合锂离子单粒子模型的正极开路电压表达式。
[0015] 技术方案中所述固相锂离子扩散方程为：
[0017] 式⑴的边界条件为：
[0020] 其中，i=p，n分别代表锂离子电池的正极和负极；
[0021] y为固相锂离子浓度；
[0022] t为时间；
[0023] Ds，i为固相锂离子扩散系数；
[0024] rf(〇，Ri)为电极内球状粒子的径向距离；
[0025] 氏为电极内球状粒子的半径；
[0026] 上为固相粒子表面的锂尚子孔壁通量。
[0027] 所述固相粒子表面的锂离子孔壁通量i的定义为：
[0029] 其中，I为锂离子电池的充放电电流；
[0030] F为法拉第常数；
[0031] Si*电极内活性区域的表面积；
[0032] 所述巴特勒-伏尔摩Butler-Volmer动力学方程为：
[0034] 其中，匕为电极反应速率常数；
[0035] cs，imax为固相最大锂离子浓度；
[0036] 为固相粒子表面锂离子浓度；
[0037] c为液相锂离子浓度；
[0038] R为普适气体常数；
[0039]T为温度；
[0040] 1为电极过电压；
[0041] 求解式（5)，得电极过电压表达式为：
[0044] 由于锂离子电池单粒子模型忽略了与液相扩散相关的反应过程，液相电压为零， 电极过电压与电极固相电压、电极开路电压之间的关系为：
[0046] 其中，（i>si为电极固相电压；
[0047] 仏为电极开路电压，其表达式是关于电极固相粒子表面锂离子浓度的非线性函 数。
[0048] 锂离子电池正极固相电压与负极固相电压之间的差值为锂离子电池端电压，为了 计算方便，引入荷电状态变量9 i:
[0050] 所述锂离子电池端电压表达式为：
[0051] v= <i>s；p-<i)s；n= (up( 0p)-un( 0 n)) + (np-nn) (9)
[0052] 式（l)-(4)、（6)、（8)-(9)为以电流为输入，端电压为输出的锂离子电池单粒子模 型数学表达式。
[0053] 技术方案中所述锂离子电池单粒子模型固相锂离子扩散方程为二阶偏微分方程， 采用三参数抛物线方法将其简化为常微分方程组和代数方程；
[0054] 采用三参数抛物线方法表示固相锂离子浓度为：
[0056] 其中，a(t)、b(t)、c(t)是待求解的系数；
[0057]式（10)满足式⑵中r= 0处的边界条件，将式（10)代入式⑴和式（3)，得：
[0060] 根据式（10)，变量固相锂离子平均浓度1(0、固相粒子表面锂离子浓度和 固相锂离子体积平均浓度通量_)可以表示为：
[0064] 由式（13)_(15)，求得a(t)、b(t)和c(t)的表达式为：
[0068] 则固相粒子锂离子平均浓度满足：
[0070] 固相锂离子体积平均浓度通量满足：
[0072] 固相粒子表面锂离子浓度为：
[0074] 式（19)-(21)为锂离子电池单粒子模型简化的固相锂离子扩散方程。
[0075] 技术方案中锂离子电池单粒子模型中的未知参数向量为0 = (Sp，Sn，Ds，p,Ds，n,kp，kn) ; 0为锂离子电池单粒子模型中的未知参数向量；
[0076] 所述采用菌群觅食优化算法辨识锂离子电池单粒子模型中的未知参数过程中，目 标函数为：
[0078] 其中，yi为锂离子电池充放电实验第i个采样点的端电压；
[0079] f(Xi，0 )为锂离子电池单粒子模型第i个采样点的端电压；
[0080] XiS充放电实验第i个采样点的输入电流；
[0081] 0 = (Sp，Sn，Ds，p,Ds，n,kp，kn)为锂离子电池单粒子模型中待辨识的参数向量；
[0082] 将目标函数作为菌群觅食优化算法的适应度函数，基于实测的锂离子电池充放电 电流和端电压数据，得到使式（22)中目标函数最小的参数向量0即为参数辨识结果。
[0083] 技术方案中所述菌群觅食优化算法模拟了细菌觅食的趋药性、聚集、复制和消散 四种智能行为；
[0084] 所述趋药性行为包括翻转行为和前进行为；翻转行为是细菌沿着任意一个新的方 向运动，前进行为是细菌沿着与上一步相同的方向运动；
[0085] -次翻转行为优化后参数向量0值为：
[0086] 0i(j+l，k，1) = 0i(j，k，1)+C⑴伞⑴(23)
[0087] 其中，0W，k，1)为第j(0 <j彡Nc)步趋药性，第k(0 <k彡Nre)步复制，第1 (0 < 1 <NJ步消散行为对应的第i个参数向量0 ;
[0088] 0i(j+l，k, 1)为一次翻转行为优化后第i个参数向量0值；
[0089] N。为趋药性行为总步数；
[0090] 为复制行为总步数；
[0091] Ned为消散行为总步数；
[0092] 巾（i)为单位长度的 随机方向向量；
[0093] C(i)为参数向量优化步长；
[0094] 所述聚集行为是菌群觅食过程中，细菌个体之间通过释放引诱剂和排斥剂来完成 菌群的聚集行为；
[0095] 细菌间聚集行为的数学表达式为：
[0097] 其中，JJ0 )为每两个细菌间的适应度函数；
[0098] S为参数变量0的个数；
[0099] dattraet为引诱剂扩散深度；
[0100] ?attrac;t为引诱剂的扩散率的大小；
[0101 ]D为待辨识参数的维数；
[0102] 9 为参数向量9在第m维空间的分量；
[0103] 化为第i个参数向量9在m维空间的分量；
[0104] h,epellant为排斥剂扩散高度；
[0105] ?Mpellant为排斥剂扩散率的大小；
[0106] 所述复制行为是一个周期的趋药性行为后，根据细菌健康函数值的大小进行复制 与淘汰，健康函数值较小的细菌开始复制行为，生成与父代完全相同的子代，没有复制的细 菌将被淘汰，复制细菌的个数与淘汰细菌的个数相等，以此来维持菌群的个数不变；
[0107] 健康函数为趋药性行为中参数向量0的适应度函数值的和：
[0109] 其中，表示第i个参数向量0的健康函数，，越大表示参数向量0优化效 果越差；
[0110] J(i，j，k，l)表示对应第i个参数向量0的适应度函数值；
[0111] 所述消散行为是指细菌被外力杀死或者被驱散到新的区域中的行为；消散行为破 坏了细菌的趋药性过程，但是，这使得细菌可能寻找到食物更加丰富的区域。
[0112] 技术方案中所述菌群觅食优化算法辨识参数的过程为：
[0113] 步骤1 :在参数优化范围内，随机给出参数向量0 = (Sp,Sn,Ds，p,Ds，n,kp,kn)的初 始值；
[0114] 步骤2 :通过细菌趋药性行为优化参数向量0，使目标函数减小：
[0115] 首先，根据式（22)计算目标函数值。考虑菌群聚集行为，根据式（24)对菌群适应 度函数进行修正，加快参数辨识速度；
[0116] 然后，参数向量0随机选择一个方向进行翻转行为优化，根据式（23)计算优化后 的参数向量0，并计算目标函数值J(0)。如果J(0)减小，则在该方向上继续前进行为优 化，直到J(0)不再减小，或参数向量0在该方向上达到了最大的前进行为优化次数，否则 参数向量0将随机选择另外一个方向进行翻转行为优化；
[0117] 步骤3 :参数向量0优化后，通过细菌复制行为对参数向量0进行复制与淘汰：
[0118] 首先，根据式（25)计算参数向量0的健康函数值，并对各参数向量0的健康函 数值按从大到小的原则进行排序；
[0119] 然后，对健康函数值低的参数向量0进行复制，将健康函数值高的参数向量0淘 汰，淘汰的参数向量0的个数为*5； =|，以保证参数向量0的总数不变；
[0120] 步骤4:参数向量0复制与淘汰后，通过细菌的消散行为将参数向量0按消散概 率ped重新随机分布到寻优区间，使算法逃逸出局部极值，从而求出全局最优点；
[0121] 步骤5 :参数向量0重新分布后，计算参数向量0对应的目标函数值，使目标函 数最小的参数向量0即为参数辨识结果。
[0122] 技术方案中所述锂离子电池单粒子模型的正极开路电压表达式的拟合步骤为：
[0123] (1)采用0. 1C电流对锂离子电池进行恒流放电实验，测得锂离子电池端电压
[0124] (2)在矩阵实验室Matlab中搭建锂离子电池单粒子模型，将正、负极开路电压的 经验公式代入到所搭建的锂离子电池单粒子模型中；采用0. 1C电流对锂离子电池单粒子 模型进行恒流放电仿真实验，得到锂离子电池单粒子模型端电压V;
[0125] (3)当锂离子电池进行小电流充放电实验时，端电压与开路电压差近似相等，而且 负极开路电压较小，则锂离子电池单粒子模型正极开路电压通过式（28)计算：
[0127] 其中，Up( 0p)是锂离子电池单粒子模型正极开路电压；
[0128] 是正极开路电压经验值；
[0129] 是通过恒流放电实验测得的锂离子电池端电压；
[0130] V是锂离子电池单粒子模型端电压；
[0131] 将拟合得到锂离子电池单粒子模型正极开路电压表达式代入到锂离子电池单粒 子模型中，得到锂离子电池单粒子模型端电压；
[0132] (4)将（3)中得到的锂离子电池单粒子模型端电压与（1)中实验测得的锂离子电 池端电压比较，若差值小于±〇.〇5V，则（3)中拟合得到的锂离子电池单粒子模型正极开路 电压表达式为所求，否则，用该表达式代替正极开路电压经验公式，返回步骤（2)，重复以上 步骤，直至端电压差小于±〇.〇5V。
[0133] 技术方案中建立锂离子电池单粒子模型的条件为：
[0134] 1)假设锂离子电池电极由多个具有相同大小和动力学特性的球状粒子组成，并且 电流通过电极时在所有活性粒子内均匀分布；
[0135] 2)假设在固相颗粒内或颗粒之间电压降为零；
[0136] 3)假设在整个锂离子电池内部液相锂离子浓度恒定，并且在时间和空间上均匀分 布；
[0137]4)忽略液相电压对电池端电压的影响；
[0138] 5)忽略锂离子电池充放电过程中产生的热量。
[0139] 与现有技术相比本发明的有益效果是：
[0140] (1)锂离子电池单粒子模型忽略液相电压对端电压的影响，采用三参数抛物线方 法将锂离子电池单粒子模型中表示固相扩散方程的偏微分方程组简化为常微分方程组和 代数方程，简化了锂离子电池单粒子模型的结构。
[0141] (2)采用菌群觅食优化算法辨识锂离子电池单粒子模型中的未知参数，辨识速度 快，得到了全局最优解。
[0142] (3)锂离子电池单粒子模型能够较准确地模拟锂离子电池的放电行为，锂离子电 池单粒子模型端电压数据与锂离子电池端电压数据可以很好地吻合，最大误差在±〇.〇7V 左右。此外，锂离子电池单粒子模型可以为锂离子电池状态估计，寿命预测，特性分析提供 理论支持。
【附图说明】
[0143] 下面结合附图对本发明作进一步的说明：
[0144] 图1为锂离子电池单粒子模型建模流程图；
[0145] 图2为锂离子电池单粒子模型示意图；
[0146] 图3为锂离子电池单粒子模型仿真模型结构框图；
[0147] 图4为菌群觅食优化算法辨识锂离子电池单粒子模型参数流程图；
[0148] 图5为适应度函数随趋药性行为步数变化曲线；
[0149] 图6为磷酸铁锂电池正极开路电压与正极荷电状态变量之间的关系曲线；
[0150] 图7为磷酸铁锂电池负极开路电压与负极荷电状态变量之间的关系曲线；
[0151] 图8为1. 6Ah磷酸铁锂电池端电压曲线与单粒子模型端电压曲线；
[0152] 图9为1. 6Ah磷酸铁锂电池单粒子模型端电压误差曲线。
【具体实施方式】
[0153] 下面结合附图对本发明作详细的描述：
[0154] 锂离子电池单粒子模型建模流程图参阅图1所示，本发明基于菌群觅食优化算法 的锂离子电池单粒子模型建模步骤如下：
[0155] 1、根据锂离子电池工作原理，建立锂离子电池单粒子模型。
[0156] 锂离子电池主要由正极、负极和隔膜组成，正极活性物质为富锂金属氧化物与相 应金氧化物的混合物质，负极活性物质是石墨或近似石墨结构的碳，隔膜是一种经特殊成 型的高分子薄膜，有微孔结构，锂离子可以自由通过，而电子不能通过。锂离子电池充电 时，正极上有锂离子生成，生成的锂离子通过电池隔膜运动到电池负极，而作为负极的碳呈 现层状结构，它有很多微孔，运动到负极的锂离子就嵌入到碳层的微孔中，嵌入的锂离子越 多，充电量越高。当对锂离子电池进行放电时，嵌在负极碳层中的锂离子析出回到正极，回 到正极的锂离子越多，放电量越高。
[0157] 锂离子电池单粒子模型是指利用一个球状粒子特性来代表整个电极特性而建立 的一种简化的锂离子电池机理模型，锂离子电池单粒子模型示意图如图2所示。建立锂离 子电池单粒子模型的条件为：
[0158] (1)假设锂离子电池电极由多个具有相同大小和动力学特性的球状粒子组成，并 且电流通过电极时在所有活性粒子内均匀分布；
[0159] (2)假设在固相颗粒内或颗粒之间电压降为零；
[0160] (3)假设在整个锂离子电池内部液相锂离子浓度恒定，并且在时间和空间上均匀 分布；
[0161] (4)忽略液相电压对电池端电压的影响；
[0162] (5)忽略锂离子电池 充放电过程中产生的热量。
[0163] 依据上述假设条件，根据锂离子电池内的机理反应过程和物理特性，锂离子电池 单粒子模型包括固相扩散方程，巴特勒-伏尔摩（Butler-Volmer)动力学方程和端电压表 达式。
[0164] 锂离子在电极固相中的扩散为非稳态扩散，根据菲克第二定律，在球形电极上，固 相锂离子扩散方程为：
[0166] 式⑴的边界条件为：
[0169] 其中，i=p，n分别代表锂离子电池的正极和负极；为固相锂离子浓度；t为时 间；Ds;i为固相锂离子扩散系数；r# (〇，Ri)为电极内球状粒子的径向距离况为电极内球 状粒子的半径；1为固相粒子表面的锂离子孔壁通量。
[0170] 固相粒子表面的锂离子孔壁通量1的定义为：
[0172] 其中，I为锂离子电池的充放电电流；F为法拉第常数A为电极内活性区域的表 面积。巴特勒-伏尔摩（Butler-Volmer)动力学方程描述固相活性粒子表面与电解液溶液 临界面处的电化学反应过程，该动力学方程为：
[0174] 其中，&为电极反应速率常数；Cs，imax为固相最大锂离子浓度；为固相粒子表面 锂离子浓度；C为液相锂离子浓度；R为普适气体常数；T为温度；n电极过电压。
[0175]求解式（5)，得电极过电压表达式为：
[0178] 由于锂离子电池单粒子模型忽略了与液相扩散相关的反应过程，液相电压为零， 电极过电压与电极固相电压、电极开路电压之间的关系为：
[0180] 其中，为电极固相电压；为电极开路电压，是关于电极固相粒子表面锂 离子浓度的非线性函数。
[0181] 锂离子电池正极固相电压与负极固相电压之间的差值为锂离子电池端电压，为了 计算方便，引入荷电状态变量9 i:
[0183] 锂离子电池端电压表达式为：
[0184] v= <i>s；p-<i)s；n= (up( 0p)-un( 0 n)) + (np-nn) (9)
[0185] 式（l)-(4)、（6)、（8)-(9)即为以电流为输入，端电压为输出的锂离子电池单粒子 模型数学表达式，锂离子电池单粒子模型仿真模型结构框图参阅图3所示。
[0186] 2、采用三参数抛物线方法简化锂离子电池单粒子模型中的固相扩散方程。
[0187] 锂离子电池单粒子模型固相扩散方程为二阶偏微分方程，采用三参数抛物线方法 将其简化为常微分方程组和代数方程。采用三参数抛物线方法表示固相锂离子浓度为：
[0189] 其中，a(t)、b(t)、c(t)是待求解的系数。
[0190] 显然，式（10)满足式⑵中r= 0处的边界条件，将式（10)代入式⑴和式（3)， 得：
[0193] 根据式（10)，变量固相锂离子平均浓度己(0、固相粒子表面锂离子浓度和固 相锂离子体积平均浓度通量砂)可以表示为：
[0197] 由式（13)-(15)，可求得3(〇、13(〇、(^)的表达式为：
[0201] 结合式（1)-⑶求解，得：
[0205] 式（19)-(21)即为锂离子电池单粒子模型简化的固相扩散方程。
[0206] 3、采用菌群觅食优化算法辨识锂离子电池单粒子模型中的未知参数。
[0207] 锂离子电池单粒子模型中的正极和负极活性区域的表面积、正极和负极固相锂 离子扩散系数、正极和负极反应速率常数无法直接测量，并且会随着锂离子电池的使用 而变化，所以需要对这些参数进行辨识。本发明采用菌群觅食优化算法辨识参数0 = (Sp，Sn，Ds，p,Ds，n,kp，kn) 〇
[0208] 细菌在觅食时希望在最短的时间内获得最多的营养。为了达到这个目标，细菌 要决定觅食策略，通过感应自身周围的化学物质浓度，寻找食物所在的区域，找到后决 定进入该区域还是继续寻找，进入后什么时候离开该区域。通过模拟细菌的觅食策略， Passino在文献Biomimicryofbacterialforagingfordistributedoptimizationand control(ControlSystems,IEEE, 2002, 22(3) : 52-67.)中提出 了菌群觅食优化算法，算法 流程图参阅图4所示。
[0209] 在菌群觅食优化算法辨识锂离子电池单粒子模型未知参数过程中，目标函数为：
[0211] 其中，yi为锂离子电池充放电实验第i个采样点的端电压；f(Xi，0)为锂离子电 池单粒子模型第i个采样点的端电压；^为充放电实验第i个采样点的输入电流；0 = (Sp，Sn，Ds，p,Ds，n,kp，kn)为锂离子电池单粒子模型中待辨识的参数向量。
[0212] 将目标函数作为菌群觅食优化算法的适应度函数，基于实测的锂离子电池端电压 和充放电电流数据，采用菌群觅食优化算法，得到使式（22)中目标函数最小的参数向量0 即为参数辨识结果。菌群觅食优化算法模拟了细菌觅食的趋药性、聚集、复制和消散四种智 能行为。
[0213] ⑴趋药性行为
[0214] 由于细菌经常生存在化学引诱剂环境中，它对环境的应激反应称为趋药性行为。 细菌的基本趋药性行为包括：翻转行为和前进行为。翻转行为是细菌沿着任意一个新的方 向运动，而前进行为是细菌沿着与上一步相同的方向运动。
[0215] -次翻转行为优化后参数向量0值为：
[0216] 0i(j+l，k，1) = 0i(j，k，1)+C⑴伞⑴(23)
[0217] 其中，0i(j，k，l)为第j(0<j<N。）步趋药性，第k(0<k<Nj步复制，第1(0 < 1 <NJ步消散行为对应的第i个参数向量0 ; 0 1 (j+1,k, 1)为一次翻转行为优化后第 i个参数向量9 ;N。为趋药性行为总步数；NM为复制行为总步数；Nm为消散行为总步数； (i)为单位长度的随机方向向量；C(i)为参数向量优化步长。
[0218] ⑵聚集行为
[0219] 在菌群觅食过程中，细菌个体之间通过释放引诱剂和排斥剂来完成菌群的聚集行 为。弓丨诱剂使细菌聚集在一起，甚至出现"抱团"现象；排斥剂使每个细菌都有一定的位置， 使其能在该位置上获取营养，以维持生存。
[0220] 细菌间聚集行为的数学表达式为：
[0222] 其中，上。（0 )为每两个细菌间的适应度函数；S为参数变量0的个数；dattMet为 引诱剂扩散深度，即释放引诱剂的量；《attM;t为引诱剂的扩散率的大小；D为寻优空间的维 数，即待辨识参数的维数；0m为参数向量0在第m维空间的分量；％为第i个参数向量0 在m维空间的分量；hMIrellant为排斥剂扩散高度，即释放排斥剂的量；《 MIrellant为排斥剂扩散 率的大小。
[0223] ⑶复制行为
[0224] 一个周期的趋药性行为后，根据细菌健康函数值的大小进行复制与淘汰，健康函 数值较小的细菌开始复制行为，生成与父代完全相同的子代，没有复制的细菌将被淘汰，复 制细菌的个数与淘汰细菌的个数相等，以此来维持细菌的个数不变。
[0225] 健康函数为趋药性行为中参数向量0的适应度函数值的和：
[0227] 其中，/Lm表示第i个参数向量9的健康函数，越大表示参数向量9优化效 果越差；J(i,j,k, 1)表示对应第i个参数向量0的适应度函数值。
[0228] ⑷消散行为
[0229] 消散行为是指细菌被外力杀死或者被驱散到新的区域中的行为。消散行为破坏了 细菌的趋药性过程，但是，这使得细菌可能寻找到食物更加丰富的区域。
[0230] 菌群觅食优化算法辨识参数的过程可以归纳为：
[0231] 步骤1在参数优化范围内，随机给出参数向量0 = (Sp,Sn,Ds，p,Ds，n,kp,kn)的初始 值；
[0232] 步骤2通过细菌趋药性行为优化参数向量0，使目标函数减小：
[0233] 首先，根据式（22)计算目标函数值。考虑菌群聚集行为，根据式（24)对菌群适应 度函数进行修正，加快参数辨识速度；
[0234] 然后，参数向量0随机选择一个方向进行翻转行为优化，根据式（23)计算优化后 的参数向量0，并计算目标函数值J(0)。如果J(0)减小，则在该方向上继续前进行为优 化，直到J(0)不再减小，或参数向量0在该方向上达到了最大的前进行为优化次数，否则 参数向量0将随机选择另外一个方向进行翻转行为优化；
[0235] 步骤3参数向量0优化后，通过细菌复制行为对参数向量0进行复制与淘汰：
[0236] 首先，根据式（25)计算参数向量0健康函数值，并对各参数向量0的健康函数 值按从大到小的原则进行排序；
[0237]然后，对健康函数值低的参数向量0进行复制，将健康函数值高的参数向量0淘 汰，淘汰的参数向量Q个数为，以保证参数向量Q的总数不变；
[0238] 步骤4参数向量0复制与淘汰后，通过细菌的消散行为将参数向量0按消散概 率ped重新随机分布到寻优区间，使算法逃逸出局部极值，从而求出全局最优点；
[0239] 步骤5参数向量0重新分布后，计算参数向量0对应的目标函数值，使目标函数 最小的参数向量0即为参数辨识结果。
[0240] 4、拟合锂离子电池单粒子模型的正极开路电压表达式。
[0241] 锂离子电池负极的活性物质为石墨或近似石墨结构的碳，结构稳定，其开路电压 与负极荷电状态变量之间的关系可以用经验公式表示。锂离子电池正极的活性物质为富锂 金属氧化物与相应金属氧化物的混合物质，并且因为锂离子电池的正极开路电压决定了电 池的开路电压差，所以需要拟合锂离子电池单粒子模型中的正极开路电压表达式。正极开 路电压Up(0p)是正极荷电状态变量0p的非线性函数，本发明通过锂离子电池小电流恒流 放电实验，基于实验测得的锂离子电池端电压与锂离子电池单粒子模型端电压之间的差， 结合正极开路电压经验公式，拟合锂离子电池单粒子模型正极开路电压表达式。
[0242] 当锂离子电池进行小电流充放电实验时，端电压与开路电压差近似相等，而且负 极开路电压较小，则锂离子电池单粒子模型正极开路电压可以通过式（26)计算：
[0244] 其中，Up(0p)是锂离子电池单粒子模型正极开路电压；是正极开路电压经 验值；la是通过恒流放电实验测得的锂离子电池端电压；V是锂离子电池单粒子模型端 电压。
[0245] 以磷酸铁锂电池为例，根据文献Modelingofacommercialgraphite/ LiFeP04cell(JournaloftheElectrochemicalSociety,2011， 158(5):A562_A571)，正、 负极开路电压经验公式为：
[0248] 其中，W吹）为负极开路电压经验值。
[0249] 锂离子电池单粒子模型正极开路电压表达式具体拟合过程为：
[0250] (1)采用0. 1C电流对锂离子电池进行恒流放电实验，测得锂离子电池端电压。
[0251] (2)在矩阵实验室（Matlab)中搭建锂离子电池单粒子模型，将正、负极开路电压 的经验公式代入到所搭建的锂离子电池单粒子模型中。采用〇. 1C电流对锂离子电池单粒 子模型进行恒流放电仿真实验，得到锂离子电池单粒子模型端电压V。
[0252] (3)根据式（26)，拟合得到锂离子电池单粒子模型正极开路电压表达式，并将其 代入到锂离子电池单粒子模型中，得到锂离子电池单粒子模型端电压。
[0253] (4)将（3)中得到的锂离子电池单粒子模型端电压与⑴中实验测得的锂离子电 池端电压比较，若差值小于±〇.〇5V，则（3)中拟合得到的锂离子电池单粒子模型正极开路 电压表达式即为所求，否则，用该表达式代替正极开路电压经验公式，返回步骤（2)，重复以 上步骤，直至端电压差小于±〇.〇5V。
[0254] 下面以一个具体实例解释说明本发明。
[0255] 本实例以正极为LiFeP04，负极为1^(：6的1.6Ah锂离子电池为研宄对象，具体实施 过程如下：
[0256] 1、根据锂离子电池内部电化学反应机理，得到锂离子电池单粒子模型的数学表达 式为式（1)-(4)、（6)、（8)-(9)。
[0257] 2、采用三参数抛物线方法化简锂离子电池单粒子模型，得到化简后的锂离子电池 单粒子模型数学表达式为式（6)、（8) _(9)、（19)-(21)。
[0258] 3、采用菌群觅食优化算法辨识锂离子电池单粒子模型中的未知参数0 = (Sp，Sn，Ds，p,Ds，n,kp，kn) 〇
[0259] 首先对磷酸铁锂电池进行0. 1C恒流放电实验得到电池端电压数据。
[0260] 然后确定锂离子电池单粒子模型中待辨识参数0 = (Sp，Sn，Ds,p，Ds,n，kp，kn)的范 围：SpG[0? 1，l]、SnG[0? 1，l]、Ds，pG[0? 5X1(T14,lXl(T13]、Ds，nG[0? 5X1(T16, 1XKT15]、 kpG[0? 1X10 -11，1X1(T10]、knG[0? 1X10 -12, 1XKT11] 〇
[0261] 其次参数初始化：设辨识参数初始值为61,;,(〇)，待辨识参数向量0的维数D= 6 ;根 据文献"菌群优化算法的研宄"（哈尔滨工业大学，2009.)取参数向量0的个数S= 20， 趋药性行为总步数凡=60,趋药性行为中前进行为优化最大步数NS= 4,复制行为总步数 Nra= 8,消散行为总步数Ned= 2,复制参数向量0组数
消散概率口《1= 〇. 25,引诱 剂扩散深度dattrart= 0.05,引诱剂的扩散率的大小《 attrart= 0. 1，排斥剂扩散高度h_llant = 0.05,排斥剂的扩散率的大小 〇repellant4〇
[0262] 最后根据图4菌群觅食优化算法辨识锂离子电池单粒子模型参数流程，结合实验 数据进行参数辨识，辨识结果为：
[0263] Sp= 0? 4945m2，Sn= 0? 4937m2，Ds，p= 8. 13Xl(T14m2 .s'Ds，n= 6. 64Xl(T16m2 .s' kp= 4. 25Xl〇-nm2.5 ?m〇r0.5 ?s'kn= 2. 09X10_12m2.5 ?m〇r0.5 ?s-1。
[0264] 参数向量适应度函数随趋药性行为步数变化曲线如图5所示，适应度函数值随着 趋药性行为步数的增加不断减小最后趋于稳定。表1为磷酸铁锂电池单粒子模型参数值。
[0265] 表 1
[0266]
[0267] 4、拟合锂离子电池单粒子模型的正极开路电压表达式。
[0268] 经过多次循环迭代和仿真实验，得到磷酸铁锂电池单粒子模型的正极开路电压表 达式为：
[0270] 磷酸铁锂电池单粒子模型正极开路电压与正极荷电状态之间的关系曲线参阅图6 所示，磷酸铁锂电池单粒子模型负极开路电压与负极荷电状态之间的关系参阅图7所示。
[0271] 用上面得到的磷酸铁锂电池单粒子模型模拟电池0. 1C恒流放电行为，端电压曲 线参阅图8所示，误差曲线参阅图9所示。参阅图8，在放电过程中，当t< 10000s时，磷酸 铁锂电池单粒子模型可以准确地模拟电池放电的压降行为，端电压误差在±0.01V以内。 在放电截止时，磷酸铁锂电池正极锂离子浓度增大，根据式（8)，正极荷电状态随之增大，根 据式（29)，正极荷电状态对锂离子电池正极开路电压的影响为指数形式，磷酸铁锂电池正 极开路电压会产生急速下降的趋势。因为正极开路电压决定电池端电压，所以放电截止时 端电压产生陡降行为，磷酸铁锂电池单粒子模型端电压误差增大，但整体上，磷酸铁锂电池 单粒子模型的放电端电压可以与磷酸铁锂电池实际放电端电压的采样点较好地吻合。
【主权项】
1. 一种锂离子电池机理建模方法，其特征在于，包括以下步骤： 步骤1 :建立锂离子电池单粒子模型； 所述锂离子电池单粒子模型描述锂离子电池内的机理反应过程和物理特性，锂离子电 池单粒子模型包括固相锂离子扩散方程、巴特勒-伏尔摩Butler-Volmer动力学方程和端 电压表达式； 步骤2:采用三参数抛物线方法简化锂离子电池单粒子模型中的固相锂离子扩散方 程； 所述锂离子电池单粒子模型中的固相锂离子扩散方程为二阶偏微分方程，采用三参数 抛物线方法将其简化为常微分方程组和代数方程； 步骤3 :采用菌群觅食优化算法辨识锂离子电池单粒子模型中的未知参数； 所述锂离子电池单粒子模型中的未知参数为：正极活性区域的表面积Sp，负极活性区 域的表面积Sn，正极固相锂离子扩散系数Ds，p，负极固相锂离子扩散系数D s，n，正极反应速率 常数kp，负极反应速率常数kn; 步骤4 :拟合锂离子单粒子模型的正极开路电压表达式。2. 根据权利要求1所述的一种锂离子电池机理建模方法，其特征在于： 所述固相锂离子扩散方程为：(1) 式（1)的边界条件为：(2) (3) 其中，i = ρ，η分别代表锂离子电池的正极和负极； Cy为固相锂离子浓度； t为时间； Dsii为固相锂离子扩散系数； rie (OJi)为电极内球状粒子的径向距离； Ri为电极内球状粒子的半径； 1为固相粒子表面的锂离子孔壁通量。 所述固相粒子表面的锂离子孔壁通量Ji的定义为：(4) 其中，I为锂离子电池的充放电电流； F为法拉第常数； S#电极内活性区域的表面积； 所述巴特勒-伏尔摩Butler-Volmer动力学方程为：其中，1^为电极反应速率常数；cEUmax为固相最大锂离子浓度； 固相粒子表面锂离子浓度； c为液相锂离子浓度； R为普适气体常数； T为温度； n i为电极过电压； 求解式（5)，得电极过电压表达式为：由于锂离 子电池单粒子模型忽略了与液相扩散相关的反应过程，液相电压为零，电极 过电压与电极固相电压、电极开路电压之间的关系为： >ι=Φ,J-UiKn (7) 其中，为电极固相电压； 仏为电极开路电压，其表达式是关于电极固相粒子表面锂离子浓度的非线性函数； 锂离子电池正极固相电压与负极固相电压之间的差值为锂离子电池端电压，为了计算 方便，引入荷电状态变量9i:(8) 所述锂离子电池端电压表达式为： ν= Φ3,Ρ_Φ3,η= (υρ( θΡ)-υη( θη)) + (ηρ-ηη) (9) 式（1)-(4)、（6)、（8)-(9)为以电流为输入，端电压为输出的锂离子电池单粒子模型数 学表达式。3.根据权利要求1所述的一种锂离子电池机理建模方法，其特征在于： 采用三参数抛物线方法表示固相锂离子浓度为：(10) 其中，a (t)、b (t)、c (t)是待求解的系数； 式（10)满足式⑵中r = 0处的边界条件，将式（10)代入式⑴和式（3)，得：根据式（10)，变量固相锂离子平均浓度己(0、固相粒子表面锂离子浓度和固相 锂离子体积平均浓度通量获t)可以表示为：由式（13)-(15)，求得a(t)、b(t)和c(t)的表达式为：则固相粒子锂离子平均浓度满足：(19) 固相锂离子体积平均浓度通量满足：(20) 固相粒子表面锂离子浓度为：(21) 式（19)-(21)为锂离子电池单粒子模型简化的固相锂离子扩散方程。4. 根据权利要求1所述的一种锂离子电池机理建模方法，其特征在于： 锂离子电池单粒子模型中的未知参数向量为Θ = (Sp，Sn，Ds，p，Ds， n，kp，kn) ; Θ为锂离子 电池单粒子模型中的未知参数向量； 所述采用菌群觅食优化算法辨识锂离子电池单粒子模型中的未知参数过程中，目标函 数为：(22) 其中，Yi为锂离子电池充放电实验第i个采样点的端电压； f(Xi，Θ )为锂离子电池单粒子模型第i个采样点的端电压； Xi为充放电实验第i个采样点的输入电流； 将目标函数作为菌群觅食优化算法的适应度函数，基于实测的锂离子电池充放电电流 和端电压数据，得到使式（22)中目标函数最小的参数向量Θ即为参数辨识结果。5. 根据权利要求4所述的一种锂离子电池机理建模方法，其特征在于： 所述菌群觅食优化算法模拟了细菌觅食的趋药性、聚集、复制和消散四种智能行为； 所述趋药性行为包括翻转行为和前进行为；翻转行为是细菌沿着任意一个新的方向运 动，前进行为是细菌沿着与上一步相同的方向运动； 一次翻转行为优化后参数向量Θ值为： Θ 1 (j+1，k，1) = Θ i (j，k，I) +C ⑴ Φ ⑴(23) 其中，9\」，1^，1)为第」_(0<_]_彡凡）步趋药性，第1^0<1^彡〇步复制，第1(0 < I < NJ步消散行为对应的第i个参数向量θ ; Θ 1 (j+1, k, 1)为一次翻转行为优化后第i个参数向量Θ值； N。为趋药性行为总步数； Nm为复制行为总步数； Ned为消散行为总步数； Φ (i)为单位长度的随机方向向量； C(i)为参数向量优化步长； 所述聚集行为是菌群觅食过程中，细菌个体之间通过释放引诱剂和排斥剂来完成菌群 的聚集行为； 细菌间聚集行为的数学表达式为：其中，U Θ )为每两个细菌间的适应度函数； S为参数变量Θ的个数； dattrac；t为引诱剂扩散深度； Wattaart为引诱剂的扩散率的大小； D为待辨识参数的维数； 为参数向量Θ在第m维空间的分量； <为第i个参数向量Θ在m维空间的分量； Iwllant为排斥剂扩散高度； ω ^llant为排斥剂扩散率的大小； 所述复制行为是一个周期的趋药性行为后，根据细菌健康函数值的大小进行复制与淘 汰，健康函数值较小的细菌开始复制行为，生成与父代完全相同的子代，没有复制的细菌将 被淘汰，复制细菌的个数与淘汰细菌的个数相等，以此来维持菌群的个数不变； 健康函数为趋药性行为中参数向量Θ的适应度函数值的和：(25) 其中，7LflM表示第i个参数向量Θ的健康函数，越大表示参数向量Θ优化效果越 差； J(i，j，k, 1)表示对应第i个参数向量Θ的适应度函数值； 所述消散行为是指细菌被外力杀死或者被驱散到新的区域中的行为；消散行为破坏了 细菌的趋药性过程，但是，这使得细菌可能寻找到食物更加丰富的区域。6.根据权利要求4所述的一种锂离子电池机理建模方法，其特征在于： 所述菌群觅食优化算法辨识参数的过程为： 步骤1 :在参数优化范围内，随机给出参数向量Θ = (sp, Sn, Ds，p, Ds，n, kp, kn)的初始值； 步骤2 :通过细菌趋药性行为优化参数向量θ，使目标函数减小： 首先，根据式（22)计算目标函数值。考虑菌群聚集行为，根据式（24)对菌群适应度函 数进行修正，加快参数辨识速度； 然后，参数向量Θ随机选择一个方向进行翻转行为优化，根据式（23)计算优化后的参 数向量Θ，并计算目标函数值J(0)。如果J(0)减小，则在该方向上继续前进行为优化， 直到J(9)不再减小，或参数向量Θ在该方向上达到了最大的前进行为优化次数，否则参 数向量Θ将随机选择另外一个方向进行翻转行为优化； 步骤3:参数向量Θ优化后，通过细菌复制行为对参数向量Θ进行复制与淘汰： 首先，根据式（25)计算参数向量Θ的健康函数值，并对各参数向量Θ的健康函数值 按从大到小的原则进行排序； 然后，对健康函数值低的参数向量Θ进行复制，将健康函数值高的参数向量Θ淘汰， 淘汰的参数向量Θ的个数为，以保证参数向量Θ的总数不变； 步骤4:参数向量Θ复制与淘汰后，通过细菌的消散行为将参数向量Θ按消散概率Ped重新随机分布到寻优区间，使算法逃逸出局部极值，从而求出全局最优点； 步骤5:参数向量Θ重新分布后，计算参数向量Θ对应的目标函数值，使目标函数最 小的参数向量Θ即为参数辨识结果。7. 根据权利要求1所述的一种锂离子电池机理建模方法，其特征在于： 所述锂离子电池单粒子模型的正极开路电压表达式的拟合步骤为： (1) 采用0.1 C电流对锂离子电池进行恒流放电实验，测得锂离子电池端电压Pcc,; (2) 在矩阵实验室Matlab中搭建锂离子电池单粒子模型，将正、负极开路电压的经验 公式代入到所搭建的锂离子电池单粒子模型中；采用0.1 C电流对锂离子电池单粒子模型 进行恒流放电仿真实验，得到锂离子电池单粒子模型端电压V ; (3) 当锂离子电池进行小电流充放电实验时，端电压与开路电压差近似相等，而且负极 开路电压较小，则锂离子电池单粒子模型正极开路电压通过式（28)计算： Ull(Oll)=U^dp)+ (Vocr-V) (28) 其中，Up( Θ p)是锂离子电池单粒子模型正极开路电压； WA)是正极开路电压经验值； Kct是通过恒流放电实验测得的锂离子电池端电压； V是锂离子电池单粒子模型端电压； 将拟合得到锂离子电池单粒子模型正极开路电压表达式代入到锂离子电池单粒子模 型中，得到锂离子电池单粒子模型端电压； (4) 将（3)中得到的锂离子电池单粒子模型端电压与（1)中实验测得的锂离子电池 端电压比较，若差值小于±〇.〇5V，则（3)中拟合得到的锂离子电池单粒子模型正极开路电 压表达式为所求，否则，用该表达式代替正极开路电压经验公式，返回步骤（2)，重复以上步 骤，直至端电压差小于±〇.〇5V。8. 根据权利要求1所述的一种锂离子电池机理建模方法，其特征在于： 建立锂离子电池单粒子模型的条件为： 1)假设锂离子电池电极由多个具有相同大小和动力学特性的球状粒子组成，并且电流 通过电极时在所有活性粒子内均匀分布； 2) 假设在固相颗粒内或颗粒之间电压降为零； 3) 假设在整个锂离子电池内部液相锂离子浓度恒定，并且在时间和空间上均匀分布； 4) 忽略液相电压对电池端电压的影响； 5) 忽略锂离子电池充放电过程中产生的热量。
【专利摘要】本发明属于电动汽车锂离子动力电池技术领域，涉及一种锂离子电池机理建模方法；克服了锂离子电池电化学模型结构复杂、参数难以辨识，经验模型精度低的缺点；包括以下步骤：1)建立锂离子电池单粒子模型；2)采用三参数抛物线方法简化锂离子电池单粒子模型中的固相扩散方程；3)采用菌群觅食优化算法辨识锂离子电池单粒子模型中的未知参数；4)拟合锂离子单粒子模型的正极开路电压表达式；本发明采用三参数抛物线方法，简化了锂离子电池单粒子模型的结构；采用菌群觅食优化算法辨识锂离子电池单粒子模型中的未知参数，辨识速度快，得到了全局最优解；本发明为锂离子电池状态估计，寿命预测，特性分析提供理论支持。
【IPC分类】H01M10/0525, G06F19/00
【公开号】CN104899439
【申请号】CN201510296543
【发明人】马彦, 殷明月, 李炳思, 周秀文, 应振华, 茹敬佩, 何鹏才, 孙延帅, 王留, 陈虹
【申请人】吉林大学
【公开日】2015年9月9日
【申请日】2015年6月2日
最新回复(0)