Закономерности движения механических систем с односторонними связями

( Motion Lanes for Mechanic Systems with One-Sided Relations
Preprint, Inst. Appl. Math., the Russian Academy of Science)

Кугушев Е.И., Сорокина О.В.
(E.I.Kugushev, O.V.Sorokina)

ИПМ им. М.В.Келдыша РАН

Москва, 2002
Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (проекты №№ 01-01-00508, 02-01-00352 и 02-07-90027)

Аннотация

Рассматриваются механические системы с голономными односторонними связями. На основе обобщенных уравнений Лагранжа 1-го и 2-го рода формулируются теоремы об изменении количества движений системы и его момента, теоремы Аппеля, Рауса, Нётер.

Abstract

Mechanics systems with one-sided relations are studied. On the base of generalized Lagrange equations of motion the theorems of pulse and pulse moment, the theorems of Appel, Routh, Noether are proved.

 

Содержание

Введение.............................................................................................................. 3

1.   Принцип Даламбера-Лагранжа...................................................................... 4

2.   Теорема Аппеля........................................................................................... 11

3.   Теоремы Карно и непрерывность кинетической энергии.......................... 12

4.   Теорема об изменении количества движения............................................. 12

5.   Теорема об изменении момента количества движения.............................. 12

6.   Циклические интегралы и теорема Рауса.................................................... 12

7.   Преобразование координат......................................................................... 12

8.   Теорема Нётер.............................................................................................. 12

9.   Несколько примеров.................................................................................... 12

ЛИТЕРАТУРА................................................................................................... 12

 


Введение

 

          В работе рассматриваются механические системы с односторонними связями. Обзор работ по данной тематике можно найти в [1,2].

          В публикации [3] был рассмотрен обобщённый принцип Даламбера-Лагранжа и уравнения движения для систем с идеальными односторонними связями в форме уравнений Лагранжа первого и второго рода с импульсными правыми частями. В данной работе, на основе этого принципа, основные законы движения механических систем с идеальными односторонними связями  (законы изменения полного импульса системы и его момента) распространяются на системы с идеальными голономными односторонними связями. Формулируются в общем виде теоремы Аппеля и Карно. Помимо этого на системы с идеальными голономными односторонними связями распространяется действие теорем Рауса и Нётер о существовании и использовании циклических интегралов.

          Опишем кратко содержание работы. В первом разделе следуя [3] мы формулируем обобщённый принцип Даламбера-Лагранжа и уравнения движения для механических систем с односторонними связями. Это нам понадобится в дальнейшем для вывода основных закономерностей движения таких систем. В локальной формулировке принцип Даламбера-Лагранжа для систем с односторонними связями известен уже давно [4]. Однако, исходя из него затруднительно получать уравнения движения справедливые на всем отрезке движения, поскольку в его формулировке присутствуют ускорения точек системы, которые могут быть неопределенными в моменты выхода на границу односторонних связей (т.е. в точках удара). В данном разделе мы даём интегральную формулировку принципа Даламбера-Лагранжа, которая позволяет избавиться от этих трудностей.

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

В следующих двух разделах рассматриваются теоремы, которые для систем с идеальными удерживающими связями относятся к числу основных законов движения. В четвёртом разделе работы для систем с идеальными односторонними связями формулируются теоремы об изменении количества движения  системы и о движении ее центра масс. В пятом разделе для систем с идеальными односторонними связями формулируются теоремы об изменении момента количества движения.

Следует заметить, что известны формулировки этих законов для систем с односторонними связями [4]. Они выводятся из локального принципа Даламбера-Лагранжа. Однако, использующиеся при этом рассуждения справедливы только для тех участков траектории, где скорость движения, по крайней мере, непрерывна (т.е. можно вычислить ускорения). Поэтому для законов сохранения в точках разрыва скоростей не гарантируется сохранение констант соответствующих интегралов движения. Приводимые в данной работе доказательства основаны на обобщенных уравнениях Лагранжа первого рода и лишены указанных недостатков.

В шестом разделе показывается, что теория Рауса о понижении порядка лагранжевой системы при наличии циклических координат справедлива и для систем с идеальными односторонними связями. Известно, что уравнения Лагранжа второго рода  инвариантны относительно замены координат. В седьмом разделе мы показываем, что это остается справедливым и при наличии идеальных односторонних связей. В восьмом разделе мы показываем, что теорема Нётер о существовании первых интегралов лагранжевых систем справедлива и при наличии идеальных односторонних связей. Как и в предыдущих разделах, наши рассуждения почти дословно повторяют традиционные.

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

 

1.     Принцип Даламбера-Лагранжа.

 

В данном разделе, следуя [3] мы сформулируем обобщенный принцип Даламбера-Лагранжа для механических систем с односторонними связями. Это нам понадобится в дальнейшем для вывода основных закономерностей движения таких систем.

