Nothing Special   »   [go: up one dir, main page]

CN110633792B - End-to-end bearing health index construction method based on convolution cyclic neural network - Google Patents

End-to-end bearing health index construction method based on convolution cyclic neural network Download PDF

Info

Publication number
CN110633792B
CN110633792B CN201911004003.1A CN201911004003A CN110633792B CN 110633792 B CN110633792 B CN 110633792B CN 201911004003 A CN201911004003 A CN 201911004003A CN 110633792 B CN110633792 B CN 110633792B
Authority
CN
China
Prior art keywords
crnn
neural network
bearing
convolutional
cyclic neural
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201911004003.1A
Other languages
Chinese (zh)
Other versions
CN110633792A (en
Inventor
徐光华
陈龙庭
张四聪
况佳臣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong 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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201911004003.1A priority Critical patent/CN110633792B/en
Publication of CN110633792A publication Critical patent/CN110633792A/en
Application granted granted Critical
Publication of CN110633792B publication Critical patent/CN110633792B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/044Recurrent networks, e.g. Hopfield networks
    • 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
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/048Activation functions
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/084Backpropagation, e.g. using gradient descent

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Biomedical Technology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

一种端到端基于卷积循环神经网络的轴承健康指标构建方法,先搭建卷积循环神经网络模型CRNN,得到能衡量轴承退化程度的健康指标CRNN‑HI;然后进行卷积循环神经网络模型CRNN模型训练;最后进行卷积循环神经网络模型CRNN评估;本发明通过整合卷积神经网络CNNs和循环神经网络RNNs的结构优势,利用RNNs编码CNNs输出的特征图的时序信息,一方面消除了CNNs不能编码时序特征和感受野小的缺陷,另一方面也消除了RNNs不能自适应地从原始数据中提取HI相关特征的缺陷,使得CRNN‑HI在轴承健康状态退化评估方面取得了较高的相关性、单调性和准确性;同时本发明构建的CRNN‑HI指标将轴承的非线性退化过程近似表征为随时间变化的线性过程,为滚动轴承的健康状态评估和退化程度的确定提供了方便。

Figure 201911004003

An end-to-end method for building a bearing health index based on a convolutional cyclic neural network. First, a convolutional cyclic neural network model CRNN is built to obtain a health index CRNN‑HI that can measure the degree of bearing degradation; then the convolutional cyclic neural network model CRNN is obtained. Model training; finally, the convolutional cyclic neural network model CRNN evaluation is performed; the present invention integrates the structural advantages of the convolutional neural network CNNs and the cyclic neural network RNNs, and uses the RNNs to encode the time series information of the feature maps output by the CNNs. The defects of encoding temporal features and small receptive field, on the other hand, also eliminates the defect that RNNs cannot adaptively extract HI-related features from the original data, so that CRNN‑HI achieves high correlation in bearing health degradation assessment , monotonicity and accuracy; at the same time, the CRNN-HI index constructed by the present invention approximates the nonlinear degradation process of the bearing as a time-varying linear process, which provides convenience for the evaluation of the health state of the rolling bearing and the determination of the degradation degree.

Figure 201911004003

Description

端到端基于卷积循环神经网络的轴承健康指标构建方法End-to-end Convolutional Recurrent Neural Network-based Bearing Health Index Construction Method

技术领域technical field

本发明属于回转机械故障诊断技术领域,具体涉及一种端到端基于卷积循环神经网络的轴承健康指标构建方法。The invention belongs to the technical field of fault diagnosis of rotating machinery, and in particular relates to an end-to-end bearing health index construction method based on a convolutional cyclic neural network.

背景技术Background technique

近年来,随着物联网(IoT)技术和网络协同制造技术的发展,使得研究者能搜集到机械设备的大量状态监测数据,这些数据对于进一步提升智能制造水平显得十分重要。滚动轴承作为回转机械动力装备中关键部件之一,其意外的故障发生很可能会造成巨大的经济损失和灾难性事故。为了解决上述问题,对回转机械轴承部件实行计划性预测性维修显得十分重要,在上述大数据背景下,许多研究学者将注意力转向了基于数据驱动型的故障预示方法研究。In recent years, with the development of Internet of Things (IoT) technology and network collaborative manufacturing technology, researchers can collect a large amount of condition monitoring data of mechanical equipment, which is very important to further improve the level of intelligent manufacturing. As one of the key components in the power equipment of rotating machinery, the unexpected failure of rolling bearing is likely to cause huge economic losses and catastrophic accidents. In order to solve the above problems, it is very important to implement planned predictive maintenance for rotary mechanical bearing components. Under the background of the above big data, many researchers have turned their attention to the research on data-driven failure prediction methods.

总的来说,数据驱动型故障预示方法主要包含3个步骤:(1)故障相关数据采集;(2)健康指标HI构建;(3)失效时间预测。其中,健康指标HI构建的好坏对于定量评估轴承退化程度的准确性产生至关重要的影响。近几年以大数据分析为中心的人工智能技术在轴承健康指标HI构建领域取得了很大的进步,其典型代表就是基于深度学习DL的方法。基于DL的HI构建方法主要分成两类,第一类是基于卷积神经网络CNNs的方法,其核心是利用CNNs的分层特征学习和分层特征表达能力,自适应地从原始数据如振动信号中学习HI相关的特征,通过网络的权重共享和局部感受野属性,浅层网络能学习到局部特征而高层网络能学习到更加全局和具有丰富语义信息的特征;有人利用一维卷积神经网络为滚动轴承构建了HI指标,并将HI中的毛刺波动现象考虑进去,在相关性和单调性方面取得了很大的性能提升。然而该方法的缺点在于,一方面直接采用分类的CNNs网络结构用于HI的回归,其回归能力比较差;另一方面,CNNs不能编码输入数据的时间序列信息,而该信息对于轴承HI的构建是很重要的。第二类是基于循环神经网络RNNs的方法,其核心是利用RNNs网络能学习长期依赖关系的记忆性去编码输入特征的时序信息,有人通过给RNNs网络输入人工自定义的特征,包括经典的时域特征、频域特征,以及时频域特征,为滚动轴承构建了健康指标RNN-HI。RNN-HI在衡量滚动轴承性能退化方面具有较高的单调性。然而,该类方法的缺陷在于需要人工自定义特征输入到网络中,即需要一定的先验知识,且整个网络不能进行端到端的优化。In general, the data-driven fault prediction method mainly includes three steps: (1) fault-related data collection; (2) health index HI construction; (3) failure time prediction. Among them, the construction of health index HI has a crucial impact on the accuracy of quantitatively evaluating the degree of bearing degradation. In recent years, artificial intelligence technology centered on big data analysis has made great progress in the field of bearing health index HI construction, and its typical representative is the method based on deep learning DL. DL-based HI construction methods are mainly divided into two categories. The first category is the method based on convolutional neural network CNNs. Through the weight sharing and local receptive field attributes of the network, the shallow network can learn local features and the high-level network can learn more global and semantically rich features; some people use one-dimensional convolutional neural networks. The HI index is constructed for rolling bearings, and the burr fluctuation phenomenon in HI is taken into account, resulting in a great performance improvement in terms of correlation and monotonicity. However, the disadvantage of this method is that, on the one hand, the classified CNNs network structure is directly used for the regression of HI, and its regression ability is relatively poor; is very important. The second type is the method based on cyclic neural network RNNs, the core of which is to use the memory of RNNs network to learn long-term dependencies to encode the timing information of input features. Some people input artificially customized features to the RNNs network, including classic time Domain features, frequency domain features, and time-frequency domain features construct a health index RNN-HI for rolling bearings. RNN-HI has high monotonicity in measuring the performance degradation of rolling bearings. However, the disadvantage of this type of method is that it needs to manually define features to be input into the network, that is, certain prior knowledge is required, and the entire network cannot be optimized end-to-end.

发明内容SUMMARY OF THE INVENTION

为了克服上述现有技术的缺点,本发明的目的在于提供了一种端到端基于卷积循环神经网络的轴承健康指标构建方法,通过整合卷积神经网络和循环神经网络的结构优势,利用RNNs编码CNNs输出的特征图的时序信息,一方面消除了CNNs不能编码时序特征和感受野小的缺陷,另一方面也消除了RNNs不能自适应地从原始数据中提取HI相关特征的缺陷。In order to overcome the above-mentioned shortcomings of the prior art, the purpose of the present invention is to provide an end-to-end method for constructing a bearing health index based on a convolutional cyclic neural network. Encoding the temporal information of the feature maps output by CNNs eliminates the defect that CNNs cannot encode temporal features and small receptive fields on the one hand, and also eliminates the defect that RNNs cannot adaptively extract HI-related features from the original data.

为了达到上述目的,本发明采取的技术方案是:In order to achieve the above object, the technical scheme adopted by the present invention is:

一种端到端基于卷积循环神经网络的轴承健康指标构建方法,包括以下步骤:An end-to-end convolutional recurrent neural network based bearing health index construction method, including the following steps:

1)搭建卷积循环神经网络模型CRNN:卷积循环神经网络模型CRNN的输入为被标准化的原始轴承振动信号数据,卷积循环神经网络模型CRNN包含两个组件,一个组件是用来提取局部特征的卷积神经网络CNNs,由交替堆叠的卷积层和池化层组成,为了加速网络的训练过程和收敛速度,在卷积神经网络CNNs中加入了BatchNorm层和残差连接结构;卷积神经网络CNNs包含4组残差单元,每一个残差单元包含有2个卷积层且每个残差单元对输入的数据均进行降采样,卷积神经网络CNNs的第一个卷积层采用预激活设计,整个卷积层中卷积核的大小均为16;另外一个组件用来提取带有时间序列信息的全局特征,并循环连接上述局部特征的循环神经网络,采用长短期记忆网络LSTM实现,长短期记忆网络LSTM在每一个时刻的输入为卷积神经网络CNNs输出的特征图上的每一个特征点,时间长度等于特征图的维度;长短期记忆网络LSTM在t时刻的输出通过以下公式计算:1) Build a convolutional cyclic neural network model CRNN: The input of the convolutional cyclic neural network model CRNN is the standardized original bearing vibration signal data. The convolutional cyclic neural network model CRNN contains two components, one component is used to extract local features The convolutional neural network CNNs consist of alternately stacked convolutional layers and pooling layers. In order to speed up the training process and convergence speed of the network, a BatchNorm layer and a residual connection structure are added to the convolutional neural network CNNs; convolutional neural network Network CNNs contain 4 groups of residual units, each residual unit contains 2 convolutional layers and each residual unit downsamples the input data, the first convolutional layer of convolutional neural network CNNs adopts pre- Activation design, the size of the convolution kernel in the entire convolution layer is 16; another component is used to extract global features with time series information, and cyclically connect the cyclic neural network of the above local features, which is implemented by long short-term memory network LSTM , the input of the long short-term memory network LSTM at each time is each feature point on the feature map output by the convolutional neural network CNNs, and the time length is equal to the dimension of the feature map; the output of the long short-term memory network LSTM at time t passes the following formula calculate:

ft=σ(Wf·[ht-1,xt]+bf) (1)f t =σ(W f ·[h t-1 ,x t ]+b f ) (1)

Figure BDA0002242201670000031
Figure BDA0002242201670000031

it=σ(Wi·[ht-1,xt]+bi) (3)i t =σ(W i ·[h t-1 ,x t ]+b i ) (3)

Figure BDA0002242201670000032
Figure BDA0002242201670000032

ot=σ(Wo·[ht-1,xt]+bo) (5)o t =σ(W o ·[h t-1 ,x t ]+b o ) (5)

ht=ot⊙tanh(Ct) (6)h t =o t ⊙tanh(C t ) (6)

其中,ht-1和xt分别代表长短期记忆网络LSTM在t-1时刻的隐层状态和t时刻的输入特征;Wf,Wi和Wo分别表示长短期记忆网络LSTM遗忘门、输入门和输出门的权重参数;bf,bi和bo分别表示长短期记忆网络LSTM遗忘门、输入门和输出门的偏置参数;ft,it和ot分别表示长短期记忆网络LSTM遗忘门、输入门和输出门的激活值;激活函数σ为sigmoid函数;符号tanh表示双曲正切函数;符号⊙和·分别表示逐元素乘法和矩阵乘法;

Figure BDA0002242201670000033
表示长短期记忆网络LSTM新学习到的特征信息,WC和bC分别为相对应的权重参数和偏置参数;Ct表示当前t时刻长短期记忆网络LSTM细胞cell的内部状态;当前时刻隐藏层状态ht包含了带时间序列信息的全局特征,随后将其通过一个全连接层的非线性映射,从而得到能衡量轴承退化程度的健康指标CRNN-HI;Among them, h t-1 and x t respectively represent the hidden layer state of the long short-term memory network LSTM at time t-1 and the input features at time t; W f , Wi and Wo represent the long short-term memory network LSTM forgetting gate, The weight parameters of the input gate and the output gate; b f , b i and b o represent the bias parameters of the long short-term memory network LSTM forgetting gate, input gate and output gate, respectively; f t , i t and o t represent the long short-term memory network, respectively The activation values of the network LSTM forget gate, input gate and output gate; the activation function σ is the sigmoid function; the symbol tanh represents the hyperbolic tangent function; the symbols ⊙ and · represent element-wise multiplication and matrix multiplication, respectively;
Figure BDA0002242201670000033
Represents the newly learned feature information of the long short-term memory network LSTM, W C and b C are the corresponding weight parameters and bias parameters respectively; C t represents the internal state of the long short-term memory network LSTM cell at the current time t; hidden at the current time The layer state h t contains global features with time series information, and then passes it through a nonlinear mapping of a fully connected layer to obtain a health index CRNN-HI that can measure the degree of bearing degradation;

2)卷积循环神经网络模型CRNN模型训练:卷积循环神经网络模型CRNN的损失函数采用均方根回归误差,训练集χ由滚动轴承的全寿命周期振动信号组成,χ={(xt,yt)|1≤t≤T},其中xt∈RN,yt∈[0,1],N表示采集到的每一个振动信号的数据样本长度,T表示滚动轴承的寿命,yt是一个实数值,代表在当前时刻t下轴承的残余寿命百分比;卷积循环神经网络模型CRNN的损失函数J通过下式计算:2) Convolutional cyclic neural network model CRNN model training: The loss function of the convolutional cyclic neural network model CRNN adopts the root mean square regression error, and the training set χ consists of the vibration signals of the rolling bearing throughout its life cycle, χ={(x t , y t )|1≤t≤T}, where x t ∈ R N , y t ∈ [0,1], N represents the data sample length of each vibration signal collected, T represents the life of the rolling bearing, y t is a Real value, representing the residual life percentage of the bearing at the current time t; the loss function J of the convolutional recurrent neural network model CRNN is calculated by the following formula:

Figure BDA0002242201670000041
Figure BDA0002242201670000041

其中,

Figure BDA0002242201670000042
是HI的预测值,当某个轴承的
Figure BDA0002242201670000043
值十分接近0的时候,就表示这个轴承处于失效状态;整个卷积循环神经网络模型CRNN采用反向传播算法BP进行训练,优化器选择为Adam,初始学习率lr设置为0.001,每个batch的样本量大小设置为32;in,
Figure BDA0002242201670000042
is the predicted value of HI, when the
Figure BDA0002242201670000043
When the value is very close to 0, it means that the bearing is in a failed state; the entire convolutional cyclic neural network model CRNN is trained by the back-propagation algorithm BP, the optimizer is selected as Adam, the initial learning rate lr is set to 0.001, and each batch of The sample size is set to 32;

3)卷积循环神经网络模型CRNN评估,挑选最优的训练过的卷积循环神经网络模型CRNN去预测待评估的滚动轴承的健康指标CRNN-HI,这里采用相关性指标Corr和单调性指标Mono定量评估训练过程中卷积循环神经网络模型CRNN在验证集上的性能,相关性指标Corr衡量了CRNN-HI衰减趋势与设备运行时间之间的相关性,其值计算如下:3) Convolutional cyclic neural network model CRNN evaluation, select the optimal trained convolutional cyclic neural network model CRNN to predict the health index CRNN-HI of the rolling bearing to be evaluated, here the correlation index Corr and the monotonicity index Mono are used to quantify To evaluate the performance of the convolutional recurrent neural network model CRNN on the validation set during the training process, the correlation index Corr measures the correlation between the CRNN-HI decay trend and the running time of the device, and its value is calculated as follows:

Figure BDA0002242201670000044
Figure BDA0002242201670000044

符号Ii和ti分别表示CRNN-HI的值和运行时刻,

Figure BDA0002242201670000045
单调性指标Mono衡量了CRNN-HI随时间递增或递减的趋势,其计算公式如下:The symbols I i and t i represent the value and running time of CRNN-HI, respectively,
Figure BDA0002242201670000045
The monotonicity indicator Mono measures the increasing or decreasing trend of CRNN-HI over time, and its calculation formula is as follows:

Figure BDA0002242201670000051
Figure BDA0002242201670000051

dI表示相邻时刻两个CRNN-HI之间的一阶差分,看出Corr和Mono的值的取值范围均在0到1之间,当最优模型参数选择完毕之后,就将采集到的待分析轴承的振动信号输入到最优模型中,从而得到最优模型的输出,及轴承当前时刻的CRNN-HI值,从而定量评估该滚动轴承的当前退化程度。dI represents the first-order difference between two CRNN-HIs at adjacent times. It can be seen that the values of Corr and Mono range from 0 to 1. When the optimal model parameters are selected, the collected The vibration signal of the bearing to be analyzed is input into the optimal model, so as to obtain the output of the optimal model and the CRNN-HI value of the bearing at the current moment, so as to quantitatively evaluate the current degradation degree of the rolling bearing.

本发明的有益效果为:The beneficial effects of the present invention are:

本发明将卷积神经网络和循环神经网络的结构优势进行结合,构建了端到端的基于卷积循环神经网络的轴承健康指标方法,克服了基于CNNs模型和基于RNNs模型的HI构建方法的缺陷,使得CRNN-HI在轴承健康状态退化评估方面取得了较高的相关性、单调性和准确性。同时,本发明构建的CRNN-HI指标将轴承的非线性退化过程近似表征为随时间变化的线性过程,为滚动轴承的健康状态评估和退化程度的确定提供了方便。The invention combines the structural advantages of the convolutional neural network and the cyclic neural network, constructs an end-to-end bearing health index method based on the convolutional cyclic neural network, and overcomes the defects of the HI construction method based on the CNNs model and the RNNs model. As a result, CRNN-HI achieves high correlation, monotonicity and accuracy in bearing health state degradation assessment. At the same time, the CRNN-HI index constructed by the present invention approximates the nonlinear degradation process of the bearing as a linear process that changes with time, which provides convenience for evaluating the health state of the rolling bearing and determining the degree of degradation.

附图说明Description of drawings

图1为本发明提出的卷积循环神经网络模型图。FIG. 1 is a diagram of a convolutional cyclic neural network model proposed by the present invention.

图2为本发明LSTM神经网络模型的结构示意图。FIG. 2 is a schematic structural diagram of the LSTM neural network model of the present invention.

图3为本发明实施例的应用效果图。FIG. 3 is an application effect diagram of an embodiment of the present invention.

具体实施方式Detailed ways

下面结合附图和实施例对本发明作详细描述。The present invention will be described in detail below with reference to the accompanying drawings and embodiments.

参照图1,一种端到端基于卷积循环神经网络的轴承健康指标构建方法,包括以下步骤:Referring to Figure 1, an end-to-end method for constructing a bearing health index based on a convolutional recurrent neural network includes the following steps:

1)搭建卷积循环神经网络模型CRNN:卷积循环神经网络模型CRNN的输入为被标准化的原始轴承振动信号数据,卷积循环神经网络模型CRNN包含两个组件,一个组件是用来提取局部特征的卷积神经网络CNNs,由交替堆叠的卷积层和池化层组成,为了加速网络的训练过程和收敛速度,在卷积神经网络CNNs中加入了BatchNorm层和残差连接结构;卷积神经网络CNNs包含4组残差单元,每一个残差单元包含有2个卷积层且每个残差单元对输入的数据均进行降采样,卷积神经网络CNNs的第一个卷积层采用预激活设计,整个卷积层中卷积核的大小均为16;另外一个组件用来提取带有时间序列信息的全局特征,并循环连接上述局部特征的循环神经网络,采用长短期记忆网络LSTM实现,参照图2。长短期记忆网络LSTM在每一个时刻的输入为卷积神经网络CNNs输出的特征图上的每一个特征点,时间长度等于特征图的维度;长短期记忆网络LSTM在t时刻的输出通过以下公式计算:1) Build a convolutional cyclic neural network model CRNN: The input of the convolutional cyclic neural network model CRNN is the standardized original bearing vibration signal data. The convolutional cyclic neural network model CRNN contains two components, one component is used to extract local features The convolutional neural network CNNs consist of alternately stacked convolutional layers and pooling layers. In order to speed up the training process and convergence speed of the network, a BatchNorm layer and a residual connection structure are added to the convolutional neural network CNNs; convolutional neural network Network CNNs contain 4 groups of residual units, each residual unit contains 2 convolutional layers and each residual unit downsamples the input data, the first convolutional layer of convolutional neural network CNNs adopts pre- Activation design, the size of the convolution kernel in the entire convolution layer is 16; another component is used to extract global features with time series information, and cyclically connect the cyclic neural network of the above local features, which is implemented by long short-term memory network LSTM , refer to Figure 2. The input of the long short-term memory network LSTM at each moment is each feature point on the feature map output by the convolutional neural network CNNs, and the time length is equal to the dimension of the feature map; the output of the long short-term memory network LSTM at time t is calculated by the following formula :

