Skip to main content

Section 12.7 Examples

Example 12.7.1.

Consider the rectangular plate 1 of Example 12.2.4 Using the same assumptions regarding the forces acting on body 1, apply the KET to find a scalar equation. How does this equation relate to the FTD?
Figure 12.7.2.
Solution.
We found the total external power \(\Pow_{\bar{1}\to 1/0}\) in Example 12.2.4.
\begin{equation*} \Pow_{\bar{1}\to 1/0} = - \tfrac{4}{3} \mu a^3 b \dte^2+ \cC \dte \end{equation*}
Furthermore the kinetic energy of body 1 is given by
\begin{equation*} \kin_{1/0} = \half \bom_{1/0} \cdot \cI_O (\bom_{1/0}) = \tfrac{2}{3} ma^2 \dte^2 \end{equation*}
Application of the KET to the motion of body 1 gives
\begin{equation*} \tfrac{d}{dt} \Big( \tfrac{2}{3} ma^2 \dte^2 \Big) = - \tfrac{4}{3} \mu a^3 b \dte^2 +\cC \dte \qquad (1) \end{equation*}
To show the relationship between this result and the FTD, recall that the KET can be written in the following form:
\begin{equation*} \{ \cD_{1/0} \} \cdot \{\cV_{1/0} \} = \{ \cA_{\bar{1}\to 1}\} \cdot \{\cV_{1/0} \} \end{equation*}
which gives
\begin{equation*} \dte \bz_0 \cdot \bD_{O, 1/0} = \dte \bz_0 \cdot \bM_{O, \bar{1}\to 1} \end{equation*}
In other words, equation (1) is nothing but the \(\bz_0\)-component of the dynamic moment equation about \(O\) premultiplied by \(\dte\text{.}\)

Example 12.7.3.

A spinning top 1\((O,\bx_1, \by_1,\bz_1)\) of mass center \(G\) and mass \(m\) is in motion in a referential 0\((O,\bx_0, \by_0,\bz_0)\) in such a way that its tip \(O\) remains stationary relative to 0. The connection between 0 and 1 is equivalent to a frictionless spherical joint. The upward vertical is directed along \(\bz_0\text{.}\) The orientation of body 1 relative to 0 is defined in terms of the three angles \((\psi,\te, \phi)\) according the following sequence of rotations:
\begin{equation*} (\bx_0, \by_0,\bz_0) \xrightarrow{\cR_{\psi , \bz_0}} (\bu, \bv,\bz_0) \xrightarrow{\cR_{\te , \bu}} (\bu, \bw,\bz_1) \xrightarrow{\cR_{\phi , \bz_1}} (\bx_1, \by_1,\bz_1) \end{equation*}
The top is axisymmetric about axis \((O,\bz_1)\text{.}\) The position of its mass center is defined by \(\br_{OG}= l \bz_1\text{.}\) Its inertia operator about \(O\) is entirely defined by its axial moment of inertia \(C= \bz_1 \cdot \cI_O (\bz_1)\) and its transverse moment of inertia \(A= \bx_1 \cdot \cI_O (\bx_1)= \by_1 \cdot \cI_O (\by_1)\text{.}\) The top is under the sole effect of gravity and contact forces, neglecting rolling and pivoting friction.
Figure 12.7.4.
Apply the Kinetic Energy Theorem to body 1 and find a first-integral of motion.
Solution.
Step 1: kinematics. With \(\bom_{1/0} = \dpsi \bz_0 + \dte \bu + \dphi \bz_1\) and \(\vel_{O\in 1/0} =\bze\text{,}\) we find the expression of the kinematic screw of the top
\begin{equation*} \{ {\cal V} _{1/0} \} = \begin{Bmatrix} \dpsi \bz_0 + \dte \bu + \dphi \bz_1 \\ \bze \end{Bmatrix} _O \end{equation*}
Step 2: kinetic energy. We find the KE of the top as a body rotating about fixed point \(O\) by representing inertia operator \(\cI_O\) on basis \(b (\bu,\bw,\bz_1)\) (using \(\bom_{1/0} \cdot\bz_1 = \dpsi \cos\te +\dphi\text{,}\) \(\bom_{1/0}\cdot \bw= \dpsi \sin\te\text{,}\) and \(\bom_{1/0}\cdot \bu= \dte\)):
\begin{align*} 2\kin_{1/0} \amp = \bom_{1/0} \cdot \cI_O (\bom_{1/0})\\ \amp = \bom_{1/0} \cdot \begin{bmatrix} A \amp 0 \amp 0 \\ 0 \amp A \amp 0\\ 0 \amp 0 \amp C \end{bmatrix}_{b} \begin{bmatrix} \dte \\ \dpsi \sin\te\\ \dpsi \cos\te +\dphi \end{bmatrix}_{b}\\ \amp = A (\dpsi^2 \sin^2\te + \dte^2) + C (\dpsi \cos\te +\dphi)^2 \end{align*}
Step 3: action screw. We find the contributions to action screw \(\{\cA_{\bar{1}\to 1}\}\)
\begin{equation*} \{\cA_{\bar{1}\to 1}\}= \begin{Bmatrix} -m g \bz_0 \\\bze \end{Bmatrix}_G + \begin{Bmatrix} \bR_O \\ \bze \end{Bmatrix}_O = \begin{Bmatrix} \bR_O -m g \bz_0 \\ mg l \sin\te \bu \end{Bmatrix}_O \end{equation*}
Step 4: We find the power \(\Pow_{\bar{1}\to 1/0}\text{:}\)
\begin{equation*} \Pow_{\bar{1}\to 1/0}= mg l \sin\te \bu \cdot \bom_{1/0} = mg l \dte \sin\te = - \frac{d}{dt}( mg l \cos\te) \end{equation*}
Step 5: We apply the KET \(\tfrac{d}{dt} \kin_{1/0} = \Pow_{\bar{1}\to 1/0}\) to find
\begin{equation*} \frac{d}{dt}\Big( A (\dpsi^2 \sin^2\te + \dte^2) + C (\dpsi \cos\te +\dphi)^2 \Big) = - \frac{d}{dt}( 2mg l \cos\te) \end{equation*}
which yields a first-integral of motion (conservation of mechanical energy):
\begin{equation*} A (\dpsi^2 \sin^2\te + \dte^2) + C (\dpsi \cos\te +\dphi)^2 + 2mg l \cos\te = Cst \end{equation*}
Note that we could have used the potential energy \(\pot_{earth \to 1 /0} = mg l \cos\te\text{.}\)

Example 12.7.5.

Figure 12.7.6 shows a truncated right-circular cone 2 of half-angle \(\theta\text{,}\) uniformly-distributed mass \(m\text{,}\) and mass center \(G\text{.}\) It is in motion relative to a horizontal support 1 in such a way that (i) one of its generatrix lines is in contact with 1, and (ii) its axis always passes through fixed point \(O\) of 1, (iii) it does not slip relative to 1.
Figure 12.7.6.
Its motion relative to 1 is entirely specified by the angles \(\psi\) and \(\phi\text{.}\) Denote by \(I_1\) the moment of inertia of 1 about axis \(OZ\text{,}\) by \(A_2\) and \(C_2\) the transverse and axial moments of inertia which characterize the inertia operator \(\cI_{O,2}\) of body 2.
  1. Assume that support 1 is stationary. Apply the KET to body 2 relative to body 1.
  2. Now assume that support 1 is in rotation about vertical axis \(OZ\) at angular velocity \(\Omega\) under the action of a motor of couple \(\cC\) mounted between 1 and the supporting frame 0. Apply the KET to 2 relative to 0, and then to the system \(\Sigma=\{1, 2\}\) relative to 0. Assume that the connection between 1 and 0 is an ideal pivot of axis \(OZ\text{.}\)
