发明名称 飞行器大迎角运动四元数模型的故障诊断和容错控制方法
摘要 本发明给出了一种飞行器大迎角运动四元数模型的故障诊断和容错控制方法,该方法按照飞行器无故障时的模型和状态组合指标判别式进行故障诊断;飞行器出现故障时首先对故障分类,再针对确定的故障类别分析飞行器全机参数和飞行动力学特性的影响,然后重新对故障发生后突变系统建模,得到飞行器出现故障时新的模型描述;容错控制器的设计方法为设计出的控制量使得无故障模型输出达到军标要求,且该控制量和飞行器故障引起的非人为附加控制量共同组成的飞行器出现故障时的控制输入使得飞行器出现故障时新的模型输出稳定;并且可以直接进行飞行器大迎角飞行故障诊断并进行容错控制。
申请公布号 CN103149927B 申请公布日期 2015.04.08
申请号 CN201310095165.7 申请日期 2013.03.24
申请人 西安费斯达自动化工程有限公司 发明人 史忠科
分类号 G05B23/02(2006.01)I;G05B13/04(2006.01)I 主分类号 G05B23/02(2006.01)I
代理机构 代理人
主权项 一种飞行器大迎角运动四元数模型的故障诊断和容错控制方法,其特点是包括以下步骤:(a) 飞行器无故障时大迎角运动的四元数非线性模型为:<img file="922429dest_path_image001.GIF" wi="480" he="428" />(1)式中:四元数<img file="363905dest_path_image002.GIF" wi="145" he="29" />,<img file="727891dest_path_image003.GIF" wi="429" he="103" />状态变量<img file="901514dest_path_image004.GIF" wi="112" he="25" />分别表示飞行高度、飞行速度、气流迎角、侧滑角、滚转角速度、俯仰角速度、偏航角速度;<img file="163999dest_path_image005.GIF" wi="49" he="22" />分别指滚转角、俯仰角、偏航角;<img file="596118dest_path_image006.GIF" wi="16" he="18" />为重力加速度;<img file="898350dest_path_image007.GIF" wi="129" he="26" />分别表示动压、机翼面积、飞行器质量、绕机体轴系X轴的转动惯量、绕机体轴系Y轴的转动惯量、绕机体轴系Z轴的转动惯量、机体轴系XZ轴的转动惯量积;<img file="570771dest_path_image008.GIF" wi="226" he="26" />分别表示纵向力、侧向力和法向力;<img file="117290dest_path_image009.GIF" wi="234" he="25" />分别为滚转力矩、偏航力矩、俯仰力矩;<img file="41514dest_path_image010.GIF" wi="14" he="20" />为控制量;通常侧滑角小于<img file="924632dest_path_image011.GIF" wi="25" he="22" />,<img file="830271dest_path_image012.GIF" wi="41" he="22" />不会为零;(b) 故障检出方法为:当状态组合指标<img file="864086dest_path_image013.GIF" wi="37" he="20" />,则判定飞行器有故障;其中:<img file="513374dest_path_image014.GIF" wi="472" he="188" />(2)式中:<img file="dest_path_image015.GIF" wi="441" he="26" />分别为根据(1)式在时间<img file="253928dest_path_image016.GIF" wi="154" he="22" />计算得到的<img file="dest_path_image017.GIF" wi="222" he="26" />值;<img file="202905dest_path_image018.GIF" wi="433" he="26" />分别为在时间<img file="724016dest_path_image016.GIF" wi="154" he="22" />实际测量及计算四元数得到的<img file="849098dest_path_image017.GIF" wi="222" he="26" />值,<img file="dest_path_image019.GIF" wi="325" he="26" />为设定的检测阈值;<img file="709737dest_path_image020.GIF" wi="150" he="49" />为符号函数;<img file="891933dest_path_image021.GIF" wi="14" he="20" />表示第<img file="900340dest_path_image021.GIF" wi="14" he="20" />个采样点,对应采样时间为<img file="829113dest_path_image022.GIF" wi="45" he="20" />;<img file="544259dest_path_image023.GIF" wi="20" he="20" />为设定误差累计个数,<img file="962602dest_path_image024.GIF" wi="16" he="18" />为采样周期;(c) 飞行器出现故障时首先对故障分类,再针对确定的故障类别分析飞行器全机参数和飞行动力学特性的影响,然后重新对故障发生后突变系统建模,得到飞行器出现故障时大迎角运动的四元数非线性模型描述为:<img file="192726dest_path_image025.GIF" wi="516" he="488" />(3)式中:<img file="984577dest_path_image026.GIF" wi="145" he="26" />别表示飞行器出现故障时的机翼面积、飞行器质量、绕机体轴系X轴的转动惯量、绕机体轴系Y轴的转动惯量、绕机体轴系Z轴的转动惯量、机体轴系XZ轴的转动惯量积;<img file="616547dest_path_image027.GIF" wi="265" he="26" />分别表示飞行器出现故障时的纵向力、侧向力和法向力;<img file="205791dest_path_image028.GIF" wi="89" he="26" />、<img file="188791dest_path_image029.GIF" wi="90" he="26" />和<img file="521683dest_path_image030.GIF" wi="92" he="26" />飞行器出现故障时的滚转力矩、偏航力矩、俯仰力矩;<img file="8159dest_path_image031.GIF" wi="77" he="52" />飞行器出现故障时的控制输入,飞行器故障引起的非人为附加控制量<img file="33884dest_path_image032.GIF" wi="29" he="26" />;(d) 容错控制器的设计方法为:设计出的控制量<img file="238600dest_path_image010.GIF" wi="14" he="20" />使得非线性模型(1)式稳定裕度达到军标要求,且该控制量<img file="372254dest_path_image010.GIF" wi="14" he="20" />和飞行器故障引起的非人为附加控制量<img file="713237dest_path_image032.GIF" wi="29" he="26" />共同组成的飞行器出现故障时的控制输入<img file="909863dest_path_image031.GIF" wi="77" he="52" />使得非线性模型(3)式输出稳定。
地址 710075 陕西省西安市高新区科技路金桥国际广场12101号