发明名称 一种超级电容充放电过程切换控制方法
摘要 一种超级电容充放电过程切换控制方法,包括如下步骤:步骤一、超级电容充放电过程动态模型设计;步骤二、超级电容的充放电模型变换设计;步骤三、超级电容的充放电切换控制器设计;步骤四,电容充放电的执行。此方法简洁、易用,可用于指导实际的超级电容充放电过程控制。
申请公布号 CN103904734B 申请公布日期 2016.07.06
申请号 CN201410085603.6 申请日期 2014.03.10
申请人 浙江工业大学 发明人 宋秀兰;俞立;孟利民;彭宏
分类号 H02J7/00(2006.01)I 主分类号 H02J7/00(2006.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵;黄美娟
主权项 一种超级电容充放电过程切换控制方法,包括如下步骤:步骤一、根据超级电容充放电机理,建立超级电容充放电过程动态模型,并通过适当假设来简化该模型;超级电容充放电过程动态模型设计以超级电容的剩余荷电状态的变化表示充放电过程,荷电状态等于剩余电量/最大储能电量,应用能量守恒原理,建立超级电容的充放电过程的动态数学模型,如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfrac><mrow><mi>d</mi><mi>Q</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mo>-</mo><msub><mi>&sigma;</mi><mi>c</mi></msub><mi>Q</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&eta;</mi><mi>c</mi></msub><mi>P</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000931790980000011.GIF" wi="1078" he="107" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfrac><mrow><mi>d</mi><mi>Q</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mo>-</mo><msub><mi>&sigma;</mi><mi>d</mi></msub><mi>Q</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><mi>P</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><msub><mi>&eta;</mi><mi>d</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000931790980000012.GIF" wi="1078" he="117" /></maths>其中,式(1)表示充电过程动态模型,式(2)表示放电过程动态模型,t表示充放电的时间(min),Q表示超级电容的荷电状态,P表示充放电功率(MW/min),σ<sub>c</sub>表示充电过程的自放电率(%/min),η<sub>c</sub>表示充电效率(%),σ<sub>d</sub>表示放电过程的自放电率(%/min),1/η<sub>d</sub>表示放电效率(%);步骤二、考虑实际超级电容的充放电的情况,对超级电容充放电模型进行变换设计;令Q<sub>max</sub>为最大贮能电量,Q<sub>ce</sub>为充电过程控制目标贮能电量;P<sub>ce</sub>对应Q<sub>ce</sub>的控制目标充电功率;Q<sub>de</sub>为放电过程控制目标贮能电量;P<sub>de</sub>对应Q<sub>de</sub>的控制目标放电功率;其中,Q<sub>ce</sub>和Q<sub>de</sub>由上一级电网充放电功率配置单元计算给定;定义充电过程状态变量<img file="FDA0000931790980000013.GIF" wi="267" he="127" />输入变量u=P‑P<sub>ce</sub>;放电过程状态变量<img file="FDA0000931790980000014.GIF" wi="270" he="127" />输入变量u=P‑P<sub>de</sub>,整理可得<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfrac><mrow><mi>d</mi><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mo>-</mo><msub><mi>&sigma;</mi><mi>c</mi></msub><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><msub><mi>&eta;</mi><mi>c</mi></msub><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><msub><mi>Q</mi><mi>max</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000931790980000015.GIF" wi="1061" he="118" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfrac><mrow><mi>d</mi><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mo>-</mo><msub><mi>&sigma;</mi><mi>d</mi></msub><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><msub><mi>&eta;</mi><mi>d</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000931790980000016.GIF" wi="1048" he="118" /></maths>步骤三、采用共同控制Lyapunov函数,构造一个解析的切换状态反馈控制律,考虑超级电容充放电变换模型式(3)和(4),令A<sub>1</sub>=‑σ<sub>c</sub>,<img file="FDA0000931790980000017.GIF" wi="214" he="119" />A<sub>2</sub>=‑σ<sub>d</sub>,<img file="FDA0000931790980000018.GIF" wi="261" he="127" />则求解如下的不等式方程组<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>A</mi><mn>1</mn></msub><mi>X</mi><mo>+</mo><msub><mi>XA</mi><mn>1</mn></msub><mo>-</mo><msubsup><mi>XB</mi><mn>1</mn><mn>2</mn></msubsup><mi>X</mi><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>A</mi><mn>2</mn></msub><mi>X</mi><mo>+</mo><msub><mi>XA</mi><mn>2</mn></msub><mo>-</mo><msubsup><mi>XB</mi><mn>2</mn><mn>2</mn></msubsup><mi>X</mi><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000931790980000019.GIF" wi="1214" he="163" /></maths>其中X是未知变量,且X&gt;0;求解式(5)得解X=X<sub>0</sub>&gt;0   (6)利用不等式公式(6)的值,构造如下充电过程输入函数u<sub>c</sub>(t)和放电过程输入函数u<sub>d</sub>(t):<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>u</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>&sigma;</mi><mi>c</mi></msub><msub><mi>Q</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>-</mo><msub><mi>&theta;</mi><mn>1</mn></msub><msqrt><mrow><msup><msub><mi>Q</mi><mi>max</mi></msub><mn>2</mn></msup><msup><msub><mi>&eta;</mi><mi>c</mi></msub><mn>2</mn></msup><mo>+</mo><mfrac><mrow><mn>4</mn><msub><mi>&theta;</mi><mn>2</mn></msub><msup><msub><mi>X</mi><mn>0</mn></msub><mn>2</mn></msup><msup><msub><mi>&eta;</mi><mi>c</mi></msub><mn>4</mn></msup></mrow><mrow><msup><msub><mi>Q</mi><mi>max</mi></msub><mn>2</mn></msup></mrow></mfrac></mrow></msqrt></mrow><msub><mi>&eta;</mi><mi>c</mi></msub></mfrac><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000931790980000021.GIF" wi="1548" he="229" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>u</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msub><mi>&sigma;</mi><mi>d</mi></msub><msub><mi>Q</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><msub><mi>&eta;</mi><mi>d</mi></msub><mo>-</mo><msub><mi>&theta;</mi><mn>1</mn></msub><msqrt><mrow><msup><msub><mi>Q</mi><mi>max</mi></msub><mn>2</mn></msup><msup><msub><mi>&eta;</mi><mi>d</mi></msub><mn>2</mn></msup><msup><msub><mi>&eta;</mi><mi>c</mi></msub><mn>2</mn></msup><mo>+</mo><mfrac><mrow><mn>4</mn><msub><mi>&theta;</mi><mn>2</mn></msub><msup><msub><mi>X</mi><mn>0</mn></msub><mn>2</mn></msup></mrow><mrow><msup><msub><mi>Q</mi><mi>max</mi></msub><mn>2</mn></msup><msup><msub><mi>&eta;</mi><mi>d</mi></msub><mn>4</mn></msup></mrow></mfrac></mrow></msqrt><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000931790980000022.GIF" wi="1525" he="159" /></maths>其中,θ<sub>1</sub>&gt;0和θ<sub>2</sub>&gt;0为可调参数,则定义如下充电过程控制器(9),其中<img file="FDA0000931790980000023.GIF" wi="265" he="127" />P(t)=u<sub>c</sub>(t)+P<sub>ce</sub>      (9)放电过程功率控制器,其中<img file="FDA0000931790980000024.GIF" wi="243" he="130" />P(t)=u<sub>d</sub>(t)+P<sub>de</sub>    (10)步骤四,电容充放电的执行,在超级电容充放电控制计算机上运行实施,可以大致分为3个阶段:4.1参数设置,包括模型参数和充放电过程目标参数;在模型导入界面中,输入超级电容充放电最大的贮能电荷量Q<sub>max</sub>;充电和放电过程自放电率(%/min)σ<sub>c</sub>和σ<sub>d</sub>;超级电容充电效率为η<sub>c</sub>;超级电容放电效率1/η<sub>d</sub>;在控制参数设置界面中,输入充电过程控制目标贮能电量Q<sub>ce</sub>和控制目标充电功率P<sub>ce</sub>;放电过程控制目标贮能电量Q<sub>de</sub>和控制目标放电功率P<sub>de</sub>;输入参数确认后,由控制计算机将设置数据送入计算机存储单元RAM中保存;4.2离线调试;点击组态界面中的“调试”按钮,超级电容充放电控制系统进入控制器调试阶段;分别考虑充电过程和放电过程,以状态变量变化1个单位量为控制目标,调整组态界面中的控制器参数θ<sub>1</sub>和θ<sub>2</sub>,观测状态变量即荷电量和输入变量即充放电功率的控制效果,由此确定一组能良好控制超级电容充放电过程的控制器参数值;参数θ<sub>1</sub>和θ<sub>2</sub>的取值规则:θ<sub>1</sub>和θ<sub>2</sub>为正实数,即θ<sub>1</sub>&gt;0和θ<sub>2</sub>&gt;0;参数θ<sub>1</sub>和θ<sub>2</sub>的取值规则:增大θ<sub>1</sub>和θ<sub>2</sub>的值将缩短充放电过程的过渡时间,但增大充放电过程的电容荷电量变化和功率值,增加对充放电过程扰动的敏感性;相反,减小θ<sub>1</sub>和θ<sub>2</sub>的值将延长充放电过程的过渡时间,但减小充放电过程的电容荷电量变化和功率值,降低对充放电过程扰动的敏感性;因此,实际调试控制器参数θ<sub>1</sub>和θ<sub>2</sub>时,应在充放电过程的过渡时间、电容荷电量变化和功率值容许的范围内综合权衡;4.3在线运行;点击组态界面“运行”按钮,启动超级电容充放电控制计算机的CPU读取超级电容充放电过程模型参数、充放电过程目标参数和最佳控制器参数,并执行“超级电容充放电过程控制程序”,通过在线测量超级电容的实际荷电量,控制充电和放电过程的充入功率和放出功率,实现超级电容充放电过程的有效控制。
地址 310014 浙江省杭州市下城区潮王路18号