В качестве аппарата, пригодного для описания движения систем с односторонними связями на любых участках траектории, предлагается использовать теорию функций с ограниченной вариацией (или, говоря иначе, мер Лебега-Стилтьеса). Остановимся на мотивах именно такого выбора. В классической механике систем с идеальными удерживающими связями движение механической системы представляет собой достаточно гладкую кривую в конфигурационном пространстве. При наложении односторонних связей скорость движения может уже претерпевать скачки. В большинстве случаев её представляют  в виде гладкой функции, к которой добавляется кусочно-постоянная функция, называемая функцией скачков. Таким образом, реальное движение систем с односторонними связями можно было бы описывать в пространстве функций, которые представимы в виде такой суммы. Однако такое пространство не является удобным для исследования свойств движения аналитическими средствами. Например, в нем затруднительно решать вопрос о корректности систем обыкновенных дифференциальных уравнений, описывающих движение. Удачным пространством было бы, например, какое-либо банахово, т.е. полное метрическое пространство функций. В самом деле, когда (при исследовании систем с удерживающими связями) мы берем в качестве пространства возможных скоростей движения банахово пространство непрерывных функций с естественной метрикой, то оказывается, что решения дифференциальных уравнений (т.е. траектории) существуют и, даже, являются гладкими.

Здесь уместно поставить вопрос о том, каково минимальное пространство функций, которое объемлет гладкие функции и функции скачков и при этом является банаховым. Естественной метрикой в пространстве функций скачков является полная вариация. И минимальным подходящим для нас банаховым пространством с этой метрикой является пространство функций с ограниченной вариацией. Именно поэтому удобным представляется выбор пространства возможных траекторий систем с односторонними связями в виде пространства таких функций, производная (скорость) которых является функцией с ограниченной вариацией. Сами траектории при этом, по крайней мере, являются абсолютно непрерывными функциями.

При таком описании движения реакции связей также описываются как меры Лебега-Стилтьеса. Уравнения движения в этом случае представляют собой обыкновенные дифференциальные уравнения, содержащие импульсные члены в правой части.

Обсудим теперь кратко вопрос о том, как понимать обыкновенные дифференциальные уравнения с импульсной правой частью. Рассмотрим сначала уравнение

                                                                           

Запишем его немного иначе

                                                    (1.1)

где , , . Пусть  и    непрерывные  вектор-функции, а  – одномерная мера Лебега-Стилтьеса.  Естественно предположить, что решение этого уравнения мы ищем в классе вектор-функций  с ограниченной вариацией. В этом случае и в левой и в правой частях уравнения стоят векторные меры Лебега-Стилтьеса, и само уравнение надо понимать как их совпадение.

          Дадим теперь строгое определение. Пусть на отрезке  задана вектор-функция   с ограниченной вариацией. Она является решением указанного уравнения с начальными условиями , если  и для любой непрерывной вектор-функции , заданной на отрезке  выполнено

                                            

Это условие эквивалентно выполнению для каждого  следующего уравнения

                                                

Интеграл здесь понимается как интеграл Стилтьеса на отрезке .

          В данной работе нас будут интересовать несколько более сложные импульсные дифференциальные уравнения. В общем виде их можно записать следующим образом.

                                                   

где , ,, , . Пусть  и    непрерывные  вектор-функции, а  – одномерная мера Лебега-Стилтьеса. Оно рассматривается совершенно аналогично.

Дадим строгое определение решения. Пусть на отрезке  задана абсолютно непрерывная вектор-функция   и функция с ограниченной вариацией . Они являются решением приведенной выше системы уравнений с начальными условиями , если   и для любого  на отрезке  выполнено

                                

В частном случае, когда , эта система может быть записана более кратко следующим образом

Перейдем теперь собственно к принципу Даламбера-Лагранжа. Рассмотрим систему из  материальных точек, перемещающихся в пространстве под действием приложенных к ним сил. Координаты точек объединим в один вектор , где . В отсутствие связей движение системы описывается вторым законом Ньютона

                   ,                                            (1.1)

где  – массы точек, а  – вектор сил, действующих на точки системы.

          Наложим на систему семейство  линейных идеальных удерживающих связей:

,                              (1.2)

где  – число связей. Запишем (1.2) в матричной форме.

                                                                                  (1.3)

где  – матрица , .

Независимость связей. Удерживающие связи независимы, если  всюду.

Наложим на нашу систему семейство односторонних голономных связей

                   ,                                                   (1.4)

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

Внутренностью конфигурационного пространства будем называть множество точек , в которых

                   ,                                                   (1.5)

Границей конфигурационного пространства или поверхностью удара будем называть множество точек  удовлетворяющих хотя бы одному уравнению

                   ,                                                   (1.6)

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

          Возможным перемещением удерживающих связей в точке  называется любой вектор , удовлетворяющий условию

.                                                                  (1.7)

Эти возможные перемещения образуют линейное пространство   возможных перемещений удерживающих связей (размерности , если связи независимы).

Возможные перемещения. Пространство возможных перемещений  в точке  обозначается . Для точек, лежащих внутри области (1.4) оно  определяется как пространство  возможных перемещений удерживающих связей. Для точек, лежащих на границе области (1.4), мы требуем также, чтобы вектор возможного перемещения  касался этой границы, т.е. чтобы вдобавок к (1.7) выполнялось следующее условие:

           для тех , для которых    (1.8)

здесь  означает обычное скалярное произведение векторов  и .

          Конус допустимых направлений. Мы также будем предполагать, что в каждой точке поверхности удара можно определить конус  допустимых направлений. Он представляет собой замыкание конуса возможных перемещений, сдвигаясь вдоль которых мы попадаем внутрь области, допустимой связями (1.4). Формально конус допустимых направлений, в добавление к условиям (1.7), определяется системой неравенств

           для тех , для которых    (1.9)

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

Мы будем рассматривать движение на некотором отрезке времени .

          Вариацией кривой  будем называть непрерывную вектор-функцию .

