CN111832216A - 基于eemd-mcnn-gru的滚动轴承剩余使用寿命预测方法 - Google Patents

基于eemd-mcnn-gru的滚动轴承剩余使用寿命预测方法 Download PDF

Info

Publication number
CN111832216A
CN111832216A CN202010290930.0A CN202010290930A CN111832216A CN 111832216 A CN111832216 A CN 111832216A CN 202010290930 A CN202010290930 A CN 202010290930A CN 111832216 A CN111832216 A CN 111832216A
Authority
CN
China
Prior art keywords
layer
convolution
data
prediction
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.)
Granted
Application number
CN202010290930.0A
Other languages
English (en)
Other versions
CN111832216B (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.)
Xinjiang University
Original Assignee
Xinjiang University
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 Xinjiang University filed Critical Xinjiang University
Priority to CN202010290930.0A priority Critical patent/CN111832216B/zh
Publication of CN111832216A publication Critical patent/CN111832216A/zh
Application granted granted Critical
Publication of CN111832216B publication Critical patent/CN111832216B/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
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/04Bearings
    • G01M13/045Acoustic or vibration analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling 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)
  • Software Systems (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computational Linguistics (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Acoustics & Sound (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明公开了一种基于EEMD‑MCNN‑GRU的滚动轴承剩余使用寿命预测方法,包括以下步骤:对滚动轴承的全生命周期振动数据进行预处理;利用集合经验模态构建轴承振动信号在不同时间尺度上的退化特征,以此构建各个滚动轴承的退化特征数据集;集成多尺度卷积神经网络特征提取层和门控制循环单元神经网络时间序列预测层搭建剩余寿命预测模型;对输入预测模型的退化特征数据集进行归一化;对归一化后的数据进行分割处理;将退化特征数据集划分为训练集和测试集,将训练数据作为整个预测网络模型的输入;使用移动平均法对预测模型输出进行平滑处理,输出最优寿命预测结果。本发明为滚动轴承的剩余使用寿命的分析提供了一种更可靠,鲁棒性、泛化性更强的预测方法。

Description

基于EEMD-MCNN-GRU的滚动轴承剩余使用寿命预测方法
技术领域
本发明涉及一种机械装备的故障预测与健康管理领域,具体来说涉及一种基于EEMD-MCNN-GRU的滚动轴承剩余使用寿命预测方法。
背景技术
设备故障预测与健康管理(Prognostics Health Management,PHM)是设备正常运行维护管理中重要的一部分,剩余使用寿命预测(Remaining Useful Life,RUL) 是实现PHM的关键技术之一。滚动轴承广泛的应用于各类机械传动部件中,是保证设备可靠运行的关键部件之一,其剩余寿命是衡量滚动轴承性能的重要指标,因此开展滚动轴承寿命预测研究具有重要的意义。设备剩余寿命预测方法发展至今,主要分为失效机理(Physics-of-Failure,PoF)分析和数据驱动 (Data-Driven,DD)这两种方法。
失效机理分析主要通过物理方法研究零部件变形、裂纹和腐蚀等退化现象, 由此建立设备的失效机理模型,结合特定设备的经验知识和相关检测方法实现设备的剩余寿命预测。WHLER早在1847年提出的S-N疲劳寿命曲线开创了疲劳破坏的经典强度理论。后续,GRIFFITH通过研究裂纹扩展规律提出了裂纹扩展能量理论,给剩余寿命预测探究提供了新的思路。在21世纪前后对疲劳断裂理论的研究不断完善,而且在考虑到腐蚀的影响后提出了腐蚀寿命预测技术。由于基于失效机理的预测方法虽然能够较为准确的预测设备的剩余使用寿命,但是由于其必须根据特定工况和设备具体分析,此外随着装备复杂性的不断增大失效机理难以获得,因此,此类方法的推广受到了一定的限制。
基于数据驱动的方法可分为统计数据驱动方法、基于浅层的机器学习 (MachineLearning,ML)的方法和基于深度学习(Deep Learning,DL)的方法。统计数据驱动方法中应用较多的方法有基于Gamma过程的方法和基于Wiener 过程的方法。统计数据驱动的寿命预测结果能反映机械产品寿命的一般规律和整体特性,但是在实际工程中,其假设的模型可能会使预测精度受到限制并且需要大量试验和数据的积累作为支撑。
相较于基于机理分析的方法和统计数据驱动方法,用浅层机器学习的方法和深度学习的方法在一定程度上克服机理模型难以建立和假设的随机退化过程带来的限制,而且对输入数据的要求也较为宽泛,可以更加全面的考虑退化因素。在通过ML进行剩余寿命预测的过程中无需过多的先验知识,可实现在线监测,得到预测结果。
浅层机器学习的方法被应用到寿命预测中也取得了一定的发展。浅层机器学习中以支持向量机(Support Vector Machine,SVM)为主。SVM首次提出是被用来解决ML的分类和回归问题,适用于分析小样本和多维化数据。发表于《.Energies》2017,11(1)中的《emaining useful life estimation of aircraft engines using a modifiedsimilarity and supporting vector machine(SVM)approach》使用 SVM和改进的相似性理论相结合应用于航空发动机的剩余寿命预测。但是SVM 存在一些如随着样本的增大,线性增加,出现过拟合、计算时间增加等固有缺陷,所以发表于《Mechanical Systems&SignalProcessing》2012,31(8)中的《Combining relevance vector machines and exponentialregression for bearing residual life estimation》使用了能较为有效克服过拟合问题的相关向量机 (Relevance vector machine,RVM)和指数回归的方法结合预测推力轴承的剩余寿命。发表于《计算机集成制造系统》2018,24(11)中的《基于支持向量数据描述的剩余寿命预测方法》使用支持向量描述(Support Vector Data Description, SVDD)的超球面与测试样本点的距离计算退化指数,并实现了轴承剩余寿命的预测。
在现代工业中监测设备运行时可采集各种类型的大量信号数据,传统的浅层机器学习要依靠专家先验知识进行特征提取,泛化性能差。深度学习的模型可以映射更为复杂的非线性关系,具有强大的特征提取能力。因此深度学习作为更加智能的方式为剩余寿命预测开拓了新方向。深度学习中经典的深度神经网络(Deep neural network,DNN)的方法、卷积神经网络(Convolutional neural network,CNN)、深度置信网络(Deep beliefnetwork,DBN)]以及递归神经网络 (Recurrent neural network,RNN)等都被应用于剩余寿命的预测。发表于《山东大学学报(工学版)》2017,47(5)中的《基于深度学习的缓变故障早期诊断及寿命预测》用DNN进行特征提取,然后进行PCA降维得到退化曲线并以此进行寿命预测。发表于《计算机集成制造系统》2019,25(07)中的《基于Autoencoder-BLSTM 的涡扇发动机剩余寿命预测》将CNN用于特征提取并直接进行了寿命预测,该文将轴承的全生命周期的振动数据进行傅里叶变换,将频域特征作为CNN的输入预测结果有效,证明了CNN在特征提取中的能力。因此,发表于《IEEE Access》2018:1-1中的《Prediction of BearingRemaining Useful Life with Deep Convolution Neural Network》将CNN用于特征提取并直接进行了轴承RUL预测,其将轴承的全生命周期的振动信号进行傅里叶变换,再将频域特征作为CNN的输入且预测结果有效。DBN是一种半监督的深度学习算法,能通过预训练能减少训练时间,并且能有效克服随机初始化权值造成的局部最优问题。由于RNN 及其改进结构能够有效处理时序问题,所以发表于《Reliability Engineering& System Safety》2019,185中的《Gated recurrent unit based recurrent neural network forremaining useful life prediction of nonlinear deterioration process》通过门控制循环单元神经网络(Gated Recurrent Unit,GRU)对航空发动机的剩余寿命进行了预测。发表于《计算机集成制造系统》2019,25(07)中的《基于Autoencoder-BLSTM 的涡扇发动机剩余寿命预测》通过结合自编码器和双向长短时记忆神经网络 (Bidirectional LongShort-Term Memory,BLSTM)实现了对航空发动机寿命的预测。通过上述研究可知,深度学习在剩余寿命预测中,成效显著、应用前景较广阔。但以往的研究中,大多数只是在单一的时间尺度上对样本数据进行特征提取,这将这导致特征的提取不完全。
发明内容
本发明所要解决的技术问题是:提供一种能有效解决现有的滚动轴承剩余使用寿命预测研究中存在的特征提取不完备、预测模型鲁棒性不强、预测精度较低问题的基于EEMD-MCNN-GRU的滚动轴承剩余使用寿命预测方法。
为解决上述技术问题,本发明所采用的技术方案为:一种基于EEMD- MCNN-GRU的滚动轴承剩余使用寿命预测方法,包括以下步骤:
S1:基于EEMD的滚动轴承退化特征数据集构建
收集滚动轴承从开始使用到发生故障的全生命周期振动数据,然后按其采样周期,将采用集合经验模态分解EEMD将原始振动数据分解为多个限制带宽的不同频段的IMF分量,如式(1)所表示,构建轴承振动信号在不同时间尺度上的退化特征,选择反映轴承退化趋势的至少前六个IMF分量作为样本,以此构建各个滚动轴承的退化特征数据集;
Figure RE-GDA0002664462640000041
式中,X(t)为原始信号,N为分解的IMFs数量,IMFi(t)为第i个固有模式函数,R(t)为余量;
S2:退化特征数据集预处理
S21:对所有输入预测模型的退化特征数据集用式(2)进行归一化,将数据压缩在[0,1]之间;
Figure RE-GDA0002664462640000042
式中,xi,j表示数据集中第i个IMF分量的第j个数据点,
Figure RE-GDA0002664462640000043
Figure RE-GDA0002664462640000044
表示第 i个IMF分量的最小值和最大值;
Figure RE-GDA0002664462640000045
表示对xi,j归一化后的值;
S22:对归一化后的数据进行分割处理以满足预测模型的输入要求;
特征集的长度为N,将特征集中的数据按照一次采样所采数据点的个数确定为分割步长L,即一次EEMD分解的长度;数据集的特征纬度为n,则模型的输入维度为(N-L,L,n)的三维张量;
S23:将退化特征数据集划分为训练集和测试集,将训练数据作为整个网络模型的输入;
S3:MCNN-GRU预测模型的搭建
S31:预测模型的多尺度卷积神经网络MCNN特征提取模块搭建;
预测模型的多尺度特征提取模块包含卷积层,池化层和连接层;
首先,由若干个包含不同尺寸卷积核的一维卷积神经网络1D-CNN组成多尺度卷积层,同时在不同尺度上进行全局的特征提取,其中每个1D-CNN连接一个平均池化层,用以缩减特征;其次,在连接层将多尺度卷积层提取的特征进行合并,以便输入后续网络结构;最后,通过1D-CNN和最大池化操作完成提取特征;
S32:GRU时序预测模块搭建
将提取的特征输入GRU预测模块进行时间序列预测,再通过全连接层输出寿命预测结果:GRU网络模型的xt,ht分别表示t时刻的GRU网络的输入和输出,其中ht由下列公式迭代计算:
zt=σ(Whzht-1+Wxzxt+bz) (6)
rt=σ(Whrht-1+Wxrxt+br) (7)
ht=tanh(rt*Whhht-1+Wxhxt+bh) (8)
Figure RE-GDA0002664462640000051
式中,zt,rt分别为更新门和重置门的输出,Whz,分别为t-1时刻输出和 t时刻输入到更新门的权重矩阵,Whr、Wxr分别为t-1时刻输出和t时刻输入到重置门的权重矩阵,这两个门输出都经过一个sigmoid激活函数,其值域为 [0,1];
Figure RE-GDA0002664462640000052
ht分别为候选隐藏状态和隐藏状态,它们分别控制上一时刻隐藏状态ht-1中信息的保留和遗忘;Whh,Wxh分别为t-1时刻输出和t时刻输入到候选隐藏状态的权重矩阵,候选隐藏状态的输出经过一个tanh函数;bz,br和bh表示更新门、重置门和候选隐藏状态的偏置;
S4:MCNN-GRU预测模型的训练
S41:初始化网络参数
初始化网络参数:数据分割步长L,特征纬度n,模型网络结构,初始网络权重θ,初始学习率η,指数衰减参数γ,防止0分母错误参数ε,训练次数E,批尺度B;
S42:确定模型的损失函数
模型的损失函数为均方误差(mean square error,MSE),训练过程中以MSE 最小为优化目标:
Figure RE-GDA0002664462640000061
式中:yi为模型的标签值,oi为模型的实际输出;
S43:对样本设定数据标签
数据标签使用理想寿命退化曲线,即将轴承从初始运行到完全故障的时间作为最终寿命然后每隔一次采样周期递减1构造一维数组,然后使用滚动轴承经验最大寿命,把所有数据集中的数据归一化到[0,1]作为其剩余寿命,这样就对每一个样本设置了对应的标签;
S44:网络权重更新
使用RMSProp算法更新网络权重,算法的更新步骤为:
Gt=γ·Gt-1+(1-γ)gt·gt (11)
Figure RE-GDA0002664462640000062
式中Gt表示t时刻的梯度平方和;γ表示指数衰减参数;gt表示t时刻的梯度;θt表示t时刻的参数;η表示学习率;ε是为了避免式(12)分母为0的防止0分母错误参数;
S5:基于移动平均法的预测结果局部波动平滑处理
将测试数据输入训练好的预测模型,得到初步预测结果,采用移动平均法,如式(13)对网络模型的初步预测结果进行平滑处理,去除局部波动对预测性能的影响,获得在滚动轴承全寿命周期上明显的退化趋势:
yt'=(yt-n+yt-n+1+yt-n+2+…+yt)/(n+1) (13)
式中yt'表示在t时刻平滑之后的值,yt表示其t时刻未平滑的数据。
作为一种优选的方案,所述步骤S4中MCNN-GRU模型的具体训练算法过程为:
算法:MCNN-GRU模型的训练算法;
输入:经EEMD分解的轴承退化振动数据集;
初始化参数:数据分割步长L,特征纬度n,模型网络结构,初始网络权重θ,初始学习率η,指数衰减参数γ,防止0分母错误参数ε,训练次数E,批尺度B;
训练过程:
a、数据标准化
b、数据分割
c、
a)For E∈N do
b)将标准化后的数据输入MCNN-GRU通过向前传播得到输出oi
c)计算损失函数
Figure RE-GDA0002664462640000071
d)使用RMSProp算法对网络权重进行更新;
e)end for
d、输出:训练完成的MCNN-GRU模型。
作为一种优选的方案,每一个所述1D-CNN包括输入层、卷积层、池化层、全连接层、输出层;
卷积层:卷积核对输入信号的局部区域进行卷积运算,如式(3),得到相应的特征,卷积核逐步在输入数据上进行卷积然后将每次卷积的结果依次排列成一个输出向量:
Figure RE-GDA0002664462640000072
式中,yl(i,k)表示第l层的第i个卷积核对第k个卷积区域的结果;
Figure RE-GDA0002664462640000081
表示第l层的第i个卷积核的第j个权值;
Figure RE-GDA0002664462640000082
表示第l层的第i个卷积核对第k个卷积区域的卷积对象;
Figure RE-GDA0002664462640000083
表示第l层的第i个卷积核的偏置量;s表示卷积尺寸;
池化层:通过池化算子对卷积结果进行池化操作,以缩减模型大小,提高计算速度及所提取特征的鲁棒性;
全连接层:将池化层的输出映射到对应的标签;该层的正向传播如下式 (4):
Figure RE-GDA0002664462640000084
式中,
Figure RE-GDA0002664462640000085
为第l层第i个神经元的输出值;
Figure RE-GDA0002664462640000086
为第l层第i个神经元与第 l+1层第j个神经元之间的权值;
Figure RE-GDA0002664462640000087
为第l层所有神经元的偏置值;W为对应的节点个数;
用多个不同尺寸卷积核的卷积层同时对原始特征信号进行卷积操作实现多尺度时间上的特征提取,以表征更加完备的特征一个多尺度1D-CNN结构的公式如(5)式,
Figure RE-GDA0002664462640000088
式中,
Figure RE-GDA0002664462640000089
表示第l层第m个卷积操作的第i个卷积核对第k个卷积区域的结果;
Figure RE-GDA00026644626400000810
表示第l层第m个卷积操作的第i个卷积核;
Figure RE-GDA00026644626400000811
表示第l层第m 个卷积操作的卷积核对应的第k个卷积区域;
Figure RE-GDA00026644626400000812
表示第l层第m个卷积操作的第i个卷积核的偏置。
本发明的有益效果是:
(1)EEMD被广泛的使用在信号处理和故障诊断中,可以自适应的在不同时间尺度上完成分解,并对故障特征反应较为敏感,本发明利用EEMD构建轴承振动信号在不同时间尺度上的退化特征,通过从不同时间尺度上对滚动轴承原始振动数据进行特征提取以提高预测模型的有效性和泛化能力。
(2)为了提取轴承在不同时间尺度上的退化特征,本发明使用多个不同卷积核尺寸的CNN同时对输入数据在不同尺度上并行提取特征,MCNN这种多尺度并行结构能使全局和局部信息同步,对RUL预测有显著贡献的特征。
(3)改进模型门控制循环单元神经网络(GRU)作为一种可有效解决RNN 梯度爆炸和梯度消失问题的变体,被成功地广泛应用于具有时序依赖关系的问题处理中,而且较另一个类似的RNN变体长短期记忆神经网络(Long Short-term Memory,LSTM)需要更少的参数,可以降低训练时的计算成本。滚动轴承的退化过程具有明显的时序性,因此将GRU应用于该问题中将有效的学习特征中的时序信息。因此,结合MCNN与GRU构建成同时具有多尺度特征提取和时序关系处理特性的网络模型,可用做剩余寿命预测的多特性预测模型。
附图说明
图1是本发明的实施例的具体流程图。
图2是MCNN-GRU网络预测模型结构示意图。
图3是多尺度一维卷积示意图。
图4是实施例所用法国FEMTO-ST研究所公开的标准数据集中7个滚动轴承振动时域图。
图5是Bearing1-1EEMD分解的前六个IMF分量。
图6是Bearing1-1分解的IMF频谱图。
图7是实施例中基于用本发明所提方法对4个测试滚动轴承RUL预测的结果。
具体实施方式
下面结合附图,详细描述本发明的具体实施方案。
本发明提出的一种基于EEMD-MCNN-GRU的滚动轴承剩余使用寿命预测方法实施例的具体流程图如附图2所示,包括以下步骤:
S1:基于EEMD的滚动轴承退化特征数据集构建
收集7个滚动轴承从开始使用到发生故障的全生命周期振动数据,然后按其采样周期,将采用EEMD将原始振动数据分解为多个限制带宽的不同频段的IMF分量,如式(1)所表示,大幅度降低模态混叠现象产生的数据波动对算法的预测性能所造成的影响。构建轴承振动信号在不同时间尺度上的退化特征,选择反映轴承退化趋势的前六个IMF分量作为样本,以此构建各个滚动轴承的退化特征数据集。
Figure RE-GDA0002664462640000101
式中,X(t)为原始信号,N为分解的IMFs数量,IMFi(t)为第i个固有模式函数,R(t)为余量。
将每个轴承的振动数据按采样样本进行EEMD分解,将7个轴承分解后前六个具有主要信息的IMF分量作为退化特征数据集。以Bearing1-1分解后的结果为例,对其分解后的前六个具有主要信息的IMF分量进行分析,如附图4 所示是前六个IMF的时域图。从附图4中可以明显地看出IMF1和IMF2存在一些噪音和冲击成分,IMF3和IMF4是次高频成分但是噪声明显减少,IMF5 和IMF6的幅值呈现明显逐渐增大趋势,能较好的表征轴承的退化过程。
取Bearing1-1退化前、中、后、期的前六个IMF分量进行傅里叶变换,其频谱图如附图5所示,a、b、c代表轴承退化的前、中、后三个时期,数字代表第几IMF个分量。可以从IMF1看出其中的高频成分不断向高频移动,说明其中的冲击和噪音在加剧。IMF3在退化开始时出现故障频率。通过分析可以看出,固有模式函数中IMF1、IMF2主要包含的是高频冲击成分,IMF5主要包含的是旋转引起的低频的工频分量。从以上分析可以看出各个阶段的集合经验模式分解是合理的,并且在保留了信号的时域和频域信息的同时表达出了信号的退化特征。
S2:退化特征数据集预处理
S21:对所有输入预测模型的退化特征数据集用式(2)进行归一化,将数据压缩在[0,1]之间,以便得到更好的训练精度。
Figure RE-GDA0002664462640000111
式中,xi,j表示数据集中第i个IMF分量的第j个数据点,
Figure RE-GDA0002664462640000112
Figure RE-GDA0002664462640000113
表示第 i个IMF分量的最小值和最大值。
Figure RE-GDA0002664462640000114
表示对xi,j归一化后的值。
S22:为了满足预测模型的输入要求,对归一化后的数据进行分割处理。
特征集的长度为N,将特征集中的数据按照一次采样所采数据点的个数确定为分割步长L,即一次EEMD分解的长度,数据集的特征纬度为n,则模型的输入维度为(N-L,L,n)的三维张量。
S23:将退化特征数据集划分为训练集和测试集,将训练数据作为整个预测网络模型的输入。
Bearing1-1、Bearing1-2和Bearing1-5的EEMD分解结果(即退化特征数据集)作为训练集,将Bearing1-3、Bearing1-4、Bearing1-6和Bearing1-7的 EEMD分解结果(即退化特征数据集)作为测试集。这样划分数据集是为了平衡样本数据在轴承突发失效和退化失效两种类型的分布,让预测网络模型能够学习到更加完善的信息。
S3:MCNN-GRU预测模型的搭建
MCNN-GRU预测模型结构主要可分为基于MCNN特征提取模块,基于 GRU时间序列预测模块和Dense全连接模块。其网络结构如附图6所示,具体描述如下:
S31:预测模型的MCNN特征提取模块搭建。
预测模型的多尺度特征提取模块包含卷积层,池化层和连接层。首先,由若干个包含不同尺寸卷积核的1D-CNN组成多尺度卷积层,同时在不同尺度上进行全局的特征提取,其中每个1D-CNN连接一个平均池化层,用以缩减特征;其次,在连接层将多尺度卷积层提取的特征进行合并,以便输入后续网络结构;最后,通过1D-CNN和最大池化操作完成提取特征。
S32:GRU时序预测模块搭建
将提取的特征输入GRU预测模块进行时间序列预测,再通过全连接层输出寿命预测结果。
S4:MCNN-GRU预测模型的训练
将Bearing1-1、Bearing1-2和Bearing1-5的EEMD分解结果(即退化特征数据集)作为训练集,输入MCNN-GRU预测模型中得到最优的网络结构参数。
S41:MCNN模型
卷积层:用多个不同尺寸卷积核的卷积层同时对原始特征信号进行卷积操作实现多尺度时间上的特征提取,如式(3),卷积核逐步在输入数据上进行卷积然后将每次卷积的结果依次排列成一个输出向量,表征更加完备的特征:
Figure RE-GDA0002664462640000121
式中,
Figure RE-GDA0002664462640000122
表示第l层第m个卷积操作的第i个卷积核对第k个卷积区域的结果;
Figure RE-GDA0002664462640000123
表示第l层第m个卷积操作的第i个卷积核;
Figure RE-GDA0002664462640000124
表示第l层第m 个卷积操作的卷积核对应的第k个卷积区域;
Figure RE-GDA0002664462640000125
表示第l层第m个卷积操作的第i个卷积核的偏置。
池化层:通过池化算子对卷积结果进行池化操作,以缩减模型大小,提高计算速度及所提取特征的鲁棒性。
全连接层:将池化层的输出映射到对应的标签。该层的正向传播如式(4):
Figure RE-GDA0002664462640000126
式中,
Figure RE-GDA0002664462640000127
为第l层第i个神经元的输出值;
Figure RE-GDA0002664462640000128
为第l层第i个神经元与第l+1层第j个神经元之间的权值;bi l为第l层所有神经元的偏置值;W为对应的节点个数。
S42:GRU时序预测
将提取的特征输入GRU预测模块进行时间序列预测,再通过全连接层输出寿命预测结果。GRU网络模型的xt,ht分别表示t时刻的GRU网络的输入和输出,其中ht由公式(5)—(8)迭代计算:
zt=σ(Whzht-1+Wxzxt+bz) (5)
rt=σ(Whrht-1+Wxrxt+br) (6)
ht=tanh(rt*Whhht-1+Wxhxt+bh) (7)
Figure RE-GDA0002664462640000131
式中,zt,rt分别为更新门和重置门的输出,Whz,分别为t-1时刻输出和 t时刻输入到更新门的权重矩阵,Whr、Wxr分别为t-1时刻输出和t时刻输入到重置门的权重矩阵,这两个门输出都经过一个sigmoid激活函数,其值域为 [0,1];
Figure RE-GDA0002664462640000132
ht分别为候选隐藏状态和隐藏状态,它们分别控制上一时刻隐藏状态ht-1中信息的保留和遗忘;Whh,Wxh分别为t-1时刻输出和t时刻输入到候选隐藏状态的权重矩阵,候选隐藏状态的输出经过一个tanh函数;bz,br和bh表示更新门、重置门和候选隐藏状态的偏置。
S42:初始化网络参数
初始化网络参数:数据分割步长L,特征纬度n,模型网络结构,初始网络权重θ,初始学习率η,指数衰减参数γ,防止0分母错误参数ε,训练次数E,批尺度B。
S43:确定模型的损失函数
模型的损失函数为均方误差(mean square error,MSE),训练过程中以MSE 最小为优化目标。
Figure RE-GDA0002664462640000141
式中:yi为模型的标签值,oi为模型的实际输出。
S44:对样本设定数据标签
数据标签使用理想寿命退化曲线,即将轴承从初始运行到完全故障的时间作为最终寿命然后每隔一次采样周期递减1构造一维数组,然后使用滚动轴承经验最大寿命2803为最大值,把所有数据集中的数据归一化到[0,1]作为其剩余寿命,这样就对每一个样本设置了对应的标签。
S45:网络权重更新
使用RMSProp算法更新网络权重,算法的更新步骤为:
Gt=γ·Gt-1+(1-γ)gt·gt (10)
Figure RE-GDA0002664462640000142
式中Gt表示t时刻的梯度平方和;γ表示指数衰减参数;gt表示t时刻的梯度;θt表示t时刻的参数;η表示学习率;为了避免式(11)分母为0,使用一个很小的正数ε作为修正。RMSProp算法使用指数衰减平均调整历史梯度的比重,能自适应调节学习率,能加速模型快速收敛。
优化后的网络结构模型参数如附表1所示,MCNN模块的第一层由3个拥有不同尺寸卷积核的并行1D-CNN实现,其通道数均为6,卷积核尺寸分别为 4、8、16,同时在不同尺度上进行全局的特征提取。其中每个1D-CNN后连接一个尺寸为3平均池化层,用以缩减特征。其次,在Concatenate层将提取的特征进行合并,以便输入后续1D-CNN网络结构,进行两次卷积和最大池化操作,第一次卷积操作的卷积核尺寸为16,通道数为8,最大池化尺寸为4,第二次卷积操作的卷积核尺寸为8,通道数为8的卷积层,最大池化尺寸为4。该模块中所有1D-CNN均采用tanh激活函数,完成特征提取操作。将提取的特征连接一个输出8个节点的GRU结构,学习其时序关系。最后通过一个输出为1激活函数为sigmoid的Dense全连接层得到RUL预测结果。
表1网络结构模型参数
Figure RE-GDA0002664462640000151
S5:基于移动平均法的预测结果局部波动平滑处理
将Bearing1-3、Bearing1-4、Bearing1-6和Bearing1-7的EEMD分解结果 (即退化特征数据集)作为测试集,输入步骤S4得到训练好的MCNN-GRU 预测模型中,得到初步预测结果。
然后,采用移动平均法,如式(12),对网络模型的初步预测结果进行平滑处理,去除局部波动对预测性能的影响,获得在滚动轴承全寿命周期上明显的退化趋势。
y't=(yt-n+yt-n+1+yt-n+2+…+yt)/(n+1) (12)
式中y't表示在t时刻平滑之后的值,yt表示其t时刻未平滑的数据。
得到4个测试集轴承的RUL最终预测结果,如附图7所示。
并将本结果与常用的几种方法对比,结果如附表2所示。
表2本发明与其它常用方法的对比结果
Figure RE-GDA0002664462640000161
将本发明所提出的方法与常用的几种方法:结合自编码器和多层感知机网络模型(Autoencoder DNN)、SVM、DNN-FSPS和传统CNN进行对比分析,结果如附表2所示。由表2所知:(1)通过对比分析本文方法(MCNN-GRU) 与传统的CNN预测方法,发现本文平均预测误差为0.013,而传统的CNN平均预测误差为0.065,相较于传统CNN,本文所提方法平均预测精度提高了约 80%,说明本文所提的多尺度结构在轴承退化特征提取中效果较好;(2)通过与其它模型对比,本文所提方法的平均预测精度得到了明显提升。
上述的实施例仅例示性说明本发明创造的原理及其功效,以及部分运用的实施例,而非用于限制本发明;应当指出,对于本领域的普通技术人员来说,在不脱离本发明创造构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。

