Physlib.ClassicalMechanics.Pendulum.SimplePendulum.Basic
The simple gravity pendulum
i. Overview
A simple gravity pendulum is a bob of mass `m` fixed to one end of a rigid massless rod of length `ℓ`, the other end of which is pinned at a pivot, swinging in a vertical plane under a uniform gravitational acceleration `g`. Its configuration is the angle `θ` of the rod from the downward vertical. The bob moves on the circle of radius `ℓ` about the pivot, so its moment of inertia about the pivot is `I = m ℓ²` and its kinetic energy is `T = ½ I θ̇²`. Swinging out to the angle `θ` raises the bob by `ℓ (1 - cos θ)` against gravity, so the potential energy is `V = m g ℓ (1 - cos θ)`, normalized to vanish at the bottom of the swing. Balancing the rate of change of the angular momentum about the pivot against the torque `-m g ℓ sin θ` of gravity gives the equation of motion `I θ̈ = -m g ℓ sin θ`, equivalently `θ̈ + (g/ℓ) sin θ = 0`. The mass drops out of the motion, which is governed by the single quantity `ω = √(g/ℓ)`, the angular frequency of the small oscillations about the bottom.
The configuration of the pendulum is genuinely an angle modulo a full turn, an element of the circle `SimplePendulum.ConfigurationSpace`. As for the harmonic oscillator, the dynamics in this file are written instead on the Euclidean lift `Time → EuclideanSpace ℝ (Fin 1)`: the angle is carried by a real number, from which the configuration is recovered by `SimplePendulum.ConfigurationSpace.ofAngle`, and the one-dimensional Euclidean space stands in for both the configuration space and its tangent space, so that the Euler–Lagrange operator of Physlib applies verbatim. Two lifts differing by `2π n` describe the same motion, as a subsequent contribution proves; the connection of the model here with the geometric configuration space is made in a later module.
ii. Key results
- `SimplePendulum` contains the input data of the problem: the mass `m` of the bob, the length `ℓ` of the rod and the gravitational acceleration `g`. - `SimplePendulum.ω` is the angular frequency `√(g/ℓ)` of the small oscillations, and `SimplePendulum.inertia` is the moment of inertia `m ℓ²` of the bob about the pivot. They are tied together by `SimplePendulum.ω_sq_mul_inertia`, the identity by which the mass cancels from the equation of motion. - `SimplePendulum.kineticEnergy`, `SimplePendulum.potentialEnergy` and `SimplePendulum.energy` are the energies, with the bounds `potentialEnergy_nonneg`, `potentialEnergy_le` and `potentialEnergy_eq_zero_iff`, the gradient `gradient_potentialEnergy` of the potential and the time derivatives `kineticEnergy_deriv`, `potentialEnergy_deriv` and `energy_deriv`. - `SimplePendulum.lagrangian` is the Lagrangian `T - V` of the pendulum, and `SimplePendulum.torque` is the torque about the pivot, the generalized force conjugate to the angle. - `SimplePendulum.EquationOfMotion` is the equation of motion `I θ̈ = τ(θ)`, with its scalar form `equationOfMotion_iff_scalar` and its independence of the mass `equationOfMotion_iff_of_eq_ω`; `SimplePendulum.IsSolution` is a smooth solution of it. - `SimplePendulum.gradLagrangian` is the variational derivative of the action, computed by `gradLagrangian_eq_eulerLagrangeOp` and `gradLagrangian_eq_torque`.
iii. Table of contents
- A. The input data - A.1. The structure of the input data - A.2. Simple inequalities for the input data - B. Frequency and moment of inertia - B.1. The angular frequency - B.2. The moment of inertia - C. The energies - C.1. The definitions of the energies - C.2. Simple equalities and bounds for the energies - C.3. Smoothness of the energies and the gradient of the potential - C.4. Time derivatives of the energies - D. The Lagrangian - D.1. The definition of the Lagrangian and equalities for it - D.2. Smoothness of the Lagrangian - D.3. Gradients of the Lagrangian - E. The torque and the equation of motion - E.1. The torque - E.2. The equation of motion - E.3. Smooth solutions - E.4. The scalar equation and independence of the mass - F. The variational derivative of the action - F.1. The definition of the variational derivative - F.2. Equality with the Euler–Lagrange operator - F.3. The variational derivative in terms of the torque
iv. References
References for the simple gravity pendulum include: - Landau & Lifshitz, Mechanics, 3rd ed., §5 and §21. - Arnold, Mathematical Methods of Classical Mechanics, 2nd ed., §4.
A. The input data
We start by defining a structure containing the input data of the simple pendulum, and proving basic properties thereof. The input data consists of the mass `m` of the bob, the length `ℓ` of the rod, and the gravitational acceleration `g`; everything else in this file is built from these three numbers.
A.1. The structure of the input data
The three numbers are carried by a structure, together with the positivity assumptions: a pendulum with a massless bob, a rod of zero length or no gravity is not a pendulum.
A.2. Simple inequalities for the input data
The positivity of the input data is used most often through the corresponding non-vanishing statements, which is the form in which the field-clearing tactics consume it.
B. Frequency and moment of inertia
Two derived quantities control the dynamics of the pendulum: the angular frequency `ω = √(g/ℓ)` of the small oscillations about the bottom of the swing, and the moment of inertia `I = m ℓ²` of the bob about the pivot.
The mass enters the equation of motion only through `I`, where it cancels against the mass in the torque of gravity; what survives is `ω`. The identity performing that cancellation is `ω_sq_mul_inertia`.
B.1. The angular frequency
Linearizing `sin θ ≈ θ` about the bottom of the swing turns the equation of motion into that of a harmonic oscillator of angular frequency `√(g/ℓ)`. The exact motion is not harmonic, but this frequency is the natural time scale of the pendulum and appears throughout its analysis.
B.2. The moment of inertia
The bob is a point mass at the fixed distance `ℓ` from the pivot, so the moment of inertia of the pendulum about the pivot is `m ℓ²`. It is the coefficient relating the angular acceleration to the torque, and so plays for the angle the role that the mass plays for a position.
C. The energies
The simple pendulum has a kinetic energy determined by the rate of change of its angle, and a potential energy determined by the height of the bob, hence by the angle itself. These combine to give the total energy of the pendulum.
Here we state and prove a number of properties of these energies, including the gradient of the potential energy, which is the object entering the equation of motion.
C.1. The definitions of the energies
We define the three energies; it is these energies which control the dynamics of the pendulum, through the Lagrangian.
C.2. Simple equalities and bounds for the energies
Besides the definitional unfoldings, the potential energy of the pendulum is non-negative and vanishes exactly at the bottom of the swing, just as the potential energy of the harmonic oscillator is non-negative and vanishes exactly at the origin. What has no harmonic-oscillator analogue is the upper bound: the potential energy of the pendulum is at most `2 m g ℓ`, its value at the top of the swing.
C.3. Smoothness of the energies and the gradient of the potential
The potential energy is a smooth function of the angle, and its gradient on the one-dimensional Euclidean lift is `m g ℓ sin θ` times the unit vector of the angular coordinate. This gradient is what the equation of motion balances against the angular acceleration, so we record it here, once. The subsection also records that, along a smooth lift of the angle, each of the three energies is a differentiable function of the time — differentiability in time along the lift, as distinct from the differentiability in the angle of the potential energy — which is the differentiability that the time derivatives of section C.4 consume.
C.4. Time derivatives of the energies
For a general smooth lift of the angle, which need not satisfy the equation of motion, we can compute the time derivatives of the energies. Each is an inner product against the angular velocity: the equation of motion will be exactly the statement that the two contributions cancel.
D. The Lagrangian
The pendulum is a conservative system, so its Lagrangian is the kinetic energy minus the potential energy, `L = ½ I θ̇² - m g ℓ (1 - cos θ)`. As for the harmonic oscillator, it is defined as a function on phase space, of the time, the angle and the angular velocity separately; that it is `T - V` along a lift of the angle is then a lemma rather than the definition.
The Lagrangian carries no explicit time dependence, the pendulum being autonomous; the time argument is kept because it is the type the Euler–Lagrange operator of Physlib expects.
D.1. The definition of the Lagrangian and equalities for it
The Lagrangian is written directly in terms of the moment of inertia and the potential energy, so that the equalities below are the two ways of reading it: expanded in the input data, and as the kinetic energy minus the potential energy along a lift of the angle.
D.2. Smoothness of the Lagrangian
The Lagrangian is a smooth function of all of its arguments jointly. This is the hypothesis that the Euler–Lagrange theorem of Physlib places on a Lagrangian, so it is recorded on the uncurried form `↿S.lagrangian`.
D.3. Gradients of the Lagrangian
The Euler–Lagrange operator is built from the two partial gradients of the Lagrangian. The gradient in the angle is minus the gradient of the potential energy, that is the torque of section E; the gradient in the angular velocity is the angular momentum `I θ̇`.
E. The torque and the equation of motion
Gravity exerts on the bob a torque `-m g ℓ sin θ` about the pivot, the generalized force conjugate to the angle, and the equation of motion balances it against the rate of change `I θ̈` of the angular momentum.
We take that pointwise relation as the definition of the equation of motion, rather than the vanishing of the variational derivative of the action, which is how the harmonic oscillator defines its own. The reason is that the variational derivative is defined to be `0` whenever no variational gradient exists, so its vanishing holds vacuously for every lift of the angle too rough to admit one; it says what it is meant to say only under a smoothness assumption. The pointwise equation is totalized too — `∂ₜ` is `fderiv`, which is `0` off differentiability — but its totalization cannot make the equation vacuously true: both sides remain genuine, and generally unequal, functions of time. A rough lift can still satisfy the equation accidentally — a discontinuous lift hopping between equilibrium angles solves it, as section E.3 explains — which is why the notion of a solution, `IsSolution`, demands smoothness as well. It is also the form in which the equation of motion is solved and used. The two agree for smooth lifts, by `equationOfMotion_iff_gradLagrangian_zero`, proved in a subsequent contribution; the present module goes as far as `gradLagrangian_eq_torque`, from which that equivalence is one rearrangement away.
E.1. The torque
The pendulum is conservative, so the generalized force conjugate to the angle is minus the gradient of the potential energy. It is a torque about the pivot rather than a force, the angle being the coordinate; this is why it is `m g ℓ sin θ` and not `m g sin θ`.
E.2. The equation of motion
The equation of motion of the simple pendulum equates the rate of change of the angular momentum about the pivot with the torque of gravity, at every instant.
E.3. Smooth solutions
A solution of the pendulum is a smooth lift satisfying the equation of motion. Smoothness is part of the definition because the bare pointwise equation, being totalized, admits unphysical solutions: a lift jumping between the equilibrium angles `0` and `π` has zero torque everywhere, and — being locally constant wherever it is differentiable at all — it has `∂ₜ θ`, and hence `∂ₜ (∂ₜ θ)`, identically zero, so it satisfies the equation even when it is nowhere continuous. Demanding smoothness excludes such junk, and is the regularity under which the variational description of the motion agrees with the pointwise one.
E.4. The scalar equation and independence of the mass
The angle is a single number, so the vector equation of motion is equivalent to the scalar equation obtained by reading off its one component. Dividing that component by the moment of inertia, using `ω_sq_mul_inertia`, cancels the mass and leaves `θ̈ + ω² sin θ = 0`: two pendulums with the same `ω = √(g/ℓ)` have exactly the same angular motions, whatever their masses.
F. The variational derivative of the action
The action of the simple pendulum is the time integral of the Lagrangian along a lift of the angle. Its variational derivative is computed here, in two steps: it is the Euler–Lagrange operator of the Lagrangian, and that operator is the torque minus the rate of change of the angular momentum.
F.1. The definition of the variational derivative
The variational derivative is that of Physlib's variational calculus, applied to the action of the pendulum. Recall that it is defined to be `0` when no variational gradient exists, so the lemmas below are stated for smooth lifts of the angle.
F.2. Equality with the Euler–Lagrange operator
For a smooth lift of the angle the variational derivative of the action is the Euler–Lagrange operator of the Lagrangian, by the general theorem `euler_lagrange_varGradient`; the hypotheses of that theorem are the smoothness of the lift and `contDiff_lagrangian`.
F.3. The variational derivative in terms of the torque
Evaluating the Euler–Lagrange operator with the gradients of section D.3 gives the variational derivative as the torque minus the rate of change of the angular momentum. Its vanishing is therefore the equation of motion of section E; that equivalence, `equationOfMotion_iff_gradLagrangian_zero`, is proved in a subsequent contribution, together with energy conservation, so that this module carries the model of the pendulum alone.
54 declarations
The Mass of a Simple Pendulum is Non-zero
For a simple pendulum with mass , the mass is non-zero, i.e., .
The length of the rod is non-zero
For a simple pendulum , the length of the rod is non-zero, i.e., .
for a simple pendulum
For a simple gravity pendulum, the gravitational acceleration is non-zero ().
Angular frequency of a simple pendulum
For a simple pendulum with a rod of length and subject to a uniform gravitational acceleration , the angular frequency is defined as This quantity represents the angular frequency of the pendulum's small oscillations about the bottom of its swing.
For a simple pendulum, the angular frequency is strictly positive ().
The angular frequency of the simple pendulum is not equal to zero, i.e., .
For a simple pendulum with rod length and gravitational acceleration , the square of its angular frequency is equal to the ratio of the gravitational acceleration to the rod length:
Inverse Square of Angular Frequency
For a simple pendulum with rod length and gravitational acceleration , the inverse of the square of its angular frequency is given by
Moment of inertia of a simple pendulum
For a simple gravity pendulum with mass and rod length , its moment of inertia about the pivot is defined as .
The Moment of Inertia
For a simple gravity pendulum , the moment of inertia about its pivot is strictly positive, i.e., .
The moment of inertia is non-zero ()
For a simple gravity pendulum, the moment of inertia about the pivot is not equal to zero ().
For a simple gravity pendulum with mass , rod length , and gravitational acceleration , let be the angular frequency and be the moment of inertia about the pivot. The product of the square of the angular frequency and the moment of inertia is equal to :
Kinetic energy of a simple pendulum
For a simple pendulum with moment of inertia , given a trajectory representing the angular position in a 1-dimensional Euclidean space, the kinetic energy at time is defined as: where denotes the time derivative of the angle at time , and is the standard inner product on the 1-dimensional Euclidean space.
Potential energy of a simple pendulum
For a simple gravity pendulum with mass , rod length , and gravitational acceleration , the potential energy at an angular position (where denotes the angle from the downward vertical) is given by This value represents the work done against gravity to raise the bob from the equilibrium position (), where the potential energy is normalized to vanish.
Total energy of a simple pendulum
For a simple pendulum, given a trajectory representing the angular position in a 1-dimensional Euclidean space, the total energy at time is defined as the sum of the kinetic energy and the potential energy : where is the kinetic energy of the trajectory at time and is the potential energy at the angular position .
Kinetic Energy of a Simple Pendulum equals
For a simple pendulum with moment of inertia , the kinetic energy of an angular trajectory at any time is given by where denotes the time derivative of the angle at time , and is the standard inner product on the 1-dimensional Euclidean space.
Potential energy of a simple pendulum is
For a simple pendulum with mass , rod length , and gravitational acceleration , the potential energy at a configuration (where represents the angle from the downward vertical) is given by:
Total Energy of a Simple Pendulum is
For a simple pendulum, given a trajectory representing the angular position in a 1-dimensional Euclidean space, the total energy at time is the sum of the kinetic energy and the potential energy : where is the kinetic energy of the trajectory at time and is the potential energy at the angular position .
The potential energy of a simple pendulum is non-negative:
For a simple gravity pendulum with mass , rod length , and gravitational acceleration , the potential energy at any angular position is non-negative: This reflects the fact that the potential energy is normalized to vanish at the bottom of the swing (), which is the lowest point of the bob's trajectory.
The potential energy of a simple pendulum is at most
For a simple gravity pendulum with mass , rod length , and gravitational acceleration , the potential energy at any angular position satisfies the inequality where . This maximum value is attained at the top of the swing.
For a simple pendulum with angle , the potential energy is equal to if and only if .
The Potential Energy of a Simple Pendulum is Smooth ()
The potential energy of a simple pendulum is a smooth function of the angle . Specifically, for any , the potential energy function is -times continuously differentiable.
The Potential Energy is Differentiable
The potential energy of a simple pendulum, where is the mass of the bob, is the length of the rod, is the gravitational acceleration, and is the angle from the downward vertical, is a differentiable function with respect to .
Gradient of the Potential Energy for a Simple Pendulum
For a simple gravity pendulum with mass , rod length , and gravitational acceleration , let be a point in the one-dimensional Euclidean configuration space representing the angular position, where denotes the angle from the downward vertical. The gradient of the potential energy at is given by where is the unit basis vector of the angular coordinate.
Smoothness of Implies Differentiability of Kinetic Energy in Time
For a simple pendulum, let be a smooth () trajectory representing the angular position as a function of time. Then the kinetic energy is differentiable with respect to time .
Smoothness of Implies Differentiability of Potential Energy in Time
For a simple pendulum with mass , rod length , and gravitational acceleration , let be a smooth () function mapping time to the angular position in the configuration space. Then the potential energy as a function of time, , is differentiable with respect to time .
Smoothness of implies differentiability of total energy in time
For a simple pendulum , let be a smooth () trajectory representing the angular position as a function of time. Then the total energy is differentiable with respect to time .
for the Simple Pendulum
Consider a simple pendulum with moment of inertia . For a smooth trajectory representing the angular position, the time derivative of the kinetic energy is given by the inner product of the angular velocity and the rate of change of the angular momentum : where and denote the first and second time derivatives of at time , respectively.
for the Simple Pendulum
For a simple pendulum , let be a smooth function representing the angular position over time. The time derivative of the potential energy evaluated along the trajectory is given by the inner product of the angular velocity and the gradient of the potential energy :
for the Simple Pendulum
For a simple pendulum with moment of inertia and potential energy , let be a smooth trajectory representing the angular position. The time derivative of the total energy is given by the inner product of the angular velocity and the sum of the term and the gradient of the potential energy : where and denote the first and second time derivatives of the trajectory at time , respectively.
Lagrangian of a simple pendulum
The Lagrangian of a simple gravity pendulum is a function of time , angular position , and angular velocity , defined as the kinetic energy minus the potential energy: where is the moment of inertia of the bob about the pivot, is the squared norm of the angular velocity vector in the one-dimensional Euclidean space, and is the potential energy at angular position . The Lagrangian is autonomous, meaning it does not explicitly depend on the time parameter .
The Lagrangian of a simple pendulum with mass , rod length , gravitational acceleration , and moment of inertia , is given by: where and are the angular position and velocity vectors in a one-dimensional Euclidean space, and denotes the single component of the position vector.
Lagrangian of a Simple Pendulum is
For a simple gravity pendulum, let be a trajectory representing the lift of the angular position. At any time , the Lagrangian evaluated at the position and angular velocity is equal to the kinetic energy of the trajectory at that time minus the potential energy at that position: where denotes the time derivative of at time .
The Simple Pendulum Lagrangian is Smooth
The Lagrangian of a simple gravity pendulum, considered as a function of time , angular position , and angular velocity , is smooth. Specifically, for any , the uncurried mapping is -times continuously differentiable over .
The Position Gradient of the Simple Pendulum Lagrangian is
For a simple gravity pendulum with mass , rod length , and gravitational acceleration , let be the Lagrangian. For any time , angular position , and angular velocity , the gradient of the Lagrangian with respect to the angular position is given by where is the angular coordinate (the first component of the vector ) and is the unit basis vector of the one-dimensional Euclidean space.
for a simple pendulum
For a simple gravity pendulum with moment of inertia , given a time , an angular position , and an angular velocity (both represented as vectors in a one-dimensional Euclidean space), the gradient of the Lagrangian with respect to the velocity is equal to the product of the moment of inertia and the velocity: This quantity represents the angular momentum of the pendulum about its pivot.
Torque of a simple pendulum
For a simple gravity pendulum, the torque (or generalized force) at an angular configuration is defined as the negative gradient of the potential energy with respect to : where is represented as an element of a one-dimensional Euclidean space (the angle ).
Torque of a Simple Pendulum
For a simple gravity pendulum with mass , rod length , and gravitational acceleration , the torque at an angular configuration (where is the angle from the downward vertical) is given by where is the unit vector of the angular coordinate.
The single component of the torque is
For a simple gravity pendulum with mass , rod length , and gravitational acceleration , let be the configuration in the one-dimensional Euclidean space representing the angle from the downward vertical. The single component of the torque is given by
Equation of motion of the simple pendulum:
Let be a function representing the angle of a simple pendulum over time. The equation of motion is the condition that at every instant , the product of the moment of inertia and the second time derivative of the angle is equal to the torque exerted by gravity: where is the gravitational torque about the pivot.
The Equation of Motion is equivalent to
For a simple gravity pendulum with moment of inertia and potential energy , a function representing the angle satisfies the equation of motion if and only if, at every instant , the sum of the product of the moment of inertia and the second time derivative of the angle, , and the gradient of the potential energy vanishes: This represents the rotational form of Newton's second law for the pendulum.
Smooth solution of the simple pendulum equation
A function representing the angle of a simple pendulum is a solution if it satisfies two conditions: 1. It is infinitely differentiable, i.e., . 2. It satisfies the equation of motion for all , where is the moment of inertia and is the gravitational torque about the pivot.
Simple Pendulum Solutions are Smooth
Let be a simple pendulum. If a function is a solution to the simple pendulum's equation of motion (as defined by `IsSolution`), then is infinitely differentiable, which is denoted as .
A solution to the simple pendulum satisfies the equation of motion
Let be a simple pendulum with mass , rod length , and gravitational acceleration . If a function is a solution for , then satisfies the equation of motion at every instant , where is the moment of inertia of the bob about the pivot and is the gravitational torque.
Equation of motion of a simple pendulum
Let be a simple pendulum with angular frequency . For a function representing the angular configuration of the pendulum, the equation of motion (defined by the balance of angular momentum and gravitational torque, ) holds if and only if at every time , the scalar equation is satisfied, where is the second time derivative of the angle and is the scalar component of the angular lift.
Equivalence of Equations of Motion for Pendulums with Equal
Consider two simple pendulums, and , and let and be their respective angular frequencies (defined by ). If , then for any function representing the angle of the pendulum over time, satisfies the equation of motion for if and only if it satisfies the equation of motion for . This implies that the angular motion of a simple pendulum is determined solely by its angular frequency , and is independent of its mass .
Variational gradient of the action
For a simple pendulum , let be a function representing the angular position of the bob over time. The variational gradient is the function mapping to the functional derivative of the action integral , where is the Lagrangian of the pendulum. If the path is not sufficiently smooth for the gradient to exist, the function returns the zero function.
The Variational Gradient of the Action Equals the Euler-Lagrange Operator for a Simple Pendulum
For a simple pendulum and a smooth trajectory of the angle , the variational gradient of the action is equal to the Euler-Lagrange operator applied to the pendulum's Lagrangian .
for the Simple Pendulum
For a simple pendulum with moment of inertia and torque function , let be an infinitely differentiable function representing the angular position over time. The variational derivative of the action is given by the difference between the torque and the product of the moment of inertia and the angular acceleration: where denotes the second time derivative of .
for the Simple Pendulum
For a simple pendulum with moment of inertia and gravitational torque , let be an infinitely differentiable () trajectory of the pendulum's angle. Then satisfies the equation of motion if and only if the variational derivative of the action vanishes: This expresses that the smooth motions of the pendulum are precisely the critical points of its action.
is a solution iff and
For a simple gravity pendulum , let be a function representing the angular position over time. Then is a solution to the equation of motion if and only if is infinitely differentiable () and the variational derivative of the action is identically zero.
Conservation of Energy for the Simple Pendulum
For a simple pendulum, let be a smooth () trajectory representing the angular position. If satisfies the equation of motion , then the time derivative of the total energy along this trajectory is zero: where is the sum of the kinetic and potential energies of the pendulum.
Energy Conservation for the Simple Pendulum:
For a simple gravity pendulum , let be a smooth trajectory representing the angle of the bob over time. If satisfies the equation of motion , then for any time , the total energy is equal to the initial energy .
Energy Conservation for Simple Pendulum Solutions ()
For a simple pendulum , let be a smooth solution to the equation of motion. Then for any time , the total energy of the pendulum at time is equal to its initial value at time : where is defined as the sum of the kinetic and potential energies of the pendulum at time .
