128
Po podstawieniu powyższych zależności do (4.33) otrzymamy
Ixx2 + Iyy2 + I-z2 - 2Dxy ■ xy - 2Dyl ■ yz - 2Dxz ■ xz = 1 (4.37)
Z teorii form kwadratowych wiadomo, że rówmanie (4.37) przedstawia elipsoidę trójosiową, jeśli tylko Ix, ly, I. są dodatnie. Elipsoida ta nazywa się elipsoidą bezwładności. Jest ona miejscem geometrycznym punktów, których odległość od jej środka jest odwrotnie proporcjonalna do pierwiastka kwadratowego z momentu bezwładności względem osi przechodzącej przez dany punkt i środek elipsoidy.
Trójosiową elipsoida ma układ trzech średnic sprzężonych i wzajemnie do siebie prostopadłych, z których jedna jest największa ze wszystkich średnic, a jedna najmniejsza. W szczególnych przypadkach niektóre średnice mogą być jednakowe, tworząc wówczas elipsoidę obrotową. Osie, na których leżą średnice sprzężone są osiami symetrii dla elipsoidy i zwą się osiami głównymi. Odpowiadające im momenty bezwładności są tzw. momentami głównymi bezwładności. Jeśli punkt O pokrywa się ze środkiem masy układu punktów materialnych, to mówimy, że odpowiadająca mu elipsoida jest centralną elipsoidą bezwładności, główne osie są centralnymi, a odpowiadające tym osiom momenty bezwładności są głównymi centralnymi momentami bezwładności.
Znajdźmy ekstrema funkcji /„ odpowiadające głównym momentom bezwładności. W tym cdi u zastosujemy metodę mnożników Lagrange’a.
Badaną funkcją jest J, = /.(a, /?, y). Dodatkowym warunkiem (warunków takich musi być jednak mniej niż niewiadomych) jest związek
(4.38)
/(a, /J, y) = cos2a + cos2fi + cos2y -1 = 0.
Jest to tzw. równanie więzów.
Warunkiem koniecznym istnienia ekstremum funkcji 7, przy ograniczeniu (4.38) jest zerowanie się wyrażeń
(4.39)
8<P _ 8<P _ 8<P da dp dy
gdzie
<P(a, fi, y) = 7,(a, fi, y) - A/(a, P, y). (4.40)
Parametr A nazywa się mnożnikiem Lagrange’a i jest wielkością stalą. Parametr ten, a także niewiadome a„ /?„ y„ przy których 7, osiąga 'ekstremum, wyznaczamy z równań (4.38) i (4.39). Zatem
<P = Ix cos2a + Iy cos2p + I. cos2y - 2 Dxy cos a cos P +
- 27)^ cos P cos y - 2Dzxcosy cos a +
+ A(cos2a + cos2/? + cos2y - 1).
Ponieważ w wyrażeniu określającym zależność <P kąty wchodzą wyłącznie poprzez ich kosinusy, więc badajmy ekstrema nie względem kątów, a względem ich kosinusów. Po wstawieniu tej wartości do (4.39) otrzymujemy
dtp
8 cos a 8<P
8 cos P 8<P
dc osy
czyli
= 27xcosa - 2DxycosP - 2D.xcosy — 2 A cos a = 0,
= 27^ cos p — 2Dyz cos y — 2Dxy cos ot — 2 A cos P = 0,
= 27. cos y - 2DjXcosa - 2DyzcosP - 2 A cos y = 0,
(7X - A) cos a - cos P - Dx. cos y = 0,
-D^ cos.a + (Iy - A) cos P - Dyz cos y = 0, (4.41)
-D:xcosa - D:ycosP + (I, - A)cosy = 0.
Wyznacznik główny tego równania musi się zerować, by równanie miało rozwiązanie niezerowe (zerowym rozwiązaniem jest cosa = cosP = cosy = 0). Wyznacznik ten ma postać
ł Ix Dxy ~DXZ
7A Xy, Iy A, Dys -Dxz, -Dyz, lz -A
= 0