# Investigation on the Effects of Structural Dynamics on Rolling Bearing Fault Diagnosis by Means of Multibody Simulation

Authors: Reza Golafshan; Georg Jacobs; Matthias Wegerhoff; Pascal Drichel; Joerg Berroth
Journal: International Journal of Rotating Machinery
 (2018)

---

## Abstract

The present study aims to combine the fields modal analysis and signal processing and to show the use of Frequency Response Function (FRF), as a vibration transfer path, in enhancing reliability and abilities of the next generation vibration-based rolling bearing condition monitoring (CM) systems in complex mechanical systems. In line with this purpose, the hereby-presented paper employs an appropriate numerical model, that is, Multibody Simulation (MBS) of a vehicle’s drivetrain as a manner for numerical modal and structural analyses. For this, first, the principles of vibration-based bearing fault detection are reviewed and presented. Following that, a summary of MBS modelling and validating strategies are given. Then, the validated MBS model is used as a case study for further investigations. The results can confirm existence of challenges in fault detection of rolling bearings, in particular in complex mechanical systems. In further discussions, the capability of FRFs in fault localization and determination of ideal sensor positions is discussed in some detail. Finally, concluding remarks and suggestions for future works are summarized.


---

## Body

## 1. Introduction

In a vibration-based fault detection and fault diagnosis system, mechanical damage is referred to changes on geometric properties of the analyzed machine. These sudden changes may generate additional dynamical forces acting in the system. Then, based on the system responses, damage/faults can be detected and localized. Since all machines vibrates even in their good operational condition due to interactions between machine components, the captured vibrations signals are required to analyze using signal postprocessing techniques. Many diagnostic procedures assume a linear behavior, and thus each individual spectral component (or a group of them) is supposed to relate to a specific dynamic force in the corresponding machine [1]. Therefore, captured time domain vibration signals need to be transformed to the frequency domain to distinguish the various sources and to detect the probable faults, accordingly [2]. However, due to the mechanical complexities and varying operating conditions, it is not always possible to capture, separate, process, and identify the faulty signals, and this could be the case in many complex systems; hence, an overview of dynamics of the analyzed structure/machine may be needed for the reliable and next generation vibration-based fault detection systems [3, 4]. It is known that the captured vibration signals contain structural properties (i.e., vibration transfer path) of the system being analyzed and they are assumed to remain same under varying operating conditions (e.g., torque loads and rotating speeds). Thus, a prestudy on dynamics of the analyzed structure/machine may be necessary to develop suitable and reliable condition monitoring (CM) and detection techniques to detect and diagnose various failures and achieve cost savings to the industry. The modal behavior of a system can be a useful insight to gain beneficial information on the system response to various possible faults and failures. According to the findings in the literature, modal domain features (i.e., natural frequencies, mode shapes, and damping ratios) are widely used in the field Structural Health Monitoring (SHM), in particular for composite and large civil structures [5]. By tracking and analyzing changes over time the presence of a failure, from crack propagation to unbalance localization, can be approved. However, for mechanical and rotating systems, it is assumed that modal features remain more or less constant in case of presence of a mechanical fault, for example, pitting in bearings. This mentioned assumption is the basis for some of the recent and nonstationary signal processing techniques [6]. By eliminating the estimated structural response characteristics, a less sensitive signal to speed and load changes can be obtained [7]. These methods are based on removing the estimated vibration transfer path effects and modal characteristics of the system from measured signals. Note that the dynamical analysis of a structure/machine can be performed using the well-known Experimental Modal Analysis (EMA) techniques to extract the exact vibration transfer path in its frequency domain. In this regard, dynamic analysis can also be performed using continuous (permanent) measurements. For example, Operational Modal Analysis (OMA), as one of the most common techniques for dynamic analyses of running machines, is used for CM of wind turbines, for example, for tower and blades [8, 9], where some important modal info and parameters can be extracted, tracked, and analyzed for a specific application.In addition to the traditional use of modal testing for design, Noise Vibration Harshness (NVH), and fault localization purposes, the Frequency Response Function (FRF) in modal analysis, as a vibration transfer path function, can be used in other applications such as fault isolation and ideal sensor placement. However, it should be pointed out that the modal analysis on a real structure/machine may not be feasible due to test environment restrictions or due to the lack of reliable mount data. It can therefore be shown that, with a reliable numerical model in hand, the corresponding dynamical analysis can be carried out numerically and reliably. The Multibody Simulation (MBS) approach has approved its usability in identification of system dynamic (modal) characteristics, as a reliable alternative for the real machine [10]. The MBS models can be used for assessments of load interactions between various components in a mechanical system [11]. These models can also represent a testing manner [12], in case there are some restrictions to performing measurements on the real machine. Recently, it is shown that the online or up-to-minute condition monitoring is also a growing application field for the so-called Digital Twin concept using MBS modelling [13, 14]. Only relying on sensor data is sometimes not sufficient. In this situation the simulation (i.e., MBS) models, provided by the Digital Twin, can regularly be modified through measured data and reused after these modifications. By using the simulation models, it is also possible to interpret the measurements in a different way, rather than just detecting deviations from the norm.For improving the performance of the machine monitoring and fault detection procedures, it can be shown that, as a complementary investigation in the scope of ideal sensing positions as well as estimation of the fault location, a model-oriented approach using MBS modelling may be a useful strategy. In addition to this, the model-oriented designing approach for CM and fault diagnosis algorithms may offer a new degree-of-freedom, in particular in the context of Blind Source Separation (BSS) family methods. The core of BSS-based signal processing method may be optimized according to insights from a comprehensive MBS model, where the excitation sources can be specified numerically and reliably. This becomes even more critical for those systems having many subsystems (i.e., components) running in varying operating conditions.Rolling bearings, as one of the most common and important components in many mechanical engineering applications (e.g., drivetrains), are always considered as a part with potentially high-risk against the failures. Bearings have very significant impact on the global vibrations of a rotary mechanical system. Furthermore, in electrified vehicles, a bearing failure may lead to increasing the rotational friction of the rotor, and thus a decrease in the performance of the electric motor. Therefore, vibration-based fault detection and diagnosis of the rolling bearings have become a vigorous area of work and have attracted more and more attention in the literature [15, 16]. These studies in condition monitoring are mainly focused on finding the best signal processing technique to detect an incipient failure. Since rolling element bearings are subjected to moving distributed radial and axial loads, any probable fault can excite the system resonances in a specific frequency range depending on neighboring components. It can be concluded that the captured bearing vibration signals carry the information about the system structure in a comparatively high frequency range, which can be interpreted as the transmission path(s) of the vibration. It is also known that the inherent frequency of a bearing is fixed resonance frequency/frequencies even if the fault repetition frequencies (bearing fundamental frequencies) tend to change with shaft speed. This can be a key factor in localizing the similar faulty bearing(s) rotating at the same speed.Many research studies are dedicated to the field of vibration-based condition monitoring (CM), diagnosis, and fault detection, mostly in oil and gas, wind energy, and railway, including tracks and vehicles, sectors, yet fewer ones are focused on specific requirements and demands for a similar system in the field of automotive engineering. Although numerous general-purpose fault detection and diagnosis algorithms and methods using signal processing techniques are reported in the literature, implementation, applicability, reliability, and ideal (optimal) sensing positions for both offline and online CM systems are still a challenge due to the involved computational complexity and expensive hardware. On-Board Diagnostic (OBD) systems in automotive engineering applications are recently employed and used for various purposes (e.g., Co2 emission and electronic components). In order to achieve a reliable and cost-effective an OBD CM and fault detection system for the drivetrain of a vehicle, the diagnostic algorithm should be simple, robust, and beneficial from ideally located sensors. Since thousands of data points are captured and need to be processed during the monitoring process, characteristics of the captured signals, number and locations of the sensor(s), and the detection algorithm itself have become critical parameters for a precision monitoring.The present study therefore aims to address and investigate some challenges and complexities in vibration-based fault detection and diagnosis for rolling bearings in automotive sector for a drivetrain of an electrified vehicle. In line with this purpose, an experimentally validated MBS model is used for the investigations. Note that a large number of manufacturers, in particular in automotive and wind energy sectors, already have and use numerical models for NVH and testing purposes, so that the modelling step in the development and design procedure of a dedicated and/or implemented condition monitoring systems would not be a time consuming step in designing the next generation monitoring and testing systems. It is worth stating here that this study should be considered as a basis for future works in combination of the fields vibration-based fault detection and numerical modelling.In the present study, the use of structural dynamics in condition monitoring and fault detection of the complex mechanical systems are studied. The MBS model of a drivetrain for an electrified vehicle created and validated in a previous study is used in numerical analyses as a case study. First, a comprehensive theoretical basis of the vibration-based fault detection for rolling bearings is given. Then, a summary of the modelling strategy as well as subsequent revisions applied to the model to create faulty components (i.e., bearings) is described. Three artificially created faulty cases using the MBS model are then used for numerical investigations on bearing fault diagnosis. In further discussions, the role of Frequency Response Function (FRF), as the transfer path function, in faulty bearing(s) localization and also in finding the ideal sensor location is studied. Finally, some concluding remarks are made about the reliability and applicability of the approach used in this study for bearing fault diagnosis and sensor placement.


