Question 9.8: For the satellite of Example 9.2, which is reproduced in Fig......

For the satellite of Example 9.2, which is reproduced in Figure 9.13, the data are as follows: N = 0.1 rad/s and θ˙=0.01 rad/s\dot{\theta}=0.01  rad / s \text {, } in the directions shown: θ=40\theta=40^{\circ} and d0=1.5 m.d_0 = 1.5  m. The length, width, and thickness of the panel are =6 m, \ell = 6  m, w = 2 m, and t = 0.025 m. The uniformly distributed mass of the panel is 50 kg. Find the angular momentum of the panel relative to the center of mass O of the satellite.

9.13
Step-by-Step
The 'Blue Check Mark' means that this solution was answered by an expert.
Learn more on how do we answer questions.

We can treat the panel as a thin parallelepiped. The panel’s xyz axes have their origin at the center of mass G of the panel and are parallel to its three edge directions. According to Figure 9.10(c), the moments of inertia relative to the xyz coordinate system are

Gx=112m(2+t2)=11250(62+0.0252)=150.0 kgm2Gy=112m(w2+t2)=11250(22+0.0252)=16.67 kgm2Gz=112m(w2+2)=11250(22+62)=166.7 kgm2Gxy=Gxz=Gyz=0\begin{aligned}& \ell_{G x}=\frac{1}{12} m\left(\ell^2+t^2\right)=\frac{1}{12} \cdot 50 \cdot\left(6^2+0.025^2\right)=150.0  kg \cdot m ^2 \\\\& \ell_{G y}=\frac{1}{12} m\left(w^2+t^2\right)=\frac{1}{12} \cdot 50 \cdot\left(2^2+0.025^2\right)=16.67  kg \cdot m ^2 \\\\& \ell_{G z}=\frac{1}{12} m\left(w^2+\ell^2\right)=\frac{1}{12} \cdot 50 \cdot\left(2^2+6^2\right)=166.7  kg \cdot m ^2 \\\\& \ell_{G x y}=\ell_{G x z}=\ell_{G y z}=0\end{aligned}                     (a)

In matrix notation,

[IG]=[150.000016.67000166.7] (kgm2)\left[I_G\right]=\left[\begin{array}{ccc}150.0 & 0 & 0 \\0 & 16.67 & 0 \\0 & 0 & 166.7\end{array}\right]  \left( kg \cdot m ^2\right)                              (b)

The unit vectors of the satellite’s xyzx^{\prime} y^{\prime} z^{\prime} system are related to those of the panel’s xyz frame. By inspection,

i^=sinθi^+cosθk^=0.6428i^+0.7660k^j^=j^k^=cosθi^+sinθk^=0.7660i^+0.6428k^\begin{aligned}& \hat{ i }^{\prime}=-\sin \theta \hat{ i }+\cos \theta \hat{ k }=-0.6428 \hat{ i }+0.7660 \hat{ k } \\& \hat{ j }^{\prime}=-\hat{ j } \\& \hat{ k }^{\prime}=\cos \theta \hat{ i }+\sin \theta \hat{ k }=0.7660 \hat{ i }+0.6428 \hat{ k }\end{aligned}                      (c)

The matrix [Q] of the transformation from xyz to xyzx^{\prime} y^{\prime} z^{\prime} comprises the direction cosines of i^,j^, and k^ : \hat{i}^{\prime}, \hat{j}^{\prime} \text {, and } \hat{ k }^{\prime} \text { : }

[Q]=[0.642800.76600100.766000.6428][ Q ]=\left[\begin{array}{ccc}-0.6428 & 0 & 0.7660 \\0 & -1 & 0 \\0.7660 & 0 & 0.6428\end{array}\right]                           (d)

In Example 9.2, we found that the absolute angular velocity of the panel, in the satellite’s xyzx^{\prime} y^{\prime} z^{\prime} frame of reference, is

