CN110765659A - 一种针对航空发动机故障诊断的支持三阶张量机建模方法 - Google Patents

一种针对航空发动机故障诊断的支持三阶张量机建模方法 Download PDF

Info

Publication number
CN110765659A
CN110765659A CN201911138245.XA CN201911138245A CN110765659A CN 110765659 A CN110765659 A CN 110765659A CN 201911138245 A CN201911138245 A CN 201911138245A CN 110765659 A CN110765659 A CN 110765659A
Authority
CN
China
Prior art keywords
tensor
data
aircraft engine
order tensor
matrix
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.)
Pending
Application number
CN201911138245.XA
Other languages
English (en)
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.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
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 Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN201911138245.XA priority Critical patent/CN110765659A/zh
Publication of CN110765659A publication Critical patent/CN110765659A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Testing Of Engines (AREA)

Abstract

本发明属于航空发动机的故障诊断技术领域,提供一种针对航空发动机故障诊断的支持三阶张量机建模方法。将航空发动机的时间序列数据构建为一种具有三阶张量结构的数据形式;在构建的三阶张量数据下,利用拉格朗日乘子方法获得具有张量內积运算的对偶支持三阶张量机模型;通过张量的Tucker分解方法来近似表示该张量內积运算,并获得具有张量Tucker分解形式下的对偶支持三阶张量机。本发明的优点:一方面它避免了数据本身在训练的过程中可能遇到的维度灾难问题,另一方面它能够保留原始航空发动机时间序列数据的大部分耦合信息,并能精准预测航空发动机在未来飞行中是否发生喘振故障。

Description

