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 to much more 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.
1. 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–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.
2. 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 (s), carrier (c), ring (r), and three planets (p). 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 (x), longitudinal (y), and torsional degrees (u).
The dynamic model for a planetary gear system.
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 φpi (Figure 2) measured relatively to the rotating coordinate system with φp1=0.
A sun-planet mesh.
The global equation of motion for the system can be simplified as [18, 19]
(1)Mq¨+ΩcGq˙+[Kb+Ke(t)-Ωc2KΩ]q=T(t)+FW(t)+Fe(t),
where q represents the vector of the degrees of freedom
(2)q={xc,yc,uc,xr,yr,ur,xs,ys,us,lllllx1,y1,u1,…,xN,yN,uN}T.Ωc 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
(3)M=[Mc00000⋯00Mr0000⋯000Ms000⋯0000M100⋯00000M20⋯0⋮⋮⋮⋮⋮⋮⋮0000000⋯MN],Mj=[mj000mj000Ijrj2],j=c,r,s,1,…,NG=[Gc00000⋯00Gr0000⋯000Gs000⋯0000G100⋯00000G20⋯0⋮⋮⋮⋮⋮⋮⋮0000000⋯GN],Gj=[0-2mj02mj00000],j=c,r,s,1,…,N,
where mj and rj are the mass and the base radius of part j,Ij is the inertial moment of part j with respect to its rotational axis.
Kb is the bearing stiffness matrices and are expressed by
(4)lKb=[Kcb000⋯00Krb00⋯000Ksb0⋯00000⋯0⋮⋮⋮⋮⋮00000⋯0],Kjb=[kjx000kjy000kju],j=c,r,s,
where kjx,kjy,kju 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
(5)KΩ=[KΩc0000⋯00KΩr000⋯000KΩs00⋯0000KΩ10⋯00000KΩ2⋯0⋮⋮⋮⋮⋮⋮0000000KΩN],KΩj=[mj000mj0000],j=c,r,s,1,…,N.
The time-varying mesh stiffness matrix is represented by Ke(t). The matrices are expressed as
(6)Ke(t)=[∑Kc1i00Kc21Kc22⋯Kc2N0∑Kr1i0Kr21Kr22⋯Kr2N00∑Ks1iKs21Ks22⋯Ks2NKc21Kr21Ks21K10⋯0Kc22Kr22Ks220K2⋯0⋮⋮⋮⋮⋮⋮0Kc2NKr2NKs2N00⋯KN],
where
(7)Kr1i=krpi(t)[sin2φrpi-cosφrpicosαrsinφrpi-cosφrpicosφrcos2φrpicosφrpisinφrpicosφrpi1],Kc1i=kpi[10-sinφpi01cosφpi-sinφpicosφpi1],Ks1i=kspi(t)[sin2φspi-cosφspisinφspi-sinφspi-cosφspisinφspicos2φspicosφspi-sinφspicosφspi1],Kc2i=kpi[-cosφpisinφpi0-sinφpi-cosφpi00-10],Ks2i=kspi(t)[sinφspisinαssinφspicosαs-sinφspi-cosφspisinαs-cosφspicosαs-cosφspi-sinαs-cosαs1],Kr2i=krpi(t)[-sinφrpisinαrsinφrpicosαrsinφrpicosφrpisinαr-cosφrpicosαr-cosφrpisinαrcosαr-1],Ki=Kc3i+Kr3i+Ks3i,Kc3i=[kpi000kpi0000]Ks3i=kspi(t)[sin2αscosαssinαs-sinαscosαssinαscos2αs-cosαs-sinαs-cosαs1],Kr3i=krpi(t)[sin2αr-cosαrsinαr-sinαr-cosαrsinαrcos2αrcosαr-sinαrcosαr1],
where kjpi(t) (j=s,r) is the meshing stiffness between the part j and the ith planet. kpi is the bear supporting stiffness of the ith planet in the transverse and longitudinal direction. αj is the pressure angle of part j. φspi=φpi-αs and φrpi=φpi+αr, φpi is illustrated in a sun-planet mesh (Figure 2). i=1,2,…,N, N is the number of the planet gears.
T(t) represents the force matrix and can be expressed by
(8)T(t)=[0,0,Tcrc,0,0,Trrr,0,0,Tsrs,0,…,0]T,
where Tj(j=c,s,r) is the external torques applied on part j. rj is the base radius of part j.
FW(t) is the impact force caused by faults, the mathematical expression of which will be given in Section 3.
Fe(t) is the excitation force induced by manufacturing errors and can be written as
(9)Fe(t)=[0,Fc,Fs,F1,F2,…,FN],
where
(10)Fc=[∑i=1Nkrpi(t)eri(t)sin(φrpi),lllllllllll-∑i=1Nkrpi(t)eri(t)cos(φrpi),-∑i=1Nkrpi(t)eri(t)]T,Fs=[∑i=1Nkspi(t)esi(t)sin(φspi),lllllllllll-∑i=1Nkspi(t)esi(t)cos(φspi),-∑i=1Nkspi(t)esi(t)]T,Fn={∑i=1N-krpi(t)eri(t)sin(αr)+kspi(t)esi(t)sin(αs),llllllllll∑i=1Nkrpi(t)eri(t)sin(αr)-kspi(t)esi(t)cos(αs)llllllllll∑i=1Nkrpi(t)eri(t)-kspi(t)esi(t)}lllllllllllllllllllllllllllllllllllln=1,2,…,Neji(t) (j=s,r) is the errors between the part j and the ith planet and can be written as(11a)esi(t)=Epis(t)-Espi(t)+espi(t)leri(t)=Erpi(t)-Epir(t)+erpi(t),
where Ejpi(t) and Epij(t) (j=s,r) 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],
(11b)Ejpi(t)=Ejsin[fmeshzjt+εj-ψji],Epij(t)=Episin[fmeshzpt+εpi±αj],
where Ej and Epi are the amplitudes of these errors, fmesh is mesh frequency, zj and zp are tooth numbers of the component j and planet, εj and εpi are phase angles, ψji is relative position angle, and the arithmetic sign of the pressure angle αj is plus for j=r and minus for j=s.
In (11a), espi and erpi are gear transmission errors. Inalpolat and Kahraman [17] believed that those errors were amplitude and frequency modulated by Ejpi(t) and Epij(t). Therefore, another equation is able to express the gear transmission errors as given by Inalpolat and Kahraman [17],
(11c)espi(t)=Aspiespisin[Bspi(t)fmesht+zsφpi+ϕspi]erpi(t)=Arpierpisin[Brpi(t)fmesht+zrφpi+γsr+ϕrpi],
where
(11d)Aspi(t)=1+βssin[fmeshzst+ϕβsi]lllllllllllllll+βpisin[fmeshzpt+ϕβpi]Bspi(t)=1+β^ssin[fmeshzst+ϕ^β^si]llllllllllllll+β^pisin[fmeshzpt+ϕ^β^pi](11e)Arpi(t)=1+βrsin(fmeshzct+ϕβri)lllllllllllllll+ϑpisin[fmeshzpt+ϕϑpi+π]Brpi(t)=1+β^rsin[fmeshzct+ϕ^β^ri]lllllllllllllll+ϑ^risin[fmeshzpt+ϕ^ϑ^pi+π],where ϕjpi(j=s,r) is the phase angle of ejpi(t), γsr is the phase angle between espi(t) and erpi(t), βs, βpi, βr, ϑpi are amplitude coefficients, ϕβsi, ϕβpi, ϕβri, ϕϑpi are the phase angles of amplitude modulation, β^s, β^pi, β^r, ϑ^ri are amplitude coefficients, and ϕβri, ϕϑpi, ϕ^β^ri, ϕ^ϑ^pi are the phase angles of frequency modulation. The research results by Conry and Seireg [21] showed that the peak-to-peak amplitude of espi(t) might experience up to 4–6% variation due to the influence of Espi(t) and Epis(t) even with the modest error magnitudes.
3. 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 (Fw(t)).
3.1. 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 schematic diagram of a cantilever beam.
The total deformation of the point B is expressed by
(12)δBi=uBsinγi+vBcosγi+wBsinγi+δf+δc,
where uB, vB, wB, δf, δc represent bending, compressive, shear, fillet-foundation, and contact deformations, respectively, which can be calculated by the method used in [12], γi (see Figure 3) defines the angle between the tooth normal force and the tooth centerline. The stiffness of one gear tooth at point B is calculated as ki=F/δBi. The total stiffness of one pair of meshing tooth is obtained by kpair=(k1×k2)/(k1+k2) according to the serial relation. The stiffness when double-teeth meshing is k=kpair1+kpair2 according to the parallel relation.
The parameters of the studied planetary gear system are listed in Table 1. Young’s modulus E=2.06×1011 Pa, Poisson’s ratio υ=0.3. Figure 4 shows the mesh stiffness of the sun gear and ring gear with planet 1 (ksp1 and krp1). The planet 1 with and without tooth breakage is considered and the mesh stiffness is seriously affected by tooth breakage.
Parameters of the planetary gear system.
Parameters
Sun gear
Ring
Carrier
Planets
Number of tooth
17
79
—
31
Module (mm)
3
3
—
3
Face width (mm)
35
35
—
35
Mass (kg)
1.75
3.552
7
0.979
Mass moment of inertia (kg/m2)
1
3
4.56
0.544
Bearing support stiffness (N/m)
KP=Ksx=Ksy=Krx=Kry=Kcx=Kcy=2×108
Torsional stiffness (N/m)
Kru=109;Ksu=Kcu=2×107
Pressure angle (°)
20°
20°
—
20°
Helix angle (°)
0°
0°
—
0°
The mesh stiffness of planet 1 with and without tooth breakage: (a) planet-sun gear pair, (b) planet-ring gear pair.
3.2. 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 FWpf(t) generated by the tooth breakage of the planet gear can be expressed as (planet 1 is assumed to be the fault gear)
(13)FWs1(t)=[0,0,Fws(t),0,Fw1(t),0,0]TFWr1(t)=[0,Fwr(t),0,0,Fw1(t),0,0]T,
where
(14)Fws(t)=[Fs1(t)sin(φsp1),-Fs1(t)cos(φsp1),-Fs1(t)]TFwr(t)=[Fr1(t)sin(φrp1),-Fr1(t)cos(φrp1),Fr1(t)]TFw1(t)=[-Fr1(t)sin(αr)+Fs1(t)sin(αs),Fr1(t)sin(αr)lllllllllllllll+Fs1(t)cos(αs),Fr1(t)-Fs1(t)]T,
where FWj1(t) represents the impact force vector of gear j (j=s,r), Fj1(t) (j=s,r) denotes the impact force and can be expressed as
(15)Fj1(t)={Cj1(t)tIbj1<t<tIfj10else,
where tIbj1 and tIfj1 represent the beginning and finishing time of the impact force of gear j. The relationship between Cr(t) and Cs(t) is expressed by Cr(t)=H(t)Cs(t). Because the impact acceleration will be generated, H(t)>1 when the fault tooth meshes with ring gear and H(t)<1 when the fault tooth meshes with sun gear.
The additional impact force FWrf(t) generated by the tooth breakage of the ring gear could be expressed by
(16)FWr1(t)=[0,Fwr1(t),Fws1(t),0,Fw1(t),0,0]TFWr2(t)=[0,Fwr2(t),Fws2(t),0,0,Fw2(t),0]TFWr3(t)=[0,Fwr3(t),Fws3(t),0,0,0,Fw3(t)]T,
where
(17)Fwri(t)=[Fri(t)sin(φrpi),-Fri(t)cos(φrpi),-Fri(t)]TFwsi(t)=[Fsi(t)sin(φspi),-Fsi(t)cos(φspi),-Fsi(t)]TFwi(t)=[-Fri(t)sin(αr)+Fsi(t)sin(αs),llllllFri(t)sin(αr)+Fsi(t)cos(αs),Fri(t)-Fsi(t)]T,
where FWri(t) represents the impact force matrix when planet i meshes with ring gear, Fsi(t) and Fri(t) are the pulsating impact forces, the mathematical expressions of which are similar to (15) and could be expressed by
(18)Fji(t)={Cji(t)tIbji<t<tIfji0else.
Similarly, there are several coefficients Hji(t) used to describe the relation among Cji(t). The starting time of Fsi(t) and Fri(t) varies with i.
The impact force generated by the tooth breakage of the sun gear is FWsf(t), the mathematical expression of which is the same as that of FWrf(t). However, the values for beginning and finishing times of Fji(t) 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 (Tm). 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 tIfji-tIbji=Tm/10 and Cji(t)=2F¯jpi=const.
4. The Dynamic Response of Planetary Gear System without Errors
Physical parameters of the studied planetary gear system are listed in Table 1. 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 (Ωc) and the load of the output part (sun gear) is 200 N·m. The rotating frequency of the carrier (fc), sun gear (fs), planets (fp), and the mesh frequency (fmesh) can be calculated by
(19)fc=Ωc2π,fs=zs+zrzsfc,fp=fszs-fc(zs+zp)zp,fmesh=fczr.
Obviously, these frequency values are fc=0.441 Hz, fs=2.49 Hz, and fp=0.68 Hz, respectively, and the mesh frequency is fmesh=34.9 Hz.
The errors espi, erpi, Epi, Es, Er are set to be zero. The dynamic equation of the model is solved by the NEWMARK method. Figure 5(a) shows the acceleration spectrum of the ring gear in longitudinal direction (without breakage). The spectrum is occupied by the mesh frequency and no 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 fmesh±n(fp+fc), fmesh±n(fs-fc) and fmesh±nfc (n=1,2,…) successively. Therefore, the mesh frequency fmesh 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 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.
5. 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
(20)a(t)=∑i=1Nwi(t)(C¯ikspi(t)δspi(t)+D¯ikrpi(t)δrpi(t)),
where C¯i and D¯i are constants, δspi(t) and δrpi(t) are defined as
(21)δspi(t)=yscosφspi-xssinφspi-xisinαs+us+ui-esi(t),δrpi(t)=yrcosφrpi-xrsinφrpi-xisinαr+ur-ui+eri(t)xj, yj, uj (j=s,r,i) can be obtained by the dynamic model.
wi(t) is a weighting function for a planet gear positioned at the angle of φpi and wi(t) can be defined as
(22)wi(t)=W(t-φpifc)Ui(t),
where W(t) is a Hanning function used to represent the influence of rolling planets on fixed accelerometer (scrolling effect) and Ui(t) is the unit step function. Both of them are given as follows:
(23)W(t)=12-12[cos(2πNfct)]Ui(t)=∑n=1∞{u[t-((n-1)N+i-1N)Tc]lllllllllll-u[t-((n-1)N+iN)Tc]},
where u(t-a)=1 for t>a and u(t-a)=0 for t<a, N is the number of planets. In this study, N is equal to 3 and i is 1, 2, and 3.
Run-out errors of the planets, sun gear, and ring gear are taken into account and espi=0.6μm, erpi=0.3μm, Epi=10μm, and Es=Er=15μm. According to Conry’s study [21], modulation coefficients are set as βs=0.033, ϕβsi=2(i-1)π/3, βpi=0.034, ϕβpi=2(i-2)π/3, ϑpi=0.033, ϕϑpi=2(i-2)π/3, βr=0.035, ϕβri=2(i-1)π/3 and all other modulation coefficients are set to be zero. The dynamic equation of the model is solved to get the δspi(t) and δrpi(t), then acceleration is calculated by (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.
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 fmesh), (c) sun gear with and without tooth breakages (the right side of fmesh), (d) ring gear with and without tooth breakages.
Figure 6(b) shows the acceleration spectrum of the system with and without tooth breakage on planet 1. Since the mesh frequency fmesh is modulated by the frequencies of weighting function (which characteristic frequencies are mfc) and the frequencies of n(fp+fc) (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 fmesh±mfc±n(fp+fc) (m,n=1,2,…). Therefore, the fault characteristic frequencies of the planet gear are fpf=fmesh±mfc±n(fp+fc) and can be simplified as fmesh±mfc±nfp (m,n=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 fmesh±nfc±mfs (m,n=1,2,…) in Figure 6(c). The main reason is that the mesh frequency fmesh is modulated by the frequencies of weighting function (mfc) and the frequencies of n(fs-fc) (caused by the influence of tooth breakage on mesh stiffness and additional impact force). Therefore, the fault characteristic frequencies of the sun gear are fmesh±mfc±n(fs-fc) and can be as fsf=fmesh±nfc±mfs (m,n=1,2,…).
Similarly, because the mesh frequency fmesh is modulated by the frequencies of the weighting function mfc and the frequencies of nfc (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 fmesh±nfc±mfc (m,n=1,2,…), as shown in Figure 6(d). Therefore, the fault characteristic frequencies of the ring gear are frf=fmesh±nfc±mfc and can be simplified as frf=fmesh±nfc (n=1,2,…).
6. 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.
Parameters of the parallel-stage gear.
Level
Stage1
Stage2
Gear
Input
Output
Input
Output
Number of tooth
73
21
66
23
The test bench of the planetary gear system.
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) and 8(c). Vibration signals are collected for a normal mode (without breakage) and the two fault modes. The sampling frequency is 16384 Hz.
The acceleration sensor and breakage tooth: (a) acceleration sensor position, (b) planet gear with tooth breakage, and (c) ring gear with tooth breakage.
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.
The acceleration spectrum of the normal system.
6.1. Planet Gear with Tooth Breakage
Figure 10 shows the acceleration spectrum of the system with and without tooth breakage on planet.
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.
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 fmesh±mfc±nfp (m,n=1,2,…), except for the frequencies of fmesh+2fc and fmesh+fc+3fp.
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 fmesh±mfc±nfp (m,n=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 fmesh±mfc±nfp (m,n=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 fpf=fmesh±mfc±nfp (m,n=1,2,…).
6.2. Ring Gear with Breakage
Figure 11 shows the acceleration spectrum of the system with and without tooth breakage on the ring gear.
The spectrum of the system with and without tooth breakage on ring gear at different loads and speeds: (a) load = 200 N·m, input speed = 2.77 rad/s, (b) load = 229 N·m, input speed = 2.77 rad/s, and (c) load = 229 N·m, input speed = 2.24 rad/s.
In Figure 11, the load and the speed of the system are 200 N·m and 2.77 rad/s, respectively; therefore, the mesh frequency is 34.9 Hz, similar with the method used in Section 6.1. The amplitudes of the acceleration spectrum are compared between the normal system and the tooth breakage system. The fault characteristic frequencies of fmesh±nfc (n=1,2,…) are further confirmed by the experiment results.
The change law of amplitudes with loads and speeds of the normal and the tooth breakage systems is same as the law described in Section 6.1. The conclusion can be obtained through Figures 11(a)–11(c).
7. Conclusions
Detection and localization of tooth fault on wind turbine planetary gear system are investigated numerically and experimentally in this paper. Through detailed comparisons on both numerical and experimental results, the fault characteristic frequencies can be obtained after considering the gear manufacturing errors.
The sidebands of the planetary gear system become much more complex when manufacturing errors are taken into consideration.
The fault characteristic frequencies of planet gears, sun gear and ring gear with tooth breakage are fpf=fmesh±mfc±nfp, fsf=fmesh±nfc±mfs and frf=fmesh±nfc, respectively.
As the fault characteristic frequencies for the sun, planet, and ring gears are different, so the tooth breakage fault could also be located by monitoring the values of these frequencies.
The above results are useful for the detection and localization of tooth faults in actual planetary gear system of wind turbines.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The research work described in this paper was supported by Natural Science Foundation of China under Grant 51335006, Tsinghua University Initiative Scientific Research Program under Grant 2011Z08137, and the State Key Laboratory of Tribology Initiative Research Program under Grant SKLT11A02.
McNamesJ.Fourier series analysis of epicyclic gearbox vibration2002124115015210.1115/1.14037352-s2.0-0036400303KusiakA.LiW.The prediction and diagnosis of wind turbine faults2011361162310.1016/j.renene.2010.05.0142-s2.0-79751491467HatchC.Improved wind turbine condition monitoring using acceleration envelopingOrbit, pp. 58–61, 2004WalfordC.Wind turbine reliability: Understanding and minimizing wind turbine operation and maintenance costsMarch 2006SAND2006-1100Sandia National LaboratoriesLeiY.KongD.LinJ.ZuoM. J.Fault detection of planetary gearboxes using new diagnostic parameters201223505560510.1088/0957-0233/23/5/0556052-s2.0-84859884263LinJ.ZuoM. J.Gearbox fault diagnosis using adaptive wavelet filter20031761259126910.1006/mssp.2002.15072-s2.0-0041880035LiuW. Y.ZhangW. H.HanJ. G.WangG. F.A new wind turbine fault diagnosis method based on the local mean decomposition20124841141510.1016/j.renene.2012.05.0182-s2.0-84862998489McFaddenP. D.SmithJ. D.An explanation for the asymmetry of the modulation sidebands about the tooth meshing frequency in epicyclic gear vibration1985199165702-s2.0-0021784348McFaddenP. D.A technique for calculating the time domain averages of the vibration of the individual planet gears and the sun gear in an epicyclic gearbox1991144116317210.1016/0022-460X(91)90739-72-s2.0-0026419440FengZ.ZuoM. J.Vibration signal models for fault diagnosis of planetary gearboxes2012331224919493910.1016/j.jsv.2012.05.0392-s2.0-84863303418ChaariF.FakhfakhT.HaddarM.Dynamic analysis of a planetary gear failure caused by tooth pitting and cracking200662737810.1361/154770206X993432-s2.0-33646769250ChaariF.FakhfakhT.HaddarM.Analytical investigation on the effect of gear teeth faults on the dynamic response of a planetary gear set200637891510.1260/0957456067786008572-s2.0-33749597764ChengZ.HuN.GuF.QinG.Pitting damage levels estimation for planetary gear sets based on model simulation and grey relational analysis20113534034172-s2.0-84855772767ChengZ.HuN. Q.ZuoM. J.FanB.Crack level estimation approach for planetary gear sets based on simulation signal and GRA2012364012076JiaS.HowardI.Comparison of localised spalling and crack damage from dynamic modelling of spur gear vibrations200620233234910.1016/j.ymssp.2005.02.0092-s2.0-27744520895LigataH.KahramanA.SinghA.An experimental study of the influence of manufacturing errors on the planetary gear stresses and planet load sharing2008130404170110.1115/1.28851942-s2.0-45149089914InalpolatM.KahramanA.A dynamic model to predict modulation sidebands of a planetary gear set having manufacturing errors2010329437139310.1016/j.jsv.2009.09.0222-s2.0-71849114668ChaariF.FakhfakhT.HbaiebR.LouatiJ.HaddarM.Influence of manufacturing errors on the dynamic behavior of planetary gears2006277-873874610.1007/s00170-004-2240-22-s2.0-29344452911LinJ.ParkerR. G.Analytical characterization of the unique properties of planetary gear free vibration1999121331632110.1115/1.28939822-s2.0-0032642480BarszczT.RandallR. B.Application of spectral kurtosis for detection of a tooth crack in the planetary gear of a wind turbine20092341352136510.1016/j.ymssp.2008.07.0192-s2.0-58949090476ConryT. F.SeiregA.A mathematical programming technique for the evaluation of load distribution and optimal modifications for gear systems1973954111511222-s2.0-001549075910.1115/1.3438259LinT. J.OuH.LiR. F.A finite element method for 3D static and dynamic contact/impact analysis of gear drives20071969–121716172810.1016/j.cma.2006.09.0142-s2.0-33845664086YaoW. X.WeiR. Z.Research on the mesh impact of spur gear199045761WangL. X.JunL. L.LeS. Z.Impact incentive research of broken teeth gear2012334264266