Claims (3)

1.一种基于EEMD-MCNN-GRU的滚动轴承剩余使用寿命预测方法,包括以下步骤:
S1:基于EEMD的滚动轴承退化特征数据集构建
收集滚动轴承从开始使用到发生故障的全生命周期振动数据,然后按其采样周期,将采用集合经验模态分解EEMD将原始振动数据分解为多个限制带宽的不同频段的IMF分量,如式(1)所表示,构建轴承振动信号在不同时间尺度上的退化特征,选择反映轴承退化趋势的至少前六个IMF分量作为样本,以此构建各个滚动轴承的退化特征数据集;
Figure FDA0002450355930000011
式中,X(t)为原始信号,N为分解的IMFs数量,IMFi(t)为第i个固有模式函数,R(t)为余量;
S2:退化特征数据集预处理
S21:对所有输入预测模型的退化特征数据集用式(2)进行归一化,将数据压缩在[0,1]之间;
Figure FDA0002450355930000012
式中,xi,j表示数据集中第i个IMF分量的第j个数据点,
Figure FDA0002450355930000013
Figure FDA0002450355930000014
表示第i个IMF分量的最小值和最大值;
Figure FDA0002450355930000015
表示对xi,j归一化后的值;
S22:对归一化后的数据进行分割处理以满足预测模型的输入要求;
特征集的长度为N,将特征集中的数据按照一次采样所采数据点的个数确定为分割步长L,即一次EEMD分解的长度;数据集的特征纬度为n,则模型的输入维度为(N-L,L,n)的三维张量;
S23:将退化特征数据集划分为训练集和测试集,将训练数据作为整个网络模型的输入;
S3:MCNN-GRU预测模型的搭建
S31:预测模型的多尺度卷积神经网络MCNN特征提取模块搭建;
预测模型的多尺度特征提取模块包含卷积层,池化层和连接层;
首先,由若干个包含不同尺寸卷积核的一维卷积神经网络1D-CNN组成多尺度卷积层,同时在不同尺度上进行全局的特征提取,其中每个1D-CNN连接一个平均池化层,用以缩减特征;其次,在连接层将多尺度卷积层提取的特征进行合并,以便输入后续网络结构;最后,通过1D-CNN和最大池化操作完成提取特征;
S32:GRU时序预测模块搭建
将提取的特征输入GRU预测模块进行时间序列预测,再通过全连接层输出寿命预测结果:GRU网络模型的xt,ht分别表示t时刻的GRU网络的输入和输出,其中ht由下列公式迭代计算:
zt=σ(Whzht-1+Wxzxt+bz) (6)
rt=σ(Whrht-1+Wxrxt+br) (7)
ht=tanh(rt*Whhht-1+Wxhxt+bh) (8)
Figure FDA0002450355930000021
式中,zt,rt分别为更新门和重置门的输出,Whz,分别为t-1时刻输出和t时刻输入到更新门的权重矩阵,Whr、Wxr分别为t-1时刻输出和t时刻输入到重置门的权重矩阵,这两个门输出都经过一个sigmoid激活函数,其值域为[0,1];
Figure FDA0002450355930000022
ht分别为候选隐藏状态和隐藏状态,它们分别控制上一时刻隐藏状态ht-1中信息的保留和遗忘;Whh,Wxh分别为t-1时刻输出和t时刻输入到候选隐藏状态的权重矩阵,候选隐藏状态的输出经过一个tanh函数;bz,br和bh表示更新门、重置门和候选隐藏状态的偏置;
S4:MCNN-GRU预测模型的训练
S41:初始化网络参数
初始化网络参数:数据分割步长L,特征纬度n,模型网络结构,初始网络权重θ,初始学习率η,指数衰减参数γ,防止0分母错误参数ε,训练次数E,批尺度B;
S42:确定模型的损失函数
模型的损失函数为均方误差(mean square error,MSE),训练过程中以MSE最小为优化目标:
Figure FDA0002450355930000031
式中:yi为模型的标签值,oi为模型的实际输出;
S43:对样本设定数据标签
数据标签使用理想寿命退化曲线,即将轴承从初始运行到完全故障的时间作为最终寿命然后每隔一次采样周期递减1构造一维数组,然后使用滚动轴承经验最大寿命,把所有数据集中的数据归一化到[0,1]作为其剩余寿命,这样就对每一个样本设置了对应的标签;
S44:网络权重更新
使用RMSProp算法更新网络权重,算法的更新步骤为:
Gt=γ·Gt-1+(1-γ)gt·gt (11)
Figure FDA0002450355930000032
式中Gt表示t时刻的梯度平方和;γ表示指数衰减参数;gt表示t时刻的梯度;θt表示t时刻的参数;η表示学习率;ε是为了避免式(12)分母为0的防止0分母错误参数;
S5:基于移动平均法的预测结果局部波动平滑处理
将测试数据输入训练好的预测模型,得到初步预测结果,采用移动平均法,如式(13)对网络模型的初步预测结果进行平滑处理,去除局部波动对预测性能的影响,获得在滚动轴承全寿命周期上明显的退化趋势:
y′t=(yt-n+yt-n+1+yt-n+2+…+yt)/(n+1) (13)
式中y′t表示在t时刻平滑之后的值,yt表示其t时刻未平滑的数据。
2.如权利要求1所述的一种基于EEMD-MCNN-GRU的滚动轴承剩余使用寿命预测方法,其特征在于:上述步骤S4中MCNN-GRU模型的具体训练算法过程为:
算法:MCNN-GRU模型的训练算法;
输入:经EEMD分解的轴承退化振动数据集;
初始化参数:数据分割步长L,特征纬度n,模型网络结构,初始网络权重θ,初始学习率η,指数衰减参数γ,防止0分母错误参数ε,训练次数E,批尺度B;
训练过程:
a、数据标准化
b、数据分割
c、
a)For E∈N do
b)将标准化后的数据输入MCNN-GRU通过向前传播得到输出oi;
c)计算损失函数
Figure FDA0002450355930000041
d)使用RMSProp算法对网络权重进行更新;
e)end for
d、输出:训练完成的MCNN-GRU模型。
3.如权利要求1或2所述的一种基于EEMD-MCNN-GRU的滚动轴承剩余使用寿命预测方法,其特征在于:每一个所述的1D-CNN都包括输入层、卷积层、池化层、全连接层、输出层;
卷积层:卷积核对输入信号的局部区域进行卷积运算,如式(3),得到相应的特征,卷积核逐步在输入数据上进行卷积然后将每次卷积的结果依次排列成一个输出向量:
Figure FDA0002450355930000051
式中,yl(i,k)表示第l层的第i个卷积核对第k个卷积区域的结果;
Figure FDA0002450355930000052
表示第l层的第i个卷积核的第j个权值;
Figure FDA0002450355930000053
表示第l层的第i个卷积核对第k个卷积区域的卷积对象;
Figure FDA0002450355930000054
表示第l层的第i个卷积核的偏置量;s表示卷积尺寸;
池化层:通过池化算子对卷积结果进行池化操作,以缩减模型大小,提高计算速度及所提取特征的鲁棒性;
全连接层:将池化层的输出映射到对应的标签;该层的正向传播如下式(4):
Figure FDA0002450355930000055
式中,
Figure FDA0002450355930000056
为第l层第i个神经元的输出值;
Figure FDA0002450355930000057
为第l层第i个神经元与第l+1层第j个神经元之间的权值;
Figure FDA0002450355930000058
为第l层所有神经元的偏置值;W为对应的节点个数;
用多个不同尺寸卷积核的卷积层同时对原始特征信号进行卷积操作实现多尺度时间上的特征提取,以表征更加完备的特征一个多尺度1D-CNN结构的公式如(5)式,
Figure FDA0002450355930000059
式中,
Figure FDA00024503559300000510
表示第l层第m个卷积操作的第i个卷积核对第k个卷积区域的结果;
Figure FDA00024503559300000511
表示第l层第m个卷积操作的第i个卷积核;
Figure FDA00024503559300000512
表示第l层第m个卷积操作的卷积核对应的第k个卷积区域;
Figure FDA00024503559300000513
表示第l层第m个卷积操作的第i个卷积核的偏置。
CN202010290930.0A 2020-04-14 2020-04-14 基于eemd-mcnn-gru的滚动轴承剩余使用寿命预测方法 Active CN111832216B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010290930.0A CN111832216B (zh) 2020-04-14 2020-04-14 基于eemd-mcnn-gru的滚动轴承剩余使用寿命预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010290930.0A CN111832216B (zh) 2020-04-14 2020-04-14 基于eemd-mcnn-gru的滚动轴承剩余使用寿命预测方法