Возможные вариации Вариацию  кривой  будем называть возможной, если  для всех , т.е. вектор  является возможным перемещением в точке . Отличие от определения возможной вариации для случая только удерживающих связей состоит в том, что во всех точках , в которых кривая  выходит на границу односторонних связей, вектор  касается границы односторонних связей.

Допустимые вариации. Вариацию  кривой  будем называть допустимой, если  для всех , т.е. вектор  является допустимым перемещением в точке . Отличие от определения возможной вариации  состоит в том, что во всех точках , в которых кривая  выходит на границу односторонних связей, вектор  лежит в конусе допустимых направлений  . Поскольку конус допустимых перемещений содержит пространство возможных перемещений, то любая возможная вариация является допустимой.

Уравнения движения мы будем строить исходя из следующих принципов.

Принцип освобождения от связей. Пусть  – траектория движения системы со связями (1.3, 1.4). Тогда  систему можно освободить от связей и добавить некую силу – реакцию связей  таким образом, что кривая  останется траекторией освобожденной системы. При этом компоненты реакции связей представляют собой меры Лебега-Стилтьеса , имеющие особенности, сосредоточенные на множестве тех моментов времени, в которые траектория  выходит на односторонние ограничения.

Траектория системы представляет собой такую дифференцируемую функцию , производная которой,  является функцией с ограниченной вариацией. При этом  будут справедливы следующие уравнения.

                   ,                          (1.10)

Идеальность связей (полная). Связи (1.3, 1.4) называются (вполне) идеальными, если для любой траектории системы ,  для любой ее допустимой вариации  интегральная элементарная работа сил реакции связей неотрицательна, т.е.

                                                     (1.11)

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

          Двусторонняя идеальность связей. Связи называются двусторонне идеальными, если для любой траектории движения системы  для любой ее возможной вариации  интегральная элементарная работа сил реакции односторонних связей равна нулю, т.е.

                                                     (1.12)

          Из полной идеальности связей вытекает их двусторонняя идеальность. В самом деле, допустим противное, что связи являются вполне, но не двусторонне идеальными. Тогда существует возможная (и, следовательно, допустимая) вариация  такая, что . Взяв возможную (и, следовательно, допустимую) вариацию  получим , что противоречит (1.11). Доказательство завершено.

Найдя  из (1.10) и подставив в (1.11) получим эквивалентную форму записи условия  идеальности связей. Для любой допустимой вариации  должно выполняться

                                                      (1.13)

Аналогично для условия двусторонней идеальности связей (1.12). Для любой возможной вариации  должно выполняться

                                   (1.14)

         

Принцип Даламбера-Лаграгжа. Пусть непрерывная кривая  удовлетворяет уравнениям связей (1.3, 1.4), а ее производная существует почти всюду и является функцией ограниченной вариации. Кривая   является траекторией движения системы с идеальными связями (1.3, 1.4) тогда и только тогда, когда для любой допустимой вариации  выполнены соотношения (1.13) и, следовательно, для любой возможной вариации, соотношения  (1.14).

          Исходя из этого принципа в [3] были получены уравнения Лагранжа первого рода.

Уравнения Лагранжа первого рода. Принцип Даламбера-Лагранжа эквивалентен следующим уравнениям Лагранжа первого рода. Пусть непрерывная кривая  удовлетворяет уравнениям связей (1.3, 1.4), а ее производная существует почти всюду и является функцией ограниченной вариации. Кривая  является траекторией движения системы с идеальными  связями (1.3, 1.4) тогда и только тогда, когда найдутся такие векторные меры  Лебега-Стилтьеса  и  , что

      (1.15)

При этом каждая мера  неотрицательна и сосредоточена на множестве моментов времени, в которые .

          Для безударных движений, когда траектории являются гладкими (скорости непрерывны по времени), эти уравнения совпадают с обычными уравнениями Лагранжа первого рода.

Матричная форма уравнений. Пусть  – диагональная матрица, с элементами  на диагонали,  – сводный вектор сил,  – транспонированная матрица  из (1.3),  – транспонированная матрица Якоби

                                                                       (1.16)

Тогда уравнения (1.15) можно записать в более компактной форме

                                                         (1.17)

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

               (1.18)

или, более кратко

                                                       (1.19)

Поскольку во все время движения выполнены уравнения удерживающих связей (1.3), то  и из (1.19) получаем

                                                                    (1.20)

Таким образом, скачки  и  не являются независимыми. Они связаны системой  линейных соотношений.

 

2.     Теорема Аппеля.

 

Для систем с удерживающими связями известна теорема Аппеля о сохранении обобщенных импульсов в касательной плоскости поверхности удара [5, 6]. Мы сформулируем ее, основываясь на уравнениях Лагранжа первого рода.

Кинетическая метрика. Кинетической метрикой будем называть Риманову метрику, порождаемую матрицей масс . В этой метрике кинетическая энергия системы  выражается формулой . Это половина квадрата нормы вектора скорости системы. Термин кинетическая метрика не является общепринятым, но мы будем им пользоваться для краткости.

Вектор импульса системы. Вектором импульса механической системы называется сводный вектор импульсов всех ее точек  , где . Кратко это можно записать следующим образом

                                                                                    (2.1)

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

Пусть система двигается по траектории  и    точка скачка вектора импульса. Будем называть ее точкой удара. Обозначим

                            (2.2)  

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

                                                            (2.3)

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

   и  ,                            (2.4)

для тех , для которых .

