喜歡這套資料就充值下載吧。。。資源目錄里展示的都可在線預覽哦。。。下載后都有，，請放心下載，，文件全都包含在內，，【有疑問咨詢QQ:1064457796 或 414951605】 =================================喜歡這套資料就充值下載吧。。。資源目錄里展示的都可在線預覽哦。。。下載后都有，，請放心下載，，文件全都包含在內，，【有疑問咨詢QQ:1064457796 或 414951605】 =================================

An order tracking technique for the gear fault diagnosis using local mean decomposition method Junsheng Cheng, Kang Zhang, Yu Yang State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body, Hunan University, Changsha, 410082, PR China College of Mechanical and Vehicle Engineering, Hunan University, Changsha, 410082, PR China article info abstract Article history: Received 17 November 2010 Received in revised form 13 December 2011 Accepted 30 April 2012 Available online 28 May 2012 Local mean decomposition (LMD) is a new self-adaptive timefrequency analysis method, which is particularly suitable for the processing of multi-component amplitude-modulated and frequency-modulated (AMFM) signals. By using LMD, any complicated signal can be decomposed into a number of product functions (PFs), each of which is the product of an envelope signal and a purely frequency modulated signal from which physically meaningful instantaneous frequencies can be obtained. Theoretically, each PF is exactly a mono-component AMFMsignal. Therefore, the procedure of LMD can beregardedastheprocess of demodulation. While fault occurs in gear, the vibration signals would exactly present AMFM characteristics. Therefore, targeting the modulation feature of gear fault vibration signal in run-ups and run- downs and the fact that fault characteristics found in gear vibration signal could often be related to revolution of the shaft in the transient process, a gear fault diagnosis method in which order tracking technique and local mean decomposition is put forward. The analysis results from the practical gearbox vibration signal demonstrate that the proposed algorithm is effective in gear fault feature extraction. 2012 Elsevier Ltd. All rights reserved. Keywords: Order tracking technique Local mean decomposition Demodulation Gear Fault diagnosis 1. Introduction Gears are the important and frequently encountered components in the rotating machines that find widespread industrial applications. Therefore, the corresponding gear fault diagnosis has been the subject of extensive research. The key step of gear fault diagnosis is the extraction of fault feature. On the one hand, the conventional gear fault diagnosis methods focus on examining the frequency spectrum analysis of vibration signal at a fixed rotation speed. Unfortunately, the information obtained thus is only partial because some faults maybe do not respond significantly at the fixed operation speed. Since faults commonly found in gear could often be related to revolution of the shaft, more comprehensive information may be acquired by measuring the gear vibration signal in the process of run-up and run-down 1. In addition, vibration signals derived fromgearinthetransientprocessthatarespeed-dependentalwaysdisplaynon-stationaryfeature.Iffrequencyspectrumanalysis is directly applied to the non-stationary vibration signal, frequency mixing would occur inevitably, which will bring undesirable effect to the fault feature extraction. In past research, order-tracking technique, which normally exploits a vibration signal supplemented with information of shaft speed of rotating machinery, has become one of the significant approaches for fault diagnosis in rotating machinery 2,3. Essentially, order-tracking technique can transform a non-stationary signal in time domain into stationary one in angular domain, which can highlight the vibration information related to rotation speed and restrain the unrelated information. Therefore, order tracking is a desirable method to extract gear fault feature in the process of run-up and run-down. Mechanism and Machine Theory 55 (2012) 6776 Corresponding author at: State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body, Hunan University, Changsha, 410082, PR China. Tel.: +86 731 88664008; fax: +86 731 88711911. E-mail address: (J. Cheng). 0094-114X/$ see front matter 2012 Elsevier Ltd. All rights reserved. doi:10.1016/j.mechmachtheory.2012.04.008 Contents lists available at SciVerse ScienceDirect Mechanism and Machine Theory journal homepage: On the other hand, while faults occur in gears, the vibration signal picked up in run-up and run-down process always present the characteristics of amplitude-modulated and frequency-modulated (AMFM). In order to extract the modulation feature of gear fault vibration signals, demodulation analysis is one of the most popular methods 4,5. However, conventional demodulation approaches such as Hilbert transform demodulation and traditional envelope analysis have their own limitations 6. These drawbacks include two aspects: (1) in practice most gear fault vibration signals are all multi-component AMFM signals. For these signals, in conventional demodulation approaches, they are usually decomposed into single component AMFM signals by band-pass filter and then demodulated to extract frequencies and amplitudes information. However, both the number of the carrier frequency components and the magnitude of the carrier frequency are hard to be determined in practice, so the selection of central frequency of band-pass filter carries great subjectivity that would bring demodulation error and make it ineffective to extract the characteristic of machinery fault vibration signal; (2) owing to the inevitable window effect of Hilbert transform, when Hilbert transform is used to extract the modulate information, the demodulation results present non- instantaneous response characteristic, that is, at the two ends of the modulated signal which has been demodulated as well as the middle part with break would produce modulation again, which makes the amplitude get fluctuation in an exponential attenuation way, and then the demodulation error would increase 7. In order to overcome the first drawback, an appropriate decomposition method should be looked for to separate multi-component signal into a number of single component AMFM signals before the envelope analysis. Since EMD (Empirical mode decomposition) could adaptively decompose a complicated multi-component signal into a sum of intrinsic mode functions (IMFs) whose instantaneous frequencies have physical significance 8,9, order tracking method based on EMD has been widely used in the gear fault diagnosis 1013. However, there still exist many deficiencies in EMD such as the end effects 14 and modes mixing 15 that are still underway. In addition, after the original signal is decomposed by EMD, the drawback produced by Hilbert transform (above mentioned) is inevitable when IMF is performed envelope analysis by Hilbert transform. Moreover, sometimes the unexplainable negative instantaneous frequency would appear when calculating instantaneous frequency by performing Hilbert transform to each IMF 16. Local mean decomposition (LMD) is a novel demodulation analysis method,which is particularlysuitable for the processing of multi-component amplitude-modulated and frequency-modulated (AMFM) signals 16. By using LMD, any complicated signal can be decomposed into a number of product functions (PFs), each of which is the product of an envelope signal (obtained directly by the decomposition) from which instantaneous amplitude of the PF can be obtained and a purely frequency modulated signal from which a well-defined instantaneous frequency could be calculated. In essence, each PF is exactly a mono-component AMFM signal. Therefore, the procedure of LMD could be, in fact, regarded as the process of demodulation. Modulation informationcan be extracted by performing spectrum analysis to the instantaneous amplitude (envelope signal, obtained directly by the decomposition) of each PF component rather than by performing Hilbert transform to the PF components. Hence, when LMD and EMD are applied to the demodulation analysis respectively, compared with EMD, the prominent advantage of LMD is to avoid the Hilberttransform. In addition, the LMD iterationprocess which uses smoothed local means and local magnitudes avoids the cubic spline approach used in EMD, which maybe bring the envelope errors and influence on the precision of the instantaneous frequency and amplitude.Moreover, compared with EMD the end effect is not obvious in LMD approachbecause of faster algorithm speed and less iterative times 17. Based upon the above analysis, order-tracking analysis and the recent development of demodulation techniques, LMD, are combined and applied to the gear fault diagnosis of various shaft speeds process. Firstly, order tracking technique is used to transformthe gear vibration signals from time domainto angular domain. Secondly,decompose the re-sampling signal of angular domain by LMD, thus s series PF components and corresponding instantaneous amplitudes and instantaneous frequencies can be obtained.Finally,spectrumanalysisiscarriedoutto theinstantaneousamplitudesof the PFcomponentcontainingdominant fault information. The analysis results from the experimental vibration signal show that the proposed method can extract fault feature of the gear effectively and classify working condition accurately. This paper is organized as follows. A theory of the LMD approach is given in Section 2. In Section 3 a gear fault diagnosis approach in which order tracking technique and LMD are combined is put forward and the practice applications of proposed method are demonstrated. In addition, the comparison between LMD-based and EMD-based method is also given in Section 3. Finally, we offer the conclusion in Section 4. 2. LMD analysis method As mentioned above, the nature of LMD is to demodulate AMFM signals. By using LMD a complicated signal can be decomposed into a set of product functions, each of which is the product of an envelope signal and a purely frequency modulated signal.Furthermore,the completedtimefrequencydistributionoftheoriginalsignalcanbeobtained. Foranysignal x(t),it canbe decomposed as follows 16: (1) Determine all local extrema n i of the original signal x(t), and then the mean value m i of two successive extrema n i and n i+1 can be calculated by m i n i n i1 2 1 All mean value m i of two successive extreme are connected by straight lines, and then local mean function m 11 (t) can be formed by using moving averaging to smooth the local means m i . 68 J. Cheng et al. / Mechanism and Machine Theory 55 (2012) 6776 (2) A corresponding envelope estimate a i is given by a i n i n i1 C12 C12 C12 C12 2 2 Similarly, the envelope estimate a i is smoothed in the same way and the corresponding envelope function a 11 (t)is formed. (3) The local mean function m 11 (t) is subtracted from the original signal x(t) and the resulting signal h 11 (t) is given by h 11 txtm 11 t 3 (4) h 11 (t) can be amplitude demodulated by dividing it by envelope function a 11 (t) s 11 th 11 t=a 11 t 4 Ideally, s 11 (t) is a purely frequency modulated signal, namely, the envelope function a 12 (t)ofs 11 (t)shouldsatisfy a 12 (t)=1.Ifa 12 (t)1, then s 11 (t) is regarded as the original signal and the above procedure needs to be repeated until a purely frequency modulated signal s 1n (t) that meets 1s 1n (t)1 is derived. In other words, envelope function a 1(n+1) (t) of the resulting s 1n (t)shouldsatisfya 1(n+1) (t)=1. Therefore h 11 txtm 11 t h 12 s 11 tm 12 t h 1n ts 1 n1 tm 1n t 8 : 5 in which, s 11 th 11 t=a 11 t s 12 th 12 t=a 12 t s 1n th 1n t=a 1n t 8 : 6 where the objective is that lim n a 1n t1 7 In practice, a variation can be determined in advance. If 1a 1(n+1) (t)1+ and 1s 1n (t)1, then iterative process would be stopped. (5) Envelope signal a 1 (t), namely, instantaneous amplitude function, can be derived by multiplying together the successive envelope estimate functions that are acquired during the iterative process described above. a 1 ta 11 ta 12 ta 1n t n q1 a 1q t 8 where q is the times of the iterative process. (6) Multiplying envelope signal a 1 (t) by the purely frequency modulated signal s 1n (t) the first product function PF 1 of the original signal can be obtained. PF 1 ta 1 ts 1n t 9 PF 1 contains the highest frequency oscillations of the original signal. Meantime, it is a mono-component AMFM signal, whose instantaneous amplitude is exactly the envelope signal a 1 (t) and instantaneous frequency is defined from the purely frequency modulated signal s 1n (t)as f 1 t 1 2 d arccos s 1n tC138 dt 10 (7) Subtract the first PF component PF 1 (t) from the original signal x(t) and we have a new signal u 1 (t), which becomes the new originalsignalandthewholeoftheaboveprocedureisrepeated,i.e.uptoktimes,untilu k becomesmonotonicfunction u 1 txtPF 1 t u 2 tu 1 tPF 2 t u k tu k1 tPF k t 8 : 11 69J. Cheng et al. / Mechanism and Machine Theory 55 (2012) 6776 Thus, the original signal x(t) was decomposed into k-product and a monotonic function u k xt X k p1 PF p t u k t 12 where p is the number of the product function. Furthermore, the corresponding complete timefrequency distribution could be obtained by assembling the instantaneous amplitude and instantaneous frequency of all PF components. 3. The gear fault diagnosis method based on order tracking technique and LMD 3.1. Order tracking analysis and the corresponding fault diagnosis method Order-trackingtechniquecouldtransformanon-stationarysignalintimedomainintoastationarysignalinangulardomainby applyingequi-angularre-samplingto vibrationsignalwithreference toshaftspeed.Furthermore, orderspectrumcan beobtained by using spectrum analysis to stationary signal in angular domain, thus the information related to rotation speed can be highlighted and the unrelated one could be restrained. Therefore, order-tracking is suitable for the vibration signal analysis of rotation machine. There are three popular techniques for producing synchronously sampled data: a traditional hardware solution, computed order tracking (COT) and order tracking based on estimation of instantaneous frequency 1820. The traditional hardware approach, which uses specialized hardware to dynamically adapt the sample rate, is only suitable for the case that rotating speed of shaft is relatively smooth, thus resulting to a high cost. The method of order tracking based on estimation of instantaneous frequency has no need for specialized hardware and thus cost is relatively low, however, it has failed to analyze multiple component signal. While in practice most gear fault vibration signals exactly present the characteristic of multi-component. Therefore, this technique has little practice significance. COT technique realized equi-angular re-sampling by software, therefore it not onlyrequires nospecialized hardware,but alsohave nolimitation for analysissignal thatmeans it is more flexibleand more accurate. Just for this reason, COT is introduced into the gear fault detection in this paper. The step of the gear fault diagnosis method based on order tracking technique and LMD can be listed as follows: (1) The vibration signals and a tachometer signal are asynchronously sampled, that is, they are sampled conventionally at equal time incrementst; (2) Calculate the time series t i corresponding to equi-angular increments by tachometer signals; (3) According to the time series t i , apply interpolation to the vibration signals, thus the synchronous sampling signal, namely, stationary signal in angular domain, can be obtained; (4) Use LMD to decompose the equi-angular re-sampling signal, thus s series PF componentsand corresponding instantaneous amplitudes and instantaneous frequencies can be acquired; (5) Apply spectrum analysis to the instantaneous amplitude of each PF component, and then we have the order spectrum. 3.2. Application Since the gear fault vibration signal in run-up and run-down process are always multiple component AMFM signals and fault feature frequency would vary with rotation speed, the fault diagnosis method in which order tracking technique and LMD are combined would be suitable for gear fault detection. To verify the effectiveness of the proposed method, the fault diagnosis method based on order tracking technique and LMD was applied to the experimental gear vibration signals analysis. An experiment has been carried out on the rotating machinery test rig that is used for modeling different gear faults 21. Here we consider three working conditions that are gear with normal condition, with cracked tooth and with broken tooth. Standard gears with teeth number z=55 and z=75 are used on input and output shafts respectively, in which the crack fault is introduced into the gear on the input shaft by cutting slot with laser in the root of tooth, and the width of the slot is 0.15 mm, as well as its depth is 0.3 mm. Therefore, the mesh order is x m =55 and the fault feature order is x c =1.Figs. 1 and 2 give the rotation speed signal r(t) picked up by a tachometer and vibration acceleration signal s(t) of the gear with crack fault collected by a piezoelectric acceleration sensor respectively, in which the sample frequency is 8192 Hz and total sample time is 20 s, and from which we know the speed of input shaft increased gradually from 150 rpm to 1410 rpm, then decreased to 820 rpm. Meantime, the amplitude of vibration acceleration signal accordingly changed, from which a section of signal s 1 (t)of5s7 s in the run-up progress is intercepted for further analysis. Fig. 3 gives the spectrum of s 1 (t)by applying spectrum analysis directly to vibration signal. For the rotation speed changes with time, the frequency mixing arises. Therefore, it is impossible to find meshing frequency and fault feature frequency in Fig. 3. As a result, actual gear working condition cannot be identified. Replace direct spectrum analysis by the order tracking method. Firstly, assume sample point per rotation is 400, namely,the maximum analysis orderis 200. Secondly,angular domainsignal j 1 () shown in Fig. 4 can be obtained by performing order re-sampling to s 1 (t), in which horizontal ordinate has changed from time to radian. Thirdly, the corresponding order spectrum of j 1 () can be calculated that is illustrated in Fig. 5, from which we can find obvious spectral peak 70 J. Cheng et al. / Mechanism and Machine Theory 55 (2012) 6776 values at order O=55 and O=110 corresponding to gear meshing order and the double. Thus it means that frequency aliasing phenomenon has been eliminated to a large degree. However, j 1 () is still a multiple component MAMF signal. Therefore, side frequency band reflecting fault feature frequency is indistinct. To extract fault characteristic effectively, apply LMD to j 1 (), thus sevenPF componentsand aresiduecan beobtainedshowninFig.6, whichmeansLMDis a demodulation progress. Therefore,it is possible to extract gear fault feature by utilizing spectrum analysis to the instantaneous amplitude of PF component containing dominant fault information. By analysis, we know that the main failure information is included in the first PF component. Therefore, Figs. 7 and 8 give instantaneous amplitude a 1 () of the first PF component PF 1 () and the corresponding order spectrum of a 1 (), from which it is clear that there are distinct spectral peak value at the 1st order (O=1) corresponding to gear fault feature order x c , which accords with the actual working condition of the gear. Figs. 9 and 10 show the rotation speed signal n(t) and the time domain waveform of vibration acceleration signal s(t) of the gear with broken tooth respectively, in which the sample rate is 8192 Hz and total sample time is 20 s. The broken tooth fault is introduced into the gear on the input shaft by cutting slot with laser in the root of tooth. Firstly, a section of signal s 1 (t)of5s7s in the run-up progress is intercepted for further analysis; secondly, assume sample point per rotation is 400; thirdly, angular domain signal j 1 () shown in Fig. 11 can be obtained by performing order re-sampling to s 1 (t); fourthly, apply LMD to j 1 (); finally, the corresponding order spectrum shown in Fig. 12 of instantaneous amplitude of t