ω=θ˙j^+Nk^=0.01j^+0.1k^ (rad/s)\omega=-\dot{\theta} \hat{ j }^{\prime}+N \hat{ k }^{\prime}=-0.01 \hat{ j }^{\prime}+0.1 \hat{ k }^{\prime}  ( rad / s )

That is,

{ω}={00.010.1} (rad/s)\left\{ \omega ^{\prime}\right\}=\left\{\begin{array}{c}0 \\-0.01 \\0.1\end{array}\right\}  ( rad / s )                    (e)

To find the absolute angular momentum {HG}\left\{ H _G^{\prime}\right\} in the satellite system requires the use of Eqn (9.39),

{HxHyHz}=[IxIxyIxzIyxIyIyzIzxIzyIz]{ωxωyωz}\left\{\begin{array}{l}H_x \\H_y \\H_z\end{array}\right\}=\left[\begin{array}{ccc}I_x & I_{x y} & I_{x z} \\I_{y x} & I_y & I_{y z} \\I_{z x} & I_{z y} & I_z\end{array}\right]\left\{\begin{array}{c}\omega_x \\\omega_y \\\omega_z\end{array}\right\}                            (9.39a)

{H}=[I]{ω}\{ H \}=[ I ]\{ \omega \}               (9.39b)

{HG}=[IG]{ω}\left\{ H _G^{\prime}\right\}=\left[ I _G^{\prime}\right]\left\{ \omega ^{\prime}\right\}                       (f)

Before doing so, we must transform the components of the moments of the inertia tensor in Eqn (b) from the unprimed system to the primed system, by means of Eqn (9.49),

[I]=[Q][I][Q]T\left[ I ^{\prime}\right]=[ Q ][ I ][ Q ]^T                      (9.49a)

[IxIxyIxzIyxIyIyzIzxIzyIz]=[Q11Q12Q13Q21Q22Q23Q31Q32Q33][IxIxyIxzIyxIyIyzIzxIzyIz][Q11Q21Q31Q12Q22Q32Q13Q32Q33]\left[\begin{array}{ccc}I_{x^{\prime}} & I_{x^{\prime} y^{\prime}} & I_{x^{\prime} z^{\prime}} \\I_{y^{\prime} x^{\prime}} & I_{y^{\prime}} & I_{y^{\prime} z^{\prime}} \\I_{z^{\prime} x^{\prime}} & I_{z^{\prime} y^{\prime}} & I_{z^{\prime}}\end{array}\right]=\left[\begin{array}{lll}Q_{11} & Q_{12} & Q_{13} \\Q_{21} & Q_{22} & Q_{23} \\Q_{31} & Q_{32} & Q_{33}\end{array}\right]\left[\begin{array}{ccc}I_x & I_{x y} & I_{x z} \\I_{y x} & I_y & I_{y z} \\I_{z x} & I_{z y} & I_z\end{array}\right]\left[\begin{array}{lll}Q_{11} & Q_{21} & Q_{31} \\Q_{12} & Q_{22} & Q_{32} \\Q_{13} & Q_{32} & Q_{33}\end{array}\right]                          (9.49b)

 

[IG]=[Q][IG][Q]T=[0.642800.76600100.766000.6428][150.000016.67000166.7][0.642800.76600100.766000.6428]\left[ I _G^{\prime}\right]=[ Q ]\left[ I _G\right][ Q ]^T=\left[\begin{array}{ccc}-0.6428 & 0 & 0.7660 \\0 & -1 & 0 \\0.7660 & 0 & 0.6428\end{array}\right]\left[\begin{array}{ccc}150.0 & 0 & 0 \\0 & 16.67 & 0 \\0 & 0 & 166.7\end{array}\right]\left[\begin{array}{ccc}-0.6428 & 0 & 0.7660 \\0 & -1 & 0 \\0.7660 & 0 & 0.6428\end{array}\right]

so that

[IG]=[159.808.205016.6708.2050156.9] (kgm2)\left[I_G^{\prime}\right]=\left[\begin{array}{ccc}159.8 & 0 & 8.205 \\0 & 16.67 & 0 \\8.205 & 0 & 156.9\end{array}\right]  \left( kg \cdot m ^2\right)                        (g)

