CN106951634B - 一种航空发动机鲁棒跟踪控制器设计方法 - Google Patents
一种航空发动机鲁棒跟踪控制器设计方法 Download PDFInfo
- Publication number
- CN106951634B CN106951634B CN201710164988.9A CN201710164988A CN106951634B CN 106951634 B CN106951634 B CN 106951634B CN 201710164988 A CN201710164988 A CN 201710164988A CN 106951634 B CN106951634 B CN 106951634B
- Authority
- CN
- China
- Prior art keywords
- model
- engine
- control system
- uncertainty
- tracking
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 239000011159 matrix material Substances 0.000 claims description 23
- 230000004044 response Effects 0.000 claims description 12
- 238000002485 combustion reaction Methods 0.000 claims description 11
- 239000000446 fuel Substances 0.000 claims description 6
- 239000000700 radioactive tracer Substances 0.000 abstract 1
- 238000004088 simulation Methods 0.000 description 11
- 230000000694 effects Effects 0.000 description 9
- 230000008569 process Effects 0.000 description 6
- 230000009471 action Effects 0.000 description 4
- 230000005764 inhibitory process Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000003416 augmentation Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013643 reference control Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
本发明公开了一种航空发动机鲁棒跟踪控制器设计方法,是面向航空发动机及其控制系统中存在不确定性和噪声等问题的Leitmann鲁棒跟踪控制方法,包括:建立航空发动机标称状态变量模型;建立含有不确定性和噪声的航空发动机控制系统状态变量模型;建立基于参考模型的航空发动机跟踪控制系统动态模型;Leitmann鲁棒跟踪控制器设计。提出“一致跟踪有界”和“一致最终跟踪有界”的跟踪控制系统性能,提出了基于Leitmann控制器的鲁棒跟踪方法,克服了航空发动机控制系统中存在的噪声和不确定性等问题,满足系统的跟踪要求,保证了被控系统在发动机工作环境、工作状态、动态性能变化时仍然具有良好鲁棒性。
Description
技术领域
本发明属于航空发动机控制技术领域,具体涉及一种航空发动机鲁棒跟踪控制器设计方法。
背景技术
近年来飞机性能指标不断提高,比如推重比、经济性和控制性能等,这对航空发动机也提出了更高的要求。为了获得更好的航空发动机控制性能,航空发动机控制系统控制量和被控量不断增加,系统也变得愈发复杂,这包括多控制输入多输出间的动态耦合,非线性变化的运行环境和运行状态(通常以转速表征)。考虑到控制器设计时,通常存在多个建模过程,如首先基于航空发动机气动热力学工作原理和试验数据建立非线性的部件级模型,随后根据该非线性部件级模型通过小扰动法、系统辨识等理论和方法建立用于控制器设计的线性模型。在这个过程中不可避免的存在建模误差。线化模型中参数的受到空气来流影响还有各个发动机之间存在的个体差异,航空发动机在运行过程中噪声额影响等。综上所及的非线性、建模误差、参数摄动、个体差异和噪声问题给航空发动机控制器设计者带来的挑战。
发动机控制系统除了面临上述来自环境和发动机自身的挑战外,随着飞机机动性要求提高,需要发动机机跟随飞行员操控提供更高动态、稳态性能的推力,意味着控制系统需要具有良好的动态指令跟踪特性
因此,针对上述航空发动机控制系统设计中面临的问题和需求,有必要发明一种新的航空发动机鲁棒跟踪控制方法,在克服非线性、建模误差、参数摄动、个体差异和噪声对控制系统影响的同时,保证被控系统具有良好的跟踪性能。
发明内容
发明目的:为了克服现有技术中存在的不足,本发明提供一种基于Leitmann方法的航空发动机鲁棒跟踪控制器设计方法,将系统中非线性、建模误差、参数摄动、个体差异视为不确定性。针对噪声、不确定性,设计鲁棒跟踪控制器,保证被控系统在飞行包线内良好的跟踪性能、抗干扰能力。
技术方案:为实现上述目的,本发明采用的技术方案为:
一种航空发动机鲁棒跟踪控制器设计方法,包括以下步骤:
步骤1)建立航空发动机标称状态变量模型;
步骤2)建立包含噪声和不确定性的航空发动机控制系统状态变量模型;
步骤3)建立基于参考模型的航空发动机跟踪控制系统动态模型;
步骤4)设计Leitmann鲁棒跟踪控制器。
进一步的,所述步骤1)中建立航空发动机状态变量模型的具体步骤如下:
步骤1.1),根据航空发动机气动热力学特性和典型部件特性数据建立发动机部件级模型,带加力双轴涡扇发动机模型主要部件包括进气道、风扇、压气机、燃烧室、高压涡轮、低压涡轮、混合室、加力燃烧室和尾喷管;
步骤1.2),根据所建发动机部件级模型,采用小扰动法和拟合法建立发动机控制系统状态方程:
yp(t)=Cxp(t)+Dup(t).
其中,xp(t)=[NLπT]T为状态量,NL和πT分别为风扇转速和涡轮落压比;up(t)=[WfbA8]T为控制量,Wfb和A8分别为主燃烧室供油量和尾喷管喉道面积;yp(t)=[NLπT]T为输出量。A,B,C,D是系统矩阵。
进一步的,所述步骤1.2)中建立发动机状态变量模型的具体步骤如下:
步骤1.2.1)采用小扰动法求出系数矩阵A,C的初始解;
步骤1.2.2)对发动机部件级模型作控制量up(t)阶跃得到其动态响应;
步骤1.2.3)根据非线性模型动态响应的稳态终值计算出矩阵B,D的各个元素。
进一步的,步骤2)中所述建立包含不确定性和噪声的航空发动机控制系统状态变量模型的具体步骤如下:
步骤2.1)建立包含噪声的发动机控制系统动态模型:
yp(t)=Cxp(t)+Dup(t)+D1wp(t).
其中,B1和D1是具有适当维数的噪声系数矩阵,wp(t)表示了系统模型中的噪声信号,满足:
||wp(t)||≤ωmax
其中,ωmax是噪声的范数上界;
步骤2.2)进一步考虑不确定性因素,建立含不确定性的发动机控制系统动态模型:
y(t)=Cxp(t)+Dup(t)+D1wp(t).
其中,ΔA是2×2维的不确定参数矩阵函数,表示了系统模型中的参数不确定性;假定所考虑的参数不确定性的表达形式为:
ΔA(xp(t),α(t),t)=α(t)·A
其中,标量α(t)是不确定性系数,且满足||α(t)||≤αmax,其中αmax是不确定性系数上界。
进一步的,步骤3)中所述建立基于参考模型的航空发动机跟踪控制系统动态模型中,参考模型的形式如下:
yss(t)=Cxss(t)+Duss(t).
其中xss(t)=[NLssπTss]T为状态量,NLss和πTss分别为参考模型的风扇转速和涡轮落压比;upss(t)=[Wfbss A8ss]T为控制量,Wfbss和A8ss分别为参考模型的主燃烧室供油量和尾喷管喉道面积;ypss(t)=[NLssπTss]T为输出量。
基于参考模型的航空发动机跟踪控制系统模型形式如下:
y(t)=Cx(t)+Du(t)+D1w(t).
其中,x(t)=xp(t)-xss(t),u(t)=up(t)-uss(t),y(t)=yp(t)-yss(t)。
进一步的,步骤4)中所述的H∞/Leitmann鲁棒控制器设计的具体步骤如下:
步骤4.1),针对基于参考模型的航空发动机跟踪控制系统动态模型所描述的被控对象,选取由p(t)和q(t)组成状态反馈控制器,其具体形式如下:
u(t)=q(t)+p(t)
步骤4.2),将基于参考模型的航空发动机跟踪控制系统动态模型所描述的被控对象,分为不含确定性的标称系统和不确定性两部分;针对标称系统,设计状态控制器q(t),具体形式如下:
q(t)=K1x(t)
则标称系统形式如下:
y(t)=Cx(t)+Du(t)+D1wp(t).
步骤4.3),定义变量ρ>0,矩阵X=XT>0和待求的矩阵W;
步骤4.4),利用Matlab中的LMI工具箱求解矩阵不等式:
X>0,
则是标称系统鲁棒跟踪控制器控制增益阵K1的具体形式如下:
K1=WX-1
步骤4.5),确定噪声的范数上界ωmax和不确定性上界αmax。
步骤4.6),定义
确定矩阵E(t)和G,满足等式:q
ΔA(t)=BE(t),
B1=BG.
定义标量v(t)和ρ(x,w,t),满足:
ν(t)=E(t)x(t)+Gw(t),
||ν(t)||=||Ex(t)+Gw(t)||
≤||E||||x(t)||+||G||||w(t)||
≤||E||max||x(t)||+||G||||wmax||=ρ(x,w,t).
则针对基于参考模型的航空发动机跟踪控制系统具体形式如下:
步骤4.7),选取常数ε,定义
μ(x,v,t)=2BTX-1x(t)ρ(x,v,t)
不确定性部分跟踪控制器为:
步骤4.8),H∞/Leitmann鲁棒跟踪控制器具体形式如下:
进一步的,步骤4.7)中所述的常数ε取0.5。ε越小对不确定性影响抑制效果越好。但是,过小的ε导致Leitmann控制器频繁切换,使得状态量波动,控制效果变差。针对发动机对象,综合考虑其动态特性,选取ε=0.5。
有益效果:本发明提供的一种航空发动机鲁棒跟踪控制器设计方法,采用以上技术方案与现有技术相比,具有以下技术效果:
(1)本发明是基于Leitmann控制方法并结合H∞控制理论,提出的Leitmann鲁棒跟踪控制器设计方法,该方法所设计控制器使得航空发动机控制系统存在噪声和不确定性时,系统状态跟踪动态指令,且跟踪误差一致有界、一致最终有界。
(2)本发明通过噪声和不确定性来综合考虑发动机及其控制系统的动态特性、工作环境、外部噪声及其建模过程中的非线性、建模误差、个体差异和参数摄动,今儿针对噪声和不确定性设计Leitmann鲁棒跟踪控制器,因而被控系统对于不确定性有良好的鲁棒性;
(3)本发明能够保证动态跟踪系统的动态和静态性能良好,适用于发动机控制系统。
附图说明
图1是采用Leitmann方法的航空发动机鲁棒跟踪控制器的发动机控制系统结构图,图中表明基于本发明的控制流程以及本发明内容中所涉及跟踪控制实现策略。
图2是航空发动机气路工作截面标识图,表明本发明所应用于的涡扇发动机结构,也即涡扇发动机气动热力学模型所涉及的发动机部件,截面标识数将作为发动机某些信号量的下标,标识该信号位置,如“A8”尾喷口喉道面积;
图3是某涡扇发动机Leitmann鲁棒控制仿真结果图,表明基于本发明的发动机控制系统中转速、压比控制的全数字仿真效果。
图4是某涡扇发动机H∞鲁棒控制仿真结果图,对比表明基于本发明的发动机控制系统具有更好的抑制不确定性能力。
图5是本发明控制器用于H=5km,Ma=0.5时,控制系统中转速、压比控制的全数字仿真效果。
图6是H∞鲁棒控制器用于H=5km,Ma=0.5时,控制系统中转速、压比控制的全数字仿真效果。
图7是本发明控制器以及H∞鲁棒控制器用于包线内4个不同工作点时,低压转子转速指令NL,r做±5%阶跃时,实际输出的波动范围大小,验证本发明控制器具有较好的鲁棒性。
图8是本发明控制器以及H∞鲁棒控制器用于包线内4个不同工作点时,涡轮落压比指令πT,r做±5%阶跃时,实际输出的波动范围大小,验证本发明控制器具有较好的鲁棒性。
具体实施方式
下面结合附图对本发明作更进一步的说明。
本发明公开了一种面向航空发动机及其控制系统中存在不确定性和噪声等问题的新的鲁棒跟踪控制器设计方法。该方法包括:建立航空发动机标称状态变量模型;建立含有不确定性和噪声的航空发动机控制系统状态变量模型;建立基于参考模型的航空发动机跟踪控制系统动态模型;Leitmann鲁棒跟踪控制器设计。本发明并创新性地提出了“一致跟踪有界”和“一致最终跟踪有界”的跟踪控制系统性能,提出了基于Leitmann控制器的鲁棒跟踪方法。本发明克服了航空发动机控制系统中存在的噪声和不确定性等问题,满足了系统的跟踪要求,使得系统满足“一致跟踪有界”和“一致最终跟踪有界”的性能指标,保证了被控系统在发动机工作环境、工作状态、动态性能变化时仍然具有良好鲁棒性。
一种航空发动机鲁棒跟踪控制器设计方法,包括以下步骤:
步骤1)建立航空发动机标称状态变量模型,具体步骤如下:
步骤1.1),根据航空发动机气动热力学特性和典型部件特性数据建立发动机部件级模型。带加力双轴涡扇发动机模型主要部件包括进气道、风扇、压气机、燃烧室、高压涡轮、低压涡轮、混合室、加力燃烧室和尾喷管;
步骤1.2),根据所建发动机部件级模型,采用小扰动法和拟合法建立发动机控制系统状态方程:
yp(t)=Cxp(t)+Dup(t).
其中,xp(t)=[NLπT]T为状态量,NL和πT分别为风扇转速和涡轮落压比;up(t)=[WfbA8]T为控制量,Wfb和A8分别为主燃烧室供油量和尾喷管喉道面积;yp(t)=[NLπT]T为输出量,A,B,C,D是系统矩阵。
其中步骤1.2)中建立发动机状态变量模型的具体步骤如下:
步骤1.2.1)采用小扰动法求出系数矩阵A,C的初始解;
步骤1.2.2)对发动机部件级模型作控制量up(t)阶跃得到其动态响应;
步骤1.2.3)根据非线性模型动态响应的稳态终值计算出矩阵B,D的各个元素。
步骤2)建立的包含噪声和不确定性的发动机控制系统动态状态变量模型,具体步骤如下:
步骤2.1)建立包含噪声的发动机控制系统动态模型:
yp(t)=Cxp(t)+Dup(t)+D1wp(t).
其中,B1和D1是具有适当维数的噪声系数矩阵,wp(t)表示了系统模型中的噪声信号,假定噪声的范数是有界且表示为:
||wp(t)||≤ωmax
其中,ωmax是噪声的范数上界。
步骤2.2)进一步考虑不确定性因素,建立含不确定性的发动机控制系统动态模型:
y(t)=Cxp(t)+Dup(t)+D1wp(t).
其中,ΔA是具有适当维数的不确定参数矩阵函数,表示了系统模型中的参数不确定性;假定所考虑的参数不确定性的表达形式为:
ΔA(xp(t),α(t),t)=α(t)·A
其中,标量α(t)是不确定性系数,且范数有界,满足||α(t)||≤αmax,其中αmax是不确定性系数上界。
步骤3)建立基于参考模型的航空发动机跟踪控制系统动态模型,具体步骤如下:
建立如下参考模型:
yss(t)=Cxss(t)+Duss(t).
其中xss(t)=[NLssπTss]T为状态量,NLss和πTss分别为参考模型的风扇转速和涡轮落压比;upss(t)=[Wfbss A8ss]T为控制量,Wfbss和A8ss分别为参考模型的主燃烧室供油量和尾喷管喉道面积;ypss(t)=[NLssπTss]T为输出量。
基于参考模型的航空发动机跟踪控制系统模型形式如下:
y(t)=Cx(t)+Du(t)+D1w(t).
其中,x(t)=xp(t)-xss(t),u(t)=up(t)-uss(t),y(t)=yp(t)-yss(t)。
步骤4)设计Leitmann鲁棒跟踪控制器,具体步骤如下:
步骤4.1),针对基于参考模型的航空发动机跟踪控制系统动态模型所描述的被控对象,选取由p(t)和q(t)组成状态反馈控制器,其具体形式如下:
u(t)=q(t)+p(t)
步骤4.2),将基于参考模型的航空发动机跟踪控制系统动态模型所描述的被控对象,分为不含确定性的标称系统和不确定性两部分。针对标称系统,设计状态控制器q(t),具体形式如下:
q(t)=K1x(t)
则标称系统形式如下:
y(t)=Cx(t)+Du(t)+D1wp(t).
步骤4.3),定义变量ρ>0,矩阵X=XT>0和W;
步骤4.4),利用Matlab中的LMI工具箱求解矩阵不等式:
X>0,
则是标称系统鲁棒跟踪控制器控制增益阵K1的具体形式如下:
K1=WX-1
步骤4.5),确定噪声的范数上界ωmax和不确定性上界αmax。
步骤4.6),定义
确定矩阵E(t)和G,满足等式:
ΔA(t)=BE(t),
B1=BG.
定义标量v(t)和ρ(x,w,t),满足:
ν(t)=E(t)x(t)+Gw(t),
||ν(t)||=||Ex(t)+Gw(t)||
≤||E||||x(t)||+||G||||w(t)||
≤||E||max||x(t)||+||G||||wmax||=ρ(x,w,t).
则针对基于参考模型的航空发动机跟踪控制系统具体形式如下:
步骤4.7),选取常数ε,定义
μ(x,v,t)=2BTX-1x(t)ρ(x,v,t)
不确定性部分跟踪控制器为:
步骤4.8),H∞/Leitmann鲁棒跟踪控制器具体形式如下:
实施例
本发明具体实施方式以某型涡扇发动机的多变量控制系统为例,图1为本发明的发动机控制系统结构图。基于涡扇发动机气动热力学的部件级模型,建立其状态变量模型,并以此描述被控涡扇发动机。利用参考模型获得基于给定输出指令的参考状态变量和参考控制变量,并与涡扇发动机状态变量模型共同构建跟踪控制系统模型。基于该模型,设计鲁棒控制器,获得控制变量传输给被控涡扇发动机,完成闭环控制。
本发明使用的某型涡扇发动机气路部件截面标识图如图2所示,主要部件包括进气道、风扇、压气机、燃烧室、高压涡轮、低压涡轮、混合室、加力燃烧室和尾喷管等。发动机在保险内某点、某工作稳定状态下的状态变量模型如下:
yp(t)=Cxp(t)+Dup(t).
本发明中采用小扰动法与拟合法相结合求取系数矩阵。该方法的原理是,首先采用小扰动法求出系数矩阵A,C的初始解;对发动机部件级模型作控制量阶跃得到其动态响应;根据非线性模型动态响应的稳态终值计算出矩阵B,D的各个元素。
结合下述实施例,在MATLAB环境下进行开展某涡扇发动机的跟踪控制系统全数字仿真,进一步验证本发明鲁棒跟踪控制器的有效性。
在发动机设计点H=0,Ma=0,Wf=2.48kg/s,A8=0.2597m2,首先利用小扰动法和拟合法结合求解在该工作点下的状态变量模型系数矩阵,各系数矩阵为:
选取不确定性系数表示为:
α(t)=0.2[1.8(-0.5+rand(1,1))+0.1sin(t)].
其中,rand(1,1)是0到1之间的随机数,sin(·)是正弦函数。那么不确定性上界为:
αmax=0.2.
选取噪声系数为:
求解线性矩阵不等式,可得到:
选取ε=0.5,那么不确定性部分控制器为:
其中,
因此,Leitmann鲁棒跟踪控制器为:
在该控制器作用下,对低压转子转速和涡轮落压比分别做±2%,±3%和±5%阶跃,转速、压比响应曲线如图3。为了验证不确定性部分控制器p(t)的控制效果,将本发明控制器中标称系统部分H∞控制器用于系统控制,同样对低压转子转速和涡轮落压比分别做±2%,±3%和±5%阶跃,转速、压比响应曲线如图4。通过对比可以发现,闭环系统应用本发明Leitmann鲁棒跟踪控制器相较于H∞控制器,转速和压比波动较小,受不确定性影响减小。因此,本发明航空发动机控制器能很好的抑制不确定性和噪声的影响。
本发明中,进一步在H=0,Ma=0.3、H=5km,Ma=0.5、H=10km,Ma=1.0进行仿真分析,验证Leitmann鲁棒跟踪控制器的鲁棒性,在此给出H=5km,Ma=0.5时的仿真结果,如图5所示。H=5km,Ma=0.5时发动机在H∞控制器作用下,仿真结果如图6。进一步对比观察本发明控制器与H∞鲁棒控制器对不确定性的抑制作用,统计两控制器在上述飞行包线内4个工作点,低压转子转速指令NL,r做±5%阶跃时,NL跟踪指令NL,r过程中的摄动量范围大小如图7所示。涡轮落压比指令πT,r做±5%阶跃时,πT跟踪指令πT,r的过程中摄动量范围大小如图8所示。其结果表明,本发明控制器在非设计点时依然具有良好的控制效果和抑制不确定性的能力,具有很好的鲁棒性。
本发明中状态变量模型是在发动机某稳态工作点对其气动热力学模型进行泰勒级数展开,忽略高阶无穷小项的影响,采用小扰动法和拟合法结合的方法求解系数矩阵而获得。考外界噪声扰动、建模误差、模型参数摄动、个体差异,在系统中引入了不确定性描述这些因素对发动机特性的影响;基于参考模型和状态变量模型建立跟踪控制系统;提出由标称系统的H∞控制器和针对不确定性的Leitmann控制器共同构成的Leitmann鲁棒跟踪控制器,保证跟踪控制系统“一致跟踪有界”以和“一致最终跟踪有界”。将该控制器应用于某型涡扇发动机转速、压比的控制系统实例中。在不确定性上界设为0.2,ε=0.5时,计算Leitmann鲁棒跟踪控制器。在此控制器作用下,令低压转子转速和涡轮落压比指令分别在包线内4个不同工作点做±2%,±3%和±5%阶跃,发动机低压转子转速、压比响应结果如图3-图8。同时,开展仅由H∞控制器作用的控制系统仿真,其系统响应结果同样绘于图3-图8中。由图可知,相较H∞控制器,本发明中Leitmann鲁棒控制器,在包线内各典型工作点处都能够更好地抑制噪声和不确定性对系统的影响,控制系统具有更好的鲁棒性。因此,本发明中基于参数不确定性的建模方法能够描述发动机及其控制系统中的参数摄动、非线性等特性,鲁棒控制方法保证了(1)航空发动机控制系统的性能基线,即对于存在噪声的发动机标称系统渐进稳定,且对噪声具有H∞抑制能力;(2)航空发动机不确定动态跟踪控制系统的一致有界和一致最终有界。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (3)
1.一种航空发动机鲁棒跟踪控制器设计方法,其特征在于:包括以下步骤:
步骤1)建立航空发动机标称状态变量模型:
步骤1.1),根据航空发动机气动热力学特性和典型部件特性数据建立发动机部件级模型,带加力双轴涡扇发动机模型主要部件包括进气道、风扇、压气机、燃烧室、高压涡轮、低压涡轮、混合室、加力燃烧室和尾喷管;
步骤1.2),根据所建发动机部件级模型,采用小扰动法和拟合法建立发动机控制系统状态方程:
yp(t)=Cxp(t)+Dup(t).
其中,xp(t)=[NL πT]T为状态量,NL和πT分别为风扇转速和涡轮落压比;up(t)=[WfbA8]T为控制量,Wfb和A8分别为主燃烧室供油量和尾喷管喉道面积;yp(t)=[NL πT]T为输出量,A,B,C,D是系统矩阵;
步骤2)建立包含噪声和不确定性的航空发动机控制系统状态变量模型:
步骤2.1)建立包含噪声的发动机控制系统动态模型:
yp(t)=Cxp(t)+Dup(t)+D1wp(t).
其中,B1和D1是具有适当维数的噪声系数矩阵,wp(t)表示了系统模型中的噪声信号,满足:
||wp(t)||≤ωmax
其中,ωmax是噪声的范数上界;
步骤2.2)进一步考虑不确定性因素,建立含不确定性的发动机控制系统动态模型:
y(t)=Cxp(t)+Dup(t)+D1wp(t).
其中,ΔA是2×2维的不确定参数矩阵函数,表示了系统模型中的参数不确定性;假定所考虑的参数不确定性的表达形式为:
ΔA(xp(t),α(t),t)=α(t)·A
其中,标量α(t)是不确定性系数,且满足||α(t)||≤αmax,其中αmax是不确定性系数上界;
步骤3)建立基于参考模型的航空发动机跟踪控制系统动态模型;
参考模型的形式如下:
yss(t)=Cxss(t)+Duss(t).
其中xss(t)=[NLss πTss]T为状态量,NLss和πTss分别为参考模型的风扇转速和涡轮落压比;uss(t)=[Wfbss A8ss]T为控制量,Wfbss和A8ss分别为参考模型的主燃烧室供油量和尾喷管喉道面积;yss(t)=[NLss πTss]T为输出量;
基于参考模型的航空发动机跟踪控制系统模型形式如下:
y(t)=Cx(t)+Du(t)+D1w(t).
其中,x(t)=xp(t)-xss(t),u(t)=up(t)-uss(t),y(t)=yp(t)-yss(t);
步骤4)设计H∞/Leitmann鲁棒跟踪控制器:
步骤4.1),针对基于参考模型的航空发动机跟踪控制系统动态模型所描述的被控对象,选取由p(t)和q(t)组成状态反馈控制器,其具体形式如下:
u(t)=q(t)+p(t)
步骤4.2),将基于参考模型的航空发动机跟踪控制系统动态模型所描述的被控对象,分为不含确定性的标称系统和不确定性两部分;针对标称系统,设计状态控制器q(t),具体形式如下:
q(t)=K1x(t)
则标称系统形式如下:
y(t)=Cx(t)+Du(t)+D1wp(t).
步骤4.3),定义变量ρ>0,矩阵X=XT>0和待求的矩阵W;
步骤4.4),利用Matlab中的LMI工具箱求解矩阵不等式:
则是标称系统鲁棒跟踪控制器控制增益阵K1的具体形式如下:
K1=WX-1
步骤4.5),确定噪声的范数上界ωmax和不确定性上界αmax;
步骤4.6),定义
确定矩阵E(t)和G,满足等式:
ΔA(t)=BE(t),
B1=BG.
定义标量v(t)和ρ(x,w,t),满足:
ν(t)=E(t)x(t)+Gw(t),
||ν(t)||=||Ex(t)+Gw(t)||
≤||E||||x(t)||+||G|| ||w(t)||
≤||E||max||x(t)||+||G|| ||wmax||=ρ(x,w,t).
则针对基于参考模型的航空发动机跟踪控制系统具体形式如下:
步骤4.7),选取常数ε,定义
μ(x,v,t)=2BTX-1x(t)ρ(x,v,t)
不确定性部分跟踪控制器为:
步骤4.8),H∞/Leitmann鲁棒跟踪控制器具体形式如下:
2.根据权利要求1所述的一种航空发动机鲁棒跟踪控制器设计方法,其特征在于:所述步骤1.2)中建立发动机部件级模型的具体步骤如下:
步骤1.2.1)采用小扰动法求出系统矩阵A,C的初始解;
步骤1.2.2)对发动机部件级模型作控制量up(t)阶跃得到其动态响应;
步骤1.2.3)根据非线性模型动态响应的稳态终值计算出系统矩阵B,D的各个元素。
3.根据权利要求1所述的一种航空发动机鲁棒跟踪控制器设计方法,其特征在于,步骤4.7)中所述的常数ε取0.5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710164988.9A CN106951634B (zh) | 2017-03-20 | 2017-03-20 | 一种航空发动机鲁棒跟踪控制器设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710164988.9A CN106951634B (zh) | 2017-03-20 | 2017-03-20 | 一种航空发动机鲁棒跟踪控制器设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106951634A CN106951634A (zh) | 2017-07-14 |
CN106951634B true CN106951634B (zh) | 2019-10-11 |
Family
ID=59472123
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710164988.9A Expired - Fee Related CN106951634B (zh) | 2017-03-20 | 2017-03-20 | 一种航空发动机鲁棒跟踪控制器设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106951634B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108828947B (zh) * | 2018-07-13 | 2021-06-01 | 南京航空航天大学 | 一种航空发动机含时滞的不确定性模糊动态模型建模方法 |
CN109031951B (zh) * | 2018-07-16 | 2020-02-18 | 南京航空航天大学 | 基于精确偏导数的航空发动机状态变量模型在线建立方法 |
CN110907120B (zh) * | 2018-09-17 | 2021-05-11 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种引射式短舱出口落压比实时反馈控制方法 |
CN109446605B (zh) * | 2018-10-16 | 2020-06-05 | 南京航空航天大学 | 涡轴发动机非线性动态逆控制方法及装置 |
CN110362960B (zh) * | 2019-08-09 | 2022-11-11 | 哈尔滨工业大学 | 基于多胞折合平衡流形展开模型的航空发动机系统辨识方法 |
CN111046568B (zh) * | 2019-12-18 | 2022-09-09 | 哈尔滨工程大学 | 一种基于航空发动机及尾喷管机电液多系统联合仿真的控制参数优化方法 |
CN111273554B (zh) * | 2020-04-04 | 2022-06-07 | 西北工业大学 | 航空发动机最大推力状态降保守性二自由度h∞控制器 |
CN111443607B (zh) * | 2020-04-04 | 2022-05-31 | 西北工业大学 | 基于健康退化的航空发动机降保守性鲁棒增益调度控制器 |
CN111459028B (zh) * | 2020-04-04 | 2022-06-07 | 西北工业大学 | 航空发动机最大推力状态降保守性二自由度μ控制器 |
CN111523184B (zh) * | 2020-04-21 | 2022-01-04 | 南京航空航天大学 | 一种航空发动机控制系统的鲁棒性分析方法及存储介质 |
CN114415506B (zh) * | 2022-01-07 | 2023-08-04 | 大连理工大学 | 航空发动机双模跟踪预测控制系统设计方法 |
CN114779637A (zh) * | 2022-04-18 | 2022-07-22 | 南京航空航天大学 | 基于变换函数的航空发动机输入饱和鲁棒控制方法 |
CN114791702A (zh) * | 2022-05-11 | 2022-07-26 | 南京航空航天大学 | 基于线性变参数模型的航空发动机滑模控制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004102938A (ja) * | 2002-09-12 | 2004-04-02 | Ricoh Co Ltd | 機構制御系設計方法、機構制御系設計システム、機構制御系設計プログラム及びそのプログラムが記録された記録媒体 |
CN102855349A (zh) * | 2012-08-06 | 2013-01-02 | 南京航空航天大学 | 航空发动机气路故障诊断的快速原型设计方法及平台 |
CN105093931A (zh) * | 2015-06-08 | 2015-11-25 | 南京航空航天大学 | 一种航空发动机非线性系统控制器设计方法 |
-
2017
- 2017-03-20 CN CN201710164988.9A patent/CN106951634B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004102938A (ja) * | 2002-09-12 | 2004-04-02 | Ricoh Co Ltd | 機構制御系設計方法、機構制御系設計システム、機構制御系設計プログラム及びそのプログラムが記録された記録媒体 |
CN102855349A (zh) * | 2012-08-06 | 2013-01-02 | 南京航空航天大学 | 航空发动机气路故障诊断的快速原型设计方法及平台 |
CN105093931A (zh) * | 2015-06-08 | 2015-11-25 | 南京航空航天大学 | 一种航空发动机非线性系统控制器设计方法 |
Non-Patent Citations (3)
Title |
---|
《基于LMI的航空发动机鲁棒H∞控制器设计》;谢光华;《航空学报》;20000331;第175-178页 * |
《带有网络时延的航空发动机分布式系统H∞鲁棒控制》;李勇;《伺服驱动系统及其应用》;20131231;第39-41,44页 * |
《航空发动机H∞状态反馈控制的鲁棒稳定界》;王曦;《航空动力学报》;20011030;第408-410页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106951634A (zh) | 2017-07-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106951634B (zh) | 一种航空发动机鲁棒跟踪控制器设计方法 | |
CN106647253B (zh) | 航空发动机分布式控制系统多性能鲁棒跟踪控制方法 | |
CN108762089B (zh) | 一种基于模型预测的航空发动机在线优化及多变量控制设计方法 | |
Garg | Aircraft turbine engine control research at NASA Glenn research center | |
Kulikov et al. | Dynamic modelling of gas turbines: identification, simulation, condition monitoring and optimal control | |
Chapman et al. | Propulsion system simulation using the toolbox for the modeling and analysis of thermodynamic systems (t mats) | |
CN108828947B (zh) | 一种航空发动机含时滞的不确定性模糊动态模型建模方法 | |
EP3073102B1 (en) | Control scheme using variable area turbine and exhaust nozzle to reduce drag | |
CN111680357B (zh) | 一种变循环发动机机载实时模型的部件级无迭代构建方法 | |
Liu et al. | Design for aircraft engine multi-objective controllers with switching characteristics | |
Mohammadian et al. | Simulation of startup operation of an industrial twin-shaft gas turbine based on geometry and control logic | |
Connolly et al. | Advanced control considerations for turbofan engine design | |
Zhang et al. | General design method of control law for adaptive cycle engine mode transition | |
CN112257256A (zh) | 一种基于稳态数据的发动机简化动态模型设计方法 | |
CN113359484B (zh) | 基于半交替优化的航空发动机模型预测控制方法及装置 | |
Albert et al. | Aerodynamic design optimization of nacelle and intake | |
Smith et al. | Optimizing aircraft performance with adaptive, integrated flight/propulsion control | |
Csank et al. | A turbine based combined cycle engine inlet model and mode transition simulation based on HiTECC tool | |
Yarlagadda | Performance analysis of J85 turbojet engine matching thrust with reduced inlet pressure to the compressor | |
Zhen et al. | Accelerated formulation of optimal control law for adaptive cycle engines: A novel design methodology | |
Tudosie et al. | Aircraft Gas-Turbine Engine with Coolant Injection for Effective Thrust Augmentation as Controlled Object | |
Rotaru et al. | Aircraft engine mathematical model–linear system approach | |
Kim | Application of machine learning and its effectiveness in performance model adaptation for a turbofan engine | |
Xia et al. | Research on power regulation schedule control system for turboprop engine | |
Lutambo et al. | Turbofan engine modelling and control design using linear quadratic regulator (LQR) |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191011 Termination date: 20210320 |
|
CF01 | Termination of patent right due to non-payment of annual fee |