一种针对航空发动机故障诊断的支持三阶张量机建模方法
技术领域
本发明属于航空发动机的故障诊断技术领域,涉及一种针对航空发动机故障诊断的支持三阶张量机建模方法,具体是针对航空发动机的喘振故障来设计一种支持三阶张量机的诊断预测模型。
背景技术
随着航空事业的不断发展,飞机的安全性和可靠性越来越引起了人们的重视。如果飞机一旦发生故障问题,轻则影响飞机的性能,重则机毁人亡。另外,值得注意的是航空发动机既是飞机的心脏,也是一个最容易发生故障的核心部件。在航空发动机的众多故障问题中,喘振故障是一种最为常见的危险因素,它不仅直接影响着航空发动机的性能,而且对飞机的安全性以及整个航空事业的发展造成了巨大的威胁。解决这个问题的有效方法主要是针对航空发动机以往的喘振故障数据提出一种数据驱动模型,从而预测航空发动机的未来是否发生喘振故障并采取相应的措施来保证飞机的安全性。
目前,预测航空发动机喘振故障问题的方法有以下几种:
1)基于故障树的分析方法。
该方法主要是从喘振故障的监测出发来模拟飞机飞行中的实时数据,并根据监控的状况采用故障树方法来分析故障的详细原因。然而,这种分析方法要求分析故障的人员必须熟悉所分析的对象系统以及熟悉地应用该分析方法,这样就造成了不同的分析人员会给出不同的故障树结果。另外对于故障树的计算也是非常的复杂,很难精确的计算发生故障的具体原因。
2)基于最小二乘支持向量机的分析方法。
该方法是利用航空发动机的气路参数,建立最小二乘支持向量机模型来对航空发动机进行状态监控。即,根据建立的模型来监控航空发动机的低压转子转速(n1),高压转子转速(n2)和尾喷管出口温度(T6),并通过预测值与真实值的相对误差率来分析喘振故障。但是,最小二乘支持向量机是把喘振的数据集作为一种向量模式来处理,这样就没有充分考虑到数据之间相互耦合,相互影响的自然关系。另外,如果数据集被强行的采用向量表示,也会导致原始数据的时序关联性遭到破坏,产生的数值误差也是难免的。
综合以上论述,本发明设计的支持三阶张量机是一种能够精确预测具有耦合时间序列数据的航空发动机喘振故障。
发明内容
本发明针对最小二乘支持向量机模型在航空发动机故障诊断中所带来的局限性问题,提供了一种支持三阶张量机模型,并获得了更好的预测精确度。由于航空发动机是一种高度复杂的气动-热力-机械系统,它所生成的时间序列数据具有很强的时序关联性,耦合性与多模态特征,因此,如何在多变的全包线环境下来预测航空发动机的喘振故障一直是一个挑战性的难题。
本发明的技术方案:
一种针对航空发动机故障诊断的支持三阶张量机建模方法,用于预测航空发动机在未来飞行时间中是否会发生喘振故障问题;首先,将航空发动机的时间序列数据构建为一种具有三阶张量结构的数据形式;然后,在构建的三阶张量数据下,利用拉格朗日乘子方法获得具有张量內积运算的对偶支持三阶张量机模型;最后,通过张量的Tucker分解方法来近似表示该张量內积运算,并获得具有张量Tucker分解形式下的对偶支持三阶张量机;
具体步骤如下:
步骤1:通过奇异谱分析的嵌入思想,把具有一维向量模式的航空发动机数据x=[x1,x2,…,xn]T构建为一种长宽高分别为I1,I2,I3的三阶张量数据
Figure BDA00022801437000000212
(1.1)首先,将航空发动机的一维向量数据x=[x1,x2,…,xn]T分割为一个如下的矩阵形式:
Figure BDA0002280143700000021
其中,x1,x2,…,xn是航空发动机在飞行过程中离散采集的n个时间序列数据,是将上述离散时间序列数据通过分割移动构成的I1×l的矩阵,n=I1l是数据的总数,l=I2+I3-1是矩阵
Figure BDA0002280143700000022
的列数,I1,I2与I3分别为张量
Figure BDA0002280143700000023
在长、宽与高等方向的总层数;
(1.2)其次,通过奇异谱分析的嵌入方法将矩阵
Figure BDA0002280143700000024
的元素嵌入到三阶张量
Figure BDA0002280143700000025
在长方向上的每一层切片中,即第i层切片是如下的矩阵形式
Figure BDA0002280143700000026
其中,
Figure BDA00022801437000000215
是矩阵
Figure BDA0002280143700000027
中的第i行的元素,Xi∷是第i层长方向上的切片矩阵;
(1.3)最后,通过定义三阶张量数据
Figure BDA0002280143700000028
中的每一个元素来获得构建的张量
Figure BDA0002280143700000029
即张量
Figure BDA00022801437000000210
的第i1i2i3元素是如下形式:
Figure BDA00022801437000000211
其中i1=1,2,…,I1;i2=1,2,…,I2;i3=1,2,…,I3
Figure BDA0002280143700000031
是矩阵
Figure BDA0002280143700000032
中的第i1行与第i2+i3-1列处的元素;
步骤2:在构建的三阶张量数据下,利用拉格朗日乘子方法获得具有张量內积运算的对偶支持三阶张量机模型,即內积形式的对偶支持三阶张量机模型如下:
Figure BDA0002280143700000033
其中,Xp与Xq分别是构造的第p个与第q个三阶张量数据,yp与yq是航空发动机数据定义的正常标签或喘振标签,αpq是第p个与第q个拉格朗日乘子,d是张量数据的总个数,c是正则化参数;
步骤3:针对具有內积形式的对偶支持三阶张量机模型,通过张量的Tucker分解方法来近似获得具有Tucker分解形式的对偶支持三阶张量机,即Tucker分解形式的对偶支持三阶张量机模型如下:
Figure BDA0002280143700000034
其中,
Figure BDA0002280143700000035
是对张量Xp(Xq)进行Tucker分解后获得的核张量Gp(Gq)中的第k1k2k3(r1r2r3)元素,符号<,>是內积运算,∑是求和运算,∏是求积运算,
Figure BDA0002280143700000036
是对张量Xp(Xq)进行Tucker分解后获得的因子矩阵
Figure BDA0002280143700000037
中的第kt(rt)列向量,n1,n2,n3(s1,s2,s3)是核张量Gp(Gq)在长,宽与高等方向的总层数,t是因子矩阵的下角标,并说明第t个因子矩阵是
本发明的有益效果:Tucker分解形式的对偶支持三阶张量机模型充分考虑了航空发动机数据之间相互耦合,相互影响的自然关系。一方面它避免了数据本身在训练的过程中可能遇到的维度灾难问题,另一方面它能够保留原始航空发动机时间序列数据的大部分耦合信息,并能精准预测航空发动机在未来飞行中是否发生喘振故障。
附图说明
图1是航空发动机在发生喘振前后,通过传感器测量获得的低压转子转速(n1)的部分数据趋势图。注:趋势图中的实心点表示航空发动机的正常数据,空心点表示航空发动机的喘振数据。
图2是航空发动机在发生喘振前后,通过传感器测量获得的高压转子转速(n2)的部分数据趋势图。注:趋势图中的实心点表示航空发动机的正常数据,空心点表示航空发动机的喘振数据。
图3是航空发动机在发生喘振前后,通过传感器测量获得的尾喷管出口温度(T6)的部分数据趋势图。注:趋势图中的实心点表示航空发动机的正常数据,空心点表示航空发动机的喘振数据。
具体实施方式
以下结合附图和技术方案,进一步说明本发明的具体实施方式。
实施例
步骤1:通过奇异谱分析的嵌入方法,把具有一维时间序列的低压转子转速(图1),高压转子转速(图2)和尾喷管出口温度(图3)等数据集构建为一种具有元素(1)形式的三阶张量数据集。
步骤2:在低压转子转速(图1),高压转子转速(图2)和尾喷管出口温度(图3)构建的三阶张量数据集下,通过拉格朗日乘子方法获得具有內积形式(2)描述的对偶支持三阶张量机模型。
步骤3:通过张量的Tucker分解方法来近似內积形式的对偶支持三阶张量机模型(2),并获得了具有Tucker分解形式(3)描述的对偶支持三阶张量机模型,从而来进一步的预测航空发动机是否发生喘振故障。
表1:航空发动机基于支持三阶张量机模型的故障诊断结果
Figure BDA0002280143700000041
实施结果
1)从表1中可以看出,在相同的测试数据集下,随着训练样本数据集的增加,具有Tucker分解形式的对偶支持三阶张量机模型对该数据集的分类精确度越来越高。
2)从表1中可以看出,在相同的训练数据集下,随着测试样本数据集的减少,具有Tucker分解形式的对偶支持三阶张量机模型对该数据集的分类精确度也越来越高。
因此,这样的结果符合支持三阶张量机模型的本质特性。同时也证明了支持三阶张量机模型对航空发动机的喘振故障问题具有精准的预测诊断能力。

