CN111859555A - 输入受限的航空发动机最大推力状态鲁棒容错控制器 - Google Patents

输入受限的航空发动机最大推力状态鲁棒容错控制器 Download PDF

Info

Publication number
CN111859555A
CN111859555A CN202010543149.XA CN202010543149A CN111859555A CN 111859555 A CN111859555 A CN 111859555A CN 202010543149 A CN202010543149 A CN 202010543149A CN 111859555 A CN111859555 A CN 111859555A
Authority
CN
China
Prior art keywords
engine
fault
input
vector
maximum thrust
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202010543149.XA
Other languages
English (en)
Inventor
张猛
缑林峰
刘志丹
蒋宗霆
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN202010543149.XA priority Critical patent/CN111859555A/zh
Publication of CN111859555A publication Critical patent/CN111859555A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明提出一种输入受限的航空发动机最大推力状态鲁棒容错控制器。最大推力状态鲁棒控制器组容错控制模块产生控制向量v并输出给输入限制模块,输入限制模块产生限制后的控制输入向量u并输出给航空发动机本体,气路部件故障诊断模块诊断发动机的气路部件故障;最大推力状态鲁棒控制器组容错控制模块计算得到适应的鲁棒控制器,并产生控制输入向量u。本发明能够在保证发动机安全工作的前提下,在发动机机气路部件故障的情况下在最大推力状态依旧对真实发动机进行良好控制,具有较强的鲁棒性,最大限度的提高发动机在最大推力状态的性能,使发动机气路部件故障时在最大推力状态不仅稳定工作,并且具有最优的性能,提高战斗机的机动性能。

Description

