151. Момент инерции относительно точки. Моментом инерции данной системы материальных частиц относительно данной точки, или полюса, называется сумма произведений масс частиц на квадраты их расстояний от взятого полюса. Пусть данный полюс $O$ служит началом координат; радиус-вектор частицы с массой $m_{v}$ и координатами $\boldsymbol{x}_{v}, \boldsymbol{y}_{v}, \boldsymbol{z}_{v}$ обозначим $r_{v}$; тогда, если число частиц равно $n$, момент инерции $J_{o}$ данной системы относительно полюса $O$ представится так:
\[
J_{O}=\sum_{v=1}^{n} m_{v} r_{v}^{2}=\sum_{v=1}^{n} m_{v}\left(x_{v}^{2}+y_{v}^{2}+z_{v}^{2}\right) .
\]

Для масс, распределённых сплошным образом, сумма в этом выражении заменится интегралом, распространённым на объём, занимаемый массами:
\[
J_{o}=\iiint r^{2} d m=\iiint \sigma r^{2} d v,
\]

где $\sigma$-плотность, а $d v$-элемент объёма. Из определения момента инерции относительно точки следует, что он представляет собой величину неотрицательную; он может обратиться в нуль только в случае единственной частицы, совпадающей с полюсом.

Найдем зависимость между моментом инерции $J_{0}$ системы относительно произвольного полюса $O$ и её моментом $J_{c}$ относительно центра масс. Рассмотрим произвольную частицу системы, имеющую массу $m_{y}$; обозначим ее радиусы-векторы, проведынные из полюсов $O$ и $C$, соответственно $r_{v}$ и $\rho_{v}$; вектор $\overline{O C}$ назовём $r_{C}$; тогда
\[
r_{v}=r_{s}+\bar{\rho}_{v},
\]

и, следовательно,
\[
J_{O}=\sum_{v} m_{v} r_{v}^{2}=\sum_{v} m_{v}\left(r_{C}+\bar{\rho}_{v}\right)^{2}=\left(\sum_{v} m_{v}\right) r_{C}^{2}+2 r_{C} \cdot \sum_{v} m_{v} \bar{\rho}_{v}+\sum_{v} m_{v} \rho_{v}^{2} .
\]

Здесь $\sum m_{v}$ есть масса всей системы частиц; мы её обозначим $M$; сумма $\sum_{\vee} m_{v} \overline{P_{v}}$ представляет собой статический момент системы относительно центра масс и потому равна нулю (§ 144); наконец, последнее слагаемое равно моменту инерции $J_{C}$ системы относительно центра масс. Игак, мы получили:
\[
J_{o}=M r_{C}^{2}+J_{C},
\]
т. е. момент инерции системы относительно любого полюса равняется ее моменту инериии относительно центра масс, сложенному с моментом инериии самого центра масс относительно взятого центра, если примем, что масса центра масс равна массе всей системы.

Выражению для $J_{C}$ легко можем дать иной вид. Для этого рассмотрим тождество
\[
\sum_{\mu=1}^{n} m_{\mu} \cdot \sum_{v=1}^{n} m_{
u} \rho_{
u}^{2}=\left(\sum_{
u=1}^{n} m_{
u} \bar{\rho}_{v}\right)^{2}+\sum_{v=1}^{n-1} \sum_{\mu=v+1}^{n} m_{v} m_{\mu}\left(\bar{\rho}_{
u}-\bar{\rho}_{\mu}\right)^{2},
\]

координаты вернин кәкого-либо из них равны $x_{a}, y_{a}, x_{b}, y_{b}, x_{c}, y_{c}$, то площадь ero равна половине абсолютной величины определителя:
\[
\left|\begin{array}{lll}
x_{a} & y_{a} & 1 \\
x_{b} & y_{b} & 1 \\
x_{c} & y_{c} & 1
\end{array}\right| .
\]

Если сечение сделано на высоте $z$, то из выражений (25.14) ясно, что плошадь треугольника, а следовательно, и площадь $P$ многоугольника, образованиого в сечении, представится целой квадратичной функцией от $z$ :
\[
P=A+B z+C z^{3} \text {. }
\]

Таким образом, теорема примера 65 применима и в рассматриваемом случае, если только формула (25.15) годится для сечения на любой высоте, т. е. если боковые рёбра идут, не прерываясь от верхнего основания до нижнего.

Приложив полученный результат к полной пирамиде, находим, что е» центр масс отстоит от основания на четверти высоты пирамиды.
150. Теоремы Паппуса-Гюльдена. Непосредственно пользуясь выражениями для координат центра масс данной крнвой или данной площади, легко убеждаемся в справедливости следуюиих двух теорем Папнуса-Гюльдена(Pappus-Guldin):
1) Если плоскую кривую повернуть на некоторый угол вокруг какой-либо оси, лежащей в плоскости кривой и ее не пересекающей, то величина площади части поверхности вращения, описанной кривою, равняется длине Фиг. 97 .
Фиг. 98. кривой, умноженной на длину пути ее центра масс $C$.
2) Если какую-либо замкнутую плоскую фигуру повернуть на некоторый угол вокруг оси, лежалей с ней в одной плоскости, но не встречающей ее контура, то объем части тела вращения, описанного фигурой, равняется произведению величины площади фигуры на длину пути ее центра масс $C$.

Действительно, площадь полоски, описанной элементом $d s$ кривой $A B$ при повороте ее на угол $\Delta$ ‘ , равна $x \Delta p \cdot d s$ (фиг. 97); следовательно, вся площадь $\mathcal{S}$, описанная кривой $A B$, если $\overparen{A B}=s$, будет иметь выражение
\[
S=\Delta \varphi \int x d s=\Delta \varphi x_{C} \cdot s .
\]

Для доказательства второй теоремы рассмотрим элемент abde плоской фигуры $A B D E$ (фиг. 98): если его площадь обозначить $d S$, то объем, им описанный при повороте фигуры на угол $\Delta^{\prime} \rho$, будет $x \Delta^{\prime} \cdot d S$; следовательно, весь объём $v$, описанный фигурой $A B D E$, если площадь ее обозначить через $S$, будет равен
\[
v=\Delta \varphi \iint x d S=\Delta \varphi x_{c} \cdot S .
\]

где последняя сумма распространяется на все члены, в которых $
u<\mu$. Если $\bar{\rho}$, попрежнему есть радиус-вектор $v$-ой частицы, шроведённый из центра масс, то из этого равенства мы получаем:
\[
J_{C}=\frac{1}{M} \sum_{
u=1}^{n-1} \sum_{\mu=v+1}^{n} m_{v} m_{\mu} \rho_{
u \mu}^{2},
\]

где $P_{
u \mu}$ — расстояние между массами $m_{
u}$ и $m_{\mu}$.

152. Момент инерции относительно оси. Моментоминерции данной системы материальных частиц относительно некоторой ос и называется сумма произведений масс частиц на квадраты их расстояний от взятой оси, т. е. если момент инерции относительно оси $и$ мы обозначим через $J_{u t}$, а расстояние частицы с массой $m_{v}$ от этой оси через $\rho_{
u}$, то
\[
J_{u u}=\sum_{v} m_{v} \rho_{v}^{2}
\]

Из этого выражения видим, что момент инерции относительно оси есть величина неотрицательная; он может обратиться в нуль лишь в том случае, когда все материальные частицы расположены на прямой, совпадающей со взятой осью. Если массы расположены непрерывно, то в выражении момента инерции вместо суммы получим интеграл
\[
J_{t u}=\iiint \rho^{2} d m=\iiint \sigma \rho^{2} d v,
\]

Иной раз приходится рассматривать моменты инерции масс, расположенных по поверхностям или по линиям (ср. § 147). Конечно, как и в главе о центре масс, такое расположение масс следует рассматривать как вспомогательное геометрическое построение или принимать, что одно или два измерения объёма, заполненного массами, настолько малы, что мы считаем себя в праве не принимать их в расчёт. Для масс, распределённых по поверхности, тройной интеграл (26.3) заменится двойным:
\[
J_{u u}=\iint_{0} \rho^{2} d \dot{S},
\]

где $\sigma_{1}$ — поверхностная плотность, а $d S$ — элемент поверхности. Для масс линейных имеем обыкновенный определённый интеграл:
\[
J_{u u}=\int \sigma_{2} \rho^{2} d s
\]

где $\sigma_{2}$ — линейная плотность, а $d s$ — элемент дуги той кривой, по которой распределены массы.

153. Зависимость между моментами инерции относительно параллельных осей. Для того чтобы найти зависимость между моментами инерции относительно параллельных осей, сравним момент инерции $J$ относительно какой-либо оси $A A^{\prime}$ с моментом инерции $J_{C}$ тех же масс относительно оси $C C^{\prime}$, параллельной $A A^{\prime}$, но проходящей через центр масс $C$ рассматриваемњх масс; расстояние между этими осями пусть будет $d$. Поместим начало координат в центре масс $C$ (фиг. 99), ось $z$