## 2. Vibration-Based Rolling Bearing Fault Diagnosis: Theoretical Background

In general, a mechanical system’s dynamic characteristics (i.e., vibration response) are specified by the corresponding transfer path and excitation sources, as governed by the mathematics of convolution. According to the previous findings in the literature [2], a typical captured vibration signal by a sensor, x(t), of a machine consists of shaft(s), electric motor, rolling bearing(s), gear stage(s), coupling(s), and so on in its operating condition with an individual focus on rolling bearing vibrations can be considered as(1)xt=ht∗s1t+s2tyt+nt,where yt, ht, and nt correspond to the bearing, the transfer path which modulates the fault(s) response to a higher frequency range (i.e., system resonances), and noise originating from measurement system, numerical errors, and/or incidental and transient excitations, and so on, respectively. While both s1t and s2t are related to uncontaminated periodic responses of all other machine parts (due to electric motor, gears, shaft misalignment and unbalances, etc.), s1t represents the unmodulated part and s2t represents the modulated part of periodic responses with the bearing signal. Also, ∗ stands for the convolution operator. Eq. (1) can then be extended as(2)xt=ht∗s1t+s2t∑n=1N∑k=-∞∞antδt-kTn+nt,where n, ant, and δt, respectively, are the number of fault(s) existing in the bearing (varying from 1 up to N = 4, as there are generally 4 types of faults in a rolling bearing), the amplitude of the vibration response to the nth fault which reflects the strength of the fault, and Dirac delta function, where Tn denotes the time period between impulses (i.e., 1/BPFO or 1/BPFI). Note that the maximum amplitude of ant can be assumed to be constant over time for outer race defects due to its stationary condition [17]. The amplitude spectrum of the vibration signal model in (2) for a system having a bearing with a single fault can therefore be derived using Fourier Transform as(3)Fxt=Xω=Fht∗s1t+s2t∑k=-∞∞atδt-kT+nt.In case the signal energy distribution is desired, the energy spectral density of ω can also be obtained as(4)Xω2=XωXω-,where Xω- denotes the complex conjugate of Xω.It is well-known that the amplitude spectrum of vibration signals may not detect and display the bearing fault frequencies. In such cases, envelope spectrum is the most common method for the fault frequency detection. Nevertheless, in harsh operating conditions, the calculated envelope spectrums are not capable of displaying the fault diagnostic information (i.e., bearing fault frequencies) due to heavy background noise and/or existence other excitation sources. As a result, some preprocessing steps, prior to the envelope analysis, are strongly recommended in the literature.Recall from convolution properties(5)Ff∗gt=FωGωand Nyquist sampling theorem(6)Fft∑k=-∞∞δt-kTr=1Tr∑k=-∞∞Fω-kωr.Using the convolution and Nyquist sampling properties, the right-hand side of (3) can be expanded as follows:(7)Xω=Fht∗s1t+Fht∗s1t∑k=-∞∞atδt-kT+Fht∗nt.Then,(8)Xω=HωS1ω︸non-relateduncontaminatedperiodiccomponentse.g. gear  stages+HωS2ω1T∑k=-∞∞Aω-kω︸Modulated  impulsivecomponents,  whichcontain  bearing  diagnosticinformation,Aω-kω+HωNω︸Noiserepresenting  random  andtransient  components.Eq. (8) reveals that Xω contains various and crucial information regarding current status of the test machine (i.e., bearing system). These three constituents can be categorized as (I) uncontaminated period components, (II) modulated impulsive components, and (III) noise. From rolling bearing fault detection point of view, the diagnostic info is associated with the term HωS2ωyω.A wide range of studies available in the literature are dedicated to signals containing multiple sources [6, 18]. Source separation techniques are proposed in order to separate and eliminate (if needed) the nonrelated periodic components. Perhaps the most common prewhitening technique for eliminating the nonrelated periodic components in rolling bearing vibration signals is autoregressive (AR) filtering, documented and suggested in [2]. Prior to further postprocessing, the other component needs to be eliminated is the noise, HωNω, representing the random and unwanted signals. This can be carried out using denoising algorithms, such as SVD-based denoising method presented in [19]. As a result, after the prewhitening and denoising processes, the denoised version of Xω, called Xdω, can be estimated as(9)Xdω=HωS2ωyω=HωS2ω1T∑k=-∞∞Aω-kω.In order to extract the modulating component,Aω-kω, in (9), demodulation should be performed. For this, an appropriate estimation of Hω is required in order to extract the fault diagnostic info. One method to estimate Hω in the desired frequency range (i.e., system resonance) is the spectral kurtosis (SK) [20]. Using the same property of convolution theorem given earlier, demodulating in the frequency domain can then be formulated as(10)S2ωAω-kω≈lowpass  filterHω^∗S2ωyω,where Hω^ is an estimation of Hω. Again, here S2 represents, for example, shaft-related components originating from unbalances, torsional vibrations, and eccentric faults. In some cases, in particular in complex mechanical systems, estimation of Hω can be a complex and time consuming step.Back to the present theoretical formulation in (9), in practice, however, demodulation on the bearing signal, yt, is carried out by a bandpass filtering (e.g., SK-based optimized filter) with a central-frequency at the system resonance, and then Hilbert Transform in the time domain as(11)At=y2t+iyH2t,where At is also known as the envelope signal and yH is the Hilbert Transform of yt as(12)yHt=1π∫-∞∞yt1t-τdτ.The so-called envelope spectrum is then be obtained by Fourier Transform of the envelope signal, At, as(13)Envyt=FAt.As is being shown, the transfer path function,Hω, plays a major role in diagnostic vibration signals captured from a mechanical system. In the following sections in the present study it is shown that the transfer path can be analyzed in more detail in order to estimate fault location and ideal sensing positions. This function can be represented by Frequency Response Function.


