CC BY
54
УДК 531.011
Вестник СПбГУ. Сер. 1. Т. 3(61). 2016. Вып. 4
ГАШЕНИЕ КОЛЕБАНИЙ ТЕЛЕЖКИ С ДВОЙНЫМ МАЯТНИКОМ С ПОМОЩЬЮ УПРАВЛЕНИЯ ЕЕ УСКОРЕНИЕМ
, Е. А. Шатров, М. П. Юшков
Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
Задача о переводе за заданное время механической системы из одного фазового состояния в другое является одной из важнейших проблем теории управления. Модельной задачей в этом случае является нахождение оптимальной управляющей силы, переводящей горизонтально движущуюся тележку с маятниками, например, из состояния покоя за заданное время на заданное расстояние в новое состояние покоя. В предыдущих работах авторов было показано, что при решении такой задачи с помощью принципа максимума Понтрягина с минимизацией функционала от квадрата управляющей силы автоматически выполняется связь высокого порядка (например, связь восьмого порядка при движении тележки с двумя маятниками). Поэтому для решения этой же задачи был использован обобщенный принцип Гаусса, что позволило найти управляющую силу в виде полинома. В настоящей статье с помощью того же принципа решается задача о гашении колебаний тележки с двойным маятником. Предлагается предварительно вместо силы искать в качестве управления ускорение тележки, а затем по найденному закону изменения оптимального ускорения тележки отыскивать непосредственно управляющую силу и движение всей механической системы. Библиогр. 4 назв. Ил. 2.
Ключевые слова: управление движением, управляющая сила, принцип максимума Понт-рягина, обобщенный принцип Гаусса.
С. А. Зегжда
Постановка задачи и размерные уравнения движения. Одной из важнейших задач теории управления является задача о переводе механической системы за заданное время из одного фазового состояния в другое. Весьма часто для решения таких задач применяются метод динамического программирования [1] и метод, опирающийся на использование принципа максимума Понтрягина [2]. Модельной задачей среди подобных задач можно считать нахождение горизонтальной оптимальной управляющей силы, гасящей колебания тележки с маятниками при переводе этой механической системы за заданное время из одного состояния покоя в другое. В ряде предыдущих работ авторов (см., например, [3]) при решении таких задач с помощью принципа максимума Понтрягина с минимизацией функционала от квадрата управляющей силы [4] показывалось, что при этом движение механической системы оказывается подчиненным выполнению неголономной связи высокого порядка (при движении тележки с двумя маятниками выполняется связь восьмого порядка). Поэтому для решения той же задачи о гашении колебаний было предложено использовать обобщенный принцип Гаусса, свойственный движению неголономных систем со связями высокого порядка. Это позволило найти управляющую силу в виде полинома, при этом движение системы оказалось более плавным, чем движение, полученное с помощью принципа максимума Понтрягина.
В предлагаемой статье развиваются идеи применения обобщенного принципа Гаусса для гашения колебаний. Рассматривается система, состоящая из тележки с двойным маятником, причем используется новый подход для отыскания управляющей силы путем предварительного определения оптимального управляющего ускорения тележки.
(¡5 Санкт-Петербургский государственный университет, 2016 БОТ: 10.21638/11701/зрЪи01.2016.418
Итак, ставится следующая задача теории управления.
Рассмотрим горизонтальное движение тележки массы т вдоль оси х. К тележке прикреплена ось двойного маятника с длинами 1а и массами та, а = 1, 2 (рис.1). Требуется сформировать такую оптимальную горизонтальную силу Е, приложенную к тележке, которая за заданное время Т переместит тележку на расстояние причем вся механическая система должна перейти из первоначального состояния покоя в новое состояние равновесия.
Рис. 1. Движение тележки с двойным маятником
Обозначим углы отклонений маятников от вертикали через и у>2. Тогда при малых колебаниях уравнения Лагранжа второго рода запишутся следующим образом:
(т-1 + т-2 + т) X - (т-1 + т2) 1\ф\ - т2Ьф2 = Е, т111(111р1 - X) + т2Ь(11ф1 + 12Ф2 - X) = -(т1 + т2) §¿1^1 , (1)
т2^2(^2^2 + ¿1^1 - X) = -т212ё¥2 ■
Из второго и третьего уравнений системы (1) следует, что определение силы Е, при которой осуществляется переход системы из состояния покоя в новое состояние покоя, сводится к определению ускорения тележки, так как затем, используя первое уравнение, можно определить и силу Е. Ниже покажем, что задача об определении ускорения тележки с двойным маятником сводится к задаче об определении ускорения тележки с двумя независимыми математическими маятниками. Этот переход осуществляется на основе определения собственных форм и собственных частот колебаний двойного маятника.
Определение собственных частот и собственных форм колебаний двойного маятника. Введем обозначения
2 §
М = ш \ + т о + т, к" = — ,
¿1
а также безразмерное перемещение тележки
х
Тогда второе и третье уравнения системы (1) перепишутся в виде
ф1 + авфъ + к2ф1 = ж ,
2 — (2) ф1 + аф>2 + к ф2 = ж .
Здесь а и ¡3 таковы:
¡2 о т2
а = — , р
¡1 Ш1 + Ш2
Определим собственные частоты и собственные формы колебаний двойного маятника. Решение однородной системы, соответствующей системе неоднородных дифференциальных уравнений (2), будем искать в виде
фа = вш(Ш + 6), а = 1, 2 ,
где П — искомая размерная собственная частота. Теперь для определения постоянных и В2 получим систему
(1 - Л2) Б1 - а^Л2 В2 =0 ,
-Л2 Б1 + (1 - аЛ2) Б2 =0 , Л
Приравнивая к нулю определитель этой системы, находим собственные частоты механической системы:
02 _ л2-2 л2 _ 1 + аТу/«2~2а + 4а/3+1
— ЛД > К —-77/-ж-' ^ —
2(а - ар)
Постоянные Па, а =1, 2, соответствующие собственному значению Л2, V =1, 2, обозначим как Оиа. Как известно, они пропорциональны алгебраическим дополнениям элементов последней линейно зависимой строки определителя, когда в него подставлена величина Л^. Их значения вычисляются по формулам
= аЛ , Д^2 = (1 - Л1), V =1, 2 . (3)
Переход от системы с двойным маятником к системе с двумя независимыми маятниками. Согласно общей теории малых колебаний, собственные векторы, задаваемые выражениями (3), позволяют связать координаты ф1 и ф2 с главными координатами £1 и £2 следующим образом:
2
фа , а =1, 2 . (4)
^=1
Подставляя (4) в уравнения (2), придем к следующей системе уравнений относительно £1, £2:
2Д
Т " (5)
Рассматривая (5) как систему двух линейных алгебраических уравнений относительно у = £1 + и х = £2 + ^2^2, получим
а + = А1 х, £2 + ^£2 = А2 х,
где
а2(Д22 - А21) л2 (ди - Д12)
Л-1 = --:-:-:-, Ло =
(6)
Ац Д22 - Д12Д21 ' Д11Д22 - Д12 Д21
Перепишем систему (6) в виде
а + А^ а =Агё, ¿1
+ £2 =А2Ш. ¿1
(7)
Умножим уравнения системы (7) соответственно на 1/А1 и 1/А2 и введем новые переменные 'ф 1, Ф2 по формулам
А1 А2
Тогда система (7) перепишется в виде
1 У , §
■1р 1 = X ,
(9)
■цф1 + —ф1=Х,
1 - ё
ТТ 'Ф2 + Т 'Ф2 = X .
А2 ¿1
Пусть теперь к той же самой тележке подвешен не двойной маятник, а два независимых маятника, имеющих длины ¡1 и ¿|. Углы отклонения этих маятников обозначим соответственно через Ф1 и '2. Тогда, если перемещение тележки характеризуется координатой х, эти углы будут удовлетворять уравнениям
(10)
¿1 ф1 + § ф1 = X, ¿2, ф2 + § Ф2 = X ■ Видим, что при выполнении соотношений
С = 0"= 1,2,
система (9) совпадает с системой (10).
Таким образом, задача об определении ускорения тележки с двойным маятником, которое входит во второе и третье уравнения системы (1), эквивалентна задаче об определении ускорения, входящего в систему (10).
Как следует из выражений (4) и (8), искомые углы у>1, ^2 двойного маятника связаны с координатами '1, Ф2 следующими соотношениями:
¥>1 = Ац£1 + А21£2 = а/ЗА1 '1 + авА2ф2 , ^2 = А12£1 + Д22£2 = (^2 ~ А+ (^2 - ^2 . ^^
Как и в статье [3], размерное управление и (в нашем случае управлением является ускорение ж), входящее в систему (10), в соответствии с применением обобщенного принципа Гаусса будем искать в виде полинома, зависящего от времени. Однако этот полином по предложению П. Е. Товстика будем представлять в виде
6
и = ж = ^ вкгк (т - г), (12)
к=1
где Т — время перемещения тележки на расстояние Б из начального состояния покоя в конечное состояние покоя. Постоянные вк определяются из граничных условий
ж(Т) = Б, ж(?) = 0 , фа(Т)=0 , фа(Т)=0 , а = 1, 2 .
Интегрируя систему (10) при найденном управлении (12), получим углы Ф1 и Ф2 как функции времени. Затем по формулам (11) найдем и углы ф1 и ф2 как функции времени. В результате, используя первое уравнение системы (1), найдем силу Е, обеспечивающую переход системы из состояния покоя в новое состояние покоя при перемещении тележки за время Т на расстояние Б. На рис. 2 приведена зависимость силы Е, выраженной в долях ш^, от безразмерного времени т = г/Т. Исходные параметры двойного маятника таковы: Ш1 = Ш2 = ш, ¡1 = ¡2. Расчеты соответствуют перемещению тележки на расстояние Б = ¡1 за время Т = 1.5 (2п/П1), где П1 = у/2- ^^Ё/к-
Литература
1. Беллман Р. Динамическое программирование. М.: Изд-во иностр. лит., 1960. 400 с.
2. Понтрягин Л. С., Болтянский В. Г., Гамкрелидзе Р. В., Мищенко Е.Ф. Математическая теория оптимальных процессов. М.: Наука, 1983. 392 с.
3. Юшков М. П., Зегжда С. А., Солтаганов Ш. Х., Пашкина А. А. О связи теории управления с неголономной механикой // Вестн. С.-Петерб. ун-та. Сер. 1. Математика. Механика. Астрономия. 2014. Т. 1(59). Вып. 4. С. 15-23.
4. Черноусько Ф. Л., Акуленко Л. Д., Соколов Б.Н. Управление колебаниями. М.: Наука, 1980. 384 с.
Статья поступила в редакцию 30 января 2016 г.
Зегжда Сергей Андреевич —доктор физико-математических наук, профессор (1935—2015) Шатров Егор Александрович — аспирант; [email protected]
Юшков Михаил Петрович —доктор физико-математических наук, профессор; [email protected] SUPPRESSION OF OSCILLATION OF A TROLLEY WITH
A DOUBLE PENDULUM BY MEANS OF CONTROL OF ITS ACCELERATION
Sergey A. Zegzhda^, Egor A. Shatrov, Mikhail P. Yushkov
St. Petersburg State University, Universitetskaya nab., 7—9, St. Petersburg, 199034, Russian Federation; [email protected], [email protected]
The problem of transition of a mechanical system from one phase state to another is one of the most important problems of the control theory. In this case, a model problem consists in finding the optimal control force which transports the trolley with pendulums moving horizontally, for example, from a state of rest to a new state of rest over a given distance during the fixed time. In their previous papers the authors have shown that when solving such a problem with the help of the Pontryagin maximum principle with minimization of the functional of the control force squared, a high-order constraint is realized automatically (for instance, an eighth-order constraint for the motion of a trolley with two pendulums). That is why, for solving the same problem the generalized Gauss principle has been used that made it possible to find the control force as a polynomial. In the present paper the problem of suppression of oscillation of a trolley with a double pendulum is solved by means of the same principle. It is offered first to find the acceleration of the trolley as a control instead of the force, and then to seek immediately the control force by the obtained law of variation of the optimal acceleration of the trolley. Refs 4. Figs 2.
Keywords: control of motion, control force, Pontryagin maximum principle, generalized Gauss principle.
References
1. Bellman R., Dynamic programming (Princeton Univ. Press, 2010, 392 p.).
2. Pontryagin L. S., Boltyanskiy V. G., Gamkrelidze R. V., Mishchenko E. F., Mathematical theory of optimal processes (Nauka, Moscow, 1983, 392p.) [in Russian].
3. Yushkov M. P., Zegzhda S. A., Soltakhanov Sh. Kh., Pashkina A. A., "On relationship between the control theory and nonholonomic mechanics", Vestnik St. Petersburg University: Mathematics 47, issue 4, 181-188 (2014).
4. Chernous'ko F. L., Akulenko L. D., Sokolov B.N., Control of oscillation (Nauka, Moscow, 1980, 384 p.) [in Russian].
Для цитирования: Зегжда С. А., Шатров Е. А., Юшков М. П. Гашение колебаний тележки с двойным маятником с помощью управления ее ускорением // Вестник Санкт-Петербургского университета. Серия 1. Математика. Механика. Астрономия. 2016. Т. 3(61). Вып. 4. С. 683—688. DOI: 10.21638/11701/spbu01.2016.418
For citation: Zegzhda S.A., Shatrov E. A., Yushkov M. P. Suppression of oscillation of a trolley with a double pendulum by means of control of its acceleration. Vestnik of Saint Petersburg University. Series 1. Mathematics. Mechanics. Astronomy, 2016, vol. 3(61), issue 4, pp. 683—688. DOI: 10.21638/11701/spbu01.2016.418
