^{1}School of Mechanical Engineering, Baoji University of Arts and Sciences, Baoji, 721016, China
Journal of Vibroengineering, Vol. 21, Issue 6, 2019, p. 17381750.
https://doi.org/10.21595/jve.2019.20471
Received 24 December 2018; received in revised form 15 March 2019; accepted 25 March 2019; published 30 September 2019
JVE Conferences
With the increasing service life of a gear transmission system, the gear teeth become constantly worn, and the gear clearance increases. The increase in the clearance will produce a series of nonlinear changes that change the stability of the system and can even result in loss of stability. In this paper, dimensionless dynamic equations of a multistage gear transmission system that contains a twostage fixedaxis gear and a onestage planetary gear were established. Planetary wear fault was simulated by changes in the gear clearance. System bifurcation diagrams with an increase in the clearance were studied. The frequency characteristics of planetary gears under different excitation frequencies and different degrees of wear were studied. The influence of planetary gear wear on the fixedaxis was discovered. The vibration mechanism and a fault diagnostic method in a multistage gear transmission system were obtained.
 By studying the timedomain, frequencydomain, phase diagram and Poincaré section of each motion state, the transition processes and frequency characteristics of each motion state were analyzed.
 Frequency characteristics of the planetary gears under different excitation frequencies and different degrees of wear were obtained.
 The influence of the planetary gear wear on the fixedaxis gear were obtained.
