Gear stress analysis: Lewis equation vs AGMA J and I factors

Gear teeth fail in two main ways: they crack at the root from repeated bending, or the tooth faces pit from repeated contact. The Lewis equation estimates the first. The AGMA 2001 method rates both, using two geometry factors, J and I, and fatigue-based allowable stresses.

Two ways a gear tooth fails

Which one limits a design depends on material and size. Hardened steel gears are usually limited by pitting; softer or smaller teeth more often by bending.

What is the Lewis bending equation?

In 1892 Wilfred Lewis treated a gear tooth as a cantilever beam with the whole load at its tip. The root bending stress is:

σ = Wt · P / (F · Y) metric: σ = Wt / (F · m · Y)

Wt is the tangential load at the pitch circle, P the diametral pitch (or m the module), F the face width, and Y the Lewis form factor, which depends only on the tooth shape. For 20° full-depth teeth (Shigley, Table 14-2):

TeethYTeethYTeethY
120.245200.320500.408
150.289260.3441000.446
170.302300.3583000.471
190.314400.389Rack0.484

Lewis is a quick sizing check, but it leaves out three things that matter: the stress concentration in the root fillet, the fact that a second tooth pair shares the load while the load is near the tip, and fatigue. It also says nothing about pitting.

How does AGMA calculate bending stress (the J factor)?

AGMA 2001 keeps Lewis's structure and adds the missing pieces:

σ = Wt · Ko · Kv · Ks · (P / F) · (Km · KB / J)

J is computed from the actual tooth (AGMA 908). The load sits at the highest point of single tooth contact, the highest point on the tooth where one pair carries the whole load, rather than at the tip. The critical section is where an inscribed Lewis parabola touches the root fillet. J then divides by a stress concentration factor for the fillet radius (Dolan and Broghamer's photoelastic results). Some AGMA 908 values for 20° full-depth gears cut with a 0.25-module tip radius:

Pinion / gear teethJ pinionJ gearI
21 / 210.330.330.078
21 / 350.340.370.091
26 / 550.370.410.101
35 / 1350.410.450.120
55 / 1350.450.470.112

J depends on the mating gear (it sets where single-tooth contact ends), on the cutter's tip radius, and on profile shift. Pinions below about 18 teeth at 20° are missing from the tables because they undercut; see undercut and profile shift.

How does AGMA calculate contact stress (the I factor)?

Contact stress comes from Hertz's theory of two curved surfaces pressed together:

σc = Cp · √( Wt · Ko · Kv · Ks · Km · Cf / (dP · F · I) )

dP is the pinion's pitch diameter. Cp is the elastic coefficient of the two materials, about 2300 √psi (191 √MPa) for steel on steel. I is the pitting geometry factor, which describes how curved the two profiles are where contact stress peaks. At the pitch point it has a simple form:

I = cos α · sin α / 2 · mG / (mG ± 1)

Here mG is the gear ratio (gear teeth over pinion teeth): + for external gears, − for internal. AGMA 908 evaluates the curvatures at the pinion's lowest point of single tooth contact instead, which gives slightly smaller and more accurate values. The minus sign for internal gears matters in planetary sets: a ring's concave teeth wrap around the planet, so I is much larger and contact stress lower than on the sun mesh.

How do you calculate a gear safety factor?

AGMA compares the stresses with fatigue strengths, not yield. For through-hardened steel the allowable stresses rise with Brinell hardness (Grade 1, 107 cycles, 99% reliability):

St = 77.3 HB + 12 800 psi Sc = 322 HB + 29 100 psi

The safety factors are SF = St / σ for bending and SH = Sc / σc for contact. Because contact stress grows with the square root of load, an SH of 1.4 is a margin of about 2 on load. Compare SF with SH2 to see which failure mode governs.

Plastic gears are rated differently. Nylon and acetal gear data from suppliers uses its own Lewis-based method with speed, temperature and lubrication factors, and is not compatible with AGMA allowables.

How are planetary gears rated differently?

Worked example

GearStudio's default set: sun 19, planets 26 (three of them), ring 71, 20 DP, 20°, 0.5 in face width, 50 in·lbf on the sun, Ko 1.25, Km 1.3, Kv 1.2, Ks 1.0. Sun and planets are 360 HB steel, the ring 150 HB.

  1. Tooth force per planet: 50 / (3 × 0.475) = 35.09 lbf, times Kγ 1.10 gives Wt = 38.60 lbf.
  2. Sun bending, J = 0.341: σ = 38.60 × 1.95 × 20 / (0.5 × 0.341) = 8.83 ksi. Allowable at 360 HB is 40.6 ksi, so SF = 4.60.
  3. Sun–planet contact, I = 0.0848, dP = 0.95 in: σc = 2291 × √(38.60 × 1.95 / (0.95 × 0.5 × 0.0848)) = 99.0 ksi. Allowable is 145.0 ksi, so SH = 1.46, or 2.14 on load.
  4. Planet bending, J = 0.368: 8.19 ksi against 70% of 40.6 = 28.4 ksi, so SF = 3.47.
  5. Ring: bending 6.15 ksi (SF 3.97); contact 53.6 ksi against 77.4 ksi, SH 1.44 (2.08 on load).

Contact governs: the smallest margin on load is 2.08, at the ring. Lewis alone would have reported a comfortable bending margin and missed that the tooth faces are the limit.

Planetary stress calculator

Ring fixed, sun driven, steel gears, standard centre distance. It uses the same rating code as GearStudio, with a 0.25-module cutter tip radius and no backlash, so J comes out a little lower than in the worked example above (which uses GearStudio's default fillet).

GearJBendingSFContactSH
Sun0.3349.02 ksi4.5099.0 ksi1.46
Planet0.3588.41 ksi3.3899.0 ksi1.46
Ring0.4706.40 ksi3.8153.6 ksi1.44
  • Lowest design factor 2.08 (ring contact), at or above the target of 2.

GearStudio runs this rating live as you design, with any materials, profile shift, fillet radius and backlash, and puts it in a printable engineering report.

Open GearStudio

This is a sizing method, not a certified AGMA or ISO 6336 rating. Check critical designs independently before manufacture.

Related guides

Undercut and profile shiftHow shift thickens the root and raises J. Planetary assembly conditionsTooth-count rules before you size for strength.