[go: up one dir, main page]

CN111456856A - Robust controller for reducing conservative maximum thrust state of aero-engine - Google Patents

Robust controller for reducing conservative maximum thrust state of aero-engine Download PDF

Info

Publication number
CN111456856A
CN111456856A CN202010261759.0A CN202010261759A CN111456856A CN 111456856 A CN111456856 A CN 111456856A CN 202010261759 A CN202010261759 A CN 202010261759A CN 111456856 A CN111456856 A CN 111456856A
Authority
CN
China
Prior art keywords
engine
maximum thrust
degradation
thrust state
state
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.)
Granted
Application number
CN202010261759.0A
Other languages
Chinese (zh)
Other versions
CN111456856B (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 CN202010261759.0A priority Critical patent/CN111456856B/en
Publication of CN111456856A publication Critical patent/CN111456856A/en
Application granted granted Critical
Publication of CN111456856B publication Critical patent/CN111456856B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Mechanical Engineering (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

The invention provides a conservative robust controller for reducing the maximum thrust state of an aircraft engine, which is characterized in that a gain scheduling controller group is improved by adding a degradation parameter estimation loop, a conservative robust controller for reducing the maximum thrust state under a certain degradation degree of the engine is added, and a resolving module for reducing the conservative robust controller group for reducing the maximum thrust state is obtained. The designed robust controller for reducing the conservative property of the maximum thrust state adopts a small perturbation uncertainty engine model, eliminates a degradation term in the uncertainty of the engine, reduces the perturbation range of the uncertainty model, and reduces the conservative property of the robust gain scheduling controller. The degradation parameter estimation loop realizes reliable estimation of degradation parameters, and gain scheduling control during engine performance degradation is realized by using the degradation parameters. The invention has strong robustness and low conservation, and improves the performance of the engine in the maximum thrust state to the maximum extent, so that the engine not only stably works in the maximum thrust state, but also improves the thrust in the maximum thrust state.

Description

航空发动机最大推力状态降保守性鲁棒控制器Conservative robust controller for maximum thrust state drop of aero-engine

技术领域technical field

本发明涉及航空发动机控制技术领域,尤其涉及一种航空发动机最大推力状态降保守性鲁棒控制器。The invention relates to the technical field of aero-engine control, in particular to an aero-engine maximum thrust state-down conservative robust controller.

背景技术Background technique

航空发动机是一个复杂的非线性动力学系统,其控制系统容易受到工作条件,发动机性能下降,环境条件变化的影响,并且很难事先知道外部干扰和测量噪声的影响。由于飞机发动机的工作过程非常复杂,难以建立准确的数学模型,所以数学模型与实际系统之间总是存在差异。因此,有必要设计一种鲁棒控制器,用于在外部干扰信号,噪声干扰,未建模的动态特性和参数变化的情况下稳定航空发动机控制系统,并具有良好的性能。Aeroengine is a complex nonlinear dynamic system, and its control system is easily affected by working conditions, engine performance degradation, changes in environmental conditions, and it is difficult to know in advance the influence of external disturbances and measurement noise. Because the working process of an aircraft engine is very complex, it is difficult to establish an accurate mathematical model, so there is always a difference between the mathematical model and the actual system. Therefore, it is necessary to design a robust controller for stabilizing the aero-engine control system with good performance in the presence of external disturbance signals, noise disturbances, unmodeled dynamic characteristics and parameter changes.

战斗机由于需要实现高机动性,发动机的最大推力状态的性能至关重要。传统的鲁棒控制器虽然可以对发动机在最大推力状态实现稳定控制,然而,它们是非常保守的,因为它们将发动机退化看作发动机模型的不确定性进行鲁棒控制器的设计。事实上,发动机的性能退化程度可以通过测量参数来估计,从而消除不确定性模型中的退化项,缩小不确定性模型的范围,降低鲁棒控制器的保守性,提高发动机在最大推力状态的性能,从而使得飞机具有更好的机动性,在战斗中具有更加明显的优势。Due to the need to achieve high maneuverability in fighter aircraft, the performance of the engine's maximum thrust state is critical. Although the traditional robust controllers can achieve stable control of the engine at the maximum thrust state, however, they are very conservative because they consider the engine degradation as the uncertainty of the engine model to design the robust controller. In fact, the degree of performance degradation of the engine can be estimated by measuring parameters, thereby eliminating the degradation term in the uncertainty model, reducing the range of the uncertainty model, reducing the conservatism of the robust controller, and improving the engine's performance in the maximum thrust state. performance, so that the aircraft has better maneuverability and has a more obvious advantage in combat.

发明内容SUMMARY OF THE INVENTION

为解决现有技术存在的问题,本发明提出一种航空发动机最大推力状态降保守性鲁棒控制器,具有强的鲁棒性并且保守性低,最大限度的提高发动机在最大推力状态的性能,使发动机在最大推力状态不仅稳定工作,并且提高发动机最大推力状态的性能,提高战斗机的机动性。In order to solve the problems existing in the prior art, the present invention proposes an aero-engine maximum thrust state reduction conservative robust controller, which has strong robustness and low conservativeness, and maximizes the performance of the engine in the maximum thrust state. Making the engine in the maximum thrust state not only works stably, but also improves the performance of the engine at the maximum thrust state and improves the maneuverability of the fighter.

本发明的技术方案为:The technical scheme of the present invention is:

所述一种航空发动机最大推力状态降保守性鲁棒控制器,其特征在于:包括最大推力状态降保守性鲁棒控制器组解算模块和退化参数估计回路;The aero-engine maximum thrust state reduction conservative robust controller is characterized in that: it includes a maximum thrust state reduction conservative robust controller group solution module and a degradation parameter estimation loop;

其中最大推力状态降保守性鲁棒控制器组解算模块、退化参数估计回路与航空发动机本体以及航空发动机上的若干传感器组成退化参数调度控制回路;Among them, the maximum thrust state reduction conservative robust controller group solution module, the degradation parameter estimation loop, the aero-engine body and several sensors on the aero-engine form the degradation parameter scheduling control loop;

所述最大推力状态降保守性鲁棒控制器组解算模块产生控制输入向量u并输出给航空发动机本体,传感器得到航空发动机测量参数y;控制输入向量u以及测量参数y共同输入到退化参数估计回路,退化参数估计回路解算得到航空发动机的退化参数h,并输出到最大推力状态降保守性鲁棒控制器组解算模块;The maximum thrust state degradation conservative robust controller group solution module generates a control input vector u and outputs it to the aero-engine body, and the sensor obtains the aero-engine measurement parameter y; the control input vector u and the measurement parameter y are jointly input to the degradation parameter estimation The loop, the degradation parameter estimation loop solves the aero-engine degradation parameter h, and outputs it to the maximum thrust state degradation conservative robust controller group solution module;

所述最大推力状态降保守性鲁棒控制器组解算模块内设计有两个最大推力状态降保守性鲁棒控制器,所述两个最大推力状态降保守性鲁棒控制器采用以下过程得到:分别在发动机正常状态h1和设定退化程度hbase处,在航空发动机最大推力状态下对包含退化参数的发动机非线性模型进行线性化得到2个线性化模型,对线性化模型加入不含发动机性能退化的摄动块得到小摄动不确定性发动机模型,并对这2个小摄动不确定性发动机模型分别设计鲁棒控制器,作为对应的最大推力状态降保守性鲁棒控制器;The maximum thrust state drop conservative robust controller group solution module is designed with two maximum thrust state drop conservative robust controllers, and the two maximum thrust state drop conservative robust controllers are obtained by using the following process : At the normal state of the engine h 1 and the set degradation degree h base , the nonlinear model of the engine including the degradation parameters is linearized under the maximum thrust state of the aero-engine to obtain two linearized models, and the linearized model is added without The engine model with small perturbation uncertainty is obtained from the perturbation block whose engine performance is degraded, and robust controllers are designed for these two engine models with small perturbation uncertainty as the corresponding maximum thrust state reduction conservative robust controller. ;

所述最大推力状态降保守性鲁棒控制器组解算模块根据输入的退化参数h,利用内部设计的两个最大推力状态降保守性鲁棒控制器计算得到适应的最大推力状态降保守性鲁棒控制器,该最大推力状态降保守性鲁棒控制器根据参考输入r和测量参数y的差值e产生控制输入向量u。The maximum thrust state reduction conservative robust controller group calculation module calculates the adapted maximum thrust state reduction conservative robustness by using two internally designed maximum thrust state reduction conservative robust controllers according to the input degradation parameter h. Rod controller, the maximum thrust state drop conservative robust controller generates a control input vector u according to the difference e between the reference input r and the measured parameter y.

进一步的,所述退化参数估计回路中包括非线性机载发动机模型和最大推力状态处卡尔曼滤波器;Further, the degradation parameter estimation loop includes a nonlinear airborne engine model and a Kalman filter at the maximum thrust state;

所述非线性机载发动机模型为带退化参数的发动机非线性模型:The nonlinear airborne engine model is an engine nonlinear model with degradation parameters:

Figure BDA0002439600450000021
Figure BDA0002439600450000021

y=g(x,u,h)y=g(x,u,h)

其中

Figure BDA0002439600450000022
为控制输入向量,
Figure BDA0002439600450000023
为状态向量,
Figure BDA0002439600450000024
为输出向量,
Figure BDA0002439600450000025
为退化参数向量,f(·)为表示系统动态的n维可微非线性向量函数,g(·)为产生系统输出的m维可微非线性向量函数;非线性机载发动机模型输入为控制输入向量u以及上一周期的退化参数h,其输出的健康稳态参考值(xaug,NOBEM,yNOBEM)作为最大推力状态处卡尔曼滤波器当前周期的估计初始值;in
Figure BDA0002439600450000022
For the control input vector,
Figure BDA0002439600450000023
is the state vector,
Figure BDA0002439600450000024
is the output vector,
Figure BDA0002439600450000025
is the degradation parameter vector, f(·) is the n-dimensional differentiable nonlinear vector function representing the system dynamics, g(·) is the m-dimensional differentiable nonlinear vector function that generates the system output; the input of the nonlinear airborne engine model is the control The input vector u and the degradation parameter h of the previous cycle, and the output healthy steady-state reference value (x aug, NOBEM , y NOBEM ) is used as the estimated initial value of the current cycle of the Kalman filter at the maximum thrust state;

所述最大推力状态处卡尔曼滤波器的输入为测量参数y以及非线性机载发动机模型输出的健康稳态参考值(xaug,NOBEM,yNOBEM),根据公式The input of the Kalman filter at the maximum thrust state is the measured parameter y and the healthy steady-state reference value (x aug, NOBEM , y NOBEM ) output by the nonlinear airborne engine model, according to the formula

Figure BDA0002439600450000031
Figure BDA0002439600450000031

计算得到当前周期的发动机的退化参数h;其中

Figure BDA0002439600450000032
K为卡尔曼滤波的增益,满足
Figure BDA0002439600450000033
P为Ricati方程
Figure BDA0002439600450000034
的解;系数Aaug和Caug根据公式Calculate the degradation parameter h of the engine in the current cycle; where
Figure BDA0002439600450000032
K is the gain of the Kalman filter, satisfying
Figure BDA0002439600450000033
P is the Ricati equation
Figure BDA0002439600450000034
The solution of ; the coefficients A aug and C aug according to the formula

Figure BDA0002439600450000035
Figure BDA0002439600450000035

确定,而A、C、L、M是将退化参数h看作发动机的控制输入,并对非线性机载发动机模型在健康稳态参考点处进行线性化得到的反映发动机性能退化的增广线性状态变量模型Determined, while A, C, L, and M are the augmented linearity reflecting the degradation of engine performance obtained by taking the degradation parameter h as the control input of the engine, and linearizing the nonlinear airborne engine model at the healthy steady-state reference point state variable model

Figure BDA0002439600450000036
的系数:
Figure BDA0002439600450000036
The coefficient of :

Figure BDA0002439600450000037
Figure BDA0002439600450000037

Figure BDA0002439600450000038
Figure BDA0002439600450000038

w为系统噪声,v为测量噪声,相应的协方差矩阵为对角阵Q和R。w is the system noise, v is the measurement noise, and the corresponding covariance matrices are the diagonal matrices Q and R.

进一步的,所述最大推力状态降保守性鲁棒控制器组解算模块根据航空发动机正常状态h1和设定退化程度hbase处的最大推力状态降保守性鲁棒控制器K、Kh_base,通过公式Further, the maximum thrust state degrades conservative robust controller group calculation module according to the normal state h 1 of the aero-engine and the maximum thrust state degrade conservative robust controllers K, K h_base at the set degradation degree h base , by formula

Figure BDA0002439600450000039
Figure BDA0002439600450000039

计算得到航空发动机当前退化状态适应的最大推力状态降保守性鲁棒控制器KhThe conservative robust controller K h is obtained by calculating the maximum thrust state adapted to the current degradation state of the aero-engine.

进一步的,所述测量参数包括进气道出口、风扇出口、压气机出口、高压涡轮后、低压涡轮后的温度和压力,风扇转速和压气机转速。Further, the measurement parameters include the temperature and pressure at the outlet of the intake duct, the outlet of the fan, the outlet of the compressor, after the high-pressure turbine and after the low-pressure turbine, the rotational speed of the fan and the rotational speed of the compressor.

有益效果beneficial effect

与现有技术相比较,本发明的航空发动机最大推力状态降保守性鲁棒控制器利用传统鲁棒控制器的设计方法,通过新增退化参数估计回路,并对增益调度控制器组进行了改进,新增了发动机一定退化程度下的最大推力状态降保守性鲁棒控制器,得到最大推力状态降保守性鲁棒控制器组解算模块。设计的最大推力状态降保守性鲁棒控制器采用小摄动不确定性发动机模型,消除了发动机不确定性中的退化项,降低了不确定模型的摄动范围,降低了鲁棒增益调度控制器的保守性。退化参数估计回路实现了退化参数的可靠估计,利用退化参数实现发动机性能退化时的增益调度控制。本发明实现发动机最大推力状态的降保守性鲁棒控制,具有强的鲁棒性并且保守性低,最大限度的提高发动机在最大推力状态的性能,使发动机在最大推力状态不仅稳定工作,并且提高发动机最大推力状态的推力,提高战斗机的机动性能。Compared with the prior art, the conservative robust controller of the aero-engine maximum thrust state reduction of the present invention utilizes the design method of the traditional robust controller, adds a degradation parameter estimation loop, and improves the gain scheduling controller group. , a conservative robust controller for the maximum thrust state drop under a certain degree of engine degradation is added, and the maximum thrust state drop conservative robust controller group solution module is obtained. The designed maximum thrust state-drop conservative robust controller adopts the engine model with small perturbation uncertainty, which eliminates the degradation term in the engine uncertainty, reduces the perturbation range of the uncertain model, and reduces the robust gain scheduling control. Conservativeness of the device. The degradation parameter estimation loop realizes the reliable estimation of the degradation parameters, and uses the degradation parameters to realize the gain scheduling control when the engine performance is degraded. The invention realizes the conservative and robust control of the maximum thrust state of the engine, has strong robustness and low conservatism, and maximizes the performance of the engine in the maximum thrust state, so that the engine not only works stably in the maximum thrust state, but also improves the performance of the engine. The thrust of the engine at the maximum thrust state improves the maneuverability of the fighter.

本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。Additional aspects and advantages of the present invention will be set forth, in part, from the following description, and in part will be apparent from the following description, or may be learned by practice of the invention.

附图说明Description of drawings

本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:The above and/or additional aspects and advantages of the present invention will become apparent and readily understood from the following description of embodiments taken in conjunction with the accompanying drawings, wherein:

图1是本发明航空发动机最大推力状态降保守性鲁棒控制器的结构简图;Fig. 1 is a schematic diagram of the structure of the aero-engine maximum thrust state drop conservative robust controller of the present invention;

图2是本实施例退化参数调度控制回路中化参数估计回路的结构示意图;FIG. 2 is a schematic structural diagram of a quantification parameter estimation loop in the degradation parameter scheduling control loop of the present embodiment;

图3是本实施例退化参数估计回路中卡尔曼滤波器的结构示意图;3 is a schematic structural diagram of a Kalman filter in the degradation parameter estimation loop of the present embodiment;

图4是发动机模型摄动结构图;Fig. 4 is the perturbation structure diagram of the engine model;

图5是退化项分离的发动机模型摄动结构图;Fig. 5 is the engine model perturbation structure diagram with the separation of the degradation term;

图6是退化后新的发动机模型摄动结构图;Fig. 6 is the perturbation structure diagram of the new engine model after degradation;

图7是不确定模型结构示意图。Figure 7 is a schematic diagram of the structure of the uncertain model.

具体实施方式Detailed ways

战斗机由于需要实现高机动性,发动机的最大推力状态的性能至关重要。传统的鲁棒控制器虽然可以对发动机在最大推力状态实现稳定控制,然而,它们是非常保守的,因为它们将发动机退化看作发动机模型的不确定性进行鲁棒控制器的设计,这严重降低了发动机的性能。针对这一问题,下面给出本发明的分析研究过程。Due to the need to achieve high maneuverability in fighter aircraft, the performance of the engine's maximum thrust state is critical. Although the traditional robust controllers can achieve stable control of the engine at the maximum thrust state, however, they are very conservative because they regard engine degradation as the uncertainty of the engine model for robust controller design, which seriously reduces the engine performance. In view of this problem, the analysis and research process of the present invention is given below.

1、发动机性能退化的估计1. Estimation of engine performance degradation

发动机性能退化是指发动机经过多次循环运行后,由于自然磨损、疲劳、积垢等原因造成的正常老化现象。此时,有些发动机的性能会慢慢偏离额定状态。以涡轮部件为例,当它与发动机一起工作多个周期时,其工作效率会缓慢下降。将高温高压气体转化为机械能的能力将会降低,发动机在一个工作点处的线性化模型也会改变。Engine performance degradation refers to the normal aging phenomenon caused by natural wear, fatigue, fouling, etc. of the engine after many cycles of operation. At this time, the performance of some engines will slowly deviate from the rated state. Take the turbine component, for example, as it works with the engine for multiple cycles, its efficiency slowly decreases. The ability to convert high temperature and high pressure gas into mechanical energy will be reduced and the linearized model of the engine at one operating point will change.

发动机性能退化的最终特征是不同转子部件的工作效率和流量的变化,风扇、压气机、主燃烧、高压涡轮和低压涡轮部件的效率系数或流量系数的变化可以表征发动机性能的退化,风扇、压气机、主燃烧室、高压涡轮和低压涡轮部件的效率系数或流量系数被称为退化参数或健康参数。The final characteristic of engine performance degradation is the change in operating efficiency and flow rate of different rotor components. Changes in the efficiency coefficient or flow coefficient of fan, compressor, main combustion, high pressure turbine and low pressure turbine components can characterize engine performance degradation. The efficiency or flow coefficients of the engine, main combustion chamber, high pressure turbine and low pressure turbine components are referred to as degradation parameters or health parameters.

基于部件法,建立带退化参数的发动机非线性模型Building a nonlinear engine model with degradation parameters based on the component method

Figure BDA0002439600450000051
Figure BDA0002439600450000051

y=g(x,u,h)y=g(x,u,h)

其中

Figure BDA0002439600450000052
为控制输入向量,
Figure BDA0002439600450000053
为状态向量,
Figure BDA0002439600450000054
为输出向量,
Figure BDA0002439600450000055
为退化参数向量,f(·)为表示系统动态的n维可微非线性向量函数,g(·)为产生系统输出的m维可微非线性向量函数。in
Figure BDA0002439600450000052
For the control input vector,
Figure BDA0002439600450000053
is the state vector,
Figure BDA0002439600450000054
is the output vector,
Figure BDA0002439600450000055
is the degradation parameter vector, f(·) is the n-dimensional differentiable nonlinear vector function representing the system dynamics, and g(·) is the m-dimensional differentiable nonlinear vector function that produces the system output.

将退化参数h看作发动机的控制输入,采用小扰动法或拟合法对发动机非线性模型在健康稳态参考点处进行线性化。The degradation parameter h is regarded as the control input of the engine, and the nonlinear model of the engine is linearized at the healthy steady-state reference point by the small disturbance method or the fitting method.

Figure BDA0002439600450000056
Figure BDA0002439600450000056

其中in

A′=A,B′=(B L),C′=C,A'=A, B'=(B L), C'=C,

D′=(D M),Δu′=(ΔuΔh)T D′=(DM),Δu′=(ΔuΔh) T

w为系统噪声,v为测量噪声,h为退化参数,Δh=h-h0;上述w与v皆为不相关的高斯白噪声,其均值均为0,协方差矩阵为对角阵Q和R,即满足条件如下:w is the system noise, v is the measurement noise, h is the degradation parameter, Δh=hh 0 ; the above-mentioned w and v are both uncorrelated Gaussian white noise, their mean values are 0, and the covariance matrices are diagonal matrices Q and R, That is, the following conditions are met:

E(w)=0E[wwT]=QE(w)=0E[ww T ]=Q

E(v)=0E[vvT]=RE(v)=0E[vv T ]=R

Δ表示该参数的变化量,h0表示发动机初始状态退化参数。Δ represents the variation of this parameter, and h 0 represents the initial state degradation parameter of the engine.

进一步得到了反映发动机性能退化的增广线性状态变量模型Further, an augmented linear state variable model reflecting engine performance degradation is obtained

Figure BDA0002439600450000061
Figure BDA0002439600450000061

其中系数矩阵可由下式得到:The coefficient matrix can be obtained by the following formula:

Figure BDA0002439600450000062
Figure BDA0002439600450000062

Figure BDA0002439600450000063
Figure BDA0002439600450000063

这些系数在发动机不同的工作状态具有不同的值。These coefficients have different values in different operating states of the engine.

实际上,退化参数很难测量,甚至不可能测量,而发动机各部分的压力、温度、转速等参数比较容易通过测量得到,通常称为“测量参数”,主要包括进气道出口、风扇出口、压气机出口、高压涡轮后、低压涡轮后的温度和压力,风扇转速和压气机转速。当发动机工作环境不发生变化时,退化参数的变化会引起被测参数的相应变化,二者之间存在气动热力学关系。因此,可以设计最优估计滤波器,通过测量参数来实现退化参数的最优估计。In fact, the degradation parameters are difficult to measure, or even impossible to measure, while the pressure, temperature, speed and other parameters of each part of the engine are relatively easy to obtain by measurement, which are usually called "measured parameters", mainly including the outlet of the intake duct, the outlet of the fan, Temperature and pressure at compressor outlet, after high pressure turbine, after low pressure turbine, fan speed and compressor speed. When the working environment of the engine does not change, the change of the degradation parameter will cause the corresponding change of the measured parameter, and there is an aero-thermodynamic relationship between the two. Therefore, an optimal estimation filter can be designed to achieve the optimal estimation of the degradation parameters by measuring the parameters.

由于发动机的性能退化过程相对较慢,可以做出以下合理假设,即Δh的变化率

Figure BDA0002439600450000064
将退化参数进一步转化为状态变量,可以得到Since the engine's performance degradation process is relatively slow, a reasonable assumption can be made that the rate of change of Δh
Figure BDA0002439600450000064
By further transforming the degradation parameters into state variables, we can get

Figure BDA0002439600450000065
Figure BDA0002439600450000065

其中in

Figure BDA0002439600450000066
Figure BDA0002439600450000066

Figure BDA0002439600450000067
Figure BDA0002439600450000067

建立的退化参数估计回路主要由两部分组成,一部分是基于性能退化的非线性机载发动机模型,另一部分是由最大推力状态处模型和稳态点对应的卡尔曼滤波器组成的最大推力状态处卡尔曼滤波器。基本工作原理是将非线性机载发动机模型的输出作为最大推力状态处卡尔曼滤波器的稳态参考值,并扩展退化参数,通过最大推力状态处卡尔曼滤波器进行在线实时估计,最后反馈给非线性机载发动机模型进行在线实时更新。实现对实际发动机的实时跟踪,建立发动机的机载自适应模型。The established degradation parameter estimation loop is mainly composed of two parts, one part is the nonlinear airborne engine model based on performance degradation, and the other part is the maximum thrust state composed of the model at the maximum thrust state and the Kalman filter corresponding to the steady state point. Kalman filter. The basic working principle is to take the output of the nonlinear airborne engine model as the steady-state reference value of the Kalman filter at the maximum thrust state, expand the degradation parameters, conduct online real-time estimation through the Kalman filter at the maximum thrust state, and finally feed back to Online real-time updating of nonlinear airborne engine models. Real-time tracking of the actual engine is realized, and an airborne adaptive model of the engine is established.

卡尔曼估计方程为:The Kalman estimator equation is:

Figure BDA0002439600450000071
Figure BDA0002439600450000071

K为卡尔曼滤波的增益,满足

Figure BDA0002439600450000072
P为Ricati方程
Figure BDA0002439600450000073
的解;利用非线性机载模型输出的健康稳态参考值(xaug,NOBEM,yNOBEM)作为式K is the gain of the Kalman filter, satisfying
Figure BDA0002439600450000072
P is the Ricati equation
Figure BDA0002439600450000073
The solution of ; use the healthy steady-state reference values (x aug, NOBEM , y NOBEM ) output by the nonlinear airborne model as the formula

Figure BDA0002439600450000074
Figure BDA0002439600450000074

的初值,可得计算公式:The initial value of , the calculation formula can be obtained:

Figure BDA0002439600450000075
Figure BDA0002439600450000075

根据该计算公式可以得到发动机的退化参数h。According to this calculation formula, the degradation parameter h of the engine can be obtained.

2.具有退化参数的不确定模型的鲁棒控制器设计2. Robust controller design for uncertain models with degenerate parameters

任何实际系统都不可避免地存在不确定性,它可以分为两类:扰动信号和模型不确定性。扰动信号包括干扰、噪声等。模型的不确定性代表了数学模型与实际对象之间的差异。Uncertainty is inevitable in any real system, and it can be divided into two categories: disturbance signal and model uncertainty. The disturbance signal includes interference, noise, and the like. Model uncertainty represents the difference between the mathematical model and the actual object.

模型不确定性可能有几个原因:线性模型中总有一些参数是有误差的;线性模型中的参数可能由于非线性或工作条件的变化而变化;建模时人为的简化;由于磨损等因素发动机性能的退化。There may be several reasons for model uncertainty: there are always some parameters in the linear model that are in error; parameters in the linear model may vary due to nonlinearity or changes in operating conditions; artificial simplification during modeling; due to factors such as wear and tear Degradation of engine performance.

不确定性可能会对控制系统的稳定性和性能产生不利影响。Uncertainty can adversely affect the stability and performance of the control system.

实际的发动机和标称模型(标称模型是一个常规的不带退化参数的发动机非线性模型)之间的误差可以表示为一个摄动块Δ。请参阅图4,在标称模型加入摄动块建立发动机不确定模型The error between the actual engine and the nominal model (the nominal model is a conventional nonlinear model of the engine without degradation parameters) can be expressed as a perturbation block Δ. Please refer to Figure 4, adding a perturbation block to the nominal model to establish an engine uncertainty model

Figure BDA0002439600450000076
Figure BDA0002439600450000076

Figure BDA0002439600450000077
Figure BDA0002439600450000077

它也可以表示为It can also be expressed as

G(s)=[I+Δ(s)]Gnom(s)G(s)=[I+Δ(s)]G nom (s)

式中G(s)为发动机的不确定模型,Gnom(s)为标称模型,Δ(s)为摄动块。where G(s) is the uncertainty model of the engine, Gnom (s) is the nominal model, and Δ(s) is the perturbation block.

摄动块Δ(s)包含性能退化,请参阅图5,可以通过测量参数进行预测。将摄动块Δ(s)分为不含发动机性能退化的摄动块Δh(s)和退化参数。请参阅图6,在标称模型加入不含发动机性能退化的摄动块Δh(s)与退化参数,将发动机不确定模型表示为The perturbation block Δ(s) contains performance degradation, see Figure 5, which can be predicted by measuring parameters. The perturbation block Δ(s) is divided into the perturbation block Δh (s) without engine performance degradation and the degradation parameters. Referring to Figure 6, the perturbation block Δh (s) and the degradation parameters without engine performance degradation are added to the nominal model, and the engine uncertainty model is expressed as

Figure BDA0002439600450000081
Figure BDA0002439600450000081

Figure BDA0002439600450000082
Figure BDA0002439600450000082

它也可以表示为G(s)=[I+Δh(s)]Gh_nom(s)It can also be expressed as G(s)=[I+ Δh (s)]G h_nom (s)

式中Δh(s)为不含发动机性能退化的摄动块,Gh_nom(s)为在发动机性能退化状态h下的新的标称模型,满足where Δ h (s) is the perturbation block without engine performance degradation, G h_nom (s) is the new nominal model under the engine performance degradation state h, satisfying

G(s)=[I+Δ(s)]Gnom(s)G(s)=[I+Δ(s)]G nom (s)

=[I+Δh(s)+h(s)]Gnom(s)=[I+ Δh (s)+h(s)] Gnom (s)

=[I+Δh(s)]Gh_nom(s)=[I+ Δh (s)]G h_nom (s)

我们可以得到,

Figure BDA0002439600450000083
we can get,
Figure BDA0002439600450000083

请参阅图7,上、下小圆区域分别代表无退化和性能退化h的发动机线性不确定模型,大圆区域代表一般鲁棒控制器设计中发动机线性不确定模型。在一般鲁棒控制器的设计中,直接将发动机的退化看作是模型中的不确定性,不改变发动机的标称模型。因此,不确定项的不确定半径必须足够大,以容纳退化发动机的不确定模型,使不确定模型的摄动半径过大。本专利针对发动机性能退化h的情况,在此状态下建立了新的标称模型,并以新的标称模型为圆心建立了不确定发动机模型。针对某一退化状态下的新的标称模型,在选择不含发动机性能退化的摄动块Δh(s)时,要选择能够覆盖发动机除退化外所有不确定性的最小摄动半径摄动块。请参阅图7,通过对发动机性能退化的估计,发动机不确定性中摄动块的摄动半径||Δh||=||Δ||-||h||<||Δ||,不确定性模型的摄动范围减小了Referring to Figure 7, the upper and lower small circle areas represent the linear uncertainty model of the engine without degradation and performance degradation h, respectively, and the large circle area represents the linear uncertainty model of the engine in the general robust controller design. In the design of general robust controllers, the degradation of the engine is directly regarded as the uncertainty in the model, and the nominal model of the engine is not changed. Therefore, the uncertainty radius of the uncertainty term must be large enough to accommodate the uncertainty model of the degraded engine, making the perturbation radius of the uncertainty model too large. In this patent, a new nominal model is established under the condition of engine performance degradation h, and an uncertain engine model is established with the new nominal model as the center of the circle. For a new nominal model in a degraded state, when selecting the perturbation block Δh (s) without engine performance degradation, the smallest perturbation radius perturbation that can cover all the uncertainties of the engine except for the degradation should be selected. piece. Referring to Fig. 7, by estimating the degradation of engine performance, the perturbation radius of the perturbation block in the engine uncertainty || Δh ||=||Δ||-||h||<||Δ||, The perturbation range of the uncertainty model is reduced

Figure BDA0002439600450000091
Figure BDA0002439600450000091

最后根据小摄动不确定模型利用传统的鲁棒控制器设计方法设计鲁棒控制器,这里设计的鲁棒控制器保守性更低。Finally, according to the small perturbation uncertainty model, the traditional robust controller design method is used to design a robust controller, and the robust controller designed here is less conservative.

3、控制器的插值3. Interpolation of the controller

这部分说明了图1中的最大推力状态降保守性鲁棒控制器组解算模块通过退化参数调度线性插值获得相应的最大推力状态降保守性鲁棒控制器的调度计算原理。This part explains the scheduling calculation principle of the maximum thrust state drop conservative robust controller group solution module in Fig. 1 to obtain the corresponding maximum thrust state drop conservative robust controller through linear interpolation of degenerate parameter scheduling.

分别在发动机最大推力状态下的正常状态和性能退化hbase状态下设计降保守性鲁棒控制器。这将产生图1中的最大推力状态降保守性鲁棒控制器组解算模块中的控制器Kh、Kh_base A conservative robust controller is designed in the normal state and the performance degraded h base state under the maximum thrust state of the engine, respectively. This will produce the maximum thrust state drop in Fig. 1 The controllers K h , K h_base in the conservative robust controller group solver module

根据公式According to the formula

Figure BDA0002439600450000092
Figure BDA0002439600450000092

计算得到航空发动机当前退化状态适应的最大推力状态降保守性鲁棒控制器Kh,并对发动机进行有效控制。The conservative robust controller K h is obtained by calculating the maximum thrust state adapted to the current degradation state of the aero-engine, and the engine is effectively controlled.

基于上述过程,下面给出本实施例中提出的一种航空发动机最大推力状态降保守性鲁棒控制器,如图1所示,主要包括最大推力状态降保守性鲁棒控制器组解算模块和退化参数估计回路。Based on the above process, an aero-engine maximum thrust state drop conservative robust controller proposed in this embodiment is given below. As shown in Figure 1, it mainly includes a maximum thrust state drop conservative robust controller group solution module and degradation parameter estimation loops.

其中最大推力状态降保守性鲁棒控制器组解算模块、退化参数估计回路与航空发动机本体以及航空发动机上的若干传感器组成退化参数调度控制回路10。The maximum thrust state reduction conservative robust controller group solution module, the degradation parameter estimation loop, the aero-engine body and several sensors on the aero-engine form the degradation parameter scheduling control loop 10 .

所述最大推力状态降保守性鲁棒控制器组解算模块产生控制输入向量u并输出给航空发动机本体,传感器得到航空发动机测量参数y;控制输入向量u以及测量参数y共同输入到退化参数估计回路,退化参数估计回路解算得到航空发动机的退化参数h,并输出到最大推力状态降保守性鲁棒控制器组解算模块。The maximum thrust state degradation conservative robust controller group solution module generates a control input vector u and outputs it to the aero-engine body, and the sensor obtains the aero-engine measurement parameter y; the control input vector u and the measurement parameter y are jointly input to the degradation parameter estimation The loop, the degradation parameter estimation loop solves to obtain the degradation parameter h of the aero-engine, and outputs it to the maximum thrust state degradation conservative robust controller group solution module.

所述最大推力状态降保守性鲁棒控制器组解算模块内设计有两个最大推力状态降保守性鲁棒控制器,所述两个最大推力状态降保守性鲁棒控制器采用以下过程得到:分别在发动机正常状态h1和设定退化程度hbase处,在航空发动机最大推力状态下对包含退化参数的发动机非线性模型进行线性化得到2个线性化模型,对线性化模型加入不含发动机性能退化的摄动块得到小摄动不确定性发动机模型,并对这2个小摄动不确定性发动机模型分别设计鲁棒控制器,作为对应的最大推力状态降保守性鲁棒控制器。所述小摄动不确定性发动机模型消除了发动机不确定性模型中的退化项,降低了不确定模型的摄动范围。The maximum thrust state drop conservative robust controller group solution module is designed with two maximum thrust state drop conservative robust controllers, and the two maximum thrust state drop conservative robust controllers are obtained by using the following process : At the normal state of the engine h 1 and the set degradation degree h base , the nonlinear model of the engine including the degradation parameters is linearized under the maximum thrust state of the aero-engine to obtain two linearized models, and the linearized model is added without The engine model with small perturbation uncertainty is obtained from the perturbation block whose engine performance is degraded, and robust controllers are designed for these two engine models with small perturbation uncertainty as the corresponding maximum thrust state reduction conservative robust controller. . The small perturbation uncertainty engine model eliminates the degradation term in the engine uncertainty model and reduces the perturbation range of the uncertainty model.

所述最大推力状态降保守性鲁棒控制器组解算模块根据输入的退化参数h,利用内部设计的两个最大推力状态降保守性鲁棒控制器计算得到适应的最大推力状态降保守性鲁棒控制器,该最大推力状态降保守性鲁棒控制器根据参考输入r和测量参数y的差值e产生控制输入向量u。The maximum thrust state reduction conservative robust controller group calculation module calculates the adapted maximum thrust state reduction conservative robustness by using two internally designed maximum thrust state reduction conservative robust controllers according to the input degradation parameter h. Rod controller, the maximum thrust state drop conservative robust controller generates a control input vector u according to the difference e between the reference input r and the measured parameter y.

优选的一种具体实现方式,可以根据输入的退化参数h插值得到的适应的最大推力状态降保守性鲁棒控制器:A preferred specific implementation manner, the adaptive maximum thrust state drop conservative robust controller can be obtained by interpolation according to the input degradation parameter h:

根据航空发动机正常状态h1和设定退化程度hbase处的最大推力状态降保守性鲁棒控制器K、Kh_base,通过公式According to the normal state h 1 of the aero-engine and the maximum thrust state at the set degradation degree h base , the conservative robust controllers K, K h_base are reduced by the formula

Figure BDA0002439600450000101
Figure BDA0002439600450000101

计算得到航空发动机当前退化状态适应的最大推力状态降保守性鲁棒控制器KhThe conservative robust controller K h is obtained by calculating the maximum thrust state adapted to the current degradation state of the aero-engine.

所述退化参数估计回路中包括非线性机载发动机模型和最大推力状态处卡尔曼滤波器;The degradation parameter estimation loop includes a nonlinear airborne engine model and a Kalman filter at the maximum thrust state;

所述非线性机载发动机模型为带退化参数的发动机非线性模型:The nonlinear airborne engine model is an engine nonlinear model with degradation parameters:

Figure BDA0002439600450000102
Figure BDA0002439600450000102

y=g(x,u,h)y=g(x,u,h)

其中

Figure BDA0002439600450000103
为控制输入向量,
Figure BDA0002439600450000104
为状态向量,
Figure BDA0002439600450000105
为输出向量,
Figure BDA0002439600450000106
为退化参数向量,f(·)为表示系统动态的n维可微非线性向量函数,g(·)为产生系统输出的m维可微非线性向量函数;非线性机载发动机模型输入为控制输入向量u以及上一周期的退化参数h,其输出的健康稳态参考值(xaug,NOBEM,yNOBEM)作为最大推力状态处卡尔曼滤波器当前周期的估计初始值。in
Figure BDA0002439600450000103
For the control input vector,
Figure BDA0002439600450000104
is the state vector,
Figure BDA0002439600450000105
is the output vector,
Figure BDA0002439600450000106
is the degradation parameter vector, f(·) is the n-dimensional differentiable nonlinear vector function representing the system dynamics, g(·) is the m-dimensional differentiable nonlinear vector function that generates the system output; the input of the nonlinear airborne engine model is the control The input vector u and the degradation parameter h of the previous cycle, and the output healthy steady-state reference value (x aug, NOBEM , y NOBEM ) are used as the estimated initial value of the current cycle of the Kalman filter at the maximum thrust state.

所述最大推力状态处卡尔曼滤波器的输入为测量参数y以及非线性机载发动机模型输出的健康稳态参考值(xaug,NOBEM,yNOBEM),根据公式The input of the Kalman filter at the maximum thrust state is the measured parameter y and the healthy steady-state reference value (x aug, NOBEM , y NOBEM ) output by the nonlinear airborne engine model, according to the formula

Figure BDA0002439600450000111
Figure BDA0002439600450000111

计算得到当前周期的发动机的退化参数h;其中

Figure BDA0002439600450000112
K为卡尔曼滤波的增益,满足
Figure BDA0002439600450000113
P为Ricati方程
Figure BDA0002439600450000114
的解;系数Aaug和Caug根据公式Calculate the degradation parameter h of the engine in the current cycle; where
Figure BDA0002439600450000112
K is the gain of the Kalman filter, satisfying
Figure BDA0002439600450000113
P is the Ricati equation
Figure BDA0002439600450000114
The solution of ; the coefficients A aug and C aug according to the formula

Figure BDA0002439600450000115
Figure BDA0002439600450000115

确定,而A、C、L、M是将退化参数h看作发动机的控制输入,并对非线性机载发动机模型在健康稳态参考点处进行线性化得到的反映发动机性能退化的增广线性状态变量模型Determined, while A, C, L, and M are the augmented linearity reflecting the degradation of engine performance obtained by taking the degradation parameter h as the control input of the engine, and linearizing the nonlinear airborne engine model at the healthy steady-state reference point state variable model

Figure BDA0002439600450000116
的系数:
Figure BDA0002439600450000116
The coefficient of :

Figure BDA0002439600450000117
Figure BDA0002439600450000117

Figure BDA0002439600450000118
Figure BDA0002439600450000118

w为系统噪声,v为测量噪声,相应的协方差矩阵为对角阵Q和R。w is the system noise, v is the measurement noise, and the corresponding covariance matrices are the diagonal matrices Q and R.

尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。Although the embodiments of the present invention have been shown and described above, it should be understood that the above embodiments are exemplary and should not be construed as limiting the present invention, and those of ordinary skill in the art will not depart from the principles and spirit of the present invention Variations, modifications, substitutions, and alterations to the above-described embodiments are possible within the scope of the present invention without departing from the scope of the present invention.

Claims (4)

1.一种航空发动机最大推力状态降保守性鲁棒控制器,其特征在于:包括最大推力状态降保守性鲁棒控制器组解算模块和退化参数估计回路;1. an aero-engine maximum thrust state drops conservative robust controller, it is characterized in that: comprise maximum thrust state drop conservative robust controller group solution module and degradation parameter estimation loop; 其中最大推力状态降保守性鲁棒控制器组解算模块、退化参数估计回路与航空发动机本体以及航空发动机上的若干传感器组成退化参数调度控制回路;Among them, the maximum thrust state reduction conservative robust controller group solution module, the degradation parameter estimation loop, the aero-engine body and several sensors on the aero-engine form the degradation parameter scheduling control loop; 所述最大推力状态降保守性鲁棒控制器组解算模块产生控制输入向量u并输出给航空发动机本体,传感器得到航空发动机测量参数y;控制输入向量u以及测量参数y共同输入到退化参数估计回路,退化参数估计回路解算得到航空发动机的退化参数h,并输出到最大推力状态降保守性鲁棒控制器组解算模块;The maximum thrust state degradation conservative robust controller group solution module generates a control input vector u and outputs it to the aero-engine body, and the sensor obtains the aero-engine measurement parameter y; the control input vector u and the measurement parameter y are jointly input to the degradation parameter estimation The loop, the degradation parameter estimation loop solves the aero-engine degradation parameter h, and outputs it to the maximum thrust state degradation conservative robust controller group solution module; 所述最大推力状态降保守性鲁棒控制器组解算模块内设计有两个最大推力状态降保守性鲁棒控制器,所述两个最大推力状态降保守性鲁棒控制器采用以下过程得到:分别在发动机正常状态h1和设定退化程度hbase处,在航空发动机最大推力状态下对包含退化参数的发动机非线性模型进行线性化得到2个线性化模型,对线性化模型加入不含发动机性能退化的摄动块得到小摄动不确定性发动机模型,并对这2个小摄动不确定性发动机模型分别设计鲁棒控制器,作为对应的最大推力状态降保守性鲁棒控制器;The maximum thrust state drop conservative robust controller group solution module is designed with two maximum thrust state drop conservative robust controllers, and the two maximum thrust state drop conservative robust controllers are obtained by using the following process : At the normal state of the engine h 1 and the set degradation degree h base , the nonlinear model of the engine including the degradation parameters is linearized under the maximum thrust state of the aero-engine to obtain two linearized models, and the linearized model is added without The engine model with small perturbation uncertainty is obtained from the perturbation block whose engine performance is degraded, and robust controllers are designed for these two engine models with small perturbation uncertainty as the corresponding maximum thrust state reduction conservative robust controller. ; 所述最大推力状态降保守性鲁棒控制器组解算模块根据输入的退化参数h,利用内部设计的两个最大推力状态降保守性鲁棒控制器计算得到适应的最大推力状态降保守性鲁棒控制器,该最大推力状态降保守性鲁棒控制器根据参考输入r和测量参数y的差值e产生控制输入向量u。The maximum thrust state reduction conservative robust controller group calculation module calculates the adapted maximum thrust state reduction conservative robustness by using two internally designed maximum thrust state reduction conservative robust controllers according to the input degradation parameter h. Rod controller, the maximum thrust state drop conservative robust controller generates a control input vector u according to the difference e between the reference input r and the measured parameter y. 2.根据权利要求1所述一种航空发动机最大推力状态降保守性鲁棒控制器,其特征在于:所述退化参数估计回路中包括非线性机载发动机模型和最大推力状态处卡尔曼滤波器;2. A kind of aero-engine maximum thrust state drop conservative robust controller according to claim 1, it is characterized in that: described degradation parameter estimation loop comprises nonlinear airborne engine model and Kalman filter at maximum thrust state ; 所述非线性机载发动机模型为带退化参数的发动机非线性模型:The nonlinear airborne engine model is an engine nonlinear model with degradation parameters:
Figure FDA0002439600440000011
Figure FDA0002439600440000011
y=g(x,u,h)y=g(x,u,h) 其中
Figure FDA0002439600440000012
为控制输入向量,
Figure FDA0002439600440000013
为状态向量,
Figure FDA0002439600440000014
为输出向量,
Figure FDA0002439600440000015
为退化参数向量,f(·)为表示系统动态的n维可微非线性向量函数,g(·)为产生系统输出的m维可微非线性向量函数;非线性机载发动机模型输入为控制输入向量u以及上一周期的退化参数h,其输出的健康稳态参考值(xaug,NOBEM,yNOBEM)作为最大推力状态处卡尔曼滤波器当前周期的估计初始值;
in
Figure FDA0002439600440000012
For the control input vector,
Figure FDA0002439600440000013
is the state vector,
Figure FDA0002439600440000014
is the output vector,
Figure FDA0002439600440000015
is the degradation parameter vector, f(·) is the n-dimensional differentiable nonlinear vector function representing the system dynamics, g(·) is the m-dimensional differentiable nonlinear vector function that generates the system output; the input of the nonlinear airborne engine model is the control The input vector u and the degradation parameter h of the previous cycle, and the output healthy steady-state reference value (x aug, NOBEM , y NOBEM ) is used as the estimated initial value of the current cycle of the Kalman filter at the maximum thrust state;
所述最大推力状态处卡尔曼滤波器的输入为测量参数y以及非线性机载发动机模型输出的健康稳态参考值(xaug,NOBEM,yNOBEM),根据公式The input of the Kalman filter at the maximum thrust state is the measured parameter y and the healthy steady-state reference value (x aug, NOBEM , y NOBEM ) output by the nonlinear airborne engine model, according to the formula
Figure FDA0002439600440000021
Figure FDA0002439600440000021
计算得到当前周期的发动机的退化参数h;其中
Figure FDA0002439600440000022
K为卡尔曼滤波的增益,满足
Figure FDA0002439600440000023
P为Ricati方程
Figure FDA0002439600440000024
的解;系数Aaug和Caug根据公式
Calculate the degradation parameter h of the engine in the current cycle; where
Figure FDA0002439600440000022
K is the gain of the Kalman filter, satisfying
Figure FDA0002439600440000023
P is the Ricati equation
Figure FDA0002439600440000024
The solution of ; the coefficients A aug and C aug according to the formula
Figure FDA0002439600440000025
Caug=(C M)
Figure FDA0002439600440000025
C aug = (CM)
确定,而A、C、L、M是将退化参数h看作发动机的控制输入,并对非线性机载发动机模型在健康稳态参考点处进行线性化得到的反映发动机性能退化的增广线性状态变量模型Determined, while A, C, L, and M are the augmented linearity reflecting the degradation of engine performance obtained by taking the degradation parameter h as the control input of the engine, and linearizing the nonlinear airborne engine model at the healthy steady-state reference point state variable model
Figure FDA0002439600440000026
Figure FDA0002439600440000026
的系数:The coefficient of :
Figure FDA0002439600440000027
Figure FDA0002439600440000027
Figure FDA0002439600440000028
Figure FDA0002439600440000028
w为系统噪声,v为测量噪声,相应的协方差矩阵为对角阵Q和R。w is the system noise, v is the measurement noise, and the corresponding covariance matrices are the diagonal matrices Q and R.
3.根据权利要求1所述一种航空发动机最大推力状态降保守性鲁棒控制器,其特征在于:所述最大推力状态降保守性鲁棒控制器组解算模块根据航空发动机正常状态h1和设定退化程度hbase处的最大推力状态降保守性鲁棒控制器K、Kh_base,通过公式3. a kind of aero-engine maximum thrust state drops conservative robust controller according to claim 1, it is characterized in that: described maximum thrust state drops conservative robust controller group solution module according to aero-engine normal state h 1 and the maximum thrust state drop conservative robust controllers K, K h_base at the set degradation degree h base , through the formula
Figure FDA0002439600440000031
Figure FDA0002439600440000031
计算得到航空发动机当前退化状态适应的最大推力状态降保守性鲁棒控制器KhThe conservative robust controller K h is obtained by calculating the maximum thrust state adapted to the current degradation state of the aero-engine.
4.根据权利要求1所述一种航空发动机最大推力状态降保守性鲁棒控制器,其特征在于:所述测量参数包括进气道出口、风扇出口、压气机出口、高压涡轮后、低压涡轮后的温度和压力,风扇转速和压气机转速。4. a kind of aero-engine maximum thrust state drop conservative robust controller according to claim 1, is characterized in that: described measurement parameter comprises inlet port outlet, fan outlet, compressor outlet, after high pressure turbine, low pressure turbine After the temperature and pressure, fan speed and compressor speed.
CN202010261759.0A 2020-04-04 2020-04-04 Robust controller for reducing conservative maximum thrust state of aero-engine Active CN111456856B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010261759.0A CN111456856B (en) 2020-04-04 2020-04-04 Robust controller for reducing conservative maximum thrust state of aero-engine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010261759.0A CN111456856B (en) 2020-04-04 2020-04-04 Robust controller for reducing conservative maximum thrust state of aero-engine

Publications (2)

Publication Number Publication Date
CN111456856A true CN111456856A (en) 2020-07-28
CN111456856B CN111456856B (en) 2023-03-24

Family

ID=71681638

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010261759.0A Active CN111456856B (en) 2020-04-04 2020-04-04 Robust controller for reducing conservative maximum thrust state of aero-engine

Country Status (1)

Country Link
CN (1) CN111456856B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112947064A (en) * 2021-01-21 2021-06-11 西北工业大学 Aero-engine maximum thrust control optimization method considering gas circuit component faults

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20010021900A1 (en) * 1998-09-28 2001-09-13 Aspen Technology, Inc. Robust steady-state target calculation for model predictive control
US6539783B1 (en) * 1998-12-28 2003-04-01 General Electric Co. Methods and apparatus for estimating engine health
US20040123600A1 (en) * 2002-11-13 2004-07-01 Brunell Brent Jerome Adaptive model-based control systems and methods for controlling a gas turbine
EP1538319A1 (en) * 2003-12-05 2005-06-08 General Electric Company Apparatus for model predictive control of aircraft gas turbine engines
CA2491208A1 (en) * 2004-03-02 2005-09-02 General Electric Company Model-based control systems and methods for gas turbine engines
JP2008180225A (en) * 2008-03-10 2008-08-07 Hitachi Ltd Engine control device
US20090198350A1 (en) * 2008-01-31 2009-08-06 Fisher-Rosemount Systems, Inc Robust adaptive model predictive controller with tuning to compensate for model mismatch
CA2697760A1 (en) * 2009-05-28 2010-11-28 General Electric Company Real-time scheduling of linear models for control and estimation
US20100318336A1 (en) * 2009-06-13 2010-12-16 Falangas Eric T Method of modeling dynamic characteristics of a flight vehicle
CA2864821A1 (en) * 2012-02-15 2013-08-22 Rolls-Royce Corporation Gas turbine engine performance seeking control
CN106647253A (en) * 2016-09-28 2017-05-10 南京航空航天大学 Aero-engine distributed control system multi-performance robust tracking control method
CN106960084A (en) * 2017-03-06 2017-07-18 西北工业大学 A kind of aero-engine limitation protector method for designing with risk assessment of transfiniting
CN108628330A (en) * 2018-05-09 2018-10-09 南京理工大学 A kind of spacecraft amplitude limit Adaptive Attitude collaboration fault tolerant control method
CN108803336A (en) * 2018-06-28 2018-11-13 南京航空航天大学 A kind of adaptive LQG/LTR controller design methods of aero-engine
CN109441644A (en) * 2018-12-11 2019-03-08 大连理工大学 A kind of fanjet stable state transition state multi-variant control method based on active disturbance rejection theory
CN110161855A (en) * 2019-05-21 2019-08-23 中国电子科技集团公司第三十八研究所 A kind of design method based on robust servo gain scheduling unmanned aerial vehicle (UAV) control device

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20010021900A1 (en) * 1998-09-28 2001-09-13 Aspen Technology, Inc. Robust steady-state target calculation for model predictive control
US6539783B1 (en) * 1998-12-28 2003-04-01 General Electric Co. Methods and apparatus for estimating engine health
US20040123600A1 (en) * 2002-11-13 2004-07-01 Brunell Brent Jerome Adaptive model-based control systems and methods for controlling a gas turbine
EP1538319A1 (en) * 2003-12-05 2005-06-08 General Electric Company Apparatus for model predictive control of aircraft gas turbine engines
CA2491208A1 (en) * 2004-03-02 2005-09-02 General Electric Company Model-based control systems and methods for gas turbine engines
US20090198350A1 (en) * 2008-01-31 2009-08-06 Fisher-Rosemount Systems, Inc Robust adaptive model predictive controller with tuning to compensate for model mismatch
JP2008180225A (en) * 2008-03-10 2008-08-07 Hitachi Ltd Engine control device
CA2697760A1 (en) * 2009-05-28 2010-11-28 General Electric Company Real-time scheduling of linear models for control and estimation
US20100318336A1 (en) * 2009-06-13 2010-12-16 Falangas Eric T Method of modeling dynamic characteristics of a flight vehicle
CA2864821A1 (en) * 2012-02-15 2013-08-22 Rolls-Royce Corporation Gas turbine engine performance seeking control
CN106647253A (en) * 2016-09-28 2017-05-10 南京航空航天大学 Aero-engine distributed control system multi-performance robust tracking control method
CN106960084A (en) * 2017-03-06 2017-07-18 西北工业大学 A kind of aero-engine limitation protector method for designing with risk assessment of transfiniting
CN108628330A (en) * 2018-05-09 2018-10-09 南京理工大学 A kind of spacecraft amplitude limit Adaptive Attitude collaboration fault tolerant control method
CN108803336A (en) * 2018-06-28 2018-11-13 南京航空航天大学 A kind of adaptive LQG/LTR controller design methods of aero-engine
CN109441644A (en) * 2018-12-11 2019-03-08 大连理工大学 A kind of fanjet stable state transition state multi-variant control method based on active disturbance rejection theory
CN110161855A (en) * 2019-05-21 2019-08-23 中国电子科技集团公司第三十八研究所 A kind of design method based on robust servo gain scheduling unmanned aerial vehicle (UAV) control device

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
LINFENG GOU等: "A Linearization Model of Turbofan Engine for Intelligent Analysis Towards Industrial Internet of Things", 《IEEE ACCESS》 *
李睿超等: "基于性能退化缓解的双发推力匹配控制", 《航空工程进展》 *
李睿超等: "涡扇发动机性能退化缓解控制与推力设定", 《航空发动机》 *
纪仓囤等: "时变纯延迟发动机模型的参数辨识及控制方法", 《科学技术与工程》 *
许秀玲等: "某型涡扇发动机控制系统传感器故障诊断研究", 《航空动力学报》 *
贾秋生等: "涡扇发动机变参数鲁棒H_∞滤波器设计", 《推进技术》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112947064A (en) * 2021-01-21 2021-06-11 西北工业大学 Aero-engine maximum thrust control optimization method considering gas circuit component faults

Also Published As

Publication number Publication date
CN111456856B (en) 2023-03-24

Similar Documents

Publication Publication Date Title
CN111273554A (en) Two-degree-of-freedom H-infinity controller for conservative state reduction of maximum thrust of aircraft engine
CN111271181A (en) Two-degree-of-freedom [ mu ] controller for conservative gain reduction scheduling of aero-engine
CN111859555A (en) Robust fault-tolerant controller for maximum thrust state of aero-engine with limited input
CN108062428B (en) Method and system for on-line component fault diagnosis of turbofan engine
CN108829928A (en) A kind of turboshaft engine self-adaptive component grade simulation model construction method
CN107357176B (en) Modeling method for test run data of aero-engine
CN103195730B (en) Compressor inlet adjustable guide vane is used to control the method for combustion turbine exhaustion temperature
CN105389427B (en) Aerial engine air passage unit failure detection method based on adaptive particle filter
CN104047748A (en) Active disturbance observation-based fuel pressure controller and control method thereof
CN111880403A (en) Fault-tolerant two-degree-of-freedom [ mu ] controller for maximum thrust state of aircraft engine
CN110516395A (en) A kind of aero-engine control method based on non-linear mould predictive
CN116151020A (en) Aeroengine health estimation method based on improved extended Kalman filtering
CN115903484A (en) Multivariate Robust Controller Optimization Method for Aeroengine Based on Cooperative Game
CN111305954A (en) Input-limited aero-engine conservative robust gain reduction scheduling controller
JP7349670B2 (en) Engine control method, engine control system, and ship
CN111456856A (en) Robust controller for reducing conservative maximum thrust state of aero-engine
CN114123238A (en) Kalman filtering control method for enabling electrolytic aluminum load to participate in power system frequency modulation
CN111852663A (en) Robust Gain Scheduling Controller with Reduced Conservatism for Variable Cycle Engines
CN120447396B (en) A data-driven closed-loop dynamic real-time operation optimization and predictive control method
CN111456857B (en) A two-degree-of-freedom H∞ controller with decreasing conservative gain scheduling for aeroengines
CN111459028A (en) Conservative two-degree-of-freedom μ controller for maximum thrust state drop of aero-engine
CN111443606A (en) Input-limited aero-engine maximum thrust state conservative robust controller
CN116698349B (en) Air intake control method, system and storage medium for engine air intake simulation
CN111443607A (en) Health degradation-based aeroengine conservative robust gain reduction scheduling controller
CN111443595B (en) Aeroengine gain scheduling controller based on health degradation

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
GR01 Patent grant
GR01 Patent grant