Circular Motion, Centrifugal Stress, and the Campbell Diagram
Table of Contents
- Rotation and Engineering: Engines, Turbines, Flywheels
- Rotational Kinematics: ω, α, θ
- Centripetal Acceleration and Centripetal Force
- The "Centrifugal Force" in Rotating Reference Frames
- Centrifugal Stress in Turbine Blades
- Flywheel Hoop Stress Design
- Gyroscopic Effects and Bearing Loads
- Campbell Diagram and Critical Speed
- Cross-Topics
1. Rotation and Engineering: Engines, Turbines, Flywheels
Rotating machinery is everywhere in engineering: jet engines, steam turbines, wind turbines, electric motors, gas compressors, automotive drivetrain components. The physics of circular motion — centripetal acceleration, rotational inertia, gyroscopic moments — creates structural loads that are unique to rotating systems and must be handled specifically in FEM analysis.
Jet engine turbine blades spin at like 10,000 RPM — what kind of stress does that actually create in the blade? I've heard it's enormous.
Enormous is right. At 10,000 RPM, a typical high-pressure turbine blade root carries centrifugal stress of 200–300 MPa from rotation alone — before adding the aerodynamic loads and thermal stresses from 1400°C gas temperature. The blade is essentially hanging by its root like a heavy weight on a very fast merry-go-round. The centrifugal force is the blade mass times the centripetal acceleration: at 10,000 RPM and a blade radius of 0.3m, the centripetal acceleration is ω²r = (1047 rad/s)² × 0.3 ≈ 328,000 m/s² — about 33,000g. That's thirty-three thousand times gravity acting outward on every particle of the blade.
2. Rotational Kinematics: ω, α, θ
Rotational motion is described by angular position $\theta$, angular velocity $\omega$, and angular acceleration $\alpha$:
For constant angular acceleration, the rotational analogs of the SUVAT equations apply. The relationship between rotational and translational quantities at radius $r$:
Converting RPM to rad/s: $\omega = 2\pi N/60$ where $N$ is in RPM.
| Machine | Operating Speed | ω (rad/s) | Tip Speed (m/s) |
|---|---|---|---|
| Wind turbine (large) | 5–15 RPM | 0.5–1.6 | 60–90 (tip) |
| Car engine (cruise) | 2000 RPM | 209 | ~10 (crankpin) |
| Electric motor (2-pole, 50Hz) | 3000 RPM | 314 | Varies |
| Centrifugal compressor | 20,000 RPM | 2094 | ~400 (impeller) |
| Gas turbine HP stage | 10,000–20,000 RPM | 1047–2094 | 300–500 |
| Dental drill | 300,000 RPM | 31,416 | ~50 |
3. Centripetal Acceleration and Centripetal Force
For uniform circular motion (constant $\omega$, constant $r$), the centripetal acceleration points inward (toward the rotation axis):
The centripetal force required to maintain circular motion:
This force must be provided by some real physical mechanism — the tension in a string for a ball on a string, the normal force from a banked road for a car in a curve, the blade root attachment for a turbine blade, the rim for a flywheel. FEM analysis of rotating components finds the stress distribution created by this centripetal force requirement.
4. The "Centrifugal Force" in Rotating Reference Frames
In an inertial (non-rotating) reference frame, centrifugal force doesn't exist — there's only the centripetal force directed inward. The centrifugal force appears when you analyze problems in a rotating reference frame (attached to the rotating body). In that frame, a fictitious outward force $F_{\text{centrifugal}} = m\omega^2 r$ seems to push everything outward.
For structural analysis of rotating components, the rotating frame is often more convenient: the blade is stationary in the rotating frame, and the centrifugal body force $\mathbf{b} = \rho\omega^2 r\hat{r}$ is applied as a distributed load. This is exactly how FEM analyses of rotating machinery are set up:
*DLOAD
BLADE_SET, CENT, 1.096e6, 0., 0., 1., 0., 0., 0.
! CENT (centrifugal) = ρω² applied as body force density
! 1.096e6 = ρ × ω² = 7900 × (1047 rad/s)²
5. Centrifugal Stress in Turbine Blades
For a rotating disk or blade attachment region, the centrifugal stress is found by integrating the centrifugal body force. For a simple uniform bar of cross-section $A$, density $\rho$, rotating at $\omega$ with tip at radius $R$:
Where $L = R - r_{\text{root}}$ is the blade length. For a typical HPT blade (Nickel superalloy, $\rho = 8400$ kg/m³, $\omega = 1047$ rad/s, $L = 0.10$ m, hub at $r = 0.25$ m):
And that's the centrifugal component alone. Add thermal stresses (50–100 MPa from temperature gradients) and aerodynamic bending stresses (30–60 MPa): total blade root stress approaches 400 MPa. The material yield at temperature (~950°C) for advanced Ni-superalloy is only ~600–700 MPa. Safety factor ≈ 1.5 — very tight.
If the stress is so high, why don't turbine blades just fly off or break during normal operation?
Because the design is exquisitely optimized — and that's where FEM earns its keep. Turbine blades are made from single-crystal nickel superalloys that have no grain boundaries to initiate creep cracks at high temperature. They're hollow with internal cooling channels that keep the blade metal below 1050°C even in 1400°C gas. The fir-tree root geometry distributes the centrifugal load over multiple lobes to minimize stress concentration. And every blade is individually inspected with X-ray CT scanning. The safety margins are intentionally small to maximize thermodynamic efficiency — any extra metal in the blade adds centrifugal load faster than it adds strength. So these blades are designed right to the limit of what's physically possible.
6. Flywheel Hoop Stress Design
A flywheel is a thin rotating disk designed to store kinetic energy. The dominant stress in a rotating disk is the hoop (circumferential) stress from centrifugal loading:
For a solid disk with $r=0$, the maximum hoop stress is at the center:
The specific energy stored per unit mass: $E/m = \frac{1}{2}\frac{I\omega^2}{m} = \frac{1}{4}\omega^2 R^2$. The maximum speed is limited by $\sigma_{\theta,\max} \leq \sigma_y$, giving:
This means flywheel energy density is determined purely by the ratio $\sigma_y/\rho$ (specific strength). Carbon fiber composite ($\sigma_y/\rho \approx 1000$ kJ/kg) outperforms steel ($\sigma_y/\rho \approx 50$ kJ/kg) by 20× — why CFRP flywheels are used in high-performance energy storage systems.
7. Gyroscopic Effects and Bearing Loads
A spinning rotor resists changes in its rotation axis orientation — gyroscopic effect. When an external torque $M$ changes the direction of angular momentum $\mathbf{L}$, the rotor precesses:
For a jet engine mounted under an aircraft wing, gyroscopic couples arise during maneuvers (pitch and yaw). These create significant bearing loads that must be included in engine mount FEM analysis. For a turbofan with rotor angular momentum $L = I\omega = 500$ kg·m²/s, a 2°/s yaw rate creates a gyroscopic moment: $M = 500 \times 0.035 = 17.5$ kN·m — entirely from gyroscopic effect, before any aerodynamic loads.
8. Campbell Diagram and Critical Speed
The Campbell diagram is the primary tool for resonance avoidance design in rotating machinery. It plots natural frequencies (horizontal lines) against rotation speed (x-axis), with engine order lines as diagonals:
Intersections of horizontal (natural frequency) lines with diagonal (engine order) lines are resonance crossings — potentially dangerous operating speeds.
Critical Speed of a Rotor-Bearing System
The critical speed of a simple rotor on two bearings is the speed at which rotor eccentricity (unbalance) excites its own natural frequency. Below critical speed: rotor deflects toward unbalance; above critical: rotor deflects away from unbalance (self-centering). Modern turbomachinery often operates supercritically — above one or more critical speeds — with carefully designed squeeze-film damper bearings to safely traverse the critical speed during run-up.
I'm doing FEM of a centrifugal compressor impeller. How do I set up the rotational analysis correctly in Abaqus?
Several steps: First, apply a steady-state spinning analysis (*STATIC with NLGEOM for stress stiffening) to get the centrifugal stress field. Use a rotating reference frame with CENTRIFUGAL load. Second, run a prestressed modal analysis (*FREQUENCY, PERTURBATION) starting from the spun-up stress state — rotation stiffens blades through geometric stiffening, so natural frequencies at speed are higher than at rest (the "spin stiffening" or centrifugal stiffening effect can raise blade frequencies by 20–50% at operating speed). Third, plot the natural frequencies vs. speed to create your Campbell diagram and check for resonance crossings at sustained operating speeds. For variable-speed machines, identify which orders could cross your frequencies in the operating range and verify those crossings have high damping or happen during transient run-through, not at steady state.
9. Cross-Topics
| Topic | Connection | Link |
|---|---|---|
| Newton's Laws | Centripetal force = Newton's 2nd Law applied to circular path | Newton's Laws of Motion |
| Simple Harmonic Motion | SHM is the projection of circular motion onto a line | Simple Harmonic Motion |
| Momentum & Impulse | Angular momentum conservation governs gyroscopic behavior | Momentum and Impulse |
| Stress & Strain | Centrifugal body force creates hoop and radial stresses in rotating disks | Stress and Strain Basics |
| Vibration & Dynamics | Rotating machinery dynamics, gyroscopic matrix, critical speed analysis | Vibration & Dynamics |