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

УДК 531.011 Вестник СПбГУ. Сер. 1. Т. 3(61). 2016. Вып. 2

НОВЫЙ ПОДХОД К НАХОЖДЕНИЮ УПРАВЛЕНИЯ, ПЕРЕВОДЯЩЕГО СИСТЕМУ

ИЗ ОДНОГО ФАЗОВОГО СОСТОЯНИЯ В ДРУГОЕ

, Е. А. Шатров, М. П. Юшков

Санкт-Петербургский государственный университет,

Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7-9

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

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

Ключевые слова: управление движением, управляющая сила, принцип максимума Понт-рягина, обобщенный принцип Гаусса, расширенная краевая задача.

С. А. Зегжда

1. Введение. Одной из центральных задач теории управления является задача об оптимальном переходе механической системы за заданное время из одного фазового состояния, в котором заданы обобщенные координаты и скорости, в другое фазовое состояние, в котором система должна иметь требуемые обобщенные координаты и скорости. В частности, если рассматриваемая механическая система в исходном и конечном состояниях покоится, говорят о гашении колебаний. Среди методов оптимального решения таких задач Н. Н. Моисеев [1] выделяет принцип максимума Понтрягина [2] и теорию локальных экстремумов Ф. Л. Черноусько [3]. К ним полезно добавить и метод динамического программирования Р. Беллмана [4].

(¡5 Санкт-Петербургский государственный университет, 2016

Модельной для исследования этой проблемы можно считать задачу о гашении колебаний горизонтально перемещающейся тележки с в маятниками, которая за заданное время должна пройти заданный путь в предположении, что в начале и в конце пути углы поворота маятников, их скорости, а также скорость тележки равны нулю. В ряде предыдущих работ авторов для решения подобных задач было предложено использовать аппарат неголономной механики со связями высокого порядка. Дело в том, что при решении такой задачи классическим методом, использующим принцип максимума Понтрягина [2] при минимизации квадрата управляющей силы [5], было замечено, что в процессе движения тележки с в маятниками непрерывно выполняется неголономная связь порядка 2в + 4. Это наталкивает на мысль о том, что для решения подобной задачи целесообразно использовать теорию движения него-лономных систем со связями высокого порядка. Применение обобщенного принципа Гаусса [6], свойственного этой теории, позволило добиться более плавного движения, чем движение, полученное при использовании принципа максимума Понтрягина. Это объясняется тем, что обобщенный принцип Гаусса позволяет находить управление в виде полинома, в то время как принцип максимума Понтрягина дает управление, содержащее гармоники с собственными частотами системы. При длительном движении это вводит систему в резонанс. Помимо этого обобщенный принцип Гаусса позволяет ставить и решать расширенную краевую задачу, когда указываются значения в начале и в конце движения не только обобщенных координат и скоростей, но и про-зводных от обобщенных координат по времени любого порядка. В частности, в случае приравнивания нулю ускорений в начале и в конце пути движения тележки удается найти управление без скачков этой функции. Эти результаты были подытожены в статье [7].

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

Отметим, что к идее отыскания управления в виде полинома, опираясь на метод интегродифференциальных соотношений, пришли также Г. В. Костин и В. В. Саурин [8, 9]. Помимо них подобными задачами активно занимается и М. А. Чуев [10].

2. Постановка задачи и классический подход к решению задачи о гашении колебаний. Рассмотрим горизонтальное движение тележки массы т, несущей два маятника, массы которых тх, т-2 и длины ¡1, ¡2 соответственно (рис.1). Для описания движения системы введем горизонтальное перемещение тележки х и углы поворота маятников ¿1 и ¿2. Решается задача о гашении колебаний, когда отыскивается горизонтальная управляющая сила Е, перемещающая за заданное время Т тележку на расстояние Б, причем механическая система должна перейти из первоначального состояния покоя в новое состояние покоя. Для этого должны выполняться следующие краевые условия:

¿х(0) = У1(Т)= 0 , ^51(0) = Т) = 0 ,

