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

шилов В. С. Динамика космических систем с тросовыми и шарнирными соединениями. Ижевск: Регулярная и хаотическая динамика, 2007. 558 с.

С. П. Безгласный, Е. В. Куркина

Самарский государственный аэрокосмический университет им. акад. С.П. Королева, кафедра теоретической механики E-mail: bezglasnsp@rambler.ru, ekaterina.kurkina@mail.ru

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

Ключевые слова: гамильтониан, неавтономная система, программное движение, функция Ляпунова.

15. Пановко Я. Г. Введение в теорию механических колебаний. М.: Наука, 1991. 256 с.

Construction and Stabilization Program Motions of Nonautonomous Hamiltonian Systems

S. P. Bezglasnyi, E. V. Kurkina

Samara State Aerospace University

named after academician S.P. Korolyov,

Chair of Theoretical Mechanics

E-mail: bezglasnsp@rambler.ru, ekaterina.kurkina@mail.ru

We consider program motion of Hamiltonian system and solve the problem of construction asymptotically stability programm motion. The programm motion can be any function. Control is received in the method and the method of limiting functions and systems. In this case we use the Lyapunov's functions having constant signs derivatives. The following examples are considered: stabilization of program motions of homogeneous rod of variable length and stabilization of program motions of mathematical pendulum variable length in the rotation plane.

Keywords: Hamiltonian, nonautonomous system, program motion, Lyapunov's function.

УДК 62.534(031)

ПОСТРОЕНИЕ И СТАБИЛИЗАЦИЯ ПРОГРАММНЫХ ДВИЖЕНИЙ НЕАВТОНОМНЫХ ГАМИЛЬТОНОВЫХ СИСТЕМ

ВВЕДЕНИЕ

Задачи об управляемых программных движениях механических систем являются актуальными и привлекают внимание многих исследователей. Построение и исследование свойств и условий устойчивости таких движений рассматривались в работах многих ученых, например [1-4].

В данной работе ставится и решается задача об определении управления, реализующего и стабилизирующего произвольные заданные движения механической гамильтоновой системы. Решение задачи сводится к исследованию нулевого решения неавтономной системы и проводится на основе прямого метода Ляпунова [5]. Метод предельных систем [6] и его модификация [7] позволяют при использовании функций Ляпунова со знакопостоянными производными строить искомое управление в замкнутой аналитической форме в классе непрерывных функций.

1. ПОСТАНОВКА ЗАДАЧИ

Рассматривается управляемая механическая система, движение которой описывается уравнениями Гамильтона:

¿Ь др /1)

ёр дН (Ь,а,р) (

— + V,

■. ёЬ дд

где д = (д1 )т — п-вектор обобщенных координат в действительном линейном пространстве

Яп с нормой ||д||, р = (р1 ,...,рп)т — п-вектор обобщенных моментов, р е Яп с нормой ||р|| и Н(Ь,д,р) — гамильтониан системы, и, V — силы управляющих воздействий.

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

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

Введем следующее определение.

Определение. Программным (желаемым) движением системы назовем пару р* = р* (Ь), д* = = д*(Ь), где д*(Ь) и р*(Ь) — ограниченные, дважды кусочно-непрерывно дифференцируемые п-мерные вектор-функции, описывающие некоторое заданное движение механической системы.

В общем случае функции д(Ь), р*(Ь) могут не являться решением системы (1). Поэтому сначала поставим и решим задачу о реализации программного движения д* (Ь), р*(Ь) системы (1).

2. ПРОГРАММНОЕ УПРАВЛЕНИЕ И УРАВНЕНИЯ В ОТКЛОНЕНИЯХ

Пусть необходимо, чтобы система совершала некоторое программное движение (д*(Ь),р*(Ь)). Добавим к правой части системы (1) управляющие силы вида ир и ур:

д = р =

дн{г,д,р)

дН(Ь, д,р)

др

дд

+ ир, + ур.

(2)

Подставляя программное движение в систему (1) с управляющими силами ири ур, получим уравнения управляемого движения вида

р* =

др дд*

(3)

+ ур.

Из системы (3) вычислим программное управление:

ир = (* -

дН(Ь, д*,р*) др%

, дН&д*,р*)

(4)

ур = р* +

дд*

Подставим (4) в систему (2):

■ дН&д,р) дН$,д*,р*)

У ~-я--г 9--я-'

др дрг

. = дН{1,д,р) +г + дН{1, д\р*)

(5)

дд

дд*