Не нарушая общности, для краткости записи, будем считать, что в точке удара все односторонние связи выходят на ограничение. Тогда плоскость удара  задается системой соотношений

                   ,                                       (2.5)

где  – матрица Якоби односторонних связей (1.16).        

Теорема Аппеля. Возьмем  – любой вектор в плоскости удара. Тогда из (2.3, 2.5) имеем . Т.е. сохраняются проекции вектора импульса системы на плоскость касательную поверхности удара. По-другому это равенство можно записать следующим образом

                                                                (2.6)

Оно означает, что в кинетической метрике сохраняются проекции скорости системы на плоскость удара.

          Введем вектора

Это градиенты односторонних связей в кинетической метрике.

          Нормальная плоскость удара. Плоскостью , нормальной к плоскости удара  будем называть линейное пространство, натянутое на вектора  в пересечении с линейным пространством  возможных перемещений удерживающих связей (1.7).

                                     (2.7)

          Нетрудно видеть, что  и  взаимно ортогональны в кинетической метрике и их линейная оболочка совпадает с . Таким образом,  есть прямая сумма касательной и нормальной плоскостей . В соответствии с этим скорости  и  однозначно раскладываются на касательную  и нормальную составляющие.

                                                        (2.8)

Причем

                                                                   (2.9)

для любого вектора  из касательной плоскости удара. Тогда из (2.6, 2.8) получаем

                                                                         (2.10)

т.е. теорема Аппеля означает, что сохраняется проекция скорости на касательную плоскость удара в кинетической метрике.

                   ,                                                   (2.11)

Отсюда следует, что

                                               (2.12)

Замена координат. Естественно, все указанные соотношения сохраняются при замене координат и, в частности, при интегрировании части связей и переходе к обобщенным координатам. Вектор импульса при этом определяется соотношением

 

3.     Теоремы Карно и непрерывность кинетической энергии.

 

В теоремах Карно (см., например,[6]) рассматривается поведение кинетической энергии в точках скачка скорости (точках удара). При этом используются следующие термины.

     Потерянные скорости – это компоненты вектора приращения скорости при ударе

                                                                              (3.1)

Из (2.11) следует, что  , т.е. что вектор потерянных скоростей лежит в нормальной плоскости удара.

          Кинетическая энергия потерянных скоростей – это кинетическая энергия системы, которую она имела бы, двигаясь с потерянной скоростью

                                                                                      (3.2)

Из (2.8, 2.11) следует, что

                                                       (3.3)

Потерянная кинетическая энергия определяется как разность кинетических энергий до и после удара

                                (3.4)

Из (2.8, 2.11) следует, что

                                                     (3.5)

          Первая теорема Карно обычно формулируется как явление, происходящее при внезапном наложении удерживающих связей. В нашей терминологии она звучит следующим образом. Если сразу после удара траектория касается границы удерживающих связей (т.е. ), то потерянная при ударе кинетическая энергия равна кинетической энергии потерянных скоростей (т.е. ). Это утверждение сразу вытекает из (3.3, 3.5).

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

Коэффициент восстановления. Нормальная составляющая скорости  не обязана сохраняться при ударе. Для определения ее изменения необходимы дополнительные предположения о природе связей наложенных на систему. Если связи не зависят от времени, одна из возможных моделей – это линейная взаимозависимость величин нормальной составляющей скорости до и после удара, задаваемая т.н. коэффициентом восстановления. Более подробную информацию об этом можно найти в [2].

          Коэффициентом восстановления  называется неотрицательная величина, определяемая следующим образом

                                                              (3.6)

Если , то  – это отношение модуля нормальной составляющей скорости  (до удара) к модулю нормальной составляющая скорости  (после удара) в кинетической метрике. Задание коэффициента восстановления позволяет определить  по , если удар не является кратным, т.е. мы вышли на ограничение ровно одной односторонней связи. В этом случае  представляет собой одномерное пространство.

Типы связей. Заметим, что, из физических соображений, при стационарности связей энергия в момент удара не должна возрастать. Поскольку активные силы не совершают мгновенной работы, то это фактически означает, что не возрастает кинетическая энергия системы. Отсюда следует, что . Если , то говорят, что в системе происходят абсолютно упругие соударения. Если , то говорят, что в системе происходят абсолютно неупругие соударения. Если же , то говорят, что в системе происходят просто неупругие соударения.

Третья теорема Карно. Рассмотрим теперь более подробно случай, когда удар не является кратным, т.е. траектория выходит на ограничение ровно одной односторонней связи. В этом случае  представляет собой одномерное пространство и поэтому.

                                                                                                (3.7)

В самом деле, если , или , то (3.7) следует из (3.6). Если  и , то из (3.6) следует, что, либо выполнено (3.7), либо . Но если бы было выполнено последнее, то это означало бы, что траектория находится за пределами области, определяемой односторонними связями – либо при подходе к точке удара, либо при отходе от нее, что невозможно.

Итак, выполнено соотношение (3.7).

          Формулировка теоремы. Третья теорема Карно формулируется следующим образом [6]. Потерянная кинетическая энергия не может быть отрицательной. Она равна доле кинетической энергии потерянных скоростей, вычисляемой по формуле

                                                                                    (3.8)

          Доказательство. Подставив (3.7) в (3.3) получим

                  

Подставив (3.7) в (3.5) получим

                  

Отсюда сразу следует (3.8). Поскольку , то и . Доказательство закончено.

