Influence of Root Crack on Time-Varying Meshing Stiffness of Wind Power Planetary Helical Gears

As a researcher focused on wind turbine drivetrain systems, I have extensively studied the dynamic behavior of planetary gear systems, particularly the effects of tooth root cracks on the time-varying meshing stiffness of helical gears. Wind turbines operate in harsh environments, such as coastal areas, grasslands, and deserts, where the planetary gear sets in the gearbox are subjected to low-speed and heavy-load conditions. The tooth root, being a stress concentration zone, is highly susceptible to crack initiation, which can lead to catastrophic failures like tooth breakage. In this article, I analyze how different types and extents of root cracks impact the meshing stiffness of planet gears in wind turbine systems, using nonlinear dynamics theory and numerical methods. The study considers both penetrating and non-penetrating cracks, modeled with parabolic fits, and evaluates their influence through stiffness degradation rates. The findings provide critical insights for fault diagnosis and durability enhancement in wind energy applications.

Planetary gear systems are integral to wind turbine gearboxes due to their high transmission ratios, stability, and load-bearing capacity. In a typical setup, the ring gear is fixed, the planet carrier rotates driven by the main shaft, and the planet gears revolve around the sun gear, enabling speed increase. The sun gear, being the smallest and fastest-moving component, experiences the most frequent meshing cycles, making it prone to root cracks. My analysis centers on a wind turbine planetary helical gear system with parameters as follows: the sun gear has 25 teeth, the planet gears have 32 teeth, a normal module of 20 mm, a pressure angle of 20 degrees, a helix angle of 7.5 degrees, and a face width of 120 mm. The material used is 17CrNiMo6, with a density of 7850 kg/m³, Young’s modulus of 2.07×10¹¹ Pa, and Poisson’s ratio of 0.3. These parameters form the basis for modeling and simulation.

To model root cracks in the planet gears, I classified them into penetrating and non-penetrating types based on their propagation paths. For non-penetrating cracks, I used parabolic fits in both the depth and width directions. The crack initiates at point A on the root fillet, with a tangent at 30 degrees to the tooth centerline. The depth direction is modeled as a parabola ABC, where A is the start point, C is the imaginary end point symmetric about the tooth axis, and B is the lowest point, with AB forming a 75-degree angle to the centerline. The actual crack depth is defined by point R on this parabola. In the width direction, another parabola RG is used, with R as the vertex and G as the endpoint. The crack extent is quantified as a percentage of the maximum possible depth or width. For depth, it is given by: $$\text{Depth Percentage} = \frac{\text{Actual Crack Depth}}{\text{Maximum Crack Depth}} \times 100\%$$ I considered four depth percentages: 10%, 20%, 50%, and 80%. For width, I defined it as the actual crack width relative to the total face width (120 mm), with percentages of 25%, 50%, 75%, and 100%. Penetrating cracks, on the other hand, are modeled with a parabolic fit in the depth direction and extend linearly across the entire face width. A total of twenty different crack configurations were analyzed, including combinations of depth and width extents, with a constant crack thickness of 0.02 mm.

The time-varying meshing stiffness is calculated using finite element contact statics analysis. I developed a three-dimensional model of the gear pair, with the planet gear’s inner surface fully constrained and the sun gear’s inner surface subjected to a tangential torque of \( T = 1.65 \times 10^5 \, \text{N·mm} \), while other degrees of freedom are restricted. The model employs hexahedral meshing, with refined grids in the contact zone and crack regions to ensure accuracy. Under load, the gears deform, and the driving gear rotates by a small angle \( \Delta \theta \), which is derived from the average deformation \( \Delta \delta \) at the inner surface and the radius \( r_1 \): $$\Delta \theta = \frac{\Delta \delta}{r_1}$$ The torsional mesh stiffness \( K_t \) is then calculated as: $$K_t = \frac{T}{\Delta \theta}$$ The mesh stiffness \( K \) relates to the torsional stiffness through the base circle radius \( r_b \): $$K = \frac{K_t}{r_b^2} = \frac{T}{r_b^2 \cdot \Delta \theta}$$ This formulation allows me to compute the stiffness over a meshing cycle, considering the helical gear’s engagement characteristics, where the contact line length varies gradually, leading to linear transitions in stiffness.