Keywords: planetary gear, nonlinear dynamical, wear, frequency characteristics, fault diagnosis.
1. Introduction
With the increasing service life of a gear transmission system, the gear teeth become continuously worn, and the gear clearance increases. Clearance is the main reason for the collision motion inside the system. It will make the motion state of the system change or even create instability. Therefore, it is very important to study the wear mechanism and the change in the system nonlinear state, as these attributes are of great significance for the design, monitoring and diagnosis of the system. When both planetary gear and fixedaxis gear exist in the system, the vibration signal of planetary gear is weaker than that of fixedaxis gear, so the wear fault of planetary gear in multistage gear transmission system is more difficult to identify.
For the study of wear failure mechanism, based on the different classification of wear, a large number of studies have established wear calculation models based on adhesion theory [13], fatigue theory [4, 5], energy wear theory [6] and material wear rate [7]. It has been studied from various angles such as tooth profile [3], gap [8], dynamic response [9], and experiment [10]. However, there are many factors affecting gear wear, and the wear mechanism is closely related to the actual working condition, so it is difficult to establish a wear prediction model completely consistent with the actual transmission system. At the same time, under the condition of multifactor coupling, wear failure has strong nonlinear characteristics. Therefore, it is particularly important to study the nonlinear characteristics of wear. Shen [11] focused on the nonlinear dynamics of a spur gear pair with slight wear fault, where the backlash, timevarying stiffness and wear fault were all included. Ding [9] and Liu [12] introduced the dynamic wear coefficient proposed by Priest [13] into the calculation of spur gear wear. Based on the multifactor excitation nonlinear dynamic model, the coupling effect between tooth wear and system dynamic response was studied. These literature all aim at the single pair gear model, and neglect the coupling characteristics between multistage gears. However, in the multigear transmission system, the fault of planetary gear is more difficult to be identified due to its weak signal. Therefore, it is necessary to further study the nonlinear dynamic characteristics of the multistage gear transmission system with planetary gear wear fault.
For the fault diagnosis of wear failure, the most common test methods are ferro spectral analysis [14, 15] and timefrequency characteristic analysis of vibration signal [16, 17]. In order to facilitate the comparison of theoretical and experimental results, the timefrequency characteristic analysis method will be used in this paper. Through the mechanism analysis of gear wear, the transmission characteristics of gear wear are identified. At present, some scholars have started relevant researches. He [18] proposed a novel nonlinear timefrequency feature based on a timefrequency manifold (TFM) technique. The effectiveness and the merits of the proposed TFM feature were confirmed by case study on gear wear diagnosis, bearing defect identification and defect severity. Liu [12] proposed a dynamic wear prediction methodology to investigate the coupling effects between surface wear and dynamics of spur gear systems. This paper will also use the nonlinear timefrequency characteristics of wear fault to present a diagnostic method of wear fault.
In this paper, a nonlinear dynamic model of multistage gear transmission system that contains a twostage fixedaxis gear and an onestage planetary gear was established. The simulation method was used to analyze the bifurcation characteristics of the system under different clearances to find out the excitation frequency range of wear and the change in the system motion state. In order to understand the frequency characteristics of the system. The transition process and frequency characteristics of each motion state under different excitation frequencies were discussed. The vibration mechanism and a diagnostic method of wear were obtained.
2. Torsional dynamic model of a multistage gear transmission system
The system studied in this paper is a test rig with a multistage gear transmission system that contains a twostage, fixedaxis gear and a onestage planetary gear, where spur gears 1 and 2 compose the 1st  stage fixedaxis gear for the input, spur gears 3 and 4 compose the 2nd  stage fixedaxis gear, and the planet carrier is for the output. The torsional dynamic model is established using the lumped mass method (Fig. 1). The model does not consider the transverse vibration displacement of the gears. The gear parameters are simulated with a spring and a damper.
Fig. 1. Torsional dynamic model of the multistage gear transmission system
In Fig. 1, ${\theta}_{s}$, ${\theta}_{c}$, ${\theta}_{pn}$, ${\theta}_{1}$, ${\theta}_{2}$, ${\theta}_{3}$, and ${\theta}_{4}$ represent the angular displacements of the sun gear, planet carrier, planetary gear $n$ ($n=$ 1, 2, 3, 4), and spur gears 1, 2, 3, and 4, respectively. Throughout this paper, the subscripts $s$, $c$, $pn$, $r$, and 1, 2, 3, and 4 denote the sun gear, planet carrier, planetary gear, ring gear, and spur gears 1, 2, 3, and 4, respectively. The quantities ${r}_{s}$, ${r}_{c}$, ${r}_{pn}$, ${r}_{1}$, ${r}_{2}$, ${r}_{3}$, and ${r}_{4}$ are the base circle radius values of the gears. The quantities ${K}_{spn}\left(t\right)$, ${K}_{rpn}\left(t\right)$, ${K}_{1}\left(t\right)$, and ${K}_{2}\left(t\right)$ denote the meshing stiffness of the sun gear with planetary gear $n$, the ring gear with planetary gear $n$, the 1ststage fixed gear, and the 2ndstage fixed gear, respectively. The quantities ${C}_{spn}$, ${C}_{rpn}$, ${C}_{1}$, and ${C}_{2}$ denote the damping of the sun gear with planetary gear $n$, the ring gear with planetary gear $n$, the 1ststage fixed gear, and the 2ndstage fixed gear, respectively. ${T}_{in}$ is the input, and ${T}_{out}$ is the output.
2.1. Motion differential equations of the system
By using the Lagrangian equation, the clearance, timevarying meshing stiffness, and composite error are considered to establish the motion differential equations of the system. After that, the motion differential equation is processed in a dimensionless way. This process has been deduced and detailed in my paper [19]. Only the final system dimensionless motion differential equations are listed here:
$=\frac{{T}_{in}{r}_{1}}{{J}_{1}{w}_{h}^{2}{b}_{1}}+\frac{{e}_{a1}}{{b}_{1}}{\mathrm{\Omega}}_{1}^{2}\mathrm{s}\mathrm{i}\mathrm{n}\left({\mathrm{\Omega}}_{1}\tau +{\phi}_{1}\right),$
${\ddot{\stackrel{}{x}}}_{2}\frac{{C}_{1}}{{m}_{e1}{w}_{h}}{\dot{\stackrel{}{x}}}_{1}\frac{{K}_{1}\left(\tau \right)}{{m}_{e2}{w}_{h}^{2}}f\left({\stackrel{}{x}}_{1}\right)+\frac{{C}_{2}}{{m}_{e3}{w}_{h}^{}}{\dot{\stackrel{}{x}}}_{2}+\frac{{K}_{2}\left(\tau \right)}{{m}_{e3}{w}_{h}^{2}}f\left({\stackrel{}{x}}_{2}\right)\sum _{n=1}^{4}\frac{{C}_{spn}}{{m}_{4s}{w}_{h}^{}}{\dot{\stackrel{}{x}}}_{spn}\sum _{n=1}^{4}\frac{{K}_{spn}\left(\tau \right)}{{m}_{4s}{w}_{h}^{2}}f\left({\stackrel{}{x}}_{spn}\right)=\frac{{e}_{a2}}{{b}_{1}}{\mathrm{\Omega}}_{1}^{2}\mathrm{s}\mathrm{i}\mathrm{n}\left({\mathrm{\Omega}}_{2}\tau +{\phi}_{2}\right),$
${\ddot{\stackrel{}{x}}}_{spn}+\frac{1}{{m}_{4s}{w}_{h}}\sum _{n=1}^{4}{C}_{spn}{\dot{\stackrel{}{x}}}_{spn}+\frac{1}{{m}_{pn}{w}_{h}}\sum _{n=1}^{4}{C}_{spn}{\dot{\stackrel{}{x}}}_{spn}+\frac{1}{{m}_{4s}{w}_{h}^{}}{C}_{spn}{\dot{\stackrel{}{x}}}_{spn}$
$+\frac{1}{{m}_{4s}{w}_{h}^{2}}\sum _{n=1}^{4}{K}_{spn}\left(\tau \right)f\left({\stackrel{}{x}}_{spn}\right)+\frac{1}{{m}_{c}{w}_{h}^{2}}\sum _{n=1}^{4}{K}_{spn}\left(\tau \right)f\left({\stackrel{}{x}}_{spn}\right)+\frac{1}{{m}_{pn}{w}_{h}^{2}}{K}_{spn}\left(\tau \right)f\left({\stackrel{}{x}}_{spn}\right)$
$\frac{1}{{m}_{pn}{w}_{h}}{C}_{rpn}{\dot{\stackrel{}{x}}}_{rpn}+\frac{1}{{m}_{c}{w}_{h}}\sum _{n=1}^{4}{C}_{rpn}{\dot{\stackrel{}{x}}}_{rpn}\frac{1}{{m}_{pn}{w}_{h}^{2}}{K}_{rpn}\left(\tau \right)f\left({\stackrel{}{x}}_{rpn}\right)$
$+\frac{1}{{m}_{c}{w}_{h}^{2}}\sum _{n=1}^{4}{K}_{rpn}\left(\tau \right)f\left({\stackrel{}{x}}_{rpn}\right)\frac{{r}_{4}{C}_{2}}{{m}_{4s}{r}_{s}{w}_{h}}{\dot{\stackrel{}{x}}}_{2}\frac{{r}_{4}{K}_{2}\left(\tau \right)}{{m}_{4s}{r}_{s}{w}_{h}^{2}}f\left({\stackrel{}{x}}_{2}\right)$
$=\frac{{r}_{c}{T}_{out}}{{J}_{ce}{w}_{h}^{2}{b}_{1}}+\frac{{e}_{aspn}}{{b}_{1}}{\mathrm{\Omega}}_{spn}^{2}\mathrm{s}\mathrm{i}\mathrm{n}\left({\mathrm{\Omega}}_{spn}^{}\tau +{\phi}_{spn}\right),$
${\ddot{\stackrel{}{x}}}_{rpn}\frac{1}{{m}_{pn}{w}_{h}}{C}_{spn}{\dot{\stackrel{}{x}}}_{spn}+\frac{1}{{m}_{c}{w}_{h}}\sum _{n=1}^{4}{C}_{spn}{\dot{\stackrel{}{x}}}_{spn}\frac{1}{{m}_{pn}{w}_{h}^{2}}{K}_{spn}\left(\tau \right)\left({\stackrel{}{x}}_{spn}\right)$
$+\frac{1}{{m}_{c}{w}_{h}^{2}}\sum _{n=1}^{4}{K}_{spn}\left(\tau \right)f\left({\stackrel{}{x}}_{spn}\right)+\frac{1}{{m}_{pn}{w}_{h}^{}}{C}_{rpn}{\dot{\stackrel{}{x}}}_{rpn}\frac{1}{{m}_{c}{w}_{h}^{}}\sum _{n=1}^{4}{C}_{rpn}{\dot{\stackrel{}{x}}}_{rpn}$
$+\frac{1}{{m}_{c}{w}_{h}^{2}}{K}_{rpn}\left(\tau \right)f\left({\stackrel{}{x}}_{rpn}\right)\frac{1}{{m}_{c}{w}_{h}^{2}}\sum _{n=1}^{4}{K}_{rpn}\left(\tau \right)f\left({\stackrel{}{x}}_{rpn}\right)$
$=\frac{{r}_{c}{T}_{out}}{{J}_{ce}{w}_{h}^{2}{b}_{1}}\frac{{e}_{arpn}}{{b}_{1}}{\mathrm{\Omega}}_{rpn}^{2}\mathrm{s}\mathrm{i}\mathrm{n}\left({\mathrm{\Omega}}_{rpn}^{}\tau +{\phi}_{rpn}\right).$
2.2. Wear fault
The gear tooth surface will wear out with an increase in the running time. The clearance between the teeth will change when wear occurs. When gear teeth are evenly worn, the tooth clearance will be larger than normal. To model wear failure, Flodin [20] simulated mild wear spur gears. Park [21] presented an approximate method to predict the surface wear of hypoid gears using surface interpolation. In this paper, the clearance was controlled by changing the size of the dimensionless composite error amplitude ${\overline{e}}_{ai}$ [22]. The bifurcation diagram of the multistage gear transmission system with the change in excitation frequency was calculated in the normal state. According to the clearance calculation formula given by the Chinese national standard GB 236390, the dimensionless clearance interval of the meshing gear in the normal state is [0.1, 1].
3. Nonlinear dynamic behavior analysis of the system with increased planetary gear wear
3.1. System bifurcation diagrams with increased planetary gear wear
By increasing the wear degree of the planetary gear 1, the influence of planetary gear wear on the excitation frequency range of the multistage gear transmission system was investigated. The minimum clearance of the multistage gear transmission system was chosen as the initial state, that is, ${\overline{e}}_{ai}=$ 0.1. At this point, the system is in a periodic motion state, and it is easy to observe the changes caused by the planetary gear wear. The variablestep RungeKutta method was used to solve the nonlinear differential Eq. (1). The bifurcation diagrams were investigated on the meshing points of the 1st  stage fixedaxis gear, and the planetary gear with the ring gear, when the fixedaxis gear clearance was constant (${\stackrel{}{e}}_{a1}={\stackrel{}{e}}_{a2}=$ 0.1) and the planetary clearance increased gradually (${\stackrel{}{e}}_{aspn}={\stackrel{}{e}}_{arpn}=$ 0.11), as shown in Fig. 2. The gear parameters of the multistage gear transmission system are shown in Table 1 and 2. The pressure angle $\alpha =$ 20°, ${T}_{in}=$ 6.5 N·m, ${T}_{out}=$ 8.5 N·m, ${b}_{i}=$ 5 μm, and the damping ratio of meshing pairs ${\xi}_{i}=$ 0.07. These parameters were derived from the multistage gear transmission system test rig. We set the sampling interval as ${\mathrm{\Omega}}_{1}=$0.01:0.01:5. Each ${\mathrm{\Omega}}_{1}$ calculates 80 rotating periods. The number of samples for each rotating period was 256 points. The bifurcation diagram started at the 25th period. The color in Fig. 2 represents 55 rotating periods. The initial velocity and displacement were set to 0. The calculation time of bifurcation diagram was 2820 s. The CPU computing time was obtained using a PC with Windows 7 operating system (Intel (E3) 3.2 GHz CPU) and 16 GB memory, and the computations were performed using MATLAB 2014a.
Table 1. Gear parameters
Gear

