CN108829928B - 一种涡轴发动机自适应部件级仿真模型构建方法 - Google Patents

一种涡轴发动机自适应部件级仿真模型构建方法 Download PDF

Info

Publication number
CN108829928B
CN108829928B CN201810432419.2A CN201810432419A CN108829928B CN 108829928 B CN108829928 B CN 108829928B CN 201810432419 A CN201810432419 A CN 201810432419A CN 108829928 B CN108829928 B CN 108829928B
Authority
CN
China
Prior art keywords
component
parameters
engine
model
performance characteristic
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
Application number
CN201810432419.2A
Other languages
English (en)
Other versions
CN108829928A (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201810432419.2A priority Critical patent/CN108829928B/zh
Publication of CN108829928A publication Critical patent/CN108829928A/zh
Application granted granted Critical
Publication of CN108829928B publication Critical patent/CN108829928B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开了一种涡轴发动机自适应部件级仿真模型构建方法,包括以下步骤:步骤A,根据涡轴发动机各部件气动热力学特性,获得涡扇发动机各个工作截面的参数,建立涡轴发动机慢车以上状态的非线性部件级动态通用模型;步骤B,设计容积跟踪滤波器,估算得到的非线性部件级动态通用模型中的涡轴发动机气路部件不可测的性能特征参数;气路部包括压气机、燃气涡轮、动力涡轮;步骤C,利用估算出的性能特征参数自动修正各气路部件的流量和效率特性图,将调整后的气路部件特性参数用于部件气动热力参数的计算,得到修正后的模型输出数据,建立慢车以上状态的自适应部件级仿真模型。本发明解决了涡轴发动机个体差异与性能蜕化引起的模型失配问题。

Description

一种涡轴发动机自适应部件级仿真模型构建方法
技术领域
本发明涉及航空发动机建模与仿真领域,尤其涉及一种涡轴发动机自适应部件级仿真模型构建方法。
背景技术
涡轴发动机工作在高温、高速、强振动、大应力的恶劣环境下,发动机工作状态经常变化,对安全性和可靠性要求较高。发动机自适应模型建模技术的研究一直是一项重要课题,发动机自适应模型能反映发动机之间的差异以及使用期内性能蜕化等因素对发动机性能的影响,是发动机自适应调整控制、保证发动机正常工作的基础。同时发动机控制系统传感器故障诊断、隔离和容错控制等技术的实现也需要准确的发动机模型作为前提条件,所以建立精确的自适应模型具有重要的理论意义和工程实用价值。对于基于模型的发动机控制和估计系统而言,考虑到发动机之间存在着个体差异、实际发动机零件的安装公差及使用期内的性能蜕化等因素的影响,对应的模型若不加以适当的调整,则不能满足在线性能寻优控制或故障诊断的精度要求,所设计的控制和估计系统性能出现不同程度的下降,无法达到设计的工作状态。
涡轴发动机的主流仿真模型有两种:非线性部件级模型和线性模型。发动机线性模型实在发动机非线性部件级模型的基础上,对模型进行局部线性化,建立状态变量模型和稳态基点模型,利用线性卡尔曼滤波器实现部件性能参数的估计与自适应。线性模型计算量较小、对资源消耗低,但这一方法在对非线性模型线性化时,不可避免地会引入二次建模误差,且线性模型对于发动机动态过程的拟合精度较低。发动机非线性部件级模型建模方法主要有转子动力学法和容积动力学法。相对于发动机线性模型,非线性部件级模型不会引入二次建模误差,对于发动机的动态过程具有较高的跟踪精度,能较准确的模拟包线内涡轴发动机的不同工况。非线性部件级模型是依据试车数据和各部件热力学原理建立的通用模型,没有考虑发动机个体差异和寿命期内性能蜕化引起的模型误差,当非线性模型适配发动机个体时,需要利用经验手动调整发动机气路部件特性图,这将带来巨大工作量。随着滤波估计技术的发展,一些非线性卡尔曼滤波器可以直接应用于非线性系统,避免模型的线性化,实现较为准确的状态估计。
发明内容
本发明的目的是提供一种涡轴发动机自适应部件级仿真模型构建方法,以解决涡轴发动机个体差异与性能蜕化引起的模型失配问题。
为实现上述目的,本发明采用的技术方案为:
一种涡轴发动机自适应部件级仿真模型构建方法,包括以下步骤:
步骤A,根据涡轴发动机各部件气动热力学特性,获得涡扇发动机各个工作截面的参数,建立涡轴发动机慢车以上状态的非线性部件级动态通用模型;
步骤B,设计容积跟踪滤波器,估算步骤A得到的非线性部件级动态通用模型中的涡轴发动机气路部件不可测的性能特征参数;气路部包括压气机、燃气涡轮、动力涡轮;
步骤C,利用步骤B估算出的性能特征参数自动修正各气路部件的流量和效率特性图,将调整后的气路部件特性参数用于部件气动热力参数的计算,得到修正后的模型输出数据,建立慢车以上状态的自适应部件级仿真模型。
所述步骤A的具体步骤如下:
步骤A1,根据涡轴发动机各部件气动热力学特性、设计点参数以及试车数据建立慢车以上状态的涡轴发动机各部件的数学模型,再根据流量连续、功率平衡及转子动力学原理建立各部件之间的共同工作方程,最后使用非线性方程数值解法迭代求解,获得发动机各个工作截面的参数,实现慢车以上状态的涡轴发动机非线性部件级动态通用模型的建立;并引入发动机气路部件性能特征参数来表征发动机个体性能差异或者使用时间带来的性能蜕化,气路部件性能特征参数选取旋转部件的效率系数SEi和流量系数SWi,定义如下
Figure BDA0001653745740000021
式中:ηi,wi为部件的实际效率和流量,而
Figure BDA0001653745740000022
为部件效率和流量的理想值;
步骤A2,选择需要使用的发动机模型工作截面的传感器量测参数,包括:动力涡轮转速NP,燃气涡轮转速NG,压气机出口总温T3,压气机出口总压P3,动力涡轮入口总温T45,动力涡轮入口总压P45,动力涡轮出口总温T5,动力涡轮出口总压P5
所述步骤B的具体步骤如下:
步骤B1,将步骤A中获得的各个工作截面的参数进行相似归一化处理;
步骤B2,利用容积跟踪滤波器估计涡轴发动机气路部件不可测的性能特征参数,获得发动机气路部件性能调整的具体数值。
所述步骤B2中利用容积跟踪滤波器计算不可测的气路部件性能特征参数的详细步骤如下:
步骤B2.1,初始化性能特征参数向量的后验估计值和后验方差矩阵。
步骤B2.2,根据上一时刻的性能特征参数后验估计和后验方差生成此时刻的性能特征参数容积点,调用非线性部件级动态通用模型并对各个性能特征参数容积点进行状态更新,计算性能特征参数一步预测的先验估计和先验方差;
步骤B2.3,根据性能特征参数先验估计和先验方差选取新的性能特征参数容积点,调用非线性部件级动态通用模型并对性能特征参数容积点进行量测更新,计算自相关矩阵和互相关矩阵,进而获得卡尔曼增益矩阵;性能特征参数容积点值加权求和可得出此时刻的先验量测,计算得到这一时刻的性能特征参数的后验估计值和后验方差矩阵;
步骤B2.4,以后时刻重复进行步骤B2.2和步骤B2.3完成性能特征参数的递推估计。
所述步骤C的具体步骤如下:
步骤C1,将所得的性能特征参数中各气路部件的效率、流量系数,输入至发动机部件级模型的对应部件中,更新气路部件的流量、效率特性图;在同样的折合转速下,保持各转子部件特性图曲线的压比坐标数值不变,将特性图中效率、流量的数值沿坐标轴方向进行缩放修正,对应的性能特征参数即为特性图中效率、流量曲线缩放的比例;
步骤C2,将调整后的气路部件特性参数用于部件气动热力参数的计算,进行部件特性图调整后的非线性部件级模型各截面参数的计算,得到修正后的模型输出信号,建立慢车以上状态的自适应部件级仿真模型。
所述步骤C1中,
对于压气机部件,流量-压比图中,特性曲线沿x轴方向进行缩放,缩放比率为对应的压气机的流量性能特征参数;流量-效率图中,首先特性曲线沿x轴方向进行缩放,缩放比率为对应的压气机的流量性能特征参数,然后曲线沿y轴方向进行缩放,缩放比率为对应的压气机的效率性能特征参数;
对于燃气涡轮、动力涡轮部件,效率-压比图中,特性曲线沿y轴方向进行缩放,缩放比率为对应的燃气涡轮、动力涡轮的效率性能特征参数;流量-压比图中,特性曲线沿y轴方向进行缩放,缩放比率为对应的燃气涡轮、动力涡轮的流量性能特征参数。
有益效果:本发明采用以上技术方案与现有技术相比,具有以下技术效果:
(1)本发明提出的涡轴发动机自适应部件级仿真模型构建方法,直接使用非线性部件级模型实现仿真模型的自适应,对实际发动机个体动态过程的输出跟踪精度高,不会因为模型的线性化过程而引入二次建模误差;
(2)本发明提出的涡轴发动机自适应部件级仿真模型构建方法,能减少因涡轴发动机个体差异和寿命期内性能蜕化引起的模型误差,具有较高的模型精度,减少目前利用经验手动调整发动机气路部件特性图所带来的巨大工作量;
(3)本发明设计的涡轴发动机自适应部件级仿真模型,能获得涡轴发动机气路部件性能特征变化情况,为涡轴发动机视情维修提供性能参考依据。
附图说明
图1是涡轴发动机气路工作截面标识图;
图2是涡轴发动机自适应部件级仿真模型原理图;
图3是容积跟踪滤波器计算流程图;
图4a-图4e是地面设计点模拟压气机性能变化的气路性能估计与涡轴发动机自适应部件级仿真模型的跟踪效果图,以及压气机部件特性图的修正;
图5a-图5e是地面设计点模拟动力涡轮性能变化的气路性能估计与涡轴发动机自适应部件级仿真模型的跟踪效果图,以及动力涡轮部件特性图的修正;
图6a-图6b是地面动态过程中,涡轴发动机模型与涡轴发动机自适应部件级仿真模型输入参数;
图7a-图7c是地面动态过程中,涡轴发动机自适应部件级仿真模型的气路性能估计结果与模型输出的跟踪效果图;
图8a-图8b是包线内高空动态过程中,涡轴发动机模型与涡轴发动机自适应部件级仿真模型输入参数;
图9a-图9c是包线内高空动态过程中,涡轴发动机自适应部件级仿真模型的气路性能估计结果与模型输出的跟踪效果图。
具体实施方式
本发明针对先进航空发动机基于模型的多变量控制与预测健康管理的需求,对现有航空发动机仿真模型进行扩展和设计开发,建立慢车以上状态自适应部件级仿真模型模型,能减少发动机个体差异和性能蜕化引起的模型误差,保证发动机模型具有较高置信度。
下面结合具体实施例及附图对本发明做更进一步的解释。
实施例
本实施例以某型涡轴发动机自适应部件级仿真模型构建为例,图2是涡轴发动机自适应部件级仿真模型原理图,该仿真模型的建立包括以下步骤:
步骤A,根据涡轴发动机各部件气动热力学特性,获得涡扇发动机各个工作截面的参数,建立涡轴发动机慢车以上状态的非线性部件级动态通用模型;的详细步骤如下:
步骤A1,某型涡轴发动机各部件及工作截面如图1所示,根据涡轴发动机部件特性、设计点参数以及试车数据建立发动机部件级模型,包括进气道、压气机、燃烧室、燃气涡轮、动力涡轮和喷管,再根据流量连续、功率平衡及转子动力学等原理建立各部件之间的共同工作方程,最后使用非线性方程数值解法迭代求解,获得发动机各个工作截面的参数。该部件法建模业内比较成熟,在此不加详述。发动机部件级模型是根据部件热力学特性和试车数据等得到的平均模型,不能准确的反映同型不同个体发动机的输出,同时随着发动机服役时间的增加,气路部件的性能也会发生不同程度的蜕化。因此,引入发动机气路部件性能特征参数来表征发动机个体性能差异或者使用时间带来的性能蜕化,气路部件性能特征参数选取旋转部件的效率系数SEi和流量系数SWi,定义如下
Figure BDA0001653745740000051
式中:下标i(i=1,2,3)表示发动机各部件编号。本发明用例的发动机共有三个旋转部件,压气机效率和流量系数为SE1,SW1,燃气涡轮效率和流量系数为SE2,SW2,动力涡轮的效率和流量系数为SE3,SW3,ηi,wi为部件的实际效率和流量,而
Figure BDA0001653745740000052
为部件效率和流量的理想值,将健康参数向量h定义为h=[SE1,SW1,SE2,SW2,SE3,SW3,SE4,SW4]T
步骤A2,考虑容积跟踪滤波器是利用发动机测量参数与模型输出量之间的残差来实现发动机部件性能特征参数的估计,因此需合理选取发动机模型输出参数。所选取的该发动机模型传感器包括:动力涡轮转速NP,燃气涡轮转速NG,压气机出口总温T3,压气机出口总压P3,动力涡轮入口总温T45,动力涡轮入口总压P45,动力涡轮出口总温T5,动力涡轮出口总压P5
步骤B,设计容积跟踪滤波器,估算涡轴发动机气路部件不可测的性能特征参数;具体步骤为:气路部件包括压气机、燃气涡轮、动力涡轮;
步骤B1,不同量测参数具有不同的物理意义,彼此的数量级相差很大,这将会带来矩阵的计算和数据存储的问题。因此,根据发动机相似准则,对步骤A中获得的各个工作截面的参数做相似归一化处理。参数相似归一化过程如下:
Figure BDA0001653745740000061
式中,下标ds表示涡轴发动机设计点参数,T2、P2为压气机进口总温和总压,N′P,N′G,T′3,P′3,T′45,P′45,T′5,P′5为对应参数的相似归一化后的值。
步骤B2,假设涡轴发动机部件级非线性数学模型如下:
Figure BDA0001653745740000062
式中,f()为涡轴发动机状态转移方程,h()为涡轴发动机量测方程,k为时间参数,ωk和νk分别为系统独立的系统噪声和量测噪声,且满足ωk~N(0,Q2),vk~N(0,R2),Q、R分别为噪声的协方差矩阵,选取Q=0.0015×I8×8,R=0.0015×I8×8。xk代表系统的状态变量,uk为系统的输入量,yk为系统的传感器量测值,变量选择为xk=[N′P,N′G,hT]T,uk=[Wfα]T,y=[N′P,N′G,T′3,P′3,T′45,P′45,T′5,P′5]T,其中Wf为燃油流量,α为总矩角。
zk为飞行条件参数向量,包含飞行高度H、马赫数Ma和进口温度T1
步骤B.2.1,初始化状态量的后验估计值
Figure BDA0001653745740000063
和后验方差矩阵P0|0
步骤B.2.2,根据Cubature准则,计算状态容积点集(Xi,k-1|k-1i),Xi,k-1|k-1为容积点的值,ωi为对应容积点的权值,具体计算公式为:
Figure BDA0001653745740000071
式中,N为待估的状态量x的维数,Sk-1|k-1=chol(Pk-1|k-1),Pk-1|k-1为前一时刻的后验估计方差阵,chol()表示对矩阵进行Cholesky分解,即
Figure BDA0001653745740000072
[1]i为集合[1]的第i列,以N=3为例,
[1]={[1,0,0]T,[0,1,0]T,[0,0,1]T,[-1,0,0]T,[0,-1,0]T,[0,0,-1]T}。
调用通用模型并对各个容积点进行状态更新,通过已更新的容积点计算一步预测的先验估计
Figure BDA0001653745740000073
和先验方差Pk|k-1,完成时间更新过程,计算公式为:
Figure BDA0001653745740000074
式中,f()即为式(3)中发动机状态转移方程。
步骤B2.3,根据一步预测的先验估计和先验方差选取新的容积点Xi,k|k-1,计算公式为:
Figure BDA0001653745740000075
式中,Sk-1|k-1=chol(Pk|k-1)。
调用通用模型并对状态量容积点进行量测更新,可得出此时刻的先验量测
Figure BDA0001653745740000076
根据容积点的更新值计算自相关矩阵Pyy,k|k-1和互相关矩阵Pxy,k|k-1,进而获得卡尔曼增益矩阵Kk,计算公式为:
Figure BDA0001653745740000081
由此时的传感器量测数据与先验量测值的差可得出先验估计的残差,并计算得到这一时刻状态量的后验估计值
Figure BDA0001653745740000082
和后验方差矩阵Pk|k,具体计算过程如下:
Figure BDA0001653745740000083
步骤B2.4,以后时刻重复进行步骤B2.2)和步骤B2.3)完成性能特征参数的递推估计。
步骤C,利用估计出的性能特征参数自动修正各气路部件的流量和效率特性图,将调整后的气路部件特性参数用于部件气动热力参数的计算,得到修正后的模型输出数据,建立慢车以上状态的自适应部件级仿真模型。详细步骤如下:
将所得的性能特征参数中各气路部件的效率、流量系数输入至发动机部件级模型的对应部件中,对原有通用模型的部件特性图的流量、效率曲线进行缩放修正,估计出的性能特征参数中各气路部件的效率、流量系数作为各气路部件特性图中效率、流量数值的缩放系数。具体计算过程如下:
Figure BDA0001653745740000084
式中,SE′i,SW′i为估计出的性能特征参数中各转子部件的效率、流量系数,η′i,w′i为部件的调整后的效率和流量。在同样的折合转速下,保持各转子部件特性图曲线的压比坐标数值不变,将特性图中效率、流量曲线坐标数值沿坐标轴方向进行缩放修正。
对于压气机部件,流量-压比图的特性曲线沿x轴方向进行缩放,缩放比率为对应估计出的压气机的流量性能特征参数SW′1;流量-效率图中,首先特性曲线沿x轴方向进行缩放,缩放比率为对应的压气机的流量性能特征参数SW′1,然后曲线沿y轴方向进行缩放,缩放比率为对应的压气机的效率性能特征参数SE′1
对于燃气涡轮、动力涡轮部件,压比-效率图的特性曲线沿y轴方向进行缩放,缩放比率为对应的燃气涡轮、动力涡轮的效率性能特征参数SE′2,SE′3;压比-流量图的特性曲线沿y轴方向进行缩放,缩放比率为对应的燃气涡轮、动力涡轮的流量性能特征参数SW′2,SW′3
将调整后的气路部件特性参数图用于部件气动热力参数的计算,进行部件特性图调整后的非线性部件级模型的计算,得到调整后的模型输出参数,实现慢车以上状态的自适应部件级仿真模型的建立。
为了验证本发明所设计的涡轴发动机自适应部件级仿真模型构建方法的有效性,在MATLAB环境下进行了如下数字仿真。
在涡轴发动机地面点处H=0m,Ma=0,Wf=0.05kg/s,α=2.6°,图4a,4b,4c给出了模拟发动机压气机效率下降3%、流量下降1%时,涡轴发动机自适应部件级仿真模型的输出参数跟踪结果和部件特性修正参数的估计结果(此处只给出了高低压转子转速的跟踪结果),涡轴发动机自适应部件级仿真模型能够很好的跟踪上发动机个体模型的输出。图4d,4e给出了在此性能变化下,压气机部件特性图的修正示意图(以SE′1=0.97,SW′1=0.99为例)。在同样的折合转速下,保持各转子部件特性图曲线的压比坐标数值不变,将特性图中效率、流量数值沿坐标轴方向进行缩放修正。在压气机部件的流量-效率图中,x轴方向缩放比值为0.99,y轴方向缩放比值为0.97。在压气机部件的流量-压比图中,特性图曲线进行x轴方向上的缩放变化,缩放比值为0.99。
在涡轴发动机地面点处,模拟发动机动力涡轮性能变化(效率下降2%、流量上升1%)时,涡轴发动机自适应部件级仿真模型的输出参数跟踪结果和部件特性修正参数的估计结果如图5a,5b,5c所示,涡轴发动机自适应部件级仿真模型能够很好的跟踪上发动机个体模型的输出。图5d,5e给出了在此性能变化下,动力涡轮部件特性参数的修正示意图(以SE′3=0.98,SW′3=1.01为例)。在同样的折合转速下,保持动力涡轮部件特性图曲线的压比坐标数值不变,将特性图中效率、流量曲线数值沿坐标轴方向进行缩放修正。在动力涡轮的压比-效率图中,y轴方向缩放比值为0.98。在动力涡轮的流量-压比图中,特性图曲线进行y轴方向上的缩放变化,缩放比值为1.01。
为了验证涡轴发动机自适应部件级仿真模型对发动机动态过程的跟踪精度,在地面点处对涡轴发动机个体模型和涡轴发动机自适应部件级仿真模型做如图6a,6b所示的Wf、α变化过程,同时模拟压气机效率下降3%,燃气涡轮、动力涡轮转速和部件特性修正参数的仿真结果如图7a-7c所示。仿真结果表明,在模拟该动态过程中,涡轴发动机自适应部件级仿真模型能够很好的跟踪上发动机模型的输出,模型最大误差不超过0.8%。
为了验证包线内不同工作点的模型跟踪精度,在高空状态(H=1000m,Ma=0)对发动机个体模型和涡轴发动机自适应部件级仿真模型做如图8a,8b所示的Wf、α变化过程,同时模拟压气机效率下降3%的性能变化,燃气涡轮、动力涡轮转速和部件特性修正参数的仿真结果如图9a-9c所示。仿真结果表明,在模拟该动态过程中,涡轴发动机自适应部件级仿真模型能够很好的跟踪上发动机模型的输出,模型最大误差不超过0.7%。可以看出在不同飞行状态的动态过程中,涡轴发动机自适应部件级仿真模型均能准确的估计出部件特性参数,使模型的输出具有较高的精度。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (4)

