CN111608808A - 输入受限航空发动机增益调度容错控制器 - Google Patents

输入受限航空发动机增益调度容错控制器 Download PDF

Info

Publication number
CN111608808A
CN111608808A CN202010544561.3A CN202010544561A CN111608808A CN 111608808 A CN111608808 A CN 111608808A CN 202010544561 A CN202010544561 A CN 202010544561A CN 111608808 A CN111608808 A CN 111608808A
Authority
CN
China
Prior art keywords
engine
fault
input
vector
tolerant
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
CN202010544561.3A
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 CN202010544561.3A priority Critical patent/CN111608808A/zh
Publication of CN111608808A publication Critical patent/CN111608808A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02CGAS-TURBINE PLANTS; AIR INTAKES FOR JET-PROPULSION PLANTS; CONTROLLING FUEL SUPPLY IN AIR-BREATHING JET-PROPULSION PLANTS
    • F02C9/00Controlling gas-turbine plants; Controlling fuel supply in air- breathing jet-propulsion plants
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05DINDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
    • F05D2260/00Function
    • F05D2260/80Diagnostics
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05DINDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
    • F05D2260/00Function
    • F05D2260/81Modelling or simulation
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05DINDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
    • F05D2270/00Control
    • F05D2270/70Type of control algorithm
    • F05D2270/706Type of control algorithm proportional-integral-differential
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05DINDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
    • F05D2270/00Control
    • F05D2270/70Type of control algorithm
    • F05D2270/71Type of control algorithm synthesized, i.e. parameter computed by a mathematical model

Landscapes

  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明提出一种输入受限航空发动机增益调度容错控制器,包括PID控制器组容错控制模块、输入限制模块和气路部件故障诊断模块;PID控制器组容错控制模块产生控制向量v并输出给输入限制模块,输入限制模块产生限制后的控制输入向量u并输出给航空发动机本体,传感器得到发动机测量参数y;控制输入向量u以及测量参数y共同输入到气路部件故障诊断模块,气路部件故障诊断模块诊断发动机的气路部件故障,解算得到发动机的健康参数h,并输出到PID控制器组容错控制模块;PID控制器组容错控制模块内设计有若干PID控制器,PID控制器组容错控制模块根据输入的健康参数h以及调度参数α,利用内部设计的若干PID控制器计算得到适应的PID控制器,该PID控制器根据参考输入r和测量参数y的差值e产生控制输入向量u。本发明能够在保证发动机安全工作的前提下,在发动机机气路部件故障的情况下依旧对真实发动机进行良好控制,保证发动机安全工作,充分发挥发动机的性能,提高飞机的安全性和性能。

Description