Number of teeth
$z$

${r}_{i}$
/ mm

Mass / g

Rotational inertia ${J}_{i}$ / (g∙m^{2})

Face width / mm

Modulus
m/mm

1

29

20.4

125

0.05

30

1.5

2

100

70.5

1224.5

6

30

1.5

3

36

25.3

224

0.14

30

1.5

4

90

63.4

1111

4

20

1.5

s

28

13

41

0.007

20

1.0

$pn$

36

17

34.6

0.01

20

1.0

$c$

–

30

848.7

0.76

20

1.0

$r$

100

47

–

–

20

1.0

Table 2. Calculation parameters
Parameters of calculation

Value

Gear clearance $b$ (μm)

5

comprehensive meshing error amplitude ${e}_{a}$ (μm)

2

Meshing pair damping ratio $\xi $

0.07

Gear contact ratio

1.68

As can be seen from Fig. 2(b), the planetary gear is in singleperiod motion state when the excitation frequency is small. When ${\mathrm{\Omega}}_{1}=$ 1, the system changes from singleperiod motion to multiperiod motion. Until ${\mathrm{\Omega}}_{1}=$ 3.5, it changes from multiperiod motion to quasiperiodic motion with large oscillations. The 1st stage fixedaxis gear is also in singleperiod motion when the excitation frequency is small (Fig. 2(a)). Affected by the synchronism of the planetary gears, the motion of the 1st stage fixedaxis gear changes from singleperiod to multiperiod, when ${\mathrm{\Omega}}_{1}=$ 1, and when ${\mathrm{\Omega}}_{1}=$ 3.5, it is also in quasiperiodic motion with large oscillations. As the wear increases (Fig. 2(d)), the excitation frequency of the quasiperiodic motion of the planetary gear advances to ${\mathrm{\Omega}}_{1}=$ 2.5, and the amplitude of the quasiperiodic motion increases. Fixedaxis gear has not been affected due to the small increase in clearance (Fig. 2(c)).
Fig. 2. Bifurcation diagram of the system with the change of excitation frequency under different gaps: a) and b) ${\stackrel{}{e}}_{aspn}={\stackrel{}{e}}_{arpn}=$ 0.1, c) and d) ${\stackrel{}{e}}_{aspn}={\stackrel{}{e}}_{arpn}=$ 0.2, e) and f) ${\stackrel{}{e}}_{aspn}={\stackrel{}{e}}_{arpn}=$ 0.4, g) and h) ${\stackrel{}{e}}_{aspn}={\stackrel{}{e}}_{arpn}=$ 0.8
a)
b)
c)
d)
e)
f)
g)
h)
As the wear continues to increase (Fig. 2(f)), the excitation frequency of planetary gear advances to ${\mathrm{\Omega}}_{1}=$ 2.1 when entering quasiperiodic motion, and quasiperiodic motion amplitude also increases. At this time, the planetary gear vibration is transmitted to the fixedaxis, so that the amplitude of the 1st stage fixedaxis gear is increased in the interval of ${\mathrm{\Omega}}_{1}$∈ [3, 4] (Fig. 2(e)). As the wear continues to increase (Fig. 2(h)), the excitation frequency of planetary gear advances to ${\mathrm{\Omega}}_{1}=$ 1.7 when entering the quasiperiodic motion, and the amplitude of quasiperiodic motion increases significantly with $\mathrm{\Omega}1\ge $4. The influence of synchronization of the planetary gear on the fixedaxis gear is more serious. The amplitude of the interval of ${\mathrm{\Omega}}_{1}$∈ [3, 4] increases continuously (Fig. 2(g)).
It can be seen from Fig. 2 that the planetary gear wear mainly affects the interval of ${\mathrm{\Omega}}_{1}\ge $ 2. In another article [23] of the author, it is found that when ${\mathrm{\Omega}}_{1}\le $ 2, it is mainly affected by the fixedaxis wear. The wear of the planetary gear causes the excitation frequency of the quasiperiodic motion of the planetary gear advances and increases the amplitude of the vibration. At this time, due to the synchronization, the changes of the planetary gear are transmitted to the fixedaxis gear, so that the amplitude of the quasiperiodic motion of the 1st stage fixedaxis gear is obviously increased in the interval of ${\mathrm{\Omega}}_{1}$∈ [3, 4].
3.2. Frequency characteristics of the planetary gears under different excitation frequencies
In order to understand the frequency characteristics of the system under different excitation frequencies and different degrees of wear, more detailed analysis of time domain, frequency domain, phase diagram and Poincaré section of each motion state are required. The large clearance state was (Fig. 2(f)) selected as the study object. At this time, the characteristics of the planetary gear and the fixedaxis gear are obvious. The dimensionless characteristic frequencies of the gears in the multistage gear transmission system are shown in Table 3. When ${\stackrel{}{e}}_{aspn}={\stackrel{}{e}}_{arpn}=$ 0.4, the dynamic characteristics of the planetary gear at different excitation frequencies were calculated, and Figs. 38 were obtained.
Table 3. Dimensionless characteristic frequencies of multistage gear transmission system (Hz)
Characteristic frequency