Абсолютно упругий удар. Остановимся более подробно на случае абсолютно упругих соударений (). В этом случае  и из (2.8) следует, что кинетическая энергия сохраняется, несмотря на скачок скорости.

Приведенные рассуждения можно резюмировать в следующем утверждении.

Непрерывность кинетической энергии. В случае абсолютно упругих соударений кинетическая энергия системы является непрерывной функцией времени. Если связи гладкие (класса два), то, как было показано в [3], скорость является суммой абсолютно непрерывной составляющей и функции скачков (непрерывная сингулярная составляющая отсутствует). Тогда из непрерывности кинетической энергии следует ее абсолютная непрерывность.

Закон сохранения энергии. Если связи стационарны, и силы потенциальны, то непрерывность кинетической энергии немедленно влечет за собой закон сохранения энергии. В самом деле, пусть  – силовая функция. Тогда полная энергия  – это непрерывная функция. На каждом открытом отрезке, не содержащем скачков она постоянна в силу обычных законов сохранения для удерживающих связей (на участках, когда траектория выходит на границу односторонних ограничений, их можно рассматривать как удерживающие). Но поскольку множество скачков не более чем счетно, то полная энергия постоянна на всем отрезке времени.

 


4.     Теорема об изменении количества движения.

 

Для систем с идеальными удерживающими связями среди основных законов динамики формулируются теоремы об изменении количества движения  системы и о движении ее центра масс (см., например, [7]). Эти теоремы верны и при наличии идеальных односторонних связей.

Количество движения системы. Это трехмерный вектор  равный сумме количеств движения всех составляющих систему точек. Количество движения материальной точки - это трехмерный вектор , где  - масса точки, а  - трехмерный вектор ее скорости.

                                                                                       (4.1)

Выше мы использовали сводный вектор координат точек системы   и сводный вектор сил , . И в этих терминах были выписаны уравнения Лагранжа первого рода (1.17). Кроме того, мы ввели вектор импульса системы  (2.1).  В этих обозначениях  , где  - матрица , составленная из  экземпляров трехмерных единичных матриц . Мы будем использовать и то и другое описание, понимая, что они эквивалентны.

Теорема об изменении количества движения. Если удерживающие и односторонние связи, наложенные на систему, идеальны и допускают поступательный сдвиг всех точек системы как твердого тела вдоль какого-нибудь направления  постоянного во времени, то проекция количества движения системы на это направление  является абсолютно непрерывной функцией и скорость ее изменения равна , т.е. равна суммарной проекции на это направление вектора активных сил.

                                                                                     (4.2)

          Эта теорема непосредственно вытекает из аналогичного утверждения для сводных векторов системы.

          Теорема об изменении вектора импульса. Если удерживающие и односторонние связи, наложенные на систему, идеальны и допускают поступательный сдвиг всех точек системы как твердого тела вдоль какого-нибудь направления  постоянного во времени, то проекция вектора импульса на это направление  является абсолютно непрерывной функцией и скорость ее изменения равна , т.е. равна проекции на это направление сводного вектора активных сил

                                                                                      (4.3)

          Доказательство. Возьмем какую-либо траекторию движения системы , . Для нее выполнены уравнения Лагранжа первого рода (1.17). Используя вектор импульса системы  их можно переписать следующим образом

                                                            (4.4)

Условие теоремы означает, что в каждой точке конфигурационного пространства вектор  является возможным перемещением, т.е.  во все время движения и  в тех точках траектории, которые расположены на границе односторонних связей, т.е. в тех точках, в которых сосредоточена мера . Отсюда следует, что во все время движения  и . Заметим, что эти равенства, как и (4.4), надо понимать как равенство мер Лебега-Стилтьеса. Домножив обе части (4.4) на , получаем

                                                                                         (4.5)

Поскольку  есть функция ограниченной вариации, то  является абсолютно непрерывной функцией и выполнено (4.3). Доказательство закончено.

                Возьмем , т.е. составим вектор  из  экземпляров вектора . Если домножить теперь (4.5) слева на , то получим (4.3). Что доказывает теорему об изменении количества движения.

          Обычными рассуждениями из теоремы об изменении количества движения выводится следующие утверждения.

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

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

          Следствие. Пусть выполняются указанные выше предположения о связях. Если суммарная проекция всех активных сил на направление  равна нулю, то центр масс системы в проекции на эту ось двигается с постоянной скоростью.

 

5.     Теорема об изменении момента количества движения.

 

Для систем с идеальными удерживающими связями помимо теорем об изменении количества движения  формулируются теоремы об изменении момента количества движения (см., например, [7]). Эти теоремы верны и при наличии идеальных односторонних связей.

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

Момент количества движения системы относительно начала координат - это трехмерный вектор  равный сумме моментов количества движения всех составляющих систему точек. Момент количества движения материальной точки  относительно начала координат - это трехмерный вектор , где  - масса точки, а  и  - ее трехмерный радиус-вектор, и вектор скорости.

                                                                        (5.1)

Момент сил системы относительно начала координат - это трехмерный вектор  равный сумме моментов всех активных сил системы. Момент силы  относительно начала координат - это трехмерный вектор , где -  трехмерный радиус-вектор точки системы, к которой приложена сила.

                                                                            (5.2)

Так же, как и выше мы будем использовать сводные вектора координат точек системы  , сил , , импульса системы  и т.п.

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

Вектор моментов сил. Аналогично вводится сводный вектор  моментов сил системы относительно начала координат. Каждая тройка его координат равна тройке координат момента силы .

В этих терминах 

,    

