Physlib

Physlib.ClassicalMechanics.RigidBody.Motion

Rigid body motion

The static `RigidBody` records a body-fixed mass distribution. To describe a rigid body *in motion* we record, in addition, the trajectory of its centre of mass in the inertial frame and the body's time-dependent orientation (a rotation about the centre of mass).

From this configuration we define the velocity of the centre of mass and the body's linear momentum. The reference point is taken to be the centre of mass, following the decomposition of a rigid motion into a translation of the centre of mass plus a rotation about it.

References

21 declarations

theorem

R(t)R(t)=IR(t) R(t)^\top = I for Rigid Body Orientation

For a rigid body motion MM in dd-dimensional space, let R(t)R(t) denote the orientation matrix of the body at time tt. Then R(t)R(t)=IR(t) R(t)^\top = I, where II is the d×dd \times d identity matrix.

definition

Velocity of the Center of Mass V(t)=ddtRcom(t)\mathbf{V}(t) = \frac{d}{dt}\mathbf{R}_{\text{com}}(t)

For a rigid body motion MM in dd-dimensional space, the velocity of its center of mass is the function mapping each time tTimet \in \text{Time} to the time derivative of the center-of-mass trajectory Rcom(t)Space d\mathbf{R}_{\text{com}}(t) \in \text{Space } d. Mathematically, this is expressed as: V(t)=ddtRcom(t)\mathbf{V}(t) = \frac{d}{dt} \mathbf{R}_{\text{com}}(t) This velocity corresponds to the translational component V\mathbf{V} in the Landau–Lifshitz decomposition of the velocity of a point on the body, v=V+Ω×r\mathbf{v} = \mathbf{V} + \boldsymbol{\Omega} \times \mathbf{r}, where Ω\boldsymbol{\Omega} is the angular velocity and r\mathbf{r} is the position relative to the center of mass.

theorem

Vcom(t)=ddtRcom(t)\mathbf{V}_{\text{com}}(t) = \frac{d}{dt} \mathbf{R}_{\text{com}}(t)

For a rigid body motion MM in dd-dimensional space, the velocity of the center of mass is equal to the time derivative of the trajectory of the center of mass. This is expressed as: V(t)=ddtRcom(t)\mathbf{V}(t) = \frac{d}{dt} \mathbf{R}_{\text{com}}(t) where V\mathbf{V} is the velocity of the center of mass and Rcom\mathbf{R}_{\text{com}} is the trajectory of the center of mass.

theorem

Rcom(t)=c    V(t)=0\mathbf{R}_{\text{com}}(t) = \mathbf{c} \implies \mathbf{V}(t) = \mathbf{0}

For a rigid body motion MM in dd-dimensional space, if the trajectory of the center of mass is constant (i.e., Rcom(t)=c\mathbf{R}_{\text{com}}(t) = \mathbf{c} for all times tt and some constant vector c\mathbf{c}), then the velocity of the center of mass is zero (V(t)=0\mathbf{V}(t) = \mathbf{0} for all tt).

definition

Linear momentum p(t)=mV(t)\mathbf{p}(t) = m \mathbf{V}(t)

For a rigid body motion MM in dd-dimensional space, the linear momentum is the function that maps each time tTimet \in \text{Time} to the vector obtained by multiplying the body's total mass mm by the velocity of its center of mass V(t)\mathbf{V}(t). Mathematically, the linear momentum p(t)\mathbf{p}(t) is defined as: p(t)=mV(t)\mathbf{p}(t) = m \mathbf{V}(t) where mm is the total mass of the rigid body and V(t)Space d\mathbf{V}(t) \in \text{Space } d is the velocity of the center of mass at time tt.

theorem

p(t)=mV(t)\mathbf{p}(t) = m \mathbf{V}(t)

For a rigid body motion MM in dd-dimensional space, the linear momentum p(t)\mathbf{p}(t) is equal to the product of the total mass mm of the rigid body and the velocity of its center of mass V(t)\mathbf{V}(t). Mathematically, this is expressed as: p(t)=mV(t)\mathbf{p}(t) = m \mathbf{V}(t) where mm is the mass of the rigid body and V(t)\mathbf{V}(t) is the center-of-mass velocity.

definition

Rigid displacement at time tt

For a rigid body motion MM in dd-dimensional space Space d\text{Space } d, the displacement at time tt is the map that carries a point yy from the body-fixed frame to its position in the inertial frame. It is defined by the formula: displacement(t,y)=R(t)(yrcm)+rtraj(t) \text{displacement}(t, y) = \mathbf{R}(t)(y - \mathbf{r}_{\text{cm}}) + \mathbf{r}_{\text{traj}}(t) where: - R(t)\mathbf{R}(t) is the rotation matrix (the `orientation` of the body at time tt), - rcm\mathbf{r}_{\text{cm}} is the fixed `centerOfMass` of the rigid body in the body frame, - rtraj(t)\mathbf{r}_{\text{traj}}(t) is the `comTrajectory` representing the position of the center of mass in the inertial frame at time tt.

