Physlib

Physlib.ClassicalMechanics.Pendulum.SimplePendulum.Geometric.PhysicalSpace

The simple pendulum in physical space

i. Overview

The trajectory module places the bob in the plane along a lifted trajectory; this module fixes a simple pendulum `S` and identifies the dynamics of `SimplePendulum.Basic`, written on the Euclidean lift of the angle, with the dynamics of the bob moving in physical space. Along a lift `θ` the bob is at `S.spaceTrajectory θ`, at the fixed distance `ℓ` from the pivot; its velocity is obtained by differentiating the position componentwise, and the square of its speed is `ℓ² θ̇²`. The chart kinetic energy `½ I θ̇²` is then exactly the kinetic energy `½ m ‖v‖²` of the bob, the chart potential energy `m g ℓ (1 - cos θ)` is exactly the gravitational potential `m g h` of the bob with the height `h` measured from the bottom of the swing, and the chart Lagrangian is the constrained Lagrangian `T - V` of a point mass moving on the circle of radius `ℓ` about the pivot. Velocities in this module are physical-space time derivatives; the geometric velocity as a tangent vector to the configuration circle is deferred, as discussed in the trajectory module.

ii. Key results

- `SimplePendulum.spaceTrajectory` : the bob's position in the plane along a lifted trajectory, with the rod-length constraint `SimplePendulum.norm_spaceTrajectory`. - `SimplePendulum.deriv_spaceTrajectory` : the bob's velocity along a lifted trajectory, with the square of its speed `SimplePendulum.norm_sq_deriv_spaceTrajectory`. - `SimplePendulum.kineticEnergy_eq_space`, `SimplePendulum.potentialEnergy_eq_height` and `SimplePendulum.lagrangian_eq_space` : the chart kinetic energy, potential energy and Lagrangian are those of the bob in physical space; `SimplePendulum.energy_eq_space` is the corresponding statement for the total energy.

iii. Table of contents

  • A. The bob's position along a lifted trajectory
  • B. The bob's velocity and speed
  • C. The energies and the Lagrangian in physical space

iv. References

- `Physlib.ClassicalMechanics.Pendulum.SimplePendulum.Basic` (the lifted dynamics: the energies and the Lagrangian on the Euclidean lift of the angle). - `Physlib.ClassicalMechanics.Pendulum.SimplePendulum.Geometric.Trajectory` (trajectories on the configuration circle and the position of the bob along them). - Landau & Lifshitz, Mechanics, 3rd ed., §5, Problems 1–3 (pendulum configurations).

A. The bob's position along a lifted trajectory

Fixing a pendulum `S`, the position of the bob in the plane along the trajectory described by a lift `θ` of the angle is obtained by composing the lifted trajectory with the map to physical space for a rod of length `S.ℓ`. Along it the bob is at `(ℓ sin θ, -ℓ cos θ)`, and since the length of the rod is positive the bob stays at distance `ℓ` from the pivot at all times.

B. The bob's velocity and speed

Differentiating the position of the bob componentwise gives its velocity in the plane: along a differentiable lift the bob's position is differentiable in time, and its velocity is `(ℓ cos θ, ℓ sin θ)` times the angular velocity — the vector tangent to the circle of radius `ℓ`, of length `ℓ |θ̇|`. The square of the speed is therefore `ℓ² θ̇²`, which is the identity behind the identification of the kinetic energies in the next section. The two chain-rule computations for the sine and the cosine of the angle read from a lift are recorded as stand-alone lemmas.

C. The energies and the Lagrangian in physical space

The identities of the previous section identify the energies of `SimplePendulum.Basic`, written on the lift of the angle, with the energies of the bob in the plane. The chart kinetic energy `½ I θ̇²` is the kinetic energy `½ m ‖v‖²` of the bob, since `I = m ℓ²` and the square of the speed is `ℓ² θ̇²`; the chart potential energy `m g ℓ (1 - cos θ)` is the gravitational potential `m g h` of the bob, the height `h` above the bottom of the swing being the vertical position of the bob plus `ℓ`. Consequently the chart Lagrangian is the constrained Lagrangian `T - V` of a point mass moving on the circle of radius `ℓ`, and the chart energy is its total energy `T + V`.

14 declarations

definition

Physical position of the bob along an angular lift θ\theta

For a simple pendulum SS with rod length ℓ\ell, this function maps a time-dependent angular lift θ:Time→R\theta : \text{Time} \to \mathbb{R} to the trajectory of the pendulum bob in 2D physical space R2\mathbb{R}^2. At any time tt, the position of the bob is given by the vector (ℓsin⁡(θ(t)),−ℓcos⁡(θ(t)))(\ell \sin(\theta(t)), -\ell \cos(\theta(t))), where the pivot is located at the origin.

theorem

The physical trajectory S.spaceTrajectory(θ)S.\text{spaceTrajectory}(\theta) equals (ℓsin⁡θ,−ℓcos⁡θ)(\ell \sin \theta, -\ell \cos \theta)