Claims (1)

1.一种针对航空发动机故障诊断的支持三阶张量机建模方法,用于预测航空发动机在未来飞行时间中是否会发生喘振故障问题;首先,将航空发动机的时间序列数据构建为一种具有三阶张量结构的数据形式;然后,在构建的三阶张量数据下,利用拉格朗日乘子方法获得具有张量内积运算的对偶支持三阶张量机模型;最后,通过张量的Tucker分解方法来近似表示该张量內积运算,并获得具有张量Tucker分解形式下的对偶支持三阶张量机;
其特征在于,具体步骤如下:
步骤1:通过奇异谱分析的嵌入思想,把具有一维向量模式的航空发动机数据x=[x1,x2,…,xn]T构建为一种长宽高分别为I1,I2,I3的三阶张量数据
Figure FDA0002280143690000011
(1.1)首先,将航空发动机的一维向量数据x=[x1,x2,…,xn]T分割为一个如下的矩阵形式:
Figure FDA0002280143690000012
其中,x1,x2,…,xn是航空发动机在飞行过程中离散采集的n个时间序列数据,
Figure FDA0002280143690000013
是将上述离散时间序列数据通过分割移动构成的I1×l的矩阵,n=I1l是数据的总数,l=I2+I3-1是矩阵
Figure FDA0002280143690000014
的列数,I1,I2与I3分别为张量在长、宽与高等方向的总层数;
(1.2)其次,通过奇异谱分析的嵌入方法将矩阵
Figure FDA0002280143690000016
的元素嵌入到三阶张量
Figure FDA0002280143690000017
在长方向上的每一层切片中,即第i层切片是如下的矩阵形式
Figure FDA0002280143690000018
其中,
Figure FDA0002280143690000019
是矩阵
Figure FDA00022801436900000110
中的第i行的元素,Xi::是第i层长方向上的切片矩阵;
(1.3)最后,通过定义三阶张量数据
Figure FDA00022801436900000111
中的每一个元素来获得构建的张量
Figure FDA00022801436900000112
即张量
Figure FDA00022801436900000113
的第i1i2i3元素是如下形式:
Figure FDA00022801436900000114
其中i1=1,2,…,I1;i2=1,2,…,I2;i3=1,2,…,I3是矩阵中的第i1行与第i2+i3-1列处的元素;
步骤2:在构建的三阶张量数据下,利用拉格朗日乘子方法获得具有张量內积运算的对偶支持三阶张量机模型,即內积形式的对偶支持三阶张量机模型如下:
Figure FDA0002280143690000021
其中,Xp与Xq分别是构造的第p个与第q个三阶张量数据,yp与yq是航空发动机数据定义的正常标签或喘振标签,αp,αq是第p个与第q个拉格朗日乘子,d是张量数据的总个数,c是正则化参数;
步骤3:针对具有内积形式的对偶支持三阶张量机模型,通过张量的Tucker分解方法来近似获得具有Tucker分解形式的对偶支持三阶张量机,即Tucker分解形式的对偶支持三阶张量机模型如下:
Figure FDA0002280143690000022
其中,
Figure FDA0002280143690000023
是对张量Xp(Xq)进行Tucker分解后获得的核张量Gp(Gq)中的第k1k2k3(r1r2r3)元素,符号<,>是內积运算,∑是求和运算,Π是求积运算,是对张量Xp(Xq)进行Tucker分解后获得的因子矩阵
Figure FDA0002280143690000025
中的第kt(rt)列向量,n1,n2,n3(s1,s2,s3)是核张量Gp(Gq)在长,宽与高等方向的总层数,t是因子矩阵
Figure FDA0002280143690000026
的下角标,并说明第t个因子矩阵是
Figure FDA0002280143690000027
CN201911138245.XA 2019-11-20 2019-11-20 一种针对航空发动机故障诊断的支持三阶张量机建模方法 Pending CN110765659A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911138245.XA CN110765659A (zh) 2019-11-20 2019-11-20 一种针对航空发动机故障诊断的支持三阶张量机建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911138245.XA CN110765659A (zh) 2019-11-20 2019-11-20 一种针对航空发动机故障诊断的支持三阶张量机建模方法