Publications (2)

Publication Number Publication Date
CN111832216A true CN111832216A (zh) 2020-10-27
CN111832216B CN111832216B (zh) 2022-07-12

Family

ID=72914089

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010290930.0A Active CN111832216B (zh) 2020-04-14 2020-04-14 基于eemd-mcnn-gru的滚动轴承剩余使用寿命预测方法

Country Status (1)

Country Link
CN (1) CN111832216B (zh)

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112541576A (zh) * 2020-12-14 2021-03-23 四川翼飞视科技有限公司 Rgb单目图像的生物活体识别神经网络及其构建方法
CN112560201A (zh) * 2020-12-24 2021-03-26 新疆大学 复杂载荷工况下风机叶片复合材料可靠性分析方法
CN112580263A (zh) * 2020-12-24 2021-03-30 湖南工业大学 基于时空特征融合的涡扇发动机剩余使用寿命预测方法
CN112785092A (zh) * 2021-03-09 2021-05-11 中铁电气化局集团有限公司 一种基于自适应深层特征提取的道岔剩余寿命预测方法
CN113051689A (zh) * 2021-04-25 2021-06-29 石家庄铁道大学 基于卷积门控循环网络的轴承剩余使用寿命预测方法
CN113052365A (zh) * 2021-02-26 2021-06-29 浙江工业大学 一种基于mswr-lrcn的旋转类机械寿命预测方法
CN113052060A (zh) * 2021-03-22 2021-06-29 六盘水师范学院 基于数据增强的轴承剩余寿命预测方法、装置及电子设备
CN113065283A (zh) * 2021-03-26 2021-07-02 深圳技术大学 一种电池寿命预测方法、系统、电子装置及存储介质
CN113096378A (zh) * 2021-02-23 2021-07-09 北京航空航天大学 基于深度集合经验模态分解的高速公路od预测方法
CN113205050A (zh) * 2021-05-09 2021-08-03 西北工业大学 一种基于gru-ae网络的舰船辐射噪声线谱提取方法
CN113449465A (zh) * 2021-06-15 2021-09-28 南京航空航天大学 一种用于滚动轴承的寿命预测方法
CN113486578A (zh) * 2021-06-28 2021-10-08 北京科技大学 一种工业过程中设备剩余寿命的预测方法
CN113807005A (zh) * 2021-08-12 2021-12-17 北京工业大学 基于改进fpa-dbn的轴承剩余寿命预测方法
CN114154743A (zh) * 2021-12-13 2022-03-08 四川大学 基于vetmrrn的空间滚动轴承剩余寿命预测方法
CN114298210A (zh) * 2021-12-24 2022-04-08 江苏国科智能电气有限公司 一种多视角融合风电齿轮箱轴承剩余寿命预测方法
CN114564877A (zh) * 2021-12-29 2022-05-31 西安交通大学 一种滚动轴承寿命预测方法、系统、设备及可读存储介质
CN114861527A (zh) * 2022-04-15 2022-08-05 南京工业大学 一种基于时间序列特征的锂电池寿命预测方法
CN115446663A (zh) * 2022-10-14 2022-12-09 常州先进制造技术研究所 基于物理指导的深度学习网络的刀具磨损状态监测方法及应用
CN115660198A (zh) * 2022-11-07 2023-01-31 中国科学院空间应用工程与技术中心 滚动轴承剩余使用寿命预测方法
CN117556261A (zh) * 2024-01-08 2024-02-13 浙江大学 一种基于mcnn的隔膜泵单向阀寿命预测方法与系统
CN117629636A (zh) * 2023-12-05 2024-03-01 哈尔滨工程大学 一种燃机滚动轴承的健康评估、故障诊断方法及系统

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100147627A1 (en) * 2008-12-16 2010-06-17 Rockwell Automation Technologies, Inc. Machine conditioning monitoring closed loop lubrication system and method
CN103048137A (zh) * 2012-12-20 2013-04-17 北京航空航天大学 一种变工况下的滚动轴承故障诊断方法
CN105590023A (zh) * 2015-12-08 2016-05-18 三峡大学 一种基于信息熵的滚动轴承性能退化模糊粒化预测方法
CN105973594A (zh) * 2016-04-25 2016-09-28 西北工业大学 一种基于连续深度置信网络的滚动轴承故障预测方法
CN106053066A (zh) * 2016-05-23 2016-10-26 华东交通大学 基于经验模态分解和逻辑回归的滚动轴承性能退化评估方法
CN108280424A (zh) * 2018-01-22 2018-07-13 北京建筑大学 一种基于稀疏编码的滚动轴承剩余寿命预测方法
CN109460618A (zh) * 2018-11-13 2019-03-12 华中科技大学 一种滚动轴承剩余寿命在线预测方法及系统
CN109472241A (zh) * 2018-11-14 2019-03-15 上海交通大学 基于支持向量回归的燃机轴承剩余使用寿命预测方法
CN109726524A (zh) * 2019-03-01 2019-05-07 哈尔滨理工大学 一种基于cnn和lstm的滚动轴承剩余使用寿命预测方法
CN110232249A (zh) * 2019-06-17 2019-09-13 中国人民解放军陆军装甲兵学院 一种滚动轴承剩余寿命预测方法
CN110597240A (zh) * 2019-10-24 2019-12-20 福州大学 一种基于深度学习的水轮发电机组故障诊断方法
CN110633792A (zh) * 2019-10-22 2019-12-31 西安交通大学 端到端基于卷积循环神经网络的轴承健康指标构建方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100147627A1 (en) * 2008-12-16 2010-06-17 Rockwell Automation Technologies, Inc. Machine conditioning monitoring closed loop lubrication system and method
CN103048137A (zh) * 2012-12-20 2013-04-17 北京航空航天大学 一种变工况下的滚动轴承故障诊断方法
CN105590023A (zh) * 2015-12-08 2016-05-18 三峡大学 一种基于信息熵的滚动轴承性能退化模糊粒化预测方法
CN105973594A (zh) * 2016-04-25 2016-09-28 西北工业大学 一种基于连续深度置信网络的滚动轴承故障预测方法
CN106053066A (zh) * 2016-05-23 2016-10-26 华东交通大学 基于经验模态分解和逻辑回归的滚动轴承性能退化评估方法
CN108280424A (zh) * 2018-01-22 2018-07-13 北京建筑大学 一种基于稀疏编码的滚动轴承剩余寿命预测方法
CN109460618A (zh) * 2018-11-13 2019-03-12 华中科技大学 一种滚动轴承剩余寿命在线预测方法及系统
CN109472241A (zh) * 2018-11-14 2019-03-15 上海交通大学 基于支持向量回归的燃机轴承剩余使用寿命预测方法
CN109726524A (zh) * 2019-03-01 2019-05-07 哈尔滨理工大学 一种基于cnn和lstm的滚动轴承剩余使用寿命预测方法
CN110232249A (zh) * 2019-06-17 2019-09-13 中国人民解放军陆军装甲兵学院 一种滚动轴承剩余寿命预测方法
CN110633792A (zh) * 2019-10-22 2019-12-31 西安交通大学 端到端基于卷积循环神经网络的轴承健康指标构建方法
CN110597240A (zh) * 2019-10-24 2019-12-20 福州大学 一种基于深度学习的水轮发电机组故障诊断方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
TIAN Z: "A neural network approach for remaining useful life prediction utilizing both failure and suspension histories", 《MECHANICAL SYSTEMS AND SIGNAL PROCESSING》 *
张焱: "融合失效样本与截尾样本的滚动轴承寿命预测", 《振动与冲击》 *

Cited By (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112541576A (zh) * 2020-12-14 2021-03-23 四川翼飞视科技有限公司 Rgb单目图像的生物活体识别神经网络及其构建方法
CN112541576B (zh) * 2020-12-14 2024-02-20 四川翼飞视科技有限公司 Rgb单目图像的生物活体识别神经网络构建方法
CN112560201A (zh) * 2020-12-24 2021-03-26 新疆大学 复杂载荷工况下风机叶片复合材料可靠性分析方法
CN112580263A (zh) * 2020-12-24 2021-03-30 湖南工业大学 基于时空特征融合的涡扇发动机剩余使用寿命预测方法
CN112560201B (zh) * 2020-12-24 2022-12-06 新疆大学 复杂载荷工况下风机叶片复合材料可靠性分析方法
CN113096378A (zh) * 2021-02-23 2021-07-09 北京航空航天大学 基于深度集合经验模态分解的高速公路od预测方法
CN113052365A (zh) * 2021-02-26 2021-06-29 浙江工业大学 一种基于mswr-lrcn的旋转类机械寿命预测方法
CN112785092B (zh) * 2021-03-09 2024-05-07 中铁电气化局集团有限公司 一种基于自适应深层特征提取的道岔剩余寿命预测方法
CN112785092A (zh) * 2021-03-09 2021-05-11 中铁电气化局集团有限公司 一种基于自适应深层特征提取的道岔剩余寿命预测方法
CN113052060A (zh) * 2021-03-22 2021-06-29 六盘水师范学院 基于数据增强的轴承剩余寿命预测方法、装置及电子设备
CN113052060B (zh) * 2021-03-22 2024-02-09 六盘水师范学院 基于数据增强的轴承剩余寿命预测方法、装置及电子设备
CN113065283A (zh) * 2021-03-26 2021-07-02 深圳技术大学 一种电池寿命预测方法、系统、电子装置及存储介质
CN113051689A (zh) * 2021-04-25 2021-06-29 石家庄铁道大学 基于卷积门控循环网络的轴承剩余使用寿命预测方法
CN113205050A (zh) * 2021-05-09 2021-08-03 西北工业大学 一种基于gru-ae网络的舰船辐射噪声线谱提取方法
CN113449465A (zh) * 2021-06-15 2021-09-28 南京航空航天大学 一种用于滚动轴承的寿命预测方法
CN113486578B (zh) * 2021-06-28 2023-10-20 北京科技大学 一种工业过程中设备剩余寿命的预测方法
CN113486578A (zh) * 2021-06-28 2021-10-08 北京科技大学 一种工业过程中设备剩余寿命的预测方法
CN113807005A (zh) * 2021-08-12 2021-12-17 北京工业大学 基于改进fpa-dbn的轴承剩余寿命预测方法
CN114154743B (zh) * 2021-12-13 2023-04-07 四川大学 基于vetmrrn的空间滚动轴承剩余寿命预测方法
CN114154743A (zh) * 2021-12-13 2022-03-08 四川大学 基于vetmrrn的空间滚动轴承剩余寿命预测方法
CN114298210A (zh) * 2021-12-24 2022-04-08 江苏国科智能电气有限公司 一种多视角融合风电齿轮箱轴承剩余寿命预测方法
CN114564877A (zh) * 2021-12-29 2022-05-31 西安交通大学 一种滚动轴承寿命预测方法、系统、设备及可读存储介质
CN114564877B (zh) * 2021-12-29 2024-05-10 西安交通大学 一种滚动轴承寿命预测方法、系统、设备及可读存储介质
CN114861527A (zh) * 2022-04-15 2022-08-05 南京工业大学 一种基于时间序列特征的锂电池寿命预测方法
CN115446663B (zh) * 2022-10-14 2023-10-20 常州先进制造技术研究所 基于物理指导的深度学习网络的刀具磨损状态监测方法及应用
CN115446663A (zh) * 2022-10-14 2022-12-09 常州先进制造技术研究所 基于物理指导的深度学习网络的刀具磨损状态监测方法及应用
CN115660198A (zh) * 2022-11-07 2023-01-31 中国科学院空间应用工程与技术中心 滚动轴承剩余使用寿命预测方法
CN117629636A (zh) * 2023-12-05 2024-03-01 哈尔滨工程大学 一种燃机滚动轴承的健康评估、故障诊断方法及系统
CN117629636B (zh) * 2023-12-05 2024-05-24 哈尔滨工程大学 一种燃机滚动轴承的健康评估、故障诊断方法及系统
CN117556261A (zh) * 2024-01-08 2024-02-13 浙江大学 一种基于mcnn的隔膜泵单向阀寿命预测方法与系统
CN117556261B (zh) * 2024-01-08 2024-05-14 浙江大学 一种基于mcnn的隔膜泵单向阀寿命预测方法与系统

Also Published As

Publication number Publication date
CN111832216B (zh) 2022-07-12

Similar Documents

Publication Publication Date Title
CN111832216B (zh) 基于eemd-mcnn-gru的滚动轴承剩余使用寿命预测方法
CN110321603B (zh) 一种用于航空发动机气路故障诊断的深度计算模型
He et al. Ensemble transfer CNNs driven by multi-channel signals for fault diagnosis of rotating machinery cross working conditions
Kumar et al. Bearing defect size assessment using wavelet transform based Deep Convolutional Neural Network (DCNN)
Cofre-Martel et al. Deep convolutional neural network-based structural damage localization and quantification using transmissibility data
Elsayed et al. Deep gated recurrent and convolutional network hybrid model for univariate time series classification
Xia et al. Multi-stage fault diagnosis framework for rolling bearing based on OHF Elman AdaBoost-Bagging algorithm
US11593611B2 (en) Neural network cooperation
CN110110809B (zh) 基于机器故障诊断的模糊自动机的构建方法
Ayodeji et al. Causal augmented ConvNet: A temporal memory dilated convolution model for long-sequence time series prediction
Luo et al. Integration of Multi-Gaussian fitting and LSTM neural networks for health monitoring of an automotive suspension component
Ji et al. Model-free fault diagnosis for autonomous underwater vehicles using sequence convolutional neural network
Shang et al. A remaining life prediction of rolling element bearings based on a bidirectional gate recurrent unit and convolution neural network
CN114429150A (zh) 一种基于改进深度子域适应网络的变工况下滚动轴承故障诊断方法及系统
Li et al. Domain adaptation remaining useful life prediction method based on AdaBN-DCNN
CN114266278B (zh) 一种基于双重注意力网络的设备剩余使用寿命预测方法
CN115659174A (zh) 基于图正则化CNN-BiLSTM的多传感器故障诊断方法、介质及设备
Wang et al. Spatiotemporal non-negative projected convolutional network with bidirectional NMF and 3DCNN for remaining useful life estimation of bearings
Wang et al. Remaining useful life prediction of bearings based on convolution attention mechanism and temporal convolution network
Zhao et al. A neural architecture search method based on gradient descent for remaining useful life estimation
CN117421667A (zh) 基于改进灰狼算法优化的Attention-CNN-LSTM工业过程故障诊断方法
Qiu et al. A piecewise method for bearing remaining useful life estimation using temporal convolutional networks
Fu et al. MCA-DTCN: A novel dual-task temporal convolutional network with multi-channel attention for first prediction time detection and remaining useful life prediction
Xu et al. Global attention mechanism based deep learning for remaining useful life prediction of aero-engine
Ayodeji et al. An empirical evaluation of attention-based multi-head models for improved turbofan engine remaining useful life prediction

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