技术领域
[0001] 本
发明涉及一种航空发动机过渡态控制规律优化的设计方法，属于航空发动机过渡态优化及控制技术领域。
背景技术
[0002] 本发明依托背景为某型小涵道比涡扇发动机的非线性部件级数学模型的过渡态控制规律优化。一般来说，在发动机控制中，主要包括稳态控制和过渡态控制。稳态控制是保证发动机在某一稳态点受到扰动时使其性能不至于
波动从而能够恢复到稳定状态，是一种小偏离控制问题。在发动机的过渡态控制，发动机部分或全部性能随着时间发生改变，通常所说的加减速性能就是发动机过渡态性能的主要表现形式，而发动机的加减速性能直接影响着飞机
起飞加速性能的好坏。为了获得良好的过渡态性能，应正确设计发动机的加减速计划，使发动机在满足约束条件下，尽可能保证其从一个工作状态到另一个工作状态的过渡态时间最短。这对于军用飞机和作战飞机而言，能够保证其快速性的作战需求；同样对于民用飞机，出于安全飞行的
角度，也要保证发动机过渡态时间最短。在发动机过渡态过程中的加减速计划的设计对过渡态性能有重要影响，因此，对发动机过渡态控制规律的研究很有必要。
[0003] 在发动机的过渡态中，加减速计划是根据经验上下平移稳态工作线且不超一些限制边界线如过喘振计划和贫油熄火线而设定的，这就导致加减速计划的设计在一定程度上留有很大的裕度，使得发动机的过渡态性能不能充分发挥。依据现有的文献，对发动机过渡态的优化方法有功率提取法，但该方法无法考虑发动机部件容积效应影响；此外还有动态稳定法，但该方法需提前已知诸多耦合因素，设计很不方便。本发明采用一种基于SQP方法，对发动机过渡态控制规律进行优化设计，使在发动机在不超过约束边界条件，获得最短的过渡态时间，改善发动机过渡态性能。同时在不超过边界限制的前提下，调整燃油流量速率限制值，设计非最短过渡态时间下的加减速计划。同时将设计出的不同过渡态时间下的加减速计划可以用到发动机过渡态闭环控制中，为其提供一定的理论依据，达到以不同的过渡态时间作为目标下的全闭环控制，实现过渡态优化与闭环控制的结合。
发明内容
[0004] 为了保证发动机过渡态不超限，同时能够满足发动机过渡态时间要求，以及针对发动机在加减速过程中开环控制问题，本发明提供了一种航空发动机过渡态控制规律优化的设计方法。
[0005] 本发明的技术方案为：
[0006] 一种航空发动机过渡态控制规律优化的设计方法，步骤如下：
[0007] S1.基于SQP
算法，对发动机的过渡态进行优化，实现沿约束边界运行的发动机过渡态控制规律设计；
[0008] 沿约束边界运行的发动机过渡态控制规律设计步骤如下：
[0009] S1.1确定优化
节点的个数以及整个优化过程时间的长短；
[0010] S1.2将发动机在闭环
控制器作用下得到的燃油流量作为参考，选取燃油流量初值，使最终确定的燃油流量初值加载到发动机模型后，发动机各项输出不超过限制值；
[0011] S1.3对SQP算法的有关优化选项最大函数计算次数MaxFunEvals、最大
迭代次数MaxIter和函数
精度TolFun参数进行设置；
[0012] S1.4确定发动机过渡态控制规律优化的边界条件，包括
压气机物理转N2的限制、
风扇物理转速N1限制、压气机出口压
力P3限制、低压
涡轮出口
温度T5限制、风扇喘振SMF喘振边界、压气机喘振SMC喘振边界、燃油流量Wf限制和过渡态燃油流量速率△Wf限制条件，并依据限制条件建立相应的约束函数，其形式如下所示：
[0013] N2,min≤N2≤N2,max
[0014] N1≤N1,max
[0015] P3≤P3,max
[0016] T5≤T5,max
[0017] SMF≥SMFmin
[0018] SMC≥SMCmin
[0019] Wf,min≤Wf≤Wf,max
[0020] ΔWf≤ΔWf,max
[0021] 其中，N2,min为压气机物理最低物理转速，N2,max为压气机最高物理转速，N1,max为风扇最低物理转速，P3,max为压气机出口压力最大值，T5,max为低压涡轮出口温度最大值，SMCmin压气机
喘振裕度最小值，SMFmin风扇喘振裕度最小值，Wf,min和Wf,max分别为燃油流量最小值和最大值，ΔWf,max为燃油流量速率最大值；
[0022] S1.5结合发动机过渡态对时间的要求，建立发动机过渡态控制规律优化的目标函数，其形式如下：
[0023] J＝100*norm(N2-N2,cmd,2)+100*norm(N1-N1,cmd,2)
[0024] 其中，N1,cmd和N2,cmd为风扇物理转速和压气机物理转速的期望值；
[0025] S2.在不超过发动机限制边界的前提下，调整发动机的燃油流量速率限制值而其他约束条件不变，由此建立不同燃油速率限制条件下的发动机过渡态控制规律；
[0026] 建立不同燃油速率限制条件下的发动机过渡态控制规律设计步骤如下：
[0027] S2.1优化节点个数以及整个优化过程时间、燃油流量初值以及各个优化选项的设定与沿约束边界下的过渡态优化保持一致；
[0028] S2.2发动机过渡态中的压气机物理转速N2的限制、风扇物理转速N1限制、压气机出口压力P3限制、低压涡轮出口温度T5限制、风扇喘振SMF喘振边界、压气机喘振SMC喘振边界和燃油流量Wf限制条件不变，在不超过过渡态燃油流量速率最大限制值ΔWf,max的前提下，对燃油流量速率限制进行改变，分别设为ΔWf,1，ΔWf,2，ΔWf,3…ΔWf,N，在不同燃油流量速率限制条件下重新进行过渡态优化，得到不同条件下的控制规律即燃油流量Wf,1，Wf,2，Wf,3…Wf,N；
[0029] S2.3将在不同燃油流量速率限制条件下优化得到的发动机燃油流量Wf,1，Wf,2，Wf,3…Wf,N进行保存，为之后的燃油流量调度做准备；
[0030] S3.完成不同燃油速率限制条件下的过渡态控制规律设计后，构建发动机过渡态时间的评估函数，从而确定不同过渡态控制规律下的过渡态时间；
[0031] 构建发动机过渡态时间的评估函数步骤如下：
[0032] S3.1产生发动机模型的输出：将发动机的燃油流量Wf加载到发动机模型中，得到发动机的输出：风扇相对物理转速N1、风扇相对换算转速N1cor、压气机相对物理转速N2、压气机相对发动机进口的换算转速N2cor；
[0033] S3.2以发动机输出中的压气机相对物理转速N2作为发动机过渡态时间的评估变量，当其稳定在其稳态值的上下浮动0.2％则默认进入稳定状态，则进入时间为发动机的过渡态时间；
[0034] S3.3将不同限制条件下的发动机过渡态控制规律中的燃油流量Wf,1，Wf,2，Wf,3…Wf,N加载到发动机模型中，得到不同限制条件下发动机模型输出N2,1，N2,2，N2,3…N2,N，计算不同限制条件下的过渡态时间Ts,1，Ts,2，Ts,3…Ts,N并保存；
[0035] S4.对不同的过渡态控制规律下的燃油流量和其对应的过渡态时间从小到大排序，并建立lookuptable插值表；根据建立的插值表，实现在线过渡态控制计划调度，即输入期望的发动机过渡态时间调度得到该过渡态时间下对应的燃油流量；
[0036] 实现在线过渡态控制计划调度的步骤如下：
[0037] S4.1将不同限制条件下燃油流量以及其对应的过渡态时间Ts,1，Ts,2，Ts,3…Ts,N按从小到大的顺序排列；
[0038] S4.2以各个过渡态时间Ts,1，Ts,2，Ts,3…Ts,N作为插值节点1，对应的燃油流量Wf,1，Wf,2，Wf,3…Wf,N作为Tabledata，由于加载到发动机模型的燃油流量是时间序列格式，将燃油流量的时间序列作为插值节点2，建立lookuptable插值表；
[0039] S4.3选择期望的过渡态时间 通过建立的插值表进行插值，得到与过渡态时间对应的一组燃油流量 实现在线过渡态时间调度；
[0040] S4.4将期望过渡态时间下的燃油流量 作为发动机的输入，保存发动机的相关输出，以保存的输出作为闭环控制器的参考指令；
[0041] S5.将期望过渡态时间下的燃油流量作为发动机输入，得到发动机的输出，相关输出参数作为发动机闭环控制的参考指令，实现发动机的全闭环控制。
[0042] 本发明的有益效果：
[0043] 本发明提出的一种航空发动机过渡态控制规律优化的设计方法能够在发动机过渡态不超限的条件下，通过优化设计出满足过渡态时间要求的过渡态控制规律，和以往的发动机过渡态沿边界运行，只能寻找出最短的过渡态时间相比，本发明在发动机过渡态时间的选择上更加灵活，可以依据过渡态时间进行在线过渡态控制计划调度；此外，航空发动机过渡态优化方法具有一定的普适性，可以推广到其他类型的航空涡轮发动机过渡态优化；同时由过渡态优化而设计出的控制规律可以作为发动机全闭环控制中参考指令，解决了以往基于加减速计划编排的过渡态开环控制的弊端，为发动机闭环控制提供了可靠的理论依据，实现了过渡态优化与全闭环控制的结合。
附图说明
[0044] 图1为航空发动机过渡态控制规律优化的设计方法结构示意图；
[0045] 图2为航空发动机过渡态控制规律优化的设计方法的整体
流程图；
[0046] 图3为沿约束边界运行的发动机过渡态控制规律设计算法流程图；
[0047] 图4为发动机在线过渡态控制计划调度流程图；
[0048] 图5为依据优化得到的燃油流量和计算出的过渡态时间而建立的插值表三维曲线图；
[0049] 图6为基于此种方法优化得到的微型
燃气涡轮发动机的过渡态控制规律。
具体实施方式
[0050] 下面结合附图对本发明作进一步说明，系统结构图如图1所示。
[0051] 如图2所示，一种航空发动机过渡态控制规律优化的设计方法，包括以下步骤：
[0052] S1.采用一种基于SQP的算法，对发动机的过渡态进行优化，实现沿约束边界运行的发动机过渡态控制规律设计；
[0053] S2.在不超过发动机限制边界的前提下，调整发动机的燃油速率限制值而其他约束条件不变，通过优化建立不同限制条件下的过渡态控制规律；
[0054] S3.完成不同限制条件下的过渡态控制规律设计后，构建发动机过渡态时间的评估函数，从而确定不同过渡态控制规律下的过渡态时间；
[0055] S4.对不同的过渡态控制规律下的燃油流量和其对应的过渡态时间从小到大排序，并建立lookuptable插值表；根据建立的插值表，实现在线过渡态控制计划调度，即输入期望的发动机过渡态时间调度得到该过渡态时间下对应的燃油流量；
[0056] S5.将期望过渡态时间下的燃油流量作为发动机输入，得到发动机的输出，相关输出参数可作为发动机闭环控制的参考指令，实现发动机的全闭环控制；
[0057] 如图3所示，沿约束边界运行的发动机过渡态控制规律设计步骤是：
[0058] S1.确定优化的过渡态时间为T＝8秒，由于模型本身计算步长为0.025秒，采取模型计算点数和优化节点数相等的原则，由此确定优化节点数为0：0.025：8，总计321个；
[0059] S2.将发动机在闭环控制器作用下得到的燃油流量大小作为参考，取优化的燃油流量初值为数值大小59.1Kg/h，长度为321的常值序列，经验证在该燃油的作用下，发动机各项输出不超过限制值，该值为可行解；
[0060] S3.对SQP方法的有关优化选项最大函数计算次数MaxFunEvals、最大迭代次数MaxIter和函数精度TolFun等参数进行设置，取MaxFunEvals＝INF，MaxIter＝INF，TolFun＝1e-8；
[0061] S4.确定发动机过渡态控制规律优化的边界条件，主要包括压气机物理转速N2的限制、风扇物理转速N1限制、压气机出口压力P3限制、低压涡轮出口温度T5限制、风扇喘振SMF喘振边界、压气机喘振SMC喘振边界、燃油流量Wf限制、过渡态燃油流量速率△Wf限制等条件，并依据限制条件建立相应的约束函数，其形式如下所示：
[0062] N2,min≤N2≤N2,max
[0063] N1≤N1,max
[0064] P3≤P3,max
[0065] T5≤T5,max
[0066] SMF≥SMFmin
[0067] SMC≥SMCmin
[0068] Wf,min≤Wf≤Wf,max
[0069] ΔWf≤ΔWf,max
[0070] 其中，N2,min为压气机物理最低物理转速，N2,max为压气机最高物理转速，N1,max为风扇最低物理转速，P3,max为压气机出口压力最大值，T5,max为低压涡轮出口温度最大值，SMCmin压气机喘振裕度最小值，SMFmin风扇喘振裕度最小值，Wf,min和Wf,max分别为燃油流量最小值和最大值，ΔWf,max为燃油流量速率最大值，取N2,min＝68.9％，N2,max＝102％，N1,max＝102％，P3,max＝1310kPa，T5,max＝873K，Wf,min＝39.6kg/h，Wf,max＝465kg/h，ΔWf,max＝2.5kg/h/25ms，SMFmin＝3％，SMCmin＝3％则约束函数表示为：
[0071] 68.9％≤N2≤102％
[0072] N1≤102％
[0073] P3≤1310
[0074] T5≤873
[0075] SMF≥3％
[0076] SMC≥3％
[0077] 39.6≤Wf≤465
[0078] ΔWf≤2.5
[0079] S5.结合发动机过渡态对时间的要求，建立发动机过渡态控制规律优化的目标函数,其形式如下：
[0080] J＝100*norm(N2-N2,cmd,2)+100*norm(N1-N1,cmd,2)
[0081] 其中N1,cmd和N2,cmd为风扇物理转速和压气机物理转速的期望值，取N2,cmd＝100，N1,cmd,＝100，则目标函数表示为
[0082] J＝100*norm(N2-100,2)+100*norm(N1-100,2)
[0083] 建立不同限制条件下的过渡态控制规律设计步骤是：
[0084] S1.发动机的优化时间和优化结点数的设定，燃油流量初值的选取以及各个优化选项的设定和沿约束边界下的过渡态优化保持一致；
[0085] S2.发动机过渡态中的压气机物理转速N2的限制、风扇物理转速N1限制、压气机出口压力P3限制、低压涡轮出口温度T5限制、风扇喘振SMF喘振边界、压气机喘振SMC喘振边界、燃油流量Wf限制条件不变，在不超过过渡态燃油流量速率最大限制值ΔWf,max＝2.5的前提下，对燃油流量速率限制进行改变，分别设为ΔWf,1＝2.2，ΔWf,2＝2.0，ΔWf,3＝1.8，ΔWf,3＝1.6，ΔWf,3＝1.4，ΔWf,3＝1.2，在不同燃油速率限制条件下重新进行过渡态优化，得到不同条件下的控制规律即燃油流量Wf,1，Wf,2，Wf,3，Wf,4，Wf,5；
[0086] S3.将在不同燃油流量速率限制条件下优化得到的发动机燃油流量Wf,1，Wf,2，Wf,3，Wf,4，Wf,5进行保存；
[0087] 构建发动机过渡态时间的评估函数步骤是：
[0088] S1产生发动机模型的输出：将优化后的发动机的燃油流量Wf加载到发动机模型中，得到发动机的输出：风扇相对物理转速N1、风扇相对换算转速N1cor、压气机相对物理转速N2、压气机相对发动机进口的换算转速N2cor等；
[0089] S2以发动机输出中的压气机相对物理转速N2作为发动机过渡态时间的评估变量，当其稳定在其稳态值的上下0.2％则默认进入稳定状态，则进入时间为发动机的过渡态时间；
[0090] S3将不同限制条件下的发动机过渡态控制规律中的燃油流量Wf,1，Wf,2，Wf,3，Wf,4，Wf,5加载到发动机模型中，得到不同限制条件下发动机模型输出N2,1，N2,2，N2,3，N2,4，N2,5，计算不同限制条件下的过渡态时间Ts,1，Ts,2，Ts,3，Ts,4，Ts,5并保存；
[0091] 如图4所示，实现在线过渡态控制计划调度的步骤是：
[0092] S1.将不同限制条件下燃油流量以及其对应的过渡态时间Ts,1，Ts,2，Ts,3…Ts,N按从小到大的顺序排列；
[0093] S2.以各个过渡态时间Ts,1，Ts,2，Ts,3…Ts,N作为插值节点1，对应的燃油流量Wf,1，Wf,2，Wf,3…Wf,N作为Tabledata，由于加载到模型的燃油流量是时间序列格式，将燃油流量的时间序列作为插值节点2，建立lookuptable插值表，如图5所示，以过渡态时间作为调度量，依据期望的过渡态时间调度出对应的燃油流量；
[0094] S3.选择期望的过渡态时间 通过建立的插值表进行插值，得到与过渡态时间对应的一组燃油流量 实现在线过渡态时间调度；
[0095] S4.将期望过渡态时间下的燃油流量 作为发动机的输入，保存发动机的相关输出，以此作为发动机闭环控制的参考指令；
[0096] 图5为依据优化得到的燃油流量和计算出的过渡态时间而建立的插值表三维曲线图；
[0097] 图6为基于此种方法优化得到的微型燃气涡轮发动机的过渡态控制规律；
[0098] 综上可见本发明提出的航空发动机过渡态控制规律化方法是有效的、可行的，并且具有普适性，能够应用到其他类型的发动机的过渡态控制规律优化中。
