CN111665050A - Rolling bearing fault diagnosis method based on clustering K-SVD algorithm - Google Patents

Rolling bearing fault diagnosis method based on clustering K-SVD algorithm Download PDF

Info

Publication number
CN111665050A
CN111665050A CN202010501122.4A CN202010501122A CN111665050A CN 111665050 A CN111665050 A CN 111665050A CN 202010501122 A CN202010501122 A CN 202010501122A CN 111665050 A CN111665050 A CN 111665050A
Authority
CN
China
Prior art keywords
rolling bearing
algorithm
dictionary
clustering
eigenmode
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
Application number
CN202010501122.4A
Other languages
Chinese (zh)
Other versions
CN111665050B (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.)
Yanshan University
Original Assignee
Yanshan 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 Yanshan University filed Critical Yanshan University
Priority to CN202010501122.4A priority Critical patent/CN111665050B/en
Publication of CN111665050A publication Critical patent/CN111665050A/en
Application granted granted Critical
Publication of CN111665050B publication Critical patent/CN111665050B/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
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/23Clustering techniques
    • G06F18/232Non-hierarchical techniques
    • G06F18/2321Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
    • G06F18/23213Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/28Determining representative reference patterns, e.g. by averaging or distorting; Generating dictionaries

Abstract

The invention discloses a rolling bearing fault diagnosis method based on a clustering K-SVD algorithm, which mainly comprises the following steps: firstly, a time-varying filtering empirical mode decomposition algorithm based on particle swarm optimization is utilized to carry out self-adaptive decomposition on an original signal to obtain each eigenmode component, and correlation kurtosis indexes (K) of each component are calculatedcr) A value; then, the maximum K is selectedcrThe eigen-mode component of the index value is used as an input sample of the clustering K-SVD algorithm to carry out dictionary learning to obtain an ultra-complete dictionary DNew(ii) a Finally, using the overcomplete dictionary DNewAnd sparse feature extraction is carried out on the original signals of the rolling bearing by combining an orthogonal matching tracking algorithm, and envelope spectrum analysis is carried out on sparse representation results to extract the fault frequency features of the rolling bearing. The method effectively solves the problems that the learning precision of the classical K-SVD algorithm on the fault impact characteristics of the rolling bearing is low and the like, and realizes weak fault diagnosis of the rolling bearingTherefore, the method has important significance.

Description

Rolling bearing fault diagnosis method based on clustering K-SVD algorithm
Technical Field
The invention relates to the technical field of equipment maintenance, in particular to a rolling bearing fault diagnosis method based on a clustering K-SVD algorithm.
Background
The rolling bearing is used as an important supporting part and widely applied to rotary machinery, but various damages inevitably happen to the rolling bearing due to severe operating environment and complex and changeable operating conditions, so that the safe and reliable operation of the whole machine is influenced. However, due to the influences of factors such as the installation position of the sensor, other rotating parts, a noisy working environment, electromagnetic interference and the like, the vibration signal acquired by the vibration sensor not only contains periodic pulse components caused by faults, but also has a large amount of noise and harmonic interference, so that effective characteristic information is submerged, and the difficulty in fault diagnosis of the rolling bearing is increased. Therefore, the research on the rolling bearing fault diagnosis technology is developed, the effective extraction of the rolling bearing fault characteristics in strong background noise is realized, the health state of the rolling bearing is accurately identified in time, and the method has important significance for guaranteeing the safe and reliable operation of equipment.
K-singular value decomposition (K-SVD for short) is a completely new non-stationary signal processing method proposed by Elad equal to 2006. The method is a classic dictionary learning algorithm, only one atom is updated in each iteration process, so that the method is more convergent, and the method trains and learns the overcomplete dictionary containing the impact components, so that the defect of poor adaptability of the fixed dictionary structure is overcome. And (4) alternately performing sparse representation of the constraint problem and a singular value decomposition updating algorithm to finally obtain the self-adaptive super-complete dictionary. The algorithm has a complete theoretical basis, and is widely applied to fault diagnosis with strong adaptability and high efficiency.
However, the vibration signal of the rolling bearing not only contains periodic impact components, but also contains a large amount of harmonic components and noise interference, so that the learned dictionary inevitably contains atoms similar to the noise or harmonic components, thereby reducing effective extraction of the impact characteristics of the rolling bearing fault and leading to failure in diagnosing the rolling bearing fault.
Disclosure of Invention
The invention aims to provide a rolling bearing fault diagnosis method based on a clustering K-SVD algorithm, which improves the extraction precision of rolling bearing fault characteristics by realizing the self-adaptive effective extraction of the impact characteristics in vibration signals and weakening the influence of noise and harmonic interference on the sparse representation of the rolling bearing fault impact characteristics.
The technical scheme adopted by the invention is as follows:
the invention provides a rolling bearing fault diagnosis method based on a clustering K-SVD algorithm, which comprises the following steps: s1, acquiring a vibration signal of the rolling bearing; s2, carrying out self-adaptive decomposition on the collected rolling bearing vibration signals by adopting a time-varying filtering empirical mode decomposition algorithm based on particle swarm optimization to obtain a plurality of eigenmode components; s3, calculating correlation kurtosis index (K) of each eigenmode componentcr) (ii) a S4, adopting clustering K-SVD algorithm to carry out correlation kurtosis index KcrPerforming dictionary learning on the largest eigenmode component to obtain an ultra-complete dictionary D containing obvious impact characteristic atomsNew(ii) a S5, adopting orthogonal matching pursuit algorithm and overcomplete dictionary DNewSparse representation is carried out on the vibration signals of the rolling bearing; s6, demodulating and analyzing the sparse representation result by using a Teager energy operator to obtain a demodulation envelope spectrum; and S7, identifying the fault frequency characteristics of the rolling bearing by utilizing the demodulation envelope spectrum, and diagnosing the fault of the rolling bearing.
Further, the step S2 specifically includes:
(1) setting the search ranges of a bandwidth threshold and a B spline order n in a time-varying filtering empirical mode decomposition algorithm as [ 0.11 ] and [ 530 ], and enabling the two parameters to jointly form a position coordinate of the particle;
(2) carrying out self-adaptive decomposition on the vibration signal by using different particle position information and adopting a time-varying filtering empirical mode decomposition algorithm to obtain a plurality of eigen-mode components, and calculating all eigen-modes corresponding to each particle position informationEnvelope entropy E of mode componentpAnd taking the minimum value of the envelope entropy as the fitness value of the particle;
envelope entropy E of signal x (i) (1, 2, …, N)PThe calculation formula is as follows:
Figure BDA0002524846880000021
where p (i) is the envelope of the analytic signal obtained by Hilbert transform of signal x (i), and N is the signal length.
(3) Updating the individual local extremum and the overall local extremum by utilizing the fitness value of each particle, and further updating the position information and the speed information of the particle swarm;
(4) when the iteration times reach the maximum set value or the error is smaller than the set threshold value, outputting the optimal fitness value and the corresponding particle position information;
(5) substituting the obtained optimal particle position information, namely the bandwidth threshold and the B spline order n, into a time-varying filtering empirical mode decomposition algorithm to obtain a final decomposition result u ═ u { for the vibration signal of the rolling bearing1,u2,…uKIn which ukDenotes the kth eigenmode component, and K denotes the number of eigenmode components.
Further, the step S3 specifically includes: for the K eigenmode components obtained in step S2, each eigenmode component u is calculatedkRelated kurtosis index KcrThe calculation formula is as follows:
Figure BDA0002524846880000031
where x represents the original signal and x represents the original signal,
Figure BDA0002524846880000032
and
Figure BDA0002524846880000033
respectively, the original signal x and the k-th eigenmode component ukN is the signal length.
Further, the step S4 specifically includes:
(1) by calculating the correlation kurtosis indexes K of K eigenmode componentscrSelecting the eigenmode component u with the maximum correlation kurtosis indexmAs input to the clustering K-SVD algorithm;
(2) setting K-SVD algorithm parameters, and adopting K-SVD algorithm to carry out eigenmode component umPerforming dictionary learning to obtain an ultra-complete dictionary D;
(3) and analyzing the overcomplete dictionary D by adopting a K-means clustering algorithm.
Further, the analyzing the overcomplete dictionary D by using the K-means clustering algorithm specifically includes:
(1) calculating each atom D in the overcomplete dictionary Dj(J ═ 1,2, …, J denotes the number of atoms in dictionary D) is used as a kurtosis index
Figure BDA0002524846880000034
Index of sum margin
Figure BDA0002524846880000035
As characterizing atom djFeature vector of
Figure BDA0002524846880000036
(2) Constructing a sample set Λ ═ v using the feature vectors of all atoms12,…,νJSetting the clustering number to be 2, and randomly selecting 2 samples from Λ as an initial mean vector (mu) as an input of a K-means clustering algorithm1,μ2};
(3) Calculating v of each input samplejAnd each mean vector mul(l is 1,2) distance
Figure BDA0002524846880000041
Determining v according to the nearest mean vectorjSo that all samples are divided into two classes C1And C2;
(4) updating the corresponding class with all samples in each classAnother mean vector; if the mean vector does not change before and after updating, stopping clustering and outputting a clustering result C ═ C1,C2}; otherwise, returning to execute the step (3);
(5) selecting one kind of atoms with larger characteristic vectors from the clustering result to reconstruct the dictionary, and obtaining an overcomplete dictionary D containing the atoms with obvious impact characteristicsNew
Further, the step S5 specifically includes:
using overcomplete dictionaries DNewSparse representation is carried out on the vibration signal x of the rolling bearing:
Figure BDA0002524846880000042
wherein d iskRepresenting an overcomplete dictionary DNewα denotes the sparse representation coefficient and r denotes the residual. the problem of solving the sparse representation coefficient α is summarized as being based on l0Solving the problem of the coefficients of the norm:
(l0):min||α||0subject to x=Dα+r
solving the l by using an orthogonal matching pursuit algorithm with the sparsity of 10And solving the problem by the norm coefficient to obtain a sparse representation result of the vibration signal of the rolling bearing.
Compared with the prior art, the invention has the following beneficial effects:
the diagnosis method provided by the invention adopts the time-varying filtering empirical mode decomposition algorithm based on particle swarm optimization to carry out self-adaptive decomposition on the vibration of the rolling bearing, so that the self-adaptive selection of algorithm parameters is realized, and the problem of mode aliasing is solved; meanwhile, the relevant kurtosis indexes are used for selecting the mode components containing rich fault information as the input of a K-SVD algorithm, dictionary learning is carried out, and the learning precision of atoms with impact morphological characteristics is improved; meanwhile, the dictionary is further analyzed by adopting a K-means clustering algorithm, atoms irrelevant to the impact form can be eliminated, the influence of noise and harmonic components on the fault feature extraction precision of the rolling bearing is avoided, and the fault diagnosis accuracy of the rolling bearing can be effectively improved.
Drawings
FIG. 1 is a schematic flow chart of a rolling bearing fault diagnosis method based on a clustering K-SVD algorithm, which is provided by the invention;
FIG. 2 is a schematic time domain waveform of a vibration signal of a rolling bearing;
FIG. 3 is a schematic diagram of a demodulated envelope spectrum of a vibration signal of a rolling bearing;
FIG. 4 is a diagram illustrating empirical mode decomposition results of time-varying filtering based on particle swarm optimization;
FIG. 5 is a schematic diagram of atomic clustering;
FIG. 6 is DNewSchematic diagram of atomic waveforms in (1);
FIG. 7 is a schematic diagram of a time domain waveform obtained by the present method;
fig. 8 is a schematic diagram of a demodulation envelope spectrum obtained by the method.
Detailed Description
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
The invention provides a rolling bearing fault diagnosis method based on a clustering K-SVD algorithm, as shown in figure 1, the diagnosis method specifically comprises the following steps:
s1, acquiring a vibration signal of the rolling bearing;
the method specifically comprises the following steps: the rolling bearing of the locomotive with the outer ring peeling fault is arranged on a vibration test bench to be tested to acquire a vibration signal, in the embodiment, the frequency conversion is set to be fr6.33Hz, with the sampling frequency set to fs12.8kHz, data length set to 8000; the time domain waveform of the vibration signal of the rolling bearing and the demodulation envelope spectrum thereof are shown in fig. 2 and 3.
S2, carrying out self-adaptive decomposition on the collected rolling bearing vibration signals by using a time-varying filtering empirical mode decomposition algorithm optimized by particle swarm to obtain a plurality of eigenmode components;
the method specifically comprises the following steps: firstly, setting initial parameters of a particle swarm algorithm, in this embodiment, the population number is set to 50, the maximum evolution algebra is set to 2000, the solving precision is set to 1e-5, and the search ranges of a bandwidth threshold and a B-spline order n in the time-varying filtering empirical mode decomposition algorithm are respectively [ 0.11]And [ 530 ]]And the two parameters jointly form the position coordinates of the particles; secondly, carrying out self-adaptive decomposition on the vibration signal by using different particle position information and adopting a time-varying filtering empirical mode decomposition algorithm to obtain a plurality of eigenmode components, and calculating the envelope entropy E of all eigenmode components corresponding to each particle position informationpWherein, in the step (A),
Figure BDA0002524846880000061
p (i) is the envelope of the analytic signal obtained by Hilbert transform of a signal x (i), N is the signal length, and the minimum value of the envelope entropy is taken as the fitness value of the particle; thirdly, updating the individual local extremum and the overall local extremum by utilizing the fitness value of each particle, and further updating the position information and the speed information of the particle swarm; then, when the iteration times reach the maximum set value or the error is smaller than the set threshold value, outputting the optimal fitness value and the corresponding particle position information; and finally, substituting the obtained optimal particle position information, namely the bandwidth threshold value of 0.43 and the B spline order n of 30 into a time-varying filtering empirical mode decomposition algorithm to obtain a final decomposition result u of the vibration signal of the rolling bearing, wherein the final decomposition result u is { u ═ u { (u) }1,u2,…,u11As shown in fig. 4.
S3, calculating correlation kurtosis index (K) of each eigenmode componentcr);
The method specifically comprises the following steps: using formulas
Figure BDA0002524846880000062
Where x represents the original signal and x represents the original signal,
Figure BDA0002524846880000063
and
Figure BDA0002524846880000064
respectively, the original signal x and the k-th eigenmode component ukThe correlation kurtosis index of each eigenmode component is calculated, and the result is shown in fig. 4.
S4, adopting clustering K-SVD algorithm to carry out correlation kurtosis index KcrPerforming dictionary learning on the largest eigenmode component to obtain an ultra-complete dictionary D containing obvious impact characteristic atomsNew
The method specifically comprises the following steps: firstly, through comparative analysis of correlation kurtosis indexes of 11 eigenmode components, a 1 st eigenmode component u is selected1As input samples of the K-SVD algorithm, wherein the parameters of the K-SVD algorithm are set as: the sample/atom length is 280 data points, and the atom number is 650; modeling the component u using the K-SVD algorithm1Performing dictionary learning to obtain an ultra-complete dictionary D; then, analyzing the overcomplete dictionary D by using a K-means clustering algorithm, calculating the kurtosis index and the margin index of each atom in the dictionary D to form a characteristic vector representing the atom, taking the characteristic vectors of all the atoms as the input of the K-means clustering algorithm, automatically dividing the atoms in the dictionary D into two types, and selecting one type of atoms with larger characteristic vectors to reconstruct the dictionary to obtain the overcomplete dictionary D containing the obviously impacted characteristic atoms as shown in figure 5New. As shown in FIG. 6, a secondary dictionary D is providedNewThe time domain waveform of 4 randomly selected atoms in the cluster is shown in fig. 6, and the clustered atom form has obvious pulse characteristics.
S5, adopting orthogonal matching pursuit algorithm and overcomplete dictionary DNewSparse representation is carried out on the vibration signals of the rolling bearing;
the method specifically comprises the following steps: using overcomplete dictionaries DNewCarrying out sparse representation on the bearing vibration by adopting an orthogonal matching pursuit algorithm with the sparsity of 1, wherein the obtained result is shown in figure 7; in comparison with FIG. 2, in the time domain waveform shown in FIG. 7The noise is obviously reduced, and the periodic impact characteristic is obvious.
S6, demodulating and analyzing the sparse representation result by using a Teager energy operator to obtain a demodulation envelope spectrum, as shown in FIG. 8.
And S7, identifying the fault frequency characteristics of the rolling bearing by utilizing the demodulation envelope spectrum, and diagnosing the fault of the rolling bearing. In the demodulation envelope spectrum, except the fault characteristic frequency f of the outer ring of the rolling bearingoBesides, a multiple harmonic component 2 × f of the fault characteristic frequency can be foundo、3×fo、4×fo、5×foThe local damage of the outer ring of the rolling bearing is predicted; and the identification precision of the fault characteristic frequency is greatly improved.
The rolling bearing fault diagnosis method based on the clustering K-SVD algorithm successfully extracts the fault characteristic frequency of the bearing, compared with the classic K-SVD method, the method utilizes the time-varying filtering empirical mode decomposition algorithm of particle swarm optimization to carry out self-adaptive decomposition on the signal, and utilizes the relevant kurtosis indexes to select the mode components containing abundant fault information, thereby weakening the influence of noise and harmonic interference on the dictionary learning precision; meanwhile, atoms irrelevant to the impact form are further removed by utilizing a K-means clustering algorithm, so that the extraction precision of the fault impact characteristics of the rolling bearing is improved, and the method has important significance for realizing weak fault diagnosis of the rolling bearing.
The above-mentioned embodiments are merely illustrative of the preferred embodiments of the present invention, and do not limit the scope of the present invention, and various modifications and improvements of the technical solution of the present invention by those skilled in the art should fall within the protection scope defined by the claims of the present invention without departing from the spirit of the present invention.