theorem

Coordinate-wise expression of rigid displacement

For a rigid body motion MM in dd-dimensional space, the kk-th coordinate of the displacement of a point yy (from the body-fixed frame to the inertial frame) at time tt is given by: (displacement(t,y))k=j=1dRkj(t)(yj(rcm)j)+(rtraj(t))k (\text{displacement}(t, y))_k = \sum_{j=1}^{d} R_{kj}(t) (y_j - (\mathbf{r}_{\text{cm}})_j) + (\mathbf{r}_{\text{traj}}(t))_k where: - Rkj(t)R_{kj}(t) is the (k,j)(k, j)-th entry of the rotation matrix R(t)\mathbf{R}(t) (representing the orientation of the body at time tt), - yjy_j is the jj-th coordinate of the point yy in the body frame, - (rcm)j(\mathbf{r}_{\text{cm}})_j is the jj-th coordinate of the fixed center of mass in the body frame, - (rtraj(t))k(\mathbf{r}_{\text{traj}}(t))_k is the kk-th coordinate of the center of mass trajectory in the inertial frame at time tt.

theorem

The displacement map is CC^\infty

For a rigid body motion MM in dd-dimensional space Space d\text{Space } d and a fixed time tt, the displacement map ydisplacement(t,y)y \mapsto \text{displacement}(t, y), which maps a point yy from the body-fixed frame to its position in the inertial frame, is infinitely differentiable (of class CC^\infty).

definition

Mass distribution of a rigid body motion at time tt

For a rigid body motion MM in dd-dimensional space Space d\text{Space } d, the mass distribution at time tt is defined as the pushforward of the body-fixed mass distribution ρ\rho along the rigid displacement map Φt:Space dSpace d\Phi_t : \text{Space } d \to \text{Space } d. Specifically, the distribution ρt\rho_t acts on a test function ff according to the formula: ρt(f)=ρ(fΦt) \rho_t(f) = \rho(f \circ \Phi_t) where Φt\Phi_t is the displacement at time tt that maps coordinates from the body-fixed frame to the inertial frame.

theorem

Conservation of total mass in rigid body motion

For a rigid body motion MM in dd-dimensional space, the total mass of the mass distribution at any time tt is equal to the total mass of the reference (body-fixed) rigid body. Mathematically, this is expressed as mass(ρt)=mass(ρ)\text{mass}(\rho_t) = \text{mass}(\rho), where ρt\rho_t is the mass distribution at time tt and ρ\rho is the original body-fixed mass distribution.

theorem

The center of mass of a moving rigid body tracks its prescribed trajectory comTrajectory(t)\text{comTrajectory}(t).

For a rigid body motion MM in dd-dimensional space with non-zero total mass, the center of mass of the body's mass distribution at any time tt is equal to the position specified by the prescribed center of mass trajectory comTrajectory(t)\text{comTrajectory}(t) at that time.

definition

Velocity of a material point yy in rigid body motion

Given a rigid body motion MM in dd-dimensional space Space d\text{Space } d and a material point yy in the body-fixed frame, the velocity vy:TimeSpace d\mathbf{v}_y: \text{Time} \to \text{Space } d is the function that assigns to each time tt the time derivative of the trajectory of yy in the inertial frame. It is defined by: vy(t)=ddtdisplacement(t,y) \mathbf{v}_y(t) = \frac{d}{dt} \text{displacement}(t, y) where displacement(t,y)\text{displacement}(t, y) is the position of the material point yy in the inertial frame at time tt.

theorem

Velocity is the time derivative of displacement: vy(t)=ddtdisplacement(t,y)\mathbf{v}_y(t) = \frac{d}{dt} \text{displacement}(t, y)

For a rigid body motion MM in dd-dimensional space Space d\text{Space } d, let yy be a material point in the body-fixed frame and tt be a time. The velocity vy(t)\mathbf{v}_y(t) is equal to the time derivative of the displacement (position in the inertial frame) of yy evaluated at time tt: vy(t)=ddsdisplacement(s,y)s=t \mathbf{v}_y(t) = \left. \frac{d}{ds} \text{displacement}(s, y) \right|_{s=t}

theorem

The ii-th component of velocity is the time derivative of the ii-th displacement coordinate

Let MM be a rigid body motion in dd-dimensional space. For a material point yy in the body-fixed frame, let r(t,y)\mathbf{r}(t, y) denote its displacement (position in the inertial frame) at time tt, and let v(t,y)\mathbf{v}(t, y) be its velocity. If the trajectory tr(t,y)t \mapsto \mathbf{r}(t, y) is differentiable at time tt, then the ii-th component of the velocity is the time derivative of the ii-th coordinate of the trajectory: vi(t,y)=ddtri(t,y) v_i(t, y) = \frac{d}{dt} r_i(t, y) where i{0,,d1}i \in \{0, \dots, d-1\}.

theorem

The velocity of the material point at the center of mass equals the center-of-mass velocity v(rcm)=V\mathbf{v}(\mathbf{r}_{\text{cm}}) = \mathbf{V}

