发明名称 一种基于角加速度传感器的多操纵面无人机闭环控制分配方法
摘要 本发明公开了一种基于角加速度传感器的多操纵面无人机闭环控制分配方法,可应用于具有冗余多操纵面无人机的飞行控制。通过无人机角速度指令和测量到的角速度信息得到期望角加速度指令,通过角加速度传感器获得无人机的角加速度,得到期望力矩增量与实际力矩增量的误差,实时地对该力矩增量误差进行控制分配解算,得到各个操纵面的偏转角度。本发明提供的方法利用了无人机角加速度传感器得到角加速度反馈信息,从而构成了闭环控制分配结构,在操纵面正常工作和发生故障时,均能达到较好的控制分配效果,并且计算简单、实现方便。
申请公布号 CN102736631A 申请公布日期 2012.10.17
申请号 CN201210191778.6 申请日期 2012.06.11
申请人 北京航空航天大学 发明人 王宏伦;盖文东;樊国玮
分类号 G05D1/10(2006.01)I;G01P15/08(2006.01)I 主分类号 G05D1/10(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 赵文利
主权项 1.一种基于角加速度传感器的多操纵面无人机闭环控制分配方法,其特征在于,具体包括以下几个步骤:步骤一:通过无人机上的角速度传感器测得当前的无人机角速度,根据期望的角速度指令,得到无人机的期望角加速度;角速度传感器测得的角速度为[p,q,r]<sup>T</sup>,其中,p表示滚转角速度,q表示俯仰角速度,r表示偏航角速度,期望的角速度指令为[p<sub>c</sub>,q<sub>c</sub>,r<sub>c</sub>]<sup>T</sup>,其中,p<sub>c</sub>表示滚转角速度指令,q<sub>c</sub>表示俯仰角速度指令,r<sub>c</sub>表示偏航角速度指令,则无人机的期望角加速度为:<maths num="0001"><![CDATA[<math><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mover><mi>p</mi><mo>&CenterDot;</mo></mover><mi>d</mi></msub></mtd></mtr><mtr><mtd><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mi>d</mi></msub></mtd></mtr><mtr><mtd><msub><mover><mi>r</mi><mo>&CenterDot;</mo></mover><mi>d</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>&omega;</mi><mi>p</mi></msub></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>&omega;</mi><mi>q</mi></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>&omega;</mi><mi>r</mi></msub></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>p</mi><mi>c</mi></msub><mo>-</mo><mi>p</mi></mtd></mtr><mtr><mtd><msub><mi>q</mi><mi>c</mi></msub><mo>-</mo><mi>q</mi></mtd></mtr><mtr><mtd><msub><mi>r</mi><mi>c</mi></msub><mo>-</mo><mi>r</mi></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>式中,ω<sub>p</sub>,ω<sub>q</sub>和ω<sub>r</sub>为滚转、俯仰和偏航通道的频带或增益,<img file="FDA00001750937300012.GIF" wi="43" he="43" />表示期望的滚转角加速度,<img file="FDA00001750937300013.GIF" wi="36" he="42" />表示期望的俯仰角加速度,<img file="FDA00001750937300014.GIF" wi="29" he="42" />表示期望的偏航角加速度;步骤二:通过无人机上的角加速度传感器测得当前的无人机角加速度<img file="FDA00001750937300015.GIF" wi="169" he="60" />其中,<img file="FDA00001750937300016.GIF" wi="25" he="40" />表示当前的滚转角加速度,<img file="FDA00001750937300017.GIF" wi="19" he="39" />表示当前的俯仰角加速度,<img file="FDA00001750937300018.GIF" wi="19" he="31" />表示当前的偏航角加速度,并与无人机的期望角加速度<img file="FDA00001750937300019.GIF" wi="193" he="60" />相减,再与转动惯量矩阵J<sub>c</sub>相乘,得到无人机的力矩增量误差ΔM<sub>e</sub>,即期望力矩增量与实际力矩增量之差:<maths num="0002"><![CDATA[<math><mrow><mi>&Delta;</mi><msub><mi>M</mi><mi>e</mi></msub><mo>=</mo><msubsup><mi>J</mi><mi>c</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mover><mi>p</mi><mo>&CenterDot;</mo></mover><mi>d</mi></msub><mo>-</mo><mover><mi>p</mi><mo>&CenterDot;</mo></mover></mtd></mtr><mtr><mtd><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mi>d</mi></msub><mo>-</mo><mover><mi>q</mi><mo>&CenterDot;</mo></mover></mtd></mtr><mtr><mtd><msub><mover><mi>r</mi><mo>&CenterDot;</mo></mover><mi>d</mi></msub><mo>-</mo><mover><mi>r</mi><mo>&CenterDot;</mo></mover></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>式中,转动惯量矩阵J<sub>c</sub>为<maths num="0003"><![CDATA[<math><mrow><msub><mi>J</mi><mi>c</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>J</mi><mi>z</mi></msub><mo>/</mo><mi>&Gamma;</mi></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>J</mi><mi>xz</mi></msub><mo>/</mo><mi>&Gamma;</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn><mo>/</mo><msub><mi>J</mi><mi>y</mi></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>J</mi><mi>xz</mi></msub><mo>/</mo><mi>&Gamma;</mi></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>J</mi><mi>x</mi></msub><mo>/</mo><mi>&Gamma;</mi></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>式中,J<sub>x</sub>表示绕机体轴x轴的转动惯量,J<sub>y</sub>表示绕机体轴y轴的转动惯量,J<sub>z</sub>表示绕机体轴z轴的转动惯量,J<sub>xz</sub>为惯性积,<img file="FDA000017509373000112.GIF" wi="265" he="51" />步骤三:采用控制分配方法对力矩增量误差ΔM<sub>e</sub>进行分配,得到不同执行机构k时刻的控制增量U<sub>v</sub>(k);步骤四:将当前时刻计算得到的执行机构控制增量U<sub>v</sub>(k)与存储的前一时刻执行机构控制输入U<sub>v</sub>(k-1)相加,得到当前时刻不同执行机构的控制输入U(k);U(k)=U(k-1)+U<sub>v</sub>(k)式中,U(k)=[u<sub>1</sub>(k),…,u<sub>m</sub>(k)]<sup>T</sup>,m表示执行机构的数量,u<sub>i</sub>(k)i=1,…,m表示第i个执行机构k时刻的控制输入;步骤五:将不同执行机构的控制输入U(k)分别接入相应的执行机构,从而得到不同操纵面的偏转角度D(k);设有m个执行机构,则这些执行机构的模型为:<img file="FDA00001750937300021.GIF" wi="607" he="299" />式中,K<sub>i</sub>表示第i个执行机构的频带特性,i=1,…,m,s表示Laplace算符,D(k)=[δ<sub>1</sub>(k),…,δ<sub>m</sub>(k)]<sup>T</sup>,m表示执行机构的数量,δ<sub>i</sub>(k)i=1,…,m表示第i个执行机构k时刻的输出,即第i个操纵面在k时刻的偏转角度;步骤六:在后续的飞行控制过程中,不断重复前述五个步骤,从而实时得到无人机不同操纵面的偏转角角度D(k);将得到的不同执行机构控制输入U(k),以及不同操纵面的偏转角角度D(k),均通过测控链路下传显示,供操控人员监控参考。
地址 100191 北京市海淀区学院路37号