УДК 514.85 Вестник СПбГУ. Математика. Механика. Астрономия. 2017. Т. 4 (62). Вып. 4
МЯО 53205
ОСОБЕННОСТИ ДВИЖЕНИЯ МАЯТНИКА
С СИНГУЛЯРНЫМ КОНФИГУРАЦИОННЫМ ПРОСТРАНСТВОМ
С. Н. Бурьян
Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
Рассматривается плоский двойной математический маятник, незакрепленный конец которого движется по эллипсу. В общем случае конфигурационное пространство маятника представляет собой две непересекающиеся кривые. Можно подобрать параметры так, что эти кривые будут пересекаться под ненулевым углом. При определенных условиях на кривизну эллипса кривые не только пересекаются, но и имеют касание первого порядка. В этом случае возникает геометрическая неопределенность — по какой ветви двигаться маятнику после прохождения особой точки.
Показывается, что для случая, когда кривые конфигурационного пространства пересекаются под ненулевым углом, множители Лагранжа по модулю стремятся к бесконечности при движении к особой точке. Наблюдаемое же движение оказывается динамически определенным. Маятник всегда переходит с одной ветви на другую при прохождении особой точки. Предложено качественное объяснение этого эффекта. Библиогр. 11 назв. Ил. 3.
Ключевые слова: конфигурационные пространства с особенностями, множители Лагран-жа, двойной маятник.
1. Введение. Механические системы с особенностями конфигурационного пространства естественно возникают в классической механике. Например, если в шарнирном механизме часть звеньев может складываться в одно звено, то его конфигурационное пространство уже не будет многообразием. Так, шарнирная рама со звеньями (1; 1; 1; 1) и одним закрепленным звеном допускает три типа движения. Его конфигурационное пространство представляет собой объединение трех попарно пересекающихся окружностей. Другим примером служит двойной плоский маятник со специальным ограничением, рассматриваемым далее в статье. На этом маятнике физически реализуется особенность касания первого порядка.
На данный момент нет общей теории для описания поведения систем с особенностями конфигурационного пространства. Постулируется, что вне особенности движение механической системы описывается классическими уравнениями Лагранжа или Гамильтона.
С точки зрения лагранжевой механики двойной маятник со связью — это механическая система с тремя голономными связями. Уравнения связей независимы вне особой точки. В особой точке, где геометрически возможен переход с одной кривой движения на другую, ранг матрицы градиентов уравнений связей падает на единицу. Это приводит к трудностям, связанными с уравнениями Лагранжа первого рода. Похожие результаты, но в неголономной механике, получены в [1] для движения балки с двумя коньками. При анализе и составлении уравнений будем следовать [2, 3]. Вопросы о реализации связей рассматриваются, например, в статье [4].
Существует несколько подходов к решению проблемы особенностей, основанных на обобщении гамильтонова формализма, который связан с различными обобщениями понятия гладкого многообразия: пространства Сикорского [5], пространства
(¡5 Санкт-Петербургский государственный университет, 2017 https://doi.org/10.21638/11701/spbu01.2017.402
Фрелихера [6], диффеология [7], диффеотопия [8] и др. Каждый из указанных методов позволяет решить отдельные вопросы механики и дифференциальной геометрии (см., например, [9-11]) и претендует на универсальность. Однако есть основания полагать, что применение перечисленных методов к задачам теоретической механики имеет существенные ограничения.
В этой работе для исследования маятника рассмотрен «предельный» лагранже-вый подход.
2. Маятник. Рассмотрим плоскую механическую систему (рис. 1): к точке А(0,0) шарнирно прикреплен стержень ¿1 = АВ длины 2, а к точке В(х\,у\) шар-нирно прикреплен стержень ¿2 = ВС длины 1. В дальнейшем будем называть ее свободным маятником. Обычные обобщенные координаты этой системы —это пара углов (1р,ф), которые отсчитывают отклонение стержней ¿1 и ¿2 от вертикали соответственно. Конфигурационное пространство — двумерный тор Т2 в четырехмерном пространстве.
Рис. 1. Свободный двойной маятник.
Пусть точка C(x2, У2) движется по линии, заданной уравнением F(x,y) = 0 (рис. 2). Будем называть эту систему маятником со связью. У такой системы остается одна степень свободы.
Зафиксируем значение |AC| = d (1 < d < 3) и направление вектора AB. Рассмотрим ДABC, длины всех сторон которого известны, следовательно он определяется с точностью до отражения относительно оси AC. Значит, есть два положения точки B, отвечающие ему: когда стержень I2 находится правее стержня li, и когда стержень I2 находится левее li.
Тогда и для каждого d (1 < d < 3) есть два различных описанных положения системы. Множеству положений маятника, при которых короткий стержень опережает длинный (или запаздывает), соответствует гладкая кривая в конфигурационном пространстве маятника. Эти две кривые будем также называть ветвями. Чтобы перейти
с одной ветви на другую, маятник должен пройти через положение, при котором его стержни лежат на одной прямой. Это возможно лишь при |АС| = 1 (сложенный маятник) или |АС| = 3 (раскрытый маятник). Далее мы будем рассматривать первый случай. Точку конфигурационного пространства, соответствующую положению маятника с |АС | = 1 естественно называть точкой ветвления — геометрически движение из нее возможно по двум различным кривым.
В общем случае ветви пересекаются под ненулевым углом. Рассмотрим окрестность точки ветвления.
Параметризация исходной системы через координаты двумерного тора не является удобной для вывода уравнений движения. Рассмотрим вместо нее полярную систему координат. Ее полюсом будет точка А, а полярной осью (АХ) —вертикальная ось с положительным направлением по направлению действия силы тяжести. Ось АУ направим ортогонально АХ (рис.2).
Рис. 2. Двойной маятник со связью.
Обозначим угол между осью АХ и прямой АС через ф. Введение этого угла, отличного от угла у, позволит выписать явную зависимость координат и скоростей точек В и С при заданной линии связи. Пусть связь ^(х, у) = 0 такова, что прямая у = tg ф х пересекает ее в одной точке (точке С). Это позволит также параметризовать кривую углом ф: уравнением связи является ! = !(ф).
Имеет смысл рассматривать !(ф) как некоторую гладкую функцию на интервале [ф1, ф2]. Например, пусть связь реализована в виде эллипса с фокусами и ^ на оси АХ (выше точки А), точка (а, 0) —центр эллипса, Г1 и Г2 —большие полуоси. Уравнение описанного эллипса примет вид
(ж2 - а)'
.+ у1= 1
^ „2
Г.
(1)
2
Тогда при ф € [— \, \] прямая у = tg фх пересекает эллипс в точке С = (х, xtg ф), где
= х(ф) =
■ +
tg2 ф
1
tg2 ф
(2)
1
1
а
2
а
2
2
2
Г
Г
Г
1
1
2
X
1
2
Г
Г
1
2
Выражение под корнем в (2) всегда положительно, т.к. |a| < ri. Остальные функции в (2) гладкие, поэтому и с1(ф) = \Jх(ф)2 + у(ф)2 = х(ф) \/l + tg2 ф — гладкая функция. Явное выражение для функции ¿(ф) не столь короткое.
Рассмотрим зависимость координат точек B и C от пары (d, ф) — полярных координат. Координаты точки C записываются просто:
(x2,y2) = (dcos ф, dsin ф).
Обозначим теперь в = ZCAB. Тогда можем записать Z(XAB) = ф + в и, значит, В = (2 cos(</>+0), 2 sin(</>+0)). Из A ABC по теореме косинусов имеем cos в = (d2+3)/4d. Тогда sin в = ±Vl0d2 - d4 - 9/(4d), где знак «±» отвечает за выбор ветви маятника. По формулам суммы косинусов и синусов получаем
, , /З + d2 л/Юd? - d4 - 9 3 + d2 л/lOd2 - d4 - 9 \ (X, , У1 ) = cos ф T -^-Sin Ф, — sin ф ± --- cos ¿ j .
(3)
Рассмотрим движение при заданной параметризации d = ¿(ф). Каждая ветвь движения параметризуется углом ф. Скорости точек в этой параметризации имеют вид
(VCx, VCy) = (x2, y2) = (d' cos ф — dsinф, d' sin ф + dcos ф)ф, (4)
• ,^2-3 %/lOd2—d4 —9 . , 9 — d4 ..
— sin©-;—heos©—-^-d=Fcos©---=Fsm©---d 5
^ 2d ^ ^ 2d? ^ ^ 2d ^ 2d2A/10d2 —d4 —9 / W
^ = =
/ ,3+d2 d2 — 3 . , Vl0d2-d4-9 , 9 — d4 „V
= cos©-;—hsm©—-^-d=Fsm©---±cos©---d (6)
^ ^ 2d ^ ^ 2d2 ^ 2d 2d2A/10d2 —d4 —9 уГ W
3. Связь скоростей и кривизны линии связи. Из формулы (3) следует, что при d = l в окрестности положения равновесия координаты точек B и C являются гладкими функциями угла ф. Рассмотрим положение равновесия. Пусть в окрестности Uo точки ф = 0 функция d^) является гладкой, d(0) = 1 и d^) = d(-ф). Тогда d'(0) = 0. Предположим, что d"(0) > 0. По лемме Адамара получаем
<ф) = 1 + ф2д(ф), д(ф) G C~(Uo).
Тогда справедливо неравенство 2д(0) = d''(0) > 0. По непрерывности, при малых ф > 0 определена функция и := \/d — 1 = ф^/д, и она является гладкой. Следовательно, можем записать VlOd2 - d4 — 9 = \Jd — 1 \/[d + 1)(9 — d2) = и{ф)А{ф), где u,A G C~(U0).
Вычислим формулу кривизны со знаком kor кривой d = d^), заданной в полярных координатах (d, ф) :
х'у" - у'х" _ d2 + 2d'2 - dd"2 = (х'2 +у'2)3/2 = (¿2 + ¿/2)3/2 '
В нашем случае d(0) = 1, d'(0) = 0, поэтому kor (0) = 1 — d''(0). Если кривизна линии d = d^) в нуле отлична от 1, то d''(0) = 0, и ветви различаются направлением скоростей (5), (6).
Рассмотрим случай kor (0) = 1. Положим d = 1 + ф"д„(ф), gn(0) = 0, n > 2 (где n — первая такая степень в разложении d в ряде Адамара). В этом случае имеем
d' фп-1С (ф)
К ' -+0, ф —^ о,
где C(ф) —некоторая гладкая функция. Следовательно, скорости по двум ветвям (5), (6) в положении равновесия совпадают. Таким образом, условие равенства скоростей описаны в геометрических терминах: достаточно выполнения d'(0) = 0,
kor (0) = 1.
Для эллипса (1) условие на единичную кривизну можно получить из стандартной параметризации x = a + ri cos в, y2 = Г2 sin в:
М©) = , 2 ■ 2п ?Г22-J7™ , МО) = Г-\= 1 ^ Г! = г22. (7)
(r2 Sin2 в + r2 cos2 в)3/2 r2
Значит, эллипс должен быть «вытянут» по оси AX.
4. Движение свободного маятника вблизи точки ветвления. Рассмотрим свободный двойной маятник со стандартными обобщенными координатами (у, ф) и изучим его поведение в окрестности точки (0, п). Массу в точке B обозначим через mi, в точке C — через m2.
Координаты точек B(xi, yi) и C(x, У2) задаются равенствами
x1 = 2 cos у, yi = 2 sin у, x2 =2 cos у + cos ф, y2 = 2 sin у + sin ф.
В положении равновесия у = 0, ф = п имеем
Ф=^У1, Ф = У1~У2- (8)
Уравнения Лагранжа второго рода в криволинейных координатах (у, ф) (величина n = mi /(mi + m.2)) запишем в виде
¿1<у + n/^cos^ — ф)ф + sin^ — ф)ф2] + g sin у = 0, (9)
/2ф + — ф)у + sin(у — ф)у2 ] + g sin ф = 0. (10)
Это система линейных уравнений относительно обобщенных ускорений ф и у, и ее определитель Д всегда положителен: Д = li/2 — n/i/2 cos(у — ф) > 0. Уравнения движения разрешимы относительно старшей производной.
Итак, зная начальные данные у* = 0, ф* = п, у*, ф*, решим задачу Коши. Пусть у(£), ф(£) —ее решение. Определим, будет ли короткий стержень опережать длинный или наоборот. Для этого перейдем в более удобную, подвижную систему отсчета, связанную со звеном AB: одну ось направим по этому звену, а другую, AW, ортогонально ей. В этой системе координат точка B неподвижна, а точка C движется по окружности. Ось AB разбивает систему координат на две части — правую и левую, которые соответствуют случаям, когда короткий стержень опережает длинный и наоборот. Переход происходит, когда маятник проходит через AB.
Новые координаты выражаются через старые как поворот на угол у:
жЛ _ Í cos у sin у\ /2 cos у\ _ Í2 yi} у—sin у cos уу у 2 sin у у ^0
ж2\ / cos у sin у\ /2 cos у + cos ф\ /2\ /cos^ — у)
у2 у у — sin у cos уу у 2 sin у + sin ф у \^т(ф — у)
Значит, за переход отвечает знак sin^ — у): если он меняет знак, то происходит переход с одной ветви движения на другую.
Рассмотрим случай, когда в положении равновесия свободный маятник имеет те же скорости, что и маятник со связью, который пришел в особую точку по эллипсу с kor (0) = 1. Слагаемые с особенностями в (5), (6) не равны нулю, и поэтому из (8) имеем ф — у = 0. Более того, знак последнего выражения зависит от выбора ветвей. Разложим синус разности углов ф и у в ряд:
sin^(t) — у(*)) = 0 + (ф — у)* + o(t).
При достаточно малых t > 0 знак последнего выражения определяется знаком ф — у, поэтому при t > 0 и t < 0 синус разности углов ф и у меняет знак — происходит смена ветвей.
5. Обобщенные реакции. Изучим вопросы применимости уравнений Лагран-жа первого рода к движению двойного маятника по эллипсу (когда конфигурационное пространство имеет особенности). Будем рассматривать случай, когда mi = m-2 = m, т. е. массы грузов одинаковы.
Выпишем уравнения голономных связей в пространстве свободной системы R4 = (z) = (xi, yi, Ж2, У2):
f1 = \{x\+y\-A)=Q-f = \({xí-x2)2+ (У1-У2)2-l) = Q-
Первое уравнение соответствует фиксации длины стержня AB, второе — фиксации длины стержня BC. Третье уравнение — условие, что точка C движется по эллипсу (1). Градиенты связей определяются формулой
'V/^ / xi yi 0 0
V/2 I = I xi — Х2 yi — y2 X2 — ж i y2 — yi vV/V V 0 0 X2 — a/r2 У2/Г2
В положении равновесия (2, 0,1, 0) ранг матрицы связей равен 2. В окрестности этой точки он равен 3, так как есть невырожденная подматрица 3 х 3.
Уравнения Лагранжа первого рода для голономных связей [2, с. 26-27] (Z — вектор действующих сил) записываются в виде
mZM = + ^ Лк V/к, м = 1..4, к =1, 2, 3; (11)
К* ^
=°- (12)
Вектор сил . имеет координаты (шц, 0, шц, 0), т. к. сила тяжести действует вдоль оси АХ. Из последних уравнений получаем линейную систему относительно множителей Лагранжа с матрицей = V/к • V/— матрицей скалярных произведений градиентов. Вектор правых частей и при этом составлен из выражений
Явное вычисление показывает, что справедливо выражение
( (X)2 + (у1)2 +
дг . . , -iJ2df>
(x1 - ¿2)2 (yl - У2)2 I . (13)
fe)2 1 (У2)2 , x2-a
r 2 ' r — ^^
Вне особой точки матрица невырожденная, так как связи независимы. В особой точке имеем V/1 = (2, 0,0, 0), V/2 = (l, 0, -1, 0) , V/3 = (0, 0, c, 0), где c = (1 -a)/(r2). Заметим также, что 1 — а = ri, если эллипс (1) при ф = 0 проходит через точку (1, 0). В данном пункте рассматривается только этот случай.
Матрица скалярных произведений градиентов в особой точке имеет вид
/4 2 0 \
Gkk*U=(2.0.i.0) = (2 2 -cl . (14)
0 -c c2
Ранг этой матрицы равен двум, поэтому не для всех векторов правых частей есть решение. Если правая часть есть U = (ui,u2,u3), то для существования решения необходимо и достаточно, чтобы выполнялось равенство
Мз , Ul п
--\-U2- — =0. (15)
c2
Заметим, что при kor (0) = 1 вектор скорости изображающей точки z = (¿1, yi, ¿2, У2) в положении равновесия w = (2, 0,1, 0) будет распределен так: Z = (0, 2v, 0, v) (см. (5), (6)), где v = ф. Подставив эти значения в формулу (13), получим вектор U0 = (4v2 + 2g, v2,v2/r2 + cg). Проверим условие (15): u3/c + u2 - u1/2 = v2(1/cr2 - 1). Получаем значение выражения cr2 = r2(1 - a)/r2 = r2/ri = 1 для эллипса с kor(0) = 1 (см. (7)). Значит, система |г=(2.о.1.о)(Л) = U0 имеет решение. Обобщенные реакции связей в этом случае определены неоднозначно.
Рассмотрим теперь движение по эллипсу с kor (0) = 1. В этом случае скорости точки B при движении (в R4) по разным ветвям при подходе к положению равновесия (2,0,1,0) (или, равносильно, ф = 0) различны. Согласно (6) получаем y2 = 2ф + C, где С = 2и'(0),и = y/d — 1. Скорости остальных координат в положении равновесия не изменятся (d(0) = 0). Скорость изображающей точки в положении равновесия теперь есть Z = (0, 2v + C, 0, v). Мы уже знаем, что при C = 0 матрица (14) имеет решение для правой части (13). Матрица (14) зависит только от c = 1/r1. Если поменять Г2, то матрица не изменится, но мы получим уже трансверсальное пересечение ветвей движения маятника. Добавка C в скорость изображающей точки добавляет к правой части (13) вектор U1 = (2vC + C2, 2vC + C2 - 2vC, 0). Для вектора U1 условие (15) выполняется только при C = 0, т. е. когда r1 = r22 , т. е. в случае касания ветвей движения, который сейчас не выполнен. Тогда, если бы система Сккф|2=(2.0.10)(Л) = U0 + U1 имела решение, то имела бы решение и Сккф|2=(2.0.10)(Л) = U1, что не выполнено.
В трансверсальном случае множители Лагранжа как решение системы (12) при ф ^ 0 по абсолютной величине стремятся к бесконечности. Действительно, решая систему (12) по правилу Крамера, замечаем, что в знаменателе будет стоять определитель системы матрицы (14), которая не вырождена в окрестности особой точки и ее ранг падает на 1 в точке ветвления, а числитель при ф ^ 0 имеет конечный ненулевой предел, поскольку каждые два столбца матрицы линейно независимы, а вектор
Uo + Ui не лежит в линейной оболочке столбцов Поскольку числитель имеет
ненулевой предел, то множители Лагранжа |ЛК| ^ œ, ф ^ 0.
Поэтому и с точки зрения предельного перехода при движении системы в транс-версальном случае к точке ветвления возникает вопрос, связанный с реализацией связи как силы реакции.
Для особого эллипса решение системы существует. В точке (2, 0,1, 0) у множителей Лагранжа есть неопределенность вида 0/0 — числитель тоже обращается в 0.
6. Уравнение движения на линии связи. Рассмотрим уравнения Лагранжа второго рода и уравнения Гамильтона, которые можно написать только вне особой точки. Изучим, что происходит при смене ветвей в точке ветвления. Угол в := BAC отсчитывается против часовой стрелки, как и угол ф. По теореме косинусов получаем
i 3 + d2
<9± = e+íd) = ± arceos-—.
4d
Выбор знака означает выбор ветви, когда длинный стержень опережает короткий или наоборот. Имеем
(x±,y±) = (li cos^ + e±(d)),1i sin^ + в±(d))). Выпишем скорости точек, взяв за параметр движения угол ф:
(x±, y±) = /1ф(1 ± в^d/)(- sin^ ± в), cos^ ± в));
V2 = (x±)2 + (y ±)2 = l 2ф2(1 ± в^d¿)2;
V2 = (x2)2 + (y2)2 = (d2 + (d; )2)ф2.
Кинетическая и потенциальная энергии запишутся в виде
ТНФ, ф) = (I(1 ± W2 + у(d2 + (dtfj) ф2 := МФ)(Ф)2; (16)
V±(ф) = —mgl 1 cos^ ± в(й)) — mgdcosф. (17)
Далее будем писать просто F(ф) и V(ф), подразумевая движение по одной из ветвей. Запишем оператор Лагранжа:
d dT dT
7гТф-ЭФ = 2тФ + Р*тФ)-
Обобщенная сила Q^ равна — V. Уравнение движения примет вид
ф 2Р{ф)
Вторую производную всегда можно выделить, поскольку Г > 0 за счет слагаемого + (с1'(р)2)ф2 в (16). Все участвующие в уравнении функции зависят лишь от координаты на линии связи.
Уравнения движения можно записать в форме Гамильтона: обобщенный импульс
2
р = дЬ/дф, тогда Н = рф — Ь = Т + У. Энергия системы имеет вид Н = ^рщ + У(Ф)-Уравнения движения запишутся следующим образом:
p =
Ф--
Р'(Ф)
Р
2 FW
Р2 - Уф(Ф);
(18)
Вне окрестности сингулярной точки на каждой взятой ветке правые части уравнений есть гладкие функции, поэтому при движении до положения ф = 0 координаты изменяются гладко.
В случае ког(0) = 1 имеем значение i1±(0) = Щf- + 2m в- Можно также проверить, что F'(0)± = ±c, где c — ненулевая константа. Так как V'(0)± = 0, то в случае перехода с одной ветви движения на другую у системы происходил бы разрыв вдоль оси p: при этом вектор поля сохранял длину, но менял бы направление на симметричное. На оси импульсов при ф = 0 вектор скорости (cip2,c2p) переходил бы в вектор (—cip2,C2p) для постоянных ci,2, т.е. вектор в фазовом пространстве должен был «вращаться».
7. Эксперимент. Было сделано несколько экспериментальных моделей двойного маятника, свободный конец которого движется по эллипсу. Связь реализовывалась в виде нити, которая крепилась в точках Fi и F2 — фокусах эллипса (они выбирались на оси AX выше точки A); нить проходила через кольцо, закрепленное в точке C (рис.3). Связь была неудерживающая. Зафиксировав угол у, точку C можно перевести с одной ветви движения на другую, меняя только угол ф. Данная система зависела еще от одного параметра — длины нити l. Если l = lo = 1 — a (см. (1)), то маятник может (геометрически) пройти положение, соответствующее ф = 0 только в «сложенном» виде. Если же длина нити l > lo, то маятник не может, оставаясь на связи, перейти с одной ветви движения на другую, так как на связи нет «сложенного» положения. Но в этом случае вблизи ф = 0 маятник может перейти с одной ветви (по связи) на другую, покинув связь.
Рис. 3. Модель двойного маятника со связью. Вестник СПбГУ. Математика. Механика. Астрономия. 2017. Т. 4(62). Вып. 4
В ходе экспериментов фокусы эллипса не менялись, но менялась длина нити. При / = /о + е (е достаточно мало) наблюдалось следующее движение: при начальных отклонениях системы |ф| < п/2 (точка C была на эллипсе) маятник вначале двигался по одной линии связи, вблизи точки (у, ф) = (0, п) он сходил со связи; как свободный маятник он производил удар о связь (при этом меняя ветвь движения, т. е. точка C проходила через ось стержня AB), а затем двигался по другой линии связи. При е ^ 0 точки схода со связи и удара стремились к точке (0, п). При е = 0 (т.е. когда конфигурационное пространство маятника стало пространством с особенностью) движение маятника выглядело симметрично: точка C всегда оставалась на эллипсе, и маятник при проходе особой точки всегда переходил с одной ветви на другую.
Литература
1. Закалюкин И. В. Особенности вырождения неголономных связей и управляемость // Труды МАИ (электр. журнал). 2010. Выпуск №39.
2. Зегжда С. А., Солтаханов Ш. Х., Юшков М. П. Уравнения движения неголономных систем и вариационные принципы механики. Новый класс задач управления. М.: ФИЗМАТЛИТ, 2005. 272 с.
3. Поляхов Н.Н., Зегжда С. А., Юшков М.П. Теоретическая механика. Л.: Изд-во Ленингр. ун-та, 1985. 536 с.; М.: Высшая школа, 2000. 592 с.
4. Журавлев В. Ф. Понятие связи в аналитической механике // Нелинейная динам. 2012. Т. 8, №4. С. 853-860.
5. Sniatycki J. Orbit of families of vector field on subcartensian space. arXiv:math/0211212 [math.DG].
6. Batubenge T.A., Sasin W. An approach to hamiltonian mechanics on glued symplectic pseudomanifolds // Demonstratio mathematica. 2008. Vol.41, N4. P. 941-960. https://doi.org/10.1515/dema-2013-0122
7. Watts J. Diffeologies, Differential Spaces, and Symplectic Geometry. Ph.D thesis. Department of Mathematics University of Toronto. 2012.
8. Виноградов А. М., Кальницкий В. С. Принцип наблюдаемости в примерах и задачах. СПб.: Изд. дом СПбГУ, 2012.
9. Кальницкий В. С. Симметрии плоской алгебры косимволов дифференциальных операторов // Зап. научн. сем. ПОМИ. 2016. Т. 443. С. 95-105.
10. Кальницкий В. С. Алгебра обобщенных полей Якоби // Зап. научн. сем. ПОМИ. 1995. Т. 231. С. 222-224.
11. Кальницкий В. С. Алгебры Якоби плоских многообразий // Зап. научн. сем. ПОМИ. 2003. Т. 299. С. 152-161.
Статья поступила в редакцию 25 апреля 2017 г.; рекомендована в печать 22 июня 2017 г. Сведения об авторе
Бурьян Сергей Николаевич — аспирант; burianserg@yandex.ru
BEHAVIOUR OF THE PENDULUM WITH A SINGULAR CONFIGURATION SPACE
Sergey N. Burian
St. Petersburg State University, Universitetskaya nab., 7-9, St. Petersburg, 199034, Russian Federation; burianserg@yandex.ru
The flat double mathematical pendulum is considered, the loose end of which moves along an ellipse. In the general case, the configuration space of this pendulum is two disjoint curves. It is possible to choose parametres so that these curves intersect transversally. The observed trajectory of motion in this case forms an angle. Moreover, there are special parameters in which the curves have a first-order tangency. In this case, there is a geometric uncertainty — how the pendulum have to move after passing a singular point.
It is shown that for the transversal case the inverse dynamic problem is not solvable, and the Lagrange multipliers tend to infinity as they move to a singular point in the configuration space. The observed motion is dynamically determined. The pendulum always moves from one branch of movement to another during the passage of a singular point. A qualitative explanation of this effect is proposed. Refs 11. Figs 3. Keywords: configuration spaces with singularities, Lagrange multipliers, double pendulum.
References
1. Zakalyukin I. V., "Singularities of degeneration of non-holonomic constraints and controllability", Trudy MAI (electronic journal), issue 39 (2010) [in Russian].
2. Zergda S. A., Soltakhanov Sh. Kh., Yushkov M.P., Equations of motion of nonholonomic systems and variational principles of mechanics. A new class of control problems (FIZMATLIT Publ., Moscow, 2005) [in Russian].
3. Polyakhov N.N., Zegzhda S. A., Yushkov M.P., Theoretical mechanics (Leningr. Univ. Press, Leningrad, 1985; Vysshaya Shkola Publ., Moscow, 2000) [in Russian].
4. Zhuravlev V. F., "The notion of connection in analytic mechanics", Nonlinear dynamics, 8(4), 853-860 (2012) [in Russian].
5. Sniatycki J., "Orbit of families of vector field on subcartensian space", arXiv:math/0211212 [math.DG].
6. Batubenge T. A., Sasin W., "An approach to hamiltonian mechanics on glued symplectic pseudomanifolds" Demonstratio mathematica 41(4), 941-960 (2008). https://doi.org/10.1515/dema-2013-0122
7. Watts J., Diffeologies, Differential Spaces, and Symplectic Geometry (Ph.D thesis, Department of Mathematics University of Toronto, 2012).
8. Vinogradov A. M., Kalnitsky V. S., The principle of observability in examples and problems (Ed. House of St. Petersburg State University, St. Petersburg, 2012) [in Russian].
9. Kalnitsky V. S., "Symmetries of a flat cosymbol algebra of differential operators", Zap. Nauchn. Sem. POMI 443, 95-105 (2016) [in Russian].
10. Kalnitsky V. S., "The algebra of generalized Jacobi fields", J. Math. Sci. 91, issue 6, 3476-3491 (1998). https://doi.org/10.1007/BF02434926
11. Kalnitsky V. S., "Jacobi Algebras of flat manifolds", J. Math. Sci. 131, issue 1, 5345-5350 (2005). https://doi.org/10.1007/s10958-005-0406-6
Для цитирования: Бурьян С. Н. Особенности движения маятника с сингулярным конфигурационным пространством // Вестник СПбГУ. Математика. Механика. Астрономия. 2017. Т. 4(62). Вып. 4. С. 541-551. https://doi.org/10.21638/11701/spbu01.2017.402
For citation: Burian S. N. Behaviour of the pendulum with a singular configuration space. Vestnik SPbSU. Mathematics. Mechanics. Astronomy, 2017, vol. 4(62), issue 4, pp. 541-551. https://doi.org/10.21638/11701/spbu01.2017.402