Claims (6)

1. A rolling bearing fault diagnosis method based on a clustering K-SVD algorithm is characterized in that: the method comprises the following steps:
s1, acquiring a vibration signal of the rolling bearing;
s2, carrying out self-adaptive decomposition on the collected rolling bearing vibration signals by adopting a time-varying filtering empirical mode decomposition algorithm based on particle swarm optimization to obtain a plurality of eigenmode components;
s3, calculating correlation kurtosis index (K) of each eigenmode componentcr);
S4, adopting clustering K-SVD algorithm to carry out correlation kurtosis index KcrPerforming dictionary learning on the largest eigenmode component to obtain an ultra-complete dictionary D containing obvious impact characteristic atomsNew
S5, adopting orthogonal matching pursuit algorithm and overcomplete dictionary DNewSparse representation is carried out on the vibration signals of the rolling bearing;
s6, demodulating and analyzing the sparse representation result by using a Teager energy operator to obtain a demodulation envelope spectrum;
and S7, identifying the fault frequency characteristics of the rolling bearing by utilizing the demodulation envelope spectrum, and diagnosing the fault of the rolling bearing.
2. The rolling bearing failure diagnosis method according to claim 1, characterized in that: the step S2 specifically includes:
(1) setting the search ranges of a bandwidth threshold and a B spline order n in a time-varying filtering empirical mode decomposition algorithm as [ 0.11 ] and [ 530 ], and enabling the two parameters to jointly form a position coordinate of the particle;
(2) carrying out self-adaptive decomposition on the vibration signal by using different particle position information and adopting a time-varying filtering empirical mode decomposition algorithm to obtain a plurality of eigenmode components, and calculating the envelope entropy E of all eigenmode components corresponding to each particle position informationpAnd taking the minimum value of the envelope entropy as the fitness value of the particle;
envelope entropy E of signal x (i) (1, 2, …, N)PThe calculation formula is as follows:
Figure FDA0002524846870000021
where p (i) is the envelope of the analytic signal obtained by Hilbert transform of signal x (i), and N is the signal length.
(3) Updating the individual local extremum and the overall local extremum by utilizing the fitness value of each particle, and further updating the position information and the speed information of the particle swarm;
(4) when the iteration times reach the maximum set value or the error is smaller than the set threshold value, outputting the optimal fitness value and the corresponding particle position information;
(5) substituting the obtained optimal particle position information, namely the bandwidth threshold and the B spline order n, into a time-varying filtering empirical mode decomposition algorithm to obtain a final decomposition result u ═ u { for the vibration signal of the rolling bearing1,u2,…uKIn which ukDenotes the kth eigenmode component, and K denotes the number of eigenmode components.
3. The rolling bearing failure diagnosis method according to claim 2, characterized in that: the step S3 specifically includes: for the K eigenmode components obtained in step S2, each eigenmode component u is calculatedkRelated kurtosis index KcrThe calculation formula is as follows:
Figure FDA0002524846870000022
where x represents the original signal and x represents the original signal,
Figure FDA0002524846870000023
and
Figure FDA0002524846870000024
respectively, the original signal x and the k-th eigenmode component ukN is the signal length.
4. The rolling bearing failure diagnosis method according to claim 3, characterized in that: the step S4 specifically includes:
(1) by calculating the correlation kurtosis indexes K of K eigenmode componentscrSelecting the eigenmode component u with the maximum correlation kurtosis indexmAs input to the clustering K-SVD algorithm;
(2) setting K-SVD algorithm parameters, and adopting K-SVD algorithm to carry out eigenmode component umPerforming dictionary learning to obtain an ultra-complete dictionary D;
(3) and analyzing the overcomplete dictionary D by adopting a K-means clustering algorithm.
5. The rolling bearing failure diagnosis method according to claim 4, characterized in that: the analyzing the overcomplete dictionary D by adopting the K-means clustering algorithm specifically comprises the following steps:
(1) calculating each atom D in the overcomplete dictionary Dj(J ═ 1,2, …, J denotes the number of atoms in dictionary D) is used as a kurtosis index
Figure FDA0002524846870000031
Index of sum margin
Figure FDA0002524846870000032
As characterizing atom djFeature vector of
Figure FDA0002524846870000033
(2) Constructing a sample set Λ ═ v using feature vectors for all atoms1,v2,…,νJSetting the clustering number to be 2, and randomly selecting 2 samples from Λ as an initial mean vector (mu) as an input of a K-means clustering algorithm1,μ2};
(3) Calculating v of each input samplejAnd each mean vector mul(l is 1,2) distance
Figure FDA0002524846870000034
Determining v according to the nearest mean vectorjSo that all samples are divided into two classes C1And C2;
(4) updating the mean vector of the corresponding category with all samples in each category; if the mean vector does not change before and after updating, stopping clustering and outputting a clustering result C ═ C1,C2}; on the contrary, the method can be used for carrying out the following steps,returning to execute the step (3);
(5) selecting one kind of atoms with larger characteristic vectors from the clustering result to reconstruct the dictionary, and obtaining an overcomplete dictionary D containing the atoms with obvious impact characteristicsNew
6. The rolling bearing failure diagnosis method according to claim 5, characterized in that: the step S5 specifically includes:
using overcomplete dictionaries DNewSparse representation is carried out on the vibration signal x of the rolling bearing:
Figure FDA0002524846870000035
wherein d iskRepresenting an overcomplete dictionary DNewα denotes the sparse representation coefficient and r denotes the residual. the problem of solving the sparse representation coefficient α is summarized as being based on l0Solving the problem of the coefficients of the norm:
(l0):min||α||0subject to x=Dα+r
solving the l by using an orthogonal matching pursuit algorithm with the sparsity of 10And solving the problem by the norm coefficient to obtain a sparse representation result of the vibration signal of the rolling bearing.
CN202010501122.4A 2020-06-04 2020-06-04 Rolling bearing fault diagnosis method based on clustering K-SVD algorithm Active CN111665050B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010501122.4A CN111665050B (en) 2020-06-04 2020-06-04 Rolling bearing fault diagnosis method based on clustering K-SVD algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010501122.4A CN111665050B (en) 2020-06-04 2020-06-04 Rolling bearing fault diagnosis method based on clustering K-SVD algorithm