направим по оси $C C^{\prime}$, а плоскость $C y z$ проведём через ось $A A^{\prime}$; обозиачим расстолния магериальной частиты $M$, от осей $C C^{\prime}$ и $A A^{\prime}$ соответственно через $\rho_{y}$ и $\rho_{v}$, а угол отрезка $C M_{v}$ с осью у назовём $\varphi_{v}$. Из треугольника $C M, B$ имеем
\[
\rho^{\prime 2}=\rho_{v}^{2}+d^{2}-2 d \rho_{v} \cos \varphi_{v}=\rho_{v}^{2}+d^{2}-2 d y_{v} .
\]

Умножив полученное выражение на массу $m_{v}$ частицы и просуммировав его по всем частицам, полуінм:
\[
\sum_{v=1}^{n} m_{v} p_{v}^{2}=\sum_{v=1}^{n} m_{v} p_{v}^{2}+d^{2} \sum_{v=1}^{n} m_{v}-2 d \sum_{v=1}^{n} m_{v} y_{v}
\]

По определению имеем
\[
\sum_{v=1}^{n} m_{v} \rho_{v}^{2}=J ; \quad \sum_{v=1}^{n} m_{v} \rho_{v}^{2}=J_{c}
\]

далее,
\[
\sum_{v=1}^{n} m_{1}=M
\]

есть масса всей системы; наконец, сумма $\sum_{v=1}^{n} m_{
u} y_{
u}$ входяцая в последнее слагаемое, равна нулю как статический момент относительно плоскости, проходящей через центр масс (§144):
\[
\sum_{v=1}^{n} m_{v} y_{v}=M y_{C}=0 .
\]

Следовательно, равенство (26.4) перепииется так:
\[
J=J_{c}+M \cdot d^{2} \text {, }
\]
т. с. момент инерии системы частии относительно некоторой оси равен моменту инерции системы относительно параллельной оси, проходячей через центр масс, сложенному с моментом инерции центра масс относительно данной оси в предположении,
$\Phi_{\text {HI: }} 99$.
что его масса равна массе всей системы. Доказанное положение носит названне теоремы Гюйгенса-Штейнера (Huygtrens-Steiner).

Из доказанной теоремы легко вывести соотношение между моментами инериин $J$ и $J_{1}$ относительно любых параллельных осей $A A^{\prime}$ и $A_{1} A_{1}^{\prime}$; действигельно, цля момента $J_{1}$ имеем по теореме (26.5)
\[
J_{1}=J_{C}+M \cdot d_{1}^{2}
\]

если $d_{1}$ ознатает расстояние между осями $A A^{\prime}$ и $C C^{\prime}$; вычтя последнее. равенство из равенства (26.5), голучим искомую зависимость:
\[
J-J_{1}=M \cdot\left(d^{2}-d_{1}^{2}\right) .
\]

154. Зависимость между моментами инериии относительно осей, проходящих через данную точку. Произведения инериии. Эллипсоид инерции. Сравним теперь между собой моменты инерции некоторой системы ча.тии относительно осей, проходящих через одну и ту же

точку. Берём эту точку за начало декартовых осей $O x y z$ (фиг. 100). Пусть направление произвольной оси $O u$ характеризуется единичным век тором
\[
u^{0}=a x^{0}+\beta y^{0}+\gamma z^{0} ;
\]

величины $\alpha, \beta, \gamma$, очевидно, являются направляющими косинусами этой оси. Радиус-вектор какой-либо частицы $M_{v}$ с массой $m_{v}$ обозначим $r_{v}$, а расстояние этой частицы от оси $u$ назовём $\rho_{.}$. Момент инерции системы частии относительно оси $O u$ имеет выражение
\[
J_{u u}=\sum_{v} m_{v} \rho_{v}^{2}
\]

но
\[
\rho_{v}^{2}=r_{v}^{2}-\left(r_{v} \cdot u^{0}\right)^{2}=\left(r_{v} \times u^{0}\right)^{2}
\]
[см. формулу (1.30) на стр. 11]. Выразив здесь векторное произведение через проекции сомножителей, получаем:

Фиг. 100.
\[
\rho_{v}^{2}=\left(y_{v} \gamma-z_{v} \beta\right)^{2}+\left(z_{v} \alpha-x_{v} \gamma\right)^{2}+\left(x_{v} \beta-y_{v} \alpha\right)^{2},
\]

или, если расположить выражение по степеням косинусов,
\[
\begin{array}{r}
\rho_{v}^{2}=\left(y_{v}^{2}+z_{v}^{2}\right) \alpha^{2}+\left(z_{v}^{2}+x_{v}^{2}\right) \beta^{2}+\left(x_{v}^{2}+y_{v}^{2}\right) \gamma^{2}-2 y_{
u} \mathcal{z}_{v} \beta \gamma-2 z_{v} x, \gamma \alpha- \\
-2 x y_{
u} \alpha .
\end{array}
\]

Умножим это равенство на $m$, и просуммируем по всем частицам системы. Введём обозначения:
\[
\begin{array}{l}
J_{x x}=\sum_{v} m_{v}\left(y_{v}^{2}+z_{\vee}^{2}\right)_{,} \quad J_{v v}=\sum_{v} m_{v}\left(z_{v}^{2}+x_{v}^{2}\right), \\
J_{z z}=\sum_{v} m_{v}\left(x_{v}^{2}+y_{v}^{2}\right), \\
J_{y z}=\sum_{\vee} m_{v} y_{v} z_{v}, \quad J_{z x}=\sum_{v} m_{,} z_{v} x_{v}, \\
J_{x y}=\sum_{v} m_{v} x_{v} y_{v} ; \\
\end{array}
\]

в результате мы получим:
\[
J_{x z}=J_{x x} \alpha^{2}+J_{y y} \beta^{2}+J_{z z} \gamma^{2}-2 J_{y z} \beta \gamma-2 J_{z x} \gamma \alpha-2 J_{x y} \alpha \beta .
\]

Величины $J_{x x} J_{y y} J_{z z}$ представляют собой моменты инерции рассматриваемой сисгемы частиц относительно осей $O x, O y$ и $O z$. Величины $J_{y z}$, $J_{z x}, J_{x y}$ соответственно носят название произведений инериии, или центробежных моментов инерции относительно осей $y$ н $z, z$ и $x$, $x$ и $y$; в обозначениях произведений инерции, когда это почему-либо удобно, мы будем ставить индексы и в обратном поряке, иными словами, $J_{x y}=J_{v x}$ и т. д. В случае непрерывно распределённых масс все суммы (26.6) заменяются интегралами по типу формулы (26.3). Полученное нами выражение (26.7) показывает, что если мы знаем три момента инерции и три произведения инерции относительно координатных осей, то сумеем вылислить моменты инерции относительно любой оси, проходящей через начало координат.

Чтобы нагляднее выразить зависимость величины момента инерции от направления оси, употребим следующий геометрический приём. Огложим от начала координат по направлению взятой оси вектор $r$, длина которого обратно пропорциональна корню квадратному из соогветственного момента инерции, т. е. пусть
\[
r=\frac{l}{\sqrt{J_{u n}}},
\]

где $l$ — некоторая постоянная, определяющая масщтаб построения, и найдём геометрическое место конца этого вектора. Обозначим координаты конца построенного нами вектора $x, y, z$; тогда
\[
x=\frac{l}{\sqrt{J_{m i}}} a, \quad y=\frac{l}{\sqrt{J_{m u}}} \beta, \quad z=\frac{l}{\sqrt{J_{u p}}} \gamma .
\]

Если теперь умножим обе части равенства (26.7) на $\frac{l^{2}}{J_{u a}}$ и воспользуемся соотношениями (26.9), то легко найдём:
\[
J_{x x} x^{2}+J_{y y} y^{2}+J_{z z} z^{2}-2 J_{y z} y z-2 J_{z x} z x-2 J_{x y} z y=l^{2} .
\]

Мы видим, что геометрическим местом конца построенного вектора служит некоторая центральная поверхность второго порядка. Как видно из формул (26.8) и (26.7), поверхность эта не может иметь бесконечно удалённых точек за исключением того случая, когда все материальные частицыі лежат на прямой, проходящей через взятый полюс. В последнем случае, как легко показать, поверхность (26.10) обращается в цилиндр вращения. При всяком другом расположении масс рассматриваемая центральная поверхность, следовательно, может быгь только эллипсоидом; поэтому она и носит название эллипсоида инердии, соответствующего взятому полюсу. Собственно говоря, данному полюсу соответствует бесчисленное множество подобных друг другу эллипсоидов инерции, отличающихся один от другого значениями постоянной $l$; иногда из семейства этих подобных погерхностей выделяется как представительница та, уравнение которой в какомлибо отношении пишется наиболее просто; так, если принять $l=1$, то уравнение эллипсоида инерции представится в таком виде:
\[
J_{x x} x^{2}+J_{v y} y^{2}+J_{z z} z^{2}-2 J_{y z} y z-2 J_{z x} z x-2 J_{x y} x y=1 .
\]

Если за оси координат возьмём направления главных диаметров, или осей поверхности (26.10), то уравнение эллипсоида напишется следующим образом:
\[
J_{x x} x^{2}+J_{y, y} y^{2}+J_{z z} z^{2}=l^{2} ;
\]

следовательно, вместо формулы (26.7) мы в этом случае будем иметь:
\[
J_{t a}=J_{x x} \alpha^{2}+J_{y
u} \beta^{2}+J_{z z} \gamma^{2} .
\]

Из возможности привести уравнение эллипсоида инерции к виду (26.12) мы заключаем, что через любой полюс всегда проходят три такие взаимно перпендикулярные оси, что для них произведения инерции обрацаются в нули. Эти оси носят название главных осей инерции для взятого полюса, а моменты инерции, им соответствующие, называются главнымимоментами инерции для данного полюса. Когда за полюс взят центр масс, то эллипсоид инерции, главные оси и главные мо-

менты инерции соогветсгвенно называют центральным эллипсоидом, главными центральными осями главными центральными моментами инерции.

Когда радиус-вектор $r$, проведённый из центра эллипсоида к точке, лежащей на его поверхности, приходит в совпадение с одной из осей эллипсоида, его модуль $r$, как известно, принимает стационарное значение, т. е. для этих положений
\[
d r=0 ;
\]

при совпадении с большей и меньшей осями трёхосного эллипсоида $r$ соответственно проходит через максимум и минимум, при совпадении со средней осью $r$ проходит через так называемый минимакс. На основании формулы (26.8) отсюда, как следствие, вытекает, что и момент инерции $J_{\text {tu }}$ относительно каждой из главных осей инерции принимает стационарное значение. При этом, когда эллипсоид инерции для взятого полюса трёосный, го все три главные момента инерции различны и один из них имеет наибольшее, а другой наименьшее значсние среди моментов инерции относительно всех других осей, проходящих через взятый полюс. Когда эллипсоид инеріии для взятого полюса представляет собой поверхность вращения, то два из главных моментов инерции становятся равными, причём попрежнему момент инерции относительно любой другой (не главной) оси, проходящей через данный полюс, содержится между двумя неравными друг другу главными моментами инерции. Накснец, когда эллипсоид инерции обращается в сферу, все моменты инерции относительно осей, проходящих через данный полюс, равны между собой.

Заметим, что не всякий эллипсоид может служить эллипсоидом инерции. Пусть главные моменты инериии так располагаются по своей всличине, что

тотда во всяком случае
\[
\begin{array}{l}
J_{x x} \leqslant J_{v y} \leqslant J_{z z} ; \\
J_{x x}+J_{y y} \geqslant J_{z z} ;
\end{array}
\]

это видно из формул (26.6), поскольку
\[
J_{x x}+J_{v y}=J_{z z}+2 \sum_{v} m_{v} z_{v}^{2}
\]

так как в уравнении (26.12) коэффициенты $J_{x x} J_{y y} J_{z z}$ обратно пропорциональны квадратам полуосей поверхности, то отсюда заключаем, что эллипсоид
\[
\frac{x^{2}}{a^{2}}+\frac{y^{2}}{b^{2}}+\frac{z^{2}}{c^{2}}=1
\]

при
\[
a \geqslant b \geqslant c
\]

может служить эллипсоидом инерции лишь при условии
\[
\frac{1}{a^{2}}+\frac{1}{b^{2}} \geqslant \frac{1}{c^{2}} \text {. }
\]
155. Радиус инерции. Гирационный эллипсоид. Выражению для момента инерции $J$ материальной системы дают иногда. следующий вид:
\[
J=M \rho^{2},
\]

где $M$ — масса всей системы. Из размерности ве.иччны $J$ нетрудно заключить, что $\rho$ представляет-собой длину. Длина эта носит название радиуса, или плеча, инерции данных масс по отношению ко взятой оси. Радиусы инершин относительно главных осей называются главным Ейли уравненге (26.10) эллипсоида инерции отнесём к главным диаметрам, дадим постоянной $l^{2}$ значение $M$ и введём главные радиусы инерции $a, b, c$, т. е. положим
\[
J_{x x}^{\prime}=M a^{2} ; \quad J_{y y}=M b^{2} ; \quad J_{z z}=M c^{2},
\]

то получим:
\[
a^{2} x^{2}+b^{2} y^{2}+c^{2} z^{2}=1 .
\]

Кроме эллипсоида инерции, роль в механике играет и другая поверхность, носящая название гирационного эллипссида. Фиг. Iui.

Получается она из эллипсоида инерции’следующим построением. Берём на эллипсонде инерции
\[
f=J_{x x} x^{2}+J_{y y} y^{2}+J_{z z} z^{2}-l^{2}=0
\]

произвольную точку $m$,. определяемую радиусом-вектором
\[
r=x x^{0}+y y^{0}+z z^{0},
\]

и строим в этой точке касательную плоскость (фиг. 101). На эту плоскость опускаем из центра $O$ эллипсоида перпендикуляр $O D=\delta$ и производим инверсию точки $D$ относительно сферы произвольного радиуса $R$, т. е. на том же луче $O D$ эткладываем вектор $\overline{O \mu}=\bar{\rho}$, по модулю равный
\[
\rho=\frac{R^{2}}{8} \text {. }
\]

Геометрическое место точки $\mu$ и будет искомой поверхностью. Для получения ее уравнения заметим, что,
\[
\delta=r \cdot \frac{\operatorname{grad} f}{|\operatorname{grad} f|}
\]

следовательно,
\[
\rho=\frac{R^{2}|\operatorname{grad} f|}{r \cdot \operatorname{grad} f}
\]

откуда
\[
\bar{\rho}=\frac{R^{2} \operatorname{grad} f}{r \cdot \operatorname{grad} f} .
\]

Выразив все входящие в это равенство векторы через их проекции, получим соотношение
\[
x_{\mu} x^{0}+y_{\mu} y^{0}+z_{\mu} z^{0}=\frac{R^{2}}{l^{2}}\left(J_{x x} x x^{0}+J_{y y} y y^{0}+J_{z z} z z^{0}\right) ;
\]

следовательно,
\[
x=\frac{l^{2} x_{\mu}}{R^{2} J_{x x}}, \quad y=\frac{l^{2} y_{\mu}}{R^{2} J_{y y}}, \quad z=\frac{l^{2} z_{\mu}}{R^{2} J_{z y}} .
\]

Вставив эти значения $\boldsymbol{x}, \boldsymbol{y}, \boldsymbol{z}$ в уравнение (26.17) и опустив в окончательном результате индекс $\mu$, мы получим уравнение искомой поверхности в виде
\[
\frac{x^{2}}{J_{x x}}+\frac{y^{2}}{J_{y y}}+\frac{z^{2}}{J_{z z}}=\frac{R^{4}}{l^{2}} .
\]

Если преобразование, посредством которого мы из эллипсоида инерции получили гирационный эллипсоид, применить к последнему, то снова получится эллипсоид инерции; вследствие этого оба эти эллипсоида называются взаимными. Если радиус инверсии $R$ взять равным единице, то эллипсоидом, взаимным с эллипсоидом (26.16), будет следующий
\[
\frac{x^{2}}{a^{2}}+\frac{y^{2}}{b^{2}}+\frac{z^{2}}{c^{2}}=1 \text {. }
\]

156. Эллиптические координаты. Прежде чем перейти к изучению распределения главных осей инерции в пространстве, познакомимся с той системой криволинейных координат, которая носит название эллиптической. Рассмотрим равенство
\[
\frac{x^{2}}{a^{2}+\lambda}+\frac{y^{2}}{b^{2}+\lambda}+\frac{z^{2}}{c^{2}+\lambda}-1=0 .
\]

Когда $\lambda$ имеет постоянное значение, а $x, y, z$ представляют собой текущие координаты, то это равенство служит уравнением некоторой центральной поверхности второго порядка. Если же дадим $x, y, z$ какиелибо постоянные значения, то оно обращается в кубичное уравнение относительно $\lambda$, имеющее, как сейчас убедимся, три действительных корня. Пусть
\[
c^{2}<b^{2}<a^{2} .
\]

Обозначим $Q(\lambda)$ левую часть уравнения (26.20):
\[
Q(\lambda)=\frac{x^{2}}{a^{2}+\lambda}+\frac{y^{2}}{b^{2}+\lambda}+\frac{z^{2}}{c^{2}+\lambda}-1 .
\]

Положим $\lambda=-a^{2}+\varepsilon^{2}$, где $\varepsilon^{2}$ — сколь угодно малая положительная величина. Тогда знак всего выражения $Q(\lambda)$ будет зависеть от знака первого члена и, следовательно, будет
\[
Q\left(-a^{2}+\varepsilon^{2}\right)>0 \text {. }
\]

Голожив, далее, $\lambda=-b^{2}-\varepsilon^{2}$, найдём, что
\[
Q\left(-b^{2}-\varepsilon^{2}\right)<0 .
\]

Из этих неравенств заключаем, что между $-a^{2}$ и $-b^{2}$ уравнение $Q(i)=0$ имеет, по крайней мере, один действительный корень; обозначим его $\lambda_{1}$. Дав теперь $\downarrow$ значения $-b^{2}+\varepsilon^{2}$ и затем $-c^{2}-\varepsilon^{2}$, убеждаемся, что
\[
Q\left(-b^{2}+\varepsilon^{2}\right)>0, \quad Q\left(-c^{2}-\varepsilon^{2}\right)<0 ;
\]

следовательно, уравнение $Q(\lambda)=0$ имеет действительный корень $\lambda_{2}$ между $-b^{2}$ и $-c^{2}$. Наконец, положив $\lambda=-c^{2}+\varepsilon^{2}$ и $\lambda=+\infty$, видим, что
\[
Q\left(-c^{2}+\varepsilon^{2}\right)>0, \quad Q(+\infty)=-1<0,
\]

т. е. третий действительный корень $\lambda_{3}$ лежит между $-c^{2}$ и $+\infty$. Итак, корни уравнения $Q(\lambda)$ можно раснолюжить в таком порядке:
\[
-a^{2}<\lambda_{1}<-b^{2}<\lambda_{2}<-c^{2}<\lambda_{3}<+\infty .
\]

Три уравнения:
\[
\begin{array}{l}
f_{1}=\frac{x^{2}}{a^{2}+\lambda_{1}}+\frac{y^{2}}{b^{2}+\lambda_{1}}+\frac{z^{2}}{c^{2}+\lambda_{1}}-1=0, \\
f_{2}=\frac{x^{2}}{a^{2}+\lambda_{2}}+\frac{y^{2}}{b^{2}+\lambda_{2}}+\frac{z^{2}}{c^{2}+\lambda_{2}}-1=0, \\
f_{3}=\frac{x^{2}}{a^{2}+\lambda_{3}}+\frac{y^{2}}{b^{2}+\lambda_{3}}+\frac{z^{2}}{c^{2}+\lambda_{3}}-1=0
\end{array}
\]

дают нам три центральные поверхности второго порядка, проходящие через данную точку. Если примем во внимание неравенства (26.21) и (26.23), то увидим, что первая из них будет двухполостным гиперболоидом, вторая — однополостным, третья — эллипсоидом. Так как между квадратами нолуосей этих поверхностей имеют место соотношения
\[
\left(a^{2}+\lambda_{i}\right)-a^{2}=\left(\dot{o}^{2}+\lambda_{i}\right)-b^{2}=\left(c^{2}+\lambda_{i}\right)-c^{2}=\lambda_{l} \quad(i=1,2,3),
\]

то все три поверхности софокусны с основным эллипсоидом:
\[
\frac{x^{2}}{a_{2}}+\frac{y^{2}}{b_{2}}+\frac{z^{2}}{c_{2}}-1=0 .
\]

Величины $\lambda_{1}, \lambda_{2}, \lambda_{3}$ можно принять за три криволинейные координаты данной точки, определяемой декартовыми координатами $x, y, z$. Всякой точке будет соответствовать своя система значений этих координат. Наоборот, давши некоторые значения величинам $\lambda_{1}, \lambda_{2}$ и $\lambda_{3}$, мы можем определить положение точки, впрочем, не однозначно, потому что поверхности $(26.24),(26.25)$ и (26.26) пересекаются не в одной, а в восьми точках. Координаты $\lambda_{1}, \lambda_{2}, \lambda_{3}$ носят название эллиптических.

Когда величина $\lambda_{1}$ проходит систему значений между $-a^{2}$ и $-b^{2}$, двухполостный гиперболоид (26.24) изменяется от плоскости $O y z$ до той части плоскости $O z x$, которая ‘лежит по положительную сторону гиперболы
\[
\frac{x^{2}}{a^{2}-b^{2}}-\frac{z^{2}}{b^{2}-c^{2}}-1=0,
\]
т. е. той части плоскости, ксэрдинаты любой точки которой обращают левую часть уравнения (26.27) в положительную величину. Когда величина $\lambda_{2}$ проходит систему значений между $-b^{2}$ и $-c^{2}$, однополостный гиперболоид (26.25) изменяется от гиперболической пластинки, лежащей в плоскости Ozx по отрицательную сторону гипербозы (26.27), до той части плоскости $O x y$, которая лежит по положительную сторону эллипса
\[
\frac{x^{2}}{a^{2}-c^{2}}+\frac{y^{2}}{b^{2}-c^{2}}-1=0 .
\]

Наконец, когда величина $\lambda_{3}$ проходит систему значений от $-c^{2}$ до $+\infty$, эллипсоид ( 26.26 ) изменяется от эллиптической пластинки, лежащей в плоскости $O x y$ и представляющей собой место точек, лежащих по положительную сторону эллипса (26.28), до сферы бесконечно болыного радиуса.

Обе кривые (26.27) и (26.28) лежат во взанмно перпендикулярных плоскостях, проходят соответственно одна через фокусы другой и представляют собой геометрические места тех точек, для которых две эллип-

тические координзты равны между собой: для гиперболы (26.27) имеем $\lambda_{1}=\lambda_{2}$, а для эллипса (26.28). $\lambda_{2}=\lambda_{3}$.

Докажем, что эллиптическая система координат представляет собой систему ортогональную. Координатные поверхности (26.24), (26.25), (26.26) встречаются ортогонально, если соблюдены условия
\[
\operatorname{grad} f_{1} \cdot \operatorname{grad} f_{2}=0, \operatorname{grad} f_{2} \cdot \operatorname{grad} f_{3}=0, \quad \operatorname{grad} f_{3} \cdot \operatorname{grad} f_{1}=0 .
\]

Если мы вычтем уравнение (26.25) из (26.24) и полученный результат сократим на $\lambda_{2}-\lambda_{1}$, то найдём:
\[
\frac{x^{2}}{\left(a^{2}+\lambda_{1}\right)\left(a^{2}+\lambda_{2}\right)}+\frac{y^{2}}{\left(b^{2}+\lambda_{1}\right)\left(b^{2}+\lambda_{2}\right)}+\frac{z^{2}}{\left(c^{2}+\lambda_{1}\right)\left(c^{2}+\lambda_{2}\right)}=0 ;
\]

а это равенство и будет первым из условий (26.29), выраженным в скалярной форме. То же можно показать относительно каждой из других двух пар координатных поверхностей.

Чтобы выразить декартовы координаты явным образом через эллиптические, снова обращаемся к уравнедниям (26.24), (26.25) и (26.26). Исключив из них сперва одну из координат, например $z$, получим:
\[
\begin{array}{l}
\frac{\left(c^{2}-a^{2}\right) x^{2}}{\left(a^{2}+\lambda_{1}\right)\left(a^{2}+\lambda_{3}\right)}+\frac{\left(c^{2}-b^{2}\right) y^{2}}{\left(b^{2}+\lambda_{1}\right)\left(b^{2}+\lambda_{3}\right)}=1, \\
\frac{\left(c^{2}-a^{2}\right) x^{2}}{\left(a^{2}+\lambda_{2}\right)\left(a^{2}+\lambda_{3}\right)}+\frac{\left(c^{2}-b^{2}\right) y^{2}}{\left(b^{2}+\lambda_{2}\right)\left(b^{2}+\lambda_{3}\right)}=1 .
\end{array}
\]

Исключив теперь отсюда $y^{2}$, получим первое из искомых выражении:
\[
x= \pm \sqrt{\frac{\left(a^{2}+\lambda_{1}\right)\left(a^{2}+\lambda_{2}\right)\left(a^{2}+\lambda_{3}\right)}{\left(a^{2}-b^{2}\right)\left(a^{2}-c^{2}\right)}},
\]

подобным образом найдём и остальные:
\[
\begin{array}{l}
y= \pm \sqrt{\frac{\left(b^{2}+\lambda_{1}\right)\left(b^{2}+\lambda_{2}\right)\left(b^{2}+\lambda_{3}\right)}{\left(b^{2}-c^{2}\right)\left(b^{2}-a^{2}\right)}}, \\
z= \pm \sqrt{\frac{\left(c^{2}+\lambda_{1}\right)\left(c^{2}+\lambda_{2}\right)\left(c^{2}+\lambda_{3}\right)}{\left(c^{3}-a^{2}\right)\left(c^{2}-b^{2}\right)}} .
\end{array}
\]

Если принять во внимание неравенства (26.21) и (26.23), то нетрудно увидеть, что подрадикальные выражения здесь всегда положительны.

В заключение выведсм некоторые вспомогательные формулы, которыми нам придётся впоследствии воспользоваться. Прежде всего заметим, что равенство (26.22) можно переписать так:
\[
Q(\lambda) \cdot\left(a^{2}+\lambda\right)\left(b^{2}+\lambda\right)\left(c^{2}+\lambda\right)=-\left(\lambda-\lambda_{1}\right)\left(\lambda-\lambda_{2}\right)\left(\lambda-\lambda_{3}\right) .
\]

Действительно, с одной стороны, корнями уравнения $Q(\lambda)=0$ служат $\lambda_{1}, \lambda_{2}$ н $\lambda_{3}$, а с цругой стороны, если все четыре члена выражения $Q(\lambda)$ привести к общему знаменателю
\[
\left(a^{2}+\lambda\right)\left(\dot{b}^{2}+\lambda\right)\left(c^{2}+\lambda\right),
\]

то коэффициент при высшей степени $\lambda$ в числителе будет равен отрицательной единице, что мы и имеем в правой части. Из равенства (26.34) внтекает, что
\[
Q(\lambda)=-\frac{\left(\lambda-i_{1}\right)\left(\lambda-\lambda_{2}\right)\left(\lambda-\lambda_{3}\right)}{\left(a^{2}+\lambda\right)\left(b^{2}+\lambda\right)\left(c^{2}+\lambda\right)} .
\]

Найдём теперь значение производной $-\frac{\partial Q}{\partial \lambda}$ для $\lambda=\lambda_{\mathbf{1}}$. Нетрудно видеть, что
\[
-\frac{\partial Q}{\partial \lambda}=\frac{\left(\lambda-\lambda_{2}\right)\left(\lambda-\lambda_{3}\right)}{\left(a^{2}+\lambda\right)\left(b^{2}+\lambda\right)\left(c^{2}+\lambda\right)}+M \cdot\left(\lambda-\lambda_{1}\right),
\]

где через $M$ обозначены члены, не обращающиеся в бесконечность для $\lambda=\lambda_{1}$. Дав $\lambda$ значение $\lambda_{1}$, находим отсюда:
\[
\left\{-\frac{\partial Q}{\partial \lambda}\right\}_{A=\lambda_{1}}=\frac{\left(\lambda_{1}-\lambda_{2}\right)\left(\lambda_{1}-\lambda_{3}\right)}{\left(a^{2}+\lambda_{1}\right)\left(b^{2}+\lambda_{1}\right)\left(c^{2}+\lambda_{1}\right)} .
\]

С другой стороны, из выражения (26.22) имеем
\[
\left\{-\frac{\partial Q}{\partial \lambda}\right\}_{h=\lambda_{1}}=\frac{x^{2}}{\left(a^{3}+\lambda_{1}\right)^{2}}+\frac{y^{2}}{\left(b^{2}+\lambda_{1}\right)^{2}}+\frac{z^{2}}{\left(c^{2}+\lambda_{1}\right)^{2}} .
\]

Введя для $\left\{-\frac{\partial Q}{\partial \lambda}\right\}_{\wedge=\lambda_{1}}$ обозначенне $L_{1}^{2}$, мы получаем, следовательно, такое соотношение:
\[
L_{1}^{2}=\frac{x^{2}}{\left(a^{2}+\lambda_{1}\right)^{2}}+\frac{y^{2}}{\left(b^{2}+\lambda_{1}\right)^{2}}+\frac{z^{2}}{\left(c^{2}+\lambda_{1}\right)^{2}}=\frac{\left(\lambda_{1}-\lambda_{2}\right)\left(\lambda_{1}-\lambda_{8}\right)}{\left(a^{2}+\lambda_{1}\right)\left(b^{2}+\lambda_{1}\right)\left(c^{2}+\lambda_{1}\right)} .
\]

К этой формуле, очевидно, могут быть присоединены ешё две:
\[
\begin{array}{l}
L_{2}^{2}=\frac{x^{2}}{\left(a^{2}+\lambda_{2}\right)^{2}}+\frac{y^{2}}{\left(b^{2}+\lambda_{3}\right)^{2}}+\frac{z^{2}}{\left(c^{2}+\lambda_{9}\right)^{2}}=\frac{\left(\lambda_{2}-\lambda_{3}\right)\left(\lambda_{2}-\lambda_{1}\right)}{\left(a^{2}+\lambda_{2}\right)\left(b^{2}+\lambda_{2}\right)\left(c^{2}+\lambda_{2}\right)}, \\
L_{3}^{2}=\frac{x^{2}}{\left(a^{2}+\lambda_{3}\right)^{2}}+\frac{y^{2}}{\left(b^{2}+\lambda_{3}\right)^{2}}+\frac{z^{2}}{\left(c^{2}+\lambda_{3}\right)^{2}}=\frac{\left(\lambda_{3}-\lambda_{1}\right)\left(\lambda_{3}-\lambda_{2}\right)}{\left(a^{2}+\lambda_{3}\right)\left(b^{2}+\lambda_{3}\right)\left(c^{2}+\lambda_{3}\right)} .
\end{array}
\]

157. Распределение главных осей инерции в пространстве. Займёмся теперь решением следующего вопроса: как по данным главным центральным моментам инерции каких-либо масс и по данным направлениям главных центральных осей инерции найти момент инерции тех же масс огносительно любой оси, а также опрәделить направления главных осей инерции для люєого полюса.

Возьмём центр масс $C$ за начало декартовых координат и совместим координатные оси с главными центральными осями инерции (фиг. 102). Пусть моменты инерции относительно осей $C x, C y, C z$ будут соответственно $J_{x x}, J_{y y}, J_{z z}$. Берём произвольный полюс $K$ с радиусом-вектором
\[
r_{k}=x_{k} x^{0}+y_{k} y^{0}+z_{k} z^{0} ;
\]

пусть через него проходит ось $K u$, определяемая единичным вектором
\[
u^{0}=\alpha x^{0}+\beta y^{0}+\gamma z^{0} .
\]

Искомый момент инерции относительно оси $K u$ обозначим $J_{u n}$, а момент инерции относительно оси $\mathrm{Cu}^{\prime}$, ей параллельной и проходящей через на-

чало координат, т. е. центр масс, пусть будет $J_{n u}^{\prime}$. Тогда по теореме Гюйгенса [см. формулу (26.5) на стр. 225)] имеем
\[
J_{u z}=J_{t u}^{\prime}+M d^{2},
\]

где $M$ — масса всей системы частиц, а $d$-расстояние между осями $K u$ и $\mathrm{Cu}^{\prime}$. Для квадрата этого расстояния получаем следуюшее выражение:
\[
\begin{aligned}
d^{2}=r_{k}^{2}-\left(r_{k} \cdot u^{
atural}\right)^{2} & =\left(r_{k} \times u^{0}\right)^{2}= \\
& =\left(y_{k} \gamma-z_{k} \beta\right)^{2}+\left(z_{k} \alpha-x_{k} \gamma\right)^{2}+\left(x_{k} \beta-y_{k} a\right)^{2},
\end{aligned}
\]

или
\[
\begin{aligned}
d^{2}=\left(y_{k}^{2}+z_{k}^{2}\right) \alpha^{2}+\left(z_{k}^{2}+x_{k}^{2}\right) \beta^{2}+ & \left(x_{k}^{2}+y_{k}^{2}\right) \gamma^{2}-2 y_{k} z_{k} \beta \gamma- \\
& -2 z_{k} x_{k} \gamma \alpha-2 x_{k} y_{k} \alpha \beta .
\end{aligned}
\]

С другой стороны, по формуле (26.13) мы имеем
\[
J_{u^{\prime} u^{\prime}}=J_{x x} a^{2}+J_{y y} \beta^{2}+J_{z z} \gamma^{2} .
\]

Введём обозначения:
\[
\left.\begin{array}{l}
J_{x x}^{(k)}=J_{x x}+M\left(y_{k}^{2}+z_{k}^{2}\right), \\
J_{y y}^{(k)}=J_{v y}+M\left(z_{k}^{2}+x_{k}^{2}\right), \\
J_{z z}^{(k)}=J_{z z}+M\left(x_{k}^{2}+y_{k}^{2}\right), \\
J_{y z}^{(k)}=M y_{k} z_{k}, \\
J_{z x}^{(k)}=M z_{k} x_{k}, \\
J_{x y}^{(k)}=M x_{k} y_{k} ;
\end{array}\right\}
\]

тогда, объединив полученные результаты, мы можем, вместо формулы (26.38), написать
\[
J_{t u}=J_{x x}^{(k)} \alpha^{2}+J_{y y}^{(k)} \beta^{2}+J_{z z}^{(k)} \gamma^{2}-2 J_{y z}^{(k) \beta \gamma}-2 J_{z \lambda}^{(k)} \gamma \alpha-2 J_{x y}^{(k)} \alpha \beta .
\]

