Flywheel (Solid Rotating Disk)

energy-storagestressmass-cost

Verified build 8 relations · 7 identities proven · 2 modeling steps · 6 parity samples

A flywheel is a battery made of motion: spin mass up to store work, let it back down to release it. Punch presses bank energy between strokes, engines smooth the gaps between cylinder firings, satellites trade spin with reaction wheels to point themselves, and grid-scale flywheel farms soak up seconds-long power surges. The catch is that the flywheel loads itself — every gram needs a centripetal pull of ρω2r\rho\omega^2 r per unit volume to stay on its circle, and that pull has to come through the disk as stress:

Ek=12Izω2σmax=3+ν8ρω2R2emax=2σy(3+ν)ρE_k = \tfrac{1}{2} I_z\,\omega^2 \qquad \sigma_{max} = \frac{3+\nu}{8}\,\rho\,\omega^2 R^2 \qquad e_{max} = \frac{2\,\sigma_y}{(3+\nu)\,\rho}

Three things to notice, one per formula:

The energy-in configuration runs the same relations the other way: name the energy you need, and the widget finds the speed that stores it — then tells you what that speed costs in stress and margin. Watch ωy\omega_y, the yield-onset speed: it is the model’s own ceiling, and it corresponds to a fixed rim speed ωyR=8σy/((3+ν)ρ)\omega_y R = \sqrt{8\sigma_y/((3+\nu)\rho)}, another pure material number.

How long does charging take? Spinning the disk up from rest with a constant drive torque TdT_d takes time, and angular impulse–momentum sets the clock — the torque’s angular impulse TdtT_d\,t must supply the angular momentum IzωI_z\,\omega the flywheel carries at speed:

tspin=IzωTdt_{spin} = \frac{I_z\,\omega}{T_d}

More torque spins it up faster; more inertia — a bigger or denser disk — or a higher target speed draws it out (swap in a lighter alloy and tspint_{spin} falls with the inertia, invariant 3 again). This is the number the chain builder’s headline example holds in tension with shaft stress: the gear ratio ahead of the flywheel sets TdT_d, so cranking the ratio up shortens tspint_{spin} while winding the driveshaft that carries TdT_d correspondingly tighter — storing energy quickly is not free. A negative TdT_d brakes rather than accelerates, so the widget withholds tspint_{spin} while the energy, stress, and margin readouts still stand — the spin-up time is the only thing that depends on how the disk is driven. (At exactly Td=0T_d = 0 the formula divides by zero, and the whole readout set refuses.)

Try it

Material

T3, bare flat sheet 0.010-0.128 in. thick, AMS 4037 / AMS-QQ-A-250/4 (MIL-HDBK-5J Table 3.2.3.0(b1), p. 3-71)

Bound properties of 2024-T3 aluminum sheet (bare)
nu0.33 1typicalmil-hdbk-5j
rho0.1 lb/inch**3typicalmil-hdbk-5j
sigma_y47 ksidesign min.mil-hdbk-5j
Inputs
N·m
Disk mass
kg
Moment of inertia (spin axis)
kg·m²
Stored kinetic energy
Specific energy
Peak stress (disk centre)
Safety factor (first yield)
Yield-onset speed
Spin-up time (from rest to ω)

3 materials in the database are not listed here: no published value in our cited sources for every property this THING needs.

Materials modeled here: 2024-T3 aluminum sheet (bare) 304 stainless steel 6061-T6 aluminum 7075-T6 aluminum AISI 1045 medium-carbon steel AISI 4340 low-alloy steel (Ni-Cr-Mo) ASTM A36 structural steel (hot-rolled) C26000 Cartridge Brass (70/30) Nylon 6/6 (PA66), unfilled Ti-6Al-4V

Governing relations

m=ρπR2tm = \rho\,\pi R^2 t

Assumes: uniform solid disk — constant thickness, no hub, no rim band, no central hole

Source: Hibbeler, R. C., Engineering Mechanics: Dynamics, 14th ed., Pearson, 2016 — §17.1 (mass moments of inertia; the uniform solid disk I = mR²/2), ch. 18 (planar kinetics, work and energy: E = ½Iω² for rotation about a fixed axis), and ch. 19 (planar kinetics, impulse and momentum: the principle of angular impulse and momentum, I_O·ω₁ + Σ∫M_O dt = I_O·ω₂, §19.1–19.2) — the basis for the constant-torque spin-up time T_d·t = I_z·ω.

