电池离线参数辨识

一、电池模型确定

    综合模型精度与工程应用复杂度分析，采用二阶RC电路模型，能在较小的计算量下得到更好的精度。等效电路如图1所示。

图1、二阶RC等效电路模型

二、HPPC测试

    在上一篇已经详细介绍过电池充放电测试的具体流程，通过HPPC实验获得电池的电压工作曲线如图2所示，电流工作曲线如图3所示。获得的HPPC试验数据用于磷酸铁锂电池二阶Thevenin等效电路模型的参数离线辨识。HPPC试验即在电池不同SOC处进行脉冲充放电试验，放电瞬间对应的电压下降值可以被用来辨识电池内部的欧姆内阻，脉冲中的短暂放电过程对应电工学上容阻环的零状态响应，放电后的静置对应电工学上容阻环的零输入响应。

图2、HPPC电压曲线

图3、HPPC电流曲线

    摘取图二电压曲线其中一部分电压响应曲线，如图4所示

图4、脉冲试验电压响应曲线

    图4中A点为加载放电电流的时刻，C点为放电结束的采样时刻，E点为端电压趋于稳定的时刻。

    AB端的电压降用来计算等效电路中的欧姆内阻R0，如式（1）所示。

                                                            （1）

    D点到E点是电池极化反应导致电压缓慢上升，对D点到E点进行指数拟合就可以求出R1、R2、C1、C2的值;D点到E点的电压响应为零输入响应，起电压表达式如式（2）所示。

                             （2）

    SOC为90%时使用matlab对D点到E点进行指数拟合，如图5所示。

图5、指数拟合

    其中函数表达式如式（3）所示

                             （3）

    式（3）对应式（2），可得出

                                                                     

    通过MATLAB拟合工具，拟合结果如图6所示：

图6、拟合结果

    其中a为U0是当前开路电压值，由实验数据可得出a=3.3258；通过公式计算可以得出电池的电容电阻；公式如图7所示

图9、计算公式

    上述求出电池SOC为90%时的参数，更改不同电压段得出其他SOC值时的参数；参数辨识结果如表1。

表1、参数辨识结果

soc	90%	80%	70%	60%	50%	40%
R0	0.0027	0.0026	0.0026	0.0026	0.0026	0.0026
R1	6.665e-04	8.362e-04	7.963e-04	7.914e-04	9.227e-04	9.716e-04
C1	1.423e+05	7.949e+04	4.112e+04	6.875e+04	4.943e+04	3.516e+04
R2	0.001013	7.812e-04	7.144e-04	4.208e-04	4.715e-04	2.133e-04
C2	6.752e+05	8.102e+05	1.325e+06	1.082e+07	7.290e+06	3.208e+07
soc	30%	20%	10%
R0	0.0024	0.0025	0.0025
R1	0.0011	0.0012	0.0017
C1	3.014e+04	3.063e+04	2.509e+04
R2	2.821e-04	2.959e-04	4.34e-04
C2	1.071e+07	1.16e+07	6.858e+06