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=121m(ℓ2+t2)=121⋅50⋅(62+0.0252)=150.0 kg⋅m2ℓGy=121m(w2+t2)=121⋅50⋅(22+0.0252)=16.67 kg⋅m2ℓGz=121m(w2+ℓ2)=121⋅50⋅(22+62)=166.7 kg⋅m2ℓGxy=ℓGxz=ℓGyz=0 (a)
In matrix notation,
[IG]=⎣⎢⎡150.000016.67000166.7⎦⎥⎤ (kg⋅m2) (b)
The unit vectors of the satellite’s x′y′z′ 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^ (c)
The matrix [Q] of the transformation from xyz to x′y′z′ comprises the direction cosines of i^′,j^′, and k^′ :
[Q]=⎣⎢⎡−0.642800.76600−100.766000.6428⎦⎥⎤ (d)
In Example 9.2, we found that the absolute angular velocity of the panel, in the satellite’s x′y′z′ frame of reference, is
ω=−θ˙j^′+Nk^′=−0.01j^′+0.1k^′ (rad/s)
That is,
{ω′}=⎩⎪⎨⎪⎧0−0.010.1⎭⎪⎬⎪⎫ (rad/s) (e)
To find the absolute angular momentum {HG′} in the satellite system requires the use of Eqn (9.39),
⎩⎪⎨⎪⎧HxHyHz⎭⎪⎬⎪⎫=⎣⎢⎡IxIyxIzxIxyIyIzyIxzIyzIz⎦⎥⎤⎩⎪⎨⎪⎧ωxωyωz⎭⎪⎬⎪⎫ (9.39a)
{H}=[I]{ω} (9.39b)
{HG′}=[IG′]{ω′} (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 (9.49a)
⎣⎢⎡Ix′Iy′x′Iz′x′Ix′y′Iy′Iz′y′Ix′z′Iy′z′Iz′⎦⎥⎤=⎣⎢⎡Q11Q21Q31Q12Q22Q32Q13Q23Q33⎦⎥⎤⎣⎢⎡IxIyxIzxIxyIyIzyIxzIyzIz⎦⎥⎤⎣⎢⎡Q11Q12Q13Q21Q22Q32Q31Q32Q33⎦⎥⎤ (9.49b)
[IG′]=[Q][IG][Q]T=⎣⎢⎡−0.642800.76600−100.766000.6428⎦⎥⎤⎣⎢⎡150.000016.67000166.7⎦⎥⎤⎣⎢⎡−0.642800.76600−100.766000.6428⎦⎥⎤
so that
[IG′]=⎣⎢⎡159.808.205016.6708.2050156.9⎦⎥⎤ (kg⋅m2) (g)
Then Eqn (f) yields
{HG′}=⎣⎢⎡159.808.205016.6708.2050156.9⎦⎥⎤⎩⎪⎨⎪⎧0−0.010.1⎭⎪⎬⎪⎫=⎩⎪⎨⎪⎧0.8205−0.166715.69⎭⎪⎬⎪⎫ (kg⋅m2/s)
or, in vector notation,
HG=0.8205i^′−0.1667j^′+15.69k^′ (kg⋅m2/s) (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 relative to the satellite’s center of mass O,
HPrel=HG+rG/P×mvG/P (9.27)
HOrel =HG+rG/O×mvG/O (i)
rG/O is the position vector from O to G,
rG/O=(dO+2ℓ)j^′=(1.5+26)j^′=4.5j^′ (m) (j)
The velocity of G relative to O, vG/O, is found from Eqn (9.2),
vB=vA+ω×RB/A (9.2)
vG/O=ωsatellite ×rG/O=Nk^′×rG/O=0.1k^′×4.5j^′=−0.45i^′ (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^′ (kg⋅m2/s) (l)
Note that we were unable to use Eqn (9.21) to find the absolute angular momentum HO because that requires knowing the absolute velocity vG, which in turn depends on the absolute velocity of O, which was not provided.
HP=HG+rG/P×mvG (9.21)