基于模型预测控制的锂电池组主动均衡控制方法与流程

本发明涉及锂电池组主动均衡器的控制方法，特别涉及一种基于模型预测控制的锂电池组主动均衡控制方法。

背景技术：

电池组的均衡是实现电池电量有效利用、保障电池组充放电安全的有效技术。由于电池在生产时个体间存在差异以及在使用时工况上的差别，导致随时间增长，各单体电池间的容量差异越来越大，在充电时如果存在某节电池电量较多，当该节电池先充满时，整个电池组就需要停止继续充电(否则已经充满的电池会因过充而爆炸)，如图1(1)所示；同理在放电过程中电量最低的电池首先耗完电荷，导致整个电池组停止继续放电(否则已经放空的电池会因过放导致永久损坏)，如图1(2)所示。上述两种情况都会因为提前终止充电或者放电，限制了电堆电荷的有效利用，因此需要均衡技术保持各电池组电量的一致。目前电动汽车主要使用被动均衡技术，其特点是电路简单可靠，通过与电池并联的开关和放电电阻将电量多的电荷耗散，但该方法造成了电荷的损耗，且均衡电流较小。为提高均衡电流和减少均衡造成的额电荷耗散，主动均衡被提出。

主动均衡根据均衡器拓扑的不同，其工作原理稍有差别，但都是通过电容或者电感作为电荷临时存储的元件，实现高效的电荷转移，本专利使用的拓扑如图2,该拓扑使用双向反激电路，如图3所示，可以实现电流的双向流动，即可以通过电堆对单体电池充电，也可以单体电池对电堆充电，其中mosfet使用凌力尔特公司的ltc3300芯片提供驱动。为解决每个反激电路应该对电池充电还是放电以达到高效均衡的问题，本发明建立了该主动均衡器的模型，使用模型预测控制，将原始优化目标问题转化为最小二乘问题，采用分支界定法以较小的计算代价求得最优的控制序列。

技术实现要素：

本发明主要提供了一种基于模型预测控制的电池组主动均衡控制方法，其提出了一种主动均衡器的模型和均衡的优化目标函数，并通过将原始问题转化为最小二乘的整数规划问题求解最优的控制序列，该方法解决了在使用直接模型预测控制使用暴力枚举法计算量大的问题，对主动均衡器控制策略具有重要意义。

本发明采用如下技术方案：一种基于模型预测控制的锂电池组主动均衡控制方法，所述锂电池组中包含多个单体电池，每个单体电池并联一个由双向反击电路构成的dc/dc转换器，包括以下步骤：

步骤一：建立主动均衡器的状态空间模型：

其中x(k)表示k时刻的系统状态变量，x＝[soc1soc2…socs]^t,soci表示第i节电池的荷电状态，s表示电堆单体电池的数目；k为常量，其中，δt为dc/dc转换器传输电荷的周期，i0为单周期内转移的电荷，c为单个电池的容量；u为均衡器的控制变量矩阵u＝[u1u2…us]^t，其中ui为第i节电池的dc/dc转换器的控制变量，定义为

y(k)为均衡器输出，即输出为所有电池均衡后的荷电状态。

步骤二：确定均衡器的控制目标函数

s.t.ui(k)∈{-1,0,1},i＝1,2,…s(3)

其中n表示模型预测控制的预测步长，为s个组成的s×1向量，表示二范数的平方；λ为权重系数，δu(l|k)＝u(l|k)-u(l-1|k)。

步骤三：将原问题化为最小二乘问题。取u(k)＝[u^t(k)…u^t(k+n-1)]^t，由公式(1)第一个式子可以得到：

x(l+1)＝x(k)+[k*b…k*b…0]*u(k)(5)

0为s×s的零矩阵；

将上式带入公式(1)的第二个式子得到

y(k)＝γx(k)+γu(k)(6)

γ包含n个s×s的单位矩阵i，即y(k)＝[y^t(l+1)…y^t(l+n)]^t，控制目标函数可以化为：

式中y^*为y对应的参考输出。将目标函数进一步写成紧凑形式

式中θ(k)、θ(k)为中间变量，

θ(k)＝((γx(k)-y^*(k))^tγ-λ^t(eu(k-1))^ts)^t

q＝γ^tγ+λs^ts

将公式(8)写成平方形式

j＝(u(k)+q^-1θ(k))^tq(u(k)+q^-1θ(k))+const(k)(9)

const(k)为与u(k)无关的部分。不考虑约束条件(3)(4)，上述目标函数得到的无约束解为

uunc(k)＝-q^-1θ(k)(10)

由定义可知q为对称正定矩阵，可以使用cholesky分解q＝h^th，取可得到

所以最优控制序列uopt(k)可表示为

步骤四：使用分支界定法搜索最优控制序列。引入球半径ρ(k)>0和球心uunc(k)，使得

当u(k)满足在该超球面内且满足约束条件(3)(4)时即是可行解，并在球半径内的可行解中搜索使得控制目标函数j最小的开关序列作为uopt(k)。

步骤五：利用步骤四uopt(k)对各个双向dc/dc转换器实施均衡控制。

本发明的有益效果是：由于本发明使用的主动均衡器的均衡电流为恒定值，因此每个反激电路的状态可以看成离散的三种状态，即对单体电池充电、放电或者静置，因此主动均衡器模型的控制变量为离散整数，在求最优的均衡控制序列过程中涉及整数规划问题，在使用模型预测控制时，只能通过暴力枚举法搜寻最优的控制序列，而本发明通过将原优化命题转化为最小二乘的整数规划为题，使用分支界定法，有效的将待求最优解约束在以无约束解为球心的超球体内，大大减少了搜索过程的计算量。