## 3. Summary of Numerical Modelling

In this section, a brief summary of numerical modelling is presented. In addition, some details of the model revision for fault detection purposes are described.In general, the dynamical systems can be presented in terms of their equation of motion (E.o.M.) as(14)Mu¨+Cu˙+Ku=f,where {u}, {f}, [M], [C], and [K] stand for the displacement vector, the force vector (including external and internal torques and moments), inertia (mass), damping, and stiffness matrices of the involved bodies, respectively. In Elastic Multibody Simulation (EMBS), the E.o.M. may be solved in an iterative way, either in time domain or frequency domain. For this, the system matrices and the force (e.g., containing the excitations) vector need to be formed first. However, prior to assemblies of the machine parts by implementation of the masses, stiffness, and the external excitations, EMA is employed in order to validate the modal properties of each elastic parts used in the model. Following this bottom-up process, the entire dynamical system (i.e., drivetrain) consisting of all validated subsystems is validated against the experiments. Some details of the numerical modelling and validation process are readily available in [21]. The previously created EMBS model of a test Battery Electric Vehicle (BEV) for which the numerical investigations are performed is shown in Figure 1.Figure 1
The validated MBS model of an electrified vehicle’s drivetrain.In line with the purpose of rolling bearing fault diagnosis, the numerical model investigated and validated in the previous paper is used in numerical analysis. As discussed in the previous section, in rolling bearings with localized defects, an excitation force may be produced when the rolling elements move over the faulty area [22]. This force may then results in sudden and sharp changes in vibration trend of the system, that is, y(t) in (1), accordingly. The magnitude and the duration of this force are related to the radial load carried by the rolling element, the defect natures (severity and age), velocity, and the material properties used in the bearing [23]. Assuming a dependency between the excitation force due to the localized bearing defect and its generated vibrations, the bearing faults are modelled in terms of their excitation forces by amplifying the magnitudes of the corresponding force element in a certain zone in the force vector in (14). A computer program developed in Fortran Programming Language is used to generate the excitation force, considering the bearing kinematics. Again, the magnitude of the excitation force depends on various parameters, and this is beyond the scopes of the present study. Some various shapes and magnitudes, however, are generated and used in numerical investigations. The developed computer program is capable of locating the given number of defective bearings at the given positions and creates the dynamic excitation forces suitable for calculations (i.e., time integration) in a Multibody Simulation platform, namely, SIMPACK. As the rolling elements pass through the faulty area in the bearing, duration of each impact by interaction between rolling elements and faulty area highly depends on the shaft speed. In general, the impulse-like force, fet, generated due to repetitive interactions of the rolling elements with the defective area, can be expressed as(15)fet=At∑k=-∞∞δt-kT,where At and T are magnitude of the force and time period between impulses (i.e., 1/BPFO), respectively. Note that T may also be considered as a time-varying parameter for nonstationary operating conditions (i.e., nonconstant shaft velocity). It is then assumed that the desired excitation force depends on the impulse-like force obtained in (15), and the radial load, Pθ. The excitation force in the radial direction, fbt, representing the fault, can be therefore obtained as(16)fbt=fetPtcos⁡θ,where Pt and θ are the radial load on the rolling element/raceway and rolling element position, respectively. Figure 2 shows three distinct forces artificially generated and applied to the test bearings in the MBS model. In the present study, it is assumed that fe1, fe2, and fe3 represent three different fault stages on outer race as a tiny crack (e.g., initial stage), medium-size fault (e.g., intermediate-stage), and a large-size fault (e.g., advanced stage), respectively.Figure 2
Illustration of three distinct impulse-like forces for the single localized bearing fault with respect to angular position of the rolling element.Three deep groove ball bearings, two of which are 6205 type and one is 6305 type, are used as the test bearings. The test bearings are labeled and shown in Figure3.Figure 3
Illustration of the faulty test bearings, namely, B1, B2, and B3.Table1 lists the corresponding geometric details of ball bearings used in the tests. The fundamental fault frequency for outer race, BPFO, where outer race is stationary and inner race rotates in shaft can be formulated as(17)BPFO=NB2fr1-Dbcos⁡βDp,where are NB, fr, Db, Dp, and β number of rolling elements, shaft rotating frequency, rolling element diameter, pitch diameter, and rolling bearing contact angle, respectively.Table 1
The geometry details of test bearings.
Bearing Bearingtype Pitchdiameter Balldiameter Numberof balls Contactangle B1 6205 type 38.9 mm 7.9 mm 9 0° B2 6205 type 38.9 mm 7.9 mm 9 0° B3 6305 type 43.4 mm 11.1 mm 7 0°In addition, four sensor locations are specified and shown in Figure4. Note that sensors measure acceleration in normal to surface direction. The related fault frequencies for the test bearings are obtained and presented in Table 2. Three case studies, namely, single faulty bearing, multiple faulty bearings rotating in different velocities (different shafts), and multiple faulty bearings rotating in the same velocity (same shaft), used in numerical investigations, are also listed and detailed in Table 2. The fault stages for all the case studies are assumed as a tiny crack (i.e., fe1) on outer races. The other fault stages (i.e., fe2 and fe3) have been investigated in the section further discussions in some detail.Table 2
Details of the numerical case studies.
Case study Test bearing(s) Fault type Fault freq. Hz Fault freq. ord #1 B1 Outer raceExc.: fe1 32 [Hz] 3.551 #2 B1 Outer raceExc.: fe1 32 [Hz] 3.551 B3 Outer raceExc.: fe1 11 [Hz] 2.61 #3 B1 Outer raceExc.: fe1 32 [Hz] 3.551 B2 Outer raceExc.: fe1 32 [Hz] 3.551Figure 4
Illustration of the sensor positions, namely, S1, S2, S3, and S4.The acceleration signals of the MBS model in normal to surface direction are captured as fault signatures. The data are collected for 2 s at 20 kHz rate in SIMPACK. In all the simulations, both drive and intermediate-stage shafts are given fixed nominal angular velocities in the steady-state area which are shown in Figure5. As seen, both shafts are subjected to some torsional vibrations due to fluctuations in torque loads from electric motor. The simulation results are exported to MATLAB for postprocessing. Note that only the second half of the simulations, that is, from 1 s to 2 s, are used for postprocessing to avoid high amount of nonstationarities during run-up.Figure 5
Angular velocities of (a) drive shaft (b) intermediate-stage shaft.
(a) (b)It is also worth mentioning that in order to avoid from unnecessary complexities, the gears mesh stiffness values are assumed to be constant over time. As a result, the termHωS1ω in (8) has not been taken into account in the present study.