Ы0)= Ы т)=0 , ¿2

х(0) = Х(0) = Х( Т)=0 , х( Т) = Б.

Рис. 1. Тележка с двумя маятниками.

Дифференциальные уравнения движения системы в линейной постановке имеют вид (т1 + Ш2 + т)Х - ш^Ьф! - тъ^фъ = ^,

¡1ф1 + ёф1 = X ,

¡2ф2 + gф2 = X ,

где g — ускорение силы тяжести.

Ясно, что первое уравнение выражает закон движения центра масс рассматриваемой механической системы:

„ т111ф1 т212ф2 ,,

Мхг = г , хг = х--—---—— , М = т\ + пго + т .

М М

В свою очередь второе и третье уравнения системы (2.2) можно записать в виде

(1 - в)ф>1 - ^аф>2 + к2ф1 = -вф1 + а(1 - 7)ф2 + к2ф2

= Т/7"

где введены обозначения

а = — , /3: 11

т1

1

М

¡1

Главные формы колебаний и собственные частоты для однородной системы, соответствующей неоднородным уравнениям (2.3), получаем из системы

(,к2 - (1 - ДО2) С1 + 7аП2С2 = 0 ,

вп2с\ + (к2 - а(1 - 7)П2) С2 =0 .

Собственные частоты можно представить в виде

% = к\3 , з = 1, 2 , Л2_2 =

1 + а - ¡3 «7-^/

Обозначая главные координаты для системы (2.3) через Х1, Х2, выразим углы отклонений маятников через них следующим образом:

1

Ф1 =

0,1 - 0,2

1

ф2 =

а7

((1 + 02) (1 - (1 - ¡3)Л\) Х1 - (1 + 01) (1 - (1 - в)Л2) Х2) .

т2 ,2 g

Здесь

Вводя еще одну главную координату х0 и безразмерное время т по формулам

хо = > т = П ^, ¡1

перепишем дифференциальные уравнения движения в главных координатах:

х0

2

х'3 + ^¿Хз = и, ] = 1, 2

Е П,

Ш-1 =

Здесь штрихами обозначаются производные по безразмерному времени т. Систему дифференциальных уравнений (2.4) будем решать при краевых условиях

х0(0) = х0(0) = х0'(0)=0 , х,(0) = х,(0)=0 , Т = Т,

Б , (2.5) х0(Т) = а=—, х'0(Т) = Хо(Т) = 0, х^Т) = х'^Т) = 0, .7 = 1,2,

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

Система (2.4) является системой дифференциальных уравнений шестого порядка, а формулами (2.5) ставятся четырнадцать краевых условий. Поэтому согласно теории, изложенной в статье [7], требуется использовать обобщенный принцип Гаусса восьмого порядка, который обеспечивает нахождение безразмерного управления и из дифференциального уравнения восьмого порядка

(8) п и = 0 .

По рекомендации П. Е. Товстика будем отыскивать управление и в виде

6

и(т) = ^ с, ¿=1

*(т) = Е с,т3 (Т - т),

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

х0(т)= / и(тх)(т - тх) йтх, )0

х^(т) = — [ и(т1)8тш3-(г -Т1)д,Т1 , ^ = 1,2-^^ 0

^^ 0

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

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

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

При такой постановке задачи вместо системы дифференциальных уравнений (2.2) естественным образом записываем систему уравнений

х = и,

¡1Ф1 + g^l = и , ¡2^2 + g^2 = и ,

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

Если теперь ввести безразмерные координаты, управление и время по формулам

хо =

Х1 = у>1, Х2 = а.^2

и

система уравнений (3.1) запишется в виде

Х1 + Х1 = иФ

х'п Л--Х2 = и* .

а

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

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

= МgUф - ш^х'/ - m2gx'J¡

х

¡1

g

и

4=

а соответствующее ей новое безразмерное управление иN, выраженное в долях Мg, запишется в виде

= тТ~ = и* С1 ~ Р ~ 7) + Р + — х2  М g а