Iz=12mR2I_z = \tfrac{1}{2} m R^2

Assumes: rigid uniform disk spinning about its axis of symmetry

Source: Hibbeler, R. C., Engineering Mechanics: Dynamics, 14th ed., Pearson, 2016 — §17.1 (mass moments of inertia; the uniform solid disk I = mR²/2), ch. 18 (planar kinetics, work and energy: E = ½Iω² for rotation about a fixed axis), and ch. 19 (planar kinetics, impulse and momentum: the principle of angular impulse and momentum, I_O·ω₁ + Σ∫M_O dt = I_O·ω₂, §19.1–19.2) — the basis for the constant-torque spin-up time T_d·t = I_z·ω.

Ek=12Izω2E_k = \tfrac{1}{2} I_z \omega^2

Assumes: rigid body in pure spin about a fixed axis — no wobble, no translation

Source: Hibbeler, R. C., Engineering Mechanics: Dynamics, 14th ed., Pearson, 2016 — §17.1 (mass moments of inertia; the uniform solid disk I = mR²/2), ch. 18 (planar kinetics, work and energy: E = ½Iω² for rotation about a fixed axis), and ch. 19 (planar kinetics, impulse and momentum: the principle of angular impulse and momentum, I_O·ω₁ + Σ∫M_O dt = I_O·ω₂, §19.1–19.2) — the basis for the constant-torque spin-up time T_d·t = I_z·ω.

e=Ekme = \frac{E_k}{m}

Assumes: the flywheel figure of merit — energy stored per kilogram carried

Source: Ashby, M. F., Materials Selection in Mechanical Design, 4th ed., Butterworth-Heinemann, 2011 — the "materials for flywheels" case study: maximum kinetic energy per unit mass is governed by the material index σ_y/ρ, independent of the disk's size.

σmax=3+ν8ρω2R2\sigma_{max} = \frac{3+\nu}{8}\,\rho\,\omega^2 R^2

Assumes: linear elastic, isotropic, uniform-thickness disk in plane stress (thin next to its radius); free outer rim, solid centre — NO central hole (even a small bore roughly doubles the peak); radial and hoop stress are equal and maximal at the centre, falling monotonically to the rim · Valid while: The centre of the disk has reached yield — beyond this speed the elastic field (and every number built on it) stops being the truth. First yield is NOT burst: a ductile disk redistributes stress and survives well past this line, but this model ends here. Thickness is no longer small next to the radius — the plane-stress thin-disk solution is fading; a disk this chunky wants the axisymmetric thick-disk treatment.

Source: Timoshenko, S. P., & Goodier, J. N., Theory of Elasticity, 3rd ed., McGraw-Hill, 1970 — ch. 4 (two-dimensional problems in polar coordinates), §32 "Rotating Disks": plane-stress field of a solid uniform disk, σ_r = σ_θ = (3+ν)ρω²R²/8 at the centre, σ_r(R) = 0 at the free rim.

SF=σyσmax\mathrm{SF} = \frac{\sigma_y}{\sigma_{max}}

Assumes: margin against first yield at the disk centre, not against burst — burst arrives later, governed by ultimate strength and crack behaviour

Source: Budynas, R. G., & Nisbett, J. K., Shigley's Mechanical Engineering Design, 10th ed., McGraw-Hill, 2015 — ch. 16 (flywheels: energy exchange and sizing) and ch. 1 (design factor against first yield).

σy=3+ν8ρωy2R2\sigma_y = \frac{3+\nu}{8}\,\rho\,\omega_y^2 R^2

Assumes: the spin speed at which the centre first reaches yield — the elastic model's own ceiling

Source: Timoshenko, S. P., & Goodier, J. N., Theory of Elasticity, 3rd ed., McGraw-Hill, 1970 — ch. 4 (two-dimensional problems in polar coordinates), §32 "Rotating Disks": plane-stress field of a solid uniform disk, σ_r = σ_θ = (3+ν)ρω²R²/8 at the centre, σ_r(R) = 0 at the free rim.

tspin=IzωTdt_{spin} = \frac{I_z\,\omega}{T_d}

