On the use of signal smoothness estimates and signal similarity in bearing fault detection

The concepts of signal smoothness and signal similarity are very common in digital signal processing. In this paper, both concepts have been explored in the ball bearing failure diagnostics, where the bearing fault, bearing inner race and outer race faults were analysed. It was found out that the signal smoothness and the signal similarity estimates might reduce computational tolerance in the vibro-acoustic signal analysis. On this basis, numerous experiments have been carried out to demonstrate usefulness of the above-mentioned estimates (as useful parameters of signal features) in the signal quality analysis.


Introduction
The rolling bearing is one of the most important and widely used components in various mechanical units and mechanisms.Unexpected fault of the item can cause significant material damage, interrupting the units work.Defective bearing often leads out of a stable rotational axis of the rotor, then a chain reaction of faults can occur throughout the unit.To avoid such accidents, it is highly important to identify the bearing defect at the earliest stage.That is why the diagnostic methods of bearing defects have been intensively investigated in recent years.Vibro-acoustic signals are often used in the condition monitoring and preventive maintenance of the rolling bearing [1].
In this study, a new fault bearing feature extraction method based on signal smoothness and signal similarity estimates is proposed.
The organization of the rest of this paper is as follows: Section 2 reviews the notions and mathematical interpretation of signal smoothness and signal similarity.The basic defect detection principles are introduced in Section 3. Section 4 gives some experimental analysis results.Finally, some conclusions are drawn in Section 5.

Signal smoothness and signal similarity: mathematical interpretation
The notion of smoothness of a digital vibro-acoustic signal = ( (0) (1) … ( − 1)) ( = 2 , ∈ ℕ) is introduced in a spectral domain, i.e. the original signal , first and foremost, is subjected to an arbitrary discrete transform (Fourier, Cosine, Walsh-Hadamard, etc.).It is well known that any discrete transform (DT) is uniquely specified by its matrix , namely: here: is the DT spectrum of and is the DT matrix.If the DT basis vectors (rows of the matrix ) are represented in a frequency order, then the spectral coefficients ( ) ( ∈ 0, 1, … , − 1 ) have tendency to decrease in absolute value, as their serial numbers increase.So, the array of coefficients | ( )| = 1, 2, … , − 1 can be approximated by a hyperbolic curve = / , where 0 and 0.
The quantity , characterizing the shape of the hyperbolic curve, i.e. the rate of "decay" of spectral coefficients is assumed to be the smoothness parameter (level, class) of the signal .
In this paper, we explore the discrete cosine transform (DCT), because this transform distinguishes itself with very good signal energy "compaction" (at lower frequencies) property [2].
For more details (computational procedures, two-dimensional case), you can refer to [3].
For assessing the fragmental similarity between and , the mean squared error (metrics) is used, i.e.
The vibro-acoustic signals (fragments) and are said to be similar if = ( , ) ≤ ( being a priori fixed positive number).
It has been proved [4] that small changes in some values of the fragment led to small changes in the smoothness parameter value , i.e.: where depends on .Consequently, the precondition for non-similarity of signal fragments and is defined by the following relationship (the necessary condition for fragmental similarity): In other words, signal fragments and cannot be similar if their smoothness parameter values differ significantly (in the mean squared error sense).

The general bearing fault detection idea
To describe the potential capability and effectiveness of application of both the signal smoothness estimates and the signal similarity to vibro-acoustic signal defect detection process, the following considerations have been taken into account: (1) The original signal (rolling bearing signal under normal conditions) = ( (0) (1) … ( − 1)) ( = 2 , ∈ ℕ) is divided into a finite number of non-overlapping fragments of equal size 2 ( < ).For each fragment , the smoothness parameter value (estimate) is found.Then, all possible fragments (of double or triple size, in comparison with ) of the same signal are processed, i.e. their smoothness estimates are computed.Finally, for each fragment (there are 2 fragments, in total), a plot of points (a "cloud") | − |, ( , ) is determined, and the centre | − | ., .( , ) of the "cloud" is found.
(2) By using all the obtained centres | − | ., .( , ) , the bearing fault detection criterion domain ℬ is formed, namely: (3) For the (rolling bearing) test signal , precisely the same computational steps are performed.The only difference, fragments are taken from the original signal and compared with larger fragments derived from the signal .If a priori prescribed percentage (usually, established experimentally) of centre points − ., . ( , ) fall into the criterion domain (rectangle) ℬ (Eq.( 6)), then the test signal indicates "no fault", and otherwise (too many centres escape the domain ℬ), the rolling bearing fault is detected.

Experimental analysis results
To corroborate the proposed theoretical ideas, a few rolling bearing (vibro-acoustic) signals of size 8192 have been analysed.The test bearing data, i.e. fault diameters (0.1778 mm), frequency (12000 Hz), shaft rotating speed of the motor (1772 rpm) and the signals themselves were obtained from the Bearing Data Center of Case Western Reserve University [5].
The rolling bearing signal under normal conditions is presented in Fig. 1, and the rolling bearing signals under ball fault, inner race fault and outer race fault conditions are shown in Fig. 2. First of all, the defect detection criterion domain ℬ was found out, that it is: ℬ = | − |, ( , ) 0.296 ≤ | − | ≤ 2.009; 0.048 ≤ ( , ) ≤ 0.101 (Fig. 3).Then for each fragment (| | = 16), of the vibro-acoustic signal , all possible pairings with double-sized fragments of the test signal (Fig. 2) have been examined.We here note, that: in the case of under ball fault conditions, 58 % of centre points escaped ℬ; in the case of under inner and outer race fault conditions, all centre points (i.e. 100 %) escaped ℬ.
Often, for practical applications, to minimize computational time expenditure, the necessary signal similarity condition (Eq.( 5)) is employed, i.e. in the general bearing fault detection steps, only pairings ↔ , with | − | ≤ ( being a priori prescribed positive number) are analysed.

Conclusions
This paper presents mathematical interpretation of fragmentary smoothness and similarity of vibro-acoustic signals.The smoothness of the signal is assumed to be a real number characterizing the downtrend of the spectral DT coefficients, as their serial numbers increase.
The main attention is paid to a newly developed vibro-acoustic signal defect detection technique based on the use of signal smoothness and signal similarity estimates.
Preliminary experimental analysis results showed that the proposed technique was sensitive enough to detect abnormal conditions in the rolling bearing signal.

Fig. 3 .
Fig. 3.The formation of the rectangular domain ℬ of signal defect detection criterion In Fig. 4-Fig.6, both the criterion domain (rectangle) ℬ and some typical situations, with earlier-mentioned centre points −