Эта формула и определяет по нацим данным нскомый момент инерции $J_{t z *}$.
Направления главных осей инерции для полюса $K$ можно найти из условия стационарности момента инерции $J_{\text {mu }}$ относительно каждой из 9тих осей. Формула (26.42) выражает момент инерции $J_{u a}$ в функции переменных $\alpha, \beta, \gamma$, но последние связаны между собой очевидным соотношением
\[
f=\alpha^{2}+\beta^{2}+\gamma^{2}-1=0 .
\]

По известной теореме математического анализа, когда исследуемая функция принимает стационарное значение, её частные производные пропорциональны частным производным от функции, выражающей связь между переменными. Поэтому в нашем случае, называя множитель пропорциональности $J$, мы имеем
\[
\frac{\partial J_{u t}}{\partial a}=J \frac{\partial f}{\partial \alpha}, \quad \frac{\partial J_{u z}}{\partial \beta}=J \frac{\partial f}{\partial \beta}, \quad \frac{\partial J_{u t}}{\partial \gamma}=J \frac{\partial f}{\partial \gamma} ;
\]

напишем эти равенства в развёрнутом виде и перенесём все члены в левые части; в результате мы получим:
\[
\left.\begin{array}{rl}
\left(J_{x x}^{k k}-J\right) \alpha-J_{x y}^{(k)} \beta-J_{x z}^{(k) \gamma} & \doteq 0, \\
-J_{y x}^{k l} \alpha+\left(J_{i y}^{(k)}-J\right) \beta-J_{y z}^{(k) \gamma} & =0, \\
-J_{z x}^{(k)} \alpha+J_{z y}^{k y} \beta+\left(J_{z z}^{(k)}-J\right) \gamma & =0 .
\end{array}\right\}
\]

Условием, при которо: эта система однородных уравнений удовлетворяется значениями $\alpha, \beta, \gamma$, не равными нулю одновременно, является равенство нулю определителя системы:
\[
\left|\begin{array}{ccc}
J_{ \pm x}^{(k y}-J & -J_{x y}^{k)} & -J_{x z}^{(k)} \\
-J_{y x}^{(k)} & J_{y y}^{(k)}-J & J_{y z}^{(k)} \\
-J_{z x}^{(k)} & -J_{z y}^{(k)} & J_{z z}^{(k)}-J
\end{array}\right|=0 .
\]
множитель $J$ явлнется его корнем. В высшей алгебре доказывается, что если коэффициенты векового уравнения действительны (что в нашем случае имеет место), то все корни его также действительны. С другой стороны, если равенства (26.43) умножим соответственно на $\alpha, \beta, \gamma$, затем сложим и сравним результат с выражением (26.42), то убедимся, что $J$ представляет собой одно из искомых стационарных значений функции $J_{u u}$. Следовательно, корни $J_{1}, J_{2}, J_{3}$ векового уравнения (26.43) служат главными моментами инерции для полюса $K$; соответствуюшие этим корням направления главных осей инерции определяются косинусами $\alpha_{v}, \beta_{v}, \gamma_{v}$, где $
u=1,2,3$; их значения найдутся из уравнений (26.43), если туда вместо $J$ вставим соответствующий корень $J_{
u}$ :

Чтобы уяснить, по какому закону меняются направления главных осей инерции для различных точек пространства, преобразуем уравнения (26.43) к другому виду. Первое из них, если вернуться к прежним обозначениям (26.40) и (26.41), можно теперь написать так:
\[
\left[J_{x x}+M\left(y_{k}^{2}+z_{k}^{2}\right)-J_{v}\right] \alpha_{\chi}=M x_{k} y_{k} \beta_{v}+M x_{k} z_{k} \gamma_{v} .
\]

Прибавим к обеим частям $M x_{k}^{i} a_{v}$; введём радиусы инерции, положив
\[
J_{x x}=M a^{2}, \quad J_{v y}=M b^{2}, \quad J_{z z}=M c^{2} ; \quad J_{v}=M \rho_{v}^{2} ;
\]

сдиничный вектор рассматриваемой главной оси инерцйи назовём
\[
u_{v}^{0}=\alpha_{v} x^{0}+\beta_{v} y^{0}+\gamma_{\gamma} z^{0},
\]

радиус-вектор точки $K$ назовём попрежнему
\[
r_{k}=x_{k} x^{0}+y_{k} y^{0}+z_{k} z^{0} ;
\]

наконец, для сокращения письма обозначим
\[
r_{k}^{2}-\rho_{i}^{2}=\lambda_{v}
\]

тогда после очевидных преобразований мы получим:
\[
a_{v}=\frac{u_{v}^{0} \cdot r_{k}}{a^{2}+\lambda_{v}} x_{k}
\]

Подобным образом найдём из остальных двух уравнений выражения
\[
\beta_{v}=\frac{u^{0} \cdot r_{k}}{b^{2}+\lambda_{v}} y_{k}, \quad \gamma_{v}=\frac{u^{0} \cdot r_{k}}{c^{2}+\lambda_{v}} z_{k} .
\]

Умножив последние три равенства соответственно на $x_{k}, y_{k}, z_{k}$, сложив и выполнив сокращение, мы получим:
\[
\frac{x_{k}^{2}}{a^{2}+\lambda_{
u}}+\frac{y_{k}^{2}}{b^{2}+\lambda_{
u}}+\frac{z_{k}^{2}}{c^{2}+\lambda_{
u}}-1=0,
\]

что мы короче запишем так:
\[
f_{v}\left(x_{k}, y_{k}, z_{k}\right)=0 .
\]

С другой стороны, умножив эти же равенства соответственно на $x^{0}, y^{0}, z^{0}$ и сложив, получим:
\[
u_{v}^{0}=u_{v}^{0} \cdot r_{k}\left(\frac{x_{k}}{a^{2}+\lambda_{v}} x^{0}+\frac{y_{k}}{b^{2}+\lambda_{v}} y^{0}+\frac{z_{k}}{c^{2}+\lambda_{v}} z^{0}\right),
\]

или
\[
\boldsymbol{u}_{v}^{0}=\boldsymbol{u}_{\vee}^{0} \cdot \boldsymbol{r}_{k} \operatorname{grad} f_{v} .
\]

Выраженне (26.48) показывает, что через взятый полюс проходят три поверхности 2-го порядка, соответствуюцие значениям $
u=1,2,3$ и софокусные с центральным гирационным эллипсоидом (26.19)
\[
\frac{x^{2}}{a^{2}}+\frac{y^{2}}{b^{2}}+\frac{z^{2}}{c^{2}}=1 .
\]

Параметры этих софокусных поверхностей связаны с главными моментами инериии для полюса $K$ соотношениями (26.45), которые можно переписать так:
\[
\lambda_{v}=r_{k}^{2}-\frac{J_{v}}{M} .
\]

Из соотношения (26.49) мы усматриваем, что направления главных осей инерции для полюса $K$ совпадают с нормалями к трём выше упомянутым софокусным поверхностям; при этом согласно формуле (26.48) нормаль к эллипсоиду служит осью наименьшего, а нормаль к двухполостному гиперболоиду осью наибольшего момента инерции ( $\S 156$ ).

Когда главные центральные моменты инерции $J_{x x}, J_{y y}, J_{z 2}$ не равны между собой, то эллипсоид инерии и и для полюса $K$ будет, вообще говоря, трёхосный; эллиясоид инерции для полюса $K$ обратится в поверхность врашения, лишь когда точка $K$ лежит на одной из кривых (26.27) или (26.28). В предельном случае $b \rightarrow c$ эги кривые пересекаются в точках оси $x$ с координатами
\[
x= \pm \sqrt{a^{2}-b^{2}},
\]

причём это единственный случай, когда рассмагриваемые кривые имеют общие точки; сами кривые при этом соответственны превращаются в отрезки оси $x$ :
\[
|x| \geqslant \sqrt{a^{2}-b^{2}} \quad \text { и }|x| \leqslant \sqrt{a^{2}-b^{2}} .
\]

Таким образом, эллипсонд инершии обращается в шар лишь дтя двух полюсов и притом только тогда, когда центральный гирационный элличсоид является яйцевидным эллипсоидом вращения (а не сфероидом) с полуосями
\[
a>b=c .
\]

Упомянутые два полюса лежат в таком случае на оси вращения гирационного эллипсоида на расстояниях от центра масс, равных (при принлтом маситтабе) $\sqrt{a^{2}-b^{2}}$.

В заключение заметим следуюшее. Как показывают формулы (26.46) и (26.47), проекции вектора $\boldsymbol{u}^{\prime}$ являются линейными функциями проекций

вектора $\boldsymbol{r}_{\boldsymbol{k}}$. Обозначим $\left\|\Delta_{\vee}\right\|$ матрицу рассматриваемого преобразования
\[
\left\|\Delta_{
u}\right\|=\left\|\begin{array}{ccc}
\frac{u_{v}^{0} \cdot r_{k}}{a^{2}+\lambda_{
u}} & 0 & 0 \\
0 & \frac{u_{v}^{0} \cdot r_{k}}{b^{2}+\lambda_{v}} & 0 \\
0 & 0 & \frac{u_{v}^{0} \cdot r_{k}}{c^{2}+\lambda}
\end{array}\right\| .
\]

Принято говорить, что вектор $\boldsymbol{u}_{v}^{0}$ получается аффинным преобразованием из вектора $r_{k}$ при помощи матришы $\left\|\Delta_{v}\right\|$; в записи это выглядит так:
\[
u^{0}=\left\|\Delta_{v}\right\| r_{k}
\]

эта векторная формула эквивалентна трём выше приведённым скалярным.

158. Момент инерции относительно плоскости. Моментом инерции данной системы материальных частии относительно плоскости называется сумма произведений масс частиц на квадраты их расстояний от рассматриваемой плоскости. На основании этого определения мы имеем следуюшие выражения для моментов инерции относительно координатных плоскостей:
\[
J_{\text {пл, } y z}=\sum_{v} m_{v} x_{v}^{2}, \quad J_{\text {пл. } z x}=\sum_{y} m_{v} y_{v}^{2}, \quad J_{\text {п., } x y}=\sum_{v} m_{v} z_{v .}^{2}
\]

В случае масс, расположенных сплошным образом, мы; вместо сумм, получаем интегралы:
\[
J_{\text {пл. } v z}=\iiint \sigma x^{2} d v \text { и т. д. }
\]

Для момента инерции огносительно плоскости можно было бы доказать теорему, аналогичную теореме Гюйгенса [формула (26.5) на стр. 255]:
\[
J=J_{C}+M d^{2} .
\]

Из выражений для моментов инерции относительно начала координат, осей координат и плоскостей координат видно, что они связаны между собой следующими соотношениями:
\[
\left.\begin{array}{rl}
J_{0} & =J_{\text {пл. } y z}+J_{\text {пл. } z x}+J_{\text {пл. } x y}=\frac{J_{x x}+J_{y y}+J_{z z}}{2}, \\
J_{x x} & =J_{\text {пл. } z x}+J_{\text {пл. } x y} \text { и } \quad \text { г. д., } \\
J_{x x}+J_{y y} & =J_{z z}+2 J_{\text {пл. } x y} \geqslant J_{z z} \text { и } \text { т. д. }
\end{array}\right\}
\]

Если массы расположены в плоскости (в плоскости $O x y$ ), то как частный случай последнен формулы мы нмеем
\[
J_{x x}+J_{y y}=J_{z z}
\]

или, как говорят, сумма двух экваториальных моментов инерции относительно взаимно перпендикулярных осей равняется полярному моменту инерции относительно точки их пересечения.
159. Вычисление моментов инерции. Как было показано в предыдущем параграфе, мы сумеем вычислить момент инершии данной системы материаль частиц относигельно данной оси, если нам известны глав-
ные центральные моменты инериии и направления главных центральных осей инерции. В некоторых случаях направления главных осей могут быть указаны непосредственно.
a) Если массы расположены симметрично относительно некоторой плоскосги (или лежат в одной плоскости), то для полюсов на этой плоскости одна из главних осей инерции направлена по перпендикуляру к плоскости симметрии; иначе говоря, всякая линия, перпендикулярная к плоскости симметрии, будет главной осью. Действительно, возьмём начало координат в произвольной точке плоскости симметрии, а эту плоскость примем за плоскость $O x y$; тогда в суммах
\[
J_{v z}=\sum_{v=1}^{n} m_{v} y_{v} z_{v}, \quad J_{z x}=\sum_{v=1}^{n} m_{v} z_{v} x_{v}
\]

каждой массе $m_{v}$ с координатами $x_{v}, y_{v}, z_{v}$ будет соответствовать, по условию, равная масса $m_{\mu}=m_{v}$ с координатами $x_{v}, y_{v}-z_{v}$, за исключением масс $m_{y}$, лежащих в самой плоскости симметрии; но для таких масс $z_{
u}=0$. Следовательно, в нашем случае
\[
J_{y z}=J_{z x}=0
\]

а потому по формуле (26.10) уравнение эллипсоида инерции для начала координат (которое помещено в произвольной точке плоскости симметрии) будет
\[
J_{x x} x^{2}+J_{y y} y^{2}+J_{z z} z^{2}-2 J_{x y} x y=l^{2},
\]
т. е. ось $O z$ служит осью этого эллипсоида или, что то же, главной осью инерции.
б) Если массы имеют три взаимноперпендикулярные плоскости симметрии, то линии встречи этих пілоскостей служат главными центральными осями инерции [ср. свойство 5), § 146].

Когда пользуются декартовой системой координат, то большей частью вычислять непосредственно приходится лишь моменты инерции относительно координатных плоскостей, а затем уже с помощью формул (26.51) получают выражения для моментов инерции относительно осей координат.

Если плотность системы масс постоянна, то ее выносят за знак интеграла и вычисляют сперва оставшееся выражение, или так называемый геометрический момент инерции (момент инерции объёма, поверхности или линии); умножив полученный результат на плотность, получают требуемый момент инерции: в отличие от геометрического его называют физическим.

Пример 69. Геометрический момент инерции площади прямоугольного треугольника относительно одного из катетов. Примем вершину $A$ за начало координат, а ось $x$ направим по катету $A C$ (фиг. 103). Тогда, если длины катетов $B C$ и $A C$ — соответственно обовначим $a$ и $b$, то уравнение прямой $A B$ будет
\[
y=\frac{a}{b} x,
\]

а потому имеем
\[
J_{x x}=\int_{0}^{b} d x \int_{0}^{y} y^{2} d y=\frac{1}{3} \int_{0}^{b} y^{3} d x=\frac{a^{3}}{3 b^{3}} \int_{0}^{b} x^{3} d x=\frac{a^{3} b}{12},
\]

или
\[
j_{x x}=\frac{S a^{2}}{6},
\]

где $S$ — площадь треугольника.
Пример 70. Геометрический момент инерции площади косоугольного треугольника относительно одной из сторон. Пусть высота треугольника $A B C$ равна $C C^{\prime}=h$ (фиг. 104 и 105 ); тогда, рассматривая данный треугольник как сумму (фиг. 104) или разность (фиг. 105) двух

Фиг. 103.
Фиг. 104.
Фиг. 105.

прямоугольных треугольников $A C C^{\prime}$ и $B C C^{\prime}$, легко найдём, если $S$ — площадь данного түеугольника, что искомый момент инерции относительно основания $A B$ равняется $\frac{S h^{2}}{6}$.

Пример 71. Главные центральные геометрические моменты ииердии плошади элдипса, лежащегов плоскости $O x y$ и заданого уравнением
\[
\frac{x^{2}}{a^{2}}+\frac{y^{2}}{b^{2}}=1 \text {. }
\]

Вычисляем:
\[
J_{x x}=4 \int_{0}^{b} y^{2} d y \int_{0}^{2} d x=\frac{4 a}{b} \int_{0}^{b} y^{2} \sqrt{b^{2}-y^{2}} d y .
\]

Вводим новую переменную $\varphi$, положив $y=b \sin \varphi$; тогда
\[
\begin{aligned}
J_{x x} & =4 a b^{2} \int_{0}^{\frac{\pi}{2}} \sin ^{2} \varphi \cos ^{2} \varphi \cdot d \varphi=a b^{3} \int_{0}^{\frac{\pi}{2}} \sin ^{2} 2 \varphi \cdot d \varphi= \\
& =\frac{a b^{3}}{2} \int_{0}^{\frac{\pi}{2}}(1-\cos 4 \varphi) d \varphi=\frac{\pi a b^{8}}{4},
\end{aligned}
\]

или, если $S$ — площадь эллипса,
\[
J_{x x}=\frac{S b^{2}}{4} .
\]

Подобным образом найдём:
\[
J_{y y}=\frac{S a^{2}}{4} .
\]

Следовательно, по формуле (26.52)
\[
J_{z z}=\frac{S\left(a^{2}+b^{2}\right)}{4} .
\]

Когда $a=b=R$, т. е. эллипс превращается в круг радиуса $R$, то
\[
J_{x x}=J_{y y}=\frac{S R^{2}}{4}, \quad J_{z z}=\frac{S R^{2}}{2} .
\]

Пример 72. Главные центральные геометрическиемоменты инерции прямоугольного параллелепипеда, стороны которогосоответственно равны $2 a, 2 b, 2 c$. Направяя оси, $x, y, z$ параллельно сторонам $2 a, 2 b, 2 c$ и называя $Q_{x}$ площадь сечения, перпендикулярного оси $x$, и $V$ — объе п параллелепипеда, имеем
\[
J_{\text {пл. } y z}=2 \int_{0}^{a} Q_{x} x^{2} d x=8 b c \int_{0}^{a} x^{2} d x=\frac{8}{3} a^{3} b c=\frac{V a^{2}}{3},
\]

