Dynamic Simulation and Diagnosis of Fatigue Pitting Fault in Helical Gears

Helical gears are critical components in various mechanical systems, including automotive transmissions, aerospace applications, and industrial machinery. They operate under high-speed conditions and are subjected to complex loads, leading to common failure modes such as tooth surface pitting. Pitting, a form of surface fatigue, occurs primarily due to repeated contact stresses and can significantly degrade gear performance and lifespan. Early detection and diagnosis of pitting faults are essential for preventing catastrophic failures and ensuring reliable operation. In this study, we focus on developing a dynamic simulation model to analyze the effects of single-tooth and multi-tooth pitting faults in helical gears, combining theoretical analysis with experimental validation. We explore the time-domain and frequency-domain characteristics of vibration signals to identify fault patterns, providing a foundation for predictive maintenance and fault diagnosis in gear transmission systems.

The dynamic behavior of helical gears under pitting faults is complex due to the helical tooth geometry, which introduces axial forces and varying contact conditions. Traditional studies often concentrate on spur gears, leaving a gap in understanding helical gear dynamics. Our approach involves creating detailed finite element and dynamic models to simulate pitting defects, followed by experimental verification. We employ envelope spectrum analysis to extract fault features from vibration signals, enabling the identification of pitting severity and the number of affected teeth. This comprehensive methodology enhances the accuracy of fault diagnosis and supports the development of health monitoring systems for helical gears in practical applications.

To establish the dynamic simulation model, we first define the geometric and material parameters of the helical gears. The gears are modeled with specific tooth numbers, module, pressure angle, helix angle, and face width, as summarized in Table 1. These parameters are derived from a typical 7-speed double-clutch automatic transmission first-stage helical gear pair, ensuring realism in our simulations. The material properties, including density, elastic modulus, and Poisson’s ratio, are based on high-strength alloy steel commonly used in automotive applications, as detailed in Table 2. The helical gears’ design involves a helix angle that promotes smooth engagement and higher load capacity compared to spur gears, but it also complicates the dynamic response due to the gradual tooth contact.

Table 1: Geometric Parameters of the Helical Gear Pair
Parameter Driving Gear Driven Gear
Number of Teeth 17 60
Normal Module (mm) 2.1 2.1
Pressure Angle (°) 17.5 17.5
Helix Angle (°) 29 29
Face Width (mm) 19.8 16.9
Table 2: Material Properties of the Helical Gears
Material Density (kg/m³) Elastic Modulus (GPa) Poisson’s Ratio
20MnCrS5 7840 210 0.278

The dynamic simulation model incorporates a multi-body dynamics approach, where the helical gears are represented as rigid bodies with rotational degrees of freedom. We apply a driving speed of 2500 rpm to the input shaft and a load torque of 840 N·m to the output shaft, replicating typical operating conditions. The contact between the helical gear teeth is modeled using a nonlinear spring-damper system, with the meshing stiffness calculated based on the potential energy method. The total meshing stiffness \( K \) for helical gears accounts for bending stiffness \( K_b \) and contact stiffness \( K_n \), derived from the following equations:

$$ \frac{1}{K_n} = \frac{4(1 – \nu^2)}{\pi E B} $$

$$ \frac{1}{K_b} = \int_0^d \frac{[(d – x) \cos \alpha – h \sin \alpha]^2}{E I_x} dx $$

$$ K = \frac{1}{\frac{1}{K_n} + \frac{1}{K_{b1}} + \frac{1}{K_{b2}}} $$

Here, \( E \) is the elastic modulus, \( \nu \) is Poisson’s ratio, \( B \) is the face width, \( I_x \) is the moment of inertia at a distance \( x \) from the base circle, \( d \) is the distance from the base circle along the tooth height, \( h \) is the distance from the contact point to the gear centerline, and \( \alpha \) is the angle between the contact force and the gear centerline. For helical gears, the integration is performed over the entire contact line, considering the helix angle effect. The meshing stiffness varies with time due to the changing number of tooth pairs in contact, which is characteristic of helical gears’ continuous engagement.

In the simulation, we define pitting faults as elliptical or rectangular surface defects on the tooth flanks, with dimensions and depths based on real-world observations. For single-tooth pitting, we consider minor pitting (e.g., 5 mm major axis, 1.5 mm minor axis) and severe pitting (e.g., 10 mm major axis, 1.5 mm minor axis), both with a depth of 0.2 mm. For multi-tooth pitting, we model defects on two or three teeth with varying sizes. The contact parameters, including meshing stiffness, collision index, damping coefficient, and penetration depth, are set as per Table 3 to ensure accurate dynamic interaction. The dynamic equations of motion are solved using numerical integration, and the resulting contact forces and vibrations are analyzed.

Table 3: Contact Parameters for Dynamic Simulation
Parameter Value
Meshing Stiffness (N/mm) 3.741 × 10⁵
Collision Index 1.5
Damping Coefficient 48
Penetration Depth (mm) 0.1

