Abstract:
A method for digital processing of vibration signals from rolling bearings in rotating machines is presented, allowing subsequent fault detection with high reliability. A linear, adaptive filter is applied to the acquired vibration signal and iteratively tuned to increase the statistical asymmetry of its output. In this process, the filter removes phase- and amplitude distortion from underlying fault impulses. Furthermore, suppression of sinusoidal disturbances is simultaneously achieved with high robustness to measurement noise. The result is a processed signal from which rolling bearing defects are more easily detected.

Description:
BACKGROUND OF THE INVENTION 
       [0001]    Rolling bearings enable low-friction movement under heavy loads in rotating machinery with applications ranging from industrial processes to wind turbines and vehicles. It is an accepted fact that even if all necessary operating conditions such as sufficient lubrication, proper alignment and a good surrounding environment are met, fatigue and mechanical wear of the contact surfaces bearings is unavoidable. Proper functionality of bearings is often critical and a breakdown can cause large economical loss or even pose a threat to personal safety. In heavy industrial processes such as steelmaking, the degradation can be quite rapid due to inherently harsh environments combined with severe loads. 
         [0002]    To avoid unpredicted failures and enable planned maintenance, scheduled or continuous condition monitoring of a bearing&#39;s state is standard practice where feasible. Vibration analysis is a common condition-monitoring technique in which the vibrations from a spinning bearing are measured with a suitable transducer and analysed on some digital processing platform, looking for signs of incipient or imminent damage. Contact surface defects create impulsive vibrations, occurring periodically as rolling elements strike developing pits or cracks in the bearing races. If sufficiently strong, these periodic impacts show up as peaks in the vibration signal&#39;s frequency spectrum which may then be matched to the bearing&#39;s characteristic frequencies, corresponding to each of the different expected faults (e.g. inner/outer race damage). The characteristic frequencies are calculated in advance from the bearing&#39;s dimensions and the axis rotation speed, see N. Tandon and A. Choudhury, “A review of vibration and acoustic measurement methods for the detection of defects in rolling element bearings,”  Tribology international,  vol. 32, pp. 469-480, August 1999. Vibration analysis can be done either by using a portable, hand-held instrument (scheduled monitoring) or with a permanently mounted system performing continuous trend monitoring. 
         [0003]    The main challenge in performing successful vibration analysis is to detect and enhance the weak, impulsive vibrations associated with bearing defects in measured signals plagued by measurement noise and disturbing vibrations. Furthermore, the unknown and non-trivial transfer function between the fault-impulse sources (i.e. rolling elements striking surface defects) and the transducer results in phase and amplitude distortion that tends to obscure the vibration patterns of interest by reducing the signal&#39;s ‘peakiness’, see S. G. Braun, “The signature analysis of sonic bearing vibrations,”  IEEE Transactions on Sonics and Ultrasonics,  vol. SU-27, pp. 317-328, November 1980. A suitable signal pre-processing method is generally required to improve the signal quality before the actual bearing condition assessment is made. Examples of such processing can be found in U.S. Pat. No. 6,868,348 B1 and U.S. Pat. No. 6,648,700 B1. 
         [0004]    The shock pulse method (SPM) is a simple, widespread technique utilizing a vibration transducer with a tuned resonance frequency to detect the presence of fault-induced impulses in a bandpass frequency range around 32 kHz. This center frequency is selected to, in most cases, produce a high signal-to-noise ratio by making the transducer sensitive to fault impulses while filtering out lower frequencies, where most disturbing vibrations reside. A related, more advanced approach is the popular high-frequency resonance technique (HFRT), also called envelope analysis, in which a bandpass-filtered vibration signal is subject to envelope detection to reveal impulsive vibrations. Despite its popularity, HFRT is known to have issues with reliability, especially if the bearing damage is well progressed, (see N. Tandon and A. Choudhury, “A review of vibration and acoustic measurement methods for the detection of defects in rolling element bearings,”  Tribology international,  vol. 32, pp. 469-480, August 1999) and requires careful selection of the bandpass filter&#39;s center frequency depending on the measurement setting, see P. D. McFadden and J. D. Smith, “Vibration monitoring of rolling element bearings by the high-frequency resonance technique—a review,”  Tribology international,  vol. 17, pp. 3-40, February 1984. Also, the non-linear operations involved in envelope detection requires a high signal-to-noise ratio (SNR) for reliable results. 
         [0005]    A different approach that avoids the use of specialized transducers or non-linear signal processing is to employ the principle of deconvolution. A linear, digital filter is used to compensate for the distorting transfer function between the fault-vibration source and the transducer so that any underlying impulses appear less distorted. Simple deconvolution examples for other applications can be found in U.S. Pat. No. 5,744,722 and U.S. Pat. No. 5,130,951. In the deconvolution technique, the measured vibration signal is modeled as a stochastic process and the transfer function as an unknown, linear system. The parameters of an adaptive deconvolution filter are then adjusted in an iterative algorithm to optimise the resulting output signal quality according to some suitable statistical measure. Traditionally, kurtosis is employed as the statistical measure to maximize for deconvolution of bearing vibration signals, see N. Sawalhi, R. B. Randall, arid H. Endo, “The enhancement of fault detection and diagnosis in rolling element bearings using minimum entropy deconvolution combined with spectral kurtosis,”  Mechanical systems and signal processing,  vol. 21, pp. 2616-2633. August 2007. Kurtosis is defined as the normalized fourth moment of a random variable and reflects a signal&#39;s peakiness. The rationale is that the filter setting that maximizes the kurtosis of the output (i.e. maximizes its peakiness) also gives optimum transfer function distortion removal under the given conditions to enhance any underlying impulsive vibrations as much as possible. 
         [0006]    While simple and intuitively attractive, deconvolution by maximizing kurtosis proves less useful in practice as the adaptive deconvolution filter is likely to amplify disturbing sinusoidal vibrations, which will therefore mask any weak impulsive vibrations and cause bearing defects to go un-detected. The reason for this is that the sinusoidal disturbances typically found in vibrations from rotating machines are not clearly distinguished from impulses in terms of their kurtosis value. While sinusoids may be suppressed with an adaptive, noise-cancelling filter prior to deconvolution, this increases riot only the implementation cost but also the risk of accidentally cancelling defect signatures, see D. Ho and R. B. Randall, “Optimisation of bearing diagnostic techniques using simulated and actual bearing fault signals,”  Mechanical signals and signal processing,  vol. 14, pp. 763-788, September 2000. 
         [0007]    More sophisticated processing techniques, such as wavelet-based methods and neural network algorithms have also been considered. While these may give decent reliability, they are far more computationally demanding and are thus only suitable when power and processing capacity is abundant. This excludes compact, low-cost or energy-efficient solutions such as handheld devices and embedded processing systems. 
       OBJECTIVE OF THE INVENTION 
       [0008]    Against the provided background, there exists a need for a reliable signal preprocessing method capable of giving fault-impulse enhancement combined with suppression of sinusoidal disturbances and robustness against noise for condition monitoring of rolling bearings. Furthermore, the growing industrial demand for compact, low-power monitoring systems places emphasis on computational simplicity without sacrificing performance. 
       SUMMARY OF THE INVENTION 
       [0009]    The present invention improves traditional deconvolution of bearing vibration signals by noting that fault impulses (signals of interest) are different from noise and disturbances (undesired signals) in terms of statistical asymmetry, an indicator of which is first defined here for a random variable. 
       Definition 1 (Statistical Asymmetry of a Random Variable) 
       [0010]    Let X be a random variable with zero mean and probability density function (pdf) ƒ X (         ). We define an indicator Φ of the statistical asymmetry of X to be 
         [0000]    
       
         
           
             
               
                 
                   
                     φ 
                      
                     
                       ( 
                       X 
                       ) 
                     
                   
                    
                   
                     = 
                     Δ 
                   
                    
                   
                     
                       
                         ∫ 
                         
                           - 
                           ∞ 
                         
                         ∞ 
                       
                        
                       
                         
                           g 
                            
                           
                             ( 
                             x 
                             ) 
                           
                         
                          
                         
                           
                             f 
                             X 
                           
                            
                           
                             ( 
                             x 
                             ) 
                           
                         
                          
                         
                             
                         
                          
                         
                            
                           x 
                         
                       
                     
                     
                       h 
                        
                       
                         ( 
                         X 
                         ) 
                       
                     
                   
                 
               
               
                 
                   ( 
                   1 
                   ) 
                 
               
             
           
         
       
     
         [0000]    where g(         ) is an arbitrary non-linear and odd function of          , and h(X) is chosen, dependent on g(x), to enforce scale-invariance so that Φ(αX) =Φ(X) for any α&gt;0. 
         [0011]    Φ(X) =0 indicates that X is statistically symmetric. By modeling a vibration signal as a discrete-time (sampled), stationary, stochastic process, statistical asymmetry of such a signal can now be defined. 
       Definition 2 (Statistical Asymmetry of a Stochastic Process) 
       [0012]    Let ν(n) be a discrete-time, stationary, stochastic process over time index n for which each sample is a unit-variance random variable X with pdf ƒ X (         ). Then, the statistical asymmetry of the process ν(n) is Φ(X). 
         [0013]    Under Definitions 1 and 2 above, fault vibrations from rolling bearings are recognized as being statistically asymmetric (Φ≠0) while linearly distorted variants, as well as disturbances such as sinusoids and measurement noise, are essentially symmetric (Φ=0), which therefore makes them distinguishable from the vibration patterns of interest. The linear deconvolution filter should be used to increase the statistical asymmetry of the measured vibration signal for the following reasons:
       1. Increasing the statistical asymmetry reduces transfer function distortion. The phase- and amplitude distortion caused by the transfer function between the fault-vibration source and the transducer tends to obscure the underlying impulses by making them appear as damped oscillations or, in other words, as more statistically symmetric vibrations. As the deconvolution filter should ‘undo’ this distortion, it should hence be tuned to increase the asymmetry of its output.   2. An algorithm increasing the statistical asymmetry is robust against noise and disturbances. Unlike kurtosis, statistical asymmetry distinguishes between the signal of interest and disturbances. Therefore, increasing asymmetry will not accidentally enhance disturbances.   3. Increasing statistical asymmetry can be achieved using low-cost algorithms. Simple embodiments are possible that require relatively small amounts of computations and hence have low power requirements.       
 
         [0017]    Statistical asymmetry as an indicator of bearing defects has been proposed in U.S. Pat. No. 5,852,793 and U.S. Pat. No. 7,505,852 B2. Deconvolution methods related to enhancing statistical asymmetry have also been evaluated for other applications, see for example J.-Y. Lee and A. K. Nandi, “Blind deconvolution of impacting signals using higher-order statistics,”  Mechanical systems and signal processing,  vol. 12, pp. 357-371, March 1998. However, the approach of increasing statistical asymmetry with an adaptive, linear filter for condition monitoring of rolling bearings has not previously been explored. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0018]      FIG. 1  is a block diagram showing the principle of condition monitoring using vibration analysis. 
           [0019]      FIG. 2  is a block diagram of the process of enhancing statistical asymmetry with a linear, adaptive filter. 
           [0020]      FIG. 3  illustrates the vibration signal (a) and corresponding histogram (b) from a functional rolling bearing. 
           [0021]      FIG. 4  illustrates the impulsive vibration source signal (a) and corresponding histogram (b) from a defective rolling bearing. 
           [0022]      FIG. 5  illustrates a measured noise- and disturbance-free vibration signal (a) and corresponding histogram (b) from a defective rolling bearing. 
           [0023]      FIG. 6  illustrates generic measurement noise (a) and corresponding histogram (b). 
           [0024]      FIG. 7  illustrates a sinusoidal disturbance (a) and corresponding histogram (b). 
           [0025]      FIG. 8  illustrates a typical measured vibration signal (a) and corresponding histogram (b) from a defective rolling bearing. 
           [0026]      FIG. 9  illustrates the output signal (a) and corresponding histogram (b) from a linear adaptive filter aimed to increase the statistical asymmetry of a measured vibration signal from a defective rolling bearing. 
       
    
    
     DESCRIPTION OF THE PREFERRED EMBODIMENT 
       [0027]    An embodiment of the invention is now described in which it is assumed that a vibration signal from a rotating machine with one or several bearings is converted to electrical form using a suitable transducer, such as an piezoelectric accelerometer, and given as input to a vibration analysis process. Other possible means of supplying such an input include reading a previously measured signal from a storage media. 
         [0028]    An overview of the steps involved in general vibration analysis is found in  FIG. 1 . illustrating the use of signal pre-processing followed by bearing condition assessment. The input vibration signal  1  from the transducer is first converted from analog to digital form. The sampling rate used in the analog-to-digital (A-to-D) conversion  2  must be greater than twice the highest vibration frequency of interest for the current application. The digital signal is then input to a signal preprocessing step  3  to suppress transfer-function distortion and disturbances that act to obscure any underlying impulsive vibrations significant for fault detection. This type of processing can for example be performed on specialized hardware such a a digital signal processor (DSP), a general-purpose microcontroller or any generic computer platform. Once processed, the resulting signal is input to a fault detection algorithm  4 , producing an estimate  5  of the bearing condition. The bearing condition may be presented in the form of a value on a continuous scale (e.g. expected remaining life time) or as one of a finite set of indicators (e.g. OK/replace). 
         [0029]    The present invention concerns the signal pre-processing step  3 , the principle of which is shown in  FIG. 2 . A digital, adaptive, linear filter  6  takes as input the signal          (n)  7  from the A-to-D converter  2 . The filter parameters are iteratively adjusted towards increasing the magnitude of the statistical asymmetry Φ of the output signal y(n)  8  using a suitable learning rule  9 . 
         [0030]    As one example of an embodiment of the invention, the adaptive filter is assumed to be of finite impulse response (FIR) type, represented at time n by the N×1 parameter vector 
         [0000]      f n           [ƒ 0   (n) ƒ 1   (n)  . . . ƒ N−1   (n)]   T ,   (2)
 
         [0000]    where ƒ i   (n)  is the i th  parameter at time n. The filter output signal sample y(n) at time n is obtained as 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       y 
                        
                       
                         ( 
                         n 
                         ) 
                       
                     
                     = 
                     
                       
                         
                           f 
                           n 
                           T 
                         
                          
                         
                           x 
                           n 
                         
                       
                       = 
                       
                         
                           ∑ 
                           
                             i 
                             = 
                             0 
                           
                           
                             N 
                             - 
                             1 
                           
                         
                          
                         
                           
                             f 
                             i 
                           
                            
                           
                             x 
                              
                             
                               ( 
                               
                                 n 
                                 - 
                                 i 
                               
                               ) 
                             
                           
                         
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   3 
                   ) 
                 
               
             
           
         
       
     
         [0000]    where T denotes vector transpose and x n, is an N× 1 vector storing the N latest samples of the input signal          (n), 
         [0000]      x(n)         [         (n)         (n− 1 ) . . .          (n−N+ 1 )] T .  (4)
 
         [0031]    With each new input sample, the corresponding output sample y(n) is calculated as (3) and the filter parameter vector f n  is updated in an iterative manner to increase the magnitude of the output signal&#39;s statistical asymmetry. A possible realization of this is to form a parameter-vector update term aimed at increasing the skewness of the output signal y(n) through a small parameter adjustment. The skewness Ψ of a zero-mean random variable X is defined as its normalized third moment, 
         [0000]    
       
         
           
             
               
                 
                   
                     Ψ 
                      
                     
                       ( 
                       X 
                       ) 
                     
                   
                    
                   
                     = 
                     Δ 
                   
                    
                   
                     
                       
                         E 
                          
                         
                           { 
                           
                             X 
                             3 
                           
                           } 
                         
                       
                       
                         
                           ( 
                           
                             E 
                              
                             
                               { 
                               
                                 X 
                                 2 
                               
                               } 
                             
                           
                           ) 
                         
                         
                           3 
                           / 
                           2 
                         
                       
                     
                     . 
                   
                 
               
               
                 
                   ( 
                   5 
                   ) 
                 
               
             
           
         
       
     
         [0032]    For a signal ν(n) modelled as a discrete-time, stationary, stochastic process where each sample is a random variable X, the skewness of ν(n) is simply defined as the skewness of X, i.e. Ψ(ν(n))=Ψ(X). 
         [0033]    The filter update term is formed from an estimate of the N×1 gradient vector ∇ Ψ (n) of the output signal skewness with respect to the current parameter vector, 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       ∇ 
                       Ψ 
                     
                      
                     
                       ( 
                       n 
                       ) 
                     
                   
                    
                   
                     = 
                     Δ 
                   
                    
                   
                     
                       
                         ∂ 
                         
                           Ψ 
                            
                           
                             ( 
                             
                               y 
                                
                               
                                 ( 
                                 n 
                                 ) 
                               
                             
                             ) 
                           
                         
                       
                       
                         ∂ 
                         
                           f 
                           n 
                         
                       
                     
                     = 
                     
                       
                         
                           [ 
                           
                             
                               
                                 ∂ 
                                 
                                   Ψ 
                                    
                                   
                                     ( 
                                     
                                       y 
                                        
                                       
                                         ( 
                                         n 
                                         ) 
                                       
                                     
                                     ) 
                                   
                                 
                               
                               
                                 ∂ 
                                 
                                   f 
                                   0 
                                   
                                     ( 
                                     n 
                                     ) 
                                   
                                 
                               
                             
                              
                             
                               
                                 ∂ 
                                 
                                   Ψ 
                                    
                                   
                                     ( 
                                     
                                       y 
                                        
                                       
                                         ( 
                                         n 
                                         ) 
                                       
                                     
                                     ) 
                                   
                                 
                               
                               
                                 ∂ 
                                 
                                   f 
                                   1 
                                   
                                     ( 
                                     n 
                                     ) 
                                   
                                 
                               
                             
                              
                             
                                 
                             
                              
                             … 
                              
                             
                                 
                             
                              
                             
                               
                                 ∂ 
                                 
                                   Ψ 
                                    
                                   
                                     ( 
                                     
                                       y 
                                        
                                       
                                         ( 
                                         n 
                                         ) 
                                       
                                     
                                     ) 
                                   
                                 
                               
                               
                                 ∂ 
                                 
                                   f 
                                   
                                     N 
                                     - 
                                     1 
                                   
                                   
                                     ( 
                                     n 
                                     ) 
                                   
                                 
                               
                             
                           
                           ] 
                         
                         T 
                       
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   6 
                   ) 
                 
               
             
           
         
       
     
         [0034]    Hence, the vector ∇ Ψ (n) gives the direction of adjustment for f n  giving the most rapid increase of the output signal skewness. An estimate {circumflex over (∇)} Ψ (n) of (6) is calculated from a finite set of signal samples and used to form the iterative filter parameter update at time n as 
         [0000]        f   n+1   =f   n +μ{circumflex over (∇)} Ψ ( n ),   (7)
 
         [0000]    where μ is a small, constant stepsize parameter. The iterative update process (7) proceeds either indefinitely or until some criterion is attained, such as a certain level of skewness in y(n). 
         [0035]    The rationale behind the use of statistical asymmetry in the present invention is illustrated in  FIGS. 3-8 , showing the characteristics of different bearing-related signals.  FIGS. 3   a  and  4   a  show simulated vibration patterns from a functional and defective hearing respectively. The corresponding histograms, or estimated pdf&#39;s, are displayed in  FIGS. 3   b  and  4   b  respectively, which also show the estimated skewness (Ψ) of the signals. Note how the impulses present in the signal from the defective bearing ( FIG. 4   a ) increases the asymmetry in the pdf resulting in a higher skewness, while the skewness of the functional bearing vibration signal ( FIG. 3   a ) is close to zero. In  FIG. 5   a,  the defective-bearing impulses from  FIG. 4   a  have been convolved with an exponentially decaying sinusoid, a typical model of a transfer function between the impulsive source and the transducer. Hence,  FIG. 5   a  represents an ideal (noise- and disturbance free) vibration signal from a faulty bearing. Note from  FIG. 5   b  how the phase- and amplitude distortion caused by the transfer function decreases the skewness to result in an almost symmetric signal. 
         [0036]    The characteristics of two common kinds of disturbances are found in  FIGS. 6 and 7 .  FIG. 6  shows white Gaussian measurement noise which is symmetric.  FIG. 7  illustrates a common disturbing vibration in the form of a sinusoid, also being essentially symmetric. 
         [0037]    The signal in  FIG. 8   a  is the sum of those from  FIGS. 5   a,    6   a  and  7   a  with different individual scaling, modelling a typical measured vibration signal from a faulty bearing consisting of weak, distorted impulsive vibrations plus noise and a strong sinusoidal disturbance. As seen from the resulting histogram in  FIG. 8   b,  this vibration pattern is practically symmetric, reflected by its low skewness. Note how the originally clear impulsive fault signatures from  FIG. 4   a  are now hard to distinguish in  FIG. 8   a.    
         [0038]    To demonstrate the functionality of the invention, the signal from  FIG. 8   a  was applied to an adaptive filter as in  FIG. 2 , subject to skewness enhancement according to (7) using an FIR filter with 64 filter parameters. The resulting output signal after 500 iterations is shown in  FIG. 9   a.  The impulsive vibration pattern becomes clearly visible and the increased asymmetry is apparent from  FIG. 9   b.  The restored impulsive pattern confirms that increasing asymmetry with an adaptive filter can successfully perform suppression of both transfer function distortion and sinusoidal disturbances in the presence of measurement noise. 
         [0039]    Alternative embodiments of the invention may be obtained by letting the filter parameter vector f n  be iteratively adjusted to increase, instead of skewness, the magnitude of any odd, normalized moment of order ≧3 of the filter output signal. The m th -order odd moment Ψ m  of a zero-mean random variable X is defined as 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         Ψ 
                         m 
                       
                        
                       
                         ( 
                         X 
                         ) 
                       
                     
                      
                     
                       = 
                       Δ 
                     
                      
                     
                       
                         E 
                          
                         
                           { 
                           
                             X 
                             m 
                           
                           } 
                         
                       
                       
                         
                           ( 
                           
                             E 
                              
                             
                               { 
                               
                                 X 
                                 2 
                               
                               } 
                             
                           
                           ) 
                         
                         
                           m 
                           / 
                           2 
                         
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   8 
                   ) 
                 
               
             
           
         
       
     
         [0000]    where m is any positive, odd integer greater ≧3. Note that the choice m=3 gives skewness as defined by (5). The filter update term is thus formed from an estimate of the gradient vector 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         ∇ 
                         
                            
                           
                             Ψ 
                             m 
                           
                            
                         
                       
                        
                       
                         ( 
                         n 
                         ) 
                       
                     
                      
                     
                       = 
                       Δ 
                     
                      
                     
                       
                         ∂ 
                         
                            
                           
                             
                               Ψ 
                               m 
                             
                              
                             
                               ( 
                               
                                 y 
                                  
                                 
                                   ( 
                                   n 
                                   ) 
                                 
                               
                               ) 
                             
                           
                            
                         
                       
                       
                         ∂ 
                         
                           f 
                           n 
                         
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   9 
                   ) 
                 
               
             
           
         
       
     
         [0000]    and the filter parameter update becomes 
         [0000]        f   n+1   =f   n μ{circumflex over (∇)} |Ψm| ( n ),   (10)
 
         [0000]    where μ it is a small, constant stepsize parameter and {circumflex over (∇)} |Ψm| (n) is an estimate of ∇ |Ψm| (n). 
         [0040]    While the present invention has been described by preferred embodiments above, it is not the intention of the applicants to restrict or limit the scope of the claims accordingly. It is recognized that numerous other variants of embodiments are possible to those skilled in the art.