串联式TBCC发动机的模态转换控制方法

技术领域

本发明涉及一种串联式涡轮基组合循环(TBCC)发动机的模态转换控制方法，属于航空宇航推进系统控制技术领域。

背景技术

涡轮基组合循环(TBCC)发动机作为高超声速飞行器的核心部件，一直是组合动力领域关注和研究的热点。对于并联式涡轮基组合循环(TBCC)发动机，作为加速器的低速通道发动机需要在传统马赫数2一级涡轮发动机的基础上，进一步扩展飞行包线至马赫数4以上，以适应高速通道双模态冲压发动机的起动范围。低速通道发动机应在整个工作过程中以及工作模式的转换过程中具有足够的推力以维持整个推进系统的推力连续。随着飞行马赫数的增加，传统马赫数2一级涡轮发动机由于进口总温的增加导致高/低压转子相对换算转速下降，在低速-高速通道模态转换点附近已接近风车状态，共同工作点远离设计点，发动机推力极低[Haney J W.Highly reusable space transportation system study[C]//AIP Conference Proceedings.American Institute of Physics,1997,387(1):1277-1284.]。在这个过程中，极有可能出现发动机安装推力小于飞行器阻力(推阻不平衡)的情况导致飞行器持续减速，无法为高速通道发动机正常起动提供合适的来流条件，或者即使高速通道发动机能够起动但仍旧无法实现推阻平衡而导致模态转换失败。这类由于推阻不平衡导致飞行器无法完成低速通道向高速通道模态转换的问题被称为并联式TBCC发动机的推力陷阱问题[郑日恒,陈操斌.涡轮基组合循环发动机推力陷阱问题解决方案[J].火箭推进,2021,47(06):21-32.]，如图1所示。

目前对串联式TBCC发动机模态转换控制方法的研究以定点式模态转换控制为主，即在某一恒定的飞行马赫数和飞行高度条件(称为模态转换设计点)下完成模态转换，保持模态转换过程中发动机总推力恒定[Chen M,Tang H,Zhang K,et al.Turbine-basedcombined cycle propulsion system integration concept design[J].Proceedings ofthe Institution of Mechanical Engineers,Part G:Journal of AerospaceEngineering,2013,227(7):1068-1089.]。然而，定点式模态转换控制下，由于串联式TBCC发动机的亚燃模态最低工作马赫数的限制，模态转换设计点仍然处于涡轮发动机工作马赫数的上限，虽然相较并联式TBCC发动机的低速-高速通道模态转换马赫数有所降低，但此时涡轮发动机仍处于推力较低的工作区间，仍容易陷入推力陷阱。另一种可取的模态转换控制方式是区间式模态转换控制，即在某一较低的飞行马赫数下开始模态转换，随着马赫数的增加，当达到亚燃模态的最低工作马赫数时，完成模态转换。日本HYPR90-C发动机在2.5～3Ma之间进行模态转换，马赫数3以上时冲压发动机单独工作[Soga Y,Kurosaki M,Tsuzuki Y,et al.Control of HYPR demonstrator engine[C]//34th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit.1998:3283.]。黄红超[黄红超,王占学,蔡元虎.基于推力连续准则的小型涡轮冲压组合发动机模态转换过程分析[J].航空动力学报,2009,24(12):2756-2762.DOI:10.13224/j.cnki.jasp.2009.12.020.]分析了不同的流量调节阀开度对模态转换过程参数变化的影响，按照推力连续的准则,确定了模态转换过程流量调节阀开度随马赫数的变化规律。相较定点式模态转换，区间式模态转换控制方式由于在较低的马赫数区间打开冲压通道，转换过程中可以获得更大的发动机的推力，为解决推力陷阱问题提供了可能性。然而，目前针对区间式模态转换的研究较少，且公开文献中研究的侧重点大多集中在模态转换前后发动机推力和空气流量光滑衔接和平稳过渡上，尚未深入研究以解决推力陷阱，增加推力为目标的模态转换控制方案，同时，缺乏区间式模态转换方式增加推力的原理分析。对于作为并联式TBCC发动机低速通道加速器使用的串联式TBCC发动机，为避免推力陷阱，同时实现高加速性，亟需以实现高推力性能为目标来制定串联式TBCC发动机的模态转换控制方案。

发明内容

本发明所要及解决的技术问题在于克服现有技术不足，提供一种串联式TBCC发动机的模态转换控制方法，对于作为并联式TBCC发动机低速通道加速器使用的串联式TBCC发动机，可有效避免推力陷阱，同时实现高加速性。

本发明具体采用以下技术方案解决上述技术问题：

一种串联式TBCC发动机的模态转换控制方法，所述串联式TBCC发动机作为并联式TBCC发动机低速通道加速器使用；对串联式TBCC发动机的模态转换进行以下分区控制：

在来流马赫数为0～第一马赫数区间内，控制所述串联式TBCC发动机以涡扇模态最大状态工作，所述第一马赫数为发动机工作在最大状态下时，冲压涵道的回流裕度

在来流马赫数为第一马赫数～第二马赫数区间内，打开低速通道冲压涵道，并令低速通道冲压涵道开度大于0且低于或等于临界开度；所述临界开度是指在当前来流马赫数下令冲压涵道气流流量与涡轮通道气流流量的合流量等于进气道所能提供的最大进气量的冲压涵道开度，所述第二马赫数为满足低速通道冲压发动机单独工作的最小马赫数；

在来流马赫数为第二马赫数～4.0区间内，控制所述串联式TBCC发动机以完全冲压模态工作。

优选地，在来流马赫数为第一马赫数～第二马赫数区间内，令低速通道冲压涵道开度始终等于临界开度。

优选地，对涡扇涵道进行H

优选地，对低速通道冲压涵道开度进行开环控制。

相比现有技术，本发明技术方案具有以下有益效果：

本发明从原理角度分析并设计了以实现高推力为目标的作为并联式TBCC发动机低速通道的串联式TBCC发动机模态转换控制方案，可以在保持涡扇各部件工作在最大状态的同时，开启冲压涵道，通过增加发动机的进气流通能力，同时降低发动机溢流阻力的方式，提高发动机最大安装推力性能并实现更好的进气道-发动机匹配工作，从而有效解决并联式TBCC发动机的推力陷阱问题；本发明还进一步提出了最优的冲压涵道开度调节方案，可最大程度地提升模态转换过程中发动机的推力，实现进气道进气量的最佳有效利用。

附图说明

图1为并联式TBCC发动机的推力陷阱问题示意图；

图2为动压为34.29kPa的飞行轨迹；

图3为常规定点式模态转换控制计划下的发动机推力F与Ma之间的关系示意图；

图4为涡扇模态下发动机关键参数随马赫数变化情况示意图；

图5为回流裕度RM随Ma的变化示意图；

图6为本发明所提出的增推控制结构示意图；

图7为H∞鲁棒控制的结构原理示意图；

图8为Ma＝2.5时RVABI开度对发动机特性的影响；

图9为不同来流马赫数条件下，RVABI开度对发动机特性的影响；

图10为RVABI开度最优调节规律；

图11为外并联增推式TBCC发动机部件级模型结构原理示意图；

图12为飞行轨迹及模态划分示意图；

图13为并联式TBCC发动机采用不同模态转换控制方案的仿真对比结果。

具体实施方式

作为并联式TBCC发动机低速通道加速器使用的串联式TBCC发动机，其任务是将并联式TBCC发动机加速至高速通道双模态超燃冲压发动机最小起动马赫数Ma＝4，在此过程中提供足够的推力避免产生推力陷阱；针对这一问题，本发明的解决思路是采用分区间的模态转换控制方案，具体如下：

对串联式TBCC发动机的模态转换进行以下分区控制：

在来流马赫数为0～第一马赫数区间内，控制所述串联式TBCC发动机以涡扇模态最大状态工作，所述第一马赫数为发动机工作在最大状态下时，冲压涵道的回流裕度

在来流马赫数为第一马赫数～第二马赫数区间内，打开低速通道冲压涵道，并令低速通道冲压涵道开度大于0且低于或等于临界开度；所述临界开度是指在当前来流马赫数下令冲压涵道气流流量与涡轮通道气流流量的合流量等于进气道所能提供的最大进气量的冲压涵道开度，所述第二马赫数为满足低速通道冲压发动机单独工作的最小马赫数；

在来流马赫数为第二马赫数～4.0区间内，控制所述串联式TBCC发动机以完全冲压模态工作。

优选地，在来流马赫数为第一马赫数～第二马赫数区间内，令低速通道冲压涵道开度始终等于临界开度。

优选地，对涡扇涵道进行H

优选地，对低速通道冲压涵道开度进行开环控制。

为了便于公众理解，下面结合附图对本发明的技术方案进行详细说明：

以TBCC推进系统为动力的高超声速飞行器一般是沿着等动压线进行加速爬升。STARRUNNER概念机飞行轨迹的动压约为95kPa，HYPR计划中飞行轨迹的动压约为20kPa，提高飞行器的动压可以对其性能进行充分的开发，但是动压的上升也会使得飞行器所承受的机械负荷以及温度载荷的上升，所以飞行器必须在合适的最大动压下进行爬升。动压q

式中，ρ

以高超声速飞行器沿动压值q

在追求发动机高推力的同时必须考虑发动机状态参数限制，串联式TBCC发动机工作过程中的主要限制包括最大机械/气动负荷限制、最大热负荷限制以及气动稳定性限制。基于现有发动机水平，须满足的限制参数值如表2所示。

表2发动机状态参数限制

沿上述q

沿飞行轨迹，发动机推力变化如图3所示：

如图3所示，随着马赫数的增加，涡扇模态下发动机推力显著减小，在模态转换点(Ma＝3,H＝20km),涡扇发动机最大推力仅为地面状态的41％。Ma>3后进入到冲压模态，发动机推力重新上升。整个过程中，在模态转换点附近推力下降严重，极易出现推阻不平衡现象。

高马赫数条件下涡扇模态推力不足的原因主要包括两个方面:一方面，随着飞行马赫数Ma增加，来流空气总温T

式中，n

由以上分析可知，串联式TBCC发动机工作过程中，在模态转换点附近发动机总推力出现明显下降，且进气道工作在亚临界状态，溢流阻力较大，导致发动机安装推力下降严重，极易出现推阻不平衡现象。因此，极有必要在高马赫数条件下开展增推控制方案研究工作。

对于一定高度和马赫数的飞行状态，串联式TBCC发动机推力公式可表示为：

F＝W

式中W

由温度，压力的静参数与迟滞参数关系式公式(4)，可得喷管出口马赫数Ma

喷管出口截面当地声速

忽略燃气流量与空气流量的差别，则有：

当喷管扩张段为全管超声速流时，喷管出口仅可能出现三种状态：出口斜激波、设计状态、出口膨胀波，而这三种状态对应的喷管出口静压p

进一步假设喷管与外界无热交换。则有T

式中，σ

由公式(9)可知，如果在涡扇模态高马赫数推力不足的工作区间打开冲压涵道，同时保证风扇进口流量W

基于以上原理分析，为得到最大输出推力，需保证风扇进口流量W

W

另外，在一定的飞行条件下，风扇进口总压p

如式(12)所示，通过闭环控制W

为了保证增推控制过程平稳，安全进行，增推控制的适用条件除了需要满足发动机最大机械/气动负荷限制、最大热负荷限制以及气动稳定性限制(见表2)外，还需要满足以下三个条件：

1)冲压涵道气流不倒流(若涡轮后的高温燃气经过冲压涵道回流到风扇进口，即形成了燃气倒流现象，将造成发动机风扇超温损坏)。2)转换过程中加力/冲压燃烧室燃烧稳定。3)流入发动机的总流量W