输入受限的航空发动机最大推力状态鲁棒容错控制器
技术领域
本发明涉及航空发动机控制技术领域,尤其涉及一种输入受限的航空发动机最大推力状态鲁棒容错控制器。
背景技术
航空发动机是一个复杂的非线性动力学系统,其控制系统容易受到工作条件,发动机性能下降,环境条件变化的影响,并且很难事先知道外部干扰和测量噪声的影响。由于飞机发动机的工作过程非常复杂,难以建立准确的数学模型,所以数学模型与实际系统之间总是存在差异。因此,有必要设计一种鲁棒控制器,用于在外部干扰信号,噪声干扰,未建模的动态特性和参数变化的情况下稳定航空发动机控制系统,并具有良好的性能。
战斗机由于需要实现高机动性,发动机的最大推力状态的性能及安全性至关重要。传统的鲁棒控制器虽然可以对发动机在最大推力状态实现稳定控制。然而,现代战机对航空发动机性能的要求不断提高,其结构也越来越复杂,并且由于发动机工作环境的恶劣多变,发动机故障约占飞机总故障的1/3。其中,气路部件故障占发动机总体故障的90%以上,其维护费用占发动机总体维护费用的60%。为了保证发动机安全工作并使故障发动机提供足够的性能来保证飞机安全飞行或具有高的机动性,必须对故障的发动机性能进行恢复,并且对发动机进行容错控制,保证控制系统正常稳定工作且具有良好的性能。因此,研究发动机气路部件故障容错控制方法具有重要意义。
传统的气路部件故障容错控制方法在航空发动机出现气路部件故障时通过修正控制规律,使得发动机的推力与油门杆始终匹配,有效的保证了发动机的推力。然而,这些设计方法并没有解决当前控制器和发动机模型不匹配从而导致控制系统性能下降甚至不稳定的问题。当发动机发生气路部件故障时,发动机在同一工作点的线性化模型也会发生较大变化。因此,根据正常状态的发动机模型设计的控制器一般无法保证气路部件故障时发动机的性能,甚至无法保证控制系统的闭环稳定。
此外,过大的控制输入会导致发动机损坏,因此我们需要考虑控制输入受限的控制器的设计。
发明内容
为解决现有技术存在的问题,本发明提出一种输入受限的航空发动机最大推力状态鲁棒容错控制器,具有较强的鲁棒性,并且在最大推力状态能够在发动机机气路部件故障的情况下依旧对真实发动机进行良好控制,保证发动机安全工作,充分发挥发动机最大推力状态的性能,提高飞机的安全性和性能,提高战斗机的机动性。并且考虑控制输入受限,保证发动机安全工作。
本发明的技术方案为:
所述一种输入受限的航空发动机最大推力状态鲁棒容错控制器,其特征在于:包括最大推力状态鲁棒控制器组容错控制模块、输入限制模块和气路部件故障诊断模块;
其中最大推力状态鲁棒控制器组容错控制模块、输入限制模块、气路部件故障诊断模块与航空发动机本体以及航空发动机上的若干传感器组成气路部件故障调度控制回路;
所述最大推力状态鲁棒控制器组容错控制模块产生控制向量v并输出给输入限制模块,输入限制模块产生限制后的控制输入向量u并输出给航空发动机本体,传感器得到航空发动机测量参数y;控制输入向量u以及测量参数y共同输入到气路部件故障诊断模块,气路部件故障诊断模块诊断发动机的气路部件故障情况得到航空发动机的健康参数h,并输出到最大推力状态鲁棒控制器组容错控制模块;
所述输入限制模块限制了控制输入向量的幅值,避免给发动机过大的控制输入导致发动机损坏;
所述最大推力状态鲁棒控制器组容错控制模块内设计有若干鲁棒控制器,所述鲁棒控制器是利用若干线性不确定性发动机模型而分别设计得到的,所述线性不确定性发动机模型是对航空发动机最大推力状态下的、不同气路部件故障下的航空发动机非线性模型进行线性化后再加入摄动块得到的;
所述最大推力状态鲁棒控制器组容错控制模块根据输入的健康参数h,利用内部设计的若干鲁棒控制器计算得到适应的鲁棒控制器,该鲁棒控制器根据参考输入r和测量参数y的差值e产生控制输入向量u。
进一步的,所述最大推力状态鲁棒控制器组容错控制模块内设计若干鲁棒控制器的过程为:在航空发动机最大推力状态对包含健康参数的发动机非线性模型进行线性化得到含有健康参数的线性化模型,通过调整健康参数的值,得到分别在发动机无气路部件故障和特定气路部件故障处的11个线性化模型再加入摄动块得到11个线性不确定性发动机模型,并对这11个线性不确定性发动机模型分别设计相应的鲁棒控制器从而组成最大推力状态鲁棒控制器组。
进一步的,所述气路部件故障诊断模块中包括非线性机载发动机模型和线性化卡尔曼滤波器;
所述非线性机载发动机模型为带健康参数的发动机非线性模型:
Figure BDA0002539641440000031
y=g(x,u,h)
其中
Figure BDA0002539641440000032
为控制输入向量,
Figure BDA0002539641440000033
为状态向量,
Figure BDA0002539641440000034
为输出向量,
Figure BDA0002539641440000035
为健康参数向量,f(·)为表示系统动态的n维可微非线性向量函数,g(·)为产生系统输出的m维可微非线性向量函数;非线性机载发动机模型输入为控制输入向量u以及上一周期的健康参数h,其输出的健康稳态参考值(xaug,NOBEM,yNOBEM)作为线性化卡尔曼滤波器当前周期的估计初始值;
所述线性化卡尔曼滤波器的输入为测量参数y以及非线性机载发动机模型输出的健康稳态参考值(xaug,NOBEM,yNOBEM),根据公式
Figure BDA0002539641440000036
计算得到当前周期的发动机的健康参数h;其中
Figure BDA0002539641440000037
K为卡尔曼滤波的增益,满足
Figure BDA0002539641440000038
P为Ricati方程
Figure BDA0002539641440000039
的解;系数Aaug和Caug根据公式
Figure BDA00025396414400000310
确定,而A、C、L、M是将健康参数h看作发动机的控制输入,并对非线性机载发动机模型在健康稳态参考点处进行线性化得到的反映发动机性能退化的增广线性状态变量模型
Figure BDA0002539641440000041
的系数:
Figure BDA0002539641440000042
Figure BDA0002539641440000043
w为系统噪声,v为测量噪声,相应的协方差矩阵为对角阵Q和R。
进一步的,所述最大推力状态鲁棒控制器组容错控制模块根据输入的健康参数h插值得到的适应的鲁棒控制器。
进一步的,所述最大推力状态鲁棒控制器组容错控制模块根据航空发动机最大推力状态对应发动机无部件故障的控制器K0,各种典型部件故障Δhbase_j的控制器
Figure BDA0002539641440000044
Δhbase_j表示向量Δh的第j个元素的值为Δhbase,其他元素的值为0,即Δhbase_j表示10种不同的部件故障,例如Δhbase_1表示风扇发生了故障且风扇的效率变化量为Δhbase。根据公式
Figure BDA0002539641440000045
计算得到航空发动机最大推力状态处发动机的当前部件故障程度(健康参数为h)下的鲁棒控制器K(式中Δhj为向量Δh的第j个元素;仅考虑||Δh||≤||Δhmax||的发动机气路部件故障情况,当||Δh||>||Δhmax||发动机已失效)。
进一步的,所述输入限制模块采用多维矩形饱和函数,控制输入向量u为:
Figure BDA0002539641440000051
Figure BDA0002539641440000052
其中v1和vm为控制向量v的元素,v1,max和vm,max为控制向量v对应元素的限幅值。
进一步的,所述测量参数包括进气道出口、风扇出口、压气机出口、高压涡轮后、低压涡轮后的温度和压力,风扇转速和压气机转速。
有益效果
与现有技术相比较,本发明的输入受限的航空发动机最大推力状态鲁棒容错控制器利用传统增益调度控制器中固有的模块,通过新增气路部件故障诊断模块,并对最大推力状态鲁棒控制器组容错控制模块进行了改进,新增了发动机不同气路部件故障下的多组鲁棒控制器。气路部件故障诊断模块通过健康参数的可靠估计实现了气路部件故障的准确判断,在保证发动机安全工作的前提下,实现发动机在最大推力状态气路部件故障时的鲁棒容错控制,具有较强的鲁棒性,最大限度的提高发动机在最大推力状态的性能,使发动机气路部件故障时在最大推力状态不仅稳定工作,并且具有最优的性能,提高战斗机的机动性能。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1是本发明输入受限的航空发动机最大推力状态鲁棒容错控制器的结构简图;
图2是本实施例气路部件故障调度控制回路中气路部件故障诊断模块的结构示意图;
图3是本实施例气路部件故障诊断模块中卡尔曼滤波器的结构示意图。
具体实施方式
战斗机由于需要实现高机动性,发动机的最大推力状态的性能至关重要。传统的鲁棒控制器虽然可以对发动机在最大推力状态实现稳定控制,然而,航空发动机在运行过程中由于自然磨损、腐蚀、积垢以及热蠕变等因素都会引起气路部件性能发生蜕化,并且当性能蜕化到一定程度会引发故障;另外,由于外物吸入引起的损伤、机械疲劳断裂等原因,也会引起气路部件故障的发生。前者故障发生的过程较为缓慢,而后者故障发生的过程是迅速的。当发动机气路部件发生故障却并未失效时,此时发动机的部分性能将会严重偏离额定状态。以涡轮部件为例,当其发生故障后,其工作效率将会下降,即将具有高温、高压的燃气转换为机械能的能力将会有所降低,但仍能为风扇或压气机部件提供相应的动力,使其工作在新的平衡状态下。此时发动机也已较大偏离原始状态。气路部件故障会导致发动机设计时所建立的非线性模型和气路部件故障时的真实发动机严重不匹配,进而导致根据该非线性模型设计的增益调度控制器无法对气路部件故障后的发动机进行良好的控制,严重降低了发动机的性能,甚至不能保证控制系统的稳定性,无法保证发动机安全工作。针对这一问题,下面给出本发明的分析研究过程。
1、发动机气路部件故障诊断
气路部件故障会导致部件对应的特性参数发生变化。发动机气路部件故障最终表征在不同转子部件的工作效率和流通量的改变上,即可以从风风扇、压气机、主燃烧、高压涡轮和低压涡轮部件的效率系数或者流量系数的变化来揭示发动机故障位置以及故障程度,风扇、压气机、主燃烧室、高压涡轮和低压涡轮部件的效率系数或流量系数被称为健康参数。
基于部件法,建立带健康参数的发动机非线性模型
Figure BDA0002539641440000061
y=g(x,u,h)
其中
Figure BDA0002539641440000062
为控制输入向量,
Figure BDA0002539641440000063
为状态向量,
Figure BDA0002539641440000064
为输出向量,
Figure BDA0002539641440000065
为健康参数向量,f(·)为表示系统动态的n维可微非线性向量函数,g(·)为产生系统输出的m维可微非线性向量函数。
将健康参数h看作发动机的控制输入,采用小扰动法或拟合法对发动机非线性模型在健康稳态参考点处进行线性化。
Figure BDA0002539641440000071
其中
A′=A,B′=(B L),C′=C,
D′=(D M),Δu′=(Δu Δh)T
w为系统噪声,v为测量噪声,h为健康参数,Δh=h-h0;上述w与v皆为不相关的高斯白噪声,其均值均为0,协方差矩阵为对角阵Q和R,即满足条件如下:
E(w)=0 E[wwT]=Q
E(v)=0 E[vvT]=R
Δ表示该参数的变化量,h0表示发动机初始状态健康参数。
进一步得到了反映发动机性能退化的增广线性状态变量模型
Figure BDA0002539641440000072
其中系数矩阵可由下式得到:
Figure BDA0002539641440000073
Figure BDA0002539641440000074
这些系数在发动机不同的工作状态具有不同的值。
实际上,健康参数很难测量,甚至不可能测量,而发动机各部分的压力、温度、转速等参数比较容易通过测量得到,通常称为“测量参数”,主要包括进气道出口、风扇出口、压气机出口、高压涡轮后、低压涡轮后的温度和压力,风扇转速和压气机转速。当发动机工作环境不发生变化时,健康参数的变化会引起被测参数的相应变化,二者之间存在气动热力学关系。因此,可以设计最优估计滤波器,通过测量参数来实现健康参数的最优估计。
对于渐变型部件故障,对应故障部件健康参数变化缓慢,因此在进行单次故障诊断的时间周期内,可以认为满足
Figure BDA0002539641440000081
对于突变型部件故障,更加关心的是发动机在故障发生后再次稳定工作时部件故障的严重程度,发动机再次稳定工作后,故障部件的健康参数变化依旧满足
Figure BDA0002539641440000082
将健康参数进一步转化为状态变量,可以得到
Figure BDA0002539641440000083
其中
Figure BDA0002539641440000084
Figure BDA0002539641440000085
建立的气路部件故障诊断模块主要由两部分组成,一部分是基于健康参数的非线性机载发动机模型,另一部分是分段线性卡尔曼滤波器。基本工作原理是将非线性机载发动机模型的输出作为分段线性卡尔曼滤波器的稳态参考值,并扩展健康参数,通过分段线性卡尔曼滤波器进行在线实时估计,最后反馈给非线性机载发动机模型进行在线实时更新,实现对实际发动机的实时跟踪。
卡尔曼估计方程为:
Figure BDA0002539641440000086
K为卡尔曼滤波的增益,满足
Figure BDA0002539641440000087
P为Ricati方程
Figure BDA0002539641440000088
的解;利用非线性机载模型输出的健康稳态参考值(xaug,NOBEM,yNOBEM)作为式
Figure BDA0002539641440000089
的初值,可得计算公式:
Figure BDA0002539641440000091
根据该计算公式可以得到发动机的健康参数h,实现发动机的气路部件故障诊断。
2.具有健康参数的不确定模型的鲁棒控制器设计
任何实际系统都不可避免地存在不确定性,它可以分为两类:扰动信号和模型不确定性。扰动信号包括干扰、噪声等。模型的不确定性代表了数学模型与实际对象之间的差异。
模型不确定性可能有几个原因:线性模型中总有一些参数是有误差的;线性模型中的参数可能由于非线性或工作条件的变化而变化;建模时人为的简化;由于磨损等因素发动机性能的退化。
不确定性可能会对控制系统的稳定性和性能产生不利影响。
实际的发动机和标称模型(标称模型是一个常规的不带健康参数的发动机非线性模型)之间的误差可以表示为一个摄动块Δ。在标称模型加入摄动块建立发动机不确定模型
Figure BDA0002539641440000092
Figure BDA0002539641440000093
它也可以表示为
G(s)=[I+Δ(s)]Gnom(s)
最后根据不确定模型利用传统的鲁棒控制器设计方法设计鲁棒控制器。
3、控制器的插值
这部分说明了图1中的最大推力状态鲁棒控制器组容错控制模块通过健康参数调度线性插值获得相应的鲁棒控制器的调度计算原理。
在发动机最大推力状态,分别在发动机正常状态和各种典型部件故障Δhbase_j状态下设计一系列线性鲁棒控制器,对发动机进行控制。这将产生图1中的最大推力状态鲁棒控制器组容错控制模块中的控制器
Figure BDA0002539641440000101
然后根据健康参数h对控制器进行插值,继而使用得到的插值控制器来控制系统。
根据发动机最大推力状态对应发动机无部件故障的控制器K0,各种典型部件故障Δhbase_j的控制器
Figure BDA0002539641440000102
Δhbase_j表示向量Δh的第j个元素的值为Δhbase,其他元素的值为0,即Δhbase_j表示10种不同的部件故障,例如Δhbase_1表示风扇发生了故障且风扇的效率变化量为Δhbase。通过线性插值可以得到发动机最大推力状态处气路部件故障h处的控制器
Figure BDA0002539641440000103
并对发动机进行有效控制。
4.系统的输入限制
请参考图1,图1中的输入限制模块是为了建模系统控制输入上的物理限制,使用了多维矩形饱和函数。限制航空发动机控制的输入,尤其是对于燃油流量输入。多维饱和函数还可以处理其他控制输入的限制,包括尾喷嘴的喉部面积。该函数是一个多维矩形饱和函数,定义为
Figure BDA0002539641440000104
其中v1和vm为控制向量v的元素,v1,max和vm,max为控制向量v对应元素的限幅值。对于所有的
Figure BDA0002539641440000106
下式给出了sat(·)
Figure BDA0002539641440000105
基于上述过程,下面给出本实施例中提出的一种输入受限的航空发动机最大推力状态鲁棒容错控制器,如图1所示,主要包括最大推力状态鲁棒控制器组容错控制模块、输入限制模块和气路部件故障诊断模块。
其中最大推力状态鲁棒控制器组容错控制模块、输入限制模块、气路部件故障诊断模块与航空发动机本体以及航空发动机上的若干传感器组成气路部件故障调度控制回路10。
所述最大推力状态鲁棒控制器组容错控制模块产生控制向量v并输出给输入限制模块,输入限制模块产生限制后的控制输入向量u并输出给航空发动机本体,传感器得到航空发动机测量参数y;控制输入向量u以及测量参数y共同输入到气路部件故障诊断模块,气路部件故障诊断模块解算得到航空发动机的健康参数h,并输出到最大推力状态鲁棒控制器组容错控制模块。
所述输入限制模块限制了控制输入向量的幅值,避免给发动机过大的控制输入导致发动机损坏。
优选的一种具体实现方式,所述输入限制模块采用多维矩形饱和函数,控制输入向量u为:
Figure BDA0002539641440000111
Figure BDA0002539641440000112
其中v1和vm为控制向量v的元素,v1,max和vm,max为控制向量v对应元素的限幅值。
所述最大推力状态鲁棒控制器组容错控制模块内设计有若干鲁棒控制器,所述鲁棒控制器是利用若干线性不确定性发动机模型而分别设计得到的,所述线性不确定性发动机模型是对航空发动机最大推力状态下的、不同气路部件故障下的航空发动机非线性模型进行线性化后再加入摄动块得到的。
优选的一种具体实现方式,可以通过以下过程得到设计若干鲁棒控制器:在航空发动机最大推力状态对包含健康参数的发动机非线性模型进行线性化得到含有健康参数的线性化模型,通过调整健康参数的值,得到分别在发动机无气路部件故障和特定气路部件故障处的11个线性化模型再加入摄动块得到11个线性不确定性发动机模型,并对这11个线性不确定性发动机模型分别设计相应的鲁棒控制器从而组成最大推力状态鲁棒控制器组。
所述最大推力状态鲁棒控制器组容错控制模块根据输入的健康参数h,利用内部设计的若干鲁棒控制器计算得到适应的鲁棒控制器,该鲁棒控制器根据参考输入r和测量参数y的差值e产生控制输入向量u。
优选的一种具体实现方式,可以根据输入的健康参数h插值得到的适应的鲁棒控制器:
根据发动机最大推力状态对应发动机无部件故障的控制器K0,各种典型部件故障Δhbase_j的控制器
Figure BDA0002539641440000121
Δhbase_j表示向量Δh的第j个元素的值为Δhbase,其他元素的值为0,即Δhbase_j表示10种不同的部件故障,例如Δhbase_1表示风扇发生了故障且风扇的效率变化量为Δhbase。通过线性插值可以得到发动机最大推力状态处气路部件故障h处的控制器
Figure BDA0002539641440000122
并对发动机进行有效控制。
所述气路部件故障诊断模块中包括非线性机载发动机模型和线性化卡尔曼滤波器。
所述非线性机载发动机模型为带健康参数的发动机非线性模型:
Figure BDA0002539641440000123
y=g(x,u,h)
其中
Figure BDA0002539641440000124
为控制输入向量,
Figure BDA0002539641440000125
为状态向量,
Figure BDA0002539641440000126
为输出向量,
Figure BDA0002539641440000127
为健康参数向量,f(·)为表示系统动态的n维可微非线性向量函数,g(·)为产生系统输出的m维可微非线性向量函数;非线性机载发动机模型输入为控制输入向量u以及上一周期的健康参数h,其输出的健康稳态参考值(xaug,NOBEM,yNOBEM)作为线性化卡尔曼滤波器当前周期的估计初始值。
所述线性化卡尔曼滤波器的输入为测量参数y以及非线性机载发动机模型输出的健康稳态参考值(xaug,NOBEM,yNOBEM),根据公式
Figure BDA0002539641440000131
计算得到当前周期的发动机的健康参数h。
其中
Figure BDA0002539641440000132
K为卡尔曼滤波的增益,满足
Figure BDA0002539641440000133
P为Ricati方程
Figure BDA0002539641440000134
的解;系数Aaug和Caug根据公式
Figure BDA0002539641440000135
确定,而A、C、L、M是将健康参数h看作发动机的控制输入,并对非线性机载发动机模型在健康稳态参考点处进行线性化得到的反映发动机性能退化的增广线性状态变量模型
Figure BDA0002539641440000136
的系数:
Figure BDA0002539641440000137
Figure BDA0002539641440000138
w为系统噪声,v为测量噪声,相应的协方差矩阵为对角阵Q和R。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (7)