输入受限航空发动机增益调度容错控制器
技术领域
本发明涉及航空发动机控制技术领域,尤其涉及一种输入受限航空发动机增益调度容错控制器。
背景技术
航空发动机是一个复杂的非线性动力学***,在范围宽广的飞行包线内工作时,发动机的工作状态随着外部条件和飞行条件的变化而不断变化。采用单一的线性控制器无法对发动机在全飞行包线内进行良好的控制,所以有人提出了非线性控制器设计方法,然而该方法不成熟而且非常复杂。针对发动机的控制,使用较多的是增益调度控制,通过对非线性模型中的一些稳定设计点进行线性化以获得对应各个点的线性化模型,然后针对各线化模型分别设计控制器,运用拟合或插值的方法将各个控制器连接起来,从而使用线性控制器对非线性发动机进行有效控制。
然而,现代战机对航空发动机性能的要求不断提高,其结构也越来越复杂,并且由于发动机工作环境的恶劣多变,发动机故障约占飞机总故障的1/3。其中,气路部件故障占发动机总体故障的90%以上,其维护费用占发动机总体维护费用的60%。为了保证发动机安全工作并使故障发动机提供足够的性能来保证飞机安全飞行或具有高的机动性,必须对故障的发动机性能进行恢复,并且对发动机进行容错控制,保证控制***正常稳定工作且具有良好的性能。因此,研究发动机气路部件故障容错控制方法具有重要意义。
传统的气路部件故障容错控制方法在航空发动机出现气路部件故障时通过修正控制规律,使得发动机的推力与油门杆始终匹配,有效的保证了发动机的推力。然而,这些设计方法并没有解决当前控制器和发动机模型不匹配从而导致控制***性能下降甚至不稳定的问题。当发动机发生气路部件故障时,发动机在同一工作点的线性化模型也会发生较大变化。因此,根据正常状态的发动机模型设计的控制器一般无法保证气路部件故障时发动机的性能,甚至无法保证控制***的闭环稳定。
此外,过大的控制输入会导致发动机损坏,因此我们需要考虑控制输入受限的控制器的设计。
发明内容
为解决现有技术存在的问题,本发明提出一种输入受限航空发动机增益调度容错控制器,能够在发动机机气路部件故障的情况下依旧对真实发动机进行良好控制,保证发动机安全工作,充分发挥发动机的性能,提高飞机的安全性和性能。并且考虑控制输入受限,保证发动机安全工作。
本发明的技术方案为:
所述一种输入受限航空发动机增益调度容错控制器,其特征在于:包括PID控制器组容错控制模块、输入限制模块和气路部件故障诊断模块;
其中PID控制器组容错控制模块、输入限制模块、气路部件故障诊断模块与航空发动机本体以及航空发动机上的若干传感器组成气路部件故障调度控制回路;
所述PID控制器组容错控制模块产生控制向量v并输出给输入限制模块,输入限制模块产生限制后的控制输入向量u并输出给航空发动机本体,传感器得到航空发动机测量参数y;控制输入向量u以及测量参数y共同输入到气路部件故障诊断模块,气路部件故障诊断模块诊断发动机的气路部件故障情况得到航空发动机的健康参数h,并输出到PID控制器组容错控制模块;
PID控制器组容错控制模块、输入限制模块、与航空发动机本体以及航空发动机上的若干传感器还组成调度参数调度控制回路;由传感器输出调度参数α至PID控制器组容错控制模块;
所述输入限制模块限制了控制输入向量的幅值,避免给发动机过大的控制输入导致发动机损坏;
所述PID控制器组容错控制模块内设计有若干PID控制器,所述PID控制器是利用若干线性化模型而分别设计得到的,所述线性化模型是对航空发动机不同设定工作点下的、不同气路部件故障下的航空发动机非线性模型进行线性化后得到的;
所述PID控制器组容错控制模块根据输入的健康参数h以及调度参数α,利用内部设计的若干PID控制器计算得到适应的PID控制器,该PID控制器根据参考输入r和测量参数y的差值e产生控制向量v。
进一步的,所述PID控制器组容错控制模块内设计若干PID控制器的过程为:在全飞行包线内根据调度参数α选取q个工作点对包含健康参数的发动机非线性模型进行线性化得到q个含有健康参数的线性化模型,通过调整健康参数的值,得到分别在发动机无气路部件故障和特定气路部件故障处的11q个线性化模型,并对这11q个线性化模型分别设计相应的PID控制器从而组成PID控制器组。
进一步的,所述气路部件故障诊断模块中包括非线性机载发动机模型和分段线性化卡尔曼滤波器;
所述非线性机载发动机模型为带健康参数的发动机非线性模型:
Figure BDA0002540106730000031
y=g(x,u,h)
其中
Figure BDA0002540106730000032
为控制输入向量,
Figure BDA0002540106730000033
为状态向量,
Figure BDA0002540106730000034
为输出向量,
Figure BDA0002540106730000035
为健康参数向量,f(·)为表示***动态的n维可微非线性向量函数,g(·)为产生***输出的m维可微非线性向量函数;非线性机载发动机模型输入为控制输入向量u以及上一周期的健康参数h,其输出的健康稳态参考值(xaug,NOBEM,yNOBEM)作为分段线性化卡尔曼滤波器当前周期的估计初始值;
所述分段线性化卡尔曼滤波器的输入为测量参数y以及非线性机载发动机模型输出的健康稳态参考值(xaug,NOBEM,yNOBEM),根据公式
Figure BDA0002540106730000036
计算得到当前周期的发动机的健康参数h;其中
Figure BDA0002540106730000037
K为卡尔曼滤波的增益,满足
Figure BDA0002540106730000038
P为Ricati方程
Figure BDA0002540106730000039
的解;系数Aaug和Caug根据公式
Figure BDA00025401067300000310
确定,而A、C、L、M是将健康参数h看作发动机的控制输入,并对非线性机载发动机模型在健康稳态参考点处进行线性化得到的反映发动机性能退化的增广线性状态变量模型
Figure BDA0002540106730000041
的系数:
Figure BDA0002540106730000042
Figure BDA0002540106730000043
w为***噪声,v为测量噪声,相应的协方差矩阵为对角阵Q和R。
进一步的,所述PID控制器组容错控制模块根据输入的健康参数h以及调度参数α插值得到的适应的PID控制器。
进一步的,所述PID控制器组容错控制模块先根据航空发动机当前的调度参数α选择前后相邻的两个设定工作点αi和αi+1,并获取两个设定工作点αi和αi+1对应发动机无部件故障的控制器Ki和Ki+1,各种典型部件故障Δhbase_j的控制器
Figure BDA0002540106730000044
Δhbase_j表示向量Δh的第j个元素的值为Δhbase,其他元素的值为0,即Δhbase_j表示10种不同的部件故障,例如Δhbase_1表示风扇发生了故障且风扇的效率变化量为Δhbase。根据公式
Figure BDA0002540106730000045
Figure BDA0002540106730000046
计算得到航空发动机所选工作点αi和αi+1处发动机的当前部件故障程度(健康参数为h)下的PID控制器Ki和Ki+1(式中Δhj为向量Δh的第j个元素;仅考虑||Δh||≤||Δhmax||的发动机气路部件故障情况,当||Δh||>||Δhmax||发动机已失效);再根据公式
Figure BDA0002540106730000047
计算得到航空发动机当前适应的容错PID控制器K(α)。
进一步的,所述输入限制模块采用多维矩形饱和函数,控制输入向量u为:
Figure BDA0002540106730000051
Figure BDA0002540106730000052
其中v1和vm为控制向量v的元素,v1,max和vm,max为控制向量v对应元素的限幅值。
进一步的,所述调度参数α包括航空发动机的风扇转速或者压气机转速。
进一步的,所述测量参数包括进气道出口、风扇出口、压气机出口、高压涡轮后、低压涡轮后的温度和压力,风扇转速和压气机转速。
有益效果
与现有技术相比较,本发明的输入受限航空发动机增益调度容错控制器利用传统增益调度控制器中固有的模块,通过新增气路部件故障诊断模块、输入限制模块,并对PID控制器组容错控制模块进行了改进,新增了发动机不同气路部件故障下的多组PID控制器。气路部件故障诊断模块通过健康参数的可靠估计实现了气路部件故障的准确判断,进而结合传统的调度参数,在保证发动机安全工作的前提下,实现发动机气路部件故障时的增益调度控制,保证发动机发生气路部件故障时仍安全工作,最大限度的提高发动机气路部件故障时增益调度的控制精度,缩短控制***的过渡时间,减小其动态偏差及静态偏差。经控制器控制非线性受控***,以使***在整个工作范围内,获得理想的动静态控制品质。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1是本发明输入受限航空发动机增益调度容错控制器的结构简图;
图2是本实施例气路部件故障调度控制回路中气路部件故障诊断模块的结构示意图;
图3是本实施例气路部件故障诊断模块中卡尔曼滤波器的结构示意图;
图4是本发明非线性发动机模型示意图。
具体实施方式
航空发动机在运行过程中由于自然磨损、腐蚀、积垢以及热蠕变等因素都会引起气路部件性能发生蜕化,并且当性能蜕化到一定程度会引发故障;另外,由于外物吸入引起的损伤、机械疲劳断裂等原因,也会引起气路部件故障的发生。前者故障发生的过程较为缓慢,而后者故障发生的过程是迅速的。当发动机气路部件发生故障却并未失效时,此时发动机的部分性能将会严重偏离额定状态。以涡轮部件为例,当其发生故障后,其工作效率将会下降,即将具有高温、高压的燃气转换为机械能的能力将会有所降低,但仍能为风扇或压气机部件提供相应的动力,使其工作在新的平衡状态下。此时发动机也已较大偏离原始状态。气路部件故障会导致发动机设计时所建立的非线性模型和气路部件故障时的真实发动机严重不匹配,进而导致根据该非线性模型设计的增益调度控制器无法对气路部件故障后的发动机进行良好的控制,严重降低了发动机的性能,甚至不能保证控制***的稳定性,无法保证发动机安全工作,并且过大的控制输入会导致发动机损坏。针对这一问题,下面给出本发明的分析研究过程。
1、发动机气路部件故障诊断
气路部件故障会导致部件对应的特性参数发生变化。发动机气路部件故障最终表征在不同转子部件的工作效率和流通量的改变上,即可以从风风扇、压气机、主燃烧、高压涡轮和低压涡轮部件的效率系数或者流量系数的变化来揭示发动机故障位置以及故障程度,风扇、压气机、主燃烧室、高压涡轮和低压涡轮部件的效率系数或流量系数被称为健康参数。
基于部件法,建立带健康参数的发动机非线性模型
Figure BDA0002540106730000071
y=g(x,u,h)
其中
Figure BDA0002540106730000072
为控制输入向量,
Figure BDA0002540106730000073
为状态向量,
Figure BDA0002540106730000074
为输出向量,
Figure BDA0002540106730000075
为健康参数向量,f(·)为表示***动态的n维可微非线性向量函数,g(·)为产生***输出的m维可微非线性向量函数。
将健康参数h看作发动机的控制输入,采用小扰动法或拟合法对发动机非线性模型在健康稳态参考点处进行线性化。
Figure BDA0002540106730000076
其中
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 BDA0002540106730000077
其中系数矩阵可由下式得到:
Figure BDA0002540106730000078
Figure BDA0002540106730000079
这些系数在发动机不同的工作状态具有不同的值。
实际上,健康参数很难测量,甚至不可能测量,而发动机各部分的压力、温度、转速等参数比较容易通过测量得到,通常称为“测量参数”,主要包括进气道出口、风扇出口、压气机出口、高压涡轮后、低压涡轮后的温度和压力,风扇转速和压气机转速。当发动机工作环境不发生变化时,健康参数的变化会引起被测参数的相应变化,二者之间存在气动热力学关系。因此,可以设计最优估计滤波器,通过测量参数来实现健康参数的最优估计。
对于渐变型部件故障,对应故障部件健康参数变化缓慢,因此在进行单次故障诊断的时间周期内,可以认为满足
Figure BDA0002540106730000081
对于突变型部件故障,更加关心的是发动机在故障发生后再次稳定工作时部件故障的严重程度,发动机再次稳定工作后,故障部件的健康参数变化依旧满足
Figure BDA0002540106730000082
将健康参数进一步转化为状态变量,可以得到
Figure BDA0002540106730000083
其中
Figure BDA0002540106730000084
Figure BDA0002540106730000085
建立的气路部件故障诊断模块主要由两部分组成,一部分是基于健康参数的非线性机载发动机模型,另一部分是分段线性卡尔曼滤波器。基本工作原理是将非线性机载发动机模型的输出作为分段线性卡尔曼滤波器的稳态参考值,并扩展健康参数,通过分段线性卡尔曼滤波器进行在线实时估计,最后反馈给非线性机载发动机模型进行在线实时更新,实现对实际发动机的实时跟踪。
卡尔曼估计方程为:
Figure BDA0002540106730000086
K为卡尔曼滤波的增益,满足
Figure BDA0002540106730000087
P为Ricati方程
Figure BDA0002540106730000088
的解;利用非线性机载模型输出的健康稳态参考值(xaug,NOBEM,yNOBEM)作为式
Figure BDA0002540106730000091
的初值,可得计算公式:
Figure BDA0002540106730000092
根据该计算公式可以得到发动机的健康参数h,实现发动机的气路部件故障诊断。
2、增益调度容错控制设计
增益调度控制的实质是设计一组线性化的控制器,然后将它们有规律地组合起来,从而能够控制非线性***。增益调度容错控制的基本原理是选择一系列的工作点,获得不同设定工作点下的、不同气路部件故障下的发动机线性化模型并分别设计对应的PID控制器得到图1中的PID控制器组。
请参阅图4,选择一组调度参数值αi,i=1,2,...,q,代表***的动态范围,并将飞行包线划分为几个子区间,并将这些点作为工作点。在工作点,有这些方程
Figure BDA0002540106730000093
Figure BDA0002540106730000094
其中
Figure BDA0002540106730000095
为所选择的第i个工作点,udi为在时刻
Figure BDA0002540106730000096
保持平衡所需的稳态控制输入,hdi为时刻
Figure BDA0002540106730000097
的健康参数。
利用小扰动法,可以得到各工况点健康参数的线性模型,并得到发动机正常状态和性能退化h状态下的线性模型。
请参阅图4,上、下实线分别表示发动机无气路部件故障和有气路部件故障h的非线性模型。一系列的黑色小圆点表示发动机不同的工作点,在每一个工作点进行线性化得到线性模型。针对发动机正常状态和不同气路部件故障状态下的线性模型,分别设计一系列PID控制器得到图1中的PID控制器组。然后,在选定的工作点之间,对控制器增益进行线性内插,使得对于所有的固定参数值,闭环***都稳定且具有良好的性能。参数α是调度参数,这里可以定义为航空发动机的风扇转速或者压气机转速,可以实时测量。控制***的另一个调度变量是反映发动机气路部件故障程度的健康参数h。工作原理是图1中的PID控制器组容错控制模块根据调度参数和健康参数进行线性插值获得相应的PID控制器来控制***。
3、控制器的插值
这部分说明了图1中的PID控制器组容错控制模块通过调度参数和健康参数调度线性插值获得相应的PID控制器的调度计算原理。
分别在发动机正常状态和各种典型部件故障Δhbase_j状态下设计一系列线性PID控制器,对每个选定的工作点αi进行控制。这将产生图1中的PID控制器组容错控制模块中的控制器
Ki
Figure BDA0002540106730000101
然后根据调度参数α和健康参数h对控制器进行插值,继而使用得到的插值控制器来控制***。
根据发动机当前的调度参数α选定周围相邻的两个工作点αi和αi+1,并获取两个设定工作点αi和αi+1对应发动机无部件故障的控制器Ki和Ki+1,各种典型部件故障Δhbase_j的控制器
Figure BDA0002540106730000102
Δhbase_j表示向量Δh的第j个元素的值为Δhbase,其他元素的值为0,即Δhbase_j表示10种不同的部件故障,例如Δhbase_1表示风扇发生了故障且风扇的效率变化量为Δhbase。通过线性插值可以得到在工作点αi处气路部件故障h处的控制器
Figure BDA0002540106730000103
同样的,可以得到在工作点αi+1处气路部件故障h处的控制器
Figure BDA0002540106730000104
我们使用分段线性插值的方法,从PID控制器集K1,K2,...,Kq中对每一对控制器之间进行线性插值。得到当前调度参数α当前退化程度h处的线性插值控制器K(α),i=1,2,...,q-1为
Figure BDA0002540106730000111
根据该公式可以得到某一工作点某一气路部件故障下相应的控制器,并对发动机进行有效控制。
4.***的输入限制
请参考图1,图1中的输入限制模块是为了建模***控制输入上的物理限制,使用了多维矩形饱和函数。限制航空发动机控制的输入,尤其是对于燃油流量输入。多维饱和函数还可以处理其他控制输入的限制,包括尾喷嘴的喉部面积。该函数是一个多维矩形饱和函数,定义为
Figure BDA0002540106730000112
其中v1和vm为控制向量v的元素,v1,max和vm,max为控制向量v对应元素的限幅值。对于所有的
Figure BDA0002540106730000113
下式给出了sat(·)
Figure BDA0002540106730000114
基于上述过程,下面给出本实施例中提出的一种输入受限航空发动机增益调度容错控制器,如图1所示,主要包括PID控制器组容错控制模块、输入限制模块和气路部件故障诊断模块。
其中PID控制器组容错控制模块、输入限制模块、气路部件故障诊断模块与航空发动机本体以及航空发动机上的若干传感器组成气路部件故障调度控制回路10。
所述PID控制器组容错控制模块产生控制向量v并输出给输入限制模块,输入限制模块产生限制后的控制输入向量u并输出给航空发动机本体,传感器得到航空发动机测量参数y;控制输入向量u以及测量参数y共同输入到气路部件故障诊断模块,气路部件故障诊断模块解算得到航空发动机的健康参数h,并输出到PID控制器组容错控制模块。
PID控制器组容错控制模块、输入限制模块、与航空发动机本体以及航空发动机上的若干传感器还组成调度参数调度控制回路20;由传感器输出调度参数α至PID控制器组容错控制模块。
所述输入限制模块限制了控制输入向量的幅值,避免给发动机过大的控制输入导致发动机损坏。
优选的一种具体实现方式,所述输入限制模块采用多维矩形饱和函数,控制输入向量u为:
Figure BDA0002540106730000121
Figure BDA0002540106730000122
其中v1和vm为控制向量v的元素,v1,max和vm,max为控制向量v对应元素的限幅值。
所述PID控制器组容错控制模块内设计有若干PID控制器,所述PID控制器是利用若干线性化模型而分别设计得到的,所述线性化模型是对航空发动机不同设定工作点下的、不同气路部件故障下的航空发动机非线性模型进行线性化后得到的。
优选的一种具体实现方式,可以通过以下过程得到设计若干PID控制器:在全飞行包线内根据调度参数α选取q个工作点对包含健康参数的发动机非线性模型进行线性化得到q个含有健康参数的线性化模型,通过调整健康参数的值,得到分别在发动机无气路部件故障和特定气路部件故障处的11q个线性化模型,并对这11q个线性化模型分别设计相应的PID控制器从而组成PID控制器组。
所述PID控制器组容错控制模块根据输入的健康参数h以及调度参数α,利用内部设计的若干PID控制器计算得到适应的PID控制器,该PID控制器根据参考输入r和测量参数y的差值e产生控制输入向量u。
优选的一种具体实现方式,可以根据输入的健康参数h以及调度参数α插值得到的适应的PID控制器:
先根据航空发动机当前的调度参数α选择前后相邻的两个设定工作点αi和αi+1,并获取两个设定工作点αi和αi+1对应发动机无部件故障的控制器Ki和Ki+1,各种典型部件故障Δhbase_j的控制器
Figure BDA0002540106730000131
Δhbase_j表示向量Δh的第j个元素的值为Δhbase,其他元素的值为0,即Δhbase_j表示10种不同的部件故障,例如Δhbase_1表示风扇发生了故障且风扇的效率变化量为Δhbase。根据公式
Figure BDA0002540106730000132
Figure BDA0002540106730000133
计算得到航空发动机所选工作点αi和αi+1处发动机的当前部件故障程度(健康参数为h)下的PID控制器Ki和Ki+1(式中Δhj为向量Δh的第j个元素;仅考虑||Δh||≤||Δhmax||的发动机气路部件故障情况,当||Δh||>||Δhmax||发动机已失效);再根据公式
Figure BDA0002540106730000134
计算得到航空发动机当前适应的容错PID控制器K(α)。
所述气路部件故障诊断模块中包括非线性机载发动机模型和分段线性化卡尔曼滤波器。
所述非线性机载发动机模型为带健康参数的发动机非线性模型:
Figure BDA0002540106730000135
y=g(x,u,h)
其中
Figure BDA0002540106730000136
为控制输入向量,
Figure BDA0002540106730000137
为状态向量,
Figure BDA0002540106730000138
为输出向量,
Figure BDA0002540106730000139
为健康参数向量,f(·)为表示***动态的n维可微非线性向量函数,g(·)为产生***输出的m维可微非线性向量函数;非线性机载发动机模型输入为控制输入向量u以及上一周期的健康参数h,其输出的健康稳态参考值(xaug,NOBEM,yNOBEM)作为分段线性化卡尔曼滤波器当前周期的估计初始值。
所述分段线性化卡尔曼滤波器的输入为测量参数y以及非线性机载发动机模型输出的健康稳态参考值(xaug,NOBEM,yNOBEM),根据公式
Figure BDA0002540106730000141
计算得到当前周期的发动机的健康参数h。
其中
Figure BDA0002540106730000142
K为卡尔曼滤波的增益,满足
Figure BDA0002540106730000143
P为Ricati方程
Figure BDA0002540106730000144
的解;系数Aaug和Caug根据公式
Figure BDA0002540106730000145
确定,而A、C、L、M是将健康参数h看作发动机的控制输入,并对非线性机载发动机模型在健康稳态参考点处进行线性化得到的反映发动机性能退化的增广线性状态变量模型
Figure BDA0002540106730000146
的系数:
Figure BDA0002540106730000147
Figure BDA0002540106730000148
w为***噪声,v为测量噪声,相应的协方差矩阵为对角阵Q和R。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (7)

