Abstract
In a multistage gear transmission system, the motion state of the system will change with the excitation frequency, and the frequency characteristics will also change accordingly. If this change is not taken into account, there is often a great deviation in identifying and judging system faults according to unified standards, especially when the system has such early undetectable fault as crack. In this paper, the dimensionless differential equations of motion of multistage gear transmission system are established. The stiffness model of gear tooth crack is established by potential energy method. The changes of the motion state of the system with the increase of excitation frequency are obtained by calculating the displacement bifurcation diagram of cracked gear tooth. The influences of crack fault on each motion state are studied by using time domain, frequency domain, phase diagram and Poincaré cross section, and the fault frequency characteristics are summarized. By comparing the theoretical and experimental data of the vibration response characteristics of the system, the motion state of the system can be effectively determined and the crack fault can be identified.
Highlights
 The stiffness model of gear tooth crack is established by potential energy method.
 The changes of the motion state of the system with the increase of excitation frequency are obtained by calculating the displacement bifurcation diagram of cracked gear tooth.
 The influences of crack fault on each motion state are studied by using time domain, frequency domain, phase diagram and Poincaré cross section, and the fault frequency characteristics are summarized.
1. Introduction
Fatigue crack has always been one of the main forms of gear failure. It is difficult to identify in the early stage, but when the crack is discovered in the middle and later stage, it may have a great impact. Therefore, how to identify the crack fault earlier and more effectively is the focus of scholars' research. Previous scholars have focused on studying the bending fatigue strength of a single gear tooth and testing the bending fatigue life from the perspective of material fracture mechanics [1]. Although these results have proved effective. However, the research results have not been popularized in the multistage gear system and lack of effective diagnostic standards. Due to the complexity of the working conditions of the multistage gear system, the double increase of manpower and material resources consumed in the experiment, and the coupling characteristics of the system, the vibration response of the crack fault of the multistage gear system is not clear. Therefore, it is particularly important to find out the vibration response of the multistage gear system crack failure, and to combine the existing theories with practical applications to propose effective diagnostic indicators.
For the research of crack failure, a large number of scholars focus on crack model [2], strength analysis [3] and life prediction [4]. Using finite element software simulation method to test the bending fatigue life of gears has become an important means to study the bending fatigue of gears at home and abroad [5]. Although scholars have continuously improved finite element software, it is still difficult to obtain accurate data for the complexity of gears only by software simulation. In order to obtain highprecision bending fatigue life data of gears, it is necessary to combine actual tests to assist in completing the bending fatigue life test of gears. The test method of comprehensive experiment and simulation has been accepted by engineers in various countries. However, finite element software is only suitable for the strength analysis of a single gear. The dynamic solution of the system is required for the mutual coupling between multistage gears.
Because research on dynamic coupling mechanism and fault mechanism of multistage gear transmission system is not perfect, in order to solve the problem of fault identification, many scholars choose to analyze the experimental signals. With the rapid development of big data, neural network, machine learning and other disciplines in the past decade, many new analysis algorithms for complex signals are proposed, such as ant colony optimization algorithm [6], structured sparsity timefrequency analysis [7], particle swarm optimization (PSO) neural network [8], support vector machine [9], multiclass relevance vector machine (mRVM) [10], convolutional neural networks (CNN) [11], rotational system noncoaxiality [12], etc.
It is also a mainstream research direction to synthesize various algorithms to improve recognition efficiency. For example, Zhixiong Li [1315] uses wavelet packet, empirical mode decomposition, Wigner distribution, AR model, blind source separation and other hybrid intelligent signal analysis methods to identify a variety of coupling faults. In addition, vibration signals collected by multichannel sensors are used to monitor the bounded component analysis and detect the gear cracks of complex gearbox [16]. Yonglong Li [17] combines wavelet packet transform with support vector machine to detect gear fault. Xinbin Li [18] proposed a method combining EMD and support vector machine to identify different bearing faults. Dikang Peng [19] use transmission error and meshing phasing to distinguish local fault types and identify faulty gears. In recent years, people have carried out indepth research on composite fault detection using intelligent methods. The intelligent classifier [20] can well identify different composite fault modes. But the existing intelligent algorithm lack of physical significance. In other words, it is difficult to understand how neural networks work [15].
One way to solve the above problems is to separate the fault by decoupling and separating the subsignals corresponding to the fault. In the past few decades, great progress has been made in the decoupling diagnosis of complex gearbox signals, including wavelet decoupling, EMD decoupling, sequential tracking decoupling, sparse decomposition decoupling and independent component analysis decoupling [21]. However, there are few researches on the correlation of fault vibration sources. In any case, components work together to complete the function of a mechanism. Inevitably, the close connection and coupling effect between components will affect their dynamic response and performance. The vibration signals excited by different components and / or faulty components are related to a certain extent. This dependency/correlation may have a significant impact on the results of fault diagnosis, leading to misdiagnosis or missed diagnosis [22]. Unfortunately, although there is often a correlation in the actual complex vibration signal, the problem has not been well studied in the fault decoupling diagnosis. The correlation of fault vibration sources will greatly increase the difficulty of fault detection. Therefore, it is very important to find out the vibration source and the correlation of faults, and excavate the vibration response.
In order to analyze the fault vibration response of a multistage system from the source, it is necessary to study from the perspective of nonlinear dynamics. By establishing a nonlinear dynamic model, the changes of motion state, vibration characteristics and fault state of the multistage system are studied step by step from the forward direction. For the nonlinear dynamics modeling of gear system, there are many researches on the nonlinear dynamics of singlestage gear. For example, Wang [23] studied the bendingtorsion coupling response of spur gear system under complex excitation, and Lin [24] calculated the transmission error using the bendingtorsion axial coupling dynamics model of transmission system. Research on multistage gear system is not perfect due to its complexity and variety. The nonlinear dynamic models that can be used as references include twoparallel shaft gear transmission system [25], twostage closedform planetary gear set [26], twostage helical planetary gear transmission system [27], multistage planetary gear system [28], multiple degrees of freedom nonlinear system [29], etc. This paper will study a multistage gear transmission system composed of a twostage fixedaxis gearbox and a onestage planetary gearbox. The nonlinear dynamic model in the author's previous research work [30] is used for simulation calculation.
In the author’s previous research, it was found that the motion state of the multistage gear transmission system is variable [31]. With the change of motion state, its vibration characteristics were also changed. The characteristics include the correlation between the gears at all stages. The characteristics of faults are different under different conditions. The change factor was not taken into account in previous studies. There may be a great deviation in identifying and judging system faults according to a unified standard. In this paper, fault characteristics under different motion states are discussed respectively. The vibration response of the system is excavated. The motion state and fault characteristics of the system are identified by comparing with the experimental data. This study has strong theoretical value and engineering application value.
2. System dynamics model
The gear transmission system shown in Fig. 1 was studied. There are two stage fixed shaft gears on the left and a planetary gear on the right. The lateral vibration displacements of the gears are not considered, only the torsional vibration is considered. Where spur gear 1, 2 compose the 1st stage fixedaxis gear for the input, spur gear 3, 4 compose the 2nd stage fixedaxis gear, the planet carrier is for the output. In Fig. 1. ${\theta}_{s}$, ${\theta}_{c}$, ${\theta}_{pn}$, ${\theta}_{1}$, ${\theta}_{2}$, ${\theta}_{3}$, ${\theta}_{4}$ represent the angular displacement of sun gear, planet carrier, planetary gear $n$ ($n=$1, 2, 3, 4), spur gear 1, 2, 3, 4, respectively. Throughout this paper, the subscripts $s$, $c$, $pn$, $r$, 1, 2, 3, 4 denote sun gear, planet carrier, planetary, ring gear and spur gear 1, 2, 3, 4. Quantities ${r}_{s}$,_{}${r}_{c}$,_{}${r}_{pn}$,_{}${r}_{1}$,_{}${r}_{2}$,_{}${r}_{3}$, ${r}_{4}$ are the base circle radius of gears. Quantities ${K}_{spn}\left(t\right)$, ${K}_{rpn}\left(t\right)$, ${K}_{1}\left(t\right)$, ${K}_{2}\left(t\right)$ denote the meshing stiffness of sun gear with planetary gear $n$, ring gear with planetary gear $n$, 1st stage fixed gear and 2nd stage fixed gear. Quantities ${C}_{spn}$, ${C}_{rpn}$, ${C}_{1}$, ${C}_{2}$ denote the damping of sun gear with planetary gear $n$, ring gear with planetary gear $n$, 1st stage fixed gear and 2nd stage fixed gear. ${T}_{in}$ is the input and ${T}_{out}$ is the output.
Fig. 1Gear transmission model [30]
2.1. Dimensionless motion differential equation
The motion differential equation of the system was established, and the dimensionless treatment was carried out to obtain Eq. (1). The specific derivation process had been described in detail in the author’s paper [30]:
where, ${x}_{i}$ is the relative displacement of each stage meshing line, ($i=$1, 2, $spn$, $rpn$), ${x}_{1}={r}_{1}{\theta}_{1}{r}_{2}{\theta}_{2}{e}_{1}\left(t\right)$, ${x}_{2}={r}_{3}{\theta}_{3}{r}_{4}{\theta}_{4}{e}_{2}\left(t\right)\text{,}$${x}_{spn}={r}_{s}{\theta}_{s}{r}_{pn}{\theta}_{pn}{r}_{c}{\theta}_{c}{e}_{spn}\left(t\right)\text{,}$${x}_{rpn}={r}_{pn}{\theta}_{pn}{r}_{c}{\theta}_{c}{e}_{rpn}\left(t\right)$; ${m}_{e1}$, ${m}_{e2}$, ${m}_{e3}$, ${m}_{e4}$ are the equivalent mass of gear, ${J}_{ce}$ is the equivalent moment of inertia of gear, ${J}_{23}={J}_{2}+{J}_{3}$, ${J}_{4s}={J}_{4}+{J}_{s}$: ${J}_{1}$, ${J}_{2}$, ${J}_{3}$, ${J}_{4}$, ${J}_{s}$, ${J}_{pn}$, ${J}_{c}$ are the moment of inertia of each gear on the shaft:
$f\left({x}_{i}\right)$ is the clearance nonlinear function, ($i=$1, 2, $spn$, $rpn$), written as:
where, ${b}_{i}$ is half of the clearance, ($i=$1, 2, $spn$, $rpn$).
The damping coefficient form:
where, ${\xi}_{1}$, ${\xi}_{2}$, ${\xi}_{spn}$, ${\xi}_{rpn}$ are damping ratio and ${m}_{1}$, ${m}_{2}$, ${m}_{3}$, ${m}_{4}$, ${m}_{s}$, ${m}_{pn}$, ${m}_{r}$ are the mass of each gear.
The comprehensive meshing error of gear pair using the 1st harmonic form of meshing function, that is:
where, ${e}_{ai}$ is the comprehensive meshing error amplitude of each gear, ($i=$1, 2, $spn$, $rpn$). ${\phi}_{i}$ is the comprehensive meshing error initial phase of each gear, ($i=$1, 2, $spn$, $rpn$); ${w}_{mi}$ is the meshing frequency of each gear, ($i=$1, 2, $spn$, $rpn$).
Define time nominal scale ${w}_{h}$, order $\tau ={w}_{h}t$, where ${w}_{h}=\sqrt{{K}_{m1}\times {m}_{e1}}$; ${\stackrel{}{x}}_{i}$ is dimensionless displacement, ${\stackrel{}{x}}_{i}={x}_{i}/{b}_{1}$, $i=$1, 2, $spn$, $rpn$; ${\mathrm{\Omega}}_{i}$ is dimensionless excitation frequency, ${\mathrm{\Omega}}_{i}={w}_{mi}/{w}_{h}$, $i=$1, 2, $spn$, $rpn$; ${\stackrel{}{e}}_{ai}$ is dimensionless comprehensive meshing error amplitude, ${\stackrel{}{e}}_{ai}={e}_{ai}/{b}_{1}$, $i=$1, 2, $spn$, $rpn$. ${\mathrm{\Omega}}_{i}$ is the dimensionless excitation frequency of the 1ststage fixedaxis gear. When ${\mathrm{\Omega}}_{i}=$1, it is the critical speed.
Fig. 2Meshing stiffness of gear tooth with crack propagation
2.2. Failure models
The stiffness of cracked teeth is calculated by potential energy method [2]. The simulated crack occurs at the tooth root of gear 1 and extends inward at an angle of 70 ° to the center line. After reaching the center line, it changes the propagation direction and propagates to the other tooth root until fracture. The gear parameters of the system are shown in Table 1. The meshing stiffness of cracked gears with crack propagation is calculated as shown in Fig. 2.
Table 1Gear parameters
Gear  Number of teeth  Mass ${m}_{i}$ / g  ${J}_{i}$ / (g·m^{2})  Face width / mm 
1  29  125  0.05  30 
2  100  1224.5  6  30 
3  36  224  0.14  30 
4  90  1111  4  20 
$s$  28  41  0.007  20 
$pn$  36  34.6  0.01  20 
$c$  848.7  0.76  20  
$r$  100  20 
3. Variation of nonlinear dynamic characteristics of crack fault
3.1. Gear displacement bifurcation diagram
In order to compare the motion characteristics of the system under normal and fault conditions, the bifurcation diagrams of the system varying with the excitation frequency under the condition of no crack and 75 % crack were calculated. The structural parameters of the system are: ${T}_{in}=$6.5 N·m, ${T}_{out}=$8.5 N·m, gear clearance ${b}_{1}=$3 μm, comprehensive meshing error amplitude ${e}_{a}=$5 μm, meshing pair damping ratio $\xi =$0.07, gear contact ratio is 1.68. As the crack fault was added to gear 1, the bifurcation diagram of gear 1 was calculated. The variable step RungeKutta method was used to solve the Eq. (1), and the bifurcation diagrams of the relative displacement of the fixed shaft gear 1 varying with the dimensionless excitation frequency ${\mathrm{\Omega}}_{1}$ were obtained, as shown in Fig. 3.
The bifurcation diagram is a twodimensional graph showing that the dynamic system varies with a certain parameter, that is, to calculate the value of the fixed point corresponding to each parameter (discrete system) or the maximum value of oscillation (continuous system). Fig. 3 shows the oscillatory change of the relative displacement of gear 1 with the increase of the dimensionless excitation frequency ${\mathrm{\Omega}}_{1}$. In Fig. 3(a), when ${\mathrm{\Omega}}_{1}$ is small, ${\stackrel{}{x}}_{1}$ is a thin curve. It shows that the value of ${\stackrel{}{x}}_{1}$ is almost the same during the movement. The vibration has only one period, and the motion state is single periodic motion. When ${\mathrm{\Omega}}_{1}$ increase to 0.5, ${\stackrel{}{x}}_{1}$ changes from a thin curve to a wide curve, and the system changes from single periodic motion to quasiperiodic motion. At this time, the trajectory of ${\stackrel{}{x}}_{1}$ is densely distributed in the area defined by the amplitude. When ${\mathrm{\Omega}}_{1}$ increase to 0.8, the resonance of the system makes the quasiperiodic motion changes into chaotic motion. The trajectory of ${\stackrel{}{x}}_{1}$ becomes random and irregular, and the oscillation range increases. Then the system bifurcates into twoperiodic motion. After a short period of chaotic motion, it bifurcates into threeperiodic motion, which divided into three branches in the range of ${\mathrm{\Omega}}_{1}\in $[1.6, 2]. Finally, the system reenters the quasiperiodic motion.
Comparing Figs. 3(a) and 3(b), I found that when the system has a fixedaxis crack fault, it affects the periodic, quasiperiodic and threeperiodic motion of ${\mathrm{\Omega}}_{1}<$3. The fault periodic motion is added into the original periodic motion. That is, a new periodic motion amplitude caused by the fault shock appears outside the original amplitude. Although the overall motion state does not change, it is obvious that the shock is more obvious in the vibration characteristics, and the fault frequency in the spectrum also increases. The motion state of the system can be roughly seen from the bifurcation diagram. In order to study the dynamic characteristics of each motion state in more details, it need to make a detailed analysis of each state.
Fig. 3Displacement bifurcation diagram of gear 1: a) normal state, and b) fixedaxis gear 75 % crack fault state
a)
b)
3.2. Nonlinear characteristics of each motion state of the system
The nonlinear characteristics of different motion states are studied respectively. The time domain, frequency domain, phase diagram and Poincaré section of the system with the increase of ${\mathrm{\Omega}}_{1}$ are investigated. Phase diagram can express the state trajectory of the dynamic system on the phase plane. Poincaré section is a cross section of the trajectory of continuous motion, and the shape of the motion is succinctly judged according to the situation that the trajectory passes through the cross section. For each motion state, a sample is selected as the representative. The variable step size RungeKutta method is used to solve the nonlinear differential Eq. (1), and Fig. 4 is obtained. The characteristic frequencies of gears at all stages of the system are shown in Table 2.
Table 2Characteristic frequencies of gears at all stages
Characteristic frequency  Dimensional frequency 
Meshing frequency of the 1ststage fixedaxis gear ${f}_{1}$  1 
Meshing frequency of the 2ndstage fixedaxis gear ${f}_{2}$  0.3599 
Meshing frequency of the planetary gear ${f}_{3}$  0.0877 
Frequency of the 1ststage fixedaxis failure ${f}_{d}$  0.0345 
Frequency of the planetary gear failure ${f}_{r}$  0.0024 
As can be seen from Fig. 4, when ${\mathrm{\Omega}}_{1}=$ 0.4, the amplitude in both time domain and frequency domain are very small. The phase diagram is a tiny ring. The Poincaré section is a point. The system moves in single period motion. When ${\mathrm{\Omega}}_{1}=$ 0.6, the amplitude in time domain and frequency domain increases. The amplitude of the circle of phase diagram increases. The Poincaré section becomes a point cluster. The frequency characteristic is the 1ststage fixedaxis meshing frequency ${f}_{1}$. At this time, the system enters a quasiperiodic motion state.
When ${\mathrm{\Omega}}_{1}=$ 0.9, the system is in chaotic motion. The circle in the phase diagram has been filled, and the amplitude increases. The point cluster of Poincaré section increases. There are many frequency components in the spectrum. In addition to ${f}_{1}$, there are also ${f}_{2}$, its double frequency $2{f}_{2}$ and modulation frequency ${f}_{1}{f}_{2}$.
Due to the appearance of other frequency components, the vibration energy is no longer concentrated, which leads to the decrease of the amplitude of ${f}_{1}$. Then the system gradually bifurcates into twoperiodic motion (${\mathrm{\Omega}}_{1}=$ 1.4). At this time, Poincaré section is two points, phase diagram returns to a twoperiod circle, and the frequency spectrum is dominated by the subharmonic resonance frequency ${f}_{1}$/2.
When ${\mathrm{\Omega}}_{1}=$ 1.5, the circle in the phase diagram is filled again, the Poincaré section becomes a point cluster, and the system experiences a short period of chaotic motion. Then it bifurcates into triple periodic motion (${\mathrm{\Omega}}_{1}=$ 2). Poincaré section bifurcates into 3 clusters. Phase diagram is a threeperiod circle. Subharmonic resonance frequency ${f}_{1}$/3 is included in the spectrum. The triple periodic motion lasts a little longer, and the subharmonic resonance frequency ${f}_{1}$/3 is found in the following experimental signals.
Fig. 4Nonlinear characteristics in each state of motion: a) Ω1= 0.4, b) Ω1= 0.6, c) Ω1= 0.9, d) Ω1= 1.4, e) Ω1= 1.5, f) Ω1= 2, g) Ω1= 2.5, and h) Ω1= 3
a)
b)
c)
d)
e)
f)
g)
h)
When ${\mathrm{\Omega}}_{1}=$2.5 and 3, the system returns to quasiperiodic motion. The Poincaré section is a limit cycle and gradually reduces to a point cluster. The phase diagram changes into a single periodic ring. The amplitude of ${f}_{1}$ is the highest in the spectrum. In addition, there is ${f}_{2}$ in the spectrum.
It can be seen that the vibration frequency characteristics of the system will change according to different excitation frequencies and motion states. If we do not explore the vibration characteristics of the system in each motion state, but simply take the vibration characteristics of a certain time as the benchmark sample for reference and comparison, the conclusion will appear a lot of deviation. The change of vibration frequency peak may not be caused by a fault but by the change of motion state. In order to study the effect of cracks on the vibration characteristics of the system, the analysis should be carried out separately for different motion states.
4. The influence of crack fault on each motion state
It can be seen from Section 3.2 that the motion states of the system in this paper are mainly divided into several categories: single periodic motion (${\mathrm{\Omega}}_{1}=$ 0.4, 0.6), chaotic motion (${\mathrm{\Omega}}_{1}=$ 0.9, 1.5), bifurcation (${\mathrm{\Omega}}_{1}=$ 1.4, 2) and quasi periodic motion (${\mathrm{\Omega}}_{1}=$ 2.5, 3). The author finds that with the increase of crack length, the change of fault characteristics under the same motion state is consistent. Therefore, only one sample is listed for each motion state to explain the changing trend of fault characteristics.
4.1. Single periodic motion
In the state of single period operation, when ${\mathrm{\Omega}}_{1}=$ 0.6, the change of the nonlinear characteristics of the system is shown in Fig. 5 when the crack increases from 0 to 100 %.
As can be seen from Fig. 5, the initial crack propagation has little effect on the motion state of the system, and there is no obvious fault feature, which is difficult to detect. It is after the crack is more than half that the fault frequency of gear appears around the first meshing frequency and increases exponentially with the increase of crack. This is why the early crack failure is hard to detect. According to the principle of gear bending fatigue fracture mechanics, scholars divide the crack bending fatigue life into three stages: low speed expansion, stable expansion and unstable expansion [5]. Nonlinear vibration characteristics of the gear are basically consistent with the three phases. It shows that the model is reasonable and the research results are universal.
Fig. 5The change of nonlinear characteristics of the system caused by cracks when Ω1= 0.6: a) Perfect, b) 10 % crack, c) 25 % crack, d) 50 % crack, e) 75 % crack, and f) Tooth missing.
a)
b)
c)
d)
e)
f)
4.2. Changes of nonlinear characteristics of the system in different motion states
The nonlinear characteristics of the system with crack fault under chaotic motion (${\mathrm{\Omega}}_{1}=$ 0.9), triple periodic motion (${\mathrm{\Omega}}_{1}=$ 2) and quasiperiodic motion (${\mathrm{\Omega}}_{1}=$ 2.5) are calculated respectively, as shown in Figs. 68. From Section 4.1, it is known that the nonlinear fault responses with cracks of 10 % and 25 % are not obvious, so they are no longer listed. Since the nonlinear characteristic diagram in normal state has been given in Fig. 4. Only the nonlinear characteristic diagrams with cracks of 50 % and 100 % in three kinds of motion states are listed.
Compared with Figs. 6 and 4(c), it can be seen that the impact of the system in the time domain is greater under chaotic motion. However, due to the irregular motion of chaotic motion itself, the phase diagram and Poincaré section have no obvious change before tooth missing. In frequency domain, the fault sideband mainly occurs at the 2ndstage meshing double frequency $2{f}_{2}$.
Comparing Figs. 7 and 8 with Figs. 4(f) and 4(g), both triple periodic motion and quasiperiodic motion operate stably without impact under normal conditions. The phase diagram is a circle. After the occurrence of crack failure, the impact increases. The impact of triple periodic motion is more obvious. The phase diagrams are filled with fault cycles. The Poincaré section of the triple periodic motion gathers inward, while the quasiperiodic motion spreads outward. In the spectrum diagram, the fault period occurs around ${f}_{1}$/3.
Thus, it can be seen that even if the same fault occurs in different motion states, the location of the fault sideband is also different.
Fig. 6The change of nonlinear characteristics of the system caused by cracks when Ω1= 0.9: a) 50 % crack, and b) tooth missing
a)
b)
Fig. 7The change of nonlinear characteristics of the system caused by cracks when Ω1= 2: a) 50 % crack, and b) tooth missing
a)
b)
Fig. 8The change of nonlinear characteristics of the system caused by cracks when Ω1= 2.5: a) 50 % crack, and b) tooth missing
a)
b)
5. Crack fault experimental analysis
Fig. 9 is a gear transmission system test rig. The axial direction of the input end of the fixedaxis gearbox was selected as the measuring point. The arrangement of the vibration accelerometer is shown in Fig. 9. The parameters are shown in Tables 1 and 2. The test rig was tested and analyzed in normal and crack states. The acceleration sensor was used for signal acquisition with a sampling frequency of 3000 Hz.
Fig. 9The multistage gear transmission system test rig: 1 – motor; 2 – torque sensor and encoder; 3 – twostage fixedaxis gearbox; 4 – radial bearing load; 5 – onestage planetary gearbox; 6 – brake
In order to investigate the variation of fault characteristics with excitation frequency, the speedup process signal of the motor was tested and analyzed. The input frequency was set to 10 Hz, 20 Hz, 30 Hz and 40 Hz respectively (corresponding to the motor speed 600 r/min, 1200 r/min, 1800 r/min, 2400 r/min). The vibration signals of the axial measuring point of the fixedaxis gearbox were collected. In order to compare with the theoretical data, the spectrum was dimensionless. The vibration signal spectrum of the system under crack fault is shown in Fig. 10.
As can be seen from Fig. 10, the vibration responses of the system are different at different rotational speeds (excitation frequencies), which are consistent with the results of theoretical analysis. Comparing the frequency diagrams of Fig. 10 with Fig. 4, the motion state of the system at different rotational speeds can be roughly judged. For example, Fig. 10(a) is similar to Fig. 4(b), the system is in single periodic motion. Fig. 10(b) is similar to Fig. 4(c), the system is in chaotic motion. Fig. 10(c) is similar to Fig. 4(f), the system is in triple periodic motion. Fig. 10(d) is similar to Fig. 4(g), the system is in quasiperiodic motion.
In engineering applications, the 1ststage fixedaxis failure frequency ${f}_{d}$ and its frequency multiplication are generally selected to monitor the 1ststage fixedaxis gear failure. As can be seen in Fig. 10, ${f}_{d}$ appears around the main peak frequency in the form of sideband. This indicates that the system contains a 1ststage fixedaxis gear fault. However, due to different excitation frequencies, the locations of fault features are different.
Fig. 10Dimensionless spectrum of crack fault at different speeds: a) motor speed 600 r/min, b) motor speed 1200 r/min, c) motor speed 1800 r/min, and d) motor speed 2400 r/min
a)
b)
c)
d)
In the state of single periodic motion, the fault characteristic frequency ${f}_{d}$ appears around ${f}_{1}$, which is consistent with the conclusion of the study in Section 4.1. That is, the characteristic of Fig. 10(a) is consistent with that of Fig. 5. During chaotic motion, ${f}_{d}$ appears around ${f}_{1}$${f}_{2}$, and the characteristic of Fig. 10(b) is consistent with that of Fig. 6. During triple periodic motion, ${f}_{d}$ appears around ${f}_{1}$/3, and the characteristic of Fig. 10(c) is consistent with that of Fig. 7. During quasiperiodic motion, ${f}_{d}$ appears around ${f}_{1}$. Thus, it can be seen that the theoretical results are almost consistent with the experimental results. The influence of excitation frequency on the motion state and fault characteristic of the system is verified.
Because of the structural diversity of multistage gear transmission system, the change with excitation frequency is different under different system structures. The relationship between crack length and vibration characteristic needs to be verified by a large number of experiments. At present, no universal comprehensive fault identification method has been constructed. In order to apply the theoretical research results to practice, the author will try to combine theoretical research results with algorithms such as machine learning and neural network. The simulation data under different conditions will be used as “training” samples to try to improve the accuracy of fault identification.
6. Conclusions
In order to deeply explore the vibration characteristics and crack fault characteristics of multistage gear transmission system under different motion states, the vibration mechanism of crack fault is explored from a nonlinear point of view. According to the testbed of multistage gear transmission system, the dynamic simulation model of the system and the crack stiffness model of gear teeth are established. Compare the nonlinear characteristics of the system under different excitation frequencies, the transition process of different motion states of the multistage gear transmission system was summarized. According to different motion states, the nonlinear dynamic changes of the system with crack propagation are studied respectively. The fault characteristics of the system under different motion states are summarized. After the dimensionless processing of the experimental signal, the motion state of the system and the location of the fault can be roughly judged from the nonlinear point of view according to the theoretical research results. This study provides a theoretical basis for the vibration and fault response of multistage gear transmission system.
References

