Next Article in Journal
Bio-Inspired Optimization-Based Path Planning Algorithms in Unmanned Aerial Vehicles: A Survey
Next Article in Special Issue
Weighted Domain Adaptation Using the Graph-Structured Dataset Representation for Machinery Fault Diagnosis under Varying Operating Conditions
Previous Article in Journal
Back-Projected Signal-Based Self-Interferometric Phase Analysis Technique for Sea Surface Observation Using a Single Scatterometer System
Previous Article in Special Issue
Remaining Useful Life Prediction Based on Adaptive SHRINKAGE Processing and Temporal Convolutional Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Acoustic-Based Rolling Bearing Fault Diagnosis Using a Co-Prime Circular Microphone Array

1
School of Mechanical Engineering, Shenyang University of Technology, Shenyang 110178, China
2
School of Mechanical Engineering, Shenyang Ligong University, Shenyang 110159, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(6), 3050; https://doi.org/10.3390/s23063050
Submission received: 9 February 2023 / Revised: 2 March 2023 / Accepted: 9 March 2023 / Published: 12 March 2023
(This article belongs to the Special Issue Fault Diagnosis and Prognosis in Rotating Machines)

Abstract

:
This study proposes a high-efficiency method using a co-prime circular microphone array (CPCMA) for the bearing fault diagnosis, and discusses the acoustic characteristics of three fault-type signals at different rotation speeds. Due to the close positions of various bearing components, radiation sounds are seriously mixed, and it is challenging to separate the fault features. Direction-of-arrival (DOA) estimation can be used to suppress noise and directionally enhance sound sources of interest; however, classical array configurations usually require a large number of microphones to achieve high accuracy. To address this, a CPCMA is introduced to raise the array’s degrees of freedom in order to reduce the dependence on the microphone numbers and computation complexity. The estimation of signal parameters via rotational invariance techniques (ESPRIT) applied to a CPCMA can quickly figure out the DOA estimation without any prior knowledge. By using the techniques above, a sound source motion-tracking diagnosis method is proposed according to the movement characteristics of impact sound sources for each fault type. Additionally, more precise frequency spectra are obtained, which are used in combination to determine the fault types and locations.

1. Introduction

Industrial equipment nowadays has become increasingly large scale and complex. On the other hand, the concept of intelligent manufacturing involves ongoing improvement of equipment performance and resource efficiency. As a result, timely and accurate machine health monitoring plays a far more important role [1]. Bearing failure, which is generally brought on by improper mounting or inappropriate lubrication, is one of the most common causes of a mechanical breakdown, and has consequently led to a remarkable amount of continuous studies [2]. Significant progress has been made in fault diagnosis based on vibration signals employing methods, including time/frequency domain analysis, modal analysis, finite element analysis, etc. [3,4,5]. Techniques such as wavelet transform (WT) and mode-decomposition methods are proposed and effectively extract the information from non-stationary, non-linear signals [6,7].
However, vibration diagnosis still has limitations and cannot work well in cases such as slow rotating and complex dynamics [8]. Acoustic-based diagnosis has gradually developed as a new effective approach. Diagnosis using acoustic emissions (AE) exploits the transient elastic waves when deformation occurs within a material, from which the energy loss could span a wide range of frequencies. AE sensors are able to capture much higher frequencies than vibration sensors, so the technique is more sensitive to early faults [9,10,11,12]. However, AE waves are greatly attenuated during propagation; therefore, just like vibration measurement, AE sensors should be placed as close as possible to the components being tested [13]. Detection based on acoustic radiation signals is another acoustic-based method, which offers the special advantage of non-contact signal acquisition, which is particularly applicable to occasions where sensor attachment is problematic. It also avoids downtime and operational inconveniences induced by sensor installation and maintenance for large equipment [14,15]. Some studies have shown that for certain statistical parameters, acoustic signals may represent clearer features for early defects than vibration signals in rotating machinery fault diagnosis [16,17]. Moreover, acoustic signals are often used in motor fault diagnosis, such as permanent-magnet demagnetization [18], commutator motor with rotor coil or gear defects [19], rotation speed identification, etc. [20]. However, the signal obtained from a single-channel microphone only provides sound pressure values; hence, acoustic measurement is highly sensitive to measuring points and the features of interest are very likely to be obscured by high-level noise, which explains the relative lack of investigation and application in this field [21,22].
Signal processing techniques based on multi-channel microphones have been presented as a solution to the issue. Spatial distribution information of the acoustic field around the sound sources can be reconstructed through the specific positional relationship between microphones, which allows for the directional enhancement and denoising of signals even from the moving sound sources. Near-field acoustic holography (NAH) and beamforming based on far-field signal models are two efficiently used detection techniques. According to Lu [23] and Hou [24], gray-level co-occurrence matrix features were retrieved from sound field images for bearing fault pattern identification. Chen adopted cyclic spectral density as the reconstruction indicator instead of the complex sound pressure to improve the efficiency of the holographic data processing and precisely identified the sources of the stationary sound field [25]. Nejade proposed the multi-reference holography processing method based on residual spectra using a simple structure model and examined its reliability with two industrial machines [26]. Ma developed a vibroacoustic transfer function using the finite element method for multiple rotating components enclosed in a small close space, which was proven to be effective by near-field experiments [27]. Deep learning methods were used in recent works to improve the accuracy of NAH [28,29]. Cabada et al. combined spectral kurtosis (SK) with beamforming to evaluate both the locations and frequencies of bearing failures using a spiral microphone array [30]. Verellen studied the beamformed acoustic data with a variety of signal-to-noise ratios as a reference for ultrasonic fault detection [31]. One of the most typical uses for multi-channel acoustic diagnosis is wayside identification, where researchers have invested large efforts to overcome signal distortion due to the Doppler effect and noise interference [32,33,34,35]. The studies mentioned above basically adopt two research ideas. One is to localize the possible position areas of the defect sound sources, and the other is to filter the signal through directional enhancement. However, microphone arrays provide the characteristic of spatial domain information, which is rarely used in the field of fault diagnosis.
The main purpose of the study is to propose a novel sound source motion-tracking approach for fault-type recognition. Spatial domain information is fully made use of, which provides a fresh perspective on fault diagnosis as an addition to the established analysis methods based on time and frequency domains. Sound sources produced by faults for various components and types are distinct from each other in obedience to the operating mechanism, which gives more logical proofs for diagnosis using motion trajectories indicated by direction-of-arrival (DOA) estimates. DOA estimates are also used for the signal beamforming, and frequency spectra with good noise reduction and obvious fault features can be extracted. The motion trajectories and the spectra work together to identify the fault types, as well as the positions; accordingly, the characteristics of the bearing fault acoustic signals are further observed and analyzed. However, precise DOA estimation of sound sources coming from large equipment often needs to be located in both the vertical and horizontal dimensions. For this reason, microphone configurations such as planar or circular arrays are customarily adopted, and a large number of microphones are used to achieve good accuracy; therefore, the practical applicability of the technique is constrained and subject to an expensive cost and high computational time. The study presents an adaptive array signal processing method that provides high-accuracy direction-of-arrival (DOA) estimation and fault feature extraction with fewer array elements, in order to reduce the reliance on the number of microphones in real applications. For the proposed method, a co-prime circular array (CPCA) configuration is attempted for signal acquisition, and the estimation of signal parameters via the rotational invariance techniques (ESPRIT) algorithm is used for the data processing. The peculiar structure of a co-prime array was first proposed and validated by Xiang and Bush in 2015, which makes it possible to construct virtual array elements. By exploiting the second-order statistics, more degrees of freedom are acquired, which overcomes the Nyquist theorem’s restriction on array apertures [36,37,38,39]. Meanwhile, when using a circular array, ESPRIT exploits the rotational invariance property to obtain paired two-dimensional DOA estimation without any prior knowledge, which avoids the secondary peak-search and accelerates the process speed [40,41,42].
In this paper, acoustic signals from a rolling bearing test rig are collected using a uniform circular microphone array (UCMA), as well as a co-prime circular microphone array (CPCMA) simultaneously, where exactly the same number of microphones are applied in both arrays. Three test conditions are investigated, including outer race fault, inner race fault, and rolling element fault. Signals between the two arrays are observed and compared after the ESPRIT processing, and the effectiveness in fault diagnosis is evaluated.

