As a graduate student specializing in mechanical engineering, my research focused on the dynamic behavior of heavy-load and high-speed involute spur gears used in aviation applications. The aviation industry is a strategic high-tech sector that directly reflects national defense capabilities and economic strength. Gear transmissions are fundamental components of aero-engines, and most operate under severe conditions involving heavy loads and high linear velocities. The relatively high sliding speeds on tooth surfaces and the substantial loads induce significant vibration, noise, and additional dynamic forces, all of which threaten the safety and stability of the entire transmission system. Long-term engineering practice has revealed that gear failures are a major contributor to engine malfunctions, significantly restricting the advancement of our aviation industry. This paper presents my comprehensive investigation into the dynamic factor used in the bending strength verification of aviation involute spur gears, aiming to clarify the mechanism of internal additional dynamic loads, compare classical calculation methods, and propose a modified formula tailored to actual dynamic responses.
1. Introduction and Research Significance
High-performance gear systems are essential for modern aviation applications, and the trend toward high reliability, efficiency, and lightweight design has intensified the challenges in gear transmission development. Aviation gear transmissions exhibit unique characteristics including elevated sliding velocities, high instantaneous tooth surface temperatures, increased probabilities of scuffing and wear, and significant vibration and noise. These factors directly impact the dynamic load, a critical parameter in gear design. My thesis addresses a specific problem: the classical formulas used for dynamic factor calculation often yield unreasonable results when applied to supercritical speed operations of heavy-load high-speed involute spur gears. This problem motivated me to conduct a detailed study on the dynamic mechanism, analyze the inherent characteristics of gear structures, and develop a new calculation method for the dynamic factor.
2. Dynamic Load Mechanism of Aviation Spur Gears
To understand the dynamic behavior of spur gears, I began with the fundamental mechanism of dynamic load generation. The dynamic load in gear meshing arises from two primary sources: manufacturing errors and time-varying mesh stiffness. These factors create internal excitation that leads to vibration and additional dynamic forces on the tooth surfaces. As the rotational speed increases, these vibrations become more pronounced, particularly when the gear’s precision is relatively low.
I used a torsional vibration model to analyze the spur gear pair, where the angular displacements of the driving and driven gears are denoted as θ₁ and θ₂. The normal dynamic load Wd can be expressed using the tooth stiffness Ki and combined error ei at the contact point:
$$W_d = \sum_i K_i (x_1 – x_2 – e_i) \tag{1}$$
Here, x₁ = r₁θ₁ and x₂ = r₂θ₂ represent the displacements along the line of action. The total mesh stiffness is,
$$K(t) = \sum_i K_i \tag{2}$$
and the equivalent mass M is given by,
$$M = \frac{I_1 I_2}{I_1 r_2^2 + I_2 r_1^2} \tag{3}$$
Considering damping, the relative torsional vibration equation becomes,
$$M\ddot{x} + 2\zeta\sqrt{MK(t)}\,\dot{x} + K(t)x = W + \sum_i K_i e_i \tag{4}$$
From this equation, the natural frequency fe of the system is obtained as,
$$f_e = \frac{1}{2\pi}\sqrt{\frac{\bar{K}(t)}{M}} \tag{5}$$
where \bar{K}(t) is the average mesh stiffness over one mesh cycle. The fundamental excitation frequency is the mesh frequency fm, defined as,
$$f_m = \frac{n z}{60} \tag{6}$$
The ratio fm/fe is the critical speed ratio, a key parameter that determines the operating regime. When this ratio equals 1, the system experiences resonance, leading to severe vibration and potential failure. The critical speed of the gear pair, based on the impact theory, is found when the effective tooth profile error equals the displacement under load during the single-tooth engagement period. This yields,
$$n_1 = \frac{30}{z_1}\sqrt{\frac{W}{2 m e}} \tag{7}$$
For a more rigorous approach, I employed the vibration theory, which considers the gear as a spring-mass system. The critical speed in terms of the mesh stiffness cγ and induced mass mred is,
$$n_E = \frac{30\times10^3}{z_1 \pi}\sqrt{\frac{c_\gamma}{m_{red}}} \tag{8}$$
This value defines the threshold between the subcritical (N<1), critical (N=1), and supercritical (N>1) operating regions.
To evaluate the dynamic load, I compared several classical standards. The ISO 6336 standard provides a widely used method that divides the operating range into different regimes. The key formulas are summarized in the following table.
| Operating range | Calculation formula in ISO 6336 |
|---|---|
| Subcritical (N ≤ 0.85) | $$K_v = 1 + N(C_{v1} B_p + C_{v2} B_f + C_{v3} B_k)$$ |
| Main resonance (0.85 < N ≤ 1.15) | $$K_v = C_{v1} B_p + C_{v2} B_f + C_{v4} B_k + 1$$ |
| Transition (1.15 < N < 1.5) | $$K_v = K_{v(N=1.5)} + \frac{K_{v(N=1.5)} – K_{v(N=1.15)}}{0.35}(N – 1.15)$$ |
| Supercritical (N ≥ 1.5) | $$K_v = C_{v5} B_p + C_{v6} B_f + C_{v7}$$ |
In my initial evaluation using data from an actual aviation gear pair, I discovered a significant discrepancy. At supercritical speeds, the ISO B method yielded a dynamic factor exceeding 3, implying a very high dynamic load that would logically lead to tooth breakage. However, the gear was functioning normally in practice. This observation confirmed that a specialized approach is necessary for high-speed heavy-load spur gears operating in the supercritical region.
3. Dynamic Analysis and Simulation of Spur Gear Pairs
The dynamic factor is directly related to the dynamic meshing force. I established a nonlinear torsional vibration model for an aviation spur gear pair to investigate the dynamic meshing force under supercritical speed. The model incorporates time-varying mesh stiffness, transmission error, and the backlash nonlinearity, which is essential for accurately capturing the dynamic behavior.
The governing differential equation using the dynamic transmission error as the degree of freedom is,
$$m_e \ddot{x}(t) + c_m \dot{x}(t) + k(t) g[x(t)] = F_m – m_e \ddot{e}(t) \tag{9}$$
where me is the equivalent mass, cm is the meshing damping, k(t) is the time-varying mesh stiffness, and e(t) is the static transmission error. The backlash function g[x(t)] is defined as,
$$g[x(t)] = \begin{cases} x(t) – b, & x(t) > b \\ 0, & |x(t)| \le b \\ x(t) + b, & x(t) < -b \end{cases} \tag{10}$$
where b is half the backlash value.
For the mesh stiffness calculation, I used the Weber-Banascheck energy method, which considers the tooth as a combination of bending, shear, and compressive deformations, plus the flexibility of the gear body and local contact deformation. The time-varying mesh stiffness was computed numerically and is shown to vary periodically due to the alternating single and double tooth engagement. To incorporate this into the dynamic model, I decomposed the stiffness into a mean and a fluctuating part and expressed it as a Fourier series:
$$k(t) = k_m + \sum_{n=1}^{N} [a_n \cos(n\omega_m t) + b_n \sin(n\omega_m t)] \tag{11}$$
The static transmission error, generated by manufacturing inaccuracies, was also modeled as a periodic function using a Fourier series representation.
I solved the dynamic equation using the fourth-order Runge-Kutta method in MATLAB. The variation of the dynamic factor Kv with rotational speed, derived from the numerical simulation of the dynamic meshing force, is presented in the figure below. This curve shows that the dynamic factor exhibits a fluctuating trend, initially increasing with speed, then decreasing and stabilizing at a lower level in the supercritical region. This is in contrast to the prediction of the ISO B method, which continues to increase.
To further validate my numerical results, I performed a dynamic simulation using ADAMS software. The impact function in ADAMS, based on the Hertzian contact theory, was used to simulate the gear contact. The contact stiffness parameter was calculated from the material properties and geometry:
$$K = \frac{4}{3} \sqrt{R}\,E^* \tag{12}$$
where R is the equivalent radius of curvature and E* is the equivalent elastic modulus. The simulation results for the contact force, after the system reached steady-state, yielded a dynamic factor of approximately 1.22 at the rated speed, which is consistent with the numerical calculation. These findings confirmed that the actual dynamic load for high-speed heavy-load spur gears at supercritical speeds is much lower than predicted by the ISO B method.
| Parameter | Value |
|---|---|
| Nonlinear exponent (Force exponent) | 1.5 |
| Damping coefficient | 50 N·s/mm |
| Static friction coefficient | 0.3 |
| Dynamic friction coefficient | 0.1 |
| Penetration depth | 0.1 mm |
4. Natural Vibration Characteristics of Aviation Spur Gears
The natural frequency of a gear structure is a critical parameter for avoiding resonance. I performed a comprehensive modal analysis of the aviation spur gear pair using the finite element method (FEM). The analysis was conducted in two stages: free modal analysis (unconstrained) and prestressed contact modal analysis, which accounts for the contact constraints and centrifugal stresses arising from high-speed rotation.
The undamped free vibration equation of a gear structure is:
$$[M]\{\ddot{X}\} + [K]\{X\} = \{0\} \tag{13}$$
The eigenvalue problem is then formulated as,
$$\left([K] – \omega_i^2 [M]\right)\{X\}_i = \{0\} \tag{14}$$
where ωi is the i-th natural frequency and {X}i is the corresponding mode shape.
I first established precise involute spur gear models in ANSYS using APDL command streams, then meshed the gear pair and established the contact pair. For the prestressed analysis, I performed a static analysis by applying the actual working speed of 54,600 r/min as an angular velocity to the driving gear, enabling the calculation of centrifugal stresses. The results showed that the prestressed natural frequencies are slightly higher than the non-prestressed ones, indicating the stiffening effect of centrifugal forces. The first 20 natural frequencies for both cases are summarized below.
| Order | Natural Frequency (Non-prestressed, Hz) | Natural Frequency (Prestressed, Hz) | Mode Shape Description |
|---|---|---|---|
| 1 | 9002.7 | 9020.6 | Rigid body rotation about Z axis |
| 2 | 9009.5 | 9162.6 | Large gear rocking around X axis |
| 3 | 9166.1 | 9828.2 | Large gear first-order nodal diameter |
| 4 | 9453.6 | 10959 | Large gear second-order nodal diameter (saddle) |
| 5 | 10296 | 15804 | Large gear bending causing pinion rocking |
| 6 | 10744 | 21392 | Pinion rocking around X axis |
| 7 | 21224 | 22156 | Pinion rocking around Y axis |
| 8 | 21396 | 23773 | Pinion first-order nodal diameter |
| 9 | 23063 | 25044 | Pinion and large gear nodal diameter |
| 10 | 24876 | 25899 | Pinion and large gear nodal diameter |
| 11 | 24946 | 26574 | Pinion and large gear nodal diameter |
| 12 | 25111 | 28922 | Large gear radial vibration |
| 13 | 25880 | 30264 | Large gear radial vibration |
| 14 | 28657 | 32409 | Large gear radial vibration |
| 15 | 28728 | 34630 | Combined nodal diameter vibration |
| 16 | 31664 | 39186 | Pinion meshing bending vibration |
| 17 | 40843 | 41137 | Combined nodal diameter vibration |
| 18 | 40892 | 41896 | Combined nodal diameter vibration |
| 19 | 41114 | 42956 | Combined nodal diameter vibration |
| 20 | 41811 | 50222 | Higher-order nodal diameter |
I also investigated the effect of key design parameters on the natural frequencies of the gear. This analysis is crucial because it provides a pathway to optimize designs to avoid resonance. I systematically varied the pressure angle, module, tooth number, addendum coefficient, and face width, and observed the following trends:
- Pressure angle: Natural frequencies increase with increases in the pressure angle.
- Module: There is a strong negative correlation between module value and natural frequencies.
- Tooth number: Natural frequencies decrease as the number of teeth increases.
- Addendum coefficient: Natural frequencies slightly decrease with increases in the addendum coefficient, leveling off near 1.
- Face width: Natural frequencies initially increase with face width, then stabilize in constrained conditions.
The mesh frequency of the gear pair at rated speed was calculated as fm = 18,200 Hz. A comparison with the natural frequencies listed above confirmed that the mesh frequency is safely separated from the system’s natural frequencies, thus avoiding resonance. This confirms the smooth dynamic behavior of the gear during operation.
5. Development of a New Dynamic Factor Formula
My investigation highlighted the limitations of existing dynamic factor calculation methods for heavy-load high-speed spur gears. The ISO B method yields excessively high and incorrect values in the supercritical region, whereas the simplified C method is not suitable for high-speed conditions. To address this, I developed a new, statistically-based formula specifically for the dynamic factor of aviation spur gears.
My new approach considers three primary influencing factors: processing accuracy, rotational speed, and mesh stiffness. I defined the dynamic factor as a function of these parameters:
$$K_v = f(C, n, k_m) = 1 + B_n + B_C + B_k \tag{15}$$
where Bn, BC, and Bk are dimensionless coefficients representing the influence of rotational speed, processing accuracy, and mesh stiffness, respectively. To determine the functional form of these coefficients, I performed extensive calculations using various methods, including the ISO B method, the ISO C method, the Russian standard, and the dynamic simulations. I generated a large dataset of dynamic factors for various gears and operating conditions, and then used the 1stOpt software with nonlinear fitting algorithms to develop the final formula.
The resulting new formula is expressed as:
$$K_v = 1 + \frac{79219\,E_1}{n^2} + \frac{3}{C} + \frac{100}{S_1} \tag{16}$$
For consistency with existing standards, I also proposed the formula within a segmented framework, using the critical speed ratio N. For the subcritical region (N ≤ 0.85), the formula is:
$$K_v = 1 + \frac{79219\,E_1}{n_1^2} + \frac{3}{C} + \frac{100}{S_1} \tag{17}$$
where E₁ is a combined amplification factor for shape and pitch errors, and S₁ is the stiffness correction factor. For the supercritical region (N ≥ 1.5):
$$K_v = 1 + \frac{79219\,E_2}{n_1^2} + \frac{3}{100} + \frac{100}{S_2} \tag{18}$$
In this formulation, E₂ and S₂ utilize specific amplification factors from the ISO B method for the supercritical region. The new formula is applicable for gears with a processing precision grade of 4–8, a tooth number range of 17–100 (up to 500 for general gears), a gear module range of 1.5–10 mm, and consistency with the actual operational state.
Verification and Comparative Analysis
To validate the new formula, I used three sets of aviation spur gear pairs from an actual aero-engine transmission system: two auxiliary drive gear pairs (No. 1 and No. 2) and one center drive gear pair (No. 3). The gear pair parameters are listed in the table below.
| Parameter | Gear Pair No.1 (Pinion/Gear) | Gear Pair No.2 (Pinion/Gear) | Gear Pair No.3 (Pinion/Gear) |
|---|---|---|---|
| Number of teeth | 13 / 16 | 26 / 32 | 20 / 26 |
| Module (mm) | 1.5 | 1.5 | 1.5 |
| Pressure angle (°) | 25 | 25 | 25 |
| Addendum coefficient | 1 | 1 | 1 |
| Face width (mm) | 4 / 3 | 6 / 3 | 4 / 3 |
| Speed (r/min) | 14,788 / 12,015 | 42,000 / 34,125 | 54,600 / 42,000 |
| Mesh stiffness (N/(mm·μm)) | 8.24 | 10.21 | — |
| Processing precision | 5 | 5 | 5 |
Gear pair No. 1 operates in the subcritical region, while gear pairs No. 2 and No. 3 operate in the supercritical region. The dynamic factors calculated using different methods—the new formula, the ISO B and C methods, and the Russian (ГОСТ) formula—are compared below.
| Gear Pair | Kv (ISO B/C) | Kv (New Formula) | Relative Error vs. ISO B/C (%) | Kv (Russian) | Relative Error vs. Russian (%) |
|---|---|---|---|---|---|
| No.1 (subcritical) | 1.2109 (B) | 1.1611 | 4.11 | — | — |
| No.2 (supercritical) | 1.4530 (C) | 1.3730 | 5.50 | — | — |
| No.3 (supercritical) | 2.9915 (B) | 1.4974 | 49.94 | 1.4558 | 2.86 |
For gear pair No. 1, operating in the subcritical region, the new formula shows good agreement with the ISO B method, with a relative error of only 4.11%. For gear pair No. 2, in the supercritical region, the dynamic factor from the new formula is 5.5% lower than the ISO C method result. For gear pair No. 3, the new formula yields a significantly lower and more realistic value than the ISO B method (49.94% difference), while showing excellent agreement with the Russian formula, with an average error of just 2.26% across the speed range. This confirms that the new formula effectively resolves the overestimation issue present in the ISO method for supercritical speed operations of heavy-load high-speed spur gears.
I further compared the new formula against existing methods under varying processing precision, mesh stiffness, and rotational speed. A comprehensive comparison of the dynamic factor values for the gear pairs across a range of processing precision grades is presented in the following table.
| Processing Precision C | Kv (ISO C) | Kv (New Formula) | Relative Error (%) |
|---|---|---|---|
| 4 | 1.2477 | 1.3291 | 6.52 |
| 5 | 1.4530 | 1.4774 | 1.68 |
| 6 | 1.6777 | 1.7907 | 6.73 |
| 7 | 2.2408 | 2.3130 | 3.22 |
| 8 | 2.6083 | 3.0441 | 16.71 |
The dynamic factor increases with precision grade, and the maximum error is 16.71% at grade 8. However, since aviation gears are typically grade 5, the new formula performs exceptionally well, confirming its practical applicability.
The comparison across varying mesh stiffness for gear pair No. 1 is summarized in the table below. Both methods show the same trend of a slight increase in the dynamic factor with increasing stiffness, with the maximum error below 4%.
| Mesh Stiffness k (N/(mm·μm)) | Kv (ISO B) | Kv (New Formula) | Relative Error (%) |
|---|---|---|---|
| 13.89 | 1.1952 | 1.1735 | 1.82 |
| 15.16 | 1.2070 | 1.1794 | 2.29 |
| 16.42 | 1.2181 | 1.1853 | 2.69 |
| 17.68 | 1.2287 | 1.1913 | 3.04 |
| 18.95 | 1.2390 | 1.1972 | 3.37 |
| 20.21 | 1.2488 | 1.2032 | 3.65 |
When the rotational speed is varied for gear pair No. 1 in the subcritical region, both the new formula and the standard methods (B and C) show that the dynamic factor increases with speed. The average relative error between the new formula and the ISO B method is less than 10%, and similarly for the ISO C method.
For the high-speed gear pair No. 3, the comparison across different speeds is presented in the following table. The ISO B method produces values nearly double those of the new formula. In contrast, the Russian standard and the new formula yield very similar results, with the maximum error being only 5.84%.
| Speed (r/min) | Kv (New) | Kv (ISO B) | Kv (Russian) | Error vs. ISO B (%) | Error vs. Russian (%) |
|---|---|---|---|---|---|
| 30000 | 1.2424 | 2.5481 | 1.2735 | 51.24 | 2.44 |
| 34000 | 1.2832 | 2.6119 | 1.3100 | 50.87 | 2.04 |
| 38000 | 1.3291 | 2.6738 | 1.3464 | 50.29 | 1.28 |
| 42000 | 1.3801 | 2.7342 | 1.3829 | 49.52 | 0.20 |
| 46000 | 1.4362 | 2.8308 | 1.4194 | 49.26 | 1.18 |
| 50000 | 1.4974 | 2.9915 | 1.4558 | 49.94 | 2.86 |
| 54000 | 1.5636 | 3.2089 | 1.4923 | 51.27 | 4.78 |
| 56000 | 1.5987 | 3.2948 | 1.5105 | 51.48 | 5.84 |
6. Conclusions
This research addressed the problem of inaccurate dynamic factor calculations for heavy-load, high-speed involute spur gears in aviation applications. The key conclusions from my work are summarized as follows:
- Dynamic Load Mechanism: The dynamic load in high-speed spur gears is primarily caused by combined excitation from time-varying mesh stiffness and gear errors. The severity of the dynamic load strongly depends on the operating speed relative to the natural frequency of the gear system. In the supercritical region, the dynamic load tends to stabilize or even decrease with increasing speed, which contradicts the predictions of some classical methods.
- Dynamic Simulation: Both numerical simulations and dynamic simulations in ADAMS demonstrated that the dynamic factor for an aviation spur gear at supercritical speeds is around 1.2 to 1.5, significantly lower than values derived from the ISO B method. This confirms the practical feasibility of high-speed gear operation and the necessity for a dedicated dynamic factor calculation method.
- Natural Frequency Analysis: Finite element modal analysis of the aviation gear pair revealed that the natural frequencies are influenced by centrifugal stresses due to high-speed rotation. The mesh frequency is below the fundamental natural frequencies of the system, thus avoiding resonance. Design parameter studies show that gear module has the most significant influence on natural frequency, offering a path for optimizing gear dynamics.
- New Dynamic Factor Formula: A new statistical formula has been developed based on a comprehensive set of data from various calculation methods and dynamic simulations. This formula accounts for the effects of processing accuracy, rotational speed, and mesh stiffness. Extensive validation using real aviation gear pairs shows that it provides accurate and reliable results in both the subcritical and supercritical regions, with errors generally less than 5% compared to validated methods.
This work provides a robust theoretical foundation and a practical calculation tool for the design of high-performance aviation spur gears. Future work should focus on experimental validation using specialized test rigs to measure actual dynamic tooth loads, enabling further refinement of the proposed formula. This will lead to more precise dynamic load predictions and more efficient, reliable gear designs for aerospace applications.