In this study, we investigate the dynamic behavior of high contact ratio spur gears, specifically cylindrical gears, under the influence of tooth surface wear. Cylindrical gears are widely used in mechanical systems due to their efficiency and reliability, but wear over time can significantly alter their performance. We focus on calculating wear-induced changes in mesh stiffness and analyzing how these changes affect the gear system’s dynamics. Our approach integrates the Archard wear model with energy-based stiffness calculations to derive a comprehensive mesh stiffness function, which is then incorporated into a single-degree-of-freedom dynamic model. We explore both light and heavy load conditions to assess the sensitivity of cylindrical gear dynamics to wear. Throughout this article, we emphasize the term “cylindrical gear” to highlight the gear type under consideration, and we present our findings using extensive tables, formulas, and analyses to ensure clarity and depth.
The importance of cylindrical gears in industrial applications cannot be overstated, especially those with high contact ratios (typically greater than 2), which offer reduced noise and improved dynamic characteristics. However, surface wear from prolonged operation can degrade these benefits, leading to unexpected vibrations, noise, and even failure. Prior research has addressed wear in standard cylindrical gears, but less attention has been given to high contact ratio variants. We aim to fill this gap by developing a detailed model that accounts for non-uniform wear patterns and their impact on mesh stiffness and dynamic response. Our work builds on existing methods but extends them to handle the unique geometry of high contact ratio cylindrical gears, providing insights that can aid in design and maintenance decisions.
We begin by establishing the dynamic model for a cylindrical gear system. The single-degree-of-freedom equation of motion is expressed as:
$$ m_e \ddot{x}(t) + c \dot{x}(t) + k(t) f(x) = F_m + F_h(t) $$
Here, \( m_e \) is the equivalent mass of the system, \( c \) is the damping coefficient, \( k(t) \) is the time-varying mesh stiffness, \( f(x) \) is the backlash function, \( F_m \) is the equivalent external excitation, and \( F_h(t) \) is the internal excitation given by \( F_h(t) = m_e f_h \omega_h^2 \cos(\omega_h t + \phi_h) \), where \( f_h \) is the amplitude, \( \omega_h \) is the frequency, and \( \phi_h \) is the phase. For cylindrical gears, this model captures the essential dynamics, but we must incorporate wear effects into \( k(t) \) to accurately predict behavior.
To simplify analysis, we non-dimensionalize the equation. Define the natural frequency \( \omega_n = \sqrt{K_m / m_e} \), where \( K_m \) is the average mesh stiffness. Let \( \tau = \omega_n t \) be the non-dimensional time, \( l \) be a nominal length, and \( \bar{x}(\tau) = x(t)/l \) be the non-dimensional displacement. The non-dimensional form becomes:
$$ \bar{x}” + 2\xi \bar{x}’ + \bar{K}(\tau) \bar{f}(\bar{x}) = \bar{F}_m + \bar{F}_h(\tau) $$
where \( \xi = c / (2 \sqrt{K_m m_e}) \), \( \bar{K}(\tau) = k(t)/K_m \), \( \bar{F}_m = F_m/(K_m l) \), \( \bar{F}_h(\tau) = \bar{\omega}_h^2 \cos(\bar{\omega}_h \tau + \psi) \) with \( \bar{\omega}_h = \omega_h / \omega_n \), and \( \bar{f}(\bar{x}) \) is the non-dimensional backlash function. This formulation allows us to study cylindrical gear dynamics across different operating conditions.
Next, we introduce wear characteristics into the model. For cylindrical gears, the time-varying mesh stiffness per tooth pair is computed using the energy method, which considers bending, shear, axial compression, contact, and fillet foundation stiffnesses. The stiffness components are derived from the gear geometry, which for high contact ratio cylindrical gears involves a modified tooth profile due to manufacturing with a hob cutter. The bending stiffness \( k_b \), shear stiffness \( k_s \), and axial stiffness \( k_a \) for a single tooth are given by:
$$ \frac{F^2}{2k_b} = \int_{x_G}^{x_E} \frac{M_1^2}{2EI_{x1}} dx_1 + \int_{x_E}^{x_D} \frac{M_2^2}{2EI_{x2}} dx_2 + \int_{x_D}^{x_C} \frac{M_3^2}{2EI_{x3}} dx_3 + \int_{x_C}^{x_K} \frac{M_4^2}{2EI_{x4}} dx_4 $$
$$ \frac{F^2}{2k_s} = \int_{x_G}^{x_E} \frac{1.2 (F \cos \beta)^2}{2GA_{x1}} dx_1 + \int_{x_E}^{x_D} \frac{1.2 (F \cos \beta)^2}{2GA_{x2}} dx_2 + \int_{x_D}^{x_C} \frac{1.2 (F \cos \beta)^2}{2GA_{x3}} dx_3 + \int_{x_C}^{x_K} \frac{1.2 (F \cos \beta)^2}{2GA_{x4}} dx_4 $$
$$ \frac{F^2}{2k_a} = \int_{x_G}^{x_E} \frac{(F \sin \beta)^2}{2EA_{x1}} dx_1 + \int_{x_E}^{x_D} \frac{(F \sin \beta)^2}{2EA_{x2}} dx_2 + \int_{x_D}^{x_C} \frac{(F \sin \beta)^2}{2EA_{x3}} dx_3 + \int_{x_C}^{x_K} \frac{(F \sin \beta)^2}{2EA_{x4}} dx_4 $$
In these integrals, \( x_G, x_E, x_D, x_C, x_K \) are coordinates along the tooth profile, \( E \) is Young’s modulus, \( G = E/(2(1+\nu)) \) is the shear modulus with Poisson’s ratio \( \nu \), \( \beta = \alpha_K – \theta_K \) where \( \alpha_K \) is the pressure angle at point K, \( \theta_K = \pi/(4z) – \text{inv} \alpha_K + \text{inv} \alpha_0 \), \( z \) is the number of teeth, and \( \alpha_0 \) is the standard pressure angle. The terms \( M_i \), \( I_{xi} \), and \( A_{xi} \) represent moment, area moment of inertia, and cross-sectional area for segments i=1 to 4. The contact stiffness \( k_h \) and fillet foundation stiffness \( k_f \) are calculated based on Hertzian contact theory and empirical formulas, respectively. The single-tooth mesh stiffness \( k \) for a cylindrical gear pair is then:
$$ k = \frac{1}{\frac{1}{k_h} + \frac{1}{k_{b1}} + \frac{1}{k_{b2}} + \frac{1}{k_{s1}} + \frac{1}{k_{s2}} + \frac{1}{k_{a1}} + \frac{1}{k_{a2}} + \frac{1}{k_{f1}} + \frac{1}{k_{f2}}} $$
where subscripts 1 and 2 denote the driving and driven cylindrical gears. This stiffness is essential for understanding how cylindrical gears transmit load under varying conditions.
To account for wear, we use the Archard wear model to calculate tooth surface wear depth. For cylindrical gears under constant load, the wear depth \( h \) at a given meshing point is:
$$ h = 2a \lambda n t \varepsilon_\alpha I_h $$
Here, \( a \) is the half-width of contact, \( \lambda \) is the sliding coefficient, \( n \) is rotational speed, \( t \) is operating time, \( \varepsilon_\alpha \) is the contact ratio, and \( I_h \) is the wear rate. The sliding coefficient and wear rate depend on gear geometry and material properties. We compute these for a cylindrical gear pair with parameters listed in Table 1. The wear distribution along the tooth profile is non-uniform, as shown in our calculations.
| Parameter | Driving Gear | Driven Gear |
|---|---|---|
| Number of Teeth, \( z \) | 32 | 25 |
| Module, \( m \) (mm) | 3.25 | 3.25 |
| Pressure Angle, \( \alpha_0 \) | 20° | 20° |
| Addendum Coefficient | 1.35 | 1.35 |
| Young’s Modulus, \( E \) (N/mm²) | 206,000 | |
| Poisson’s Ratio, \( \nu \) | 0.3 | |
| Face Width, \( L \) (mm) | 20 | |
| Driven Gear Speed (rpm) | 600 | |
| Torque (N·m) | 50 | |
Our wear calculations reveal that for cylindrical gears, wear depth varies with pressure angle, reaching a minimum at the pitch point. After \( N = 10^6 \) meshing cycles, the wear pattern indicates that initial wear has minimal impact on stiffness, but as cycles increase to \( 10^7 \), stiffness reductions become noticeable. This is critical for cylindrical gear longevity and performance.
We now derive the total mesh stiffness \( K(t) \) for the cylindrical gear pair, considering wear. For a high contact ratio cylindrical gear with \( \varepsilon = 2.23 \), the number of tooth pairs in contact alternates between 2 and 3. The total stiffness is the sum of individual tooth stiffnesses over the meshing period \( T \). We express it as a piecewise function:
$$ K(t) =
\begin{cases}
k(Y(t)) + k(Y(t) + T) + k(Y(t) + 2T), & 0 \leq Y(t) \leq (\varepsilon – 2)T \\
k(Y(t)) + k(Y(t) + T), & (\varepsilon – 2)T \leq Y(t) \leq T
\end{cases} $$
where \( Y(t) = \text{mod}(t, T) \), and \( k(t) \) is the single-tooth stiffness modified by wear. We fit \( k(t) \) to a continuous function for ease of computation. The effect of wear on \( K(t) \) is summarized in Table 2, which shows stiffness reduction percentages for different wear cycles.
| Meshing Cycles, \( N \) | Stiffness Reduction (%) | Notes |
|---|---|---|
| 0 (no wear) | 0 | Baseline for cylindrical gear |
| 2 × 10⁶ | 0.23 | Minor effect on cylindrical gear dynamics |
| 6 × 10⁶ | 0.72 | Moderate change in cylindrical gear response |
| 10⁷ | 1.26 | Significant impact on cylindrical gear behavior |
This stiffness function is inserted into the dynamic model to analyze system response. We solve the non-dimensional equation numerically for both light and heavy load conditions, with parameters \( \bar{f}_h = 0.1 \), \( \bar{\omega}_h = 1 \), \( \bar{b} = 0.5 \), and \( \xi = 0.05 \). For light load, we set \( \bar{F}_m = 0.04 \); for heavy load, \( \bar{F}_m = 0.08 \). The results illustrate how wear influences cylindrical gear dynamics.
Under light load without wear, the cylindrical gear system exhibits periodic motion, specifically a 5-times quasi-periodic response, as seen in time histories, Poincaré maps, and frequency spectra. However, with wear introduced at \( N = 2 \times 10^6 \) cycles, the system transitions to chaotic motion, characterized by irregular time series, scattered Poincaré points, and broadband noise in spectra. At \( N = 6 \times 10^6 \), the cylindrical gear shows near-periodic behavior, indicating a partial return to order, but at \( N = 10^7 \), chaos re-emerges. This complex progression—quasi-period → chaos → near-period → chaos—highlights the sensitivity of cylindrical gear dynamics to wear under light loads. The changes are attributed to stiffness variations that alter resonance conditions and nonlinear interactions.
In contrast, under heavy load, the cylindrical gear system remains predominantly periodic regardless of wear level. Even at \( N = 10^7 \) cycles, the dynamic characteristics—time history, Poincaré map, and spectrum—are nearly identical to the unworn case. This insensitivity stems from the dominant external load, which suppresses wear-induced stiffness changes, allowing the cylindrical gear to maintain stable operation. This finding is crucial for applications where cylindrical gears operate under high loads, as wear may have less immediate impact on performance.
To further quantify these effects, we analyze key dynamic indicators. Table 3 summarizes the response types for cylindrical gears under different wear and load conditions, based on our simulations.
| Load Condition | Wear Cycles, \( N \) | Response Type | Implication for Cylindrical Gear |
|---|---|---|---|
| Light Load (\( \bar{F}_m = 0.04 \)) | 0 | Quasi-Periodic (5×) | Stable cylindrical gear operation |
| 2 × 10⁶ | Chaotic | Increased vibration in cylindrical gear | |
| 6 × 10⁶ | Near-Periodic | Partial recovery in cylindrical gear dynamics | |
| 10⁷ | Chaotic | Severe instability in cylindrical gear | |
| Heavy Load (\( \bar{F}_m = 0.08 \)) | 0 | Periodic (2×) | Robust cylindrical gear performance |
| 2 × 10⁶ | Periodic (2×) | Minimal change in cylindrical gear behavior | |
| 6 × 10⁶ | Periodic (2×) | Cylindrical gear maintains stability | |
| 10⁷ | Periodic (2×) | Wear has negligible effect on cylindrical gear |
These results underscore the importance of considering wear in cylindrical gear design, especially for light-load applications where dynamics are more susceptible to degradation. The cylindrical gear’s high contact ratio amplifies these effects due to the alternating contact patterns, making wear management essential.
We also explore the mathematical underpinnings of wear progression. The Archard model can be extended to account for dynamic loads in cylindrical gears. The wear depth over time is governed by:
$$ \frac{dh}{dt} = 2a \lambda n \varepsilon_\alpha I_h $$
Integrating this for constant parameters gives the linear wear accumulation, but for cylindrical gears under varying conditions, numerical methods are required. We approximate the wear rate \( I_h \) using material properties and lubrication states, which for steel cylindrical gears typically ranges from \( 10^{-15} \) to \( 10^{-13} \) m²/N. This influences long-term predictions for cylindrical gear life.
Furthermore, the mesh stiffness function \( K(t) \) can be expressed in a Fourier series to facilitate analytical solutions. For a cylindrical gear with wear, we write:
$$ K(t) = K_m + \sum_{i=1}^{\infty} [A_i \cos(i\omega_m t) + B_i \sin(i\omega_m t)] $$
where \( \omega_m \) is the meshing frequency, and coefficients \( A_i, B_i \) depend on wear depth. This representation helps in studying parametric excitations in cylindrical gear systems. Wear reduces \( K_m \) and alters higher harmonics, potentially triggering instabilities.
In summary, our study provides a comprehensive framework for analyzing cylindrical gear dynamics with surface wear. We have derived formulas for stiffness and wear, presented tabular data, and simulated dynamic responses. The key takeaways are that cylindrical gears under light loads experience complex dynamic shifts due to wear, while under heavy loads, they remain stable. This has practical implications for industries relying on cylindrical gears, such as automotive and aerospace, where predictive maintenance can be informed by wear monitoring. Future work could incorporate thermal effects, lubrication variations, and multi-degree-of-freedom models to enhance accuracy for cylindrical gear applications. Throughout, we emphasize that cylindrical gear performance is intricately linked to wear evolution, necessitating continued research and innovation.
To conclude, cylindrical gears are fundamental components in machinery, and understanding their wear-induced dynamics is vital for reliability. Our integrated approach of modeling stiffness and wear offers a pathway to optimized cylindrical gear design and operation. We hope this contribution aids engineers in developing more durable and efficient cylindrical gear systems.