1.一种输入受限航空发动机增益调度容错控制器,其特征在于:包括PID控制器组容错控制模块、输入限制模块和气路部件故障诊断模块;
其中PID控制器组容错控制模块、输入限制模块、气路部件故障诊断模块与航空发动机本体以及航空发动机上的若干传感器组成气路部件故障调度控制回路;
所述PID控制器组容错控制模块产生控制向量v并输出给输入限制模块,输入限制模块产生限制后的控制输入向量u并输出给航空发动机本体,传感器得到航空发动机测量参数y;控制输入向量u以及测量参数y共同输入到气路部件故障诊断模块,气路部件故障诊断模块诊断发动机的气路部件故障情况得到航空发动机的健康参数h,并输出到PID控制器组容错控制模块;
PID控制器组容错控制模块、输入限制模块与航空发动机本体以及航空发动机上的若干传感器还组成调度参数调度控制回路;由传感器输出调度参数α至PID控制器组容错控制模块;
所述输入限制模块限制了控制输入向量的幅值,避免给发动机过大的控制输入导致发动机损坏;
所述PID控制器组容错控制模块内设计有若干PID控制器,所述PID控制器是利用若干线性化模型而分别设计得到的,所述线性化模型是对航空发动机不同设定工作点下的、不同气路部件故障下的航空发动机非线性模型进行线性化后得到的;
所述PID控制器组容错控制模块根据输入的健康参数h以及调度参数α,利用内部设计的若干PID控制器计算得到适应的PID控制器,该PID控制器根据参考输入r和测量参数y的差值e产生控制输入向量u。
2.根据权利要求1所述一种输入受限航空发动机增益调度容错控制器,其特征在于:所述PID控制器组容错控制模块内设计若干PID控制器的过程为:在全飞行包线内根据调度参数α选取q个工作点对包含健康参数的发动机非线性模型进行线性化得到q个含有健康参数的线性化模型,通过调整健康参数的值,得到分别在发动机无气路部件故障和特定气路部件故障处的11q个线性化模型,并对这11q个线性化模型分别设计相应的PID控制器从而组成PID控制器组。
3.根据权利要求1或2所述一种输入受限航空发动机增益调度容错控制器,其特征在于:所述PID控制器组容错控制模块根据输入的健康参数h以及调度参数α插值得到适应的PID控制器。
根据权利要求3所述一种输入受限航空发动机增益调度容错控制器,其特征在于:所述PID控制器组容错控制模块先根据航空发动机当前的调度参数α选择前后相邻的两个设定工作点αi和αi+1,并获取两个设定工作点αi和αi+1对应发动机无部件故障的控制器Ki和Ki+1,各种典型部件故障Δhbase_j的控制器
Figure FDA0002540106720000021
Δhbase_j表示向量Δh的第j个元素的值为Δhbase,其他元素的值为0,即Δhbase_j表示10种不同的部件故障,例如Δhbase_1表示风扇发生了故障且风扇的效率变化量为Δhbase。根据公式
Figure FDA0002540106720000022
Figure FDA0002540106720000023
计算得到航空发动机所选工作点αi和αi+1处发动机的当前部件故障程度(健康参数为h)下的PID控制器Ki和Ki+1(式中Δhj为向量Δh的第j个元素;仅考虑||Δh||≤||Δhmax||的发动机气路部件故障情况,当||Δh||>||Δhmax||发动机已失效);再根据公式
Figure FDA0002540106720000024
计算得到航空发动机当前适应的容错PID控制器K(α)。
4.根据权利要求1所述一种输入受限航空发动机增益调度容错控制器,其特征在于:所述气路部件故障诊断模块中包括非线性机载发动机模型和分段线性化卡尔曼滤波器;
所述非线性机载发动机模型为带健康参数的发动机非线性模型:
Figure FDA0002540106720000025
y=g(x,u,h)
其中
Figure FDA0002540106720000031
为控制输入向量,
Figure FDA0002540106720000032
为状态向量,
Figure FDA0002540106720000033
为输出向量,
Figure FDA0002540106720000034
为健康参数向量,f(·)为表示***动态的n维可微非线性向量函数,g(·)为产生***输出的m维可微非线性向量函数;非线性机载发动机模型输入为控制输入向量u以及上一周期的健康参数h,其输出的健康稳态参考值(xaug,NOBEM,yNOBEM)作为分段线性化卡尔曼滤波器当前周期的估计初始值;
所述分段线性化卡尔曼滤波器的输入为测量参数y以及非线性机载发动机模型输出的健康稳态参考值(xaug,NOBEM,yNOBEM),根据公式
Figure FDA0002540106720000035
计算得到当前周期的发动机的健康参数h;其中
Figure FDA0002540106720000036
K为卡尔曼滤波的增益,满足
Figure FDA0002540106720000037
P为Ricati方程
Figure FDA0002540106720000038
的解;系数Aaug和Caug根据公式
Figure FDA0002540106720000039
Caug=(C M)
确定,而A、C、L、M是将健康参数h看作发动机的控制输入,并对非线性机载发动机模型在健康稳态参考点处进行线性化得到的反映发动机性能退化的增广线性状态变量模型
Figure FDA00025401067200000310
的系数:
Figure FDA00025401067200000311
Figure FDA00025401067200000312
w为***噪声,v为测量噪声,相应的协方差矩阵为对角阵Q和R。
5.根据权利要求1所述一种输入受限航空发动机增益调度容错控制器,其特征在于:
所述输入限制模块采用多维矩形饱和函数,控制输入向量u为:
Figure FDA0002540106720000041
Figure FDA0002540106720000042
其中v1和vm为控制向量v的元素,v1,max和vm,max为控制向量v对应元素的限幅值。
6.根据权利要求1所述一种输入受限航空发动机增益调度容错控制器,其特征在于:所述调度参数α包括航空发动机的风扇转速或者压气机转速。
7.根据权利要求1所述一种输入受限航空发动机增益调度容错控制器,其特征在于:所述测量参数包括进气道出口、风扇出口、压气机出口、高压涡轮后、低压涡轮后的温度和压力,风扇转速和压气机转速。
CN202010544561.3A 2020-06-15 2020-06-15 输入受限航空发动机增益调度容错控制器 Pending CN111608808A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010544561.3A CN111608808A (zh) 2020-06-15 2020-06-15 输入受限航空发动机增益调度容错控制器

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010544561.3A CN111608808A (zh) 2020-06-15 2020-06-15 输入受限航空发动机增益调度容错控制器