For a simple pendulum SS with rod length ℓ\ell, let θ:Time→R\theta: \text{Time} \to \mathbb{R} be an angular lift representing the angle of the pendulum over time. The physical trajectory of the bob in the 2D plane, denoted by S.spaceTrajectory(θ)S.\text{spaceTrajectory}(\theta), is given by the function t↦(ℓsin⁡(θ(t)),−ℓcos⁡(θ(t))),t \mapsto (\ell \sin(\theta(t)), -\ell \cos(\theta(t))), where the origin is the pivot of the pendulum.

theorem

The horizontal position of the pendulum bob is ℓsin⁡(θ(t))\ell \sin(\theta(t))

For a simple pendulum SS with rod length ℓ\ell, let θ(t)\theta(t) be the time-dependent angular lift. The horizontal coordinate (the component with index 00) of the pendulum bob's position in physical space at time tt is given by ℓsin⁡(θ(t))\ell \sin(\theta(t)).

theorem

Vertical position of the bob is −ℓcos⁡θ-\ell \cos \theta

Let SS be a simple pendulum with rod length ℓ\ell. For any time-dependent angular lift θ:Time→R\theta : \text{Time} \to \mathbb{R}, the vertical position (the second coordinate) of the pendulum bob in 2D physical space at time tt is given by y(t)=−ℓcos⁡(θ(t))y(t) = -\ell \cos(\theta(t)) where the pivot is located at the origin (0,0)(0,0) and the second axis points upwards.

theorem

The distance of the pendulum bob from the pivot is ℓ\ell (i.e., ∥S.spaceTrajectory θ(t)∥=ℓ\| S.\text{spaceTrajectory } \theta(t) \| = \ell)

For a simple pendulum SS with rod length ℓ\ell, let θ:Time→R\theta : \text{Time} \to \mathbb{R} be a time-dependent angular lift. At any time tt, the Euclidean norm of the bob's position in physical space, denoted as ∥S.spaceTrajectory θ(t)∥\| S.\text{spaceTrajectory } \theta(t) \|, is equal to the length of the rod ℓ\ell. This reflects the physical constraint that the bob remains at a fixed distance ℓ\ell from the pivot.

theorem

ddtsin⁡(θ)=cos⁡(θ)θ˙\frac{d}{dt} \sin(\theta) = \cos(\theta) \dot{\theta}

Let θ:Time→R\theta : \text{Time} \to \mathbb{R} be a differentiable angular lift of a simple pendulum. At any time tt, the time derivative of the sine of the angle is the product of the cosine of the angle and the angular velocity: ddtsin⁡(θ(t))=cos⁡(θ(t))⋅θ˙(t)\frac{d}{dt} \sin(\theta(t)) = \cos(\theta(t)) \cdot \dot{\theta}(t) where θ˙(t)\dot{\theta}(t) denotes the time derivative of the angle θ\theta at time tt.

theorem

ddtcos⁡(θ(t))=−sin⁡(θ(t))θ˙(t)\frac{d}{dt} \cos(\theta(t)) = -\sin(\theta(t)) \dot{\theta}(t)

Let θ:Time→R\theta : \text{Time} \to \mathbb{R} be a differentiable function representing the angular lift of a pendulum. For any time tt, the time derivative of the cosine of the angle, cos⁡(θ(t))\cos(\theta(t)), is equal to the negative of the sine of the angle multiplied by the angular velocity: ddtcos⁡(θ(t))=−sin⁡(θ(t))θ˙(t),\frac{d}{dt} \cos(\theta(t)) = -\sin(\theta(t)) \dot{\theta}(t), where θ˙(t)\dot{\theta}(t) denotes the time derivative of θ\theta at time tt.

theorem

Differentiability of the Pendulum Bob's Physical Trajectory

Let SS be a simple pendulum. If the angular lift θ:Time→R\theta : \text{Time} \to \mathbb{R} is a differentiable function of time, then the physical trajectory of the pendulum bob in 2D space, denoted by S.spaceTrajectory(θ)S.\text{spaceTrajectory}(\theta) (which maps time tt to the position (ℓsin⁡(θ(t)),−ℓcos⁡(θ(t)))(\ell \sin(\theta(t)), -\ell \cos(\theta(t)))), is also differentiable with respect to time.

theorem

Velocity of the Pendulum Bob in Physical Space

Let SS be a simple pendulum with rod length ℓ\ell. If the angular lift θ:Time→R\theta : \text{Time} \to \mathbb{R} is a differentiable function of time, then the velocity of the pendulum bob in 2D physical space at time tt is the time derivative of its position S.spaceTrajectory(θ)S.\text{spaceTrajectory}(\theta), given by the vector: ddt(S.spaceTrajectory(θ))(t)=(ℓcos⁡(θ(t))θ˙(t),ℓsin⁡(θ(t))θ˙(t))\frac{d}{dt} (S.\text{spaceTrajectory}(\theta))(t) = \left( \ell \cos(\theta(t)) \dot{\theta}(t), \ell \sin(\theta(t)) \dot{\theta}(t) \right) where θ˙(t)\dot{\theta}(t) denotes the angular velocity (the time derivative of θ\theta) at time tt.

