Physlib

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

theorem

The Mass mm of a Simple Pendulum is Non-zero

For a simple pendulum with mass mm, the mass is non-zero, i.e., m0m \neq 0.

theorem

The length \ell of the rod is non-zero

For a simple pendulum SS, the length \ell of the rod is non-zero, i.e., S.0S.\ell \neq 0.

theorem

g0g \neq 0 for a simple pendulum

For a simple gravity pendulum, the gravitational acceleration gg is non-zero (g0g \neq 0).

definition

Angular frequency ω=g/\omega = \sqrt{g/\ell} of a simple pendulum

For a simple pendulum with a rod of length \ell and subject to a uniform gravitational acceleration gg, the angular frequency ω\omega is defined as ω=g\omega = \sqrt{\frac{g}{\ell}} This quantity represents the angular frequency of the pendulum's small oscillations about the bottom of its swing.

theorem

ω>0\omega > 0

For a simple pendulum, the angular frequency ω\omega is strictly positive (ω>0\omega > 0).

theorem

ω0\omega \neq 0

The angular frequency ω\omega of the simple pendulum is not equal to zero, i.e., ω0\omega \neq 0.

theorem

ω2=g/\omega^2 = g / \ell

For a simple pendulum with rod length \ell and gravitational acceleration gg, the square of its angular frequency ω\omega is equal to the ratio of the gravitational acceleration to the rod length: ω2=g\omega^2 = \frac{g}{\ell}

theorem

Inverse Square of Angular Frequency ω2=/g\omega^{-2} = \ell/g

For a simple pendulum with rod length \ell and gravitational acceleration gg, the inverse of the square of its angular frequency ω\omega is given by (ω2)1=g (\omega^2)^{-1} = \frac{\ell}{g}

definition

Moment of inertia II of a simple pendulum

For a simple gravity pendulum with mass mm and rod length \ell, its moment of inertia II about the pivot is defined as I=m2I = m \ell^2.

theorem

The Moment of Inertia I>0I > 0

For a simple gravity pendulum SS, the moment of inertia II about its pivot is strictly positive, i.e., I>0I > 0.

theorem

The moment of inertia II is non-zero (I0I \neq 0)

For a simple gravity pendulum, the moment of inertia II about the pivot is not equal to zero (I0I \neq 0).

theorem

ω2I=mg\omega^2 I = m g \ell

For a simple gravity pendulum with mass mm, rod length \ell, and gravitational acceleration gg, let ω\omega be the angular frequency and II 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 mgm g \ell: ω2I=mg\omega^2 I = m g \ell

definition

Kinetic energy T=12Iθ˙2T = \frac{1}{2} I \|\dot{\theta}\|^2 of a simple pendulum

For a simple pendulum with moment of inertia II, given a trajectory θ:TimeR1\theta : \text{Time} \to \mathbb{R}^1 representing the angular position in a 1-dimensional Euclidean space, the kinetic energy at time tt is defined as: T(t)=12Iθ˙(t),θ˙(t)T(t) = \frac{1}{2} I \langle \dot{\theta}(t), \dot{\theta}(t) \rangle where θ˙(t)\dot{\theta}(t) denotes the time derivative of the angle θ\theta at time tt, and ,\langle \cdot, \cdot \rangle is the standard inner product on the 1-dimensional Euclidean space.

definition

Potential energy of a simple pendulum V(θ)=mg(1cosθ)V(\theta) = m g \ell (1 - \cos \theta)

For a simple gravity pendulum with mass m m , rod length \ell , and gravitational acceleration g g , the potential energy at an angular position x x (where x0 x_0 denotes the angle θ \theta from the downward vertical) is given by V(θ)=mg(1cosθ) V(\theta) = m g \ell (1 - \cos \theta) This value represents the work done against gravity to raise the bob from the equilibrium position (θ=0 \theta = 0 ), where the potential energy is normalized to vanish.

definition

Total energy of a simple pendulum E=T+VE = T + V