Assumes: a constant drive torque T_d accelerates the disk from rest (ω = 0) up to the spin speed ω; angular impulse–momentum about the spin axis: the angular impulse ∫T dt equals the change in angular momentum I_z·Δω, so a constant T_d gives T_d·t = I_z·ω and t = I_z·ω / T_d; the drive torque alone sets the time — bearing drag, windage, and the disk's own stress field do not enter it (an idealization; real spin-up against losses takes longer) · Valid while: A non-positive drive torque never brings the disk up to speed, so the spin-up time is withheld. For a NEGATIVE (braking) torque only the spin-up time goes — the flywheel's stored energy, stress, and yield margin do not depend on how it is driven, so those readouts still stand. At EXACTLY zero torque the spin-up formula divides by T_d = 0, an undefined value, and the whole evaluation refuses.

Source: Hibbeler, R. C., Engineering Mechanics: Dynamics, 14th ed., Pearson, 2016 — §17.1 (mass moments of inertia; the uniform solid disk I = mR²/2), ch. 18 (planar kinetics, work and energy: E = ½Iω² for rotation about a fixed axis), and ch. 19 (planar kinetics, impulse and momentum: the principle of angular impulse and momentum, I_O·ω₁ + Σ∫M_O dt = I_O·ω₂, §19.1–19.2) — the basis for the constant-torque spin-up time T_d·t = I_z·ω.

Derivation

Steps marked modeling step are where physics enters by citation — every other line is machine-proven to follow from them, and the cited models are independently re-derived in the test pipeline where possible. See what is and isn't machine-verified.

σr=ω2ρ(R2r2)(ν+3)8\sigma_{r} = \frac{\omega^{2} \rho \left(R^{2} - r^{2}\right) \left(\nu + 3\right)}{8}

