Detection and Localization of Tooth Breakage Fault on Wind Turbine Planetary Gear System considering Gear Manufacturing Errors

Sidebands of vibration spectrum are sensitive to the fault degree and have been proved to be useful for tooth fault detection and localization. However, the amplitude and frequency modulation due to manufacturing errors (which are inevitable in actual planetary gear system) lead tomuchmore complex sidebands.Thus, in the paper, a lumped parameter model for a typical planetary gear system with various types of errors is established. In the model, the influences of tooth faults on time-varying mesh stiffness and tooth impact force are derived analytically. Numerical methods are then utilized to obtain the response spectra of the system with tooth faults with and without errors. Three system components (including sun, planet, and ring gears) with tooth faults are considered in the discussion, respectively. Through detailed comparisons of spectral sidebands, fault characteristic frequencies of the system are acquired. Dynamic experiments on a planetary gear-box test rig are carried out to verify the simulation results and these results are of great significances for the detection and localization of tooth faults in wind turbines.


Introduction
Planetary gear systems have been widely used in wind power systems because of the advantages of compact structure, large carrying capacity, and high transmission efficiency [1].In recent years, the tooth faults that have occurred in planetary gear systems have brought numerous troubles to wind power plants [2][3][4].Therefore, fault diagnosis of the planetary gear system is of great significant to the safe operation of wind turbine.
Tooth pitting, spalling, cracking, and breakage are some of the common fault modes that have occurred in wind turbine planetary gear systems.Detection and localization of those faults are of great significance.Vibration based diagnosis is one of the most effective and frequently used health monitoring technologies.In recent years, much attention has been paid on analyzing the vibration signal in time and frequency domains.To name a few, Lei et al. [5] extracted two diagnostic parameters based on the examination of vibration characteristics of a planetary gearbox in both time and frequency domains.Experiments revealed that the proposed diagnostic parameters performed better than other parameters.Lin and Zuo [6] introduced a method based on the Morlet wavelet time-frequency analysis.The method was found more effective in detection of tooth cracks than two other types of discrete wavelet transform.Liu et al. [7] carried out fault diagnosis based on the local mean decomposition (LMD) method.The method was suitable for obtaining instantaneous frequencies in frequency domains.
The above-mentioned studies are helpful for the fault detection of planetary gear system.In addition, the sidebands of vibration signals also attract the attention of many researchers, because the sidebands contain abundant information which is useful for the detection and localization of faults in planetary gear systems.McFadden and Smith [8] found that the modulation sidebands of the planetary gear system are typically asymmetric and the frequency with largest amplitude is the multiple of the carrier rotation frequency.Based on the research, McFadden [9] presented a method to calculate the time domain average of the vibration signal and the method for the separation of planet vibration and sun gear vibration is extracted.Feng and Zuo [10] deduced the relationship between fault characteristic frequencies and rotating frequencies of the components in a planetary gear system.The complicity of sidebands and fault characteristic frequencies were verified by simulations and experiments.Chaari et al. [11] modeled a planetary gear system which contained the fault modes of tooth pitting and cracking.The dynamic responses of the fault systems and a healthy gear system were compared in frequency domains.Chaari et al. [12] also studied the influence of tooth pitting and cracking on the sun gear vibration spectrum based on a planetary gear model.Sidebands were found to be induced because of amplitude modulation.Cheng et al. [13,14] estimated the degree of gear pitting and cracking by comparing the dynamic response of a fault system with a health system.A model containing tooth spalling and cracking was developed by Jia and Howard [15].His research indicated that the amplitude and phase modulation of the vibration signal could indicate the difference between tooth spalling and cracking effectively.
However, various gear manufacturing errors were not considered in the above-mentioned investigations.Typical gear errors, such as eccentricities, pitch-line run-out errors, tooth spacing, and indexing errors, are inevitable in actual planetary gear system.Ligata et al. [16] pointed out that the aforementioned errors of a planetary gear set might cause modulations to the sidebands of the system.Inalpolat and Kahraman [17] developed a nonlinear time-varying dynamic model to predict modulated sidebands of a planetary gear system, and both the numerical and experimental results showed that the sidebands were greatly changed due to the amplitude and frequency modulations caused by gear errors.As the gear errors can influence the sidebands significantly [16,17], the fault characteristics might be greatly changed after considering the gear errors.Therefore, it is necessary to take the influence of errors into consideration to increase the accuracy and precision of the tooth faults diagnosis.
In this paper, a lumped parameter model of a typical planetary gear system is established.The influences of typical tooth faults on the time-varying mesh stiffness and tooth impact forces are derived analytically.Numerical methods are utilized to obtain the response spectra of the system with tooth faults with and without gear errors.Three system components (including sun, planet, and ring gears) with tooth faults are considered, respectively, in the discussion.Through detailed comparisons on spectral sidebands, fault characteristic frequencies of the system are acquired.Dynamic experiments on a planetary gear-box test rig are carried out to verify the simulation results.

