Управление системой нелинейно связанных перевернутых маятников Текст научной статьи по специальности «Математика»

Journal of Siberian Federal University. Engineering & Technologies, 2018, 11(3), 280-290

yflK 517.9

Nonlinear Related System of Inverted Pentulums Control

Mikhail E. Semenova, Mikhail A. Popovb and Olesya I. Kanishchevaa*

aMilitary Education and Research Centre of Military-Air Forces

«Military-Air Academy Named After Professor N.E. Zhukovsky and Yu.A. Gagarin» 54a Starykh Bol'shevikov Str., Voronezh, 394064, Russia bVoronezh State Technical University 84 20 let Oktyabrya Str., Voronezh, 394006, Russia

Received 05.12.2017, received in revised form 04.01.2018, accepted 21.04.2018

In this paper we propose a linear and nonlinear mathematical model of linked inverse pendulums. We investigate dynamics of this mechanical system and determined the stability parameters. After that we presented results of experiments for various system configurations. In conclusion we constructed stability zones in the parameter space for linear and nonlinear systems. The system is controlled by feedback. The introduced nonlinear spring stiffness is part of the dynamic control. There are three stationary points in the phase space, however, only one of them has real coordinates. As a result of the study, it was shown that a complex unstable system consisting of oscillators with a nonlinear coupling can be described by a fairly simple system of equations, and its stabilization, under certain conditions, is possible with the help of a fairly simple control of the periodic feedback function.

Keywords: inverted pendulum, linked oscillators, stabilization, controlling.

Citation: Semenov M.E., Popov M.A., Kanishcheva O.I. Nonlinear related system of inverted pentulums control, J. Sib. Fed. Univ. Eng. technol., 2018, 11(3), 280-290. DOI: 10.17516/1999-494X-0040.

© Siberian Federal University. All rights reserved

Corresponding author E-mail address: mkl150@mail.ru, soeltic@gmail.com, oleka_olesya@mail.ru

Управление системой нелинейно связанных перевернутых маятников

М.Е. Семёнов3, М.А. Попов®, О.И. Канищеваа

аВоенный учебно-научный центр Военно-воздушных сил

«Военно-воздушная академия имени профессора Н.Е. Жуковского и Ю.А. Гагарина» Россия, 371600, Воронеж, ул. Старых Большевиков, 54а бВоронежский государственный технический университет Россия, 394064, Воронеж, ул. 20-лет Октября, 84

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

Ключевые слова: обратный маятник, связанные осцилляторы, стабилизация, управление.

Введение

Теория колебаний нелинейных систем широко применяется при моделировании различных физических явлений: колебания молекул в жидкостях и твердых телах, электрических цепях, состоящих из нескольких взаимосвязанных контуров и т.д. [1] В таких системах реализуются разнообразные дисперсионные зависимости, на основе которых исследуется распространение волн в нелинейных средах. Подавляющее большинство моделей таких систем основывается на законах движения простейших связанных осцилляторов и их цепочек, динамика которых формализуется посредством как линейных, так и нелинейных уравнений. В большинстве подобных задач рассматриваются колебания маятников, у которых нулевое положение является устойчивым [2-4]. Подробный обзор последних результатов в этой области приведен, например, в [5].

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

ятника с помощью вертикальных осцилляций точки крепления хорошо изучена и сводится к уравнению Матье [8]. Теоретическое объяснение этого явления было сделано Стефенсоном в 1908 г. [9]. Физическое объяснение динамической стабилизации перевернутого маятника вертикальными осцилляциями точки крепления (с помощью управления по обратной связи) было предложено академиком П.Л. Капицей в 1951 г., выполнившим детальное экспериментальное исследование этого явления [10]. Также отметим работы, в которых стабилизация верхнего неустойчивого положения осуществлялась посредством гистерезисного управления [11-14].

Постановка задачи

В настоящей работе рассматривается система, состоящая из двух перевернутых маятников массами m¡ и m2, связанных между собой пружиной с жесткостью k. Основная задача работы заключается в формировании такого закона управления горизонтальным движением точки крепления, при котором фазовые координаты системы остаются в ограниченной области фазового пространства. Иными словами, речь идет об условиях, накладываемых на параметры системы и управление, обеспечивающих диссипативность системы. Рассмотрим решение данной задачи для различных вариантов управления рассматриваемой системой (управление может быть приложено как к одному маятнику системы, так и к системе в цел м) и в случае линейной и нелинейной жесткости системы.

