Spiral bevel gears are widely used in helicopter reducers, vehicle drivelines, and construction machinery because their high overlap, smooth motion, and high load capacity allow them to transmit torque between intersecting shafts under severe duty. Their dynamic behavior is not determined by tooth geometry alone. Mesh stiffness, support flexibility, transmitted load, and operating speed interact to produce vibration, dynamic mesh force, and impact-like effects at tooth engagement and disengagement.
Most mesh models represent time-varying stiffness as a function of rotation angle. That captures the nominal meshing cycle, but it does not account for the fact that a loaded tooth pair enters contact earlier and leaves contact later than the theoretical contact segment. The resulting extended contact raises the effective contact ratio, and the increase depends on load. The analysis here treats this effect explicitly for a spiral bevel gear pair and compares steady operation with a rapid speed ramp.
Contact Beyond the Theoretical Mesh Segment
In the normal section at mid-face width, the theoretical contact path of a spiral bevel gear pair runs between points B and A. At zero load, contact begins at B and ends at A. Under load, elastic compression of the teeth allows a pair to enter before B and leave after A, so the actual contact path extends from B′ to A′. The separation between an entering or leaving tooth pair and the theoretical contact point is defined as the gear teeth separation distance (GTSD), denoted ξa and ξb at the two ends.
Let S be the instantaneous contact point and |BS| its position along the line of action. The position is periodic in pinion rotation:
\[
|BS|=\mathrm{mod}(r_{b1}\theta_1,|B’A’|)
\]
Within the theoretical segment BA, the relative mesh deflection is
\[
\lambda_0=(x_p-x_g)a_1+(y_p-y_g)a_2+(z_p-z_g)a_3+r_{b1}\theta_p-r_{b2}\theta_g
\]
where x, y, and z are translations, θ is angular displacement, a1, a2, and a3 are direction cosines of the line of action, and rb1 and rb2 are the base radii at the mesh point. Outside the theoretical segment, only the part of the normal approach that closes the GTSD contributes:
\[
\lambda=
\begin{cases}
\max(\delta_0-\xi_a,0), & \text{before }B\text{ (early engagement)},\\[4pt]
\max(\delta_0-\xi_b,0), & \text{after }A\text{ (delayed release)}.
\end{cases}
\]
Here δ0 is the reference normal approach along the line of action. This piecewise description makes the effective contact ratio load-dependent: as the transmitted load rises, the tooth pair penetrates farther into the extended zones, and the actual contact interval grows. The mesh stiffness therefore depends not only on rotation angle but also on the load state.
Eight-Degree-of-Freedom Mesh Model
A lumped-parameter model represents the pinion and gear as rigid bodies connected by a mesh spring-damper and supported by linear stiffness and damping. Each member has three translations and one rotation, giving eight degrees of freedom. The two shaft axes intersect at the origin, with the pinion axis along X and the gear axis along Z. Tooth friction and torsional swing are neglected, so the model isolates the influence of mesh stiffness, extended contact, and speed variation.
The equations of motion are assembled in matrix form:
\[
\mathbf M\ddot{\mathbf X}+\mathbf C\dot{\mathbf X}+\mathbf K\mathbf X=\mathbf F,
\qquad
\mathbf X=[x_p,y_p,z_p,\theta_p,x_g,y_g,z_g,\theta_g]^{T}.
\]
The normal mesh force is
\[
F_n=k_n(t)\lambda+c_n\dot\lambda,
\]
where kn(t) is the time-varying mesh stiffness and cn is the mesh damping. The normal force is projected onto the coordinate directions through the line-of-action direction cosines. For the torsional coordinates,
\[
I_p\ddot\theta_p=T_p-F_n r_{b1},
\qquad
I_g\ddot\theta_g=-T_g+F_n r_{b2},
\]
where Ip and Ig are moments of inertia, Tp is input torque, and Tg is load torque. The eight second-order equations are reduced to sixteen first-order equations and integrated with MATLAB’s ode45 solver.
Case Definition and Steady-State Response
The case is a first-stage reduction pair representative of an electrically driven helicopter transmission. The main gear parameters are listed in Table 1.
| Parameter | Pinion | Gear |
|---|---|---|
| Number of teeth | 23 | 75 |
| Module (mm) | 3.25 | 3.25 |
| Normal pressure angle | 20° | 20° |
| Midpoint spiral angle | 35° | 35° |
| Pitch cone angle | 17.05° | 72.95° |
At a pinion speed of 5000 r/min and a load of 300 N·m, the mesh frequency is
\[
f_m=\frac{5000\times 23}{60}\approx 1916\ \mathrm{Hz}.
\]
With initial displacement terms set to zero, the response is periodic in displacement, velocity, and normal mesh force. Impact-like peaks appear at the start and end of each mesh cycle, where tooth pairs enter or leave the extended contact zone and the mesh stiffness changes abruptly. Among the translational coordinates, the axial z displacement is largest, consistent with the axial force component produced by the spiral bevel geometry. The pinion and gear displacement and velocity traces follow the same variation trends.
The dynamic mesh force remains periodic and no tooth separation occurs, so the transmission stays in continuous contact. In the frequency domain, the main amplitudes occur at the mesh frequency of 1916 Hz and its multiples. The largest amplitude appears at five times the mesh frequency, approximately 9.58 kHz.
Mesh-Force Spectrum and Dynamic Load Factor
The dynamic load factor compares the maximum dynamic mesh force with the maximum static mesh force:
\[
K_v=\frac{\max(DF_m(t))}{\max(SF_m(t))}.
\]
A discrete speed sweep from 0 to 5000 r/min shows that Kv rises rapidly at low speed and then fluctuates, forming several pronounced peaks and valleys. These peaks identify speed bands of elevated dynamic amplification rather than a single critical speed. For design purposes, the existence of multiple high-response bands means that a gear pair cannot be assessed at one nominal operating point alone.
| Operating case | Input condition | Observed response |
|---|---|---|
| Steady state | Pinion 5000 r/min; load 300 N·m | Periodic displacement, velocity, and mesh force; entry and exit shocks; no tooth separation; spectrum dominated by 1916 Hz and harmonics, with maximum at five times mesh frequency |
| Speed ramp | Pinion 0 to 5000 r/min in 0.5 s; load 300 N·m | Mean angular speed increases linearly; larger dynamic-force fluctuations than steady state; no tooth separation |
Variable-Speed Run-Up
For the transient case, the pinion speed increases uniformly from 0 to 5000 r/min in 0.5 s at a load of 300 N·m. The mean angular velocity increases linearly over the ramp, with an average angular acceleration of approximately \(1.05\times10^3\ \mathrm{rad/s^2}\). This confirms that the model can represent a spiral bevel gear transmission during speed change, not only at a fixed operating point.
Compared with steady operation, the variable-speed transient produces significantly larger amplitude fluctuations in the dynamic mesh force. The tooth pair still does not separate, so contact is maintained throughout the ramp. The larger excursions are consistent with sweeping the mesh excitation through speed-dependent amplification bands while the load-dependent contact ratio changes continuously. A variable-speed drive therefore cannot be evaluated only at its maximum speed or at a few steady operating points; the run-up path itself affects the dynamic force history.
Design Implications and Model Scope
Two design consequences follow from the analysis. First, the effective contact ratio of a spiral bevel gear pair is load-dependent because tooth deflection extends contact beyond the theoretical segment. A nominal contact ratio calculated from unloaded geometry does not fully describe the mesh excitation under load. Second, rapid speed changes produce larger dynamic mesh-force fluctuations than steady operation at the same final speed. Speed schedules, control ramps, and transient torque management therefore matter for vibration and load control.
The model is suited to comparative studies of mesh stiffness, support stiffness, load level, and run-up rate. It does not include tooth-surface friction, torsional swing, thermal effects, or detailed contact pressure. It therefore provides system-level vibration and force trends rather than local stress or fatigue-life prediction. The numerical results apply to the stated geometry, load, and speed range; other designs require reanalysis with their own parameters.
Conclusions
Load-induced extended contact changes the mesh excitation of spiral bevel gear pairs. In steady operation at 5000 r/min and 300 N·m, the response is periodic, with impact-like effects at mesh entry and exit and a dynamic-force spectrum dominated by the 1916 Hz mesh frequency and its harmonics; the largest amplitude occurs at five times the mesh frequency. A 0.5 s ramp from 0 to 5000 r/min produces larger dynamic-force fluctuations but no loss of tooth contact. The results support load- and speed-dependent dynamic assessment of spiral bevel gear transmissions, especially for variable-speed drivelines in which transient operation can be more demanding than a fixed design point.