1. The modeling step: every annulus of a spinning disk needs a net inward pull of ρω²r per unit volume to stay on its circle — the disk is its own load. Radial equilibrium plus strain compatibility (plane stress, thin uniform disk) is a solvable ODE in σ_r; the two integration constants are fixed by σ_r staying finite at the centre and the rim being free, σ_r(R) = 0. This is the classical rotating-disk field. (The build's test suite re-substitutes it into the equilibrium ODE with a computer algebra system as an independent check.) — elasticity: radial equilibrium + compatibility, plane stress modeling step

σt=ω2ρ(R2(ν+3)r2(3ν+1))8\sigma_{t} = \frac{\omega^{2} \rho \left(R^{2} \left(\nu + 3\right) - r^{2} \left(3 \nu + 1\right)\right)}{8}

2. The same solution's hoop component. Both stresses are tensile everywhere — spinning stretches the disk in every in-plane direction at once, which is why a flywheel cracks radially when it lets go. — elasticity: hoop component of the same field modeling step

σr+σt=ω2r2ρ(1ν)4- \sigma_{r} + \sigma_{t} = \frac{\omega^{2} r^{2} \rho \left(1 - \nu\right)}{4}

3. Subtract: the hoop stress leads the radial stress by (1−ν)ρω²r²/4 ≥ 0 everywhere, so the hoop direction always governs, and the gap closes to zero at the centre where the two components meet. — algebra on the field

σmax=R2ω2ρ(ν+3)8\sigma_{max} = \frac{R^{2} \omega^{2} \rho \left(\nu + 3\right)}{8}

4. Evaluate either component at r = 0: the peak sits at the centre, where σ_r = σ_θ = (3+ν)ρω²R²/8. Read it as ρ(ωR)² times a number near ½ — the stress is set by the rim speed ωR, squared. Note what moved compared to the torsion shaft: there the stress ignored the material; here density IS the load, so the stress moves when the material does. — evaluate the field at the centre

Iz=R2mdisk2I_{z} = \frac{R^{2} m_{disk}}{2}

5. The energy side needs the spin inertia: I = ∫r² dm over the uniform disk gives ∫₀ᴿ r²·ρt·2πr dr = ρtπR⁴/2 = mR²/2. The r² weighting means the outer half of the radius does almost all the storing — which is why real flywheels are rim-heavy. (The test suite re-runs this integral too.) — mass moment of inertia by integration

Ek=R2mdiskω24E_{k} = \frac{R^{2} m_{disk} \omega^{2}}{4}

6. E = ½I_zω² in rim-speed form: a uniform disk stores ¼·m·(ωR)² — exactly half of what the same mass would carry moving in a straight line at rim speed. Energy, like stress, is a rim-speed-squared quantity. The collision is coming. — rotational kinetic energy

em=2σmaxρ(ν+3)e_{m} = \frac{2 \sigma_{max}}{\rho \left(\nu + 3\right)}

7. The collision: eliminate ω between energy and stress. Radius and thickness cancel completely — the energy per kilogram is locked to the peak stress by e = 2σ/((3+ν)ρ). No geometry can beat this. Run the disk up to its strength limit σ_y and every flywheel of a given material, any size, stores the same e_max = 2σ_y/((3+ν)ρ): energy per mass is a pure material index, σ_y/ρ. This is why grid-storage flywheels are carbon-fibre, not steel. — eliminate ω — the material merit index

Rωy=22σyρν+3R \omega_{y} = \frac{2 \sqrt{2} \sqrt{\sigma_{y}}}{\sqrt{\rho} \sqrt{\nu + 3}}

8. The same elimination read as a speed limit: yield onset happens at a fixed RIM SPEED, ω_y R = √(8σ_y/((3+ν)ρ)) — a material constant with units of m/s. A big slow disk and a small fast one hit the same wall. Spin speed in rpm is negotiable; rim speed is not. — yield-onset rim speed

tspin=IzωTdt_{spin} = \frac{I_{z} \omega}{T_{d}}

9. How long to get there. A constant drive torque T_d delivers angular impulse T_d·t; angular impulse–momentum sets that equal to the change in angular momentum I_z·ω when the disk is spun up from rest, so t = I_z·ω / T_d. More inertia or a higher target speed costs time; more torque buys it back. This is the lever a drivetrain pulls: the gear ratio ahead of the flywheel sets T_d, so cranking the ratio up spins the disk up faster — while winding the shaft that delivers T_d correspondingly tighter. Storing energy quickly is not free, and the chain-builder's headline example lets you feel that trade. — angular impulse–momentum, constant torque from rest

How it fails

The widget’s margin is first yield at the disk centre, under a steady spin, for a solid uniform disk. Real flywheels die harder and uglier:

  • Torsional Oscillator (Disk on a Shaft)

    A disk on an elastic shaft is a torsion pendulum: twist it and let go and it rings at one natural frequency, ω_n = √(k_t/J_d). The pitch is set entirely by the shaft's stiffness and the disk's inertia — not by how hard you twist it — while the stress it survives is set by the amplitude.

    • dynamics
    • stress
    • torque-power
  • DC Motor (Permanent Magnet)

    The machine that turns current into torque — and its whole personality is one straight line. At fixed voltage a PM DC motor trades speed for torque along T = T_stall(1 − ω/ω₀): two datasheet numbers pin every operating point, and the peak power hides at half the no-load speed.

    • torque-power
    • kinematics
  • Four-Bar Linkage (Position)

    Four pinned links — ground, crank, coupler, rocker — and the oldest mechanism in the book. Spin the crank and the rocker answers through pure geometry. Every position has TWO valid assemblies (open and crossed): the first THING in the catalog where one input has two honest answers, and the widget lets you pick the circuit.

    • kinematics
  • Slider-Crank (Exact Kinematics and Gas Torque)

    Crank, connecting rod, piston — the four-bar linkage with one pivot pushed to infinity, and the heart of every reciprocating engine and pump. Spin the crank at a fixed speed and the piston's position, velocity, and acceleration follow by pure geometry and two derivatives; push on the piston with a gas force and the connecting-rod obliquity turns it into a crank torque that swings from zero to a peak and back every revolution. The classic two-term r/l approximation rides alongside the exact form so you can watch it drift.

    • kinematics
    • torque-power
  • Rotating Disk with a Central Bore

    Drill the smallest possible shaft hole through a spinning disk and the peak stress exactly doubles — not "roughly increases": doubles, in the limit of a vanishing bore. The solid flywheel's optimistic numbers meet the hole every real rotor needs.

    • stress
    • energy-storage
    • mass-cost
  • Cantilever Beam (End Load)

    A beam fixed at one end, loaded at the other — the fruit-fly of structures. One widget shows why stiffness (E) and strength (σ_y) are independent axes: swap steel for titanium and deflection goes UP while the safety factor also goes up.

    • stress
    • mass-cost

Chains with

Outputs whose SI dimension and quantity kind match another THING's input — the only wires the planner's connectionLegal accepts (invariant 2, computed at build time, not hand-listed). Wire these on the chaining demo.

+ 11 more THINGs its outputs can legally feed (showing the first 8 in course order).

Sources