发明名称 实现串联锂离子电池组均衡充电的方法
摘要 本发明涉及一种实现串联锂离子电池组均衡充电的方法,方法具体如下:首先,根据锂离子电池充放电的特性,建立锂离子单体电池的物理模型并确定其各个参数;其次,根据单体电池物理模型,在恒流充电过程中采用安时法结合扩展的卡尔曼滤波法计算各时刻单体电池的剩余容量和路端电压;然后,对电池组恒流充电至各单体电池端电压离散时,根据此时端电压求出其剩余容量和剩余容量差;最后,对电池组用小电流进行全程剩余容量调整,直到各单体电池剩余容量均衡。本发明与传统串联锂离子电池组均充方法相比,用小电流对电池组实行全程均充,克服传统方法只在充电末期用大电流均充导致元器件功耗大的缺点,且可靠易行。
申请公布号 CN104065117A 申请公布日期 2014.09.24
申请号 CN201410239931.7 申请日期 2014.05.30
申请人 浙江大学 发明人 宋晨路;姜威;刘涌
分类号 H02J7/00(2006.01)I 主分类号 H02J7/00(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 韩介梅
主权项 实现串联锂离子电池组均衡充电的方法,其特征是步骤如下:1)根据锂离子电池充放电的特性,建立锂离子单体电池的物理模型,该物理模型由理想电压源OCV、电池欧姆内阻R<sub>0</sub>以及R<sub>1</sub>C<sub>1</sub>组成的并联电路、R<sub>2</sub>C<sub>2</sub>组成的并联电路串联而成,设:模型中电池欧姆内阻R<sub>0</sub>、两个电池极化电阻R<sub>1</sub>和R<sub>2</sub>、两个极化电容C<sub>1</sub>和C<sub>2</sub>,两个时间常数τ<sub>1</sub>和τ<sub>2</sub>;根据电池脉冲放电法,通过电池内阻压降法计算电池的欧姆内阻:<img file="FDA0000513932330000011.GIF" wi="229" he="126" />通过锂离子单体电池模型中零输入响应和零状态响应的两个方程如(1)、(2)式所示,用最小二乘拟合法分别计算电池模型中的R<sub>1</sub>C<sub>1</sub>并联回路和R<sub>2</sub>C<sub>2</sub>并联回路的时间常数τ<sub>1</sub>、τ<sub>2</sub>和极化电阻R<sub>1</sub>、R<sub>2</sub>,零输入响应:V<sub>out</sub>=OCV‑U<sub>1</sub>exp(‑t/τ<sub>1</sub>)‑U<sub>2</sub>exp(‑t/τ<sub>2</sub>)     (1)零状态响应:V<sub>out</sub>=U'‑IR<sub>1</sub>(1‑exp(‑t/τ<sub>1</sub>))‑IR<sub>2</sub>(1‑exp(‑t/τ<sub>2</sub>))    (2)其中,U'为脉冲平稳后的路端电压,U<sub>1</sub>、U<sub>2</sub>分别是R<sub>1</sub>C<sub>1</sub>并联回路和R<sub>2</sub>C<sub>2</sub>并联回路的电压,t是脉冲放电的时间,I是单体电池的脉冲电流,V<sub>out</sub>为单体电池的路端电压;根据式(3)状态方程和式(4)量测方程采用安时法结合扩展的卡尔曼滤波法计算充电过程中各个单体电池每个时刻的电池剩余容量S<sub>k</sub>以及相对应的单体电池端电压V<sub>outk</sub><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>S</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>U</mi><mrow><mn>1</mn><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>U</mi><mrow><mn>2</mn><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>&Delta;t</mi><mo>/</mo><msub><mi>&tau;</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>&Delta;t</mi><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>*</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>S</mi><mi>k</mi></msub></mtd></mtr><mtr><mtd><msub><mi>U</mi><mrow><mn>1</mn><mo>,</mo><mi>k</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>U</mi><mrow><mn>2</mn><mo>,</mo><mi>k</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mfrac><mi>&eta;&Delta;t</mi><mi>Q</mi></mfrac></mtd></mtr><mtr><mtd><msub><mi>R</mi><mn>1</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>&Delta;t</mi><mo>/</mo><msub><mi>&tau;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>R</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>&Delta;t</mi><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>*</mo><msub><mi>i</mi><mi>k</mi></msub><mo>+</mo><msub><mi>&omega;</mi><mi>k</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000513932330000012.GIF" wi="1712" he="322" /></maths>V<sub>outk</sub>=OCV<sub>k</sub>+i<sub>k</sub>R<sub>0</sub>+U<sub>1,k</sub>+U<sub>2,k</sub>+ν<sub>k</sub>             (4)式中,S<sub>k</sub>、V<sub>outk</sub>分别为锂离子单体电池采样时刻点k的电池剩余容量和路端电压,η为库仑系数,i<sub>k</sub>、OCV<sub>k</sub>、U<sub>1,k</sub>和U<sub>2,k</sub>分别是采样时刻点k的电流、电池开路电压、R<sub>1</sub>C<sub>1</sub>并联电路电压和R<sub>2</sub>C<sub>2</sub>并联电路电压,Q是电池标称容量,Δt为采样周期,ω<sub>k</sub>为状态方程噪声系数,ν<sub>k</sub>为量测方程噪声系数;将状态方程(3)和量测方程(4)改成x<sub>k+1</sub>=f(x<sub>k</sub>)+ω<sub>k</sub>和y<sub>k</sub>=g(x<sub>k</sub>)+ν<sub>k</sub>,令x<sub>k</sub>=(S<sub>k</sub>,U<sub>1,k</sub>,U<sub>2,k</sub>)<sup>T</sup>                 (5)<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>A</mi><mi>k</mi></msub><mo>=</mo><mfrac><mrow><mo>&PartialD;</mo><mi>f</mi></mrow><mrow><mo>&PartialD;</mo><mi>x</mi></mrow></mfrac><msub><mo>|</mo><mrow><mi>x</mi><mo>=</mo><msub><mi>x</mi><mi>k</mi></msub></mrow></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>&Delta;t</mi><mo>/</mo><msub><mi>&tau;</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>&Delta;t</mi><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000513932330000021.GIF" wi="1306" he="233" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>C</mi><mi>k</mi></msub><mo>=</mo><mfrac><mrow><mo>&PartialD;</mo><mi>g</mi></mrow><mrow><mo>&PartialD;</mo><mi>x</mi></mrow></mfrac><msub><mo>|</mo><mrow><mi>x</mi><mo>=</mo><msubsup><mi>x</mi><mi>k</mi><mo>-</mo></msubsup></mrow></msub><mo>=</mo><mrow><mo>(</mo><mfrac><mrow><mi>dOC</mi><msub><mi>V</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>S</mi><mo>)</mo></mrow></mrow><mi>dS</mi></mfrac><msub><mo>|</mo><mrow><mi>S</mi><mo>=</mo><msubsup><mi>S</mi><mi>k</mi><mo>-</mo></msubsup></mrow></msub><mo>,</mo><mo>-</mo><mn>1</mn><mo>,</mo><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000513932330000022.GIF" wi="1307" he="157" /></maths>其中<img file="FDA0000513932330000023.GIF" wi="60" he="66" />为采样时刻点k之前的x<sub>k</sub>的估计值;根据式(8)至式(12)循环计算:预测计算、预报误差协方差、滤波增益、滤波计算以及预报误差协方差,得到恒流充电时单体电池各个时刻路端电压V<sub>outk</sub>所对应的电池剩余容量S<sub>k</sub>,预测计算:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>x</mi><mi>k</mi><mo>-</mo></msubsup><mo>=</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000513932330000024.GIF" wi="1095" he="74" /></maths>预报误差协方差:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>D</mi><mrow><mo>(</mo><msubsup><mi>x</mi><mi>k</mi><mo>-</mo></msubsup><mo>)</mo></mrow><mo>=</mo><msub><mi>A</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mi>D</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><msubsup><mi>A</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>T</mi></msubsup><mo>+</mo><mi>D</mi><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000513932330000025.GIF" wi="934" he="77" /></maths>滤波增益:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>L</mi><mi>k</mi></msub><mo>=</mo><mi>D</mi><mrow><mo>(</mo><msubsup><mi>x</mi><mi>k</mi><mo>-</mo></msubsup><mo>)</mo></mrow><msubsup><mi>C</mi><mi>k</mi><mi>T</mi></msubsup><msup><mrow><mo>(</mo><msub><mi>C</mi><mi>k</mi></msub><mi>D</mi><mrow><mo>(</mo><msubsup><mi>x</mi><mi>k</mi><mo>-</mo></msubsup><mo>)</mo></mrow><msubsup><mi>C</mi><mi>k</mi><mi>T</mi></msubsup><mo>+</mo><mi>D</mi><mrow><mo>(</mo><mi>&upsi;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000513932330000026.GIF" wi="1087" he="78" /></maths>滤波计算:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>x</mi><mi>k</mi></msub><mo>=</mo><msubsup><mi>x</mi><mi>k</mi><mo>-</mo></msubsup><mo>+</mo><msub><mi>L</mi><mi>k</mi></msub><mo>[</mo><msub><mi>y</mi><mi>k</mi></msub><mo>-</mo><mi>g</mi><mrow><mo>(</mo><msubsup><mi>x</mi><mi>k</mi><mo>-</mo></msubsup><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000513932330000027.GIF" wi="1077" he="71" /></maths>滤波误差协方差:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>D</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mi>E</mi><mo>-</mo><msub><mi>L</mi><mi>k</mi></msub><msub><mi>C</mi><mi>k</mi></msub><mo>)</mo></mrow><mi>D</mi><mrow><mo>(</mo><msubsup><mi>x</mi><mi>k</mi><mo>-</mo></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000513932330000028.GIF" wi="932" he="74" /></maths>其中,E为单位矩阵,D(ω)和D(υ)分别是状态方程噪声系数ω和量测方程噪声系数υ的方差;2)将n个串联锂离子电池组在电池管理系统BMS上进行恒流充电,充电过程中,当在充电电压特征曲线上出现n个单体电池的端电压离散时,根据所测的各单体电池的端电压值V<sub>out,i</sub>求出各个单体电池的剩余容量S<sub>i</sub>,i=1~n,以及所有单体电池的平均剩余容量S=(S<sub>1</sub>+S<sub>2</sub>+…+S<sub>n</sub>)/n;3)计算出各个单体电池的剩余容量S<sub>i</sub>与所有单体电池平均剩余容量S之差ΔQ<sub>i</sub>=S<sub>i</sub>‑S;4)利用电池管理系统BMS双向DC/DC转换电路对ΔQ<sub>i</sub>>0的单体电池进行放电,对ΔQ<sub>i</sub><0的单体电池进行充电,直到各个单体电池的剩余容量S<sub>i</sub>相等,即ΔQ<sub>i</sub>=0。
地址 310027 浙江省杭州市西湖区浙大路38号