Publications (2)

Publication Number Publication Date
CN111665050A true CN111665050A (en) 2020-09-15
CN111665050B CN111665050B (en) 2021-07-27

Family

ID=72386040

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010501122.4A Active CN111665050B (en) 2020-06-04 2020-06-04 Rolling bearing fault diagnosis method based on clustering K-SVD algorithm

Country Status (1)

Country Link
CN (1) CN111665050B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113295410A (en) * 2021-05-14 2021-08-24 上海交通大学 Bearing fault diagnosis method under variable rotating speed working condition
CN113418705A (en) * 2021-07-23 2021-09-21 燕山大学 Early fault diagnosis method for rolling bearing
CN113743240A (en) * 2021-08-12 2021-12-03 武汉钢铁有限公司 Fault detection method, device and computer readable storage medium

Citations (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6178261B1 (en) * 1997-08-05 2001-01-23 The Regents Of The University Of Michigan Method and system for extracting features in a pattern recognition system
WO2006106508A2 (en) * 2005-04-04 2006-10-12 Technion Research & Development Foundation Ltd. System and method for designing of dictionaries for sparse representation
CN102269644A (en) * 2010-06-07 2011-12-07 北京化工大学 Diagnosis method for impact type failure between rolling bearing and gear based on optimal self-adaptive wavelet filter
CN102721545A (en) * 2012-05-25 2012-10-10 北京交通大学 Rolling bearing failure diagnostic method based on multi-characteristic parameter
CN102735447A (en) * 2012-06-29 2012-10-17 西安交通大学 Method for quantitatively identifying performance degradation degree of rolling bearing
CN102778355A (en) * 2012-08-07 2012-11-14 北京交通大学 Rolling bearing state identification method based on empirical mode decomposition (EMD) and principal component analysis (PCA)
CN102797671A (en) * 2011-05-25 2012-11-28 中国石油大学(北京) Fault detection method and device of reciprocating compressor
CN102854015A (en) * 2012-10-15 2013-01-02 哈尔滨理工大学 Diagnosis method for fault position and performance degradation degree of rolling bearing
US8863044B1 (en) * 2013-09-06 2014-10-14 International Business Machines Corporation Layout assessment method and system
CN106441897A (en) * 2016-10-27 2017-02-22 石家庄铁道大学 Morphological component bearing failure diagnosis method based on dictionary study
CN108171263A (en) * 2017-12-26 2018-06-15 合肥工业大学 Based on the Fault Diagnosis of Roller Bearings for improving variation mode decomposition and extreme learning machine
CN108388908A (en) * 2018-01-18 2018-08-10 华南理工大学 Rolling bearing impulse fault diagnostic method based on K-SVD and sliding window noise reduction
CN108507789A (en) * 2018-05-31 2018-09-07 燕山大学 The sparse diagnostic method of rolling bearing fault based on the weak orthogonal matching pursuit of mean random
CN108896306A (en) * 2018-03-26 2018-11-27 四川大学 Method for Bearing Fault Diagnosis based on adaptive atom dictionary OMP
CN109635334A (en) * 2018-11-12 2019-04-16 武汉科技大学 Fault Diagnosis of Roller Bearings, system and medium based on particle group optimizing
CN109707615A (en) * 2019-02-26 2019-05-03 东北石油大学 Reciprocating compressor method for diagnosing faults based on fine multi-fractal
CN109813547A (en) * 2019-01-08 2019-05-28 华南理工大学 Rotating machinery local type method for diagnosing faults based on sparse decomposition optimization algorithm
CN110765965A (en) * 2019-10-30 2020-02-07 兰州理工大学 Quick dictionary learning algorithm for sparse representation of mechanical vibration signals
CN110987435A (en) * 2019-12-18 2020-04-10 北京化工大学 Bearing fault diagnosis method based on improved K-SVD dictionary learning and compressed sensing

Patent Citations (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6178261B1 (en) * 1997-08-05 2001-01-23 The Regents Of The University Of Michigan Method and system for extracting features in a pattern recognition system
WO2006106508A2 (en) * 2005-04-04 2006-10-12 Technion Research & Development Foundation Ltd. System and method for designing of dictionaries for sparse representation
CN102269644A (en) * 2010-06-07 2011-12-07 北京化工大学 Diagnosis method for impact type failure between rolling bearing and gear based on optimal self-adaptive wavelet filter
CN102797671A (en) * 2011-05-25 2012-11-28 中国石油大学(北京) Fault detection method and device of reciprocating compressor
CN102721545A (en) * 2012-05-25 2012-10-10 北京交通大学 Rolling bearing failure diagnostic method based on multi-characteristic parameter
CN102735447A (en) * 2012-06-29 2012-10-17 西安交通大学 Method for quantitatively identifying performance degradation degree of rolling bearing
CN102778355A (en) * 2012-08-07 2012-11-14 北京交通大学 Rolling bearing state identification method based on empirical mode decomposition (EMD) and principal component analysis (PCA)
CN102854015A (en) * 2012-10-15 2013-01-02 哈尔滨理工大学 Diagnosis method for fault position and performance degradation degree of rolling bearing
US8863044B1 (en) * 2013-09-06 2014-10-14 International Business Machines Corporation Layout assessment method and system
CN106441897A (en) * 2016-10-27 2017-02-22 石家庄铁道大学 Morphological component bearing failure diagnosis method based on dictionary study
CN108171263A (en) * 2017-12-26 2018-06-15 合肥工业大学 Based on the Fault Diagnosis of Roller Bearings for improving variation mode decomposition and extreme learning machine
CN108388908A (en) * 2018-01-18 2018-08-10 华南理工大学 Rolling bearing impulse fault diagnostic method based on K-SVD and sliding window noise reduction
CN108896306A (en) * 2018-03-26 2018-11-27 四川大学 Method for Bearing Fault Diagnosis based on adaptive atom dictionary OMP
CN108507789A (en) * 2018-05-31 2018-09-07 燕山大学 The sparse diagnostic method of rolling bearing fault based on the weak orthogonal matching pursuit of mean random
CN109635334A (en) * 2018-11-12 2019-04-16 武汉科技大学 Fault Diagnosis of Roller Bearings, system and medium based on particle group optimizing
CN109813547A (en) * 2019-01-08 2019-05-28 华南理工大学 Rotating machinery local type method for diagnosing faults based on sparse decomposition optimization algorithm
CN109707615A (en) * 2019-02-26 2019-05-03 东北石油大学 Reciprocating compressor method for diagnosing faults based on fine multi-fractal
CN110765965A (en) * 2019-10-30 2020-02-07 兰州理工大学 Quick dictionary learning algorithm for sparse representation of mechanical vibration signals
CN110987435A (en) * 2019-12-18 2020-04-10 北京化工大学 Bearing fault diagnosis method based on improved K-SVD dictionary learning and compressed sensing

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
FATAOHOU 等: ""Weak fault feature extraction of rolling bearings based on globally optimized sparse coding and approximate SVD"", 《MECHANICAL SYSTEMS AND SIGNAL PROCESSING》 *
S. WU, X. CHEN, Z. ZHAO AND R. LIU: ""Data-driven discriminative K-SVD for bearing fault diagnosis"", 《2017 PROGNOSTICS AND SYSTEM HEALTH MANAGEMENT CONFERENCE》 *
刘敏: ""基于K-SVD的字典学习算法及应用研究"", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *
李继猛 等: ""采用粒子群算法的冲击信号自适应单稳态随机共振检测方法"", 《机械工程学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113295410A (en) * 2021-05-14 2021-08-24 上海交通大学 Bearing fault diagnosis method under variable rotating speed working condition
CN113418705A (en) * 2021-07-23 2021-09-21 燕山大学 Early fault diagnosis method for rolling bearing
CN113743240A (en) * 2021-08-12 2021-12-03 武汉钢铁有限公司 Fault detection method, device and computer readable storage medium

Also Published As

Publication number Publication date
CN111665050B (en) 2021-07-27

Similar Documents

Publication Publication Date Title
CN111665050B (en) Rolling bearing fault diagnosis method based on clustering K-SVD algorithm
Wang et al. A novel weighted sparse representation classification strategy based on dictionary learning for rotating machinery
Gao et al. Rolling bearing compound fault diagnosis based on parameter optimization MCKD and convolutional neural network
CN110543860B (en) Mechanical fault diagnosis method and system based on TJM (machine learning model) transfer learning
CN110991295B (en) Self-adaptive fault diagnosis method based on one-dimensional convolutional neural network
CN108106830A (en) A kind of Variable Speed Rotating Machinery method for diagnosing faults based on time-frequency spectrum segmentation
CN112257530B (en) Rolling bearing fault diagnosis method based on blind signal separation and support vector machine
CN111562108A (en) Rolling bearing intelligent fault diagnosis method based on CNN and FCMC
CN112084974A (en) Multi-label rolling bearing fault diagnosis method based on meta-learning
CN113176092B (en) Motor bearing fault diagnosis method based on data fusion and improved experience wavelet transform
CN111191740B (en) Fault diagnosis method for rolling bearing
CN112906644A (en) Mechanical fault intelligent diagnosis method based on deep migration learning
CN112036301A (en) Driving motor fault diagnosis model construction method based on intra-class feature transfer learning and multi-source information fusion
CN109374293B (en) Gear fault diagnosis method
Jin et al. A multi-scale convolutional neural network for bearing compound fault diagnosis under various noise conditions
CN111189638B (en) HMM and QPSO optimization algorithm-based bearing fault degree identification method
CN111753891A (en) Rolling bearing fault diagnosis method based on unsupervised feature learning
CN112881942A (en) Abnormal current diagnosis method and system based on wavelet decomposition and empirical mode decomposition
CN109117896B (en) Rolling bearing fault feature extraction method based on KSVD dictionary learning
CN114112400A (en) Mechanical bearing fault diagnosis method based on multi-angle information fusion
CN110782041A (en) Structural modal parameter identification method based on machine learning
CN110231165B (en) Mechanical equipment fault diagnosis method based on expectation difference constraint confidence network
CN110044619B (en) Multi-fault feature identification method based on sparse multi-cycle group lasso
CN111025100A (en) Transformer ultrahigh frequency partial discharge signal mode identification method and device
CN115905832A (en) Transformer fault diagnosis method based on vibration signal clustering identification

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