Analysis of Nasa Bearing Dataset of the University of Cincinnati by Means of Hjorth’S Parameters
Total Page:16
File Type:pdf, Size:1020Kb
View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Archivio istituzionale della ricerca - Università di Modena e Reggio Emilia ANALYSIS OF NASA BEARING DATASET OF THE UNIVERSITY OF CINCINNATI BY MEANS OF HJORTH’S PARAMETERS Cavalaglio Camargo Molano, J., Strozzi, M., Rubini, R., Cocconcelli, M. DISMI, University of Modena and Reggio Emilia Via G. Amendola 2 - Pad. Morselli, 42122, Reggio Emilia, Italy [email protected] ABSTRACT Hjorth’s parameters are statistical time-domain parameters used in signal processing and introduced by Bo Hjorth in 1970. These parameters are Activity, Mobility and Complexity. They are related to the variance of the signal and of its subsequent derivatives. They are commonly used in the analysis of electroencephalography (EEG) signals for feature extraction, but also in the tactile signal analysis in robotic area. In this paper, Hjorth’s parameters are applied to vibration signals for fault detection in ball bearings. In particular, an open-access dataset is used, those data cover a complete lifetime of a set of four bearings. They were run to failure under constant load and running conditions. The robustness of the dataset has been proven by a recent analysis proposed by Gousseau et al., who assessed that only two of the three datasets available are consistent and can be used for bearing diagnosis. Hjorth’s parameters clearly identify the time when a damage occurs and its evolution in the time-domain. An intuitive explanation is given in the text to justify the good results obtained in the detection of fault on the outer race, inner race and rolling elements of the bearing. Hjorth’s parameters can be easily computed with a minimum of computational resources, suggesting their use in real-time application for condition monitoring of ball bearings. KEYWORDS: Hjorth’s parameter, bearing, diagnostics 1. INTRODUCTION With respect to the last decades, today the condition monitoring and diagnostics of mechanical components is a need for every maintenance center in Industry, despite the product developed or services offered by the company. A sudden break in production line has costs of missing product that abundantly exceed the costs of the component itself. As a consequence, the condition monitoring techniques have been developed so far in the literature, for the most common components such as gearboxes [1], ball bearings [2] and rotating shafts [3]. This paper focuses on the condition monitoring of ball bearing, i.e. on the early detection of incipient damage in the components the bearing is made of. In particular, a ball bearing is made of an inner race integral with the rotating shaft, an outer race integral with the fixed frame, rolling elements that link together the two races and a cage that keep the rolling elements equispaced in their seat. The most used diagnostics technique consists in the signal processing of the vibration signal acquired by an accelerometer placed near the bearing [4]. Incipient fault induces cyclic impacts of the rolling elements that could be collected and analyzed. The periodicity of impacts is highlighted in frequency domain, despite the need of a Signal-to-Noise ratio enhancement to reduce mechanical or electric noise [5]. So far [6], thousand of papers on condition monitoring appeared in Literature proposing from simple statistical parameters [7] to complex cyclostationary indexes [8]. These papers have their foundations on the mechanical behaviour of a faulted bearing and its dynamic analysis, or on the characteristics of the expected vibration signal [9]. More recently, machine learning techniques came into play powerfully [10]. Machine learning techniques take into account raw data or a set of features from both healthy and faulted bearings, determining a similarity between a test dataset and the reference ones. As drawback, these techniques require historical data, that are not available so frequently. The model-based and data-driven approaches could be combined to offer a quick diagnostics service and reliable maintenance of the algorithms [11]. This paper focuses on the application of simple statistical values for the diagnostics of ball bearing: the Hjiorth’s parameters. They are statistical time-domain parameters used in signal processing and introduced by Hjorth and Elema-Schönander in 1970 [12]. These parameters are Activity, Mobility and Complexity. They are related to the variance of the signal and of its subsequent derivatives. They are commonly used in the analysis of electroencephalography (EEG) signals for feature extraction [13], but also in the tactile signal analysis in robotic area [14]. In this paper, Hjorth’s parameters are applied to vibration signals for fault detection in ball bearings. In particular, an open-access dataset is used, those data cover a complete lifetime of a set of four bearings [15]. They were run to failure under constant load and running conditions. The results shows that Hjort’s parameter could early detect an incipient fault on different parts of a ball bearing. The small computational efforts needed, suggests their use in real-time application for condition monitoring of ball bearings or as a feature array for machine learning techniques. The paper is structured as follows: Section2 introduces the mathematical definition of the Hjorth’s parameters both in time and frequency domain. Section3 shows the results on open-access experimental dataset and Section 4 draws conclusions. 2. MATHEMATICAL THEORY OF HJORTH’S PARAMETERS In this section, the mathematical theory of Hjorth’s parameters is developed. These parameters are initially defined in the time domain, but could be transferred to the frequency (spectral) domain also. Finally, a dynamic relationship correlating the parameters both in the time and in the frequency domain is obtained. 2.1. Time domain analysis Let us consider a generic vibration signal y = f (t) (e.g. displacement, velocity or acceleration) defined in the time domain t. The variance s 2(y) of a vibration signal y represents a measure of the dispersion or variability of the signal about its mean value, in the form [16]: s 2(y) = y2 − (y)2 (1) where y2 is the mean square value (i.e. the mean power) of the signal [16]: T 1 Z y2 = lim y2(t)dt (2) T!¥ T 0 y is the mean value of the signal [16]: T 1 Z y = lim y(t)dt (3) T!¥ T 0 and T is the signal acquisition period. It should be noted that, for a vibration signal y presenting mean value y equal to zero, the variance s 2(y) is equal to the mean power y2 of the signal, i.e. to the square of the root mean square value RMS(y) of the signal itself, which is defined as [16]: v u T u 1 Z RMS(y) = lim u y2(t)dt (4) T!¥ tT 0 Moreover, the square root of the variance s 2(y) is called standard deviation s(y) of the signal y, in the form [16]: q s(y) = y2 − (y)2 (5) After introducing the previous statistical quantities, now we can start defining Hjorth’s parameters in the time domain. Hjorth’s parameters [12] are statistic time-domain parameters adopted in the vibration signal processing that are related to the definition of variance s 2(y) and are calculated on the basis of the vibration signal amplitude y, its time-derivative (velocity in the following)y ˙ = dy=dt and its two-times derivative (acceleration in the following) y¨ = d2y=dt2 [17]. The first Hjorth’s parameter is called Activity. The Activity of a vibration signal is defined as the variance of the vibration signal amplitude y: Activity(y) = s 2( f (t)) [y] (6) where the dimensions of the Activity parameter are the same of the vibration signal considered. The second Hjorth’s parameter is called Mobility. The Mobility of a vibration signal is given by the square root of the ratio between the Activity of the vibration signal velocityy ˙ and the Activity of the vibration signal amplitude y (i.e. the square root of the ratio between the variances of vibration signal velocityy ˙ and amplitude y): s s Activity(y˙) s 2(d f (t)=dt) Mobility(y) = = [t−1] (7) Activity(y) s 2( f (t)) where the dimensions of the Mobility parameter are expressed as a ratio per time unit. The last Hjorth’s parameter is called Complexity. The Complexity of a vibration signal is denoted by the ratio between the Mobility of the vibration signal velocityy ˙ and the Mobility of the vibration signal amplitude y (i.e., the ratio between the square root of the ratio between the variances of vibration signal accelerationy ¨ and velocity y˙ and the square root of the ratio between the variances of vibration signal velocityy ˙ and amplitude y): r s 2(d2 f (t)=dt2) Mobility(y˙) s 2(d f (t)=dt) Complexity(y) = = r (8) Mobility(y) s 2(d f (t)=dt) s 2( f (t)) It must be underlined that the Complexity is a dimensionless parameter. Moreover, from (6-8), it should be observed that the Activity, Mobility and Complexity parameters are time-related, as the Activity is associated to the variance of the vibration signal amplitude y, the Mobility is related to the variance of the vibration signal velocityy ˙ and the Complexity is connected to the variance of the vibration signal accelerationy ¨ in the time domain t. 2.2. Frequency domain analysis The generic vibration signal y = f (t) previously defined in the time domain t can be transferred to the frequency (spectral) domain w, in the form y = F(w), by means of the Fourier transform [16]: +¥ Z F(w) = f (t)e− jwt dt (9) −¥ where the corresponding inverse Fourier transform is: +¥ 1 Z f (t) = F(w)e− jwt dw (10) 2p −¥ and the complex conjugate of the Fourier transform is: +¥ Z F∗(w) = f ∗(t)e jwt dt (11) −¥ Let us introduce the power spectrum function S(w), which is defined as the product of the Fourier transform F(w) and its complex conjugate function F∗(w): +¥ Z S(w) = F(w) · F∗(w) = f (t) · f ∗(t)dt (12) −¥ where the contribution of the phase jwt in the frequency domain w is not taken into account.