1. Introduction
Models of machine failures enable the generation of data used for the training of various data analysis systems [
1], which play important roles in maintenance scheduling. This training generally covers data pre-processing, feature selection, and data analysis. Generating synthetic data is far cheaper, faster, and more convenient than collecting real data from individual machine failures. Such models of machine damage are typically constructed on the basis of the physical nature of the damage [
2,
3] (physical models), or on the basis of data generated by the object, so-called data-driven models [
4,
5]. Theoretically, the physical model should generate data corresponding to any possible machine damage; however, they are burdened with a relatively high level of approximations and assumptions. Additionally, any real material variations (such as manufacturing variations) or operational variations must be modeled separately, resulting in a large set of model variables. Data-driven models, on the other hand, are inherently associated with true operational conditions of machinery, but their range of fault representation (e.g., failure mode, size of the damage, fault severity) is limited to the current true machine state and is, therefore, burdened with undetermined uncertainty beyond it. An interesting approach to rolling-element bearing (REB) signal modelling is given in [
6], where the authors claim that lack of data on failure models within data-driven models could be solved by data fusion techniques. The authors propose to combine a complete dataset representing a machine healthy state with synthetically generated failure mode data. As shown, this approach could provide the data necessary to train models in a wide variety of failure classifications. Nevertheless, the failure mode data used in multi-body techniques and data-driven techniques (namely, neural networks (NN) and support vector machines (SVM)) are in favor within feature extractions, classification, diagnostics, and analysis of remaining useful life (RUL).
The basic impulsive character of faulty rolling-element bearing is addressed in, among others, [
7], as well as by McFadden and Smith [
8,
9]. The following research examines the influence of manufacturing errors of fault signature [
10], and the effect of instantaneous load and speed with respect to rolling-element bearing defect size [
11]. A comparative study of a simulated physical model with a finite-element model of faulty REB is presented in [
12]. Recently, Wrzochal and Adamczyk [
13] provide a detailed mathematical study on different REB modeling techniques, namely, the basic model with four degrees of freedom, a model with changing defect topography, a model with outer ring deformation, a model with waviness of rings, and, finally, a dynamic model with a variable viscosity damping coefficient. In their study, they tackle the problem of modeling simultaneous deformations within the ball–raceway contact zone, lubrication, friction process, and the geometrical structure of working surfaces. As concluded in [
14,
15,
16,
17], some of these problems, considered analytically in [
13] might be solved by the numerical techniques shown in [
18], especially regarding the modeling of dynamic conditions originating from different manufacturing, operational conditions, and assembly imperfections. Two year later, Borghesani et al. supplemented analytical REB models with analytical consideration in an assessment of true- vs. pseudo-cyclostationarity in REB signals [
19].
The current paper proposes a hybrid model, which combines phenomenological features of an REB-induced vibration signal with its behavioral features. This model partially takes advantage of the analytical consideration presented in [
20], where the REB signal model for constant speed presented in Equation (1) is extended to a more general, transient state, presented in Equation (2):
where
represents average machine speed between (
k − 1)-th vs.
k-th REB impact,
is the average angular period of pulses,
n(
t) is the environmental noise represented by a repetition of the impulse response of the structure
s(
t), triggered, and A-modulated by subsequent REB impacts. The standard variable T represents expected period time, while
represents jitter (in both time and angle domains). As shown in the paper, the accepted scenario eventually follows the concept of the reversed squared envelope spectrum presented in [
20], as it enables the generation of distinctive envelope components for REB vibration signals under simultaneous variable speed regime and jitter. While [
20] discusses the general difference in true bearing characteristic frequencies and analytical values with respect to transient states, the current paper discusses the influence of the speed profile generation method on the smearing of the envelope spectrum in the order domain.
The paper is organized as follows:
Section 2 shows a step-by-step construction of a generalized synthetic vibration signal from a real single pulse to the entire model.
Section 3 illustrates the frequency and envelope analysis of the generated signal.
2. Signal Generation Methodology
2.1. Overview
The presented method enables the simulation of rotary machine faults under non-stationary conditions. In this paper, a local fault of a rolling-element bearing is considered. Formally, the presented model could be classified as a hybrid parametric model, because it combines the parametric phenomenological model with the parametric behavioral model.
The
first part of the model—the phenomenological part—in contrast to the physical nature of failure modes used in some models [
2,
3], does not operate on primary object properties (mass, damping, and stiffness), but takes advantage of failure signal patterns, which are associated with particular failure modes. These patterns are based on the kinetostatic analysis of a particular drive train, as well as on well-accepted knowledge on vibration analysis and failure modes [
21].
The
second part of the model—the behavioral part—is, in a sense, a data-driven method, but modified. In this part, the real signal recorded from a machine of interest is used to determine characteristics of random structural vibrations for failure simulation. In this way, experimental observation of the behavior of the machine is used to approximate the real signal. As this observation is made without precise apparatus, which give detailed information about root processes (just a simple measurement of system response is required), this part of the model is classified as a behavioral part. The scope of observations considered in the presented model includes instantaneous speed characteristics, REB slip, amplitude–frequency characteristics of structural vibrations, phase–frequency characteristics of structural vibrations, and generalized amplitude and frequency modulations of 2nd order cyclostationary components [
22], extended to a cyclo-non-stationary regime [
23].
Figure 1 illustrates the general methodology of the generation of a hybrid vibration signal simulating a REB fault. The proposed method has three general stages. In the first stage, a real vibration pulse is collected from the mechanical object of interest. In stage 2, this signal is used to generate a set of pulses, which occur at pre-calculated timestamps. Finally, the amplitude of individual pulses is scaled according to the real amplitude–frequency characteristics of the object as a function of instantaneous machine speed. In this way, the synthetic signal enables the generation of a faulty REB signal, without an actual REB fault.
Section 2.2,
Section 2.3 and
Section 2.4 describe the consecutive steps of the methodology in detail.
2.2. Base Pulse
First, the housing of the REB of interest needs to be provided, as illustrated in
Figure 2. With the use of some kind of metal object, vibrations of the housing were induced. This response is similar to a system response generated by a faulty bearing, but this action does not require the bearing to be faulty. During the experiment, a standard bearing housing on the AVM test bench was selected.
This operation was repeated multiple times. As illustrated in
Figure 3, in practice, some pulses generate saturation of acquisition path, while others use a relatively small channel range. From a set of pulses, a single pulse is selected, as illustrated in
Figure 4. This repetition enables the selection of such a pulse, which covers a majority of the physical channel range, yet it does not saturate it. In the selected example, the pulse covers over 90% of the physical input range of a 24-bit ADC (analog-to-digital) converter.
Figure 4 shows a time plot of the pulse. As visually assessed, the full damping of the pulse lasts ca. 0.5 s. The next question is how much could it be shortened to still represent structural vibrations of the investigated object and be sufficient for signal generation. This shortening determines boundary condition of the proposed algorithm.
Figure 5 illustrates four versions of this pulse: 0.5 s (full length), 0.2 s, 0.05 s, and 0.01 s. Individual pulses are presented in the same figure. Listed duration times correspond to REB characteristic frequencies from 1 Hz up to 100 Hz.
The applicability of individual pulses is verified by the shape of their amplitude–frequency characteristics of the approximated frequency response function. From each shortened version of the base pulse, a corresponding full-resolution, one-sided, scaled spectrum is plotted.
Figure 6 shows spectra for five different duration times of pulses: 2 s, 1 s, 0.2 s, 0.05 s, and 0.01 s, which correspond to full resolutions equal to 2 Hz, 5 Hz, 20 Hz, and 100 Hz, respectively. As shown, for all investigated duration times, the amplitude–frequency characteristics show similar dominant resonant frequencies of the structure, which are ca. 1.7–2 kHz and 5–8 kHz. This means that the presented technique enables the modelling of faulty REB signals, for which the characteristic frequencies (ball-passing frequency of the outer race (BPFO), ball-passing frequency of the inner race (BPFI), fundamental train frequency (FTF), and ball spin frequency (BSF)) are all in the range 1–100 Hz.
After initial selection and verification, a single base pulse was selected. This pulse is shown in full-scale in
Figure 7, and in a zoomed version in
Figure 8. This step is referred to as Stage 1 in
Figure 1. The base pulse was used in two ways. During the construction of the hybrid vibration signal, individual pulses were generated from the SHORTENED base pulse, taking into account the duration of each consecutive pulse. Next, the FULL LENGTH pulse was used to generate high-resolution amplitude–frequency characteristics, followed by calculations of the instantaneous amplitude of pulses on the basis of machine instantaneous speed.
The base pulse illustrated in
Figure 7 and
Figure 8 reaches roughly ±50 [g] in amplitude, with the no. of unique points exceeding 5.6 k values (for a 24-bit ADC data acquisition unit (in MATLAB, use length(unique(x)), where x represents the pulse raw waveform)).
In the studied example, the signal is modelled with sampling frequency 25 kHz and duration 0.5 s.
2.3. Generalized Frequency Modulation
The method presented in this paper is oriented toward modeling signals of non-stationary machine operational conditions, which requires instantaneous speed data as input, such as the one presented in
Figure 9. Machine speed profiles can be synthetically generated [
24], or they can be extracted from machine process data with additional up-sampling.
The nominal speed selected was 1000 RPM (rotations per minute), which is ca. 16.667 Hz. In the studied example, the arbitrary speed fluctuation was selected to be +/−5.5 Hz (referred to as “
FM_factor”). The characteristic bearing BPFO order was arbitrarily selected to be equal to 4.88 (data coming from the author’s experience in industry), which results in the nominal frequency of pulses equal to 81.33 Hz. Taking into account the machine instantaneous speed data, and the nominal frequency of bearing pulses, the instantaneous frequency of pulses is calculated, and illustrated in
Figure 10.
For each pulse in the signal, the time of the pulse was calculated taking the base pulse and trimming it to the duration corresponding to instantaneous frequency of this pulse. Additionally, bearing jitter was selected to be 2% [
25], and added to the signal. The subsequent pulses were concentrated. This step is referred as Stage 2 in
Figure 1.
Figure 11 illustrates the constant amplitude pulses generated from the base pulse, the duration of which are a function of the machine instantaneous speed.
In next step, the amplitude of the pulses were modified.
2.4. Generalized Amplitude Modulation
The amplitude of individual pulses were scaled according to the amplitude–frequency characteristics of the frequency response function of the object and the machine operational speed, as illustrated in
Figure 12. For the considered machine, the speed oscillates in the range of 11–22 Hz, so only this range of amplitude–frequency characteristics of simplified frequency response function (FRF) is of interest. As the resolution of this characteristic is 2 Hz (due to 0.5 s pulse length), there are only 7 points in the 11–22 Hz range. On the other hand, the resolution of the instantaneous speed is much higher (25,000 pints). For this reason, a large portion of FRF profile data needed to be interpolated [
24].
Generally, the presented scaling methodology represents the change in power of the vibration signal generated by a rotary machinery as its speed changes. In particular, the scaling signal amplitude, as a function of instantaneous frequency data, and the FRF profile represent variable input forces of the system, which result in the variable system response for the time-invariant FRF.
2.5. Resultant Hybrid Signal
Figure 13 represents resultant hybrid vibration signal and corresponding machine instantaneous speed. As observed, frequency of pulses changes proportionally to the speed, which results in frequency signal modulation (FM). Simultaneously, the amplitude of pulses changes as a function of speed and FRF, resulting in profiled amplitude modulation (AM). For this reason, such data is classified as a generalized AM–FM signal.
The data presented in
Figure 13 corresponds to Stage 3 in
Figure 1. The following part of the paper shows spectral analysis of the generated hybrid signal.
4. Conclusions
The presented technique shows the concept of a hybrid model of vibration signal generated by a faulty rolling-element bearing. The model uses real pulse response data for two purposes; namely, as a base pulse, and as an amplitude ratio profile calculated from the amplitude–frequency characteristics of the frequency response function of the structure. The instantaneous machine speed data are used as a separate input entry, and this can be real historical data, or a smooth simulated profile. The presented model is verified for signals of length 1–10 s, sampling frequency 25 kHz, and a rolling-element bearing nominal characteristic frequency from 10 Hz to 100 Hz. Due to its hybrid nature, and the simplicity of generation, the presented model serves as an attractive alternative to other physical or data-driven models from the literature discussed in the paper. The proposed method enables the generation of vibration signals corresponding to defects of the real selected bearing mounted on the machinery, without constructing its physical model, and without collecting data from individual failure modes. Finally, the paper shows that the methods of analysis of cyclo-nonstationary signal components, such as those generated by a faulty REB, could perform better if these components are treated as individual, consecutive, time-invariant signal fragments under a global time-variant regime.