Why collapse the nutation train to a single stage
Nutation gearing extracts a large reduction from a wobbling motion rather than from a conventional epicyclic arrangement. A bevel gear whose axis is tilted by a small, constant nutation angle ε relative to the input shaft is rolled against a stationary mating gear; one full wobble of the tilted gear advances the mesh by only the tooth difference between the pair, so the ratio is set by that difference and not by the absolute tooth counts. Because the mechanism is a variant of small-tooth-difference planetary gearing, a compact diameter can carry a large ratio, which is attractive wherever radial space is the binding constraint — most notably in electric submersible screw-pump oil-lifting systems, where the reducer must fit inside the wellbore.
The bevel-gear nutation reducers built so far have mostly been two-stage systems. They reach very high ratios, but they pay for it with a larger part count, two simultaneous gear pairs whose meshing behaviour must be phased correctly, and dynamic characteristics that are harder to keep stable. A single-stage layout using double-circular-arc spiral bevel gears removes one of those meshing pairs while retaining the ratio and load capacity. The design target examined here was a radial envelope of no more than 108 mm, a transmitted power of 20 kW and a ratio of roughly 10.
Kinematics and mechanical layout
The classical illustration is a coin spinning on a table: before it settles, its axis sweeps a cone, and that slow precession is the nutation motion. Replace the coin and the table with an external bevel gear and a fixed internal bevel gear, and the kinematic pair is complete. The external gear rotates about its own axis, which is inclined to the input axis by the constant angle ε; the internal gear is axially and rotationally fixed. Describing the geometry through an imaginary crown gear that meshes with both members — its pitch angle being π/2 — the pitch cone angles satisfy
$$\delta_1+\delta_2+\varepsilon=\pi$$
where δ1 and δ2 are the pitch cone angles of the external and internal gears. The transmission ratio follows from the pitch-cone relation and from the tooth counts z1, z2:
$$i=\frac{\sin\delta_1}{\sin(\delta_1+\varepsilon)-\sin\delta_1}=\frac{z_1}{z_1-z_2}$$
Two consequences follow directly. First, as ε shrinks the ratio grows without bound, so a small nutation angle is the route to high reduction. Second, for a given ratio the two tooth counts are not independent; only one of them is a free design quantity. With the internal gear having more teeth than the external one, the output turns in the opposite sense to the input.
The hardware is correspondingly simple. An inclined shaft, machined integrally with the input shaft at the nutation angle ε, carries the external bevel gear, which is free to rotate on it. That gear is driven into mesh with the stationary internal gear. The spin of the nutating gear about the inclined axis is then passed to the output shaft through a universal joint, and the reduction appears between input and output. Since the ratio is generated by the tooth difference, keeping that difference small — below four teeth — is important: an excessive difference increases the unbalanced wobble excitation that the rotating tilted gear imposes on the supports.
Tooth geometry of the double-circular-arc spiral bevel gear
The teeth are generated by sweeping a double-circular-arc normal basic profile along the spiral tooth trace, so each tooth flank carries a convex and a concave arc segment. Design parameters are defined at the small end of the tooth, which therefore serves as the reference section: face width b, pitch cone distance R, normal module at the small end mn, helix angle β, mean pitch diameters dm1 and dm2, and small-end pitch diameters d1 and d2. The basic rack profile follows GB/T 12759-1991 for double-circular-arc cylindrical gears, and the load-capacity calculation follows GB/T 13799-1992, adapted to the bevel geometry through the cone angles.

Optimization formulation
The design problem was posed as volume minimisation. The overall size of the reducer is governed almost entirely by the two bevel gears, so the objective V(X) sums the material envelope of the external and internal gear bodies, each treated as a conical shell of face width b whose radial extent follows from its tooth count, module, helix angle and pitch cone angle.
The nutation angle ε and the two pitch cone angles are fixed in advance — they are linked through the kinematic relation and the required ratio — and do not enter the optimisation. With the ratio fixed, z1 and z2 are coupled, so only z1 is carried as a variable. The design vector is therefore
$$\mathbf{X}=[z_1,\;m_n,\;\beta,\;b]^{\mathsf T}=[X_1,\;X_2,\;X_3,\;X_4]^{\mathsf T}$$
and the problem takes the standard form
$$\min_{\mathbf X} V(\mathbf X)\quad\text{subject to}\quad g_j(\mathbf X)\le 0,\qquad j=1,\dots,12$$
Twelve constraints were imposed, grouped as follows.
| Group | Constraint | Limit |
|---|---|---|
| Root bending fatigue (g1) | σF computed from torque, application and dynamic factors, load-sharing factors K1 and KF2, and the integer part of the contact ratio | σF ≤ [σF] |
| Surface contact fatigue (g2) | σH from the same load terms with KH2 and the contact-strength material, ratio, helix and arc-length factors | σH ≤ [σH] |
| Tooth count (g3) | Lower bound on the external gear to preserve contact ratio and smoothness | z1 > 5 |
| Tooth difference (g4, g5) | Internal gear larger than external, difference kept small to limit wobble excitation | z2 > z1; z2 − z1 < 4 |
| Radial envelope (g6) | Large-end pitch diameter of the internal gear, allowing for housing wall thickness | z1mn + 2b sin δ2 ≤ 90 mm |
| Face width coefficient (g7, g8) | Ratio of face width to small-end pitch diameter | 1/5 ≤ b/d1 ≤ 1/3 |
| Helix angle (g9, g10) | Band chosen for transmission smoothness | 20° < β < 40° |
| Contact ratio (g11, g12) | Integer part of the overlap ratio, which enters the load-sharing terms | 1 ≤ ⌊b sin β/(π mn)⌋ ≤ 6 |
The integer part of the contact ratio used in the stress expressions is
$$\mu_\varepsilon=\left\lfloor\frac{b\sin\beta}{\pi m_n}\right\rfloor$$
with the raw overlap b sin β/(π mn) computed from the geometry. Because the model mixes integer tooth counts with continuous dimensions, it is a discrete optimisation problem; it was solved with LINGO.
Optimized gear parameters
| Parameter | Symbol | Value |
|---|---|---|
| External gear tooth count | z1 | 20 |
| Internal gear tooth count | z2 | 22 |
| Small-end normal module | mn | 2 mm |
| Helix angle | β | 35° |
| Face width | b | 12.2 mm |
The solution converges on a two-tooth difference, which satisfies the wobble-excitation limit, and on a helix angle in the middle of the permitted band. Derived checks confirm that the active constraints are respected: the small-end pitch diameter of the external gear is z1mn = 40 mm, giving a face width coefficient b/d1 ≈ 0.305 inside the 0.2–0.333 window, and the overlap b sin β/(π mn) ≈ 1.11, whose integer part of 1 lies within the required range. A three-dimensional model of both members was then built from these parameters using the double-circular-arc tooth surface equations.
| Operating quantity | Value |
|---|---|
| Rated power | 20 kW |
| Input speed | 3000 r/min |
| Input torque, derived from power and speed | ≈63.7 N·m |
| Ratio magnitude (z1 = 20, z2 = 22) | 10, output opposite to input |
| Output torque at ratio 10 | ≈637 N·m, matching the 636 N·m load applied in the analysis |
| Radial envelope target | ≤108 mm |
Finite element model and loading
The assembled gear pair was analysed with an explicit solver, treating the mesh as a transient contact problem rather than a static one, so that the rolling of the tilted gear through engagement could be followed in time.
| Item | Setting |
|---|---|
| Material | Structural steel |
| Elastic modulus | 210 GPa |
| Poisson ratio | 0.3 |
| Density | 7850 kg/m³ |
| Element type | 10-node quadratic tetrahedra, free meshing |
| Global seed size, teeth / gear blank | 1 mm / 4 mm |
| Contact definition | Hard contact, finite sliding |
| Boundary conditions | External gear hinged on the inclined shaft, five degrees of freedom removed and rotation about that axis free; all degrees of freedom of the internal gear fixed |
| Input speed | 3000 r/min |
| Output load | 636 N·m |
| Solution control | 1 s total, 100 increments |
| Evaluation window | Last 50 increments (0.5–1 s), after the mesh has settled |
Contact behaviour and stress distribution
Contact area per increment was extracted over the settled portion of the run. It peaks at 98 mm² and averages 42.8 mm². The magnitude is a direct consequence of the tooth form and the phasing: several teeth of the pair are in mesh simultaneously, and each engaged tooth carries two contact traces, one on the convex flank and one on the concave flank. The load is therefore distributed over a comparatively large area rather than concentrated on a single line of contact, which is the mechanism behind the smooth transmission observed for this gear type.
The Mises stress results show a consistent asymmetry. At every evaluated instant the maximum stress in the external gear exceeds that in the internal gear. The peak locations are concentrated at the tooth tip near the small end of the external gear, where the tooth is thinnest and the section changes abruptly; the geometric discontinuity, not the nominal contact pressure, drives the local value. Away from that corner the stress distribution over both members is broadly uniform, which indicates that the optimized proportions do not create a systematic overload zone along the tooth length.
Design implications and open questions
The single-stage configuration delivers the intended combination: a ratio of about 10 from a two-tooth difference, 20 kW transmitted through a gear pair whose largest pitch diameter stays near 90 mm, and a part count well below that of a two-stage train. For downhole applications, where the reducer diameter sets the maximum tool size, that trade is the decisive one.
The stress concentration at the small-end tip is the clearest target for the next iteration. Tip relief, edge chamfering, or a local correction of the double-circular-arc profile would move the peak away from the thinnest section without altering the ratio or the envelope. Attention to the small-end geometry is more productive here than increasing the module, which would grow the diameter against the radial constraint.
The scope of the reported analysis should be read narrowly. It covers one operating point — full load at 3000 r/min — with an elastic material model, hard contact and ideally rigid supports, over one second of simulated meshing. Lubricant film formation, friction and heat generation are outside the model, as are manufacturing errors and misalignment of the inclined shaft, all of which affect where contact actually lands on the flank. The wobble-induced excitation that the small tooth difference is meant to control has been bounded only through the geometric constraint, not evaluated dynamically. Efficiency measurement, thermal behaviour and bench testing of a prototype remain open, and the stress result should be treated as a screening indication for profile modification rather than as a validated fatigue prediction.
