CN108960328A - 一种基于汉克尔矩阵的诊断汽车轮毂轴承故障的方法 - Google Patents
一种基于汉克尔矩阵的诊断汽车轮毂轴承故障的方法 Download PDFInfo
- Publication number
- CN108960328A CN108960328A CN201810724155.8A CN201810724155A CN108960328A CN 108960328 A CN108960328 A CN 108960328A CN 201810724155 A CN201810724155 A CN 201810724155A CN 108960328 A CN108960328 A CN 108960328A
- Authority
- CN
- China
- Prior art keywords
- fault
- sampling
- fault mode
- neural networks
- convolutional 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.)
- Withdrawn
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/29—Graphical models, e.g. Bayesian networks
- G06F18/295—Markov models or related models, e.g. semi-Markov models; Markov random fields; Networks embedding Markov models
-
- 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
-
- 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/2415—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on parametric or probabilistic models, e.g. based on likelihood ratio or false acceptance rate versus a false rejection rate
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Probability & Statistics with Applications (AREA)
- Character Discrimination (AREA)
Abstract
本发明涉及一种基于汉克尔矩阵的诊断汽车轮毂轴承故障的方法,其方案为:首先,通过构造汉克尔矩阵实现对原始振动信号的二维矩阵表示,提高信号的相关性,其次,建立已知故障模式的卷积神经网络增强隐马尔科夫模型;最后,利用该模型诊断未知故障类型的故障模式。卷积神经网络是一种数据驱动的特征学习方法,可以从原始数据中利用卷积神经网络模型对汉克尔矩阵表示的二维信号进行卷积和下采样操作,完整地保留表征信号特征的部分,并降低信号中的高维干扰成分,解决故障信号的特征自动学习问题;最后利用自学习得到的特征建立已知故障模式的隐马尔科夫模型。本发明方法一方面利用卷积神经网络自动学习特征,降低数据维度,以获得数据的分布式特征表示;另一方面,利用隐马尔科夫模型动态序列建模能力和时序模式分类能力判断故障类型。
Description
技术领域
本发明属于轴承故障诊断领域,具体是指一种基于汉克尔矩阵的诊断汽车轮毂轴承故障的方法,特别是一种基于汉克尔矩阵的卷积神经网络增强隐马尔科夫模型诊断汽车轮毂轴承故障的方法。
背景技术
汽车轮毂轴承是汽车底盘***的关键零部件之一,具有承载车体重量和引导轮毂转向的作用。汽车轮毂轴承通常安装在悬架与轮毂之间,用于连接轮毂和轮胎,其运行状态直接影响汽车行驶时的舒适和安全。一旦发生故障,就会引发振动加剧和噪声增大,导致行车过程中强烈的颠簸感和尖锐的噪声,存在巨大的安全隐患。因此,对汽车轮毂轴承的工作状态进行实时监测,准确诊断出其故障,对于汽车行驶安全具有重要意义。
汽车轮毂轴承发生故障时,其振动信号往往会表现出周期性的冲击特征。然而,由于故障引发的异响(通常表现为“嗡嗡声”、“吱吱声”和“咯噔声”等),以及环境噪声的污染,其故障特征提取困难。因此,迫切需要寻找一种新的汽车轮毂轴承诊断方法,解决复杂、微弱故障的特征提取问题,实现轴承故障的可靠识别。
卷积神经网络(Convoutional neural network,CNN),本质上是一种人工神经网络(Artificial neural network,ANN),它可以实现原始振动信号的分布式特征表达,通过逐层构建多层神经网络,并利用卷积(convolution)和下采样(sub-sampling)操作,保留信号中的特征成分,减少高维干扰,最后由全连接层重新生成数据特征,从而实现特征自动学习。隐马尔可夫模型(Hidden Markov Model,HMM),凭借其动态过程时间序列建模能力,广泛应用于语音识别、计算金融、时序分析以及机器翻译等领域。在利用隐马尔可夫模型进行故障模式分类时,振动信号即为模型的观测序列,通过振动信号分割成多个样本输入到模型中去,完成训练过程。然而,原始振动信号维度大、信息分散,直接影响到故障诊断的成功与否。
发明内容
本发明的目的是为了克服现有技术存在的缺点和不足,而提供一种基于汉克尔矩阵的诊断汽车轮毂轴承故障的方法,该方法基于汉克尔矩阵的卷积神经网络增强隐马尔科夫模型来诊断汽车轮毂轴承故障,取得优秀的多故障分类结果。
为实现上述目的,本发明的技术方案是包括有:
S1:通过构造汉克尔矩阵实现对原始振动信号的二维矩阵表示,提高信号的相关性,
S2:建立已知故障模式的卷积神经网络增强隐马尔科夫模型;
S3:利用该模型诊断未知故障类型的故障模式,确定故障类型,
最终实现故障诊断的目的。
进一步设置是所述的步骤S1为:对各类已知的故障类型进行故障模式预定义:建立涵盖各种故障类型的一个完备事件组并定义故障模式,收集所述各类已知故障类型的原始振动信号,该原始振动信号为一维信号,其序列xs=[x1,x2…,xN];其汉克尔矩阵H表示为:
式(1)中,n是嵌入维数,N-n+1是每一个子序列的数据长度。
进一步设置是所述的步骤S2为:
a.基于汉克尔矩阵H,将其作为卷积神经网络的输入图片即:将其输入卷积层可得:
式(2)中上标S0为对H进行第0次下采样操作得到的特征图;为经过第一次卷积操作后得到的第k张特征图;为第一个卷积层中的第k个卷积核;为第一个卷积层中的第k个偏置,f(·)为输出的激活函数,得到经过第一次卷积操作之后的特征表示;
b.经过卷积运算后,得到了原始图片的特征表示,继而对各张卷积后的特征图采用下采样操作降低数据的维度和计算复杂度,得到:
式(3)中,表示的是经过第一次下采样操作后得到的k张特征图片,下采样操作后,特征图片数与上层一致;是该下采样层的取用的步长;为该下采样层的第k个偏置;down(·)为下采样函数,由此,得到经过第一次下采样操作之后的特征表示;
c.继续进行第二次卷积操作,其输入图片为第一次下采样操作得到多张特征图片,有:
式(4)中,Mk表示包含在第一下采样层中的特征图片集合;表示包含在Mk集合中的第i张特征图片;表示对包含在Mk中的特征图片求卷积时的第k个卷积核函数;表示在第二卷积层中的第k个偏置;
同样地,第二次下采样操作后,得到:
式(5)中,表示的是经过第二次卷积操作后得到的k张特征图片,下采样操作后,特征图片数与上层一致;是该下采样层的取用的步长;为该下采样层的第k个偏置;down(·)为下采样函数;
d.全连接层将二维矩阵连接为一维向量,将其作为隐马尔科夫模型的输入,建立卷积神经网络增强的隐马尔科夫模型;
e.对于一个隐状态序列为S={S1,S2,…,Sn,…,SN}和观测序列O={O1,O2,…,Ot,…,OT}的隐马尔科夫模型,定义如下:
h={π,A,B} (6)
式(6)中,π={π1,π2,…,πn,…,πN})为隐状态的初始概率向量,πn=P(c1=Sn),(1≤n≤N c1为t=1时刻,隐状态为Sn的概率;A是隐状态间的状态转移概率矩阵;B是观测值概率矩阵,即隐状态为Sn时观测值为Ot的概率,采用混合高斯模型来拟合各个隐状态下的观测值概率密度函数:
式(7)中,gmm(·)是M个混合高斯模型的加权和,为高斯分布的概率密度函数,μm,σm分别表示第m个高斯模型的均值和方差,权重wm需满足:
对于一个给定的隐马尔科夫模型,状态转移矩阵A和发射矩阵B也都是已知的,每个隐状态Sn由隐马尔科夫模型中M个高斯分布分配不同的权重来构建;
最后,将卷积神经网络训练全连接层输出的一维向量作为观测序列输入隐马尔科夫模型中训练,可以得到已知故障模式的卷积神经网络增强的隐马尔科夫模型。
进一步设置是所述的步骤S3为:对于需要监测的汽车轮毂轴承未知故障类型的振动信号利用式(1)构建二维矩阵,将其作为测试样本,输入至已知汽车轮毂轴承故障模式的卷积神经网络增强隐马尔科夫模型中,确定故障类型,最终诊断出轴承故障。
本申请的创新原理和有益效果是:本申请提出了一种基于汉克尔矩阵的卷积神经网络增强隐马尔科夫模型诊断汽车轮毂轴承故障的方法。首先,通过构造汉克尔矩阵实现对原始振动信号的二维矩阵表示,提高信号的相关性,其次,建立已知故障模式的卷积神经网络增强隐马尔科夫模型;最后,利用该模型诊断未知故障类型的故障模式。该方法对汽车轮毂轴承取得了非常满意的多故障分类结果,有关这方面研究,目前尚无报道。
本申请的技术方案可以从原始数据中利用卷积神经网络模型对汉克尔矩阵表示的二维信号进行卷积和下采样操作,完整地保留表征信号特征的部分,并降低信号中的高维干扰成分,解决故障信号的特征自动学习问题;最后利用自学习得到的特征建立已知故障模式的隐马尔科夫模型。本发明方法一方面利用卷积神经网络自动学习特征,降低数据维度,以获得数据的分布式特征表示;另一方面,利用隐马尔科夫模型动态序列建模能力和时序模式分类能力判断故障类型。本申请的具体效果见实施例图文数据。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,根据这些附图获得其他的附图仍属于本发明的范畴。
图1本发明的原理流程图;
图2本发明的卷积神经网络增强隐马尔科夫模型轴承故障诊断流程图;
图3本发明的卷积神经网络结构图;
图4本发明的隐马尔科夫模型图;
图5本发明实施例案例一的汽车轮毂轴承故障信号波形图;
图6本发明实施例案例一的汽车轮毂轴承故障特征学习结果图;
图7本发明实施例案例一的卷积神经网络--隐马尔科夫模型故障分类结果图;
图8本发明实施例案例一的两种模型故障分类对比折线图;
图9本发明实施例案例二的汽车轮毂轴承故障信号波形图;
图10本发明实施例案例二的汽车轮毂轴承故障特征学习结果图;
图11本发明实施例案例二的卷积神经网络--隐马尔科夫模型故障分类结果图;
图12本发明实施例案例二的两种模型故障分类对比折线图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述。
为本发明实施例中,包括:
S1:通过构造汉克尔矩阵实现对原始振动信号的二维矩阵表示,提高信号的相关性;
S2:建立已知故障模式的卷积神经网络增强隐马尔科夫模型;
S3:最后,利用该模型诊断未知故障类型的故障模式,确定故障类型,最终实现故障诊断的目的。诊断流程如图1-2所示,包括以下步骤:
1、构造汉克尔矩阵实现对原始振动信号的二维矩阵表示。
对各类已知的故障类型进行故障模式预定义:建立涵盖各种故障类型的一个完备事件组(需要监测的未知故障类型属于该完备事件组)并定义故障模式。为此,需要收集各类故障类型的原始振动信号(一维信号)。
汉克尔矩阵常用于实现对一维信号的二维矩阵表示,可提高原始一维信号的相关性,常用于在非平稳信号的分解,具有良好的去噪效果及弱信号的提取能力。对于一个原始振动信号(一维信号)序列xs=[x1,x2…,xN],其汉克尔矩阵H表示为:
式(1)中,n是嵌入维数,N-n+1是每一个子序列的数据长度。显然,汉克尔矩阵H的大小就是n×(N-n+1)。本发明方法在汽车轮毂轴承故障诊断实践中发现:嵌入维数n=50效果好。
2、建立已知故障模式的卷积神经网络增强隐马尔科夫模型
(1)利用卷积神经网络对故障信号特征进行自动学习
利用卷积神经网络模型对二维矩阵进行卷积和下采样操作,完整地保留表征信号特征的部分,并降低信号中的高维干扰成分,解决故障信号的特征自动学习问题。
卷积神经网络最初用于处理图像领域,是一种具有特有操作的深度学习模型,即卷积层、下采样层和全连接层,如图3所示。一个卷积神经网络模型具有三大特性:局部视野(local field),下采样(subsampling)和权值共享(weight sharing)。卷积神经网络正是通过对输入图片不断重复卷积、下采样操作,完整地保留表征信号特征的部分,并降低信号中的高维干扰成分。因此,对已知故障模式的信号构建汉克尔矩阵H后,将其作为卷积神经网络的输入图片即:将其输入卷积层可得:
式(2)中,为了方便描述卷积神经网络的卷积、下采样这两个不断重复的操作,给出了上标S0,可以将其理解为对H进行第0次下采样操作得到的特征图;为经过第一次卷积操作后得到的第k张特征图,每张图片经过卷积操作后,都会表征为不同的特征(在该卷积层中表现为多张特征图的形式);为第一个卷积层(也就是经过第一次卷积操作得到的特征图片所在的层)中的第k个卷积核;为第一个卷积层中的第k个偏置。f(·)为输出的激活函数。由此,得到经过第一次卷积操作之后的特征表示。
经过卷积运算后,得到了原始图片的特征表示,继而对各张卷积后的特征图采用下采样操作降低数据的维度和计算复杂度,得到:
式(3)中,表示的是经过第一次下采样操作后得到的k张特征图片,下采样操作后,特征图片数与上层一致;是该下采样层的取用的步长;为该下采样层的第k个偏置;down(·)为下采样函数。由此,得到经过第一次下采样操作之后的特征表示。
继续进行第二次卷积操作,其输入图片为第一次下采样操作得到多张特征图片,有:
式(4)中,Mk表示包含在第一下采样层中的特征图片集合;表示包含在Mk集合中的第i张特征图片;表示对包含在Mk中的特征图片求卷积时的第k个卷积核函数;表示在第二卷积层中的第k个偏置。
同样地,第二次下采样操作后,得到:
式(5)中,表示的是经过第二次卷积操作后得到的k张特征图片,下采样操作后,特征图片数与上层一致;是该下采样层的取用的步长;为该下采样层的第k个偏置;down(·)为下采样函数。
最后,全连接层将二维矩阵连接为一维向量,将其作为隐马尔科夫模型的输入,可建立卷积神经网络增强的隐马尔科夫模型。
(2)建立已知故障模式的卷积神经网络增强隐马尔科夫模型
隐马尔科夫模型是一个双重随机过程,即其隐状态序列间的状态转移和各个隐状态下的观测值都是随机的。对于一个隐状态序列为S={S1,S2,…,Sn,…,SN}和观测序列为O={O1,O2,…,Ot,…,OT}的隐马尔科夫模型,定义如下:
h={π,A,B} (6)
式(6)中,π={π1,π2,…,πn,…,πN})为隐状态的初始概率向量,πn=P(c1=Sn),(1≤n≤N c1为t=1时刻,隐状态为Sn的概率;A是隐状态间的状态转移概率矩阵;B是观测值概率矩阵,即隐状态为Sn时观测值为Ot的概率,可采用混合高斯模型(Gaussian mixturemodel,GMM)来拟合各个隐状态下的观测值概率密度函数:
式(7)中,gmm(·)是M个混合高斯模型的加权和,为高斯分布的概率密度函数,μm,σm分别表示第m个高斯模型的均值和方差,权重wm需满足:
对于一个给定的隐马尔科夫模型,状态转移矩阵A和发射矩阵B也都是已知的。此外,每个隐状态Sn由隐马尔科夫模型中M个高斯分布分配不同的权重来构建。本发明方法在汽车轮毂轴承故障诊断中实际中发现:混合高斯模型M=3效果好。
最后,将卷积神经网络训练全连接层输出的一维向量作为观测序列输入隐马尔科夫模型中训练,可以得到已知故障模式的卷积神经网络增强的隐马尔科夫模型。
3、利用基于汉克尔矩阵的卷积神经网络增强隐马尔科夫模型诊断未知故障类型的故障模式
对于需要监测的汽车轮毂轴承未知故障类型的振动信号利用式(1)构建二维矩阵,将其作为测试样本,输入至已知汽车轮毂轴承故障模式的卷积神经网络增强隐马尔科夫模型中,确定故障类型,最终诊断出轴承故障。
实施例案例一:
本案例中取用模拟某汽车轮毂轴承故障的振动数据,其轴承型号为:NSK/38BWD12,滚动体列数为双列,滚动体个数为30个,轴承的厚度为36mm,轴承的内径为38mm,轴承的外径为72mm。轴承的故障类型为,(1)正常轴承,预定义为故障模式1;(2)0.1mm的内圈故障,预定义为故障模式2;(3)0.1mm的滚动体故障,预定义为故障模式2;(4)0.1mm的外圈故障,预定义为故障模式4;(5)0.3mm的内圈故障,预定义为故障模式5;(6)0.3mm的滚动体故障,预定义为故障模式6;(7)0.3mm的外圈故障,预定义为故障模式7;(8)0.5mm的内圈故障,预定义为故障模式8;(9)0.5mm的滚动体故障,预定义为故障模式9;(10)0.5mm的外圈故障,预定义为故障模式10;(11)0.7mm的内圈故障,预定义为故障模式11;(12)0.7mm的滚动体故障,预定义为故障模式12。轴承转动频率为21.67Hz,采样频率12k Hz。
原始振动信号波形如图5所示:(a)故障模式1;(b)故障模式2;(c)故障模式3;(d)故障模式4;(e)故障模式5;(f)故障模式6;(g)故障模式7;(h)故障模式8;(i)故障模式9;(j)故障模式10;(k)故障模式11;(l)故障模式12。由图可知,每一类故障信号除了振幅方面有一些差异,很难观察到故障特征。因此,有必要寻找一种数据驱动的特征自动学习方法,实现每一类故障模式特征的分布式表达。
为了实现故障数据的自动特征学习,在采样点中选择200000点数作为预处理的原始数据,每段截取100个数据点构造汉克尔矩阵,嵌入维数n=50。因此,输入图片的尺寸是50×50。每种故障模式训练集选取120000(1200×100)个数据点,测试集选取80000(800×100)个数据点,即每种故障模式分别有1200个训练样本和800个测试样本。
通过旋转三维散点图,可以完整全面地观察每一个故障模式的学***均分类准确率。为进一步提高故障模式的分类准确率,利用隐马尔科夫模型对自学习的特征进行模式识别,最终实现故障分类。
图7给出了采用卷积神经网络-隐马尔科夫模型后,训练样本和测试样本的模式识别结果。如图7(a)所示,对于所有的训练样本,分类平均准确率可达到100%。图7(b)所示为测试样本的结果,故障模式3被误判为故障模式6的概率为4%,故障模式4被误判为故障模式10的概率为4%,故障模式5被误判为故障模式9的概率为5%,故障模式6被误判为故障模式10的概率为5%,故障模式8被误判为故障模式10的概率为1%,故障模式9被误判为故障模式4和故障模式10的概率为5%和1%,故障模式11被误判为故障模式8的概率为3%。尽管如此,测试样本的平均分类准确率仍然高达98.125%,比单独的卷积神经网络模型的平均分类准确率要高出很多。
为进一步说明方法的有效性,给出了单独的卷积神经网络和卷积神经网络-隐马尔科夫模型对所有故障模式的分类准确率的对比折现,如图8所示。图8(a)为训练样本的分类准确率对比,可以看出,尽管两个模型的平均分类准确率都在86%以上,但从每一类故障来看,单独的卷积神经网络模型对于各个故障模式的分类准确率参差不齐,例如,故障模式5只有87%,故障模式8为89%,故障模式6和9为90%,故障模式10为94%,故障模式4为95%,故障模式7,12为97%,故障模式11为99%,故障模式1,2为100%。而利用卷积神经网络-隐马尔科夫模型,每一类都达到了100%,说明了该模型在对每一类故障模式的分类准确率上具有更强的稳定性。如图8(b)所示,在测试样本中,卷积神经网络-隐马尔科夫模型的优势更加明显,对每一类故障的分类准确率均高于94%。而卷积神经网络模型对于每一类故障模式的分类准确率则存在较大波动,例如,故障模式8只有86%,故障模式6只有84%,故障模式10只有83%,故障模式5为77%,而故障模式9仅为23%。
实施案例二:
取某汽车传动***轮毂轴承故障数据,其轴承型号为:SKF/DAC30600037,滚动体列数为双列,滚动体个数为28个,轴承的厚度为37mm,轴承的内径为30mm,轴承的外径为60mm。轴承的故障类型为,(1)正常轴承,预定义为故障模式1;(2)0.2mm的内圈故障,预定义为故障模式2;(3)0.2mm的外圈故障,预定义为故障模式3;(4)0.2mm的滚动体故障,预定义为故障模式4;(5)0.3mm的内圈故障,预定义为故障模式5;(6)0.3mm的外圈故障,预定义为故障模式6。轴承转动频率为25Hz,采样频率12k Hz。
原始振动信号波形如图9所示:(a)故障模式1;(b)故障模式2;(c)故障模式3;(d)故障模式4;(e)故障模式5;(f)故障模式6。由图可知,每一类故障信号除了振幅方面有一些差异,很难观察到故障特征。因此,有必要寻找一种数据驱动的特征自动学习方法,实现每一类故障模式特征的分布式表达。
为了实现故障数据的自动特征学习,在采样点中选择240000点数作为预处理的原始数据,每段截取100个数据点构造汉克尔矩阵,嵌入维数n=50。因此,输入图片的尺寸是50×50。每种故障模式训练集选取180000(1800×100)个数据点,测试集选取60000(600×100)个数据点,即每种故障模式分别有1800个训练样本和600个测试样本。
通过旋转三维散点图,可以完整全面地观察每一个故障模式的学***均分类准确率。为进一步提高故障模式的分类准确率,利用隐马尔科夫模型对自学习的特征进行模式识别,最终实现故障分类。
图11给出了采用卷积神经网络-隐马尔科夫模型后,训练样本和测试样本的模式识别结果。如图11(a)所示,对于所有的训练样本,分类平均准确率可达到100%。图11(b)所示为测试样本的结果,故障模式2被误判为故障模式5的概率为1%,故障模式3被误判为故障模式1和故障模式4的概率为1%和5%,故障模式4被误判为故障模式2的概率为3%。尽管如此,测试样本的平均分类准确率仍然高达98.3%,比单独的卷积神经网络模型的平均分类准确率要高出很多。
为进一步说明方法的有效性,给出了单独的卷积神经网络和卷积神经网络-隐马尔科夫模型对所有故障模式的分类准确率的对比折现,如图12所示。图12(a)为训练样本的分类准确率对比,可以看出,尽管两个模型的平均分类准确率都在88%以上,但从每一类故障来看,单独的卷积神经网络模型对于各个故障模式的分类准确率参差不齐,例如,故障模式2只有79%,故障模式4为87%,故障模式6为88%,故障模式3为95%,故障模式5为97%,故障模式1为100%。而利用卷积神经网络-隐马尔科夫模型s模型,每一类都达到了100%,说明了该模型在对每一类故障模式的分类准确率上具有更强的稳定性。如图12(b)所示,在测试样本中,卷积神经网络-隐马尔科夫模型的优势更加明显,对每一类故障的分类准确率均高于94%。而单独的卷积神经网络对于每一类故障模式的分类准确率则存在较大波动,例如,故障模式2只有88%,故障模式4为79%,故障模式6为67%,而故障模式5仅为19%。
实施例不应视为对发明的限制,但任何基于本发明的精神所作的改进,都应在本发明的保护范围之内。
发明的优点是:
本发明方法通过构造汉克尔矩阵实现对原始振动信号的二维矩阵表示,提高信号的相关性,继而一方面利用卷积神经网络自动学习特征,降低数据维度,以获得数据的分布式特征表示;另一方面,利用隐马尔科夫模型动态序列建模能力和时序模式分类能力进行模式识别,最终可靠地检测出轴承故障类型。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,所述的程序可以存储于一计算机可读取存储介质中,所述的存储介质,如ROM/RAM、磁盘、光盘等。
以上所揭露的仅为本发明较佳实施例而已,当然不能以此来限定本发明之权利范围,因此依本发明权利要求所作的等同变化,仍属本发明所涵盖的范围。
Claims (4)
1.一种基于汉克尔矩阵的诊断汽车轮毂轴承故障的方法,其特征在于包括:
S1:通过构造汉克尔矩阵实现对原始振动信号的二维矩阵表示,提高信号的相关性,
S2:建立已知故障模式的卷积神经网络增强隐马尔科夫模型;
S3:利用该模型诊断未知故障类型的故障模式,确定故障类型,最终实现故障诊断的目的。
2.根据权利要求1所述的方法,其特征在于:所述的步骤S1为:对各类已知的故障类型进行故障模式预定义:建立涵盖各种故障类型的一个完备事件组并定义故障模式,收集所述各类已知故障类型的原始振动信号,该原始振动信号为一维信号,其序列xs=[x1,x2…,xN];其汉克尔矩阵H表示为:
式(1)中,n是嵌入维数,N-n+1是每一个子序列的数据长度。
3.根据权利要求2所述的方法,其特征在于:所述的步骤S2为:
a.基于汉克尔矩阵H,将其作为卷积神经网络的输入图片即:将其输入卷积层可得:
式(2)中上标S0为对H进行第0次下采样操作得到的特征图;为经过第一次卷积操作后得到的第k张特征图;为第一个卷积层中的第k个卷积核;为第一个卷积层中的第k个偏置,f(·)为输出的激活函数,得到经过第一次卷积操作之后的特征表示;
b.经过卷积运算后,得到了原始图片的特征表示,继而对各张卷积后的特征图采用下采样操作降低数据的维度和计算复杂度,得到:
式(3)中,表示的是经过第一次下采样操作后得到的k张特征图片,下采样操作后,特征图片数与上层一致;是该下采样层的取用的步长;为该下采样层的第k个偏置;down(·)为下采样函数,由此,得到经过第一次下采样操作之后的特征表示;
c.继续进行第二次卷积操作,其输入图片为第一次下采样操作得到多张特征图片,有:
式(4)中,Mk表示包含在第一下采样层中的特征图片集合;表示包含在Mk集合中的第i张特征图片;表示对包含在Mk中的特征图片求卷积时的第k个卷积核函数;表示在第二卷积层中的第k个偏置;
同样地,第二次下采样操作后,得到:
式(5)中,表示的是经过第二次卷积操作后得到的k张特征图片,下采样操作后,特征图片数与上层一致;是该下采样层的取用的步长;为该下采样层的第k个偏置;down(·)为下采样函数;
d.全连接层将二维矩阵连接为一维向量,将其作为隐马尔科夫模型的输入,建立卷积神经网络增强的隐马尔科夫模型;
e.对于一个隐状态序列为S={S1,S2,…,Sn,…,SN}和观测序列O={O1,O2,…,Ot,…,OT}的隐马尔科夫模型,定义如下:
h={π,A,B} (6)
式(6)中,π={π1,π2,…,πn,…,πN})为隐状态的初始概率向量,πn=P(c1=Sn),(1≤n≤Nc1为t=1时刻,隐状态为Sn的概率;A是隐状态间的状态转移概率矩阵;B是观测值概率矩阵,即隐状态为Sn时观测值为Ot的概率,采用混合高斯模型来拟合各个隐状态下的观测值概率密度函数:
式(7)中,gmm(·)是M个混合高斯模型的加权和,为高斯分布的概率密度函数,μm,σm分别表示第m个高斯模型的均值和方差,权重wm需满足:
对于一个给定的隐马尔科夫模型,状态转移矩阵A和发射矩阵B也都是已知的,每个隐状态Sn由隐马尔科夫模型中M个高斯分布分配不同的权重来构建;
最后,将卷积神经网络训练全连接层输出的一维向量作为观测序列输入隐马尔科夫模型中训练,可以得到已知故障模式的卷积神经网络增强的隐马尔科夫模型。
4.根据权利要求3所述的方法,其特征在于:所述的步骤S3为:
对于需要监测的汽车轮毂轴承未知故障类型的振动信号利用式(1)构建二维矩阵,将其作为测试样本,输入至已知汽车轮毂轴承故障模式的卷积神经网络增强隐马尔科夫模型中,确定故障类型,最终诊断出轴承故障。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810724155.8A CN108960328A (zh) | 2018-07-04 | 2018-07-04 | 一种基于汉克尔矩阵的诊断汽车轮毂轴承故障的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810724155.8A CN108960328A (zh) | 2018-07-04 | 2018-07-04 | 一种基于汉克尔矩阵的诊断汽车轮毂轴承故障的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108960328A true CN108960328A (zh) | 2018-12-07 |
Family
ID=64485662
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810724155.8A Withdrawn CN108960328A (zh) | 2018-07-04 | 2018-07-04 | 一种基于汉克尔矩阵的诊断汽车轮毂轴承故障的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108960328A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110555273A (zh) * | 2019-09-05 | 2019-12-10 | 苏州大学 | 基于隐马尔科夫模型和迁移学习的轴承寿命预测方法 |
CN110647830A (zh) * | 2019-09-12 | 2020-01-03 | 华中科技大学 | 基于卷积神经网络和高斯混合模型的轴承故障诊断方法 |
CN110705581A (zh) * | 2019-07-18 | 2020-01-17 | 浙江运达风电股份有限公司 | 一种基于改进的隐马尔可夫模型的变桨轴承故障识别方法 |
CN111678699A (zh) * | 2020-06-18 | 2020-09-18 | 山东大学 | 一种面向滚动轴承早期故障监测与诊断方法及*** |
CN112561306A (zh) * | 2020-12-11 | 2021-03-26 | 领伟创新智能***(浙江)有限公司 | 一种基于Hankel矩阵的滚动轴承健康状态评价方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5577166A (en) * | 1991-07-25 | 1996-11-19 | Hitachi, Ltd. | Method and apparatus for classifying patterns by use of neural network |
CN105205475A (zh) * | 2015-10-20 | 2015-12-30 | 北京工业大学 | 一种动态手势识别方法 |
CN106649201A (zh) * | 2016-09-27 | 2017-05-10 | 厦门大学 | 一种基于指数信号的范德蒙分解的数据补全方法 |
CN107817106A (zh) * | 2017-10-11 | 2018-03-20 | 温州大学 | 基于贝叶斯残余变换‑奇异值分解‑高斯混合隐马尔科夫模型框架的滚动轴承故障诊断方法 |
-
2018
- 2018-07-04 CN CN201810724155.8A patent/CN108960328A/zh not_active Withdrawn
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5577166A (en) * | 1991-07-25 | 1996-11-19 | Hitachi, Ltd. | Method and apparatus for classifying patterns by use of neural network |
CN105205475A (zh) * | 2015-10-20 | 2015-12-30 | 北京工业大学 | 一种动态手势识别方法 |
CN106649201A (zh) * | 2016-09-27 | 2017-05-10 | 厦门大学 | 一种基于指数信号的范德蒙分解的数据补全方法 |
CN107817106A (zh) * | 2017-10-11 | 2018-03-20 | 温州大学 | 基于贝叶斯残余变换‑奇异值分解‑高斯混合隐马尔科夫模型框架的滚动轴承故障诊断方法 |
Non-Patent Citations (1)
Title |
---|
SHUHUI WANG 等: "Convolutional neural network-based hidden Markov models for rolling element bearing fault identification", 《KNOWLEDGE-BASED SYSTEMS》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110705581A (zh) * | 2019-07-18 | 2020-01-17 | 浙江运达风电股份有限公司 | 一种基于改进的隐马尔可夫模型的变桨轴承故障识别方法 |
CN110555273A (zh) * | 2019-09-05 | 2019-12-10 | 苏州大学 | 基于隐马尔科夫模型和迁移学习的轴承寿命预测方法 |
CN110647830A (zh) * | 2019-09-12 | 2020-01-03 | 华中科技大学 | 基于卷积神经网络和高斯混合模型的轴承故障诊断方法 |
CN110647830B (zh) * | 2019-09-12 | 2021-12-03 | 华中科技大学 | 基于卷积神经网络和高斯混合模型的轴承故障诊断方法 |
CN111678699A (zh) * | 2020-06-18 | 2020-09-18 | 山东大学 | 一种面向滚动轴承早期故障监测与诊断方法及*** |
CN111678699B (zh) * | 2020-06-18 | 2021-06-04 | 山东大学 | 一种面向滚动轴承早期故障监测与诊断方法及*** |
CN112561306A (zh) * | 2020-12-11 | 2021-03-26 | 领伟创新智能***(浙江)有限公司 | 一种基于Hankel矩阵的滚动轴承健康状态评价方法 |
CN112561306B (zh) * | 2020-12-11 | 2023-12-08 | 领伟创新智能***(浙江)有限公司 | 一种基于Hankel矩阵的滚动轴承健康状态评价方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108960328A (zh) | 一种基于汉克尔矩阵的诊断汽车轮毂轴承故障的方法 | |
CN109117987B (zh) | 一种基于深度学习的个性化交通事故风险预测推荐方法 | |
CN104792530B (zh) | 一种基于SDA和Softmax回归的深度学习滚动轴承故障诊断方法 | |
CN110334764B (zh) | 基于集成深度自编码器的旋转机械智能故障诊断方法 | |
CN106895975B (zh) | 基于Stacked SAE深度神经网络的轴承故障诊断方法 | |
CN109781411B (zh) | 一种结合改进稀疏滤波器与kelm的轴承故障诊断方法 | |
CN103868690B (zh) | 基于多种特征提取和选择的滚动轴承状态自动预警方法 | |
CN112906644B (zh) | 基于深度迁移学习的机械故障智能诊断方法 | |
CN103544392B (zh) | 基于深度学习的医学气体识别方法 | |
CN103728551B (zh) | 一种基于级联集成分类器的模拟电路故障诊断方法 | |
CN106323636A (zh) | 栈式稀疏自动编码深度神经网络的机械故障程度特征自适应提取与诊断方法 | |
CN108256556A (zh) | 基于深度信念网络的风力发电机组齿轮箱故障诊断方法 | |
CN106407649A (zh) | 基于时间递归神经网络的微震信号到时自动拾取方法 | |
Shi et al. | Designing a lightweight 1D convolutional neural network with Bayesian optimization for wheel flat detection using carbody accelerations | |
CN112284735B (zh) | 基于一维卷积与动态路由的多传感器滚动轴承故障诊断 | |
CN105679313A (zh) | 一种音频识别报警***及方法 | |
CN106202489A (zh) | 一种基于大数据的农业病虫害智能诊断*** | |
CN108694408B (zh) | 一种基于深度稀疏滤波卷积神经网络的驾驶行为识别方法 | |
CN115758212A (zh) | 一种基于并行网络和迁移学习的机械设备故障诊断方法 | |
CN112329520B (zh) | 基于生成对抗学习的货车轴承故障识别方法 | |
CN115905976B (zh) | 基于注意力机制的Highway Bi-LSTM轴承故障诊断方法、***及设备 | |
CN110991471A (zh) | 一种高速列车牵引***故障诊断方法 | |
CN115221958A (zh) | 一种基于复合权重的机械设备选择性迁移故障诊断方法 | |
CN115356599B (zh) | 一种多模态城市电网故障诊断方法及*** | |
CN117112991A (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 | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20181207 |
|
WW01 | Invention patent application withdrawn after publication |