CN109583386B - 一种智能旋转机械故障深度网络特征辨识方法 - Google Patents

一种智能旋转机械故障深度网络特征辨识方法 Download PDF

Info

Publication number
CN109583386B
CN109583386B CN201811458792.1A CN201811458792A CN109583386B CN 109583386 B CN109583386 B CN 109583386B CN 201811458792 A CN201811458792 A CN 201811458792A CN 109583386 B CN109583386 B CN 109583386B
Authority
CN
China
Prior art keywords
fault
rotary machine
matrix
reconstructed
sequence
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
CN201811458792.1A
Other languages
English (en)
Other versions
CN109583386A (zh
Inventor
刘辉
龙治豪
李燕飞
段铸
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Central South University
Original Assignee
Central South 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 Central South University filed Critical Central South University
Priority to CN201811458792.1A priority Critical patent/CN109583386B/zh
Publication of CN109583386A publication Critical patent/CN109583386A/zh
Application granted granted Critical
Publication of CN109583386B publication Critical patent/CN109583386B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/08—Feature extraction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00—Testing of machine parts
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/24—Classification techniques
    • G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/12—Classification; Matching

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • Signal Processing (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种智能旋转机械故障深度网络特征辨识方法，通过在列车滚动轴承待检测旋转机械部位设置振动传感器，收集滚动轴承工作时的原始振动序列，然后通过奇异谱分析方法将原始振动序列进行分解重构，并提取重构振动序列的均方根值、标准差、偏度指标和峰值，利用支持向量机训练得到的旋转机械故障位置诊断模型来判断故障位置，然后将重构振动序列进行集合经验模态分解，计算分解后的一组固有模态分量各自的排列熵值，将排列熵值的排列组合作为检测特征，利用支持向量机训练得到的旋转机械故障类型诊断模型来判断故障类型。本发明能更及时地检测出旋转机械的故障位置和故障类型，提高故障诊断的准确性和可靠性。

Description

一种智能旋转机械故障深度网络特征辨识方法
技术领域
本发明涉及机械系统故障辨识领域，特别是指一种智能旋转机械故障深度网络特征辨识方法。
背景技术
随着高速铁路技术的不断进步和智能高铁计划的提出，高铁运营安全越来越受到人们的关注。高速铁路组成复杂，而旋转机械在其中占据了重要地位，例如转向架电机，牵引制动装置都包含了大量旋转机械。不过旋转机械在长期使用过程中，极其容易出现不同程度的损耗而产生各种故障，如果发现不及时导致机械故障积累，轻则出现晚点事故导致经济损失，重则产生安全隐患造成安全事故。
现在旋转机械的诊断方法以采集振动信号为主，对其进行分析诊断。但是振动信号在故障早期异常变化极其微小，常规诊断方法往往效果不佳。现有的诊断算法中同时针对诊断部位和诊断类型进行检测，导致训练样本维度高，而且对原始数据的前处理较为简单，效果一般，故而导致最终模型检测精度不高。
因此，为了进一步提高高速铁路的稳定运营和安全运行，需要更好更及时的检测出旋转机械的故障部位和故障类型，提高故障诊断的准确性和可靠性，一种智能可靠的机械故障诊断方法急需被提出。
发明内容
本发明提供一种智能旋转机械故障深度网络特征辨识方法，通过对振动信号进行特征提取，基于深度网络对特征进行辨识诊断旋转机械是否发生故障、故障位置以及故障类型，提高故障诊断的准确性和可靠性。
为实现上述技术目的，本发明采用如下技术方案：
一种智能旋转机械故障深度网络特征辨识方法，包括以下步骤：
步骤1，获取旋转机械的原始振动信号；
步骤2，利用奇异谱分析方法，对原始振动信号进行分解重构得到重构振动序列；
步骤3，计算重构振动序列的均方根值RMSE、标准差SDE、偏度指标Skewness和峰值M；
步骤4，将重构振动序列进行集合经验模态分解，得到一组固有模态分量C1(t),C2(t),…Cn(t)，并计算得到各个固有模态分量Ci(t)的排列熵值Hi；
步骤5，构建训练样本，训练诊断模型；
将旋转机械的原始振动信号对应的均方根值RMSE、标准差SDE、偏度指标Skewness、峰值M和机械故障位置，作为旋转机械的第一训练样本；
将旋转机械的原始振动信号对应的各个固有模态分量Ci(t)的排列熵值Hi和机械故障类型，作为旋转机械的第二训练样本；
以旋转机械的第一训练样本为输入、故障位置为输出，训练第一支持向量机，得到旋转机械故障位置诊断模型；
以旋转机械的第二训练样本为输入、故障类型为输出，训练第二支持向量机，得到旋转机械故障类型诊断模型；
步骤6，实时采集待检测旋转机械的原始振动信号，按步骤2-3获取待检测旋转机械的原始振动信号所对应的均方根值RMSE、标准差SDE、偏度指标Skewness和峰值M，并输入到旋转机械故障位置诊断模型，旋转机械故障位置诊断模型输出待检测旋转机械的故障位置；
当旋转机械故障位置诊断模型输出的故障位置为无故障，则结束故障检测；否则执行步骤7；
步骤7，按步骤4获取待检测旋转机械的重构振动信号所对应的各个固有模态分量Ci(t)的排列熵值Hi，并输入到旋转机械故障类型诊断模型，旋转机械故障类型诊断模型输出待检测旋转机械的故障类型。
进一步地，所述利用奇异谱分析方法，对原始振动信号进行分解重构得到重构振动序列的具体过程为：
步骤A1，嵌入操作；
将原始振动序列X＝(x1，x2，x3…，xN)变换为Z＝[Z1，Z2，Z3…，ZK]，其中，Zi＝(xi，xi+1，xi+2,…，xi+L-1)，K＝N-L+1，L表示嵌入维数，且L∈[2,N]；最终映射的结果为轨迹矩阵Z：
Figure GDA0002542929020000021
其中，xi表示原始振动信号，1≤i≤N，N为每T个时间间隔的采样周期内采集的样本数目，T时间内采集的原始信号序列汇总为原始振动序列X＝(x1，x2，x3…，xN)；K为轨迹矩阵Z的列数；
步骤A2，将轨迹矩阵Z进行奇异值分解；
设矩阵H＝ZZT，且λ1,λ2,…λL为矩阵H的特征值，λ1≥λ2≥…≥λL≥0，其中ZT是轨迹矩阵Z的转置矩阵，通过奇异值分解的方法，将轨迹矩阵Z进行如下变换从而得到Zf；
Zf＝Z1+Z2+…+Zd，
Figure GDA0002542929020000031
其中
Figure GDA0002542929020000032
为轨迹矩阵Z的奇异值，Ui和Vi分别表示轨迹矩阵Z的左右特征向量，其中d＝L*＝min{L,K}；
步骤A3，将分解后的矩阵Zf进行分组；
将区间{1,2,…,d}分割为非连接的子集I1,I2,…Im，则轨迹矩阵Z可表示为：
Figure GDA0002542929020000033
步骤A4，将分组后的轨迹矩阵Z重构为长度为N的重构振动序列Xc；
轨迹矩阵Z是L×K的矩阵，L*＝min(L,K)，K*＝max(L,K)，N＝L+K-1，当L＜K时，
Figure GDA0002542929020000034
否则
Figure GDA0002542929020000035
重构振动序列为Xc＝(xc1,xc2,…,xcn)，且xcn为：
Figure GDA0002542929020000036
进一步地，所述均方根值RMSE、标准差SDE、偏度指标Skewness和峰值M的计算公式为：
Figure GDA0002542929020000037
Figure GDA0002542929020000041
Figure GDA0002542929020000042
M＝max|xcn(i)|；
其中，xcn(i)表示重构振动序列中的第i个元素，
Figure GDA0002542929020000043
表示重构振动序列中元素的均值。
进一步地，所述将重构振动序列进行集合经验模态分解的具体过程为：
步骤B1，求得重构振动序列Xc＝(xc1,xc2,…,xcn)中的极大值和极小值，采用三次样条函数进行插值拟合上包络线bmax(t)和下包络线bmin(t)；
步骤B2，计算上下包络线的平均值m(t)，提取h(t)＝x(t)-m(t)，判断h(t)是否满足固有模态分量条件，如果满足，则设h(t)为第一个固有模态分量；否则将h(t)作为新的重构振动序列Xc；
步骤B3，重复步骤B1和B2，直到经过n次筛选后的差值hn(t)满足固有模态分解条件，则将该差值hn(t)设为第一个固有模态分量，记作C1(t)＝hn(t)；
步骤B4，将重构振动序列Xc转化为剩余序列R1(t)＝Xc-C1(t)；
步骤B5，将剩余序列R1(t)重复步骤B1-B4，得到其余的固有模态分量，直到余量Rn(t)为单调函数时终止循环，得到一组固有模态分量C1(t),C2(t),…Cn(t)。
进一步地，所述各个固有模态分量Ci(t)的排列熵值Hi的具体过程为：
步骤C1，对固有模态分量C1(t)进行相空间重构，得到矩阵
Figure GDA0002542929020000044
其中j＝1,2,...,W，d为嵌入维数，τ为延迟时间，W为重构分量的个数；
步骤C2，将矩阵Y中每一行{x(j),x(j+τ),…,x(j+(d-1)τ}按照升序重新排列，将每一行重新表示为{x(j1),x(j2+τ),…,x(jd+(d-1)τ}，从而每一行得到一组符号序列：
S(l)＝{j1,j2,…,jd}；
式中l＝1,2,...,W；
步骤C3，计算每一组符号序列出现的概率，得到P1,P2,…PW，那么固有模态分量C1(t)的排列熵值为：
Figure GDA0002542929020000051
步骤C4，对剩余的固有模态分量重复步骤C1-C3，计算出所有固有模态分量的排列熵值，按照顺序排列为H1,H2,…Hn。
进一步地，当旋转机械故障类型诊断模型输出的故障类型为无故障，则重复步骤7；当重复步骤7达到预设次数时，输出的故障类型仍为无故障，则结束故障检测。
进一步地，所述故障位置包括无故障、外圈故障、内圈故障、滚动体故障和保持架故障。
进一步地，所述故障类型包括无故障、疲劳剥落、磨损、腐蚀和胶合。
进一步地，所述训练样本集包括低速运行、中速运行和高速运行三种工况下的旋转机械。
进一步地，所述旋转机械为滚动轴承。
有益效果
1、对监测对象即旋转机械所执行的故障位置检测和故障类型检测分开使用两个支持向量机进行，简化单个支持向量机的训练样本的特征数量，有利于提高对支持向量机的训练率，进而提高对旋转机械故障检测的准确率；
2、利用奇异谱分析方法将原始振动信号分解重构为重构振动序列，起到主特征提取和降噪的作用，提高后续对旋转机械故障检测的准确率；
3、利用集合经验模态分解将主成分，即重构振动序列，分解为一组固有模态分量，计算每个固有模态分量的排列熵值，从而通过位置诊断模型根据该组排列熵值来区分不同的故障类型，该排列熵值便于计算、易于区分、标签空间广，并且可以提取并放大重构振动序列中的微小变化，可以有效区分多种故障类型；
4、提取的特征，包括均方根值RMSE、标准差SDE、偏度指标Skewness、峰值M，固有模态分量的排列熵，属于时频特征，对于旋转机械的非平稳非线性信号的描述性更好，即对更能体现旋转机械的非平稳非线性信号所具有的特征，从而更进一步提高对旋转机械故障检测的准确率。
附图说明
图1为本发明智能旋转机械故障深度网络特征辨识方法的整体流程图。
具体实施方式
本发明实施例提供的智能旋转机械故障深度网络特征辨识方法，通过在列车滚动轴承待检测旋转机械部位设置振动传感器，收集滚动轴承工作时的原始振动序列，然后通过奇异谱分析方法将原始振动序列进行分解重构，并提取重构振动序列的均方根值、标准差、偏度指标和峰值，利用支持向量机训练得到的旋转机械故障位置诊断模型来判断故障位置，然后将重构振动序列进行集合经验模态分解，计算分解后的一组固有模态分量各自的排列熵值，将排列熵值的排列组合作为检测特征，利用支持向量机训练得到的旋转机械故障类型诊断模型来判断故障类型。
在本实施例中旋转机械以滚动轴承为例，如图1所示，本发明方法包括如下步骤：
步骤1，获取训练样本集。
在实验环境下，通过布置在轴承检测处的振动传感器，分别采集滚动轴承的原始振动信号xi，其中1≤i≤N，N为每T个时间间隔的采样周期内采集的样本数目，这样T时间内采集的原始信号序列汇总为原始振动序列X＝(x1，x2，x3…，xN)。
在滚动轴承正常运行于低速、中速和高速三种工况，每次对滚动轴承的内圈、外圈、滚动体和保持架四个部位均正常运作时的原始振动序列各采集100组，由于速度越快发生故障的概率和危害越大，因此在本实施例采集滚动轴承的原始振动信号训练模型时，低速运行取1次数据、中速运行取2次数据，高速运行取3次数据，从而提高训练得到的故障位置诊断模型和故障类型诊断模型对故障诊断的准确性。四个部位分别发生故障时(即某个部位故障、另外三个部位正常)的原始振动序列各采集100组，如此对故障位置诊断模型和故障类型诊断模型的训练效果更好。
步骤2，利用奇异谱分析方法，对原始振动信号进行分解重构得到重构振动序列，具体步骤为：
步骤A1，嵌入操作；
将原始振动信号X＝(x1，x2，x3…，xN)变换为Z＝[Z1，Z2，Z3…，ZK]，其中，Zi＝(xi，xi+1，xi+2,…，xi+L-1)，K＝N-L+1，L表示嵌入维数，且L∈[2,N]；最终映射的结果为轨迹矩阵Z：
Figure GDA0002542929020000071
步骤A2，将轨迹矩阵Z进行奇异值分解；
设矩阵H＝ZZT，且λ1,λ2,…λL为矩阵H的特征值，λ1≥λ2≥…≥λL≥0，其中ZT是轨迹矩阵Z的转置矩阵，通过奇异值分解的方法，将轨迹矩阵Z进行如下变换从而得到Zf；
Zf＝Z1+Z2+…+Zd，
Figure GDA0002542929020000072
其中
Figure GDA0002542929020000073
为轨迹矩阵Z的奇异值，Ui和Vi分别表示轨迹矩阵Z的左右特征向量，其中d＝L*＝min{L,K}；
步骤A3，将分解后的矩阵Zf进行分组；
将区间{1,2,…,d}分割为非连接的子集I1,I2,…Im，则轨迹矩阵Z可表示为：
Figure GDA0002542929020000074
步骤A4，将分组后的轨迹矩阵Z重构为长度为N的重构振动序列Xc；
轨迹矩阵Z是L×K的矩阵，L*＝min(L,K)，K*＝max(L,K)，N＝L+K-1，当L＜K时，
Figure GDA0002542929020000075
否则
Figure GDA0002542929020000076
重构振动序列为Xc＝(xc1,xc2,…,xcn)，且xcn为：
Figure GDA0002542929020000077
利用奇异谱分析方法将原始振动信号分解重构为重构振动序列，起到主特征提取和降噪的作用，提高后续对旋转机械故障检测的准确率。
步骤3，计算重构振动序列Xc的均方根值RMSE、标准差SDE、偏度指标Skewness和峰值M，其计算公式为：
Figure GDA0002542929020000081
Figure GDA0002542929020000082
Figure GDA0002542929020000083
M＝max|xcn(i)|；
其中，xcn(i)表示重构振动序列中的第i个元素，
Figure GDA0002542929020000084
表示重构振动序列中元素的均值。
这些提取到的重构振动序列Xc的均方根值RMSE、标准差SDE、偏度指标Skewness和峰值M属于时频特征，对于旋转机械的非平稳非线性信号的描述性更好，即对更能体现旋转机械的非平稳非线性信号所具有的特征，从而更进一步提高对旋转机械故障检测的准确率。
步骤4，将重构振动序列进行集合经验模态分解，得到一组固有模态分量C1(t),C2(t),…Cn(t)，并计算得到各个固有模态分量Ci(t)的排列熵值Hi。
其中，对将重构振动序列Xc进行集合经验模态分解的具体过程为：
步骤B1，求得重构振动序列Xc＝(xc1,xc2,…,xcn)中的极大值和极小值，采用三次样条函数进行插值拟合上包络线bmax(t)和下包络线bmin(t)；
步骤B2，计算上下包络线的平均值m(t)，提取h(t)＝x(t)-m(t)，判断h(t)是否满足固有模态分量条件，如果满足，则设h(t)为第一个固有模态分量；否则将h(t)作为新的重构振动序列Xc；
在本实施例中，选用标准差SDE作为筛选停止标准，当h(t)的标准差SDE为0.2～0.3时，则认为该h(t)符合条件；
步骤B3，重复步骤B1和B2，直到经过n次筛选后的差值hn(t)满足固有模态分解条件，则将该差值hn(t)设为第一个固有模态分量，记作C1(t)＝hn(t)；
步骤B4，得到C1(t)后，将重构振动序列Xc转化为剩余序列R1(t)＝Xc-C1(t)；
步骤B5，将剩余序列R1(t)重复步骤B1-B4，得到其余的固有模态分量，直到余量Rn(t)为单调函数时终止循环，得到一组固有模态分量C1(t),C2(t),…Cn(t)。
另外，计算各个固有模态分量Ci(t)的排列熵值Hi的具体过程为：
步骤C1，对固有模态分量C1(t)进行相空间重构，得到矩阵
Figure GDA0002542929020000091
其中j＝1,2,…,W，d为嵌入维数，τ为延迟时间，W为重构分量的个数；
步骤C2，将矩阵Y中每一行{x(j),x(j+τ),…,x(j+(d-1)τ}按照升序重新排列，那么每一行都可以重新表示为{x(j1),x(j2+τ),…,x(jd+(d-1)τ}，从而每一行都可以得到一组符号序列：
S(l)＝{j1,j2,…,jd}；
式中l＝1,2,…,W；
步骤C3，计算每一组符号序列出现的概率，得到P1,P2,…PW，那么固有模态分量C1(t)的排列熵值为：
Figure GDA0002542929020000092
步骤C4，对剩余的固有模态分量重复步骤C1-C3，计算出所有固有模态分量的排列熵值，按照顺序排列为H1,H2,…Hn。
利用集合经验模态分解将主成分，即重构振动序列，分解为一组固有模态分量，计算每个固有模态分量的排列熵值，从而通过位置诊断模型根据该组排列熵值来区分不同的故障类型，该排列熵值便于计算、易于区分；而且，利用集合经验模态分解可得到大量模态分量，从而有大量对应的排列熵值，将所有的排列熵值排列组合可得到大量的排列熵，这样繁多的排列组合便构成了宽广的标签空间，并且可以提取并放大重构振动序列中的微小变化，可以有效区分多种故障类型；
步骤5，构建训练样本，训练诊断模型；
将旋转机械的原始振动信号对应的均方根值RMSE、标准差SDE、偏度指标Skewness、峰值M和机械故障位置，作为旋转机械的第一训练样本；
将旋转机械的原始振动信号对应的各个固有模态分量Ci(t)的排列熵值Hi和机械故障类型，作为旋转机械的第二训练样本。
以旋转机械的第一训练样本为输入、故障位置为输出，训练第一支持向量机，得到旋转机械故障位置诊断模型。其中，故障位置包括无故障、外圈故障、内圈故障、滚动体故障和保持架故障，在第一支持向量机模型中分别表示为0、1、2、3、4。
以旋转机械的第二训练样本为输入、故障类型为输出，训练第二支持向量机，得到旋转机械故障类型诊断模型。其中，故障类型包括无故障、疲劳剥落、磨损、腐蚀和胶合，在第二支持向量机模型中分别表示为0、1、2、3、4。
步骤6，在线检测。
利用布置在待检测滚动轴承处的振动传感器，实时采集待检测滚动轴承的原始振动信号，与训练旋转机械故障位置诊断模型时获取第一训练样本的均方根值RMSE、标准差SDE、偏度指标Skewness和峰值M的方法相同，获取待检测滚动轴承的原始振动信号的均方根值RMSE、标准差SDE、偏度指标Skewness峰值M，并将该四个值输入到旋转机械故障位置诊断模型，旋转机械故障位置诊断模型输出待检测旋转机械的故障位置。
当旋转机械故障位置诊断模型输出的故障位置为无故障，则结束故障检测；否则执行步骤7；
步骤7，与训练旋转机械故障类型诊断模型时获取第二训练样本的的方法相同，按步骤4获取待检测旋转机械的重构振动信号所对应的各个固有模态分量Ci(t)的排列熵值Hi，并将排列熵值组成的序列H1,H2,…Hn输入到旋转机械故障类型诊断模型，旋转机械故障类型诊断模型输出待检测旋转机械的故障类型。
当旋转机械故障类型诊断模型输出的故障类型为无故障，则重复步骤7；当重复步骤7达到预设次数3时，输出的故障类型仍为无故障，则结束故障检测。
当旋转机械故障位置诊断模型输出故障位置、旋转机械故障类型诊断模型输出故障类型，将故障位置和故障类型发送给人机交互终端或者诊断结果响应设备，以进一步采取故障排除措施。
至此，本领域技术人员应认识到，虽然本文已详尽示出和描述了本发明的示例性实施例，但是，在不脱离本发明精神和范围的情况下，仍可根据本发明公开的内容直接确定或推导出符合本发明原理的许多其他变形或修改。因此，本发明的范围应该被理解和认定为覆盖了所有这些其他变形或修改。

Claims (10)

1.一种智能旋转机械故障深度网络特征辨识方法，其特征在于，包括以下步骤：
步骤1，获取旋转机械的原始振动信号；
步骤2，利用奇异谱分析方法，对原始振动信号进行分解重构得到重构振动序列；
步骤3，计算重构振动序列的均方根值RMSE、标准差SDE、偏度指标Skewness和峰值M；
步骤4，将重构振动序列进行集合经验模态分解，得到一组固有模态分量C1(t),C2(t),…Cn(t)，并计算得到各个固有模态分量Ci(t)的排列熵值Hi；
步骤5，构建训练样本，训练诊断模型；
将旋转机械的原始振动信号对应的均方根值RMSE、标准差SDE、偏度指标Skewness、峰值M和机械故障位置，作为旋转机械的第一训练样本；
将旋转机械的原始振动信号对应的各个固有模态分量Ci(t)的排列熵值Hi和机械故障类型，作为旋转机械的第二训练样本；
以旋转机械的第一训练样本为输入、故障位置为输出，训练第一支持向量机，得到旋转机械故障位置诊断模型；
以旋转机械的第二训练样本为输入、故障类型为输出，训练第二支持向量机，得到旋转机械故障类型诊断模型；
步骤6，实时采集待检测旋转机械的原始振动信号，按步骤2-3获取待检测旋转机械的原始振动信号所对应的均方根值RMSE、标准差SDE、偏度指标Skewness和峰值M，并输入到旋转机械故障位置诊断模型，旋转机械故障位置诊断模型输出待检测旋转机械的故障位置；
当旋转机械故障位置诊断模型输出的故障位置为无故障，则结束故障检测；否则执行步骤7；
步骤7，按步骤4获取待检测旋转机械的重构振动信号所对应的各个固有模态分量Ci(t)的排列熵值Hi，并输入到旋转机械故障类型诊断模型，旋转机械故障类型诊断模型输出待检测旋转机械的故障类型。
2.根据权利要求1所述的方法，其特征在于，所述利用奇异谱分析方法，对原始振动信号进行分解重构得到重构振动序列的具体过程为：
步骤A1，嵌入操作；
将原始振动序列X＝(x1，x2，x3…，xN)变换为Z＝[Z1，Z2，Z3…，ZK]，其中，Zi＝(xi，xi+1，xi+2,…，xi+L-1)，K＝N-L+1，L表示嵌入维数，且L∈[2,N]；最终映射的结果为轨迹矩阵Z：
Figure FDA0002542929010000021
其中，xi表示原始振动信号，1≤i≤N，N为每T个时间间隔的采样周期内采集的样本数目，T时间内采集的原始信号序列汇总为原始振动序列X＝(x1，x2，x3…，xN)；K为轨迹矩阵Z的列数；
步骤A2，将轨迹矩阵Z进行奇异值分解；
设矩阵H＝ZZT，且λ1,λ2,…λL为矩阵H的特征值，λ1≥λ2≥…≥λL≥0，其中ZT是轨迹矩阵Z的转置矩阵，通过奇异值分解的方法，将轨迹矩阵Z进行如下变换从而得到Zf；
Zf＝Z1+Z2+…+Zd，
Figure FDA0002542929010000022
其中
Figure FDA0002542929010000023
为轨迹矩阵Z的奇异值，Ui和Vi分别表示轨迹矩阵Z的左右特征向量，其中d＝L*＝min{L,K}；
步骤A3，将分解后的矩阵Zf进行分组；
将区间{1,2,…,d}分割为非连接的子集I1,I2,…Im，则轨迹矩阵Z可表示为：
Figure FDA0002542929010000024
步骤A4，将分组后的轨迹矩阵Z重构为长度为N的重构振动序列Xc；
轨迹矩阵Z是L×K的矩阵，L*＝min(L,K)，K*＝max(L,K)，N＝L+K-1，当L＜K时，
Figure FDA0002542929010000025
否则
Figure FDA0002542929010000026
重构振动序列为Xc＝(xc1,xc2,…,xcn)，且xcn为：
Figure FDA0002542929010000031
3.根据权利要求2所述的方法，其特征在于，所述均方根值RMSE、标准差SDE、偏度指标Skewness和峰值M的计算公式为：
Figure FDA0002542929010000032
Figure FDA0002542929010000033
Figure FDA0002542929010000034
M＝max|xcn(i)|；
其中，xcn(i)表示重构振动序列中的第i个元素，
Figure FDA0002542929010000035
表示重构振动序列中元素的均值。
4.根据权利要求2所述的方法，其特征在于，所述将重构振动序列进行集合经验模态分解的具体过程为：
步骤B1，求得重构振动序列Xc＝(xc1,xc2,…,xcn)中的极大值和极小值，采用三次样条函数进行插值拟合上包络线bmax(t)和下包络线bmin(t)；
步骤B2，计算上下包络线的平均值m(t)，提取h(t)＝x(t)-m(t)，判断h(t)是否满足固有模态分量条件，如果满足，则设h(t)为第一个固有模态分量；否则将h(t)作为新的重构振动序列Xc；
步骤B3，重复步骤B1和B2，直到经过n次筛选后的差值hn(t)满足固有模态分解条件，则将该差值hn(t)设为第一个固有模态分量，记作C1(t)＝hn(t)；
步骤B4，将重构振动序列Xc转化为剩余序列R1(t)＝Xc-C1(t)；
步骤B5，将剩余序列R1(t)重复步骤B1-B4，得到其余的固有模态分量，直到余量Rn(t)为单调函数时终止循环，得到一组固有模态分量C1(t),C2(t),…Cn(t)。
5.根据权利要求4所述的方法，其特征在于，所述各个固有模态分量Ci(t)的排列熵值Hi的具体过程为：
步骤C1，对固有模态分量C1(t)进行相空间重构，得到矩阵
Figure FDA0002542929010000041
其中j＝1,2,...,W，d为嵌入维数，τ为延迟时间，W为重构分量的个数；
步骤C2，将矩阵Y中每一行{x(j),x(j+τ),…,x(j+(d-1)τ}按照升序重新排列，将每一行重新表示为{x(j1),x(j2+τ),…,x(jd+(d-1)τ}，从而每一行得到一组符号序列：
S(l)＝{j1,j2,…,jd}；
式中l＝1,2,...,W；
步骤C3，计算每一组符号序列出现的概率，得到P1,P2,…PW，那么固有模态分量C1(t)的排列熵值为：
Figure FDA0002542929010000042
步骤C4，对剩余的固有模态分量重复步骤C1-C3，计算出所有固有模态分量的排列熵值，按照顺序排列为H1,H2,…Hn。
6.根据权利要求1所述的方法，其特征在于，当旋转机械故障类型诊断模型输出的故障类型为无故障，则重复步骤7；当重复步骤7达到预设次数时，输出的故障类型仍为无故障，则结束故障检测。
7.根据权利要求1所述的方法，其特征在于，所述故障位置包括无故障、外圈故障、内圈故障、滚动体故障和保持架故障。
8.根据权利要求1所述的方法，其特征在于，所述故障类型包括无故障、疲劳剥落、磨损、腐蚀和胶合。
9.根据权利要求1所述的方法，其特征在于，所述训练样本集包括低速运行、中速运行和高速运行三种工况下的旋转机械。
10.根据权利要求1所述的方法，其特征在于，所述旋转机械为滚动轴承。
CN201811458792.1A 2018-11-30 2018-11-30 一种智能旋转机械故障深度网络特征辨识方法 Active CN109583386B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811458792.1A CN109583386B (zh) 2018-11-30 2018-11-30 一种智能旋转机械故障深度网络特征辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811458792.1A CN109583386B (zh) 2018-11-30 2018-11-30 一种智能旋转机械故障深度网络特征辨识方法

Publications (2)

Publication Number Publication Date
CN109583386A CN109583386A (zh) 2019-04-05
CN109583386B true CN109583386B (zh) 2020-08-25

Family

ID=65926345

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811458792.1A Active CN109583386B (zh) 2018-11-30 2018-11-30 一种智能旋转机械故障深度网络特征辨识方法

Country Status (1)

Country Link
CN (1) CN109583386B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110160813A (zh) * 2019-05-22 2019-08-23 湃方科技(北京)有限责任公司 旋转设备状态监测方法及装置
CN110108474B (zh) * 2019-06-04 2020-05-15 山东大学 一种旋转机械运行稳定性在线监测与评估方法及系统
CN110232415B (zh) * 2019-06-12 2020-10-30 中南大学 一种基于生物信息特征的列车转向架故障识别方法
CN110361505B (zh) * 2019-07-25 2021-06-22 中南大学 一种车外大气污染环境下列车乘员健康预警系统的方法
CN110781854B (zh) * 2019-11-04 2022-07-26 上海联影智能医疗科技有限公司 故障检测模型的训练方法和机电设备的故障检测方法
CN111310711B (zh) * 2020-03-03 2022-10-21 广东工业大学 基于二维奇异谱分析融合emd的人脸图像识别方法及系统
CN111444988B (zh) * 2020-05-11 2021-09-17 北华大学 一种滚动轴承故障诊断系统
CN112149845A (zh) * 2020-09-23 2020-12-29 山东通维信息工程有限公司 基于大数据和机器学习的智能运维方法
CN112699601B (zh) * 2020-12-28 2022-05-31 电子科技大学 一种传感器网络数据的空时重构方法
CN113065418A (zh) * 2021-03-18 2021-07-02 上海电力大学 一种基于ssa-wdcnn的滚动轴承故障诊断方法
CN114088400B (zh) * 2021-11-01 2024-04-09 中国人民解放军92728部队 一种基于包络排列熵的滚动轴承故障诊断方法
CN115571196A (zh) * 2022-11-17 2023-01-06 北京乐码仕智能科技有限公司 一种基于规则引擎的正线故障智能处置方法及系统

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101750210B (zh) * 2009-12-24 2012-05-30 重庆大学 基于正交局部保持映射（olpp）特征约简的故障诊断方法
CN104833534A (zh) * 2015-04-21 2015-08-12 广州市地下铁道总公司 一种基于多源信息融合的列车走行系故障诊断装置及方法
CN105275833B (zh) * 2015-10-30 2017-01-11 北京航空航天大学 一种基于CEEMD-STFT时频信息熵和multi-SVM的离心泵故障诊断方法
CN108254179A (zh) * 2017-08-08 2018-07-06 常州路航轨道交通科技有限公司 一种基于meemd排列熵的高速列车轮对轴承故障诊断方法
CN108332970A (zh) * 2017-11-17 2018-07-27 中国铁路总公司 一种基于ls-svm和d-s证据理论的轴承故障诊断方法

Also Published As

Publication number Publication date
CN109583386A (zh) 2019-04-05

Similar Documents

Publication Publication Date Title
CN109583386B (zh) 一种智能旋转机械故障深度网络特征辨识方法
CN111353482B (zh) 一种基于lstm的疲劳因子隐性异常检测及故障诊断方法
CN111089726B (zh) 一种基于最优维数奇异谱分解的滚动轴承故障诊断方法
CN107451004B (zh) 一种基于定性趋势分析的道岔故障诊断方法
CN100437072C (zh) 离心泵振动故障融合诊断方法及振动信号采集装置
Zhou et al. Mechanical equipment fault diagnosis based on redundant second generation wavelet packet transform
CN107657250B (zh) 轴承故障检测及定位方法及检测定位模型实现系统和方法
Dou et al. A rule-based intelligent method for fault diagnosis of rotating machinery
CN105973595A (zh) 一种滚动轴承故障的诊断方法
CN109946080B (zh) 一种基于嵌入式循环网络的机械设备健康状态识别方法
CN106226074A (zh) 基于卷积神经网络和小波灰度图的旋转机械故障诊断方法
CN109255333B (zh) 一种大型风电机组滚动轴承故障混合诊断方法
CN105738722A (zh) 一种航空机电作动器故障诊断方法
CN112182912B (zh) 基于概率描述与谱分析的制造装备主轴轴承健康评估方法
CN111931625A (zh) 基于非对称损失神经网络的产品关键零件剩余寿命预测方法
CN112116029A (zh) 一种多尺度结构与特征融合的齿轮箱智能故障诊断方法
CN107451624B (zh) 基于正交化局部连接网络的机械装备健康状态识别方法
CN110207974A (zh) 基于振动信号时频能量分布特征的断路器故障识别方法
Wang et al. Fault diagnosis of industrial robots based on multi-sensor information fusion and 1D convolutional neural network
CN115586023A (zh) 一种轨道车辆传动系统故障诊断方法及系统
CN113485863B (zh) 基于改进生成对抗网络生成异构不平衡故障样本的方法
CN110530631B (zh) 一种基于混合分类器的齿轮单类型故障检测方法
Yu et al. Rolling bearing fault feature extraction and diagnosis method based on MODWPT and DBN
CN110160781B (zh) 一种旋转机械故障分类的测试集重构和预测方法
CN113516159A (zh) 轨道车辆走行部小齿轮轴裂纹故障诊断方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant