Torsional Oscillator (Disk on a Shaft)

dynamicsstresstorque-power

Verified build 11 relations · 5 identities proven · 1 modeling step · 3 parity samples

Clamp one end of a shaft, hang a disk on the other, give it a twist, and let go: it rings. The disk winds the shaft up, the shaft winds it back, and the two trade energy back and forth at one fixed rate — the system’s natural frequency. This is the torsional cousin of a mass bobbing on a spring, and it is the single most important idea in machine dynamics, because every real driveline, crankshaft, and torsion-bar is one of these hiding in plain sight.

ωn=ktJd,kt=GJpL,f=ωn2π,T=1f\omega_n = \sqrt{\frac{k_t}{J_d}}, \qquad k_t = \frac{G J_p}{L}, \qquad f = \frac{\omega_n}{2\pi}, \qquad T = \frac{1}{f}

The shaft is a torsion spring of stiffness kt=GJp/Lk_t = G J_p/L — literally the inverse of the twist the shaft-in-torsion page computes, θ=TL/GJp\theta = TL/GJ_p. The disk is the inertia Jd=12mdR2J_d = \tfrac12 m_d R^2 — literally what the flywheel page computes. Bolt the two together and you get an oscillator. Three things to notice:

One material, one simplification. Here a single material sets both the shaft and the disk, so you cannot make only the disk denser in the widget — increasing density raises JdJ_d (lowering the pitch) and changes GG at the same time. The counter-intuitive “a heavier disk rings lower” story still lives in the readouts: watch JdJ_d climb with density. Independent disk-and-shaft materials are a Phase 3 capability (multi-material binding slots); this page states the shared-material assumption up front rather than fake it.

The chain teaser (coming in Phase 4). This page is built to be chained. Its two governing inputs are exactly two other pages’ outputs: the shaft stiffness that torsion-shaft embodies and the disk inertia that flywheel-disk produces. The peak restoring torque TdynT_{dyn} it reports back is an honest torque you could feed into the shaft page to size it. Wire motor → shaft → flywheel and you have simulated a real driveline’s natural frequency, every number carrying its citation. The chain-builder that lets you draw those wires arrives in Phase 4; for now the ports are here, waiting.

A note on units: the natural frequency ff is measured in hertz and the spin-like ωn\omega_n in radians per second, and although both are “per second” the site treats them as different quantity kinds — you cannot silently wire a 35 Hz ring into a port expecting a 35 rad/s spin. The factor of 2π2\pi between them is always explicit. That distinction is exactly what the dimensional type system exists to protect.

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)
G4 Msitypicalmil-hdbk-5j
sigma_y47 ksidesign min.mil-hdbk-5j
rho0.1 lb/inch**3typicalmil-hdbk-5j
Inputs
N·m
Polar second moment of area (shaft)
m⁴
Disk mass
kg
Disk mass moment of inertia (spin axis)
kg·m²
Shaft mass moment of inertia
kg·m²
Natural angular frequency
Natural frequency
Hz
Period
Peak restoring torque (at amplitude)
N·m
Peak surface shear stress (at amplitude)
Safety factor (shear yield at amplitude)
Static twist under applied torque

4 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) Ti-6Al-4V

Governing relations

Jp=πd432J_p = \frac{\pi d^4}{32}

Assumes: solid circular shaft cross-section

Source: Gere, J. M., & Goodno, B. J., Mechanics of Materials, 9th ed., Cengage, 2018 — §3.3–3.5 (circular shafts in torsion: the torsion formula τ = 16T/πd³, the angle of twist θ = TL/GJ, and hence the torsional stiffness k_t = GJ/L).

md=ρπR2tdm_d = \rho\,\pi R^2 t_d

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 = ½ m R² about its spin axis) and ch. 22 (undamped free vibration: a single mass on a spring oscillates at ω_n = √(k/m), the translational analogue of the torsional system here).

Jd=12mdR2J_d = \tfrac{1}{2} m_d R^2

Assumes: rigid uniform disk, polar mass moment about its own spin axis (the shaft axis)

Source: Hibbeler, R. C., Engineering Mechanics: Dynamics, 14th ed., Pearson, 2016 — §17.1 (mass moments of inertia; the uniform solid disk I = ½ m R² about its spin axis) and ch. 22 (undamped free vibration: a single mass on a spring oscillates at ω_n = √(k/m), the translational analogue of the torsional system here).

Js=ρπd432LJ_s = \rho\,\frac{\pi d^4}{32}\,L

Assumes: polar mass moment of the solid shaft about its own axis (ρ·J_p·L), used only to check the lumped-parameter assumption — it is NOT added to J_d in this model

Source: Hibbeler, R. C., Engineering Mechanics: Dynamics, 14th ed., Pearson, 2016 — §17.1 (mass moments of inertia; the uniform solid disk I = ½ m R² about its spin axis) and ch. 22 (undamped free vibration: a single mass on a spring oscillates at ω_n = √(k/m), the translational analogue of the torsional system here).