где  - матрица , составленная из  экземпляров трехмерных единичных матриц .

Сводным векторным произведением векторов ,  будем называть вектор , определяемый по следующему правилу. Пусть , ,  - трехмерные вектора, составленные из троек координат векторов ,  и , как и выше. Тогда , где "" означает обычное векторное произведение. Сводное векторное произведение будем обозначать знаком "". Таким образом, .

Заметим, что для сводного векторного произведения справедливы те же утверждения, что и для обычного векторного произведения в трехмерном пространстве, в частности, верна формула смешанного произведения: .

Во введенных обозначениях можно записать

          ,                                                     (5.3)

          Теорема об изменении момента количества движения. Если удерживающие и односторонние связи, наложенные на систему, идеальны и допускают поворот всех точек системы как твердого тела вокруг какой-нибудь постоянной оси с направляющим вектором   и проходящей через начало координат, то проекция момента количества движения системы на эту ось  является абсолютно непрерывной функцией и скорость ее изменения равна , т.е. равна суммарной проекции на это направление векторов моментов активных сил.

                                                                                    (5.4)

     Эта теорема непосредственно вытекает из аналогичного утверждения для сводных векторов системы.

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

  Пусть у нас есть поворот в исходном трехмерном пространстве всех точек системы как твердого тела вокруг какой-нибудь постоянной оси с направляющим вектором   и проходящей через начало координат. Ему соответствует поворот в сводном пространстве координат точек системы  с направляющим вектором , координаты которого составлены из  экземпляров координат вектора .

          Теорема об изменении вектора моментов импульса. Если удерживающие и односторонние связи, наложенные на систему, идеальны и допускают поворот всех точек системы как твердого тела вокруг какой-нибудь постоянной оси с направляющим вектором   и проходящей через начало координат, то проекция сводного вектора моментов импульса системы на эту ось  является абсолютно непрерывной функцией и скорость ее изменения равна , т.е. равна проекции на это направление сводного вектора моментов активных сил.

                                                                            (5.5)

          Доказательство. Возьмем какую-либо траекторию движения системы , . Для нее выполнены уравнения Лагранжа первого рода (4.4).

                                                            (5.6)

Условие теоремы означает, что в каждой точке конфигурационного пространства возможным перемещением является вектор , т.е.  во все время движения и  в тех точках траектории, которые расположены на границе односторонних связей, т.е. в тех точках, в которых сосредоточена мера . Отсюда следует, что во все время движения  и . Эти равенства надо понимать как равенства мер Лебега-Стилтьеса. Домножив обе части (5.6) на , получаем

                  

или

                                                               (5.7)

Поскольку  - абсолютно непрерывна, то применима формула Лейбница для мер Лебега-Стилтьеса (см. [3]). Поэтому . Но  и  диагональная матрица масс размером , у которой на диагонали расположены тройки одинаковых чисел (массы соответствующих точек системы). Поэтому  и, значит . Тогда из (5.7)

                  

или

                  

Поскольку  есть функция ограниченной вариации, то  является абсолютно непрерывной функцией и выполнено (5.5). Доказательство закончено.

          Обычными рассуждениями из теоремы об изменении момента количества движения выводится следующее утверждение.

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

 

 

6.     Циклические интегралы и теорема Рауса.

 

В этом разделе мы покажем, что теория Рауса игнорирования циклических координат справедлива и при наличии идеальных односторонних связей.

В [3] мы вывели уравнения Лагранжа второго рода для систем с голономными удерживающими и односторонними связями. Пусть  обобщенные координаты нашей системы тогда

,                                    (6.1)

Мы считаем, что функции  имеют, по крайней мере, второй класс гладкости. Тогда траектория движения  - это абсолютно непрерывная вектор-функция, производная которой,  является функцией ограниченной вариации.

Односторонние связи в обобщенных координатах задаются системой неравенств

          ,                                                   (6.2)

Если силы потенциальны и имеют силовую функцию , то введем, как обычно, функцию Лагранжа , где  - кинетическая энергия системы. Тогда уравнения движения будут иметь следующий вид

                                     (6.3)

Здесь везде  - неотрицательные меры Лебега-Стилтьеса, каждая из которых сосредоточена на множестве моментов времени, в которые траектория движения выходит на границу соответствующего одностороннего ограничения .

                Заметим, что по форме уравнения (6.3) совпадают с уравнениями Аппеля для систем с непроинтегрированными связями.

          Далее рассуждаем по обычной схеме.

Циклические координаты. Обобщенная координата  называется циклической, если она не входит явно в функцию Лагранжа  и в неравенства односторонних связей, т.е.

                   ,   ,                                       (6.4)

Не циклические координаты называются позиционными.

Циклические интегралы. В силу уравнений (6.3) и условий (6.4) для циклической координаты  выполнено . Значит, величина   остается постоянной, т.е. является первым интегралом уравнений движения, который называется  циклическим интегралом.

Функция Рауса. Пусть координаты  являются циклическими. Им соответствуют циклические интегралы

          ,                                                  (6.5)

где  - константы интегралов. Циклические координаты в эти уравнения не входят. Будем считать, что уравнения (6.5) функционально независимы и, следовательно, разрешимы относительно циклических скоростей, и мы можем их найти как функции позиционных координат , позиционных скоростей , констант циклических интегралов  и времени:

                   ,                                 (6.6)

Составим функцию Рауса

                                                              (6.7)