аналогично
\[
J_{\text {п. } z x}=\frac{V b^{2}}{3}, \quad J_{\text {пл. } x y}=\frac{V c^{2}}{3} .
\]

Огсюда по формуле (26.51) находим:
\[
J_{x x}=\frac{V\left(b^{2}+c^{2}\right)}{3}, \quad J_{y y}=\frac{V\left(c^{2}+a^{2}\right)}{3}, \quad J_{z z}=\frac{V\left(a^{2}+b^{2}\right)}{3} .
\]

Для куба $a=b=c$, и, следовательно,
\[
J_{x x}=J_{y y}=J_{z z}=\frac{2 V a^{2}}{3} .
\]

Чтобы получить геометрический момент инерции стержня длины $2 a$ относительно оси $z$, ему перпендикулярной и проходящей через его середину, можно взять третью из формул (26.53) и положить в ней $b=0$, а вместо $V$ поставить $2 a$; тогда мы получим:
\[
J_{z z}=\frac{2 a^{3}}{3} \text {. }
\]

Пример 73. Главные центральные геометрические моменты инерции круговогоцилиндра радиуса $R$ и высоты $H$. Поместим начало координат в центре цилиндра и ось $\boldsymbol{z}$ направим по его оси. Так как элемент объема в цилиндрических координатах равен $d \rho \cdot \rho d \varphi \cdot d r z$, то
\[
J_{x z}=\int_{-\frac{H}{2}}^{\frac{H}{2}} d z \int_{\gamma}^{2 \pi} d \varphi \int_{\gamma}^{R} \rho^{\mathrm{s}} \cdot d \rho=\frac{\pi R^{4} H}{2}=\frac{V R^{2}}{2},
\]

где $V$-объем цилиндра. Далее
\[
J_{x x}=J_{\text {пл. } x y}+J_{\text {пл. } z x}, \quad J_{y y}=J_{\text {пл. } x y}+J_{\text {пл. } y z} ;
\]

отсюда, так как по соображениям симметрии $J_{x x}=J_{y y}$, имеем
\[
J_{x x}=\frac{J_{x x}+J_{y y}}{2}=J_{\text {пл. } x y}+\frac{J_{\text {пл. } z x}+J_{\text {пл. } y z}}{2}=J_{\text {пл. } x y}+\frac{J_{z x}}{2} ;
\]

но
\[
J_{\text {เл. } x y}=\int_{-\frac{H}{2}}^{\frac{H}{2}} d z \int_{0}^{2 \pi} d \varphi \int_{0}^{R} z^{2} \cdot d \rho=\frac{V H^{2}}{12} ;
\]

окончательно получаем:
\[
J_{x x}=J_{y v}=\frac{V}{12}\left(H^{2}+6 R^{2}\right) .
\]

Пример 74. Главные центральные геометрические моменты инерции эллипсоида
\[
\frac{x^{2}}{a^{2}}+\frac{y^{2}}{b^{2}}+\frac{z^{2}}{c^{2}}=1 .
\]

Вычисляем сперва момент инериии относительно плоскости $O y z$; назвав $Q_{x}$ плошадь сечения, параллельного плоскости $O y z$ и отстоящего от неё на $x$, и $V-$ объём эллипсоида, находим:
\[
J_{\text {пл. } y z}=2 \int_{0}^{a} Q_{x} x^{2} d x
\]

по формуле (25.13) на стр. 251 имеем
\[
Q_{x}=\pi b c \cdot\left(1-\frac{x^{2}}{a^{2}}\right) ;
\]

подставив это выражение в предыдущую формулу, получаем:
\[
J_{\text {ил. } \dot{y} z}=\frac{4 \pi a^{3} b c}{15}=\frac{V a^{2}}{5} .
\]

Подобным образом находим:
\[
J_{\text {пл. } z x}=\frac{V b^{2}}{5}, \quad J_{\text {п.. } x y}=\frac{V c^{2}}{5} .
\]

Следовательно,
\[
J_{x x}=\frac{V\left(b^{2}+c^{2}\right)}{5}, \quad J_{v y}=\frac{V\left(c^{2}+a^{2}\right)}{5}, \quad J_{z z}=\frac{V\left(a^{2}+b^{2}\right)}{5} .
\]

Если $a=b=c=R$, т. е. эллипсоид обращается в шар радиуса $R$, то
\[
J_{x x}=J_{v y}=J_{z z}=\frac{2 V R^{2}}{5} .
\]

Пример 75. Дан однородный параллелепипед со сторонами $2 a, 2 a, 2 b$ (фиг. 106); исходя из результатов примера 72 , определить для середины $K$ ребра $A B$ главные геометрические моменты инерции и направления главных осей инерции. Назвав $V$ объём параллелепипеда, имеем
\[
J_{x x}=J_{y y}=\frac{V\left(a^{2}+b^{2}\right)}{3}, \quad J_{z z}=\frac{2 V a^{3}}{3} ;
\]
Фиг. 106.
далее, по формулам (26.40) и (26.41), заменив массу $M$ объёмом $V$, находим:
\[
\begin{array}{l}
J_{x x}^{(k)}=J_{y y}^{(k)}=\frac{V\left(4 a^{2}+b^{2}\right)}{3}, \quad j_{z z}^{(k)}=\frac{4 V a^{2}}{3}, \\
J_{y z}^{(k)}=J_{z x}^{(k)}=0 ; \quad J_{x y}^{(k)}=V a^{2} .
\end{array}
\]

Составляем уравнения (26.43), определяюшие направляющие косинусы $a_{v}, \beta_{v}, \dot{\gamma}$ главных осей инерции полюса $K$; положив для краткости письма $\delta_{v}=\frac{V_{\tau_{v}}}{3}$, имеем
\[
\left.\begin{array}{rrr}
\left(4 a^{2}+b^{2}-\tau_{v}\right) a_{v}- & 0 a^{2} \beta_{v}+ & 0 \cdot \gamma_{v}=0 \\
-3 a^{2} \alpha_{v}+\left(4 a^{2}+b^{2}-\tau_{v}\right) \beta_{
u}+ & 0 \cdot \gamma_{v}=0 \\
0 \cdot \alpha_{v}+ & 0 \cdot \beta_{
u}+\left(2 a^{2}-\tau_{v}\right) \gamma_{v}=0
\end{array}\right\}
\]

Вековое уравнение после развертывания определителя примет вид

Отсюда находим:
\[
\begin{array}{c}
\left(4 a^{2}+b^{2}-\tau_{\downarrow}\right)^{2}\left(2 a^{2}-\tau_{\downarrow}\right)-9 a^{4}\left(2 a^{2}-\tau_{v}\right)=0 . \\
\tau_{1}=2 a^{2}, \quad \tau_{2}=7 a^{2}+b^{2}, \quad \tau_{3}=a^{2}+b^{2} .
\end{array}
\]

Следовательно, главные момеить инерции для полюса $K$ равны
\[
J_{1}=\frac{2 V a^{2}}{3}, \quad J_{2}=\frac{V\left(7 a^{2}+b^{2}\right)}{3}, \quad J_{3}=\frac{V\left(a^{2}+b^{2}\right)}{3} .
\]

Косинусы $a_{v}, \beta_{v}, \gamma_{v}$ находим из уравнений (26.54), последовательно полагая $y=1,2,3$. Если ранг определителя системы $n$ однородных уравнений с $n$ неизвестными равен $n-1$ (что имеет место в нашем случае), то, как известно, нетривиальные решения (т. е. не равные нулю одновременно) иаходим, написав, что неизвестные пропорциональны соответственным адъюнктам определителя системы. На этом основании мы находим:
\[
\frac{\alpha_{1}}{0}=\frac{\beta_{1}}{0}=\frac{\gamma_{1}}{1}, \quad \frac{\alpha_{2}}{1}=\frac{\beta_{2}}{-1}=\frac{\gamma_{2}}{0}, \quad \frac{\alpha_{8}}{1}=\frac{\beta_{B}}{1}=\frac{r_{3}}{0} .
\]

Если, кроме того, мы примем во внимание, что сумма квадратов косинусов равна единице, то окончательно получия:
\[
\begin{array}{lll}
a_{1}=0, & \beta_{1}=0, & \gamma_{1}=1 ; \\
a_{2}=\frac{\sqrt{2}}{2}, & \beta_{2}=\frac{-\sqrt{2}}{2}, & \gamma_{2}=0 ; \\
a_{3}=\frac{\sqrt{2}}{2}, & \beta_{3}=\frac{\sqrt{2}}{2}, & \gamma_{3}=0 .
\end{array}
\]

Разумеется, в данном примере направления главных осей инериии для полюса $K$ можно было бы также найги непосредственно из соображений симметрии.