Dynamic Model of Planetary Gear System
Figure 1 shows the schematic diagram of the dynamic model for a planetary gear system.The system composes of sun gear (), carrier (), ring (), and three planets ().Bearings are considered as linear springs and gear meshes are considered as linear springs acting on the lines of action.Each component has three degrees of freedom, transverse (), longitudinal (), and torsional degrees ().The reference coordinate system rotates with the carrier and the origin of the coordinate coincides with the center of the carrier.The circumferential positions of planets are specified by the fixed angles    (Figure 2) measured relatively to the rotating coordinate system with   1 = 0.
The global equation of motion for the system can be simplified as [18,19] where q represents the vector of the degrees of freedom Ω  is the rotating speed of the carrier.Due to the unstable wind conditions, frequent low speed start and emergency brakes, the rotating speed is often nonstationary.However, if the concerned time period is small, the work condition of the wind turbine could be considered as quasi-stable [20].Thus, the premise of constant rotating speed is used in this paper.M, G represent the mass and gyroscopic matrices, respectively, and are expressed by where   and   are the mass and the base radius of part ,   is the inertial moment of part  with respect to its rotational axis.
K  is the bearing stiffness matrices and are expressed by where   ,   ,   represent the bear supporting stiffness long the three degrees of freedom.
The centripetal stiffness matrix due to the carrier rotation is denoted by K Ω .Consider ,  = , , , 1, . . ., . ( The time-varying mesh stiffness matrix is represented by K  ().The matrices are expressed as where where    () ( = , ) is the meshing stiffness between the part  and the th planet.   is the bear supporting stiffness of the th planet in the transverse and longitudinal direction.  is the pressure angle of part .   =   −   and    =   +   ,   is illustrated in a sun-planet mesh (Figure 2). = 1, 2, . . ., ,  is the number of the planet gears.T() represents the force matrix and can be expressed by where   ( = , , ) is the external torques applied on part .  is the base radius of part .F  () is the impact force caused by faults, the mathematical expression of which will be given in Section 3.
F  () is the excitation force induced by manufacturing errors and can be written as where () ( = , ) is the errors between the part  and the th planet and can be written as where    () and     () ( = , ) are errors associated with rotational frequencies of gears.A unified mathematical formula which is able to express those errors (gear eccentricities, pitch-line run-out errors, tooth spacing, and indexing errors) was extracted by Inalpolat and Kahraman [17], where   and    are the amplitudes of these errors,  mesh is mesh frequency,   and   are tooth numbers of the component  and planet,   and    are phase angles,   is relative position angle, and the arithmetic sign of the pressure angle   is plus for  =  and minus for  = .In (11a),    and    are gear transmission errors.Inalpolat and Kahraman [17] believed that those errors were amplitude and frequency modulated by    () and     ().Therefore, another equation is able to express the gear transmission errors as given by Inalpolat and Kahraman [17], where  [21] showed that the peak-to-peak amplitude of    () might experience up to 4-6% variation due to the influence of    () and     () even with the modest error magnitudes.

Modeling of Tooth Faults
In this section, the typical tooth fault of tooth breakage is modeled to illustrate the influence of tooth faults on gear mesh.The tooth breakage can be modeled by the reduction of the time-varying mesh stiffness and the generation of additional impact forces (F  ()).

The Influence of Tooth Breakage on Time-Varying Stiffness.
As shown in Figure 3, a tooth is often considered as a cantilever beam with varying section when calculating deformation.Deformation consists of five parts: bending, compressive, shear, fillet-foundation, and contact deformations.
The total deformation of the point  is expressed by where   , V  ,   ,   ,   represent bending, compressive, shear, fillet-foundation, and contact deformations, respectively, which can be calculated by the method used in [12],   (see Figure 3) defines the angle between the tooth normal force and the tooth centerline.The stiffness of one gear tooth at point  is calculated as   = /  .The total stiffness of one pair of meshing tooth is obtained by according to the serial relation.The stiffness when doubleteeth meshing is  =  pair1 +  pair2 according to the parallel relation.
The parameters of the studied planetary gear system are listed in Table 1.Young's modulus  = 2.06×10 11 Pa, Poisson's ratio  = 0.3.Figure 4 shows the mesh stiffness of the sun gear and ring gear with planet 1 (  1 and   1 ).The planet 1 with and without tooth breakage is considered and the mesh stiffness is seriously affected by tooth breakage.

The Impact Force
Generated by Tooth Breakage.The impact force is generated during the approaching and separation of the meshing teeth [22].As one gear tooth is broken, the normal gear meshing process is interrupted.When the next gear tooth is going to mesh, additional impacts would occur between the driven and driving gears.The dynamic response of the system will be influenced by the additional impact force.Thus, only the change of mesh stiffness could not reflect the effect of tooth breakage upon the dynamic response of the planetary gear system.In this section, the additional tooth impact forces are modeled.
The additional impact force F   () generated by the tooth breakage of the planet gear can be expressed as (planet 1 is assumed to be the fault gear)  where where F  1 () represents the impact force vector of gear  ( = , ),  1 () ( = , ) denotes the impact force and can be expressed as where   1 and   1 represent the beginning and finishing time of the impact force of gear .The relationship between   () and   () is expressed by   () = ()  ().Because the impact acceleration will be generated, () > 1 when the fault tooth meshes with ring gear and () < 1 when the fault tooth meshes with sun gear.The additional impact force F   () generated by the tooth breakage of the ring gear could be expressed by where where F  () represents the impact force matrix when planet  meshes with ring gear,   () and   () are the pulsating impact forces, the mathematical expressions of which are similar to (15) and could be expressed by Without breakage Breakage

Amplitude
Frequency (Hz) Without breakage Breakage The impact force generated by the tooth breakage of the sun gear is F   (), the mathematical expression of which is the same as that of F   ().However, the values for beginning and finishing times of   () are different.
The research results by Yao and Wei [23] showed that the action time of impact force is usually about 5-10% of one mesh cycle (  ).Wang et al. [24] pointed out that the maximal impact force induced by tooth breakage is about twice the average mesh force.Thus, we can assume    −   =   /10 and   () = 2   = const.

The Dynamic Response of Planetary Gear System without Errors
Physical parameters of the studied planetary gear system are listed in Table   Without breakage Breakage    sidebands are found around the mesh frequency.Then tooth breakage is added to the planet gear, sun gear, and ring gear successively, and the acceleration spectrum (the ring gear in longitudinal direction) is showed in Figures 5(b)-5(d).
Compared with the spectrum of the normal system (without breakage), modulation sidebands appeared in the spectrum of the breakage system and the peaks of the sidebands are at the frequencies of  mesh ± (  +   ),  mesh ± (  −   ) and  mesh ±   ( = 1, 2, . ..) successively.Therefore, the mesh frequency  mesh is modulated by those frequencies, which are caused by the influence of tooth breakage on mesh stiffness and additional impact force.So, the sidebands of the planetary gear system are comparatively simple if the errors are not included in the system.

The Dynamic Response of Planetary Gear System with Errors
In order to display the modulation effect of errors more reasonably, a formula was given by Inalpolat and Kahraman [17] to express the acceleration response of the system where   and   are constants,    () and    () are defined as ,   ,   ( = , , ) can be obtained by the dynamic model.  () is a weighting function for a planet gear positioned at the angle of    and   () can be defined as where () is a Hanning function used to represent the influence of rolling planets on fixed accelerometer (scrolling  (20).Figure 6(a) shows the system acceleration spectrum of a normal system (without breakage).Compared with Figure 5(a), complex sidebands appeared around the mesh frequency due to the existence of errors.They are symmetrically arranged on both sides of the mesh frequency.The sidebands of the system have undergone tremendous changes when errors are taken into consideration.The phenomenon is bound to influence the fault characteristic frequencies of the system components.
Figure 6(b) shows the acceleration spectrum of the system with and without tooth breakage on planet 1.Since the mesh frequency  mesh is modulated by the frequencies of weighting function (which characteristic frequencies are   ) and the frequencies of (  +   ) (caused by the influence of tooth breakage on mesh stiffness and additional impact force), the amplitudes of the tooth breakage system are much larger than those of the normal (without breakage) system at the frequencies of  mesh ±   ± (  +   ) (,  = 1, 2, . ..).Therefore, the fault characteristic frequencies of the planet gear are   =  mesh ±   ± (  +   ) and can be simplified as  mesh ±   ±   (,  = 1, 2, . ..).Likewise, the amplitudes of the sun gear tooth breakage system are found to be larger than those of the normal system (without breakage) at the frequencies of  mesh ±   ±   (,  = 1, 2, . ..) in Figure 6(c).The main reason is that the mesh frequency  mesh is modulated by the frequencies of weighting function (  ) and the frequencies of (  −   ) (caused by the influence of tooth breakage on mesh stiffness and additional impact force).Therefore, the fault characteristic frequencies of the sun gear are  mesh ±   ± (  −   ) and can be as   =  mesh ±   ±   (,  = 1, 2, . ..).
Similarly, because the mesh frequency  mesh is modulated by the frequencies of the weighting function   and the frequencies of   (caused by the influence of tooth breakage on mesh stiffness and additional impact force), the amplitudes of the ring gear tooth breakage system are larger than those of the normal system (without breakage) at frequencies of  mesh ±   ±   (,  = 1, 2, . ..), as shown in Figure 6(d).Therefore, the fault characteristic frequencies of the ring gear are   =  mesh ±   ±   and can be simplified as   =  mesh ±   ( = 1, 2, . ..).

Experimental Validations
Figure 7 shows the test bench of a planetary gear system.The drive motor is controlled by frequency converter and provides power for the system.The speed of the system is slowed down by the speed reducer and speeded up by the speed multiplier.A loading motor continuously pumps high-pressure oil into the loader to induce a workload.The structure and the size of the speed multiplier are proportional to the gearbox of a real wind turbine and it consists of a planetary and a parallel gear boxes.The parallel gear box is composed of two pairs of parallel-stage gearing.The parameters of the planetary gear system are listed in Table 1 and those of the parallel-stage are listed in Table 2.The structure style and parameter setting of the speed reducer are the same as those of the speed multiplier.
As shown in Figure 8(a), an acceleration sensor is mounted on the planetary gear box of the speed multiplier.The tooth breakage of planet gear and ring gear are added to the planetary gear system separately, as shown in Figures 8(b (without breakage) and the two fault modes.The sampling frequency is 16384 Hz. Figure 9 shows the acceleration spectrum of the normal system (without breakage).The drive motor speed is 156.24 rad/s and the input speed of the speed multiplier is 2.77 rad/s.Therefore, the mesh frequency is 34.9 Hz. Figure 9 shows that sidebands appeared around the mesh frequency.The results indicate the existence of errors in the system.6.1.Planet Gear with Tooth Breakage.Figure 10 shows the acceleration spectrum of the system with and without tooth breakage on planet.
In Figure 10(a), the load and the speed of the system are 200 N⋅m and 2.74 rad/s, respectively.Therefore, the mesh frequency is 34.87 Hz.Compared with the amplitudes of the normal system (without breakage), those of the tooth breakage system change significantly.In order to make it more clear, the right side of Figure 10(a) is zoomed in, as shown in Figure 10(b).All amplitudes increase at the frequencies of  mesh ±   ±   (,  = 1, 2, . ..), except for the frequencies of  mesh + 2  and  mesh +   + 3  .
Then, the load is increased to 229 N⋅m and the speed remains the same.The differences of amplitudes between the normal system and the tooth breakage system are increased at the frequencies of  mesh ±   ±   (,  = 1, 2, . ..), as shown in Figure 10(c).
Finally, the load is kept unchanged at a value of 229 N⋅m and the input speed is decreased to 2.22 rad/s.Therefore, the mesh frequency is changed to 27.86 Hz.The amplitudes of the tooth breakage system are also larger than those of the normal system at the frequencies of  mesh ±  ±  (,  = 1, 2, . ..) significantly, as shown in Figure 10(d).
Figure 10 verified the correctness of the simulation results in Section 5.The fault characteristic frequencies of the planet gear with tooth breakage are   =  mesh ±   ±   (,  = 1, 2, . ..).

Ring Gear with Breakage.
Figure 11 shows the acceleration spectrum of the system with and without tooth breakage on the ring gear.
In Figure 11, the load and the speed of the system are 200 N⋅m and 2.77 rad/s, respectively; therefore, the mesh

Figure 1 :
Figure 1: The dynamic model for a planetary gear system.

Figure 3 :
Figure 3: The schematic diagram of a cantilever beam.

Figure 4 :
Figure 4: The mesh stiffness of planet 1 with and without tooth breakage: (a) planet-sun gear pair, (b) planet-ring gear pair.

Figure 5 :
Figure 5: The acceleration spectrum of the ring gear in longitudinal direction in the system with and without tooth breakages (without errors): (a) without breakage, (b) planet gear with and without tooth breakages, (c) sun gear with and without tooth breakages, (d) ring gear with and without tooth breakages.

Figure 6 :
Figure 6: The spectrum of the system with and without tooth breakages (with errors): (a) without breakage, (b) planet gear with and without tooth breakages (the right side of  mesh ), (c) sun gear with and without tooth breakages (the right side of  mesh ), (d) ring gear with and without tooth breakages.

Figure 7 :Figure 8 :
Figure 7: The test bench of the planetary gear system.

Figure 9 :
Figure 9: The acceleration spectrum of the normal system.

Figure 10 :
Figure 10: The spectrum of the system with and without tooth breakage on planet at different speeds and loads: (a) load = 200 N⋅m, input speed = 2.74 rad/s, (b) zoom the right side of (a), (c) load = 229 N⋅m, input speed = 2.74 rad/s, and (d) load = 229 N⋅m, input speed = 2.22 rad/s.

Figure 11 :
Figure7shows the test bench of a planetary gear system.The drive motor is controlled by frequency converter and provides power for the system.The speed of the system is slowed down by the speed reducer and speeded up by the speed multiplier.A loading motor continuously pumps high-pressure oil into the loader to induce a workload.The structure and the size of the speed multiplier are proportional to the gearbox of a real wind turbine and it consists of a planetary and a parallel gear boxes.The parallel gear box is composed of two pairs of parallel-stage gearing.The parameters of the planetary gear system are listed in Table1and those of the parallel-stage are listed in Table2.The structure style and parameter setting of the speed reducer are the same as those of the speed multiplier.As shown in Figure8(a), an acceleration sensor is mounted on the planetary gear box of the speed multiplier.The tooth breakage of planet gear and ring gear are added to the planetary gear system separately, as shown in Figures8(b) and 8(c).Vibration signals are collected for a normal mode

Table 1 :
Parameters of the planetary gear system.  =   =   =   =   =   =   = 2 × 10 8 −   (  +   )   ,  mesh =     .Obviously, these frequency values are   = 0.441 Hz,   = 2.49 Hz, and   = 0.68 Hz, respectively, and the mesh frequency is  mesh = 34.9Hz.The errors    ,    ,    ,   ,   are set to be zero.The dynamic equation of the model is solved by the NEWMARK method.Figure5(a) shows the acceleration spectrum of the ring gear in longitudinal direction (without breakage).The spectrum is occupied by the mesh frequency and no . The ring gear is fixed and there are three planets rotating around the sun gear.The speed of the input part (carrier) is 2.77 rad/s (Ω  ) and the load of the output part (sun gear) is 200 N⋅m.The rotating frequency of the carrier (  ), sun gear (  ), planets (  ), and the mesh frequency ( mesh ) can be calculated by

Table 2 :
Parameters of the parallel-stage gear.