For a simple pendulum, given a trajectory θ:TimeR1\theta : \text{Time} \to \mathbb{R}^1 representing the angular position in a 1-dimensional Euclidean space, the total energy EE at time tt is defined as the sum of the kinetic energy TT and the potential energy VV: E(t)=T(t)+V(t) E(t) = T(t) + V(t) where T(t)T(t) is the kinetic energy of the trajectory θ\theta at time tt and V(t)V(t) is the potential energy at the angular position θ(t)\theta(t).

theorem

Kinetic Energy of a Simple Pendulum equals 12Iθ˙,θ˙\frac{1}{2} I \langle \dot{\theta}, \dot{\theta} \rangle

For a simple pendulum with moment of inertia II, the kinetic energy TT of an angular trajectory θ:TimeR1\theta: \text{Time} \to \mathbb{R}^1 at any time tt is given by T(θ)(t)=12Iθ˙(t),θ˙(t)T(\theta)(t) = \frac{1}{2} I \langle \dot{\theta}(t), \dot{\theta}(t) \rangle where θ˙(t)\dot{\theta}(t) denotes the time derivative of the angle at time tt, and ,\langle \cdot, \cdot \rangle is the standard inner product on the 1-dimensional Euclidean space.

theorem

Potential energy of a simple pendulum is mg(1cosθ)m g \ell (1 - \cos \theta)

For a simple pendulum with mass mm, rod length \ell, and gravitational acceleration gg, the potential energy at a configuration xR1x \in \mathbb{R}^1 (where x0x_0 represents the angle θ\theta from the downward vertical) is given by: V(x)=mg(1cosθ)V(x) = m g \ell (1 - \cos \theta)

theorem

Total Energy of a Simple Pendulum is E=T+VE = T + V

For a simple pendulum, given a trajectory θ:TimeR1\theta : \text{Time} \to \mathbb{R}^1 representing the angular position in a 1-dimensional Euclidean space, the total energy EE at time tt is the sum of the kinetic energy TT and the potential energy VV: E(θ)(t)=T(θ)(t)+V(θ(t)) E(\theta)(t) = T(\theta)(t) + V(\theta(t)) where T(θ)(t)T(\theta)(t) is the kinetic energy of the trajectory θ\theta at time tt and V(θ(t))V(\theta(t)) is the potential energy at the angular position θ(t)\theta(t).

theorem

The potential energy of a simple pendulum is non-negative: V(θ)0V(\theta) \ge 0

For a simple gravity pendulum with mass mm, rod length \ell, and gravitational acceleration gg, the potential energy V(θ)V(\theta) at any angular position θ\theta is non-negative: V(θ)0 V(\theta) \geq 0 This reflects the fact that the potential energy is normalized to vanish at the bottom of the swing (θ=0\theta = 0), which is the lowest point of the bob's trajectory.

theorem

The potential energy of a simple pendulum is at most 2mg2mg\ell

For a simple gravity pendulum with mass mm, rod length \ell, and gravitational acceleration gg, the potential energy V(θ)V(\theta) at any angular position θ\theta satisfies the inequality V(θ)2mg V(\theta) \leq 2 m g \ell where V(θ)=mg(1cosθ)V(\theta) = m g \ell (1 - \cos \theta). This maximum value is attained at the top of the swing.

theorem

V(θ)=0    cosθ=1V(\theta) = 0 \iff \cos \theta = 1

For a simple pendulum with angle θ\theta, the potential energy V(θ)V(\theta) is equal to 00 if and only if cosθ=1\cos \theta = 1.

theorem

The Potential Energy of a Simple Pendulum is Smooth (CC^\infty)

The potential energy V(θ)=mg(1cosθ)V(\theta) = m g \ell (1 - \cos \theta) of a simple pendulum is a smooth function of the angle θ\theta. Specifically, for any nN{}n \in \mathbb{N} \cup \{\infty\}, the potential energy function is nn-times continuously differentiable.

theorem

The Potential Energy V(θ)V(\theta) is Differentiable

The potential energy V(θ)=mg(1cosθ)V(\theta) = m g \ell (1 - \cos \theta) of a simple pendulum, where mm is the mass of the bob, \ell is the length of the rod, gg is the gravitational acceleration, and θ\theta is the angle from the downward vertical, is a differentiable function with respect to θ\theta.

