CC BY
44
Физика
УДК 531.31
Исследование неголономности некоторых гамильтоновых
полей
И. Е. Каспирович, В. А. Попова, В. И. Санюк
Кафедра теоретической физики и механики Российский университет дружбы народов ул. Миклухо-Маклая, д. 6, Москва, Россия, 117198
В классической механике понятие неголономности применяется, как правило, лишь к связям, наложенным на систему. При этом динамической системе с наложенной кинетической неголономной связью можно сопоставить векторное поле. Одной из характеристик такого поля является степень неголономности, которая определяет свойства геометрии данного поля. Однако использование этой характеристики в геометрии векторных полей ограничивалось полями в евклидовом пространстве. В данной статье предложено обобщение понятия степени неголономности на поля, определённые в неевклидовых пространствах. Для этого степень неголономности рассматривается как трёхлинейная форма. Коэффициенты этой формы, очевидно, связаны с компонентами метрического тензора пространства, в котором определено векторное поле. Соответственно, обобщение метрического тензора на случай неевклидового пространства порождает обобщения коэффициентов трёхлинейной формы, которые, в свою очередь, обобщают понятие степени неголономности. В качестве примера в данной статье проводится анализ него-лономности гамильтоновых векторных полей. Также ставится вопрос о возможности применения данного метода и о существовании механической трактовки полученных результатов.
Ключевые слова: степень неголономности, неголономные связи, гамильтоново поле, интегрируемость дифференциальных форм, форма Пфаффа, ковариант Фробениуса.
1. Степень неголономности в неевклидовом пространстве
В геометрии векторных полей поле £ называется голономным, если существует семейство поверхностей, ортогональных данному полю:
£ = , №1 + Ь&Х2 + & ёжз = 0. (1)
Иными словами, поле голономно, если соответствующее ему уравнение Пфаффа (1) интегрируемо, то есть представляет собой полный дифференциал некоторой функции. Последним определением удобнее воспользоваться для обобщения на многомерный случай. Имеет место некоторый критерий голономности: согласно теореме Якоби [1], поле голономно только тогда, когда оно ортогонально своему ротору, что позволяет ввести некоторую величину, характеризующую голо-номность поля
Р =(£ го1(0). (2)
Эту величину иногда называют степенью неголономности. Формула (2) применима для векторного поля в трёхмерном евклидовом пространстве, однако она редко находит применение в работах по неголономной механике в силу отсутствия содержательных примеров с одной кинематической связью.
Для начала необходимо обобщить определение степени неголономности для векторных полей в многомерном евклидовом пространстве. Для этого потребуется исследовать интегрируемость линейных дифференциальных форм в многомерном случае. Рассмотрим форму следующего вида, которую будем называть в
Статья поступила в редакцию 9 апреля 2015 г.
дальнейшем пфаффовой или формой Пфаффа
Ш = £хёЖ1 + &&Х2 +-----+ ^п&Хп. (3)
Уравнение ш = 0, соответственно, называется уравнением Пфаффа. Для его интегрируемости, как показал Фробениус, необходимо чтобы все коэффициенты соответствующей билинейной формы Ф были равны нулю, то есть
Ф = ёш(^) — 5ш(с1) = - Р^&х^х^ — 5х1&хз). (4)
д.6 д&
Здесь коэффициенты Р™ = —— — —— образуют 2-валентный кососимметричный
ОХ1 ох^
тензор, а форма Ф называется билинейным ковариантом Фробениуса. Таким образом, чтобы система (3) была вполне интегрируема, все Fj^ должны быть равны нулю. Это утверждение носит название теоремы Фробениуса. Для перехода от внешних форм к векторным полям необходимо рассмотреть некоторое обобщённое условие (4), которое будет являться критерием неголономности.
Действительно, рассмотрим следующее выражение, состоящее из коэффициентов коварианта Фробениуса:
= Рк] + + (5)
Как легко показать, оно преобразуется по тензорному закону, и, если в явном виде расписать все коэффициенты, то будет видно, что для трёхмерного случая это выражение будет в точности совпадать с формулой для степени неголоном-ности (2). Запишем (5) в следующем виде:
Щк = 1,3, к = 1, 2,...,п, (6)
где квадратные скобки означают альтернирование, то есть сумму по всем перестановкам индексов с учётом чётности перестановок. Этот тензор называется объектом неголономности, а равенство нулю всех его компонент и есть необходимое и достаточное условие неголономности векторного поля в многомерном евклидовом пространстве.
Далее рассмотрим обобщение формул (1), (2), (5) на случай полей в неевклидовых пространствах. Начнём с трёхмерного случая. Формально, в формуле (2) для степени голономности стоит скалярное произведение поля на свой ротор, поэтому, на первый взгляд, для обобщения достаточно вместо евклидовой метрики использовать, например, риманову. Но при замене координат получившаяся формула будет вести себя странно, так как обычный, полярный вектор п умножается на псевдовектор. Обойти это можно следующим образом: для удобства запишем формулу (2) в виде трёхлинейной формы, получим:
Р = Ш^к £гдэ ^к. (7)
Задача сводится к выявлению связи между коэффициентами формы (7) и компонентами метрического тензора. К слову, эти коэффициенты изменяются по трёхковариантному закону, то есть являются тензором. Также видно, что при = наша формула для степени неголономности переходит в обычный вид (2) для евклидова пространства. Стоит отметить, что предполагаемое обобщение не потребует замены частных производных на ковариантные [2], так как кососимметричные коэффициенты F^j взаимно аннулируют все неевклидовы дополнения, оставляя лишь евклидовы частные производные. Теперь запишем формулу (2) в виде скалярного произведения р = д^^г^, где г\ = £íjkдj£к.
Сравнивая коэффициенты при соответствующих членах, получим следующее выражение для трёхиндексных коэффициентов в (7):
= 9И . (8)
Теперь распишем в виде трёхлинейной формы формулу для неголономности полей в многомерном пространстве:
Пг^к = Шцг^дг^г + Ши + Шцк£г 9г£к +----+
+----+ ^г] г£гд] & + ^гЦ^гд] ^ + £гд] Ск +----+
+----+ ШкЫ^к дк + Шкк]^к дк ^ + Шккк £к дк £к. (9)
Заметим, что объект неголономности, определённый таким образом, не всегда преобразуется по тензорному закону, он является тензором, только если векторное поле определено в евклидовом пространстве с декартовыми координатами. Таким образом, если поле является голономным в одной системе координат, то оно голономно и в другой системе координат. Элементы с повторяющимися индексами, которые были тождественно равны нулю в евклидовом случае, могут стать ненулевыми. Таким образом, для анализа неголономности полей, заданных в неевклидовом пространстве, необходимо учитывать все компоненты объекта неголономности.
2. Применение в механике и некоторые свойства
Для начала укажем некоторые свойства объекта неголономности.
Теорема 1. Пусть а(ж) и Ъ(ж) векторные поля в и-,мерном евклидовом пространстве, причём аг = /(х)Ьг. Здесь х = (х\,... ,хп), а / — произвольная гладкая функция. Тогда объекты неголономности этих векторных полей связаны нелинейным соотношением:
^'к (а) = / 2(ж)Пг^к (Ъ).
Доказательство. Для доказательства распишем объект неголономности в явном виде:
и- к(а) = а-(— - —+ а-( — - —^ + ак( — - —^ гз г \дxj дхк) ° \дхк дхг) \дхг дх^)
= /(*)М/(,)£ -/СО£ + *^ -ь,+ +
+, , (1) ^ -, ^+Ь1?Ж - ^
(х)
Все слагаемые, содержащие производные типа —-, взаимно уничтожаются,
охг
поэтому в формуле остаётся лишь объект неголономности поля Ь, умноженный на квадрат функции /(ж). □
В частном случае, когда /(ж) = 1/д/(Ъ, Ъ), а является просто нормированным полем. Тогда, применяя полученное выше утверждение, получим:
Пг ок(а) = (ъ Ъ) Пок(Ъ).
Таким образом, теорему Якоби можно обобщить на любое ненормированное векторное поле, поделив степень неголономности на квадрат модуля векторного поля.
Теперь рассмотрим объект неголономности от суммы двух полей.
Теорема 2. Пусть a(x) и Ъ(ж) — векторные поля в n-мерном евклидовом пространстве, .здесь х = (х\,... ,хп), тогда объект неголономности суммы векторных полей определяется по формуле:
üijk(a + Ъ) = üijk(a) + üijk(b) + a\idjbk¡ + b^djak¡.
В трёхмерном случае формула приобретает вид:
^ijk(а + 6) = ^íjk(а) + ^íjk(b) + (a, rotó) + (b, rota).
Доказательство. Для доказательства распишем объект неголономности в явном виде:
Г. I , u\ í , í. \(dak daj dbk dbj \
üijk (a + b) = (üí + bí)[---+ ---+
\OXj OXk OXj OXk J
+ (а - + b ) (- + ^L _ +ía + b )(_ + _ ^Л
3 3 \dxk dxi dxk dxi) \dxi dxj dxi dxj)
Раскрывая скобки и группируя слагаемые, получим в точности формулу для неголономности суммы полей. Таким образом, утверждение доказано. □
Понятие неголономности встречается в аналитической механике при рассмотрении свойств связей, наложенных на динамическую систему. Если связь выражена неинтегрируемым дифференциальным выражением, то она называется неголономной. В случае, если связь является линейной формой по скоростям, то критерий неголономности такой связи совпадает с критерием неголономности векторного поля, соответствующего связи. Одним из методов решения задач неголономной механики является переход к квазикоординатам. Уравнение Лагранжа в квазиикоординатах ^k, как известно [3], имеет вид:
d dL* dL* dL*
+ ^kij -^-Kj = 0,
di дтгk d^k dñi
где L* — лагранжиан, зависящий уже от квазикоординат, а ^kíj — трёхиндексные символы Больцмана-Гамеля, определяемые следующим образом:
/ dais daii \
^ = bsk Ц — - —j .
Здесь bsk компоненты обратной матрицы к а^k. Видно, что символы Больцмана-Гамеля можно переписать через введённые ранее коэффициенты билинейной формы Фробениуса Fij, тем самым можно проследить некоторую связь с объектом неголономности.
3. Вычисление объектов неголономности для гамильтоновых полей
Как известно [4], гамильтоново поле определяется как поле изменения координат фазового пространства во времени. Оно имеют следующий вид: £ = {¿1, ¿2,..., ¿2 п}, где Хг — координаты фазового пространства, а производные определяются симплектической структурой следующим образом:
^г — ¿г — Н} — Qz ' ,
где Н — гамильтониан системы.
Известно [5], что гамильтоново векторное поле является объектом симплекти-ческой геометрии. Фазовое пространство четырёхмерно, поэтому воспользуемся обобщённой формулой (9). Сначала вычислим компоненты Шг^к, воспользовавшись формулой (8), в которой gгj определим как метрические коэффициенты симплектической формы в четырёхмерном случае. Следовательно, ненулевыми коэффициентами будут, например, 12, 21, ш323, ш332,.... Заметим, что в этом случае объект неголономности не всегда кососимметричен, однако ввиду того, что слагаемые в каждом элементе с перестановкой индекса меняют знак, нам достаточно выписать лишь по одному элементу из всевозможных перестановок.
о , 2 9 ¿1 А , (9 6 да 2
0123 = ч-Ъь) -—
\дС2 дСз)1
0124 = -^ аа - Кг) + ЧК - д<2)
0 = да 3
0134 = - Ч асз- ^3
0234=- Из- 5с4}+61 ^
^ а 9 (2.
( д^Л
^\ди дС1 у
*С2 ^С3) .
(10)
Подставляя в формулу определение гамильтонова поля через канонические уравнения, получим следующие формулы:
0
123
_дн ( д2н д2н \ дн / д2н д2н \
0
124
др 1\др2дд 1 др 18дд 1\др 1др2 дд 18,
дн ( д2н д2н \ дн / д2н д2н \
др 2\др 2дд 1 др 18 д2 ) д д2\др 1др 2 дд 18 ® /'
дн ( д2н д2н \ дн / д2н д2н
др 2\др1др 2 дд 18 д2) 8 д2\др2д д1 др 18 д2
_ дн / д2н д2н \ дн / д2н д2н
134 дд 1\8р2дд 1 др 1дд2) др 1\др1др2 дд 1дд2
234
(11)
Выше были рассмотрены только компоненты с неповторяющимися индексами. Однако, как отмечалось ранее, для полноценного анализа необходимо учесть все ненулевые элементы объекта неголономности, так как в них будут встречаться все коэффициенты Рг^, следовательно, необходимым и достаточным условием для голономности гамильтонова поля будет равенство нулю всех этих коэффициентов,
т.е. ^ = 0, VI ,3 = 1,..., п. (12)
Заметим, что для произвольного негамильтонова векторного поля в симплек-тическом пространстве это не так. Гамильтонианы, как правило, квадратичны по обобщённым импульсам, следовательно, в формулах для объекта неголономно-сти будут появляться линейные, квадратичные и формы третьей степени по импульсам. Тогда равенство нулю всех коэффициентов Fгj автоматически означает равенство нулю всех компонент объекта неголономности. Выпишем их в явном виде:
Р12 = Р34 =
8 2н
д 2н
34
др 18 д2 др 2дд 1
Р24 =
д2н д2н
+
др 2
д д1
^14 = Е23 =
23
Р13 =
д2н др2
д2н
др1др2
д 2н
+
д2н
дд 18д2,
+
дд2
4. Заключение
Анализируя получившийся результат, приходим к выводу о том, что даже на первый взгляд, поля простых динамических систем, таких как свободное движение на плоскости, будут неголономны. Одно из немногих голономных гамильто-новых полей будет поле порождаемое гамильтонианом
Приведённые выше соображения могут рассматриваться как некоторое подтверждение высказывания одного из пионеров неголономной динамики Воронца о неголономности гамильтоновых динамических систем [6].
1. Аминов Ю. А. Геометрия векторного поля. — М.: Наука, 1990. — 215 с.
2. Дубровин Б. А, Новиков С. П., Фоменко А. Т. Современная геометрия: Методы и приложения. — М.: Наука, 1986. — 760 с.
3. Неймарк Ю. И., Фуфаев Н. А. Динамика неголономных систем. — М.: Наука, 1967. — 521 с.
4. Трещев Д. В. Гамильтонова механика // Лекционные курсы НОЦ. — 2006. — № 4. — С. 64.
5. Татаринов Я. В. Лекции по классической динамике. — М.: Изд-во Моск. ун-та, 1984. — 296 с.
6. Сумбатов А. С. Интегралы, линейные относительно скоростей. Обобщения теоремы Якоби // Итоги науки и техники. Общая механика. — 1979. — Т. 5. —
UDC 531.31
Analysis of Nonholonomicity Value of Some Hamiltonian
Fields
I. E. Kaspirovich, V. A. Popova, V. I. Sanyuk
In classical mechanics such notion as nonholonomicity is applied only to constraints put on a dynamical system. Besides, Pfaffian nonholonomic constraints might be associated with vector fields. The Nonholonomicity value is one of the principal characteristics of such fields, which determines properties of geometry of these vector fields. However, the application of this characteristic to the geometry of vector fields was restricted only to fields in Euclidean spaces. Some generalization of nonholonomicity value of vector fields in non-Euclidean spaces is proposed in this paper. For this purpose the nonholonomicity value is considered as a trilinear form. It is obvious that the coefficients of this form are connected with the components of the metric tensor of the space, where a vector field is defined. So generalization of metric tensor to non-Euclidean spaces generates the generalization of the coefficients of trilinear form, which in its turn generates the generalization of nonholonomicity value. As an example, the nonholonomicity values of Hamiltonian vector fields in sympletic spaces are analyzed in this article. Also it is important to find out whether a mechanical interpretation of the obtained results exists and can we actually apply this method to Hamiltonian fields.
Key words and phrases: nonholonomicity value, nonholonomic constraints, Hamiltonian vector field, integrability of differential forms, Pfaffian form, Frobenius theorem.
Литература
С. 3-57.
Department of Theoretical Physics and Mechanics
Peoples' Friendship University of Russia 6, Miklukho-Maklaya str., Moscow, Russia, 117198
References
1. Y. A. Aminov, Geomerty of Vector Fields, Nauka, Moscow, 1990, in Russian.
2. B. A. Dubrovin, S. P. Novikov, A. T. Fomenko, Modern Geometry: Methods and Applications, Nauka, Moscow, 1986, in Russian.
3. Y. I. Neymark, N. A. Fufaev, Dynamics of Nonholonomic Systems, Nauka, Moscow, 1967, in Russian.
4. D. V. Tretschev, Hamiltonian mechanics, Lectures of the Educational Center of Steklov Mathematical Institute (4), in Russian.
5. Y. V. Tatarinov, Lectures on Classical Dynamics, Moscow University Publishing, Moscow, 1984, in Russian.
6. A. S. Sumbatov, Linear by Velocities Integrals. The Generalizations of Jacobi's Theorem, The Results of Science and Technology. General Mechanics 5 (1979) 357, in Russian.
