CN112327611A - 变循环发动机增益调度二自由度h∞控制器 - Google Patents

变循环发动机增益调度二自由度h∞控制器 Download PDF

Info

Publication number
CN112327611A
CN112327611A CN202010542137.5A CN202010542137A CN112327611A CN 112327611 A CN112327611 A CN 112327611A CN 202010542137 A CN202010542137 A CN 202010542137A CN 112327611 A CN112327611 A CN 112327611A
Authority
CN
China
Prior art keywords
degree
freedom
engine
infinity
scheduling
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
CN202010542137.5A
Other languages
English (en)
Other versions
CN112327611B (zh
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 CN202010542137.5A priority Critical patent/CN112327611B/zh
Publication of CN112327611A publication Critical patent/CN112327611A/zh
Application granted granted Critical
Publication of CN112327611B publication Critical patent/CN112327611B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G05B13/045Adaptive 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 using a perturbation signal
    • 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

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

本发明提出一种变循环发动机增益调度二自由度H∞控制器。二自由度H∞控制器组调度计算模块产生控制输入向量u并输出给发动机本体,由传感器输出调度参数α和模式选择活门打开程度msv至二自由度H∞控制器组调度计算模块;二自由度H∞控制器组调度计算模块根据输入的调度参数α和模式选择活门打开程度msv,利用内部设计的若干二自由度H∞控制器计算得到适应的二自由度H∞控制器,该二自由度H∞控制器根据参考输入r和测量参数y的差值e产生控制输入向量u。本发明能够在全飞行包线内对变循环发动机在不同的工作模式下进行良好控制,可以同时保证系统的鲁棒稳定性和鲁棒性能,保证控制系统的稳定性并充分发挥变循环发动机的性能。

Description

变循环发动机增益调度二自由度H∞控制器
技术领域
本发明涉及变循环航空发动机控制技术领域,尤其涉及一种变循环发动机增益调度二自由度H∞控制器。
背景技术
航空发动机是一个复杂的非线性动力学系统,在范围宽广的飞行包线内工作时,发动机的工作状态随着外部条件和飞行条件的变化而不断变化。针对航空发动机的强非线性和模型的不确定性,现有技术中有提出鲁棒增益调度控制方法,将发动机划分为一系列工作点,并在每一个工作点设计鲁棒控制器,最终采用增益调度的方法选择合适的鲁棒控制器对发动机进行控制。
上述航空发动机鲁棒增益调度控制方法可以对航空发动机进行控制。然而,现代战争要求先进战斗机具备长航程亚声速巡航的能力,同时在作战时又要具备快速反应能力,未来航空发动机将向长巡航里程、高推重比、宽工作范围三个方向不断发展。通过研究常规发动机速度特性,研究者发现超声速状态下涡喷发动机具有较高的单位推力和较低的单位燃油消耗率而亚声速状态下大涵道比涡扇发动机具有较低的单位燃油消耗率。考虑现代战争对战斗机推进系统的性能要求,涡扇发动机更加适合亚声速飞行,而涡喷发动机更适合超声速飞行。因此,便有了性能更好的变循环发动机。在发动机不同的工作状态下,通过采用调节特征部件的几何形状、物理位置或尺寸大小等不同的技术手段,将涡扇和涡喷两种不同的航空发动机的性能优势集中一体,从而保证变循环发动机在亚声速巡航状态下以涡扇发动机类似构型工作,从而获得较高的经济性,在超声速作战状态下以涡喷发动机类似构型工作,从而获得持续可靠的高单位推力,达到了将涡扇、涡喷发动机的性能优势融为一体的目的,使变循环发动机在发动机工作全过程中均具有优良的性能。变循环发动机在不同的工作模式下,可以调整自身热力循环状态,除了需要具备特殊的特征部件外,还需要与之适应、性能优良的控制系统。相比传统航空发动机,变循环发动机工作包线宽广、工作模式复杂、可调部件众多,控制系统中控制变量更多、变量间耦合性更强。传统的鲁棒增益调度控制无法适用于变循环发动机,另外,针对传统单自由度控制器无法同时兼顾航空发动机控制系统的鲁棒稳定性和鲁棒性能,二自由度H∞控制器设计方法在传统H∞控制器的基础上加入前置滤波器和反馈控制器,通过调整反馈控制器使扰动抑制的能力达到最佳,在此基础上调整前置滤波器使系统的指令跟踪能力达到最佳。因此,研究变循环发动机增益调度二自由度H∞控制方法具有重要意义。
发明内容
为解决现有技术存在的问题,本发明提出一种变循环发动机增益调度二自由度H∞控制器,能够在全飞行包线内对变循环发动机在不同的工作模式下进行良好控制,同时兼顾变循环发动机控制系统的鲁棒稳定性和鲁棒性能,保证控制系统的稳定性并充分发挥变循环发动机的性能。
本发明的技术方案为:
所述一种变循环发动机增益调度二自由度H∞控制器,其特征在于:包括二自由度H∞控制器组调度计算模块;
其中二自由度H∞控制器组调度计算模块与变循环发动机本体以及变循环发动机上的若干传感器组成调度控制回路;由传感器输出调度参数α和模式选择活门打开程度msv至二自由度H∞控制器组调度计算模块;
所述二自由度H∞控制器组调度计算模块产生控制输入向量u并输出给变循环发动机本体,传感器得到变循环发动机测量参数y;
所述二自由度H∞控制器组调度计算模块内设计有若干二自由度H∞控制器,所述二自由度H∞控制器包含前置滤波器和反馈控制器,对若干不确定性发动机模型利用H∞回路成形方法分别设计得到;
所述线性不确定性发动机模型是对变循环发动机不同模式选择活门MSV打开程度、不同调度参数下的非线性发动机模型进行线性化后再加入摄动块得到的;
所述二自由度H∞控制器组调度计算模块根据输入的调度参数α和模式选择活门打开程度msv,利用内部设计的若干二自由度H∞控制器计算得到适应的二自由度H∞控制器,该二自由度H∞控制器根据参考输入r和测量参数y的差值e产生控制输入向量u。
进一步的,所述二自由度H∞控制器组调度计算模块内设计若干二自由度H∞控制器的过程为:在全飞行包线内根据调度参数α和模式选择活门打开程度msv,选取q*J个工作点对非线性发动机模型进行线性化得到q*J个线性化模型再加入摄动块得到q*J个线性不确定性发动机模型,并对这q*J个线性不确定性发动机模型分别设计相应的二自由度H∞控制器从而组成二自由度H∞控制器组。
所述非线性发动机模型为:
Figure BDA0002539278950000031
y=g(x,u,msv)
其中
Figure BDA0002539278950000032
为控制输入向量,
Figure BDA0002539278950000033
为状态向量,
Figure BDA0002539278950000034
为输出向量,
Figure BDA0002539278950000035
为模式选择活门MSV打开程度,f(·)为表示系统动态的n维可微非线性向量函数,g(·)为产生系统输出的m维可微非线性向量函数。
所述线性化过程为:从模式选择活门MSV完全闭合msv0到完全打开msvJ等间距的选取模式选择活门MSV打开程度msvj,j=1,2,...,J,在每一个固定的msvj,j=1,2,...,J,选择一组调度参数值αi,i=1,2,...,q,代表系统的动态范围,并将飞行包线划分为几个子区间,并将这些点作为工作点,在工作点,采用小扰动法或拟合法对发动机非线性模型在稳态参考点处进行线性化,得到
Figure BDA0002539278950000036
其中Δ表示该参数的变化量,系数矩阵可由下式得到:
Figure BDA0002539278950000037
Figure BDA0002539278950000038
这些系数在发动机不同的工作状态具有不同的值。
所述线性不确定性发动机模型得到方法为:得到的各个工作点的线性化模型为标称模型,线性状态空间方程
Figure BDA0002539278950000041
也可以表示为
Figure BDA0002539278950000042
实际的发动机和标称模型之间的误差可以表示为一个摄动块Δ。在标称模型加入摄动块建立发动机不确定模型
Figure BDA0002539278950000043
进一步的,所述二自由度H∞控制器组调度计算模块根据输入的调度参数α和模式选择活门打开程度msv插值得到的适应的二自由度H∞控制器。
进一步的,所述二自由度H∞控制器组调度计算模块根据变循环发动机当前的调度参数α和模式选择活门打开程度msv选择相邻的四个设定工作点xi,j、xi,j+1、xi+1,j和xi+1,j+1(xi,j表示调度参数为αi,模式选择活门打开程度为msvj时的工作点),并获取四个设定工作点对应的线性控制器Ki,j、Ki,j+1、Ki+1,j和Ki+1,j+1(i=1,2,...,q,j=1,2,...,J)。
根据公式
Figure BDA0002539278950000044
计算得到变循环发动机当前适应的二自由度H∞控制器K(α,msv)。
进一步的,所述调度参数α包括变循环发动机的风扇转速或者压气机转速。
进一步的,所述测量参数包括进气道出口、风扇出口、压气机出口、高压涡轮后、低压涡轮后的温度和压力,风扇转速和压气机转速。
有益效果
与现有技术相比较,本发明的变循环发动机增益调度二自由度H∞控制器利用传统鲁棒增益调度控制器中固有的调度参数,通过新增模式选择活门打开程度msv这一新的调度参数,并对二自由度H∞控制器组调度计算模块进行了改进,新增了变循环发动机不同工作模式下的多组二自由度H∞控制器。结合传统的调度参数,实现变循环发动机的增益调度控制,能够在全飞行包线内对变循环发动机在不同的工作模式下进行良好控制,可以同时保证系统的鲁棒稳定性和鲁棒性能,保证控制系统的稳定性并充分发挥变循环发动机的性能。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1是本发明变循环发动机增益调度二自由度H∞控制器的结构简图;
图2是本发明变循环发动机结构示意图;
图3是本发明变循环发动机可调参数示意图;
图4是本发明变循环发动机可调部件示意图;
图5是本发明变循环发动机双外涵模式气流分布图;
图6是本发明变循环发动机单外涵模式气流分布图;
图7是本发明变循环发动机非线性模型线性化示意图;
图8是二自由度控制器的闭环系统结构图;
图9是二自由度控制器的闭环系统标准结构图;
图10是带有干扰和噪声的闭环结构图;
图11是指定L的期望奇异值示意图。
具体实施方式
航空发动机是一个复杂的非线性动力学系统,在范围宽广的飞行包线内工作时,发动机的工作状态随着外部条件和飞行条件的变化而不断变化,具有强非线性和模型的不确定性。针对发动机的控制,使用较多的是鲁棒增益调度控制。然而,相比传统航空发动机,变循环发动机工作包线宽广、工作模式复杂、可调部件众多,控制系统中控制变量更多、变量间耦合性更强。传统的鲁棒增益调度控制无法适用于变循环发动机,另外,传统单自由度控制器无法同时兼顾航空发动机控制系统的鲁棒稳定性和鲁棒性能。针对这一问题,下面给出本发明的分析研究过程。
1、变循环发动机工作原理
本发明以带有核心驱动风扇级(CDFS)的双外涵变循环发动机为主要研究对象,其主要结构如图2所示,包含的主要部件有进气道、风扇、核心驱动风扇级、高压压气机、燃烧室、高压涡轮、低压涡轮、混合室、加力燃烧室、尾喷管。相比普通双轴涡扇发动机,其显著的结构特点是在风扇和高压压气机之间增加了CDFS,同时在风扇和CDFS后分别设置副外涵和主外涵。在变循环发动机的不同工作状态下,通过改变CDFS的导叶角度,可以大幅度调节发动机外涵道和核心机的空气流量,进而调节发动机内外涵空气流量、涵道比和增压比等循环参数,使发动机的热力循环调节更加灵活。
相比普通双轴涡扇发动机,变循环发动机具有更多的可调部件。带有CDFS部件的变循环发动机主要有8个可调部件,具体如图3所示,可调部件示意图如图4所示。
与传统发动机相比,变循环发动机的性能优势主要体现在由于其可调部件增加,通过改变可调部件参数,调节发动机在工作过程中的气动热力循环,在保证推力基本不变时显著降低单位燃油消耗率,大大提高发动机的经济效益,同时可调部件的增加,使控制系统的调节过程更加灵活,风扇、压气机等部件的稳定裕度将大大提高。
变循环发动机具有单/双外涵两种典型的工作模式,它们通过模式选择活门MSV、FVABI、RVABI等可变阀门实现切换。当MSV完全打开时,气流经风扇后被分为两部分,一股气流流入副外涵,这部分气流最终在主外涵出口截面与主外涵气流有效掺混,流入总外涵。另一股气流流入CDFS,这股气流部分经RVABI被引入总外涵,其余气流将流入核心机。由于尾端涵道和RVABI的存在,总外涵气流在出口处会分为两部分,一股气流直接通过尾端涵道流入尾喷管,另一股气流则会进入混合室,与通过核心机的气流进行掺混后经加力燃烧室燃烧后,流入尾喷管,具体气流分布如图5所示。上述工作过程中,主外涵和副外涵均有气流通过,故被命名为双外涵模式。
当模式选择活门MSV完全关闭时,流经风扇的气流全部流入CDFS,风扇以压气机模式工作,副外涵不再有气流通过,这一过程被命名为单外涵工作模式,其具体气流分布如图6所示。
变循环发动机在不同的工作模式下进行切换时,内部热力循环状态会随之发生改变。为保证发动机能够持续保持稳定可靠地工作,平稳实现单双外涵模式转换,在模式切换过程中应该满足以下基本条件:
(1)风扇进口流量基本保持不变;
(2)风扇增压比基本保持不变;
(3)核心驱动风扇级的增压比随切换过程平稳变化;
(4)涵道比随MSV位移的改变而平稳变化;
(5)保证回流裕度始终大于0,即不存在气流绕CDFS的倒流;
(6)避免发生持续的超温、超转现象,避免喘振现象。
为了满足上述条件,在调节MSV位移时,应该配合调节其他可调部件参数,模式选择活门MSV打开程度可以表征变循环发动机的工作模式。目前已证实切实可行的模式切换调节策略是:在单外涵至双外涵的模式切换过程中,通过调节MSV位移,增大副外涵进口截面面积,为避免风扇压比的大幅度降低,需要配合减小CDFS进口导流叶片角度αi,同时减小可调涡轮导向器角度αt。由双外涵至单外涵的模式切换过程,调节策略相反。变循环发动机在不同工作模式工作时,为获得理想的涵道比同时保证气流不发生喘振或其他非正常工作状态,需要调节CDFS导流叶片角度αi以改变内涵空气流量,使之与发动机工作状态匹配。
2、变循环发动机增益调度控制设计
增益调度控制的实质是设计一组线性化的控制器,然后将它们有规律地组合起来,从而能够控制非线性系统。变循环发动机增益调度控制的基本原理是选择一系列的模式选择活门MSV打开程度,分别获得不同设定工作点下的发动机线性化模型并分别设计对应的二自由度H∞控制器得到图1中的二自由度H∞控制器组。
非线性模型是得到线性化模型的基础,基于部件法,建立变循环发动机非线性模型
Figure BDA0002539278950000071
y=g(x,u,msv)
其中
Figure BDA0002539278950000072
为控制输入向量,
Figure BDA0002539278950000073
为状态向量,
Figure BDA0002539278950000074
为输出向量,
Figure BDA0002539278950000075
为模式选择活门MSV打开程度,f(·)为表示系统动态的n维可微非线性向量函数,g(·)为产生系统输出的m维可微非线性向量函数。
请参阅图7,从模式选择活门MSV完全闭合msv0到完全打开msvJ等间距的选取模式选择活门MSV打开程度msvj,j=1,2,...,J,在每一个固定的msvj,j=1,2,...,J,选择一组调度参数值αi,i=1,2,...,q,代表系统的动态范围,并将飞行包线划分为几个子区间,并将这些点作为工作点,在工作点,采用小扰动法或拟合法对发动机非线性模型在稳态参考点处进行线性化,得到
Figure BDA0002539278950000081
其中Δ表示该参数的变化量,系数矩阵可由下式得到:
Figure BDA0002539278950000082
Figure BDA0002539278950000083
这些系数在发动机不同的工作状态具有不同的值。
请参阅图7,上、下实线分别表示模式选择活门MSV完全打开和模式选择活门MSV完全闭合时发动机的非线性模型。一系列的黑色小圆点表示发动机不同的工作点,在每一个工作点进行线性化得到线性模型。针对所有的线性模型,分别设计一系列二自由度H∞控制器得到图1中的二自由度H∞控制器组。然后,在选定的工作点之间,对控制器增益进行线性内插,使得对于所有的固定参数值,闭环系统都稳定且具有良好的性能。参数α是调度参数,这里可以定义为变循环发动机的风扇转速或者压气机转速,可以实时测量。控制系统的另一个调度变量是改变发动机工作模式的模式选择活门MSV的打开程度msv。工作原理是图1中的二自由度H∞控制器组调度计算模块根据调度参数α和模式选择活门打开程度msv进行线性插值获得相应的二自由度H∞控制器来控制系统。
3、不确定模型的二自由度H∞控制器设计
任何实际系统都不可避免地存在不确定性,它可以分为两类:扰动信号和模型不确定性。扰动信号包括干扰、噪声等。模型的不确定性代表了数学模型与实际对象之间的差异。
模型不确定性可能有几个原因:线性模型中总有一些参数是有误差的;线性模型中的参数可能由于非线性或工作条件的变化而变化;建模时人为的简化;由于磨损等因素发动机性能的退化。
不确定性可能会对控制系统的稳定性和性能产生不利影响。
上面得到的各个工作点的线性化模型为标称模型,线性状态空间方程
Figure BDA0002539278950000091
也可以表示为
Figure BDA0002539278950000092
实际的发动机和标称模型之间的误差可以表示为一个摄动块Δ。在标称模型加入摄动块建立发动机不确定模型
Figure BDA0002539278950000093
最后根据不确定模型利用传统的二自由度H∞控制器设计方法设计二自由度H∞控制器。
4、二自由度H∞控制器设计
具有二自由度控制器的闭环系统的框图如图8所示。该系统有一个参考输入(r),输出干扰(d)和两个输出误差(z1)和(z2)。系统M0是闭环系统应该匹配的理想模型。在这种结构中,除了内部稳定性要求之外,两个信号e和u将被最小化。信号e显示系统输出和参考模型输出之间的差异。u是控制信号,也与摄动中的鲁棒稳定性有关。在图8中,包含两个加权函数以反映这两个惩罚信号之间的特征。
所谓二自由度控制,就是对使目标值跟踪特性为最优的参数和使外扰抑制特性为最优的参数分别独立进行整定,使用反馈控制器(Ky)来实现内部稳定性、鲁棒稳定性和干扰抑制等,并且在前馈路径上设计另一个控制器(Kr)以满足跟踪需求,最小化整个系统的输出与参考模型M的输出之间的差异,使两特性同时达到最优。
图8的结构可以通过定义w=r,
Figure BDA0002539278950000094
来重新排列为图9的标准结构。控制器K由用于干扰衰减的反馈控制器Ky和前置滤波器Kr组成,以实现期望的闭环性能,并表示为
K=[Kr Ky]
多变量传递函数的
Figure BDA0002539278950000101
回路整形设计的思想是通过设置前置或后置补偿器对开环对象进行扩展,以使开环频率响应的奇异值有期望的形状。如图10中所示的结构,系统G和控制器K的互连是由参考命令r、输入干扰di、输出干扰do和测量噪声n驱动的,y为要控制的输出,u为控制信号。
对于输入灵敏度函数Si=(I+KG)-1、输出灵敏度函数So=(I+GK)-1和输出互补灵敏度函数To=GK(I+GK)-1,有如下关系:
Figure BDA0002539278950000102
这些关系决定了几个闭环目标:
1.对于输入干扰的衰减,使
Figure BDA0002539278950000103
很小。
2.对于输出干扰衰,使
Figure BDA0002539278950000104
很小。
3.用于噪声抑制,使
Figure BDA0002539278950000105
很小。
4.为好的参考跟踪,使
Figure BDA0002539278950000106
在经典的环路整形中,所整形的是开环传递函数L=GK幅值的大小,如图11所示。
采用
Figure BDA0002539278950000107
回路整形设计方法设计的控制器阶数高,限制了控制器的实时性,难以实现。利用绝对误差逼近法对设计的鲁棒控制器进行适当的降阶,得到降阶控制器Kr(s),即使下式最小
||K(s)-Kr(s)||
5、控制器的插值
这部分说明了图1中的二自由度H∞控制器组调度计算模块通过调度参数和健康参数调度线性插值获得相应的二自由度H∞控制器的调度计算原理。
根据变循环发动机当前的调度参数α和模式选择活门打开程度msv选择相邻的四个设定工作点xi,j、xi,j+1、xi+1,j和xi+1,j+1(xi,j表示调度参数为αi,模式选择活门打开程度为msvj时的工作点),并获取四个设定工作点对应的线性控制器Ki,j、Ki,j+1、Ki+1,j和Ki+1,j+1(i=1,2,...,q,j=1,2,...,J)。根据公式
Figure BDA0002539278950000111
计算得到变循环发动机当前适应的二自由度H∞控制器K(α,msv),并对发动机进行有效控制。
基于上述过程,下面给出本实施例中提出的一种变循环发动机增益调度二自由度H∞控制器,如图1所示,主要包括二自由度H∞控制器组调度计算模块。
其中二自由度H∞控制器组调度计算模块与变循环发动机本体以及变循环发动机上的若干传感器组成多维调度控制回路10。
所述二自由度H∞控制器组调度计算模块产生控制输入向量u并输出给变循环发动机本体,传感器得到变循环发动机测量参数y;传感器同时得到发动机调度参数α和模式选择活门打开程度msv,并输出到二自由度H∞控制器组调度计算模块。
所述二自由度H∞控制器组调度计算模块内设计有若干二自由度H∞控制器,所述二自由度H∞控制器包含前置滤波器和反馈控制器,对若干不确定性发动机模型利用H∞回路成形方法分别设计得到。
所述线性不确定性发动机模型是对变循环发动机不同模式选择活门MSV打开程度、不同调度参数下的非线性发动机模型进行线性化后再加入摄动块得到的。
优选的一种具体实现方式,可以通过以下过程得到设计若干二自由度H∞控制器:在全飞行包线内根据调度参数α和模式选择活门打开程度msv,选取q*J个工作点对非线性发动机模型进行线性化得到q*J个线性化模型再加入摄动块得到q*J个线性不确定性发动机模型,并对这q*J个线性不确定性发动机模型分别设计相应的二自由度H∞控制器从而组成二自由度H∞控制器组。
所述二自由度H∞控制器组调度计算模块根据输入的调度参数α和模式选择活门打开程度msv,利用内部设计的若干二自由度H∞控制器计算得到适应的二自由度H∞控制器,该二自由度H∞控制器根据参考输入r和测量参数y的差值e产生控制输入向量u。
优选的一种具体实现方式,可以根据输入的调度参数α和模式选择活门打开程度msv插值得到的适应的二自由度H∞控制器:
根据变循环发动机当前的调度参数α和模式选择活门打开程度msv选择相邻的四个设定工作点xi,j、xi,j+1、xi+1,j和xi+1,j+1(xi,j表示调度参数为αi,模式选择活门打开程度为msvj时的工作点),并获取四个设定工作点对应的线性控制器Ki,j、Ki,j+1、Ki+1,j和Ki+1,j+1(i=1,2,...,q,j=1,2,...,J)。根据公式
Figure BDA0002539278950000121
计算得到变循环发动机当前适应的二自由度H∞控制器K(α,msv)。
所述非线性发动机模型为:
Figure BDA0002539278950000122
y=g(x,u,msv)
其中
Figure BDA0002539278950000123
为控制输入向量,
Figure BDA0002539278950000124
为状态向量,
Figure BDA0002539278950000125
为输出向量,
Figure BDA0002539278950000126
为模式选择活门MSV打开程度,f(·)为表示系统动态的n维可微非线性向量函数,g(·)为产生系统输出的m维可微非线性向量函数。
所述一系列线性化模型为,从模式选择活门MSV完全闭合msv0到完全打开msvJ等间距的选取模式选择活门MSV打开程度msvj,j=1,2,...,J,在每一个固定的msvj,j=1,2,...,J,选择一组调度参数值αi,i=1,2,...,q,代表系统的动态范围,并将飞行包线划分为几个子区间,并将这些点作为工作点,在工作点,采用小扰动法或拟合法对发动机非线性模型在稳态参考点处进行线性化,得到
Figure BDA0002539278950000127
其中Δ表示该参数的变化量,系数矩阵可由下式得到:
Figure BDA0002539278950000128
Figure BDA0002539278950000129
这些系数在发动机不同的工作状态具有不同的值。
所述线性不确定性发动机模型得到方法为:得到的各个工作点的线性化模型为标称模型,线性状态空间方程
Figure BDA0002539278950000131
也可以表示为
Figure BDA0002539278950000132
实际的发动机和标称模型之间的误差可以表示为一个摄动块Δ。在标称模型加入摄动块建立发动机不确定模型
Figure BDA0002539278950000133
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (7)

1.一种变循环发动机增益调度二自由度H∞控制器,其特征在于:包括二自由度H∞控制器组调度计算模块;
其中二自由度H∞控制器组调度计算模块与变循环发动机本体以及变循环发动机上的若干传感器组成调度控制回路;由传感器输出调度参数α和模式选择活门打开程度msv至二自由度H∞控制器组调度计算模块;
所述二自由度H∞控制器组调度计算模块产生控制输入向量u并输出给变循环发动机本体,传感器得到变循环发动机测量参数y;
所述二自由度H∞控制器组调度计算模块内设计有若干二自由度H∞控制器,所述二自由度H∞控制器包含前置滤波器和反馈控制器,对若干不确定性发动机模型利用H∞回路成形方法分别设计得到;
所述线性不确定性发动机模型是对变循环发动机不同模式选择活门MSV打开程度、不同调度参数下的非线性发动机模型进行线性化后再加入摄动块得到的;
所述二自由度H∞控制器组调度计算模块根据输入的调度参数α和模式选择活门打开程度msv,利用内部设计的若干二自由度H∞控制器计算得到适应的二自由度H∞控制器,该二自由度H∞控制器根据参考输入r和测量参数y的差值e产生控制输入向量u。
2.根据权利要求1所述一种变循环发动机增益调度二自由度H∞控制器,其特征在于:所述二自由度H∞控制器组调度计算模块内设计若干二自由度H∞控制器的过程为:在全飞行包线内根据调度参数α和模式选择活门打开程度msv,选取q*J个工作点对非线性发动机模型进行线性化得到q*J个线性化模型再加入摄动块得到q*J个线性不确定性发动机模型,并对这q*J个线性不确定性发动机模型分别设计相应的二自由度H∞控制器从而组成二自由度H∞控制器组。
3.根据权利要求1或2所述一种变循环发动机增益调度二自由度H∞控制器,其特征在于:所述二自由度H∞控制器组调度计算模块根据输入的调度参数α和模式选择活门打开程度msv插值得到的适应的二自由度H∞控制器。
4.根据权利要求3所述一种变循环发动机增益调度二自由度H∞控制器,其特征在于:所述二自由度H∞控制器组调度计算模块根据变循环发动机当前的调度参数α和模式选择活门打开程度msv选择相邻的四个设定工作点xi,j、xi,j+1、xi+1,j和xi+1,j+1(xi,j表示调度参数为αi,模式选择活门打开程度为msvj时的工作点),并获取四个设定工作点对应的线性控制器Ki,j、Ki,j+1、Ki+1,j和Ki+1,j+1(i=1,2,...,q,j=1,2,...,J)。根据公式
Figure RE-FDA0002874479100000021
计算得到变循环发动机当前适应的二自由度H∞控制器K(α,msv)。
5.根据权利要求1所述一种变循环发动机增益调度二自由度H∞控制器,其特征在于:所述非线性发动机模型为:
Figure RE-FDA0002874479100000022
y=g(x,u,msv)
其中
Figure RE-FDA0002874479100000023
为控制输入向量,
Figure RE-FDA0002874479100000024
为状态向量,
Figure RE-FDA0002874479100000025
为输出向量,
Figure RE-FDA0002874479100000026
为模式选择活门MSV打开程度,f(·)为表示系统动态的n维可微非线性向量函数,g(·)为产生系统输出的m维可微非线性向量函数。
所述线性化过程为:从模式选择活门MSV完全闭合msv0到完全打开msvJ等间距的选取模式选择活门MSV打开程度msvj,j=1,2,...,J,在每一个固定的msvj,j=1,2,...,J,选择一组调度参数值αi,i=1,2,...,q,代表系统的动态范围,并将飞行包线划分为几个子区间,并将这些点作为工作点,在工作点,采用小扰动法或拟合法对发动机非线性模型在稳态参考点处进行线性化,得到
Figure RE-FDA0002874479100000027
其中Δ表示该参数的变化量,系数矩阵可由下式得到:
Figure RE-FDA0002874479100000031
Figure RE-FDA0002874479100000032
这些系数在发动机不同的工作状态具有不同的值。
所述线性不确定性发动机模型得到方法为:得到的各个工作点的线性化模型为标称模型,线性状态空间方程
Figure RE-FDA0002874479100000033
也可以表示为
Figure RE-FDA0002874479100000034
实际的发动机和标称模型之间的误差可以表示为一个摄动块Δ。在标称模型加入摄动块建立发动机不确定模型
Figure RE-FDA0002874479100000035
6.根据权利要求1所述一种变循环发动机增益调度二自由度H∞控制器,其特征在于:所述调度参数α包括变循环发动机的风扇转速或者压气机转速。
7.根据权利要求1所述一种变循环发动机增益调度二自由度H∞控制器,其特征在于:所述测量参数包括进气道出口、风扇出口、压气机出口、高压涡轮后、低压涡轮后的温度和压力,风扇转速和压气机转速。
CN202010542137.5A 2020-06-15 2020-06-15 变循环发动机增益调度二自由度h∞控制器 Active CN112327611B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010542137.5A CN112327611B (zh) 2020-06-15 2020-06-15 变循环发动机增益调度二自由度h∞控制器

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010542137.5A CN112327611B (zh) 2020-06-15 2020-06-15 变循环发动机增益调度二自由度h∞控制器

Publications (2)

Publication Number Publication Date
CN112327611A true CN112327611A (zh) 2021-02-05
CN112327611B CN112327611B (zh) 2022-06-10

Family

ID=74302999

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010542137.5A Active CN112327611B (zh) 2020-06-15 2020-06-15 变循环发动机增益调度二自由度h∞控制器

Country Status (1)

Country Link
CN (1) CN112327611B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114526164A (zh) * 2022-04-24 2022-05-24 中国航发四川燃气涡轮研究院 一种适用于双工作模式核心机的过渡态性能建模方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105631140A (zh) * 2015-12-30 2016-06-01 中国航空工业集团公司沈阳发动机设计研究所 一种变循环发动机稳态性能分析及优化方法
CN110083869A (zh) * 2019-03-27 2019-08-02 南京航空航天大学 一种评估模式变换对涡喷/涡扇变循环发动机稳定裕度影响的计算方法
CN111273554A (zh) * 2020-04-04 2020-06-12 西北工业大学 航空发动机最大推力状态降保守性二自由度h∞控制器
CN111271181A (zh) * 2020-04-04 2020-06-12 西北工业大学 航空发动机降保守性增益调度二自由度μ控制器

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105631140A (zh) * 2015-12-30 2016-06-01 中国航空工业集团公司沈阳发动机设计研究所 一种变循环发动机稳态性能分析及优化方法
CN110083869A (zh) * 2019-03-27 2019-08-02 南京航空航天大学 一种评估模式变换对涡喷/涡扇变循环发动机稳定裕度影响的计算方法
CN111273554A (zh) * 2020-04-04 2020-06-12 西北工业大学 航空发动机最大推力状态降保守性二自由度h∞控制器
CN111271181A (zh) * 2020-04-04 2020-06-12 西北工业大学 航空发动机降保守性增益调度二自由度μ控制器

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
ZHIDAN LIU等: "Design of Gain-Scheduling Robust Controller for Aircraft Engine", 《2019 CHINESE CONTROL CONFERENCE (CCC)》, 17 October 2019 (2019-10-17), pages 870 - 875 *
李嘉: "变循环航空发动机自适应控制技术研究", 《中国博士学位论文全文数据库(电子期刊)》, no. 2, 15 February 2019 (2019-02-15), pages 031 - 31 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114526164A (zh) * 2022-04-24 2022-05-24 中国航发四川燃气涡轮研究院 一种适用于双工作模式核心机的过渡态性能建模方法
CN114526164B (zh) * 2022-04-24 2022-07-26 中国航发四川燃气涡轮研究院 一种适用于双工作模式核心机的过渡态性能建模方法

Also Published As

Publication number Publication date
CN112327611B (zh) 2022-06-10

Similar Documents

Publication Publication Date Title
CN111271181B (zh) 航空发动机降保守性增益调度二自由度μ控制器
CN111273554B (zh) 航空发动机最大推力状态降保守性二自由度h∞控制器
CN111856918A (zh) 变循环发动机增益调度控制器
US5305599A (en) Pressure-ratio control of gas turbine engine
CN112613119B (zh) 一种多几何参数可调的进/排/发一体化航空推进系统建模方法
CN111856927B (zh) 变循环发动机增益调度二自由度μ控制器
CN111046568B (zh) 一种基于航空发动机及尾喷管机电液多系统联合仿真的控制参数优化方法
US5211007A (en) Method of pressure-ratio control of gas turbine engine
CN112085145B (zh) 自适应循环航空发动机模式转换过程控制规律设计方法
CN111680357B (zh) 一种变循环发动机机载实时模型的部件级无迭代构建方法
CN112327611B (zh) 变循环发动机增益调度二自由度h∞控制器
CN111852663A (zh) 变循环发动机降保守性鲁棒增益调度控制器
CN111880403A (zh) 航空发动机最大推力状态容错二自由度μ控制器
US9822731B2 (en) Control scheme using variable area turbine and exhaust nozzle to reduce drag
CN111856928B (zh) 变循环发动机鲁棒增益调度控制器
CN112327602A (zh) 变循环发动机气路部件故障增益调度容错控制器
CN111305954B (zh) 输入受限的航空发动机降保守性鲁棒增益调度控制器
CN111456857B (zh) 航空发动机降保守性增益调度二自由度h∞控制器
CN112711278B (zh) 变循环发动机模态转换恒定流量控制方法
CN115217635A (zh) 一种涡扇发动机全包线自适应加速控制方法
CN112904716A (zh) 一种变循环发动机最低油耗控制优化方法
CN114637210A (zh) 航空发动机多回路综合控制方法及系统
CN111852662A (zh) 航空发动机最大推力状态容错二自由度h∞控制器
CN111830827A (zh) 航空发动机容错增益调度二自由度μ控制器
CN115680903B (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