## 4. Results

In this section, the captured vibration signals of simulations are analyzed and investigated in both time and the frequency domains. The aim is to show and discuss the advantages, disadvantages, and limitations of the well-known envelope analysis in bearings fault detection and diagnosis. In what follows, the details of these simulations are presented for three different case studies.


### 4.1. Case Study 1: Single Faulty Bearing

This case study investigates the dynamical behavior of the MBS model having one single faulty bearing on the drive shaft. Figure6 illustrates the artificially generated excitation force applied to B1 in time domain and its representation in frequency domain. As stated earlier, this excitation force is generated using the excitation function, fe1, to analyze the vibration behavior of the model, as well as the real drivetrain, in a condition with a bearing having an initial localized fault such as a tiny crack. The magnitude of force acting on the outer race of B1 with a certain repetition frequency (i.e., BPFO) is assumed to be 120 N in this case study.Figure 6
Illustration of the applied excitation force to B1 in (a) time domain and (b) frequency domain.
(a) (b)Figure7 shows the signals captured from four different locations in the time domain as well as their spectrums up to the Nyquist frequency. As expected, the excitation force excites the system resonances in a specific frequency range. However, by analyzing the force spectrum in Figure 6(b) and the obtained spectrums in Figure 7, it is revealed that small amount of energy in a high frequency range strongly excites the system resonances. This behavior may be explained by transfer path of the vibration from B1 to the sensing locations. In the next section, it is shown that using this valuable insight, an ideal sensor location (between four sensing locations) can be selected. Again, the results presented here suggest that S3 can yield better diagnostic results in terms of the amplitude of faulty signal.Figure 7
Case study 1: the captured time domain vibration signals from S1 to S4 in (a) to (d), respectively, and their corresponding spectrums in (e) to (h).The corresponding envelope spectrums of time domain signals in Figure7 are illustrated, respectively, in Figure 8. Since the variations in shaft velocities are very low, that is, shown in Figure 5, the smearing effect in envelope spectrums is negligible. However, as expected, the fault frequencies are highly modulated with torsional vibrations, and it is more obvious in signal captured from S3. It is seen that all envelope and envelope order spectrums can reveal the fault frequency. Note that the results presented in Figure 8 make these findings more obvious where the diagnostic information obtained in envelope spectrums is very sharp and obvious for S3, and the fault frequencies have small magnitudes in S1.Figure 8
Envelope spectrums (a–d) and envelope order spectrums (e–h) for the time domain signals in Figure7, respectively, from (a) to (d).In this case, it is seen that the envelope analysis, briefly reviewed and explained in Section2, is capable of detecting the fault frequency, taking into account the fact that the magnitudes are greatly influenced by related transfer path and system structure. This suggests that, using a proper FRF function, it is possible to select an ideal sensing location for monitoring and test purposes.


### 4.2. Case Study 2: Multiple Faulty Bearings on Different Shafts

