Bearing fault signal analysis is an important means of bearing fault diagnosis. To effectively eliminate noise in a fault signal, an adaptive multiscale combined morphological filter is proposed based on the theory of mathematical morphology. Both simulation and experimental results show that the adaptive multiscale combined morphological filter can remove noise more thoroughly and retain details of the fault signal better than the dual-tree complex wavelet filter, traditional morphological filter, adaptive singular value decomposition method (ASVD), and improved switching Kalman filter (ISKF). The adaptive multiscale combined morphological filter considers both positive and negative impulses in the signal; therefore, it has strong adaptability to complex noise in the environment, making it an effective new method for bearing fault diagnosis.
National Natural Science Foundation of China513054541. Introduction
Bearings are important components in rotating mechanical equipment, and their operating condition affects the overall working state of mechanical equipment. A bearing is prone to defects after prolonged use, and these defects can be useful during bearing fault diagnosis [1–4]. When a bearing fails, its vibration signal is composed of a modulation signal, periodic pulse signal, and noise signal. The periodic pulse signal is caused by the bearing fault, and the bearing fault type can be determined by analyzing the characteristic frequency of the periodic pulse signal. Due to the complex working environment of the bearing, the vibration signals of the bearing collected by relevant equipment show obvious nonstationary and nonlinear characteristics [5–9].
The mathematical morphology filter, developed on the basis of mathematical morphology transformation, is an effective method of analyzing nonlinear signals. The mathematical morphology filter selects appropriate structuring elements to decompose the signal into several components, according to the geometric characteristics of the signal. Even if there is strong noise or distortion in the signal, the basic morphological characteristics of the signal can be preserved after the signal is filtered. Compared with other signal analysis methods, the mathematical morphology filter, which has been widely used in pattern recognition, image processing, computer vision, power signals, ECG and EEG signal processing, mechanical equipment fault diagnosis, and other fields, has the advantages of easy implementation, fast algorithms, and minimal calculations [10–14]. Because the collected vibration signal of the bearing inevitably contains noise, the frequencies of the noise and fault signal will overlap in the frequency spectrum. A traditional morphological filter uses a single-scale operation, so its noise-filtering effects are general; in addition, it may filter out useful fault signals, which may prevent a true reflection of the bearing fault signals [15–18].
An adaptive multiscale combined morphological filter is developed based on the concept of multiscale overall filtering to overcome the shortcomings of the traditional morphological filter. Both simulation and experimental results show that the adaptive multiscale combined morphological filter has a better denoising effect and can retain useful signals better than the traditional filter.
2. Summary of Mathematical Morphology Filtering Theory2.1. Basic Morphology Operators
Dilation and erosion are the two basic operators of mathematical morphology. The collected vibration signal of the bearing is a one-dimensional dispersed signal, so this paper only introduces a one-dimensional dispersed gray-value morphological transformation.
fn is defined as a disperse one-dimensional input signal, and its domain is F=0,1,2⋯,N−1. gm is the structuring element, and its domain is G=0,1,2⋯,M−1; in addition, M≪N.
The dilation of fn on gm is
(1)f⊕gn=maxfn−m+gmm∈G.
The erosion of fn on gm is
(2)fΘgn=minm∈Gfn+m−gm.
The dilation operation is equivalent to the maximum value filtering of the disperse function in the sliding filter window (structuring element), which widens the peak value as well as increases the valley value of the signal. The erosion operation is equivalent to the minimum value filtering of the disperse function in the sliding filter window (structuring element), which widens the valley value as well as increases the peak value of the signal. Dilation and erosion, as two basic operations, can be used to form opening operations and closing operations, and their expressions are as follows:
The opening operation
(3)f∘gn=fΘg⊕gn.
The closing operation
(4)f·gn=f⊕gΘgn.
From Equation (3) and Equation (4), it can be seen that the opening operation means the signal is eroded first and then dilated by the structuring element, and the closing operation means the signal is dilated first and then eroded by the structuring element. The opening operation can eliminate the isolated point and weaken the positive pulse interference in the signal to smoothen the signal waveform. The closing operation can repair the holes and gaps, fill the grooves to make them flat, and suppress the negative pulse interference in the signal waveform.
Selection of the structuring element used in morphological transformation is also very important because different structuring elements have different effects on signal processing. Several common types of structuring elements are flat, circular, cosine, triangular, and curved. The flat structuring element has the following advantages: the amplitude is zero, the calculation is simple and fast, the signal processing efficiency is high, and the processing effect can be guaranteed. Therefore, in this paper, a flat structuring element is selected for the morphological filtering of the bearing fault signal.
2.2. Traditional Morphological Filter
Because the dilation, erosion, opening, and closing operations can only filter one of the positive and negative impulse noise in the signal, in many cases, the opening and closing operations are combined to construct the morphological closing-opening filter and morphological opening-closing filter, which are as shown in Equation (5) and Equation (6), respectively.
(5)FCOfn=f·g∘gn.(6)FOCfn=f∘g·gn.
Morphological closing-opening and opening-closing filters combine two kinds of operations (opening and closing), which can weaken both positive and negative impulse noise in the signal, but they cannot avoid the statistical bias of the filtering results caused by the dilatability of the closing operation and the erodibility of the opening operation. To achieve a better filtering effect, the combined morphological filter is usually constructed by combining morphological closing-opening and opening-closing filters, as shown in Equation (7).
(7)FCO−OCfn=12FCOfn+FOCfn.
2.3. Multiscale Morphological Filter
Originally, multiscale morphology was established to enrich the practicability of morphology in the field of image processing and shorten the operation time when the size of the structuring element is large, which is based on decomposing morphological structuring elements. When applied to the research of a one-dimensional vibration signal, the structuring element of the multiscale morphological filter can be obtained by corresponding time dilations of the unit structuring element itself.
If the structuring element of scale i is ig and the unit structuring element is g, then ig can be obtained by i−1 times dilation of g itself, as shown in Equation (8).
(8)ig=g⊕g⊕⋯⊕g⏟i−1times.
The dilation and erosion operations of signal f by multiscale structuring elements are as shown in Equation (9) and Equation (10), respectively:
(9)f⊕i+1gn=f⊕ign⊕gn,(10)fΘi+1gn=fΘignΘgn.
Similarly, the multiscale morphological opening, closing, opening-closing, and closing-opening operations of signals are as shown in Equations (11) to (14):
(11)f∘gin=fΘign⊕ign,(12)f·gin=f⊕ignΘign,(13)FCOfn=f·ign∘ign,(14)FOCfn=f∘ign·ign.
The environment and working conditions of the bearing are often complex and changeable; consequently, noise in the collected signal of the bearing is also complex and diverse; thus, it is difficult to achieve an ideal filtering effect by using a single-scale structuring element. To enable the morphological filter to deal with noise in various situations, multiple groups of structuring elements with different scales are adopted to carry out morphological closing-opening and opening-closing combined filtering for a signal based on the concept of overall multiscale filtering. In addition, structuring elements with different scales are used successively in the same group. The final scale signal is obtained by adaptive weighting combination for each group of operation results. The principle is as follows:
(15)F1COfn=f·g11∘g12n,F1OCfn=f∘g11·g12n,F2COfn=f·g21∘g22n,F2OCfn=f∘g21·g22n,⋯⋯FiCOfn=f·gi1∘gi2n,FiOCfn=f∘gi1·gi2n.
In Equation (15), g11⊆g12,g21⊆g22,andgi1⊆gi2.
The final signal filtered by k groups of structuring elements is as shown in Equation (16):
(16)Ff=∑i=1kλi1FiCOfn+λi2FiOCfn,where λi1 represents the weight of the results of the morphological closing-opening operation on the signal by the ith group of structuring elements in the final signal. λi2 represents the weight of the results of the morphological opening-closing operation on the signal by the ith group of structuring elements in the final signal. λi1 and λi2 can be calculated by Equations (17) to (20):
(17)δi1=∑n=1NFiCOfn−fn,(18)δi2=∑n=1NFiOCfn−fn,(19)λi1=δi1∑i=1kδi1+δi2,(20)λi2=δi2∑i=1kδi1+δi2.
It can be seen that the morphological closing-opening (opening-closing) operation results from large-scale structuring elements have a large proportion of weight in the final signal, and the morphological closing-opening (opening-closing) operation results from small-scale structuring elements have a small proportion of weight, which considers denoising and protecting details under the premise of giving first priority to denoising.
The process of the proposed morphological filter is as follows:
Read the signal fn to be processed
Calculate the value of d (the number of sampling points in a pulse period) according to the sampling frequency and pulse frequency (or theoretical fault frequency)
Obtain the width 0.6d of the flat structuring element
According to Equation (15), use 2i (the value of i is optional) groups of different scale structuring elements to carry out morphological closing-opening and morphological opening-closing operations for f1n and obtain the final signal f1n according to Equations (16) to (20)
Implement a fast Fourier transform of f1n to obtain its spectrum F1n
Output f1n and F1n
3. Analysis of the Filtering Effect for the Simulation Signal
To test the effect of the proposed adaptive multiscale combined morphological filter, a simulation signal is created for filtering research, which is shown in Equation (21):
(21)xt=x1t+xnt.
In Equation (21), x1t is a periodic pulse attenuation signal with a frequency of 60 Hz and a weekly attenuation function of e−800·t·cos2π·800·t; xnt is a white Gaussian noise signal. The sampling frequency of the signal is 6000 Hz, and the sampling time is 1 s.
The amplitude of the flat structuring element is 0. According to [19], the width of the structuring element was selected as 0.6d (d is the number of sampling points in a pulse period). Since the sampling frequency is 6000 Hz and the pulse frequency is 60 Hz, it can be seen that d=100, so the width of the structuring element is 60.
The waveform and frequency spectrum of the simulated signal are shown in Figure 1. Clearly, there is considerable noise in the simulation signal; the signal-to-noise ratio is low; the pulse signal is completely submerged by noise; the frequency spectrum is relatively messy; neither the waveform nor the frequency spectrum can reflect the periodic characteristics of the pulse signal.
Simulated signal and its frequency spectrum. (a) The waveform of the simulated signal. (b) The frequency spectrum of the simulated signal.
The waveform and frequency spectrum of the simulated signal denoised by the dual-tree complex wavelet are shown in Figure 2. It can be seen obviously that noise in the signal is suppressed to a certain extent by the dual-tree complex wavelet, but considerable noise still remains. The characteristic frequency of the pulse signal can be reflected in the frequency spectrum, but there is still considerable interference. Moreover, the periodic characteristic of the signal is not obvious.
Waveform and frequency spectrum of the simulated signal denoised by the dual-tree complex wavelet. (a) Waveform. (b) Frequency spectrum.
The waveform and frequency spectrum of the simulated signal denoised by the traditional morphological filter with small-scale structuring elements are shown in Figure 3. It can be seen obviously that the denoising effect is better than that of the dual-tree complex wavelet, but a certain amount of noise remains. The characteristic frequency of the pulse signal can be reflected in the frequency spectrum, but the periodic characteristic of the signal is still not sufficiently obvious because of the interference of the left noise.
Waveform and frequency spectrum of simulated signal denoised by traditional morphological filter with small-scale structuring elements. (a) Waveform. (b) Frequency spectrum.
The waveform and frequency spectrum of the simulated signal denoised by the traditional morphological filter with large-scale structuring elements are shown in Figure 4. It can be seen obviously that some useful signals in the waveform are filtered, resulting in distortion of the signal, and the frequency spectrum only reflects the primary frequency of the pulse signal.
Waveform and frequency spectrum of simulated signal denoised by traditional morphological filter with large-scale structuring elements. (a) Waveform. (b) Frequency spectrum.
The waveform and frequency spectrum of the simulated signal denoised by the adaptive multiscale combined morphological filter are shown in Figure 5. It can be seen obviously that not only a great quantity of noise in the signal is effectively filtered but also the distortion of the effective signal is avoided. Only a miniscule quantity of clutter remains in the frequency spectrum, which reflects both the primary and double frequency of the pulse signal. The denoising effect is obviously better than that of the dual-tree complex wavelet and traditional morphological filter.
Waveform and frequency spectrum of the simulated signal denoised by the adaptive multiscale combined morphological filter. (a) Waveform. (b) Frequency spectrum.
4. Analysis of the Filtering Effect for the Bearing Fault Signal
The bearing fault signal is collected through an experiment involving a two-stage transmission gearbox whose structure is shown in Figure 6. In the first-stage transmission, gear A on the input shaft (I) and gear B on the intermediate shaft (II) mesh with each other. The tooth numbers of gear A and gear B are 23 and 47, respectively. In the second-stage transmission, gear C on the intermediate shaft (II) and gear D on the output shaft (III) mesh with each other. The tooth numbers of gear C and gear D are 19 and 78, respectively. During the experiment, the outer ring of bearing R3 on the intermediate shaft (II) is set to be defective, and the acceleration sensor is set at the bearing block of R3 to collect the fault signal. The speed of the input shaft (I) is 1520 r/min, the sampling frequency is 10 kHz, and the number of sampling points is 5000. It can be calculated that the characteristic frequency of the bearing R3 outer ring defect signal is 126.4 Hz.
Structure diagram of the tested gearbox.
The width of the structuring element is chosen as 0.6d (d is the number of sampling points in a pulse period). Since the sampling frequency is 10 kHz and the pulse frequency is 126.4 Hz, it can be seen that d=79, so the width of the structuring element is 47.
The waveform and frequency spectrum of bearing R3 with an outer ring defect obtained by the experiment are shown in Figure 7.
Waveform and frequency spectrum of bearing R3 with an outer ring defect. (a) Waveform. (b) Frequency spectrum.
The waveform and frequency spectrum of bearing R3 with an outer ring defect signal denoised by a dual-tree complex wavelet are shown in Figure 8. It can be clearly seen that the dual-tree complex wavelet can filter out most of the noise in the signal, but the filtering is insufficient. The frequency spectrum reflects the primary frequency (126 Hz) and double frequency (252 Hz) of the pulse signal; however, there is excessive clutter in the frequency spectrum, which easily produces interference.
Waveform and frequency spectrum of bearing R3 with an outer ring defect signal denoised by a dual-tree complex wavelet. (a) Waveform. (b) Frequency spectrum.
The waveform and frequency spectrum of bearing R3 with an outer ring defect signal denoised by a traditional morphological filter with small-scale structuring elements are shown in Figure 9. It can be seen obviously that the noise is less than that by the dual-tree complex wavelet; the frequency spectrum obviously highlights the primary frequency (126 Hz) and double frequency (252 Hz) of the pulse signal, and the clutter is further reduced as well.
Waveform and frequency spectrum of bearing R3 with an outer ring defect signal denoised by a traditional morphological filter with small-scale structuring elements. (a) Waveform. (b) Frequency spectrum.
The waveform and frequency spectrum of bearing R3 with an outer ring defect signal denoised by a traditional morphological filter with large-scale structuring elements are shown in Figure 10. It can be seen obviously that the signal is distorted, with a loss of partial effective information; the frequency spectrum only reflects the primary frequency (126 Hz) of the pulse signal.
Waveform and frequency spectrum of bearing R3 with an outer ring defect signal denoised by a traditional morphological filter with large-scale structuring elements. (a) Waveform. (b) Frequency spectrum.
The waveform and frequency spectrum of bearing R3 with an outer ring defect signal denoised by an adaptive multiscale combined morphological filter are shown in Figure 11. It can be seen obviously that the signal details are well preserved while the noise is effectively filtered out; the frequency spectrum clearly reflects the primary frequency (126 Hz) and double frequency (252 Hz) of the pulse signal; furthermore, the clutter is extremely small.
Waveform and frequency spectrum of bearing R3 with an outer ring defect signal denoised by an adaptive multiscale combined morphological filter. (a) Waveform. (b) Frequency spectrum.
In recent years, Qin et al. [20] proposed an adaptive singular value decomposition method (ASVD) for gearbox fault signal filtering. Cui et al. [21] proposed an improved switching Kalman filter (ISKF) to filter bearing fault signals. We use these two methods to denoise the outer ring defect signal of bearing R3, and the results are shown in Figures 12 and 13. It can be seen obviously that although these two methods can filter out most of the noise in the signal, yet they are both inferior to our proposed filter. The frequency spectrum reflects the primary frequency (126 Hz) and double frequency (252 Hz) of the pulse signal, but there is still much interference in the high-frequency part. By comparing with these two methods, the novelty of our proposed adaptive multiscale combined morphological filter is reflected.
Waveform and frequency spectrum of bearing R3 with an outer ring defect signal denoised by ASVD. (a) Waveform. (b) Frequency spectrum.
Waveform and frequency spectrum of bearing R3 with an outer ring defect signal denoised by ISKF. (a) Waveform. (b) Frequency spectrum.
5. Conclusions
To address the insufficient denoising effects of traditional filters, an adaptive multiscale combined morphological filter is proposed, which is composed of morphological closing-opening and opening-closing operations. The advantages of the proposed filter are verified by simulations and experiments, as specifically manifested in the following three aspects:
The morphological filter proposed in this paper uses multiple groups of structuring elements with different scales to carry out morphological closing-opening and opening-closing combined filtering for signals. In addition, structuring elements with different scales are used successively in the same group. The final scale signal is obtained by a self-adaptive weighting combination of each group of operation results. The morphological closing-opening (opening-closing) operation results by large-scale structuring elements have a large proportion of weight in the final signal, and the morphological closing-opening (opening-closing) operation results by small-scale structuring elements have a small proportion of weight, which can take into account the denoising and protecting details of the signal under the premise of giving first priority to denoising
The proposed filter can filter out noise more thoroughly and extract pulse signals more effectively than dual-tree complex wavelets, ASVD and ISKF. Compared with the traditional morphological filter with small-scale structuring elements, which has an incomplete denoising effect, the proposed filter can filter out noise more effectively; compared with the traditional morphological filter with large-scale structuring elements, which loses part of the useful signal, the proposed filter can preserve the useful signal much better
The proposed filter takes into account both positive and negative pulses in the signal; therefore, it has strong adaptability to complex and variable noise in actual environments. Thus, the filter provides a new effective method for rotating machinery fault diagnosis
Data Availability
The data used to support this research can be obtained from the corresponding author on request.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
Acknowledgments
The research has been supported by the National Natural Science Foundation of China (grant number 51305454).
BlochI.Duality vs. adjunction for fuzzy mathematical morphology and general form of fuzzy erosions and dilations2009160131858186710.1016/j.fss.2009.01.0062-s2.0-67349086131ChenH.JiangB.LuN.MaoZ.Deep PCA based real-time incipient fault detection and diagnosis methodology for electrical drive in high-speed trains20186764819483010.1109/TVT.2018.28185382-s2.0-85044321735GuX.ChenC.Rolling bearing fault signal extraction based on stochastic resonance-based denoising and VMD2017201712359587110.1155/2017/35958712-s2.0-85042224572HaralickR. M.SternbergS. R.ZhuangX.Image analysis using mathematical morphology19879453255010.1109/tpami.1987.47679412-s2.0-002338318721869411RodríguezP. H.AlonsoJ. B.FerrerM. A.TraviesoC. M.Application of the Teager-Kaiser energy operator in bearing fault diagnosis201352227828410.1016/j.isatra.2012.12.0062-s2.0-8487523561323352553JayaswalP.WadhwaniA. K.MulchandaniK. B.Machine fault signature analysis200820081058398210.1155/2008/5839822-s2.0-48349102089LiC.LiangM.ZhangY.HouS.Multi-scale autocorrelation via morphological wavelet slices for rolling element bearing fault diagnosis20123142844610.1016/j.ymssp.2012.03.0122-s2.0-84861591054LiC.LiangM.Continuous-scale mathematical morphology-based optimal scale band demodulation of impulsive feature for bearing defect diagnosis2012331265864587910.1016/j.jsv.2012.07.0452-s2.0-84866620291NamdariF.SalehiM.High-speed protection scheme based on initial current traveling wave for transmission lines employing mathematical morphology201732124625310.1109/TPWRD.2016.25713412-s2.0-85010447439PetersR. A.A new algorithm for image noise reduction using mathematical morphology19954555456810.1109/83.3824912-s2.0-002930961018290006RajA. S.MuraliN.Early classification of bearing faults using morphological operators and fuzzy inference201360256757410.1109/tie.2012.21882592-s2.0-84866598130ShihF. Y.ChengS.Adaptive mathematical morphology for edge linking20041671-492110.1016/j.ins.2003.07.0202-s2.0-9544254865ShuaiJ.ShenC.ZhuZ.Adaptive morphological feature extraction and support vector regressive classification for bearing fault diagnosis2017201710238418410.1155/2017/23841842-s2.0-85016583915WeiY.PeiL.Detection of power quality disturbances based on mathematical morphology (MM) filter20024101317GongX.QiaoW.Bearing fault diagnosis for direct-drive wind turbines via current-demodulated signals20136083419342810.1109/tie.2013.22388712-s2.0-84876205701YuH.LiH.-r.TianZ.-k.WangW.-g.Rolling bearing degradation state identification based on LPP optimized by GA2016201610928109810.1155/2016/92810982-s2.0-84983757037ZhangJ.YaoH.RizzoniG.Fault diagnosis for electric drive systems of electrified vehicles based on structural analysis20176621027103910.1109/TVT.2016.25566912-s2.0-85015025531ZhangL.XuJ.YangJ.YangD.WangD.Multiscale morphology analysis and its application to fault diagnosis200822359761010.1016/j.ymssp.2007.09.0102-s2.0-38349039525NikolaouN. G.AntoniadisI. A.Application of morphological operators as envelope extractors for impulsive- type periodic signals20031761147116210.1006/mssp.2002.15762-s2.0-0042381324QinY.ZhangQ. L.ZhaoY.Fault diagnosis method for planetary gearboxes based on adaptive SVD20183717122127CUIL.XinW.HuaqingW.YonggangX.JianyuZ.Feature extraction of bearing fault based on improved switching Kalman filter2019557445110.3901/JME.2019.07.0442-s2.0-85067652789