CN114841078A - 一种轴承剩余使用寿命的预测方法及系统 - Google Patents
一种轴承剩余使用寿命的预测方法及系统 Download PDFInfo
- Publication number
- CN114841078A CN114841078A CN202210586316.8A CN202210586316A CN114841078A CN 114841078 A CN114841078 A CN 114841078A CN 202210586316 A CN202210586316 A CN 202210586316A CN 114841078 A CN114841078 A CN 114841078A
- Authority
- CN
- China
- Prior art keywords
- bearing
- unit
- service life
- health
- predicting
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000036541 health Effects 0.000 claims abstract description 81
- 230000015654 memory Effects 0.000 claims abstract description 59
- 238000012549 training Methods 0.000 claims abstract description 38
- 239000000203 mixture Substances 0.000 claims abstract description 26
- 230000007246 mechanism Effects 0.000 claims abstract description 14
- 238000009826 distribution Methods 0.000 claims description 35
- 239000013598 vector Substances 0.000 claims description 35
- 239000011159 matrix material Substances 0.000 claims description 30
- 238000010276 construction Methods 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 9
- 230000002787 reinforcement Effects 0.000 claims 1
- 238000006731 degradation reaction Methods 0.000 description 13
- 230000015556 catabolic process Effects 0.000 description 12
- 210000004027 cell Anatomy 0.000 description 11
- 238000013528 artificial neural network Methods 0.000 description 7
- 238000012360 testing method Methods 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000001965 increasing effect Effects 0.000 description 3
- 230000004913 activation Effects 0.000 description 2
- 239000000306 component Substances 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000007787 long-term memory Effects 0.000 description 2
- 238000005096 rolling process Methods 0.000 description 2
- 208000000044 Amnesia Diseases 0.000 description 1
- 208000026139 Memory disease Diseases 0.000 description 1
- 238000005299 abrasion Methods 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 210000003850 cellular structure Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000006984 memory degeneration Effects 0.000 description 1
- 208000023060 memory loss Diseases 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 230000000306 recurrent effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
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
- G06N3/044—Recurrent networks, e.g. Hopfield networks
-
- 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
- G06N3/048—Activation functions
-
- 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
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
-
- 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)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Mathematical Physics (AREA)
- Computational Linguistics (AREA)
- Molecular Biology (AREA)
- Computing Systems (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)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明涉及一种轴承剩余使用寿命的预测方法及系统，方法包括：采集轴承生命周期内的原始振动信号；利用高斯混合模型、KL散度和无监督健康指标方法，根据所述原始振动信号生成健康指标集；基于注意力机制构建强化记忆门控循环单元；所述强化记忆门控循环单元包括当前时刻之前的两个时刻的状态信息；根据所述健康指标集对所述强化记忆门控循环单元进行训练和预测，得到预测结果；根据预测结果计算轴承剩余使用寿命。本发明提出的强化记忆门控循环单元能够对构建的历史趋势信息进行充分学习，提高了网络对剩余使用寿命预测的精度和鲁棒性。
Description
技术领域
本发明涉及轴承寿命预测技术领域，特别是涉及一种轴承剩余使用寿命的预测方法及系统。
背景技术
设备的安全可靠运行是保证工厂生产活动的核心内容。为了确保设备的运行状态和制定合理的维修计划，研究具有高精度和高鲁棒性的寿命预测方法非常重要。滚动轴承作为核心零部件之一，被广泛的应用于飞机、高铁等。然而轴承常常处于恶劣的工作环境下，这容易使轴承发生磨损、点蚀等失效故障，进而影响设备的正常运行。因此有必要对滚动轴承的剩余使用寿命进行高精度的预测。
发明内容
为了克服现有技术的不足，本发明的目的是提供一种轴承剩余使用寿命的预测方法及系统。
为实现上述目的，本发明提供了如下方案：
一种轴承剩余使用寿命的预测方法，包括：
采集轴承生命周期内的原始振动信号；
利用高斯混合模型、KL散度和无监督健康指标方法，根据所述原始振动信号生成健康指标集；
基于注意力机制构建强化记忆门控循环单元；所述强化记忆门控循环单元包括当前时刻之前的两个时刻的状态信息；
根据所述健康指标集对所述强化记忆门控循环单元进行训练和预测，得到预测结果；
根据预测结果计算轴承剩余使用寿命。
优选地，所述利用高斯混合模型、KL散度和无监督健康指标方法，根据所述原始振动信号生成健康指标集，包括：
将所述原始振动信号的数据集中的前H个样本作为健康样本，通过所述高斯混合模型估计所述健康样本的基准分布；H为设定健康样本阈值；
根据所述高斯混合模型估计所述原始振动信号的数据集中的剩余样本的概率分布；
利用所述KL散度对所述基准分布和所述概率分布进行计算，得到KLD值；所述KLD值用于表征轴承的健康状态；
根据所述KLD值构建所述健康指标集。
优选地，所述根据所述健康指标集对所述强化记忆门控循环单元进行训练和预测，得到预测结果，包括：
对所述健康指标集进行划分，得到训练集；
获取预设的细胞单元数k，并根据所述细胞单元数和所述训练集构建输入矩阵；
将所述输入矩阵的前k个行向量作为所述强化记忆门控循环单元的输入，将所述输入矩阵的最后一个行向量作为所述强化记忆门控循环单元的输出，并根据损失函数对所述强化记忆门控循环单元进行训练；
将所述输入矩阵的最后k个行向量输入至训练好的强化记忆门控循环单元，得到所述预测结果；所述预测结果包括至少一个健康指标向量。
优选地，所述根据预测结果计算轴承剩余使用寿命，包括：
当所述健康指标向量超过预设失效阈值时，通过预设采样时间和记录间隔计算轴承剩余使用寿命。
一种轴承剩余使用寿命的预测系统，包括：
采集模块，用于采集轴承生命周期内的原始振动信号；
指标生成模块，用于利用高斯混合模型、KL散度和无监督健康指标方法，根据所述原始振动信号生成健康指标集；
模型构建模块，用于基于注意力机制构建强化记忆门控循环单元；所述强化记忆门控循环单元包括当前时刻之前的两个时刻的状态信息；
预测模块，用于根据所述健康指标集对所述强化记忆门控循环单元进行训练和预测，得到预测结果；
计算模块，用于根据预测结果计算轴承剩余使用寿命。
优选地，所述指标生成模块具体包括：
基准单元，用于将所述原始振动信号的数据集中的前H个样本作为健康样本，通过所述高斯混合模型估计所述健康样本的基准分布；H为设定健康样本阈值；
概率单元，用于根据所述高斯混合模型估计所述原始振动信号的数据集中的剩余样本的概率分布；
状态计算单元，用于利用所述KL散度对所述基准分布和所述概率分布进行计算，得到KLD值；所述KLD值用于表征轴承的健康状态；
指标集构建单元，用于根据所述KLD值构建所述健康指标集。
优选地，所述预测模块具体包括：
划分单元，用于对所述健康指标集进行划分，得到训练集；
矩阵构建单元，用于获取预设的细胞单元数k，并根据所述细胞单元数和所述训练集构建输入矩阵；
训练单元，用于将所述输入矩阵的前k个行向量作为所述强化记忆门控循环单元的输入，将所述输入矩阵的最后一个行向量作为所述强化记忆门控循环单元的输出，并根据损失函数对所述强化记忆门控循环单元进行训练；
预测单元，用于将所述输入矩阵的最后k个行向量输入至训练好的强化记忆门控循环单元，得到所述预测结果；所述预测结果包括至少一个健康指标向量。
优选地，所述计算模块具体包括：
寿命计算单元，用于当所述健康指标向量超过预设失效阈值时，通过预设采样时间和记录间隔计算轴承剩余使用寿命。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明提供了一种轴承剩余使用寿命的预测方法及系统，方法包括：采集轴承生命周期内的原始振动信号；利用高斯混合模型、KL散度和无监督健康指标方法，根据所述原始振动信号生成健康指标集；基于注意力机制构建强化记忆门控循环单元；所述强化记忆门控循环单元包括当前时刻之前的两个时刻的状态信息；根据所述健康指标集对所述强化记忆门控循环单元进行训练和预测，得到预测结果；根据预测结果计算轴承剩余使用寿命。本发明结合艾宾浩斯遗忘曲线和门控循环单元神经网络提出了强化记忆门控循环单元能够对构建的历史趋势信息进行充分学习，提高了网络对剩余使用寿命预测的精度和鲁棒性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明提供的实施例中的预测方法的流程图；
图2为本发明提供的实施例中的GRU网络结构示意图；
图3为本发明提供的实施例中的寿命预测流程图；
图4为本发明提供的实施例中的RMGRU网络结构图示意图；
图5为本发明提供的实施例中的PRONOSTIA试验台；
图6为本发明提供的实施例中的第一轴承的KLD-based HIs示意图；
图7为本发明提供的实施例中的第二轴承的KLD-based HIs示意图；
图8为本发明提供的实施例中的第一轴承的HI曲线；
图9为本发明提供的实施例中的第二轴承的HI曲线；
图10为本发明提供的实施例中的不同预测模型的评价指标示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
在本文中提及“实施例”意味着，结合实施例描述的特定特征、结构或特性可以包含在本申请的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例，也不是与其它实施例互斥的独立的或备选的实施例。本领域技术人员显式地和隐式地理解的是，本文所描述的实施例可以与其它实施例相结合。
本申请的说明书和权利要求书及所述附图中的术语“第一”、“第二”、“第三”和“第四”等是用于区别不同对象，而不是用于描述特定顺序。此外，术语“包括”和“具有”以及它们任何变形，意图在于覆盖不排他的包含。例如包含了一系列步骤、过程、方法等没有限定于已列出的步骤，而是可选地还包括没有列出的步骤，或可选地还包括对于这些过程、方法、产品或设备固有的其它步骤元。
本发明的目的是提供一种轴承剩余使用寿命的预测方法及系统，能够提高对剩余使用寿命预测的精度和鲁棒性。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明提供的实施例中的预测方法的流程图，如图1所示，本发明提供了一种轴承剩余使用寿命的预测方法，包括：
步骤100：采集轴承生命周期内的原始振动信号；
步骤200：利用高斯混合模型、KL散度和无监督健康指标方法，根据所述原始振动信号生成健康指标集；
步骤300：基于注意力机制构建强化记忆门控循环单元；所述强化记忆门控循环单元包括当前时刻之前的两个时刻的状态信息；
步骤400：根据所述健康指标集对所述强化记忆门控循环单元进行训练和预测，得到预测结果；
步骤500：根据预测结果计算轴承剩余使用寿命。
本实施例中为了减少人为干扰和专家知识需求，利用高斯混合模型(Gaussianmixture model，GMM)和Kullback-Leibler散度(Kullback-Leibler divergence，KLD)构造了无监督健康指标(health indicator，HI)，命名为KLD-based HI。针对当前的寿命预测方法，主要是通过改进循环神经网络的细胞结构，并未对历史状态信息进行充分的利用。因此本实施例结合艾宾浩斯遗忘曲线和门控循环单元神经网络(Gated RecurrentUnit，GRU)提出了强化记忆门控循环单元神经网络(reinforced memory GRU，RMGRU)。提出的RMGRU能够对构建的KLD-based HI的历史趋势信息进行充分学习，以此提高网络对剩余使用寿命预测(remaininguseful life，RUL)的精度和鲁棒性。
其中，门控循环单元(GRU)作为长短时记忆网络(LSTM)的一种变体，GRU的性能与LSTM相似，但它的参数更少，计算成本更低。GRU将LSTM中的输入门和遗忘门合并成一个单独的更新门。同时，它还混合了LSTM的细胞状态和隐藏状态。因此，一个GRU单元由两个门组成，即复位门和更新门。更新门决定有多少历史信息被引入到当前状态。复位门根据历史状态控制当前输入的信息。复位门的值越小，添加到当前状态的信息就越少。GRU网络结构如图2所示，其中σ和tanh分别表示s形激活函数和双曲正切激活函数。对于单个GRU单元，t时刻输出隐藏状态计算为：
zt＝σ(Wz[xt,ht-1]+bz)
rt＝σ(Wr[xt,ht-1]+br)
其中xt表示在t时刻的输入；ht-1和ht分别表示前一时刻和当前时刻的隐藏状态；
表示暂时的隐藏状态；Wz,Wr,Wh表示权重矩阵；bz，br，bh表示偏置矩阵；重置门的输出为rt，更新门的输出为zt，⊙表示点积。
图3为本发明提供的实施例中的寿命预测流程图，如图3所示，本实施例包括KLD-based HI构建和RUL预测。首先，利用高斯混合模型和Kullback-Leibler散度从原始振动信号中提取具有明显退化趋势的无监督健康指标；然后设计了RMGRU，并将其应用于KLD-based HI值的逐步预测。当预测的HI值超过预定的失效阈值，立即计算轴承的RUL。
优选地，所述步骤200具体包括：
将所述原始振动信号的数据集中的前H个样本作为健康样本，通过所述高斯混合模型估计所述健康样本的基准分布；H为设定健康样本阈值；
根据所述高斯混合模型估计所述原始振动信号的数据集中的剩余样本的概率分布；
利用所述KL散度对所述基准分布和所述概率分布进行计算，得到KLD值；所述KLD值用于表征轴承的健康状态；
根据所述KLD值构建所述健康指标集。
本实施例中KLD-based HI构建过程如下：
以所有健康样本的概率分布p0为基准。当轴承健康时，所采集样本的概率分布几乎类似于p0。一旦轴承发生退化，所采集样本的概率分布与p0不同，这种差异一般会随着退化程度的增加而增大。因此，p0与监测信号的概率分布之间的这个差值可以用来表示轴承的健康状况，可以通过KLD进行计算。基于两个单变量高斯分布的KLD，本实施例提出了KLD-based HI，其构造过程总结如下：
1)采集轴承全生命周期振动数据集。X＝[x1，x2，…，xK]表示获取的数据集，其中xi＝(x1，x2，…，xk)，K和k表示样本的数量和样本的维度。
2)数据集X前H个样本x1，x2，…，xH作为健康样本。通过GMM估计健康样本的基准分布p0。其中GMM的高斯分量设置为4。
3)采用高斯分量为4的GMM模型估计剩余样本xH+1，…，xK的概率分布p＝[pH+1，…，pK]。
4)构建KLD-based HI。使用KLD，分别计算基准分布p0与样本的概率分布p＝[pH+1，…，pK]之间的KLD值。该计算得到的KLD值能够用于表征轴承的健康状态。
优选地，所述步骤400具体包括：
对所述健康指标集进行划分，得到训练集；
获取预设的细胞单元数k，并根据所述细胞单元数和所述训练集构建输入矩阵；
将所述输入矩阵的前k个行向量作为所述强化记忆门控循环单元的输入，将所述输入矩阵的最后一个行向量作为所述强化记忆门控循环单元的输出，并根据损失函数对所述强化记忆门控循环单元进行训练；
将所述输入矩阵的最后k个行向量输入至训练好的强化记忆门控循环单元，得到所述预测结果；所述预测结果包括至少一个健康指标向量。
优选地，所述步骤500具体包括：
当所述健康指标向量超过预设失效阈值时，通过预设采样时间和记录间隔计算轴承剩余使用寿命。
本实施例中在预测轴承RUL时，RNN除了使用当前状态信息外，还将使用历史趋势信息。这与人类记忆的机制相似。因此，RNN的历史记忆能力在处理长序列时非常重要。现有的基于RNN的模型只是改进了RNN的细胞结构，而对历史趋势信息的利用不足。艾宾浩斯(H.Ebbinghaus)在1885年揭示的遗忘曲线表明，记忆在学习后立即开始丧失，在记住后的短时间内遗忘率最大。在遗忘规律的启发下，构建RMGRU网络，增强趋势特征的学习能力。RMGRU网络结构如图4所示。在该图中，x和h分别表示输入数据和隐藏状态，S表示注意力机制，st是使用注意力机制后的输出。
为了提高长时记忆能力，获得更鲁棒的预测结果，包含所学知识的隐藏状态需要进一步记忆。根据遗忘规律，对所学知识的及时巩固可以减缓记忆的丢失，因此通过结合前两个时刻的状态信息，可以增强RNN的长期记忆能力。由于前两次直接使用历史状态信息会导致信息冗余，增加网络的计算负担，因此采用注意力机制强调健康指标的趋势特征信息，为其分配计算资源。如图4所示，加入ht-2和ht-1，然后利用注意机制从前两个时刻的状态信息中获取重要的趋势特征。由图3可知，t时刻的输出隐藏状态被更新为：
s(ht-1,ht-2)＝VTtanh(Wsht-1+Usht-2)
zt＝σ(Az[xt,At-1]+dz)
rt＝σ(Ar[xt,At-1]+dr)
由上可知，隐藏状态比其他经典RNN网络包含更多的趋势信息。因此，RMGRU网络具有较强的预测能力，尤其是在少样本的情况下。
本实施例利用所提出KLD-based HI和RMGRU神经网络，可以对轴承的剩余使用寿命进行预测。本实施例提出了一种新的时间序列预测方法，具体步骤如下：
1)采集轴承生命周期内的原始振动信号。采样时间为T，记录间隔为Ts。在测试轴承的生命周期中，假设获得了n个样本。
2)通过高斯混合模型对原始振动数据的分布参数进行估计，得到振动信号的分布参数即均值和方差，然后通过KLD-based HI构建方法生成健康指标，并表示为H＝[h1,h2,…,hn]T。
3)划分训练集和测试集。对于生成的健康指标H，前m个健康指标作为训练集Strain，剩余的健康指标作为测试集Stest，即有：Strain＝[h1,h2,…,hm]T和Stest＝[hm+1,…,hn]。
4)设置输入的细胞单元数k，通过训练集Strain和k构建输入矩阵U：
其中，
Ui＝[hi hi+1 … hm-k+i+1]
5)将U的前k个行向量作为RMGRU神经网络的输入，将U的最后一个行向量设置为输出。定义RMGRU的映射函数为f，通过最小化损失函数L更新RMGRU的参数W和b:
min[L(W,b)]＝f(U1,U2,…,Uk)-Uk+1 2
6)训练RMGRU后，将最后的k个向量输入训练后的神经网络，预测下一个健康指标向量，其逐步预测的步骤如下：
7)当预测的健康指标hm-k+i+1(hm-k+i+1∈Ui)的值超过预先设定的失效阈值时，可以通过(i-m)×(T+Ts)计算出轴承剩余使用寿命。此外，失效点由实验结果所得，失效点的健康指标值即为失效阈值。
本实施例采用实验方法对上述方法进行验证，实验数据来自PRONOSTIA实验平台，结构如图5所示。本实验是专门为研究轴承故障的检测、诊断和预测方法而设计的。该实验系统可以在几小时内对轴承进行加速退化试验。实验中加速度传感器固定在轴承的外圈上，捕捉振动信号。在本次实验验证中，对轴承1_1进行剩余使用寿命预测。该轴承运行工况为内圈转速为1800转/分，载荷为4kN。采样频率为25.6kHz。采样每10秒重复一次，每次0.1s。
采用以上基于KLD-based HI构建方法，选取轴承1_1(第一轴承)和轴承1_3(第二轴承)的前200个样本作为健康样本，利用它们计算基线分布p0。接下来，构建轴承1_1和轴承1_3的KLD-based HI向量，分别如图6和图7所示。从图6和图7可以看出，生成的HI能够很好地反映两个轴承的健康阶段和退化趋势。当轴承处于健康阶段时KLD-based HI值非常接近于0，并略有波动。第一个退化点由三次标准差原则确定，即当某一点的KLD-based HI值超过健康信号的三次标准差时，将该点定义为第一个退化点。然后随着轴承的退化，KLD-based HI值和波动均逐渐增大。最后HI值急剧上升，表明轴承开始失效。KLD-based HI值都先保持稳定，然后逐渐增加，最后在轴承失效时急剧上升。因此，构建的KLD-based HI能够反映轴承的真实健康状况。通过一系列实验，发现当HI值达到3时轴承发生初始失效，因此将其设为失效阈值。
通过构建的KLD-based HI和所提出的RMGRU网络，对轴承RUL进行预测。所提出的RUL预测方法更关注第一个退化点到初始故障点之间的退化过程，因为它包含了大部分的退化信息。然而，为了充分了解退化信息，还需要将健康阶段的一些HI要素输入预测模型。轴承1_1选取从1001到2550的1550个元素，用前1450个元素预测后100个元素的变化趋势。同样，对于轴承1_3，从第901元素到第2070元素共选取1170个元素，用前1070个元素预测后100个元素。利用所提出的RUL预测方法，得到了轴承1_1和轴承1_3的预测结果。图8绘制轴承1_1的训练、预测和实际HI曲线，图9绘制轴承1_3的HI曲线。从图8和图9可以看出，训练和预测规则在几乎所有时刻都非常接近真实的HI曲线，训练和预测HI曲线的退化趋势也接近于实际的HI曲线。在图8中，HI元素不超过预设故障阈值的个数为93，在图9中为95。由此可知轴承1_1的RUL值为930s，轴承1_3的RUL值为950s。该实验结果表明了提出的RMGRU能够较好的预测轴承的剩余使用寿命。为了进一步证明提出的RMGRU的优越性，本实施例对比了长短时记忆网络(long-short termmemory，LSTM)、门控循环单元(gatedrecurrentunit，GRU)、门控双注意力单元(gated dual attention unit，GDAU)、宏微观注意力LSTM网络(macroscopic-microscopic attention LSTM，MMALSTM)以及JANET网络，并通过MAE、NRMSE和Score函数进行定量评价，其结果如图10所示。
可见，RMGRU的MAE和NRMSE最小，RMGRU的Score最高，因此RMGRU网络比经典预测网络具有更好的预测性能。原因在于，提出的RMGRU能够及时重新学习遗忘的趋势特征，即RMGRU能够比当前基于RNN的网络记住更多有用的退化趋势信息。对比实验结果也表明，对利用历史状态信息的再次使用比单纯改善细胞结构更有效。因此，提出的基于KLD-basedHI和RMGRU的寿命预测方法能够很好地预测轴承的剩余使用寿命。
对应上述方法，本实施例还提供了一种轴承剩余使用寿命的预测系统，包括：
采集模块，用于采集轴承生命周期内的原始振动信号；
指标生成模块，用于利用高斯混合模型、KL散度和无监督健康指标方法，根据所述原始振动信号生成健康指标集；
模型构建模块，用于基于注意力机制构建强化记忆门控循环单元；所述强化记忆门控循环单元包括当前时刻之前的两个时刻的状态信息；
预测模块，用于根据所述健康指标集对所述强化记忆门控循环单元进行训练和预测，得到预测结果；
计算模块，用于根据预测结果计算轴承剩余使用寿命。
优选地，所述指标生成模块具体包括：
基准单元，用于将所述原始振动信号的数据集中的前H个样本作为健康样本，通过所述高斯混合模型估计所述健康样本的基准分布；H为设定健康样本阈值；
概率单元，用于根据所述高斯混合模型估计所述原始振动信号的数据集中的剩余样本的概率分布；
状态计算单元，用于利用所述KL散度对所述基准分布和所述概率分布进行计算，得到KLD值；所述KLD值用于表征轴承的健康状态；
指标集构建单元，用于根据所述KLD值构建所述健康指标集。
优选地，所述预测模块具体包括：
划分单元，用于对所述健康指标集进行划分，得到训练集；
矩阵构建单元，用于获取预设的细胞单元数k，并根据所述细胞单元数和所述训练集构建输入矩阵；
训练单元，用于将所述输入矩阵的前k个行向量作为所述强化记忆门控循环单元的输入，将所述输入矩阵的最后一个行向量作为所述强化记忆门控循环单元的输出，并根据损失函数对所述强化记忆门控循环单元进行训练；
预测单元，用于将所述输入矩阵的最后k个行向量输入至训练好的强化记忆门控循环单元，得到所述预测结果；所述预测结果包括至少一个健康指标向量。
优选地，所述计算模块具体包括：
寿命计算单元，用于当所述健康指标向量超过预设失效阈值时，通过预设采样时间和记录间隔计算轴承剩余使用寿命。
本发明的有益效果如下：
本发明结合艾宾浩斯遗忘曲线和门控循环单元神经网络提出了强化记忆门控循环单元能够对构建的历史趋势信息进行充分学习，提高了网络对剩余使用寿命预测的精度和鲁棒性。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统小而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种轴承剩余使用寿命的预测方法，其特征在于，包括：
采集轴承生命周期内的原始振动信号；
利用高斯混合模型、KL散度和无监督健康指标方法，根据所述原始振动信号生成健康指标集；
基于注意力机制构建强化记忆门控循环单元；所述强化记忆门控循环单元包括当前时刻之前的两个时刻的状态信息；
根据所述健康指标集对所述强化记忆门控循环单元进行训练和预测，得到预测结果；
根据预测结果计算轴承剩余使用寿命。
2.根据权利要求1所述的轴承剩余使用寿命的预测方法，其特征在于，所述利用高斯混合模型、KL散度和无监督健康指标方法，根据所述原始振动信号生成健康指标集，包括：
将所述原始振动信号的数据集中的前H个样本作为健康样本，通过所述高斯混合模型估计所述健康样本的基准分布；H为设定健康样本阈值；
根据所述高斯混合模型估计所述原始振动信号的数据集中的剩余样本的概率分布；
利用所述KL散度对所述基准分布和所述概率分布进行计算，得到KLD值；所述KLD值用于表征轴承的健康状态；
根据所述KLD值构建所述健康指标集。
3.根据权利要求1所述的轴承剩余使用寿命的预测方法，其特征在于，所述根据所述健康指标集对所述强化记忆门控循环单元进行训练和预测，得到预测结果，包括：
对所述健康指标集进行划分，得到训练集；
获取预设的细胞单元数k，并根据所述细胞单元数和所述训练集构建输入矩阵；
将所述输入矩阵的前k个行向量作为所述强化记忆门控循环单元的输入，将所述输入矩阵的最后一个行向量作为所述强化记忆门控循环单元的输出，并根据损失函数对所述强化记忆门控循环单元进行训练；
将所述输入矩阵的最后k个行向量输入至训练好的强化记忆门控循环单元，得到所述预测结果；所述预测结果包括至少一个健康指标向量。
4.根据权利要求3所述的轴承剩余使用寿命的预测方法，其特征在于，所述根据预测结果计算轴承剩余使用寿命，包括：
当所述健康指标向量超过预设失效阈值时，通过预设采样时间和记录间隔计算轴承剩余使用寿命。
5.一种轴承剩余使用寿命的预测系统，其特征在于，包括：
采集模块，用于采集轴承生命周期内的原始振动信号；
指标生成模块，用于利用高斯混合模型、KL散度和无监督健康指标方法，根据所述原始振动信号生成健康指标集；
模型构建模块，用于基于注意力机制构建强化记忆门控循环单元；所述强化记忆门控循环单元包括当前时刻之前的两个时刻的状态信息；
预测模块，用于根据所述健康指标集对所述强化记忆门控循环单元进行训练和预测，得到预测结果；
计算模块，用于根据预测结果计算轴承剩余使用寿命。
6.根据权利要求5所述的轴承剩余使用寿命的预测系统，其特征在于，所述指标生成模块具体包括：
基准单元，用于将所述原始振动信号的数据集中的前H个样本作为健康样本，通过所述高斯混合模型估计所述健康样本的基准分布；H为设定健康样本阈值；
概率单元，用于根据所述高斯混合模型估计所述原始振动信号的数据集中的剩余样本的概率分布；
状态计算单元，用于利用所述KL散度对所述基准分布和所述概率分布进行计算，得到KLD值；所述KLD值用于表征轴承的健康状态；
指标集构建单元，用于根据所述KLD值构建所述健康指标集。
7.根据权利要求5所述的轴承剩余使用寿命的预测系统，其特征在于，所述预测模块具体包括：
划分单元，用于对所述健康指标集进行划分，得到训练集；
矩阵构建单元，用于获取预设的细胞单元数k，并根据所述细胞单元数和所述训练集构建输入矩阵；
训练单元，用于将所述输入矩阵的前k个行向量作为所述强化记忆门控循环单元的输入，将所述输入矩阵的最后一个行向量作为所述强化记忆门控循环单元的输出，并根据损失函数对所述强化记忆门控循环单元进行训练；
预测单元，用于将所述输入矩阵的最后k个行向量输入至训练好的强化记忆门控循环单元，得到所述预测结果；所述预测结果包括至少一个健康指标向量。
8.根据权利要求7所述的轴承剩余使用寿命的预测系统，其特征在于，所述计算模块具体包括：
寿命计算单元，用于当所述健康指标向量超过预设失效阈值时，通过预设采样时间和记录间隔计算轴承剩余使用寿命。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210586316.8A CN114841078B (zh) | 2022-05-26 | 2022-05-26 | 一种轴承剩余使用寿命的预测方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210586316.8A CN114841078B (zh) | 2022-05-26 | 2022-05-26 | 一种轴承剩余使用寿命的预测方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114841078A true CN114841078A (zh) | 2022-08-02 |
| CN114841078B CN114841078B (zh) | 2024-09-06 |
Family
ID=82571991
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210586316.8A Active CN114841078B (zh) | 2022-05-26 | 2022-05-26 | 一种轴承剩余使用寿命的预测方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114841078B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115329283A (zh) * | 2022-10-12 | 2022-11-11 | 南通翔润机电有限公司 | 一种启动马达高强度换向器寿命预测方法 |
| CN116124460A (zh) * | 2022-12-26 | 2023-05-16 | 江西理工大学 | 一种基于健康指标构建的轴承寿命预测方法及系统 |
| CN117664576A (zh) * | 2023-11-29 | 2024-03-08 | 东北大学 | 一种基于改进优化器深度学习模型的轴承健康评估方法 |
| CN118503928A (zh) * | 2024-07-17 | 2024-08-16 | 广东石油化工学院 | 一种轴承剩余寿命预测方法、装置、设备及存储介质 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4128147A (en) * | 1976-11-30 | 1978-12-05 | Antonio Lafuente Ruberte | Vehicle speed reducer or eddy current brake |
| CN110866314A (zh) * | 2019-10-22 | 2020-03-06 | 东南大学 | 多层双向门控循环单元网络的旋转机械剩余寿命预测方法 |
| CN113051689A (zh) * | 2021-04-25 | 2021-06-29 | 石家庄铁道大学 | 基于卷积门控循环网络的轴承剩余使用寿命预测方法 |
-
2022
- 2022-05-26 CN CN202210586316.8A patent/CN114841078B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4128147A (en) * | 1976-11-30 | 1978-12-05 | Antonio Lafuente Ruberte | Vehicle speed reducer or eddy current brake |
| CN110866314A (zh) * | 2019-10-22 | 2020-03-06 | 东南大学 | 多层双向门控循环单元网络的旋转机械剩余寿命预测方法 |
| CN113051689A (zh) * | 2021-04-25 | 2021-06-29 | 石家庄铁道大学 | 基于卷积门控循环网络的轴承剩余使用寿命预测方法 |
Non-Patent Citations (2)
| Title |
|---|
| JIANGHONG ZHOU ET AL.: "Remaining useful life prediction of bearings by a new reinforced memory GRU network", ADVANCED ENGINEERING INFORMATICS, vol. 53, 1 July 2022 (2022-07-01), pages 1 - 11 * |
| 秦毅 等: "基于流形学习和最小二乘支持向量机的滚动轴承退化趋势预测", 振动与冲击, vol. 34, no. 09, 15 May 2015 (2015-05-15), pages 149 - 153 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115329283A (zh) * | 2022-10-12 | 2022-11-11 | 南通翔润机电有限公司 | 一种启动马达高强度换向器寿命预测方法 |
| CN115329283B (zh) * | 2022-10-12 | 2023-03-28 | 南通翔润机电有限公司 | 一种启动马达高强度换向器寿命预测方法 |
| CN116124460A (zh) * | 2022-12-26 | 2023-05-16 | 江西理工大学 | 一种基于健康指标构建的轴承寿命预测方法及系统 |
| CN116124460B (zh) * | 2022-12-26 | 2024-01-30 | 江西理工大学 | 一种基于健康指标构建的轴承寿命预测方法及系统 |
| CN117664576A (zh) * | 2023-11-29 | 2024-03-08 | 东北大学 | 一种基于改进优化器深度学习模型的轴承健康评估方法 |
| CN118503928A (zh) * | 2024-07-17 | 2024-08-16 | 广东石油化工学院 | 一种轴承剩余寿命预测方法、装置、设备及存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114841078B (zh) | 2024-09-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114841078A (zh) | 一种轴承剩余使用寿命的预测方法及系统 | |
| CN109472110A (zh) | 一种基于lstm网络和arima模型的航空发动机剩余使用寿命预测方法 | |
| CN107144430B (zh) | 一种基于增量学习的轴承故障诊断方法 | |
| CN104077445B (zh) | 基于模糊理论的加速寿命试验统计分析方法 | |
| CN113420509A (zh) | 一种风电机组状态评价方法、装置及存储介质 | |
| US20220122103A1 (en) | Customized product performance prediction method based on heterogeneous data difference compensation fusion | |
| CN111611744B (zh) | 基于循环卷积网络及变分推理的滚动轴承寿命预测方法 | |
| Khawaja et al. | Reasoning about uncertainty in prognosis: a confidence prediction neural network approach | |
| CN111783242B (zh) | 一种基于rvm-kf的滚动轴承剩余寿命预测方法及装置 | |
| CN114266201B (zh) | 一种基于深度学习的自注意力电梯困人预测方法 | |
| CN115422687A (zh) | 一种滚动轴承的寿命预测方法 | |
| CN116992638A (zh) | 基于相似性下失效阈值确定的滚动轴承rul预测方法 | |
| CN116070508A (zh) | 一种桥梁健康状态预测预警方法 | |
| CN113673774A (zh) | 基于自编码器和时序卷积网络的航空发动机余寿预测方法 | |
| CN117454324A (zh) | 风速风向仪故障预警方法、系统、电子设备及存储介质 | |
| CN112560252A (zh) | 一种航空发动机剩余寿命预测方法 | |
| CN116822722A (zh) | 水位预测方法、系统、装置、电子设备及介质 | |
| CN116579233A (zh) | 一种机械设备的剩余寿命预测方法 | |
| CN115898927A (zh) | 一种基于机器学习的聚变装置分子泵故障诊断与预测方法 | |
| CN115296298A (zh) | 一种风电场功率预测方法 | |
| CN115459982A (zh) | 一种电力网络虚假数据注入攻击检测方法 | |
| CN114995164A (zh) | 一种基于物联网的新能源汽车安全预警方法及装置 | |
| CN109886386B (zh) | 唤醒模型的确定方法及装置 | |
| CN114595525A (zh) | 一种基于线性回归和tcn的滚动轴承退化趋势预测方法 | |
| CN113962020A (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 |