theorem

Gradient of the Potential Energy VV for a Simple Pendulum

For a simple gravity pendulum with mass mm, rod length \ell, and gravitational acceleration gg, let xx be a point in the one-dimensional Euclidean configuration space representing the angular position, where x0x_0 denotes the angle θ\theta from the downward vertical. The gradient of the potential energy VV at xx is given by V(x)=(mgsin(x0))e0, \nabla V(x) = (m g \ell \sin(x_0)) \mathbf{e}_0, where e0\mathbf{e}_0 is the unit basis vector of the angular coordinate.

theorem

Smoothness of θ(t)\theta(t) Implies Differentiability of Kinetic Energy T(θ(t))T(\theta(t)) in Time

For a simple pendulum, let θ(t)\theta(t) be a smooth (CC^\infty) trajectory representing the angular position as a function of time. Then the kinetic energy T(t)=12Iθ˙(t)2T(t) = \frac{1}{2} I \|\dot{\theta}(t)\|^2 is differentiable with respect to time tt.

theorem

Smoothness of θ(t)\theta(t) Implies Differentiability of Potential Energy V(θ(t))V(\theta(t)) in Time

For a simple pendulum with mass mm, rod length \ell, and gravitational acceleration gg, let θ(t)\theta(t) be a smooth (CC^\infty) function mapping time to the angular position in the configuration space. Then the potential energy as a function of time, V(θ(t))=mg(1cosθ(t))V(\theta(t)) = m g \ell (1 - \cos \theta(t)), is differentiable with respect to time tt.

theorem

Smoothness of θ(t)\theta(t) implies differentiability of total energy E(t)E(t) in time

For a simple pendulum SS, let θ:TimeR1\theta : \text{Time} \to \mathbb{R}^1 be a smooth (CC^\infty) trajectory representing the angular position as a function of time. Then the total energy E(t)=T(t)+V(t)E(t) = T(t) + V(t) is differentiable with respect to time tt.

theorem

ddtT=θ˙,Iθ¨\frac{d}{dt} T = \langle \dot{\theta}, I \ddot{\theta} \rangle for the Simple Pendulum

Consider a simple pendulum with moment of inertia II. For a smooth trajectory θ:TimeR1\theta : \text{Time} \to \mathbb{R}^1 representing the angular position, the time derivative of the kinetic energy TT is given by the inner product of the angular velocity θ˙\dot{\theta} and the rate of change of the angular momentum Iθ¨I \ddot{\theta}: ddtT(t)=θ˙(t),Iθ¨(t) \frac{d}{dt} T(t) = \langle \dot{\theta}(t), I \ddot{\theta}(t) \rangle where θ˙(t)\dot{\theta}(t) and θ¨(t)\ddot{\theta}(t) denote the first and second time derivatives of θ\theta at time tt, respectively.

theorem

tV(θ)=tθ,V(θ)\partial_t V(\theta) = \langle \partial_t \theta, \nabla V(\theta) \rangle for the Simple Pendulum

For a simple pendulum SS, let θ:TimeR1\theta : \text{Time} \to \mathbb{R}^1 be a smooth function representing the angular position over time. The time derivative of the potential energy VV evaluated along the trajectory θ(t)\theta(t) is given by the inner product of the angular velocity tθ(t)\partial_t \theta(t) and the gradient of the potential energy V(θ(t))\nabla V(\theta(t)): t(V(θ(t)))=tθ(t),V(θ(t)) \partial_t (V(\theta(t))) = \langle \partial_t \theta(t), \nabla V(\theta(t)) \rangle

theorem

dEdt=θ˙,Iθ¨+V(θ)\frac{dE}{dt} = \langle \dot{\theta}, I \ddot{\theta} + \nabla V(\theta) \rangle for the Simple Pendulum

