基于贝塞尔曲线的涡扇发动机加速过程控制律设计方法技术

本发明专利技术公开了一种基于贝塞尔曲线的涡扇发动机加速过程控制律设计方法，包括以下步骤：1、确定涡扇发动机在加速过程的起点时刻以及涡扇发动机各控制量在起点时刻和终点时刻的状态，并初始化涡扇发动机各控制量在时域空间的贝塞尔曲线阶数；2、以缩短控制量从初始值到目标值的耗时为目标，并结合相应约束条件，优化求解涡扇发动机各控制量在时域空间的贝塞尔曲线；3、改变涡扇发动机各控制量的贝塞尔曲线阶数，转至步骤2；4、反复迭代多次后，选择使得所述优化目标值最小的迭代步中所得到的涡扇发动机各控制量的最优贝塞尔曲线。相比现有技术，本发明专利技术从全局优化的角度出发，可以有效提高发动机在响应过程中的推力响应速率。

【技术实现步骤摘要】

本专利技术属于航空宇航推进理论与工程中的系统控制与仿真
，具体涉及一种基于贝塞尔曲线的涡扇发动机加速过程控制律设计方法。
技术介绍
在航空发动机控制系统设计中，对加速性能的要求越来越高，航空发动机是非常复杂的热力机械系统，其加速时的工作状态是强非线性，因而加速控制是航空发动机过渡态控制中最复杂的控制问题。在应急控制中，如紧急爬升、飞机复飞、追击、摆脱敌机等情况，发动机的加速性能无疑是一项极为重要指标。同时，它是研究其它过渡态控制，如进出加力、减速过程、起动过程的基础。因此，对于航空发动机加速过程的研究是航空发动机控制技术的一项重要研究内容。近年来，非线性规划理论开始应用于航空发动机优化控制，并显现出其在求解复杂非线性对象优化问题方面的极大优势。这类算法可以直接对非线性对象进行优化控制计算，并且可以在优化控制过程中充分考虑各种约束条件，因此，很多学者在该方面做了大量的工作，成功地将罚函数法、约束变尺度法、遗传算法、二次序列规划法、可行二次序列规划法、功率提取法和动稳态法等方法应用于加速优化控制中，并取得良好效果。其研究的主要思路是，先把连续的优化问离散化，按时间顺序，在每一个离散点，建立能过提高加速性能的目标函数和能够保证发动机运行时的约束方程，运用优化算法对其进行求解，因此该方法存在一个致命缺陷，它只能寻求发动机在每个离散点取得最优值，而忽略全局离散点相互之间的联系，必然导致优化结果只能使发动机在某个阶段性能得到最优，然而，由于发动机的强非线性，加速过程中在某个阶段加速最快并不代表整个过程的加速时间最短。
技术实现思路
本专利技术所要解决的技术问题在于克服现有技术不足，提供一种基于贝塞尔曲线的涡扇发动机加速过程控制律设计方法，能够避免现有技术在优化过程单纯只寻求某个阶段性能最优的缺点，该方法从全局优化的角度出发，可以有效提高发动机在响应过程中的推力响应速率。本专利技术具体采用以下技术方案解决上述技术问题：基于贝塞尔曲线的涡扇发动机加速过程控制律设计方法，包括以下步骤：步骤1、确定涡扇发动机在加速过程的起点时刻以及涡扇发动机各控制量在起点时刻和终点时刻的状态，并初始化涡扇发动机各控制量在时域空间的贝塞尔曲线阶数；步骤2、根据以下优化目标并结合包括涡扇发动机的风扇和压气机的喘振裕度和最高转子转速的限制、涡轮部件高温限制、富油熄火限制、控制量控制范围限制、最大控制量变化率限制的约束条件，优化求解涡扇发动机各控制量在时域空间的贝塞尔曲线，得到涡扇发动机各控制量在当前贝塞尔曲线阶数下的最优贝塞尔曲线：minJ＝max(tu1(nu1),tu2(nu2),…,tuM(nuM))其中，tuj(nuj)表示涡扇发动机第j个控制量uj的nuj阶贝塞尔曲线起点到终点的时间，j＝1,2,…,M，M为涡扇发动机控制量的个数；步骤3、改变涡扇发动机各控制量的贝塞尔曲线阶数，转至步骤2；步骤4、反复迭代多次后，选择使得所述优化目标值最小的迭代步中所得到的涡扇发动机各控制量的最优贝塞尔曲线，即为涡扇发动机各控制量在加速过程的最优控制曲线。优选地，所述约束条件被N点离散化为以下约束表达式：s.t.1-Nf[k+1](u[k])Nf,max≥0;1-Nc[k+1](u[k])Nc,max≥0;Smf[k+1](u[k])Smf,min-1≥0;Smc[k+1](u[k])Smc,min-1≥0;1-T4[k+1](u[k])T4,max≥0;1-FAR4[k+1](u[k])FAR4,max≥0;umin≤u[k]≤umax;Δu[k]≤Δumax;k=1,2,...,N;tuj[i]≤tuj[i+1];i=0,1,...,nuj-1;j=1,2,...,M]]>其中，Nf[k+1](u[k])、Nc[k+1](u[k])、Smf[k+1](u[k])、Smc[k+1](u[k])、T4[k+1](u[k])、FAR4[k+1](u[k])分别表示第k+1个离散点处的风扇转速、压气机转子转速、风扇喘振裕度、压气机喘振裕度、涡轮前温度、燃烧室燃气比，其为第k个离散点处的控制量u[k]的函数；Nf,max、Nc,max、Smf,min、Smc,min、T4,max、FAR4,max分别为所允许的风扇最大转速、压气机转子最大转速、最小风扇喘振裕度、最小压气机喘振裕度、最大涡轮前温度、最大燃烧室燃气比；umin、umax分别为预设的最小、最大控制量；Δu[k]为第k个离散点处的控制量变化率；Δumax为控制量变化率上限；M为涡扇发动机控制量的个数；tuj[i]表示涡扇发动机第j个控制量uj的nuj阶贝塞尔曲线的起点到第i个控制点的时间。优选地，使用FSQP算法求解优化求解涡扇发动机各控制量在时域空间的贝塞尔曲线。优选地，在3～7的范围内改变涡扇发动机各控制量的贝塞尔曲线阶数。相比现有技术，本专利技术具有以下有益效果：本专利技术技术方案是基于Bezier(贝塞尔)控制量曲线的优化，以缩短加速过程中发动机控制量最后稳定不变的时间为优化目标，可以有效地克服常规优化控制方法只考虑某个阶段最优的缺点，从而提高发动机推力的整体响应速度。附图说明图1是具体实施例中所针对的发动机界面图；图2是推力F的响应曲线；图3是燃油Wfb的优化结果曲线；图4是尾喷管喉道面积A8的优化结果曲线；图5是涡轮前温度T4的响应曲线；图6是压气机转子转速Nc的响应曲线；图7是风扇喘振裕度Smf的响应曲线；图8是风扇转子转速Nf的响应曲线；图9是压气机喘振裕度Smc的响应曲线；图10是主燃烧室的油气比FAR4的响应曲线；图2～图10中的Wfb、A8、F、和T4作了归一化处理。具体实施方式下面结合附图对本专利技术的技术方案进行详细说明：针对现有技术在优化过程单纯只寻求某个阶段性能最优的缺点，本专利技术的思路是利用Bezier曲线设计涡扇发动机控制量在加速过程中的曲线，将发动机本文档来自技高网...