Здесь всюду вместо циклических скоростей подставлены их выражения (6.6).

          Теорема Рауса. Пусть  - траектория движения нашей системы, тогда для позиционных координат выполнены уравнения Лагранжа второго рода, в которых в качестве функции Лагранжа берется функция Рауса (6.7).

                               (6.8)

          Доказательство. Для краткости записи введем вектор циклических скоростей . Тогда

                                               (6.9)

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

                                        (6.10)

Дифференцируя правую часть (6.9) получаем

                    (6.11)

Возьмем любую точку  на траектории движения, тогда выполнены (6.5) и, следовательно,

.

Поэтому (6.11) в этой точке выглядит следующим образом

                                        (6.12)

Формы (6.10) и (6.12) совпадают во всех точках траектории движения, поэтому в этих точках выполняются соотношения

                    ,   ,   ,                     (6.13)

Подставляя (6.13) в уравнения (6.3) получаем для позиционных координат уравнения (6.8). Доказательство закончено.

 

7.     Преобразование координат.

 

Известно, что уравнения Лагранжа второго рода  инвариантны относительно замены координат. В данном разделе мы покажем, что это остается справедливым и при наличии голономных односторонних связей.

Рассмотрим голономную механическую систему, которая в локальных координатах  описывается лагранжианом . Пусть на систему наложены односторонние голономные связи

                   ,                                           (7.1)

Пусть ,  - траектория движения нашей системы. Тогда для нее выполнены уравнения Лагранжа второго рода [3]. А именно, найдутся такие неотрицательные меры Лебега-Стилтьеса , , сосредоточенные в тех точках, где  (т.е. там, где траектория выходит на границу односторонних связей (7.1)), что будут выполнены уравнения

                                     (7.2)

Пусть мы переходим от старых координат  к новым , причем функции перехода  и  принадлежат второму классу гладкости и преобразование невырожденное, т.е.

                                                                                                                               (7.3)

В тех точках, где существует , также существует и , причём выполнены соотношения

                   и                                                                                          (7.4)

а также

                   и                                                 (7.5)

Рассматривая  и  как пары независимых переменных, выводим из (7.4), что

   и                                                                                          (7.6)

а также

   и                                                                          (7.7)

Для краткости записи введем векторные обозначения , . Тогда (7.2) можно переписать следующим образом

                                                                                                          (7.8)

Кроме того, обозначим

,          , 

                                                                                         (7.9)

Ясно, что

                                                                                                            (7.10)

Наша цель  - показать, что, для , будут выполнены уравнения

                                                                                                       (7.11)

Из определения (7.9) функции  имеем

                                                                                                                                                (7.12)

Пользуясь (7.4) получаем

                                                                                                        (7.13)

где знак  означает транспонирование матриц.

Поскольку на траекториях движения  - абсолютно непрерывная функция, то на траектории  применима формула Лейбница для мер Лебега-Стилтьеса [3]. Тогда из (7.13)

                                                                                  (7.14)

и, кроме того, используя (7.7) получаем

                                        (7.15)

Также на этой траектории можем написать

                                                       (7.16)

Подставив (7.12-7.16) в (7.8) получим

                               

Поскольку преобразование координат невырождено (7.3), то получаем отсюда требуемые уравнения (7.11).

 


8.     Теорема Нётер.

 

В этом разделе мы покажем, что теорема Нётер о существовании первых интегралов лагранжевых систем справедлива и при наличии односторонних связей. Как и в предыдущих разделах, наши рассуждения почти дословно повторяют традиционные (см., например, [8]). Мы будем проводить их в локальных координатах, понимая, что значение нётерова интеграла не зависит от выбора последних [9].

Рассмотрим голономную механическую систему, которая в локальных координатах  описывается лагранжианом . Пусть на систему наложены односторонние голономные связи

                   ,                                           (8.1)

          Рассмотрим однопараметрическую группу  диффеоморфизмов нашего конфигурационного пространства. Локально действие групп из этого семейства можно описать векторным полем . Пусть  - групповой параметр, а  соответствующее ему преобразование. Будем считать, что значение  соответствует тождественному преобразованию. Действие группы  можно описать функционально . При этом функции  удовлетворяют следующим уравнениям

                                                                                          (8.2)

Мы будем требовать, чтобы в рассматриваемой координатной окрестности действие группы было свободным, т.е., чтобы .

Пусть ,  - траектория движения нашей системы. Тогда для нее выполнены уравнения Лагранжа второго рода [3]. А именно, найдутся такие неотрицательные меры Лебега-Стилтьеса , , сосредоточенные в тех точках, где  (т.е. там, где траектория выходит на границу односторонних связей (8.1)), что

                                     (8.3)

Пусть  любая гладкая кривая. Применив к ней преобразования , мы получим однопараметрическое семейство кривых  . При этом . Обозначим . Тогда .

Инвариантность лагранжиана. Будем говорить, что лагранжиан  инвариантен относительно действия группы , если для любой гладкой кривой  и любых  и   выполнено

                                          (8.4)

Инвариантность односторонних связей. Будем говорить, что односторонние голономные связи (8.1) инвариантны относительно действия группы , если в любой момент времени инвариантна их граница, т.е. выполнено

, для любых  таких, что

Теорема Нётер. Пусть лагранжиан нашей системы и односторонние связи инвариантны относительно гладкого (класса два) свободного действия группы , определённой выше, тогда для любой траектории движения системы  выполнено

                                                                                                                     (8.5)

