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.
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).
Critical Speed Calculator
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 Supported | 48 | 1.0x (baseline) |
| Fixed-Fixed | 192 | 4.0x stiffer |
| Fixed-Pinned | 107 | 2.2x stiffer |
| Cantilever | 3 | 0.06x (flexible) |
Operating Guidelines
| Speed Range | Status | Action |
|---|---|---|
| Below 75% Nc | SAFE | Normal subcritical operation |
| 75% - 80% Nc | CAUTION | Increased vibration expected |
| 80% - 120% Nc | DANGER | Avoid steady operation |
| 120% - 125% Nc | CAUTION | Transition zone |
| Above 125% Nc | SAFE | Supercritical operation |