1.一种涡轴发动机自适应部件级仿真模型构建方法,其特征在于:包括以下步骤:
步骤A,根据涡轴发动机各部件气动热力学特性,获得涡轴发动机各个工作截面的参数,建立涡轴发动机慢车以上状态的非线性部件级动态通用模型;
步骤B,设计容积跟踪滤波器,估算步骤A得到的非线性部件级动态通用模型中的涡轴发动机气路部件不可测的性能特征参数;气路部包括压气机、燃气涡轮、动力涡轮;利用容积跟踪滤波器计算不可测的气路部件性能特征参数的详细步骤如下:
步骤B2.1,初始化性能特征参数向量的后验估计值和后验方差矩阵;
步骤B2.2,根据上一时刻的性能特征参数后验估计和后验方差生成此时刻的性能特征参数容积点,调用非线性部件级动态通用模型并对各个性能特征参数容积点进行状态更新,计算性能特征参数一步预测的先验估计和先验方差;
步骤B2.3,根据性能特征参数先验估计和先验方差选取新的性能特征参数容积点,调用非线性部件级动态通用模型并对性能特征参数容积点进行量测更新,计算自相关矩阵和互相关矩阵,进而获得卡尔曼增益矩阵;性能特征参数容积点值加权求和可得出此时刻的先验量测,计算得到这一时刻的性能特征参数的后验估计值和后验方差矩阵;
步骤B2.4,以后时刻重复进行步骤B2.2和步骤B2.3完成性能特征参数的递推估计;
步骤C,利用步骤B估算出的性能特征参数自动修正各气路部件的流量和效率特性图,将调整后的气路部件特性参数用于部件气动热力参数的计算,得到修正后的模型输出数据,建立慢车以上状态的自适应部件级仿真模型;具体步骤如下:
步骤C1,将所得的性能特征参数中各气路部件的效率、流量系数,输入至发动机部件级模型的对应部件中,更新气路部件的流量、效率特性图;在同样的折合转速下,保持各转子部件特性图曲线的压比坐标数值不变,将特性图中效率、流量的数值沿坐标轴方向进行缩放修正,对应的性能特征参数即为特性图中效率、流量曲线缩放的比例;
对于压气机部件,流量-压比图中,特性曲线沿x轴方向进行缩放,缩放比率为对应的压气机的流量性能特征参数;流量-效率图中,首先特性曲线沿x轴方向进行缩放,缩放比率为对应的压气机的流量性能特征参数,然后曲线沿y轴方向进行缩放,缩放比率为对应的压气机的效率性能特征参数;
对于燃气涡轮、动力涡轮部件,效率-压比图中,特性曲线沿y轴方向进行缩放,缩放比率为对应的燃气涡轮、动力涡轮的效率性能特征参数;流量-压比图中,特性曲线沿y轴方向进行缩放,缩放比率为对应的燃气涡轮、动力涡轮的流量性能特征参数;
步骤C2,将调整后的气路部件特性参数用于部件气动热力参数的计算,进行部件特性图调整后的非线性部件级模型各截面参数的计算,得到修正后的模型输出信号,建立慢车以上状态的自适应部件级仿真模型。
2.根据权利要求1所述的涡轴发动机自适应部件级仿真模型构建方法,其特征在于:所述步骤A的具体步骤如下:
步骤A1,根据涡轴发动机各部件气动热力学特性、设计点参数以及试车数据建立慢车以上状态的涡轴发动机各部件的数学模型,再根据流量连续、功率平衡及转子动力学原理建立各部件之间的共同工作方程,最后使用非线性方程数值解法迭代求解,获得发动机各个工作截面的参数,实现慢车以上状态的涡轴发动机非线性部件级动态通用模型的建立;并引入发动机气路部件性能特征参数来表征发动机个体性能差异或者使用时间带来的性能蜕化,气路部件性能特征参数选取旋转部件的效率系数SEi和流量系数SWi,定义如下
Figure FDA0002572397500000021
式中:ηi,wi为部件的实际效率和流量,而
Figure FDA0002572397500000022
为部件效率和流量的理想值;
步骤A2,选择需要使用的发动机模型工作截面的传感器量测参数,包括:动力涡轮转速NP,燃气涡轮转速NG,压气机出口总温T3,压气机出口总压P3,动力涡轮入口总温T45,动力涡轮入口总压P45,动力涡轮出口总温T5,动力涡轮出口总压P5
3.根据权利要求1所述的涡轴发动机自适应部件级仿真模型构建方法,其特征在于:所述步骤B的具体步骤如下:
步骤B1,将步骤A中获得的各个工作截面的参数进行相似归一化处理;
步骤B2,利用容积跟踪滤波器估计涡轴发动机气路部件不可测的性能特征参数,获得发动机气路部件性能调整的具体数值。
4.根据权利要求3所述的涡轴发动机自适应部件级仿真模型构建方法,其特征在于:所述步骤B1中,相似归一化过程如下:
Figure FDA0002572397500000031
Figure FDA0002572397500000032
Figure FDA0002572397500000033
Figure FDA0002572397500000034
式中,下标ds表示涡轴发动机设计点参数,T2、P2为压气机进口总温和总压,NP为动力涡轮转速,NG为燃气涡轮转速,T3为压气机出口总温,P3为压气机出口总压,T45为动力涡轮入口总温,P45为动力涡轮入口总压,T5为动力涡轮出口总温,P5动为力涡轮出口总压,N′P,N′G,T′3,P′3,T′45,P′45,T′5,P′5为对应参数的相似归一化后的值。
CN201810432419.2A 2018-05-08 2018-05-08 一种涡轴发动机自适应部件级仿真模型构建方法 Active CN108829928B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810432419.2A CN108829928B (zh) 2018-05-08 2018-05-08 一种涡轴发动机自适应部件级仿真模型构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810432419.2A CN108829928B (zh) 2018-05-08 2018-05-08 一种涡轴发动机自适应部件级仿真模型构建方法