4. Построение аналитического решения. Итак, получили безразмерное управление. Однако оказывается, что искомое ускорение тележки в безразмерных переменных при заданных длинах и массах маятников зависит от одного безразмерного параметра К = Т/Т1, т.е. от отношения безразмерного времени перемещения тележки к безразмерному периоду колебаний наиболее длинного маятника. В работе [11] показано, что при любом числе маятников имеется счетное число особых значений параметра К, при приближении к которым ускорение тележки неограниченно возрастает.

Определить управление движением тележки, которое является аналитической функцией времени при всех значениях параметра К, предлагается следующим образом. Строится новое решение рассматриваемой задачи. При его построении дополнительно предполагается, что в начале и в конце пути производная от ускорения тележки по времени равна нулю. У нового решения задачи также будут существовать особые значения параметра К, но они будут отличны от особых значений исходного решения. Пусть управление и1(т) соответствует первому решению, а и2(т) —второму. Тогда при любых значениях параметра ¡1 функция

и

будет решением рассматриваемой задачи. Избежать вычисления особых значений решений и1 и и2 и построить аналитическое решение, непрерывно зависящее от параметра К, позволяет определение параметра ¡ из условия минимальности интеграла от квадрата функции и(т) за время перемещения Т. Это решение соответствует следующему значению параметра ¡ :

З1 - З2

где

^ 31 + З3 - 232 '

31 = и2(т) йт , З2 = и1(т )и2 (т) йт , 3з = и2(т) ¿т. ио Jо и о

5. Результаты численных расчетов. Результаты численных расчетов, полученные по старому и по новому подходам для движения тележки с маятниками, сравнивались между собой. Вычислялись как функции времени управляющая сила Е и углы поворотов маятников. Сила выражалась в долях Мg, где М — масса всей системы, а углы поворотов ф1 и ф2 — в градусах. Как отмечалось выше, данные функции, определенные как функции безразмерного времени т = Ь/Т, существенно зависят от параметра К, равного отношению безразмерного времени перемещения Т к периоду Т1 первой формы колебаний механической системы. Интересным является случай, когда величина К находится в промежутке от единицы до двух [11]. С учетом этого при расчетах полагалось, что К=1.54. Так как отыскиваемые функции пропорциональны перемещению тележки Б за время Т, данное перемещение, заданное в долях длины

-0.005 --0.01 -0.015

Рис. 2. Управляющая сила для тележки с одним маятником.

Рис. 3. Колебания маятника на тележке.

Рис. 4. Управляющая сила для тележки с двумя маятниками (случай т^/тх = 1/8).

Рис. 5. Колебания первого маятника.

Рис. 6. Колебания второго маятника.

Ме 0.015

0.01 0.005

Рис. 7. Управляющая сила для тележки с двумя маятниками (случай т2 /т\ = 1/64).

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

Рис. 2 и 3 соответствуют расчету, когда на тележке укреплен один маятник, при этом выбирались следующие параметры системы:

ш 1 ¡

— = -, А =1.54, 5'=-.

Ш1 2 5

Рис. 4-7 соответствуют движению тележки с двумя маятниками. В этом случае имеем четыре независимых параметра: К, отношения массы тележки т к массе т1 первого маятника и к массе т2 второго маятника, а также отношение длины ¡2 второго маятника к длине ¡1 первого. Рис. 4-6 соответствуют параметрам системы

£ = 1, ^ = - = -, К = 1.54, 5=^1, ¡1 4 т1 8 т1 2 5

а для рис. 7 принято

Ш2 _ т1 64

