Метод, изложенный в предыдущем параграфе, зависит существенно от приведения неголономных систем к голономным путем введения дополнительных сил, обусловленных кинематическими связями. Практически это часто наиболее удобно выполняется тем, что составляют уравнения движения для каждого тела системы в отдельности. Кроме того, часто бывает выгодно пользоваться системой координат, движущейся как относительно пространства, так и относительно тела. Поэтому мы хотим сейчас составить уравнения движения твердого тела относительно системы осей, имеющих начало в центре тяжести тела и движущихся произвольным образом в пространстве ${ }^{1}$.

Пусть $G$ – центр тяжести тела, а $G x y z$ – подвижная система осей. Обозначим через $u, v, w$ компоненты скорости центра тяжести относительно подвижных осей и через $\vartheta_{1}, \vartheta_{2}, \vartheta_{3}$ – компоненты угловой скорости системы Gxyz относительно этих же осей. Компоненты угловой скорости тела относительно подвижных осей мы обозначим через $\omega_{1}, \omega_{2}, \omega_{3}$. Согласно $\S 64$ центр тяжести $G$ движется как свободная материальная точка, в которой сконцентрирована вся масса $M$ тела и к которой приложены все действующие на тело внешние силы. (При этом следует учитывать все силы реакций за исключением внутренних реакций между отдельными точками тела.) Допустим, что эти внешние силы имеют относительно подвижных осей $G x y z$ компоненты $X, Y, Z$.

Центр тяжести $G$ имеет по оси $G x$ компонент скорости $u$ и компонент ускорения, равный $\dot{u}-v \vartheta_{3}+w \vartheta_{2}$ (§17). Поэтому имеет место уравнение:
\[
M\left(\dot{u}-v \vartheta_{3}+w \vartheta_{2}\right)=X,
\]

которое может быть переписано в виде:
\[
\frac{d}{d t}\left(\frac{\partial T}{\partial u}\right)-\vartheta_{3} \frac{\partial T}{\partial v}+\vartheta_{2} \frac{\partial T}{\partial w}=X,
\]

где $T$ – кинетическая энергия тела, выраженная через $u, v, w, \omega_{1}, \omega_{2}$, $\omega_{3}$. Аналогичные уравнения имеют место и для движения центра тяжести в направлении осей $G y$ и $G z$.

Рассмотрим относительное движение тела вокруг центра тяжести $G$, которое согласно $\S 64$ не зависит от движения $G$. Согласно $\S 62,63$ момент количества движения тела относительно $G x$ равен $\frac{\partial T}{\partial \omega_{1}}$. Поэтому приращение момента количества движения тела относительно неподвижной в пространстве оси, мгновенно совпадающей с $G x$, равно:
\[
\frac{d}{d t}\left(\frac{\partial T}{\partial \omega_{1}}\right)-\vartheta_{3} \frac{\partial T}{\partial \omega_{2}}+\vartheta_{2} \frac{\partial T}{\partial \omega_{3}} .
\]
${ }^{1}$ При пользовании этим методом оси координат выбирают обычно таким образом, чтобы моменты инерции и девиации относительно них оставались постоянными. Однако это условие несущественно.

Если $L, M, N$ суть моменты внешних сил относительно осей $G x y z$, то ( $\S 40$ ) имеет место уравнение:
\[
\frac{d}{d t}\left(\frac{\partial T}{\partial \omega_{1}}\right)-\vartheta_{3} \frac{\partial T}{\partial \omega_{2}}+\vartheta_{2} \frac{\partial T}{\partial \omega_{3}}=L
\]

и два аналогичных уравнения.
Следовательно, движение тела определяется шестью уравнениями:
\[
\begin{array}{l}
\frac{d}{d t}\left(\frac{\partial T}{\partial u}\right)-\vartheta_{3} \frac{\partial T}{\partial v}+\vartheta_{2} \frac{\partial T}{\partial w}=X, \frac{d}{d t}\left(\frac{\partial T}{\partial \omega_{1}}\right)-\vartheta_{3} \frac{\partial T}{\partial \omega_{2}}+\vartheta_{2} \frac{\partial T}{\partial \omega_{3}}=L, \\
\frac{d}{d t}\left(\frac{\partial T}{\partial v}\right)-\vartheta_{1} \frac{\partial T}{\partial w}+\vartheta_{3} \frac{\partial T}{\partial u}=Y, \frac{d}{d t}\left(\frac{\partial T}{\partial \omega_{2}}\right)-\vartheta_{1} \frac{\partial T}{\partial \omega_{3}}+\vartheta_{3} \frac{\partial T}{\partial \omega_{1}}=M \\
\frac{d}{d t}\left(\frac{\partial T}{\partial w}\right)-\vartheta_{2} \frac{\partial T}{\partial u}+\vartheta_{1} \frac{\partial T}{\partial v}=Z, \frac{d}{d t}\left(\frac{\partial T}{\partial \omega_{3}}\right)-\vartheta_{2} \frac{\partial T}{\partial \omega_{1}}+\vartheta_{1} \frac{\partial T}{\partial \omega_{2}}=N
\end{array}
\]

Нетрудно видеть, что эти уравнения представляют собой не что иное, как уравнения Лагранжа в квазикоординатах и поэтому они могут быть получены при помощи теоремы $\S 30$.
ЗАДАчА 1. Начало координат подвижной системы не закреплено неподвижно в теле, а имеет относительно подвижных осей компоненты скорости $u_{1}, u_{2}, u_{3}$. Величины $\vartheta_{1}, \vartheta_{2}, \vartheta_{3}$ суть компоненты угловой скорости подвижных осей относительно этих же осей, $v_{1}, v_{2}, v_{3}$ суть компоненты скорости относительно подвижных осей той точки тела, которая в данное мгновение совпадает с началом координат, и $\omega_{1}, \omega_{2}, \omega_{3}$ суть компоненты угловой скорости тела относительно подвижных же осей. Погазать, что уравнения движения могут быть написаны в виде:
\[
\begin{array}{c}
\frac{d}{d t}\left(\frac{\partial T}{\partial v_{1}}\right)-\vartheta_{3} \frac{\partial T}{\partial v_{2}}+\vartheta_{2} \frac{\partial T}{\partial v_{3}}=X, \\
\frac{d}{d t}\left(\frac{\partial T}{\partial v_{2}}\right)-\vartheta_{1} \frac{\partial T}{\partial v_{3}}+\vartheta_{3} \frac{\partial T}{\partial v_{1}}=Y, \\
\frac{d}{d t}\left(\frac{\partial T}{\partial v_{3}}\right)-\vartheta_{2} \frac{\partial T}{\partial v_{1}}+\vartheta_{1} \frac{\partial T}{\partial v_{2}}=Z, \\
\frac{d}{d t}\left(\frac{\partial T}{\partial \omega_{1}}\right)-u_{3} \frac{\partial T}{\partial v_{2}}+u_{2} \frac{\partial T}{\partial v_{3}}-\vartheta_{3} \frac{\partial T}{\partial \omega_{2}}+\vartheta_{2} \frac{\partial T}{\partial \omega_{3}}=L, \\
\frac{d}{d t}\left(\frac{\partial T}{\partial \omega_{2}}\right)-u_{1} \frac{\partial T}{\partial v_{3}}+u_{3} \frac{\partial T}{\partial v_{1}}-\vartheta_{1} \frac{\partial T}{\partial \omega_{3}}+\vartheta_{3} \frac{\partial T}{\partial \omega_{1}}=M, \\
\frac{d}{d t}\left(\frac{\partial T}{\partial \omega_{3}}\right)-u_{2} \frac{\partial T}{\partial v_{1}}+u_{1} \frac{\partial T}{\partial v_{2}}-\vartheta_{2} \frac{\partial T}{\partial \omega_{1}}+\vartheta_{1} \frac{\partial T}{\partial \omega_{2}}=N,
\end{array}
\]

где $X, Y, Z, L, M, N$ означают компоненты и моменты внешних сил относительно подвижных осей.