9.1 Transmission Shaft Design, Critical Speeds & Deflection

9.1 Transmission Shaft Design, Critical Speeds & Deflection

Transmission shafts serve as the foundational rotating conduits for mechanical power across gearboxes, drivetrains, turbines, and industrial machinery. Sizing a reliable shaft requires simultaneously satisfying multiaxial fatigue endurance criteria, static yield thresholds, strict elastic deflection and angular slope tolerances at bearing and gear interfaces, and dynamic clearance from lateral whirling critical speeds.

---

1. Shaft Loading Mechanics: Cyclic Bending & Torsion

A rotating shaft supporting power transmission elements (spur gears, helical gears, pulleys, sprockets) is subjected to three-dimensional force systems that induce simultaneous bending moments, torsional moments, and transverse or axial shear loads.

          [Bearing A]         [Gear / Pulley]         [Bearing B]
               |                     |                     |
         +=====|=====================O=====================|=====+
         |     |                   / | \                   |     |
         +=====|=====================O=====================|=====+
             ^                     |   |                     ^
             |                     v   v                     |
           R_Ay                   F_t  F_r                 R_By

  Rotational Speed: ω (rpm)
  Bending Moment: M (reverses sign every 180° rotation -> completely reversed fatigue)
  Torque: T (steady transmission -> constant mean shear stress)

Fully Reversed Bending Stress

When a shaft rotates at angular velocity $\omega$ under a stationary transverse load (such as gear pitch-line forces or gravity), any material element located on the outer surface rotates from maximum tension at the bottom outer fiber to maximum compression at the top outer fiber. Thus, steady external transverse forces produce completely reversed cyclic bending stresses:

$$M_a = \sqrt{M_x^2 + M_y^2}, \quad M_m = 0$$ $$\sigma_a = K_f \frac{32 M_a}{\pi d^3}, \quad \sigma_m = K_f \frac{32 M_m}{\pi d^3} = 0$$

Where $K_f$ is the fatigue stress concentration factor for bending at geometric transitions (shoulders, keyways, retaining ring grooves).

Torsional Shear Stress

In steady-state power transmission, the delivered motor torque is predominantly constant:

$$T_m = T = \frac{P}{\omega} = \frac{63025 H}{n}, \quad T_a \approx 0$$ $$\tau_m = K_{fs} \frac{16 T_m}{\pi d^3}, \quad \tau_a = K_{fs} \frac{16 T_a}{\pi d^3} \approx 0$$

Where $K_{fs}$ is the torsional fatigue stress concentration factor. For reciprocating machinery (internal combustion engines, piston compressors), an alternating torque component $T_a$ exists and must be evaluated alongside $T_m$.

---

2. Multiaxial Fatigue Failure Criteria: The ASME Shaft Equation

Because the bending normal stress $\sigma$ and torsional shear stress $\tau$ act simultaneously on the shaft surface, they produce a biaxial plane stress state. Applying the Distortion Energy (DE / von Mises) theory separates the stress state into equivalent alternating and mean von Mises stresses:

$$\sigma_a' = \sqrt{\sigma_a^2 + 3 \tau_a^2} = \frac{16}{\pi d^3} \sqrt{4 (K_f M_a)^2 + 3 (K_{fs} T_a)^2}$$ $$\sigma_m' = \sqrt{\sigma_m^2 + 3 \tau_m^2} = \frac{16}{\pi d^3} \sqrt{4 (K_f M_m)^2 + 3 (K_{fs} T_m)^2}$$

DE-Goodman Criterion

Coupling the equivalent stresses with the Modified Goodman fatigue failure line yields:

$$\frac{\sigma_a'}{S_e} + \frac{\sigma_m'}{S_{ut}} = \frac{1}{n}$$

Substituting the expressions for $\sigma_a'$ and $\sigma_m'$ and solving explicitly for the shaft diameter $d$ gives the DE-Goodman shaft equation (the form used in Shigley; the withdrawn ANSI/ASME B106.1M shaft standard used the elliptic criterion instead):

$$d = \left[ \frac{16 n}{\pi} \left\{ \frac{1}{S_e}\left(4(K_f M_a)^2 + 3(K_{fs} T_a)^2\right)^{1/2} + \frac{1}{S_{ut}}\left(4(K_f M_m)^2 + 3(K_{fs} T_m)^2\right)^{1/2} \right\} \right]^{1/3}$$

For the common industrial case of purely reversed bending ($M_m = 0$) and purely steady torque ($T_a = 0$), this simplifies to:

$$d = \left[ \frac{16 n}{\pi} \left( \frac{2 K_f M_a}{S_e} + \frac{\sqrt{3} K_{fs} T_m}{S_{ut}} \right) \right]^{1/3}$$

DE-Gerber and Static Yield Check

The DE-Gerber criterion replaces the linear Goodman boundary with a parabolic relation:

$$\frac{n \sigma_a'}{S_e} + \left( \frac{n \sigma_m'}{S_{ut}} \right)^2 = 1$$