【技术保护点】
基于贝塞尔曲线的涡扇发动机加速过程控制律设计方法，其特征在于，包括以下步骤：步骤1、确定涡扇发动机在加速过程的起点时刻以及涡扇发动机各控制量在起点时刻和终点时刻的状态，并初始化涡扇发动机各控制量在时域空间的贝塞尔曲线阶数；步骤2、根据以下优化目标并结合包括涡扇发动机的风扇和压气机的喘振裕度和最高转子转速的限制、涡轮部件高温限制、富油熄火限制、控制量控制范围限制、最大控制量变化率限制的约束条件，优化求解涡扇发动机各控制量在时域空间的贝塞尔曲线，得到涡扇发动机各控制量在当前贝塞尔曲线阶数下的最优贝塞尔曲线：min J＝max(tu1(nu1),tu2(nu2),…,tuM(nuM))其中，tuj(nuj)表示涡扇发动机第j个控制量uj的nuj阶贝塞尔曲线起点到终点的时间，j＝1,2,…,M，M为涡扇发动机控制量的个数；步骤3、改变涡扇发动机各控制量的贝塞尔曲线阶数，转至步骤2；步骤4、反复迭代多次后，选择使得所述优化目标值最小的迭代步中所得到的涡扇发动机各控制量的最优贝塞尔曲线，即为涡扇发动机各控制量在加速过程的最优控制曲线。

【技术特征摘要】
1.基于贝塞尔曲线的涡扇发动机加速过程控制律设计方法，其特征在于，包括
以下步骤：
步骤1、确定涡扇发动机在加速过程的起点时刻以及涡扇发动机各控制量在起点
时刻和终点时刻的状态，并初始化涡扇发动机各控制量在时域空间的贝塞尔曲线
阶数；
步骤2、根据以下优化目标并结合包括涡扇发动机的风扇和压气机的喘振裕度和
最高转子转速的限制、涡轮部件高温限制、富油熄火限制、控制量控制范围限制、
最大控制量变化率限制的约束条件，优化求解涡扇发动机各控制量在时域空间的
贝塞尔曲线，得到涡扇发动机各控制量在当前贝塞尔曲线阶数下的最优贝塞尔曲
线：
minJ＝max(tu1(nu1),tu2(nu2),…,tuM(nuM))
其中，tuj(nuj)表示涡扇发动机第j个控制量uj的nuj阶贝塞尔曲线起点到终点的
时间，j＝1,2,…,M，M为涡扇发动机控制量的个数；
步骤3、改变涡扇发动机各控制量的贝塞尔曲线阶数，转至步骤2；
步骤4、反复迭代多次后，选择使得所述优化目标值最小的迭代步中所得到的涡
扇发动机各控制量的最优贝塞尔曲线，即为涡扇发动机各控制量在加速过程的最
优控制曲线。
2.如权利要求1所述涡扇发动机加速过程控制律设计方法，其特征在于，所述
约束条件被N点离散化为以下约束表达式：
s.t.1-Nf[k+1](u[k])Nf,max≥0;1-Nc[k+1](u[k])Nc,max≥0;Smf[k+1](u[k])Smf,min-1≥0;Smc[k+1](u[k])Smc,min-1≥0;1-T4[k+1](u[k])N4,...

【专利技术属性】
技术研发人员：郑前钢，缪丽祯，张海波，孙丰勇，李永进，叶志锋，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏;32