Spiral bevel gears carry torque between intersecting shafts with a high contact ratio, large load capacity and comparatively low noise, which is why they dominate the reduction stages of helicopter and automotive drivetrains and appear widely in construction machinery. Because the gear pair sets the boundary conditions for the rest of the drivetrain, its dynamic behaviour — vibration displacement, transmission error and dynamic mesh force — largely determines system stability. The analysis below examines how contact that extends beyond the theoretical line of action, caused by tooth deflection under load, alters that behaviour at constant speed and during a speed ramp.
Why the theoretical line of action is not sufficient
Lumped-parameter formulations of spiral bevel gears commonly couple time-varying mesh stiffness with backlash, friction, support stiffness and body displacement, and reproduce measured vibration trends with reasonable fidelity. In most of these models, mesh stiffness is written as a function of the instantaneous rotation angle alone, and meshing is assumed to begin and end at the nominal entry and exit points of the theoretical contact path. That assumption is inconsistent with a loaded gear pair. Mounting error and elastic deformation of the teeth let a tooth pair engage before the theoretical entry point and remain in contact after the theoretical exit point, so the actual contact zone is longer than the nominal one. The phenomenon is referred to here as extended contact, or out-of-mesh engagement.
Two consequences follow. First, the real contact ratio becomes load dependent: the heavier the transmitted torque, the greater the extension, and the more tooth pairs share the load at any instant. Second, the mesh stiffness discontinuities that a conventional model places at the nominal entry and exit instants are displaced in time and reduced in amplitude. Since the dynamic mesh force is driven by the rate of change of stiffness, a model that ignores this effect misrepresents both the magnitude and the phase of the excitation, and therefore the predicted force amplitude and its spectral content.
Geometry of extended contact and the separation distance
Consider a single tooth pair viewed in the normal plane at mid-face width, with A and B the theoretical exit and entry points. Under no load the teeth do not deform, the contact path keeps its nominal length and the pair transmits the theoretical ratio. Once torque is applied, the driving flank deflects, tooth pair I enters contact ahead of the theoretical entry point and tooth pair III leaves contact after the theoretical exit point, lengthening the effective contact path. The gap that separates a tooth pair about to enter, or just having left, from its mating flank in the undeformed state is defined as the gear teeth separation distance (GTSD), measured along the mesh line and denoted ξa at entry and ξb at exit.
Position along the contact path is tracked through the intersection point S of the pinion profile with the path B’A’, whose distance from the datum is a periodic function of pinion rotation:
$$|BS|=\mathrm{mod}\left(r_{b1}\theta_1,\ |B’A’|\right)$$
where mod is the remainder function and rb1 the equivalent base radius of the pinion. Inside the theoretical contact zone the deformation along the mesh line follows the relative displacement of the two bodies:
$$\lambda_0=(x_p-x_g)a_1+(y_p-y_g)a_2+(z_p-z_g)a_3+\theta_p r_{b1}-\theta_g r_{b2}$$
with a1, a2, a3 the direction cosines of the common normal in the fixed frame. Outside that zone, deformation is limited by the separation distance, so the tooth pair can only carry load once the approach exceeds the local gap:
$$\lambda=\begin{cases}\max\left(\delta_0-\xi_a,\ 0\right), & 0<|BS|<|B’B|\\[4pt] \max\left(\delta_0-\xi_b,\ 0\right), & |B’A|<|BS|<|A’A|\end{cases}$$
Here δ0 is the loaded approach of the reference pair. The formulation is continuous: as the load rises, δ0 grows, more of the extension interval satisfies δ0 > ξ, and the number of simultaneously loaded pairs increases accordingly.
Eight-degree-of-freedom dynamic model
The gear pair is represented by a lumped-parameter model with eight degrees of freedom — three translational and one rotational per body — taking the intersection of the two shaft axes as origin, the pinion axis as X and the gear axis as Z. Tooth friction and torsional pendulum motion are excluded. The mesh normal is projected onto the eight coordinates to give the vector V, whose rotational entries are −rb1 and +rb2, with α the normal pressure angle, δp the pinion pitch cone angle and β the mid-face spiral angle. The normal mesh force and its components are
$$F_n=k_n(t)\lambda+c_n\dot{\lambda}$$
$$F_x=-F_n\left(\sin\alpha\cos\delta_p+\cos\alpha\sin\beta\sin\delta_p\right)$$
$$F_y=F_n\left(\sin\alpha\sin\delta_p-\cos\alpha\sin\beta\cos\delta_p\right)$$
$$F_z=F_n\cos\alpha\cos\beta$$
where kn(t) is the time-varying mesh stiffness and cn the mesh damping. The equations of motion for the pinion (p) and gear (g) take the form
$$\mathbf{M}\ddot{\mathbf{X}}+\mathbf{C}\dot{\mathbf{X}}+\mathbf{K}\mathbf{X}=\mathbf{F}$$
with translational rows such as mpẍp + cxpẋp + kxpxp = −Fx and torsional rows Ipθ̈p = Tp − Fnrb1, Igθ̈g = −Tg + Fnrb2; M, C and K collect the mass and inertia, damping and support stiffness terms, and F the external load vector. Introducing the sixteen state variables (xp, yp, zp, θp, xg, yg, zg, θg and their first derivatives) reduces the system to sixteen first-order differential equations, integrated with the MATLAB ode45 solver from an initial state in which all displacement terms are zero.
Gear pair and operating conditions
The case examined is the first reduction stage of an electrically driven helicopter transmission. The geometric data are listed in Table 1; the two pitch cone angles sum to the 90° shaft angle.
| Parameter | Pinion | Gear |
|---|---|---|
| Number of teeth | 23 | 75 |
| Normal module (mm) | 3.25 | 3.25 |
| Normal pressure angle | 20° | 20° |
| Mid-face spiral angle | 35° | 35° |
| Pitch cone angle | 17.05° | 72.95° |
| Case | Pinion speed | Load torque on gear |
|---|---|---|
| Steady state | 5000 r/min constant | 300 N·m |
| Speed sweep | 0–5000 r/min, discrete steps | 300 N·m |
| Transient | 0 → 5000 r/min in 0.5 s, uniform acceleration | 300 N·m |
Steady-state response at 5000 r/min
At constant pinion speed the vibration displacement, vibration velocity and dynamic mesh force all vary periodically. The largest impact excursions occur at the beginning and at the end of each mesh cycle, where the number of loaded tooth pairs changes. Axial vibration in the Z direction dominates the translational response, and the pinion and gear traces follow closely matching trends in both displacement and velocity, as expected from a single mesh force acting on two bodies of comparable support stiffness. The transmission error remains bounded and the dynamic mesh force never falls to zero, so the teeth do not separate and the drive stays in continuous contact at this torque level.
In the frequency domain, the dynamic mesh force is concentrated at the mesh frequency and its harmonics. With 23 pinion teeth at 5000 r/min, the mesh frequency is
$$f_m=\frac{z_1 n_1}{60}=\frac{23\times 5000}{60}\approx 1916\ \text{Hz}$$
which matches the dominant line in the computed spectrum. The amplitude peak is not at the fundamental but at the fifth harmonic, near 9.6 kHz, a feature that follows from the shape of the stiffness variation across the extended contact interval rather than from the mean force level. For noise and fatigue assessment this means that a single-frequency check at 1916 Hz would underestimate the excitation: the harmonic content up to roughly the fifth order has to be included, and structural resonances of the housing in that band deserve attention.
Dynamic load factor over the speed range
To characterise speed sensitivity, the pinion speed was swept from 0 to 5000 r/min in discrete steps and the dynamic load factor was evaluated as
$$K_v=\frac{\max\left(DF_m(t)\right)}{\max\left(SF_m(t)\right)}$$
where DFm and SFm are the dynamic and static mesh forces. Kv rises steeply at low speed, then oscillates, forming a sequence of well-defined peaks and valleys rather than a monotonic trend. The peaks mark speeds at which the excitation harmonics coincide with a system natural frequency, and the valleys mark the speeds between them. Operating points close to a peak should be avoided in service, or damped by design, because the dynamic force can substantially exceed the static value even though the transmitted torque is unchanged.
Transient response during uniform acceleration
For the variable-speed case, the pinion is accelerated uniformly from rest to 5000 r/min in 0.5 s with the gear load held at 300 N·m, corresponding to a mean angular acceleration of about 1047 rad/s². The computed mean angular velocity increases linearly with time, which confirms that the model, including the extended-contact stiffness, remains numerically well behaved and physically consistent when the rotation angle — and hence the mesh stiffness — changes continuously with time.
The dynamic mesh force under the ramp shows noticeably larger amplitude fluctuations than the steady-state result at the same final speed, and the fluctuation persists throughout the acceleration rather than settling. No loss of contact occurs, so the amplification originates from the interaction between the accelerating mesh frequency and the modal response of the system: each instant of the ramp excites a different harmonic, and the transient has no time to decay before the mesh frequency moves on. The practical implication for variable-speed drives is that a design verified only at rated speed can still experience significantly higher dynamic tooth loads during run-up and coast-down, and that acceleration rate is a design variable in its own right.
Assessment
Treating the mesh as confined to the theoretical line of action is a simplification that costs accuracy in two specific ways. It removes the load dependence of the contact ratio, and it misplaces the stiffness transitions that generate the impact at the boundaries of each mesh cycle. Restoring extended contact into an eight-degree-of-freedom formulation of a spiral bevel gear pair reproduces the periodic shocks at the start and end of the mesh cycle, the dominance of axial vibration, the concentration of dynamic force at the mesh frequency and its harmonics with a maximum at the fifth order, and the peak–valley structure of the dynamic load factor between 0 and 5000 r/min. The same model remains stable through a 0.5 s ramp to rated speed and predicts that acceleration raises the fluctuation amplitude of the dynamic mesh force without causing tooth separation. For spiral bevel gears in helicopter and vehicle transmissions, both the harmonic content of the mesh force and the speed range traversed during transients therefore belong in the design envelope, not only the steady-state load at rated speed.