Dimensional frequency

Meshing frequency of 1ststage fixedaxis gear ${f}_{1}$

1

Meshing frequency of 2ndstage fixedaxis gear ${f}_{2}$

0.3599

Meshing frequency of planetary gear ${\mathrm{f}}_{3}$

0.0877

Fig. 3. Dynamic characteristics of the planetary gear ${\mathrm{\Omega}}_{1}=$ 1
Fig. 4. Dynamic characteristics of the planetary gear ${\mathrm{\Omega}}_{1}=$ 1.7
It can be seen from Fig. 3 that the system is in a singleperiod motion when the excitation frequency is small, and the Poincaré section is a point group. Due to the sun gear is connected to the 2nd stage fixedaxis gear, the planetary gear vibration signal is weak compared with the 2nd stage fixedaxis gear. The main peaks of the vibration spectrum are ${f}_{2}$ and ${f}_{2}\pm n{f}_{3}$. As the excitation frequency increases, the system leaves the singleperiod motion and enters a quasiperiodic motion (Fig. 5) through chaotic when the interval of ${\mathrm{\Omega}}_{1}$∈ [1.7, 2.1] (Fig. 4). Planetary gear vibration is intensified, and the amplitude of planetary gear meshing frequency becomes the main peak. After the quasiperiodic motion is stable, the amplitude of the planetary gear meshing frequency far exceeds that of the 2nd stage fixedaxis meshing frequency, and the main peaks at this time are ${f}_{3}$ and its higher harmonics. Subsequently, the system undergoes loop surface doubling (Fig. 6), and obvious nonlinear characteristics appear in the spectrum, which is represented by the high amplitude of the subharmonic ${f}_{3}$/2. After the system is stable, it maintains the quasiperiodic motion again (Fig. 7). The main peaks of the spectrum are still ${f}_{3}$ and its higher harmonics, but the amplitude of ${f}_{3}$ increases with the increase of excitation frequency (Fig. 8).
Fig. 5. Dynamic characteristics of the planetary gear ${\mathrm{\Omega}}_{1}=$ 2.4
Fig. 6. Dynamic characteristics of the planetary gear ${\mathrm{\Omega}}_{1}=$ 2.5
Fig. 7. Dynamic characteristics of the planetary gear ${\mathrm{\Omega}}_{1}=$ 2.9
Fig. 8. Dynamic characteristics of the planetary gear ${\mathrm{\Omega}}_{1}=$ 4.4
When ${\mathrm{\Omega}}_{1}\ge $ 1.7, the vibration characteristics of the planetary gears are more prominent than the fixedaxis gear. The increase of the planetary gear clearance has little effect on low frequencies, but has a great influence on high frequencies, which greatly increases the amplitude of ${f}_{r}$. The planetary vibrations are only sensitive to high excitation frequencies, and their characteristics are not visible at low excitation frequencies. However, such high excitation frequencies are often not achieved in engineering applications. Therefore, planetary gear wear is relatively the least identifiable fault in the system.
3.3. Frequency characteristics of the planetary gears under different degrees of wear
The interval where the vibration characteristics of the planetary gear are obvious is selected as the research object, that is, when ${\mathrm{\Omega}}_{1}=$ 3. The frequency characteristics of the planetary gear changed with the degree of planetary wear were investigated, as shown in Fig. 9.
It can be seen from Fig. 9 that as the clearance increases, the planetary gear keeps periodic motion, and the main peaks of the frequency are always ${f}_{3}$ and its higher harmonics, but the amplitude of the frequency increases significantly. When the clearance is small, the planetary gear is in singleperiod motion, and the vibration amplitude of the meshing frequency is also small (Fig. 9(a)). The clearance gradually increases. When ${\stackrel{}{e}}_{aspn}={\stackrel{}{e}}_{arpn}=$ 0.20.6, the amplitude increases rapidly (Figs. 9(b)(d)). At this time, the wear enters the development period, and the slight increase of wear will be reflected in the spectrum. When the wear is ${\stackrel{}{e}}_{aspn}={\stackrel{}{e}}_{arpn}=$ 0.61, the amplitude growth becomes slower (Figs. 9(e) and (f)).
Fig. 9. Frequency characteristics of the planetary gears: a) ${\stackrel{}{e}}_{aspn}={\stackrel{}{e}}_{arpn}=$ 0.1, b) ${\stackrel{}{e}}_{aspn}={\stackrel{}{e}}_{arpn}=$ 0.2, c) ${\stackrel{}{e}}_{aspn}={\stackrel{}{e}}_{arpn}=$ 0.4, d) ${\stackrel{}{e}}_{aspn}={\stackrel{}{e}}_{arpn}=$ 0.6, e) ${\stackrel{}{e}}_{aspn}={\stackrel{}{e}}_{arpn}=$ 0.8, f) ${\stackrel{}{e}}_{aspn}={\stackrel{}{e}}_{arpn}=$ 1
a)
b)
c)
d)
e)
f)
3.4. The influence of the planetary gear wear on the fixedaxis gear
The above research results only appear in the experimental simulation. In practical engineering applications, the planetary gear signals are often coupled with the fixedaxis gear signals to become complex mixed signals, which making the planetary gear signals difficult to identify. The 1st stage fixedaxis signal is closer to the input end, so the signal interference is small, and the vibration frequency is high. In the extracted experimental signals, the information components are many and clear. Therefore, frequency characteristics on the meshing point of the 1st stage fixedaxis gear are studied to investigate whether the planetary gear wear signals can affect the fixedaxis ones. The frequency characteristics on the 1st stage fixedaxis meshing point were investigated under different degrees of wear and different excitation frequencies. The results are shown in Figs. 10, 11.
Fig. 10. Frequency spectrum diagram of the 1st stage fixedaxis gear, when ${\stackrel{}{e}}_{aspn}={\stackrel{}{e}}_{arpn}=$ 0.2: a) ${\mathrm{\Omega}}_{1}=$1, b) ${\mathrm{\Omega}}_{1}=$ 2, c) ${\mathrm{\Omega}}_{1}=$ 3, d) ${\mathrm{\Omega}}_{1}=$ 4
a)
b)
c)
d)
Fig. 11. Frequency spectrum diagram of the 1st stage fixedaxis gear, when ${\stackrel{}{e}}_{aspn}={\stackrel{}{e}}_{arpn}=$ 0.4: a) ${\mathrm{\Omega}}_{1}=$ 1, b) ${\mathrm{\Omega}}_{1}=$ 2, c) ${\mathrm{\Omega}}_{1}=$ 3, d) ${\mathrm{\Omega}}_{1}=$ 4
a)
b)
c)
d)
It can be seen from Fig. 2(g) that the excitation frequency range of the 1st stage fixedaxis gear affected by the wear of the planetary gear is mainly [3, 4]. Comparing Figs. 10(a), (b) and 11(a), (b), it can be seen that when ${\mathrm{\Omega}}_{1}=$ 1 and 2, the spectrum has no change. In Figs. 10(c) and (d), the amplitude of the planetary gear meshing frequency ${f}_{3}$ is masked due to the slight wear. But with the wear increase, ${f}_{3}$ appears in Figs. 11(c) and d. In Fig. 11(d), the amplitude of ${f}_{3}$ has exceeded ${f}_{2}$, indicating that the planetary gear wear characteristics can be reflected in the fixedaxis spectrum. But the planetary gear characteristics are only obvious when the wear progresses to the medium term and the excitation frequency is sufficiently large.
4. Test signals identification and diagnosis
Fig. 12 shows a multistage gear transmission test rig consisting of a twostage fixedaxis gear and an onestage planetary gear.
Fig. 12. The test rig of multistage gear transmission system: 1 – motor, 2 – torque sensor and encoder, 3 – two stage fixedaxis gearbox, 4 – radial load of bearing, 5 – one stage planetary gearbox, 6 – brake
The planetary gear 1 in the test rig is the wear gear (see Fig. 13). It was found in 3.1 that the excitation frequency range affected by the planetary gear wear was [3, 4]. The closer the excitation frequency is to the region, the more obvious the vibration characteristics of the planetary gear are. In order to identify changes in the planetary gear wear signals, the excitation frequency was varied by adjusting the motor speed of the test rig. Since the maximum speed of the test rig is 3000 rpm, the signals were tested and analyzed with the motor speed of 1800 rpm (corresponding to the excitation frequency ${\mathrm{\Omega}}_{1}=$ 1) and 2400 rpm (corresponding to the excitation frequency ${\mathrm{\Omega}}_{1}=$ 1.5). In order to compare the measured spectrums with the dimensionless simulation results, the dimensionless spectrums of the vertical direction of the planetary gearbox at different speed were obtained in Fig. 14. The dimensionless characteristic frequencies of each stage gears of the multistage gear transmission system is shown in Table 3.
Fig. 13. Wear planetary gear
a)
b)
As can be seen from Fig. 14(b), planetary gear meshing frequency ${f}_{3}$ is not significant when the excitation frequency is low. As the excitation frequency increases (Fig. 14(d)), planetary gear meshing frequency ${f}_{3}$ amplitude increases. This conclusion is consistent with the numerical simulation result, indicating that the planetary gear has obvious vibration characteristics at high excitation frequency. It can be seen from the numerical simulation results that the increase of the planetary gear wear will increase the amplitude of the meshing frequency of the planetary gear and its higher harmonics but will not generate other side frequency components. Through the variation of the amplitude of the planetary gear meshing frequency, the numerical simulation results were used to identify the planetary gear wear.
Fig. 14. Comparison of spectrum of the planetary gearbox at different speeds: a) motor speed 1800 rpm timedomain, b) motor speed 1800 rpm frequencydomain, c) motor speed 2400 rpm timedomain, d) motor speed 2400 rpm frequencydomain
a)
b)
c)
d)
Comparing the theoretical simulation results with the experimental results, the main frequency components are compared in detail. A comparison is drawn between Fig. 10(a) and Fig. 14(b) in which the excitation frequency is the same as ${\mathrm{\Omega}}_{1}=$ 1. Because of the different amplitudes of the two, for the sake of comparison, the magnitude of the dimensionless spectrum of the theoretical simulation is scaled down, as shown in Fig. 15.
Fig. 15. Comparison of main frequency between theoretical and experimental results
As can be seen from Fig. 15, peak values of ${f}_{1}$, ${f}_{2}$, ${f}_{3}$, 2${f}_{2}$ and ${f}_{1}$${f}_{2}$ are all found in both theoretical and experimental conclusions. The difference is that there is a lack of peaks associated with ${f}_{d}$ in the theoretical value. This is because the established stiffness and error models are functions of the meshing frequency, ignoring the motor frequency ${f}_{d}$. In the ratio of the energy of the frequency components, ${f}_{1}$ in the theoretical value is much larger than other frequency components. In the experimental values, the frequency components are proportionally balanced, and when the rotational speed is increased, the ${f}_{1}$ peak is highlighted (Fig. 14(d)).
5. Conclusions
In this paper, the nonlinear dynamic characteristics of the wear of planetary gear in a multistage gear transmission system were studied. Planetary gear wear was simulated through the change of the planetary gear clearance, and the bifurcation diagrams of the fixedaxis gear and planetary gear under different clearances were analyzed. The changes of the system's motion state were investigated. The vibration mechanism of the system with wear failure was summarized. The study found that the wear of the planetary gear mainly affected the excitation frequency range of [3, 4]. As the excitation frequency increased, the motion state of the system changed from periodic motion to quasiperiodic motion. The wear caused the excitation frequency of the system to enter the quasiperiodic motion to advanced. By studying the timedomain, frequencydomain, phase diagram and Poincaré section of each motion state of the fixedaxis gear and planetary gear under different clearances, the transition processes and frequency characteristics of each motion state were analyzed. Frequency characteristics of wear were obtained.
The vibration signals of the experimental rig were tested at different excitation frequencies. When the rotational speed was increased, the amplitude of the meshing frequency of the planetary gear increased. This result is the same as the theoretical study.
Acknowledgements
This research is supported by Baoji Municipal Science and Technology Bureau (2017JH211), China.
References
 Kuang J. H., Lin A. D. The effect of tooth wear on the vibration spectrum of a spur gear pair. Journal of Vibration and Acoustics, Vol. 123, Issue 3, 2001, p. 311317. [Publisher]
 Yuksel C., Kahraman A. Dynamic tooth loads of planetary gear sets having tooth profile wear. Mechanism and Machine Theory, Vol. 39, Issue 7, 2004, p. 695715. [Publisher]
 Wojnarowski J., Onishchenko V. Tooth wear effects on spur gear dynamics. Mechanism and Machine Theory, Vol. 38, Issue 2, 2003, p. 161178. [Publisher]
 Li S., Kahraman A. A micropitting model for spur gear contacts. International Journal of Fatigue, Vol. 59, Issue 2, 2014, p. 224233. [Publisher]
 Brand J. A., Martins R., Seabra J. H. O., et al. An approach to the simulation of concurrent gear micropitting and mild wear. Wear, Vol. 324, 2015, p. 6473. [Publisher]
 Zhang Y. C., Tian H. Q., Tang J. Y., et al. Study on sliding wear of high pair based on friction work principles. China Mechanical Engineering, Vol. 3, 2010, p. 344347. [CrossRef]
 He R. G., Jian Q. Y., Yao Y. F. Numerical simulation of tooth wearing for involute helical cylindrical gears. Lubrication Engineering, Vol. 32, Issue 3, 2007, p. 8891. [CrossRef]
 Liu Z. X., Liu Z. S., Zhao J. M., Zhang G. H. Study on interactions between tooth backlash and journal bearing clearance nonlinearity in spur gear pair system. Mechanism and Machine Theory, Vol. 107, 2017, p. 229245. [Publisher]
 Ding H. L., Kahraman A. Interactions between nonlinear spur gear dynamics and surface wear. Journal of Sound and Vibration, Vol. 307, Issue 3, 2007, p. 662679. [Publisher]
 Onishchenko V. Investigation of tooth wears from scuffing of heavy duty machine spur gears. Mechanism and Machine Theory, Vol. 83, 2015, p. 3855. [Publisher]
 Shen Y. J., Yang S. P., Xing H. J., Wang X. Y. Nonlinear dynamics of a spur gear pair with fault. Key Engineering Materials, Vol. 84, Issue 353, 2007, p. 11771180. [Publisher]
 Liu X., Yang Y., Zhang J. Investigation on coupling effects between surface wear and dynamics in a spur gear system. Tribology International, Vol. 101, 2016, p. 383394. [Publisher]
 Priest M., Taylor C. M. Automobile engine tribologyapproaching the surface. Wear, Vol. 241, Issue 2, 2000, p. 193203. [Publisher]
 Feng S., Fan B., Mao J., et al. Prediction on wear of a spur gearbox by online wear debris concentration monitoring. Wear, Vol. 336, 2015, p. 18. [Publisher]
 Wu H., Kwok N. M., Liu S., et al. A prototype of online extraction and threedimensional characterisation of wear particle features from video sequence. Wear, Vol. 368, 2016, p. 314325. [Publisher]
 Hu C., Smith W. A., Randall R. B., et al. Development of a gear vibration indicator and its application in gear wear monitoring. Mechanical Systems and Signal Processing, Vol. 76, Issue 77, 2016, p. 319336. [Publisher]
 Zhang R., Gu F., Mansa H., et al. Gear wear monitoring by modulation signal bispectrum based on motor current signal analysis. Mechanical Systems and Signal Processing, Vol. 94, 2017, p. 202213. [Publisher]
 He Q. B. Timefrequency manifold for nonlinear feature extraction in machinery fault diagnosis. Mechanical Systems and Signal Processing, Vol. 35, Issues 12, 2013, p. 200218. [Publisher]
 Wang X. Stability research of the planetary multistage gear transmission system with coupling faults. Journal of Sound and Vibration, Vol. 434, 2018, p. 6377. [Publisher]
 Flodin A., Andersson S. Simulation of mild wear in spur gears. Wear, Vol. 207, Issue 2, 1997, p. 1623. [Publisher]
 Park D., Kolivand M., Kahraman A. An approximate method to predict surface wear of hypoid gears using surface interpolation. Mechanism and Machine Theory, Vol. 71, 2014, p. 6478. [Publisher]
 Wang Y. G., Zheng H. Q., Yang T. Q., Guan Z. Z., Yang J. Nonlinear dynamics behavior of gear system with fault parameters. Journal of Vibration, Measurement and Diagnosis, Vol. 31, Issues 52011, 570, p. 573662. [CrossRef]
 Wang X., Xu Y. X., Wu B. L. Study on failure characteristics of gearbox transmission system with coupling faults. Journal of Vibration and Shock, Vol. 36, Issue 12, 2017, p. 217223. [CrossRef]