Publications (1)

Publication Number Publication Date
CN110765659A true CN110765659A (zh) 2020-02-07

Family

ID=69338688

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911138245.XA Pending CN110765659A (zh) 2019-11-20 2019-11-20 一种针对航空发动机故障诊断的支持三阶张量机建模方法

Country Status (1)

Country Link
CN (1) CN110765659A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112431753A (zh) * 2021-01-25 2021-03-02 赛腾机电科技(常州)有限公司 轴向柱塞泵松靴故障的多重量化诊断方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105654110A (zh) * 2015-12-04 2016-06-08 深圳先进技术研究院 一种张量模式下的有监督学习优化方法及系统
CN109359668A (zh) * 2018-09-07 2019-02-19 南京航空航天大学 一种航空发动机多故障并发诊断方法
CN109829243A (zh) * 2019-02-25 2019-05-31 大连理工大学 一种组合式航空发动机高压齿轮泵建模仿真及故障注入方法
US20190347575A1 (en) * 2018-05-08 2019-11-14 International Business Machines Corporation Simulating quantum circuits on a computer using hierarchical storage
CN110457979A (zh) * 2018-05-08 2019-11-15 西北农林科技大学 基于张量Tucker分解模糊控制的柴油机故障诊断方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105654110A (zh) * 2015-12-04 2016-06-08 深圳先进技术研究院 一种张量模式下的有监督学习优化方法及系统
US20190347575A1 (en) * 2018-05-08 2019-11-14 International Business Machines Corporation Simulating quantum circuits on a computer using hierarchical storage
CN110457979A (zh) * 2018-05-08 2019-11-15 西北农林科技大学 基于张量Tucker分解模糊控制的柴油机故障诊断方法
CN109359668A (zh) * 2018-09-07 2019-02-19 南京航空航天大学 一种航空发动机多故障并发诊断方法
CN109829243A (zh) * 2019-02-25 2019-05-31 大连理工大学 一种组合式航空发动机高压齿轮泵建模仿真及故障注入方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
吴振宇 等: "基于Tucker分解的半监督支持张量机", 《计算机科学》 *
许小伟 等: "基于张量Tucker分解的发动机故障诊断", 《中国机械工程》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112431753A (zh) * 2021-01-25 2021-03-02 赛腾机电科技(常州)有限公司 轴向柱塞泵松靴故障的多重量化诊断方法