不倒流是增推控制过程中打开冲压涵道的基本要求，定义冲压涵道的回流裕度RM如式(13)所示。为了保证不倒流，需要保证进气道出口总压P

回流裕度RM主要受发动机工作状态和飞行马赫数Ma的影响，当发动机工作在最大状态时，RM的大小主要由飞行马赫数决定；不同马赫数下发动机RM如图5所示。

由图5可以看出，当飞行马赫数小于2.2，回流裕度RM<2％，此时若打开冲压涵道会发生燃气倒流现象；当飞行马赫数大于2.2时，RM≥2％，此时满足打开冲压涵道的条件。随着飞行马赫数的增加，RM不断增加,说明飞行马赫数越高，冲压涵道越不容易产生回流。

为了实现上述增推控制方案，在满足增推控制方案适用条件的工作区间，设计了如图6所示的增推控制结构。该增推控制结构分为两部分，包括涡扇涵道的H

H

其中，ω是扰动信号，需要设计合适的控制器K(s)抑制其对输出y的影响，保证被控对象P(s)稳定。为了实现无静差跟踪，W

为了保证发动机全程安全稳定运行，控制结构中加入了发动机限制保护回路，包括n

采用图6所示增推控制结构，在等动压线上不同来流马赫数时，仿真后可变面积涵道引射器(RVABI)不同开度下对发动机性能的影响，RVABI开度α定义如式(14)所示。