The analyses carried out for the case study 1 are repeated here for the case study 2, but this time, as shown in Figure9, two excitation forces using the excitation function, fe1, having magnitudes of 120 N and 200 N are applied, respectively, to B1 and B3. Again, from frequency domain point of view, both excitations carry a high level of energy in low frequency band, but low level of energy in high frequency range (i.e., higher than 1 kHz).Figure 9
Illustration of the excitation forces to B1 and B3 in the time domain in (a) and (b), and in the frequency domain in (c) and (d), respectively.Figure10 shows the captured vibration signals, which carry fault diagnostic information for two distinct bearings. The great influence of the sensor location on the fault detectability feature of the signals is very obvious. In fact, the amplitude of the impact trains, as the symptom of existence of bearing faults, in signal captured using S3 is approximately three times greater than that the one obtained using S4 and much sharper than other two sensors. This is an important result, showing that the performance of a CM system is also highly dependent on the sensor location.Figure 10
Case study 2: the captured time domain vibration signals from S1 to S4 in (a) to (d), respectively, and their corresponding spectrums in (e) to (h).The comments made before about the location of sensor on detecting feature of the envelope spectrum are also valid here. Once again, the highest magnitudes in the fault frequencies appear in the envelope spectrums obtained from S3. However, it is quite obvious that the diagnostic results for B3 are not as good as those of for B1.Similar to the previous case study, the envelope order spectrums using the related rpm info are computed and presented in Figure11. The resampling process from time domain to angular domain in order tracking can eliminate the modulation effect of torsional vibrations by removing the side-bands around the fault frequency and its higher harmonics.Figure 11
Envelope order spectrums for B1 (a–d) and for B3 (e–h) for the time domain signals in Figure10, respectively, from (a) to (d).


### 4.3. Case Study 3: Multiple Faulty Bearings on the Same Shaft

In previous case studies, the performance of the envelope analysis is assessed using two distinct bearings. It is shown that the envelope spectrum can yield reliable diagnostic results for cases having well separated fault frequencies. Here, again a case study using multiple faulty bearings, but on the same shaft, is investigated.Figure12 shows the force excitations applied to B1 and B2 having, respectively, magnitudes of 120 N and 80 N. Also, Figure 13 shows the time and frequency domains of the vibration responses obtained from four different locations. Again, similar to the previous cases, in contrast to the applied low external energy level in high frequency range to the system by these forces, inspection of the frequency spectrums in Figure 13 reveals that the system is more sensitive in higher frequencies due to its dynamical characteristics.Figure 12
Illustration of the excitation forces to B1 and B2 in the time domain in (a) and (b), and in the frequency domain in (c) and (d), respectively.Figure 13
Case study 3: the captured time domain vibration signals from S1 to S4 in (a) to (d), respectively, and their corresponding spectrums in (e) to (h).Figure14 shows the envelope spectrums and envelope order spectrums of the time domain signals in Figure 13. As seen, due to overlying of the diagnostic information of B1 and B2, it is not possible to detect and diagnose both the bearings by relying merely on signal frequency domain information. This is an illustration of limitation of envelope spectrum in revealing that the fault frequencies of the same bearings rotate on the same rotating shaft. In the following section, in further discussion, it is shown that this can be explained by the combination of the structural dynamics and frequency domain signal processing.Figure 14
Envelope spectrums (a–d) and envelope order spectrums (e–h) for the time domain signals in Figure13, respectively, from (a) to (d).According to the results calculated and shown in Figure14, it is also can be seen that the signals captured using S3 and S4 share more or less same level of fault magnitudes. This can be explained by existence of a fault in B2, which excites the system in a different frequency range; thus S4 and S3 may have a nearly similar transfer path to the faulty area.Again, the effects of torsional vibrations, as the termS2 in (10), are very clear in terms of side-bands around the bearing fault frequency.


## 4.1. Case Study 1: Single Faulty Bearing

This case study investigates the dynamical behavior of the MBS model having one single faulty bearing on the drive shaft. Figure6 illustrates the artificially generated excitation force applied to B1 in time domain and its representation in frequency domain. As stated earlier, this excitation force is generated using the excitation function, fe1, to analyze the vibration behavior of the model, as well as the real drivetrain, in a condition with a bearing having an initial localized fault such as a tiny crack. The magnitude of force acting on the outer race of B1 with a certain repetition frequency (i.e., BPFO) is assumed to be 120 N in this case study.Figure 6
Illustration of the applied excitation force to B1 in (a) time domain and (b) frequency domain.
(a) (b)Figure7 shows the signals captured from four different locations in the time domain as well as their spectrums up to the Nyquist frequency. As expected, the excitation force excites the system resonances in a specific frequency range. However, by analyzing the force spectrum in Figure 6(b) and the obtained spectrums in Figure 7, it is revealed that small amount of energy in a high frequency range strongly excites the system resonances. This behavior may be explained by transfer path of the vibration from B1 to the sensing locations. In the next section, it is shown that using this valuable insight, an ideal sensor location (between four sensing locations) can be selected. Again, the results presented here suggest that S3 can yield better diagnostic results in terms of the amplitude of faulty signal.Figure 7
Case study 1: the captured time domain vibration signals from S1 to S4 in (a) to (d), respectively, and their corresponding spectrums in (e) to (h).The corresponding envelope spectrums of time domain signals in Figure7 are illustrated, respectively, in Figure 8. Since the variations in shaft velocities are very low, that is, shown in Figure 5, the smearing effect in envelope spectrums is negligible. However, as expected, the fault frequencies are highly modulated with torsional vibrations, and it is more obvious in signal captured from S3. It is seen that all envelope and envelope order spectrums can reveal the fault frequency. Note that the results presented in Figure 8 make these findings more obvious where the diagnostic information obtained in envelope spectrums is very sharp and obvious for S3, and the fault frequencies have small magnitudes in S1.Figure 8
Envelope spectrums (a–d) and envelope order spectrums (e–h) for the time domain signals in Figure7, respectively, from (a) to (d).In this case, it is seen that the envelope analysis, briefly reviewed and explained in Section2, is capable of detecting the fault frequency, taking into account the fact that the magnitudes are greatly influenced by related transfer path and system structure. This suggests that, using a proper FRF function, it is possible to select an ideal sensing location for monitoring and test purposes.


## 4.2. Case Study 2: Multiple Faulty Bearings on Different Shafts