For a rigid body motion MM in dd-dimensional space, the velocity vrcm\mathbf{v}_{\mathbf{r}_{\text{cm}}} of the material point located at the body's center of mass rcm\mathbf{r}_{\text{cm}} is equal to the center-of-mass velocity V\mathbf{V} of the motion. Mathematically, this is expressed as: vrcm(t)=V(t)\mathbf{v}_{\mathbf{r}_{\text{cm}}}(t) = \mathbf{V}(t) for all times tt, where V(t)\mathbf{V}(t) is the time derivative of the center-of-mass trajectory in the inertial frame.

theorem

In pure translation, vy=V\mathbf{v}_y = \mathbf{V} for any material point yy

For a rigid body motion MM in dd-dimensional space, if the orientation is constant over time—meaning there exists a rotation matrix RR in the special orthogonal group SO(d)SO(d) such that M.orientation(t)=RM.\text{orientation}(t) = R for all tt—then for any material point yy, its velocity vy\mathbf{v}_y is equal to the velocity of the center of mass V\mathbf{V}. Mathematically, if M.orientation(t)=RM.\text{orientation}(t) = R, then vy(t)=V(t)\mathbf{v}_y(t) = \mathbf{V}(t).

theorem

Decomposition of rigid body velocity: v=R˙(yc)+V\mathbf{v} = \dot{\mathbf{R}}(\mathbf{y} - \mathbf{c}) + \mathbf{V}

For a rigid body motion MM in dd-dimensional space and a material point yy in the body-fixed frame, if the orientation R(t)\mathbf{R}(t) and the center-of-mass trajectory are differentiable, then the ii-th component of the velocity of point yy at time tt is given by the decomposition: (vy(t))i=(ddtR(t)(yc))i+(V(t))i (\mathbf{v}_y(t))_i = \left( \frac{d}{dt} \mathbf{R}(t) \cdot (\mathbf{y} - \mathbf{c}) \right)_i + (\mathbf{V}(t))_i where R(t)\mathbf{R}(t) is the rotation matrix (the first component of the orientation) at time tt, c\mathbf{c} is the center of mass of the rigid body in the body-fixed frame, and V(t)\mathbf{V}(t) is the velocity of the center of mass in the inertial frame. In this expression, yc\mathbf{y} - \mathbf{c} represents the position of the point relative to the center of mass in the body frame, and the term ddtR(t)(yc)\frac{d}{dt} \mathbf{R}(t) \cdot (\mathbf{y} - \mathbf{c}) represents the velocity due to the body's rotation.

theorem

Closed-Form Velocity Formula: v=R˙(yc)+V\mathbf{v} = \dot{\mathbf{R}}(\mathbf{y} - \mathbf{c}) + \mathbf{V}

For a rigid body motion MM in dd-dimensional space, the closed-form velocity at time tt for a material point yy (given in the body-fixed frame) is defined as: vclosed(t,y)=(ddtR(t))(yc)+V(t)\mathbf{v}_{\text{closed}}(t, y) = \left( \frac{d}{dt} \mathbf{R}(t) \right) (\mathbf{y} - \mathbf{c}) + \mathbf{V}(t) where: - R(t)\mathbf{R}(t) is the rotation matrix representing the body's orientation at time tt. - c\mathbf{c} is the center of mass of the rigid body in the body-fixed frame. - V(t)\mathbf{V}(t) is the velocity of the center of mass in the inertial frame at time tt. - ddtR(t)(yc)\frac{d}{dt} \mathbf{R}(t) (\mathbf{y} - \mathbf{c}) is the matrix-vector product of the time derivative of the orientation and the position of the point relative to the center of mass.

theorem

vclosed(t,y)=vy(t)\mathbf{v}_{\text{closed}}(t, y) = \mathbf{v}_y(t) for Differentiable Rigid Body Motion

For a rigid body motion MM in dd-dimensional space, let R(t)\mathbf{R}(t) be the rotation matrix representing the body's orientation and xcom(t)\mathbf{x}_{\text{com}}(t) be the trajectory of the center of mass in the inertial frame. If both R\mathbf{R} and xcom\mathbf{x}_{\text{com}} are differentiable at time tt, then for any material point yy in the body-fixed frame, the closed-form velocity vclosed(t,y)\mathbf{v}_{\text{closed}}(t, y) (defined by the rotation and center-of-mass velocity) is equal to the honest point velocity vy(t)\mathbf{v}_y(t) (defined as the time derivative of the point's position in the inertial frame).

theorem

The squared speed yv(t,y)2y \mapsto \|v(t, y)\|^2 of a rigid body point is smooth.

Let MM be a rigid body motion in dd-dimensional space and tt be a fixed time. The function that maps a point yRdy \in \mathbb{R}^d to its squared speed v(t,y),v(t,y)\langle v(t, y), v(t, y) \rangle (where v(t,y)v(t, y) is the velocity of the point at time tt in closed form) is infinitely differentiable (smooth).