Shaft Critical Speed Calculator

Calculate the critical speed of rotating shafts using Rayleigh's energy method and Dunkerley's equation for systems with multiple disks. Includes support configuration analysis, operating speed margin checker, and mode shape visualization for pumps, motors, and turbines.

How It Works - Shaft Critical Speed Theory

Understanding Shaft Whirl

When a rotating shaft operates at certain speeds, it can experience dangerous lateral vibrations called whirl. This occurs when the rotational speed matches one of the shaft's natural frequencies, causing resonance.

  • Synchronous Whirl: The shaft deflection rotates at the same speed as the shaft - most common type
  • Subsynchronous Whirl: Deflection rotates slower than shaft speed - often caused by fluid film effects in bearings
  • Supersynchronous Whirl: Deflection rotates faster than shaft - less common, seen in certain rotor configurations

At critical speed, even small mass imbalances cause the shaft to bow outward, creating a self-reinforcing vibration that can lead to catastrophic failure.

Rayleigh Method (Upper Bound)

The Rayleigh method uses energy principles to estimate the first natural frequency. It equates maximum potential energy (at maximum deflection) with maximum kinetic energy (at equilibrium position):

omega = sqrt(g * Sum(m_i * y_i) / Sum(m_i * y_i^2))

  • Requires knowledge of static deflection at each mass location
  • Provides an upper bound estimate - actual critical speed is at or below this value
  • More accurate when the assumed deflection curve closely matches the actual first mode shape
  • Best results when using deflections from gravity loading or first mode approximation

Dunkerley's Method (Lower Bound)

Dunkerley's equation combines individual critical speeds to estimate the system critical speed:

1/Nc^2 = 1/Nc_shaft^2 + Sum(1/Nc_i^2)

  • Each Nc_i is the critical speed if only that mass were present
  • Provides a lower bound estimate - actual critical speed is at or above this value
  • Very useful when individual component critical speeds are known
  • Accuracy improves when masses are well-separated along the shaft

Support Conditions and Their Effects

The boundary conditions at shaft ends dramatically affect stiffness and critical speed:

  • Simply Supported (Pinned-Pinned): Both ends free to rotate but not translate. Common for between-bearing rotors. K = 48 for center load.
  • Fixed-Fixed: Both ends fully constrained. Higher stiffness, higher critical speed. K = 192 for center load.
  • Cantilever (Fixed-Free): One end fixed, other free. Overhung designs like pump impellers. K = 3 for end load - lowest stiffness.
  • Fixed-Pinned: One end fixed, other pinned. Intermediate case. K = 107 for center load.

Operating Below vs Above Critical Speed

  • Subcritical Operation (below 75% of Nc):
    • Shaft acts as a stiff rotor - deflection follows heavy side
    • Unbalance response increases as speed approaches critical
    • Common for most industrial machinery
    • Simpler balancing requirements
  • Transition Zone (75-125% of Nc):
    • DANGER ZONE - Avoid steady-state operation
    • Maximum vibration amplitude occurs here
    • Must pass through quickly during startup/shutdown
    • High stress, potential for damage
  • Supercritical Operation (above 125% of Nc):
    • Shaft becomes flexible rotor - deflection moves to opposite side of heavy spot
    • Self-centering effect reduces vibration
    • Used in high-speed turbomachinery, centrifuges
    • Requires careful design and balancing

Higher Mode Critical Speeds

Beyond the first critical speed, shafts have additional mode shapes and critical speeds:

  • 1st Mode: Single bow - shaft bends in one half-wave
  • 2nd Mode: S-shape - shaft has one node point where deflection is zero
  • 3rd Mode: Double S-shape - two node points

For simply supported shafts, higher mode frequencies are approximately: 2nd = 4x, 3rd = 9x, 4th = 16x the first mode (n^2 relationship).

Shaft Deflection Mode Shapes
1st MODE (Fundamental) y_max omega Single bow deflection Maximum amplitude at center No node points 1st Critical Speed Nc1 = (pi^2/L^2)*sqrt(EI/rho*A) 2nd MODE Node S-shaped deflection One node at center Opposite phase each side 2nd Critical Speed Nc2 = 4 x Nc1 (approx) 1st Mode 2nd Mode Node Point Undeflected

Critical Speed Calculator

Simply Supported
Bearings at both ends
K = 48
Fixed-Fixed
Both ends rigidly mounted
K = 192
Cantilever
Overhung / one end free
K = 3
Fixed-Pinned
One fixed, one pinned
K = 107
Dunkerley's Method: Add multiple disks along the shaft. Each disk's contribution is combined to find the system critical speed (lower bound estimate).
Operating Speed Margin Checker: Enter your calculated critical speed and operating speed to verify safe separation margins.
Leave blank for estimate
-- RPM
1st Critical Speed
--
Subcritical Zone Danger Zone Supercritical Zone
Nc
0 RPM Nc 2x Nc
Natural Frequency
--
Angular Frequency
--
Static Deflection
--
Shaft Stiffness
--
Max Subcritical
--
Min Supercritical
--
2nd Mode (Est.)
--
Method
--

Key Formulas

Single Mass (Center Load):

Nc = (60/2*pi) * sqrt(K*E*I / (m*L^3))

Rayleigh Method:

omega = sqrt(g*Sum(m*y) / Sum(m*y^2))

Dunkerley's Equation:

1/Nc^2 = 1/Nc_s^2 + Sum(1/Nc_i^2)

Support Configuration Factors

Configuration K Factor Relative Stiffness
Simply Supported481.0x (baseline)
Fixed-Fixed1924.0x stiffer
Fixed-Pinned1072.2x stiffer
Cantilever30.06x (flexible)

Operating Guidelines

Speed Range Status Action
Below 75% NcSAFENormal subcritical operation
75% - 80% NcCAUTIONIncreased vibration expected
80% - 120% NcDANGERAvoid steady operation
120% - 125% NcCAUTIONTransition zone
Above 125% NcSAFESupercritical operation