CN107144430A - 一种基于增量学习的轴承故障诊断方法 - Google Patents
一种基于增量学习的轴承故障诊断方法 Download PDFInfo
- Publication number
- CN107144430A CN107144430A CN201710499001.9A CN201710499001A CN107144430A CN 107144430 A CN107144430 A CN 107144430A CN 201710499001 A CN201710499001 A CN 201710499001A CN 107144430 A CN107144430 A CN 107144430A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- munderover
- mtd
- mtr
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000003745 diagnosis Methods 0.000 title claims abstract description 25
- 238000012549 training Methods 0.000 claims abstract description 39
- 238000012360 testing method Methods 0.000 claims abstract description 28
- 230000001186 cumulative effect Effects 0.000 claims abstract description 17
- 238000009826 distribution Methods 0.000 claims abstract description 17
- 238000000605 extraction Methods 0.000 claims abstract description 5
- 230000006870 function Effects 0.000 claims description 7
- 230000011218 segmentation Effects 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims 1
- 230000000717 retained effect Effects 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 3
- 238000010801 machine learning Methods 0.000 description 3
- 238000013473 artificial intelligence Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000003909 pattern recognition Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012706 support-vector machine Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
Classifications
-
- 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/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
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种基于增量学习的轴承故障诊断方法,通过对轴承单一故障模式下振动信号的采集、特征提取,然后用提取的特征训练一个SVDD模型,得到相应的一系列模型参数值αk_n,n=1,2,…,N;通过每个训练样本到超球体中心的距离dk_n,得到累积概率分布Fk(d);对于轴承测试样本y,计算其到每类故障模式分类超球体中心的距离dk_y,并代入累积概率分布Fk(d)中,计算出概率Pk,然后找到的最大值,其对应的故障模式则是测试样本所属的故障模式。本发明在有新的故障知识出现时,只需对新故障模式进行训练,这样可以保留之前学习的结果并且从不断出现的新数据中学习,可以提高效率和节省计算资源。
Description
技术领域
本发明属于轴承故障诊断技术领域,更为具体地讲,涉及一种基于增量学习的轴承故障诊断方法。
背景技术
轴承是旋转机械中应用最广泛的一种部件,而且作为旋转机械关键部件之一,其工作状态的好坏会影响整个设备的运行状态。轴承一旦发生故障,必将导致旋转机械结构失效,从而带来经济损失,严重还会引发安全事故。因此,对轴承进行故障诊断具有重要的工程意义。
轴承故障诊断通常被看作模式识别问题,并且随着计算机软硬件技术的发展,故障诊断越来越趋向于以人工智能技术为核心的智能诊断。而人工智能技术主要采用基于数据驱动型的机器学习方法,譬如神经网络、支持向量机和k-近邻法等。传统的机器学习采用批量学习方法,即在进行学习之前,必须将所有的数据准备好。而且当一组新数据进来后,往往需要舍弃之前的学习成果,根据新增加的数据和过去所有数据一起重新学习。显然,如果在新数据增加后,放弃以前的学习成果,重新学习全部数据,势必会消耗大量的时间和计算资源,也不适合在线故障监测。
发明内容
本发明的目的在于克服现有技术的不足,从增量学习角度出发,提出一种基于增量学习的轴承故障诊断方法,保留之前学习的结果并且从不断出现的新数据中学习,以提高效率和节省计算资源。
为实现上述发明目的,本发明基于增量学习的轴承故障诊断方法,其特征在于,包括以下步骤:
(1)、轴承振动信号采集
使用振动数据采集仪采集轴承单一故障模式下的振动信号,得到轴承在每一类故障模式下的故障数据dk,其中,k表示故障模式编号;
(2)、振动信号分割以及特征提取
首先对故障数据dk进行分割,得到N个数据段,然后提取每个数据段的特征,并将每个数据段的特征信息作为一个训练样本sk_n,n=1,2,…,N,最后,将所有的训练样本进行归一化处理,所有归一化后的训练样本构成的集合为Sk;
(3)、利用样本集合Sk训练支持向量数据描述(Support Vector DataDescription,简称SVDD)模型,训练SVDD模型等效于求解公式(1),求解结束后可以得到一系列模型参数值αk_n,n=1,2,…,N,即其数目等于训练样本个数N并且一一对应;
其中,κ为给定的核函数,xk_i、xk_j分别为第i、j个训练样本,C是给定的惩罚参数;
(4)、SVDD模型为一个超球体,根据公式(2)计算每个训练样本xk_n到第k类超球体中心的距离dk_n:
(5)、根据距离dk_n,n=1,2,…,N,利用核密度估计方法估计距离的概率密度,并进一步转换得到距离的累积概率分布Fk(d),其中,d表示距离;
(6)、当轴承出现新的故障模式时,重复步骤(1)~(5),得到一个新的累积概率分布,这样故障模式有K个,即有K个累积概率分布Fk(d),k=1,2,…,K;
(7)、轴承的故障诊断过程包含以下步骤:
7.1)、按照步骤(1)和步骤(2)的方法,获取轴承测试样本y;
7.2)、按照步骤(4)的方法,根据公式(3)计算测试样本y距第k类故障模式分类超球体中心的距离dk_y:
7.3)、将距离dk_y作为距离代入第k个累积概率分布Fk(d)中,按照公式(4)计算测试样本属于第k类故障模式的概率Pk:
Pk=1-Fk(dk_y),k=1,2,…,K (4)
7.4)、在K个概率Pk,k=1,2,…,K中找到的最大值,其对应的故障模式则是测试样本所属的故障模式(类别)。
本发明的目的是这样实现的。
本发明基于增量学习的轴承故障诊断方法,通过对轴承单一故障模式下振动信号的采集、特征提取,然后用提取的特征训练一个SVDD模型,得到相应的一系列模型参数值αk_n,n=1,2,…,N,即其数目等于训练样本个数N;通过每个训练样本到超球体中心的距离dk_n,得到累积概率分布Fk(d);对于轴承测试样本y,计算其到每类故障模式分类超球体中心的距离dk_y,并代入累积概率分布Fk(d)中,计算出概率Pk,在K个概率Pk,k=1,2,…,K中找到的最大值,其对应的故障模式则是测试样本所属的故障模式。本发明在有新的故障出现时,只需要对新故障模式的累积概率分布Fk(d)进行训练即可,这样可以保留之前学习的结果并且从不断出现的新数据中学习,类似人类渐进获取知识的过程,可以提高效率和节省计算资源。
附图说明
图1是本发明基于增量学习的轴承故障诊断方法一种具体实施方式流程图;
图2是训练样本获取过程示意图;
图3是轴承诊断模型构建时间;
图4是轴承诊断准确率。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述,以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
图1是本发明基于增量学习的轴承故障诊断方法一种具体实施方式流程图。
在本实施例中,如图1所示,本发明基于增量学习的轴承故障诊断方法基于增量学习的轴承故障诊断方法包括以下步骤:
步骤S1:轴承振动信号采集
对于每一类故障模式,使用振动数据采集仪采集轴承单一故障模式下的振动信号,得到轴承在每一类故障模式下的故障数据dk,其中,k表示故障模式编号。
步骤S2:振动信号特征提取
首先对故障数据dk进行分割,得到N个数据段,然后提取每个数据段的特征,并将每个数据段的特征信息作为一个训练样本sk_n,n=1,2,…,N,最后,将所有的训练样本进行归一化处理,所有归一化后的训练样本构成的集合为Sk。
在本实施例中,振动数据采集仪采用三个加速器对振动信号进行采集,三个加速度分别垂向安装在驱动端、风扇端和支撑台。在本实施例中,步骤S2具体包括:
步骤S21:将每个加速器同时采集的振动数据(故障数据)进行分割,分割长度为t个数据点,在本实施例中,t=100,得到得到N个数据段;
步骤S22:计算每个数据段的均值、方差、均方根和峭度值四个状态特征变量,然后拼接成在一起构成一个特征单元;
步骤S23:将三个加速器采集振动数据相同数据段的特征单元再进行拼接作为数据段的特征信息,然后再将其作为一个训练样本sk_n,n=1,2,…,N。这样得到每个训练样本都包含3个均值、3个方差、3个均方根和3个峭度值共12个特征值{ck_m_n,m=1,2,…,12};
步骤S24:将N个训练样本的12个特征值ck_m_n,按照公式(5)归一到[0,1]区间,所有归一化后的训练样本构成的集合为Sk。
其中,是样本第m个特征归一化后的值,ck_m_max和ck_m_min分别表示N个训练样本第m个特征中的最大值和最小值。
步骤S3:利用样本集合Sk训练支持向量数据描述(Support Vector DataDescription,简称SVDD)模型,训练SVDD模型等效于求解公式(1),求解结束后可以得到一系列模型参数值αk_n,n=1,2,…,N,即其数目等于训练样本个数N并且一一对应;
其中,κ为给定的核函数,xk_i、xk_j分别为第i、j个训练样本,C是给定的惩罚参数。
在本实施例中,所述核函数κ为高斯径向基函数:
其中,γ=1/12。
在本实施例中,所述惩罚参数C取值0.5。
步骤S4:SVDD模型为一个超球体,根据公式(2)计算每个训练样本xk_n到第k类超球体中心的距离dk_n:
步骤S5:根据距离dk_n,n=1,2,…,N,利用核密度估计方法估计距离的概率密度,并进一步转换得到距离的累积概率分布Fk(d),其中,d表示距离;
步骤S6:当轴承出现新的故障模式时,重复步骤步骤S1~S6,得到一个新的累积概率分布,这样故障模式有K个,即有K个累积概率分布Fk(d),k=1,2,…,K;
步骤S7:轴承的故障诊断过程包含以下步骤:
步骤S7.1:按照步骤S1和步骤S2的方法,获取轴承测试样本y;
步骤S7.2:按照步骤S4的方法,根据公式(3)计算测试样本y距第k类故障模式分类超球体中心的距离dk_y:
步骤S7.3:将距离dk_y作为距离代入第k个累积概率分布Fk(d)中,按照公式(4)计算测试样本属于第k类故障模式的概率Pk:
Pk=1-Fk(dk_y),k=1,2,…,K (4)
步骤S7.4:在K个概率Pk,k=1,2,…,K中找到的最大值,其对应的故障模式则是测试样本所属的故障模式(类别)。
在本实施例中,以凯斯西储大学的轴承实验数据为例,具体过程如下:
该轴承故障诊断试验台由电机、转矩传感器/编码器、测力计和控制电子设备等组成。两个测试轴承(分别为驱动端轴承与风扇端轴承)支撑电机轴,并通过电火花加工技术加工了单点故障。故障严重程度(直径)包括0.007英寸、0.014英寸,0.021英寸,故障类型包含内圈故障、外圈故障和滚动体故障。由于轴承外圈固定,外圈故障还包括故障点所处位置。实验中所用轴承型号为6205-2RS JEM SKF。
当电机转速为1797转/分时,由三个加速器(分别垂向安装在驱动端、风扇端和支撑台)采集轴承不同故障类型及故障严重程度下的振动数据用于诊断分析,采样频率为12kHz,每种故障模式采样时间大约10秒。
表1
表1是用于诊断分析的测试轴承故障模式列表。
针对表1中每一类故障模式,样本数为1200,总共25类,一共的样本数为30000个。
本实施例中,将其中每类的960个样本,一共24000个样本作为训练样本,来训练SVDD模型以及构建概率模型。剩余的6000个样本作为测试样本用于诊断过程,以验证本文方法的有效性。同时,实验过程中,为了模拟类别增量过程,从2类样本开始,往后依次增加1类样本,测试时则将这些类的测试集全部测试。比如,刚开始的2类样本得到一个二分类器,测试则将这两类的测试集测试分类器的性能,当第3类样本加进来时,训练则可以得到三分类器,测试则将第3类样本的测试集加上前2类的测试集测试分类器的性能,以此类推。
此外,还将本发明与传统的类别增量算法CIL(class incremental leaning来自论文:ZHANG B F,SU J S,XU X.A class-incremental learning method for multi-class support vector machines in text classification[C]//Proceedings of the5th International Conference on Machine Learning and Cybernetics,Dalian,2006:2581-2585.)和传统的批量学习方法v-SVC进行比较,实验结果见图3和图4。
从轴承故障诊断的实验结果可以看出,本发明在模型构建时间上比其他两种方法低,模型构建效率更高,大概是CIL的10倍,v-SVC的138倍。此外,本发明在准确率方面平均大约比v-SVC低7%,但比同样是类别增量学习方法的CIL高18%。由此可见,本发明应用于轴承故障模式识别,在模型构建效率以及准确率方面也具有一定优点。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (3)
1.一种基于增量学习的轴承故障诊断方法,其特征在于,包括以下步骤:
(1)、轴承振动信号采集
使用振动数据采集仪采集轴承单一故障模式下的振动信号,得到轴承在每一类故障模式下的故障数据dk,其中,k表示故障模式编号;
(2)、振动信号特征提取
首先对故障数据dk进行分割,得到N个数据段,然后提取每个数据段的特征,并将每个数据段的特征信息作为一个训练样本sk_n,n=1,2,…,N,最后,将所有的训练样本进行归一化处理,所有归一化后的训练样本构成的集合为Sk;
(3)、利用样本集合Sk训练支持向量数据描述(Support Vector Data Description,简称SVDD)模型,训练SVDD模型等效于求解公式(1),求解结束后可以得到一系列参数值αk_n,n=1,2,…,N,即其数目等于训练样本个数N并且一一对应;
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>min</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>j</mi>
</mrow>
</msub>
<mi>&kappa;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>i</mi>
</mrow>
</msub>
<mi>&kappa;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mi>C</mi>
<mo>,</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中,xk_i、xk_j分别为第i、j个训练样本,C是给定的惩罚参数;
(4)、SVDD模型为一个超球体,根据公式(2)计算每个训练样本xk_n到第k类超球体中心的距离dk_n:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>d</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>&kappa;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>n</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>2</mn>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>i</mi>
</mrow>
</msub>
<mi>&kappa;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>n</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>c</mi>
<mi>o</mi>
<mi>n</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>n</mi>
<mi>s</mi>
<mi>t</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>j</mi>
</mrow>
</msub>
<mi>&kappa;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
(5)、根据距离dk_n,n=1,2,…,N,利用核密度估计方法估计距离的概率密度,并进一步转换得到距离的累积概率分布Fk(d),其中,d表示距离;
(6)、当轴承出现新的故障模式时,重复步骤(1)~(5),得到一个新的累积概率分布,这样故障模式有K个,即有K个累积概率分布Fk(d),k=1,2,…,K;
(7)、轴承的故障诊断过程包含以下步骤:
7.1)、按照步骤(1)和步骤(2)的方法,获取轴承测试样本y;
7.2)、按照步骤(4)的方法,根据公式(3)计算测试样本y距第k类故障模式分类超球体中心的距离dk_y:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>d</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>y</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>&kappa;</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>2</mn>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>i</mi>
</mrow>
</msub>
<mi>&kappa;</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>c</mi>
<mi>o</mi>
<mi>n</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>n</mi>
<mi>s</mi>
<mi>t</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>i</mi>
</mrow>
</msub>
<mi>&kappa;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>K</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
7.3)、将距离dk_y作为距离代入第k个累积概率分布Fk(d)中,按照公式(4)计算测试样本属于第k类故障模式的概率Pk:
Pk=1-Fk(dk_y),k=1,2,…,K (4)
7.4)、在K个概率Pk,k=1,2,…,K中找到的最大值,其对应的故障模式则是测试样本所属的故障模式(类别)。
2.根据权利要求1所述的轴承故障诊断方法,其特征在于,振动数据采集仪采用三个加速器对振动信号进行采集,三个加速度分别垂向安装在驱动端、风扇端和支撑台;
步骤(2)具体为:
步骤2.1)、将每个加速器同时采集的振动数据(故障数据)进行分割,分割长度为t个数据点,在本实施例中,t=100,得到N个数据段;
步骤2.2)、计算每个数据段的均值、方差、均方根和峭度值四个状态特征变量,然后拼接成在一起构成一个特征单元;
步骤2.3)、将三个加速器采集振动数据相同数据段的特征单元再进行拼接作为数据段的特征信息,然后再将其作为一个训练样本sk_n,n=1,2,…,N。这样得到每个训练样本都包含3个均值、3个方差、3个均方根和3个峭度值共12个特征值{ck_m_n,m=1,2,…,12};
步骤2.4)、将N个训练样本的12个特征值ck_m_n,m=1,2,…,12。按照公式(5)归一到[0,1]区间,所有归一化后的训练样本构成的集合为Sk。
<mrow>
<msub>
<mover>
<mi>c</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>m</mi>
<mo>_</mo>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>c</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>m</mi>
<mo>_</mo>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>c</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>m</mi>
<mo>_</mo>
<mi>min</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>c</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>m</mi>
<mo>_</mo>
<mi>max</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>c</mi>
<mrow>
<mi>k</mi>
<mo>_</mo>
<mi>m</mi>
<mo>_</mo>
<mi>min</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,是样本第m个特征归一化后的值,ck_m_max和ck_m_min分别表示N个训练样本第m个特征中的最大值和最小值。
3.根据权利要求1所述的轴承故障诊断方法,其特征在于:
在本实施例中,所述核函数κ为高斯径向基函数:
<mrow>
<mi>&kappa;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>&gamma;</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,γ=1/12。
在本实施例中,所述惩罚参数C取值0.5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710499001.9A CN107144430B (zh) | 2017-06-27 | 2017-06-27 | 一种基于增量学习的轴承故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710499001.9A CN107144430B (zh) | 2017-06-27 | 2017-06-27 | 一种基于增量学习的轴承故障诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107144430A true CN107144430A (zh) | 2017-09-08 |
CN107144430B CN107144430B (zh) | 2019-02-01 |
Family
ID=59785620
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710499001.9A Expired - Fee Related CN107144430B (zh) | 2017-06-27 | 2017-06-27 | 一种基于增量学习的轴承故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107144430B (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107748209A (zh) * | 2017-12-04 | 2018-03-02 | 中国水利水电科学研究院 | 一种结构损伤检测方法 |
CN107782551A (zh) * | 2017-10-30 | 2018-03-09 | 电子科技大学 | 一种机械部件损伤程度的评估方法 |
CN108229071A (zh) * | 2018-02-09 | 2018-06-29 | 杭州安脉盛智能技术有限公司 | 基于ar模型与svdd算法的刀具性能退化评估方法及系统 |
CN108376264A (zh) * | 2018-02-26 | 2018-08-07 | 上海理工大学 | 一种基于支持向量机增量学习的冷水机组故障诊断方法 |
CN108844612A (zh) * | 2018-08-27 | 2018-11-20 | 重庆大学 | 一种基于数理统计概率模型的变压器内部故障识别方法 |
CN109029699A (zh) * | 2018-06-12 | 2018-12-18 | 国网四川省电力公司乐山供电公司 | 一种变压器振动在线异常检测方法 |
CN109299705A (zh) * | 2018-10-24 | 2019-02-01 | 电子科技大学 | 基于一维深度残差卷积神经网络的旋转机械故障诊断方法 |
CN109359699A (zh) * | 2018-11-02 | 2019-02-19 | 广东工业大学 | 一种基于Alpha稳定分布的工业机组故障诊断方法 |
CN109828168A (zh) * | 2019-01-31 | 2019-05-31 | 福州大学 | 基于核密度估计的变换器故障诊断方法 |
CN109974782A (zh) * | 2019-04-10 | 2019-07-05 | 郑州轻工业学院 | 基于大数据敏感特征优化选取的设备故障预警方法及系统 |
CN110070060A (zh) * | 2019-04-26 | 2019-07-30 | 天津开发区精诺瀚海数据科技有限公司 | 一种轴承设备的故障诊断方法 |
CN110647924A (zh) * | 2019-09-05 | 2020-01-03 | 国网河北省电力有限公司电力科学研究院 | 一种基于支持向量描述与k近邻算法的gis设备状态评估方法 |
CN111401136A (zh) * | 2020-02-22 | 2020-07-10 | 上海交通大学 | 一种柱塞泵空化程度检测方法、装置及终端 |
CN111947927A (zh) * | 2020-07-16 | 2020-11-17 | 南京航空航天大学 | 一种基于色度理论的滚动轴承故障检测方法 |
CN113865873A (zh) * | 2021-12-02 | 2021-12-31 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | 一种轴承组的健康监测方法 |
CN116805157A (zh) * | 2023-08-25 | 2023-09-26 | 中国人民解放军国防科技大学 | 无人集群自主动态评估方法及装置 |
CN116839909A (zh) * | 2023-06-30 | 2023-10-03 | 大连华冶联自动化有限公司 | 一种基于增量学习和规则推理的风机轴承故障诊断方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4730484A (en) * | 1986-06-20 | 1988-03-15 | Robert Olschefski | Missing bearing detector |
CN102072829A (zh) * | 2010-11-04 | 2011-05-25 | 同济大学 | 一种面向钢铁连铸设备的故障预测方法及装置 |
US20130060417A1 (en) * | 2010-05-18 | 2013-03-07 | Navistar Canada, Inc. | Turbo-charger bearing monitor |
CN103645249A (zh) * | 2013-11-27 | 2014-03-19 | 国网黑龙江省电力有限公司 | 基于精简集下采样不均衡svm变压器在线故障检测方法 |
KR101482511B1 (ko) * | 2013-09-04 | 2015-01-19 | 주식회사 우진 | 위상 지연과 데이터 분포 형상지수를 이용한 베어링 결함 진단 시스템 및 그 진단 방법 |
CN104462846A (zh) * | 2014-12-22 | 2015-03-25 | 山东鲁能软件技术有限公司 | 一种基于支持向量机的设备故障智能诊断方法 |
CN104729853A (zh) * | 2015-04-10 | 2015-06-24 | 华东交通大学 | 一种滚动轴承性能退化评估装置及方法 |
CN105528504A (zh) * | 2016-03-01 | 2016-04-27 | 哈尔滨理工大学 | 基于cfoa-mkhsvm的滚动轴承健康状态评估方法 |
CN106323635A (zh) * | 2016-08-04 | 2017-01-11 | 中国航空综合技术研究所 | 一种滚动轴承故障在线检测与状态评估方法 |
-
2017
- 2017-06-27 CN CN201710499001.9A patent/CN107144430B/zh not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4730484A (en) * | 1986-06-20 | 1988-03-15 | Robert Olschefski | Missing bearing detector |
US20130060417A1 (en) * | 2010-05-18 | 2013-03-07 | Navistar Canada, Inc. | Turbo-charger bearing monitor |
CN102072829A (zh) * | 2010-11-04 | 2011-05-25 | 同济大学 | 一种面向钢铁连铸设备的故障预测方法及装置 |
KR101482511B1 (ko) * | 2013-09-04 | 2015-01-19 | 주식회사 우진 | 위상 지연과 데이터 분포 형상지수를 이용한 베어링 결함 진단 시스템 및 그 진단 방법 |
CN103645249A (zh) * | 2013-11-27 | 2014-03-19 | 国网黑龙江省电力有限公司 | 基于精简集下采样不均衡svm变压器在线故障检测方法 |
CN104462846A (zh) * | 2014-12-22 | 2015-03-25 | 山东鲁能软件技术有限公司 | 一种基于支持向量机的设备故障智能诊断方法 |
CN104729853A (zh) * | 2015-04-10 | 2015-06-24 | 华东交通大学 | 一种滚动轴承性能退化评估装置及方法 |
CN105528504A (zh) * | 2016-03-01 | 2016-04-27 | 哈尔滨理工大学 | 基于cfoa-mkhsvm的滚动轴承健康状态评估方法 |
CN106323635A (zh) * | 2016-08-04 | 2017-01-11 | 中国航空综合技术研究所 | 一种滚动轴承故障在线检测与状态评估方法 |
Non-Patent Citations (3)
Title |
---|
TAO XIN-MIN ETC: "A Novel Model of one-class Bearing Fault Detection using SVDD and Genetic Algorithm", 《2007 2ND IEEE CONFERENCE ON INDUSTRIAL ELECTRONICS AND APPLICATIONS》 * |
姚亚夫 等: "基于瞬时能量熵和SVM的滚动轴承故障诊断", 《电子测量与仪器学报》 * |
陈伟 等: "利用小波包和SVDD的分拣机轴承故障诊断", 《振动、测试与诊断》 * |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107782551A (zh) * | 2017-10-30 | 2018-03-09 | 电子科技大学 | 一种机械部件损伤程度的评估方法 |
CN107782551B (zh) * | 2017-10-30 | 2020-10-23 | 电子科技大学 | 一种机械部件损伤程度的评估方法 |
CN107748209B (zh) * | 2017-12-04 | 2020-05-15 | 中国水利水电科学研究院 | 一种结构损伤检测方法 |
CN107748209A (zh) * | 2017-12-04 | 2018-03-02 | 中国水利水电科学研究院 | 一种结构损伤检测方法 |
CN108229071A (zh) * | 2018-02-09 | 2018-06-29 | 杭州安脉盛智能技术有限公司 | 基于ar模型与svdd算法的刀具性能退化评估方法及系统 |
CN108376264A (zh) * | 2018-02-26 | 2018-08-07 | 上海理工大学 | 一种基于支持向量机增量学习的冷水机组故障诊断方法 |
CN109029699A (zh) * | 2018-06-12 | 2018-12-18 | 国网四川省电力公司乐山供电公司 | 一种变压器振动在线异常检测方法 |
CN109029699B (zh) * | 2018-06-12 | 2020-09-11 | 国网四川省电力公司乐山供电公司 | 一种变压器振动在线异常检测方法 |
CN108844612A (zh) * | 2018-08-27 | 2018-11-20 | 重庆大学 | 一种基于数理统计概率模型的变压器内部故障识别方法 |
CN109299705A (zh) * | 2018-10-24 | 2019-02-01 | 电子科技大学 | 基于一维深度残差卷积神经网络的旋转机械故障诊断方法 |
CN109299705B (zh) * | 2018-10-24 | 2021-08-20 | 电子科技大学 | 基于一维深度残差卷积神经网络的旋转机械故障诊断方法 |
CN109359699A (zh) * | 2018-11-02 | 2019-02-19 | 广东工业大学 | 一种基于Alpha稳定分布的工业机组故障诊断方法 |
CN109828168A (zh) * | 2019-01-31 | 2019-05-31 | 福州大学 | 基于核密度估计的变换器故障诊断方法 |
CN109974782A (zh) * | 2019-04-10 | 2019-07-05 | 郑州轻工业学院 | 基于大数据敏感特征优化选取的设备故障预警方法及系统 |
CN109974782B (zh) * | 2019-04-10 | 2021-03-02 | 郑州轻工业学院 | 基于大数据敏感特征优化选取的设备故障预警方法及系统 |
CN110070060B (zh) * | 2019-04-26 | 2021-06-04 | 天津开发区精诺瀚海数据科技有限公司 | 一种轴承设备的故障诊断方法 |
CN110070060A (zh) * | 2019-04-26 | 2019-07-30 | 天津开发区精诺瀚海数据科技有限公司 | 一种轴承设备的故障诊断方法 |
CN110647924A (zh) * | 2019-09-05 | 2020-01-03 | 国网河北省电力有限公司电力科学研究院 | 一种基于支持向量描述与k近邻算法的gis设备状态评估方法 |
CN110647924B (zh) * | 2019-09-05 | 2022-04-22 | 国网河北省电力有限公司电力科学研究院 | 一种基于支持向量描述与k近邻算法的gis设备状态评估方法 |
CN111401136A (zh) * | 2020-02-22 | 2020-07-10 | 上海交通大学 | 一种柱塞泵空化程度检测方法、装置及终端 |
CN111947927A (zh) * | 2020-07-16 | 2020-11-17 | 南京航空航天大学 | 一种基于色度理论的滚动轴承故障检测方法 |
CN113865873A (zh) * | 2021-12-02 | 2021-12-31 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | 一种轴承组的健康监测方法 |
CN113865873B (zh) * | 2021-12-02 | 2022-02-18 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | 一种轴承组的健康监测方法 |
CN116839909A (zh) * | 2023-06-30 | 2023-10-03 | 大连华冶联自动化有限公司 | 一种基于增量学习和规则推理的风机轴承故障诊断方法 |
CN116839909B (zh) * | 2023-06-30 | 2024-12-13 | 大连华冶联自动化有限公司 | 一种基于增量学习和规则推理的风机轴承故障诊断方法 |
CN116805157A (zh) * | 2023-08-25 | 2023-09-26 | 中国人民解放军国防科技大学 | 无人集群自主动态评估方法及装置 |
CN116805157B (zh) * | 2023-08-25 | 2023-11-17 | 中国人民解放军国防科技大学 | 无人集群自主动态评估方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107144430B (zh) | 2019-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107144430A (zh) | 一种基于增量学习的轴承故障诊断方法 | |
CN105910827B (zh) | 基于判别性卷积特征学习的感应电机故障诊断方法 | |
CN110375987B (zh) | 一种基于深度森林机械轴承故障检测方法 | |
CN109297689B (zh) | 一种引入权重因子的大型水力机械智能诊断方法 | |
CN110334764A (zh) | 基于集成深度自编码器的旋转机械智能故障诊断方法 | |
CN110110768A (zh) | 基于并行特征学习和多分类器的滚动轴承故障诊断方法 | |
CN109946080B (zh) | 一种基于嵌入式循环网络的机械设备健康状态识别方法 | |
CN106017876A (zh) | 基于等权局部特征稀疏滤波网络的轮对轴承故障诊断方法 | |
CN110516305A (zh) | 基于注意机制元学习模型的小样本下故障智能诊断方法 | |
CN110672343A (zh) | 基于多注意力卷积神经网络的旋转机械故障诊断方法 | |
CN106441888A (zh) | 一种高速列车滚动轴承故障诊断方法 | |
CN113139278A (zh) | 一种基于深度学习的轴承剩余寿命预测方法及系统 | |
CN110297479A (zh) | 一种基于卷积神经网络信息融合的水电机组故障诊断方法 | |
CN115859077A (zh) | 一种变工况下多特征融合的电机小样本故障诊断方法 | |
CN112633098A (zh) | 一种旋转机械故障诊断方法、系统及存储介质 | |
CN115060494A (zh) | 一种滚动轴承的故障诊断方法 | |
CN113158984A (zh) | 基于复Morlet小波和轻量级卷积网络的轴承故障诊断方法 | |
CN115452381A (zh) | 基于多重模态分解和长短期记忆网络的强噪声干扰下的滚动轴承故障诊断方法 | |
CN115358259A (zh) | 一种基于自学习的无监督跨工况轴承故障诊断方法 | |
Tao et al. | Fault diagnosis of rolling bearing using deep belief networks | |
CN115587290A (zh) | 基于变分自编码生成对抗网络的航空发动机故障诊断方法 | |
CN115993248A (zh) | 基于脉冲神经网络的轴承故障检测方法、系统、设备和存储介质 | |
CN114048787A (zh) | 一种基于Attention CNN模型的轴承故障实时智能诊断方法与系统 | |
Wang et al. | A deep neural network based on kernel function and auto-encoder for bearing fault diagnosis | |
Wang et al. | Multiscale deep attention Q network: A new deep reinforcement learning method for imbalanced fault diagnosis in gearboxes |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190201 Termination date: 20210627 |
|
CF01 | Termination of patent right due to non-payment of annual fee |