My results show that for a healthy gear pair without cracks, the time-varying meshing stiffness exhibits cyclic variations, with higher stiffness during double-tooth contact and lower stiffness during single-tooth contact. The overlap ratio for this helical gear system is 1.86, resulting in alternating single and double-tooth engagement. The stiffness transitions linearly between these states due to the progressive nature of helical gear contact. When root cracks are introduced, the stiffness decreases, with the most significant reduction occurring during single-tooth engagement. For penetrating cracks, as the depth increases, the stiffness degradation becomes more pronounced. For instance, at 10% depth, the average stiffness decreases by 4.31% compared to the uncracked case; at 20%, it drops by 5.72%; at 50%, by 11.25%; and at 80%, by 18.83%. This indicates that deeper cracks in planet gears lead to a more severe loss of stiffness, compromising system integrity.

Non-penetrating cracks also reduce stiffness, but to a lesser extent than penetrating ones. For example, at a fixed width of 100%, a 10% depth crack causes a 10.69% reduction in average stiffness, while deeper cracks show diminishing returns in degradation rate. Similarly, for a fixed depth of 80%, increasing the width from 25% to 100% results in stiffness reductions of 1.79%, 3.89%, 8.09%, and 17.69%, respectively. This highlights that width extension has a more dramatic impact on stiffness than depth extension for non-penetrating cracks. To quantify these effects, I used the stiffness degradation rate \( H \), defined as: $$H = \frac{K_P – K_C}{K_P} \times 100\%$$ where \( K_P \) is the average stiffness for the uncracked case and \( K_C \) is for the cracked case. The following tables summarize the degradation rates for various crack scenarios, emphasizing the critical role of planet gears in these dynamics.

Table 1: Stiffness Degradation Rates for Penetrating Cracks in Planet Gears
Depth Percentage Average Stiffness (10^8 N/mm) Degradation Rate (%)
0% (Uncracked) 3.314 0
10% 3.216 2.93
20% 3.093 6.67
50% 2.861 13.64
80% 2.511 24.22
Table 2: Stiffness Degradation Rates for Non-Penetrating Cracks in Planet Gears (100% Width)
Depth Percentage Average Stiffness (10^8 N/mm) Degradation Rate (%)
0% (Uncracked) 3.314 0
10% 3.056 7.78
20% 2.941 11.25
50% 2.808 15.27
80% 2.727 17.69
Table 3: Stiffness Degradation Rates for Non-Penetrating Cracks in Planet Gears (80% Depth)
Width Percentage Average Stiffness (10^8 N/mm) Degradation Rate (%)
0% (Uncracked) 3.314 0
25% 3.254 1.79
50% 3.185 3.89
75% 3.046 8.09
100% 2.727 17.69

In discussing these results, I emphasize that the behavior of planet gears under crack conditions is crucial for understanding overall system dynamics. The helical nature of the gears causes stiffness to vary cyclically, and cracks exacerbate this by introducing localized weaknesses. Penetrating cracks, due to their full-width extension, cause more significant stiffness reductions than non-penetrating ones, as they affect a larger portion of the tooth. For non-penetrating cracks, the initial depth increase leads to a high degradation rate, but further depth extension has a reduced effect. Conversely, width extension in non-penetrating cracks starts with a low degradation rate that accelerates as the crack propagates. This nonlinear behavior underscores the importance of monitoring both depth and width in fault detection for planet gears. The stiffness degradation rate serves as a reliable indicator for assessing crack severity, with higher rates signaling more critical damage that could lead to failure in wind turbine operations.

In conclusion, my analysis demonstrates that root cracks in planetary helical gears significantly alter the time-varying meshing stiffness, with penetrating cracks causing more severe degradation than non-penetrating ones. The stiffness reduction is most pronounced during single-tooth engagement, and for non-penetrating cracks, width extension has a greater impact than depth extension. The degradation rates provide a quantitative measure for evaluating crack effects, highlighting the vulnerability of planet gears in wind turbine systems. These insights can inform maintenance strategies and design improvements to enhance the reliability and lifespan of wind energy infrastructure. Future work could explore dynamic responses and crack propagation under varying operational conditions to further optimize planetary gear performance.

Scroll to Top