Расчеты, как и следовало ожидать, показали, что чем больше масса тележки по отношению к массе маятников, тем ближе результаты, полученные по первому и по второму подходам. Поэтому полагалось, что масса тележки в два раза меньше массы первого основного маятника. Длина же второго маятника и его масса по отношению к первому принимались меньшими соответственно в четыре и в восемь раз. Перемещение тележки равно 0.2 ¡1, где ¡1 —длина основного маятника. Видим, что применение нового подхода дает уменьшение колебаний маятников и по амплитуде, и по частоте. Из сравнения рис. 4 и 2 следует, что при новом подходе добавление второго маятника незначительно повлияло на управляющую силу, в то время как при обычном подходе она существенно изменилась. Интересным является и следующий результат. Пусть в задаче с двумя маятниками масса второго меньше первого не в восемь, а в шестьдесят четыре раза. Сравнивая рис. 7 и 4 видим, что при новом подходе такое существенное уменьшение массы второго маятника мало повлияло на управляющую силу, в то время как при обычном подходе она изменилась заметно.

Литература

1. Моисеев Н. Н., Иванилов Ю. П., Столярова Е. М. Методы оптимизации. М.: Наука, 1978. 352 с.

2. Понтрягин Л. С., Болтянский В. Г., Гамкрелидзе Р. В., Мищенко Е.Ф. Математическая теория оптимальных процессов. М.: Наука, 1983. 392с.

3. Черноусько Ф. Л., Баничук Н. В. Вариационные задачи механики и управления. Численные методы. М.: Наука, 1973. 238 с.

4. Беллман Р. Динамическое програмирование. М.: Изд-во иностр. лит., 1960. 400 с.

5. Черноусько Ф. Л., Акуленко Л. Д., Соколов Б. Н. Управление колебаниями. М.: Наука, 1980. 384 с.

6. Поляхов Н.Н., Зегжда С. А., Юшков М.П. Обобщение принципа Гаусса на случай неголо-номных систем высших порядков // Докл. АН СССР. 1983. Т. 269, №6. С. 1328-1330.

7. Юшков М. П., Зегжда С. А., Солтаханов Ш. Х., Пашкина А. А. О связи теории управления с неголономной механикой // Вестн. С.-Петерб. ун-та. Сер. 1. 2014. Т. 1(59), вып. 4. С. 609-618.

8. Костин Г. В., Саурин В. В. Интегродифференциальный подход к решению задач линейной теории упругости // Доклады академии наук. 2005. Т. 404, №5. С. 628-631.

9. Костин Г. В., Саурин В. В. Моделирование и оптимизация движений упругих систем методом интегродифференциальных соотношений // Доклады академии наук. 2006. Т. 408, № 6. С. 750-753.

10. Чуев М.А. Дифференциальные уравнения програмных движений механической системы // Изв. РАН. 2008. №1. С. 179-192.

11. Гаврилов Д.Н., Зегжда С. А. Гашение колебаний упругого тела при его перемещении // Вестн. С.-Петерб. ун-та. Сер. 1. 2012. Вып. 3. С. 73-83.

Статья поступила в редакцию 8 декабря 2015 г. Сведения об авторах

Зегжда Сергей Андреевич —доктор физико-математических наук, профессор (1935—2015) Шатров Егор Александрович — аспирант; egorshatroff@yandex.ru

Юшков Михаил Петрович —доктор физико-математических наук, профессор; yushkovmp@mail.ru

A NEW APPROACH TO FINDING THE CONTROL TRANSPORTING A SYSTEM FROM ONE PHASE STATE TO ANOTHER

Sergey A. Zegzhda*, Egor A. Shatrov, Mikhail P. Yushkov

St. Petersburg State University, Universitetskaya nab., 7-9, St. Petersburg, 199034, Russian Federation; egorshatroff@yandex.ru, yushkovmp@mail.ru