ωn2=ktJd=GJpLJd\omega_n^2 = \frac{k_t}{J_d} = \frac{G J_p}{L\,J_d}

Assumes: simple harmonic motion of a single disk on a massless torsion spring — the lumped one-degree-of-freedom model (physics enters here by citation; no ODE is integrated); small amplitude, so the shaft stays linear-elastic and k_t = G J_p/L is constant · Valid while: The shaft's own rotational inertia is no longer negligible next to the disk's (J_shaft/J_d exceeds 0.1). The lumped model treats the shaft as a massless torsion spring; a heavy shaft stores kinetic energy too, which lowers the true frequency (a first correction adds about one-third of J_shaft to J_d), so ω_n here reads high.

Source: Timoshenko, S. P., Young, D. H., & Weaver, W., Vibration Problems in Engineering, 4th ed., Wiley, 1974 — the torsional pendulum / disk-on-shaft system: a rotor of mass moment J on a shaft of torsional stiffness k_t oscillates in simple harmonic motion at ω_n = √(k_t/J).

f=ωn2πf = \frac{\omega_n}{2\pi}

Assumes: frequency in cycles per second; ω_n is the same motion measured in radians per second

Source: Budynas, R. G., & Nisbett, J. K., Shigley's Mechanical Engineering Design, 10th ed., McGraw-Hill, 2015 — ch. 7 (shafts and torsional systems), §3-12 (power transmission), and §5-4 (maximum-shear-stress theory, shear yield at S_sy = 0.5 S_y).

T=1f=2πJdktT = \frac{1}{f} = 2\pi\sqrt{\frac{J_d}{k_t}}

Assumes: the time for one full oscillation

Source: Budynas, R. G., & Nisbett, J. K., Shigley's Mechanical Engineering Design, 10th ed., McGraw-Hill, 2015 — ch. 7 (shafts and torsional systems), §3-12 (power transmission), and §5-4 (maximum-shear-stress theory, shear yield at S_sy = 0.5 S_y).

Tdyn=ktΘ=GJpLΘT_{dyn} = k_t\,\Theta = \frac{G J_p}{L}\,\Theta

Assumes: peak torque the shaft carries at the extremes of the swing, where twist equals the amplitude Θ; linear torsion spring

Source: Gere, J. M., & Goodno, B. J., Mechanics of Materials, 9th ed., Cengage, 2018 — §3.3–3.5 (circular shafts in torsion: the torsion formula τ = 16T/πd³, the angle of twist θ = TL/GJ, and hence the torsional stiffness k_t = GJ/L).

τmax=16Tdynπd3\tau_{max} = \frac{16\,T_{dyn}}{\pi d^3}

Assumes: linear-elastic torsion, shear stress maximal at the shaft surface (same field as a statically twisted shaft, evaluated at twist Θ); pure torsion — no bending, no stress concentrations (keyways, fillets, press-fits) · Valid while: Surface shear stress at the swing extremes exceeds the shear yield strength (σ_y/2 by the maximum-shear-stress criterion). The shaft yields during each cycle, the spring rate k_t stops being constant, and every elastic number here — starting with the frequency — drifts. This is also where torsional fatigue lives.

Source: Gere, J. M., & Goodno, B. J., Mechanics of Materials, 9th ed., Cengage, 2018 — §3.3–3.5 (circular shafts in torsion: the torsion formula τ = 16T/πd³, the angle of twist θ = TL/GJ, and hence the torsional stiffness k_t = GJ/L).

SF=σy/2τmax\mathrm{SF} = \frac{\sigma_y/2}{\tau_{max}}

Assumes: margin against shear yield at the oscillation amplitude (Tresca; conservative next to the distortion-energy value 0.577 σ_y)

Source: Budynas, R. G., & Nisbett, J. K., Shigley's Mechanical Engineering Design, 10th ed., McGraw-Hill, 2015 — ch. 7 (shafts and torsional systems), §3-12 (power transmission), and §5-4 (maximum-shear-stress theory, shear yield at S_sy = 0.5 S_y).

θst=TappLGJp\theta_{st} = \frac{T_{app}\,L}{G J_p}

Assumes: the same shaft acting as a static torsion spring — the twist a steady torque T_app would hold (the reciprocal of k_t, θ_st = T_app/k_t), a bridge to the shaft-in-torsion page · Valid while: The steady torque's surface shear stress (16 T_app/πd³) exceeds the shear yield strength (σ_y/2, maximum-shear-stress criterion): the shaft would yield before holding this twist, so the elastic θ_st shown here is fictional past this point — the real shaft takes a permanent set instead.

Source: Gere, J. M., & Goodno, B. J., Mechanics of Materials, 9th ed., Cengage, 2018 — §3.3–3.5 (circular shafts in torsion: the torsion formula τ = 16T/πd³, the angle of twist θ = TL/GJ, and hence the torsional stiffness k_t = GJ/L).

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.