The analyses carried out for the case study 1 are repeated here for the case study 2, but this time, as shown in Figure9, two excitation forces using the excitation function, fe1, having magnitudes of 120 N and 200 N are applied, respectively, to B1 and B3. Again, from frequency domain point of view, both excitations carry a high level of energy in low frequency band, but low level of energy in high frequency range (i.e., higher than 1 kHz).Figure 9
Illustration of the excitation forces to B1 and B3 in the time domain in (a) and (b), and in the frequency domain in (c) and (d), respectively.Figure10 shows the captured vibration signals, which carry fault diagnostic information for two distinct bearings. The great influence of the sensor location on the fault detectability feature of the signals is very obvious. In fact, the amplitude of the impact trains, as the symptom of existence of bearing faults, in signal captured using S3 is approximately three times greater than that the one obtained using S4 and much sharper than other two sensors. This is an important result, showing that the performance of a CM system is also highly dependent on the sensor location.Figure 10
Case study 2: the captured time domain vibration signals from S1 to S4 in (a) to (d), respectively, and their corresponding spectrums in (e) to (h).The comments made before about the location of sensor on detecting feature of the envelope spectrum are also valid here. Once again, the highest magnitudes in the fault frequencies appear in the envelope spectrums obtained from S3. However, it is quite obvious that the diagnostic results for B3 are not as good as those of for B1.Similar to the previous case study, the envelope order spectrums using the related rpm info are computed and presented in Figure11. The resampling process from time domain to angular domain in order tracking can eliminate the modulation effect of torsional vibrations by removing the side-bands around the fault frequency and its higher harmonics.Figure 11
Envelope order spectrums for B1 (a–d) and for B3 (e–h) for the time domain signals in Figure10, respectively, from (a) to (d).


## 4.3. Case Study 3: Multiple Faulty Bearings on the Same Shaft

In previous case studies, the performance of the envelope analysis is assessed using two distinct bearings. It is shown that the envelope spectrum can yield reliable diagnostic results for cases having well separated fault frequencies. Here, again a case study using multiple faulty bearings, but on the same shaft, is investigated.Figure12 shows the force excitations applied to B1 and B2 having, respectively, magnitudes of 120 N and 80 N. Also, Figure 13 shows the time and frequency domains of the vibration responses obtained from four different locations. Again, similar to the previous cases, in contrast to the applied low external energy level in high frequency range to the system by these forces, inspection of the frequency spectrums in Figure 13 reveals that the system is more sensitive in higher frequencies due to its dynamical characteristics.Figure 12
Illustration of the excitation forces to B1 and B2 in the time domain in (a) and (b), and in the frequency domain in (c) and (d), respectively.Figure 13
Case study 3: the captured time domain vibration signals from S1 to S4 in (a) to (d), respectively, and their corresponding spectrums in (e) to (h).Figure14 shows the envelope spectrums and envelope order spectrums of the time domain signals in Figure 13. As seen, due to overlying of the diagnostic information of B1 and B2, it is not possible to detect and diagnose both the bearings by relying merely on signal frequency domain information. This is an illustration of limitation of envelope spectrum in revealing that the fault frequencies of the same bearings rotate on the same rotating shaft. In the following section, in further discussion, it is shown that this can be explained by the combination of the structural dynamics and frequency domain signal processing.Figure 14
Envelope spectrums (a–d) and envelope order spectrums (e–h) for the time domain signals in Figure13, respectively, from (a) to (d).According to the results calculated and shown in Figure14, it is also can be seen that the signals captured using S3 and S4 share more or less same level of fault magnitudes. This can be explained by existence of a fault in B2, which excites the system in a different frequency range; thus S4 and S3 may have a nearly similar transfer path to the faulty area.Again, the effects of torsional vibrations, as the termS2 in (10), are very clear in terms of side-bands around the bearing fault frequency.


## 5. Further Discussions

### 5.1. Ideal Sensor Placement

As briefly discussed in introduction, in reality, the number of used sensors is limited due to various reasons, and this may result in some restrictions. As a result, optimization of the location of sensor(s) becomes a necessary and crucial step in practice for a given problem. The problem of determining the ideal or optimal sensor placements has already received much attention. Most of them share a common basis and deal with a specific application other than condition monitoring and fault detection such as modal testing and control engineering. According to the findings in the previous section, this section is concerned with the specific question of ideal sensor placement for the purpose of bearing fault detection. It should be noted that here, unlike a classical optimization problem, the sensing positions are compared and investigated in terms of their fault detectability.As discussed and shown in theory, the transfer path of the bearing vibration signal plays a major role in the bearing fault detection and has its effects on the performance of a CM system. These effects can be analyzed in two different aspects: (I) the influence of transfer path on diagnostic information, and (II) estimation of transfer path in the desired frequency range for the bandpass filtering prior to envelope analysis.Figure15 illustrates the FRFs obtained from four different locations, namely, S1, S2, S3, and S4, based on two distinct excitation points, on B1 and B3. It is well-known that the FRFs can represent the transfer path in a mechanical system. With reference to the earlier results in previous section, it is expected that the FRFs obtained from S3 have the highest peak in a frequency range of 4.5 kHz to 6.5 kHz, where the fault excites the system resonances. In fact, close inspection of results in Figure 15, in particular comparing envelope spectrums in Figures 8 and 11, reveals that the FRFs, as a transfer path, have a great influence on the amplitudes of generated impact trains, and thus detectability of the corresponding fault frequencies. It is clearly seen that the modulus of the FRF from B1 to S3 in resonance range is approximately ten times greater than other FRFs. This suggests that S3 can be the ideal location among other sensing locations investigated in this study. According to the FRFs in Figure 15(b), similar increasing trend, but this time with a less difference compared to other FRFs, can also be seen for B3.Figure 15
Simulated FRFs obtained in four different sensor locations, exciting from (a) B1 and (b) B3.
(a) (b)In case the experimental FRFs are desired for a similar purpose, with a reasonable assumption, a nearby point can be selected for excitation. In reality, to obtain the experimental FRFs, it is not always possible to excite the system from a point on (or close to) the bearing.


### 5.2. Estimation of Bearing Fault Location

