[‘] Функция $L^{*}$, которую автор называет приведенной функцией Лагранжа, называется также функцией Рауса.
Легко показать, что эта функция равна разности кинетических энергий $T^{\prime}$ и $T^{\prime \prime}$
\[
L^{*}=T^{\prime}-T^{\prime \prime},
\]
где
\[
T^{\prime}=\frac{1}{2} \sum_{i, k=1}^{m} a_{i k} \dot{q}_{i} q_{k} \quad \text { и } \quad T^{\prime \prime}=\frac{1}{2} \sum_{i, k=m+1}^{n} a_{i k} \dot{q}_{i} \dot{q}_{k} .
\]
Выразим явно частные производные $\frac{\partial T}{\partial \dot{q}}$ :
\[
\begin{array}{l}
\frac{\partial T}{\partial \dot{q}_{1}}=a_{11} \dot{q}_{1}+a_{1 s} \dot{q}_{2}+\ldots+a_{1 m} \dot{q}_{m}+a_{1 m+1} \dot{q}_{m+1}+\ldots+a_{1 n} \dot{q}_{n}, \\
\text {. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . } \\
\frac{\partial T}{\partial \dot{q}_{m}}=a_{m 1} \dot{q}_{1}+a_{m s} \dot{q}_{2}+\ldots+a_{m m} \dot{q}_{m}+a_{m m+1} \dot{q}_{m+1}+\ldots+a_{m n} \dot{q}_{n}, \\
\text {. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . } \\
\frac{\partial T}{\partial \dot{q}_{n}}=a_{n 1} \dot{q}_{1}+a_{n 2} \dot{q}_{2}+\ldots+a_{n n} \dot{q}_{m}+a_{n m+1} \dot{q}_{m+1}+\ldots+a_{n n} \dot{q}_{n} . \\
\end{array}
\]
Далее имеем
\[
\begin{aligned}
*=\frac{1}{2} \sum_{i=1}^{h} \frac{\partial T}{\partial \dot{q}_{i}} \dot{q}_{i}-\sum_{i=1}^{m} \frac{\partial T}{\partial \dot{q}_{i}} \dot{q}_{i}= & \\
& =\frac{1}{2} \sum_{i=m+1}^{n} \frac{\partial T}{\partial \dot{q}_{i}} \dot{q}_{i}-\frac{1}{2} \sum_{i=1}^{m} \frac{\partial T}{\partial \dot{q}_{i}} \dot{q}_{i} .
\end{aligned}
\]
Умножая $\frac{\partial T}{\partial \dot{q}}$ на соответствующие значения $\dot{q}$ и вычисляя суммы, найдем, что обе суммы будут содержать в себе одинаковые слагаемые, которые и исчезнут после приведения подобных членов:
\[
\left\{=\frac{1}{2} \sum_{i, k=m+1}^{n} a_{i k} \dot{q}_{i} \dot{q}_{k}-\frac{1}{2} \sum_{i, k=1}^{m} a_{i k} \dot{q}_{i} \dot{q}_{k} .\right.
\]
Если координаты $\vec{q}_{1} \ldots \bar{q}_{m}$ – циклические, то
\[
\dot{q}_{i}=\gamma_{i}+\lambda_{i} \quad(i=1,2, \ldots, m) .
\]
Заменяя в последней сумме $\dot{q}_{i}$ их значениями, найдем
\[
2^{*}=\frac{1}{2} \sum_{i, k=m+1}^{n} a_{i k} \dot{q}_{i} \dot{q}_{k}-\frac{1}{2} \sum_{k, i=1}^{m} a_{i k}\left(\gamma_{i}+\lambda_{i}\right)\left(\gamma_{k}+\dot{\lambda}_{k}\right) .
\]
[] Утверждение автора о том, что при связях, не зависящих от времени, коэффициенты $b_{j 0} \equiv 0$, в такой категорической форме ничем не оправдано и в некоторых случаях может оказаться неверным. Можно привести примеры неголономных систем, уравнения связей которых не содержат явно времени, а все коэффициенты $b_{j h}$ либо постоянны, либо являются функциями координат, причем $b_{j 0}
eq 0$.
В качестве простейшего примера можно привести трехколесный велосипед (трехколесное шасси самолета) (фиг. 22). Напишем для его переднего
Фиг. 22.
колеса уравнение связей, полагая, что и задние, и передние колеса катятся по плоскости без скольжения. Пусть $\dot{\psi}$ есть угловая скорость рамы $O A, \theta-$ абсолютная угловая скорость вилки $A C$ около вертикальной оси, а $\omega=\dot{\varphi}-$ относительная угловая скорость переднего колеса по отношению к вилке $A C$. Обозначим через $v_{0}$ скорость средней точки $O$ задней оси, направленную во все время движения по прямой $O A$, если скольжение отсутствует. Введем также единичные векторы $\boldsymbol{i}, \boldsymbol{j}, \boldsymbol{k}$, направленные соответственно по прямой $O A$, по оси переднего колеса и по нормали к плоскости, по которой катится велосипед.
Можно написать для скорости точки $C$ следующие два равенства:
\[
\begin{array}{l}
\boldsymbol{v}_{C}=\boldsymbol{\theta} \times \boldsymbol{r}, \\
\boldsymbol{v}_{C}=\boldsymbol{v}_{\boldsymbol{0}}+\dot{\psi} \boldsymbol{k} \times \overrightarrow{O A}+\dot{\theta}_{\boldsymbol{k}} \times \overrightarrow{A C},
\end{array}
\]
где $r$ есть радиус-вектор, идущий из точки касания переднего колеса с плоскостью в центр колеса. Исключая $v_{C}$ и проектируя полученное равенство на ось $O A$ и на перпендикулярное к ней направление, получим два уравнения:
\[
\begin{array}{l}
\omega r \cos \alpha=v_{0}+\dot{\theta} \lambda \sin \alpha, \\
\omega r \sin \alpha=\dot{\omega} l-\dot{\theta} \lambda \cos \alpha ;
\end{array}
\]
здесь можно положить $\omega=\dot{\varphi}$ и $\dot{\theta}=\dot{\phi}+\dot{\alpha}$. Пока мы рассматриваем уравнения в написанной здесь форме, утверждение авторов остается в силе.
Представим себе теперь, что угловые скорости $\alpha, \dot{\theta}, \dot{\psi}$ малы по сравнению с $\omega$, а угол $\alpha$ также есть малая величина; пусть, кроме того, скорость $v_{0}$ поддерживается постоянной. Тогда с точностью до малых второго порядка можно написать
\[
\omega r=v_{0}, \quad \dot{\psi}(l-\lambda)-\dot{\alpha} \lambda=v_{0} \cdot \alpha .
\]
Второе из равенств представляет собой уравнение кинематической связи между переменным $\alpha, \dot{\alpha}$ и $\dot{\psi}$; здесь коэффициенты при $\dot{\psi}$ и $\dot{\alpha}$ не зависят от времени, а правая часть равенства все же не равна нулю тождественно.
[3] Уравнения движения неголономных систем обычно пишут в иной форме, пользуясь неопределенными множителями (см., например, Суслов, Теоретическая механика). Пусть общее уравнение механики приведено к виду (79)
\[
\sum_{h=1}^{n} \tau_{h} \hat{\delta} q_{h}=\sum_{h=1}^{h} Q_{h} \hat{\delta} q_{h}
\]
а уравнения кинематических связей имеют вид
\[
\begin{array}{l}
\text {. . . . . . . . . . . . . . . . } \\
a_{m 1} \dot{q}_{1}+a_{m 2} \dot{q}_{2}+\ldots+a_{m n} \dot{q}_{n}=a_{m} . \\
\end{array}
\]
Тогда соотношения между вариациями $\hat{\partial} q_{h}$ будут
\[
a_{r 1} \hat{\imath} q_{1}+a_{r 2} \hat{\jmath} q_{2}+\ldots+a_{r n^{\hat{\jmath}}} \hat{\jmath} q_{n}=0 \quad(r=1,2, \ldots, m) .
\]
Умножая последиие уравнения соответственно на неопределенные множители $\lambda_{1}, \lambda_{2}, \ldots, \lambda_{m}$ и вычитая сумму из (79), получим
\[
\sum_{h=1}^{n}\left(\tau_{h}-Q_{h}-\sum_{r=1}^{m} a_{r h} \lambda_{r}\right) \hat{o} q_{h}=0 .
\]
Выбирая теперь $\lambda_{1}, \ldots, \lambda_{m}$ так, чтобы множители при $\delta q_{1}, \ldots, \delta q_{m}$ равнялись нулю, найдем, что множители и при остальных $n-m$ вариациях $\hat{\delta} q_{m+1}, \ldots, \hat{\delta} q_{n}$ также должны равняться нулю, так как эти вариации совершенно произвольны.
Таким образом, получаем следующие $n$ уравнений движения:
\[
\tau_{h}=Q_{h}+\sum_{r=1}^{m} a_{r h \lambda_{r}} \quad(h=1,2, \ldots, n) ;
\]
вместе с $m$ уравнениями кинематических связей они дают систему $(n+m)$ уравнений с $(n+m)$ неизвестными $q_{1}, \ldots, q_{n}, \lambda_{-1}, \ldots, \lambda_{m}$. Уравнения, которые автор называет уравнениями Маджи, можно получить, исключая из полученной системы уравнений множите.и $\lambda_{1}, \ldots, \lambda_{m}$. Результат можно представить в сжатой форме, если приравнять нулю определители порядка $(m+1)$ матрицы
\[
\left\|\begin{array}{llr}
\tau_{1}-Q_{1}, & \tau_{2}-Q_{2}, \ldots, \tau_{n}-Q_{n} \\
a_{11}, & a_{12}, & \ldots \\
\cdots & \ldots & a_{1 n} \\
a_{m 1}, & a_{m 2}, & a_{m n}
\end{array}\right\| .
\]
Число независимых уравнений будет равно $(n-m)$.
[4] В истории науки можно привести много примеров, когда одновременно над одним и тем же вопросом работали ученые разных стран и приходили к одним и тем же результатам независимо один от другого. Впервые С. А. Чаплыгин в сообщении на заседании физического отделения Общества любителей естествознания 25 октября 1895 г. изложил метод получения уравнений движения неголономных систем без неопределенных множителей.
Уравнения движения были приведены им как раз к той форме, о которой упоминают авторы в тексте, т. е. С. А. Чаплыгин выделил члены неголономности, хотя и не придал этому особого значения.
Это сообщение с дополнениями появилось в печати в 1897 г. в Трудах Отделения физических наук Общества любителей естествознания, т. IX. Работы С. А. Чаплыгина по неголономным системам собраны и изданы в небольшом сборнике \”Исследования по динамике неголономных систем“, 1949 г.
Уравнения Чаплыгина легко получить из матрицы, приведенной в примечании на стр. 423 . Пусть уравнения связей приведены к виду
\[
\begin{array}{l}
\dot{q}_{1}=a_{1}, m+1 \dot{q}_{m+1}+\ldots+a_{1 n} \dot{q}_{n}, \\
\dot{q}_{m}=a_{m}, \cdots \dot{q}_{m+1}+\cdots+\dot{a}_{m n} \dot{q}_{n} .
\end{array}
\]
Тогда только что упомянутая матрица, из которой получаются все уравнения движения, упростится и ее можно написать в следующей форме:
Приравняв нулю первый из определителей ( $m+1$ )-го порядка, получим
Предположим теперь, что первые $m$ координат $q_{1}, \ldots, q_{m}$ не входят ни в выражение кинетической энергии, ни в коэффициенты $a_{r s}$; ограничиваясь рассмотрением консервативной системы, допустим, что силовая функция $U$ также не зависит от $q_{1}, q_{2}, \ldots, q_{m}$.
Так как
\[
\tau=\frac{d}{d t} \frac{\partial T}{\partial \dot{q}}-\frac{\partial T}{\partial q}, \frac{\partial \dot{q}_{k}}{\partial \dot{q}_{m+i}}=a_{k}, m+1,
\]
то уравнению (а) можно придать следующий вид:
\[
\begin{aligned}
\frac{d}{d t}\left(\frac{\partial T}{\partial \dot{q}_{m+1}}+\sum_{k=1}^{m} \frac{\partial T}{\partial \dot{q}_{k}} \frac{\partial \dot{q}_{k}}{\partial \dot{q}_{m+1}}\right)-\sum \frac{\partial T}{\partial \dot{q}_{k}} \frac{d a_{k}, m+1}{d t}- \\
-\frac{\partial T}{\partial q_{m+1}}=\frac{\partial U}{\partial q_{m+1}} .
\end{aligned}
\]
Обозначая через $\bar{T}$ выражение кинетической энергии, получаемое после замены $\dot{q}_{1}, \ldots, \dot{q}_{m}$ через $\dot{q}_{m+1}, \ldots, \dot{q}_{n}$ и замечая, что
\[
\frac{\partial \bar{T}}{\partial q_{m+1}}=\frac{\partial T}{\partial q_{m+1}}+\sum \frac{\partial T}{\partial \dot{q}_{k}} \frac{\partial \dot{q}_{k}}{\partial q_{m+1}},
\]
найдем
\[
\begin{array}{c}
\frac{d}{d t}\left(\frac{\partial \bar{T}}{\partial \dot{q}_{m+1}}\right)-\frac{\partial \bar{T}}{\partial q_{m+1}}+ \\
+\sum_{k=1}^{m} \frac{\partial T}{\partial \dot{q}_{k}}\left\{\sum_{s=1}^{n-m}\left(\frac{\partial a_{k}, m+s}{\partial q_{m+1}}-\frac{\partial a_{k}, m+1}{\partial q_{m+s}}\right) \dot{q}_{m+s}\right\}=\frac{\partial U}{\partial q_{m+1}} .
\end{array}
\]
Такой же вид будут иметь и остальные уравнения движения. Заметим, что переменные $q_{1}, \ldots, q_{m}$ не входят в эти уравнения по предположению, а обобщенные скорости $\dot{q}_{1}, \ldots, \dot{q}_{m}$ могут быть исключены, т. е. система приведена таким образом к наименьшему возможному числу уравнений. Уравнения движения в форме (б) и были получены С. А. Чаплыгиным.
[] Автор имеет в виду в данном случае не только механические явления в узком смысле слова; методы изучения устойчивости одинаковы как для механических явлений, так и для процессов в более широком понимании. В последнем случае выбирается некоторая величина или несколько величин, изменение которых с течением времени характеризует изучаемый процесс. Эти величины с таким же правом могут быть приняты за побобщенные координаты“, как углы, прямолинейные отрезки или дуги кривых служат обобщенными координатами какой-либо механической системы. Например, говоря о колебаниях электрических или электромеханических, принимают за одну из координат количество электричества $q$, протекающего сквозь сечение проводника; в этом случае $i=d q / d t$ есть обобщенная скорость или сила тока. Для таких систем справедливы уравнения Лагранжа, и потому общие выводы, которые получены для механических явлений, будут справедливы и по отношению к процессам в более широком смысле слова.
Очень хорошее изложение вопроса о применении уравнений Лагранжа к электрическим системам можно найти в книге В. Ф. Миткевича, физические основы электротехники“, 1928.
[‘] Эта теорема есть частный случай первой теоремы А. М. Ляпунова об устойчивости. Для доказательства ее необходимо привлечь рассуждения, примененные Ляпуновым при изложении им второго метода. См. А. М. Ляпу нов, Общая задача об устойчивости движения, 1950, стр. 77 и сл.
[‘] К рассуждениям автора необходимы пояснения. По первой из теорем Ляпунова об устойчивости (см. [‘]) функция $H(x \mid t)$ должна быть знакоопределенной и, кроме того, должно быть
\[
H(x, t)=\text { const. }
\]
Если $H(x \mid t)$ – знакоопределенная функция, то при $x_{s}=0$ и произвольном $t$ она равна нулю, а при всех остальных значениях $x_{s}$ сохраняет один и тот же знак, пока $\left|x_{s}\right| \leqslant h$. Отсюда и следует, что функция $H(x, t)$ для $x_{s}=0$ имеет максимум или минимум. В случае Дирихле можно положить
\[
H=T+V=\text { const, }
\]
т. е. $H$ есть полная энергия системы. Так как $T+V$ есть знакоопределенная функция, а $T \supseteqq 0$, то и
\[
H=T+V=\text { const }>0 .
\]
Но в силу того, что $V$ не зависит от $\dot{q}_{s}$, должно быть
\[
V>0 \text {, }
\]
т. е. потенциальная энергия в положении равновесия должна иметь минимум.
Отсюда следует, что теорема Дирлхле есть просто частный случай первой теоремы Ляпунова, а потому нельзя и противопоставлять их друг другу
[‘] С этими рассуждениями авторов нельзя согласиться, так как устойчивость решений системы (16′) еще не влечет за собой устойчивости полной системы (16’) и (16\”).
Раус рассматривал вопрос не вполне так, как излагают авторы, а основная теорема Рауса об устойчивости движения голономной консервативной системы есть частный случай теоремы Ляпунова об устойчивости движения.
Разделение понятия устойчивости на пустойчивость по Дирихле“ и „устойчивость по Раусу“ ничем исторически не оправдано, так как и Дирихле и Раус не давали точного определения этого понятия. Впервые Н. Е. Жуковский обратил внимание на то, что задачу об устойчивости движения консервативной системы можно ставить иначе, чем это сделано у Рауса, и только А. М. Ляпунов дал окончательное, общепринятое теперь определение понятия об устойчивости движения.
[] На случай кратных корней характеристического уравнения впервые обратил внимание Лагранж в своей „Аналитической механике“. Излагая теорию колебаний, он указывает (Ж. .Т агран ж, Аналитическая механика, т. I, 1950, Динамика, отдел шестой, § 1, п. 7, стр. 452 и сл.), что координаты системы будут оставаться малыми только при условии, что все корни характеристического уравнения вещественны, положительны и неравны между собой. Нетрудно показать на частном примере, что последнее требование излишне. Рассмотрим колебание материальной точки, лежащей в горизонтальной плоскости и притягиваемой к неподвижному центру с силой, пропорциональной расстоянию. Уравнения движения точки будут
\[
\ddot{x}+k^{2} x=0, \ddot{y}+k^{2} y=0 .
\]
Характеристическое уравнение в данном случае имеет вид
\[
\left(-\mu^{2}+k^{2}\right)\left(-\mu^{2}+k^{2}\right)=0
\]
и допускает кратные корни
\[
\mu_{1}^{2}=k^{2}, \mu_{2}^{2}=k^{2} \text {. }
\]
Общее решение можно написать в форме
\[
x=A \cos (k t+\alpha), y=B \cos (k t+\beta),
\]
где $A, B, \alpha, \beta$ – произвольные постоянные; здесь время $t$ входит только под знаком синуса или косинуса и координаты $x$ и $y$ остаются малыми, если малы $A$ и $B$, определённые через начальные данные. Лагранж считал, что в случае кратных корней время $t$ войдет в общее решение не только под знаком синуса и косинуса, а потому координаты $x$ и $y$ будут неограниченно возрастать, т. е. для устойчивости равновесия необходимо, чтобы все корни характеристического уравнения были различными. Ошибка Лагранжа была предметом специальных исследований Вейерштрасса (Gesammelte Werke, т. I). Из.тожение этого вопроса можно найти в книге академика Крылова А. Н. „О некоторых дифференциальных уравнениях математической физики, имеющих применение в технических вопросах“, изд. 2 (1932) или „Вибрация судов“ (Собрание сочинений, т. X, Изд. Академии паук СССР, (1948)).
[10] Научное наследство А. М. Ляпунова так велико, что, говоря о „теореме Ляпунова“, надо ее цитировать, так как число этих теорем огромно. К сожалению, автор не указывает точно, где опубликована приведенная теорема, но из теорем Ляпунова можно вывести совершенно определенное суждение о значении метода малых колебаний. Сошлемся на теоремы (I, II, III) А. М. Ляпунова, помещенные на стр. 82-94 его классического сочинения \”Общая задача об устойчивости движения “, М., 1950 г. На стр. 137 того же сочинения находим: „Из предыдущего анализа следует, что в большинстве случаев вопрос об устойчивости разрешается уже исследованием одного первого приближения, ибо только в тех случаях исследование это не дает ответа на вопрос, когда определяющее уравнение, не имея корней с положительными веществен-
ными частями, имеет корни, вещественные части которых суть нули“. Из этого следует, что указанный автором настоящей книги случай неустойчивости был исследован А. М. Ляпуновым за десять тет до того. Ляпунов указал и приемы, которым надо следовать при исследовании устойчивости \”в „особых случаях“, т. е. когда корни характеристического уравнения имеют действительные части, равные нулю.
[11] Ссылка автора на теорему Ляпунова ошибочна, а его точка зрения на значение метода малых колебаний при рассмотрении частных практических вопросов может ввести читателя в заблуждение. Метод малых колебаний приводит к исчерпывающему ответу, если все корни характеристического уравнения имеют действительные отрицательные части или в том случае, когда хотя бы один из них имеет положительную вещественную часть. Если же имеются корни, действительные части которых равны нулю, то нельзя судить об устойчивости и неустойчивости по первому приближению, так как все будет зависеть от членов более высокого порядка в уравнениях возмущенного движения. Если все кории чисто мнимые, то требуется дополнительное исследование. Обычно это встречается при исследовании устойчивости консервативных систем, но в этих случаях можно вывести необходимое заключение из анализа интеграла энергии. Если в рассмотрение входят диссинативные силы, что обычно и бывает при решении технических проблем, то можно и́стребовать, чтобы все корни характеристического уравнения имели отрицательные действительные части. В тех случаях, когда все же нельзя удовлетворить этому условию и когда входит, например, один нулевой корень, следует обратиться к исследованиям \”особых случаев“ Ляпунова или изменить постановку задачи, что иногда бывает возможно.
[12] Исследование устойчивости периодических решений впервые было выполнено А. М. Ляпуновым в его диссертации „Общая задача об устойчивости движения“ (гл. III).