CN110513199B - 一种航空发动机推力衰退缓解的神经网络控制方法 - Google Patents
一种航空发动机推力衰退缓解的神经网络控制方法 Download PDFInfo
- Publication number
- CN110513199B CN110513199B CN201910664441.4A CN201910664441A CN110513199B CN 110513199 B CN110513199 B CN 110513199B CN 201910664441 A CN201910664441 A CN 201910664441A CN 110513199 B CN110513199 B CN 110513199B
- Authority
- CN
- China
- Prior art keywords
- engine
- thrust
- neural network
- performance
- recession
- 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.)
- Active
Links
- 238000013528 artificial neural network Methods 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000012937 correction Methods 0.000 claims abstract description 35
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 18
- 238000005457 optimization Methods 0.000 claims abstract description 16
- 230000000116 mitigating effect Effects 0.000 claims description 21
- 239000000446 fuel Substances 0.000 claims description 18
- 230000015556 catabolic process Effects 0.000 claims description 7
- 238000006731 degradation reaction Methods 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 7
- 238000011084 recovery Methods 0.000 claims description 7
- 238000012549 training Methods 0.000 claims description 7
- 238000010586 diagram Methods 0.000 claims description 6
- 238000013507 mapping Methods 0.000 claims 1
- 230000035772 mutation Effects 0.000 abstract description 7
- 238000004088 simulation Methods 0.000 description 7
- 239000000295 fuel oil Substances 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000012821 model calculation Methods 0.000 description 2
- 238000002485 combustion reaction Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02C—GAS-TURBINE PLANTS; AIR INTAKES FOR JET-PROPULSION PLANTS; CONTROLLING FUEL SUPPLY IN AIR-BREATHING JET-PROPULSION PLANTS
- F02C9/00—Controlling gas-turbine plants; Controlling fuel supply in air- breathing jet-propulsion plants
- F02C9/48—Control of fuel supply conjointly with another control of the plant
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05D—INDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
- F05D2270/00—Control
- F05D2270/70—Type of control algorithm
- F05D2270/709—Type of control algorithm with neural networks
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种航空发动机推力衰退缓解的神经网络控制方法,该方法包括:采用神经网络学习的NARMA‑L2模型建立内环控制系统;结合内环控制系统,设计基于变增量LP优化算法的外环指令修正模块,得到在发动机发生部件性能突变的情况下具有自调整能力的推力衰退缓解神经网络控制器。本发明解决了常规多变量控制器在发动机发生部件性能突变的情况下整机推力水平下降的问题,适用于在一定飞行包线内不同工作点的发动机推力衰退缓解控制,对于在不超温、不超转的情况下缓解由发动机部件性能突变引起的整机推力损失、提高发动机整机性能表现有着积极促进的作用。
Description
技术领域
本发明属于航空发动机推力缓解控制技术领域,尤其涉及一种航空发动机推力衰退缓解的神经网络控制方法。
背景技术
随着现代航空技术的飞速发展,飞行器需要性能更强、可靠性更高的航空发动机。由于发动机本身复杂的结构和其所处恶劣的工作环境,当发动机运行一段时间后,发动机的各个部件由于腐蚀、外物损伤等原因会发生不同程度的性能突变。这些原因使得实际发动机与同型号额定发动机之间会产生偏差,主要体现在部件的流通能力和效率的降低。为了确保飞行器的安全飞行、延长发动机使用寿命,需要对发生部件性能突变的发动机进行推力恢复,因此,研究航空发动机推力衰退缓解控制方法十分必要。
一种有效的推力缓解控制方法是通过内环控制转速,外环修正指令,进而达到推力缓解的目的。航空发动机是一个复杂的热力学系统,具有很强的不确定性和时变性,这使得内环控制器的设计变得困难。神经网络由于其良好的非线性逼近和泛化能力被广泛应用,带有反馈线性化的非线性自回归滑动平均(Nonlinear Auto regressive MovingAverage with Feedback Linearization,NARMA-L2)控制器是一种有效的人工神经网络控制器架构。在一定条件下,非线性系统的输入输出关系可以由NARMA-L2模型辨识得到,并且可以通过简单的数学变换得到控制律。然而由于NARMA-L2模型存在建模误差和训练误差,使得所设计的控制器性能在包线内应用受到影响。因此,本发明尝试提出一种具有在线修正能力的航空发动机性能衰退缓解控制器(OC-NARMA-L2),该方法利用梯度下降法在线修正控制器的神经网络参数,使其具有自适应的特性。
发明内容
针对上述技术问题,本发明提供一种航空发动机推力衰退缓解的神经网络控制方法,针对常规多变量控制器在发动机部件性能突变的情况下推力水平下降的问题,将基于变增量LP优化算法的外环指令修正模块同内环控制系统结合,得到具有自调整能力的推力衰退缓解神经网络控制器。该方法可在不超温、不超转的情况下缓解由发动机部件性能突变引起的整机推力损失,提高发动机整机性能表现。
技术方案:为实现上述目的,本发明采用的技术方案为:
一种航空发动机推力衰退缓解的神经网络控制方法,其特征在于,包括以下步骤:
步骤A)建立发动机高压转子转速、压比与尾喷管喉道面积之间的物理变量插值模型,以及利用主燃油量与高压转子转速之间的关系,采用神经网络建立性能衰退缓解的内环控制NARMA-L2模型;
步骤B)结合NARMA-L2模型,设计基于变增量LP算法的外环指令修正模块,开发在发动机发生部件性能突变的情况下具有自调整能力的推力衰退缓解神经网络控制器。
步骤C)发动机发生部件性能突变时,通过外环加入指令修正模块,对发动机指令进行修正,利用推力衰退缓解神经网络控制器追踪修正后的指令进而达到推力恢复的目的。
进一步的,所述步骤A)中利用主燃油量与高压转子转速之间的关系,采用神经网络建立性能衰退缓解的内环控制NARMA-L2模型具体步骤如下:
步骤A1),选取发动机主燃油量作为输入量u[k],高压转子转速作为输出量y[k],并根据所选输入输出量构建NARMA-L2模型为:
y[k]=f0(y[k-1],y[k-2],…,y[k-n],u[k-1],u[k-2],…,u[k-n])+g0(y[k-1],y[k-2],…,y[k-n],u[k-1],u[k-2],…,u[k-n])u[k]
其中,f0和g0可以通过训练神经网络组f0 *和g0 *近似得到,而转速控制律为:
其中,y*[k]为发动机期望输出;并根据发动机整机落压比指令、实际输出落压比和输出转速插值换算得到尾喷管喉道面积控制量;
步骤A2),利用期望输出和实际输出间的误差ec(k)对神经网络组f0 *和g0 *的拓扑结构参数进行在线修正:
W(k+1)=W(k)+αΔW(k)
V(k+1)=V(k)+αΔV(k)
其中,W和V分别为神经网络组的权值和阈值;α为学习率,其大小影响参数的调整速度,由多次试取法确定;ΔW(k)和ΔV(k)为雅可比矩阵。拓扑向量中的各元素在k时刻的更新如下:
wi(k+1)=wi(k)-αΔwi(k)
vi(k+1)=vi(k)-αΔvi(k)
其中,T={y(k),y(k-1),y(k-2),...,y(k-n),u(k-1),u(k-2),…,u(k-n)}。
进一步的,所述步骤B)中结合NARMA-L2模型,设计基于变增量LP算法的外环指令修正模块,得到在发动机发生部件性能突变的情况下具有自调整能力的推力衰退缓解神经网络控制器具体步骤如下:
步骤B1),构建基于变增量LP优化的航空发动机推力缓解神经网络控制系统结构图,包括NARMA-L2模型、内环神经网络控制器、发动机模型、性能状态滤波器、外环指令修正模块;
步骤B2),根据飞行条件及相关的控制计划确定发动机初始指令,利用发动机模型输出和限制条件设计基于变增量LP算法的外环指令修正模块,建立在发动机发生部件性能突变的情况下具有自调整能力的推力衰退缓解神经网络控制器。
进一步的,所述步骤B2)根据发动机模型输出和限制条件设计基于变增量LP算法的外环指令修正模块,建立在发动机发生部件性能突变的情况下具有自调整能力的推力衰退缓解神经网络控制器具体步骤如下:
步骤B2.1),对不可测推力和边界限制条件在工作点附近小范围进行线性化,提取性能指标、约束量和控制量之间的线性关系PSM,具体公式为:
其中,ΔSmf和ΔSmc分别为风扇和压气机的喘振裕度增量,ΔnL和ΔnH分别为低压、高压转子转速增量,ΔFn为推力增量,ΔTt6为涡轮后温度增量,ΔWf为燃油量增量,ΔA8为尾喷口面积增量,P矩阵元素可从发动机稳态模型中通过扰动法得出,具体公式为:
Pij={Mi(u0(j)+Δuj,...)-Mi(u0)}/Δuj
i=1,2,...,6,j=1,2
其中,Mi(·)是通过模型计算得到的输出矩阵;u0是控制量的初始值;Δuj是第j个控制量增量;
步骤B2.2),对于步骤B2.1)中的控制量增量Δuj,可由相应控制量uj乘以增量系数得出,增量系数的选取规则为:
其中,k为迭代次数;ΔFn为推力增量。根据步骤B2.1)中提取的PSM(推进系统矩阵,Propulsion System Matrix),将指令修正转换为LP优化问题,求解推力缓解所需指令的修正量。
进一步的,所述步骤B2.2)中将指令修正转换为LP优化问题,求解推力缓解所需指令的修正量具体步骤如下:
步骤B2.2.1),在发动机不超转、不超温、喘振裕度和控制量满足范围要求、推力不超过期望值的约束下寻找发动机最大推力,具体公式为:
其中,Fn为发动机推力,Wf为燃油量,A8为尾喷口面积,u=[Wf,A8]T为发动机控制向量,umin和umax分别为向量下界和上界,nLmin和nLmax分别为低压转子转速最低和最高值,nHmin和nHmax分别为高压转子转速最低和最高值,Tt6min和Tt6max分别为涡轮后温度最低和最高值,Smf和Smc分别为风扇和压气机喘振裕度,Fnor为实时模型计算出未发生突变故障的发动机推力期望值;
步骤B2.2.2),根据控制量约束与限制约束,利用步骤B2.1)中提取的PSM,将步骤B2.2.1)中的公式改写为:
同理求得最大推力性能指标为:
ΔFn=p61ΔWf+p62ΔA8
步骤B2.2.3),求解LP问题,经反复迭代至满足精度要求或达到最大迭代次数限制后得到全局最优解Wf和A8,并将最优解输入到发动机中计算得到对应的输出信号nH和EPR,将二者与推力恢复前指令信号作差,得到的ΔnHr和ΔEPRr作为转速指令和整机落压比指令的修正值输入到内环神经网络控制器,内环神经网络控制器跟踪修正后的指令,缓解发动机整机推力损失。
进一步的,步骤C)中,发动机发生部件性能突变时,先将发动机转速指令下降1%,然后加入外环修指令。
有益效果:本发明提出了基于变增量LP算法的航空发动机推力缓解神经网络控制器,在多变量控制器基础上设计指令修正外环回路,其中转速控制通过神经网络学习的NARMA-L2模型逆实现。当发动机性能发生突变时,首先将高压转子转速指令降低1%确保发动机安全工作,之后利用变增量LP优化算法在保证发动机不超限条件下得到高压转子转速指令和整机落压比指令的修正量,从而使发动机在性能退化情况下尽可能继续提供所需推力。
附图说明
图1是本发明基于变增量LP优化的航空发动机推力缓解神经网络控制系统结构图。
图2是内环NARMA-L2控制系统结构图。
图3是变增量LP算法结构图。
图4是H=0,Ma=0,风扇发生性能突变时供油量比较图。
图5是H=0,Ma=0,风扇发生性能突变时高压转子转速比较图。
图6是H=0,Ma=0,风扇发生性能突变时推力比较图。
图7是H=8000m,Ma=1.2,压气机发生性能突变时供油量比较图。
图8是H=8000m,Ma=1.2,压气机发生性能突变时高压转子转速比较图。
图9是H=8000m,Ma=1.2,压气机发生性能突变时推力比较图。
具体实施方式
下面结合附图对本发明的具体实施方式作更进一步的说明。
本发明说明的一种航空发动机推力衰退缓解的神经网络控制方法,包括以下步骤:
步骤A)建立发动机高压转子转速、压比与尾喷管喉道面积之间的物理变量插值模型,以及利用主燃油量与高压转子转速之间的关系,采用神经网络建立性能衰退缓解的内环控制NARMA-L2模型;
步骤B)结合NARMA-L2模型,设计基于变增量LP算法的外环指令修正模块,开发在发动机发生部件性能突变的情况下具有自调整能力的推力衰退缓解神经网络控制器。
步骤A1),选取发动机主燃油量作为输入量u[k],高压转子转速作为输出量y[k],并根据所选输入输出量构建NARMA-L2模型为:
y[k]=f0(y[k-1],y[k-2],…,y[k-n],u[k-1],u[k-2],…,u[k-n])+g0(y[k-1],y[k-2],…,y[k-n],u[k-1],u[k-2],…,u[k-n])u[k]
其中,f0和g0可以通过训练神经网络组f0 *和g0 *近似得到,而转速控制律为:
其中,y*[k]为发动机期望输出;并根据发动机整机落压比指令、实际输出落压比和输出转速插值换算得到尾喷管喉道面积控制量;
步骤A2),利用期望输出和实际输出间的误差ec(k)对神经网络组f0 *和g0 *的拓扑结构参数进行在线修正:
W(k+1)=W(k)+αΔW(k)
V(k+1)=V(k)+αΔV(k)
其中,W和V分别为神经网络组的权值和阈值;α为学习率,其大小影响参数的调整速度,由多次试取法确定;ΔW(k)和ΔV(k)为雅可比矩阵。拓扑向量中的各元素在k时刻的更新如下:
wi(k+1)=wi(k)-αΔwi(k)
vi(k+1)=vi(k)-αΔvi(k)
其中,T={y(k),y(k-1),y(k-2),...,y(k-n),u(k-1),u(k-2),…,u(k-n)}。
步骤B1),构建基于变增量LP优化的航空发动机推力缓解神经网络控制系统结构图,包括NARMA-L2模型、内环神经网络控制器、发动机模型、性能状态滤波器、外环指令修正模块;
步骤B2),根据飞行条件及相关的控制计划确定发动机初始指令,利用发动机模型输出和限制条件设计基于变增量LP算法的外环指令修正模块,建立在发动机发生部件性能突变的情况下具有自调整能力的推力衰退缓解神经网络控制器。
步骤B2.1),对不可测推力和边界限制条件在工作点附近小范围进行线性化,提取性能指标、约束量和控制量之间的线性关系PSM,具体公式为:
其中,ΔSmf和ΔSmc分别为风扇和压气机的喘振裕度增量,ΔnL和ΔnH分别为低压、高压转子转速增量,ΔFn为推力增量,ΔTt6为涡轮后温度增量,ΔWf为燃油量增量,ΔA8为尾喷口面积增量,P矩阵元素可从发动机稳态模型中通过扰动法得出,具体公式为:
Pij={Mi(u0(j)+Δuj,...)-Mi(u0)}/Δuj
i=1,2,...,6,j=1,2
其中,Mi(·)是通过模型计算得到的输出矩阵;u0是控制量的初始值;Δuj是第j个控制量增量;
步骤B2.2),对于步骤B2.1)中的控制量增量Δuj,可由相应控制量uj乘以增量系数得出,增量系数的选取规则为:
其中,k为迭代次数;ΔFn为推力增量。根据步骤B2.1)中提取的PSM,将指令修正转换为LP优化问题,求解推力缓解所需指令的修正量。
步骤B2.2.1),在发动机不超转、不超温、喘振裕度和控制量满足范围要求、推力不超过期望值的约束下寻找发动机最大推力,具体公式为:
其中,Fn为发动机推力,Wf为燃油量,A8为尾喷口面积,u=[Wf,A8]T为发动机控制向量,umin和umax分别为向量下界和上界,nLmin和nLmax分别为低压转子转速最低和最高值,nHmin和nHmax分别为高压转子转速最低和最高值,Tt6min和Tt6max分别为涡轮后温度最低和最高值,Smf和Smc分别为风扇和压气机喘振裕度,Fnor为实时模型计算出未发生突变故障的发动机推力期望值;
步骤B2.2.2),根据控制量约束与限制约束,利用步骤B2.1)中提取的PSM,将步骤B2.2.1)中的公式改写为:
同理求得最大推力性能指标为:
ΔFn=p61ΔWf+p62ΔA8
步骤B2.2.3),求解LP问题,经反复迭代至满足精度要求或达到最大迭代次数限制后得到全局最优解Wf和A8,并将最优解输入到发动机中计算得到对应的输出信号nH和EPR,将二者与推力恢复前指令信号作差,得到的ΔnHr和ΔEPRr作为转速指令和整机落压比指令的修正值输入到内环神经网络控制器,内环神经网络控制器跟踪修正后的指令,缓解发动机整机推力损失。
为了验证本发明所设计的航空发动机推力衰退缓解的神经网络控制方法的有效性,在MATLAB环境下进行了一定包线内推力衰退缓解控制的数字仿真。
本发明采用某型双转子小涵道比涡扇发动机的非线性部件模型作为被控对象。该模型通过面向对象编程的方法构建,包括了进气道,风扇,压气机,燃烧室,涡轮和尾喷管等航空发动机的重要部件,并且很容易在MATLAB环境中调用。
在进行推力缓解控制之前,首先需要训练得到内环神经网络控制器。其中,控制器中的NARMA-L2模型的两个神经网络均为单隐含层,且输入层、隐含层和输出层的神经元个数分别为5、10和1。训练首先需要利用航空发动机模型在H=0,Ma=0的飞行条件下,通过随机的输入信号生成包含10000组训练样本的主燃油量-高压转子转速数据集,其中随机输入信号Wf的范围为:96.79%到100%,模型输出信号nH的范围为:99.42%到100%。训练完成后,为了进一步得到具有自适应能力的控制器,对神经网络参数进行在线修正,其中学习率α=0.0022,至此内环神经网络控制器设计完毕,其结构如图2所示。
在内环控制器的基础上加入外环指令修正,得到完整的推力衰退缓解控制器并针对该控制器进行仿真验证,其中外环指令修正的变增量LP算法结构如图3所示。仿真共有2组,包含包线内2个飞行条件,分别为H=0,Ma=0以及H=8000m,Ma=1.2,其中第1组仿真高压转子转速指令nHr初始值为1,在10s后发动机风扇发生性能突变,其效率降低5%,在20s后开始修正指令,并且在30s时仿真结束;第2组仿真nHr初始值为0.95,在10s后发动机风扇和压气机发生性能突变,二者流量和效率均降低1%,在20s后开始修正指令,同样在30s时仿真结束。
图4-6展示了H=0,Ma=0,nHr=1稳定工作点处风扇发生性能突变时所选变量变化情况。图中的3条线分别对应一下三种情况:①发动机正常工作情况下,常规的内环控制方式下主燃油、高压转子转速和推力的变化;②发动机发生性能突变的情况下,常规的内环控制方式下主燃油、高压转子转速和推力的变化;③发动机发生突变故障的情况下,推力缓解控制方式下主燃油、高压转子转速和推力的变化。
从图4中可以看出性能突变发生后,常规转速控制通过降低燃油量保证高压转子转速继续追踪原指令,这导致发动机推力发生衰退。本发明提出的推力衰退缓解控制方法在突变发生后先将nHr降低1%以确保发动机安全工作,并在第20s将修正后的指令输入神经网络控制器进行推力恢复,但从图5中可以看出,由于此时的高压转子转速最大值nHmax=1.018,导致在变增量LP优化算法寻优过程中nH达到边界而寻优停止。从图6中可以发现,虽然推力缓解控制下的推力并未恢复到无突变时的推力,但相比于常规控制,依然在保证发动机不超限的情况下使推力得到提升。
图7-9展示了H=8000,Ma=1.2,nHr=0.95稳定工作点处压气机发生性能突变时所选变量变化情况。随着高度和马赫数的提升,发动机特性相比于地面点已发生较大变化。得益于内环神经网络控制器的自适应特性,使得性能突变后控制器依然可以很好的跟踪转速指令,但此时推力大小相比于发动机正常工作时有所下降。当加入外环指令修正回路后,为了保证发动机安全工作,转速指令信号先下调1%,此时控制器很快追踪上新的指令信号,之后利用变增量LP算法得到推力缓解所需的指令信号增量,在第20s将迭代获得的ΔnHr和ΔEPRr同时作用于内环控制器,从图9中可以看出,本文设计的推力缓解控制方法相比于常规内环控制器可以有效地对推力进行恢复。
本发明针对航空发动机性能突变导致部件性能蜕化引起的推力下降的问题,提出了一种基于变增量LP优化的航空发动机推力衰退缓解的神经网络控制方法。该方法依据发动机高压转子转速和整机落压比指令,通过发动机多变量控制器计算出燃油流量和尾喷口喉道面积控制变量,其中内环控制通过神经网络学习的NARMA-L2模型逆实现;在多变量控制器基础上设计指令修正外环回路,利用变增量LP优化方法,并根据发动机部件性能突变引起的推力衰退程度,调整发动机指令信号以缓解整机推力损失。仿真结果验证了设计的部件性能突变条件下的推力缓解控制方法的有效性。
需要指出的是,以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化和替换,都应涵盖在本发明的保护范围内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (6)
1.一种航空发动机推力衰退缓解的神经网络控制方法,其特征在于,包括以下步骤:
步骤A)建立发动机高压转子转速、压比与尾喷管喉道面积之间的物理变量插值模型,以及利用主燃油量与高压转子转速之间的关系,采用神经网络建立性能衰退缓解的内环控制NARMA-L2模型;
步骤B)结合NARMA-L2模型,设计基于变增量LP算法的外环指令修正模块,开发在发动机发生部件性能突变的情况下具有自调整能力的推力衰退缓解神经网络控制器;
步骤C)发动机发生部件性能突变时,通过外环加入指令修正模块,对发动机指令进行修正,利用推力衰退缓解神经网络控制器追踪修正后的指令进而达到推力恢复的目的。
2.根据权利要求1所述的一种航空发动机推力衰退缓解的神经网络控制方法,其特征在于:所述步骤A)中利用主燃油量与高压转子转速之间的关系,采用神经网络建立性能衰退缓解的内环控制NARMA-L2模型具体步骤如下:
步骤A1),选取发动机主燃油量作为输入量u[k],高压转子转速作为输出量y[k],并根据所选输入输出量构建NARMA-L2模型为:
y[k]=f0(y[k-1],y[k-2],…,y[k-n],u[k-1],u[k-2],…,u[k-n])+g0(y[k-1],y[k-2],…,y[k-n],u[k-1],u[k-2],…,u[k-n])u[k]
其中,f0和g0是NARMA-L2模型中的两个非线性映射,通过训练神经网络组f0 *和g0 *近似得到,而转速控制律为:
其中,y*[k]为发动机期望输出;并根据发动机整机落压比指令、实际输出落压比和输出转速插值换算得到尾喷管喉道面积控制量;
步骤A2),利用期望输出和实际输出间的误差ec(k)对神经网络组f0 *和g0 *的拓扑结构参数进行在线修正:
W(k+1)=W(k)+αΔW(k)
V(k+1)=V(k)+αΔV(k)
其中,W和V分别为神经网络组的权值和阈值;α为学习率,其大小影响参数的调整速度,由多次试取法确定;ΔW(k)和ΔV(k)为雅可比矩阵;拓扑向量中的各元素在k时刻的更新如下:
wi(k+1)=wi(k)-αΔwi(k)
vi(k+1)=vi(k)-αΔvi(k)
其中,T={y(k),y(k-1),y(k-2),...,y(k-n),u(k-1),u(k-2),…,u(k-n)}。
3.根据权利要求1所述的一种航空发动机推力衰退缓解的神经网络控制方法,其特征在于:所述步骤B)中结合NARMA-L2模型,设计基于变增量LP算法的外环指令修正模块,得到在发动机发生部件性能突变的情况下具有自调整能力的推力衰退缓解神经网络控制器具体步骤如下:
步骤B1),构建基于变增量LP优化的航空发动机推力缓解神经网络控制系统结构图,包括内环神经网络控制器、发动机模型、性能状态滤波器、外环指令修正模块;
步骤B2),根据飞行条件及相关的控制计划确定发动机初始指令,利用发动机模型输出和限制条件设计基于变增量LP算法的外环指令修正模块,建立在发动机发生部件性能突变的情况下具有自调整能力的推力衰退缓解神经网络控制器。
4.根据权利要求3所述的一种航空发动机推力衰退缓解的神经网络控制方法,其特征在于:所述步骤B2)根据发动机模型输出和限制条件设计基于变增量LP算法的外环指令修正模块,建立在发动机发生部件性能突变的情况下具有自调整能力的推力衰退缓解神经网络控制器具体步骤如下:
步骤B2.1),对不可测推力和边界限制条件在工作点附近小范围进行线性化,提取性能指标、约束量和控制量之间的线性关系PSM,具体公式为:
其中,ΔSmf和ΔSmc分别为风扇和压气机的喘振裕度增量,ΔnL和ΔnH分别为低压、高压转子转速增量,ΔFn为推力增量,ΔTt6为涡轮后温度增量,ΔWf为燃油量增量,ΔA8为尾喷口面积增量,P矩阵元素可从发动机稳态模型中通过扰动法得出,具体公式为:
Pij={Mi(u0(j)+Δuj,...)-Mi(u0)}/Δuj
i=1,2,...,6,j=1,2
其中,Mi(·)是通过发动机稳态模型计算得到的输出矩阵;u0是控制量的初始值;Δuj是第j个控制量增量;
步骤B2.2),对于步骤B2.1)中的控制量增量Δuj,可由相应控制量uj乘以增量系数得出,增量系数的选取规则为:
其中,k为迭代次数;ΔFn为推力增量;根据步骤B2.1)中提取的PSM,将指令修正转换为LP优化问题,求解推力缓解所需指令的修正量。
5.根据权利要求4所述的一种航空发动机推力衰退缓解的神经网络控制方法,其特征在于:所述步骤B2.2)中将指令修正转换为LP优化问题,求解推力缓解所需指令的修正量具体步骤如下:
步骤B2.2.1),在发动机不超转、不超温、喘振裕度和控制量满足范围要求、推力不超过期望值的约束下寻找发动机最大推力,具体公式为:
其中,Fn为发动机推力,Wf为燃油量,A8为尾喷口面积,u=[Wf,A8]T为发动机控制向量,umin和umax分别为向量下界和上界,nLmin和nLmax分别为低压转子转速最低和最高值,nHmin和nHmax分别为高压转子转速最低和最高值,Tt6min和Tt6max分别为涡轮后温度最低和最高值,Smf和Smc分别为风扇和压气机喘振裕度,Fnor为实时模型计算出未发生突变故障的发动机推力期望值;
步骤B2.2.2),根据控制量约束与限制约束,利用步骤B2.1)中提取的PSM,将步骤B2.2.1)中的公式改写为:
同理求得最大推力性能指标为:
ΔFn=p61ΔWf+p62ΔA8
步骤B2.2.3),求解LP问题,经反复迭代至满足精度要求或达到最大迭代次数限制后得到全局最优解Wf和A8,并将最优解输入到发动机中计算得到对应的输出信号nH和EPR,将二者与推力恢复前指令信号作差,得到的ΔnHr和ΔEPRr作为转速指令和整机落压比指令的修正值输入到内环神经网络控制器,内环神经网络控制器跟踪修正后的指令,缓解发动机整机推力损失。
6.根据权利要求1所述的一种航空发动机推力衰退缓解的神经网络控制方法,其特征在于:步骤C)中,发动机发生部件性能突变时,先将发动机转速指令下降1%,然后加入外环修指令。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910664441.4A CN110513199B (zh) | 2019-07-23 | 2019-07-23 | 一种航空发动机推力衰退缓解的神经网络控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910664441.4A CN110513199B (zh) | 2019-07-23 | 2019-07-23 | 一种航空发动机推力衰退缓解的神经网络控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110513199A CN110513199A (zh) | 2019-11-29 |
CN110513199B true CN110513199B (zh) | 2021-06-22 |
Family
ID=68623326
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910664441.4A Active CN110513199B (zh) | 2019-07-23 | 2019-07-23 | 一种航空发动机推力衰退缓解的神经网络控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110513199B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111258296B (zh) * | 2020-01-16 | 2021-06-01 | 西安交通大学 | 航空发动机控制系统的漂移类故障检测方法 |
CN111255574A (zh) * | 2020-03-09 | 2020-06-09 | 南京航空航天大学 | 一种航空发动机进口畸变下推力衰退自主控制方法 |
CN111767610B (zh) * | 2020-05-22 | 2022-07-15 | 成都飞机工业(集团)有限责任公司 | 一种基于线化增量方程的飞机飞行性能计算方法 |
CN112286047B (zh) * | 2020-09-28 | 2022-04-08 | 南京航空航天大学 | 基于神经网络的narma-l2多变量控制方法 |
CN113027614A (zh) * | 2021-04-07 | 2021-06-25 | 西北工业大学 | 基于改进模拟退火算法的航空发动机最大推力控制优化方法 |
CN113969842B (zh) * | 2021-10-22 | 2022-11-18 | 清华大学 | 一种航空发动机推力高动态响应控制方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5091843A (en) * | 1988-12-20 | 1992-02-25 | Allied-Signal, Inc. | Nonlinear multivariable control system |
US8306791B2 (en) * | 2009-12-21 | 2012-11-06 | United Technologies Corporation | Method and system for modeling the performance of a gas turbine engine |
CN104389685B (zh) * | 2014-11-24 | 2016-02-24 | 西北工业大学 | 一种航空发动机自适应寿命延长控制系统的设计方法 |
US9879593B2 (en) * | 2015-12-03 | 2018-01-30 | GM Global Technology Operations LLC | System and method for adaptively learning values and controlling a turbocharger of an engine based on the values |
CN106055770B (zh) * | 2016-05-26 | 2019-02-26 | 南京航空航天大学 | 一种基于滑模理论的航空发动机气路故障诊断方法 |
US10309330B2 (en) * | 2016-10-27 | 2019-06-04 | Rolls-Royce Corporation | Model reference adaptive controller |
CN108416086B (zh) * | 2018-01-25 | 2021-07-16 | 大连理工大学 | 一种基于深度学习的航空发动机模型自适应修正方法 |
CN108829928B (zh) * | 2018-05-08 | 2020-10-27 | 南京航空航天大学 | 一种涡轴发动机自适应部件级仿真模型构建方法 |
CN109162813B (zh) * | 2018-10-09 | 2019-10-11 | 南京航空航天大学 | 一种基于迭代学习修正的航空发动机智能转速控制方法 |
CN109611217B (zh) * | 2018-11-07 | 2020-12-11 | 大连理工大学 | 一种航空发动机过渡态控制规律优化的设计方法 |
CN109630281B (zh) * | 2019-01-10 | 2020-12-11 | 大连理工大学 | 一种基于误差区间观测器的航空发动机主动容错控制方法 |
-
2019
- 2019-07-23 CN CN201910664441.4A patent/CN110513199B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110513199A (zh) | 2019-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110513199B (zh) | 一种航空发动机推力衰退缓解的神经网络控制方法 | |
CN109162813B (zh) | 一种基于迭代学习修正的航空发动机智能转速控制方法 | |
WO2019144337A1 (zh) | 一种基于深度学习算法的航空发动机全包线模型自适应修正方法 | |
CN108829928B (zh) | 一种涡轴发动机自适应部件级仿真模型构建方法 | |
CN111042928B (zh) | 一种基于动态神经网络的变循环发动机智能控制方法 | |
CN111255574A (zh) | 一种航空发动机进口畸变下推力衰退自主控制方法 | |
CN105676640B (zh) | 基于贝塞尔曲线的涡扇发动机加速过程控制律设计方法 | |
CN110221537A (zh) | 控制方法、装置及航空发动机限制保护控制方法、装置 | |
CN109871653B (zh) | 航空发动机数学模型部件特性修正方法 | |
CN110579962A (zh) | 基于神经网络的涡扇发动机推力预测方法及控制器 | |
Marchukov et al. | Improving of the working process of axial compressors of gas turbine engines by using an optimization method | |
CN111680357A (zh) | 一种变循环发动机机载实时模型的部件级无迭代构建方法 | |
CN112284752A (zh) | 一种基于改进状态跟踪滤波器的变循环发动机解析余度估计方法 | |
CN114237029B (zh) | 基于主动切换逻辑的涡扇发动机加速过程控制方法及装置 | |
CN110985216B (zh) | 一种含在线修正的航空发动机智能多变量控制方法 | |
CN113705120A (zh) | 一种调水工程泵后阀最优调控方案制定方法 | |
CN113642271A (zh) | 一种基于模型的航空发动机性能恢复控制方法及装置 | |
CN112947064A (zh) | 考虑气路部件故障的航空发动机最大推力控制优化方法 | |
CN114912191A (zh) | 一种导叶可调的涡扇发动机部件级机载动态模型设计方法 | |
CN113759727A (zh) | 航空发动机多变量控制器的综合优化设计方法 | |
CN112943479A (zh) | 基于改进单纯形法的航空发动机加速过程最优控制方法 | |
CN111734533A (zh) | 一种基于涡扇发动机的模型预测方法及系统 | |
CN112943453A (zh) | 气路部件故障下基于iga的发动机最大推力控制优化方法 | |
CN111856929A (zh) | 航空发动机容错增益调度二自由度h∞控制器 | |
CN112523874B (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 |