发明名称 一种六旋翼飞行器有一个电机停机时的安全着陆方法
摘要 一种六旋翼飞行器有一个电机停机时的安全着陆方法,该方法有三大步骤:步骤一:六旋翼飞行器姿态和高度控制器设计;步骤二:六旋翼飞行器所有电机正常工作时控制分配器设计;步骤三:当六旋翼飞行器有一个电机停机时的控制重分配。六旋翼飞行器有一个电机停机后其全状态是不可控的,但是当放弃偏航通道后,其俯仰和滚转通道依然可控。本发明即在此基础上采用控制重分配的方式设计六旋翼飞行器有一个电机停机的安全着陆控制器。通过Matlab仿真和真实的实验可以验证本发明是可行的。它在飞行控制技术领域里有较好的实用价值和应用前景。
申请公布号 CN102902277B 申请公布日期 2015.09.30
申请号 CN201210398628.2 申请日期 2012.10.19
申请人 北京航空航天大学 发明人 全权;杜光勋;杨斌先;蔡开元
分类号 G05D1/10(2006.01)I 主分类号 G05D1/10(2006.01)I
代理机构 北京慧泉知识产权代理有限公司 11232 代理人 王顺荣;唐爱华
主权项 一种六旋翼飞行器有一个电机停机时的安全着陆方法,其特征在于:该方法具体步骤如下:步骤一:六旋翼飞行器姿态和高度控制器设计六旋翼飞行器姿态控制器采用如下控制器:τ<sub>a</sub>=‑K<sub>Ω</sub>Ω‑K<sub>q</sub>q   (5)其中,K<sub>Ω</sub>∈R<sup>3×3</sup>为对称正定阵,Ω为六旋翼飞行器绕机体轴的角速度,K<sub>q</sub>∈R是正系数,四元数(q<sub>0</sub>,q<sup>T</sup>)<sup>T</sup>,其中q<sub>0</sub>∈R,q∈R<sup>3</sup>表示六旋翼的方向,用欧拉角表示的姿态角x=(ψ,θ,φ)<sup>T</sup>与四元数的对应关系如下:q<sub>0</sub>=cos(φ/2)cos(θ/2)cos(ψ/2)+sin(φ/2)sin(θ/2)sin(ψ/2)q<sub>1</sub>=sin(φ/2)cos(θ/2)cos(ψ/2)‑cos(φ/2)sin(θ/2)sin(ψ/2)                                                                  (6)q<sub>2</sub>=cos(φ/2)sin(θ/2)cos(ψ/2)+sin(φ/2)cos(θ/2)sin(ψ/2)q<sub>3</sub>=cos(φ/2)cos(θ/2)sin(ψ/2)+sin(φ/2)sin(θ/2)cos(ψ/2)六旋翼飞行器的高度控制采用的控制器为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>F</mi><mo>=</mo><mfrac><mi>m</mi><mrow><mi>cos</mi><mi></mi><mi>&phi;</mi><mi>cos</mi><mi>&theta;</mi></mrow></mfrac><mrow><mo>(</mo><mo>-</mo><msub><mi>K</mi><mi>P</mi></msub><mrow><mo>(</mo><mi>h</mi><mo>-</mo><msub><mi>h</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>K</mi><mi>D</mi></msub><msub><mi>V</mi><mi>h</mi></msub><mo>-</mo><mi>g</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000722893470000011.GIF" wi="1280" he="130" /></maths>其中,m为飞行器质量,K<sub>P</sub>,K<sub>D</sub>为正系数,h为飞行器高度,h<sub>0</sub>为期望的高度,V<sub>h</sub>为垂直速度,g为重力加速度;步骤二:六旋翼飞行器所有电机正常工作时控制分配器设计虚拟控制量u和旋翼升力f之间的映射关系由下列等式表示:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>u</mi><mo>=</mo><mi>Wf</mi><mo>,</mo><mi>W</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>b</mi><mi>f</mi><mi>T</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>db</mi><mi>l</mi><mi>T</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>db</mi><mi>m</mi><mi>T</mi></msubsup></mtd></mtr><mtr><mtd><msub><mi>k</mi><mi>u</mi></msub><msubsup><mi>b</mi><mi>n</mi><mi>T</mi></msubsup></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000722893470000012.GIF" wi="1071" he="315" /></maths>其中,d>0为电机到质心的距离;k<sub>u</sub>为一个与空气密度、旋翼的尺寸、形状、扭转角以及一些其他因素有关的参数,在此设为常值;其中W称为控制效率矩阵,表示如下:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>b</mi><mi>l</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><mn>0</mn><mo>,</mo><mo>-</mo><mfrac><msqrt><mn>3</mn></msqrt><mn>2</mn></mfrac><mo>,</mo><mo>-</mo><mfrac><msqrt><mn>3</mn></msqrt><mn>2</mn></mfrac><mo>,</mo><mn>0</mn><mo>,</mo><mfrac><msqrt><mn>3</mn></msqrt><mn>2</mn></mfrac><mo>,</mo><mfrac><msqrt><mn>3</mn></msqrt><mn>2</mn></mfrac><mo>)</mo></mrow><mi>T</mi></msup></mrow>]]></math><img file="FDA0000722893470000021.GIF" wi="700" he="186" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>b</mi><mi>m</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mo>-</mo><mn>1</mn><mo>,</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>)</mo></mrow><mi>T</mi></msup></mrow>]]></math><img file="FDA0000722893470000022.GIF" wi="584" he="156" /></maths>b<sub>n</sub>=(‑1,1,‑1,1,‑1,1)<sup>T</sup>b<sub>f</sub>=(1,1,1,1,1,1)<sup>T</sup>其中,W称为控制效率矩阵,通过伪逆矩阵法得到期望的旋翼升力,其计算公式如下:f=W<sup>T</sup>(WW<sup>T</sup>)<sup>‑1</sup>u   (9)步骤三:当六旋翼飞行器有一个电机停机时的控制重分配当六旋翼飞行器有一个电机i,i∈{1,2,3,4,5,6}停机后,由公式(8)导出虚拟控制量u和电机升力f<sub>[i]</sub>的映射关系如下:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>u</mi><mo>=</mo><msub><mi>W</mi><mrow><mo>[</mo><mi>i</mi><mo>]</mo></mrow></msub><msub><mi>f</mi><mrow><mo>[</mo><mi>i</mi><mo>]</mo></mrow></msub><mo>,</mo><msub><mi>W</mi><mrow><mo>[</mo><mi>i</mi><mo>]</mo></mrow></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>b</mi><mrow><mi>f</mi><mo>[</mo><mi>i</mi><mo>]</mo></mrow><mi>T</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>db</mi><mrow><mi>l</mi><mo>[</mo><mi>i</mi><mo>]</mo></mrow><mi>T</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>db</mi><mrow><mi>m</mi><mo>[</mo><mi>i</mi><mo>]</mo></mrow><mi>T</mi></msubsup></mtd></mtr><mtr><mtd><msub><mi>k</mi><mi>u</mi></msub><msubsup><mi>b</mi><mrow><mi>n</mi><mo>[</mo><mi>i</mi><mo>]</mo></mrow><mi>T</mi></msubsup></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000722893470000023.GIF" wi="1126" he="330" /></maths>其中,b<sub>f[i]</sub>,b<sub>l[i]</sub>,b<sub>m[i]</sub>,b<sub>n[i]</sub>和f<sub>[i]</sub>分别为从b<sub>f</sub>,b<sub>l</sub>,b<sub>m</sub>,b<sub>n</sub>和f中去除第i个元素后的向量;其中,W<sub>[i]</sub>为当电机i停机时的控制效率矩阵,由伪逆矩阵法可得从虚拟控制量u到旋翼升力f<sub>[i]</sub>的分配关系如下:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>f</mi><mrow><mo>[</mo><mi>i</mi><mo>]</mo></mrow></msub><mo>=</mo><msubsup><mi>W</mi><mrow><mo>[</mo><mi>i</mi><mo>]</mo></mrow><mi>T</mi></msubsup><msup><mrow><mo>(</mo><msub><mi>W</mi><mrow><mo>[</mo><mi>i</mi><mo>]</mo></mrow></msub><msubsup><mi>W</mi><mrow><mo>[</mo><mi>i</mi><mo>]</mo></mrow><mi>T</mi></msubsup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>u</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000722893470000024.GIF" wi="1079" he="89" /></maths>然而,当六旋翼飞行器有一个电机停机时全状态是不可控的,但是当放弃偏航通道,并且满足总推力<img file="FDA0000722893470000025.GIF" wi="274" he="132" />其中K为每个电机提供的最大升力,其俯仰和滚转通道依然可控,从而控制其安全着陆;所以由公式(11)得控制重分配为<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>f</mi><mrow><mo>[</mo><mi>i</mi><mo>]</mo></mrow></msub><mo>=</mo><msubsup><mi>W</mi><mrow><mo>[</mo><mi>i</mi><mo>]</mo></mrow><mi>T</mi></msubsup><msup><mrow><mo>(</mo><msub><mi>W</mi><mrow><mo>[</mo><mi>i</mi><mo>]</mo></mrow></msub><msubsup><mi>W</mi><mrow><mo>[</mo><mi>i</mi><mo>]</mo></mrow><mi>T</mi></msubsup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><mi>u</mi><mo>*</mo></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000722893470000026.GIF" wi="1089" he="93" /></maths>其中,u<sup>*</sup>=(F<sub>v</sub>,L<sub>v</sub>,M<sub>v</sub>,0)<sup>T</sup>;L<sub>v</sub>,M<sub>v</sub>和F<sub>v</sub>分别由公式(5)和(7)产生;其中F<sub>v</sub>和τ<sub>av</sub>=(L<sub>v</sub>,M<sub>v</sub>,N<sub>v</sub>)<sup>T</sup>分别为六旋翼期望的总推力和转动力矩。
地址 100191 北京市海淀区学院路37号