Случай управления одним маятником

Рассмотрим систему с линейной жесткостью и управлением, приложенным к одному из маятников. Ее физическая модель представлена на рис. 1.

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

i Jcq> с = тс glc sin фс + kh2 (sin ф2 - sin фс)-и \j2 ф 2 = т2 gl2 sin ф2 - kh2 (sin ф2 - sin фс)

Рис. Fig.

1. Физическая модель связаниых маятников

1.Physicalmodelof coupled pendulums

где ф1, ф2 - углы отклонения маятников от вертикального положения. Предполагая отклонения от нулевого положениямалыми, линеаризуемуравнениядинамики.Тогда уравнения (1) примут вид

ф 1 =

^ mlgll kh2 ^

ф a =

mogl2 kh J

2 Л

V J 2

2 /

kh m,l, ..

Ф1+-J- ф a--p >"

Ji Ji

kh2 Ф 2+—Г Ф1

J n.

(2)

Для длпьнвпшдго рассмот^дия систем;^0 (2) удобно придставииь д форде четырех урав-иений пьрвдгв НРавеи^дил перемееные: 2 = ЧВ, 22 = Ч, 2 = Ф2, 2 = с^с2 и обозначим

кк2 теI

В = ——, Л = — В. С учетом введенных обозначений динамика системы будет описываться

следующими уравнениями:

z2 = A1 z1 + B1 z3 + z

(3)

+Z4 = B2Z1 + A2 Z3

или, is матричной форме,

f 0 1

и^де; Q =

z = Qz + bu, 0 0) 00)

(4)

4 0 B1 0 0 0 0 1

0 4 0

V B2

b =

1

0

V0/

. СИоамируим алноритм управленияпо принципу о^аепой

связи: м _ . _. arctan(z x с)' где ^ _ парамета, характеризующий «интенси-нонтьые управ-

ления. Отметим, что при отсутствии упраалиния (Л = 0) поведение системы будет определяться собственными числамваатрицы системы. Эта матрицт имеет четыре собственных числа:

ч*=±1

A1 + A2 + а (A -A2f + 4 / B2

4,4 =±1

. Два из них

вещественные (одно положительно, второе отрицательно), а два других - чисто мнимые. Очевидно, что фазовое пространство делится на два подпространства. Колебания в одном из них (которое соответствует комплексным собственным числам) будут соответствовать особой точке центр, а во втором подпространстве особая точка - седло [8].

Таким образом, задача стабилизации нулевого положения, соответствующего вертикальному положению обоих маятников, заключается в подавлении эффекта положительного собственного числа [14]. Такая система детально описана в работе [15], в которой было показано, что при достаточно большом значении параметра Б система (2) стабилизируема, - т.е. фазовые координаты остаются в ограниченной области фазового пространства. В настоящей рабо-

те рассматривается пружина с нелинейной жесткостью, е именно предполагается изменение жесткости по Зимину

к = к0 (1 + х2), (5 )

где х = к ( - р2) или, во введеныых выше обозннчениях, х = X=z1 - аз). Таким образом, жесткость пружины, зависящая от взаимных отклонений маятников, является частью динамического управления. Перепишем систему (2), подставив в нее (5):

( к) (1 + (Й ( - Z3 ))+ )( 0« f1 + (h (( - к ))2 ) h

Jj

J

Jl

К f1 ^ (h ( - .23 )) (m2gl2 k) -1 + ( ( - Z3 ))2 ) h

JJ

J

J,

mgl k0 hg

Сделаем замену переменных —j— = a, ■—— = p. Догда системе (T) примет вид

(6)

z2 = а Л-Pi (( - к ) - ) () - к К +г

73 =Z4

zA = a—z3 +132 ( - z3) + p2/2 ( - z3 )