ft=σ(Wf·[ht-1,xt]+bf) (1)f t =σ(W f ·[h t-1 ,x t ]+b f ) (1)

Figure BDA0002242201670000061
Figure BDA0002242201670000061

it=σ(Wi·[ht-1,xt]+bi) (3)i t =σ(W i ·[h t-1 ,x t ]+b i ) (3)

Figure BDA0002242201670000062
Figure BDA0002242201670000062

ot=σ(Wo·[ht-1,xt]+bo) (5)o t =σ(W o ·[h t-1 ,x t ]+b o ) (5)

ht=ot⊙tanh(Ct) (6)h t =o t ⊙tanh(C t ) (6)

其中,ht-1和xt分别代表长短期记忆网络LSTM在t-1时刻的隐层状态和t时刻的输入特征;Wf,Wi和Wo分别表示长短期记忆网络LSTM遗忘门、输入门和输出门的权重参数;bf,bi和bo分别表示长短期记忆网络LSTM遗忘门、输入门和输出门的偏置参数;ft,it和ot分别表示长短期记忆网络LSTM遗忘门、输入门和输出门的激活值;激活函数σ为sigmoid函数;符号tanh表示双曲正切函数;符号⊙和·分别表示逐元素乘法和矩阵乘法;

Figure BDA0002242201670000063
表示长短期记忆网络LSTM新学习到的特征信息,WC和bC分别为相对应的权重参数和偏置参数;Ct表示当前t时刻长短期记忆网络LSTM细胞cell的内部状态;当前时刻隐藏层状态ht包含了带时间序列信息的全局特征,随后将其通过一个全连接层的非线性映射,从而得到能衡量轴承退化程度的健康指标CRNN-HI;Among them, h t-1 and x t respectively represent the hidden layer state of the long short-term memory network LSTM at time t-1 and the input features at time t; W f , Wi and Wo represent the long short-term memory network LSTM forgetting gate, The weight parameters of the input gate and the output gate; b f , b i and b o represent the bias parameters of the long short-term memory network LSTM forgetting gate, input gate and output gate, respectively; f t , i t and o t represent the long short-term memory network, respectively The activation values of the network LSTM forget gate, input gate and output gate; the activation function σ is the sigmoid function; the symbol tanh represents the hyperbolic tangent function; the symbols ⊙ and · represent element-wise multiplication and matrix multiplication, respectively;
Figure BDA0002242201670000063
Represents the newly learned feature information of the long short-term memory network LSTM, W C and b C are the corresponding weight parameters and bias parameters respectively; C t represents the internal state of the long short-term memory network LSTM cell at the current time t; hidden at the current time The layer state h t contains global features with time series information, and then passes it through a nonlinear mapping of a fully connected layer to obtain a health index CRNN-HI that can measure the degree of bearing degradation;

2)卷积循环神经网络模型CRNN模型训练:卷积循环神经网络模型CRNN的损失函数采用均方根回归误差,训练集χ由滚动轴承的全寿命周期振动信号组成,χ={(xt,yt)|1≤t≤T},其中xt∈RN,yt∈[0,1],N表示采集到的每一个振动信号的数据样本长度,T表示滚动轴承的寿命,yt是一个实数值,代表在当前时刻t下轴承的残余寿命百分比;卷积循环神经网络模型CRNN的损失函数J通过下式计算:2) Convolutional cyclic neural network model CRNN model training: The loss function of the convolutional cyclic neural network model CRNN adopts the root mean square regression error, and the training set χ consists of the vibration signals of the rolling bearing throughout its life cycle, χ={(x t , y t )|1≤t≤T}, where x t ∈ R N , y t ∈ [0,1], N represents the data sample length of each vibration signal collected, T represents the life of the rolling bearing, y t is a Real value, representing the residual life percentage of the bearing at the current time t; the loss function J of the convolutional recurrent neural network model CRNN is calculated by the following formula:

Figure BDA0002242201670000071
Figure BDA0002242201670000071

其中,

Figure BDA0002242201670000072
是HI的预测值,当某个轴承的
Figure BDA0002242201670000073
值十分接近0的时候,就表示这个轴承处于失效状态;整个卷积循环神经网络模型CRNN采用反向传播算法BP进行训练,优化器选择为Adam,初始学习率lr设置为0.001,每个batch的样本量大小设置为32;in,
Figure BDA0002242201670000072
is the predicted value of HI, when the
Figure BDA0002242201670000073
When the value is very close to 0, it means that the bearing is in a failed state; the entire convolutional cyclic neural network model CRNN is trained by the back propagation algorithm BP, the optimizer is selected as Adam, the initial learning rate lr is set to 0.001, and the The sample size is set to 32;

3)卷积循环神经网络模型CRNN评估,挑选最优的训练过的卷积循环神经网络模型CRNN去预测待评估的滚动轴承的健康指标CRNN-HI,这里采用相关性指标Corr和单调性指标Mono定量评估训练过程中卷积循环神经网络模型CRNN在验证集上的性能,相关性指标Corr衡量了CRNN-HI衰减趋势与设备运行时间之间的相关性,其值计算如下:3) Convolutional cyclic neural network model CRNN evaluation, select the optimal trained convolutional cyclic neural network model CRNN to predict the health index CRNN-HI of the rolling bearing to be evaluated, here the correlation index Corr and the monotonicity index Mono are used to quantify To evaluate the performance of the convolutional recurrent neural network model CRNN on the validation set during the training process, the correlation index Corr measures the correlation between the CRNN-HI decay trend and the running time of the device, and its value is calculated as follows:

Figure BDA0002242201670000074
Figure BDA0002242201670000074

符号Ii和ti分别表示CRNN-HI的值和运行时刻,

Figure BDA0002242201670000081
单调性指标Mono衡量了CRNN-HI随时间递增或递减的趋势,其计算公式如下:The symbols I i and t i represent the value and running time of CRNN-HI, respectively,
Figure BDA0002242201670000081
The monotonicity indicator Mono measures the increasing or decreasing trend of CRNN-HI over time, and its calculation formula is as follows:

Figure BDA0002242201670000082
Figure BDA0002242201670000082

dI表示相邻时刻两个CRNN-HI之间的一阶差分,看出Corr和Mono的值的取值范围均在0到1之间,当最优模型参数选择完毕之后,就将采集到的待分析轴承的振动信号输入到最优模型中,从而得到最优模型的输出,及轴承当前时刻的CRNN-HI值,从而定量评估该滚动轴承的当前退化程度。dI represents the first-order difference between two CRNN-HIs at adjacent times. It can be seen that the values of Corr and Mono range from 0 to 1. When the optimal model parameters are selected, the collected The vibration signal of the bearing to be analyzed is input into the optimal model, so as to obtain the output of the optimal model and the CRNN-HI value of the bearing at the current moment, so as to quantitatively evaluate the current degradation degree of the rolling bearing.

本发明已经成功应用在滚动轴承的健康状态评估中,参照图3,该图是本发明在两种运行工况条件下对4个轴承全寿命周期振动数据的健康状态的预测结果,横坐标是运行时间,纵坐标是预测的CRNN-HI的值,由图3可见,CRNN-HI在4个轴承的退化趋势预测上均表现出了较高的相关性和单调性,且将滚动轴承的非线性退化过程表征为了随时间变化的线性过程,即使轴承工作在不同的运行工况条件下。The present invention has been successfully applied in the health state assessment of rolling bearings. Referring to FIG. 3 , the figure is the prediction result of the present invention on the health state of the vibration data of four bearings in the whole life cycle under two operating conditions. The abscissa is the operation Time, the ordinate is the predicted value of CRNN-HI. It can be seen from Figure 3 that CRNN-HI shows high correlation and monotonicity in the prediction of the degradation trend of the four bearings, and degrades the nonlinearity of rolling bearings. The process characterization is a linear process over time, even if the bearing operates under different operating conditions.

以上实例只为说明本发明的技术构思和特点,并不能以此限制本发明的保护范围。对于本领域的技术人员来说,凡是根据本发明精神实质所做的等效变换或修饰改进,都应涵盖在本发明的保护范围之内。The above examples are only intended to illustrate the technical concept and characteristics of the present invention, and do not limit the protection scope of the present invention. For those skilled in the art, all equivalent transformations or modifications made according to the spirit of the present invention should be covered within the protection scope of the present invention.

Claims (1)

