In the field of power transmission, the spur and pinion gear pair serves as a fundamental component, finding extensive application in demanding sectors such as new energy vehicles and aerospace machinery. The accurate prediction of mesh stiffness, a primary source of internal excitation, is paramount for enhancing transmission accuracy, optimizing gear design, and mitigating vibration and noise. Traditional methods for calculating mesh stiffness, including the potential energy method, finite element analysis (FEM), and hybrid approaches, have provided significant insights. However, these methods predominantly operate within a quasi-static framework, often neglecting the profound influence of operational kinematics. Specifically, the driving speed and its concomitant centrifugal effect are frequently overlooked, despite their critical role in the dynamic meshing process of high-speed spur and pinion systems.
This omission leads to an incomplete understanding of gear dynamics under realistic operating conditions. At elevated rotational speeds, significant centrifugal forces develop, altering the deformation response and effective stiffness of the gear teeth. This phenomenon, known as centrifugal stiffening or hardening, can substantially impact the system’s natural frequencies and dynamic mesh stiffness (DMS), thereby influencing its vibration signature and noise emission. Consequently, developing a dynamic model that intrinsically incorporates the effects of driving speed and centrifugal force is essential for a more authentic analysis of spur and pinion gear performance.
To address this challenge, we propose an original computational algorithm grounded in the Euler-Bernoulli beam theory and Hamilton’s principle. This framework formulates the equations of motion for a rotating flexible spur gear, explicitly integrating the kinetic energy contributions from translation and rotation, the strain energy, and the work done by centrifugal forces. The resulting model allows for the direct calculation of dynamic tooth deflection under the combined action of mesh forces and speed-dependent centrifugal loading. Subsequently, the time-varying dynamic mesh stiffness is derived. This paper details the derivation of the governing equations, presents the novel solution algorithm, validates the approach against established methods, and extensively investigates the influence of driving speed, centrifugal effects, and material properties on the dynamic mesh stiffness of spur and pinion gears.
1. Formulation of Motion Equations for a Rotating Flexible Spur Gear
The core of the proposed method lies in modeling an individual spur gear tooth as a rotating cantilever beam attached to a rigid hub. The centrifugal field arising from the gear’s rotation significantly affects its deformation state. As illustrated in the modeling framework, the gear exhibits two distinct states: an undeformed reference state under static load (dashed line) and an expanded, deformed state under centrifugal force \(F_v\) (solid line). With increasing driving speed \(\dot{\theta}\), the effective centrifugal pull causes the gear material to displace outward, influencing the boundary conditions and stiffness for tooth bending.
Consider a point on the flexible tooth. Its position vector \(\mathbf{P_T}\) in the inertial frame, accounting for both rigid-body rotation and elastic deformation \((u, v)\) in the axial and transverse directions (defined in a rotating coordinate system attached to the gear), is given by:
$$\mathbf{P_T} = \begin{bmatrix} x – u_r + u \cos\theta – v \sin\theta \\ z_{i,x} + v_r + u \sin\theta + v \cos\theta \end{bmatrix}$$
where \(u_r, v_r\) are displacements due to pure rotation, and \(\theta\) is the angular position.
The corresponding velocity vector \(\dot{\mathbf{P_T}}\) is derived by differentiation. The kinetic energy \(T_i\) of the gear element is then:
$$T_i = \frac{1}{2} \rho \int \dot{\mathbf{P_T}} \cdot \dot{\mathbf{P_T}} \, dV$$
where \(\rho\) is material density and \(dV\) is the volume element.
The nonlinear axial strain \(\varepsilon_{xx}\), considering moderate rotations, is:
$$\varepsilon_{xx} = u’ + \frac{1}{2}\left[(u’)^2 + (v’)^2\right]$$
where the prime (\(‘\)) denotes differentiation with respect to the axial coordinate \(x\). The potential (strain) energy \(U_i\) is:
$$U_i = \frac{1}{2} \int E \varepsilon_{xx}^2 \, dV$$
where \(E\) is the Young’s modulus.
Applying Hamilton’s principle,
$$\delta \int_{t_1}^{t_2} (U – T) \, dt = 0$$
yields the discretized equations of motion for the finite element model of the spur gear tooth:
$$(\mathbf{M_t} + \mathbf{M_r})\ddot{\mathbf{X}} + (\mathbf{C_r} + \mathbf{C_p})\dot{\mathbf{X}} + (\mathbf{K_e} + \mathbf{K_v})\mathbf{X} = \mathbf{F} + \mathbf{F_v} \tag{1}$$
where:
- \(\mathbf{M_t}, \mathbf{M_r}\): Translational and rotational mass matrices.
- \(\mathbf{C_r}, \mathbf{C_p}\): Rayleigh damping and gyroscopic damping matrices.
- \(\mathbf{K_e}, \mathbf{K_v}\): Elastic stiffness and centrifugal stiffness matrices.
- \(\mathbf{F}, \mathbf{F_v}\): Mesh force vector and centrifugal force vector.
- \(\mathbf{X}\): Nodal displacement vector.
The centrifugal stiffness matrix \(\mathbf{K_v}\) is symmetric and proportional to the square of the driving speed \(\dot{\theta}^2\), representing the centrifugal hardening effect. The centrifugal force vector \(\mathbf{F_v}\) for an element of length \(\Delta L\) is:
$$\mathbf{F_v} = -\frac{\rho A \Delta L \dot{\theta}^2}{2} \left( L + \frac{\Delta L}{3}, 0, 0, L + \frac{2\Delta L}{3}, 0, 0 \right)^T \tag{2}$$
This formulation captures the dynamic coupling between the spur gear’s flexibility, its rotational inertia, and the centrifugal force field, providing a foundation for accurate dynamic mesh stiffness calculation.
2. Computational Algorithm for Dynamic Mesh Stiffness Under Centrifugal Effects
The dynamic mesh stiffness (DMS) is not a static property but evolves during the meshing cycle due to changing contact conditions, dynamic excitation, and speed effects. Our algorithm calculates this by solving the dynamic response of the spur and pinion teeth as they engage under load and rotation. The process simulates the load application moving from the start to the end of the single-tooth contact path.
The external force vector \(\mathbf{F}_i\) is applied only at the instantaneous contact point \(i\) along the tooth profile:
$$\mathbf{F}_i = [0, 0, 0, \dots, F_i \sin(\beta_i), F_i \cos(\beta_i), M_i, \dots, 0, 0, 0]^T \tag{3}$$
where \(F_i\) is the mesh force, \(\beta_i\) is the pressure angle at point \(i\), and \(M_i\) is the associated moment. The pressure angle is geometry-dependent:
$$\beta_i = \arccos\left(\frac{R_{bp}}{\sqrt{x_i^2 + z_{i,x}^2}}\right) – \arctan\left(\frac{z_{i,x}}{x_i}\right)$$
where \(R_{bp}\) is the base circle radius of the spur or pinion gear.
The engagement velocity \(v_i\) at the contact point, crucial for determining dynamic interaction timing, is:
$$v_i = \dot{\theta} \sqrt{x_i^2 + z_{i,x}^2} \tag{4}$$
To solve the dynamic equation (1) for the displacement \(\mathbf{X}_i\) under moving load and centrifugal effects, we employ the Newmark-\(\beta\) time integration scheme. The time step \(\Delta t_i\) for the load moving between two successive contact points \(i\) and \(i+1\) is determined by their average engagement velocity:
$$\Delta t_i = \frac{2 \sqrt{\Delta x_i^2 + \Delta y_i^2}}{v_i + v_{i+1}} \tag{5}$$
This adaptive step ensures the solution captures the transient response as the load traverses the tooth.
The dynamic single-tooth stiffness \(k_{pi}\) for the pinion (and similarly \(k_{gi}\) for the gear) at contact point \(i\) is calculated from the dynamic deflection \((\Delta x_{i,x}, \Delta x_{i,y})\) extracted from \(\mathbf{X}_i\):
$$k_{pi} = \frac{F_i}{\Delta x_{i,x} \cos(\pi/2 – \beta_i) + \Delta x_{i,y} \cos \beta_i} \tag{6}$$
For a single tooth pair in contact, the effective dynamic mesh stiffness \(k_{ms}\) is the series combination:
$$k_{ms} = \frac{k_{pi} \cdot k_{gi}}{k_{pi} + k_{gi}} \tag{7}$$
In the double-tooth contact region, the total mesh stiffness is the sum of the stiffnesses of the two engaged spur and pinion tooth pairs. The algorithm iterates this process over the entire path of contact, incorporating the updated centrifugal force \(\mathbf{F_v}(\dot{\theta})\) at each step, thereby yielding the time-varying DMS profile influenced by centrifugal effects.
3. Model Validation and Analysis of Centrifugal Effects
To validate the proposed algorithm, we compare its results with those from a commercial finite element analysis (FEA) software (ANSYS APDL) under both quasi-static and dynamic conditions. The parameters for the example spur and pinion gear pair are listed in Table 1.
| Parameter | Pinion | Gear (Spur) |
|---|---|---|
| Number of Teeth | 27 | 41 |
| Module (mm) | 2.5 | |
| Pressure Angle (°) | 20 | |
| Face Width (mm) | 10 | |
| Young’s Modulus, \(E\) (GPa) | 207 | |
| Poisson’s Ratio | 0.3 | |
| Density, \(\rho\) (kg/m³) | 7850 | |
First, a quasi-static comparison was made by setting an extremely low driving speed (\(\dot{\theta}_p = 0.01\) rpm) to minimize centrifugal effects. The calculated Single Tooth Dynamic Stiffness (STDS) from our algorithm showed excellent agreement with the static stiffness obtained from ANSYS, validating the underlying static stiffness model.
Second, a dynamic comparison was conducted at a significant speed (\(\dot{\theta}_p = 300\) rpm). As shown in Figure 4 (result chart), our algorithm predicts a dynamic stiffness that oscillates around the static stiffness curve. The amplitude of this oscillation is larger in our results compared to the ANSYS dynamic analysis. This discrepancy is attributed to two factors: 1) The proposed algorithm fully accounts for the centrifugal stiffening effect via the \(\mathbf{K_v}\) matrix, which increases the effective restoring force, and 2) inherent differences in dynamic formulation and integration between the custom algorithm and the commercial FEM solver. The observed oscillations stem from dynamic excitation; the prior tooth deflection does not fully recover before the next engagement point is loaded, leading to a cumulative dynamic effect that our model captures.
We further investigated the impact of increasing driving speed. Figure 6 illustrates the STDS for both the spur gear and pinion at various speeds (300, 3000, 5000, 9000 rpm). Key observations are:
- The mean value of STDS increases with speed due to centrifugal hardening.
- The amplitude of fluctuation around the mean also increases with speed.
- The fluctuation is most pronounced near the initial point of contact (root of the driving pinion, tip of the driven spur gear), where kinetic energy and dynamic excitation are highest, and attenuates towards the end of the contact path.
The overall dynamic mesh stiffness over a full engagement cycle is presented in Figure 7. It conclusively demonstrates that as the driving speed rises:
- The mean dynamic mesh stiffness increases nonlinearly.
- The fluctuation amplitude within the meshing cycle grows.
- The effect is more pronounced in the double-tooth contact region because two pairs of teeth contribute to and are affected by the centrifugal field, amplifying the dynamic interaction compared to the single-tooth contact region.
This confirms a strong nonlinear relationship between centrifugal effect and dynamic mesh stiffness for the spur and pinion system.
4. Influence of Gear Material and Dynamic Frequency Analysis
The material properties of the spur and pinion gears significantly modulate how centrifugal effects influence dynamic behavior. We analyze four common gear materials: Hard Aluminum Alloy, Cast Iron, Carbon Fiber Nylon, and Ceramics. Their properties are summarized in Table 2.
| Material | Elastic Modulus, \(E\) (GPa) | Density, \(\rho\) (kg/m³) | Specific Modulus, \(E/\rho\) (m) |
|---|---|---|---|
| Hard Aluminum Alloy | 70 | 2700 | ~25,900 |
| Cast Iron | 207 | 7890 | ~26,200 |
| Carbon Fiber Nylon | 230 | 1760 | ~130,700 |
| Ceramics | 410 | 3150 | ~130,200 |
The growth rate and volatility rate (normalized fluctuation amplitude) of DMS with increasing speed for these materials are plotted in Figure 10. Two critical insights emerge:
- Elastic Modulus Dominance: Materials with a lower elastic modulus (like Aluminum) exhibit a much larger percentage increase in DMS with speed. The centrifugal stress \(\sigma_c \propto \rho \omega^2 R^2\) causes a strain; the resulting stiffening \(\Delta k \propto E \cdot \epsilon_c\) has a weaker modulus-dependence than the static stiffness \(k_{static} \propto E\), leading to a higher growth rate \((\Delta k / k_{static})\) for low-\(E\) materials.
- Density Effect: For materials with similar elastic moduli (e.g., Cast Iron and Carbon Fiber Nylon), the one with higher density (Cast Iron) shows a greater DMS growth rate and volatility. This is because the centrifugal force and the associated \(\mathbf{K_v}\) matrix scale with \(\rho \dot{\theta}^2\).
Thus, for high-speed spur and pinion applications, materials with high specific modulus (like Carbon Fiber Nylon) are advantageous as they better resist centrifugal distortion and exhibit more stable DMS.
Natural Frequency and Campbell Diagram Analysis: The centrifugal effect also alters the system’s vibrational characteristics. We analyzed the Campbell diagrams (natural frequency vs. driving speed) for the spur gear and pinion under two inertia configurations: considering only the translational mass matrix (\(\mathbf{M_t}\)), and considering both translational and rotational mass matrices (\(\mathbf{M_t}+\mathbf{M_r}\)).
Key findings from the Campbell diagrams (Figures 8 & 9):
- With only \(\mathbf{M_t}\), certain modes (e.g., the first and third for the pinion) remain constant with speed, while others (like the second and fifth for the gear) show a strong increase due to centrifugal stiffening.
- When \(\mathbf{M_t}+\mathbf{M_r}\) is included, the dynamics become richer. A distinct frequency “veering” phenomenon is observed near 12,790 rpm, indicating strong coupling between different rotational modes induced by gyroscopic and centrifugal effects.
- The inclusion of rotational inertia is critical for accurate high-speed dynamic prediction, as its effect becomes increasingly significant with speed.
The accuracy of our model in predicting natural frequencies was validated against FEA at low speed, with relative errors below 5%, as shown in a comparison table. The impact of the mass matrix formulation on DMS growth rate calculation is quantified in Table 4, confirming that the full inertia formulation (\(\mathbf{M_t}+\mathbf{M_r}\)) yields more accurate results, especially at speeds exceeding 6000 rpm.
| Mass Matrix | Material | Growth Rate @ 2000 rpm | Growth Rate @ 6000 rpm | Growth Rate @ 10000 rpm |
|---|---|---|---|---|
| \(\mathbf{M_t}\) only | Hard Aluminum Alloy | 14.31% | 50.24% | 70.10% |
| \(\mathbf{M_t}+\mathbf{M_r}\) | Hard Aluminum Alloy | 15.71% | 53.85% | 86.50% |
| \(\mathbf{M_t}\) only | Cast Iron | 12.78% | 46.77% | 65.57% |
| \(\mathbf{M_t}+\mathbf{M_r}\) | Cast Iron | 13.45% | 49.93% | 69.02% |
| \(\mathbf{M_t}\) only | Ceramics | 3.39% | 21.73% | 38.26% |
| \(\mathbf{M_t}+\mathbf{M_r}\) | Ceramics | 3.83% | 24.60% | 47.10% |
| \(\mathbf{M_t}\) only | Carbon Fiber Nylon | 3.30% | 20.34% | 36.56% |
| \(\mathbf{M_t}+\mathbf{M_r}\) | Carbon Fiber Nylon | 3.30% | 21.74% | 37.71% |
5. Conclusion
This study has developed and validated a novel computational framework for determining the dynamic mesh stiffness of spur and pinion gears, with explicit consideration of driving speed and the resulting centrifugal effects. The governing equations, derived from Hamilton’s principle for a rotating Euler-Bernoulli beam, integrate translational and rotational inertia, elastic deformation, and centrifugal stiffening in a consistent dynamic formulation.
The proposed algorithm successfully calculates the time-varying dynamic mesh stiffness by solving for the tooth’s dynamic response under a moving mesh load within a centrifugal force field. Validation and numerical analyses lead to the following principal conclusions:
- The dynamic mesh stiffness is not static but oscillates around the quasi-static stiffness value. Both the mean value and the amplitude of oscillation of the DMS increase nonlinearly with driving speed due to the centrifugal hardening effect and intensified dynamic excitation.
- The centrifugal effect has a more pronounced impact on the mesh stiffness in the double-tooth contact region compared to the single-tooth contact region, owing to the involvement of two tooth pairs.
- The proposed model provides a more authentic prediction of DMS under high-speed conditions than traditional quasi-static methods. It reveals complex dynamic phenomena such as frequency veering in Campbell diagrams, which are critical for avoiding resonant conditions in high-speed spur and pinion drives.
- Material selection profoundly influences the sensitivity of DMS to speed. A low elastic modulus leads to a high percentage increase in DMS with speed, while a high density amplifies both the growth rate and volatility of DMS. Materials with a high specific modulus (e.g., Carbon Fiber Nylon) offer superior stability for high-speed applications.
- For accurate high-speed dynamic analysis of spur and pinion gears, it is essential to include both the translational and rotational inertia matrices in the model, as their coupling with centrifugal effects becomes significant at elevated speeds.
This work provides a valuable reference for the dynamic analysis, design optimization, and vibration/noise reduction of spur and pinion gear systems operating under significant centrifugal fields, particularly in advanced automotive and aerospace transmissions.