For a simple pendulum with moment of inertia II and potential energy VV, let θ:TimeR1\theta : \text{Time} \to \mathbb{R}^1 be a smooth trajectory representing the angular position. The time derivative of the total energy E=T+VE = T + V is given by the inner product of the angular velocity θ˙\dot{\theta} and the sum of the term Iθ¨I \ddot{\theta} and the gradient of the potential energy V\nabla V: ddtE(t)=θ˙(t),Iθ¨(t)+V(θ(t)) \frac{d}{dt} E(t) = \langle \dot{\theta}(t), I \ddot{\theta}(t) + \nabla V(\theta(t)) \rangle where θ˙(t)\dot{\theta}(t) and θ¨(t)\ddot{\theta}(t) denote the first and second time derivatives of the trajectory θ\theta at time tt, respectively.

definition

Lagrangian L(t,x,v)=12Iv2V(x)L(t, x, v) = \frac{1}{2} I \|v\|^2 - V(x) of a simple pendulum

The Lagrangian L(t,x,v)L(t, x, v) of a simple gravity pendulum is a function of time tt, angular position xx, and angular velocity vv, defined as the kinetic energy minus the potential energy: L(t,x,v)=12Iv,vV(x) L(t, x, v) = \frac{1}{2} I \langle v, v \rangle - V(x) where II is the moment of inertia of the bob about the pivot, v,v\langle v, v \rangle is the squared norm of the angular velocity vector vv in the one-dimensional Euclidean space, and V(x)V(x) is the potential energy at angular position xx. The Lagrangian is autonomous, meaning it does not explicitly depend on the time parameter tt.

theorem

L(t,x,v)=12Iv2mg(1cosx0)L(t, x, v) = \frac{1}{2} I \|v\|^2 - m g \ell (1 - \cos x_0)

The Lagrangian L(t,x,v)L(t, x, v) of a simple pendulum with mass mm, rod length \ell, gravitational acceleration gg, and moment of inertia II, is given by: L(t,x,v)=12Iv,vmg(1cosx0) L(t, x, v) = \frac{1}{2} I \langle v, v \rangle - m g \ell (1 - \cos x_0) where xx and vv are the angular position and velocity vectors in a one-dimensional Euclidean space, and x0x_0 denotes the single component of the position vector.

theorem

Lagrangian of a Simple Pendulum is L=TVL = T - V

For a simple gravity pendulum, let θ:TimeR1\theta: \text{Time} \to \mathbb{R}^1 be a trajectory representing the lift of the angular position. At any time tt, the Lagrangian LL evaluated at the position θ(t)\theta(t) and angular velocity θ˙(t)\dot{\theta}(t) is equal to the kinetic energy TT of the trajectory at that time minus the potential energy VV at that position: L(t,θ(t),θ˙(t))=T(θ,t)V(θ(t)) L(t, \theta(t), \dot{\theta}(t)) = T(\theta, t) - V(\theta(t)) where θ˙(t)\dot{\theta}(t) denotes the time derivative of θ\theta at time tt.

theorem

The Simple Pendulum Lagrangian is CC^\infty Smooth

The Lagrangian L(t,x,v)L(t, x, v) of a simple gravity pendulum, considered as a function of time tTimet \in \text{Time}, angular position xR1x \in \mathbb{R}^1, and angular velocity vR1v \in \mathbb{R}^1, is CC^\infty smooth. Specifically, for any n{0,1,,}n \in \{0, 1, \dots, \infty\}, the uncurried mapping (t,x,v)L(t,x,v)(t, x, v) \mapsto L(t, x, v) is nn-times continuously differentiable over R\mathbb{R}.

theorem

The Position Gradient of the Simple Pendulum Lagrangian is (mgsinθ)e0-(m g \ell \sin \theta) \mathbf{e}_0

For a simple gravity pendulum with mass mm, rod length \ell, and gravitational acceleration gg, let L(t,x,v)L(t, x, v) be the Lagrangian. For any time tt, angular position xR1x \in \mathbb{R}^1, and angular velocity vR1v \in \mathbb{R}^1, the gradient of the Lagrangian with respect to the angular position xx is given by xL(t,x,v)=(mgsinθ)e0 \nabla_x L(t, x, v) = -(m g \ell \sin \theta) \mathbf{e}_0 where θ=x0\theta = x_0 is the angular coordinate (the first component of the vector xx) and e0\mathbf{e}_0 is the unit basis vector of the one-dimensional Euclidean space.