其中，A

首先选取Ma＝2.5飞行条件，仿真RVABI不同开度下对发动机性能参数的影响，如图8所示。随着α的放大，冲压涵道空气流量W

α＝α

此时继续增大α，将导致涡扇发动机空气流量下降，进气道工作在超临界状态，总压恢复系数σ下降，进而导致推力F下降，耗油率S

图9仿真了不同来流马赫数条件下，RVABI开度对发动机特性的影响。如图9所示，在不同的来流马赫数下，本发明提出的增推控制方案均可以增大发动机推力F，同时，降低耗油率S

由图9中的(b)、(c)可以看出，不同飞行马赫数下均存在一个RVABI的临界开度α

以上述的串联式TBCC发动机作为低速通道加速器，进一步开发了并联式TBCC发动机部件级模型，如图11所示。

低速通道采用上文建立的串联式TBCC发动机模型，高速通道为双模态超燃冲压发动机，高速通道亚燃模态下需满足公式(17)尾喷管流量平衡方程；动态过程中，需满足燃烧室容积动力学方程，如公式(18)所示。超燃模态下由于高速通道内气流全程超声速，上游气流不会受到下游气流条件的影响，因此不需要迭代计算。

式中，下标“s”指双模态超燃冲压发动机截面参数，W

沿等动压线飞行轨迹，外并联增推式TBCC发动机的工作模态划分如图12所示。低速通道串联式TBCC发动机向高速通道双模态超燃冲压发动机转换的马赫数为Ma＝4。工作模态随着来流马赫数提高可分为以下几个阶段：(1)Ma＝0～2.2区间内，处于涡扇模态，以涡扇模态最大状态工作；(2)Ma＝2.2～3.0区间内，打开低速通道冲压涵道，定义此模态为增推模态，低速通道冲压涵道开度按图10的最优调节规律变化；(3)Ma3.0～4.0区间内，满足低速通道冲压发动机单独工作来流条件，串联式TBCC发动机进入完全冲压模态；(4)Ma4.0以上，串联式TBCC发动机停止工作，进入高速通道双模态超燃冲压工作模态。

沿上述飞行轨迹，针对并联式TBCC发动机，分别采用低速通道的常规控制方案和本发明控制方案，推力、空气流量、耗油率变化的对比如图13所示：

由图13中的(a)可以看出，采用常规的控制方案，当飞行马赫数Ma>2.7时，由于发动机安装推力小于飞行器阻力，出现推阻不平衡现象致使飞行器持续减速，无法为冲压发动机正常起动提供合适的来流条件，陷入推力陷阱。当采用增推控制计划后，发动机推力最大增加量ΔF