CN113158564B - 基于数据驱动的航空发动机状态变量模型建立方法 - Google Patents
基于数据驱动的航空发动机状态变量模型建立方法 Download PDFInfo
- Publication number
- CN113158564B CN113158564B CN202110419235.4A CN202110419235A CN113158564B CN 113158564 B CN113158564 B CN 113158564B CN 202110419235 A CN202110419235 A CN 202110419235A CN 113158564 B CN113158564 B CN 113158564B
- Authority
- CN
- China
- Prior art keywords
- state variable
- model
- aeroengine
- data
- variable model
- 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
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000012549 training Methods 0.000 claims abstract description 13
- 238000013178 mathematical model Methods 0.000 claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 238000003062 neural network model Methods 0.000 claims abstract description 7
- 238000001914 filtration Methods 0.000 claims description 9
- 230000006870 function Effects 0.000 claims description 7
- 238000013528 artificial neural network Methods 0.000 claims description 5
- 238000009795 derivation Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000004913 activation Effects 0.000 claims description 2
- 239000000446 fuel Substances 0.000 description 9
- 230000008569 process Effects 0.000 description 9
- 230000000052 comparative effect Effects 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000000354 decomposition reaction Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000015556 catabolic process Effects 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
- G06F2218/06—Denoising by applying a scale-space analysis, e.g. using wavelet analysis
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biophysics (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Geometry (AREA)
- Medical Informatics (AREA)
- Computer Hardware Design (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于数据驱动的航空发动机状态变量模型建立方法,包括:根据所要建立的航空发动机状态变量模型的结构,采集所述航空发动机状态变量模型的输入、输出、状态变量数据;对所采集的数据进行去噪后,以其作为训练数据对神经网络模型进行训练,得到航空发动机智能网络模型;基于所建立的航空发动机智能网络模型,根据所要建立的航空发动机状态变量模型参数的偏导数表达,以链式求导法则进行输出量到输入量的偏导数计算,从而建立基于数据驱动的航空发动机状态变量模型。本发明可以直接根据发动机输入和输出数据,实时获取发动机个体的高精度状态变量模型,避免了对发动机部件级数学模型的依赖,实现基于发动机个体数据的准确建模。
Description
技术领域
本发明属于航空宇航推进理论与工程中的系统仿真与控制领域,具体涉及一种航空发动机状态变量模型建立方法。
背景技术
状态变量模型作为线性化数学模型的一种,是采用现代控制理论的知识对航空发动机进行多变量性能分析和设计的基础。同时,在发动机模型预测控制、性能寻优控制、健康评估中,采用状态变量模型进行发动机性能在线预测,能够提高控制系统的实时性。鉴于发动机的非线性,能实时准确反映发动机性能的状态变量模型,则可以提高控制和评估的精度。
航空发动机状态变量模型主要有两种建模方法,即拟合法和偏导数法。两种方法本质上都属于小扰动法,通过在部件级模型的稳态工作点上对输入量或者控制量分别进行小扰动来获得模型参数。拟合法可以采用群智能算法、最小二乘算法等,通过最小化模型输出和仿真数据之间的偏差来确定模型参数。偏导数法则是在部件级模型上,对系统控制变量、状态变量分别做小扰动,通过差分近似偏导数,进而得出系统的状态变量模型。这两种方法都基于发动机的部件级数学模型展开,以便于在同一个工作点上,分别进行各个变量的扰动。
因此,上述两种方法的实施的前提是具有高精度的部件级模型,而建立高精度的部件级模型一直是航空发动机控制领域的难点问题。航空发动机部件级数学模型是由多个部件模型、共同工作方程以及方程求解过程组成的非常复杂的系统。由于建模中的假设,使得部件级数学模型的精度难以保证,且这种模型多基于额定状态的部件特性建立,不能反映发动机个体差异和性能退化,基于这样的部件模型建立的状态变量模型也不能反映发动机个体特性。同时这两种状态变量建模方法在参数拟合和偏导数计算过程中耗时较多,实时性较差,只能在有限的稳态工作点上进行建模。最新研究的基于精确偏导数的状态变量模型建模方法,本质上同样是一种小扰动法,同样基于发动机部件级模型展开,扰动控制量和状态量,这就增加了状态变量模型建模的前期工作,且状态变量建模过程需要部件级模型参与计算,计算量大,模型的精度也会受到部件级模型精度的影响。
由于发动机个体差异带来的控制性能差异,也要对发动机个体进行建模,分析控制器参数调整方法。通过状态变量模型能够方便计算发动机极点位置和静态增益,对于控制器参数调整具有非常重要的作用。如果在对发动机个体建立部件级模型之后再进行状态变量模型的建立,则要耗费大量的人力和时间。同时航空发动机模型预测控制、性能寻优控制、健康评估等技术也需要对发动机个体进行状态变量模型在线建立。而对发动机个体进行建模的最方便的方法就是基于数据驱动的方法,基于发动机输入和输出数据建立智能网络模型,但是一直以来,智能网络模型都用来建立非线性模型,实现某一特定用途,比如传感器故障诊断、推力估计等,不能满足线性化建模的需求。
发明内容
本发明所要解决的技术问题在于克服现有技术不足,提供一种基于数据驱动的航空发动机状态变量模型建立方法,可以直接根据发动机输入和输出数据,实时获取发动机个体的高精度状态变量模型,避免了现有状态变量建模方法对发动机部件级数学模型的依赖,实现基于发动机个体数据的准确建模。
本发明具体采用以下技术方案解决上述技术问题:
基于数据驱动的航空发动机状态变量模型建立方法,包括:
步骤A、根据所要建立的航空发动机状态变量模型的结构,采集所述航空发动机状态变量模型的输入、输出、状态变量数据;
步骤B、对所采集的数据进行去噪后,以其作为训练数据对神经网络模型进行训练,得到航空发动机智能网络模型;所述航空发动机智能网络模型的输入为当前时刻的航空发动机状态变量模型的输入和状态变量,所述航空发动机智能网络模型的输出为下一时刻的航空发动机状态变量模型的状态变量以及当前时刻的航空发动机状态变量模型的输出;
步骤C、基于所建立的航空发动机智能网络模型,根据所要建立的航空发动机状态变量模型参数的偏导数表达,以链式求导法则进行输出量到输入量的偏导数计算,得到航空发动机状态变量模型的参数,从而建立基于数据驱动的航空发动机状态变量模型。
优选地,所述神经网络模型为在线训练极限学习机模型OS-ELM。
更进一步优选地,所述OS-ELM为三层前向神经网络结构,隐含层激活函数为sigmoid函数。
优选地,使用以下方法对所采集的数据进行去噪:先使用小波变换的残差直方图剔除数据中的野点,然后对剔除野点后的数据进行小波滤波。
更进一步优选地,所述小波变换和小波滤波均使用尺度为7的‘db8’小波。
优选地,所要建立的航空发动机状态变量模型的结构按照以下方法确定:
先将发动机离散形式的非线性工作特性表示为:
式中,x(k)、u(k)、y(k)分别为航空发动机状态变量模型的状态变量、输入量、输出量;
然后将发动机在k时刻工作点(xk,0,uk,0,yk,0)处增量形式的线性化数学模型表示为:
其中,Δ为增量符号,Δx(k)=x(k)-xk,0,Δu(k)=u(k)-uk,0,Δy(k)=y(k)-yk,0。
相比现有技术,本发明技术方案具有以下有益效果:
(1)具有准确描述发动机非线性变参数特征的能力:传统的拟合法基于一组动态响应数据,只能建立一个定常的状态变量模型,忽略了不同转速下发动机的非线性变参数特性。本发明采用基于数据驱动的方法,不同转速下的智能网络模型参数不同,依据智能网络模型参数建立的状态变量模型不同,充分体现发动机的非线性变参数特性。
(2)具有不依赖于部件级模型的特点:本发明实现了脱离部件级数学模型的状态变量建模,利用智能网络模型进行偏导数计算,避免了复杂的部件级模型建立过程,避免了部件级模型精度对状态变量模型的影响;
(3)具有建模效率高的特点:本发明基于数据驱动方法进行建模,无需进行任何扰动,在每一个仿真步长均能实时获得当前时刻的状态变量模型,
(4)具有精度高可移植性强的特点:本发明基于在线训练的方法建立发动机智能模型,模型完全依赖于对象的输入输出数据,不受发动机制造安装公差和性能退化的影响,模型精度高,状态变量模型参数计算过程不受工作机理约束,适用于各类复杂的非线性系统。
附图说明
图1为原始燃油信号小波变换残差直方图;
图2为剔除野点后燃油信号小波变换残差直方图;
图3(a)为小波滤波前的部分燃油数据;
图3(b)为小波滤波后的部分燃油数据;
图4为基于数据驱动的航空发动机数学模型结构图;
图5(a)为ELM网络n1输出及其误差的对比图;
图5(b)为ELM网络T6输出及其误差的对比图;
图5(c)为ELM网络F输出及其误差的对比图;
图6(a)为2步预测时状态变量模型n1输出及其误差的对比图;
图6(b)为2步预测时状态变量模型n2输出及其误差的对比图;
图6(c)为2步预测时状态变量模型P3输出及其误差的对比图;
图6(d)为2步预测时状态变量模型P6输出及其误差的对比图;
图6(e)为2步预测时状态变量模型T6输出及其误差的对比图;
图6(f)为2步预测时状态变量模型F输出及其误差的对比图;
图7(a)为5步预测时状态变量模型n1输出及其误差的对比图;
图7(b)为5步预测时状态变量模型n2输出及其误差的对比图;
图7(c)为5步预测时状态变量模型P3输出及其误差的对比图;
图7(d)为5步预测时状态变量模型P6输出及其误差的对比图;
图7(e)为5步预测时状态变量模型T6输出及其误差的对比图;
图7(f)为5步预测时状态变量模型F输出及其误差的对比图。
具体实施方式
针对传统航空发动机状态变量建模依赖于部件级模型及难以在线实时建模的问题,本发明根据控制系统需求,基于人工神经网络建立航空发动机的智能模型,基于建立的智能模型采取偏导数法进行状态变量模型参数计算,实现基于数据的发动机状态变量参数辨识。
具体地,本发明基于数据驱动的航空发动机状态变量模型建立方法,包括:
步骤A、根据所要建立的航空发动机状态变量模型的结构,采集所述航空发动机状态变量模型的输入、输出、状态变量数据;
步骤B、对所采集的数据进行去噪后,以其作为训练数据对神经网络模型进行训练,得到航空发动机智能网络模型;所述航空发动机智能网络模型的输入为当前时刻的航空发动机状态变量模型的输入和状态变量,所述航空发动机智能网络模型的输出为下一时刻的航空发动机状态变量模型的状态变量以及当前时刻的航空发动机状态变量模型的输出;
步骤C、基于所建立的航空发动机智能网络模型,根据所要建立的航空发动机状态变量模型参数的偏导数表达,以链式求导法则进行输出量到输入量的偏导数计算,得到航空发动机状态变量模型的参数,从而建立基于数据驱动的航空发动机状态变量模型。
为了便于公众理解,下面以某型双转子混合排气加力式涡扇发动机状态变量建模为例,并结合附图来对本发明的技术方案进行详细说明:
首先将发动机的非线性模型表示成离散状态空间形式:
式中,输入量u(k)=[Wf(k),A8(k)]T为主燃油流量和尾喷口喉道面积,输出量y(k)=[P3(k),T6(k),P6(k),F(k)]T,依次为压气机出口总压、低压涡轮出口总温、低压涡轮出口总压和推力,状态量x(k)=[n1(k),n2(k)]T为低压转速、高压转速,上标T为转置运算符号。
将发动机在k时刻工作点(xk,0,uk,0,yk,0)处增量形式的线性化数学模型表示为:
其中,Δ为增量符号,Δx(k)=x(k)-xk,0,Δu(k)=u(k)-uk,0,Δy(k)=y(k)-yk,0。
对采集到的状态变量模型输入、输出、状态变量数据进行检查,处理无效数据:
本次建模所用数据来源于某型涡扇发动机的试车数据,为避免数据中野点对建模产生影响,利用小波的局部奇异性估算能力,采用小波变换的残差直方图来确定野点判别的残差阈值。
记原始信号S经过L层小波分解后的残差为
Re=S-A (3)
其中,A为经过L层分解后的低频重构信号。
当残差较小时,低频重构信号趋近于原始信号,当残差较大时,则两者之间有高频噪声信号,当残差过大时,则可能是野点。以某次试车的燃油数据为例,采用尺度为7的‘db8’小波对信号进行变换,得到分解后的残差直方图如图1所示,其横坐标为残差,纵坐标为对应区间数据的个数。本次测试共采集数据111406组。如图1可见,绝大部分残差都处于直方图区域中间位置,具体残差分布见表1。如表1可见,处于残差[-0.721%,1.39%]区间的数据占总数据的98.54%,处于残差区间[-2.83%,3.49%]区间的数据占总数据的99.94%,在残差大于3.49%和小于-2.83%区间数据分布不具有随残差增大而减小的趋势,为此将野点判别残差阈值设定为-2.83%与3.49%,超出此范围的则认为是野点,将其用左右两侧数据的平均值来取代,对于连续的野点,则通过前后合理范围数据插值替代。可以得到剔除野点后数据小波分解后的残差直方图如图2所示,超过残差阈值范围的数据点仅剩下两个,表明有效剔除了大部分的野点。对应剩余野点,可再次进行剔除。
表1燃油信号小波变换残差分布统计
对剔除无效数据后的数据后进行滤波,消除噪声。同样采用尺度为7的‘db8’小波对信号滤波,同样以燃油为例,滤波前后部分数据对比结果如图3(a)和图3(b)。由图3可见,经过滤波后,数据中的噪声明显降低,有助于采用基于数据驱动的方法建立模型。
本实施例将滤波后的数据归一化到0.5到1.5之间。
数据归一化方法如下:
其中,xmax为变量x的最大可能测量值、xmin为变量x的最小可能测量值,为归一化后的数据。
本实施例中的智能网络模型采用在线训练极限学习机(OS-ELM),智能网络模型的输入为状态变量模型k时刻的输入量和状态量,智能网络模型的输出为状态变量模型k+1时刻的状态量和k时刻的输出量。设计基于OS-ELM的发动机数学模型结构如图4所示,采用三层前向神经网络结构,模型输入记为XELM(k)=[Wf(k) A8(k) n1(k) n2(k)]T,为k时刻发动机的主燃油流量Wf(k)、尾喷口喉道面积A8(k)、低压转速n1(k)、高压转速n2(k),输出记为YELM(k)=[n1(k+1),n2(k+2),P3(k),T6(k),P6(k),F(k)]T,为k+1时刻的低压转速n1(k+1)、高压转速n2(k+1)和k时刻的压气机出口总压P3(k)、涡轮出口总温T6(k)、涡轮出口总压P6(k)和推力F(k),隐含层节点数取为85个,隐含层激活函数取为sigmoid函数。
上述智能网络模型的具体训练过程如下:
随机生成OS-ELM输入层到隐含层的连接权值,基于第1组输入计算隐含层的输出O(1),则有
其中,O(1)∈R85×1,WELM∈R85×4为输入层到隐含层的连接权值,bELM∈R85×1为隐含层偏置,f()为隐含层激励函数。
记OS-ELM网络隐含层到输出层的连接权值为β,则网络输出
其中,β(1)∈R85×6,β(1)=[β1(1) β2(1) β3(1) β4(1) β5(1) β6(1)]T
记T为目标输出,则TT(1)=[n1(2),n2(2),P3(1),T6(1),P6(1),F(1)],则隐含层到输出层的连接权值
β(1)=(O(1)O(1)T)-1O(1)TT(1) (7)
基于递推最小二乘算法更新OS-ELM网络隐含层权值。
若OS-ELM网络的第k+1组输出数据某一变量与发动机输出数据相比误差处于0.15%<|ei|<4%,i=1,2,…,6,则网络权值需要调整。0.15%到4%为预设的权值更新阈值,小于0.15%则认为输出精度满足要求,大于4%则认为数据处于非正常误差范围,为避免错误数据对网络造成影响,不对权值进行调整。在权值更新范围内,依据递推最小二乘算法,有
其中β(k)=[β1(k) β2(k) β3(k) β4(k) β5(k) β6(k)],βi(k)∈R85×1,i=1,2…,6。
以低压转速、涡轮出口温度和推力为例,在图5(a)~图5(c)中给出了OS-ELM输出与发动机输出之间的对比,图中下标“ELM”代表OS-ELM神经网络输出,下标“engine”代表滤波后的发动机输出,下标“error”代表相对误差。图5(a)为n1转速输出及其误差的对比图,图5(b)为涡轮出口总温T6输出及其误差的对比图,图5(c)为推力F输出及其误差的对比图,图中n1、T6和F均为相对于设计点的百分比数据,误差为相对发动机输出的百分比误差。可以看出,网络对低压转速的建模效果最好,最大动态误差小于0.12%,对涡轮出口总温和推力的建模效果相近,最大误差均在3%左右,所有数据的相对误差大部分时刻可以控制在0.5%以内,相比于部件级模型,取得了较高的建模精度。
确定状态变量模型参数的偏导数表达方法。
对于所研究的涡扇发动机,其状态变量模型中的矩阵可记为
其中
基于所建立的OS-ELM网络,以链式求导法则进行偏导数计算。
OS-ELM网络采用线性输出节点,网络输出为:
n1(k+1)=OT(k)β1(k),n2(k+1)=OT(k)β2(k),P3(k)=OT(k)β3(k),
T6(k)=OT(k)β4(k),P6(k)=OT(k)β5(k),F(k)=OT(k)β6(k) (10)
则有
其中,WELM1j为输入到隐含层连接权值WELM中第1列的第j个元素,WELM第1列对应第一个输入的连接权。
按照a11的计算过程,完成k时刻所有状态变量矩阵参数的计算,构建出k时刻的状态变量模型。
为了验证本发明的有效性,选择不同时刻建立的状态变量模型进行仿真,对比状态变量模型在未来时刻输出与神经网络模型输出之间的误差。
对k时刻建立的发动机状态变量模型,在未来m时刻的输出
其中,um为从k时刻起,未来m个时刻的发动机输入,xm为从k时刻起,未来m个时刻的状态变量,ym为从k时刻起,未来m个时刻的状态变量模型输出。
以m=2和m=5为例,对比状态变量模型输出和滤波后的发动机输出之间的误差,m=2时结果如图6(a)~图6(f)所示,m=5时结果如图7(a)~图7(f)所示,图中下标“SVM”代表状态变量模型输出,下标“ELM”,代表OS-ELM网络输出。图6(a)为2步预测时状态变量模型n1输出及其误差的对比图,图6(b)为2步预测时状态变量模型n2输出及其误差的对比图,图6(c)为2步预测时状态变量模型P3输出及其误差的对比图,图6(d)为2步预测时状态变量模型P6输出及其误差的对比图,图6(e)为2步预测时状态变量模型T6输出及其误差的对比图,图6(f)为2步预测时状态变量模型F输出及其误差的对比图。图7(a)为5步预测时状态变量模型n1输出及其误差的对比图,图7(b)为5步预测时状态变量模型n2输出及其误差的对比图,图7(c)为5步预测时状态变量模型P3输出及其误差的对比图,图7(d)为5步预测时状态变量模型P6输出及其误差的对比图,图7(e)为5步预测时状态变量模型T6输出及其误差的对比图,图7(f)为5步预测时状态变量模型F输出及其误差的对比图。
由图6(a)~图6(f)可以看出,在状态变量模型两步预测中,转速的预测误差最小,低压转速最大预测误差小于0.04%,高压转速最大预测误差小于0.01%,压气机出口总压、涡轮出口总温和总压的最大误差也都小于1%,推力的预测误差偏低,但最大误差也可以控制在2%以内,稳态预测误差均趋近于零。由图7(a)~图7(f)可以看出,预测五步时,转速的最大预测误差仍然可以控制在0.1%以内,压气机出口总压、涡轮出口总温和总压的最大误差也都小于2.5%,推力的预测误差偏低,最大动态预测误差在5%以内,且稳态时误差趋近于0。可见基于数据驱动的状态变量模型取得了较好的建模效果。
为了验证本发明建立状态变量模型过程的实时性,对OS-ELM建模过程和偏导数建模过程进行统一计时,仿真运行在CPU为Intel(R)Xeon(R)E5-1620 v33.5GHz,内存16G的平台。经验证,基于111406组数据建立的111405组状态变量模型共耗时300.7935s,每个状态变量模型平均耗时2.7ms,低于20ms的采样步长,满足实时性的需求。
Claims (5)
1.基于数据驱动的航空发动机状态变量模型建立方法,其特征在于,包括:
步骤A、根据所要建立的航空发动机状态变量模型的结构,采集所述航空发动机状态变量模型的输入、输出、状态变量数据;所要建立的航空发动机状态变量模型的结构按照以下方法确定:
先将发动机离散形式的非线性工作特性表示为:
式中,x(k)、u(k)、y(k)分别为航空发动机状态变量模型的状态变量、输入量、输出量;
然后将发动机在k时刻工作点(xk,0,uk,0,yk,0)处增量形式的线性化数学模型表示为:
其中,Δ为增量符号,Δx(k)=x(k)-xk,0,Δu(k)=u(k)-uk,0,Δy(k)=y(k)-yk,0,A、B、C、D为航空发动机状态变量模型中的矩阵;
步骤B、对所采集的数据进行去噪后,以其作为训练数据对神经网络模型进行训练,得到航空发动机智能网络模型;所述航空发动机智能网络模型的输入为当前时刻的航空发动机状态变量模型的输入和状态变量,所述航空发动机智能网络模型的输出为下一时刻的航空发动机状态变量模型的状态变量以及当前时刻的航空发动机状态变量模型的输出;
步骤C、基于所建立的航空发动机智能网络模型,根据所要建立的航空发动机状态变量模型参数的偏导数表达,以链式求导法则进行输出量到输入量的偏导数计算,得到航空发动机状态变量模型的参数,从而建立基于数据驱动的航空发动机状态变量模型。
2.如权利要求1所述基于数据驱动的航空发动机状态变量模型建立方法,其特征在于,所述神经网络模型为在线训练极限学习机模型OS-ELM。
3.如权利要求2所述基于数据驱动的航空发动机状态变量模型建立方法,其特征在于,所述OS-ELM为三层前向神经网络结构,隐含层激活函数为sigmoid函数。
4.如权利要求1所述基于数据驱动的航空发动机状态变量模型建立方法,其特征在于,使用以下方法对所采集的数据进行去噪:先使用小波变换的残差直方图剔除数据中的野点,然后对剔除野点后的数据进行小波滤波。
5.如权利要求4所述基于数据驱动的航空发动机状态变量模型建立方法,其特征在于,所述小波变换和小波滤波均使用尺度为7的‘db8’小波。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110419235.4A CN113158564B (zh) | 2021-04-19 | 2021-04-19 | 基于数据驱动的航空发动机状态变量模型建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110419235.4A CN113158564B (zh) | 2021-04-19 | 2021-04-19 | 基于数据驱动的航空发动机状态变量模型建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113158564A CN113158564A (zh) | 2021-07-23 |
CN113158564B true CN113158564B (zh) | 2024-03-29 |
Family
ID=76868733
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110419235.4A Active CN113158564B (zh) | 2021-04-19 | 2021-04-19 | 基于数据驱动的航空发动机状态变量模型建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113158564B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114723184B (zh) * | 2022-06-08 | 2022-10-11 | 广东数字生态科技有限责任公司 | 一种基于视觉感知的风力发电机测量方法、装置以及设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2336918A1 (en) * | 2009-12-21 | 2011-06-22 | United Technologies Corporation | Method and system for modeling the performance of a gas turbine engine |
CN107357176A (zh) * | 2017-08-05 | 2017-11-17 | 中国科学院工程热物理研究所 | 一种航空发动机试车数据建模方法 |
CN109031951A (zh) * | 2018-07-16 | 2018-12-18 | 南京航空航天大学 | 基于精确偏导数的航空发动机状态变量模型在线建立方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7216071B2 (en) * | 2002-04-23 | 2007-05-08 | United Technologies Corporation | Hybrid gas turbine engine state variable model |
-
2021
- 2021-04-19 CN CN202110419235.4A patent/CN113158564B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2336918A1 (en) * | 2009-12-21 | 2011-06-22 | United Technologies Corporation | Method and system for modeling the performance of a gas turbine engine |
CN107357176A (zh) * | 2017-08-05 | 2017-11-17 | 中国科学院工程热物理研究所 | 一种航空发动机试车数据建模方法 |
CN109031951A (zh) * | 2018-07-16 | 2018-12-18 | 南京航空航天大学 | 基于精确偏导数的航空发动机状态变量模型在线建立方法 |
Non-Patent Citations (4)
Title |
---|
Data-driven neural network methodology to remaining life predictions for aircraft actuator components;C.S. Byington 等;2004 IEEE Aerospace Conference Proceedings (IEEE Cat. No.04TH8720);20040313;第3581-3589页 * |
From Model, Signal to Knowledge: A Data-Driven Perspective of Fault Detection and Diagnosis;Xuewu Dai 等;IEEE Transactions on Industrial Informatics;20130130;第9卷(第4期);第2226-2238页 * |
基于改进拟合法的涡扇发动机状态变量模型建立方法;胡宇 等;推进技术;20130313;第34卷(第3期);第405-410页 * |
航空发动机智能鲁棒控制研究;李秋红;中国博士学位论文全文数据库 工程科技Ⅱ辑;20120715;第C031-21页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113158564A (zh) | 2021-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108647428B (zh) | 一种涡扇发动机自适应部件级仿真模型构建方法 | |
CN108256173B (zh) | 一种航空发动机动态过程的气路故障诊断方法及系统 | |
CN110502840B (zh) | 航空发动机气路参数在线预测方法 | |
CN108573116B (zh) | 一种基于长短时记忆网络的航空发动机过渡态推力估计方法 | |
CN109635318B (zh) | 一种航空发动机传感器智能解析余度设计方法 | |
CN111042928B (zh) | 一种基于动态神经网络的变循环发动机智能控制方法 | |
DE102008055593A1 (de) | Verfahren und Systeme zur neuronalen Netzwerkmodellierung von Turbinenkomponenten | |
CN108062428B (zh) | 一种涡扇发动机在线部件故障诊断方法及系统 | |
CN111860791A (zh) | 基于相似变换的航空发动机推力估计方法及装置 | |
CN112729857B (zh) | 航空发动机健康参数估计方法及航空发动机自适应模型 | |
CN109800449B (zh) | 一种基于神经网络的航空发动机压缩部件特性修正方法 | |
CN107357176B (zh) | 一种航空发动机试车数据建模方法 | |
CN111324110A (zh) | 一种基于多块收缩自动编码器的发酵过程故障监测方法 | |
CN113158564B (zh) | 基于数据驱动的航空发动机状态变量模型建立方法 | |
CN111651940A (zh) | 基于pso-narx的涡轴发动机稳态模型辨识方法 | |
CN112149233A (zh) | 基于回声状态网络的航空发动机动态推力估计方法 | |
CN109976311A (zh) | 一阶定值控制系统传感器故障诊断方法及诊断系统 | |
CN112378619A (zh) | 具有ReMD-OSELM的FER-FSE在风洞试验冲压阶段总压实时建模中的应用 | |
CN115688609A (zh) | 一种用于航空发动机的智能推力预测及实时预警方法 | |
CN112364446B (zh) | 基于ec-rbelm算法的发动机整机性能衰减预测方法 | |
CN116085245A (zh) | 一种基于os-elm的压缩机性能在线预测方法及系统 | |
Vladov et al. | Neural Network Modeling of Helicopters Turboshaft Engines at Flight Modes Using an Approach Based on" Black Box" Models. | |
CN109635358A (zh) | 一种基于滑动窗口多尺度主元分析的机组故障检测方法 | |
CN113048012B (zh) | 基于混合高斯模型的风电机组偏航角度识别方法和装置 | |
Liu et al. | The research of aero-engine digital twin model based on model-driven and data-driven fusion method |
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 |