Publications (1)

Publication Number Publication Date
CN111608808A true CN111608808A (zh) 2020-09-01

Family

ID=72198907

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010544561.3A Pending CN111608808A (zh) 2020-06-15 2020-06-15 输入受限航空发动机增益调度容错控制器

Country Status (1)

Country Link
CN (1) CN111608808A (zh)

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 西北工业大学 考虑气路部件故障的航空发动机加速过程最优控制方法
CN112904881A (zh) * 2021-01-19 2021-06-04 杭州电子科技大学 一种高超声速飞行器动态增益调度控制器设计方法
CN112949160A (zh) * 2021-01-21 2021-06-11 西北工业大学 气路部件故障下基于ipso的发动机加速过程最优控制方法
CN112949161A (zh) * 2021-01-21 2021-06-11 西北工业大学 气路部件故障下基于iga的发动机最低油耗控制优化方法
CN112947064A (zh) * 2021-01-21 2021-06-11 西北工业大学 考虑气路部件故障的航空发动机最大推力控制优化方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109358491A (zh) * 2018-10-25 2019-02-19 中国电子科技集团公司第五十四研究所 基于卡尔曼滤波的模糊pid故障确定容错控制方法
CN110118128A (zh) * 2019-05-28 2019-08-13 南京航空航天大学 微型燃气轮机传感器故障诊断和容错控制方法
US20190383221A1 (en) * 2018-06-15 2019-12-19 Dalian University Of Technology Design method of aero-engine on-line optimization and multivariable control based on model prediction
CN111271181A (zh) * 2020-04-04 2020-06-12 西北工业大学 航空发动机降保守性增益调度二自由度μ控制器

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20190383221A1 (en) * 2018-06-15 2019-12-19 Dalian University Of Technology Design method of aero-engine on-line optimization and multivariable control based on model prediction
CN109358491A (zh) * 2018-10-25 2019-02-19 中国电子科技集团公司第五十四研究所 基于卡尔曼滤波的模糊pid故障确定容错控制方法
CN110118128A (zh) * 2019-05-28 2019-08-13 南京航空航天大学 微型燃气轮机传感器故障诊断和容错控制方法
CN111271181A (zh) * 2020-04-04 2020-06-12 西北工业大学 航空发动机降保守性增益调度二自由度μ控制器

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
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 西北工业大学 考虑气路部件故障的航空发动机加速过程最优控制方法
CN112904881A (zh) * 2021-01-19 2021-06-04 杭州电子科技大学 一种高超声速飞行器动态增益调度控制器设计方法
CN112949160A (zh) * 2021-01-21 2021-06-11 西北工业大学 气路部件故障下基于ipso的发动机加速过程最优控制方法
CN112949161A (zh) * 2021-01-21 2021-06-11 西北工业大学 气路部件故障下基于iga的发动机最低油耗控制优化方法
CN112947064A (zh) * 2021-01-21 2021-06-11 西北工业大学 考虑气路部件故障的航空发动机最大推力控制优化方法