kt=GJpLk_{t} = \frac{G J_{p}}{L}

1. Start with the shaft as a torsion spring. The shaft-in-torsion page gives the twist under a torque, θ = T L/(G J_p); invert it and the shaft's torsional stiffness — the torque it takes to wind it up one radian — is k_t = T/θ = G J_p/L. Fatter shaft (d⁴), stiffer spring; longer shaft, softer spring. — torsion: stiffness is the inverse of compliance

ωn2=ktJd\omega_{n}^{2} = \frac{k_{t}}{J_{d}}

2. The modeling step, where the physics enters by citation. Displace the disk by a small angle φ from rest; the shaft pulls it back with a torque −k_t φ, and Newton's law for rotation is J_d φ'' = −k_t φ. That is the equation of simple harmonic motion, whose solution φ(t) = Θ sin(ω_n t) oscillates at ω_n² = k_t/J_d. We do NOT integrate this ODE in the build — the SHM result is taken from the vibration literature and verified only as a consistent definition (invariant 5: SymPy checks algebra, not physics). — dynamics: single-DOF simple harmonic motion (cited) modeling step

ωn=GJpJdL\omega_{n} = \frac{\sqrt{G} \sqrt{J_{p}}}{\sqrt{J_{d}} \sqrt{L}}

3. Substitute k_t = G J_p/L. Notice what is NOT here: the amplitude Θ has cancelled. A big swing and a tiny swing ring at exactly the same pitch — the isochronism that made the pendulum clock possible. Stiffness raises the pitch, inertia lowers it; the drive that set the swing going does not. — substitute the shaft stiffness

f=ωn2πf = \frac{\omega_{n}}{2 \pi}

4. Cycles per second is radians per second over 2π. This 2π is the whole reason f and ω_n are different quantity kinds even though they share units of 1/s: on this site an f-port can never be wired into an ω-port without the explicit factor, so a 35 Hz ring is never mistaken for a 35 rad/s spin. The period is its reciprocal, T = 1/f = 2π√(J_d/k_t). — definition of frequency

τmax=GΘd2L\tau_{max} = \frac{G \Theta d}{2 L}

5. Now the amplitude comes back — in the stress. At the extremes of the swing the shaft is twisted by Θ, so its peak restoring torque is T_dyn = k_t Θ and the surface shear is τ_max = 16 T_dyn/(π d³), which collapses to G d Θ/(2L) — exactly the shaft-in-torsion stress evaluated at twist Θ. The frequency ignores the amplitude; the stress is proportional to it. Ring it twice as hard and the pitch holds while the stress doubles. — surface shear at the swing extreme

ωn2=Gd416LR4ρtd\omega_{n}^{2} = \frac{G d^{4}}{16 L R^{4} \rho t_{d}}

6. Expand J_p and J_d and the material collapses to a single group: ω_n² = (G/ρ)·d⁴/(16 L R⁴ t_d). Frequency depends on the material ONLY through G/ρ — the square of the shear-wave speed. That ratio barely moves across metals, so a steel disk-shaft and an aluminium one of the same geometry ring at nearly the same pitch; swap in a polymer (much lower G/ρ) and the pitch drops hard. Strength (σ_y) sets what amplitude it survives, not what note it plays. — eliminate geometry — the material index G/ρ

How it fails

The widget reports one clean natural frequency and the shear stress at a given amplitude. Real torsional systems fail in ways that clean single number quietly hides:

None of these are edge cases; they are the reason torsional-vibration analysis is its own discipline. The natural frequency is where that analysis starts, not where it ends.

  • 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
  • Flywheel (Solid Rotating Disk)

    The machine that stores work as spin — and loads itself doing it. Centrifugal self-loading grows with ρω²R², so the energy a flywheel can hold per kilogram is capped not by its size but by one material index: strength over density.

    • energy-storage
    • stress
    • mass-cost
  • 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
  • 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
  • Impact Loading (Falling Mass, Energy Method)

    Drop a mass onto an elastic member and the peak stress is not W/A — it is n times the static stress, where the impact factor n = 1 + √(1 + 2h/δ_st). A suddenly-applied load (h = 0) already doubles the stress; a real drop multiplies it many times over. Stiffer members take HIGHER impact stress, because a smaller static deflection means a larger n.

    • dynamics
    • stress
  • Shaft Critical Speed (Rayleigh + Dunkerley)

    Spin a shaft fast enough and it whips: at its critical speed the rotor whirls in resonance with its own static sag, ω_c = √(g/δ_st). Gravity sets the sag but cancels out of the answer — the critical speed is pure stiffness over inertia. Dunkerley's estimate folds in the shaft's own mass and is provably never higher than Rayleigh's.

    • dynamics
    • stress

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.

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

Sources