Then Eqn (f) yields

{HG}=[159.808.205016.6708.2050156.9]{00.010.1}={0.82050.166715.69} (kgm2/s)\left\{ H _G^{\prime}\right\}=\left[\begin{array}{ccc}159.8 & 0 & 8.205 \\0 & 16.67 & 0 \\8.205 & 0 & 156.9\end{array}\right]\left\{\begin{array}{c}0 \\-0.01 \\0.1\end{array}\right\}=\left\{\begin{array}{c}0.8205 \\-0.1667 \\15.69\end{array}\right\}  \left( kg \cdot m ^2 / s \right)

or, in vector notation,

HG=0.8205i^0.1667j^+15.69k^ (kgm2/s)H _G=0.8205 \hat{i} ^{\prime}-0.1667 \hat{ j }^{\prime}+15.69 \hat{ k }^{\prime}  \left( kg \cdot m ^2 / s \right)                     (h)

This is the absolute angular momentum of the panel about its own center of mass G, and it is used in Eqn (9.27) to calculate the angular momentum HOrel H _{O_{\text {rel }}} relative to the satellite’s center of mass O,

HPrel=HG+rG/P×mvG/PH _{P_{ rel }}= H _G+ r _{G / P} \times m v _{G / P}                       (9.27)

HOrel =HG+rG/O×mvG/OH _{O_{\text {rel }}}= H _G+ r _{G / O} \times m v _{G / O}                       (i)

rG/Or _{G / O} is the position vector from O to G,

rG/O=(dO+2)j^=(1.5+62)j^=4.5j^ (m)r _{G / O}=\left(d_O+\cfrac{\ell}{2}\right) \hat{ j }^{\prime}=\left(1.5+\cfrac{6}{2}\right) \hat{ j }^{\prime}=4.5 \hat{ j }^{\prime}  ( m )                        (j)

The velocity of G relative to O, vG/Ov _{G / O}, is found from Eqn (9.2),

vB=vA+ω×RB/Av _B= v _A+ \omega \times R _{B / A}                    (9.2)

vG/O=ωsatellite ×rG/O=Nk^×rG/O=0.1k^×4.5j^=0.45i^  (m/s)v _{G / O}= \omega _{\text {satellite }} \times r _{G / O}=N \hat{ k }^{\prime} \times r _{G / O}=0.1 \hat{ k }^{\prime} \times 4.5 \hat{ j }^{\prime}=-0.45 \hat{ i }^{\prime}   ( m / s )          (k)

Substituting Eqns (h), (j) and (k) into Eqn (i) finally yields

HOrel =(0.8205i^0.1667j^+15.69k^)+4.5j×[50(0.45i^)]=0.8205i^0.1667j^+116.9k^ (kgm2/s)\begin{aligned}H _{O_{\text {rel }}} & =\left(0.8205 \hat{ i }^{\prime}-0.1667 \hat{ j }^{\prime}+15.69 \hat{ k }^{\prime}\right)+4.5 j ^{\prime} \times\left[50\left(-0.45 \hat{i}^{\prime}\right)\right] \\\\& = \boxed{ 0.8205 \hat{i}^{\prime}-0.1667 \hat{ j }^{\prime}+116.9 \hat{ k }^{\prime}  \left( kg \cdot m ^2 / s \right)}\end{aligned}                    (l)

Note that we were unable to use Eqn (9.21) to find the absolute angular momentum HOH_O because that requires knowing the absolute velocity vG,v_G, which in turn depends on the absolute velocity of O, which was not provided.

HP=HG+rG/P×mvGH _P= H _G+ r _{G / P} \times m v _G                        (9.21)

9.10c

Related Answered Questions

Question: 9.12

Verified Answer:

Since the comoving frame is rigidly attached to th...
Question: 9.10

Verified Answer:

From Example 9.5, \left[ I _A\right]=\left[...