Abstract
In the present study, in order to detect the fault of the gearmeshs, two engaged gears based on research department of a major automotive company have been modeled. First off, by using the CATIA software the fault was induced to the output gear. Then, the faulty gearmesh and nonfaulty gearmesh is modeled to find the fault pattern to predict and estimate the failure of the gearmesh. The induced defect is according to the frequently practical fault that takes place to the teeth of gears. In order to record the acceleration signals to calculate the decomposition algorithm, mount the accelerometer on accessible place of the output shaft to recognize the pattern. Then, for more realistic simulation, noise is added to the output signal. At the first step by means of Butterworth low pass digital, the noise has to be removed from signals after that by using the Empirical Mode Decomposition (EMD), signals have decomposed into the Instinct Mode Function (IMF) and every IMF were tested by using the Instantaneous Frequency (IF) in way of Hillbert Transform (HT). For this purpose a code was developed in MATLAB software. Then, in order to detect the presence of the fault the frequency spectrum of IMF's are created and defect is detected in gearmesh frequency of the spectrum.
Keywords:
Signal Processing; Condition monitoring; Fault detection and EMD; Hillbert Transform (HT); Helical Gears
1 INTRODUCTION
In industrial sites and instruments there are always unavoidably faults in machinery. Therefore, because of fault existence in equipment, they do not work properly. In the other words, faults do not allow any part of the system work efficiently and will cause the irreparable damages to machineries. For this reason the condition monitoring process had invented to pave the way of maintaining industrial instruments at the healthy condition and keep them work properly. In maintenance part there are always three strategic and particular ways to keep instruments work properly with high efficiency. These strategic and particular ways include run to break, timebased maintenance and the more complex way conditionbased maintenance. First off, the simplest way is run to break which has nothing to discuss. In the second way, i.e., timebased maintenance, a company may appoint inspector to set the interval to check up the tools and may be failure occurred at the middle of interval time and sometimes it is crucial to stop the operation to run checkup. The third way is the real time update and complex way to run maintenance which checks the machinery every second without abruptly disruption in operating system. Conditionbased maintenance frequently has four basic ways which are consist of vibration analysis, lubricant analysis, performance analysis and thermography. In vibrational analysis conditionbased monitoring, every machinery or industrial tools have same pattern of vibration. These kinds of prevalent vibrations which are like mechanical signature make inspector capable to specifically recognize the health condition. For a century researchers make their researches in order to damp these types of vibrations and still they are trying to damp it, because this kind of vibrations may leave unwanted effect on parts, but nowadays because of information that these types of vibrations carry with themselves, they become more important and reliable in comparison to the past and also some researchers try to use it as a means for fault detection. So, they can be so useful in a way they have healthy information of the system. It should be noted that these vibrations are still unwanted and in some case dangerous to some systems.
In the signal processing method, a wavelet transform is able to decompose all kind of signals with respect to the mother function but because of the shape of the mother function only some sort of the signals could be decomposed Li Li et al (2007Li Li, Liangsheng Qu, Xianghui Liao, (2007) Haar wavelet for machine fault diagnosis, Mechanical Systems and Signal Processing 21: 17731786.). Some researchers make their effort to detect the faults by means of combined timefrequency entropy base on the Hilberttransform Dejie et al (2007Dejie Yu, Yu Yang, Junsheng Cheng, (2007) Application of TimeFrequency Entropy Method Based on HilbertHuang Transform to Gear Fault Diagnosis, Measurement 40: 823830.). Cheng et al (2008Junsheng Cheng, Yu Yang, Dejie Yu, (2008) Application of frequency family separation method based upon EMD and local Hilbert energy spectrum method to gear fault diagnosis, Mechanical Systems and Signal Processing 43: 712723.) used Empirical Mode Decomposition (EMD) and local Hilbert energy to detect the faults. Some authors like Chuan and Liang (2012Chuan Li, Ming Liang, (2012) TimeFrequency Signal Analysis for Gearbox Fault Diagnosis Using a Generalized Synchrosqueezing Transform, Mechanical Systems and Signal Processing 26: 205217.), Junsheng Cheng et al. (2010Junsheng Cheng, Kang Zhang, Yu Yang, (2012) An Order Tracking Technique for the Gear Fault Diagnosis Using Local Mean Decomposition Method, Mechanism and Machine Theory 55: 6776.) tried to detect the defects by using the timefrequency analysis. Chuan and Ming (2012Chuan Li, Ming Liang, (2012) TimeFrequency Signal Analysis for Gearbox Fault Diagnosis Using a Generalized Synchrosqueezing Transform, Mechanical Systems and Signal Processing 26: 205217.) introduced and used Generalized Synchrosqueezing Transform (GST) transform to detect the gearbox fault. That method of analyzing works properly at variable rotating speed of run up and down conditions, so it is suitable for amplitude and frequency modulated (AMFM) signals. Cheng et al (2010Junsheng Cheng, Yu Yang, Dejie Yu, (2010) The Envelope Order Spectrum Based on Generalized Demodulation TimeFrequency Analysis and Its Application to Gear Fault Diagnosis, Mechanical Systems and Signal Processing 24: 508521.) used both generalized demodulation timefrequency analysis and order spectrum to detect the fault in the gear matting systems. After decomposition process was done and then each envelope was resampled in an angle domain and finally the spectrum analysis applied to each resampled signal to make fault become detected. Cheng et al (2012Junsheng Cheng, Yi Yang, Yu Yang, (2012) A Rotating Machinery Fault Diagnosis Method Based on Local Mean Decomposition, Digital Signal Processing 22: 356366.) investigated envelope order spectrum based on Generalized Demodulation Method (GDM) to trace of faults in operating systems. They first decompose multicomponent signal in monocomponent signal and then they reveal fault with the aid of envelope spectrum timefrequency analysis. Cheng et al (2012Junsheng Cheng, Yi Yang, Yu Yang, (2012) A Rotating Machinery Fault Diagnosis Method Based on Local Mean Decomposition, Digital Signal Processing 22: 356366.) introduced the Local Mean Decomposition (LMD) to detect the faults. The LMD algorithm is a selfadaptive algorithm which is suitable for decompose multicomponent AMFM signal in monocomponent signal which is introduced as Productive Function (PF). Each PF's includes the meaningful frequency information of signal. Wang and Heyns (2011K.S. Wang, P.S. Heyns, (2011) An Empirical ReSampling Method on Intrinsic Mode Function to Deal with Speed Variation in Machine Fault Diagnostics, Applied Soft Computing 11: 50155027.) studied resampling method of an Intrinsic Mode Function (IMF). They showed that resampling process had an advantage to reconstruct zero mean oscillating and symmetric signals for frequency analysis of nonstationary signals. Li et al (2011Hui Li, Yuping Zhang, Haiqi Zheng, (2011) Application of Hermitian Wavelet to Crack Fault Detection in Gearbox, Mechanical Systems and Signal Processing 25: 13531363.) developed wavelet as a tool to decompose signals in different scales. They used Hermitian wavelet mother function which was derived from Gaussian function. In this article the wavelet tool applied to extract the intrinsic mode to detect the crack in variable speed. Ricci and Pennacchi (2011Roberto Ricci, Paolo Pennacchi, (2011) Diagnostics of Gear Faults Based on EMD and Automatic Selection of Intrinsic Mode Functions, Mechanical Systems and Signal Processing 25: 821838.) indicated the Hilbert Huang transform to diagnosis process. They defined the merit index to show which IMF contains fault information and by means of merit index they could choose which IMF was suitable, for frequency spectrum analysis to reveal faults. Yu et al (2007Dejie Yu, Yu Yang, Junsheng Cheng, (2007) Application of TimeFrequency Entropy Method Based on HilbertHuang Transform to Gear Fault Diagnosis, Measurement 40: 823830.) defined the mathematical expression as timefrequency entropy. In order to gear fault diagnosis, they used timefrequency entropy based on Hilbert Huang transform. Loutridis (2006S.J. Loutridis, (2006) Instantaneous energy density as a feature for gear fault detection, Mechanical Systems and Signal Processing 20: 12391253) introduced the energy density to detect the faults. In his work the feature of highly value of energy of fault signal was concerned and with aid of wavelet the fault was diagnosed. Zheng et al. (2002H. Zheng, Z. Li, X. Chen, (2002) Gear Fault Diagnosis Based on continuous Wavelet Transform, Mechanical Systems and Signal Processing 16(23): 447457) could found a new approach with aim of wavelet transform. In their work the useful concept of Time Average Wavelet Spectrum (TAWS) based on Morlet family of wavelet was used. Finally, in order to fault detection two ways which were Spectrum Comparison Method (SCM) and Feature Energy Method (FEM) were used. Chen et al. (2002D. Chen, W. J. Wang, (2002) Classification of wavelet map Patterns using MultiLayer neural networks For gear Fault Detection, Mechanical Systems and Signal Processing 16(4): 695704) by means of wavelet tool, decomposed the signal into the patterns and then by using the multilayer neural network pattern recognition instead of human directly control, the fault was recognized. Samanata (2004B. Samanata, (2004) Gear Fault Detection Using Artificial Neural Networks And Support Vector Machines With Genetic Algorithms, Mechanical Systems and Signal Processing 18: 625644) used Artificial Neural Network (ANN), Genetic Algorithm (GA) and feature of pattern recognition of artificial process to recognize the fault. Wang and Mcfadden (1996W. J. WANG, P. D. Mcfadden, (1996) Application Of Wavelets To Gearbox Vibration Signals For Fault Detection, Journal of Sound and Vibration 192(5): 816_828) were able to detect the fault of helicopter gearbox by using the wavelet method. Litak (2013Grzegorz Litak, Krzysztof Kecik, Rafal Rusinek, (2013) Cutting force response in milling of Inconel Analysis by wavelet and HilbertHuang Transforms, Latin American Journal of Solids and Structures 10: 133140.) tried to study the response of milling process with variable depth. They used the continue wavelet map, which was composed of convolution mean force response with modulated Morelet wavelet with Gaussian function and Hilbert Haung transform. They made their calculation with aim of MATLAB software. Rubio (2014Lourdes Rubio, Belen MunozAbella, Patricia Rubio, Laura Montero, (2014) Quasistatic numerical study of the breathing mechanism of an elliptical crack in an unbalanced rotating shaft, Latin American Journal of Solids and Structures 11: 23332350.) in order to keep a shaft safe from abruptly crack failure, they studied the effect of the crack breathing with existence of eccentric mass in different position. They used their result at the trend of crack propagation and evaluation of Stress Intensity Factor (SIF).
In some experimental, analytical and simulating works which reviewed above, the experimental specifications which were not really practical, was imposed theoretically ideal condition on experiment. Prior noted point left inaccuracy effect on whole result and conclusion in which no gearmating transmitting load, truly positioning sensor and practical rotation speed over really commercially nonstandard gearmesh were used. In other words, the aim of gearmating process is transmitting torque and rotation speed and all the industrial gears are not based on standard regulation and are based on the R&D section of its company. In order to reduce vibration and avoidance of fault occur in each gear which works at specific speed under specific load, they have specific tip relief and specific tooth profile. Each geometrical situation of practical gearbox may allow attaching sensor on specific bearing at specific location which may cause the modulation on vibrating signal and make it harder to decompose.
In this research the study was done in semipractical situation, because of the nature of some practical factors which are not the kind of simulation factor by the software. The selected gear was based on conventional gear in automobile industrial and also rotation speed was based on real epidemic linear vehicle speed with real torque transmitting based on engine generated torque. In order to have more practical work, the location of sensor was chosen based on accessibility in gear box.
2 EMPIRICAL MODE DECOMPOSITION ALGORITHM
In the past, every customary vibration signals which were acquired for fault detection by every tools were analyzed by means of Fourier analysis, wavelet transformation, etc. In all these types of processes it was assumed that signals are based on linear stationary signal. However, against all assumptions in real world, the generated vibrational signals are nonstationary and have transient nature in time frame Loutridis (2006S.J. Loutridis, (2006) Instantaneous energy density as a feature for gear fault detection, Mechanical Systems and Signal Processing 20: 12391253).
In the presented method, it assumed to all signals were generated by number of IMFs which is oscillated in time platform. Each generated IMF has number of extremums which allow linearly and nonlinearly to vibrate. Every point of oscillated signal represents simple vibration with same number of extremum and zero crossing for every IMF. In addition, each oscillated IMF is symmetric with respect to local average which makes symmetric signals with regard to horizontal local average axis and it is possible for each point of signal to be member of another IMF. For original signal being composed of some IMF will results in time superposition which makes a complicated signal decomposed in more physical frequency meaning Huang and Shen, (2005Norden E Huang, Samuel S P Shen, (2005) HilbertHuang Transform and Its Applications, World Scientific Publishing Co. Pte. Ltd Volume 5).
Each signal has its own local scale which defines the inherent physical meaning of what was done during data acquisition. Considering to oscillation in local scale is one of the advantages of the EMD algorithm, which allow decomposing the signal with respect to the frequency and rather than basic signal frequency, higher magnificent frequency result can be concluded. In this algorithm each IMF represent the oscillating mode and constructed by means of following instruction:

