发明名称 一种用于飞机电刹车无刷直流电机的灰色滑模控制方法
摘要 本发明提供了一种用于飞机电刹车无刷直流电机的灰色滑模控制方法,采集电机转速和两组电机绕组的电流反馈信号,经压力观测器得电机转速微分和压力观测值微分,计算观测的刹车压力后建立灰色模型,计算刹车压力误差并建立滑模线,计算角速度误差信号后进行速度环PID控制,最终计算电流误差信号,进行均流环PID控制。本发明提高了刹车系统的可靠性和鲁棒性,减少了外界环境对压力传感器的干扰。
申请公布号 CN103223938B 申请公布日期 2015.07.15
申请号 CN201310163415.6 申请日期 2013.05.06
申请人 西北工业大学 发明人 林辉;戴志勇;齐蓉;孙欣欣;郭皓
分类号 B60T13/74(2006.01)I;B60T8/17(2006.01)I 主分类号 B60T13/74(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 顾潮琪
主权项 一种用于飞机电刹车无刷直流电机的灰色滑模控制方法,其特征在于包括下述步骤:第一步,控制器分别采集电机转速ω<sub>r</sub>、电机绕组1电流反馈信号i<sub>1</sub>及电机绕组2电流反馈信号i<sub>2</sub>;第二步,经压力观测器:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mover><mi>&omega;</mi><mover><mo>^</mo><mo>&CenterDot;</mo></mover></mover><mi>r</mi></msub><mo>=</mo><mfrac><msub><mi>C</mi><mrow><mi>T</mi><mn>1</mn></mrow></msub><mi>J</mi></mfrac><msub><mi>i</mi><mn>1</mn></msub><mo>+</mo><mfrac><msub><mi>C</mi><mrow><mi>T</mi><mn>2</mn></mrow></msub><mi>J</mi></mfrac><msub><mi>i</mi><mn>2</mn></msub><mo>-</mo><mfrac><mn>1</mn><mi>J</mi></mfrac><mover><mi>P</mi><mo>^</mo></mover><mo>-</mo><mfrac><mi>B</mi><mi>J</mi></mfrac><msub><mover><mi>&omega;</mi><mo>^</mo></mover><mi>r</mi></msub><mo>+</mo><msub><mi>c</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>r</mi></msub><mo>-</mo><msub><mover><mi>&omega;</mi><mo>^</mo></mover><mi>r</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mover><mi>P</mi><mover><mo>^</mo><mo>&CenterDot;</mo></mover></mover><mo>=</mo><msub><mi>c</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>r</mi></msub><mo>-</mo><msub><mover><mi>&omega;</mi><mo>^</mo></mover><mi>r</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000603367920000011.GIF" wi="871" he="214" /></maths>得电机转速微分<img file="FDA0000603367920000012.GIF" wi="64" he="82" />和压力观测值微分<img file="FDA0000603367920000013.GIF" wi="75" he="77" />其中,C<sub>T1</sub>为电机绕组1转矩常数,C<sub>T2</sub>为电机绕组2转矩常数,J为负载转动惯量,<img file="FDA0000603367920000014.GIF" wi="46" he="65" />为压力观测值,<img file="FDA0000603367920000015.GIF" wi="62" he="68" />为电机转速观测值,c<sub>1</sub>为观测器系数,0<c<sub>1</sub><100000,c<sub>2</sub>为观测器系数,‑50000<c<sub>2</sub><0;第三步,计算观测的刹车压力<img file="FDA0000603367920000016.GIF" wi="202" he="135" />第四步,建立灰色模型GM(1,1):x<sup>(0)</sup>(k)+az<sup>(1)</sup>(k)=b,系统相应的白化方程为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfrac><mrow><msup><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mi>t</mi></mrow><mi>dt</mi></mfrac><mo>+</mo><msup><mi>ax</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>b</mi><mo>,</mo></mrow>]]></math><img file="FDA0000603367920000017.GIF" wi="430" he="145" /></maths>其中,‑a为发展系数,b为灰色作用量,连续采集t‑4,...t,时刻计算的观测的刹车压力<img file="FDA0000603367920000018.GIF" wi="353" he="83" />构成非负序列<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><mi>X</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup><mo>=</mo><mrow><mo>(</mo><mover><mi>P</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mn>4</mn><mo>)</mo></mrow><mo>,</mo><mover><mi>P</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mn>3</mn><mo>)</mo></mrow><mo>,</mo><mover><mi>P</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo><mover><mi>P</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mover><mi>P</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000603367920000019.GIF" wi="974" he="86" /></maths>X<sup>(0)</sup>的1‑AGO序列为:X<sup>(1)</sup>=(x<sup>(1)</sup>(1),x<sup>(1)</sup>(2),x<sup>(1)</sup>(3),x<sup>(1)</sup>(4),x<sup>(1)</sup>(5));其中,<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msup><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>t</mi><mo>-</mo><mn>4</mn></mrow><mrow><mi>t</mi><mo>-</mo><mn>5</mn><mo>+</mo><mi>k</mi></mrow></munderover><msup><mi>&theta;</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>,</mo><mi>k</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mo>;</mo></mrow>]]></math><img file="FDA00006033679200000113.GIF" wi="722" he="143" /></maths>第五步,计算X<sup>(1)</sup>的近邻均值生成序列:Z<sup>(1)</sup>=(z<sup>(1)</sup>(2),z<sup>(1)</sup>(3),z<sup>(1)</sup>(4),z<sup>(1)</sup>(5)),其中,<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msup><mi>z</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><msup><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msup><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo><mi>k</mi><mo>=</mo><mn>2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mo>;</mo></mrow>]]></math><img file="FDA00006033679200000114.GIF" wi="925" he="137" /></maths>第六步,计算中间矩阵Y,B<sup>T</sup>Y=[θ(t‑3),θ(t‑2),θ(t‑1),θ(t)]<sup>T</sup><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msup><mi>B</mi><mi>T</mi></msup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mrow><mo>-</mo><mi>z</mi></mrow><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mtd><mtd><msup><mrow><mo>-</mo><mi>z</mi></mrow><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mn>3</mn></mrow></mtd><mtd><msup><mrow><mo>-</mo><mi>z</mi></mrow><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mtd><mtd><msup><mrow><mo>-</mo><mi>z</mi></mrow><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA00006033679200000112.GIF" wi="901" he="162" /></maths>则灰色系统GM(1,1)的参数为:[a,b]<sup>T</sup>=(B<sup>T</sup>B)<sup>‑1</sup>B<sup>T</sup>Y;第七步,计算刹车系统t+1时刻预测压力:<img file="FDA0000603367920000021.GIF" wi="607" he="123" />第八步,控制器接收给定刹车压力P<sub>gv</sub>;第九步,计算刹车压力误差<img file="FDA0000603367920000022.GIF" wi="266" he="87" />第十步,建立滑模线<img file="FDA0000603367920000023.GIF" wi="217" he="62" />其中,<img file="FDA0000603367920000024.GIF" wi="40" he="62" />为对压力误差的一次微分;c为误差系数,0<c<10;第十一步,采集电机角速度ω,计算角速度给定ω<sub>gv</sub>,其中:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mi>gv</mi></msub><mo>=</mo><mo>[</mo><mi>Jc</mi><msub><mover><mi>P</mi><mo>&CenterDot;</mo></mover><mi>gv</mi></msub><mo>-</mo><mi>Jc&omega;</mi><mo>+</mo><mi>J</mi><msub><mover><mi>P</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mi>gv</mi></msub><mo>+</mo><mi>B&omega;</mi><mo>+</mo><msub><mi>T</mi><mi>L</mi></msub><mo>+</mo><mi>J&xi;sgn</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>+</mo><mi>Jks</mi><mo>]</mo><mo>/</mo><msub><mi>K</mi><mi>T</mi></msub></mrow>]]></math><img file="FDA0000603367920000025.GIF" wi="1158" he="92" /></maths>式中:J为负载转动惯量,<img file="FDA0000603367920000026.GIF" wi="73" he="87" />为对给定压力的一次微分,<img file="FDA0000603367920000027.GIF" wi="68" he="90" />为给定压力的二次微分,B为粘滞系数,T<sub>L</sub>为负载转矩,K<sub>T</sub>为绕组平均转矩系数,ξ>0,k>0;第十二步,计算角速度误差信号e<sub>ω</sub>=ω<sub>gv</sub>‑ω;第十三步,进行速度环PID控制,PID输出作为电流环给定<img file="FDA0000603367920000028.GIF" wi="63" he="69" />第十四步,计算均流电流给定<img file="FDA0000603367920000029.GIF" wi="118" he="85" /><maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mover><mi>i</mi><mo>^</mo></mover><mn>1</mn></msub><mo>=</mo><mover><mi>i</mi><mo>^</mo></mover><mo>/</mo><mn>2</mn></mtd></mtr><mtr><mtd><msub><mover><mi>i</mi><mo>^</mo></mover><mn>2</mn></msub><mo>=</mo><mover><mi>i</mi><mo>^</mo></mover><mo>/</mo><mn>2</mn></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA00006033679200000210.GIF" wi="230" he="181" /></maths>第十五步,计算电流误差信号<img file="FDA00006033679200000211.GIF" wi="207" he="82" />及<img file="FDA00006033679200000212.GIF" wi="239" he="98" />第十六步,进行均流环PID控制,驱动双余度无刷直流电机工作;转至第一步。
地址 710072 陕西省铜川市友谊西路127号