CN111271181A - 航空发动机降保守性增益调度二自由度μ控制器 - Google Patents
航空发动机降保守性增益调度二自由度μ控制器 Download PDFInfo
- Publication number
- CN111271181A CN111271181A CN202010261763.7A CN202010261763A CN111271181A CN 111271181 A CN111271181 A CN 111271181A CN 202010261763 A CN202010261763 A CN 202010261763A CN 111271181 A CN111271181 A CN 111271181A
- Authority
- CN
- China
- Prior art keywords
- engine
- degree
- freedom
- conservative
- controller
- 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
Links
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02C—GAS-TURBINE PLANTS; AIR INTAKES FOR JET-PROPULSION PLANTS; CONTROLLING FUEL SUPPLY IN AIR-BREATHING JET-PROPULSION PLANTS
- F02C9/00—Controlling gas-turbine plants; Controlling fuel supply in air- breathing jet-propulsion plants
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive 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/042—Adaptive 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
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling 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
本发明提出一种航空发动机降保守性增益调度二自由度μ控制器,包括降保守性二自由度μ控制器组解算模块和退化参数估计回路。本发明通过新增退化参数估计回路,并对增益调度控制器组进行改进得到降保守性二自由度μ控制器组解算模块。本发明中小摄动不确定性发动机模型消除了发动机不确定性模型中的退化项,降低了不确定模型的摄动范围,降低了鲁棒增益调度控制器的保守性。退化参数估计回路实现了退化参数的可靠估计,进而结合传统的调度参数,实现发动机性能退化时的增益调度控制,最大限度的提高发动机性能退化时增益调度的控制精度,缩短控制系统过渡时间,减小动态偏差及静态偏差,不仅具有较强的鲁棒性而且保守性低,充分发挥发动机的性能。
Description
技术领域
本发明涉及航空发动机控制技术领域,尤其涉及一种航空发动机降保守性增益调度二自由度μ控制器。
背景技术
航空发动机是一个复杂的非线性动力学系统,在范围宽广的飞行包线内工作时,发动机的工作状态随着外部条件和飞行条件的变化而不断变化。针对航空发动机的强非线性和模型的不确定性,现有技术中有提出鲁棒增益调度控制方法,将发动机划分为一系列工作点,并在每一个工作点设计鲁棒控制器,最终采用增益调度的方法选择合适的鲁棒控制器对发动机进行控制。
上述航空发动机鲁棒增益调度控制方法可以对航空发动机进行控制。然而,它们是非常保守的,因为它们将发动机退化看作发动机模型的不确定性进行鲁棒控制器的设计。事实上,发动机的性能退化程度可以通过测量参数来估计,从而消除不确定性模型中的退化项,缩小不确定性模型的范围,降低鲁棒增益调度控制器的保守性,提高发动机的性能。
另外,传统单自由度控制器无法同时兼顾航空发动机控制系统的鲁棒稳定性和鲁棒性能,而二自由度H∞控制器设计方法由于没有考虑摄动结构会使得设计保守性增强。我们引入自由度μ控制器设计方法对航空发动机进行鲁棒控制器设计。在传统μ控制器的基础上加入前置滤波器和反馈控制器,通过调整反馈控制器使扰动抑制的能力达到最佳,在此基础上调整前置滤波器使系统的指令跟踪能力达到最佳。
发明内容
为解决现有技术存在的问题,本发明提出一种航空发动机降保守性增益调度二自由度μ控制器,同时兼顾航空发动机控制系统的鲁棒稳定性和鲁棒性能,且保守性低,将发动机的性能退化程度通过测量参数来估计,从而消除不确定性模型中的退化项,缩小不确定性模型的范围,降低鲁棒增益调度控制器的保守性,能够在发动机性能发生退化的情况下依旧对真实发动机进行良好控制,充分发挥发动机的性能,提高飞机全寿命效能。
本发明的技术方案为:
所述一种航空发动机降保守性增益调度二自由度μ控制器,其特征在于:包括降保守性二自由度μ控制器组解算模块和退化参数估计回路;
其中降保守性二自由度μ控制器组解算模块、退化参数估计回路与航空发动机本体以及航空发动机上的若干传感器组成退化参数调度控制回路;
所述降保守性二自由度μ控制器组解算模块产生控制输入向量u并输出给航空发动机本体,传感器得到航空发动机测量参数y;控制输入向量u以及测量参数y共同输入到退化参数估计回路,退化参数估计回路解算得到航空发动机的退化参数h,并输出到降保守性二自由度μ控制器组解算模块;
降保守性二自由度μ控制器组解算模块与航空发动机本体以及航空发动机上的若干传感器还组成调度参数调度控制回路;由传感器输出调度参数α至降保守性二自由度μ控制器组解算模块;
所述降保守性二自由度μ控制器组解算模块内设计有若干降保守性二自由度μ控制器;所述降保守性二自由度μ控制器包含前置滤波器和反馈控制器,对若干小摄动不确定性发动机模型利用μ综合方法分别设计得到;
所述小摄动不确定性发动机模型是对航空发动机不同设定工作点下的、包含退化参数的航空发动机非线性模型进行线性化后再加入不含发动机性能退化的摄动块得到的;针对某一退化状态下的航空发动机非线性模型,所加入的不含发动机性能退化的摄动块为能够覆盖航空发动机除退化外所有不确定性的最小摄动半径摄动块;
所述降保守性二自由度μ控制器组解算模块根据输入的退化参数h以及调度参数α,利用内部设计的若干降保守性二自由度μ控制器计算得到适应的降保守性二自由度μ控制器,该降保守性二自由度μ控制器根据参考输入r和测量参数y的差值e产生控制输入向量u。
进一步的,所述降保守性二自由度μ控制器组解算模块内设计若干降保守性二自由度μ控制器的过程为:分别在发动机正常状态h1和设定退化程度hbase处,在全飞行包线内根据调度参数α选取n个工作点对包含退化参数的发动机非线性模型进行线性化得到2n个线性化模型,对线性化模型加入不含发动机性能退化的摄动块得到小摄动不确定性发动机模型,并对这2n个小摄动不确定性发动机模型分别设计相应的降保守性二自由度μ控制器从而组成降保守性二自由度μ控制器组。
进一步的,所述退化参数估计回路中包括非线性机载发动机模型和分段线性化卡尔曼滤波器;
所述非线性机载发动机模型为带退化参数的发动机非线性模型:
y=g(x,u,h)
其中为控制输入向量,为状态向量,为输出向量,为退化参数向量,f(·)为表示系统动态的n维可微非线性向量函数,g(·)为产生系统输出的m维可微非线性向量函数;非线性机载发动机模型输入为控制输入向量u以及上一周期的退化参数h,其输出的健康稳态参考值(xaug,NOBEM,yNOBEM)作为分段线性化卡尔曼滤波器当前周期的估计初始值;
所述分段线性化卡尔曼滤波器的输入为测量参数y以及非线性机载发动机模型输出的健康稳态参考值(xaug,NOBEM,yNOBEM),根据公式
确定,而A、C、L、M是将退化参数h看作发动机的控制输入,并对非线性机载发动机模型在健康稳态参考点处进行线性化得到的反映发动机性能退化的增广线性状态变量模型
w为系统噪声,v为测量噪声,相应的协方差矩阵为对角阵Q和R。
进一步的,所述降保守性二自由度μ控制器组解算模块根据输入的退化参数h以及调度参数α插值得到的适应的降保守性二自由度μ控制器。
进一步的,所述降保守性二自由度μ控制器组解算模块先根据航空发动机当前的调度参数α选择前后相邻的两个设定工作点αi和αi+1,并获取两个设定工作点αi和αi+1对应发动机正常状态h1和设定退化程度hbase处的降保守性二自由度μ控制器Ki、Ki+1和根据公式
计算得到考虑航空发动机当前退化参数h后,在两个设定工作点αi和αi+1下的降保守性二自由度μ控制器Ki和Ki+1;再根据公式
计算得到航空发动机当前适应的降保守性二自由度μ控制器K(α)。
进一步的,所述调度参数α包括航空发动机的风扇转速或者压气机转速。
进一步的,所述测量参数包括进气道出口、风扇出口、压气机出口、高压涡轮后、低压涡轮后的温度和压力,风扇转速和压气机转速。
有益效果
与现有技术相比较,本发明的航空发动机降保守性增益调度二自由度μ控制器利用传统增益调度控制器中固有的模块,通过新增退化参数估计回路,并对增益调度控制器组进行了改进,采用发动机一定退化程度下的一组降保守性二自由度μ控制器,得到降保守性二自由度μ控制器组解算模块。设计的小摄动不确定性发动机模型消除了发动机不确定性模型中的退化项,降低了不确定模型的摄动范围,降低了鲁棒增益调度控制器的保守性。退化参数估计回路实现了退化参数的可靠估计,进而结合传统的调度参数,实现发动机性能退化时的增益调度控制,最大限度的提高发动机性能退化时增益调度的控制精度,缩短控制系统的过渡时间,减小其动态偏差及静态偏差,不仅具有较强的鲁棒性而且保守性低,充分发挥发动机的性能。经控制器控制非线性受控系统,以使系统在整个工作范围内,获得理想的动静态控制品质。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1是本发明航空发动机降保守性增益调度二自由度μ控制器的结构简图;
图2是本发明航空发动机降保守性增益调度二自由度μ控制器中降保守性二自由度μ控制器组的结构示意图;
图3是本实施例退化参数调度控制回路中退化参数估计回路的结构示意图;
图4是本实施例退化参数估计回路中卡尔曼滤波器的结构示意图;
图5是发动机模型摄动结构图;
图6是退化项分离的发动机模型摄动结构图;
图7是退化后新的发动机模型摄动结构图;
图8是不确定模型结构示意图;
图9是发动机非线性模型不确定模型示意图;
图10是二自由度控制器的闭环系统结构图;
图11是二自由度控制器的闭环系统标准结构图;
图12标准M-Δ结构图;
图13是闭环系统的结构图;
图14是μ-综合结构图;
图15是控制器插值示意图。
具体实施方式
航空发动机具有强非线性和模型的不确定性,传统的鲁棒增益调度控制直接将发动机退化看作发动机模型的不确定性进行鲁棒控制器的设计,极大的增大了控制器的保守性,严重降低了发动机的性能;另外,传统单自由度控制器无法同时兼顾航空发动机控制系统的鲁棒稳定性和鲁棒性能,而二自由度H∞控制器设计方法由于没有考虑摄动结构会使得设计保守性增强。针对这一问题,下面给出本发明的分析研究过程。
1、发动机性能退化的估计
发动机性能退化是指发动机经过多次循环运行后,由于自然磨损、疲劳、积垢等原因造成的正常老化现象。此时,有些发动机的性能会慢慢偏离额定状态。以涡轮部件为例,当它与发动机一起工作多个周期时,其工作效率会缓慢下降。将高温高压气体转化为机械能的能力将会降低,发动机在一个工作点处的线性化模型也会改变。
发动机性能退化的最终特征是不同转子部件的工作效率和流量的变化,风扇、压气机、主燃烧、高压涡轮和低压涡轮部件的效率系数或流量系数的变化可以表征发动机性能的退化,风扇、压气机、主燃烧室、高压涡轮和低压涡轮部件的效率系数或流量系数被称为退化参数或健康参数。
基于部件法,建立带退化参数的发动机非线性模型
y=g(x,u,h)
将退化参数h看作发动机的控制输入,采用小扰动法或拟合法对发动机非线性模型在健康稳态参考点处进行线性化。
其中
A′=A,B′=(BL),C′=C,
D′=(DM),Δ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表示发动机初始状态退化参数。
进一步得到了反映发动机性能退化的增广线性状态变量模型
其中系数矩阵可由下式得到:
这些系数在发动机不同的工作状态具有不同的值。
实际上,退化参数很难测量,甚至不可能测量,而发动机各部分的压力、温度、转速等参数比较容易通过测量得到,通常称为“测量参数”,主要包括进气道出口、风扇出口、压气机出口、高压涡轮后、低压涡轮后的温度和压力,风扇转速和压气机转速。当发动机工作环境不发生变化时,退化参数的变化会引起被测参数的相应变化,二者之间存在气动热力学关系。因此,可以设计最优估计滤波器,通过测量参数来实现退化参数的最优估计。
其中
建立的退化参数估计回路主要由两部分组成,一部分是基于性能退化的非线性机载发动机模型,另一部分是由分段线性化模型和稳态点对应的卡尔曼滤波器组成的分段线性卡尔曼滤波器。基本工作原理是将非线性机载发动机模型的输出作为分段线性卡尔曼滤波器的稳态参考值,并扩展退化参数,通过分段线性卡尔曼滤波器进行在线实时估计,最后反馈给非线性机载发动机模型进行在线实时更新。实现对实际发动机的实时跟踪,建立发动机的机载自适应模型。
卡尔曼估计方程为:
的初值,可得计算公式:
根据该计算公式可以得到发动机的退化参数h。
2.具有退化参数的不确定模型的鲁棒控制器设计
任何实际系统都不可避免地存在不确定性,它可以分为两类:扰动信号和模型不确定性。扰动信号包括干扰、噪声等。模型的不确定性代表了数学模型与实际对象之间的差异。
模型不确定性可能有几个原因:线性模型中总有一些参数是有误差的;线性模型中的参数可能由于非线性或工作条件的变化而变化;建模时人为的简化;由于磨损等因素发动机性能的退化。
不确定性可能会对控制系统的稳定性和性能产生不利影响。
实际的发动机和标称模型(标称模型是一个常规的不带退化参数的发动机非线性模型)之间的误差可以表示为一个摄动块Δ。请参阅图5,在标称模型加入摄动块建立发动机不确定模型
它也可以表示为
G(s)=[I+Δ(s)]Gnom(s)
式中G(s)为发动机的不确定模型,Gnom(s)为标称模型,Δ(s)为摄动块。
摄动块Δ(s)包含性能退化,请参阅图6,可以通过测量参数进行预测。将摄动块Δ(s)分为不含发动机性能退化的摄动块Δh(s)和退化参数。请参阅图7,在标称模型加入不含发动机性能退化的摄动块Δh(s)与退化参数,将发动机不确定模型表示为
它也可以表示为
G(s)=[I+Δh(s)]Gh_nom(s)
式中Δh(s)为不含发动机性能退化的摄动块,Gh_nom(s)为在发动机性能退化状态h下的新的标称模型,满足
G(s)=[I+Δ(s)]Gnom(s)
=[I+Δh(s)+h(s)]Gnom(s)
=[I+Δh(s)]Gh_nom(s)
请参阅图8,上、下小圆区域分别代表无退化和性能退化h的发动机线性不确定模型,大圆区域代表一般鲁棒控制器设计中发动机线性不确定模型。在一般鲁棒控制器的设计中,直接将发动机的退化看作是模型中的不确定性,不改变发动机的标称模型。因此,不确定项的不确定半径必须足够大,以容纳退化发动机的不确定模型,使不确定模型的摄动半径过大。本专利针对发动机性能退化h的情况,在此状态下建立了新的标称模型,并以新的标称模型为圆心建立了不确定发动机模型。针对某一退化状态下的新的标称模型,在选择不含发动机性能退化的摄动块Δh(s)时,要选择能够覆盖发动机除退化外所有不确定性的最小摄动半径摄动块。请参阅图8,通过对发动机性能退化的估计,发动机不确定性中摄动块的摄动半径||Δh||=||Δ||-||h||<||Δ||,不确定性模型的摄动范围减小了
最后根据小摄动不确定模型利用μ综合方法设计鲁棒控制器,这里设计的鲁棒控制器保守性更低。
3、具有退化参数的增益调度控制设计
增益调度控制的实质是设计一组线性化的控制器,然后将它们有规律地组合起来,从而能够控制非线性系统。具有退化参数的增益调度控制的基本原理是选择一系列的工作点,获得在正常状态和某些性能退化状态的发动机线性化模型并分别设计对应的降保守性二自由度μ控制器得到图1中的降保守性二自由度μ控制器组。
请参阅图9,选择一组调度参数值αi,i=1,2,...,n,代表系统的动态范围,并将飞行包线划分为几个子区间,并将这些点作为工作点。在工作点,有这些方程
利用小扰动法,可以得到各工况点退化参数的线性模型,并得到发动机正常状态和性能退化h状态下的线性模型。
请参阅图9,上、下实线分别表示发动机无退化和有性能退化h的非线性模型。一系列的黑色小圆点表示发动机不同的工作点,在每一个工作点进行线性化得到线性模型。上、下一系列的小虚线圆分别表示没有退化和有退化下的一系列不含退化项的小摄动范围,大虚线圆表示含退化项的大范围摄动。针对发动机正常状态和退化状态下的小摄动不确定线性模型,分别设计一系列降保守性二自由度μ控制器得到图1中的降保守性二自由度μ控制器组。然后,在选定的工作点之间,对控制器增益进行线性内插,使得对于所有的固定参数值,闭环系统都具有良好的性能。参数α是调度参数,这里可以定义为航空发动机的风扇转速或者压气机转速,可以实时测量。控制系统的另一个调度变量是反映发动机性能退化的退化参数h。工作原理是图1中的降保守性二自由度μ控制器组解算模块根据调度参数和退化参数进行线性插值获得相应的降保守性二自由度μ控制器来控制系统。
4.二自由度μ控制器设计
具有二自由度控制器的闭环系统的框图如图10所示。该系统有一个参考输入(r),输出干扰(d)和两个输出误差(z1)和(z2)。系统M0是闭环系统应该匹配的理想模型。在这种结构中,除了内部稳定性要求之外,两个信号e和u将被最小化。信号e显示系统输出和参考模型输出之间的差异。u是控制信号,也与摄动中的鲁棒稳定性有关。在图10中,包含两个加权函数以反映这两个惩罚信号之间的特征。
所谓二自由度控制,就是对使目标值跟踪特性为最优的参数和使外扰抑制特性为最优的参数分别独立进行整定,使用反馈控制器(Ky)来实现内部稳定性、鲁棒稳定性和干扰抑制等,并且在前馈路径上设计另一个控制器(Kr)以满足跟踪需求,最小化整个系统的输出与参考模型M的输出之间的差异,使两特性同时达到最优。
K=[Kr Ky]
二自由度μ控制器设计与二自由度H∞控制器设计的区别在于后者是求解控制器K使得系统的范数达到最小,前者是求解控制器K使得系统的结构奇异值达到最小,如果不确定性系统的摄动是非结构化的,范数最小和结构奇异值最小是等价的,如果不确定性系统的摄动是结构化的,通过使系统范数最小由于没有考虑摄动的已知的对角结构这一条件,因此设计出来的控制器偏保守,无法充分发挥控制系统的鲁棒稳定性和鲁棒性能,而通过使系统的结构奇异值最小设计的控制器K克服了这一不足。
在航空发动机二自由度μ控制器设计问题中,所考虑的不确定性将包括非结构不确定性(如未建模动态)和参数变化。所有这些不确定部分仍然可以从动态中取出,整个系统可以重新排列在(上)线性分式变换F(M,Δ)的标准结构中,不确定性块和标称模型之间的关系可以概括为标准结构,如图12所示,以表示不确定性如何影响正在研究的控制系统的输入/输出关系。不确定块Δ将具有如下的一般形式:
如果没有Δ∈Δ使得det(I-MΔ)=0,则μΔ(M):=0。
当M是如图12所示的互连传递矩阵时,关于Δ的结构奇异值定义为
相应地,不确定性集合可以被定义为
当不确定性结构固定时,为了简洁,我们可以省略μΔ(M)的下标Δ。
结构奇异值的倒数表示频率相关的稳定裕度。下面给出了关于结构不确定性的鲁棒稳定性结果。
很明显,如果不确定性存在于单位球BΔ中,那么鲁棒稳定性条件是μΔ(M(s))≤1。μΔ(M(s))取决于频率,并且在实际应用中在合理范围内以“每个”频率计算。
考虑图13所示的线性分式变换(LFT)中的控制问题。
由P表示的系统是开环连接并表示所有已知元素,包括标称系统模型和性能加权函数以及不确定性加权函数。块Δ是来自集合Δ的不确定元素,其将所有假定的模型不确定性参数化。控制器由K表示。对P的输入是三组信号:由不确定性引起的输入uΔ,参考输入和干扰w以及控制u。产生三组输出:由不确定性,误差或受控输出z和测量值y而输出yΔ。
LFT描述了要控制的一组系统
设计目标是确定一个控制器K,使标称系统稳定;对于所有Δ∈Δ,maxωσ[Δ(jω)]≤1,闭环系统是稳定的并满足
||FU[FL(P,K),Δ]||∞<1
对于给定的任意K,可以使用对线性分式变换FL(P,K)的鲁棒性能测试来测试该性能标准。应该对扩展的不确定结构进行鲁棒性能测试
其中ΔF是一个虚构的非结构化摄动块。具有控制器K的系统可以实现可靠的性能,当且仅当
如图14所示。最后可以通过D-K迭代法进行μ-综合得到控制器K。
5、控制器的插值
这部分说明了图1中的降保守性二自由度μ控制器组解算模块通过调度参数和退化参数调度线性插值获得相应的降保守性二自由度μ控制器的调度计算原理。
分别在发动机正常状态和性能退化hbase状态下设计一系列降保守性二自由度μ控制器,对每个选定的工作点αi进行控制。这将产生图1中的降保守性二自由度μ控制器组解算模块中的控制器
然后根据调度参数α和退化参数h对控制器进行插值,继而使用得到的插值控制器来控制系统。
请参阅图15,根据发动机当前的调度参数α选定周围相邻的两个工作点αi和αi+1,根据发动机在选定的工作点αi的实际退化程度,在发动机性能退化h处的控制器Ki,使用所选的工作点αi处发动机正常状态和性能退化h-base状态的控制器Ki和通过线性插值得到
同样的,可以得到在工作点αi+1处实际退化h处的控制器
我们使用分段线性插值的方法,从降保守性二自由度μ控制器集K1,K2,…,Kn中对每一对控制器之间进行线性插值。得到当前调度参数α当前退化程度h处的降保守性二自由度μ控制器K(α),i=1,2,…,n-1为
根据该公式可以得到某一调度参数某一退化参数下相应的控制器,并对发动机进行有效控制。
基于上述过程,下面给出本实施例中提出的一种航空发动机降保守性增益调度二自由度μ控制器,如图1所示,主要包括降保守性二自由度μ控制器组解算模块和退化参数估计回路。
其中降保守性二自由度μ控制器组解算模块、退化参数估计回路与航空发动机本体以及航空发动机上的若干传感器组成退化参数调度控制回路10。
所述降保守性二自由度μ控制器组解算模块产生控制输入向量u并输出给航空发动机本体,传感器得到航空发动机测量参数y;控制输入向量u以及测量参数y共同输入到退化参数估计回路,退化参数估计回路解算得到航空发动机的退化参数h,并输出到降保守性二自由度μ控制器组解算模块。
降保守性二自由度μ控制器组解算模块与航空发动机本体以及航空发动机上的若干传感器还组成调度参数调度控制回路20;由传感器输出调度参数α至降保守性二自由度μ控制器组解算模块。
所述降保守性二自由度μ控制器组解算模块内设计有若干降保守性二自由度μ控制器;所述降保守性二自由度μ控制器包含前置滤波器和反馈控制器,对若干小摄动不确定性发动机模型利用μ综合方法分别设计得到。
所述小摄动不确定性发动机模型是对航空发动机不同设定工作点下的、包含退化参数的航空发动机非线性模型进行线性化后再加入不含发动机性能退化的摄动块得到的;针对某一退化状态下的航空发动机非线性模型,所加入的不含发动机性能退化的摄动块为能够覆盖航空发动机除退化外所有不确定性的最小摄动半径摄动块。
优选的一种具体实现方式,可以通过以下过程得到设计若干降保守性二自由度μ控制器:分别在发动机正常状态h1和设定退化程度hbase处,在全飞行包线内根据调度参数α选取n个工作点对包含退化参数的发动机非线性模型进行线性化得到2n个线性化模型,对线性化模型加入不含发动机性能退化的摄动块得到小摄动不确定性发动机模型,并对这2n个小摄动不确定性发动机模型分别设计相应的降保守性二自由度μ控制器从而组成降保守性二自由度μ控制器组。
所述降保守性二自由度μ控制器组解算模块根据输入的退化参数h以及调度参数α,利用内部设计的若干降保守性二自由度μ控制器计算得到适应的降保守性二自由度μ控制器,该降保守性二自由度μ控制器根据参考输入r和测量参数y的差值e产生控制输入向量u。
优选的一种具体实现方式,可以根据输入的退化参数h以及调度参数α插值得到的适应的降保守性二自由度μ控制器:
先根据航空发动机当前的调度参数α选择前后相邻的两个设定工作点αi和αi+1,并获取两个设定工作点αi和αi+1对应发动机正常状态h1和设定退化程度hbase处的降保守性二自由度μ控制器Ki、Ki+1和根据公式
计算得到考虑航空发动机当前退化参数h后,在两个设定工作点αi和αi+1下的降保守性二自由度μ控制器Ki和Ki+1;再根据公式
计算得到航空发动机当前适应的降保守性二自由度μ控制器K(α)。
所述退化参数估计回路中包括非线性机载发动机模型和分段线性化卡尔曼滤波器;
所述非线性机载发动机模型为带退化参数的发动机非线性模型:
y=g(x,u,h)
其中为控制输入向量,为状态向量,为输出向量,为退化参数向量,f(·)为表示系统动态的n维可微非线性向量函数,g(·)为产生系统输出的m维可微非线性向量函数;非线性机载发动机模型输入为控制输入向量u以及上一周期的退化参数h,其输出的健康稳态参考值(xaug,NOBEM,yNOBEM)作为分段线性化卡尔曼滤波器当前周期的估计初始值;
所述分段线性化卡尔曼滤波器的输入为测量参数y以及非线性机载发动机模型输出的健康稳态参考值(xaug,NOBEM,yNOBEM),根据公式
确定,而A、C、L、M是将退化参数h看作发动机的控制输入,并对非线性机载发动机模型在健康稳态参考点处进行线性化得到的反映发动机性能退化的增广线性状态变量模型
的系数:
w为系统噪声,v为测量噪声,相应的协方差矩阵为对角阵Q和R。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (7)
1.一种航空发动机降保守性增益调度二自由度μ控制器,其特征在于:包括降保守性二自由度μ控制器组解算模块和退化参数估计回路;
其中降保守性二自由度μ控制器组解算模块、退化参数估计回路与航空发动机本体以及航空发动机上的若干传感器组成退化参数调度控制回路;
所述降保守性二自由度μ控制器组解算模块产生控制输入向量u并输出给航空发动机本体,传感器得到航空发动机测量参数y;控制输入向量u以及测量参数y共同输入到退化参数估计回路,退化参数估计回路解算得到航空发动机的退化参数h,并输出到降保守性二自由度μ控制器组解算模块;
降保守性二自由度μ控制器组解算模块与航空发动机本体以及航空发动机上的若干传感器还组成调度参数调度控制回路;由传感器输出调度参数α至降保守性二自由度μ控制器组解算模块;
所述降保守性二自由度μ控制器组解算模块内设计有若干降保守性二自由度μ控制器;所述降保守性二自由度μ控制器包含前置滤波器和反馈控制器,对若干小摄动不确定性发动机模型利用μ综合方法分别设计得到;
所述小摄动不确定性发动机模型是对航空发动机不同设定工作点下的、包含退化参数的航空发动机非线性模型进行线性化后再加入不含发动机性能退化的摄动块得到的;针对某一退化状态下的航空发动机非线性模型,所加入的不含发动机性能退化的摄动块为能够覆盖航空发动机除退化外所有不确定性的最小摄动半径摄动块;
所述降保守性二自由度μ控制器组解算模块根据输入的退化参数h以及调度参数α,利用内部设计的若干降保守性二自由度μ控制器计算得到适应的降保守性二自由度μ控制器,该降保守性二自由度μ控制器根据参考输入r和测量参数y的差值e产生控制输入向量u。
2.根据权利要求1所述一种航空发动机降保守性增益调度二自由度μ控制器,其特征在于:所述降保守性二自由度μ控制器组解算模块内设计若干降保守性二自由度μ控制器的过程为:分别在发动机正常状态h1和设定退化程度hbase处,在全飞行包线内根据调度参数α选取n个工作点对包含退化参数的发动机非线性模型进行线性化得到2n个线性化模型,对线性化模型加入不含发动机性能退化的摄动块得到小摄动不确定性发动机模型,并对这2n个小摄动不确定性发动机模型分别设计相应的降保守性二自由度μ控制器从而组成降保守性二自由度μ控制器组。
3.根据权利要求1所述一种航空发动机降保守性增益调度二自由度μ控制器,其特征在于:所述退化参数估计回路中包括非线性机载发动机模型和分段线性化卡尔曼滤波器;
所述非线性机载发动机模型为带退化参数的发动机非线性模型:
y=g(x,u,h)
其中为控制输入向量,为状态向量,为输出向量,为退化参数向量,f(·)为表示系统动态的n维可微非线性向量函数,g(·)为产生系统输出的m维可微非线性向量函数;非线性机载发动机模型输入为控制输入向量u以及上一周期的退化参数h,其输出的健康稳态参考值(xaug,NOBEM,yNOBEM)作为分段线性化卡尔曼滤波器当前周期的估计初始值;
所述分段线性化卡尔曼滤波器的输入为测量参数y以及非线性机载发动机模型输出的健康稳态参考值(xaug,NOBEM,yNOBEM),根据公式
确定,而A、C、L、M是将退化参数h看作发动机的控制输入,并对非线性机载发动机模型在健康稳态参考点处进行线性化得到的反映发动机性能退化的增广线性状态变量模型
的系数:
w为系统噪声,v为测量噪声,相应的协方差矩阵为对角阵Q和R。
4.根据权利要求2所述一种航空发动机降保守性增益调度二自由度μ控制器,其特征在于:所述降保守性二自由度μ控制器组解算模块根据输入的退化参数h以及调度参数α插值得到的适应的降保守性二自由度μ控制器。
6.根据权利要求1所述一种航空发动机降保守性增益调度二自由度μ控制器,其特征在于:所述调度参数α包括航空发动机的风扇转速或者压气机转速。
7.根据权利要求1所述一种航空发动机降保守性增益调度二自由度μ控制器,其特征在于:所述测量参数包括进气道出口、风扇出口、压气机出口、高压涡轮后、低压涡轮后的温度和压力,风扇转速和压气机转速。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010261763.7A CN111271181B (zh) | 2020-04-04 | 2020-04-04 | 航空发动机降保守性增益调度二自由度μ控制器 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010261763.7A CN111271181B (zh) | 2020-04-04 | 2020-04-04 | 航空发动机降保守性增益调度二自由度μ控制器 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111271181A true CN111271181A (zh) | 2020-06-12 |
CN111271181B CN111271181B (zh) | 2022-08-19 |
Family
ID=70999739
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010261763.7A Active CN111271181B (zh) | 2020-04-04 | 2020-04-04 | 航空发动机降保守性增益调度二自由度μ控制器 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111271181B (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111608808A (zh) * | 2020-06-15 | 2020-09-01 | 西北工业大学 | 输入受限航空发动机增益调度容错控制器 |
CN111830827A (zh) * | 2020-06-15 | 2020-10-27 | 西北工业大学 | 航空发动机容错增益调度二自由度μ控制器 |
CN111856918A (zh) * | 2020-06-15 | 2020-10-30 | 西北工业大学 | 变循环发动机增益调度控制器 |
CN111856928A (zh) * | 2020-06-15 | 2020-10-30 | 西北工业大学 | 变循环发动机鲁棒增益调度控制器 |
CN111852663A (zh) * | 2020-06-15 | 2020-10-30 | 西北工业大学 | 变循环发动机降保守性鲁棒增益调度控制器 |
CN111856927A (zh) * | 2020-06-15 | 2020-10-30 | 西北工业大学 | 变循环发动机增益调度二自由度μ控制器 |
CN111856919A (zh) * | 2020-06-15 | 2020-10-30 | 西北工业大学 | 航空发动机气路部件故障增益调度容错控制器 |
CN111852662A (zh) * | 2020-06-15 | 2020-10-30 | 西北工业大学 | 航空发动机最大推力状态容错二自由度h∞控制器 |
CN111859555A (zh) * | 2020-06-15 | 2020-10-30 | 西北工业大学 | 输入受限的航空发动机最大推力状态鲁棒容错控制器 |
CN111880403A (zh) * | 2020-06-15 | 2020-11-03 | 西北工业大学 | 航空发动机最大推力状态容错二自由度μ控制器 |
CN112327611A (zh) * | 2020-06-15 | 2021-02-05 | 西北工业大学 | 变循环发动机增益调度二自由度h∞控制器 |
CN112327602A (zh) * | 2020-06-15 | 2021-02-05 | 西北工业大学 | 变循环发动机气路部件故障增益调度容错控制器 |
CN112346336A (zh) * | 2020-06-15 | 2021-02-09 | 西北工业大学 | 航空发动机气路部件故障鲁棒增益调度容错控制器 |
CN112377311A (zh) * | 2020-06-15 | 2021-02-19 | 西北工业大学 | 输入受限的航空发动机鲁棒增益调度容错控制器 |
CN112879165A (zh) * | 2021-01-14 | 2021-06-01 | 西北工业大学 | 考虑气路部件故障的航空发动机加速过程最优控制方法 |
CN112947064A (zh) * | 2021-01-21 | 2021-06-11 | 西北工业大学 | 考虑气路部件故障的航空发动机最大推力控制优化方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1666480A (zh) * | 2002-06-27 | 2005-09-07 | 罗伯特·博世有限公司 | 用于数据传输的方法和装置 |
CN104736372A (zh) * | 2012-10-23 | 2015-06-24 | 罗伯特·博世有限公司 | 用于运行至少一个驾驶员协助系统的方法以及用于至少一个驾驶员协助系统的运行装置 |
EP3346344A1 (en) * | 2017-01-04 | 2018-07-11 | General Electric Company | Method for loop gain sizing of gas turbines |
CN109441644A (zh) * | 2018-12-11 | 2019-03-08 | 大连理工大学 | 一种基于自抗扰理论的涡扇发动机稳态过渡态多变量控制方法 |
CN110008600A (zh) * | 2019-04-09 | 2019-07-12 | 上海理工大学 | 车辆稳定性控制器性能保守性的设计方法 |
CN110513198A (zh) * | 2019-08-13 | 2019-11-29 | 大连理工大学 | 一种涡扇发动机控制系统主动容错控制方法 |
-
2020
- 2020-04-04 CN CN202010261763.7A patent/CN111271181B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1666480A (zh) * | 2002-06-27 | 2005-09-07 | 罗伯特·博世有限公司 | 用于数据传输的方法和装置 |
CN104736372A (zh) * | 2012-10-23 | 2015-06-24 | 罗伯特·博世有限公司 | 用于运行至少一个驾驶员协助系统的方法以及用于至少一个驾驶员协助系统的运行装置 |
EP3346344A1 (en) * | 2017-01-04 | 2018-07-11 | General Electric Company | Method for loop gain sizing of gas turbines |
JP2018138784A (ja) * | 2017-01-04 | 2018-09-06 | ゼネラル・エレクトリック・カンパニイ | ガスタービンのループ利得の調節のための方法 |
CN109441644A (zh) * | 2018-12-11 | 2019-03-08 | 大连理工大学 | 一种基于自抗扰理论的涡扇发动机稳态过渡态多变量控制方法 |
CN110008600A (zh) * | 2019-04-09 | 2019-07-12 | 上海理工大学 | 车辆稳定性控制器性能保守性的设计方法 |
CN110513198A (zh) * | 2019-08-13 | 2019-11-29 | 大连理工大学 | 一种涡扇发动机控制系统主动容错控制方法 |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112327611A (zh) * | 2020-06-15 | 2021-02-05 | 西北工业大学 | 变循环发动机增益调度二自由度h∞控制器 |
CN111856927B (zh) * | 2020-06-15 | 2022-06-10 | 西北工业大学 | 变循环发动机增益调度二自由度μ控制器 |
CN111856918A (zh) * | 2020-06-15 | 2020-10-30 | 西北工业大学 | 变循环发动机增益调度控制器 |
CN111856928A (zh) * | 2020-06-15 | 2020-10-30 | 西北工业大学 | 变循环发动机鲁棒增益调度控制器 |
CN111852663A (zh) * | 2020-06-15 | 2020-10-30 | 西北工业大学 | 变循环发动机降保守性鲁棒增益调度控制器 |
CN111856927A (zh) * | 2020-06-15 | 2020-10-30 | 西北工业大学 | 变循环发动机增益调度二自由度μ控制器 |
CN111856919A (zh) * | 2020-06-15 | 2020-10-30 | 西北工业大学 | 航空发动机气路部件故障增益调度容错控制器 |
CN111608808A (zh) * | 2020-06-15 | 2020-09-01 | 西北工业大学 | 输入受限航空发动机增益调度容错控制器 |
CN111830827B (zh) * | 2020-06-15 | 2022-06-10 | 西北工业大学 | 航空发动机容错增益调度二自由度μ控制器 |
CN111830827A (zh) * | 2020-06-15 | 2020-10-27 | 西北工业大学 | 航空发动机容错增益调度二自由度μ控制器 |
CN111852662A (zh) * | 2020-06-15 | 2020-10-30 | 西北工业大学 | 航空发动机最大推力状态容错二自由度h∞控制器 |
CN112327602A (zh) * | 2020-06-15 | 2021-02-05 | 西北工业大学 | 变循环发动机气路部件故障增益调度容错控制器 |
CN112346336A (zh) * | 2020-06-15 | 2021-02-09 | 西北工业大学 | 航空发动机气路部件故障鲁棒增益调度容错控制器 |
CN112377311A (zh) * | 2020-06-15 | 2021-02-19 | 西北工业大学 | 输入受限的航空发动机鲁棒增益调度容错控制器 |
CN111859555A (zh) * | 2020-06-15 | 2020-10-30 | 西北工业大学 | 输入受限的航空发动机最大推力状态鲁棒容错控制器 |
CN111856928B (zh) * | 2020-06-15 | 2022-06-10 | 西北工业大学 | 变循环发动机鲁棒增益调度控制器 |
CN112327611B (zh) * | 2020-06-15 | 2022-06-10 | 西北工业大学 | 变循环发动机增益调度二自由度h∞控制器 |
CN111880403A (zh) * | 2020-06-15 | 2020-11-03 | 西北工业大学 | 航空发动机最大推力状态容错二自由度μ控制器 |
CN112879165A (zh) * | 2021-01-14 | 2021-06-01 | 西北工业大学 | 考虑气路部件故障的航空发动机加速过程最优控制方法 |
CN112947064A (zh) * | 2021-01-21 | 2021-06-11 | 西北工业大学 | 考虑气路部件故障的航空发动机最大推力控制优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111271181B (zh) | 2022-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111271181B (zh) | 航空发动机降保守性增益调度二自由度μ控制器 | |
CN111273554B (zh) | 航空发动机最大推力状态降保守性二自由度h∞控制器 | |
CN109162813B (zh) | 一种基于迭代学习修正的航空发动机智能转速控制方法 | |
Zhou et al. | An improved multivariable generalized predictive control algorithm for direct performance control of gas turbine engine | |
CN111880403A (zh) | 航空发动机最大推力状态容错二自由度μ控制器 | |
CN110647052B (zh) | 一种变循环发动机模式切换自适应身份证模型构建方法 | |
CN111859555A (zh) | 输入受限的航空发动机最大推力状态鲁棒容错控制器 | |
CN111608808A (zh) | 输入受限航空发动机增益调度容错控制器 | |
CN110516395A (zh) | 一种基于非线性模型预测的航空发动机控制方法 | |
CN111305954B (zh) | 输入受限的航空发动机降保守性鲁棒增益调度控制器 | |
CN111158386B (zh) | 一种干扰自主免疫的无人机姿态控制系统及实现方法 | |
CN111856919A (zh) | 航空发动机气路部件故障增益调度容错控制器 | |
CN111459028B (zh) | 航空发动机最大推力状态降保守性二自由度μ控制器 | |
CN111456857B (zh) | 航空发动机降保守性增益调度二自由度h∞控制器 | |
CN111856927B (zh) | 变循环发动机增益调度二自由度μ控制器 | |
CN111852663A (zh) | 变循环发动机降保守性鲁棒增益调度控制器 | |
CN111830827B (zh) | 航空发动机容错增益调度二自由度μ控制器 | |
CN111456856B (zh) | 航空发动机最大推力状态降保守性鲁棒控制器 | |
CN111443607B (zh) | 基于健康退化的航空发动机降保守性鲁棒增益调度控制器 | |
CN111852662A (zh) | 航空发动机最大推力状态容错二自由度h∞控制器 | |
CN111443595B (zh) | 基于健康退化的航空发动机增益调度控制器 | |
CN112327611B (zh) | 变循环发动机增益调度二自由度h∞控制器 | |
Litt et al. | Evaluation of an outer loop retrofit architecture for intelligent turbofan engine thrust control | |
CN111443606B (zh) | 输入受限的航空发动机最大推力状态降保守性鲁棒控制器 | |
CN112346336A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |