CN114201989A - 一种采用卷积神经网络和双向长短时记忆网络的交流电机轴承故障诊断方法 - Google Patents

一种采用卷积神经网络和双向长短时记忆网络的交流电机轴承故障诊断方法 Download PDF

Info

Publication number
CN114201989A
CN114201989A CN202111430349.5A CN202111430349A CN114201989A CN 114201989 A CN114201989 A CN 114201989A CN 202111430349 A CN202111430349 A CN 202111430349A CN 114201989 A CN114201989 A CN 114201989A
Authority
CN
China
Prior art keywords
time
layer
output
neural network
input
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.)
Pending
Application number
CN202111430349.5A
Other languages
English (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.)
Jiangsu University
Original Assignee
Jiangsu 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 Jiangsu University filed Critical Jiangsu University
Priority to CN202111430349.5A priority Critical patent/CN114201989A/zh
Publication of CN114201989A publication Critical patent/CN114201989A/zh
Pending legal-status Critical Current

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
    • G01M13/04—Bearings
    • G01M13/045—Acoustic or vibration analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/04—Architecture, e.g. interconnection topology
    • G06N3/045—Combinations of networks
    • 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)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Evolutionary Computation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Signal Processing (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Acoustics & Sound (AREA)
  • Control Of Electric Motors In General (AREA)

Abstract

本发明公开了一种采用卷积神经网络和双向长短时记忆网络的交流电机轴承故障诊断方法。该方法充分利用卷积神经网络的参数共享、局部感知和降采样等优点，有效的从原始电流数据中提取到相应地空间特征，避免了复杂的特征提取过程。然后将提取的电流数据空间特征输入双向长短时记忆网络，捕获电流数据空间特征的时序信息，进一步提高滚动轴承故障诊断准确率。本发明利用电机驱动过程获取定子电流作为故障信号，具有非侵入式且易于形成闭环“驱动—诊断”的融合系统，可有效降低监控成本；所提轴承故障诊断方法结合卷积神经网络和双向长短时记忆神经网络特性，在特征提取的深度和复杂程度上都具有良好的性能，能够实现故障特征的精准有效提取。

Description

一种采用卷积神经网络和双向长短时记忆网络的交流电机轴 承故障诊断方法
技术领域
本发明涉及交流电机状态监测领域，具体地说是一种采用电流信号分析的交流电机轴承故障诊断方法。
背景技术
与传统采用振动信号分析的轴承故障诊断方法相比，受成本和空间约束，有些工作环境无法安装额外的硬件系统来满足交流电机轴承故障诊断的需求。此外，电机周围存在的较多其他振动噪声干扰源，严重影响振动信号分析诊断轴承故障的效果。通过电机驱动控制系统获取故障信号，实现轴承状态监测是解决这一问题的重要途径。
采用驱动控制系统的轴承故障诊断方法以电流、电压、转矩和转速等不同信号作为故障诊断依据。其中，电压和转矩信号通过电机系统采集较为困难，使用这两种信号诊断轴承故障的研究比较少见。采用电机电流信号特征分析诊断轴承故障的研究最多，近些年也受到越来越多的关注。
受环境噪声、基频分量和逆变器等因素干扰，电流信号分析诊断轴承故障属于强噪声背景下微弱故障特征提取问题。传统的故障诊断方法主要依靠先进信号处理技术，如去噪技术(维纳滤波器、谱减法、信号时移、线性预测等)、高分辨率谱分析(Root-MUSIC、ESPIRT、增强双谱、调制双谱等)、解调变换(派克变换、希尔伯特变换、能量算子等)、时频分析(短时傅里叶变换、维格纳分布、小波包变换、集合经验模态分解等)、时频域统计特征(独立主元分析、四元数旋转、拟合优度测试、伽玛分布等)等，从频域或者时频域提取电流信号的轴承故障特征进行识别。然而，这些方法精度低、适应性差且严重依赖人工经验。
利用不同信号处理技术，得到电流信号中表征滚动轴承运行状态的特征量后，并与支持向量机、自适应神经模糊系统、蚁群聚类算法等机器学习算法相结合，也可实现轴承故障智能识别。但是，此类方法严重依赖信号处理技术的特征提取能力，另外，浅层机器学习泛化能力差等问题，也限制了其诊断效果。
深度学习算法直接将原始电流信号作为输入，通过多层模型逐层学习数据特征，达到故障特征有效表达，为轴承故障智能诊断提供了有益思路。这类方法无需进行额外的去噪、特征提取等操作，自动挖掘代表性信息和敏感特征，有效克服了当前诊断方法的缺陷。近年来，一维/二维卷积神经网络等被应用于采用电流信号分析的轴承故障诊断领域，诊断性能优越。虽然卷积神经网络应用在轴承故障识别方面取得了较好的效果，但对于时间序列数据样本而言，其并不擅长处理。然而，当滚动轴承发生故障时，采集的电流信号是一种具有时间连续性的数据，充分利用电流信号的时空相关性是实现滚动轴承故障识别的关键。
发明内容
针对现有技术中存在的上述不足之处，本发明要解决的技术问题是提供一种采用卷积神经网络和双向长短时记忆网络的交流电机轴承故障诊断方法，在使用卷积神经网络从原始电流数据样本中提取得到相应地空间特征信息基础上，结合双向长短时记忆网络，捕获电流数据空间特征的时序信息，进而实现交流电机轴承故障诊断。
本发明为实现上述目的所采用的技术方案是：一种采用卷积神经网络和双向长短时记忆网络的交流电机轴承故障诊断方法，包括以下步骤：
步骤1，获取交流电机定子电流信号：设定采样频率fs和采样时间T，以一定时间间隔从交流电机的数据采集系统中获取三相定子电流信号，将其作为待处理的数据；
步骤2，定子电流数据样本点划分：以定子电流信号相邻波谷对应的时间为一个时间周期，截取原始定子电流信号中m个时间周期的数据点作为一个样本，将其作为交流电机轴承故障诊断模型的样本数据将其作为交流电机轴承故障诊断模型的样本数据；
步骤3，创建数据集：根据步骤2对交流电机轴承故障类型对应的样本数据进行编码，制作标签，作为交流电机轴承故障诊断模型的数据集；
步骤4，构建卷积神经网络：将训练数据输入卷积神经网络，其包括两层卷积层和一层池化层，通过卷积池化操作提取定子电流信号中有区分性的空间特征信息；
步骤5，构建双向长短时记忆神经网络：将提取的定子电流数据空间信息输入双向长短时记忆网络，捕获定子电流数据空间特征的时序信息；
步骤6，实现滚动轴承故障类型识别：将提取的定子电流数据时空信息输入到全连接层和输出层，选择合适的损失函数、学习率和训练迭代次数，对模型进行训练，得到训练好的具有交流电机轴承故障状态识别能力的诊断模型。
进一步，所述步骤1中，在满足奈奎斯特采样定理的条件下，为了包含更多的故障信息，采样频率fs设置为：
fs>10f
其中，f为定子电流信号频率。
进一步，所述步骤2中，以定子电流信号相邻波谷对应的时间为一个时间周期，截取原始定子电流信号中连续m个时间周期的数据点作为一个样本
Figure BDA0003380008400000031
Figure BDA0003380008400000032
其中
Figure BDA0003380008400000033
表示连续m个时间周期内a相、b相、c相电流数据，
Figure BDA0003380008400000034
Figure BDA0003380008400000035
Nsub为m个时间周期内单相电流数据的长度，将三相电流视为三个通道作为网络的输入，交流电机轴承故障诊断模型的输入数据样本可表示为
Figure BDA0003380008400000036
Figure BDA0003380008400000037
p＝1,2,…Nseg，Nseg为输入数据样本长度。
进一步，所述步骤3的具体过程为：
步骤3-1，确定标签的编码方式：采用one-hot编码；
步骤3-2，将无故障交流电机对应的定子电流数据和轴承故障交流电机对应的定子电流数据按照故障类型制作标签，同时以6:2:2的数据比例将定子电流数据分为训练集、验证集和测试集，作为交流电机轴承故障诊断模型的数据集；
进一步，所述步骤4中，卷积神经网络包括两层卷积层和一层池化层，具体计算过程为：
步骤4-1，第一层：卷积层，输入数据样本Ip∈I，卷积核的数量为FN，长度为FL，步长为1，则卷积层对输入Ip在第j个卷积核处的值
Figure BDA0003380008400000041
可通过以下公式进行计算：
Figure BDA0003380008400000042
Figure BDA0003380008400000043
其中，Hjc为第j个卷积核第c通道，
Figure BDA0003380008400000044
*T为矩阵*的转置，
Figure BDA0003380008400000045
为第j个卷积核第q个卷积区域值，
Figure BDA0003380008400000046
为第j个卷积核在通道c的偏移向量，fjc为第j个卷积核第c通道的激活函数，q:q+FL-1为第j个卷积核第c通道第q个卷积区域，
Figure BDA0003380008400000047
为第j个卷积核第c通道第q个卷积区域中的元素。
输入Ip在第一层卷积层的输出特征向量
Figure BDA0003380008400000048
数学表达式为：
Figure BDA0003380008400000049
步骤4-2，第二层，卷积层：输入第一层输出的特征向量
Figure BDA00033800084000000410
卷积核与第一层卷积层具有相同的数量FN和长度FL，步长为2，则第二层卷积层的输出特征向量
Figure BDA00033800084000000411
可以简述为：
Figure BDA00033800084000000412
步骤4-3，第三层，池化层：输入为第二层输出的特征向量
Figure BDA00033800084000000413
采用最大池化操作，池化区域宽度为K，则池化操作产生过程如下所示：
Figure BDA00033800084000000414
Figure BDA00033800084000000415
其中，
Figure BDA00033800084000000416
为输入
Figure BDA00033800084000000417
的池化层输出，具有Np维度。
对于输入
Figure BDA00033800084000000418
整个池化层的输出
Figure BDA00033800084000000419
可以简化表述为：
Figure BDA00033800084000000420
进一步，所述步骤5中，通过双向长短时记忆神经网络，捕捉三相定子电流空间特征向量的时序信息，对于输入
Figure BDA00033800084000000421
在t∈[1,2,…,Np]时刻，长短时记忆神经网络内部输入门、遗忘门、输出门计算方法可以简化描述为:
Figure BDA0003380008400000051
其中，it、ft、ot分别为t时刻输入门、遗忘门和输出门向量，
Figure BDA0003380008400000052
为t时刻的输入向量，ht为t时刻隐藏层的输出向量，ht-1为t-1时刻隐藏层状态的输出向量，ct为t时刻记忆细胞状态，σ表示logistic sigmoid函数，Wzi、Whi和Wci为输入门的权重矩阵，Wzf、Whf、Wcf为遗忘门的权重矩阵，Wzc和Whc为记忆细胞的权重矩阵，Wzo、Who和Wco为输出门的权重矩阵，bi、bf、bc和bo分别为输入门、遗忘门、记忆细胞和输出门的偏置向量，
Figure BDA0003380008400000053
进一步，所述步骤5中，所述长短时记忆神经网络输出的特征向量为
Figure BDA0003380008400000054
Figure BDA0003380008400000055
L/2为长短时记忆神经网络的隐藏层神经元数量，
Figure BDA0003380008400000056
为t时刻前向长短时记忆神经网络隐藏层的输出特征向量，
Figure BDA0003380008400000057
为t时刻后向长短时记忆神经网络隐藏层的输出特征向量，双向长短时记忆神经网络在t时刻根据
Figure BDA0003380008400000058
和
Figure BDA0003380008400000059
进行矩阵拼接并输出，双向长短时记忆神经网络在t时刻叠加后的隐藏层输出
Figure BDA00033800084000000510
可以表述为：
Figure BDA00033800084000000511
Figure BDA00033800084000000512
Figure BDA00033800084000000513
其中，
Figure BDA00033800084000000514
和
Figure BDA00033800084000000515
为输入层-隐藏层的权重矩阵，
Figure BDA00033800084000000516
和
Figure BDA00033800084000000517
表示隐藏层-隐藏层的权重矩阵，
Figure BDA00033800084000000518
和
Figure BDA00033800084000000519
表示隐藏层-输出层的权重矩阵，
Figure BDA00033800084000000520
和
Figure BDA00033800084000000521
表示隐藏层偏置向量，by为输出层偏置向量，H为隐藏层激活函数。
对于输入
Figure BDA00033800084000000522
双向长短时记忆网络的输出
Figure BDA00033800084000000523
可以简述为：
Figure BDA0003380008400000061
进一步，所述步骤6的具体过程为：
步骤6-1，展平层：将提取的定子电流数据时空特征向量
Figure BDA0003380008400000062
进一步输入至全连接层实现“展平”操作，将所有的时空特征向量
Figure BDA0003380008400000063
首尾连接，整合为一维向量
Figure BDA0003380008400000064
Nz＝L×FN，同时采用随机失活(dropout)技巧操作进行处理，降低神经元减的相互依赖性，假设神经元失活概率为p，则实现过程可以简述为：
r5～Bernoulli(p) (13)
Figure BDA0003380008400000065
其中，r5是采用Bernoulli函数以概率p，随机生成的0、1向量，
Figure BDA0003380008400000066
为随机失活后的一维特征向量。
步骤6-2，全连接层：输入为展平后随机失活后的一维特征向量
Figure BDA0003380008400000067
使用ReLU激活函数进行处理，神经元个数为Nd，则输出
Figure BDA0003380008400000068
可以描述为：
Figure BDA0003380008400000069
其中，
Figure BDA00033800084000000610
和
Figure BDA00033800084000000611
为全连接层的权重矩阵和偏置向量。
步骤6-3，输出层：输入为上一层的输出
Figure BDA00033800084000000612
神经元数量与轴承故障类别数目保持一致，通过Softmax函数，直接完成轴承故障类型识别
Figure BDA00033800084000000613
其中，
Figure BDA00033800084000000614
和
Figure BDA00033800084000000615
为输出层的权重矩阵和偏置向量。
步骤6-4，设置学习率lr，迭代次数变量为Num，损失函数为交叉熵损失函数，对网络模型进行训练，当损失函数小于0.01或者等于迭代次数时，保存此时的模型参数，得到训练好的具有交流电机轴承故障状态识别能力的诊断模型。
本发明提出的一种采用卷积神经网络和双向长短时记忆网络的交流电机轴承故障诊断方法，利用交流电机定子电流信号作为故障信号，结合卷积神经网络和双向长短时记忆神经网络特性，在特征提取的深度和复杂程度上都具有良好的性能，达到滚动轴承故障特征的精准有效提取，实现了交流电机滚动轴承运行状态的非侵入式监测。具体表现在：
1.本发明采用交流电机定子电流信号替代传统的振动信号作为故障信号，不仅可以降低监控系统建设成本，而且有利于构建故障诊断与驱动控制为一体的电驱动系统，促进电驱动系统向智能执行机构转变。
2.本发明充分利用卷积神经网路和双向长短时记忆神经网络特性，提取定子电流信号中有关滚动轴承故障的空间和时序信息，实现故障特征的精准有效提取，显著提高故障识别精度。
3.本发明所提供的交流电机轴承故障诊断方法，仅需采集电气故障信号，并将特征提取和故障识别融合在一起，实现了端到端的滚动轴承运行状态监测。
附图说明
图1为本发明所述采用卷积神经网络和双向长短时记忆网络的交流电机轴承故障诊断方法流程图。
图2为图1中卷积神经网络和双向长短时记忆网络模型结构示意图。
图3为图2中长短时记忆网络的单元结构框图。
图4为不同高斯噪音背景下卷积神经网络和双向长短时记忆网络训练准确率和损失曲线。(a)无白噪声(b)加入-5db白噪声。
图5为不同高斯噪音背景下卷积神经网络训练准确率和损失曲线。(a)无白噪声(b)加入-5db白噪声。
图6为T-SNE特征可视化分布。(a)原始电流信号输入(b)最大池化层输出(c)双向长短时记忆网络层输出(d)分类层输出。
具体实施方式
本发明公开了一种采用卷积神经网络和双向长短时记忆网络的交流电机轴承故障诊断方法，包括如下步骤：给定交流电机供电频率和负载转矩，设置采样频率fs和采样时间T，以一定时间间隔从交流电机的数据采集系统中获取三相定子电流信号(ia,ib,ic)；以定子电流信号相邻波谷对应的时间为一个时间周期，截取原始定子电流信号中连续m个时间周期的数据点作为一个样本数据
Figure BDA0003380008400000081
Figure BDA0003380008400000082
采用one-hot编码，将无故障交流电机对应的定子电流数据和轴承故障交流电机对应的定子电流数据按照故障类型制作标签，同时以6:2:2的数据比例将定子电流数据分为训练集、验证集和测试集，作为交流电机轴承故障诊断模型的数据集；构建卷积神经网络，包括两层卷积层和一层池化层，通过卷积池化操作提取定子电流信号中有区分性的空间特征信息；构建双向长短时记忆神经网络：将提取的定子电流数据空间信息输入双向长短时记忆网络，捕获定子电流数据空间特征的时序信息；将提取的定子电流数据时空信息输入到全连接层和输出层，选择合适的损失函数、学习率和训练迭代次数，对模型进行训练，得到训练好的具有交流电机轴承故障状态识别能力的诊断模型。
下面结合附图，对本发明实施实例中的技术方案做进一步的详细说明。
具体实施步骤1：
步骤1-1，给定交流电机供电频率和负载转矩，设置采样频率fs和采样时间T，以一定时间间隔从交流电机的数据采集系统中获取三相定子电流信号(ia,ib,ic)；
步骤1-2，以定子电流相邻波谷对应时间为一个时间周期，截取原始定子电流中连续m个时间周期的数据点作为一个样本数据
Figure BDA0003380008400000083
Figure BDA0003380008400000084
Nsub为m个时间周期内单相电流数据的长度，将三相电流视为三个通道作为网络的输入，交流电机轴承故障诊断模型的输入数据样本可表示为
Figure BDA0003380008400000085
Figure BDA0003380008400000086
p＝1,2,…Nseg，Nseg为输入数据样本长度；
步骤1-3，采用one-hot编码，将无故障交流电机对应的定子电流数据和轴承故障交流电机对应的定子电流数据按照故障类型制作标签，同时以6:2:2的数据比例将定子电流数据分为训练集、验证集和测试集，作为交流电机轴承故障诊断模型的数据集；
具体实施步骤2：
步骤2-1，第一层卷积层，输入数据样本Ip，卷积核的数量为FN，长度为FL，步长为1，则卷积层对输入Ip在第j个卷积核处的值
Figure BDA0003380008400000091
可通过以下公式得到：
Figure BDA0003380008400000092
Figure BDA0003380008400000093
其中，Hjc为第j个卷积核第c通道，
Figure BDA0003380008400000094
*T为矩阵*的转置，
Figure BDA0003380008400000095
为第j个卷积核第q个卷积区域值，
Figure BDA0003380008400000096
为第j个卷积核在通道c的偏移向量，fjc为第j个卷积核第c通道的激活函数，q:q+FL-1为第j个卷积核第c通道第q个卷积区域，
Figure BDA0003380008400000097
为第j个卷积核第c通道第q个卷积区域中的元素。
输入Ip在第一层卷积层的输出特征向量
Figure BDA0003380008400000098
数学表达式为：
Figure BDA0003380008400000099
步骤2-2，第二层卷积层：输入第一层输出的特征向量
Figure BDA00033800084000000910
卷积核与第一层卷积层具有相同的参数配置，则第二层卷积层的输出特征向量
Figure BDA00033800084000000911
可以简述为：
Figure BDA00033800084000000912
步骤2-3，第三层池化层：输入为第二层输出的特征向量
Figure BDA00033800084000000913
采用最大池化操作，池化区域宽度为K，则池化操作产生过程如下所示：
Figure BDA00033800084000000914
Figure BDA00033800084000000915
v∈[1,2,…,Np]
其中，
Figure BDA00033800084000000916
为输入
Figure BDA00033800084000000917
的池化层输出，具有Np维度。
对于输入
Figure BDA00033800084000000918
整个池化层的输出
Figure BDA00033800084000000919
可以简化表述为：
Figure BDA00033800084000000920
步骤2-4，第四层双向长短时记忆神经网络层，将提取的定子电流数据空间信息
Figure BDA00033800084000000921
输入双向长短时记忆网络，捕获定子电流数据空间特征的时序信息
Figure BDA00033800084000000922
Figure BDA00033800084000000923
得到双向长短时记忆网络输出
Figure BDA00033800084000000924
Figure BDA00033800084000000925
步骤2-5：第五层展平层，将提取的定子电流数据时空特征向量
Figure BDA0003380008400000101
进一步输入至全连接层实现“展平”操作，将所有的时空特征向量
Figure BDA0003380008400000102
首尾连接，整合为一维向量
Figure BDA0003380008400000103
Nz＝L×FN，同时采用随机失活(dropout)技巧操作进行处理，降低神经元减的相互依赖性，假设神经元失活概率为p，则实现过程可以简述为：
r5～Bernoulli(p)
Figure BDA0003380008400000104
其中，r5是采用Bernoulli函数以概率p，随机生成的0、1向量，
Figure BDA0003380008400000105
为随机失活后的一维特征向量；
步骤2-6，第六层全连接层，输入为展平后随机失活后的一维特征向量
Figure BDA0003380008400000106
使用ReLU激活函数进行处理，神经元个数为Nd，则输出
Figure BDA0003380008400000107
由以下公式得到：
Figure BDA0003380008400000108
其中，
Figure BDA0003380008400000109
和
Figure BDA00033800084000001010
为全连接层的权重矩阵和偏置向量。
步骤2-7，第七层输出层，输入为上一层的输出
Figure BDA00033800084000001011
神经元数量与轴承故障类别数目保持一致，通过Softmax函数，直接完成轴承故障类型识别
Figure BDA00033800084000001012
其中，
Figure BDA00033800084000001013
和
Figure BDA00033800084000001014
为输出层的权重矩阵和偏置向量。
具体实施步骤3：
步骤3-1，设置学习率lr、迭代次数变量Num，损失函数等，以训练集为数据源，对所提轴承故障诊断网络模型进行逐层训练，当损失函数小于0.01或者等于迭代次数是，保存此时的模型参数；
步骤3-2，将验证集数据输入到训练后的轴承故障诊断网络模型中，将输出的结果与验证集数据标签对比计算得到识别准确率，若准确率≥95％，则得到训练好的网络模型，若准确率<95％，则调整网络模型的超参数，再次训练网络模型并进行验证，直到得到训练好的轴承故障诊断网络模型；
步骤3-3，将测试集输入到训练后的轴承故障诊断网络模型进行故障诊断，并统计测试结果；
结合案例分析
实验电机为永磁电机MS-180ST-M19015B-43P0，其所用轴承型号为6209z，轴承外圈故障通过丝锥钻孔而获得。数据采集系统主要有电流传感器、NI数据采集卡以及存储终端组成。采用的电流传感器型号为JLB-21霍尔电流传感器，采集卡为NI USB-6366型数据采集卡，通过在联想笔记本电脑上安装SignalExpress软件完成数据采集。网络训练和测试是在Windows10操作系统、英特尔酷睿i9-10900KCPU、英伟达RTX-4000GPU上进行。其中，电机供电频率为20Hz，负载扭矩为0Nm，三相定子电流信号采样频率为20kHz。表1显示了不同轴承状态的分类编号。
表1
轴承状态 故障尺寸(mm) 标签
无故障轴承 0 C0
外圈故障轴承1 4 C1
外圈故障轴承2 7 C2
在创建数据样本阶段，使用256个连续的三相定子电流数据点组成一个样本，一共生成了120个样本，训练集72个，验证集24个，测试集24个；在搭建网络模型中，第一层和第二次卷积核的数量FN都为16，卷积核的大小FL都为3，卷积步长分别为1和2，第三层最大池化层的池化区域宽度K为2，第四层长短记忆网络的隐含层单元数量为16，第五层随机失活概率p为0.5，第六层全连接层的输出为64，第七层输出层的输出为3，分别代表三种轴承状态。所提网络框架的其他主要参数如表2所示。
表2
操作 卷积核 步长 通道数 激活函数 输出尺寸
输入层 -- -- -- -- 256×3
卷积层1 3*1 1 16 Relu 254×16
卷积层2 3*1 2 16 Relu 126×16
最大池化层 2*1 4 -- -- 32×16
双向长短记忆网络 -- -- 16×2 -- 32×32
展平层 -- -- -- -- 1024
随机失活 -- -- -- -- 512
全连接层 -- -- 64 Relu 64
输出层 -- -- 3 Softmax 3
实验结果
通过添加不同的高斯白噪声，比较了所提网络框架和只使用卷积神经网络，表3显示了不同高斯白噪声情况下的实验结果对比，图4和图5为所提网络框架和只使用卷积神经网络的准确率和损失曲线。这两种方法都是采用三相电流信号进行电机轴承故障诊断。对比实验结果表明，本专利提出的网络框架在不同的高斯白噪声背景下的预测准确率均为100％。
表3
白噪音(SNR) 卷积神经网络-双向长短时记忆网络 卷积神经网络
-5db 100％ 95.83％
无白噪音 100％ 100％
图6显示了所提网络框架在不用网络层的T-SNE特征可视化分布，对比可以看出，随着网络层数的增加，原始电流信号从杂乱无序分布到不同故障郑磊之间的距离更加明显，有效实现了不同轴承故障状态的分类。

Claims (8)

1.一种采用卷积神经网络和双向长短时记忆网络的交流电机轴承故障诊断方法，其特征在于，包括以下步骤：
步骤1，获取交流电机定子电流信号：设定采样频率fs和采样时间T，以一定时间间隔从交流电机的数据采集系统中获取三相定子电流信号，将其作为待处理的数据；
步骤2，定子电流数据样本点划分：以定子电流信号相邻波谷对应的时间为一个时间周期，截取原始定子电流信号中m个时间周期的数据点作为一个样本，将其作为交流电机轴承故障诊断模型的样本数据将其作为交流电机轴承故障诊断模型的样本数据；
步骤3，创建数据集：根据步骤2对交流电机轴承故障类型对应的样本数据进行编码，制作标签，作为交流电机轴承故障诊断模型的数据集；
步骤4，构建卷积神经网络：将训练数据输入卷积神经网络，其包括两层卷积层和一层池化层，通过卷积池化操作提取定子电流信号中有区分性的空间特征信息；
步骤5，构建双向长短时记忆神经网络：将提取的定子电流数据空间信息输入双向长短时记忆网络，捕获定子电流数据空间特征的时序信息；
步骤6，实现滚动轴承故障类型识别：将提取的定子电流数据时空信息输入到全连接层和输出层，选择合适的损失函数、学习率和训练迭代次数，对模型进行训练，得到训练好的具有交流电机轴承故障状态识别能力的诊断模型。
2.根据权利要求1所述的方法，其特征在于，所述步骤1中，在满足奈奎斯特采样定理的条件下，为了包含更多的故障信息，采样频率fs设置为：
fs>10f
其中，f为定子电流信号基频分量。
3.根据权利要求1所述的方法，其特征在于，所述步骤2中，以定子电流信号相邻波谷对应的时间为一个时间周期，截取原始定子电流信号中连续m个时间周期的数据点作为一个样本
Figure FDA0003380008390000011
其中
Figure FDA0003380008390000012
表示连续m个时间周期内a相、b相、c相电流数据，
Figure FDA0003380008390000021
Nsub为m个时间周期内单相电流数据的长度，将三相电流视为三个通道作为网络的输入，交流电机轴承故障诊断模型的输入数据样本可表示为
Figure FDA0003380008390000022
Nseg为输入数据样本长度。
4.根据权利要求1所述的方法，其特征在于，所述步骤3的具体过程为：
步骤3-1，确定标签的编码方式：采用one-hot编码；
步骤3-2，将无故障交流电机对应的定子电流数据和轴承故障交流电机对应的定子电流数据按照故障类型制作标签，同时以6:2:2的数据比例将定子电流数据分为训练集、验证集和测试集，作为交流电机轴承故障诊断模型的数据集；
5.根据权利要求1所述的方法，其特征在于，所述步骤4中，卷积神经网络包括两层卷积层和一层池化层，具体计算过程为：
步骤4-1，第一层：卷积层，输入数据样本Ip∈I，卷积核的数量为FN，长度为FL，步长为1，则卷积层对输入Ip在第j个卷积核处的值
Figure FDA0003380008390000023
可通过以下公式进行计算：
Figure FDA0003380008390000024
Figure FDA0003380008390000025
其中，Hjc为第j个卷积核第c通道，
Figure FDA0003380008390000026
*T为矩阵*的转置，
Figure FDA0003380008390000027
为第j个卷积核第q个卷积区域值，
Figure FDA0003380008390000028
为第j个卷积核在通道c的偏移向量，fjc为第j个卷积核第c通道的激活函数，q:q+FL-1为第j个卷积核第c通道第q个卷积区域，
Figure FDA0003380008390000029
为第j个卷积核第c通道第q个卷积区域中的元素。
输入Ip在第一层卷积层的输出特征向量
Figure FDA00033800083900000210
数学表达式为：
Figure FDA00033800083900000211
步骤4-2，第二层，卷积层：输入第一层输出的特征向量
Figure FDA00033800083900000212
卷积核与第一层卷积层具有相同的数量FN和长度FL，步长为2，则第二层卷积层的输出特征向量
Figure FDA00033800083900000213
可以简述为：
Figure FDA0003380008390000031
步骤4-3，第三层，池化层：输入为第二层输出的特征向量
Figure FDA0003380008390000032
采用最大池化操作，池化区域宽度为K，则池化操作产生过程如下所示：
Figure FDA0003380008390000033
Figure FDA0003380008390000034
其中，
Figure FDA0003380008390000035
为输入
Figure FDA0003380008390000036
的池化层输出，具有Np维度。
对于输入
Figure FDA0003380008390000037
整个池化层的输出
Figure FDA0003380008390000038
可以简化表述为：
Figure FDA0003380008390000039
6.根据权利要求1所述的方法，其特征在于，所述步骤5中，通过双向长短时记忆神经网络，捕捉三相定子电流空间特征向量的时序信息，对于输入
Figure FDA00033800083900000310
在t∈[1,2,···,Np]时刻，长短时记忆神经网络内部输入门、遗忘门、输出门计算方法可以简化描述为:
Figure FDA00033800083900000311
其中，it、ft、ot分别为t时刻输入门、遗忘门和输出门向量，
Figure FDA00033800083900000312
为t时刻的输入向量，ht为t时刻隐藏层的输出向量，ht-1为t-1时刻隐藏层的输出向量，ct为t时刻记忆细胞状态，σ表示logistic sigmoid函数，Wzi、Whi和Wci为输入门的权重矩阵，Wzf、Whf、Wcf为遗忘门的权重矩阵，Wzc和Whc为记忆细胞的权重矩阵，Wzo、Who和Wco为输出门的权重矩阵，bi、bf、bc和bo分别为输入门、遗忘门、记忆细胞和输出门的偏置向量，
Figure FDA00033800083900000313
7.根据权利要求1所述的方法，其特征在于，所述步骤5中，所述长短时记忆神经网络输出的特征向量为
Figure FDA00033800083900000314
Figure FDA00033800083900000315
L/2为长短时记忆神经网络的隐藏层神经元数量，
Figure FDA00033800083900000316
为t时刻前向长短时记忆神经网络隐藏层的输出特征向量，
Figure FDA0003380008390000041
为t时刻后向长短时记忆神经网络隐藏层的输出特征向量，双向长短时记忆神经网络在t时刻根据
Figure FDA0003380008390000042
和
Figure FDA0003380008390000043
进行矩阵拼接并输出，双向长短时记忆神经网络在t时刻叠加后的隐藏层输出
Figure FDA0003380008390000044
可以表述为：
Figure FDA0003380008390000045
Figure FDA0003380008390000046
Figure FDA0003380008390000047
其中，
Figure FDA0003380008390000048
和
Figure FDA0003380008390000049
为输入层-隐藏层的权重矩阵，
Figure FDA00033800083900000410
和
Figure FDA00033800083900000411
表示隐藏层-隐藏层的权重矩阵，
Figure FDA00033800083900000412
和
Figure FDA00033800083900000413
表示隐藏层-输出层的权重矩阵，
Figure FDA00033800083900000414
和
Figure FDA00033800083900000415
表示隐藏层偏置向量，by为输出层偏置向量，H为隐藏层激活函数。
对于输入
Figure FDA00033800083900000416
双向长短时记忆网络的输出
Figure FDA00033800083900000417
可以简述为：
Figure FDA00033800083900000418
8.根据权利要求1所述的方法，其特征在于，所述步骤6的具体过程为：
步骤6-1，展平层：将提取的定子电流数据时空特征向量
Figure FDA00033800083900000419
进一步输入至全连接层实现“展平”操作，将所有的时空特征向量
Figure FDA00033800083900000420
首尾连接，整合为一维向量
Figure FDA00033800083900000421
Nz＝L×FN，同时采用随机失活(dropout)技巧操作进行处理，降低神经元减的相互依赖性，假设神经元失活概率为p，则实现过程可以简述为：
r5～Bernoulli(p) (13)
Figure FDA00033800083900000422
其中，r5是采用Bernoulli函数以概率p，随机生成的0、1向量，
Figure FDA00033800083900000423
为随机失活后的一维特征向量；
步骤6-2，全连接层：输入为展平后随机失活后的一维特征向量
Figure FDA00033800083900000424
使用ReLU激活函数进行处理，神经元个数为Nd，则输出
Figure FDA00033800083900000425
可以描述为：
Figure FDA0003380008390000051
其中，
Figure FDA0003380008390000052
和
Figure FDA0003380008390000053
为全连接层的权重矩阵和偏置向量；
步骤6-3，输出层：输入为上一层的输出
Figure FDA0003380008390000054
神经元数量与轴承故障类别数目保持一致，通过Softmax函数，直接完成轴承故障类型识别
Figure FDA0003380008390000055
其中，
Figure FDA0003380008390000056
和
Figure FDA0003380008390000057
为输出层的权重矩阵和偏置向量；
步骤6-4，设置学习率lr，迭代次数变量为Num，损失函数为交叉熵损失函数，对网络模型进行训练，当损失函数小于0.01或者等于迭代次数时，保存此时的模型参数，得到训练好的具有交流电机轴承故障状态识别能力的诊断模型。
CN202111430349.5A 2021-11-29 2021-11-29 一种采用卷积神经网络和双向长短时记忆网络的交流电机轴承故障诊断方法 Pending CN114201989A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111430349.5A CN114201989A (zh) 2021-11-29 2021-11-29 一种采用卷积神经网络和双向长短时记忆网络的交流电机轴承故障诊断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111430349.5A CN114201989A (zh) 2021-11-29 2021-11-29 一种采用卷积神经网络和双向长短时记忆网络的交流电机轴承故障诊断方法

Publications (1)

Publication Number Publication Date
CN114201989A true CN114201989A (zh) 2022-03-18

Family

ID=80649333

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111430349.5A Pending CN114201989A (zh) 2021-11-29 2021-11-29 一种采用卷积神经网络和双向长短时记忆网络的交流电机轴承故障诊断方法

Country Status (1)

Country Link
CN (1) CN114201989A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114893428A (zh) * 2022-04-02 2022-08-12 上海大学 一种通风机精细预诊与精准运维装置及系统
CN115034268A (zh) * 2022-06-17 2022-09-09 辽宁大学 一种双向lstm融合多尺度卷积的轴承亚健康识别算法
CN115542172A (zh) * 2022-12-01 2022-12-30 湖北工业大学 动力电池故障检测方法、系统、装置及存储介质
CN117031107A (zh) * 2023-08-09 2023-11-10 江苏大学 一种采用主成分分析的交流电机无速度传感器状态监测方法
CN117792190A (zh) * 2023-12-27 2024-03-29 浙江通势达电动技术有限公司 一种高精度电机驱动方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160266208A1 (en) * 2015-03-09 2016-09-15 Rolls-Royce Plc Fault detection and diagnosis in an induction motor
CN110954326A (zh) * 2019-12-17 2020-04-03 北京化工大学 一种自动学习特征表达的滚动轴承在线故障诊断方法
CN112162196A (zh) * 2020-08-27 2021-01-01 浙江工业大学 一种基于图注意力网络的电机故障诊断方法
CN113176092A (zh) * 2021-04-25 2021-07-27 江苏科技大学 基于数据融合和改进经验小波变换的电机轴承故障诊断方法
CN113255209A (zh) * 2021-04-27 2021-08-13 燕山大学 一种齿轮箱轴承剩余寿命的预测方法
CN113339204A (zh) * 2021-05-01 2021-09-03 北京工业大学 一种基于混合神经网络的风力发电机故障识别方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160266208A1 (en) * 2015-03-09 2016-09-15 Rolls-Royce Plc Fault detection and diagnosis in an induction motor
CN110954326A (zh) * 2019-12-17 2020-04-03 北京化工大学 一种自动学习特征表达的滚动轴承在线故障诊断方法
CN112162196A (zh) * 2020-08-27 2021-01-01 浙江工业大学 一种基于图注意力网络的电机故障诊断方法
CN113176092A (zh) * 2021-04-25 2021-07-27 江苏科技大学 基于数据融合和改进经验小波变换的电机轴承故障诊断方法
CN113255209A (zh) * 2021-04-27 2021-08-13 燕山大学 一种齿轮箱轴承剩余寿命的预测方法
CN113339204A (zh) * 2021-05-01 2021-09-03 北京工业大学 一种基于混合神经网络的风力发电机故障识别方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
吴楠 等: "基于一维CNN 与BiLSTM 结合的轴承故障诊断", 《组合机床与自动化加工技术》, no. 9, pages 39 - 45 *
闫书豪 等: "基于一维WConv－BiLSTM的轴承故障诊断算法", 《电子科技》, vol. 34, no. 4, pages 75 - 81 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114893428A (zh) * 2022-04-02 2022-08-12 上海大学 一种通风机精细预诊与精准运维装置及系统
CN114893428B (zh) * 2022-04-02 2023-02-21 上海大学 一种通风机精细预诊与精准运维装置及系统
CN115034268A (zh) * 2022-06-17 2022-09-09 辽宁大学 一种双向lstm融合多尺度卷积的轴承亚健康识别算法
CN115542172A (zh) * 2022-12-01 2022-12-30 湖北工业大学 动力电池故障检测方法、系统、装置及存储介质
CN117031107A (zh) * 2023-08-09 2023-11-10 江苏大学 一种采用主成分分析的交流电机无速度传感器状态监测方法
CN117031107B (zh) * 2023-08-09 2024-03-19 江苏大学 一种采用主成分分析的交流电机无速度传感器状态监测方法
CN117792190A (zh) * 2023-12-27 2024-03-29 浙江通势达电动技术有限公司 一种高精度电机驱动方法

Similar Documents

Publication Publication Date Title
CN114201989A (zh) 一种采用卷积神经网络和双向长短时记忆网络的交流电机轴承故障诊断方法
CN108830127A (zh) 一种基于深度卷积神经网络结构的旋转机械故障特征智能诊断方法
CN110161343A (zh) 一种智能列车外部受电装置的非侵入式实时动态监测方法
CN111751133B (zh) 一种基于非局部均值嵌入的深度卷积神经网络模型的智能故障诊断方法
Liu et al. Intelligent fault diagnosis under small sample size conditions via Bidirectional InfoMax GAN with unsupervised representation learning
CN110657984A (zh) 一种基于强化胶囊网络的行星齿轮箱故障诊断方法
Liu et al. Multi-feature fusion for fault diagnosis of rotating machinery based on convolutional neural network
CN111538947A (zh) 风力发电机轴承故障分类模型的构建方法
CN111783684A (zh) 一种基于Bi-GRU-SA的家庭用电负荷分解方法及系统
CN115905976B (zh) 基于注意力机制的Highway Bi-LSTM轴承故障诊断方法、系统及设备
CN115314287B (zh) 一种基于深度聚类的对抗异常检测系统
Mohamed et al. Induction motor broken rotor bar faults diagnosis using ANFIS-based DWT
Poh et al. LSTM and HMM comparison for home activity anomaly detection
Zhao et al. A novel deep fuzzy clustering neural network model and its application in rolling bearing fault recognition
CN111310719B (zh) 一种未知辐射源个体识别及检测的方法
CN117272230A (zh) 一种基于多任务学习模型的非侵入式负荷监测方法和系统
Li et al. A novel unsupervised anomaly detection method for rotating machinery based on memory augmented temporal convolutional autoencoder
CN111795819A (zh) 一种融合振动与电流信号协同学习的齿轮箱故障诊断方法
Xu et al. Loop optimization noise-reduced LSTM based classifier for PD detection
Sun et al. Intelligent fault diagnosis scheme for rolling bearing based on domain adaptation in one dimensional feature matching
CN111031064A (zh) 一种检测电网假数据注入攻击的方法
CN114781458A (zh) 基于自适应小波和注意力机制的时频记忆神经网络的配电网初期故障识别方法
Rauscher et al. Deep learning and data augmentation for partial discharge detection in electrical machines
CN114462306A (zh) 基于变权重时域卷积网络的非侵入式用电负荷分解方法
Xuejiao et al. A method for rolling bearing life state recognition by combining health indicator and anti-noise deep residual shrinkage network

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20220318