Daniewicz S. Conception and Development of Improved Analytical Prediction Models for Fatigue Induced Tooth Breakage Due to Cyclic Bending in Spur Gear Teeth. ProQuest Dissertations and Theses Global: ProQuest LLC, 1991.

Liang X. H., Ming J. Z., Mayank P. Analytically evaluating the influence of crack on the mesh stiffness of a planetary gear set. Mechanism and Machine Theory, Vol. 76, 2014, p. 2038.

Silori P., Shaikh A., Kumar K. C. N., Tandon T. Finite Element Analysis of Traction Gear Using ANSYS. Materials Today: Proceedings, Vol. 2, Issues 45, 2015, p. 22362245.

Patil V., Chouhan V., Pandya Y. Geometrical complexity and crack trajectory based fatigue life prediction for a spur gear having tooth root crack. Engineering Failure Analysis, Vol. 105, 2019, p. 444465.

Wang Z. Y., Xing Z. G., Wang H. D., Li G. L., Liu K. J., Xing Z. Research status of test method for bending fatigue life of heavy duty gear. Materials Review, Vol. 32, Issue 17, 2018, p. 30513059.

Deng W., Xu J. J., Zhao H. M. An improved ant colony optimization algorithm based on hybrid strategies for scheduling problem. IEEE Access, Vol. 7, 2019, p. 2028120292.

Sun R. B., Yang Z. B., Chen X. F., Tian S. H., Xie Y. Gear fault diagnosis based on the structured sparsity timefrequency analysis. Mechanical Systems and Signal Processing, Vol. 1021, 2018, p. 346363.

Pan H., Ma Q. Research on gearbox fault diagnosis method based on adjustinglearningrate PSO neural network. Journal of Donghua University, Vol. 23, 2006, p. 2932.

Wu F., Meng G. Compound rub malfunctions feature extraction based on fullspectrum cascade analysis and SVM. Mechanical Systems and Signal Processing, Vol. 20, 2006, p. 20072021.

Lei Y. G., Liu Z. Y., Wu X. H., Li N. P., Chen W., Lin J. Health condition identification of multistage planetary gearboxes using a mRVMbased method. Mechanical Systems and Signal Processing, Vol. 60, Issue 61, 2015, p. 289300.

Kolar D., Lisjak D., Pająk M., Pavkovi´ D. Fault diagnosis of rotary machines using deep convolutional neural network with wide three axis vibration signal input. Sensors, Vol. 20, Issue 14, 2020, p. 40174019.

Pająk M., Muślewski Ł., Landowski B., Grządziela A. Fuzzy Identification of The Reliability State of The Mine Detecting Ship Propulsion System. Polish Maritime Research, Vol. 26, Issue 1, 2019, p. 5564.

Li Z. X., Yan X. P., Yuan C. Q., Zhao J. B., Peng Z. X. A fault diagnosis approach for gears using multidimensional features and intelligent classifier. Noise and Vibration Worldwide, Vol. 41, Issue 10, 2010, p. 7688.

Li Z. X., Yan X. P., Yuan C. Q., Peng Z. X., Li L. Virtual prototype and experimental research on gear multifault diagnosis using waveletautoregressive model and principal component analysis method. Mechanical Systems and Signal Processing, Vol. 25, Issue 7, 2011, p. 25892607.

Li Z. X., Peng Z. X. A new nonlinear blind source separation method with chaos indicators for decoupling diagnosis of hybrid failures: A marine propulsion gearbox case with a large speed variation. Chaos, Solitons and Fractals, Vol. 89, 2016, p. 2739.

Li Z. X., Yan X. P., Wang X. P., Peng Z. X. Detection of gear cracks in a complex gearbox of wind turbines using supervised bounded component analysis of vibration signals collected from multichannel sensors. Journal of Sound and Vibration, Vol. 371, Issue 9, 2016, p. 406433.

Li Y., Shao R., Cao J. A new and effective method of gear fault diagnosis using wavelet packet transform combined with support vector machine. Journal of Northwestern Polytechnical University, Vol. 28, 2010, p. 530535.

Li X., Chen Y., Zhang S. Hybrid fault diagnosis algorithm based on fusion decision of multiple LSSVM classifiers. Journal of Sound and Vibration, Vol. 32, 2013, p. 159164.

Peng D. K., Smith W. A., Borghesani P., Randall R. B., Peng Z. X. Comprehensive planet gear diagnostics: Use of transmission error and mesh phasing to distinguish localised fault types and identify faulty gears. Mechanical Systems and Signal Processing, Vol. 12715, 2019, p. 531550.

Lei Y., He Z., Zi Y. Application of a novel hybrid intelligent method to compound fault diagnosis of locomotive roller bearings. Journal of Vibration and Acoustics, Vol. 130, 2008, p. 034501.

Li Z. X., Jiang Y., Hu C., Peng Z. Recent progress on decoupling diagnosis of hybrid failures in gear transmission systems using vibration sensor signal: A review. Measurement, Vol. 90, 2016, p. 419.

Li Z. X., Ge S. R., Zhu H. Key issues in the wear fault monitoring and diagnosis for critical components of coal cutters under deep coal seam. Tribology, Vol. 34, 2014, p. 729730.

Wang J. G., Zhang J., Yao Z. Y., Yang X. F., Sun R., Zhao Y. X. Nonlinear characteristics of a multidegreeoffreedom spur gear system with bendingtorsional coupling vibration. Mechanical Systems and Signal Processing, Vol. 12115, 2019, p. 810827.