theorem

vL=Iv\nabla_v L = I v for a simple pendulum

For a simple gravity pendulum with moment of inertia II, given a time tt, an angular position xx, and an angular velocity vv (both represented as vectors in a one-dimensional Euclidean space), the gradient of the Lagrangian LL with respect to the velocity vv is equal to the product of the moment of inertia and the velocity: vL(t,x,v)=Iv. \nabla_v L(t, x, v) = I v. This quantity represents the angular momentum of the pendulum about its pivot.

definition

Torque of a simple pendulum τ=V\tau = -\nabla V

For a simple gravity pendulum, the torque (or generalized force) τ\tau at an angular configuration xx is defined as the negative gradient of the potential energy VV with respect to xx: τ(x)=V(x) \tau(x) = -\nabla V(x) where xx is represented as an element of a one-dimensional Euclidean space (the angle θ\theta).

theorem

Torque of a Simple Pendulum τ=(mgsinθ)e0\tau = -(m g \ell \sin \theta) \mathbf{e}_0

For a simple gravity pendulum with mass mm, rod length \ell, and gravitational acceleration gg, the torque τ\tau at an angular configuration xx (where x0x_0 is the angle θ\theta from the downward vertical) is given by τ(x)=(mgsinx0)e0\tau(x) = -(m g \ell \sin x_0) \mathbf{e}_0 where e0\mathbf{e}_0 is the unit vector of the angular coordinate.

theorem

The single component of the torque is τ0=mgsinθ\tau_0 = -m g \ell \sin \theta

For a simple gravity pendulum with mass mm, rod length \ell, and gravitational acceleration gg, let xx be the configuration in the one-dimensional Euclidean space representing the angle θ=x0\theta = x_0 from the downward vertical. The single component of the torque τ(x)\tau(x) is given by τ(x)0=mgsinθ\tau(x)_0 = -m g \ell \sin \theta

definition

Equation of motion of the simple pendulum: Iθ¨=τ(θ)I \ddot{\theta} = \tau(\theta)

Let θ:TimeR1\theta : \text{Time} \to \mathbb{R}^1 be a function representing the angle of a simple pendulum over time. The equation of motion is the condition that at every instant tt, the product of the moment of inertia II and the second time derivative of the angle θ¨(t)\ddot{\theta}(t) is equal to the torque τ(θ(t))\tau(\theta(t)) exerted by gravity: Iθ¨(t)=τ(θ(t)) I \ddot{\theta}(t) = \tau(\theta(t)) where τ(θ)=mgsinθ\tau(\theta) = -m g \ell \sin \theta is the gravitational torque about the pivot.

theorem

The Equation of Motion is equivalent to Iθ¨+V=0I \ddot{\theta} + \nabla V = 0

For a simple gravity pendulum with moment of inertia II and potential energy VV, a function θ:TimeR1\theta : \text{Time} \to \mathbb{R}^1 representing the angle satisfies the equation of motion if and only if, at every instant tt, the sum of the product of the moment of inertia and the second time derivative of the angle, Iθ¨(t)I \ddot{\theta}(t), and the gradient of the potential energy V(θ(t))\nabla V(\theta(t)) vanishes: Iθ¨(t)+V(θ(t))=0 I \ddot{\theta}(t) + \nabla V(\theta(t)) = 0 This represents the rotational form of Newton's second law for the pendulum.

definition

Smooth solution of the simple pendulum equation Iθ¨=τ(θ)I \ddot{\theta} = \tau(\theta)

A function θ:TimeR1\theta : \text{Time} \to \mathbb{R}^1 representing the angle of a simple pendulum SS is a solution if it satisfies two conditions: 1. It is infinitely differentiable, i.e., θC(Time,R1)\theta \in C^\infty(\text{Time}, \mathbb{R}^1). 2. It satisfies the equation of motion Iθ¨(t)=τ(θ(t))I \ddot{\theta}(t) = \tau(\theta(t)) for all tTimet \in \text{Time}, where II is the moment of inertia and τ(θ)=mgsinθ\tau(\theta) = -m g \ell \sin \theta is the gravitational torque about the pivot.