2. Theoretical Background

2.1. Co-Prime Circular Array

Co-prime array configuration was at first applied in linear arrays. Let M and N be a pair of co-prime numbers and M < N, then consider two sparse uniform linear sub-arrays with M and N elements, respectively. Overlap the two sub-arrays, and an M + N − 1 element co-prime linear array can be created, which was validated by experiments to have comparable degrees of freedom to a uniform linear array of MN elements [36,37]. The co-prime circular array model is set up based on the co-prime linear array theory. Consider two uniform circular arrays (UCAs) with the same radius as two sub-arrays, in which there are M and N elements, respectively. To form a CPCA, interleave the two sub-arrays. The first element of the new array is shared by the two sub-arrays, as illustrated in Figure 1, where the corresponding array elements of sub-arrays are severally represented by hollow asterisks and filled circles. The number of CPCA array elements is M + N − 1.
Figure 2 shows the CPCA signal model. Assume a CPCA of radius R in the x y plane with its elements clockwise distributed over the circumference. Let the origin O of the coordinate system be the center of the array. A narrowband plane wave comes from the sound source S, and S’ is the projection point of S. Take the radius overlaid with the x-axis as the reference line passing through Element 1. The angle measured down from the z-axis to OS is the elevation angle θ of the sound arrival direction, and the angle measured from the x-axis to OS’ is the azimuth angle ϕ . γ k is the angle between the first and the kth element.
The array response vector a ( θ , ϕ ) can be formulated as:
a ( θ , ϕ ) = [ exp ( j k R sin θ cos ( ϕ γ 0 ) exp ( j k R sin θ cos ( ϕ γ 1 ) exp ( j k R sin θ cos ( ϕ γ M + N 2 ) ]
When there are P uncorrelated sources, we have:
A = [ a ( θ 1 , ϕ 1 ) , , a ( θ P , ϕ P ) ]
The received signal then can be written as:
x ( t ) = A s ( t ) + n ( t )
where s ( t ) and n ( t ) , respectively, stand for the matrices of the incoming source signal and noise signal.
Then, the covariance matrix of the signal can be obtained as:
R xx = E [ x ( t ) x H ( t ) ]
For an M + N−1 element UCA, the angle γ k = 2 π i / ( m + n 1 ) , i = 0 , 1 , , m + n 2 is easily calculated; however, arc lengths between elements in a CPCA are uneven. To solve this, vectorize R xx as:
v e c ( R xx ) = [ a * ( θ 1 , ϕ 1 ) a ( θ 1 , ϕ 1 ) , , a * ( θ P , ϕ P ) a ( θ P , ϕ P ) ] p + σ n 2 v e c ( I )
where vector p denotes the power squared of the P sources, σ n 2 is the noise power, and I is the identity matrix. The vectorization uses difference co-array (DCA) to reconstitute a virtual UCA with MN − 1 uniformly distributed elements, which notably improves the degrees of freedom compared to an M + N − 1 physical element UCA.

2.2. UCA-ESPRIT

ESPRIT was initially presented based on uniform linear arrays (ULA). ESPRIT exploits the rotational invariance of the signal subspace; accordingly, the channel manifold vectors should follow the Vandermonde structure, while the UCA signal model does not meet the condition. Transformation based on phase mode excitation converts a UCA manifold vector into a ULA-like Vandermonde matrix by using the recursive nature of Bessel functions. Consider F r H as the beamformer to map the UCA manifold vector a ( θ , ϕ ) onto the transformed beamspace, whose manifold is defined as b ( θ , ϕ ) . The mapping can be expressed as:
b ( θ , ϕ ) = F r H a ( θ , ϕ )
To keep F r H orthogonal, a rotation angle α is brought in to construct a matrix W that has the structure of centro-Hermitian rows, which is expressed as:
W = 1 K [ v ( α K ) , v ( α 0 ) , , v ( α K ) ]
where α i = 2 π i / K , i [ K , K ] . K denotes the highest excited mode passing through the array apertures. v ( ϕ ) is the vector to convert the azimuth variation into the Vandermonde-like structure, which is given as:
v ( ϕ ) = [ e - j K ϕ , , e - j ϕ , e - j 0 , e j ϕ , , e j K ϕ ] T
Elevation information of the manifold also needs to be converted into a symmetric amplitude taper, by giving the matrix C V V H that derives from the Bessel function. For a UCA with N elements, we have:
C V = d i a g { j K , , j 1 , j 0 , j 1 , , j K }
V = N [ w K , w 0 , , w K ]
where w i , i [ K , K ] is the normalized weight vector for the ith phase mode, which is expressed as:
w i = 1 N [ 1 , e - j 2 π i / N , , e - j 2 π i ( N - 1 ) / N ] T
Combined with those aforementioned, the beamformer F r H can be given as:
F r H = W H C V V H
and the received signal can be converted as:
y ( t ) = F r H x ( t )
By considering the phase delay between two adjacent array elements as an artificially defined delay matrix, the array can be split into identical sub-arrays. According to the rotational invariance, the sub-arrays share the same signal subspace, which can be obtained by eigenvalue decomposition of the array signal covariance matrix. Combine the sub-array signals into a matrix z ( t ) , the signal subspace S can be calculated by applying eigenvalue decomposition to R zz = E [ z ( t ) z H ( t ) ] and the eigenvalues of the phase delay matrix Φ can be picked out as well, which gives the paired DOA information θ and ϕ . For P sound sources, there is:
Φ = diag ( sin θ 1 e j ϕ 1 , sin θ 2 e j ϕ 2 , , sin θ P e j ϕ P )

3. Experimental Investigation

3.1. Experimental Setup

The experiment is conducted on a QPZZ-II test rig, as shown in Figure 3. The rotor apparatus is driven by an ac motor. The rotation speed is controlled by the generator frequency, which ranges from 75 rpm to 1450 rpm. The bearings (HRB N205EM CHINA) used in the tests are shown in Figure 4. Three test conditions are investigated, including outer race fault, inner race fault, and rolling element fault. When conducting the outer race fault, the outer race breach is mounted at the 5 o’clock position. For each trial, the apparatus operates under constant speeds at 290 rpm and 725 rpm, which corresponds to the generator frequency of 10 Hz and 25 Hz, and the shaft rotational speed of 4.83 Hz and 12.08 Hz. The outer race keeps fixed. Table 1 and Table 2 represent the bearing configuration and the fault feature frequencies.
The microphone array is placed at a 170 mm distance in front of the apparatus, as is shown in Figure 5. The diameter of the circumference is 200 mm. Figure 6 shows the microphone arrangement. A total of nine microphones are used in the measurement. The microphone designated Number 0 is situated at the center to capture mono signals. Numbers 2, 3, 4, 6, 7, and 8 comprise a 6-element UCA, as well as a 3-element co-prime circular sub-array, while Numbers 1, 3, 5, and 7 form another 4-element co-prime circular sub-array, consequently both UCMA and CPCMA have 6 array elements. The setup provides simultaneous signal acquisition of UCMA and CPCMA, which minimizes the effect of environmental variations on the signals collected by both. The data collector is NI PXIe-1082 and the type of microphone is G.R.A.S. 40 ph. For each trial, the signals are recorded at a sample rate of 10.24 kHz during 10 s.

3.2. Data Analysis Procedure

Firstly, the motion regularity of each fault type needs to be analyzed. The outer race fault produces sounds when balls are passing through the defect, which are always from the same direction. Inner race fault and rolling element fault generate moving sound sources that uniformly rotate with the shaft. Due to the existence of moving sources, the DOA estimation algorithm cannot be applied directly to a long-lasting data section. The collected data need to be chopped into segments in keeping with the real-time location of the defects.
However, it cannot be guaranteed that every split segment contains information of interest. When there is no collision of defects in the segment, unconcerned sound source positions are likely to be located instead. To address the issue, elevation angles are checked to screen unwanted DOA estimates out, for the impact sounds always come from the circumferences where the defects locate, and the elevation angles from the same circumference have similar values. Azimuth angles then are feasible to be observed to speculate about the fault types. Furthermore, spatial filtering is applied on the basis of the screened DOA estimates, and the resulting frequency spectra also help to determine. The signal processing procedure is carried out as follows:
  • Separate the collected signals into UCMA and CPCMA sets.
  • Figure out the rotation angles between every collision for each defect type. Since the outer race is fixed, the azimuth angle values for the outer race defect should always be the same. To calculate the angular interval of the impact sounds for the inner race defect, we have:
Δ ϕ I D = 2 π / ( f B P F I / f i )
where f i is the rotational speed of both the shaft and the inner race.
The rolling element fault has the peculiarity that impact sounds are triggered both on the inner and outer races, and both are moving sound sources, for which we have:
Δ ϕ R D O = π 2 d B / ( π L O )
Δ ϕ R D I = π 2 d B / ( π d I )
where Δ ϕ R D O and Δ ϕ R D I denote the rotation angles between defect collisions on the outer and inner races. d B is the nominal diameter of balls, L O is the outer ring land diameter, and d I is the inside diameter of the inner race bore.
An adequate interval angle Δ ϕ L should be determined to divide the data into segments. It is necessary that Δ ϕ L is smaller than Δ ϕ I D , Δ ϕ R D O , and Δ ϕ R D I , so that the positions of impact sounds can be closely tracked. From Δ ϕ L , the number of snapshots used to split data can be calculated as:
N L = Δ ϕ L L / ( 2 π f i )
where L denotes the number of snapshots per second, which depends on the sampling rate of signal acquisition.
3.
Apply ESPRIT to UCMA and CPCMA signal segments to obtain DOA estimates of both. Elevation angles for each fault type should be worked out based on the distance between the defect and the reference element of the microphone array. When an elevation angle estimate has a large difference from the calculated value in a segment, it means the segment gives invalid information and the location should be excluded.
4.
Compare the azimuth angular intervals of the screened DOA estimates with Δ ϕ I D , Δ ϕ R D O , and Δ ϕ R D I to observe whether it conforms to the moving regularity of impact sounds for a specific bearing fault type.
5.
Apply beamforming for the further determination of fault types.
The flow chart in Figure 7 represents the data processing and analysis procedure.

3.3. Results and Discussion

3.3.1. Outer Race Defect

Table 3 and Table 4 each list a group of screened DOA estimates severally from UCMA and CPCMA signals under the outer race fault condition at 290 rpm. To see clearly, values listed are all expressed in degrees. As mentioned in Section 3.1 and Section 3.2, the breach set on the outer race is placed at a 150-degree angle direction relative to the reference array element, and the elevation angle is calculated as 7.79 degrees. It can be seen that CPCMA has much closer estimates than UCMA. The root-mean-square errors (RMSEs) for each group of values are listed in Table 5.
It is noticed that negative values are occasionally found in the azimuth angle estimates, as has been shown in Table 3. That is because the incident sound waves are presumed to be plane waves in the far-field model, where aliasing in signals may induce a 180-degree phase ambiguity and mislead the estimation results.
Figure 8 shows the spectral envelopes of the non-processed mono signal, the UCMA-ESPRIT beamformed signal, and the CPCMA-ESPRIT beamformed signal. The fault feature frequencies of the outer race defect are marked.
The mono signal is captured by the microphone at the center of the array right in front of the shaft; therefore fault features are greatly masked by the modulation of the shaft rotational frequency f i . The UCMA signal after beamforming presents slightly clearer fault features, but still the modulation of f i causes confusion that cannot be ignored. The beamformed CPCMA signal makes a huge improvement and practically eliminates the interference.

3.3.2. Inner Race Defect

The elevation angle for the inner race defect is calculated as 5.29 degrees. For the inner race fault, Table 6 and Table 7 each list a group of 10 sequential screened DOA estimates at 290 rpm, respectively, from UCMA and CPCMA signals. The trajectory of rotation can be found from the azimuth angle estimates, and in each group, a full rotation is approximately covered. As noted above, values with phase ambiguity likewise come up, which have been marked in the tables.
Table 8 shows RMSE for each group of the estimates. Since the azimuth angle is not a fixed value, the angular interval, which is figured out as 46.42 degrees, is chosen as the observed value in the calculation of RMSE. The results indicate that CPCMA performs better than UCMA at locating, moreover, phase ambiguity appears mainly in UCMA signals, which manifests the same conclusion. However, judging from the elevation angle estimates, both UCMA and CPCMA show a drop in location accuracy when compared with the condition of the outer race fault. It is likely to be induced by the movement of the sound sources.
Figure 9 shows a comparison of the spectral envelopes among the three signal processing approaches.
Similar to the outer race fault circumstance, the fault feature eigenfrequency in the mono signal is overwhelmed by the f i modulation, and the harmonics are barely perceptible. Processed UCMA signal roughly extracts fault feature frequencies, yet the interference is still strong and the relative relationship between harmonic amplitudes appears disordered. The issues are all well resolved in the processed CPCMA signal.

3.3.3. Rolling Element Defect

The rolling element defect successively collides with the outer and the inner races, generating fault locating with two different elevation angles. Based on this, DOA estimates after screening are given in Table 9, Table 10, Table 11 and Table 12. Each table lists 10 segments at 290 rpm. Table 13 and Table 14 give RMSEs in the ball pass location of the outer and inner races, where calculated values of the elevation angles remain the same as before, and angular intervals are figured out as 29.03 degrees for the outer race and 54 degrees for the inner race.
It can be found from the DOA estimates of UCMA that elevation angle values are highly dispersed, and the phase ambiguity occurs more frequently, which shows a further decrease in location accuracy compared with the inner race fault condition. It is probably because the interval time between impacts is so close, which aggravates the aliasing of sound waves and makes locating more difficult. The RMSE values confirm the observation that location errors of UCMA are much larger than that under the inner race fault condition. DOA estimation precision of CPCMA also has a slight decline, which is manifested as a small rise in phase ambiguity occurrence frequency. Regardless of this, CPCMA still maintains an accuracy numerically no less than before.
Figure 10 shows the spectral envelopes of the three signal processing methods, where the rolling element defect frequency (RDF) is equivalent to twice the ball spin frequency (BSF). The spectra of the rolling element fault contain significant fundamental train frequency (FTF) components. Because of the intense interference of its modulation, fault feature frequencies are almost completely submerged and hardly distinguished. UCMA has enhanced the fault features, but due to the inaccuracy of elevation angle estimates, f i is enhanced at the same time, of which the modulation constitutes a new disturbance. In the CPCMA signal, the modulation of both FTF and f i are effectively suppressed as a result of the precise location, and harmonics up to the fourth order are clearly extracted.

3.3.4. Experimental Results at a Different Rotation Speed

In order to compare the diagnostic efficacy at different rotation speeds, experimental results under the outer and inner race fault conditions at 725 rpm are presented to show the location and fault feature extraction effects of both static and moving impact sound sources. For the outer race fault signal, Table 15, Table 16 and Table 17 give the DOA estimates and RMSE, and Figure 11 shows the spectral envelopes of each processing method. Additionally, Table 18, Table 19 and Table 20 and Figure 12 present the results under the inner race fault condition.
The comparison shows that at a higher speed, the characteristics of both defect and interference sound sources are enhanced, which is mainly evidenced by the amplitude increase in the high-order harmonics of the fault feature frequencies, as well as the f i modulation. From the DOA estimates, elevation angle location accuracies for both fault types are all improved at a higher speed, which is caused by the enhancement of the direct sounds from defects to microphones. However, as can be seen from DOA estimates for the outer race fault, both UCMA and CPCMA show an overall deviation in the azimuth angle localization. This is probably because azimuth location is highly influenced by the reflected sounds, and higher sound pressure of the sources leads to stronger reflections and coherence, which results in location bias. On the other hand, the enhancement of the impact sounds has reduced the dispersion in DOA estimates. It can be presumed that the same is true for the inner race fault, but compared with UCMA, the CPCMA signal maintains almost as good accuracy in angular interval estimates as at 290 rpm.

4. Conclusions

In this paper, a novel diagnosis method based on sound source motion tracking is proposed, and acoustic signals of three different bearing fault types at two rotation speeds are acquired to validate its effectiveness. Compared with vibration signals, acoustic signals have significantly more noise and aliasing problems. Due to the close positions of each bearing component, eigenfrequency modulation of the fault-free components, such as the shaft, has a great effect even on the high-order harmonics and in many cases overwhelms the fault feature frequencies. Sounds of the defects can be directionally enhanced by locating the collisions using a microphone array. However, UCMA requires a large number of microphones to achieve a certain location accuracy. CPCMA takes advantage of the non-uniform differential relationship between array elements to create a virtual array through the vectorization process, which considerably raises the degree of freedom of the array and lessens the dependence on the number of microphones needed for high-precision. The superiority of CPCMA performance with the same number of microphones has been verified by experiments in the study. When applying more microphones, the performance gap between the two array configurations will be larger in accordance with CPCA properties. On this basis, using ESPRIT as the signal processing method, the paired elevation and azimuth angles can be quickly estimated, which can largely improve the process speed for array signals. The CPCMA-ESPRIT method realizes the low-cost and high-efficiency sound source location of defects, and satisfies the requirement for intelligent analysis in practical engineering.

Author Contributions

Conceptualization, C.L.; experimental design and implementation, C.L.; data analysis, C.L.; thesis writing, C.L.; supervision, C.C.; funding acquisition, C.C.; facilities and laboratory space management, X.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 51675350, No. 51705337).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhou, J.; Li, P.G.; Zhou, Y.H. Toward New-Generation Intelligent Manufacturing. Engineering 2018, 4, 11–20. [Google Scholar] [CrossRef]
  2. Fu, C.; Sinou, J.J.; Zhu, W. A state-of-the-art review on uncertainty analysis of rotor systems. Mech. Syst. Signal Process. 2023, 183, 109619. [Google Scholar] [CrossRef]
  3. Feng, Z.; Liang, M.; Chu, F. Recent advances in time–frequency analysis methods for machinery fault diagnosis: A review with application examples. Mech. Syst. Signal Process. 2013, 38, 165–205. [Google Scholar] [CrossRef]
  4. Berntsen, J.; Brandt, A.; Gryllias, K. Enhanced demodulation band selection based on Operational Modal Analysis (OMA) for bearing diagnostics. Mech. Syst. Signal Process. 2022, 181, 109300. [Google Scholar] [CrossRef]
  5. Singh, S.; Köpke, U.G.; Howard, C.Q. Analyses of contact forces and vibration response for a defective rolling element bearing using an explicit dynamics finite element model. J. Sound Vib. 2014, 333, 5356–5377. [Google Scholar] [CrossRef]
  6. Kedadouche, M.; Thomas, M.; Tahan, A. A comparative study between Empirical Wavelet Transforms and Empirical Mode Decomposition Methods: Application to bearing defect diagnosis. Mech. Syst. Signal Process. 2016, 81, 88–107. [Google Scholar] [CrossRef]
  7. Barbosh, M.; Singh, P.; Sadhu, A. Empirical mode decomposition and its variants: A review with applications in structural health monitoring. Smart Mater. Struct. 2020, 29, 093001. [Google Scholar] [CrossRef]
  8. Kim, Y.H.; Tan, A.C.C.; Kosse, V. Condition monitoring of low-speed bearings—A review. Aust. J. Mech. Eng. 2008, 6, 61–68. [Google Scholar] [CrossRef]
  9. Bechhoefer, E.; Qu, Y.; Zhu, J. Signal processing techniques to improve an acoustic emissions sensor. In Proceedings of the Annual Conference of the PHM Society, New Orleans, LA, USA, 14–17 October 2013; p. 5. [Google Scholar]
  10. Qu, Y.; He, D.; Yoon, J. Gearbox tooth cut fault diagnostics using acoustic emission and vibration sensors—A comparative study. Sensors 2014, 14, 1372–1393. [Google Scholar] [CrossRef]
  11. He, M.; He, D.; Bechhoefer, E. Using deep learning based approaches for bearing fault diagnosis with AE sensors. In Proceedings of the Annual Conference of the PHM Society, Denver, CO, USA, 3–6 October 2016; p. 8. [Google Scholar]
  12. Li, R.; He, D.; Bechhoefer, E. Investigation on fault detection for split torque gearbox using acoustic emission and vibration signals. In Proceedings of the Annual Conference of the PHM Society, San Diego, CA, USA, 27 September–1 October 2009; p. 1. [Google Scholar]
  13. Leaman, F.; Vicuña, C.M.; Clausen, E. A review of gear fault diagnosis of planetary gearboxes using acoustic emissions. Acoust. Aust. 2021, 49, 265–272. [Google Scholar] [CrossRef]
  14. Li, X.P.; Shin, W.Y.; Vartuli, J. Detection of water-ice transition using a lead zirconate titanate/brass transducer. J. Appl. Phys. 2002, 92, 106–111. [Google Scholar] [CrossRef]
  15. Pichika, S.N.; Yadav, R.; Rajasekharan, S.G. Optimal sensor placement for identifying multi-component failures in a wind turbine gearbox using integrated condition monitoring scheme. Appl. Acoust. 2022, 187, 108505. [Google Scholar] [CrossRef]
  16. Jena, D.P.; Panigrahi, S.N. Automatic gear and bearing fault localization using vibration and acoustic signals. Appl. Acoust. 2015, 98, 20–33. [Google Scholar] [CrossRef]
  17. Hemmati, F.; Orfali, W.; Gadala, M.S. Roller bearing acoustic signature extraction by wavelet packet transform, applications in fault detection and size estimation. Appl. Acoust. 2016, 104, 101–118. [Google Scholar] [CrossRef]
  18. Zhu, M.; Hu, W.; Kar, N.C. Acoustic noise-based uniform permanent-magnet demagnetization detection in SPMSM for high-performance PMSM drive. IEEE Trans. Transp. Electrif. 2017, 4, 303–313. [Google Scholar] [CrossRef]
  19. Glowacz, A. Acoustic-based fault diagnosis of commutator motor. Electronics 2018, 7, 299. [Google Scholar] [CrossRef] [Green Version]
  20. Paramo-Balsa, P.; Roldan-Fernandez, J.M.; Burgos-Payan, M. A Low-Cost Non-Intrusive Method for In-Field Motor Speed Measurement Based on a Smartphone. Sensors 2021, 21, 4317. [Google Scholar] [CrossRef]
  21. Wang, C.; Hu, F.; He, Q. De-noising of wayside acoustic signal from train bearings based on variable digital filtering. Appl. Acoust. 2014, 83, 127–140. [Google Scholar] [CrossRef]
  22. Lu, S.; Zheng, P.; Liu, Y. Sound-aided vibration weak signal enhancement for bearing fault detection by using adaptive stochastic resonance. J. Sound Vib. 2019, 449, 18–29. [Google Scholar] [CrossRef]
  23. Lu, W.; Jiang, W.; Wu, H. A fault diagnosis scheme of rolling element bearing based on near-field acoustic holography and gray level co-occurrence matrix. J. Sound Vib. 2012, 331, 3663–3674. [Google Scholar] [CrossRef]
  24. Hou, J.; Jiang, W.; Lu, W. Application of a near-field acoustic holography-based diagnosis technique in gearbox fault diagnosis. J. Vib. Control 2013, 19, 3–13. [Google Scholar] [CrossRef]
  25. Chen, Z.M.; Zhu, H.C.; Peng, M. Identification and localization of the sources of cyclostationary sound fields. Appl. Acoust. 2015, 87, 64–71. [Google Scholar] [CrossRef]
  26. Nejade, A. On holographic reconstruction without reference: Application to machinery noise and comparison with multi-reference holography approach. Appl. Acoust. 2017, 116, 348–356. [Google Scholar] [CrossRef]
  27. Ma, Y.; Chin, C.S. Near-field vibro-acoustic transfer function prediction of small close fit enclosure with multiple rotating components. Appl. Acoust. 2019, 149, 129–142. [Google Scholar] [CrossRef]
  28. Olivieri, M.; Pezzoli, M.; Antonacci, F. A physics-informed neural network approach for nearfield acoustic holography. Sensors 2021, 21, 7834. [Google Scholar] [CrossRef]
  29. Wang, J.; Zhang, Z.; Li, Z. Research on joint training strategy for 3D convolutional neural network based near-field acoustical holography with optimized hyperparameters. Measurement 2022, 202, 111790. [Google Scholar] [CrossRef]
  30. Cabada, E.C.; Leclere, Q.; Antoni, J. Fault detection in rotating machines with beamforming: Spatial visualization of diagnosis features. Mech. Syst. Signal Process. 2017, 97, 33–43. [Google Scholar] [CrossRef]
  31. Verellen, T.; Verbelen, F.; Stockman, K. Beamforming applied to ultrasound analysis in detection of bearing defects. Sensors 2021, 21, 6803. [Google Scholar] [CrossRef]
  32. He, Q.; Wang, J.; Hu, F. Wayside acoustic diagnosis of defective train bearings based on signal resampling and information enhancement. J. Sound Vib. 2013, 332, 5635–5649. [Google Scholar] [CrossRef]
  33. Zhang, D.; Entezami, M.; Stewart, E. Adaptive fault feature extraction from wayside acoustic signals from train bearings. J. Sound Vib. 2018, 425, 221–238. [Google Scholar] [CrossRef]
  34. Huang, H.; Liu, F.; Geng, L. Fault diagnosis accuracy improvement using wayside rectangular microphone array for health monitoring of railway-vehicle wheel bearing. IEEE Access 2019, 7, 87410–87424. [Google Scholar] [CrossRef]
  35. He, S.; Hu, D.; Yu, G. Trackside acoustic detection of axle bearing fault using wavelet domain moving beamforming method. Appl. Acoust. 2022, 195, 108851. [Google Scholar] [CrossRef]
  36. Bush, D.; Xiang, N. Broadband implementation of coprime linear microphone arrays for direction of arrival estimation. J. Acoust. Soc. Am. 2015, 138, 447–456. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Bush, D.; Xiang, N. A model-based Bayesian framework for sound source enumeration and direction of arrival estimation using a coprime microphone array. J. Acoust. Soc. Am. 2018, 143, 3934–3945. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Dong, F.; Jiang, Y.; Yan, Y. Direction-of-arrival tracking using a co-prime microphone array: A particle filter perspective. Appl. Acoust. 2020, 170, 107499. [Google Scholar] [CrossRef]
  39. Nannuru, S.; Gerstoft, P.; Ping, G. Sparse planar arrays for azimuth and elevation using experimental data. J. Acoust. Soc. Am. 2021, 149, 167–178. [Google Scholar] [CrossRef] [PubMed]
  40. Roy, R.; Paulraj, A.; Kailath, T. ESPRIT—A subspace rotation approach to estimation of parameters of cisoids in noise. IEEE Trans. Acoust. Speech Signal Process. 1986, 34, 1340–1342. [Google Scholar] [CrossRef]
  41. Roy, R.; Kailath, T. ESPRIT-estimation of signal parameters via rotational invariance techniques. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 984–995. [Google Scholar] [CrossRef] [Green Version]
  42. Xu, Z.; Wu, S.; Yu, Z. A robust direction of arrival estimation method for uniform circular array. Sensors 2019, 19, 4427. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Co−prime circular array (CPCA) configuration.
Figure 1. Co−prime circular array (CPCA) configuration.
Sensors 23 03050 g001
Figure 2. CPCA signal model.
Figure 2. CPCA signal model.
Sensors 23 03050 g002
Figure 3. Test rig.
Figure 3. Test rig.
Sensors 23 03050 g003
Figure 4. Bearing fault conditions: (a) outer race fault; (b) inner race fault; (c) roller fault.
Figure 4. Bearing fault conditions: (a) outer race fault; (b) inner race fault; (c) roller fault.
Sensors 23 03050 g004
Figure 5. Experimental setup.
Figure 5. Experimental setup.
Sensors 23 03050 g005
Figure 6. Microphone arrangement.
Figure 6. Microphone arrangement.
Sensors 23 03050 g006
Figure 7. Signal processing and analysis procedure.
Figure 7. Signal processing and analysis procedure.
Sensors 23 03050 g007
Figure 8. Envelopes under outer race fault condition at 290 rpm of: (a) non−processed mono signal; (b) UCMA−ESPRIT beamformed signal; (c) CPCMA−ESPRIT beamformed signal.
Figure 8. Envelopes under outer race fault condition at 290 rpm of: (a) non−processed mono signal; (b) UCMA−ESPRIT beamformed signal; (c) CPCMA−ESPRIT beamformed signal.
Sensors 23 03050 g008
Figure 9. Envelopes under inner race fault condition at 290 rpm of: (a) non−processed mono signal; (b) UCMA−ESPRIT beamformed signal; (c) CPCMA−ESPRIT beamformed signal.
Figure 9. Envelopes under inner race fault condition at 290 rpm of: (a) non−processed mono signal; (b) UCMA−ESPRIT beamformed signal; (c) CPCMA−ESPRIT beamformed signal.
Sensors 23 03050 g009
Figure 10. Envelopes under rolling element fault condition at 290 rpm of: (a) non−processed mono signal; (b) UCMA−ESPRIT beamformed signal; (c) CPCMA−ESPRIT beamformed signal.
Figure 10. Envelopes under rolling element fault condition at 290 rpm of: (a) non−processed mono signal; (b) UCMA−ESPRIT beamformed signal; (c) CPCMA−ESPRIT beamformed signal.
Sensors 23 03050 g010
Figure 11. Envelopes under outer race fault condition at 725 rpm of: (a) non−processed mono signal; (b) UCMA−ESPRIT beamformed signal; (c) CPCMA−ESPRIT beamformed signal.
Figure 11. Envelopes under outer race fault condition at 725 rpm of: (a) non−processed mono signal; (b) UCMA−ESPRIT beamformed signal; (c) CPCMA−ESPRIT beamformed signal.
Sensors 23 03050 g011
Figure 12. Envelopes under inner race fault condition at 725 rpm of: (a) non−processed mono signal; (b) UCMA−ESPRIT beamformed signal; (c) CPCMA−ESPRIT beamformed signal.
Figure 12. Envelopes under inner race fault condition at 725 rpm of: (a) non−processed mono signal; (b) UCMA−ESPRIT beamformed signal; (c) CPCMA−ESPRIT beamformed signal.
Sensors 23 03050 g012
Table 1. Bearing Configuration.
Table 1. Bearing Configuration.
Pitch DiameterBall DiameterBall NumberContact Angle
39 mm7.5 mm13
Table 2. Bearing Fault Feature Frequencies at 290 rpm.
Table 2. Bearing Fault Feature Frequencies at 290 rpm.
BPFO 1BPFI 2BSF 3FTF 4
25.375 Hz37.46 Hz12.102 Hz1.952 Hz
1 Ball pass frequency outer race; 2 ball pass frequency inner race; 3 ball spin frequency; 4 fundamental train frequency.
Table 3. DOA estimation by UCMA-ESPRIT for outer race fault signal at 290 rpm.
Table 3. DOA estimation by UCMA-ESPRIT for outer race fault signal at 290 rpm.
Elevation Angle (°)7.087.167.197.047.218.388.10
Azimuth Angle (°)155.29145.56157.75145.95−31.17145.83−32.81
Table 4. DOA estimation by CPCMA-ESPRIT for outer race fault signal at 290 rpm.
Table 4. DOA estimation by CPCMA-ESPRIT for outer race fault signal at 290 rpm.
Elevation Angle (°)7.987.567.707.897.657.957.54
Azimuth Angle (°)145.43150.45150.36151.62151.95152.10149.99
Table 5. Root-mean-square error (RMSE) in outer race defect locatied by both arrays at 290 rpm.
Table 5. Root-mean-square error (RMSE) in outer race defect locatied by both arrays at 290 rpm.
UCMACPCMA
Elevation AngleAzimuth AngleElevation AngleAzimuth Angle
0.60994.88200.17522.1399
Table 6. DOA estimation by UCMA-ESPRIT for inner race fault signal at 290 rpm.
Table 6. DOA estimation by UCMA-ESPRIT for inner race fault signal at 290 rpm.
Elevation Angle (°)7.696.335.897.956.79
Azimuth Angle (°)−126.25 1−81.67 1−37.1710.2559.89
Elevation Angle (°)4.314.734.064.587.24
Azimuth Angle (°)105.16 1150.68−161.60 1−119.72 1−78.96 1
1 Phase ambiguity exists.
Table 7. DOA estimation by CPCMA-ESPRIT for inner race fault signal at 290 rpm.
Table 7. DOA estimation by CPCMA-ESPRIT for inner race fault signal at 290 rpm.
Elevation Angle (°)7.026.315.174.905.80
Azimuth Angle (°)−129.24−82.83−36.458.6558.31
Elevation Angle (°)5.854.775.285.297.14
Azimuth Angle (°)106.19152.59−162.07−119.56−72.64
Table 8. RMSE in inner race defect located by both arrays at 290 rpm.
Table 8. RMSE in inner race defect located by both arrays at 290 rpm.
UCMACPCMA
Elevation AngleAngle IntervalElevation AngleAngle Interval
1.59472.88280.84891.8818
Table 9. DOA estimation of nearby outer race by UCMA-ESPRIT for rolling element fault signal at 290 rpm.
Table 9. DOA estimation of nearby outer race by UCMA-ESPRIT for rolling element fault signal at 290 rpm.
Elevation Angle (°)6.406.178.928.016.71
Azimuth Angle (°)43.30 1−102.87108.13 1−56.52163.85 1
Elevation Angle (°)10.497.987.048.9811.27
Azimuth Angle (°)17.4949.8979.61119.71145.14
1 Phase ambiguity exists.
Table 10. DOA estimation of nearby inner race by UCMA-ESPRIT for rolling element fault signal at 290 rpm.
Table 10. DOA estimation of nearby inner race by UCMA-ESPRIT for rolling element fault signal at 290 rpm.
Elevation Angle (°)4.704.224.427.456.35
Azimuth Angle (°)−108.02−56.66−5.86−134.11 1−71.14 1
Elevation Angle (°)4.626.346.395.566.33
Azimuth Angle (°)151.99−154.2593.62 1−32.7012.80
1 Phase ambiguity exists.
Table 11. DOA estimation of nearby outer race by CPCMA-ESPRIT for rolling element fault signal at 290 rpm.
Table 11. DOA estimation of nearby outer race by CPCMA-ESPRIT for rolling element fault signal at 290 rpm.
Elevation Angle (°)6.516.527.597.597.15
Azimuth Angle (°)−142.16−110.93−80.06−54.18−20.03
Elevation Angle (°)8.757.907.397.267.95
Azimuth Angle (°)8.3239.5470.03−83.00 1125.72
1 Phase ambiguity exists.
Table 12. DOA estimation of nearby inner race by CPCMA-ESPRIT for rolling element fault signal at 290 rpm.
Table 12. DOA estimation of nearby inner race by CPCMA-ESPRIT for rolling element fault signal at 290 rpm.
Elevation Angle (°)4.784.914.695.585.68
Azimuth Angle (°)−113.39−60.28−5.5150.04105.66
Elevation Angle (°)5.665.534.805.645.35
Azimuth Angle (°)156.8530.42 1−94.03−39.8714.83
1 Phase ambiguity exists.
Table 13. RMSE in ball pass location outer race (BPLO) by both arrays at 290 rpm.
Table 13. RMSE in ball pass location outer race (BPLO) by both arrays at 290 rpm.
UCMACPCMA
Elevation AngleAngle IntervalElevation AngleAngle Interval
1.68717.53690.71702.4965
Table 14. RMSE in ball pass location inner race (BPLI) by both arrays at 290 rpm.
Table 14. RMSE in ball pass location inner race (BPLI) by both arrays at 290 rpm.
UCMACPCMA
Elevation AngleAngle IntervalElevation AngleAngle Interval
1.09345.75570.39511.3909
Table 15. DOA estimation by UCMA-ESPRIT for outer race fault signal at 725 rpm.
Table 15. DOA estimation by UCMA-ESPRIT for outer race fault signal at 725 rpm.
Elevation Angle (°)7.657.607.627.647.627.647.57
Azimuth Angle (°)145.23145.50145.22145.32145.47145.87145.20
Table 16. DOA estimation by CPCMA-ESPRIT for outer race fault signal at 725 rpm.
Table 16. DOA estimation by CPCMA-ESPRIT for outer race fault signal at 725 rpm.
Elevation Angle (°)7.697.767.697.847.707.657.80
Azimuth Angle (°)145.96145.87153.88146.33145.83146.36146.00
Table 17. Root-mean-square error (RMSE) in outer race defect located by both arrays at 725 rpm.
Table 17. Root-mean-square error (RMSE) in outer race defect located by both arrays at 725 rpm.
UCMACPCMA
Elevation AngleAzimuth AngleElevation AngleAzimuth Angle
0.17194.60390.08553.9377
Table 18. DOA estimation by UCMA-ESPRIT for inner race fault signal at 725 rpm.
Table 18. DOA estimation by UCMA-ESPRIT for inner race fault signal at 725 rpm.
Elevation Angle (°)4.554.525.726.054.81
Azimuth Angle (°)−79.95137.79 1−173.15 154.3098.68
Elevation Angle (°)6.176.155.705.506.06
Azimuth Angle (°)146.10−163.41−115.82−68.55−19.52
1 Phase ambiguity exists.
Table 19. DOA estimation by CPCMA-ESPRIT for inner race fault signal at 725 rpm.
Table 19. DOA estimation by CPCMA-ESPRIT for inner race fault signal at 725 rpm.
Elevation Angle (°)4.774.794.875.124.95
Azimuth Angle (°)−89.59−42.012.7249.9293.11
Elevation Angle (°)5.644.955.575.285.55
Azimuth Angle (°)−38.45 19.74 1−121.43−75.22−27.71
1 Phase ambiguity exists.
Table 20. RMSE in inner race defect located by both arrays at 725 rpm.
Table 20. RMSE in inner race defect located by both arrays at 725 rpm.
UCMACPCMA
Elevation AngleAngle IntervalElevation AngleAngle Interval
0.66733.55910.35011.8104
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, C.; Chen, C.; Gu, X. Acoustic-Based Rolling Bearing Fault Diagnosis Using a Co-Prime Circular Microphone Array. Sensors 2023, 23, 3050. https://doi.org/10.3390/s23063050

AMA Style

Li C, Chen C, Gu X. Acoustic-Based Rolling Bearing Fault Diagnosis Using a Co-Prime Circular Microphone Array. Sensors. 2023; 23(6):3050. https://doi.org/10.3390/s23063050

Chicago/Turabian Style

Li, Chi, Changzheng Chen, and Xiaojiao Gu. 2023. "Acoustic-Based Rolling Bearing Fault Diagnosis Using a Co-Prime Circular Microphone Array" Sensors 23, no. 6: 3050. https://doi.org/10.3390/s23063050

APA Style

Li, C., Chen, C., & Gu, X. (2023). Acoustic-Based Rolling Bearing Fault Diagnosis Using a Co-Prime Circular Microphone Array. Sensors, 23(6), 3050. https://doi.org/10.3390/s23063050

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop