• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 11
    • 下载费用:30 金币  

    重庆时时彩下载安卓: 电动汽车用锂离子电容器荷电状态估计方法.pdf

    关 键 词:
    电动汽车 锂离子 电容器 状态 估计 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201610808361.8

    申请日:

    2016.09.08

    公开号:

    CN106354964A

    公开日:

    2017.01.25

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G06F 17/50申请日:20160908|||公开
    IPC分类号: G06F17/50 主分类号: G06F17/50
    申请人: 厦门理工学院
    发明人: 卢向军; 周振威; 谢安; 张勇
    地址: 361000 福建省厦门市集美区理工路600号
    优先权:
    专利代理机构: 厦门智慧呈睿知识产权代理事务所(普通合伙) 35222 代理人: 杨唯
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201610808361.8

    授权公告号:

    |||

    法律状态公告日:

    2017.02.22|||2017.01.25

    法律状态类型:

    实质审查的生效|||公开

    摘要

    一种电动汽车用锂离子电容器荷电状态估计方法,包括:S1,建立简化的锂离子电容器的动态二阶模型S2,基于所述动态二阶模型,分别建立状态方程和观测方程;S3,对锂离子电容器参数进行辨识获;S4,将建立简化的锂离子电容器二阶模型带入人工免疫粒子滤波进行SOC估计,获取锂离子电容器荷电状态。

    权利要求书

    1.一种电动汽车用锂离子电容器荷电状态估计方法,包括以下步骤:
    S1,建立简化的锂离子电容器的动态二阶模型:UL=Uocv-I·R0-Up1-Up2 (1),其中,
    I锂离子电容器的负载电流,UL为端电压,Uocv为开路电压,R0为欧姆内阻,Up1、Up2为极化电
    压;
    S2,基于所述动态二阶模型,分别建立状态方程(2)和观测方程(3):
    <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> <mrow> <msub> <msup> <mi>U</mi> <mrow> <mi>p</mi> <mn>1</mn> </mrow> </msup> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <msup> <mi>U</mi> <mrow> <mi>p</mi> <mn>2</mn> </mrow> </msup> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow> </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> <mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&Delta;</mi> <mi>t</mi> <mo>/</mo> <msub> <mi>&tau;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&Delta;</mi> <mi>t</mi> <mo>/</mo> <msub> <mi>&tau;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mfenced open = '[' close = ']'> <mtable> <mtr> <mtd> <msub> <mi>S</mi> <mi>k</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <msup> <mi>U</mi> <mrow> <mi>p</mi> <mn>1</mn> </mrow> </msup> <mi>k</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <msup> <mi>U</mi> <mrow> <mi>p</mi> <mn>2</mn> </mrow> </msup> <mi>k</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfenced open = '[' close = ']'> <mtable> <mtr> <mtd> <mo>-</mo> <mfrac> <mrow> <mi>&eta;</mi> <mi>&Delta;</mi> <mi>t</mi> </mrow> <mi>Q</mi> </mfrac> </mtd> </mtr> <mtr> <mtd> <msub> <mi>R</mi> <mrow> <mi>p</mi> <mn>1</mn> </mrow> </msub> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mrow> <mo>-</mo> <mi>&Delta;</mi> <mi>t</mi> <mo>/</mo> <msub> <mi>&tau;</mi> <mn>1</mn> </msub> </mrow> <mo>)</mo> </mrow> <mo>)</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>R</mi> <mrow> <mi>p</mi> <mn>2</mn> </mrow> </msub> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mrow> <mo>-</mo> <mi>&Delta;</mi> <mi>t</mi> <mo>/</mo> <msub> <mi>&tau;</mi> <mn>2</mn> </msub> </mrow> <mo>)</mo> </mrow> <mo>)</mo> </mtd> </mtr> </mtable> </mfenced> <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>2</mn> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
    <mrow> <msub> <mi>U</mi> <mi>k</mi> </msub> <mo>=</mo> <msub> <mi>U</mi> <mrow> <mi>o</mi> <mi>c</mi> <mi>v</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>S</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>i</mi> <mi>k</mi> </msub> <msub> <mi>R</mi> <mn>0</mn> </msub> <mo>-</mo> <msubsup> <mi>U</mi> <mi>k</mi> <mrow> <mi>p</mi> <mn>1</mn> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>U</mi> <mi>k</mi> <mrow> <mi>p</mi> <mn>2</mn> </mrow> </msubsup> <mo>+</mo> <msub> <mi>v</mi> <mi>k</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
    其中,Q为锂离子电容器实际容量,Δt为采样周期,η为充放电效率,ik为k时刻的电流,
    放电时为负,充电时为正,Up1k表示k时刻Rp1上的极化电压估计值,Up2k表示k时刻Rp2上的极
    化电压估计值,τ1、τ2为RC时间常数,τ1=Rp1Cp1、τ2=Rp2Cp2,Sk表示k时刻的SOC值,ωk为k时刻
    的过程噪声;
    S3,对锂离子电容器参数进行辨识,从而获取Q、η、Uocv、τ1、τ2、Rp1、Rp2、Cp1和Cp2;
    S4,将建立简化的锂离子电容器二阶模型带入人工免疫粒子滤波进行SOC估计,获取锂
    离子电容器荷电状态。
    2.根据权利要求1所述的电动汽车用锂离子电容器荷电状态估计方法,其特征在于:在
    步骤S3中,所述对锂离子电容器参数进行辨识,从而获取Q、η的步骤包括:进行一定温度下、
    电流倍率条件下的锂离子电容器恒流充放电实验,得出锂离子电容器的充放电曲线,对放
    电曲线进行按时间积分得出Q、η。
    3.根据权利要求1所述的电动汽车用锂离子电容器荷电状态估计方法,其特征在于:在
    步骤S3中,所述对锂离子电容器参数进行辨识,从而获取Uocv的步骤包括:在一定温度条件
    开展间隔放电实验以得出SOC-OCV曲线,对放电电流进行按时间积分得到锂离子电容器实
    际电量,然后折算出每次测量锂离子电容器端电压时对应的SOC值,最后得出一定温度条件
    下的SOC-OCV参数。
    4.根据权利要求1所述的电动汽车用锂离子电容器荷电状态估计方法,其特征在于:在
    步骤S3中,所述对锂离子电容器参数进行辨识,从而获取τ1、τ2的步骤包括:使用Matlab软件
    以U01、U02、τ1、τ2作为待测参数,以U01exp(-t/τ1)+U02exp(-t/τ2)作为目标式进行实验数据曲
    线拟合,得到时间常数τ1和τ2。
    5.根据权利要求1所述的电动汽车用锂离子电容器荷电状态估计方法,其特征在于:在
    步骤S3中,所述对锂离子电容器参数进行辨识获,从而获取Rp1、Rp2、Cp1和Cp2的步骤包括:将
    得出的时间常数τ1和τ2和代入中,并使用Matlab
    软件将Rp1和Rp2作为待定参数,进行实验数据曲线拟合,得到极化电阻参数Rp1和Rp2,再由公
    式τ1=Rp1Cp1、τ2=Rp2Cp2计算出极化电容参数Cp1和Cp2。
    6.根据权利要求1所述的电动汽车用锂离子电容器荷电状态估计方法,其特征在于:所
    述将建立简化的锂离子电容器二阶模型带入人工免疫粒子滤波进行SOC估计的步骤包括:
    S41,首先进行初始化,确定锂离子电容器的初始状态,包括极化电压初始值和SOC初始
    值;
    S42,依据初始概率密度函数p(x0)产生粒子集所得粒子的初始权重
    进行下一步时刻k=1的估算;
    S43,在k时刻,利用锂离子电容器状态方程式(2)和k-1时刻的状态值得出粒子的值
    并用锂离子电容器观测方程式(3)计算出对应的
    S44,根据公式(4)计算粒子权值,并根据公式(5)对权值进行归一化处理,以得到每个
    粒子的权值其中,yk为真值,为计算值,
    <mrow> <msubsup> <mi>w</mi> <mi>k</mi> <mi>i</mi> </msubsup> <mo>=</mo> <mfrac> <mn>1</mn> <msqrt> <mrow> <mn>2</mn> <mi>&pi;</mi> <mi>&sigma;</mi> </mrow> </msqrt> </mfrac> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msubsup> <mi>y</mi> <mi>k</mi> <mi>i</mi> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mrow> <mn>2</mn> <msup> <mi>&sigma;</mi> <mn>2</mn> </msup> </mrow> </mfrac> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
    <mrow> <msubsup> <mover> <mi>w</mi> <mo>~</mo> </mover> <mi>k</mi> <mi>i</mi> </msubsup> <mo>=</mo> <msubsup> <mi>w</mi> <mi>k</mi> <mi>i</mi> </msubsup> <mo>/</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msubsup> <mi>w</mi> <mi>k</mi> <mi>i</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
    S45,根据人工免疫算法将求得的粒子集作为初始抗体群,并根据公式
    得出每个抗体的亲和力fitk(i),并对亲和力进行排序,排除排斥力不大
    于0.0001中的1个抗体,选择排在前面的N个抗体来更新记忆单元的数据;
    S46,将记忆单位中的N个抗体作为新一代的粒子,重返上述步骤S44-S45以得到新粒子
    集并对权值进行归一化处理,得到
    S47,若Neff<Nthres则进行随机重采样,产生1个服从[0,1]均匀分布的数ui,如果
    则第m个粒子重采样结果为并将重采样后的粒子权
    值进行归一化;如果Neff>Nthres,则转到下一步骤S48;
    S48,进行状态估计,输出k时刻估计值粒子权值更新为
    S49,返回步骤S43,k加1继续循环,直到k=M时结束。

    说明书

    电动汽车用锂离子电容器荷电状态估计方法

    技术领域

    本发明涉及一种电动汽车用锂离子电容器荷电状态估计方法。

    背景技术

    锂离子电容器是近十年来国内外发展起来的一种介于超级电容器与化学电池二
    者之间的新型储能元件。它具备超级电容器那样的放电功率,也具备化学电池储备电荷的
    能力。与超级电容器相比,锂离子电容器具备更大的电容量、较高的能量;与化学电池相比,
    锂离子电容器具有较高的比功率、极长的使用寿命(充放电循环次数可达十万次以上)且不
    用维护和无环境污染等。因此锂离子电容器是一种高效、实用、环保的能量存储装置,其优
    越的性能得到各方的重视,目前发展十分迅速。随着电极制作工艺的不断进步,锂离子电容
    器的比能量在不断提高且成本在不断下降,因此高比能锂离子电容器不仅具有良好的功率
    特性,而且其能量特性也日益接近铅酸电池、镍氢电池、锂离子电池等传统储能器件。锂离
    子电容器与高能量电池组成的混合电源应用在纯电动汽车时,可作为汽车的辅助电源与动
    力电池配合使用,减少大电流充放电对电池的损害,延长电池的使用寿命。同时锂离子电容
    器的使用可以减少车内用于电制动、电转向等子系统的布线,使车辆的稳定性得到提高。

    荷电状态(SOC)是锂离子电容器应用在混合电源电动汽车储能管理系统中最重要
    的参数之一。随时了解锂离子电容器的SOC值,可为储能系统均衡技术提供基础数据,保证
    混合电源系统整体充电与放电性能,延长混合电源系统使用寿命。但锂离子电容器的SOC值
    是无法实现直接测量。目前,卡尔曼滤波是一种较成熟的预测储能器件SOC的方法,但是由
    锂离子电容器与高能量电池组成的混合电源系统是非常复杂的非线性系统,在处理非线性
    问题时,卡尔曼滤波具有一定的局限性。

    粒子滤波法对系统状态噪声和观测噪声无限制,可适用于任何非线性系统,精度
    可以逼近最优估计,是一种很有效的非线性滤波技术。同时,为了增加粒子多样性,釆用人
    工免疫算法来对粒子滤波进行优化,可对锂离子电容器SOC的估计更加的准确,以满足实际
    需求。

    发明内容

    本发明的目的在于克服现有技术的缺点,提供一种电动汽车用锂离子电容器荷电
    状态估计方法。

    为解决上述技术问题,本发明采用了以下技术措施:

    本发明提供一种电动汽车用锂离子电容器荷电状态估计方法,包括以下步骤:

    S1,建立简化的锂离子电容器的动态二阶模型:UL=Uocv-I·R0-Up1-Up2 (1),其
    中,I锂离子电容器的负载电流,UL为端电压,Uocv为开路电压,R0为欧姆内阻,Up1、Up2为极化
    电压;

    S2,基于所述动态二阶模型,分别建立状态方程和观测方程:



    <mrow> <msub> <mi>U</mi> <mi>k</mi> </msub> <mo>=</mo> <msub> <mi>U</mi> <mrow> <mi>o</mi> <mi>c</mi> <mi>v</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>S</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>i</mi> <mi>k</mi> </msub> <msub> <mi>R</mi> <mn>0</mn> </msub> <mo>-</mo> <msubsup> <mi>U</mi> <mi>k</mi> <mrow> <mi>p</mi> <mn>1</mn> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>U</mi> <mi>k</mi> <mrow> <mi>p</mi> <mn>2</mn> </mrow> </msubsup> <mo>+</mo> <msub> <mi>v</mi> <mi>k</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> <mo>,</mo> </mrow>

    其中,Q为锂离子电容器实际容量,Δt为采样周期,η为充放电效率,ik为k时刻的
    电流,放电时为负,充电时为正,Up1k表示k时刻Rp1上的极化电压估计值,Up2k表示k时刻Rp2上
    的极化电压估计值,τ1、τ2为RC时间常数,τ1=Rp1Cp1、τ2=Rp2Cp2,Sk表示k时刻的SOC值,ωk为k
    时刻的过程噪声;

    S3,对锂离子电容器参数进行辨识获,从而获取Q、η、Uocv、τ1、τ2、Rp1、Rp2、Cp1和Cp2;

    S4,将建立简化的锂离子电容器二阶模型带入人工免疫粒子滤波进行SOC估计,获
    取锂离子电容器荷电状态。

    作为进一步改进的,在步骤S3中,所述对锂离子电容器参数进行辨识,从而获取Q、
    η的步骤包括:进行一定温度下、电流倍率条件下的锂离子电容器恒流充放电实验,得出锂
    离子电容器的充放电曲线,对放电曲线进行按时间积分得出Q、η。

    作为进一步改进的,在步骤S3中,所述对锂离子电容器参数进行辨识,从而获取
    Uocv的步骤包括:在一定温度条件开展间隔放电实验以得出SOC-OCV曲线,对放电电流进行
    按时间积分得到锂离子电容器实际电量,然后折算出每次测量锂离子电容器端电压时对应
    的SOC值,最后得出一定温度条件下的SOC-OCV参数。

    作为进一步改进的,在步骤S3中,所述对锂离子电容器参数进行辨识,从而获取
    τ1、τ2的步骤包括:使用Matlab软件以U01、U02、τ1、τ2作为待测参数,以U01exp(-t/τ1)+U02exp
    (-t/τ2)作为目标式进行实验数据曲线拟合,得到时间常数τ1和τ2。

    作为进一步改进的,在步骤S3中,所述对锂离子电容器参数进行辨识获,从而获取Rp1、Rp2、
    Cp1和Cp2的步骤包括:将得出的时间常数τ1和τ2和代入
    中,并使用Matlab软件将Rp1和Rp2作为待定参数,进行实验数据曲线拟合,得到极化电阻参
    数Rp1和Rp2,再由公式τ1=Rp1Cp1、τ2=Rp2Cp2计算出极化电容参数Cp1和Cp2。

    作为进一步改进的,所述将建立简化的锂离子电容器二阶模型带入人工免疫粒子
    滤波进行SOC估计的步骤包括:

    S41,首先进行初始化,确定锂离子电容器的初始状态,包括极化电压初始值和SOC
    初始值;

    S42,依据初始概率密度函数p(x0)产生粒子集所得粒子的初始权重
    进行下一步时刻k=1的估算;

    S43,在k时刻,利用锂离子电容器状态方程式(2)和k-1时刻的状态值得出粒子的
    值并用锂离子电容器观测方程式(3)计算出对应的

    S44,根据公式(4)计算粒子权值,并根据公式(5)对权值进行归一化处理,以得到
    每个粒子的权值其中,yk为真值,为计算值,

    <mrow> <msubsup> <mi>w</mi> <mi>k</mi> <mi>i</mi> </msubsup> <mo>=</mo> <mfrac> <mn>1</mn> <msqrt> <mrow> <mn>2</mn> <mi>&pi;</mi> <mi>&sigma;</mi> </mrow> </msqrt> </mfrac> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msubsup> <mi>y</mi> <mi>k</mi> <mi>i</mi> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mrow> <mn>2</mn> <msup> <mi>&sigma;</mi> <mn>2</mn> </msup> </mrow> </mfrac> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>

    <mrow> <msubsup> <mover> <mi>w</mi> <mo>~</mo> </mover> <mi>k</mi> <mi>i</mi> </msubsup> <mo>=</mo> <msubsup> <mi>w</mi> <mi>k</mi> <mi>i</mi> </msubsup> <mo>/</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msubsup> <mi>w</mi> <mi>k</mi> <mi>i</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>

    S45,根据人工免疫算法将求得的粒子集作为初始抗体群,并根据公式
    得出每个抗体的亲和力fitk(i),并对亲和力进行排序,排除排斥力不大
    于0.0001中的1个抗体,选择排在前面的N个抗体来更新记忆单元的数据;

    S46,将记忆单位中的N个抗体作为新一代的粒子,重返上述步骤S44-S45以得到新
    粒子集并对权值进行归一化处理,得到

    S47,若Neff<Nthres则进行随机重采样,产生1个服从[0,1]均匀分布的数ui,如果
    则第m个粒子重采样结果为并将重采样后的粒子权
    值进行归一化;如果Neff>Nthres,则转到下一步骤S48;

    S48,进行状态估计,输出k时刻估计值粒子权值更新为

    S49,返回步骤S43,k加1继续循环,直到k=M时结束。

    本发明提供的电动汽车用锂离子电容器荷电状态估计方法具有以下优点:其一,
    基于简化的锂离子电容器二阶模型和参数辨识方法,准确描述锂离子电容器工作特性,用
    于锂离子电容器在储能系统电路中的设计;其二,准确预测混合电池电动汽车用锂离子电
    容器辅助电源的荷电状态,实现锂离子电容器辅助储能系统的状态监测和健康管理。

    附图说明

    图1为简化的锂离子电容器二阶模型。

    图2为模型参数辨识流程。

    图3基于人工免疫粒子滤波的SOC估计流程图。

    具体实施方式

    下面结合附图与具体实施方式对本发明作进一步详细描述。

    改进锂离子电容器二阶模型:

    请参照图1,建立简化的锂离子电容器动态二阶模型,来表征锂离子电容器端电压
    和SOC、内阻等的动态关系。图中I锂离子电容器的负载电流,UL为端电压,Uocv为开路电压;R0
    为欧姆内阻。Rp1、Rp2为极化电阻,Cp1、Cp2为极化电容,用来描述锂离子电容器的极化电压。E
    为理想电压源,电容Cb表示电流累积产生的开路电压,E和Cb一起描述开路电压Uocv的变化。
    根据电路模型,锂离子电容器端电压与电流的关系如式(1)所示:

    UL=Uocv-I·R0-Up1-Up2 (1)

    建立状态方程和观测方程:

    基于简化的锂离子电容器二阶模型,建立离散化后的状态方程如(2)式所示,状态
    量为SOC和两个极化电压组成的矩阵。(2)式中Q为锂离子电容器实际容量;Δt为采样周期;
    η为充放电效率;ik为k时刻的电流,放电时为负,充电时为正;Up1k表示k时刻Rp1上的极化电
    压估计值;Up2k表示k时刻Rp2上的极化电压估计值;τ1、τ2为RC时间常数,τ1=Rp1Cp1、τ2=
    Rp2Cp2;Sk表示k时刻的SOC值;ωk为k时刻的过程噪声。

    锂离子电容器端电压表达式离散化得到观测方程如(3)式所示,观测值为锂离子
    电容器端电压(yk=Uk),vk为k时刻的观测噪声。


    <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> <mrow> <msub> <msup> <mi>U</mi> <mrow> <mi>p</mi> <mn>1</mn> </mrow> </msup> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <msup> <mi>U</mi> <mrow> <mi>p</mi> <mn>2</mn> </mrow> </msup> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow> </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> <mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&Delta;</mi> <mi>t</mi> <mo>/</mo> <msub> <mi>&tau;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&Delta;</mi> <mi>t</mi> <mo>/</mo> <msub> <mi>&tau;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>S</mi> <mi>k</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <msup> <mi>U</mi> <mrow> <mi>p</mi> <mn>1</mn> </mrow> </msup> <mi>k</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <msup> <mi>U</mi> <mrow> <mi>p</mi> <mn>2</mn> </mrow> </msup> <mi>k</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <mfrac> <mrow> <mi>&eta;</mi> <mi>&Delta;</mi> <mi>t</mi> </mrow> <mi>Q</mi> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mrow> <mi>p</mi> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>exp</mi> <mo>(</mo> <mrow> <mo>-</mo> <mi>&Delta;</mi> <mi>t</mi> <mo>/</mo> <msub> <mi>&tau;</mi> <mn>1</mn> </msub> </mrow> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mrow> <mi>p</mi> <mn>2</mn> </mrow> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>exp</mi> <mo>(</mo> <mrow> <mo>-</mo> <mi>&Delta;</mi> <mi>t</mi> <mo>/</mo> <msub> <mi>&tau;</mi> <mn>2</mn> </msub> </mrow> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <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>2</mn> <mo>)</mo> </mrow> </mrow>

    <mrow> <msub> <mi>U</mi> <mi>k</mi> </msub> <mo>=</mo> <msub> <mi>U</mi> <mrow> <mi>o</mi> <mi>c</mi> <mi>v</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>S</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>i</mi> <mi>k</mi> </msub> <msub> <mi>R</mi> <mn>0</mn> </msub> <mo>-</mo> <msubsup> <mi>U</mi> <mi>k</mi> <mrow> <mi>p</mi> <mn>1</mn> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>U</mi> <mi>k</mi> <mrow> <mi>p</mi> <mn>2</mn> </mrow> </msubsup> <mo>+</mo> <msub> <mi>v</mi> <mi>k</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>

    锂离子电容器参数辨识方法

    1.容量、充放电效率辨识:

    进行一定温度下、电流倍率条件下的锂离子电容器恒流充放电实验,得出锂离子
    电容器的充放电曲线,对放电曲线进行按时积分得出的锂离子电容器容量、充放电效率。

    2.Uocv辨识:

    在一定温度条件开展间隔放电实验以得出SOC-OCV曲线。先将锂离子电容器充满
    电,静置6小时后测出锂离子电容器端电压得出SOC=1时的Uocv值。在一定放电电流条件下,
    放电数秒后再静置6小时测锂离子电容器端电压,重复该过程(放电数秒,静置6小时后测锂
    离子电容器端电压)直到SOC=0。对放电电流进行按时积分得到锂离子电容器实际电量,然
    后折算出每次测量锂离子电容器端电压时对应的SOC值,最后得出一定温度条件下的SOC-
    OCV参数。

    3.电路参数辨识:

    模型中电路参数辨识参考《功率辅助型混合动力汽车用动力电池测试手册》中的
    HPPC(混合动力脉冲能力特性)方法。首先将锂离子电容器充满电,之后在SOC为1处以1C电
    流进行脉冲放电和充电,然后静置15分钟,再以1C恒流放电到SOC=0.9处,静置30分钟后,
    再进行HPPC充放电实验,重复该过程直到SOC=0。对上述HPPC实验数据进行分析,具体模型
    参数辨识流程如图2所示。

    请参照图2,图2中ΔU是负载加载瞬间电压变化值;I锂离子电容器的负载电流;Δ
    Uocv为HPPC充、放电脉冲引起的锂离子电容器端电压的改变;ΔQ为放电电量;U01和U02分别
    为初始极化电压;U′为放电开始时刻电压垂直下降后的起点值。锂离子电容器放电结束后
    端电压会缓慢上升,这是由于极化电容对极化电阻放电的过程,是RC并联回路的零输入响
    应,根据这个特点,使用Matlab软件以U01、U02、τ1、τ2作为待测参数,以U01exp(-t/τ1)+U02exp
    (-t/τ2)作为目标式进行实验数据曲线拟合,得到时间常数τ1和τ2。

    此外锂离子电容器放电时端电压缓慢下降,这是放电电流对极化电容充电的过
    程,是RC并联回路的零状态响应。根据这个特点,将得出的时间常数和代入
    式中,使用Matlab软件将Rp1和Rp2作为待定参数,
    进行实验数据曲线拟合,得到极化电阻参数Rp1和Rp2。由τ1=Rp1Cp1、τ2=Rp2Cp2计算出极化电
    容参数Cp1和Cp2。

    基于人工免疫粒子滤波的SOC估计:

    请参照图3,将建立简化的锂离子电容器二阶模型带入人工免疫粒子滤波进行SOC
    估计,具体流程如图3所示。

    1、首先进行初始化,例如初始迭代变量k=0,迭代次数为M,粒子数为N,噪声值,有
    效粒子阈值等。确定锂离子电容器的初始状态,包括极化电压初始值和SOC初始值。

    2、依据初始概率密度函数p(x0)产生粒子集所得粒子的初始权重
    进行下一步时刻k=1的估算。

    3、在k时刻,利用锂离子电容器状态方程式(2)和k-1时刻的状态值得出粒子的值
    并用锂离子电容器观测方程式(3)计算出对应的

    4、粒子权值用公式(4)计算,以得到每个粒子的权值式(4)中yk为真值,为计
    算值。用(5)式对权值进行归一化处理。

    <mrow> <msubsup> <mi>w</mi> <mi>k</mi> <mi>i</mi> </msubsup> <mo>=</mo> <mfrac> <mn>1</mn> <msqrt> <mrow> <mn>2</mn> <mi>&pi;</mi> <mi>&sigma;</mi> </mrow> </msqrt> </mfrac> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>k</mi> </msub> <mo>-</mo> <msubsup> <mi>y</mi> <mi>k</mi> <mi>i</mi> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mrow> <mn>2</mn> <msup> <mi>&sigma;</mi> <mn>2</mn> </msup> </mrow> </mfrac> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>

    <mrow> <msubsup> <mover> <mi>w</mi> <mo>~</mo> </mover> <mi>k</mi> <mi>i</mi> </msubsup> <mo>=</mo> <msubsup> <mi>w</mi> <mi>k</mi> <mi>i</mi> </msubsup> <mo>/</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msubsup> <mi>w</mi> <mi>k</mi> <mi>i</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>

    5、人工免疫算法

    (a)根据人工免疫算法特点,将求得的粒子集作为初始抗体群,并根据
    公式(6)得出每个抗体的亲和力fitk(i)。fitk(i)表示k时刻第i个抗体的亲和力,值越大表
    明抗原和抗体匹配的越差。

    <mrow> <msub> <mi>fit</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>1</mn> <mo>-</mo> <msubsup> <mover> <mi>&omega;</mi> <mo>~</mo> </mover> <mi>k</mi> <mi>i</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>

    (b)抗体克隆。根据人工免疫算法特点,根据每一个抗体的亲和力对其进行克隆,
    克隆数目按式(7)确定。式(7)中,CLnumk(i)表示k时刻第i个抗体的克隆数,round()表示向
    最近的整数取整??寺∧康脑谟谝种魄缀土Υ蟮目固?,增加亲和力小的抗体,使得整个体系
    能快速的收敛于全局最优解。

    <mrow> <msub> <mi>CLnum</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>r</mi> <mi>o</mi> <mi>u</mi> <mi>n</mi> <mi>d</mi> <mrow> <mo>(</mo> <mi>N</mi> <mo>&times;</mo> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mo>(</mo> <mfrac> <mrow> <mi>&pi;</mi> <mo>&CenterDot;</mo> <msub> <mi>fit</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> </mrow> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <mo>)</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>

    (c)变异。根据人工免疫算法特点,按式(8)对所有抗体进行变异。式(8)中randn表
    示随机抽取的满足N(0,l)分布的随机数。变异遵循原则为抗体亲和力小的变异量小,抗体
    亲和力大的变异量大。

    <mrow> <msubsup> <mi>x</mi> <mi>k</mi> <mi>i</mi> </msubsup> <mo>=</mo> <msubsup> <mi>x</mi> <mi>k</mi> <mi>i</mi> </msubsup> <mo>+</mo> <msub> <mi>fit</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>&times;</mo> <mi>r</mi> <mi>a</mi> <mi>n</mi> <mi>d</mi> <mi>n</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>

    (d)选优。根据人工免疫算法特点,计算变异后所有抗体排斥力offk(i,j)和亲和
    力fitk(i)。对亲和力进行排序,并排除排斥力不大于0.0001中的1个抗体,选在排出前面的
    N个抗体来更新记忆单元的数据。

    6、粒子权值重构。将记忆单位中的N个抗体当做新一代的粒子,重返上述步骤3、步
    骤4以得到新粒子集并对权值进行归一化处理,得到

    7、若Neff<Nthres则进行随机重采样,产生1个服从[0,1]均匀分布的数ui。如果
    则第m个粒子重采样结果为并将重采样后的粒子权
    值进行归一化。如果Neff>Nthres,则转到下一步骤8。

    8、进行状态估计,输出k时刻估计值粒子权值更新为

    9、返回步骤3,k加1继续循环,只到k=M时结束。

    以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精
    神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明?;さ姆段е?。

    关于本文
    本文标题:电动汽车用锂离子电容器荷电状态估计方法.pdf
    链接地址://www.4mum.com.cn/p-6027188.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 体育彩票31选7开奖 北京赛车587 新11选5开奖直播 腾讯二八杠 茗彩彩票游戏 江苏快三遗漏数据统计 865棋牌 棋牌游戏大厅 福彩2017131期号码预测 体彩20选5最新开奖结果查询结果 甘肃十一选五任五遗漏 14场胜负彩澳客网 腾讯欢乐捕鱼cd-key兑换码怎么兑换 官方华东15选5走势图 3d组选365前后 幸运28属于合法彩票吗 七乐彩开奖号码结果查询