Gerber models the mean of experimental failure data for ductile steels more closely, whereas Goodman is a conservative lower envelope that is widely used for design. In all designs, a first-cycle yield check (Langer / von Mises static line) must confirm that transient startup torques do not cause localized plastic deformation:

$$\sigma_{\max}' = \sqrt{(\sigma_a + \sigma_m)^2 + 3(\tau_a + \tau_m)^2} \le \frac{S_y}{n_y}$$

---

3. Shaft Deflection & Angular Slope Limits

High-strength alloy steels increase fatigue allowable stresses ($S_e, S_{ut}$) without altering Young's modulus ($E \approx 207\text{ GPa}$ for all carbon and low-alloy steels). Consequently, shafts designed purely for fatigue strength often exhibit excessive elastic deflections that degrade gear mesh kinematics and destroy bearings.

(Typical ranges from Shigley Table 7-2; bearing catalogs give the exact misalignment limits for each bearing series. Helical, bevel and crowned gears have their own manufacturer limits.)

---

4. Shaft Whirling & Critical Speeds

Dynamic unbalance occurs because the mass center of any real rotor is displaced from its geometric axis of rotation by a small radial eccentricity $e$. As the shaft rotates at angular velocity $\omega$, centrifugal forces bow the shaft outward against its elastic flexural stiffness $k$:

         Elastic Centerline (Restoring Force: k*r)
                   |
                   v
       [ Bearing ] o---------*---------*---------o [ Bearing ]
                             |         |
                         Geometric   Center
                          Center    of Mass
                            (r)     (r + e)

  Centrifugal Force: F_c = m * ω^2 * (r + e)
  Equating F_c = k * r yields whirl amplitude: r = e * (ω / ω_n)^2 / [1 - (ω / ω_n)^2]

When operating speed $\omega$ matches the lateral natural frequency $\omega_n = \sqrt{k/m}$, resonance occurs and the theoretical whirl amplitude approaches infinity. This natural frequency is the critical whirling speed $\omega_{cr}$.

Rayleigh-Ritz Method for Lumped Masses

For a shaft supporting multiple lumped masses $w_i = m_i g$ with static deflections $y_i$ caused by self-weight acting in the direction of vibration:

$$\omega_{cr} = \sqrt{\frac{g \sum_{i=1}^n w_i y_i}{\sum_{i=1}^n w_i y_i^2}}$$

For a single concentrated center load on a simple span, this reduces to $\omega_{cr} = \sqrt{g / y_{\text{static}}}$.

Dunkerley's Equation

Dunkerley's method provides a lower-bound (conservative) estimate of the fundamental critical speed by combining the bare shaft natural frequency $\omega_s$ with the critical speeds associated with each individual rotor mass $\omega_i$ acting alone on a massless shaft:

$$\frac{1}{\omega_{cr}^2} = \frac{1}{\omega_s^2} + \frac{1}{\omega_1^2} + \frac{1}{\omega_2^2} + \dots + \frac{1}{\omega_n^2}$$

Sub-Critical vs. Super-Critical Operation

• Sub-Critical Operation ($\omega < 0.75\omega_{cr}$): Standard for machine tool spindles, gearboxes, and industrial drives. The shaft remains rigid, with whirl amplitude remaining below eccentricity $e$.
• Super-Critical Operation ($\omega > 1.25\omega_{cr}$ to $1.4\omega_{cr}$): Common in high-speed turbomachinery (steam turbines, centrifuges). Above critical resonance, a $180^\circ$ phase lag occurs between unbalance and deflection; the rotor self-centers about its center of mass ($r \to -e$). Startup systems must accelerate rapidly through the critical speed with adequate external damping to prevent catastrophic deflection accumulation.

---

5. Complete Worked Shaft Design & Critical Speed Analysis

Problem Statement

A solid drive shaft fabricated from AISI 1045 cold-drawn steel ($S_{ut} = 630\text{ MPa}$, $S_y = 530\text{ MPa}$, $E = 207\text{ GPa}$) is supported by deep-groove ball bearings across a simple span of $L = 600\text{ mm}$. The shaft transmits $P = 30\text{ kW}$ at $n = 1200\text{ rpm}$ under uniform, non-pulsating torque ($T_a = 0$).

A spur gear mounted at mid-span ($x = 300\text{ mm}$) introduces a steady torque and a completely reversed bending moment of $M_a = 420\text{ N}\cdot\text{m}$ ($M_m = 0$). The total lumped rotor mass at mid-span is $m = 32\text{ kg}$.

At the gear mounting profile keyway, fatigue stress concentration factors are $K_f = 2.0$ (bending) and $K_{fs} = 1.6$ (torsion). The endurance limit corrected for size, surface finish, and reliability ($99\%$) is $S_e = 210\text{ MPa}$. The target design safety factor is $n = 2.0$.

Step 1: Calculate Transmitted Torque

