CN115130514A - Construction method and system for health index of engineering equipment - Google Patents
Construction method and system for health index of engineering equipment Download PDFInfo
- Publication number
- CN115130514A CN115130514A CN202210767222.0A CN202210767222A CN115130514A CN 115130514 A CN115130514 A CN 115130514A CN 202210767222 A CN202210767222 A CN 202210767222A CN 115130514 A CN115130514 A CN 115130514A
- Authority
- CN
- China
- Prior art keywords
- continuous
- network
- frequency domain
- time
- engineering equipment
- 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
Images
Classifications
-
- 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
- G06N3/088—Non-supervised learning, e.g. competitive learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Development Economics (AREA)
- General Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- General Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明涉及一种工程设备健康指标构建方法及系统,属于可靠性工程领域。方法包括:获取工程设备测试数据并计算其时频域信息;将时频域信息输入至连续深度置信网络,借助于对比散度算法对其进行训练;采用训练好的连续深度置信网络提取深层次故障特征;将深层次故障特征作为SOM神经网络的输入进行无监督学习;通过计算训练好的SOM神经网络的输入层的状态向量与全部单元的权重向量间的欧氏距离,构建出工程设备的健康指标。本发明方法及系统采用连续深度置信网络提取蕴含在时频域背后的深层次故障特征,利用SOM神经网络在确保原始特征拓扑结构前提下构建设备的健康指标,提高了健康指标构建的准确性。
The invention relates to a method and a system for constructing a health index of engineering equipment, and belongs to the field of reliability engineering. The method includes: acquiring test data of engineering equipment and calculating its time-frequency domain information; inputting the time-frequency domain information into a continuous depth belief network, and training it by means of a contrastive divergence algorithm; using the trained continuous depth belief network to extract the deep level Fault features; take the deep fault features as the input of the SOM neural network for unsupervised learning; by calculating the Euclidean distance between the state vector of the input layer of the trained SOM neural network and the weight vector of all units, the engineering equipment is constructed. health indicators. The method and system of the present invention use a continuous deep confidence network to extract the deep-level fault features contained in the time-frequency domain, and use the SOM neural network to construct the health index of the equipment on the premise of ensuring the original feature topology structure, thereby improving the accuracy of the construction of the health index.
Description
技术领域technical field
本发明涉及可靠性工程技术领域,特别是涉及一种工程设备健康指标构建方法及系统。The invention relates to the technical field of reliability engineering, in particular to a method and system for constructing health indicators of engineering equipment.
背景技术Background technique
随着传感监测技术的发展与电子器件工艺的改进,可在设备受到较小影响的情形下收集到越来越多的测试数据,再加上设备可靠性的日益提升,设备的测试数据数量将会激增。特别地,在当今物联网与工业4.0时代,工程设备获取的测试数据会呈现海量化、非线性与动态性等特征。当采用传统方法确定设备的健康指标时,难以提取出设备的深层次故障特征,由此导致健康指标的准确性会受到一定程度影响。为了充分挖掘海量数据蕴含的隐藏特征,Hinton提出了一种分析海量数据的有效工具——深度学习,此后一直为广大学者与工程技术人员所推崇。深度学习已成功应用于许多工程领域,如图像处理、语音识别、机器翻译等,同时基于该技术的故障特征提取与健康指标构建受到了重点关注,也取得了一系列成果。With the development of sensor monitoring technology and the improvement of electronic device technology, more and more test data can be collected under the condition that the equipment is less affected. will proliferate. In particular, in today's Internet of Things and Industry 4.0 era, the test data acquired by engineering equipment will present the characteristics of mass, nonlinearity and dynamics. When using the traditional method to determine the health index of the equipment, it is difficult to extract the deep-level fault features of the equipment, so the accuracy of the health index will be affected to a certain extent. In order to fully tap the hidden features contained in massive data, Hinton proposed an effective tool for analyzing massive data—deep learning, which has been highly praised by scholars and engineers since then. Deep learning has been successfully applied in many engineering fields, such as image processing, speech recognition, machine translation, etc. At the same time, the extraction of fault features and the construction of health indicators based on this technology have received a lot of attention, and a series of results have been achieved.
由于受限玻尔兹曼机(Restricted Boltzmann Machine,RBM)的优良特性与特殊结构,可用于解决各种分类、降维、回归等问题。作为一种双向型概率图模型,受限玻尔兹曼机为一个具有两层结构的神经网络,主要由可视层与隐含层组成,满足层内无连接,层间全连接的要求。受限玻尔兹曼机优势在于能够以概率的形式定量表示出一组输入和输出的分布,在处理二进制数据具有一定的优势,但在面对连续数据建模时,模型训练效果并非令人满意,给具有连续性测试数据的工程设备故障特征提取提出了挑战。此外,拓扑结构是健康指标需要考虑的另外一类重要因素,单纯采用深度学习难以确保原始特征的拓扑结构。一旦拓扑结构发生改变,在降维过程中将会损失更多有价值的信息。因而,亟需引入其他降维方法,方便在保障拓扑结构的前提下确定出健康指标。Due to the excellent characteristics and special structure of Restricted Boltzmann Machine (RBM), it can be used to solve various classification, dimensionality reduction, regression and other problems. As a bidirectional probabilistic graphical model, the restricted Boltzmann machine is a neural network with a two-layer structure, which is mainly composed of a visible layer and a hidden layer, which meets the requirements of no connection within the layer and full connection between layers. The advantage of restricted Boltzmann machine is that it can quantitatively represent the distribution of a set of inputs and outputs in the form of probability. It has certain advantages in processing binary data, but when faced with continuous data modeling, the model training effect is not impressive Satisfied, it presents a challenge to the fault feature extraction of engineering equipment with continuous test data. In addition, topology is another important factor that needs to be considered in health indicators, and it is difficult to ensure the topology of original features using deep learning alone. Once the topology is changed, more valuable information will be lost in the process of dimensionality reduction. Therefore, it is urgent to introduce other dimensionality reduction methods to facilitate the determination of health indicators on the premise of ensuring the topology structure.
发明内容SUMMARY OF THE INVENTION
本发明的目的是提供一种工程设备健康指标构建方法及系统,采用连续深度置信网络提取蕴含在时频域背后的深层次故障特征,利用SOM神经网络在确保原始特征拓扑结构前提下构建设备的健康指标,以提高健康指标构建的准确性。The purpose of the present invention is to provide a method and system for constructing a health index of engineering equipment, adopting a continuous deep confidence network to extract the deep-level fault features contained in the time-frequency domain, and using the SOM neural network to construct the equipment's Health indicators to improve the accuracy of health indicator construction.
为实现上述目的,本发明提供了如下方案:For achieving the above object, the present invention provides the following scheme:
一种工程设备健康指标构建方法,包括:A method for constructing health indicators of engineering equipment, comprising:
获取工程设备的测试数据并计算所述测试数据的时频域信息;Obtain the test data of the engineering equipment and calculate the time-frequency domain information of the test data;
将所述时频域信息输入至连续深度置信网络,借助于对比散度算法对所述连续深度置信网络进行训练,生成训练好的连续深度置信网络;The time-frequency domain information is input into the continuous depth belief network, and the continuous depth belief network is trained by means of a contrastive divergence algorithm to generate a trained continuous depth belief network;
采用训练好的连续深度置信网络提取所述工程设备的深层次故障特征;Use the trained continuous deep belief network to extract the deep fault features of the engineering equipment;
将所述深层次故障特征作为SOM神经网络的输入进行无监督学习,生成训练好的SOM神经网络;The deep-level fault feature is used as the input of the SOM neural network for unsupervised learning to generate a trained SOM neural network;
通过计算所述训练好的SOM神经网络的输入层的状态向量与全部单元的权重向量间的欧氏距离,构建出所述工程设备的健康指标;所述健康指标直接反映所述工程设备的真实健康状态。By calculating the Euclidean distance between the state vector of the input layer of the trained SOM neural network and the weight vectors of all units, the health index of the engineering equipment is constructed; the health index directly reflects the actual situation of the engineering equipment health status.
可选地,所述获取工程设备的测试数据并计算所述测试数据的时频域信息,具体包括:Optionally, obtaining the test data of the engineering equipment and calculating the time-frequency domain information of the test data specifically includes:
获取工程设备的测试数据;Obtain test data of engineering equipment;
提取所述测试数据的时域特征;所述时域特征包括均值、均方根、方根幅值、平均绝对幅度、偏度、峭度、方差、最大值、最小值、峰值、波形因子、峰值因子、脉冲因子、裕度因子、偏度因子、峭度因子;Extract the time domain features of the test data; the time domain features include mean, root mean square, root square amplitude, mean absolute amplitude, skewness, kurtosis, variance, maximum value, minimum value, peak value, shape factor, Crest factor, impulse factor, margin factor, skewness factor, kurtosis factor;
提取所述测试数据的频域特征;所述频域特征包括中心频率、均方根频率、频率方差;Extract the frequency domain features of the test data; the frequency domain features include center frequency, root mean square frequency, and frequency variance;
将提取出的所述时域特征和所述频域特征构成所述工程设备的时频域信息。The time-domain features and the frequency-domain features extracted form the time-frequency domain information of the engineering equipment.
可选地,所述将所述时频域信息输入至连续深度置信网络,借助于对比散度算法对所述连续深度置信网络进行训练,生成训练好的连续深度置信网络,具体包括:Optionally, the time-frequency domain information is input into a continuous depth belief network, and the continuous depth belief network is trained by means of a contrastive divergence algorithm to generate a trained continuous depth belief network, specifically including:
构建由多个连续受限玻尔兹曼机堆叠而成的连续深度置信网络;所述连续受限玻尔兹曼机包含可视层、隐含层和高斯噪声单元;constructing a continuous deep belief network formed by stacking a plurality of continuous restricted Boltzmann machines; the continuous restricted Boltzmann machine includes a visible layer, a hidden layer and a Gaussian noise unit;
将所述时频域信息输入至所述连续深度置信网络,借助于对比散度算法更新所述连续深度置信网络的连接权重和sigmoid函数斜率控制参数,完成连续深度置信网络参数训练,生成训练好的连续深度置信网络。Input the time-frequency domain information into the continuous depth belief network, update the connection weight and sigmoid function slope control parameters of the continuous depth belief network with the help of the contrast divergence algorithm, complete the continuous depth belief network parameter training, and generate a trained A continuous deep belief network for .
可选地,所述采用训练好的连续深度置信网络提取所述工程设备的深层次故障特征,具体包括:Optionally, the use of a trained continuous deep belief network to extract deep-level fault features of the engineering equipment specifically includes:
将所述测试数据的时频域信息X=[x1,x2,…,xn]作为所述训练好的连续深度置信网络的输入,将所述训练好的连续深度置信网络的输出Y=[y1,y2,…,ym]作为所述工程设备的深层次故障特征;其中,xi(0≤i≤n)表示连续深度置信网络第i个输入神经元的特征序列,n为输入特征的总数;yi(0≤i≤m)表示连续深度置信网络第i个输出神经元的特征序列,m为连续深度置信网络输出神经元总数。The time-frequency domain information X=[x 1 , x 2 , . . . , x n ] of the test data is used as the input of the trained continuous deep belief network, and the output Y of the trained continuous deep belief network =[y 1 , y 2 ,...,y m ] as the deep fault feature of the engineering equipment; where x i (0≤i≤n) represents the feature sequence of the ith input neuron of the continuous deep belief network, n is the total number of input features; y i (0≤i≤m) represents the feature sequence of the ith output neuron of the continuous deep belief network, and m is the total number of output neurons of the continuous deep belief network.
可选地,所述将所述深层次故障特征作为SOM神经网络的输入进行无监督学习,生成训练好的SOM神经网络,具体包括:Optionally, the described deep-level fault feature is used as the input of the SOM neural network to carry out unsupervised learning, and the trained SOM neural network is generated, specifically including:
将所述深层次故障特征Y=[y1,y2,…,ym]中的一行作为所述SOM神经网络的输入向量D,根据所述输入向量D更新最佳匹配单元以及拓扑邻域的权重向量,完成所述SOM神经网络的无监督学习,生成训练好的SOM神经网络。Taking a row of the deep fault feature Y=[y 1 , y 2 ,..., y m ] as the input vector D of the SOM neural network, and updating the best matching unit and topological neighborhood according to the input vector D The weight vector of the SOM neural network is completed, and the unsupervised learning of the SOM neural network is completed, and the trained SOM neural network is generated.
一种工程设备健康指标构建系统,包括:A system for constructing health indicators of engineering equipment, comprising:
时频域信息提取模块,用于获取工程设备的测试数据并计算所述测试数据的时频域信息;a time-frequency domain information extraction module, used for acquiring the test data of the engineering equipment and calculating the time-frequency domain information of the test data;
连续深度置信网络训练模块,用于将所述时频域信息输入至连续深度置信网络,借助于对比散度算法对所述连续深度置信网络进行训练,生成训练好的连续深度置信网络;The continuous depth belief network training module is used to input the time-frequency domain information into the continuous depth belief network, train the continuous depth belief network by means of a contrastive divergence algorithm, and generate a trained continuous depth belief network;
深层次故障特征提取模块,用于采用训练好的连续深度置信网络提取所述工程设备的深层次故障特征;The deep-level fault feature extraction module is used for extracting the deep-level fault features of the engineering equipment by using the trained continuous deep belief network;
SOM神经网络训练模块,用于将所述深层次故障特征作为SOM神经网络的输入进行无监督学习,生成训练好的SOM神经网络;The SOM neural network training module is used to perform unsupervised learning with the deep-level fault features as the input of the SOM neural network to generate a trained SOM neural network;
工程设备健康指标构建模块,用于通过计算所述训练好的SOM神经网络的输入层的状态向量与全部单元的权重向量间的欧氏距离,构建出所述工程设备的健康指标;所述健康指标直接反映所述工程设备的真实健康状态。The engineering equipment health index building module is used to construct the health index of the engineering equipment by calculating the Euclidean distance between the state vector of the input layer of the trained SOM neural network and the weight vectors of all units; The indicators directly reflect the real health status of the engineering equipment.
可选地,所述时频域信息提取模块具体包括:Optionally, the time-frequency domain information extraction module specifically includes:
测试数据获取单元,用于获取工程设备的测试数据;The test data acquisition unit is used to acquire the test data of the engineering equipment;
时域特征提取单元,用于提取所述测试数据的时域特征;所述时域特征包括均值、均方根、方根幅值、平均绝对幅度、偏度、峭度、方差、最大值、最小值、峰值、波形因子、峰值因子、脉冲因子、裕度因子、偏度因子、峭度因子;A time-domain feature extraction unit for extracting time-domain features of the test data; the time-domain features include mean, root mean square, root square amplitude, mean absolute amplitude, skewness, kurtosis, variance, maximum value, Minimum value, peak value, shape factor, crest factor, pulse factor, margin factor, skewness factor, kurtosis factor;
频域特征提取单元,用于提取所述测试数据的频域特征;所述频域特征包括中心频率、均方根频率、频率方差;a frequency domain feature extraction unit, used for extracting the frequency domain feature of the test data; the frequency domain feature includes a center frequency, a root mean square frequency, and a frequency variance;
时频域信息构建单元,用于将提取出的所述时域特征和所述频域特征构成所述工程设备的时频域信息。The time-frequency domain information construction unit is configured to form the time-frequency domain information of the engineering equipment from the extracted time-domain features and the frequency-domain features.
可选地,所述连续深度置信网络训练模块具体包括:Optionally, the continuous deep belief network training module specifically includes:
连续深度置信网络构建单元,用于构建由多个连续受限玻尔兹曼机堆叠而成的连续深度置信网络;所述连续受限玻尔兹曼机包含可视层、隐含层和高斯噪声单元;A continuous deep belief network building unit for constructing a continuous deep belief network composed of multiple continuous restricted Boltzmann machines stacked; the continuous restricted Boltzmann machine includes a visible layer, a hidden layer and a Gaussian noise unit;
连续深度置信网络训练单元,用于将所述时频域信息输入至所述连续深度置信网络,借助于对比散度算法更新所述连续深度置信网络的连接权重和sigmoid函数斜率控制参数,完成连续深度置信网络参数训练,生成训练好的连续深度置信网络。The continuous depth belief network training unit is used to input the time-frequency domain information into the continuous depth belief network, and update the connection weight and the sigmoid function slope control parameter of the continuous depth belief network by means of a contrastive divergence algorithm to complete the continuous depth belief network. Deep belief network parameter training to generate a trained continuous deep belief network.
可选地,所述深层次故障特征提取模块具体包括:Optionally, the deep-level fault feature extraction module specifically includes:
深层次故障特征提取单元,用于将所述测试数据的时频域信息X=[x1,x2,…,xn]作为所述训练好的连续深度置信网络的输入,将所述训练好的连续深度置信网络的输出Y=[y1,y2,…,ym]作为所述工程设备的深层次故障特征;其中,xi(0≤i≤n)表示连续深度置信网络第i个输入神经元的特征序列,n为输入特征的总数;yi(0≤i≤m)表示连续深度置信网络第i个输出神经元的特征序列,m为连续深度置信网络输出神经元总数。A deep-level fault feature extraction unit, configured to use the time-frequency domain information X=[x 1 , x 2 , . . . , x n ] of the test data as the input of the trained continuous deep belief network, The output of a good continuous deep belief network Y=[y 1 , y 2 ,...,y m ] is used as the deep fault feature of the engineering equipment; wherein, x i (0≤i≤n) represents the first continuous deep belief network The feature sequence of i input neurons, n is the total number of input features; y i (0≤i≤m) represents the feature sequence of the ith output neuron of the continuous deep belief network, m is the total number of output neurons of the continuous deep belief network .
可选地,所述SOM神经网络训练模块具体包括:Optionally, the SOM neural network training module specifically includes:
SOM神经网络训练单元,用于将所述深层次故障特征Y=[y1,y2,…,ym]中的一行作为所述SOM神经网络的输入向量D,根据所述输入向量D更新最佳匹配单元以及拓扑邻域的权重向量,完成所述SOM神经网络的无监督学习,生成训练好的SOM神经网络。The SOM neural network training unit is configured to use a row of the deep fault features Y=[y 1 , y 2 ,..., y m ] as the input vector D of the SOM neural network, and update according to the input vector D The weight vector of the best matching unit and topological neighborhood completes the unsupervised learning of the SOM neural network, and generates a trained SOM neural network.
根据本发明提供的具体实施例,本发明公开了以下技术效果:According to the specific embodiments provided by the present invention, the present invention discloses the following technical effects:
本发明提供了一种工程设备健康指标构建方法及系统,所述方法包括:获取工程设备的测试数据并计算所述测试数据的时频域信息;将所述时频域信息输入至连续深度置信网络,借助于对比散度算法对所述连续深度置信网络进行训练,生成训练好的连续深度置信网络;采用训练好的连续深度置信网络提取所述工程设备的深层次故障特征;将所述深层次故障特征作为SOM神经网络的输入进行无监督学习,生成训练好的SOM神经网络;通过计算所述训练好的SOM神经网络的输入层的状态向量与全部单元的权重向量间的欧氏距离,构建出所述工程设备的健康指标;所述健康指标直接反映所述工程设备的真实健康状态。本发明方法及系统采用连续深度置信网络提取蕴含在时频域背后的深层次故障特征,利用SOM神经网络在确保原始特征拓扑结构前提下构建设备的健康指标,提高了健康指标构建的准确性。The present invention provides a method and system for constructing health indicators of engineering equipment. The method includes: acquiring test data of engineering equipment and calculating time-frequency domain information of the test data; inputting the time-frequency domain information into continuous depth confidence network, train the continuous depth belief network by means of the contrast divergence algorithm to generate a trained continuous depth belief network; use the trained continuous depth belief network to extract the deep fault features of the engineering equipment; The hierarchical fault feature is used as the input of the SOM neural network for unsupervised learning to generate a trained SOM neural network; by calculating the Euclidean distance between the state vector of the input layer of the trained SOM neural network and the weight vector of all units, The health index of the engineering equipment is constructed; the health index directly reflects the real health state of the engineering equipment. The method and system of the present invention use a continuous deep confidence network to extract the deep fault features contained in the time-frequency domain, and use the SOM neural network to construct the health index of the equipment on the premise of ensuring the original feature topology structure, thereby improving the accuracy of the construction of the health index.
附图说明Description of drawings
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the accompanying drawings required in the embodiments will be briefly introduced below. Obviously, the drawings in the following description are only some of the present invention. In the embodiments, for those of ordinary skill in the art, other drawings can also be obtained according to these drawings without creative labor.
图1为本发明一种工程设备健康指标构建方法的流程图;Fig. 1 is the flow chart of a kind of construction method of engineering equipment health index of the present invention;
图2为本发明实施例提供的轴承实验平台的示意图;2 is a schematic diagram of a bearing experimental platform provided by an embodiment of the present invention;
图3为本发明实施例提供的轴承在全寿命周期内的典型水平振动信号示意图;3 is a schematic diagram of a typical horizontal vibration signal of a bearing provided in an embodiment of the present invention in a full life cycle;
图4为本发明实施例提供的三种运行条件中第一种运行条件下轴承的健康指标示意图;FIG. 4 is a schematic diagram of the health index of the bearing under the first operating condition among the three operating conditions provided by the embodiment of the present invention;
图5为本发明实施例提供的三种运行条件中第二种运行条件下轴承的健康指标示意图;FIG. 5 is a schematic diagram of the health index of the bearing under the second operating condition among the three operating conditions provided by the embodiment of the present invention;
图6为本发明实施例提供的三种运行条件中第三种运行条件下轴承的健康指标示意图。FIG. 6 is a schematic diagram of the health index of the bearing under the third operating condition among the three operating conditions provided by the embodiment of the present invention.
具体实施方式Detailed ways
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, but not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.
本发明针对大数据背景下工程设备,提出了一种具有普适性的故障特征提取与健康指标构建方法及系统,采用连续深度置信网络提取蕴含在时频域背后的深层次故障特征,利用SOM神经网络在确保原始特征拓扑结构前提下构建设备的健康指标,以改善健康指标的特性,提高健康指标构建的准确性。Aiming at engineering equipment under the background of big data, the present invention proposes a universal fault feature extraction and health index construction method and system. The continuous deep confidence network is used to extract the deep fault features contained in the time-frequency domain, and the SOM is used. The neural network constructs the health index of the device under the premise of ensuring the original feature topology structure, so as to improve the characteristics of the health index and improve the accuracy of the construction of the health index.
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。In order to make the above objects, features and advantages of the present invention more clearly understood, the present invention will be described in further detail below with reference to the accompanying drawings and specific embodiments.
图1为本发明一种工程设备健康指标构建方法的流程图。参见图1,本发明一种工程设备健康指标构建方法包括:FIG. 1 is a flowchart of a method for constructing a health index of engineering equipment according to the present invention. Referring to Fig. 1, a construction method of an engineering equipment health index of the present invention includes:
步骤1:获取工程设备的测试数据并计算所述测试数据的时频域信息。Step 1: Acquire test data of engineering equipment and calculate time-frequency domain information of the test data.
本发明利用信号处理技术对所收集的海量测试数据进行分析,确定相应的时域与频域特征,典型的时域特征主要有均值、均方根、方根幅值、平均绝对幅度、偏度、峭度、方差、最大值、最小值、峰值、波形因子、峰值因子、脉冲因子、裕度因子、偏度因子、峭度因子;频域特征主要包括中心频率、均方根频率、频率方差等。The present invention uses signal processing technology to analyze the collected mass test data, and determines the corresponding time domain and frequency domain characteristics. Typical time domain characteristics mainly include mean value, root mean square, root square amplitude, mean absolute amplitude and skewness. , kurtosis, variance, maximum value, minimum value, peak value, shape factor, crest factor, impulse factor, margin factor, skewness factor, kurtosis factor; frequency domain features mainly include center frequency, root mean square frequency, frequency variance Wait.
因此,所述步骤1获取工程设备的测试数据并计算所述测试数据的时频域信息,具体包括:Therefore, the
步骤1.1:获取工程设备的测试数据。Step 1.1: Obtain test data of engineering equipment.
所述工程设备包括轴承等机械设备,以及芯片、三极管、MOS管等电子设备。The engineering equipment includes mechanical equipment such as bearings, and electronic equipment such as chips, triodes, and MOS tubes.
步骤1.2:提取所述测试数据的时域特征;所述时域特征包括均值、均方根、方根幅值、平均绝对幅度、偏度、峭度、方差、最大值、最小值、峰值、波形因子、峰值因子、脉冲因子、裕度因子、偏度因子、峭度因子。Step 1.2: Extract the time domain features of the test data; the time domain features include mean, root mean square, root square amplitude, mean absolute amplitude, skewness, kurtosis, variance, maximum value, minimum value, peak value, Shape factor, crest factor, impulse factor, margin factor, skewness factor, kurtosis factor.
各时域特征的计算公式如下:The calculation formula of each time domain feature is as follows:
最大值xmax=maxx(t) (8)max x max = maxx(t) (8)
最小值xmin=minx(t) (9)Minimum x min = minx(t) (9)
峰值xpeak=max|x(t)| (10)Peak x peak =max|x(t)| (10)
其中,xi表示工程设备在单次采样间隔内的第i个测试数据;N为单次采样间隔内的测试数据数目,x(t)表示工程设备在单次采样间隔内所有数据集合。Among them, x i represents the ith test data of the engineering equipment in a single sampling interval; N is the number of test data in a single sampling interval, and x(t) represents all the data sets of the engineering equipment in a single sampling interval.
步骤1.3:提取所述测试数据的频域特征;所述频域特征主要包括中心频率、均方根频率、频率方差等。为方便理解频域特征,以中心频率、均方根频率、频率方差为例,详细阐述其计算公式:Step 1.3: Extract the frequency domain features of the test data; the frequency domain features mainly include center frequency, root mean square frequency, frequency variance, and the like. In order to facilitate the understanding of the frequency domain characteristics, the calculation formulas are described in detail by taking the center frequency, root mean square frequency and frequency variance as examples:
将采集的工程设备测试数据进行傅里叶变换,能够从时域信号中提取出所包含的频谱信息,即频率值与相应的峰值/幅值,频域特征可以基于频谱信息计算得到。其中,fi为第i时刻频谱对应的频率值,pi表示第i时刻频域对应的幅值,fc为平均频率。其他频域特征可根据现有文献,并结合上述频域特征确定。Performing Fourier transform on the collected engineering equipment test data can extract the spectral information contained in the time domain signal, that is, the frequency value and the corresponding peak/amplitude value, and the frequency domain feature can be calculated based on the spectral information. Among them, f i is the frequency value corresponding to the spectrum at the i -th time, pi is the amplitude corresponding to the frequency domain at the i-th time, and f c is the average frequency. Other frequency domain features can be determined according to the existing literature and in combination with the above frequency domain features.
步骤1.4:将提取出的所述时域特征和所述频域特征构成所述工程设备的时频域信息。Step 1.4: The time-frequency domain information of the engineering equipment is constituted by the extracted time-domain features and the frequency-domain features.
将每个监测时刻的时频域特征组成一个特征向量,对于多个监测时刻,不同特征向量按照时间顺序组成一个特征矩阵,即构成所述工程设备的时频域信息,作为连续深度置信网络的输入。The time-frequency domain features of each monitoring moment are formed into a eigenvector, and for multiple monitoring moments, different eigenvectors form a feature matrix in time order, that is, the time-frequency domain information of the engineering equipment, as the continuous depth confidence network. enter.
步骤2:将所述时频域信息输入至连续深度置信网络,借助于对比散度算法对所述连续深度置信网络进行训练,生成训练好的连续深度置信网络。Step 2: Input the time-frequency domain information into a continuous deep belief network, train the continuous deep belief network by means of a contrastive divergence algorithm, and generate a trained continuous deep belief network.
本发明将所提取的时频域信息输入至连续深度置信网络(Continuous DeepBelief Network,CDBN),借助于对比散度算法对网络模型进行训练,训练完成后可提取出工程设备的深层次故障特征。In the present invention, the extracted time-frequency domain information is input into a continuous deep belief network (Continuous Deep Belief Network, CDBN), the network model is trained by means of a contrastive divergence algorithm, and the deep-level fault features of the engineering equipment can be extracted after the training is completed.
所述步骤2将所述时频域信息输入至连续深度置信网络,借助于对比散度算法对所述连续深度置信网络进行训练,生成训练好的连续深度置信网络,具体包括:In the
步骤2.1:构建由多个连续受限玻尔兹曼机堆叠而成的连续深度置信网络;所述连续受限玻尔兹曼机包含可视层、隐含层和高斯噪声单元。Step 2.1: Construct a continuous deep belief network formed by stacking multiple continuous restricted Boltzmann machines; the continuous restricted Boltzmann machine includes a visible layer, a hidden layer and a Gaussian noise unit.
本发明构建的连续深度置信网络是由多个连续受限玻尔兹曼机(ContinuousRestricted Boltzmann Machine,CRBM)堆叠而成的神经网络,连续受限玻尔兹曼机是在传统受限玻尔兹曼机的基础上,通过在采样Sigmoid单元中增加一个零均值的高斯噪声构成,主要包含可视层、隐含层和高斯噪声单元。令sj表示以若干神经元状态集合{si}为输入的神经元j输出状态,则sj可表示为:The continuous depth belief network constructed by the present invention is a neural network formed by stacking a plurality of Continuous Restricted Boltzmann Machines (CRBM). On the basis of the Mann machine, it is formed by adding a zero-mean Gaussian noise to the sampling Sigmoid unit, which mainly includes the visible layer, the hidden layer and the Gaussian noise unit. Let s j denote the output state of neuron j with several neuron state sets {s i } as input, then s j can be expressed as:
其中,Wij为连接权重,Nj(0,1)表示均值为0、方差为1的高斯白噪声。σ为常值项,与Nj(0,1)共同组成了噪声输入部分,即nj=σ·Nj(0,1),相应的概率密度密度p(nj)为:Among them, W ij is the connection weight, and N j (0,1) represents Gaussian white noise with
φj(xj)为更为一般的sigmoid函数,且满足xj为sigmoid函数作用之前的状态,满足θL与θH分别表示φj(xj)函数的下渐近线与上渐近线,cj为一个控制sigmoid函数斜率的参数,称为sigmoid函数斜率控制参数。φ j (x j ) is a more general sigmoid function that satisfies x j is the state before the sigmoid function acts, satisfying θ L and θ H represent the lower asymptote and the upper asymptote of the φ j (x j ) function, respectively, and c j is a parameter that controls the slope of the sigmoid function, which is called the sigmoid function slope control parameter.
步骤2.2:将所述时频域信息输入至所述连续深度置信网络,借助于对比散度算法更新所述连续深度置信网络的连接权重和sigmoid函数斜率控制参数,完成连续深度置信网络参数训练,生成训练好的连续深度置信网络。Step 2.2: Input the time-frequency domain information into the continuous depth belief network, update the connection weight of the continuous depth belief network and the sigmoid function slope control parameter by means of the contrast divergence algorithm, and complete the continuous depth belief network parameter training, Generate a trained continuous deep belief network.
步骤2.1中的{si}表示的是连续受限玻尔兹曼机中输入层每个神经元,为了方便描述连续受限玻尔兹曼机的运算机制,所以统一用{si}表示,当时域频域特征(即所述时频域信息)输入至连续深度置信网络后,即第一个连续受限玻尔兹曼机输入层直接接收该时域频域特征,因此这里的{si}需要用时域频域特征进行代替,但对于后面的连续受限玻尔兹曼机,{si}需要用上一个连续受限玻尔兹曼机的输出进行代替。{s i } in step 2.1 represents each neuron in the input layer of the continuous restricted Boltzmann machine. For the convenience of describing the operation mechanism of the continuous restricted Boltzmann machine, it is uniformly represented by {s i } , after the time-frequency domain features (that is, the time-frequency domain information) are input to the continuous deep confidence network, that is, the first continuous restricted Boltzmann machine input layer directly receives the time-frequency domain features, so here { s i } needs to be replaced by time domain and frequency domain features, but for the following continuous restricted Boltzmann machine, {s i } needs to be replaced by the output of the previous continuous restricted Boltzmann machine.
当可视层的神经元输入特定状态,根据式(20)可得到隐含层神经元的状态,基于隐含层神经元的状态可重构出可视层状态,利用同样方式可确定出重构后隐含层神经元的状态,按照此方式进行类推,经过无穷多次重构,可实现模型参数的训练。但这种方式同样需要多次Gibbs采样(吉布斯采样),耗时较长,不利于模型实时求解。本发明采用一种最小化对比散度算法用于模型参数训练,显著降低了计算量,为模型参数高效训练提供了便捷途径。连续深度置信网络权重与sigmoid函数斜率控制参数可根据下式进行更新:When the neurons in the visual layer input a specific state, the state of the neurons in the hidden layer can be obtained according to formula (20), and the state of the visual layer can be reconstructed based on the state of the neurons in the hidden layer. The state of the neurons in the hidden layer after the structure is constructed, and the analogy is carried out in this way. After infinite reconstructions, the training of the model parameters can be realized. However, this method also requires multiple Gibbs sampling (Gibbs sampling), which takes a long time and is not conducive to the real-time solution of the model. The invention adopts a minimization contrast divergence algorithm for model parameter training, which significantly reduces the amount of calculation and provides a convenient way for efficient training of model parameters. The continuous depth belief network weight and sigmoid function slope control parameters can be updated according to the following formula:
式中,εw与εc分别表示网络连接权重Wij和sigmoid函数斜率控制参数cj的学习率。sisj为神经元i与神经元j当前状态的乘积。分别神经元i与神经元j一步采样状态值。对于单个连续受限玻尔兹曼机而言,当隐含层神经元总数小于可视层神经元总数时,能够有效降低输入状态的维数。ΔWij表示经过一步采样后Wij的变化量,Δcj表示经过一步采样后cj的变化量,经过多步采样后,可实现参数更新。具体地,通过将ΔWij和Δcj分别与上一次训练的参数Wij和cj相加来进行参数更新,从而完成连续深度置信网络参数训练,生成训练好的连续深度置信网络。In the formula, εw and εc represent the network connection weight W ij and the learning rate of the sigmoid function slope control parameter c j , respectively. s i s j is the product of the current state of neuron i and neuron j. The state values of neuron i and neuron j are sampled in one step, respectively. For a single continuous restricted Boltzmann machine, when the total number of neurons in the hidden layer is less than the total number of neurons in the visual layer, the dimension of the input state can be effectively reduced. ΔW ij represents the change of W ij after one-step sampling, Δc j represents the change of c j after one-step sampling, and parameter update can be realized after multi-step sampling. Specifically, the parameters are updated by adding ΔW ij and Δc j to the parameters W ij and c j of the previous training respectively, so as to complete the continuous deep belief network parameter training and generate a trained continuous deep belief network.
步骤3:采用训练好的连续深度置信网络提取所述工程设备的深层次故障特征。Step 3: Using the trained continuous deep belief network to extract the deep fault features of the engineering equipment.
假定退化设备测试信号的原始特征,即提取的所述测试数据的时频域信息可表示为X=[x1,x2,…,xn],其中,xi(0≤i≤n)表示连续深度置信网络第i个输入神经元的特征序列,n为原始特征的总数,将这些原始特征当作连续深度置信网络的输入,采用上述对比散度算法逐层确定每个连续受限玻尔兹曼机的模型参数,基于式(20)可计算出连续深度置信网络各层神经元的状态,即通过多重连续受限玻尔兹曼机可提取出反映设备健康状态的深层次故障特征。因而,连续深度置信网络的输出可视为输入数据的深层次故障特征,可表示为Y=[y1,y2,…,ym],其中,yi(0≤i≤m)表示连续深度置信网络第i个输出神经元的特征序列,m为连续深度置信网络输出神经元总数。Assuming the original characteristics of the test signal of the degraded device, that is, the extracted time-frequency domain information of the test data can be expressed as X=[x 1 , x 2 , . . . , x n ], where x i (0≤i≤n) Represents the feature sequence of the ith input neuron of the continuous deep belief network, n is the total number of original features, these original features are regarded as the input of the continuous deep belief network, and the above contrastive divergence algorithm is used to determine each continuous restricted glass layer by layer. The model parameters of the Boltzmann machine can be calculated based on the formula (20) to calculate the state of each layer of neurons in the continuous deep belief network, that is, the deep fault features reflecting the health status of the equipment can be extracted through multiple continuous restricted Boltzmann machines. . Therefore, the output of the continuous deep belief network can be regarded as the deep fault feature of the input data, which can be expressed as Y=[y 1 , y 2 ,...,y m ], where y i (0≤i≤m) represents the continuous The feature sequence of the ith output neuron of the deep belief network, m is the total number of output neurons of the continuous deep belief network.
测试数据的时频域信息X=[x1,x2,…,xn]中包含了全部数据的时域与频域特征,i表示的特征序号,xi(0≤i≤n)第i个特征的序列值,因为在多个时间点采样,所以是一个序列,每个特征直接对应连续深度置信网络输入的一个神经元。X中的每一行代表了某时刻时域和频域特征组成的向量,每一列代表的是某个特征在不同时刻的值,所以X输入到连续深度置信网络中后,X的每一行对应了连续深度置信网络第一个连续受限玻尔兹曼机的输入神经元,即与{si}是相对应的,而对于后面的连续受限玻尔兹曼机,{si}仅表示上一层连续受限玻尔兹曼机的输出。The time-frequency domain information X=[x 1 , x 2 ,..., x n ] of the test data contains the time domain and frequency domain features of all the data, i represents the feature sequence number, x i (0≤i≤n)th The sequence values of i features are a sequence because they are sampled at multiple time points, and each feature directly corresponds to a neuron input to the continuous deep belief network. Each row in X represents a vector composed of time-domain and frequency-domain features at a certain time, and each column represents the value of a feature at different times, so after X is input into the continuous deep belief network, each row of X corresponds to The input neuron of the first continuous restricted Boltzmann machine of the continuous deep belief network, that is, corresponding to {s i }, and for the following continuous restricted Boltzmann machine, {s i } only represents The output of the continuous restricted Boltzmann machine of the previous layer.
步骤4:将所述深层次故障特征作为SOM神经网络的输入进行无监督学习,生成训练好的SOM神经网络。Step 4: Use the deep-level fault features as the input of the SOM neural network to perform unsupervised learning to generate a trained SOM neural network.
本发明将连续深度置信网络所提取的故障特征作为SOM(Self-organizing map,自组织映射)神经网络的输入以实现其无监督学习,通过计算状态向量与全部单元的权重向量间的欧氏距离,定量刻画设备的健康状态。In the present invention, the fault feature extracted by the continuous deep confidence network is used as the input of the SOM (Self-organizing map, self-organizing map) neural network to realize its unsupervised learning, and the Euclidean distance between the state vector and the weight vector of all units is calculated by calculating the Euclidean distance. , quantitatively characterize the health status of the equipment.
本发明构建的SOM神经网络主要包含两部分:输入层与输出层。输入层主要接收实际中高维非线性数据,也称为竞争层;输出层由规则节点网格所构成。输出层第l个神经元与输入神经元之间的权重向量可表示为Wl=[Wl1,Wl2,…,Wlm],m为SOM神经网络输入向量的维数。对于步骤3提取的深层次故障特征Y的某一行,即特定时刻来自连续深度置信网络的深度特征可看作SOM神经网络输入向量,可记为D。根据所述输入向量D更新最佳匹配单元以及拓扑邻域的权重向量,完成所述SOM神经网络的无监督学习,从而生成训练好的SOM神经网络。The SOM neural network constructed by the present invention mainly includes two parts: an input layer and an output layer. The input layer mainly receives actual medium and high-dimensional nonlinear data, also known as the competition layer; the output layer is composed of regular node grids. The weight vector between the lth neuron of the output layer and the input neuron can be expressed as W l =[W l1 ,W l2 ,...,W lm ], where m is the dimension of the input vector of the SOM neural network. For a row of the deep fault feature Y extracted in
学习过程主要是计算输入向量和权重向量之间的欧式距离,以确定每步训练中最接近输入向量的神经元。一般将最接近的神经元称为最佳匹配单元(BestMatchingUnit,BMU),基于获得的BMU,能够调整更新BMU以及拓扑邻域的权重向量,更新结果可表示为:The learning process is mainly to calculate the Euclidean distance between the input vector and the weight vector to determine the neuron closest to the input vector in each step of training. Generally, the closest neuron is called the Best Matching Unit (BMU). Based on the obtained BMU, the weight vector of the BMU and the topological neighborhood can be adjusted and updated. The update result can be expressed as:
式中,p为迭代步长的索引,Wl(p)和Wl(p+1)分别为p步和p+1步的权重向量,D为输入向量,ε(p)为学习率,表示以BMU为中心的拓扑邻域核函数,WBMU为BMU的权重向量。In the formula, p is the index of the iteration step, W l (p) and W l (p+1) are the weight vectors of p steps and p+1 steps respectively, D is the input vector, ε(p) is the learning rate, Represents the topological neighborhood kernel function centered on the BMU, and W BMU is the weight vector of the BMU.
步骤5:通过计算所述训练好的SOM神经网络的输入层的状态向量与全部单元的权重向量间的欧氏距离,构建出所述工程设备的健康指标。Step 5: Construct the health index of the engineering equipment by calculating the Euclidean distance between the state vector of the input layer of the trained SOM neural network and the weight vector of all units.
在进行健康评估过程中,利用健康状态数据对SOM结构进行训练,训练完成后将输入层的状态向量与全部单元的权重向量进行比较,用于评价设备的健康状态。具体地,本发明选择最小量化误差作为健康指标,即通过计算状态向量与全部单元的权重向量间的欧氏距离,定量刻画设备的健康状态,即:In the process of health assessment, the SOM structure is trained by using the health state data. After the training is completed, the state vector of the input layer is compared with the weight vector of all units to evaluate the health state of the equipment. Specifically, the present invention selects the minimum quantization error as the health index, that is, by calculating the Euclidean distance between the state vector and the weight vectors of all units, quantitatively characterizes the health state of the device, that is:
Z=min||D-Wl|| (25)Z=min||DW l || (25)
由于不同时刻输入向量对应了不同的健康指标,利用式(25)可得到一系列健康指标序列Z=[Z(t1),Z(t2),…,Z(tN)],其中Z(ti)(1≤i≤N)表示ti时刻的健康指标。需要说明的是,SOM神经网络能够将多维数据转换为一维的健康指标,该健康指标直接反映了设备的真实健康状态。基于该健康指标可以进行故障诊断与寿命预测。Since the input vectors at different times correspond to different health indicators, a series of health indicator sequences Z=[Z(t 1 ), Z(t 2 ), ..., Z(t N )] can be obtained by using formula (25), where Z(t N )] (t i ) (1≤i≤N) represents the health index at time t i . It should be noted that the SOM neural network can convert multi-dimensional data into one-dimensional health indicators, which directly reflect the real health status of the equipment. Based on this health index, fault diagnosis and life prediction can be performed.
本发明提取的健康指标序列Z=[Z(t1),Z(t2),…,Z(tN)]是一个随时间变化的量,若该指标随时间有上升的趋势,意味着设备出现了退化趋势,表明健康状态会越来越差。该健康指标是借助于连续深度置信网络和SOM神经网络融合设备的所有时域频域特征,所提取出的一个综合指标,因此能够对原始特征进行抽象化表达。下面轴承的具体实施例可反映这一个过程。The health index sequence Z=[Z(t 1 ), Z(t 2 ),..., Z(t N )] extracted by the present invention is a quantity that changes with time. The equipment is showing a trend of degradation, indicating that the health status is getting worse and worse. The health index is a comprehensive index extracted by integrating all the time-domain and frequency-domain features of the device with the continuous deep belief network and the SOM neural network, so the original features can be abstractly expressed. The following specific examples of bearings may reflect this process.
下面采用公开的轴承实验数据对所述方法的有效性与优越性进行验证。该轴承实验数据组来源于西安交通大学设计科学与基础部件研究所与长兴昇阳科技有限公司共同建立的机械装备健康监测联合实验室。轴承实验平台主要由交流电机、电机速度控制器、支撑轴、两个支撑轴承以及液压加载装置组成,如图2所示。实验台的传感装置为垂直安装在实验轴承外圈的两个加速度计,主要功能是采集水平与垂直方向的振动信号。单次采样频率为25.6kHz,采样持续时间为1.28s,同时相邻两次采样间隔为1min。The effectiveness and superiority of the method are verified by using the published bearing experimental data below. The bearing experimental data set comes from the joint laboratory of mechanical equipment health monitoring jointly established by the Institute of Design Science and Basic Components of Xi'an Jiaotong University and Changxing Shengyang Technology Co., Ltd. The bearing experimental platform is mainly composed of an AC motor, a motor speed controller, a support shaft, two support bearings and a hydraulic loading device, as shown in Figure 2. The sensing devices of the experimental bench are two accelerometers vertically installed on the outer ring of the experimental bearing. The main function is to collect vibration signals in the horizontal and vertical directions. The single sampling frequency is 25.6kHz, the sampling duration is 1.28s, and the interval between two adjacent samplings is 1min.
该轴承加速退化实验是在三种实验条件下进行的,每种实验条件下测试5组轴承,因而数据组主要包含15组轴承的全寿命周期振动测试数据,表1罗列了实验轴承运行条件与寿命具体信息。The accelerated degradation test of the bearing was carried out under three experimental conditions, and 5 sets of bearings were tested under each experimental condition. Therefore, the data set mainly contains the whole life cycle vibration test data of 15 sets of bearings. Table 1 lists the operating conditions of the experimental bearings and Lifetime specific information.
表1实验轴承运行条件与寿命具体信息Table 1 Specific information on the operating conditions and life of the experimental bearing
由于实验的载荷是施加于水平方向的,安装于水平方向的加速度计能够获取更多的轴承退化信息,因此这里选择水平振动信号对测试轴承的RUL进行自适应预测。图3描绘了轴承在全寿命周期内的典型水平振动信号。从图3中能够直观看出,轴承的健康状态通常可分为两个阶段:稳定运行阶段与退化阶段。在稳定运行阶段,轴承的振动信号围绕0值小幅度上下波动,不存在任何退化趋势,因而轴承在此阶段发生退化的情形可忽略不计;然而在退化阶段,振动信号的幅值随着时间推移具有明显上升趋势,意味着在此阶段能够获取极其丰富的轴承退化信息。Since the load of the experiment is applied in the horizontal direction, the accelerometer installed in the horizontal direction can obtain more bearing degradation information, so the horizontal vibration signal is selected here to adaptively predict the RUL of the tested bearing. Figure 3 depicts a typical horizontal vibration signal over the life of a bearing. It can be seen intuitively from Figure 3 that the health state of the bearing can usually be divided into two stages: the stable operation stage and the degradation stage. In the stable operation stage, the vibration signal of the bearing fluctuates up and down in a small range around the 0 value, and there is no degradation trend, so the degradation of the bearing at this stage can be ignored; however, in the degradation stage, the amplitude of the vibration signal changes with time. There is a clear upward trend, which means that extremely rich bearing degradation information can be obtained at this stage.
为充分利用振动信号中的特征信息,这里选择步骤1.2所述的16个时域特征和13个频域特征(包括步骤1.3所述的中心频率、均方根频率、频率方差以及其余10个本领域熟知的其他频域特征),并将其输入值至包含2层CRBM结构的CDBN,以提取振动信号中蕴含的深层次故障特征。CDBN和SOM神经网络具体结构和初始参数由交叉验证确定,如表2与表3所列。In order to make full use of the feature information in the vibration signal, the 16 time-domain features and 13 frequency-domain features described in step 1.2 (including the center frequency, root mean square frequency, frequency variance and the remaining 10 parameters described in step 1.3) are selected here. Other frequency domain features well-known in the field), and input the value to the CDBN containing the 2-layer CRBM structure to extract the deep fault features contained in the vibration signal. The specific structures and initial parameters of CDBN and SOM neural networks are determined by cross-validation, as listed in Table 2 and Table 3.
表2 CDBN的具体结构以及初始参数Table 2 The specific structure and initial parameters of CDBN
表3 SOM的具体结构以及初始参数Table 3 Specific structure and initial parameters of SOM
基于CDBN提取出反映轴承健康状态的隐含深层故障特征后,可利用SOM神经网络确定三种运行条件下轴承的健康指标,分别如图4、图5、图6所示。因而,29个特征将会简化为一维特征。图4表明轴承的健康指标趋势与振动信号的幅度趋势一致,在一定程度上证实了轴承的两阶段健康状态。因此,本发明所述方法在提取轴承的深层次故障特征和确定健康指标方面是合理且有效的。After extracting the hidden deep fault features reflecting the bearing health status based on CDBN, the SOM neural network can be used to determine the bearing health indicators under three operating conditions, as shown in Figure 4, Figure 5, and Figure 6, respectively. Thus, the 29 features will be reduced to one-dimensional features. Figure 4 shows that the trend of the health index of the bearing is consistent with the amplitude trend of the vibration signal, which confirms the two-stage health state of the bearing to a certain extent. Therefore, the method described in the present invention is reasonable and effective in extracting deep fault features of bearings and determining health indicators.
为进一步体现本发明提出的融合了连续深度置信网络与SOM神经网络方法的优越性,这里分别采用传统的降维方法作为比较方法来确定健康指标,如PCA(PrincipleComponentAnalysis,主成分分析)、SOM神经网络和连续深度置信网络(CDBN)。同时,选择健康指标趋势因子来评估健康指标的质量。一般地,趋势因子越接近1,则趋势性越强。表4罗列了所提方法与三种对比方法下各轴承健康指标的趋势因子结果。In order to further reflect the superiority of the method of integrating the continuous depth belief network and the SOM neural network proposed by the present invention, the traditional dimensionality reduction method is used as a comparison method to determine the health indicators, such as PCA (Principle Component Analysis, principal component analysis), SOM neural network. Networks and Continuous Deep Belief Networks (CDBN). At the same time, the health indicator trend factor is selected to evaluate the quality of the health indicator. In general, the closer the trend factor is to 1, the stronger the trend. Table 4 lists the trend factor results of each bearing health index under the proposed method and the three comparison methods.
表4四种方法下健康指标性能比较Table 4 Comparison of health index performance under four methods
从表4数据可以发现,本发明方法计算的健康指标趋势因子大于其他三种方法下健康指标的趋势因子,意味着采用本发明方法确定的健康指标优于三种传统的降维方法。此外,基于本发明方法所确定大多数轴承健康指标的趋势因子均大于0.5,为健康指标的退化建模与RUL(Remaining Useful Life,剩余使用寿命)预测提供了技术支撑。From the data in Table 4, it can be found that the trend factor of the health index calculated by the method of the present invention is greater than the trend factor of the health index under the other three methods, which means that the health index determined by the method of the present invention is better than the three traditional dimensionality reduction methods. In addition, the trend factors of most bearing health indicators determined based on the method of the present invention are all greater than 0.5, which provides technical support for the degradation modeling of health indicators and RUL (Remaining Useful Life, remaining useful life) prediction.
本发明提出的一种新型工程设备健康指标构建方法,对于大数据背景下的工程设备,能够提取出蕴含在时频域背后的深层次故障特征,在确保原始特征拓扑结构的前提下构建设备的健康指标。本发明方法同现有技术相比,充分提取了蕴含在时频域信息背后的深层次故障特征,同时确保了降维过程中原始特征的拓扑结构,进一步发展了故障特征提取方法,丰富了设备健康指标构建理论体系,有效改善了健康指标的性能特征。The method for constructing a new type of engineering equipment health index proposed by the present invention can extract the deep-level fault features contained in the time-frequency domain for engineering equipment under the background of big data, and construct the equipment's health indicators on the premise of ensuring the original feature topology. health indicators. Compared with the prior art, the method of the invention fully extracts the deep-level fault features contained in the time-frequency domain information, and at the same time ensures the topology structure of the original features in the process of dimensionality reduction, further develops the fault feature extraction method, and enriches the equipment. Health indicators build a theoretical system, which effectively improves the performance characteristics of health indicators.
基于本发明提供的一种工程设备健康指标构建方法,本发明还提供一种工程设备健康指标构建系统,所述系统包括:Based on the construction method of engineering equipment health indicators provided by the present invention, the present invention also provides a construction system for engineering equipment health indicators, the system comprising:
时频域信息提取模块,用于获取工程设备的测试数据并计算所述测试数据的时频域信息;a time-frequency domain information extraction module, used for acquiring the test data of the engineering equipment and calculating the time-frequency domain information of the test data;
连续深度置信网络训练模块,用于将所述时频域信息输入至连续深度置信网络,借助于对比散度算法对所述连续深度置信网络进行训练,生成训练好的连续深度置信网络;The continuous depth belief network training module is used to input the time-frequency domain information into the continuous depth belief network, train the continuous depth belief network by means of a contrastive divergence algorithm, and generate a trained continuous depth belief network;
深层次故障特征提取模块,用于采用训练好的连续深度置信网络提取所述工程设备的深层次故障特征;The deep-level fault feature extraction module is used for extracting the deep-level fault features of the engineering equipment by using the trained continuous deep belief network;
SOM神经网络训练模块,用于将所述深层次故障特征作为SOM神经网络的输入进行无监督学习,生成训练好的SOM神经网络;The SOM neural network training module is used to perform unsupervised learning with the deep-level fault features as the input of the SOM neural network to generate a trained SOM neural network;
工程设备健康指标构建模块,用于通过计算所述训练好的SOM神经网络的输入层的状态向量与全部单元的权重向量间的欧氏距离,构建出所述工程设备的健康指标;所述健康指标直接反映所述工程设备的真实健康状态。The engineering equipment health index building module is used to construct the health index of the engineering equipment by calculating the Euclidean distance between the state vector of the input layer of the trained SOM neural network and the weight vectors of all units; The indicators directly reflect the real health status of the engineering equipment.
所述时频域信息提取模块具体包括:The time-frequency domain information extraction module specifically includes:
测试数据获取单元,用于获取工程设备的测试数据;The test data acquisition unit is used to acquire the test data of the engineering equipment;
时域特征提取单元,用于提取所述测试数据的时域特征;所述时域特征包括均值、均方根、方根幅值、平均绝对幅度、偏度、峭度、方差、最大值、最小值、峰值、波形因子、峰值因子、脉冲因子、裕度因子、偏度因子、峭度因子;A time-domain feature extraction unit for extracting time-domain features of the test data; the time-domain features include mean, root mean square, root square amplitude, mean absolute amplitude, skewness, kurtosis, variance, maximum value, Minimum value, peak value, shape factor, crest factor, pulse factor, margin factor, skewness factor, kurtosis factor;
频域特征提取单元,用于提取所述测试数据的频域特征;所述频域特征包括中心频率、均方根频率、频率方差;a frequency domain feature extraction unit, used for extracting the frequency domain feature of the test data; the frequency domain feature includes a center frequency, a root mean square frequency, and a frequency variance;
时频域信息构建单元,用于将提取出的所述时域特征和所述频域特征构成所述工程设备的时频域信息。The time-frequency domain information construction unit is configured to form the time-frequency domain information of the engineering equipment from the extracted time-domain features and the frequency-domain features.
所述连续深度置信网络训练模块具体包括:The continuous deep belief network training module specifically includes:
连续深度置信网络构建单元,用于构建由多个连续受限玻尔兹曼机堆叠而成的连续深度置信网络;所述连续受限玻尔兹曼机包含可视层、隐含层和高斯噪声单元;A continuous deep belief network building unit for constructing a continuous deep belief network composed of multiple continuous restricted Boltzmann machines stacked; the continuous restricted Boltzmann machine includes a visible layer, a hidden layer and a Gaussian noise unit;
连续深度置信网络训练单元,用于将所述时频域信息输入至所述连续深度置信网络,借助于对比散度算法更新所述连续深度置信网络的连接权重和sigmoid函数斜率控制参数,完成连续深度置信网络参数训练,生成训练好的连续深度置信网络。The continuous depth belief network training unit is used to input the time-frequency domain information into the continuous depth belief network, and update the connection weight and the sigmoid function slope control parameter of the continuous depth belief network by means of a contrastive divergence algorithm to complete the continuous depth belief network. Deep belief network parameter training to generate a trained continuous deep belief network.
所述深层次故障特征提取模块具体包括:The deep-level fault feature extraction module specifically includes:
深层次故障特征提取单元,用于将所述测试数据的时频域信息X=[x1,x2,…,xn]作为所述训练好的连续深度置信网络的输入,将所述训练好的连续深度置信网络的输出Y=[y1,y2,…,ym]作为所述工程设备的深层次故障特征;其中,xi(0≤i≤n)表示连续深度置信网络第i个输入神经元的特征序列,n为输入特征的总数;yi(0≤i≤m)表示连续深度置信网络第i个输出神经元的特征序列,m为连续深度置信网络输出神经元总数。A deep-level fault feature extraction unit, configured to use the time-frequency domain information X=[x 1 , x 2 , . . . , x n ] of the test data as the input of the trained continuous deep belief network, The output of a good continuous deep belief network Y=[y 1 , y 2 ,...,y m ] is used as the deep fault feature of the engineering equipment; wherein, x i (0≤i≤n) represents the first continuous deep belief network The feature sequence of i input neurons, n is the total number of input features; y i (0≤i≤m) represents the feature sequence of the ith output neuron of the continuous deep belief network, m is the total number of output neurons of the continuous deep belief network .
所述SOM神经网络训练模块具体包括:The SOM neural network training module specifically includes:
SOM神经网络训练单元,用于将所述深层次故障特征Y=[y1,y2,…,ym]中的一行作为所述SOM神经网络的输入向量D,根据所述输入向量D更新最佳匹配单元以及拓扑邻域的权重向量,完成所述SOM神经网络的无监督学习,生成训练好的SOM神经网络。The SOM neural network training unit is configured to use a row of the deep fault features Y=[y 1 , y 2 ,..., y m ] as the input vector D of the SOM neural network, and update according to the input vector D The weight vector of the best matching unit and topological neighborhood completes the unsupervised learning of the SOM neural network, and generates a trained SOM neural network.
本发明提供的一种新型工程设备健康指标构建方法及系统,充分提取了蕴含在时频域背后的深层次故障特征,同时确保了降维过程中原始特征的拓扑结构。首先,收集工程设备的测试数据并计算出相应的时频域信息;然后,将时频域信息输入至连续深度置信网络,借助于对比散度算法对网络模型进行训练,训练完成后可提取出工程设备的深层次故障特征,将该故障特征作为SOM神经网络的输入以实现其无监督学习,通过计算状态向量与全部单元的权重向量间的欧氏距离,定量刻画设备的健康状态。本发明所述方法进一步发展了故障特征提取方法,丰富了设备健康指标构建理论体系,有效改善了健康指标的性能特征,有助于后续的故障诊断与寿命预测,具有广阔的应用空间。The invention provides a new method and system for constructing health indicators of engineering equipment, which fully extracts the deep-level fault features contained in the time-frequency domain, and at the same time ensures the topology structure of the original features in the process of dimensionality reduction. First, the test data of the engineering equipment is collected and the corresponding time-frequency domain information is calculated; then, the time-frequency domain information is input into the continuous deep confidence network, and the network model is trained with the help of the contrastive divergence algorithm. The deep fault feature of engineering equipment is used as the input of SOM neural network to realize its unsupervised learning. By calculating the Euclidean distance between the state vector and the weight vector of all units, the health state of the equipment is quantitatively described. The method of the invention further develops the fault feature extraction method, enriches the theoretical system of equipment health index construction, effectively improves the performance characteristics of the health index, helps subsequent fault diagnosis and life prediction, and has broad application space.
对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。For the system disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant part can be referred to the description of the method. For the system disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant part can be referred to the description of the method.
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。In this paper, specific examples are used to illustrate the principles and implementations of the present invention. The descriptions of the above embodiments are only used to help understand the methods and core ideas of the present invention; meanwhile, for those skilled in the art, according to the present invention There will be changes in the specific implementation and application scope. In conclusion, the contents of this specification should not be construed as limiting the present invention.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210767222.0A CN115130514B (en) | 2022-06-30 | 2022-06-30 | A method and system for constructing engineering equipment health indicators |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210767222.0A CN115130514B (en) | 2022-06-30 | 2022-06-30 | A method and system for constructing engineering equipment health indicators |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115130514A true CN115130514A (en) | 2022-09-30 |
CN115130514B CN115130514B (en) | 2025-04-01 |
Family
ID=83382699
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210767222.0A Active CN115130514B (en) | 2022-06-30 | 2022-06-30 | A method and system for constructing engineering equipment health indicators |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115130514B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117196267A (en) * | 2023-11-07 | 2023-12-08 | 浙江正泰中自控制工程有限公司 | Intelligent inspection method and system based on equipment health analysis technology |
CN117407824A (en) * | 2023-12-14 | 2024-01-16 | 四川蜀能电科能源技术有限公司 | Health detection method, equipment and medium of power time synchronization device |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105973594A (en) * | 2016-04-25 | 2016-09-28 | 西北工业大学 | Rolling bearing fault prediction method based on continuous deep belief network |
WO2018076475A1 (en) * | 2016-10-26 | 2018-05-03 | 广东产品质量监督检验研究院 | Photovoltaic assembly accelerated degradation model established based on deep approach of learning, and method for predicting photovoltaic assembly lifetime |
CN109716346A (en) * | 2016-07-18 | 2019-05-03 | 河谷生物组学有限责任公司 | Distributed machines learning system, device and method |
CN111737911A (en) * | 2020-06-15 | 2020-10-02 | 大连理工大学 | A Bearing Recession Trend Prediction Method Based on Deep Belief Network and Extreme Learning Machine |
CN112613227A (en) * | 2020-12-15 | 2021-04-06 | 大连理工大学 | Model for predicting remaining service life of aero-engine based on hybrid machine learning |
CN112926273A (en) * | 2021-04-13 | 2021-06-08 | 中国人民解放军火箭军工程大学 | Method for predicting residual life of multivariate degradation equipment |
-
2022
- 2022-06-30 CN CN202210767222.0A patent/CN115130514B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105973594A (en) * | 2016-04-25 | 2016-09-28 | 西北工业大学 | Rolling bearing fault prediction method based on continuous deep belief network |
CN109716346A (en) * | 2016-07-18 | 2019-05-03 | 河谷生物组学有限责任公司 | Distributed machines learning system, device and method |
WO2018076475A1 (en) * | 2016-10-26 | 2018-05-03 | 广东产品质量监督检验研究院 | Photovoltaic assembly accelerated degradation model established based on deep approach of learning, and method for predicting photovoltaic assembly lifetime |
CN111737911A (en) * | 2020-06-15 | 2020-10-02 | 大连理工大学 | A Bearing Recession Trend Prediction Method Based on Deep Belief Network and Extreme Learning Machine |
CN112613227A (en) * | 2020-12-15 | 2021-04-06 | 大连理工大学 | Model for predicting remaining service life of aero-engine based on hybrid machine learning |
CN112926273A (en) * | 2021-04-13 | 2021-06-08 | 中国人民解放军火箭军工程大学 | Method for predicting residual life of multivariate degradation equipment |
Non-Patent Citations (6)
Title |
---|
YUBIN PAN等: "A hyb rid DBN-SOM-PF- based prognostic approach of remaining useful life for wind turbi ne gearbox", 《RENEWABLE ENERGY 》, 16 January 2020 (2020-01-16) * |
刘浩;熊炘;王小静;郭家宇;沈杰希: "基于自组织映射与受限玻尔兹曼机的滚动轴承健康评估", 机械传动, vol. 41, no. 6, 31 December 2017 (2017-12-31) * |
张全德;陈果;林桐;欧阳文理;滕春禹;王洪伟: "基于自组织神经网络的滚动轴承状态评估方法", 中国机械工程, vol. 28, no. 5, 31 March 2017 (2017-03-31) * |
赵光权;葛强强;刘小勇;彭喜元: "基于DBN的故障特征提取及诊断方法研究", 仪器仪表学报, vol. 37, no. 9, 30 September 2016 (2016-09-30) * |
邹涵: "基于深度学习的船舶导航系统健康管理", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》, 15 May 2021 (2021-05-15), pages 036 - 154 * |
郑晓飞;郭创;姚斌;冯华鑫: "基于深度学习的航空传感器故障诊断方法", 计算机工程, vol. 43, no. 7, 31 July 2017 (2017-07-31) * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117196267A (en) * | 2023-11-07 | 2023-12-08 | 浙江正泰中自控制工程有限公司 | Intelligent inspection method and system based on equipment health analysis technology |
CN117407824A (en) * | 2023-12-14 | 2024-01-16 | 四川蜀能电科能源技术有限公司 | Health detection method, equipment and medium of power time synchronization device |
CN117407824B (en) * | 2023-12-14 | 2024-02-27 | 四川蜀能电科能源技术有限公司 | Health detection method, equipment and medium of power time synchronization device |
Also Published As
Publication number | Publication date |
---|---|
CN115130514B (en) | 2025-04-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111307453B (en) | Transmission system fault diagnosis method based on multi-information fusion | |
Liang et al. | Unsupervised fault diagnosis of wind turbine bearing via a deep residual deformable convolution network based on subdomain adaptation under time-varying speeds | |
CN110516305B (en) | Intelligent fault diagnosis method under small sample based on attention mechanism meta-learning model | |
CN110491416A (en) | It is a kind of based on the call voice sentiment analysis of LSTM and SAE and recognition methods | |
CN118549823B (en) | Lithium battery electrical performance testing method and system | |
CN111461176A (en) | Multi-mode fusion method, device, medium and equipment based on normalized mutual information | |
CN112289391B (en) | Anode aluminum foil performance prediction system based on machine learning | |
CN113869208A (en) | Rolling bearing fault diagnosis method based on SA-ACWGAN-GP | |
CN104198184A (en) | Bearing fault diagnostic method based on second generation wavelet transform and BP neural network | |
CN110657984A (en) | A Fault Diagnosis Method of Planetary Gearbox Based on Enhanced Capsule Network | |
CN109948207A (en) | A method for predicting assembly errors of aero-engine high-pressure rotors | |
CN116502175A (en) | Method, device and storage medium for diagnosing fault of graph neural network | |
CN115130514A (en) | Construction method and system for health index of engineering equipment | |
CN113987910A (en) | Method and device for identifying load of residents by coupling neural network and dynamic time planning | |
Xu et al. | Zero-shot compound fault diagnosis method based on semantic learning and discriminative features | |
CN110991471B (en) | Fault diagnosis method for high-speed train traction system | |
CN116226646A (en) | Method, system, equipment and medium for predicting bearing health state and remaining life | |
CN117909881A (en) | Fault diagnosis method and device for multi-source data fusion pumping unit | |
CN114897138B (en) | System fault diagnosis method based on attention mechanism and deep residual network | |
CN112465054A (en) | Multivariate time series data classification method based on FCN | |
CN115048856A (en) | Method for predicting residual life of rolling bearing based on MS-ALSTM | |
Liang et al. | A novel intelligent fault diagnosis method for rolling bearings based on compressed sensing and stacked multi-granularity convolution denoising auto-encoder | |
CN115132324A (en) | Mental health prediction method and device, electronic device, storage medium | |
CN118690276A (en) | A high-speed rail wheel-rail sensing adversarial learning damage identification method | |
CN113052060A (en) | Bearing residual life prediction method and device based on data enhancement and electronic equipment |
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 |