In my work, I studied the dynamic behavior of a helical gear transmission system used in a high-speed train gearbox. I focused on two common local faults: tooth root cracks and tooth surface spalling. I considered the main internal excitations of the helical gear pair, including time-varying mesh stiffness, backlash, mesh damping, and transmission error. I built a six-degree-of-freedom nonlinear dynamic model, derived the time-varying mesh stiffness of a healthy helical gear pair, extended the stiffness formulation to cracked and spalled helical gear teeth, and then solved the dynamic response under different speeds and fault levels. The main goal of my study was to establish a clear link between fault parameters and fault features so that condition monitoring and fault diagnosis of helical gear systems can be supported by a reliable dynamic model.
Nonlinear factors in the helical gear system
The helical gear pair is not a linear system. Its vibration is driven by several internal factors that vary with meshing position and time. I treated time-varying mesh stiffness as the most important factor because it directly reflects tooth contact variation, crack growth, and spalling damage. I also included backlash because it controls the transition among normal meshing, single-sided impact, and double-sided impact. I included mesh damping because it dissipates vibration energy and influences the amplitude of the dynamic mesh force. I included transmission error because it represents the difference between the ideal and actual meshing positions.
The dynamic transmission error of the helical gear pair can be written as
$$
\delta(t) = (y_1 – y_2)\cos\beta_b + (z_1 – z_2)\sin\beta_b + r_{b1}\theta_1 – r_{b2}\theta_2 – e(t)
$$
where \(y_1, y_2\) are the transverse displacements, \(z_1, z_2\) are the axial displacements, \(\theta_1, \theta_2\) are the torsional displacements, \(r_{b1}, r_{b2}\) are the base circle radii, \(\beta_b\) is the base helix angle, and \(e(t)\) is the transmission error. The backlash function is
$$
f(\delta) =
\begin{cases}
\delta – b, & \delta \gt b \\
0, & -b \le \delta \le b \\
\delta + b, & \delta \lt -b
\end{cases}
$$
where \(b\) is half of the backlash. When \(f(\delta) = 0\), the helical gear pair is in the loss-of-contact state. When \(f(\delta) \gt 0\) or \(f(\delta) \lt 0\), the teeth are in driving-side or back-side contact. The dynamic mesh force is
$$
F_d = k(t) f(\delta) + c_m \dot{\delta}
$$
The mesh damping coefficient is estimated from the mean mesh stiffness and the equivalent inertia of the helical gear pair:
$$
c_m = 2\xi_g \sqrt{\frac{k_m r_{b1}^2 r_{b2}^2 I_1 I_2}{r_{b1}^2 I_1 + r_{b2}^2 I_2}}
$$
where \(\xi_g\) is the mesh damping ratio, \(k_m\) is the mean mesh stiffness, and \(I_1, I_2\) are the rotary inertias. The transmission error is represented by a Fourier series. For the dynamic analysis, I retained the first harmonic term:
$$
e(t) = e_0 + e_1 \cos(\omega_m t + \varphi_1)
$$
where \(e_0\) is the mean transmission error, \(e_1\) is the first harmonic amplitude, \(\omega_m\) is the mesh frequency, and \(\varphi_1\) is the phase angle.
Six-degree-of-freedom dynamic model of the helical gear system
I used the lumped-mass method and Newton’s second law to build a six-degree-of-freedom helical gear model. The model includes the transverse, axial, and torsional motions of the driving and driven gears. I made the following assumptions: friction on the tooth surface was neglected, the shaft mass was lumped into the gears, the gear bodies were treated as rigid, and the gear pair interaction was represented by a stiffness-damping element. The equations of motion are
$$
\begin{aligned}
I_1 \ddot{\theta}_1 + r_{b1} F_d &= T_1 \\
I_2 \ddot{\theta}_2 – r_{b2} F_d &= -T_2 \\
m_1 \ddot{y}_1 + c_{1y}\dot{y}_1 + k_{1y} y_1 &= F_y \\
m_2 \ddot{y}_2 + c_{2y}\dot{y}_2 + k_{2y} y_2 &= -F_y \\
m_1 \ddot{z}_1 + c_{1z}\dot{z}_1 + k_{1z} z_1 &= F_z \\
m_2 \ddot{z}_2 + c_{2z}\dot{z}_2 + k_{2z} z_2 &= -F_z
\end{aligned}
$$
where \(m_1, m_2\) are the gear masses, \(I_1, I_2\) are the rotary inertias, \(T_1, T_2\) are the torques, \(k_{1y}, k_{2y}, k_{1z}, k_{2z}\) are the support stiffnesses, and \(c_{1y}, c_{2y}, c_{1z}, c_{2z}\) are the support damping coefficients. The dynamic mesh force components along the transverse and axial directions are
$$
F_y = F_d \cos\beta_b, \qquad F_z = F_d \sin\beta_b
$$
To improve numerical stability, I introduced a dimensionless time \(\tau = \omega_n t\), a dimensionless mesh frequency ratio \(\Omega_h = \omega_m/\omega_n\), and dimensionless displacements \(Y_1 = y_1/b\), \(Y_2 = y_2/b\), \(Z_1 = z_1/b\), \(Z_2 = z_2/b\), and \(\lambda = \delta/b\). The dimensionless state equations were then solved by the standard fourth-order Runge–Kutta method.
| Parameter | Symbol | Driving gear | Driven gear |
|---|---|---|---|
| Number of teeth | \(z\) | 29 | 69 |
| Normal module (mm) | \(m_n\) | 7 | 7 |
| Normal pressure angle (deg) | \(\alpha_n\) | 26 | 26 |
| Helix angle (deg) | \(\beta\) | 20 | 20 |
| Face width (mm) | \(B\) | 75 | 70 |
| Mass (kg) | \(m\) | 9.65 | 56.92 |
| Support stiffness (N/m) | \(k_y, k_z\) | \(2.5 \times 10^7\) | \(2.5 \times 10^7\) |
| Support damping (N s/m) | \(c_y, c_z\) | \(1 \times 10^4\) | \(1 \times 10^4\) |
| Backlash half-width (m) | \(b\) | \(40 \times 10^{-6}\) | |
| Transmission error amplitude (m) | \(e_1\) | \(10 \times 10^{-6}\) | |
| Mesh damping ratio | \(\xi_g\) | 0.03 | |
Time-varying mesh stiffness of a healthy helical gear pair
The mesh stiffness of a helical gear pair changes along the contact line because the number of engaged tooth pairs changes and because each tooth pair moves from root to tip. I used slice theory and the potential energy method to calculate the time-varying mesh stiffness. In slice theory, the helical gear tooth is divided into a series of thin slices along the face width. Each slice is treated as a spur gear tooth with a small thickness. The total mesh stiffness is obtained by summing the stiffness of all slices that are in contact at a given meshing position.
For one slice, the tooth is treated as a cantilever beam fixed at the root circle. The potential energy includes bending energy, shear energy, and axial compression energy:
$$
U_a = \int_0^d \frac{F_a^2}{2 E A_x} dx, \qquad
U_b = \int_0^d \frac{F_b^2}{2 E I_x} dx, \qquad
U_s = \int_0^d \frac{1.2 F_b^2}{2 G A_x} dx
$$
where \(F_a\) is the axial force, \(F_b\) is the tangential force, \(E\) is Young’s modulus, \(G\) is the shear modulus, \(A_x\) is the cross-sectional area, and \(I_x\) is the area moment of inertia. The corresponding stiffness components are
$$
\frac{1}{k_a} = \int_0^d \frac{\sin^2\alpha_1}{E A_x} dx
$$
$$
\frac{1}{k_b} = \int_0^d \frac{[(d-x)\cos\alpha_1 – h\sin\alpha_1]^2}{E I_x} dx
$$
$$
\frac{1}{k_s} = \int_0^d \frac{1.2\cos^2\alpha_1}{G A_x} dx
$$
where \(\alpha_1\) is the mesh angle, \(d\) is the horizontal distance from the mesh point to the root, and \(h\) is the distance from the mesh point to the tooth centerline. The tooth profile was not simplified as a straight line. I used the actual fillet transition curve and the involute curve to build an accurate tooth profile. This step is important because the root circle and the base circle do not always coincide in a real helical gear. Ignoring the transition curve introduces a noticeable error in the mesh stiffness.
The single tooth pair stiffness is obtained by combining the bending, shear, axial compression, fillet foundation, and Hertzian contact stiffnesses in series:
$$
\frac{1}{k} =
\frac{1}{k_{b1}} + \frac{1}{k_{s1}} + \frac{1}{k_{a1}} +
\frac{1}{k_{f1}} + \frac{1}{k_h} +
\frac{1}{k_{f2}} + \frac{1}{k_{a2}} + \frac{1}{k_{s2}} + \frac{1}{k_{b2}}
$$
The Hertzian contact stiffness is
$$
k_h = \frac{E}{4(1-\nu^2)}
$$
For the helical gear pair, the total mesh stiffness at a given meshing position is
$$
K_z = \sum_{i=1}^{n} K_i
$$
where \(n\) is the number of tooth pairs in contact. Because the helical gear has an overlap ratio greater than two but less than three, the contact alternates between two-tooth and three-tooth meshing regions. This alternation produces a periodic variation in the total mesh stiffness.
I validated the analytical mesh stiffness model by the finite element method. The average relative error was 2.58%, and the maximum relative error was 4.52%. The calculated stiffness was slightly lower than the finite element result because the slice model does not fully account for shear interaction between adjacent slices. Nevertheless, the error is acceptable for dynamic analysis. I also compared linear interpolation and Fourier series fitting for the discrete stiffness data. The linear interpolation method gave a much smaller fitting error. The Fourier series method produced visible errors near stiffness jumps, with a maximum error of about 1.21%. Therefore, I used linear interpolation in the subsequent dynamic simulations.
| Method | Average error | Maximum error | Fitting quality |
|---|---|---|---|
| Potential energy plus slice theory versus finite element | 2.58% | 4.52% | Acceptable |
| Linear interpolation | Near zero | Very small | Excellent |
| Fourier series, 8th order | Moderate | 1.21% | Poor near jumps |
Tooth root crack model for the helical gear
I modeled a root crack that does not propagate uniformly along the face width. In the early and middle stages of crack growth, the crack depth varies from one end of the tooth to the other, and the crack may not pass through the entire face width. I therefore used a non-uniform crack depth function. For a through crack, the depth is
$$
q(z) = q_0 + (q_e – q_0)\frac{z^2}{L^2}
$$
For a non-through crack, the depth is
$$
q(z) =
\begin{cases}
q_0 – \frac{q_0}{L_c}z, & 0 \le z \le L_c \\
0, & L_c \lt z \le L
\end{cases}
$$
where \(q_0\) is the crack depth at the front face, \(q_e\) is the crack depth at the back face, \(L\) is the face width, and \(L_c\) is the crack length along the face width. When the crack is present, the cross-sectional area and area moment of inertia of the cracked slice are reduced. I divided the cracked slice into two cases. In the first case, the crack does not pass through the foundation. In the second case, the crack passes through the foundation. The reduced area moment of inertia and area are
$$
dI_{x,\text{crack}} =
\begin{cases}
\frac{1}{12} h_{c1}^3 \Delta y, & d_0 \gt 0, \ x \lt g_c \\
\frac{1}{12} (h_x + h_{c1})^3 \Delta y, & d_0 \gt 0, \ x \ge g_c
\end{cases}
$$
$$
dA_{x,\text{crack}} =
\begin{cases}
h_{c1} \Delta y, & d_0 \gt 0, \ x \lt g_c \\
(h_x + h_{c1}) \Delta y, & d_0 \gt 0, \ x \ge g_c
\end{cases}
$$
where \(h_{c1}\) is the reduced tooth height caused by the crack, \(d_0\) is the crack tip position relative to the root, and \(g_c\) is the transition position. These expressions replace the healthy cross-section terms in the bending and shear stiffness integrals. The axial compression stiffness, foundation stiffness, and contact stiffness are not changed by the crack because the crack mainly affects the bending and shear load paths.
For a cracked helical gear slice, the bending and shear stiffnesses are integrated over the cracked geometry. The total cracked mesh stiffness is then obtained by summing the stiffnesses of all slices along the face width, including the cracked slices and the healthy slices.
| Crack case | Front depth \(q_0\) (mm) | Back depth \(q_e\) (mm) | Crack length \(L_c\) (mm) | Crack percentage |
|---|---|---|---|---|
| Case 1 | 0 | 0 | 0 | 0% |
| Case 2 | 4 | 0 | 45 | 13% |
| Case 3 | 6 | 0 | 60 | 27% |
| Case 4 | 8 | 0 | 75 | 44% |
| Case 5 | 10 | 3 | 75 | 60% |
| Case 6 | 12 | 6 | 75 | 78% |
I defined the stiffness reduction percentage as
$$
P = \frac{k_0 – k_c}{k_0} \times 100\%
$$
where \(k_0\) is the healthy mesh stiffness and \(k_c\) is the cracked mesh stiffness. As the crack grows, the mesh stiffness decreases, and the reduction becomes larger. For the same crack, the stiffness reduction first increases and then decreases along the meshing path. The largest reduction occurs in the double-tooth meshing region near the tooth tip. This region is important because the local stiffness is already lower than in the three-tooth region, and the cracked slices are fully engaged there. In the most severe case I studied, the maximum stiffness reduction reached about 10.85%. This result indicates that for a cracked helical gear, the double-tooth region near the tip should be carefully monitored.
Tooth surface spalling model for the helical gear
I modeled spalling as a rectangular pit on the tooth surface. The spall is described by four parameters: spall length \(l_s\), spall width \(w_s\), spall depth \(h_s\), and spall starting position \(x_{\text{start}}\). The length and width are parallel to the helix line and the transverse tooth profile, respectively. I assumed that the spall is symmetric about the face-width centerline, which is a reasonable approximation for typical pitting and spalling near the pitch line. The spall affects the mesh stiffness in two ways. First, it reduces the cross-sectional area and area moment of inertia of the affected slices. Second, it reduces the effective contact line length when the meshing point enters the spall region.
For a spalled slice, the cross-sectional area and area moment of inertia are
$$
dA_x =
\begin{cases}
h_x \Delta y, & \text{no spall} \\
(h_x – h_s)\Delta y, & \text{spall region}
\end{cases}
$$
$$
dI_x =
\begin{cases}
\frac{1}{12} h_x^3 \Delta y, & \text{no spall} \\
\frac{1}{12}(h_x – h_s)^3 \Delta y + A_s \delta_{xs}^2, & \text{spall region}
\end{cases}
$$
where \(A_s\) is the remaining area of the spalled slice and \(\delta_{xs}\) is the shift of the neutral axis caused by the loss of material. The neutral axis shift is
$$
\delta_{xs} = \frac{0.5 h_s (h_x – h_s)}{h_x}
$$
The contact line length in the spall region is
$$
dL_x =
\begin{cases}
0, & x_{\text{start}} \le x \le x_{\text{end}} \\
\Delta y \cos\beta_b, & \text{otherwise}
\end{cases}
$$
The contact stiffness of a spalled slice is therefore
$$
k_h = \frac{E dL_x}{4(1-\nu^2)}
$$
When the meshing point is outside the spall region, the spall does not affect the contact stiffness. When the meshing point enters the spall region, the contact stiffness drops sharply. The bending, shear, and axial compression stiffnesses are also reduced because the cross-section is smaller. However, the contact stiffness reduction is usually more localized and more abrupt. I summed the slice stiffnesses to obtain the total spalled mesh stiffness.
| Spall parameter | Symbol | Effect on mesh stiffness | Relative importance |
|---|---|---|---|
| Spall depth | \(h_s\) | Changes the local cross-section and inertia | Moderate |
| Spall width | \(w_s\) | Changes the spall-affected width | Moderate |
| Spall starting position | \(x_{\text{start}}\) | Shifts the affected meshing interval | Moderate |
| Spall length | \(l_s\) | Strongly changes the affected face width and contact loss | Strongest |
Among the four spall parameters, the spall length had the strongest influence on the mesh stiffness. As the spall length increases, the affected face-width region becomes wider, and the contact line loss becomes more severe. The mesh stiffness reduction forms a local dip. Along the meshing path, the stiffness first decreases rapidly, then decreases slowly, then increases rapidly, and finally returns to the healthy value. The spall depth and spall width do not change the extent of the spall-affected meshing interval as much as the spall length does. Moving the spall starting position toward the tooth tip shifts the affected interval and reduces its width.
Dynamic behavior of the healthy helical gear system
I first solved the healthy helical gear system to establish a baseline. I used the linearly interpolated mesh stiffness and the standard fourth-order Runge–Kutta method. I treated the driving gear speed as the bifurcation parameter and studied the speed range from 1000 r/min to 10000 r/min. The healthy system showed rich nonlinear behavior. At low speeds, the system was in period-one motion. Near 4700 r/min, a jump appeared in the bifurcation diagram. After the jump, the system remained period-one, but the vibration amplitude changed. As the speed increased further, the system entered period-two motion, then chaotic motion, then period-four motion, and then chaotic motion again. Several periodic windows appeared inside the chaotic range, including period-five motion.
| Speed range (r/min) | Dominant motion | Poincaré map feature | Frequency feature |
|---|---|---|---|
| 1000–6300 | Period-one and some quasi-periodic intervals | One point or a small point cluster | Mesh frequency and harmonics |
| 6300–7250 | Period-two and chaos | Two points or many points | Half mesh frequency or broadband |
| 7250–7750 | Period-two to period-four | Two to four points | Subharmonic mesh frequency |
| 7750–10000 | Chaos with periodic windows | Many points or structured sets | Broadband with narrow peaks |
The healthy system showed a transition from no impact to single-sided impact and then to double-sided impact as speed increased. I evaluated the dynamic mesh force. The maximum mesh force remained positive for all speeds. The minimum mesh force changed from positive to negative as speed increased. When both maximum and minimum mesh forces were positive, the helical gear pair operated without impact. When the minimum mesh force became negative while the maximum remained positive, the system entered double-sided impact. The boundary between these states was near 4700 r/min.
Effect of a root crack on the helical gear dynamics
I introduced the cracked mesh stiffness into the dynamic model and solved the system response over the same speed range. The crack fault changed the nonlinear behavior of the helical gear system. At low speeds, the healthy system was period-one, but the cracked system became period-three. The time-domain signal showed periodic impacts with a period equal to the rotation period of the cracked gear. The frequency spectrum showed fault characteristic frequencies around the mesh frequency and its harmonics. The fault characteristic frequency was equal to the rotation frequency of the cracked gear.
At medium speeds, the cracked system became quasi-periodic. The phase trajectory changed from a single closed curve to multiple closed curves. The Poincaré map changed from a finite set of points to a point cluster. The time-domain impact features became less obvious because the impact energy was distributed over a wider frequency band. At high speeds, the cracked system entered chaos at a lower speed than the healthy system. The maximum Lyapunov exponent increased and became positive in some speed ranges. This indicates that the crack reduced the stability of the helical gear system.
| Speed condition | Healthy helical gear | Cracked helical gear | Main change |
|---|---|---|---|
| Low speed, e.g. 3500 r/min | Period-one | Period-three | Periodic impacts and fault frequency |
| Medium speed, e.g. 6500 r/min | Period-two | Quasi-periodic | Point cluster and sidebands |
| High speed, e.g. 9900 r/min | Period-three or periodic window | Chaos | Broadband spectrum and positive Lyapunov exponent |
The crack evolution study showed that the vibration amplitude increases with crack severity. At 5400 r/min, the system remained relatively stable when the crack percentage was below about 22%. When the crack percentage exceeded 22%, the number of periodic motions increased, and the vibration amplitude increased rapidly. The time-domain impact sequence became more regular, and the fault frequency amplitude in the spectrum increased. Therefore, for this helical gear system, a crack percentage of about 22% can be regarded as an early warning threshold at 5400 r/min.
At 7600 r/min, the healthy system was period-four. As the crack grew, the system changed from period-four to quasi-period-four and then to chaos. The phase trajectory expanded from two closed curves to two broad bands. The Poincaré map changed from four points to four point clusters and then to four structured line sets. The spectrum developed a continuous component. The fault feature frequency became harder to identify at high crack levels because the broadband component masked the discrete peaks.
Effect of tooth surface spalling on the helical gear dynamics
I also introduced the spalled mesh stiffness into the dynamic model. The spalling fault produced effects similar to the crack fault, but the severity was generally lower for the same fault level. At low speed, spalling changed the system from period-one to period-two. The time-domain signal showed periodic impacts with a period equal to the rotation period of the spalled gear. The spectrum showed sidebands around the mesh frequency. At medium speed, spalling changed period-four motion into quasi-period-four motion. At high speed, spalling did not change the chaotic state, but it increased the amplitude of some subharmonic components.
The maximum Lyapunov exponent increased when spalling was present, but the increase was smaller than that caused by a crack. This suggests that a root crack has a stronger influence on the stability of the helical gear system than a surface spall of comparable size. The dynamic mesh force analysis showed that spalling did not change the impact state boundaries significantly. At speeds below about 2500 r/min, spalling increased the maximum mesh force. At higher speeds, the maximum mesh force remained close to the healthy value.
| Speed condition | Healthy helical gear | Spalled helical gear | Main change |
|---|---|---|---|
| Low speed, e.g. 3000 r/min | Period-one | Period-two | Periodic impacts and sidebands |
| Medium speed, e.g. 7500 r/min | Period-four | Quasi-period-four | Point clusters and broadened spectral peaks |
| High speed, e.g. 9500 r/min | Chaos | Chaos | Increased subharmonic amplitude |
I studied spall evolution by increasing the spall length from 0 to 70 mm. At 6500 r/min, the healthy system was period-two. As the spall length increased, the system changed to quasi-period-two and then to multi-periodic motion. The phase trajectory expanded from two closed curves to two broad bands. The Poincaré map changed from two points to two point clusters and then to two line-like structures. The time-domain impact sequence became more obvious, and the fault frequency amplitude in the spectrum increased. These features made the spalling fault easier to identify at low and medium speeds.
At 8500 r/min, the healthy system was period-five. As the spall length increased, the system entered quasi-period-five, then chaos at about 43 mm spall length, and then alternated between quasi-periodic and chaotic states. The phase trajectory expanded into broad bands. The Poincaré map changed from five points to five point clusters and then to five line-like structures. The spectrum developed a continuous component at large spall lengths. The fault frequency became less identifiable when the spall length was very large because the chaotic broadband component dominated the spectrum.
Speed dependence and fault feature identifiability
The dynamic response of the faulty helical gear system depends strongly on speed. At low speeds, the fault features are clearer in both the time domain and the frequency domain. The time-domain signal shows periodic impacts, and the frequency spectrum shows discrete fault characteristic frequencies. At high speeds, the nonlinearity becomes stronger. The system tends to enter quasi-periodic or chaotic motion. The time-domain impact sequence becomes irregular, and the spectrum contains a continuous component. As a result, the fault characteristic frequency is harder to identify at high speeds. Therefore, low-speed operating conditions are more favorable for early fault detection in helical gear systems.
I also observed that the maximum Lyapunov exponent is a useful indicator of stability. For the healthy helical gear system, the maximum Lyapunov exponent was usually negative. For the cracked and spalled systems, it increased. In some high-speed cases, it became positive, indicating chaotic motion. The crack caused a larger increase in the maximum Lyapunov exponent than the spall. This means that a root crack is more likely to destabilize the helical gear system than a surface spall of similar size.
| Fault type | Low speed effect | Medium speed effect | High speed effect | Stability impact |
|---|---|---|---|---|
| Root crack | Period-one to period-three | Periodic to quasi-periodic | Periodic to chaos at lower speed | Strong |
| Surface spall | Period-one to period-two | Period-four to quasi-period-four | Chaos remains chaos | Moderate |
Time-domain statistical indicators for fault diagnosis
To evaluate the health condition of the helical gear system, I calculated several time-domain statistical indicators from the simulated vibration signals. I used kurtosis, root mean square, waveform factor, peak factor, and impulse factor. The formulas are
$$
K = \frac{\frac{1}{N}\sum_{i=1}^N (x_i – \bar{x})^4}{\left[\frac{1}{N}\sum_{i=1}^N (x_i – \bar{x})^2\right]^2}
$$
$$
RMS = \sqrt{\frac{1}{N}\sum_{i=1}^N (x_i – \bar{x})^2}
$$
$$
W_f = \frac{RMS}{\frac{1}{N}\sum_{i=1}^N |x_i|}
$$
$$
P_f = \frac{x_p}{RMS}
$$
$$
I_f = \frac{x_p}{\frac{1}{N}\sum_{i=1}^N |x_i|}
$$
where \(x_i\) is the vibration signal, \(x_p\) is the peak value, \(N\) is the number of samples, and \(\bar{x}\) is the mean value. I used two methods. In method one, I subtracted the healthy signal from the faulty signal and then calculated the statistical indicators of the residual signal. In method two, I calculated the indicators of the healthy and faulty signals separately and then computed the percentage change:
$$
R_{Xi} = \frac{X_i – X_0}{X_0} \times 100\%
$$
where \(X_i\) is the indicator for the \(i\)-th fault level and \(X_0\) is the indicator for the healthy gear. The results showed that method one gave better sensitivity. For the crack fault, the waveform factor and peak factor increased monotonically with crack severity. The other indicators, such as kurtosis and RMS, did not show a clear monotonic trend. For the spalling fault, the waveform factor and peak factor also increased monotonically with spall length. The increase became faster when the spall length exceeded about 10 mm. Therefore, the waveform factor and peak factor calculated by method one can be used as effective indicators of fault severity in helical gear systems.
| Indicator | Crack sensitivity | Spall sensitivity | Monotonic trend | Recommended use |
|---|---|---|---|---|
| Kurtosis | Low | Low to moderate | Not clear | Supplementary |
| RMS | Low | Low | Not clear | Not preferred |
| Waveform factor | High | High | Yes | Recommended |
| Peak factor | High | High | Yes | Recommended |
| Impulse factor | Low | Low | Not clear | Not preferred |
Summary of findings
My study shows that the dynamic behavior of a helical gear system is strongly affected by root cracks and surface spalling. The mesh stiffness of a healthy helical gear pair varies periodically because of the change in the number of contacting tooth pairs and the movement of the contact line. An accurate tooth profile, including the root transition curve, is necessary for correct stiffness calculation. The slice theory and potential energy method provide a fast and reasonably accurate way to calculate the time-varying mesh stiffness of a helical gear pair. The finite element method confirmed that the analytical model is reliable, with an average error of 2.58% and a maximum error of 4.52%.
A root crack reduces the mesh stiffness. The reduction depends on the crack depth, crack length, and crack position. For a non-uniform crack, the largest stiffness reduction occurs near the tooth tip in the double-tooth meshing region. A surface spall also reduces the mesh stiffness, but the effect is more localized. Among the spall parameters, the spall length has the strongest influence. The spall length changes both the affected face-width region and the contact line length.
The dynamic response of the helical gear system is highly nonlinear. The healthy system can transition from period-one motion to period-two, period-four, period-five, quasi-periodic, and chaotic motion as speed increases. A root crack and a surface spall both change these transitions. At low speeds, the faults introduce period-three or period-two motion and produce clear impact features in the time domain. At high speeds, the faults promote quasi-periodic or chaotic motion and produce broadband spectral components. The fault features are more identifiable at low speeds. The maximum Lyapunov exponent increases when a fault is present, and a root crack causes a larger increase than a surface spall.
The time-domain statistical analysis shows that the waveform factor and peak factor are the most sensitive indicators for both crack and spall faults. They increase monotonically with fault severity when the residual signal method is used. These two indicators can be used for condition monitoring and fault severity assessment of helical gear systems. The results of my study provide a theoretical basis for fault diagnosis and health management of helical gear transmission systems, especially for high-speed train gearboxes operating under variable speed conditions.
Practical implications for helical gear condition monitoring
The practical implications of my study are important for the safe operation of helical gear systems. First, the mesh stiffness model should include the actual tooth profile and the transition curve. A simplified root model may underestimate or overestimate the stiffness reduction caused by a crack. Second, the double-tooth meshing region near the tooth tip should be monitored closely because it experiences the largest stiffness reduction when a root crack is present. Third, low-speed operating conditions should be used for early fault detection whenever possible. At low speeds, the fault-induced impacts are clearer in the time domain, and the fault characteristic frequencies are more distinct in the frequency domain. Fourth, the maximum Lyapunov exponent can be used as an auxiliary indicator of stability. A transition from negative to positive values indicates that the helical gear system has entered a chaotic state.
In future work, I plan to improve the mesh stiffness model by considering the coupling between adjacent slices. I also plan to include friction and lubrication effects, because these factors influence the contact condition and the spalling evolution. Finally, I intend to build a full helical gear-rotor-bearing model to represent the whole transmission system more accurately. These improvements will make the dynamic simulation of faulty helical gear systems more realistic and will further support fault diagnosis and condition monitoring in practical engineering applications.