Similar Documents

Publication Publication Date Title
Lee et al. Real-time anomaly detection framework using a support vector regression for the safety monitoring of commercial aircraft
WO2023123593A1 (zh) 基于变分模态分解和残差网络的航空轴承故障诊断方法
CN107103658B (zh) 基于深度学习与高斯分布的航空发动机气路异常检测方法
JP4875661B2 (ja) 航空機の健全性診断装置及び方法並びにプログラム
US20230214691A1 (en) Fatigue Crack Growth Prediction
CN113469060A (zh) 多传感器融合卷积神经网络航空发动机轴承故障诊断方法
Barad et al. Neural network approach for a combined performance and mechanical health monitoring of a gas turbine engine
CN106169001B (zh) 一种基于气路性能参数图形匹配的航空发动机整机剩余寿命预测方法
US20170024649A1 (en) Anomaly detection system and method for industrial asset
US20170193372A1 (en) Health Management Using Distances for Segmented Time Series
CN113053171B (zh) 一种民机系统风险预警方法及系统
CN102246110A (zh) 用于监测航空发动机的数据标准化方法
CN113168313A (zh) 用于监控飞行器发动机的计算环境系统
CN105956514A (zh) 一种振动数据驱动的直升机旋翼异常检测方法
Tongguang et al. An early fault diagnosis method based on the optimization of a variational modal decomposition and convolutional neural network for aeronautical hydraulic pipe clamps
CN111241694A (zh) 一种基于大数据处理的飞机机群健康评估方法,设备及可读存储介质
Liu et al. Bayesian network method for fault diagnosis of civil aircraft environment control system
Puranik et al. Utilizing energy metrics and clustering techniques to identify anomalous general aviation operations
CN106067032A (zh) 用于数据异常监测的传感器选择方法的评价方法
Pan et al. A method for mechanical fault recognition with unseen classes via unsupervised convolutional adversarial auto-encoder
Pan et al. A meta network pruning framework for remaining useful life prediction of rocket engine bearings with temporal distribution discrepancy
EP3596439A1 (en) Fatigue crack growth prediction
CN110765659A (zh) 一种针对航空发动机故障诊断的支持三阶张量机建模方法
Huang et al. A modified fusion model-based/data-driven model for sensor fault diagnosis and performance degradation estimation of aero-engine
Bisanti et al. Digital twins for aircraft maintenance and operation: A systematic literature review and an IoT-enabled modular architecture

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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20200207

WD01 Invention patent application deemed withdrawn after publication