附图说明

图1为无均衡电池组放电示意图；

图2为主动均衡原理图；

图3为反激电路工作过程(t＝1,s＝12)；

图4为均衡器最优控制信号输出；

图5为soc最大差值5％下的放电均衡过程；

图6为soc最大差值5％下的充电均衡过程；

图7为可行解随着均衡过程进行的变化情况。

具体实施方式

本发明结合具体的试验验证过程对基于模型预测控制的锂电池组主动均衡控制策略作进一步的说明。

一种基于模型预测控制的锂电池组主动均衡控制方法，所述锂电池组中包含多个单体电池，每个单体电池并联一个由双向反击电路构成的dc/dc转换器，包括以下步骤：

步骤一：建立主动均衡器的状态空间模型：

均衡器工作示意图如图2所示，对于单体电池celli，i＝1,2…s

·当对celli放电时，celli失去电荷(s-1)*i0，cellj≠i得到电荷i0；

·当对celli充电时，celli得到电荷(s-1)*i0，cellj≠i失去电荷i0；

·当celli不充电不放电时，celli不会造成电荷改变；

i0表示单次转换的常量电荷数，s表示电堆电池的数目。因此可以将均衡器表示成：

其中x表示系统状态变量，取x＝[soc1soc2…socs]^t，soci表示第i节电池的荷电状态，假设已获得x(k)；表示单次转换导致soc的变化，为常量，δt为dc/dc传输电荷一次的时间，为5.417us，i0为单次传输的能量0.417uc，c为单体电池的容量，为2.9ah；u为均衡器的控制变量矩阵u＝[u1u2…us]^t，其中ui为第i节电池的dc/dc转换器的控制变量，定义为：

y(k)为均衡器输出，即输出为所有电池均衡后的荷电状态。

步骤二：确定均衡器控制目标函数。均衡器主要目的是转移电荷，使各电池间的电量差异达到最小，同时对于每一个双向dc/dc，其电流方向不能突变，即不能直接从放电状态变为充电状态或者直接从充电状态转为放电状态，并且应尽量保持dc/dc状态不变，以减少因为状态转换而导致的能量损耗。考虑直接模型预测控制，目标函数可以写为：

s.t.ui(k)∈{-1,0,1},i＝1,2,…s(3)

其中n表示直接模型预测控制的预测步长，为s个组成的s×1向量，表示二范数的平方；(2)中前半部分衡量均衡的准确性，即使各电池的电量尽量达到一致，后半部分衡量均衡器动作的频繁程度，λ为权重系数，用于调整两者的关系，δu(l|k)＝u(l|k)-u(l-1|k)。因为不允许均衡器从放电直接转换到充电或者反过来，因此有约束条件(4)。

步骤三：将原问题化为最小二乘问题。取u(k)＝[u^t(k)…u^t(k+n-1)]^t，由公示(1)中第一个式子可以得到：

x(l+1)＝x(k)+[k*b…k*b…0]*u(k)(5)

其中l＝k,…,k+n-1，0为s×s的零矩阵。

将上式带入(1)的第二个式子得到：

y(k)＝γx(k)+γu(k)(6)

γ包含n个s×s的单位矩阵i，即y(k)＝[y^t(l+1)…y^t(l+n)]^t，将(6)带入目标函数(2)可以化为：

式中y^*为y对应的参考输出。将上式进一步写成紧凑形式：

式中θ(k)、θ(k)为中间变量，

θ(k)＝((γx(k)-y^*(k))^tγ-λ^t(eu(k-1))^ts)^t

q＝γ^tυ+λs^ts

将(8)写成平方形式：

j＝(u(k)+q^-1θ(k))^tq(u(k)+q^-1θ(k))+const(k)(9)

const(k)为与u(k)无关的部分，即其值随k变化，但与优化u(k)无关。不考虑约束条件(3)(4)，上述代价函数的无约束最小二乘解为：

uunc(k)＝-q^-1θ(k)(10)

由定义可知，q为对称正定矩阵，可以cholesky分解为q＝h^th，h为下三角阵，取可得到：

所以最优控制序列uopt(k)可表示为：

步骤四：使用分支界定法搜索最优控制序列。引入球半径ρ(k)>0和球心uunc(k)，使得：

当u(k)满足在该超球面内且满足约束条件(3)(4)时即是可行解，并在球半径内的可行解中搜索使得代价函数j最小的开关序列作为uopt(k)。因为h是下三角矩阵，(13)可以写成：

其中表示的第i项，ui表示u(k)的第i项，h(i,j)是h的第(i，j)项。因此可以通过顺序计算u(k)的项，只有满足(14)的控制序列是可行解。在每步需要确定初始化的球半径，通过将上一步得到的最优控制序列做一次移位，得到：

则初始化球半径可以表示为：

最后可以通过下面的分支定界算法求解uopt(k):

步骤五：将步骤四获得的uopt(k)作用于各双向dc/dc转换器，实现对电池组的均衡。图4为均衡器最优输出控制序列，满足约束条件(4)；均衡效果如图5所示，在soc初始差值为5％时，通过均衡器对电堆均衡，各个电池的soc趋于一致，耗时约400s。同样的，在充电时，均衡器也能有效的使各电池的soc趋于一致。如果采用暴力枚举，需要枚举的次数为3⁹＝19683，图7为可行解的数目，可以看出最大可行解不超过2000个，因此本专利提出的方法可以节省约90％的计算量。