1.一种输入受限的航空发动机最大推力状态鲁棒容错控制器,其特征在于:包括最大推力状态鲁棒控制器组容错控制模块、输入限制模块和气路部件故障诊断模块;
其中最大推力状态鲁棒控制器组容错控制模块、输入限制模块、气路部件故障诊断模块与航空发动机本体以及航空发动机上的若干传感器组成气路部件故障调度控制回路;
所述最大推力状态鲁棒控制器组容错控制模块产生控制向量v并输出给输入限制模块,输入限制模块产生限制后的控制输入向量u并输出给航空发动机本体,传感器得到航空发动机测量参数y;控制输入向量u以及测量参数y共同输入到气路部件故障诊断模块,气路部件故障诊断模块解算得到航空发动机的健康参数h,并输出到最大推力状态鲁棒控制器组容错控制模块;
所述输入限制模块限制了控制输入向量的幅值,避免给发动机过大的控制输入导致发动机损坏;
所述最大推力状态鲁棒控制器组容错控制模块内设计有若干鲁棒控制器,所述鲁棒控制器是利用若干线性不确定性发动机模型而分别设计得到的,所述线性不确定性发动机模型是对航空发动机最大推力状态下的、不同气路部件故障下的航空发动机非线性模型进行线性化后再加入摄动块得到的;
所述最大推力状态鲁棒控制器组容错控制模块根据输入的健康参数h,利用内部设计的若干鲁棒控制器计算得到适应的鲁棒控制器,该鲁棒控制器根据参考输入r和测量参数y的差值e产生控制输入向量u。
2.根据权利要求1所述一种输入受限的航空发动机最大推力状态鲁棒容错控制器,其特征在于:所述最大推力状态鲁棒控制器组容错控制模块内设计若干鲁棒控制器的过程为:在航空发动机最大推力状态对包含健康参数的发动机非线性模型进行线性化得到含有健康参数的线性化模型,通过调整健康参数的值,得到分别在发动机无气路部件故障和特定气路部件故障处的11个线性化模型再加入摄动块得到11个线性不确定性发动机模型,并对这11个线性不确定性发动机模型分别设计相应的鲁棒控制器从而组成最大推力状态鲁棒控制器组。
3.根据权利要求1或2所述一种输入受限的航空发动机最大推力状态鲁棒容错控制器,其特征在于:所述最大推力状态鲁棒控制器组容错控制模块根据输入的健康参数h插值得到适应的鲁棒控制器。
4.根据权利要求3所述一种输入受限的航空发动机最大推力状态鲁棒容错控制器,其特征在于:根据发动机最大推力状态对应发动机无部件故障的控制器K0,各种典型部件故障Δhbase_j的控制器
Figure RE-FDA0002687562880000021
Δhbase_j表示向量Δh的第j个元素的值为Δhbase,其他元素的值为0,即Δhbase_j表示10种不同的部件故障,例如Δhbase_1表示风扇发生了故障且风扇的效率变化量为Δhbase。通过线性插值可以得到发动机最大推力状态处气路部件故障h处的控制器
Figure RE-FDA0002687562880000022
并对发动机进行有效控制。
5.根据权利要求1所述一种输入受限的航空发动机最大推力状态鲁棒容错控制器,其特征在于:所述气路部件故障诊断模块中包括非线性机载发动机模型和线性化卡尔曼滤波器;
所述非线性机载发动机模型为带健康参数的发动机非线性模型:
Figure RE-FDA0002687562880000023
y=g(x,u,h)
其中
Figure RE-FDA0002687562880000024
为控制输入向量,
Figure RE-FDA0002687562880000025
为状态向量,
Figure RE-FDA0002687562880000026
为输出向量,
Figure RE-FDA0002687562880000027
为健康参数向量,f(·)为表示系统动态的n维可微非线性向量函数,g(·)为产生系统输出的m维可微非线性向量函数;非线性机载发动机模型输入为控制输入向量u以及上一周期的健康参数h,其输出的健康稳态参考值(xaug,NOBEM,yNOBEM)作为线性化卡尔曼滤波器当前周期的估计初始值;
所述线性化卡尔曼滤波器的输入为测量参数y以及非线性机载发动机模型输出的健康稳态参考值(xaug,NOBEM,yNOBEM),根据公式
Figure RE-FDA0002687562880000031
计算得到当前周期的发动机的健康参数h;其中
Figure RE-FDA0002687562880000032
K为卡尔曼滤波的增益,满足
Figure RE-FDA0002687562880000033
P为Ricati方程
Figure RE-FDA0002687562880000034
的解;系数Aaug和Caug根据公式
Figure RE-FDA0002687562880000035
确定,而A、C、L、M是将健康参数h看作发动机的控制输入,并对非线性机载发动机模型在健康稳态参考点处进行线性化得到的反映发动机性能退化的增广线性状态变量模型
Figure RE-FDA0002687562880000036
的系数:
Figure RE-FDA0002687562880000037
Figure RE-FDA0002687562880000038
w为系统噪声,v为测量噪声,相应的协方差矩阵为对角阵Q和R。
6.根据权利要求1所述一种输入受限航空发动机增益调度容错控制器,其特征在于:所述输入限制模块采用多维矩形饱和函数,控制输入向量u为:
Figure RE-FDA0002687562880000041
Figure RE-FDA0002687562880000042
其中v1和vm为控制向量v的元素,v1,max和vm,max为控制向量v对应元素的限幅值。
7.根据权利要求1所述一种输入受限的航空发动机最大推力状态鲁棒容错控制器,其特征在于:所述测量参数包括进气道出口、风扇出口、压气机出口、高压涡轮后、低压涡轮后的温度和压力,风扇转速和压气机转速。
CN202010543149.XA 2020-06-15 2020-06-15 输入受限的航空发动机最大推力状态鲁棒容错控制器 Pending CN111859555A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010543149.XA CN111859555A (zh) 2020-06-15 2020-06-15 输入受限的航空发动机最大推力状态鲁棒容错控制器

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010543149.XA CN111859555A (zh) 2020-06-15 2020-06-15 输入受限的航空发动机最大推力状态鲁棒容错控制器

