Heavy-lift rotor unmanned aerial vehicles depend on a compact, high-torque transmission path from the engine or motor to each rotor shaft. The spiral bevel gear is a common final stage in that path because its curved, oblique teeth engage progressively rather than along the full face at once. This gives a high contact ratio, smooth angular transmission, and good bending capacity under heavy rotor loads. The same geometry also makes the gear sensitive to mesh impact, load transients, and resonance: the teeth enter and leave contact with a changing stiffness, and a sudden change in rotor load can produce a force pulse that is much larger than the steady torque.
For a heavy-lift rotor drive, two different dynamic questions therefore matter. The first is transient contact behavior: how the mesh force responds to speed ramps and load steps. The second is structural dynamics: whether the gear’s natural frequencies can be excited by shaft rotation, tooth mesh, or control-induced torque ripple. A contact simulation and a modal analysis address these questions at different levels, and together they define a preliminary operating envelope for the gear pair.
Meshing impact as a contact-stiffness problem
At the instant a tooth pair enters mesh, the surfaces deform elastically and the normal relative velocity across the contact is not zero. The resulting velocity mismatch produces a short impact pulse. In a rotor transmission, this pulse is superimposed on the steady tooth load and on the periodic stiffness variation caused by alternating single- and double-tooth contact. A second source of impact is external: rapid acceleration, gust response, or a sudden change in rotor load can apply a torque step to the driven gear.
A penalty-contact formulation based on Hertzian stiffness is often used to represent this behavior in a rigid-body dynamics solver. Let \(l\) be the instantaneous contact distance and \(l_1\) the separation at contact onset, so that \(\delta=l_1-l\) is the penetration. The normal contact force can be written as
$$
F=\begin{cases}
0, & \delta\le 0,\\[4pt]
\max\left[0,\;k\delta^{e}-c(\delta)\dot{\delta}\right], & \delta>0,
\end{cases}
$$
where \(k\) is the contact stiffness, \(e\) is the force exponent, and \(c(\delta)\) is a damping term that rises from zero at contact onset to a maximum value \(C_{\max}\) at a specified penetration \(d\). The stiffness for two elastic bodies in contact is obtained from
$$
k=\frac{4}{3}\,
\frac{\left(\dfrac{R_1R_2}{R_1+R_2}\right)^{1/2}}
{\dfrac{1-\nu_1^2}{E_1}+\dfrac{1-\nu_2^2}{E_2}},
$$
where \(R_1\) and \(R_2\) are the equivalent radii of the pinion and wheel at the mesh, \(E_1\) and \(E_2\) are their elastic moduli, and \(\nu_1\) and \(\nu_2\) are their Poisson ratios. For a carburized steel pair, the model used \(E_1=E_2=2.0675\times10^5\ \text{MPa}\) and \(\nu_1=\nu_2=0.3\), giving a contact stiffness of \(6.73\times10^5\ \text{N/mm}\). The remaining contact inputs are summarized in Table 1.
Table 1. Contact-model inputs for the spiral bevel gear pair.
| Parameter | Value |
|---|---|
| Elastic modulus, \(E_1=E_2\) | \(2.0675\times10^5\ \text{MPa}\) |
| Poisson ratio, \(\nu_1=\nu_2\) | 0.3 |
| Contact stiffness, \(k\) | \(6.73\times10^5\ \text{N/mm}\) |
| Force exponent, \(e\) | 1.5 |
| Maximum damping, \(C_{\max}\) | 45 N·s/mm |
| Full-damping penetration, \(d\) | 0.1 mm |
| Static / kinetic friction coefficient | 0.08 / 0.06 |
| Pinion speed stages | 1000, 2000, 3000 r/min |
| Constant resisting torque on driven gear | 100 N·mm |
| Simulation time / steps | 0.3 s / 5000 |
The contact simulation used the Adams solver within UG. The pinion speed was prescribed as a sequence of ramps: from rest to \(6000^\circ/\text{s}\) over 0.05 s, held to 0.10 s; then to \(12{,}000^\circ/\text{s}\) by 0.15 s, held to 0.20 s; then to \(18{,}000^\circ/\text{s}\) by 0.25 s. These plateaus correspond to 1000, 2000, and 3000 r/min. The driven gear carried a constant resisting torque of 100 N·mm in the contact model.
Speed ramps: impact grows with speed
The pinion followed the prescribed ramps and settled at the three speed plateaus. The driven gear reached mean angular velocities of 5999.97, 12,000.06, and 17,999.97\(^\circ/\text{s}\), respectively. The deviations from the commanded values are negligible, which confirms that the contact and constraint definitions preserve the intended transmission ratio.
The instantaneous output speed was not perfectly constant. As the mesh alternates between single- and double-tooth contact, the mesh stiffness changes periodically, and each transition introduces a local impact. The driven gear therefore oscillated around its mean speed. With increasing speed, both the amplitude and the frequency of that oscillation grew, and the contact-force trace showed a corresponding increase in the impact band. A start-up spike appeared when the gear first engaged under load, after which the force settled into a periodic pattern; at each higher speed stage, the force fluctuation became more pronounced.
This trend is important for a rotor drive because a high shaft speed does not simply scale the steady torque. It also raises the rate at which tooth pairs enter and leave contact, so the impact component becomes a larger fraction of the tooth load. Sustained operation near the highest modeled speed is therefore the least favorable condition for fatigue initiation at the tooth root and for surface distress. The practical implication is to limit continuous speed and to avoid unnecessarily fast acceleration ramps, especially when the rotor is already carrying a high thrust load.
Load steps: heavier load can smooth the mesh, but unloading is critical
The second load case held the speed constant and stepped the resisting torque through 30, 70, and 100 N·mm. At the lower load, the effective contact ratio was below its fully loaded value, so the mesh force fluctuated over a wider band after each engagement. As the load increased, more tooth surface came into effective contact, the mesh stiffness variation was reduced, and the force trace became smoother. In other words, a heavier steady load can make the mesh appear more stable, because the teeth are pressed into a broader contact pattern and the single/double-contact transition is less abrupt.
That result should not be read as a recommendation to operate at the highest possible load. In a rotor UAV, the sequence is more important than the steady-state point. Takeoff is dominated by a rapid increase in rotor speed and thrust, so the gear sees a speed-induced impact before the load has fully developed. Once the vehicle is airborne and fully loaded, the mesh can be relatively stable. The critical event is often the opposite transition: when the vehicle releases a payload or changes mission mode, the rotor load drops suddenly, and the gear experiences a rapid unloading pulse. The contact force can spike as the tooth pair recovers from deflection, and a severe pulse may initiate root cracking or tooth fracture.
For design and operation, the payload-release case should therefore be treated as a distinct load case rather than as a minor variation of cruise. The initial takeoff mass, the payload distribution, and the release sequence all affect the magnitude of that unloading transient. Adjusting flight attitude or scheduling the release so that the rotor retains a controlled minimum load can reduce the shock. The gear pair should also be checked for the combination of high speed and sudden load reduction, not only for steady high torque.
Modal analysis and resonance margins
Contact simulation describes forced response, but it does not by itself identify the frequencies at which the gear structure can amplify vibration. A linear modal analysis provides that information. The equation of motion for a damped mechanical system is
$$
\mathbf{M}\ddot{\mathbf{x}}+\mathbf{C}\dot{\mathbf{x}}+\mathbf{K}\mathbf{x}=\mathbf{F}(t),
$$
where \(\mathbf{M}\), \(\mathbf{C}\), and \(\mathbf{K}\) are the mass, damping, and stiffness matrices, and \(\mathbf{F}(t)\) is the external load vector. For free vibration with damping neglected, the equation reduces to
$$
\mathbf{M}\ddot{\mathbf{x}}+\mathbf{K}\mathbf{x}=0.
$$
Assuming harmonic motion, \(\mathbf{x}=\boldsymbol{\phi}\sin(\omega t)\), gives the eigenvalue problem
$$
\left(\mathbf{K}-\omega^2\mathbf{M}\right)\boldsymbol{\phi}=0.
$$
The natural frequencies and mode shapes are obtained by solving this eigenvalue problem. Because the lowest modes usually produce the largest relative deformation and are the most likely to be excited, the first six modes of the pinion are sufficient for a preliminary resonance check. The finite-element model used a tetrahedral mesh so that mass and stiffness were distributed through the tooth and rim geometry. The resulting frequencies are listed in Table 2.
Table 2. First six natural frequencies of the pinion and their dominant motion.
| Mode | Natural frequency (Hz) | Dominant motion |
|---|---|---|
| 1 | \(3.5803\times10^{-3}\) | Near-zero circumferential response |
| 2 | 16,286.0 | Rocking about the gear Y-axis |
| 3 | 16,445.0 | Rocking about the gear Y-axis |
| 4 | 20,928.0 | Bending |
| 5 | 20,977.0 | Torsion |
| 6 | 23,233.0 | Bending |
The first computed frequency is essentially zero. In an unconstrained modal solution this represents a near-rigid circumferential response rather than an elastic resonance, so it is not treated as an operating limit. The first elastic mode is the second mode at 16,286 Hz. Its deformation is concentrated at the large end near the tooth root, with the two sides of the gear rocking about the Y-axis. The third mode is only 159 Hz higher and has a similar rocking character. The fourth and sixth modes are bending modes with maximum deformation at the outer rim, while the fifth mode is torsional and reaches its maximum at the large-end root under axial twist.
For the rotor drive, the lowest elastic frequency of 16.3 kHz is far above the shaft rotation frequency: at 3000 r/min, the shaft turns at 50 Hz. The more relevant excitation is the tooth-mesh frequency and its harmonics, which depend on the number of teeth and the actual rotor speed. The closeness of the second and third modes means that a resonance check should not focus only on the lowest value; both modes should be compared with the mesh excitation spectrum. If a harmonic or a control-induced torque ripple falls near 16.3–16.4 kHz, the resulting dynamic load can be amplified even when the steady contact stress appears acceptable.
Operating envelope and design checks
The two analyses lead to a consistent set of design checks for a heavy-lift rotor transmission using spiral bevel gears:
- Limit sustained speed and acceleration. Higher speed increases both the amplitude and frequency of mesh-force fluctuation. Fast ramps add a start-up impact before the load has stabilized.
- Treat payload release as a design load case. A high steady load can smooth the mesh, but a sudden load drop can produce a severe unloading impact. The release sequence and the retained rotor load should be controlled.
- Keep mesh excitation away from the first elastic modes. The lowest elastic mode at 16,286 Hz and the nearby mode at 16,445 Hz should be checked against the tooth-mesh harmonics of the actual gear set.
- Use contact and modal models for different purposes. The penalty-contact simulation is a screening tool for impact trends and speed/load sensitivity; the modal model identifies frequency margins and deformation patterns. Neither replaces strain measurement, fatigue testing, or teardown inspection of the assembled rotor drive.
For a heavy-lift rotor UAV, the practical outcome is a preliminary operating envelope rather than a single safe speed. The gear should be operated with controlled acceleration, a defined minimum rotor load during payload release, and a validated separation between mesh excitation and the 16.3–16.4 kHz modal pair. Within that envelope, the spiral bevel gear can retain its main advantage: high torque capacity in a compact, right-angle rotor drive.