Publications (2)

Publication Number Publication Date
CN108829928A CN108829928A (zh) 2018-11-16
CN108829928B true CN108829928B (zh) 2020-10-27

Family

ID=64148418

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810432419.2A Active CN108829928B (zh) 2018-05-08 2018-05-08 一种涡轴发动机自适应部件级仿真模型构建方法

Country Status (1)

Country Link
CN (1) CN108829928B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109460628B (zh) * 2018-12-14 2023-05-26 中国航空工业集团公司西安飞机设计研究所 一种进气道与发动机共同工作的流量匹配评估方法
CN109614722B (zh) * 2018-12-14 2023-12-29 中国人民解放军陆军航空兵学院 基于模糊逻辑的涡轴发动机全状态参数的建模方法
CN109670244B (zh) * 2018-12-20 2023-06-09 中国人民解放军海军航空大学 一种涡轴发动机翻修后燃气涡轮导向器面积调整方法
CN109871653B (zh) * 2019-03-15 2020-09-01 南京航空航天大学 航空发动机数学模型部件特性修正方法
CN110118128B (zh) * 2019-05-28 2020-07-24 南京航空航天大学 微型燃气轮机传感器故障诊断和容错控制方法
CN110348078B (zh) * 2019-06-25 2022-04-08 南京航空航天大学 一种涡轴发动机容积动力学结合热惯性效应的建模方法
CN110513199B (zh) * 2019-07-23 2021-06-22 南京航空航天大学 一种航空发动机推力衰退缓解的神经网络控制方法
CN110442956B (zh) * 2019-07-31 2023-01-17 中国航发沈阳发动机研究所 一种燃气轮机部件级仿真方法
CN111339686A (zh) * 2020-03-26 2020-06-26 西北工业大学 一种基于试验数据的涡轴发动机逆向建模方法
CN114169084B (zh) * 2021-10-20 2023-09-05 中国航发四川燃气涡轮研究院 基于核心机特性试验数据修正核心机性能数学模型的方法
CN116577993B (zh) * 2023-06-19 2024-01-26 南京航空航天大学 涡轴发动机神经网络自适应控制方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104102769A (zh) * 2014-06-27 2014-10-15 南京航空航天大学 基于人工智能的涡轴发动机实时部件级模型建立方法
CN106126902A (zh) * 2016-06-21 2016-11-16 上海电机学院 基于粒子群优化算法辨识的燃气轮机部件特性线修正方法
GB2547704A (en) * 2016-02-29 2017-08-30 Delphi Automotive Systems Lux Method to determine the intake air temperature in an engine system
RU2649715C1 (ru) * 2016-12-06 2018-04-04 Федеральное государственное унитарное предприятие "Центральный институт авиационного моторостроения имени П.И. Баранова" Способ полетной диагностики авиационного турбореактивного двухконтурного двигателя со смешением потоков

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105389427B (zh) * 2015-10-28 2018-08-03 南京航空航天大学 基于自适应粒子滤波的航空发动机气路部件故障检测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104102769A (zh) * 2014-06-27 2014-10-15 南京航空航天大学 基于人工智能的涡轴发动机实时部件级模型建立方法
GB2547704A (en) * 2016-02-29 2017-08-30 Delphi Automotive Systems Lux Method to determine the intake air temperature in an engine system
CN106126902A (zh) * 2016-06-21 2016-11-16 上海电机学院 基于粒子群优化算法辨识的燃气轮机部件特性线修正方法
RU2649715C1 (ru) * 2016-12-06 2018-04-04 Федеральное государственное унитарное предприятие "Центральный институт авиационного моторостроения имени П.И. Баранова" Способ полетной диагностики авиационного турбореактивного двухконтурного двигателя со смешением потоков

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于典型特性图的燃气涡轮系统建模仿真研究;韦福;《计算机仿真》;20070131;第24卷(第1期);58-62 *
燃气轮机加级方案研究及气动设计;许天帮;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑(月刊)》;20120515(第 05 期);C039-176 *