Solution.
a. If body 1 is stationary, we can apply the KET to body 2 relative to body 1:
\begin{equation*} \tfrac{d}{dt}\kin_{2/1} = \Pow_{\bar{2}\to 2/1} \end{equation*}
To find power \(\Pow_{\bar{2}\to 2/1}\) we first note that \(\{\cA_{\bar{2}\to 2} \} = \{ \cA^g_{\text{Earth} \to 2}\} + \{\cA^c_{1\to 2} \}\text{:}\)
\begin{equation*} \Pow_{\bar{2}\to 2/1} = \underbrace{-m g \bz_1 \cdot \vel_{G/1}}_{0} + \underbrace{\int_{l_A}^{l_B} (N\bz_1 + F_u \bu + F_v \bv ) \cdot \vel_{Q\in 2/1} \, dl }_{0} = 0 \end{equation*}
where we have accounted for contact forces distributed along the contact line \(AB\) (with \(\bv = \bz_1 \times \bu\)), and the no-slip kinematic condition \(\vel_{Q\in 2/1} = \bze\text{.}\)
It results that kinetic energy \(\kin_{2/1}\) remains constant during the motion (thus yielding a first-integral of motion). To find \(\kin_{2/1}\text{,}\) we first find angular velocity \(\bom_{2/1}= \dpsi \bz_1 + \dphi \bz_2\text{.}\) Imposing the no-slip condition at \(Q\) gives the condition \(\dpsi = - \dphi \sin\te\) and \(\bom_{2/1}= -\dpsi\cot\te \bu\text{.}\)
The kinetic energy \(\kin_{2/1}\) is then given by
\begin{equation*} \kin_{2/1} = \half \bom_{2/1}\cdot \cI_{O,2} (\bom_{2/1}) = \half (A_2 + C_2 \cot^2\te) \dpsi^2 \cos^2\te = Cst \end{equation*}
This implies that \(\dot{\psi} =\) constant (or \(\dot{\phi} =\) constant).
b. If body 1 is in rotation, we cannot apply the KET to 2 relative to 1 which is no longer Newtonian. We find power \(\Pow_{\bar{2}\to 2/0}\) as in question a):
\begin{align*} \tfrac{d}{dt}\kin_{2/0} \amp =\Pow_{\bar{2}\to 2/0} = \underbrace{-m g \bz_1 \cdot \vel_{G/0}}_{0} + \int_{l_A}^{l_B} (N\bz_1 + F_u \bu + F_v \bv ) \cdot \vel_{Q\in 2/0} \, dl \\\\ \amp = \int_{l_A}^{l_B} F_v l \Omega \, dl \neq 0 \end{align*}
where we have used \(\vel_{Q\in 2/0} = \vel_{Q\in 2/1} +\vel_{Q\in 1/0}=\vel_{Q\in 1/0}= \Omega l \bv\) (denoting \(l= |OQ|\)). We end up with an integral expression whose value cannot be explicitly written in terms of the motion variables \((\Omega, \psi, \phi)\text{.}\)
The situation is different if we apply the KET to the system \(\Sigma=\{1,2\}\) relative to 0:
\begin{equation*} \tfrac{d}{dt}(\kin_{1/0}+\kin_{2/0}) = \underbrace{\Pow_{\text{Earth} \to \Sigma /0}^g}_{0}+ \underbrace{\Pow^c_{0 \to 1 /0}}_{0}+ \underbrace{\Pow_{\text{Motor} \to 1 /0}}_{\cC \Omega} + \underbrace{\Pow^c_{1\leftrightarrow 2}}_{0} = \cC \Omega \end{equation*}
Note here that we have to account for the internal interaction between 1 and 2 in the term \(\Pow^c_{1\leftrightarrow 2}\) which is in fact identical to \(\Pow^c_{1\to 2/1} =0\text{.}\)
The kinetic energy of the system is given by
\begin{align*} \kin_{\Si/0} \amp = \kin_{1/0}+\kin_{2/0}\\ \amp = \half I_1 \Omega^2 + \half A_2 (\dpsi+\Om)^2 \cos^2\te + \half C_2 [(\dpsi+\Om) \sin\te + \dphi]^2 \end{align*}
with \(\dphi= -\dpsi /\sin\te\text{.}\)