Similar Documents

Publication Publication Date Title
CN111608808A (zh) 输入受限航空发动机增益调度容错控制器
CN111859555A (zh) 输入受限的航空发动机最大推力状态鲁棒容错控制器
CN111271181B (zh) 航空发动机降保守性增益调度二自由度μ控制器
CN111856919A (zh) 航空发动机气路部件故障增益调度容错控制器
US8315741B2 (en) High fidelity integrated heat transfer and clearance in component-level dynamic turbine system control
EP2154588B1 (en) Control system
US8849542B2 (en) Real time linearization of a component-level gas turbine engine model for model-based control
CN111880403A (zh) 航空发动机最大推力状态容错二自由度μ控制器
US20160069277A1 (en) Turboshaft engine control
CN110502840B (zh) 航空发动机气路参数在线预测方法
CN109557815B (zh) 一种航空发动机的传感器故障调节方法
CN112729857B (zh) 航空发动机健康参数估计方法及航空发动机自适应模型
CN112284752A (zh) 一种基于改进状态跟踪滤波器的变循环发动机解析余度估计方法
CN112327602A (zh) 变循环发动机气路部件故障增益调度容错控制器
CN111856929B (zh) 航空发动机容错增益调度二自由度h∞控制器
CN112947064A (zh) 考虑气路部件故障的航空发动机最大推力控制优化方法
CN111305954A (zh) 输入受限的航空发动机降保守性鲁棒增益调度控制器
CN111830827B (zh) 航空发动机容错增益调度二自由度μ控制器
CN112346336A (zh) 航空发动机气路部件故障鲁棒增益调度容错控制器
Smith et al. Optimizing aircraft performance with adaptive, integrated flight/propulsion control
CN111852663A (zh) 变循环发动机降保守性鲁棒增益调度控制器
CN111852662A (zh) 航空发动机最大推力状态容错二自由度h∞控制器
CN112377311A (zh) 输入受限的航空发动机鲁棒增益调度容错控制器
CN114047692B (zh) 一种涡扇发动机鲁棒容错抗干扰模型参考动态输出反馈控制方法
CN112948961A (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: 20200901

WD01 Invention patent application deemed withdrawn after publication