In all data the number of extermums and zero crossing must be equal or difference in one extermum or one zero crossing.

In each point mean value of defined data with respect to the envelope of local maximum and envelope of the local minimum must be zero.
The simple harmonic functions of former decomposition method represent the constant amplitude and constant frequency modes. On the contrary, IMFs include variable amplitude and variable frequency. As it seems in Figure 1, all local extermums make envelope by green color and set of local maximums form one cluster which connected to each other by means of third orders spline. Then, by using this procedure again, set of local minimum connected to each other in another same independent cluster by aim of third orders spline. Now both of upper and lower sides of data envelope were built as illustrated in Figure 1 and by calculating the difference between upper and lower envelope first and higher frequency IMF can be extracted.
Original signal with upper and lower push and first step of mean value Huang and Shen, (2005Norden E Huang, Samuel S P Shen, (2005) HilbertHuang Transform and Its Applications, World Scientific Publishing Co. Pte. Ltd Volume 5).
Let m _{1} be the difference between two envelopes so h _{1} can be produced by using the simple subtraction which is shown in Equation 1.
The produced h _{1} and initial signal is shown in Figure 2. Practically the generated extermums which obtained after running one step of this procedure make reveal appropriate mode extermums which initially were lost into the original signals. In the other words, repeating the sifting process make more improvement with low amplitude signals and it makes signal more symmetric and eliminates the dominant low frequency mode.
Original signal with first step difference with mean value Huang and Shen, (2005Norden E Huang, Samuel S P Shen, (2005) HilbertHuang Transform and Its Applications, World Scientific Publishing Co. Pte. Ltd Volume 5).
If h _{1} does not part of IMF, the sifting process has to proceed for k step till the h _{1} _{k} meet the criteria which define the step to stop. After that h _{1} _{k} has become the first point of the IMF and can be represented by c _{1} as first IMF or c _{1} as i'th IMF. As it is shown in Figure 3, the first IMF produced after 12 iterations of the sifting process.
The result of first step sifting process Huang and Shen, (2005Norden E Huang, Samuel S P Shen, (2005) HilbertHuang Transform and Its Applications, World Scientific Publishing Co. Pte. Ltd Volume 5).
In 1998 the first stopping criteria was proposed by Haung, Huang and Shen, (2005Norden E Huang, Samuel S P Shen, (2005) HilbertHuang Transform and Its Applications, World Scientific Publishing Co. Pte. Ltd Volume 5). They were proposed to use normalized square difference for two sequential steps to stop the sifting process. This criteria as it is shown in Equation 4, was based on convergence Cauchy test. In the year of 2003, according to the Haung researches the value of SD_{k} was proposed to be 0.0002. If the calculated value of SD_{k} was greater than 0.0002, the sifting process went to the next level. Otherwise, it stopped and original signal had to subtract by created IMF in order to generate the residual signal for other IMF. As it was shown in Equation 3, the residual signal created with respect to the extracted IMF. For more illustration, original signals and residual signal after extraction of the first IMF is shown in Figure 4. Finally, Equation 6 is applied to take summation of all IMF's and final residual signal to reconstruct the original signal.
Original signal with first step residual Huang and Shen, (2005Norden E Huang, Samuel S P Shen, (2005) HilbertHuang Transform and Its Applications, World Scientific Publishing Co. Pte. Ltd Volume 5).
2.1 Numerical Test
In order to demonstrate the power of the EMD algorithm by means of MATLAB software, Set of artificial sinusoidal combined signal is created which is based on Equation 7. This set of data is included different frequency sine wave which is varied from 20 Hz to 400 Hz and also has different amplitudes to examine the EMD algorithm. The artificial signal is shown in Figure 5 which can be interpreted in terms of time in second.
After decomposition process which is taken a few times, the artificial signal is decomposed into six IMF which are in descending order of frequency variable. In Table 1, the number of iteration for each IMF is shown. Table 1 shows that after how many times iterations the signals were decomposed into six IMF's. It is important because it shows how fast and simple EMD can decompose the signals. Besides, it shows that how fast basic IMFs are decomposed (first, second and third IMF's) they take respectively 4, 1 and 4 iterations to be decomposed in comparison to other IMF's which are shown the interaction of basic IMF's.
The first triplex cluster decomposed waves include three basic frequency spectrum. The first cluster is shown in Figure 6. In this figure, frequency from upper graph to lower graph varies from 400, 200 and 40 Hz, respectively. It is clear that EMD algorithm can decompose the original signal into the basic elements of itself properly. However, these elements have some interaction with each other which are decomposed in another triplex cluster of IMFs and one residual signal. The second cluster and residual signal are shown in Figure 7. In fact, second cluster IMFs are the excess decomposition of EMD algorithm and technically these are not the basic components of original signal. It should be noted that second cluster and residual signal have meaningless frequency. They are just the result of basic element interaction. It should be mentioned that the residual signal shows the general attitude of signal which is more related to the strongest frequency component of signal. For instance, the sine wave at frequency of 400 Hz and amplitude of 2 is dominant component which defines general trend of the signal.
First cluster decomposed IMF's of numerical test in different frequencies a: 400Hz, b: 200Hz and c: 40Hz.
2.2 Noise Effects on Signal Decomposition
In real application of decomposition method, noise and its effect should be taken into account as epidemic existing factor in signal processing which may cause different kinds of problem. Equation 8 determines that the ratio of standard deviation of noise and signals is 0.1 in the same words it's actually define the power or frequency content of the noise. So, By using the equation 8, the added noise is defined as a high frequency component of signals that was generated while it was measuring. In the other words, data acquisition tool induce the high frequency component into the signal. It happens because measuring tools are the main cause of noise generation. In order to examine the power of EMD algorithm, first of all, a randomly definite noise based on Equation 8 is added to original signal. In order to have more realistic simulation the value of added noise is based on measuring data apparatus. The conventional data acquisition tools gather data with same amount of adding noise. The usual and industrial data acquisition tools corporations claim that the ratio of standard deviation of pure signal over noise is equal to 0.1. Then, noise is added to signal which is showed in Figure 8 and is going to separate by EMD algorithm. Therefore, the behavior of EMD algorithm in the presence of noise is determined. The noise which is added to the signal is generated by using Rand function in MATLAB software.
a: numerical test signal, b: noise, c: noisy numerical test signal and d: magnified comparison of noisy and pure numerical test signal.
In the first graph (part a) of Figure 8, the generated signal is shown. Besides, in part b, c and d of Figure 8, added noise, noisy signal or impure signal and zoomed pure signal versus noisy signal are shown respectively. The added noise has to be a zero mean graph which makes the noise symmetric with respect to pure signal. Finally, in the lower graph the original pure signal in the presence of the noisy signal is shown in order to the effect of noise in miniature scale become more determined.
In Figure 9, the first six IMF is shown. As it is shown, every six IMF is useless noise. They occur because of frequency order of EMD algorithm which is decomposing signal from higher frequency to lower. Added noise becomes one of the signal components and it has to be decomposed like the basic frequency part of signal. Noise always has higher frequency rather than basic frequency of signal. It also indicates that EMD algorithm treat signal as a denoising low pass filter like Butterworth low pass filters with different cut off frequency per each IMF.
Noise effect on EMD, all first IMF to sixth IMF's are noise in a, b, c, d, e and f, respectively.
Figure 10 demonstrates important component of signal which is consisting of four graphs. The graph a from Figure 10 shows incomplete part of higher component frequency of signal. This imperfection is done because of interaction with induced noise. In graph (a) of Figure 11 by means of FFT algorithm, the frequency spectrum of C _{7} IMF is shown. In this spectrum, the presence of higher component of original signal is illustrated and noise effect of signal is detected by some other peak frequency in the vicinity of 200Hz. In graph b from Figure 11, there are some evidence which indicate the presence of 200Hz frequency component and a little of 100Hz frequency component which are formed IMF C _{8} in graph b of Figure 10. In the other words, the higher frequency portion of signal is composed of incomplete IMF C _{7} and incomplete IMF C _{8}. However, some missed part of IMF C _{8} can be completed by IMF C _{7}. By looking at graph c from Figure 11, it is obvious that IMF C _{9} in Figure 10 is the second higher frequency part of original signal which is consists of only 100Hz frequency part. However, some part of this component is missed and by reconsidering graph b from Figure 11, cause of missing some part of IMF C _{9} becomes clear. Also, in graph b from Figure 11 there are two frequency peaks which are 200Hz and 100Hz and little missing part of IMF C _{9} is in the context of IMF C _{8}. At last the third part of frequency component is properly decomposed in IMF C _{10} which is shown in graph d from Figure 10 and its frequency spectrum is also shown in Figure 11.
3 BASIC CHARACTERIZATION OF SIMULATION
In order to run more reliable simulation which is more adapted to actual world, the used gear mesh is according to real productive automobile corporation standard which is now available in the market. In order to simulate the gear mesh in high velocity with respect to what research have done till now, the fourth gear engaged is chosen. In Figure 12 (Upper side), the input shaft and its installed gear is shown. In this figure, the input shaft which is the initial receiver of engine torque and angular velocity is shown. In Figure 12 (Upper side), first, second, rear and fourth gears from left to right also are shown, respectively.
Locations of gears and bearings on input shaft (upper side) and output shaft (lower side).
Figure 12 (lower side) depicts the output shaft. In this figure, the third and fourth gear from left to right is shown, respectively. There are two places for bearing to locate. One of them is hard access and another one is easy access. In order to have a better simulation the acceleration sensor in MSC Visual Nastran software is attached with respect to the accessibility of bearing and also, for output shaft the data acquisition is attached on accessible bearing.
The position of gears over the input and output shaft and the position of accessible and inaccessible bearings are shown Figure 12. Accessible bearing on input shaft and accessible bearing on output shaft are at the distance of 110.725 mm and 112.486 mm, respectively from fourth gears. In order to constraint the two shafts; four bearings are applied in the geometrical models which are placed at the real position of the bearings with the equal length to real bearings. Each shaft has two ball bearings which are placed at the end and beginning of the shaft to make the shaft under constraint and all of the ball bearings have low friction coefficient. The friction of the bearings is calculated by means of wet friction factor which is according to the SKF catalog. The friction coefficient is set to 0.002 and the clearance is set to 0.05 mm between ball bearings and outside of the shaft.
The technical properties' of couple gears are shown in Table 2. According to this table, the factor of this engagement with considers the number of tooth of input and output fourth gears, is 1.05.
In order to have correct RPM and correct input and output torque values which had to be exerted on input and output shaft respectively, the practical engine map is useful tool. Engine map includes the generated engine output torque versus revolutionary rate of engine; the fourth gear conventional velocity lay between 60Kmph to 80Kmph. So, first of all the rotating velocity of input gear has to be determined. For this purpose the CARSIM software is chosen to simulate the really feature of practical vehicle to calculate the engine RPM at 80Kmph and 60Kmph velocity of the vehicle and the result of this simulation is shown in Table 3. The chosen RPM is 2674 rpm which is equal to 44.57 Hz.
In order to determine the input and output exerted torque by using the engine map, the output torque of engine which is exerted on input shaft is read as 144 N.m. The engine map which is used for this simulation is also shown in Figure 13. For determination of output torque on output shaft with respect to the first assumption that is constant linear velocity, the force equilibrium condition is considered and output torque by using the gearbox ratio at fourth gear and input torque is calculated 136.98 N.m.
In Figure 14, the position of sensor, the length of bearings, the point where the torque exerted (input and output torque respectively on input and output gearbox) and also the scheme of assembly are shown. In practical mechanism, all ball bearings have to hold moving components at their places and allow them to rotate freely without friction forces or torques. Therefore, in the procedure of simulation the fixed constraint put on the bearings which is prevented the bearing from vary and also it should be noted that, between all six components there are impenetrate constraint which force whole six geometrical model not to penetrate each other. In order to have more effective control over impulse between all models, restitution coefficient is set to 0.8. All four ball bearings were fixed at each place and into the CARSIM software each shaft and its bearings were defined to be impenetrable to each other. This being impenetrable defined a rule which allows software to make the collision between each shaft and specific bearings with define restitution and friction factor.
3.1 Mesh Independency Analysis
In order to have best feasible size of mesh and to reach the minimum possible element size, at the first step above of all, simple static structure model with 100 N.m and different element size is run. As it is shown in Figure 15, the place where the engine exerts, 100 N.m torque is located. The only axial rotations of six degree of freedom are opened to the model and the output place of output shaft is fixed. Finally, for the boundary conditions at right hand bearings, the axial rotation of shaft is the only freedom degree of the model.
For the first step of mesh independency, the geometrical model is meshed at 10 mm lengths element which is generated 4954 and 18540 for number of nodes and number of elements of output shaft, respectively. Because it is supposed that fault is occurred in input shaft, the only input shaft is meshed in order to reduce the time and calculation source. The amount of torque which is exerted to the mesh model is 100 N.m. In Figure 16, the maximum principle stress of input shaft is shown at 0.25 mm mesh length. In the graph of Figure 17, the maximum principle stress versus the element length is plotted and the length which maximum stress is converged is chosen to remesh the output and input shaft again and make both of them ready to dynamic simulation. The chosen length for element length is 0.5 mm which is created the 559534 and 2726092 number of nodes and elements, respectively. The meshed geometrical model is also shown in Figure 18.
4 DYNAMIC SIMULATION OF ENGAGED GEARS
In this step the geometrical model is imported to the MSC Visual Nastran software and the constraint is added to the model to complete the simulation process. At the first step, the gear without fault is imported to reach the conventional no fault vibration as reference pattern. As it mentioned, the input gear is meshed with 49424 elements and 8290 nodes of tetrahedral geometrical shape. Then, the output gear is imported and meshed with 43846 elements and 74048 nodes. Then, for two gears in order to having no penetration in each other, the collision between them is defined. The next step is the import of four bearings. The four bearings are imported and in order to make sure of being bearings immobilized the fix constrain make them freeze. Fix constraints do not mean that the rotations of two shafts are ignored, but it means the all bearing in their position are fixed and they do not have any motion. The collisions between each bearings and relevant shaft are defined by means of constitution factor and friction factor of bearings which are 0.8 and 0.02, respectively. Then, 144 and 136.98 N.m torques are added to input and output shaft, respectively. In order to make the system rotation the motor on input shaft is used to create the angular velocity of 16045.2 Deg/s which is equal to 44.57 Hz. By using the meter tools in MSC Nastran software angular position, linear translation and linear acceleration of output shaft are recorded by data acquisition frequency of 524288 Hz. The pure signal of angular position, lateral position of axis and linear acceleration of output shaft is shown in Figure 19.
In order to have more real signals, the noise should be added to the signal. Therefore, based on Equation 8, based on the sensor production company catalogue, the noise is added to the pure signal. According to the sensor production company catalogue, the ratio of standard deviation of recorded noise per standard deviation of pure signals is less than 0.1. So, by using Equation 8, the noise is added to the pure signal and make it impure and ready to analyze. For more illustration, the impure signal is shown in Figure 20 at two graphs. The first upper graph is impure signal and for comparison in more detail the pure signal and impure signal is depicted at lower graph in high level zoomed view.
a: Noisy no faulty signal and b: noisy no faulty signal in comparison of pure no faulty signal of input shaft.
4.1 Fault Analysis in One Gear Tooth
One of the conventional fault that take place in gear mesh is one teeth removed suddenly. Fluctuating transmitting torque and angular velocity is one of the major cause which make one or more teeth removed or can cause crack initiation. So as the first fault, one teeth of gear teeth is removed in order to make the input gear imperfect. As it is shown in Figure 21, in order to induced the fault in gear teeth the pocket command in CATIA software is used and one of the teeth arbitrary is removed. After meshing process the entire geometrical model contains 42891 elements and 72457 nodes.
The entire simulation steps for this kind of fault are the same as no fault simulation. In Figure 22, pure signal of angular position, lateral position of axis and linear acceleration of input shaft of one teeth fault are shown.
In Figure 23, the impure faulty signal is shown at upper section and for more comparison impure and pure fault signal is also shown in lower section.
a: Noisy faulty signal and b: noisy faulty signal in comparison of pure faulty signal of faulty input shaft.
4.2 Fault Analysis of Partially One Gear Tooth
Like previous section, the other possibility of fault occurrence is partially tooth removed on effect of abruptly change in transmitting torque and angular velocity in which particular conditions always leads to have critical condition. In order to make fault on gear, afourth of one teeth helical length is removed as it mentioned by using the pocket command in CATIA software. The length of removing teeth is 3.9 mm which is afourth length of whole helix length in nofault gear. In Figure 24, the partially removed teeth is shown.
Just like pervious section, the entire simulation steps are done with this kind of fault and in Figure 25, pure signal of angular position, lateral position of axis and linear acceleration of input shaft of a fourth one teeth fault is shown, respectively. In Figure 26, the impure signal of a fourth teeth removed is shown at upper section and for more comparison impure and pure faulty signal is shown in lower section.
a: Noisy partial faulty signal and b: noisy partial faulty signal in comparison of pure partial faulty signal of faulty input shaft.
5 SIMULATION RESULTS
In the previous section, the three sets of pair gears were simulated and the noise was added to all data with respect to the Equation 8. In this section, at the first step the added noise should be removed. For this purpose, the 5th kind of Butterworth low pass filter is used with cut off frequency of 5 KHz and sampling frequency of 524288 Hz is applied. The low pass filter parameter is calculated by try and error process. In order to avoid frequency damage to data, the cut off frequency value is chosen in a way which is far away from frequency information. For more illustration, linear frequency response in decibel unit and linear phase response in radian unit is depicted in Figure 27. In the gear fault detection the important and usable frequency data is formed at the gearmesh frequency of defected gear, which is production of teeth number with rotation frequency. The gearmesh frequency is about 1800 Hz. So, in order to avoid missing important information the cut off frequency of low pass filter is selected as 5 KHz.
5.1 Analysis of No Fault Gear Engaged Simulation
For no faulty gear analysis, the first step should be denoising generated data. So, as it was mentioned in order to denoise the data the low pass 5th Butterworth filter with 5 KHz cut off frequency was applied to data and this result is shown in Figure 28 and for more illustration the figure have two segments include the unzoomed and zoomed view.
In order to reveal the fault at the second step, by using the FFT algorithm the frequency spectrum of undecomposed no faulty gearmesh data is obtained. The frequency spectrum of no faulty gearmesh is shown in Figure 29. It is obvious that there is no sign of peak frequency at 1738.23 Hz. So, in order to dig deeper and get more details from denoised data, the decomposition process should be run.
On the other hand, to obtain more resolution detail of data, the EMD process is employed to decompose the denoised signal. After decomposition, the nine IMF's are obtained with different range of frequencies on horizontal axis. After taking the FFT frequency spectrum of each IMF's, the range of frequency of each IMF's is shown in Figure 30.
As it is shown in Figure 30, the second and fairly third IMFs are capable to have fault frequency peak. So, in order to find fault information the first category of four IMF's which has higher frequency range is shown in Figure 31 and frequency spectrum of these four IMFs is depicted in Figure 32.
FFT spectrum of a: first IMF, b: second IMF, c: third IMF and d: fourth IMF of no faulty signal.
As it seems neither in frequency spectrum of IMF's of first category of four IMF's, nor in the other IMF's and residual signal, there is no sign of frequency peak at 1738.23Hz and it means there is no fault in gear mating.
5.2 Fault Analysis of One Tooth from One Gear
In this section like the previous one, analysis of one gear tooth faulty signal is discussed to find the trace of fault in frequency information. At the first step, using butter worth low pass filter the denoising process is done with 5 KHz cut off frequency and its result is also shown in Figure 33 in two moods from top to bottom, i.e., the zoomed and unzoomed view, respectively.
In Figure 34, in order to obtain the frequency information the frequency spectrum is shown. However, as it seems there is no dominant frequency at vicinity of 1738 Hz. So, this kind of frequency peak is not reliable to make sure of fault detection. Then, in order to reach more resolution and more details by using the EMD algorithm the original denoised signal is decomposed to nine IMF's.
In Figure 35, first category of four IMF's are shown which have the frequency range between 500 Hz to 4500 Hz. In this IMF's there is the sideband frequency of gear mating which is generated by the effect of impacting of teeth when they come to mesh.
As it is shown in Figure 36, the frequency spectrum of second IMF, there is peak frequency with magnitude of 1738 Hz which is the gearmesh frequency of output gear. This peak of frequency is the indication of fault existing in gear meshing. This peak is indicated at output gearmesh frequency. So, this fault is belonging to the output gear. The input gear has 41 teeth and the output gearmesh frequency is 1827.37 Hz and at 1837.37 Hz there is no reliable peak frequency.
FFT spectrum of a: first IMF, b: second IMF, c: third IMF and d: fourth IMF of one teeth faulty signal.
5.3 Fault Analysis of Partially One Teeth of Gear
In Figure 37, the denoised signal is shown at first upper segment of figure and the second IMF which is more important from others is depicted at middle segment of Figure 37 and at last there is frequency spectrum of second IMF at bottom segment of Figure 37.
a: filtered signal, b: second IMF and c: FFT spectrum of second IMF of partial faulty signal.
Just like whole tooth fault, at partially teeth fault there is frequency peak at output gearmesh frequency which indicates the present fault in the output and because of this frequency peak is the dominant peak like the other peak, it is reliable sign of fault. It should be noted that the dominancy of whole teeth fault is more recognizable. So, with respect to the partially teeth fault the whole teeth fault is the prominent fault in frequency spectrum.
6 CONCLUSION
In the present work, the fault detection of two mating gears was discussed and the fault was detected by means of signal processing method. First of all, the real and practical gear which is used in the industry was provided using CATIA software and then, by means of MSC Nastran software the no faulty gearmeshing simulation in practical speed of gear working and practical gear torque transmission was done and for more realization the noise was added to all generated signals. Then, the generated signal was processed using MATLAB software to reach the pattern of FFT spectrum. In the processing stage, the Butterworth low pass filter with cut off frequency of 5 KHz was applied to signal to denoise the entire nuisance signal and it was completed the signals purification process. Then, the EMD algorithm was decomposed the signal into the IMF's and each frequency IMF spectrum were shown the pattern of no faulty gearmeshing process and these results were saved as healthy pattern reference. In next stage of this research, in order to reach the faults pattern references, the entire previous stage was implemented in the presence of fault which is consist of two faults. First, the whole tooth broken and then, the partial tooth broken and these kinds of faults could be revealed by using EMD algorithm and FFT spectrum. After the simulation and signal processing were done, the following results are reported as follows:

The location of sensor was chosen at distance of 112 mm from input gear which may cause the more complicated signal with regard to shaft properties. This altercation in the location affect the signal and make it AMFM signal and reduced the strength or power of signals because of long distance to the acquisition location and it also may leave resonance effect on signals and randomly changes the frequency spectrum.

If it is just necessary to the denoise signal without decomposition, the EMD algorithm can denoise signal, but if it is necessary to denoise and then decompose the signals, the EMD algorithm could not be a useful tool.

In nonstationary signals EMD algorithm in comparison with wavelet which is based on constant sifting function, has an advantage to decompose the signal into basic intrinsic functions which are lay into the signal.

In the existence of noise, the EMD algorithm has bad behaviour in a way that the noise are decomposed into initials IMF's and absolutely the noise should be removed before decomposition process. The noise could leave a bad effect in decomposition process.

Detection of one tooth broken fault has high resolution with respect to the partial tooth broken.
Reference
 B. Samanata, (2004) Gear Fault Detection Using Artificial Neural Networks And Support Vector Machines With Genetic Algorithms, Mechanical Systems and Signal Processing 18: 625644
 Chuan Li, Ming Liang, (2012) TimeFrequency Signal Analysis for Gearbox Fault Diagnosis Using a Generalized Synchrosqueezing Transform, Mechanical Systems and Signal Processing 26: 205217.
 D. Chen, W. J. Wang, (2002) Classification of wavelet map Patterns using MultiLayer neural networks For gear Fault Detection, Mechanical Systems and Signal Processing 16(4): 695704
 Dejie Yu, Yu Yang, Junsheng Cheng, (2007) Application of TimeFrequency Entropy Method Based on HilbertHuang Transform to Gear Fault Diagnosis, Measurement 40: 823830.
 Grzegorz Litak, Krzysztof Kecik, Rafal Rusinek, (2013) Cutting force response in milling of Inconel Analysis by wavelet and HilbertHuang Transforms, Latin American Journal of Solids and Structures 10: 133140.
 H. Zheng, Z. Li, X. Chen, (2002) Gear Fault Diagnosis Based on continuous Wavelet Transform, Mechanical Systems and Signal Processing 16(23): 447457
 Hui Li, Yuping Zhang, Haiqi Zheng, (2011) Application of Hermitian Wavelet to Crack Fault Detection in Gearbox, Mechanical Systems and Signal Processing 25: 13531363.
 Junsheng Cheng, Kang Zhang, Yu Yang, (2012) An Order Tracking Technique for the Gear Fault Diagnosis Using Local Mean Decomposition Method, Mechanism and Machine Theory 55: 6776.
 Junsheng Cheng, Yi Yang, Yu Yang, (2012) A Rotating Machinery Fault Diagnosis Method Based on Local Mean Decomposition, Digital Signal Processing 22: 356366.
 Junsheng Cheng, Yu Yang, Dejie Yu, (2008) Application of frequency family separation method based upon EMD and local Hilbert energy spectrum method to gear fault diagnosis, Mechanical Systems and Signal Processing 43: 712723.
 Junsheng Cheng, Yu Yang, Dejie Yu, (2010) The Envelope Order Spectrum Based on Generalized Demodulation TimeFrequency Analysis and Its Application to Gear Fault Diagnosis, Mechanical Systems and Signal Processing 24: 508521.
 K.S. Wang, P.S. Heyns, (2011) An Empirical ReSampling Method on Intrinsic Mode Function to Deal with Speed Variation in Machine Fault Diagnostics, Applied Soft Computing 11: 50155027.
 Li Li, Liangsheng Qu, Xianghui Liao, (2007) Haar wavelet for machine fault diagnosis, Mechanical Systems and Signal Processing 21: 17731786.
 Lourdes Rubio, Belen MunozAbella, Patricia Rubio, Laura Montero, (2014) Quasistatic numerical study of the breathing mechanism of an elliptical crack in an unbalanced rotating shaft, Latin American Journal of Solids and Structures 11: 23332350.
 Norden E Huang, Samuel S P Shen, (2005) HilbertHuang Transform and Its Applications, World Scientific Publishing Co. Pte. Ltd Volume 5
 Roberto Ricci, Paolo Pennacchi, (2011) Diagnostics of Gear Faults Based on EMD and Automatic Selection of Intrinsic Mode Functions, Mechanical Systems and Signal Processing 25: 821838.
 S.J. Loutridis, (2006) Instantaneous energy density as a feature for gear fault detection, Mechanical Systems and Signal Processing 20: 12391253
 W. J. WANG, P. D. Mcfadden, (1996) Application Of Wavelets To Gearbox Vibration Signals For Fault Detection, Journal of Sound and Vibration 192(5): 816_828
Publication Dates

Publication in this collection
Nov 2016
History

Received
13 July 2015 
Reviewed
08 May 2016 
Accepted
31 May 2016