As shown in case study 3, the envelope spectrum is unable to detect the fault frequencies of multiple same-type faulty bearings rotating at the same angular velocity. It is shown here that using the structural information obtained from dynamics of the system, it becomes possible to estimate the probable fault location(s).Figure16(a) presents an overlaid plot for power spectrums for faulty B1, from S3, using fe1 function, each of which with different amount of force magnitude. As seen, all the fault stages in the same bearing excite the same frequency range, but with different magnitudes. In addition, Figure 16(b) presents the power spectrum for case study 3, where both B1 and B2 are faulty. These results in Figure 16 confirm that a localized fault on a specific component in a total system may excite the resonances depending on neighboring components and system dynamics.Figure 16
The power spectrum of the captured signal from S3 for (a) single faulty bearing (b) multiple faulty bearings.
(a) (b)Similar to the findings in Figure15 for B1, while a dominant frequency range (i.e., between 4.5 kHz and 6.5 kHz) can be seen for different cases having only one faulty bearing (i.e., B1), some other additional dominant frequency ranges appear in the power spectrum for a case having multiple faulty bearings (i.e., B1 and B2). Again, this is an important insight, in particular for the cases where the faulty bearings share same geometry and rotate at the same angular velocity, where envelope spectrum may not be able to diagnose both bearings. To have a better understanding of system resonance areas and the effects of fault excitations due to the localized bearing faults, Figure 17 illustrates the Campbell diagrams, representing both time and frequency domains for case studies 1 and 3. According to the results in Figure 17(a), although the faulty bearing, B1, may excite the total system in a broadband frequency range (with a specific repetition rate), this excitation reaches to its highest amplitude within a frequency range between 4.5 kHz and 6.5 kHz for case study 1. Similar insights are obtained in Figure 17(b) for case study 3, but this time two similar faulty bearings rotating at the same speed, B1 and B2, produce highest excitations in distinct frequency ranges. The comparison of results of both case studies here may confirm the existence of another bearing fault with the same repetition frequency, which may lead to a conclusion of existence a fault on B2. Note that, in practice, mechanical faults may lead to excitation of resonances in different frequency bands; however, the numerically generated FRFs excited from a point on the B2, which are shown in Figure 18, suggest that this specific system shows its highest sensitivities to the excitations from B2 in two frequency ranges, one around 4 kHz and one around 8 kHz, for S3 and S4. Similar results are previously shown in Figure 14 for case study 3.Figure 17
Time-frequency representations of the captured signals from S3 for (a) case study 1: single faulty bearing and (b) case study 3: multiple faulty bearings.
(a) (b)Figure 18
Simulated FRFs obtained in four different sensor locations exciting from B2.Until now, initial fault stages on bearings and their effects on vibration trend of the system are investigated. Here, using other two excitation functions, shown in Figure2, two other fault stages are numerically analyzed in terms of their envelope spectrums. Figure 19 shows an overlaid plot for two envelope spectrums obtained from S3, as a desired sensor position. This is a comparison of a semilocalized fault and a semidistributed fault. As seen, while the fault frequency, BPFO, and its harmonics are visible and sharp for a medium-size fault, they become subtle and difficult to detect for a large-size bearing fault with a more distributed faulty area. This is an expected result, where in reality, when bearing faults grow, they may lose their sharpness by passing the roller elements over time. This example once again illustrates the great importance of detection of bearing faults as early as possible to avoid resulting complications.Figure 19
Envelope spectrums for two various excitations.


## 5.1. Ideal Sensor Placement

As briefly discussed in introduction, in reality, the number of used sensors is limited due to various reasons, and this may result in some restrictions. As a result, optimization of the location of sensor(s) becomes a necessary and crucial step in practice for a given problem. The problem of determining the ideal or optimal sensor placements has already received much attention. Most of them share a common basis and deal with a specific application other than condition monitoring and fault detection such as modal testing and control engineering. According to the findings in the previous section, this section is concerned with the specific question of ideal sensor placement for the purpose of bearing fault detection. It should be noted that here, unlike a classical optimization problem, the sensing positions are compared and investigated in terms of their fault detectability.As discussed and shown in theory, the transfer path of the bearing vibration signal plays a major role in the bearing fault detection and has its effects on the performance of a CM system. These effects can be analyzed in two different aspects: (I) the influence of transfer path on diagnostic information, and (II) estimation of transfer path in the desired frequency range for the bandpass filtering prior to envelope analysis.Figure15 illustrates the FRFs obtained from four different locations, namely, S1, S2, S3, and S4, based on two distinct excitation points, on B1 and B3. It is well-known that the FRFs can represent the transfer path in a mechanical system. With reference to the earlier results in previous section, it is expected that the FRFs obtained from S3 have the highest peak in a frequency range of 4.5 kHz to 6.5 kHz, where the fault excites the system resonances. In fact, close inspection of results in Figure 15, in particular comparing envelope spectrums in Figures 8 and 11, reveals that the FRFs, as a transfer path, have a great influence on the amplitudes of generated impact trains, and thus detectability of the corresponding fault frequencies. It is clearly seen that the modulus of the FRF from B1 to S3 in resonance range is approximately ten times greater than other FRFs. This suggests that S3 can be the ideal location among other sensing locations investigated in this study. According to the FRFs in Figure 15(b), similar increasing trend, but this time with a less difference compared to other FRFs, can also be seen for B3.Figure 15
Simulated FRFs obtained in four different sensor locations, exciting from (a) B1 and (b) B3.
(a) (b)In case the experimental FRFs are desired for a similar purpose, with a reasonable assumption, a nearby point can be selected for excitation. In reality, to obtain the experimental FRFs, it is not always possible to excite the system from a point on (or close to) the bearing.


## 5.2. Estimation of Bearing Fault Location

As shown in case study 3, the envelope spectrum is unable to detect the fault frequencies of multiple same-type faulty bearings rotating at the same angular velocity. It is shown here that using the structural information obtained from dynamics of the system, it becomes possible to estimate the probable fault location(s).Figure16(a) presents an overlaid plot for power spectrums for faulty B1, from S3, using fe1 function, each of which with different amount of force magnitude. As seen, all the fault stages in the same bearing excite the same frequency range, but with different magnitudes. In addition, Figure 16(b) presents the power spectrum for case study 3, where both B1 and B2 are faulty. These results in Figure 16 confirm that a localized fault on a specific component in a total system may excite the resonances depending on neighboring components and system dynamics.Figure 16
The power spectrum of the captured signal from S3 for (a) single faulty bearing (b) multiple faulty bearings.
(a) (b)Similar to the findings in Figure15 for B1, while a dominant frequency range (i.e., between 4.5 kHz and 6.5 kHz) can be seen for different cases having only one faulty bearing (i.e., B1), some other additional dominant frequency ranges appear in the power spectrum for a case having multiple faulty bearings (i.e., B1 and B2). Again, this is an important insight, in particular for the cases where the faulty bearings share same geometry and rotate at the same angular velocity, where envelope spectrum may not be able to diagnose both bearings. To have a better understanding of system resonance areas and the effects of fault excitations due to the localized bearing faults, Figure 17 illustrates the Campbell diagrams, representing both time and frequency domains for case studies 1 and 3. According to the results in Figure 17(a), although the faulty bearing, B1, may excite the total system in a broadband frequency range (with a specific repetition rate), this excitation reaches to its highest amplitude within a frequency range between 4.5 kHz and 6.5 kHz for case study 1. Similar insights are obtained in Figure 17(b) for case study 3, but this time two similar faulty bearings rotating at the same speed, B1 and B2, produce highest excitations in distinct frequency ranges. The comparison of results of both case studies here may confirm the existence of another bearing fault with the same repetition frequency, which may lead to a conclusion of existence a fault on B2. Note that, in practice, mechanical faults may lead to excitation of resonances in different frequency bands; however, the numerically generated FRFs excited from a point on the B2, which are shown in Figure 18, suggest that this specific system shows its highest sensitivities to the excitations from B2 in two frequency ranges, one around 4 kHz and one around 8 kHz, for S3 and S4. Similar results are previously shown in Figure 14 for case study 3.Figure 17
Time-frequency representations of the captured signals from S3 for (a) case study 1: single faulty bearing and (b) case study 3: multiple faulty bearings.
(a) (b)Figure 18
Simulated FRFs obtained in four different sensor locations exciting from B2.Until now, initial fault stages on bearings and their effects on vibration trend of the system are investigated. Here, using other two excitation functions, shown in Figure2, two other fault stages are numerically analyzed in terms of their envelope spectrums. Figure 19 shows an overlaid plot for two envelope spectrums obtained from S3, as a desired sensor position. This is a comparison of a semilocalized fault and a semidistributed fault. As seen, while the fault frequency, BPFO, and its harmonics are visible and sharp for a medium-size fault, they become subtle and difficult to detect for a large-size bearing fault with a more distributed faulty area. This is an expected result, where in reality, when bearing faults grow, they may lose their sharpness by passing the roller elements over time. This example once again illustrates the great importance of detection of bearing faults as early as possible to avoid resulting complications.Figure 19
Envelope spectrums for two various excitations.


## 6. Concluding Remarks and Future Works

### 6.1. Concluding Remarks

(i) In the present paper, in accordance with the purpose of improving capabilities and reliabilities of the present vibration-based condition monitoring and fault detection systems for complex mechanical machines/structures in general, and for drivetrains in automotive engineering in particular, a new approach based on combination of numerical modelling and signal processing techniques is developed and investigated in some detail.(ii) For the objectives of this study, an experimentally validated MBS model of the drivetrain of an electrified vehicle is used as a tool for numerical investigations as well as for illustrating the relations between dynamics of the system and rolling bearing fault signals.(iii) Three various case studies are defined first and investigated in terms of their vibration responses. It is shown and seen that using only frequency domain information of the measured data may lead to some missed detection, for example, in monitoring of the nearby components with similar geometric details.(iv) The benefits and limitations of the well-known envelope analysis in automotive engineering applications are highlighted and discussed. A combination of the fields structural dynamics and signal processing, to overcome the available difficulties in finding a desired sensor placement and fault location(s) for the purpose of condition monitoring, is developed. It is shown that using FRFs, as the transfer path functions, delivers a new degree-of-freedom for CM systems. These insights may even be beneficial and useful for further process such as source separation purposes.(v) Finally, the results from this study show that the approach of analyzing and fault identifying the complex systems using a corresponding MBS model can also be employed as a quality control tool for manufacturers for testing purposes.


### 6.2. Future Works

In terms of validation and more detailed investigations, the proposed and used methodology needs to be studied and mature further for other complex systems with different faulty cases (e.g., gear stages). In this regard, the present study should be considered as a basis for more comprehensive studies on combination of the fields vibration-based fault diagnosis and numerical modelling. Future works will study a mechanical system having a wide range of mechanical faults.


## 6.1. Concluding Remarks

(i) In the present paper, in accordance with the purpose of improving capabilities and reliabilities of the present vibration-based condition monitoring and fault detection systems for complex mechanical machines/structures in general, and for drivetrains in automotive engineering in particular, a new approach based on combination of numerical modelling and signal processing techniques is developed and investigated in some detail.(ii) For the objectives of this study, an experimentally validated MBS model of the drivetrain of an electrified vehicle is used as a tool for numerical investigations as well as for illustrating the relations between dynamics of the system and rolling bearing fault signals.(iii) Three various case studies are defined first and investigated in terms of their vibration responses. It is shown and seen that using only frequency domain information of the measured data may lead to some missed detection, for example, in monitoring of the nearby components with similar geometric details.(iv) The benefits and limitations of the well-known envelope analysis in automotive engineering applications are highlighted and discussed. A combination of the fields structural dynamics and signal processing, to overcome the available difficulties in finding a desired sensor placement and fault location(s) for the purpose of condition monitoring, is developed. It is shown that using FRFs, as the transfer path functions, delivers a new degree-of-freedom for CM systems. These insights may even be beneficial and useful for further process such as source separation purposes.(v) Finally, the results from this study show that the approach of analyzing and fault identifying the complex systems using a corresponding MBS model can also be employed as a quality control tool for manufacturers for testing purposes.


## 6.2. Future Works

In terms of validation and more detailed investigations, the proposed and used methodology needs to be studied and mature further for other complex systems with different faulty cases (e.g., gear stages). In this regard, the present study should be considered as a basis for more comprehensive studies on combination of the fields vibration-based fault diagnosis and numerical modelling. Future works will study a mechanical system having a wide range of mechanical faults.



---

*Source: 5159189-2018-03-13.xml*