For single-tooth pitting faults, the time-domain analysis of contact forces reveals distinct patterns. Under normal conditions, the contact force remains relatively stable with minimal fluctuations. However, with minor pitting, the amplitude increases by approximately 2%, and with severe pitting, a 3% increase is observed along with periodic impacts at intervals of 0.024 seconds, corresponding to the rotational period of the faulty gear. In the frequency domain, we compute the envelope spectrum using the Hilbert transform, defined as:

$$ y(t) = \frac{1}{\pi} \int_{-\infty}^{\infty} \frac{x(\tau)}{t – \tau} d\tau $$

where \( x(\tau) \) is the original signal. The rotational frequency \( f_r \) and meshing frequency \( f_m \) are calculated as:

$$ f_r = \frac{n}{60} $$

$$ f_m = f_r z_i $$

with \( n = 2500 \) rpm and \( z_i \) as the number of teeth (17 for driving gear, 60 for driven gear), giving \( f_r = 41.67 \) Hz and \( f_m = 708.33 \) Hz. For normal helical gears, the spectrum shows peaks only at \( f_m \) and \( 2f_m \). With minor pitting, additional peaks appear at \( f_r \), \( 2f_r \), and \( 3f_r \) in the low-frequency range, and sidebands around \( f_m \) and \( 2f_m \) at intervals of \( f_r \). Severe pitting intensifies these effects, with higher amplitudes and more sidebands, indicating increased modulation due to the fault.

Multi-tooth pitting faults exhibit more pronounced dynamic responses. In the time domain, dual-tooth pitting results in a 9% increase in contact force amplitude and prominent periodic impacts at the rotational period. For triple-tooth pitting, the force distribution becomes sparse, with larger and more frequent impacts. Frequency-domain analysis shows that dual-tooth pitting enhances the low-frequency peaks at \( f_r \), \( 2f_r \), and \( 3f_r \) by 46%, and the meshing frequency peaks by 5%, with numerous sidebands. Triple-tooth pitting further amplifies low-frequency components by 150% and introduces abundant high-amplitude sidebands between 0 and 800 Hz. This progression highlights that as the number of pitted teeth increases, the sideband count and amplitude in the low-frequency region rise significantly, providing a clear indicator of fault severity in helical gears.

To validate the simulation results, we conduct experimental tests on a gear fatigue rig equipped with helical gears identical to those modeled. The setup includes a two-motor test bench where the input motor drives the gear pair at 2500 rpm, and the output motor applies a load torque of 840 N·m. Vibration acceleration signals are acquired using tri-axial accelerometers mounted on the input and output shaft end-covers, with a sampling frequency of 12,000 Hz. The experimental procedure involves running the gears through their lifecycle and capturing data at various stages of pitting progression.

In the experimental analysis for single-tooth pitting, the time-domain vibration signals show that normal gears exhibit steady acceleration without impacts. Minor pitting causes a slight increase in amplitude and sparse distribution, while severe pitting leads to a 10% amplitude rise and distinct periodic impacts at 0.024-second intervals. The frequency-domain envelope spectra confirm the simulation findings: normal gears display peaks only at \( f_m \) and \( 2f_m \); minor pitting introduces low-frequency peaks at \( f_r \), \( 2f_r \), and \( 3f_r \), along with sidebands around meshing frequencies; severe pitting increases the number and amplitude of these sidebands by 15%. This consistency between simulation and experiment underscores the reliability of our dynamic model for diagnosing single-tooth pitting in helical gears.

For multi-tooth pitting experiments, dual-tooth pitting results in denser impact patterns in the time domain, with periodic impacts at the rotational period, and a 12% increase in low-frequency peak amplitudes. Triple-tooth pitting causes a 100% surge in vibration amplitude, with sparse, large impacts. In the frequency domain, dual-tooth pitting raises the meshing frequency peaks by 11% and sideband counts, while triple-tooth pitting maximizes low-frequency amplitudes and sideband proliferation between 0 and 800 Hz. These experimental observations align closely with the simulation predictions, demonstrating that our model effectively captures the dynamic characteristics of multi-tooth pitting faults in helical gears.

The implications of this study are significant for the maintenance and design of helical gear systems. By correlating time-domain and frequency-domain features with pitting severity and extent, we establish diagnostic criteria: for single-tooth pitting, the appearance of rotational frequency harmonics and meshing frequency sidebands indicates fault presence, with impact periodicity and sideband amplitude correlating with severity. For multi-tooth pitting, the increase in low-frequency sideband count and amplitude serves as a key indicator. This approach enables early fault detection, reducing downtime and preventing failures in critical applications such as automotive transmissions and industrial machinery.

In conclusion, our dynamic simulation model for helical gears provides a robust framework for analyzing pitting faults. Through detailed modeling of contact mechanics and experimental validation, we demonstrate the ability to diagnose single-tooth and multi-tooth pitting based on vibration signal characteristics. The use of envelope spectrum analysis enhances fault detection sensitivity, making it a valuable tool for predictive maintenance. Future work could explore the effects of other fault types, such as cracks or wear, and integrate machine learning techniques for automated diagnosis. Overall, this research advances the understanding of helical gear dynamics and contributes to improved reliability in gear transmission systems.

Scroll to Top