Publications (1)

Publication Number Publication Date
CN111859555A true CN111859555A (zh) 2020-10-30

Family

ID=72986649

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010543149.XA Pending CN111859555A (zh) 2020-06-15 2020-06-15 输入受限的航空发动机最大推力状态鲁棒容错控制器

Country Status (1)

Country Link
CN (1) CN111859555A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112464387A (zh) * 2021-01-26 2021-03-09 中国人民解放军国防科技大学 一种喉栓式变推力固体发动机推力匹配设计方法
CN112879165A (zh) * 2021-01-14 2021-06-01 西北工业大学 考虑气路部件故障的航空发动机加速过程最优控制方法
CN112943453A (zh) * 2021-01-21 2021-06-11 西北工业大学 气路部件故障下基于iga的发动机最大推力控制优化方法
CN112947064A (zh) * 2021-01-21 2021-06-11 西北工业大学 考虑气路部件故障的航空发动机最大推力控制优化方法
CN112949161A (zh) * 2021-01-21 2021-06-11 西北工业大学 气路部件故障下基于iga的发动机最低油耗控制优化方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000267705A (ja) * 1999-03-15 2000-09-29 Oce Technol Bv モジュラーエンジン用の時間スケジュール方法及びスケジュラ
US20110288836A1 (en) * 2008-11-28 2011-11-24 Snecma Detection of anomalies in an aircraft engine
CN107942653A (zh) * 2017-10-30 2018-04-20 南京航空航天大学 航空电动燃油泵流量控制系统传感器故障鲁棒容错方法
CN110377043A (zh) * 2019-07-13 2019-10-25 西北工业大学 一种基于h∞回路成形算法的小型固定翼无人机姿态控制方法
CN110513198A (zh) * 2019-08-13 2019-11-29 大连理工大学 一种涡扇发动机控制系统主动容错控制方法
CN111273554A (zh) * 2020-04-04 2020-06-12 西北工业大学 航空发动机最大推力状态降保守性二自由度h∞控制器
CN111271181A (zh) * 2020-04-04 2020-06-12 西北工业大学 航空发动机降保守性增益调度二自由度μ控制器

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000267705A (ja) * 1999-03-15 2000-09-29 Oce Technol Bv モジュラーエンジン用の時間スケジュール方法及びスケジュラ
US20110288836A1 (en) * 2008-11-28 2011-11-24 Snecma Detection of anomalies in an aircraft engine
CN107942653A (zh) * 2017-10-30 2018-04-20 南京航空航天大学 航空电动燃油泵流量控制系统传感器故障鲁棒容错方法
CN110377043A (zh) * 2019-07-13 2019-10-25 西北工业大学 一种基于h∞回路成形算法的小型固定翼无人机姿态控制方法
CN110513198A (zh) * 2019-08-13 2019-11-29 大连理工大学 一种涡扇发动机控制系统主动容错控制方法
CN111273554A (zh) * 2020-04-04 2020-06-12 西北工业大学 航空发动机最大推力状态降保守性二自由度h∞控制器
CN111271181A (zh) * 2020-04-04 2020-06-12 西北工业大学 航空发动机降保守性增益调度二自由度μ控制器

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
LINFENG GOU .ETAL: "Aeroengine Robust Gain-Scheduling Control Based on Performance Degradation", 《IEEE ACCESS》 *
MEHRDAD PAKMEHR .ETAL: "Gain Scheduled Control of Gas Turbine Engines Stability and Verification", 《JOURNAL OF ENGINEERING FOR GAS TURBINES AND POWER》 *
贾秋生 等: "基于混合区域极点配置的航空发动机全包线鲁棒变参数控制器设计", 《推进技术》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112879165A (zh) * 2021-01-14 2021-06-01 西北工业大学 考虑气路部件故障的航空发动机加速过程最优控制方法
CN112943453A (zh) * 2021-01-21 2021-06-11 西北工业大学 气路部件故障下基于iga的发动机最大推力控制优化方法
CN112947064A (zh) * 2021-01-21 2021-06-11 西北工业大学 考虑气路部件故障的航空发动机最大推力控制优化方法
CN112949161A (zh) * 2021-01-21 2021-06-11 西北工业大学 气路部件故障下基于iga的发动机最低油耗控制优化方法
CN112464387A (zh) * 2021-01-26 2021-03-09 中国人民解放军国防科技大学 一种喉栓式变推力固体发动机推力匹配设计方法