Lin T. J., He Z. Y. Analytical method for coupled transmission error of helical gear system with machining errors, assembly errors and tooth modifications. Mechanical Systems and Signal Processing, Vol. 91, 2017, p. 167182.

Li W., Sun J. D., Yu J. P. Analysis of dynamic characteristics of a multistage gear transmission system. Journal of Vibration and Control, Vol. 25, Issue 10, 2019, p. 16531662.

Zhang L. N., Wang Y., Wu K., Sheng R.Y., Huang Q. L. Dynamic modeling and vibration characteristics of a twostage closedform planetary gear train. Mechanism and Machine Theory, Vol. 97, 2016, p. 1228.

Sun W., Li X., Wei J., Zhang A. Q., Ding X., Hu X. L. A study on loadsharing structure of multistage planetary transmission system. Journal of Mechanical Science and Technology, Vol. 29, Issue 4, 2015, p. 15011511.

Sun W., Ding X., Wei J., Hu X. L., Wang Q.G. An analyzing method of coupled modes in multistage planetary gear system. International Journal of Precision Engineering and Manufacturing, Vol. 15, Issue 11, 2014, p. 23572366.

Sheng L. C., Li W., Wang Y. Q., Fan M. B., Yang X. F. Nonlinear dynamic analysis and chaos control of multifreedom semidirect gear drive system in coal cutters. Mechanical Systems and Signal Processing, Vol. 1161, 2019, p. 6277.

Wang X. Stability research of the planetary multistage gear transmission system with coupling faults. Journal of Sound and Vibration, Vol. 434, 2018, p. 6377.

Wang X. A study on coupling faults’ characteristics of fixedaxis gear crack and planetary gear wear. Shock and Vibration, Vol. 2018, 2018, p. 4692796.
About this article
This research is supported by Natural Science Research Project of Shaanxi Provincial Department of Education (Program No. 19JK0912) and Natural Science Basic Research Program of Shaanxi (Program No. 2019JQ898), China.