Полученные уравнения (5) имеют решение, соответствующее выбранному программному движению (*(Ь),р*(Ь). Программные движения могут быть неустойчивыми, т.е. нереализуемыми на практике. Исследуем его на устойчивость и стабилизируем в случае необходимости до асимптотической устойчивости.

Сведем решение задачи о стабилизации программных движений к задаче стабилизации нулевого решения неавтономной гамильтоновой системы. Это позволит применить к задаче о стабилизации программных движений методы и результаты, разработанные для исследования устойчивости и стабилизации нулевого положения равновесия неавтономных систем [7,8]. Введем отклонения по правилу

{х = д - д*(Ь),

У = р - р*(Ь).

В силу линейности замены и линейности оператора дифференцирования структура уравнений Гамильтона при переходе к уравнениям в отклонениях не изменится.

Выпишем уравнения движения механической системы в отклонениях в виде системы уравнений Гамильтона. Уравнения (5) примут вид

x =

У =

dH(t, x + q* (t),y + p* (t)) dH(t, q* (t),p* (t))

dy dy

dH(t, x + q* (t),y + p* (t)) dH(t, q* (t),p* (t))

(6)

dx

+

dx

В дальнейшем будем рассматривать классические натуральные системы, для которых гамильтониан имеет вид [9]

1

(7)

H(t,p, q) = ^pTAp + U(t, q),

где А — матрица коэффициентов квадратичной формы кинетической энергии. С учетом гамильтониана (7) уравнения возмущенного движения в отклонениях (6) запишутся в виде

x = (A - A*)p* + Ay,

1 T dA {У = ~2У di"

rp dA 1 >-p P djy+2P

dA*

dq*

dA dq

* dU*

dq

d_U_

dq

(8)

Обозначим

*TdA r>u \ 1 *T

P gj=B{t,q), -p

dA* dA

dq* dq

* i p +

dU* dU

dq* dq

= Co(t, q).

Тогда систему (8) представим в виде

x = (A - A*)p* + Ay,

1 T dA {У = -2У а-У-Ву + Со.

(9)

3. ПОСТРОЕНИЕ СТАБИЛИЗИРУЮЩЕГО УПРАВЛЕНИЯ

Пусть C — положительно определенная, неисчезающая, ограниченная матрица:

соE < C < c\E (0 < с0 < c\ - const),

где E — единичная матрица.

Рассмотрим положительно определенную по отклонениям x и у, допускающую бесконечно малый высший предел, функцию Ляпунова [10]:

У(х, у, Ь) = -утАу + хтСх. Полная производная функции Ляпунова (10) по времени будет иметь вид

V - (—У • + (—V х+ —

\ду) У \dx) Х dt1

(10)

(11)

или

V = -y'

2 dt

У + \хТ^х + УТАСо + хтСАу - \уТА

dA

УтгУ

dq

+

+ (1ут^у)ТАу + (1ут^у)Т[{А ~ А*)р*] + {Сх)Т [(А " А*)р*] ■

Поставим и решим задачу о стабилизации программного движения д* (Ь), р*(Ь).

Определение. Управляющие воздействия ия = ия(х,у,Ь) и Vя = Vя(х,у,Ь) называются стабилизирующими, если они обеспечивают асимптотическую устойчивость невозмущенного движения х = у = 0 системы (9).

Постановка задачи о стабилизации: определить стабилизирующие управления ия = ия(х,у,Ь) и Vя = Vя(х,у,Ь), при добавлении которых в систему (9) ее решение х = у = 0 становится асимптотически устойчивым.

Выберем стабилизирующее управление в виде

и3 = -(А - А*)р*, у3 = -Со - Сх - А-1Бу,

(12)

где Б = Б(Ь) — симметричная, положительно определенная матрица размерности п х п. Добавив силы (12) в систему (9), имеем уравнения управляемых движений:

X = Ау,

1 дА

У = ~2уТ~д^У ~ ВУ~Сх~ А~гВу. Тогда полная производная (11) от функции Ляпунова в силу системы (13) имеет вид

(13)

У = -у1

1 дА

ЛВ~2т

У + \хТ^х + УТЛСо + хтСАу - Х-утА

дА

+

т

+ Ау + утАи = -ут

АВ-1-^ 2 дЬ

1 т дС т у + -хт—х + утАС0+

+хтСАу - т^УТА

дА У1ГУ

дд

1 дА 2У~д.х''

+ ( ^у^у ) Ау - утАС0 - утАСу - утВу =

= -У

1 дА

лв-ш

1 т дС 1 т У + ^х т—х — —у А у 2 дЬ 2

дА

1 дА т

+ ( 2У~д^У ' У~У У'

Перепишем:

где обозначено

у = -у1

1 дА

лв-2 т+Е>

1 т дС

N = —у А 2

дА дд

1 дА

2У~дхУ АУ

слагаемые старше второго порядка малости, ими можно пренебречь. В итоге имеем

V = -у1

АВ-1д4+В 2 дЬ

1 1 дС У+2Х ~тх-

(14)

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

Утверждение 1. Пусть матрицы А и С не зависят от Ь, т. е. выполнены условия

Тогда управление (12) стабилизирует программное движение д* = д* (Ь), р* = р* (Ь) при условии, что Б + АВ — положительно определенная. Утверждение 2. Пусть выполнены условия

дА дС

Тогда управление (12) стабилизирует программное движение д* = д* (Ь), р* = р* (Ь) при условии, что Б + АВ — положительно определенная. Утверждение 3. Пусть выполнены условия

дА дС

т*»' Ж = 0-

Тогда управление (12) стабилизирует программное движение д* = д* (Ь), р* = р* (Ь) при условии,

что

В>1-д4-АВ.

~ 2 дЬ

Утверждение 4. Пусть выполнены условия

дА дС

(15)

Тогда управление (13) стабилизирует программное движение д* = д* (Ь), р* = р*(Ь) при условии (15).

Представленные утверждения отличаются тем, что в утверждениях 2 и 4 производная функции Ляпунова есть отрицательно определенная по всем переменным, и наличие асимптотической устойчивости нулевого решения справедливо на основании теоремы Ляпунова об асимптотической устойчивости [5]. В утверждениях 1 и 3 производная (14) при управлении (12) является отрицательно определенной функцией только по переменным у , т. е. знакопостоянной. На основе теоремы 3 из [7] имеем асимптотическую устойчивость нулевого решения системы (9), т.е. программного решения д*(Ь),р*(Ь) исходной управляемой системы (1).

O

4. СТАБИЛИЗАЦИЯ ПРОГРАММНЫХ ДВИЖЕНИИ ОДНОРОДНОГО СТЕРЖНЯ ПЕРЕМЕННОЙ ДЛИНЫ

Рассмотрим однородный тяжелый стержень переменной длины массой m = 1, движущийся без трения в плоскости Oxy, которая вращается с постоянной угловой скоростью и вокруг неподвижной вертикальной оси Oy, лежащей в этой плоскости (рис. 1). Пусть £, п — координаты центра масс стержня, в — угол отклонения стержня от вертикали. Длина стержня изменяется по закону k = k(t) = а + bcost. Система имеет три степени свободы.

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

CJ2

У"

ю

Рис. 1

Гамильтониан системы примет вид

1

#(£,т/,0) = ^(С2 + V2 + к2 (t)e2 +LU2 к2 (t) cos2 в +UJ2 С2)- дг).

Выберем программное движение так, чтобы центр масс стержня двигался по окружности радиуса R =1 с центром в точке O(2, 2), и при этом стержень вращался в плоскости Oxy вокруг центра масс с постоянной угловой скоростью = const. Программное движение имеет вид

= cos t + 2, П* = sin t + 2, в* = u0t.

Уравнения в отклонениях (9) запишутся в виде

x i = yi,

X 2 = У2,

X 3 = Уз/к2, 2

У1 = U2X1, У2 =0,

y3 = —и2к2 sin x3 cos(x3 + 2u0t).

Выбрав управление в виде (12) согласно утверждению 3

ui = —и2х1 — c11x1 — d11 y1, U2 = —C22 X2 — d22 У2,

U3 = и2к2 sinX3 cos(x3 + 2q*) — C33X3 — к2dззyз,

получим уравнения движения стабилизированной системы:

x i = yi, x 2 = y2,

= yz/k2,

yi = -ciixi - diiyi, y2 = -022x2 - d22y2, Уз = -033x3 - k2d33y3.

5. СТАБИЛИЗАЦИЯ ПРОГРАММНЫХ ДВИЖЕНИЙ МАТЕМАТИЧЕСКОГО МАЯТНИКА ПЕРЕМЕННОЙ ДЛИНЫ

O

Рассмотрим маятник переменной длины на вращающемся основании. Пусть плоскость Oxy с подвижной системой координат Oxy вращается с постоянной угловой скоростью ш вокруг оси Oy. В точке B на подвижной оси Ox прикреплен математический маятник массой m переменной длины l = l(t), совершающий относительные движения в плоскости Oxy (рис. 2).

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

Гамильтониан системы примет вид

H(l,6) = ml2 + ml262 - у[/2 + 62l2 +w2(Co +2Co/sin6' + /2sin2 6)]-mgl cos 6.

B

Рис. 2

Выбрав программное движение в виде

l* = l0 + a sin t, в* = wot,

получаем следующие уравнения в отклонениях:

x i = yi/m,

y2 + mw0 (l0 + a sin t)2 %2 = —7-—;—i—, . ло--^o,

yi =

m(xi + l0 + a sin t)2

(y2 + mmojlo + asint)2)2 uj2 . . 2

-7-:-:—-rz--1--{{xi + lo + asmt) sm (X2+uJot) +

m(x2 + l0 + a sin t)3 m

+ £0 sin(x2 + w0t)) + mgcos(x2 + w01) - mw2(l0 + a sin t)-w2

--(£0 sin ujot + (lo + a sin t) sin2 u>ot) — mg eos u>ot,

w2£0(xi +10 + a sin t) . , w2 . . ,2 . ,

y2 =-cos(a:2 + uiot) + -—(xi + lo + a sm t) sm2(x2 + (¿ot) —

2m

w2

m

— mg(x 1 +lo + a sin t) sin(x2 + ooot) — -— (2Co (lo + a sin t) eos uJot+ + (l0 + a sint)2 sin2w0t) + mg(l0 + a sint) sinw0t.

На основе утверждения 1 управления

ui = 0,

U2 =

p*

p2

mq2 mq*2' ~pf , Pf

Vi = -7ч +

mq *3 m(q* + x i )3

w2 w2 H--(Co sin+ q{ sin2 - mgcosq?--(Co sin(g| + ж2) +

m

m

+ (q* + xi) sin (q* + x2)) - mgcos(q* + x2) - oi ^i - mdi^i, w2 w2

v2 = — (2Coq{ eos ql + qf sin 2^) + rngq\ sinq\ - — (2Co(q{ +^i)cos(g| + x2)+ + (q* + xi)2 sin2(q* + x2)) + mg(q* + xi) sin(q* + x2) - 022x2 - mqfd22y2.

Изв. Capar, ун-та. Нов. сер. Z011. Т. 11. Сер. Математика. Механика. Информатика, вып. 3, ч. Z стабилизируют выбранное программное движение. При этом управляемая система имеет вид

X i = yi/m

X 2 =

У2

yi =

m(xi + l0 + a sin t)2 У2 + 2y2mcvo (lo + a sin t)2

- ciixi - mdiiyi,

m(x2 + l0 + a sin t)3 УУ2 = -C22X2 - m(xi + lo + a sin t)2d22У2.

Проиллюстрируем полученные результаты графическим представлением численного интегрирования построенной управляемой системы. Интегрирование проведено при следующих значениях параметров системы: m = 1 кг, ш = 1 рад/с, l0 = 2 м, a = 0.5 м, ш0 = 5 рад/с, c11 = 10, c22 = 10, d11 = 10, d22 = 10 и следующих начальных отклонениях: x1 (0) = 0.1 м, x2(0) = 0.1 рад, y1(0) = 0.1 кгм/с, У2 (0) = 0.1 кг рад/с.

На рис. 3, а изображены отклонения x1, x2, а на рис. 3, б — отклонения y1, y2. Представленные графики демонстрируют асимптотическую сходимость отклонений, тем самым имеем стабилизированные произвольно выбранные программные движения маятника.

Результаты, полученные в работе, развивают и обобщают результаты из [8,10].

а б

Рис. 3

Библиографический список

1. Летов А. М. Динамика полета и управление. М.: Наука, 1969. 359 с.

2. Галиуллин А. С., Мухаметзянов И. А., Мухарля-мов Р. Г., Фурасов В. Д. Построение систем программного движения. М.: Наука, 1971. 352 с.

3. Зубов В. И. Проблема устойчивости процессов управления. Л.: Судостроение, 1980. 375 с.

4. Афанасьев В.Н., Колмановский В. Б., Носов В. Р. Математическая теория конструирования систем управления. М.: Высш. шк., 1989. 447 с.

5. Малкин И. Г. Теория устойчивости движения. М.: Наука, 1966. 530 с.

6. Artstein Z. Topological dynamics of an ordinary

equations // J. Differ. Equat. 1977. Vol. 23. P. 216-223.

7. Андреев А. С. Об устойчивости и неустойчивости нулевого решения неавтономной системы // ПММ. 1984. Вып. 2. С. 40-45.

8. Bezglasnyi S. P. The stabilization of equilibrium state of nonlinear hamiltonian systems // Seminarberichte aus dem Fachberrich Mathematik. FernUnivetsitat in Hagen. 2001. Bd. 71. P. 45-53.

9. Маркеев А. П. Теоретическая механика. М.: Физма-лит, 1990. 414 с.

10. Bezglasnyi S. P. The stabilization of program motions of controlled nonlinear mechanical systems // Korean J. Comput. and Appl. Math. 2004. Vol. 14, № 1. P. 251-266.