Similar Documents

Publication Publication Date Title
CN111859555A (zh) 输入受限的航空发动机最大推力状态鲁棒容错控制器
CN111271181B (zh) 航空发动机降保守性增益调度二自由度μ控制器
CN111608808A (zh) 输入受限航空发动机增益调度容错控制器
CN111273554B (zh) 航空发动机最大推力状态降保守性二自由度h∞控制器
CN111880403A (zh) 航空发动机最大推力状态容错二自由度μ控制器
CN108829928B (zh) 一种涡轴发动机自适应部件级仿真模型构建方法
CN111856919A (zh) 航空发动机气路部件故障增益调度容错控制器
CN110502840B (zh) 航空发动机气路参数在线预测方法
CN112729857B (zh) 航空发动机健康参数估计方法及航空发动机自适应模型
CN108062428B (zh) 一种涡扇发动机在线部件故障诊断方法及系统
WO2014004494A1 (en) Real time linearization of a component-level gas turbine engine model for model-based control
CN112284752A (zh) 一种基于改进状态跟踪滤波器的变循环发动机解析余度估计方法
CN111305954B (zh) 输入受限的航空发动机降保守性鲁棒增益调度控制器
CN111856929B (zh) 航空发动机容错增益调度二自由度h∞控制器
CN112947064A (zh) 考虑气路部件故障的航空发动机最大推力控制优化方法
CN111830827B (zh) 航空发动机容错增益调度二自由度μ控制器
CN112327602A (zh) 变循环发动机气路部件故障增益调度容错控制器
CN112346336A (zh) 航空发动机气路部件故障鲁棒增益调度容错控制器
CN111852662A (zh) 航空发动机最大推力状态容错二自由度h∞控制器
CN111852663A (zh) 变循环发动机降保守性鲁棒增益调度控制器
CN112360634A (zh) 航空发动机最大推力状态鲁棒容错控制器
CN112377311A (zh) 输入受限的航空发动机鲁棒增益调度容错控制器
CN111456856A (zh) 航空发动机最大推力状态降保守性鲁棒控制器
CN111456857B (zh) 航空发动机降保守性增益调度二自由度h∞控制器
CN112879165A (zh) 考虑气路部件故障的航空发动机加速过程最优控制方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20201030

WD01 Invention patent application deemed withdrawn after publication