Как известно, динимика лпобот сиитемысвязанг с ее поведе:нием а окреетностз стноитнмр-ных решоний [Т6[. 13 рмсс 15(^в:][)]м^аете^1н случме стационарные тзрки ынтедтлсются кал решения следующейсистемыалгебраическихуравнений:

z2 — О

а

1 zi -Pi(-i -k)-Piг ( -k3)3 + u —(

z4 — О

а2Z3 + Р2 (l - k ) - Р2) ( - k ) — 0

(7 )

Тогда в фазовом пространстве существует три точки равновесия, из которых вещественные координатыимееттолькоодна:

P:

1

(

Z 6 Т\к2

-6Г

л

z? = О

оцосДРСг^ЗДёгj -6иГ2Г3 -оф! )-Г/И(х 1(3.

Г?

2 3

6 ot2 |31/г2Г2Г3

Z4 = О

Г1 =а1а2 -а1Р2 -а2Рь Г2 = а2Ь1 + а1Ь2,

Г3 "Г52 + ^3а®р(й6Г+ (( + 27а22И2и2Г2) .

Исследуем устойчивость этого единственного стационарного решения. Для этого применим первый метод Ляпунова [17] по линейному приближению. Линеаризуем исходную систему в ощэестности стационарногорешения и найдем собственные числа матрицы линеаризованной системы:

с1еЮ (в-Я0ч) = 0 с (8)

д

где © =—, г = (г2,23,). В рассматриваемом случае матрица линеаризованной системы будев нмете вид

0 =

О 10 0

0=-^= + 1X CX Р^ + с^х сд

0 0 0 1

Р2Ц 0 а2 — P^JJ. 0

(9)

2 —2—

где |и = 1 + 3/- (-—-Г3) , х =—т-2' П и^ ^ лоются соотношением Р1. Характери-

7= 1 + ) )

стическое утавнение(8) можно перепиеатьв следчющем виде:

X4 + у^3 и у2Х2 -4- -У У4 = <0.

Подставим вычисленные ранее координаты точки Р1 и найдем собственные значения Х1>2д4. Во избежание громоздких выкладок дальнейшие вычисления приведем для случая, в котором: величина управления Б = 10; начальная жесткость пружины к0 = 14.5; высота крепления пружины h = 0.1; длины и массы осцилляторов 11 = 0.1,12 = 0.15, т1 = 0.1, т2 = 0.2. Получим собственные значения в точке Л: ^ = -1.0127 + /0.7334, = -1.0127 - /0.7334, Х3 = -0.6582, = -0.2453. Поскольку все показатели Л < 0, то при данных параметрах система асимптотически устойчива [16].

Динамика системы при значениях параметров, описанных выше, и начальных условиях: ^[0]=-0.001, г2[0]=-0.001, ¿3[0] = 0.002, г4[0]=-0.002,изображенанарис.2, 3.

Как видим, стационарное решение, определяемое координатами точки Р1, устойчиво по Ляпунову. Исследуем зависимость устойчивости стационарного решения от параметров к0 и Б. Для этого зададим диапазоны изменения интенсивности управления Б е [1;50] и параметра жесткости пружины к0 е [1;20], а значения остальных параметров зафиксируем. Как известно из [16], стационарное решение устойчиво, если старший показатель Ляпунова отрицателен. Т.е. область устойчивости в пространстве параметров к0, Б определяется как решение неравенства Л = Лшах(к0, Б) < 0. Результаты представлены на рис. 4.

На графике рис. 4 отмечены области устойчивости системы в пространстве параметров. В областях, помеченных синим цветом, показатель Ляпунова Л > 0 - это значит, что система неустойчива. В белых же областях, наоборот, система асимптотически устойчива.

ZI

0.010520! O.OIC'ÜCf O.ilOiSCi IIHSIII 0.010520!

EI

-О 07411c -0.57411l

-0.074114

-O 174114

Рис. 2. Графики координат: z1 (слева) и z3 (справа) Fig. 2. Velocity curves: z1 (left) and z3 (right)

0.*10 4. О-1"

1 4 4 10" 1 *10" i H 10"

SOO 1000

tOO 400 ООО 300

Рис. 3. Графики скоростей: z2 (слева) и z4 (cnpaiBa) Fig.3. Speedgraphs: z2 (left)and z4 (right)

:

3D (

l>

40

20 0

Рис. 4. Распределение значений Re[max(X)] в зависимости от значений k0 и D Fig. 4. Values distribution of Re[max(X)] in terms of values of k0 and D

Случ ай управления обопмимаятниками

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

Вданном случае в четвервнмураенении системы(4) добавится переменная упуккления и. Как ре в предыдр щотсхдчае, имеетет топтаееднастационамнея тонко:

P

1

( бЩа2Р, + з/12г

18Г32 h2

# h2 Г

3

z2 = 0

а1а2Р1Р2h2 (3-^18Г3 - 6мГ^Г3) - -^18а11

18а2 h2Г32Г3

Z 4 = 0

Матрицалинеаризованнойсистемы имеетвид

© =

где ц = 1 + 3// 1 (zx ~zx) , у= =

0 10 0

с^-Р^ + сд XX Р^ + ^Х сд

0 0 0 1

Р2м- + ciX c2X a2 "Р2Ц + сзХ XX -2d

(10)

Zj и z3 определяются соотношением Pi. Характери-

71(1 + (г хеУ

стическое уравнениеэквивалентно поха'ченномуранее

(з4 + У!;^3+у.А.2 -I- у3з + ++ = о- .

Найдем собгтвенные значения =ы2,з,- Во чзбежатие громоздкик выкладок, даньней-ньиевыаисления приведем для случая, в котором начальная жесткость пружины к0 = 100, остальные параметры такие оке, как и в предыдущем случае. Получим собственные значения: ^ к -О.РОРН о 314.9709, и = -0.8794 - /10.Р009,13 у -0.5И05, Ч = -0.32Гб. Тоща пегазатели

z1 =

z, =

Рис. 5. Физическая модель свянанныхмаятников Fig. 5. Coupled pendulums physic almodel

Ляпунова: Л; = Л2 = ReX2 = -0.8794, Л3 = Х3 = -0.5705, Л4 = Х4 = -0.3231. Поскольку все показатели Л < 0, то при данных выбранных параметрах система устойчива [16]. Результаты моделирования описанной системы представлены на рис. 6, 7.

Как видно из графиков (рис. 6-7), при моделировании с заданными начальными условиями система диссипативна. Исследуем зависимость устойчивости стационарного решения от параметров к0 и Б. Для этого зададим диапазоны изменения интенсивности управления Б е [1;50] и параметра жесткости пружины к0 е [1;20], а значения остальных параметров зафиксируем. Результатыданы на рис.8.

На графике рис. 8 отмечены области устойчивости системы в пространстве параметров. В областях, помеченных синим цветом, показатель Ляпунова Л > 0 - это значит, что система неустойчива. Вбелых жеобластях,наоборот, системаустойчива.

Заключение

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

Рис. 6. Графики координат: z1 (слева) иг3 (справа) Fig. 6. Coordinate graphs: z1 (left) andz3 (right)

Рис. 7. Графики скоростей: z2 (слева) и z4(cnpaBa) Fig. 7. Velocity curves: z2 (left) and z4 (right)

5 1с 15 2d

Рис. 8. Распределение значений Re[тах(Х)] в зависимости отзначений к0иО Fig. 8. Values distribution ofR.e[max(?b)]interms of values of k0 andD

Работа выполнена при поддержке РФФИ, гранты № 17-01-00251, 16-08-00312.

Список литературы

[1] Андронов А.А., Витт А.А., Хайкин С.Э. Теория колебаний. М.: Наука, 1981. 568 с. [Andronov A.A., Witt A.A., Khaikin S.E. Theory of oscillations. Moscow, Nauka, 1981. 568 p. (in Russian)]

[2] Баутин Н.Н. О числе предельных циклов, появляющихся при изменении коэффициентов из состояния равновесия типа фокуса или центра, Математический сборник, 1952, 30(72), 181-196. [Bautin N.N. On the number of limit cycles appearing when the coefficients change from an equilibrium state of the focus or center type, Matematicheskii sbornik, 1952, 30(72), 181-196. (in Russian)]

[3] Трубецков Д.И., Рожнев А.Г. Линейные колебания и волны. М.: Издательство физико-математической литературы, 2001. 416 с. [Trubetskov D.I., Rozhnev A.G. Linear oscillations and waves. Moscow, Publishing House of Physical and Mathematical Literature, 2001. 416 p. (in Russian)]

[4] Магнус К. Колебания: Введение в исследование колебательных систем. Пер. с нем. М.: Мир, 1982. 304 с. [Magnus K. Oscillations: Introduction to research of oscillatory systems. Translation from German. Moscow, Mir, 1982. 304 p. (in Russian)]

[5] Осипов Г.В. Синхронизация в неоднородных сетях осцилляторов. Нижний Новгород, 2014. 135 с. [Osipov G.V. Synchronization in heterogeneous networks of oscillators. Nizhny Novgorod, 2014. 135 p. (in Russian)]

[6] Колмогоров А.Н. О сохранении условно периодических движений при малом изменении функции Гамильтона, Доклады АН СССР, 1954, 98(4), 527-530. [Kolmogorov A.N. On the conservation of conditionally periodic motions for a small change in the Hamiltonian function, Doklady Akad. Nauk USSR, 1954, 98(4), 527-530. (in Russian)]

[7] Бутиков Е.И. Стабилизация перевернутого маятника (60 лет маятнику Капицы), Компьютерные инструменты в образовании, 2010, 5, 3951. [Butikov E.I. Stabilization of the inverted pendulum (60 years of Kapitza's pendulum), Computer tools in education, 2010, 5, 39-51. (in Russian)]

[8] Неймарк Ю.И., Коган Н.Я., Савельев В.П. Динамические модели теории управления. М.: Наука, 1985. 400 с. [Neimark Yu.I., Kogan N.Ya., Savelyev V.P. Dynamic models of control theory. Moscow, Nauka, 1985. 400 p. (in Russian)]

[9] Stephenson A. «On an induced stability», Phil. Mag, 1908, 15(233).

[10] Капица П.Л. Маятник с вибрирующим подвесом, УФН, 1951, 44, 720. [Kapitsa P.L. A pendulum with a vibrating suspension, UFN, 1951, 44, 7-20. (in Russian)]

[11] Матвеев М.Г., Семенов М.Е., Шевлякова Д.В., Канищева О.И. Зоны устойчивости и периодические решения перевернутого маятника с гистерезисным управлением, Мехатрони-ка, Автоматизация, Управление, 2012, 11, 814. [Matveev M.G., Semenov M.E., Shevlyakova D.V., Kanishcheva O.I. Zone of Stability and Periodic Solutions of the Inverted Pendulum with Hysteretic Control, Mechatronics, Automation, Control, 2012, 11, 8-14. (in Russian)]

[12] Семенов М.Е., Хатиф З., Решетова О.О., Демчук А.А., Мелешенко П.А. Модель динамики обратного маятника с гистерезисным управлением, Вестник Воронежского государственного университета. Серия: Физика. Математика, 2016, 4, 165-177. [Semenov M.E., Hatif Z., Reshetova O.O., Demchuk A.A., Meleshenko P.A. The model of the dynamics with the inverted pendulum hysteretic control, Proceedings of Voronezh State University. Series: Physics. Mathematics,

2016, 4, 165-177. (in Russian)]

[13] Семенов М.Е., Матвеев М.Г., Лебедев Г.Н., Соловьев А.М. Стабилизация обратного гибкого маятника с гистерезисными свойствами, Мехатроника, Автоматизация, Управление,

2017, 8, 516-525. [Semenov M.E., Matveev M.G., Lebedev G.N., Solovyev A.M. Stabilization of a Flexible Inverted Pendulum with the Hysteretic Properties, Mechatronics, Automation, Control, 2017, 8, 516-525 (in Russian)]

[14] Solovyov A.M., Semenov M.E., Meleshenko P.A., Reshetova O.O., Popov M.A., Kabulova E.G. Hysteretic nonlinearity and unbounded solutions in oscillating systems, Procedia Engineering, 2017, 201, 549-555.

[15] Семенов М.Е., Соловьев А.М., Попов М.А. Стабилизация неустойчивых объектов: связанные осцилляторы, Труды МАИ, 2017, 93. [Semenov M.E., Solovyov A.M., Popov M.A. Stabilization of unstable objects: coupled oscillators, Trudy MAI, 2017, 93 (in Russian)]

[16] Кузнецов С.П. Динамический хаос. М.: Физматлит, 2006. 356 с. [Kuznetsov S.P. Dynamic chaos. Moscow, Fizmatlit, 2006. 356 p. (in Russian)]

[17] Мирошник И.В. Теория автоматического управления. Нелинейные и оптимальные системы. СПб.: Питер, 2006. 272 с. [Miroshnik I.V. Theory of automatic control. Nonlinear and optimal systems. St. Petersburg, Piter, 2006. 272 p. (in Russian)]