In their previous papers the authors have considered a possibility of application of the theory of motion for nonholonomic systems with high-order constraints to solving one of the main problems of the control theory. This is a problem of transporting a mechanical system with the finite number of degrees of freedom from a given phase state to another given phase state during a fixed time. It was shown that when solving such a problem using the Pontryagin maximum principle with minimization of the integral of the control force squared, a nonholonomic high-order constraint is realized continuously during the motion of the system. But in this case, one can also apply a generalized Gauss principle, which is commonly used in the motion of nonholonomic systems with high-order constraints. It is essential that the latter principle makes it possible to find the control as a polynomial, while the use of the Pontryagin maximum principle yields the control containing harmonics with natural frequencies of the system. The latter fact determines increasing the amplitude of oscillation of the system if the time of motion is long. Besides this, a generalized Gauss principle allows us to formulate and solve extended boundary problems in which along with the conditions for generalized coordinates and velocities at the beginning and at the end of motion, the values of any-order derivatives of the coordinates are introduced at the same time instants. This makes it possible to find the control without jumps at the beginning and at the end of motion. The theory presented has been demonstrated when solving the problem of the control of horizontal motion of a trolley with pendulums. A similar problem can be considered as a model, since when the parameters are chosen correspondingly it becomes equivalent to the problem of suppression of oscillations of a given elastic body some cross-section of which should move by a given distance in a fixed time. Equivalence of these problems widens essentially a range of possible applications of the problem of a trolley with pendulums.

Earlier solving the problem has been reduced to the choice of a horizontal force that is a solution of the problem formulated. In the present paper it is offered to seek an acceleration of a trolley, with which it moves by a given distance in a fixed time, as a time function but not a force applied to the trolley, the velocities and accelerations being equal to zero at the beginning and at the end of motion. In this new problem the rotation angles of pendulums are the principal coordinates. This makes it possible to find a sought acceleration of a trolley on the basis of a generalized Gauss principle according to the technique developed before. Knowing the motion of a trolley and pendulums it is easy to determine the required control force. The results of numerical calculation are presented. Refs 11. Figs 7.

Keywords: control of motion, control force, Pontryagin maximum principle, generalized Gauss principle, extended boundary problem.

References

1. Moiseev N. N., Ivanilov Yu.P., Stolyarova E. M., Optimization methods (Nauka, Moscow, 1978, 352p.) [In Russian].

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. Chernous'ko F. L., Banichuk N. V., Variational problems in mechanics and control. Numerical methods (Nauka, Moscow, 1973, 238p.) [In Russian].

4. Bellman R., Dynamic programming (Princeton Univ. Press, 2010, 392p.).

5. Chernous'ko F. L., Akulenko L. D., Sokolov B.N., Control of oscillation (Nauka, Moscow, 1980, 384 p.) [In Russian].

6. Polyakhov N. N., Zegzhda S. A., Yushkov M. P., "Generalization of the Gauss principle to the case of nonholonomic high-order systems", Dokl. AN. SSSR 269(6), 1328-1330 (1983) [In Russian].

7. Yushkov M.P., Zegzhda S. A., Soltakhanov Sh.Kh., Pashkina A. A., "On relationship between the control theory and nonholonomic mechanics", Vestnik Saint-Petersburg University. Mathematics 47, Issue 4, 181-188 (2014).

8. Kostin G. V., Saurin V. V., "Integrodifferential approach to solving problems of linear elasticity theory", Doklady Physics 50(10), 535-538 (2005).

9. Kostin G. V., Saurin V. V., "Modeling and optimization of elastic system motion by the method of integrodifferential relations", Doklady Mathematics 73(3), 469-472 (2006).

10. Chuev M.A., "Differential equations for program motion of a mechanical system", Mech. Solids 43(1), 153-164 (2008).

11. Zegzhda S. A., Gavrilov D. N., "Suppression of vibration of an elastic body when it moves", Vestn. St.Petersb. un-ta. Ser. 1, Issue 3, 73-83 (2012) [In Russian].