1. An end-to-end convolution cyclic neural network-based bearing health index construction method is characterized by comprising the following steps:
1) constructing a convolution cyclic neural network model CRNN: the input of the convolution cyclic neural network model CRNN is standardized original bearing vibration signal data, the convolution cyclic neural network model CRNN comprises two components, one component is a convolution neural network CNNs used for extracting local characteristics and consists of convolution layers and pooling layers which are stacked alternately, and in order to accelerate the training process and the convergence speed of the network, a BatchNorm layer and a residual error connecting structure are added into the convolution neural network CNNs; the convolutional neural networks CNNs comprise 4 groups of residual error units, each residual error unit comprises 2 convolutional layers, each residual error unit performs down-sampling on input data, the first convolutional layer of the convolutional neural networks CNNs adopts a pre-activation design, and the sizes of convolutional cores in the whole convolutional layer are all 16; the other component is used for extracting global characteristics with time series information and circularly connecting the cyclic neural network of the local characteristics, and is realized by adopting a long-short term memory network LSTM, the input of the long-short term memory network LSTM at each moment is each characteristic point on a characteristic diagram output by the convolutional neural network CNNs, and the time length is equal to the dimension of the characteristic diagram; the output of the long-short term memory network LSTM at the time t is calculated by the following formula:
ft=σ(Wf·[ht-1,xt]+bf) (1)
Figure FDA0002242201660000011
it=σ(Wi·[ht-1,xt]+bi) (3)
Figure FDA0002242201660000012
ot=σ(Wo·[ht-1,xt]+bo) (5)
ht=ot⊙tanh(Ct) (6)
wherein h ist-1And xtRespectively representing the hidden layer state of the long-short term memory network LSTM at the time t-1 and the input characteristics at the time t; wf,WiAnd WoRespectively representing the weight parameters of a long-short term memory network LSTM forgetting gate, an input gate and an output gate; bf,biAnd boRespectively representing the bias parameters of a long-short term memory network LSTM forgetting gate, an input gate and an output gate; f. oft,itAnd otRespectively representing the activation values of a long-short term memory network LSTM forgetting gate, an input gate and an output gate; the activation function sigma is a sigmoid function; the notation tanh denotes a hyperbolic tangent function; the symbols & · denote element-by-element multiplication and matrix multiplication, respectively;
Figure FDA0002242201660000021
representing newly learned feature information, W, of long-and short-term memory networks LSTMCAnd bCRespectively corresponding weight parameters and bias parameters; ctRepresenting the internal state of the LSTM cell at the current time t; hidden layer state h at the current momenttThe method comprises the steps that global characteristics with time series information are included, and then the global characteristics are subjected to nonlinear mapping through a full connection layer, so that a health index CRNN-HI capable of measuring the degradation degree of a bearing is obtained;
2) training a convolution cyclic neural network model (CRNN): the loss function of the convolution cyclic neural network model CRNN adopts a root mean square regression error, a training set chi is composed of a full life cycle vibration signal of a rolling bearing, and chi { (x)t,yt) L 1 is less than or equal to T is less than or equal to T, wherein xt∈RN,yt∈[0,1]N represents the length of a data sample of each vibration signal acquired, T represents the life of the rolling bearing, ytIs a real number representing the percentage of residual life of the bearing at the current time t; the loss function J of the convolutional recurrent neural network model CRNN is calculated by:
Figure FDA0002242201660000022
wherein,
Figure FDA0002242201660000023
is the predicted value of HI for a bearing
Figure FDA0002242201660000024
A value very close to 0 indicates that the bearing is in a failure state; the whole convolution cyclic neural network model CRNN is trained by adopting a back propagation algorithm BP, an optimizer selects Adam, the initial learning rate lr is set to be 0.001, and the sample size of each batch is set to be 32;
3) evaluating a convolution cycle neural network model (CRNN), selecting an optimal trained convolution cycle neural network model (CRNN) to predict a health index (CRNN-HI) of a rolling bearing to be evaluated, wherein a correlation index (Corr) and a monotonicity index (Mono) are adopted to quantitatively evaluate the performance of the convolution cycle neural network model (CRNN) on a verification set in a training process, the correlation index (Corr) measures the correlation between the attenuation trend of the CRNN-HI and the equipment running time, and the value calculation is as follows:
Figure FDA0002242201660000031
symbol IiAnd tiRespectively representing the value of CRNN-HI and the running time,
Figure FDA0002242201660000032
the monotonicity index Mono measures the increasing or decreasing trend of the CRNN-HI along with time, and the calculation formula is as follows:
Figure FDA0002242201660000033
and dI represents a first-order difference between two CRNN-HI at adjacent moments, the value ranges of the values of Corr and Mono are both between 0 and 1, and after the optimal model parameters are selected, the acquired vibration signals of the bearing to be analyzed are input into the optimal model, so that the output of the optimal model and the CRNN-HI value of the bearing at the current moment are obtained, and the current degradation degree of the rolling bearing is quantitatively evaluated.
CN201911004003.1A 2019-10-22 2019-10-22 End-to-end bearing health index construction method based on convolution cyclic neural network Active CN110633792B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911004003.1A CN110633792B (en) 2019-10-22 2019-10-22 End-to-end bearing health index construction method based on convolution cyclic neural network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911004003.1A CN110633792B (en) 2019-10-22 2019-10-22 End-to-end bearing health index construction method based on convolution cyclic neural network

Publications (2)

Publication Number Publication Date
CN110633792A CN110633792A (en) 2019-12-31
CN110633792B true CN110633792B (en) 2022-03-22

Family

ID=68976959

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911004003.1A Active CN110633792B (en) 2019-10-22 2019-10-22 End-to-end bearing health index construction method based on convolution cyclic neural network

Country Status (1)