$$\omega = \frac{1200 \times 2\pi}{60} = 125.66\text{ rad/s}$$ $$T_m = \frac{P}{\omega} = \frac{30000\text{ W}}{125.66\text{ rad/s}} = 238.73\text{ N}\cdot\text{m}$$

Step 2: Determine Minimum Diameter via the DE-Goodman Equation

Applying the simplified DE-Goodman shaft equation for $M_m = 0$ and $T_a = 0$:

$$d = \left[ \frac{16 n}{\pi} \left( \frac{2 K_f M_a}{S_e} + \frac{\sqrt{3} K_{fs} T_m}{S_{ut}} \right) \right]^{1/3}$$

Evaluate the individual terms:

$$\frac{2 K_f M_a}{S_e} = \frac{2 \times 2.0 \times 420\text{ N}\cdot\text{m}}{210 \times 10^6\text{ Pa}} = \frac{1680}{210 \times 10^6} = 8.00 \times 10^{-6}\text{ m}^3$$ $$\frac{\sqrt{3} K_{fs} T_m}{S_{ut}} = \frac{1.732 \times 1.6 \times 238.73\text{ N}\cdot\text{m}}{630 \times 10^6\text{ Pa}} = \frac{661.60}{630 \times 10^6} = 1.050 \times 10^{-6}\text{ m}^3$$ $$\text{Sum} = (8.00 + 1.050) \times 10^{-6} = 9.050 \times 10^{-6}\text{ m}^3$$ $$\frac{16 n}{\pi} = \frac{16 \times 2.0}{\pi} = 10.186$$ $$d^3 = 10.186 \times 9.050 \times 10^{-6} = 9.218 \times 10^{-5}\text{ m}^3$$ $$d = (9.218 \times 10^{-5})^{1/3} = 0.04517\text{ m} = 45.2\text{ mm}$$

Select standard metric shaft diameter: $d = 50\text{ mm} = 0.050\text{ m}$.

Step 3: First-Cycle Static Yield Verification at $d = 50\text{ mm}$

$$\sigma_a = \frac{32 K_f M_a}{\pi d^3} = \frac{32 \times 2.0 \times 420}{\pi (0.050)^3} = 68.45\text{ MPa}$$ $$\tau_m = \frac{16 K_{fs} T_m}{\pi d^3} = \frac{16 \times 1.6 \times 238.73}{\pi (0.050)^3} = 15.56\text{ MPa}$$ $$\sigma_{\max}' = \sqrt{\sigma_a^2 + 3 \tau_m^2} = \sqrt{(68.45)^2 + 3(15.56)^2} = \sqrt{4685.4 + 726.3} = 73.56\text{ MPa}$$ $$n_y = \frac{S_y}{\sigma_{\max}'} = \frac{530\text{ MPa}}{73.56\text{ MPa}} = 7.21 \gg 1.0 \quad \text{(Safe against static yielding)}$$

Step 4: Critical Whirling Speed & Bearing Slope Verification

Calculate shaft cross-sectional area moment of inertia:

$$I = \frac{\pi d^4}{64} = \frac{\pi (0.050)^4}{64} = 3.068 \times 10^{-7}\text{ m}^4$$ $$E I = (207 \times 10^9\text{ Pa}) \times (3.068 \times 10^{-7}\text{ m}^4) = 63507\text{ N}\cdot\text{m}^2$$

Rotor weight acting at mid-span:

$$W = m g = 32\text{ kg} \times 9.81\text{ m/s}^2 = 313.92\text{ N}$$

Mid-span static deflection under rotor weight:

$$y_{\text{static}} = \frac{W L^3}{48 E I} = \frac{313.92 \times (0.600)^3}{48 \times 63507} = \frac{67.807}{3048336} = 2.224 \times 10^{-5}\text{ m} = 0.02224\text{ mm}$$

Fundamental critical whirling speed:

$$\omega_{cr} = \sqrt{\frac{g}{y_{\text{static}}}} = \sqrt{\frac{9.81\text{ m/s}^2}{2.224 \times 10^{-5}\text{ m}}} = \sqrt{441097} = 664.15\text{ rad/s}$$ $$N_{cr} = \frac{664.15 \times 60}{2\pi} = 6342\text{ rpm}$$ $$\frac{N_{\text{op}}}{N_{cr}} = \frac{1200\text{ rpm}}{6342\text{ rpm}} = 0.189 \quad (18.9\% \text{ of critical speed})$$

The operating speed is far below the $75\%$ threshold, confirming rigid sub-critical operation free of whirling resonance.

Check bearing seat slope $\theta$ at the supports:

$$\theta_{\text{bearing}} = \frac{W L^2}{16 E I} = \frac{313.92 \times (0.600)^2}{16 \times 63507} = \frac{113.01}{1016112} = 1.112 \times 10^{-4}\text{ rad} = 0.000111\text{ rad}$$

Because $0.000111\text{ rad} \ll 0.001\text{ rad}$ (the allowable alignment limit for deep-groove ball bearings), the shaft satisfies all stiffness and operational criteria.