theorem

Simple Pendulum Solutions are Smooth

Let SS be a simple pendulum. If a function θ:TimeR1\theta : \text{Time} \to \mathbb{R}^1 is a solution to the simple pendulum's equation of motion (as defined by `IsSolution`), then θ\theta is infinitely differentiable, which is denoted as θC(Time,R1)\theta \in C^\infty(\text{Time}, \mathbb{R}^1).

theorem

A solution to the simple pendulum satisfies the equation of motion Iθ¨=τ(θ)I \ddot{\theta} = \tau(\theta)

Let SS be a simple pendulum with mass mm, rod length \ell, and gravitational acceleration gg. If a function θ:TimeR1\theta : \text{Time} \to \mathbb{R}^1 is a solution for SS, then θ\theta satisfies the equation of motion Iθ¨(t)=τ(θ(t)) I \ddot{\theta}(t) = \tau(\theta(t)) at every instant tt, where I=m2I = m \ell^2 is the moment of inertia of the bob about the pivot and τ(θ)=mgsinθ\tau(\theta) = -m g \ell \sin \theta is the gravitational torque.

theorem

Equation of motion of a simple pendulum     θ¨+ω2sinθ=0\iff \ddot{\theta} + \omega^2 \sin \theta = 0

Let SS be a simple pendulum with angular frequency ω=g/\omega = \sqrt{g/\ell}. For a function θ:TimeR1\theta : \text{Time} \to \mathbb{R}^1 representing the angular configuration of the pendulum, the equation of motion (defined by the balance of angular momentum and gravitational torque, Iθ¨=τ(θ)I \ddot{\theta} = \tau(\theta)) holds if and only if at every time tt, the scalar equation θ¨(t)+ω2sin(θ(t))=0 \ddot{\theta}(t) + \omega^2 \sin(\theta(t)) = 0 is satisfied, where θ¨(t)\ddot{\theta}(t) is the second time derivative of the angle and θ(t)\theta(t) is the scalar component of the angular lift.

theorem

Equivalence of Equations of Motion for Pendulums with Equal ω\omega