Country Link
CN (1) CN110633792B (en)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111291673A (en) * 2020-02-04 2020-06-16 贵州省人民医院 A method, device, server and storage medium for predicting classification of time series signals
CN111832216B (en) * 2020-04-14 2022-07-12 新疆大学 Prediction method of remaining service life of rolling bearing based on EEMD-MCNN-GRU
CN111680661A (en) * 2020-06-19 2020-09-18 哈尔滨工业大学 A performance degradation tracking method for mechanical rotating parts based on multi-feature fusion
CN111856287B (en) * 2020-07-17 2021-07-13 上海交通大学 Lithium battery health state detection method based on stacked residual causal convolutional neural network
CN112033463B (en) * 2020-09-02 2022-09-06 哈尔滨工程大学 Nuclear power equipment state evaluation and prediction integrated method and system
CN112347898B (en) * 2020-11-03 2024-04-09 重庆大学 Rolling bearing health index construction method based on DCAE neural network
CN112729825A (en) * 2020-12-18 2021-04-30 中车永济电机有限公司 Method for constructing bearing fault diagnosis model based on convolution cyclic neural network
CN112992119B (en) * 2021-01-14 2024-05-03 安徽大学 Accent classification method based on deep neural network and model thereof
CN112766182B (en) * 2021-01-22 2024-02-02 浙江工业大学 Broaching tool wear state identification method
CN113283155B (en) * 2021-03-22 2024-04-09 中国农业科学院农业资源与农业区划研究所 Near-surface air temperature estimation method, system, storage medium and equipment
CN113553903B (en) * 2021-06-11 2022-11-22 东南大学 A Deep Temporal Convolutional Network for Health Assessment of Rotating Machinery
CN113505535B (en) * 2021-07-08 2023-11-10 重庆大学 Motor life prediction method based on gating self-adaptive hierarchical attention unit network
EP4124529B1 (en) 2021-07-29 2023-08-23 Volvo Truck Corporation Method for detecting a failure in a vehicle pneumatic system
CN114036691B (en) * 2021-11-26 2024-06-21 江苏科技大学 Bearing degradation evaluation method and system based on recursion complex network
CN114925723B (en) * 2022-05-06 2023-04-07 盐城工学院 Method for predicting residual service life of rolling bearing by adopting encoder and decoder
CN115033984B (en) * 2022-05-26 2024-09-24 西北工业大学 LSTM-based comprehensive modularized avionics system health assessment method
WO2024108577A1 (en) * 2022-11-25 2024-05-30 深圳先进技术研究院 Method for interpreting pancreatitis diagnosis model by using layer-wise relevance propagation

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106934126A (en) * 2017-02-28 2017-07-07 西安交通大学 Component of machine health indicator building method based on Recognition with Recurrent Neural Network fusion
CN108960077A (en) * 2018-06-12 2018-12-07 南京航空航天大学 A kind of intelligent failure diagnosis method based on Recognition with Recurrent Neural Network
CN109522600A (en) * 2018-10-16 2019-03-26 浙江大学 Complex equipment remaining life prediction technique based on combined depth neural network
CN109726524A (en) * 2019-03-01 2019-05-07 哈尔滨理工大学 A prediction method for the remaining service life of rolling bearings based on CNN and LSTM
CN110108456A (en) * 2019-04-16 2019-08-09 东南大学 A kind of rotating machinery health evaluating method of depth convolutional neural networks

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10410116B2 (en) * 2014-03-11 2019-09-10 SparkCognition, Inc. System and method for calculating remaining useful time of objects
WO2018226437A1 (en) * 2017-06-05 2018-12-13 Adasky, Ltd. Shutterless far infrared (fir) camera for automotive safety and driving systems

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106934126A (en) * 2017-02-28 2017-07-07 西安交通大学 Component of machine health indicator building method based on Recognition with Recurrent Neural Network fusion
CN108960077A (en) * 2018-06-12 2018-12-07 南京航空航天大学 A kind of intelligent failure diagnosis method based on Recognition with Recurrent Neural Network
CN109522600A (en) * 2018-10-16 2019-03-26 浙江大学 Complex equipment remaining life prediction technique based on combined depth neural network
CN109726524A (en) * 2019-03-01 2019-05-07 哈尔滨理工大学 A prediction method for the remaining service life of rolling bearings based on CNN and LSTM
CN110108456A (en) * 2019-04-16 2019-08-09 东南大学 A kind of rotating machinery health evaluating method of depth convolutional neural networks

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"Deep Recurrent Convolutional Neural Network for Remaining Useful Life Prediction";Meng Ma 等;《2019 IEEE International Conference on Prognostics and Health Management (ICPHM)》;20190829;全文 *
"基于全矢卷积循环神经网络的滚动轴承诊断研究";谢远东;《中国优秀博硕士学位论文全文数据库(硕士)-工程科技Ⅱ辑》;20190815;第2019年卷(第8期);C029-133 *

Also Published As

Publication number Publication date
CN110633792A (en) 2019-12-31

Similar Documents

Publication Publication Date Title
CN110633792B (en) End-to-end bearing health index construction method based on convolution cyclic neural network
CN111709448B (en) A Mechanical Fault Diagnosis Method Based on Transition Relational Network
CN110866314B (en) Rotating Machinery Remaining Lifetime Prediction Method Based on Multilayer Bidirectionally Gated Recurrent Unit Networks
CN106779069A (en) A kind of abnormal electricity consumption detection method based on neutral net
CN112149316A (en) Prediction method of remaining life of aero-engine based on improved CNN model
CN114429153B (en) Gearbox incremental fault diagnosis method and system based on lifelong learning
CN112651167A (en) Semi-supervised rolling bearing fault diagnosis method based on graph neural network
CN110108456A (en) A kind of rotating machinery health evaluating method of depth convolutional neural networks
CN110689171A (en) A method for predicting the state of health of steam turbines based on E-LSTM
CN110657984A (en) A Fault Diagnosis Method of Planetary Gearbox Based on Enhanced Capsule Network
CN112488235A (en) Elevator time sequence data abnormity diagnosis method based on deep learning
CN111459144A (en) A Fault Prediction Method of Aircraft Flight Control System Based on Deep Recurrent Neural Network
CN111507046B (en) Method and system for predicting remaining service life of electric gate valve
CN114282443A (en) Remaining service life prediction method based on MLP-LSTM supervised joint model
CN113391211B (en) A method for predicting the remaining life of lithium batteries under the condition of small samples
CN115330082A (en) PM2.5 concentration prediction method of LSTM network based on attention mechanism
CN113435321A (en) Method, system and equipment for evaluating state of main shaft bearing and readable storage medium
CN118310746A (en) An unsupervised generator bearing fault detection method based on variational autoencoder
CN117763488A (en) Bearing residual service life prediction method based on operation phase division
CN112101659A (en) Complex equipment service life prediction method based on stacking denoising autoencoder
CN114942140A (en) Rolling bearing fault diagnosis method based on multi-input parallel graph convolution neural network
CN118378167B (en) Transmission chain fault diagnosis method based on zero sequence current correlation entropy reinforcement sparse learning
CN115438897A (en) Industrial process product quality prediction method based on BLSTM neural network
CN111738398B (en) A Novel Deep Feature Learning Method for Planetary Gear Fault Diagnosis
CN118378524A (en) Water pump unit health trend prediction and model training method thereof

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