theorem

The square of the pendulum bob's speed is ℓ2θ˙2\ell^2 \dot{\theta}^2

For a simple pendulum SS with rod length ℓ\ell, let θ:Time→R\theta: \text{Time} \to \mathbb{R} be a differentiable angular trajectory. At any time tt, the square of the speed of the pendulum bob in physical space—defined as the squared norm of the time derivative of its position vector r(t)\mathbf{r}(t)—is equal to the product of the square of the rod length and the square of the angular velocity θ˙(t)\dot{\theta}(t): ∥∂t(S.spaceTrajectory θ)(t)∥2=ℓ2(dθdt)2 \|\partial_t (S.\text{spaceTrajectory } \theta)(t)\|^2 = \ell^2 \left( \frac{d\theta}{dt} \right)^2 where S.spaceTrajectory θS.\text{spaceTrajectory } \theta is the function mapping the angular lift to the bob's coordinates (ℓsin⁡θ,−ℓcos⁡θ)(\ell \sin \theta, -\ell \cos \theta) in R2\mathbb{R}^2.

theorem

The Chart Kinetic Energy of a Simple Pendulum Equals 12m∥v∥2\frac{1}{2} m \|\mathbf{v}\|^2 in Physical Space

Consider a simple pendulum SS with mass mm. Let θ:Time→R\theta: \text{Time} \to \mathbb{R} be a differentiable angular trajectory. At any time tt, the chart kinetic energy of the pendulum T(θ,t)T(\theta, t) is equal to the kinetic energy of the bob in physical space, given by: T(θ,t)=12m∥ddtr(t)∥2T(\theta, t) = \frac{1}{2} m \left\| \frac{d}{dt} \mathbf{r}(t) \right\|^2 where r(t)\mathbf{r}(t) is the position of the pendulum bob in 2D physical space at time tt (defined as S.spaceTrajectory(θ)(t)S.\text{spaceTrajectory}(\theta)(t)), and ∥ddtr(t)∥\left\| \frac{d}{dt} \mathbf{r}(t) \right\| is the speed of the bob.

theorem

The potential energy of a simple pendulum V(θ)=mg(y+ℓ)V(\theta) = m g (y + \ell)

For a simple pendulum SS with mass mm, gravitational acceleration gg, and rod length ℓ\ell, let θ(t)\theta(t) be a time-dependent angular trajectory. The chart potential energy V(θ(t))V(\theta(t)) is equal to the gravitational potential energy in physical space, given by V(θ(t))=mg(y(t)+ℓ) V(\theta(t)) = m g (y(t) + \ell) where y(t)y(t) is the vertical component (the second coordinate) of the bob's position (x(t),y(t))(x(t), y(t)) in the plane. Here, y(t)+ℓy(t) + \ell represents the height of the bob above the bottom of its swing.

theorem

The Chart Lagrangian of a Simple Pendulum Equals T−VT - V in Physical Space

Consider a simple pendulum SS with mass mm, rod length ℓ\ell, and gravitational acceleration gg. Let θ:Time→R\theta: \text{Time} \to \mathbb{R} be a differentiable angular trajectory. At any time tt, the chart Lagrangian L(t,θ(t),θ˙(t))L(t, \theta(t), \dot{\theta}(t)) of the pendulum is equal to the kinetic energy of the bob minus its gravitational potential energy in 2D physical space: L(t,θ(t),θ˙(t))=12m∥v(t)∥2−mg(y(t)+ℓ)L(t, \theta(t), \dot{\theta}(t)) = \frac{1}{2} m \|\mathbf{v}(t)\|^2 - m g (y(t) + \ell) where v(t)=ddtr(t)\mathbf{v}(t) = \frac{d}{dt} \mathbf{r}(t) is the velocity of the bob in physical space, r(t)\mathbf{r}(t) is its position vector, and y(t)y(t) is its vertical coordinate (the second component of the position vector). The term y(t)+ℓy(t) + \ell represents the height of the bob above the bottom of its swing.

theorem

The Chart Energy of a Simple Pendulum Equals T+VT + V in Physical Space

Consider a simple pendulum SS with mass mm, rod length ℓ\ell, and gravitational acceleration gg. Let θ:Time→R\theta: \text{Time} \to \mathbb{R} be a differentiable angular trajectory. At any time tt, the chart energy E(t)E(t) of the pendulum is equal to the total energy of the bob in physical space, which is the sum of its kinetic energy and its gravitational potential energy: E(t)=12m∥v(t)∥2+mg(y(t)+ℓ)E(t) = \frac{1}{2} m \|\mathbf{v}(t)\|^2 + m g (y(t) + \ell) where v(t)=ddtr(t)\mathbf{v}(t) = \frac{d}{dt} \mathbf{r}(t) is the velocity of the bob in 2D physical space, r(t)\mathbf{r}(t) is its position vector, and y(t)y(t) is the vertical component (the second coordinate) of its position. The term y(t)+ℓy(t) + \ell represents the height of the bob above the bottom of its swing.