Rotational Dynamics – Understanding The Angular Velocity

angular velocityreference framesrigid-body-dynamicsrotational-dynamicsrotational-kinematics

Is the angular velocity of a rigid body about any point the same as that about the axis of rotation. Also, can we even define angular terms (Angular Velocity, Angular Acceleration, etc) about any axis, other than the axis of rotation?

Best Answer

Given some central point $\vec c_1$, the velocity of this point $\vec v_{c_1}$ with respect to some frame of reference, and an angular velocity $\vec\omega$, the velocity of some other point $\vec x$ on a rigid body is $\vec v_x = \vec v_{c_1} + \vec\omega \times (\vec x - \vec c_1)$.

What if some other point $\vec c_2$ is chosen as the central point? The expression for the velocity of the point $\vec x$ becomes $\vec v_x = \vec v_{c_2} + \vec\omega \times (\vec x - \vec c_2)$, where $\vec v_{c_2} = \vec v_{c_1} + \vec\omega \times (\vec c_2 - \vec c_1)$. The angular velocity does not change. It is the same regardless of which point one chooses as the central point. Angular velocity is a free vector.


Update, because the above apparently isn't satisfactory to some

A rigid body is an object for which a frame of reference exists such that the location of every point on the rigid body is constant from the perspective of this frame. In other words, $(\frac {d\boldsymbol x}{dt})_F = 0$ for every point $x$ in the rigid body, with the derivative taken from the perspective of the fixed frame $F$.

Suppose at some time $t$ you know the location $\boldsymbol X(\boldsymbol c,t)$ of some fixed point $\boldsymbol c$ on the rigid body in some other frame of reference $I$. The location of the point $\boldsymbol x$ in this other frame is related to the location of the point $\boldsymbol c$ via $$\boldsymbol X(\boldsymbol x,t) = \boldsymbol X(\boldsymbol c,t) + \boldsymbol{\mathrm R}_{F\to I}(t)(\boldsymbol x - \boldsymbol c) \tag{1}$$ where $\boldsymbol{\mathrm R}_{F\to I}(t)$ is the rotation matrix that transforms coordinates from frame $F$ to frame $I$.

Suppose you know some other point $\boldsymbol c'$, also fixed with respect to the rigid body. The location of the point $\boldsymbol c'$ in the non-fixed frame is $$\boldsymbol X(\boldsymbol c',t) = \boldsymbol X(\boldsymbol c,t) + \boldsymbol{\mathrm R}_{F\to I}(t)(\boldsymbol c' - \boldsymbol c)\tag2$$ The location of the point $\boldsymbol x$ in the non-fixed frame can be re-expressed to be in terms of $\boldsymbol c'$: $$\boldsymbol X(\boldsymbol x,t) = \boldsymbol X(\boldsymbol c',t) + \boldsymbol{\mathrm R}_{F\to I}(t)(\boldsymbol x - \boldsymbol c') \tag{3}$$ Differentiating equations (1) and (3) with respect to time yields $$\dot {\boldsymbol X}(\boldsymbol x,t) = \dot {\boldsymbol X}(\boldsymbol c,t) + \dot{\boldsymbol{\mathrm R}}_{F\to I}(t)(\boldsymbol x - \boldsymbol c) = \dot{\boldsymbol X}(\boldsymbol c',t) + \dot {\boldsymbol{\mathrm R}}_{F\to I}(t)(\boldsymbol x - \boldsymbol c') \tag4$$ The time derivative of a transformation matrix $\boldsymbol{\mathrm R}(t)$ from one Cartesian frame to another can be written as a product of that matrix and a skew symmetric matrix: $$\dot{\boldsymbol{\mathrm R}}(t) = \boldsymbol{\mathrm R}(t) \boldsymbol{\mathrm S}(t)$$ This is valid for a Euclidean space of any dimensionality. In three dimensional space (and three dimensional space only), such a skew symmetric matrix can be mapped to and from a three dimensional pseudo vector: $$\dot{\boldsymbol{\mathrm R}}(t) = \boldsymbol{\mathrm R}(t) \operatorname{Sk}(\boldsymbol\omega(t))$$ Rewriting equation (4) in terms of this pseudo vector, $$\dot {\boldsymbol X}(\boldsymbol x,t) = \dot {\boldsymbol X}(\boldsymbol c,t) + \boldsymbol{\mathrm R}_{F\to I}(t)\left(\boldsymbol\omega\times(\boldsymbol x - \boldsymbol c)\right) = \dot{\boldsymbol X}(\boldsymbol c',t) + \boldsymbol{\mathrm R}_{F\to I}(t)\left(\boldsymbol\omega\times(\boldsymbol x - \boldsymbol c')\right)\tag5$$ Without loss of generality, one can choose a non-fixed frame that is instantaneously co-aligned with the fixed frame at time $t$. With this, equation (5) reduces to $$\dot {\boldsymbol X}(\boldsymbol x,t) = \dot {\boldsymbol X}(\boldsymbol c,t) + \boldsymbol\omega\times(\boldsymbol x - \boldsymbol c) = \dot{\boldsymbol X}(\boldsymbol c',t) + \boldsymbol\omega\times(\boldsymbol x - \boldsymbol c')\tag6$$ This is of course identical to what I originally wrote.

That angular velocity expressed as a skew symmetric matrix is the same regardless of which point one chooses as the origin applies in all Euclidean spaces in which time is the independent parameter of motion. Angular velocity is related to the time derivative of the transformation matrix, and changing origins doesn't change the transformation matrix in an affine transformation (e.g., equations (1) to (3) are affine transformations).

That angular velocity expressed as a pseudo vector is the same regardless of which point one chooses as the origin applies only in three dimensional Euclidean spaces in which time is the independent parameter of motion. That specialization is rather important because we apparently live in a universe that locally appears to be a three dimensional Euclidean space with time as the independent parameter of motion. In other words, we live in a universe where Newtonian mechanics is locally valid. This is the context in which this question was asked and in which I wrote this answer.