т.е. система допускает первый интеграл движения

                                                                                                                                                          (8.6)

          Доказательство. Возьмем любую точку тректории . В соответствии с теоремой о выпрямлении для обыкновенных дифференциальных уравнений [10], в некоторой окрестности этой точки систему (8.2) можно "выпрямить". Это означает, что существует такое (класса два) преобразование от старых координат  к новым . При котором уравнения действия группы (8.2) примут простейший вид

                                                                                              (8.7)

где  - постоянный единичный вектор.

          Несложно убедиться в том, что

                                                                   (8.8)

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

                                                                        (8.9)

          Инвариантность лагранжиана и односторонних связей при действии группы  сохранится и в новой системе координат. Однако в новой системе координат эти условия инвариантности становятся очень простыми. Они означают, что координата  является циклической, и выполняются условия теоремы Рауса (п. 6). Значит величина

                                                                                      (8.10)

является первым интегралом движения. Подставив сюда (8.8, 8.9) и сравнив с (8.6) убеждаемся, что . Теорема доказана.

 

9.     Несколько примеров.

 

В этом разделе, не давая подробного анализа, мы приведем несколько примеров механических систем с односторонними связями, для которых применимы результаты, изложенные в данной работе.

Пример 1.  Плоское тяжелое тело произвольной формы в вертикальной плоскости. Движение ограничено горизонтальной прямой. Это система с тремя степенями свободы. Здесь существует циклический интеграл сохранения горизонтальной составляющей количества движения системы. Исключение циклической координаты по Раусу приводит к системе с двумя степенями свободы, а именно, к плоскому биллиарду в силовом поле. Вид поля и форма области биллиарда определяется формой тела.

Пример 2. Плоское тело произвольной формы в горизонтальной плоскости. Движение ограничено круговой областью. Это система с тремя степенями свободы. Здесь существует циклический интеграл сохранения вертикальной составляющей момента количества движения системы. Исключение циклической координаты по Раусу приводит к системе с двумя степенями свободы, а именно, к плоскому биллиарду в силовом поле. Вид поля и форма области биллиарда определяется формой тела.

Пример 3. Плоское тело, свободно вращающееся вокруг вертикальной оси. В полости тела, имеющей любую форму, свободно двигается материальная точка. Здесь существует циклический интеграл сохранения вертикальной составляющей момента количества движения системы. Исключение циклической координаты по Раусу приводит к системе с двумя степенями свободы, а именно, к плоскому биллиарду в силовом поле, вид которого (как и  форма области биллиарда) определяется формой полости. Если полость заменить отрезком кривой, по которой может перемещаться точка, то задача будет сводиться уже к одномерному биллиарду в силовом поле на отрезке.

Пример 4.  Свободный волчок Лагранжа, поверхность которого является поверхностью вращения вокруг оси динамической симметрии. Движение тела ограничено горизонтальной плоскостью (). Это система с шестью степенями свободы. Здесь существует два циклических интеграла сохранения горизонтальной составляющей количества движения системы. И два циклических интеграла, имеющих место для волчка Лагранжа. Исключение циклических координат по Раусу приводит к системе с двумя степенями свободы, а именно, к плоскому биллиарду в силовом поле, вид которого (как и форма области биллиарда) определяется формой тела.

Пример 5.  Волчок Лагранжа с вертикальным люфтом точки подвеса. Отличие этой системы от обычного волчка Лагранжа состоит в том, что точка подвеса твердого тела может свободно перемещаться в пределах вертикального отрезка . Здесь существует два циклических интеграла, имеющих место для волчка Лагранжа. Исключение циклических координат по Раусу приводит к системе с двумя степенями свободы, а именно, к плоскому биллиарду в силовом поле.

Пример 6. Волчок Лагранжа с радиальным люфтом точки подвеса. Отличие этой системы от обычного волчка Лагранжа состоит в том, что точка подвеса твердого тела может свободно перемещаться в пределах некоторого отрезка расположенного на оси динамической симметрии твердого тела. Здесь существует два циклических интеграла, имеющих место для волчка Лагранжа. Исключение циклических координат по Раусу приводит к системе с двумя степенями свободы, а именно, к плоскому биллиарду в силовом поле.

         

ЛИТЕРАТУРА

 

1.        Козлов В.В., Трещев Д.В. Биллиарды. Генетическое введение в динамику систем с ударами. М., МГУ, 1991 - 168с.

2.        Иванов А.П. Динамика систем с механическими соударениями. М., "Международная программа образования", 1997 - 336с.

3.        Кугушев Е.И., Сорокина О.В. Принцип Даламбера-Лагранжа в механических системах с односторонними связями. Препринт ИПМ им. М.В. Келдыша РАН, 2002 - 32с.

4.        Жуковский Н.Е. Теоретическая механика. М., Гостехиздат, 1950 - 811с.

5.        Аппель П. Теоретическая механика. т. 2, М., Физматгиз, 1960 - 487с.

6.        Маркеев А.П. Теоретическая механика. "Регулярная и хаотическая динамика", 1999 - 569с.

7.        Березкин Е.Н. Курс теоретической механики. М., МГУ, 1974 - 646с.

8.        Голубев Ю.Ф. Основы теоретической механики. М., МГУ, 2000 - 719с.

9.           Арнольд В.И. Математические методы классической механики. М., Эдиториал УРСС, 2000 - 408с.

10.    Арнольд В.И. Обыкновенные дифференциальные уравнения. М., Наука, 1984 - 272с.