Damped pendulum Created 2026-10-06 Updated 2026-10-14
A damped pendulum is a simple pendulum whose mechanical energy is removed by a resisting torque. Here the model is undriven, with a fixed pivot, a point mass at distance , uniform gravitational acceleration , and resisting torque with constant . The angle is measured from the downward vertical in radians. An ideal viscous rotary dashpot realizes this torque law; tangential linear drag on the bob gives the same equation with . The moment of inertia is , so
Here is physical time, is the undamped small-amplitude angular frequency, and has units of inverse time. In the normalized equation , dots mean differentiation with respect to , not . The parameter is dimensionless. This separates the chosen damping law from the small-angle approximation: retaining allows finite angles even when the resisting torque is linear in angular velocity.
Figure 1.
Proposed viscous-damper experiment and illustrative small-angle decay: angle versus time and successive peak ratios, not measured data
. Details.
To motivate the model, release a pendulum without a push and record how its swing shrinks. A published example by Gintautas and Hubler uses a rotary sensor built from a computer mouse to measure the frequency and damping coefficient of a physical pendulum. [ref] The following is a proposed experiment with a chosen viscous rotary dashpot; the plotted observations are predictions, not a reproduction of their measurements.
Mount a compact bob on a light rigid rod with a fixed pivot, couple the pivot to an adjustable viscous rotary dashpot, and measure from the pivot to the bob's centre. First record decay with the dashpot disconnected, then repeat with it connected at several settings; this establishes how much decay is already due to other losses. Keep the bob, rod, camera, and release angle unchanged between these trials. The simple model requires rod and damper-rotor moments of inertia to be negligible compared with . Otherwise measure the total moment of inertia and the distance to the assembly's centre of mass: the physical pendulum equation becomes , where is the total moving mass. A damper that also exerts an elastic restoring torque would require another term.
Place a camera perpendicular to the swing plane, include the pivot and bob centre in every frame, and determine their coordinates relative to the resting vertical. For horizontal and downward separations , use ; record the frame timestamps rather than assuming a frame rate. Release from rest at about , record through several complete swings, and repeat the release. For every same-side turning point, tabulate its time and angular amplitude . Measure the oscillation period from differences ; comparing several periods exposes drift. Interpolating around each observed extremum reduces frame discretization error. Check that late peaks still exceed the angular resolution; taking the natural logarithm of tiny, uncertain peaks amplifies their relative errors.
In the small-angle approximation, the damped pendulum becomes a damped harmonic oscillator. For and release from rest at , its exact linear-model solution is
The decaying sinusoidal envelope has coefficient ; because release is from rest, its successive same-side extrema occur at and have . Thus the envelope and the sequence of actual peak heights are different descriptions of the same trace. The linear damping regimes and this solution are derived in Richard Fitzpatrick, Damped Oscillatory Motion, equations 82–93.
The experiment tests two predictions independently: nearly equal same-side peak spacing and a straight line when plotting against . Its predicted slope is , giving
Use several peaks and repeated trials rather than a single ratio, and compare the frequency identity with independently measured and . For the illustrative choices , , , and , the linear model gives , , and . The first two positive peaks after release are therefore approximately and . The undamped small-amplitude oscillation period is about : substantial loss per swing need not imply a large change in period.
Repeat at a larger release angle to distinguish two limitations. The relative error in replacing by is about for small angles: about at , but about at when compared with . Finite-amplitude swing spacing can therefore change as the swing shrinks even if is constant. Separately, friction at the pivot or quadratic drag can invalidate the assumed linear resisting torque; exponential peak decay should be checked, not imposed on all recorded motion. Peters discusses why real pendulum damping can require nonlinear terms. [ref] Repeatability, residual pivot motion, changes in damper temperature, and finite camera resolution are useful checks on this proposed measurement.
The normalized equilibrium points are . Even multiples of are attracting: their linearization has eigenvalues , giving a stable focus for , repeated negative eigenvalues at , and a stable node for . These are the underdamped, critically damped, and overdamped small-displacement regimes. Odd multiples of have eigenvalues of opposite signs and are saddle equilibria. The energy dissipation of a damped pendulum explains why trajectories in the phase portrait spiral inward near a stable focus instead of closing into the recurring orbit of an undamped simple pendulum.
Moment of inertia Created 2026-09-24 Updated 2026-10-14
The moment of inertia about a specified axis measures how a mass distribution resists changes of angular velocity about that axis. For point masses at perpendicular distances , or a continuous distribution with mass density ,
Its SI unit is . The distance is to the entire axis, not to one point on it: moving a point mass along the axis leaves its contribution unchanged. The moment of inertia depends on both the distribution and the chosen axis; it is not determined by total mass alone. About an axis fixed in an inertial frame, a rigid body with larger moment of inertia has smaller angular acceleration under the same net axial torque, and greater rotational kinetic energy at the same angular speed.[ref]
For a simple model, put four equal point masses of total mass on light rigid spokes at radius . Their moment of inertia is . Doubling quadruples the moment of inertia without changing . The illustrated placements and give and . Actual spokes and finite-sized weights contribute their own moments of inertia; the point mass model neglects them. Repositioning weights during rotation also breaks the fixed-distribution assumption, so it requires a separate work and angular momentum calculation.
Figure 1.
Same total mass, different radial placement and moment of inertia
. Original model diagrams, not experimental data. The upper pair has four equal point masses on negligible-mass spokes, at and from a perpendicular central axis. The lower pair compares a uniform disc with an ideal thin ring of the same mass and outer radius . The displayed moments of inertia determine rotational kinetic energy at a shared angular speed. The visible thickness of the ring and size of the point markers are only drawing conventions. Details.
The moment of inertia of a uniform disc follows by dividing it into concentric annuli: and . For an ideal thin ring every mass element is at , so . A ring of finite radial thickness instead uses the moment of inertia of a uniform annulus, , with inner radius . A moment of inertia is thus a compact summary of radial placement rather than a complete description of shape. For several joined pieces, their moments of inertia add only when each is calculated about the same axis.[ref]
A physical pendulum makes the moment of inertia measurable through timing. Martín Monteiro, Cecilia Stari, Cecilia Cabeza and Arturo C. Martí used a metal bar with suspension holes apart and an attached smartphone. They report bar length , bar mass and phone mass . Changing the suspension hole changes the pivot axis. Their readouts were phone gyroscope angular velocity records, analyzed by autocorrelation or a fit to a sinusoid, and an independent photogate; the logged records for the sinusoid fits lasted . They observed agreement between the three timing methods and the geometry-based model, including a minimum of oscillation period as the suspension point varied. These are reported apparatus details and observations, not measurements made for this article.[ref]
To infer a moment of inertia from such a measurement, determine the total mass including every attached part, locate the assembly's centre of mass , and measure its distance below the horizontal pivot . Record successive crossings in the same direction, or the elapsed time for complete oscillations, giving oscillation period . The physical pendulum model with a fixed pivot, a uniform gravitational field, negligible friction and the small-angle approximation gives
The crossing times, masses and geometric distances are the measured inputs; is an inference conditional on this model, while an integral over the assumed shape supplies an independent model prediction. The parallel axis theorem converts the inferred pivot moment of inertia to one about a parallel axis through . The controls and interpretation of that comparison are described there. An axial measurement determines one scalar moment of inertia, not the complete inertia tensor needed for general rotation.
Parallel axis theorem Created 2026-09-24 Updated 2026-10-14
The parallel axis theorem transfers a moment of inertia from an axis through the centre of mass to a parallel axis through :
Here is the total mass and is the shortest perpendicular distance between the two axes. Both moments of inertia concern the same mass distribution and the same axis direction. The axis for must pass through the centre of mass; may refer to an axis outside the body. No uniform-density or planar-shape assumption is needed. Within any family of parallel axes, the centre of mass axis minimizes the moment of inertia.[ref]
A short proof isolates why the centre of mass matters. Choose parallel axes perpendicular to an plane, put at the origin, and place at . Then
because in centre of mass coordinates. Individual mass elements can move closer to the new axis, but their weighted cross terms cancel. For arbitrary parallel axes and , use , where each distance is measured from the parallel centre of mass axis; adding times the squared separation of and is generally incorrect. A displacement along the axis has . Changing the axis direction instead requires the inertia tensor or, in its applicable planar setting, the perpendicular axis theorem.
For a uniform disc, mass and radius , the central perpendicular moment of inertia is . Shifting that axis to the rim gives and , three times . In the computed , example, , and . At angular speed the rim-pivoted disc has rotational kinetic energy . Its centre of mass now moves, so already includes that translation; adding a separate centre of mass translation term to it would count that contribution twice.
Figure 1.
Parallel axis theorem: a uniform disc's central and rim axes
. Original geometric and numerical model, not measured data. Both axes are perpendicular to the disc; the separation is . The bar chart separates from the additional in , normalized by . All parts of the disc remain in the calculation, including those nearer to the shifted axis. Details.
The physical pendulum measurement described under moment of inertia provides a practical check: infer at several pivot positions, then subtract to obtain . Keep all pivot axes parallel and measure to the centre of mass of the entire moving assembly. The same assembly should yield the same within measurement precision. If a smartphone is attached, its mass and its own moment of inertia belong in both the geometric prediction and the centre of mass calculation. The published bar-and-phone study explicitly adds both contributions using the parallel axis theorem.[ref]
For a careful replication, release small swings without a push, repeat timings at several small amplitudes, and check that the extrapolated oscillation period is stable as amplitude decreases. Time many complete oscillations and repeat the run; two crossings of a photogate per cycle must not be mistaken for two cycles. Inspect decay of amplitude and compare early and late timings to detect significant pivot friction or air drag. Mark the bar and phone mount to detect relative slip, tighten the attachment, and check that the pivot support does not rock or migrate between holes. Repeat with an independent photogate or video readout to check the phone's timing. These are proposed controls for the replication, not additional results reported by the authors. The ideal model excludes appreciable damping, large swings and flexing; subtraction of also magnifies relative errors when is much smaller than .
The oscillation period need not increase monotonically with pivot distance: the shift raises the moment of inertia but also increases the gravitational restoring torque. For a computed slender uniform bar model of length , and
This expression has its minimum at . With , and stipulated , the prediction is ; it is not a recorded timing or a prediction for the bar with a phone attached. The competing terms explain why a timing experiment can reveal a mass distribution rather than merely its weight.
The moment of inertia of the uniform disc about its central perpendicular axis is
The parallel axis theorem gives . For angle from the downward resting position, the gravitational torque is , so . Linearizing at gives . Thus the uniform disc pivoted at its rim is a physical pendulum with period
The uniform rod has linear mass density . Integrating squared distance from the pivot gives its moment of inertia about the axis perpendicular to the plane of motion:
Measure from the downward vertical. The center of mass is at distance from the pivot, so gravity produces the restoring torque . The fixed-axis rotational equation of motion is . Therefore the physical pendulum obeys
The pivot force has zero torque about the pivot.
Both instruments regulate a sequence of time intervals with an oscillator. A restoring action opposes displacement, while inertia carries the moving part through equilibrium. Their repeated swings control an escapement, which releases the gear train in steps and supplies enough energy to offset frictional losses.
For a simple pendulum of length ,
At small amplitude, , giving simple harmonic motion with
A real clock pendulum is a physical pendulum. If its pivot moment of inertia is and its centre of gravity is a distance below the pivot, its small-swing period is . Gravity supplies the restoring torque. A controlled swing amplitude and stable dimensions are needed: large swings, thermal expansion and changing effective gravity alter the rate.
A watch's balance wheel turns back and forth about its spindle. A torsion spring supplies an approximately linear restoring torque , so
The restoring mechanism is elastic rather than gravitational. This makes the balance-spring system suitable for a portable watch, although disturbances, imperfect balance, changing spring properties and escapement interactions still affect the rate.
Figure 1.
Gravitational and elastic restoring torques in clock oscillators
. Original schematic comparing a pendulum with a balance wheel and spiral spring. Equations show ideal small-amplitude models, not the detailed construction of a particular historical clock. Details.
Near independence of period from amplitude is isochronism. It is approximate for the simple pendulum: its leading correction is with in radians. Keeping the oscillations repeatable therefore matters as well as selecting the nominal period.
Let place the centre of mass vertically below a horizontal pivot axis. The restoring torque and rotational kinetic energy are and , so
For small angles, , and the equation reduces to simple harmonic motion. Thus the physical pendulum has
This assumes negligible friction and a horizontal axis, with . For a tilted axis the projected gravitational restoring torque changes; an axis through the centre of mass has no gravitational restoring torque.
For the general Lagrangian just obtained, the Euler-Lagrange equation gives
Using the small-angle approximation , this is simple harmonic motion with . Thus
The endpoint mass cancels in the light-rod limit. A concrete comparison shows why the finite rod cannot silently be included: the same specified endpoint with gives , distinct from the negligible-rod result. This is a physical pendulum when the rod contributes, and a simple pendulum in the light-rod limit. Both expressions require small oscillation amplitude; the full sine equation remains nonlinear.
The parallel axis theorem gives for parallel axes separated by , with the reference axis through the centre of mass. For a plane lamina, the perpendicular axis theorem gives for mutually perpendicular axes meeting at one point, with in the plane.
Treat the removed material as negative mass. With the hole's centre at , its moment of inertia about is . The remaining moment of inertia and centre of mass are
For a uniform original disk, . Write and measure from the downward equilibrium direction of its centre of mass. The physical pendulum then has Lagrangian
The Euler-Lagrange equation is . For small oscillations, , giving
Now add the removed small disk on the opposite side, with its centre a distance from so that it is internally tangent to the large disk. The total mass returns to , while its first moment about is . Thus
Replacing both the inertia and gravitational restoring moment gives
The same first-moment method applies to a general perforated-disc pendulum. The added disk is attached to the lamina and shares its rotation; treating it as a massless patch or updating only the inertia would omit its restoring torque.
Figure 1.
Removing an off-centre disk and reattaching it on the opposite side changes both the pendulum inertia and its centre of mass
. Details.
Take the positive coordinate from toward . The centre of mass of the remaining material has coordinate
Thus . If is the small angular displacement from the stable hanging orientation, gravity gives torque about the suspension point. The physical pendulum equation is
Using the small-angle approximation gives simple harmonic motion with
Both the reduced mass and displaced centre of mass matter; changing only the moment of inertia would give an incorrect frequency.
Inertia is resistance to a change of motion. Translational inertia is measured by mass in Newton's second law; rotational inertia about a specified fixed axis is measured by the moment of inertia
where is perpendicular distance from that axis. For a thin uniform rod, take axes perpendicular to the rod. With linear mass density ,
The second result uses the parallel axis theorem; direct integration over gives the same value. The axis specification is essential: an ideal thin rod has zero moment of inertia about its own longitudinal axis.
During the short elastic collision, neglect the impulse from gravity. The pivot can exert an external impulse, so total linear momentum of rod plus particle need not be conserved. Its torque about the pivot is zero, so conservation of angular momentum about that point applies. Let be the particle's signed outgoing horizontal velocity, and the rod's outgoing angular velocity, positive towards the incident motion. Then
Substituting into the kinetic energy equation gives
Discarding the zero-impulse branch , the actual collision gives
The rod's lower end initially moves forward at , so the bodies separate. Afterwards the free particle undergoes projectile motion. With the pivot as origin and downward, its trajectory is , , where at the collision. The rod instead swings about the fixed pivot in a vertical plane as a physical pendulum.
The rod's initial rotational kinetic energy is . Reaching the horizontal raises its centre of mass by , increasing gravitational potential energy by . Consequently,
For the threshold case, measure the rod's angle from the downward vertical and reset at its first horizontal turning point, so and . The restoring torque about the pivot gives the exact equation of motion
At the horizontal, , so the rod immediately swings back. Conservation of mechanical energy yields
On the initial descent its explicit time dependence is specified by the definite integral
It continues through the vertical to , reverses, and repeats without damping. An alternative expression valid throughout the oscillation uses the Jacobi elliptic sine. Put , , and let denote the complete elliptic integral of the first kind. Then
This is not simple harmonic motion at an amplitude of . Its oscillation period exceeds its own small-amplitude value by the factor .
A simple pendulum of length obeys . At the same amplitude , its oscillation period is , so
The rod swings more slowly. Its equivalent simple pendulum length is , not the distance to its centre of mass. The same oscillation period ratio holds in the common small-amplitude limit, but replacing the actual 90-degree motion by that limit would underestimate both oscillation periods.
Axis theorems. The parallel axis theorem states
where is the moment of inertia about an axis through the centre of mass, is about a parallel axis, is the body's total mass, and is the perpendicular separation of the axes. Expanding proves the result because at the centre of mass. The perpendicular axis theorem applies to a planar lamina in the plane:
All three axes meet at the same point, lie in the lamina and is normal to it; the identity follows from being the squared distance to the axis. It is not a theorem for an arbitrary three-dimensional body.
Cylinder integrals. Put the cylinder's midpoint at and its symmetry axis along . Its uniform mass density is . The transverse moment of inertia of a uniform solid cylinder about the axis is
Equivalently each thin circular slice has diametral moment of inertia , by the perpendicular axis theorem and symmetry, and shifting it by adds . The parallel axis theorem then gives the transverse end-axis moment of inertia:
The axial moment of inertia of a uniform solid cylinder is . Equating it to gives
Here the source's principal axis means the cylinder's symmetry axis; its transverse centre axes are also principal axes. This equality describes a short cylinder rather than the thin-rod limit.
Hammer dynamics and forces. Treat the compact head as a point mass and neglect its intrinsic moment of inertia. Let be the displacement from downward vertical, and let point along the rod from pivot to head, with in the direction of increasing . Define
The two external weights and act at distances and from the pivot. The frictionless pivot supplies a reaction but no torque about its axis. The total restoring torque and rotational equation of motion are
For completeness, Newton's second law for the centre of mass gives both components of the pivot force:
The negative radial component points towards the pivot. These follow from the centre of mass acceleration and the total weight . They specify the external forces even when the oscillation is not small.
With the small-angle approximation, , this is simple harmonic motion. The physical pendulum has small-amplitude oscillation period
For release from a finite amplitude , conservation of mechanical energy instead gives . Integrating a quarter-cycle gives the exact oscillation period
where is the complete elliptic integral of the first kind; recovers .
For and , compare with the simple pendulum value :
For small oscillations, the hammer's oscillation period is about shorter. The extra rod weight adds restoring torque proportionately faster than it adds pivot moment of inertia.
In the small-angle approximation the physical pendulum is a harmonic oscillator with . Its classical energy is . Equating this to gives
An oscillating bob is not a single stationary energy eigenstate. An ideal coherent state has a Poisson distribution of oscillator occupations centred on this mean, with relative spread , exceptionally small here. This quantum occupation scale of a classical oscillator fixes the excitation scale without requiring the real, damped clock to be a pure coherent state.
The moment of inertia about an axis is . The parallel axis theorem states , where a is the perpendicular separation from a parallel axis through the centre of mass, and M is total mass.
Since , the density coefficient has units . The density-power cylindrical moment of inertia follows from cylindrical shells:
For an offset pivot, the physical pendulum equation is . At small amplitude,
For the minimum-period offset of a physical pendulum, minimizing is equivalent to minimizing . Its derivative vanishes at
The second derivative is positive. The period diverges as near zero and grows as for large offset, giving one minimum. This assumes the stated ideal pivot can pass inside the body; an externally constrained pivot would restrict the allowable a. [ref]
Keeping one resonant harmonic and treating orbital eccentricity as fixed reduces mean-motion resonance dynamics to a physical pendulum equation. For the inner first-order term , leading Lagrange planetary equations give . Around its stable phase, the small-amplitude libration frequency is . Evolution of eccentricity, other harmonics, and approach to a separatrix change this approximation.
A uniform rigid rod of length and mass , pivoted at one end with a point mass at the other, has pivot moment of inertia . Its gravitational potential energy is . The planar Euler-Lagrange equation is , giving the displayed small-amplitude oscillation period. The limit gives a simple pendulum, whereas gives the uniform-rod physical pendulum. Thinness does not specify a rod's mass, so both contributions must be stated or explicitly neglected.
Stokes flow Created 2026-09-24 Updated 2026-10-14
Stokes flow is the limit of fluid mechanics in which inertia is negligible: pressure gradients, viscosity and any body forces balance at each instant. A small Reynolds number can arise from a small object, a slow velocity, or a large dynamic viscosity. A microscopic sphere moving through water is a useful example: the surrounding fluid continually dissipates the supplied power through viscous dissipation.
Figure 1.
Exact Stokes streamlines around a stationary no-slip sphere in a uniform incident flow
. Details.
The picture is a meridional slice of the three-dimensional axisymmetric flow around a sphere, with the sphere held fixed and the distant fluid moving rightward at speed . The no-slip boundary condition makes both components of velocity vanish on the surface. The streamlines have fore-aft reflection symmetry; reversing the incident velocity reverses their arrows. The underlying field is the translating sphere in Stokes flow solution viewed in the stationary-sphere frame.
For an incompressible flow of a Newtonian fluid with constant dynamic viscosity , the steady Stokes equation is
Here is pressure, is velocity, and is body force per unit volume. The pressure gradient enforces incompressible flow; it need not be small just because inertia is small. At speed and length , the motion-induced pressure gradient has typical scale , corresponding to a disturbance pressure difference of order . See Shreyas Mandre's fluid dynamics notes, sections 9.1.1–9.1.2.
The controlling Reynolds number is
where is density and is kinematic viscosity. The condition suppresses advection of momentum on the length . Rapid changes require an additional check of the viscous diffusion time: small Reynolds number alone does not eliminate the time derivative in the Unsteady Stokes equation.
Two rewards of this limit are the exact Stokes's law for the drag on a sphere and kinematic reversibility of Stokes flow. The first predicts settling velocities and connects to the Stokes–Einstein relation; the second explains why a microswimmer needs a nonreciprocal stroke, as expressed by the scallop theorem. Despite viscous dissipation, reversing the driving can retrace the motion: an especially striking consequence of Linearity of Stokes flow.
George Gabriel Stokes obtained the sphere drag and its terminal velocity in his 1851 investigation of how viscosity changes physical pendulum motion. He also used the slow fall of small water droplets to discuss clouds. The paper was read on 9 December 1850; 1851 is its publication date. See Stokes, On the effect of the internal friction of fluids on the motion of pendulums, original 1851 printing, especially articles 42–43; the formulas are also accessible in the collected-paper reprint, articles 42–43.