Also Published As

Publication number Publication date
CN108829928A (zh) 2018-11-16

Similar Documents

Publication Publication Date Title
CN108829928B (zh) 一种涡轴发动机自适应部件级仿真模型构建方法
CN108647428B (zh) 一种涡扇发动机自适应部件级仿真模型构建方法
CN108416086B (zh) 一种基于深度学习的航空发动机模型自适应修正方法
WO2019144337A1 (zh) 一种基于深度学习算法的航空发动机全包线模型自适应修正方法
CN110502840B (zh) 航空发动机气路参数在线预测方法
CN106874569B (zh) 一种变几何分轴式燃气轮机实时仿真建模方法
CN109162813B (zh) 一种基于迭代学习修正的航空发动机智能转速控制方法
CN109612738B (zh) 一种涡扇发动机的气路性能改进的分布式滤波估计方法
EP1103926A2 (en) Methods and apparatus for model-based diagnostics
CN109031951B (zh) 基于精确偏导数的航空发动机状态变量模型在线建立方法
CN108828947B (zh) 一种航空发动机含时滞的不确定性模糊动态模型建模方法
CN109800449B (zh) 一种基于神经网络的航空发动机压缩部件特性修正方法
CN110647052A (zh) 一种变循环发动机模式切换自适应身份证模型构建方法
CN109871653B (zh) 航空发动机数学模型部件特性修正方法
CN105389427A (zh) 基于自适应粒子滤波的航空发动机气路部件故障检测方法
CN111880403A (zh) 航空发动机最大推力状态容错二自由度μ控制器
CN111859555A (zh) 输入受限的航空发动机最大推力状态鲁棒容错控制器
CN112284752A (zh) 一种基于改进状态跟踪滤波器的变循环发动机解析余度估计方法
CN111666648A (zh) 一种航空发动机动态特性模拟方法
CN115903484A (zh) 基于合作博弈的航空发动机多变量鲁棒控制器优化方法
CN114154234A (zh) 一种航空发动机建模方法、系统、存储介质
Chen et al. Aero-engine real-time models and their applications
JP5845705B2 (ja) ガスタービン性能推定装置
CN112257256A (zh) 一种基于稳态数据的发动机简化动态模型设计方法
Huang et al. Gas path deterioration observation based on stochastic dynamics for reliability assessment of aeroengines

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