发明名称 一种估计电动汽车的动力电池组的荷电状态的方法
摘要 本发明涉及电动汽车的车载电池管理系统的设计领域,尤其涉及一种对电动汽车的动力电池组的荷电状态进行估计的方法。为提高动力电池的荷电状态的估计精度及估计结果的可靠性,本发明提出一种估计电动汽车的动力电池组的荷电状态的方法,采集动力电池组的端电压和充放电电流,建立动力电池组的戴维宁模型、双极化模型和3阶RC网络模型,分别配合状态观测器对动力电池组的荷电状态进行估计得到<img file="DDA0000836872660000011.GIF" wi="235" he="79" />和<img file="DDA0000836872660000012.GIF" wi="116" he="78" />对<img file="DDA0000836872660000013.GIF" wi="241" he="78" />和<img file="DDA0000836872660000014.GIF" wi="86" he="78" />进行加权计算得到动力电池组的荷电状态的估计值z<sub>融,k</sub>,<img file="DDA0000836872660000015.GIF" wi="806" he="79" />k时刻的加权系数w<sub>1</sub>(k)+w<sub>2</sub>(k)+w<sub>3</sub>(k)=1。采用该估计方法估计电动汽车的动力电池组的荷电状态,估计精度较高,估计结果稳定、可靠。
申请公布号 CN105425154A 申请公布日期 2016.03.23
申请号 CN201510732371.3 申请日期 2015.11.02
申请人 北京理工大学 发明人 熊瑞;穆浩;何洪文;孙逢春
分类号 G01R31/36(2006.01)I 主分类号 G01R31/36(2006.01)I
代理机构 北京市邦道律师事务所 11437 代理人 薛艳;王丽影
主权项 一种估计电动汽车的动力电池组的荷电状态的方法,其特征在于,该方法包括如下步骤:步骤1、在所述动力电池组充放电的过程中,对动力电池的端电压U<sub>t</sub>和充放电电流I<sub>L</sub>进行采样,且采样时间间隔为Δt;步骤2、分别建立所述动力电池组的戴维宁模型、双极化模型和3阶网络RC模型,并分别辨识出所述动力电池组的戴维宁模型、双极化模型及3阶RC网络模型的模型参数;步骤3、估计所述动力电池组的荷电状态:首先,根据所述动力电池组的等效电路模型建立所述动力电池组的状态方程和测量方程:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mi>k</mi></msub><mo>=</mo><msub><mi>Ax</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>Bu</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>G&Phi;</mi><mi>k</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>y</mi><mi>k</mi></msub><mo>=</mo><msub><mi>Cx</mi><mi>k</mi></msub><mo>+</mo><msub><mi>Du</mi><mi>k</mi></msub><mo>+</mo><msub><mi>H&Phi;</mi><mi>k</mi></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000836872630000011.GIF" wi="670" he="189" /></maths>其中,x<sub>k</sub>为所述动力电池组在k时刻的系统估计状态,y<sub>k</sub>为所述动力电池组在k时刻的观测矩阵,u<sub>k</sub>为所述动力电池组在k时刻的系统输入矩阵,Φ<sub>k</sub>为所述动力电池组的噪声向量,当所述动力电池组的等效电路模型为戴维宁模型时,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mo>-</mo><mn>1</mn><mo>/</mo><msub><mi>C</mi><mrow><mi>D</mi><mn>1</mn></mrow></msub><msub><mi>R</mi><mrow><mi>D</mi><mn>1</mn></mrow></msub></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000836872630000012.GIF" wi="550" he="206" /></maths>B=[1/C<sub>D1</sub> 1/C<sub>Nom</sub>]<sup>T</sup>,C=[1 dU<sub>oc</sub>/dz],D=R<sub>0</sub>,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>G</mi><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><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000836872630000021.GIF" wi="390" he="181" /></maths>H=[0 0 1],x<sub>k</sub>=[U<sub>D1,k</sub>,z<sub>融,k</sub>]T,U<sub>D1,k</sub>为所述动力电池组在k时刻的极化电压,z<sub>融,k</sub>为所述动力电池组在k时刻的荷电状态的估计值;当所述动力电池组的等效电路模型为双极化模型时,<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mo>-</mo><mn>1</mn><mo>/</mo><msub><mi>C</mi><mrow><mi>D</mi><mn>1</mn></mrow></msub><msub><mi>R</mi><mrow><mi>D</mi><mn>1</mn></mrow></msub></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mn>1</mn><mo>/</mo><msub><mi>C</mi><mrow><mi>D</mi><mn>2</mn></mrow></msub><msub><mi>R</mi><mrow><mi>D</mi><mn>2</mn></mrow></msub></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000836872630000022.GIF" wi="1026" he="390" /></maths>B=[1/C<sub>D1</sub> 1/C<sub>D2</sub> 1/C<sub>Nom</sub>]<sup>T</sup>,C=[1 2 dU<sub>oc</sub>/dz],D=R<sub>0</sub>,<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>G</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000836872630000023.GIF" wi="549" he="302" /></maths>H=[0 0 0 1],x<sub>k</sub>=[U<sub>D1,k</sub>,U<sub>D2,k</sub> z<sub>融,k</sub>]<sup>T</sup>,U<sub>D1,k</sub>和U<sub>D2,k</sub>为所述动力电池组在k时刻的极化电压,z<sub>融,k</sub>为所述动力电池组在k时刻的荷电状态的估计值;当所述动力电池组的等效电路模型为3阶RC网络模型时,<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mo>-</mo><mn>1</mn><mo>/</mo><msub><mi>C</mi><mrow><mi>D</mi><mn>1</mn></mrow></msub><msub><mi>R</mi><mrow><mi>D</mi><mn>1</mn></mrow></msub></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mn>1</mn><mo>/</mo><msub><mi>C</mi><mrow><mi>D</mi><mn>2</mn></mrow></msub><msub><mi>R</mi><mrow><mi>D</mi><mn>2</mn></mrow></msub></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mn>1</mn><mo>/</mo><msub><mi>C</mi><mrow><mi>D</mi><mn>3</mn></mrow></msub><msub><mi>R</mi><mrow><mi>D</mi><mn>3</mn></mrow></msub></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000836872630000031.GIF" wi="1422" he="487" /></maths>B=[1/C<sub>D1</sub> 1/C<sub>D2</sub> 1/C<sub>D3</sub> 1/C<sub>Nom</sub>]<sup>T</sup>,C=[1 1 1 dU<sub>oc</sub>/dz],D=R<sub>0</sub>,<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>G</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000836872630000032.GIF" wi="653" he="399" /></maths>H=[0 0 0 0 1],x<sub>k</sub>=[U<sub>D1,k</sub>,U<sub>D2,k</sub>,U<sub>D3,k</sub>,z<sub>融,k</sub>]<sup>T</sup>,U<sub>D1,k</sub>、U<sub>D2,k</sub>和U<sub>D3,k</sub>为所述动力电池组在k时刻的极化电压,z<sub>融,k</sub>为所述动力电池组在k时刻的荷电状态的估计值;其中,C<sub>D</sub>为所述动力电池组的极化电容,R<sub>D</sub>为所述动力电池组的极化电阻,R<sub>0</sub>为所述动力电池组的直流内阻;C<sub>Nom</sub>为所述动力电池组的额定容量,dU<sub>oc</sub>/dz为所述动力电池组的开路电压U<sub>oc</sub>对其荷电状态的一阶导数;然后,采用三个状态观测器分别与所述动力电池组的戴维宁模型、双极化模型以及3阶RC网络模型配合对所述动力电池组的荷电状态和端电压进行估计,并对所述状态观测器进行时间更新和测量更新得到所述动力电池组的荷电状态在k时刻的估计值<img file="FDA0000836872630000041.GIF" wi="230" he="78" />和<img file="FDA0000836872630000042.GIF" wi="117" he="78" />所述动力电池组在k时刻的端电压的估计值<img file="FDA0000836872630000043.GIF" wi="291" he="87" />和<img file="FDA0000836872630000044.GIF" wi="143" he="86" />最后,对<img file="FDA0000836872630000045.GIF" wi="234" he="78" />和<img file="FDA0000836872630000046.GIF" wi="100" he="78" />进行加权计算得到所述动力电池组的荷电状态在k时刻的估计值z<sub>融,k</sub>,且<img file="FDA0000836872630000047.GIF" wi="990" he="87" />其中,w<sub>1</sub>(k)、w<sub>2</sub>(k)和w<sub>3</sub>(k)为在k时刻的加权系数,分别根据所述动力电池组在k时刻的端电压的估计值<img file="FDA0000836872630000048.GIF" wi="326" he="86" />和<img file="FDA0000836872630000049.GIF" wi="118" he="86" />计算得出,且w<sub>1</sub>(k)+w<sub>2</sub>(k)+w<sub>3</sub>(k)=1。
地址 100081 北京市海淀区中关村南大街5号