Consider two simple pendulums, SS and SS', and let ωS\omega_S and ωS\omega_{S'} be their respective angular frequencies (defined by ω=g/\omega = \sqrt{g/\ell}). If ωS=ωS\omega_{S'} = \omega_S, then for any function θ:TimeR1\theta: \text{Time} \to \mathbb{R}^1 representing the angle of the pendulum over time, θ\theta satisfies the equation of motion for SS' if and only if it satisfies the equation of motion for SS. This implies that the angular motion of a simple pendulum is determined solely by its angular frequency ω\omega, and is independent of its mass mm.

definition

Variational gradient δSδθ\frac{\delta \mathcal{S}}{\delta \theta} of the action

For a simple pendulum SS, let θ:TimeR1\theta: \text{Time} \to \mathbb{R}^1 be a function representing the angular position of the bob over time. The variational gradient is the function mapping tTimet \in \text{Time} to the functional derivative δSδθ(t)\frac{\delta \mathcal{S}}{\delta \theta}(t) of the action integral S[θ]=L(t,θ(t),θ˙(t))dt\mathcal{S}[\theta] = \int L(t, \theta(t), \dot{\theta}(t)) \, dt, where LL is the Lagrangian of the pendulum. If the path θ\theta is not sufficiently smooth for the gradient to exist, the function returns the zero function.

theorem

The Variational Gradient of the Action Equals the Euler-Lagrange Operator for a Simple Pendulum

For a simple pendulum SS and a smooth trajectory of the angle θ:TimeR1\theta: \text{Time} \to \mathbb{R}^1, the variational gradient of the action δSδθ\frac{\delta \mathcal{S}}{\delta \theta} is equal to the Euler-Lagrange operator Lθddt(Lθ˙)\frac{\partial L}{\partial \theta} - \frac{d}{dt} \left( \frac{\partial L}{\partial \dot{\theta}} \right) applied to the pendulum's Lagrangian LL.

theorem

δSδθ=τ(θ)Iθ¨\frac{\delta \mathcal{S}}{\delta \theta} = \tau(\theta) - I \ddot{\theta} for the Simple Pendulum

For a simple pendulum with moment of inertia II and torque function τ\tau, let θ:TimeR1\theta: \text{Time} \to \mathbb{R}^1 be an infinitely differentiable function representing the angular position over time. The variational derivative of the action δSδθ\frac{\delta \mathcal{S}}{\delta \theta} is given by the difference between the torque and the product of the moment of inertia and the angular acceleration: (δSδθ)(t)=τ(θ(t))Iθ¨(t) \left( \frac{\delta \mathcal{S}}{\delta \theta} \right)(t) = \tau(\theta(t)) - I \ddot{\theta}(t) where θ¨\ddot{\theta} denotes the second time derivative of θ\theta.

theorem

Iθ¨=τ(θ)    δSδθ=0I \ddot{\theta} = \tau(\theta) \iff \frac{\delta \mathcal{S}}{\delta \theta} = 0 for the Simple Pendulum

For a simple pendulum SS with moment of inertia II and gravitational torque τ(θ)=mgsinθ\tau(\theta) = -m g \ell \sin \theta, let θ:TimeR1\theta: \text{Time} \to \mathbb{R}^1 be an infinitely differentiable (CC^\infty) trajectory of the pendulum's angle. Then θ\theta satisfies the equation of motion Iθ¨=τ(θ) I \ddot{\theta} = \tau(\theta) if and only if the variational derivative of the action δSδθ\frac{\delta \mathcal{S}}{\delta \theta} vanishes: δSδθ=0. \frac{\delta \mathcal{S}}{\delta \theta} = 0. This expresses that the smooth motions of the pendulum are precisely the critical points of its action.

theorem

θ\theta is a solution iff θC\theta \in C^\infty and δSδθ=0\frac{\delta \mathcal{S}}{\delta \theta} = 0

For a simple gravity pendulum SS, let θ:TimeR1\theta : \text{Time} \to \mathbb{R}^1 be a function representing the angular position over time. Then θ\theta is a solution to the equation of motion if and only if θ\theta is infinitely differentiable (θC\theta \in C^\infty) and the variational derivative of the action δSδθ\frac{\delta \mathcal{S}}{\delta \theta} is identically zero.

theorem

Conservation of Energy for the Simple Pendulum

For a simple pendulum, let θ:TimeR1\theta : \text{Time} \to \mathbb{R}^1 be a smooth (CC^\infty) trajectory representing the angular position. If θ\theta satisfies the equation of motion Iθ¨=τ(θ)I \ddot{\theta} = \tau(\theta), then the time derivative of the total energy EE along this trajectory is zero: t(E(θ))=0 \partial_t (E(\theta)) = 0 where E(θ)E(\theta) is the sum of the kinetic and potential energies of the pendulum.

theorem

Energy Conservation for the Simple Pendulum: E(t)=E(0)E(t) = E(0)

For a simple gravity pendulum SS, let θ:TimeR1\theta : \text{Time} \to \mathbb{R}^1 be a smooth trajectory representing the angle of the bob over time. If θ\theta satisfies the equation of motion Iθ¨=τ(θ)I \ddot{\theta} = \tau(\theta), then for any time tt, the total energy E(t)E(t) is equal to the initial energy E(0)E(0).

theorem

Energy Conservation for Simple Pendulum Solutions (E(t)=E(0)E(t) = E(0))

For a simple pendulum SS, let θ:TimeR1\theta: \text{Time} \to \mathbb{R}^1 be a smooth solution to the equation of motion. Then for any time tt, the total energy EE of the pendulum at time tt is equal to its initial value at time t=0t = 0: E(t)=E(0) E(t) = E(0) where E(t)E(t) is defined as the sum of the kinetic and potential energies of the pendulum at time tt.