Вестник РУДН Серия Математика. Информатика. Физика. № 4. 2008. С. 37-41
УДК 531.36, 590.86
Исследование устойчивости интегрального многообразия в случае уравнений возмущения связей с постоянными коэффициентами
Е. А. Горшков
Кафедра теоретической механики Российский университет дружбы народов ул. Миклухо-Маклая, 6, Москва, Россия, 117198
Предлагается метод построения математической модели динамики управляемой механической системы. Строится алгоритм модификации уравнений динамики, которая позволяет решить задачу стабилизации связей и обеспечить требуемую точность численного решения соответствующей системы дифференциально-алгебраических уравнений, описывающих наложенные на систему связи, её кинематику и динамику. Предлагаемый метод может быть использован для исследования динамики систем различной физической природы.
Ключевые слова: устойчивость системы, возмущение связей, уравнения возмущения, интегральное многообразие, динамика системы, стабилизация связей, диссипативная функция, асимптотическая устойчивость.
1. Введение
Современные системы управления имеют достаточно сложную структуру, содержащую элементы различной физической природы. Исследование динамики таких систем в первую очередь требует построения общей математической модели, описывающей динамику всей системы. Кинематические и динамические аналогии позволяют использовать для моделирования кинематики и динамики управляемых систем методы классической механики [1]. Основываясь на известных принципах, уравнения динамики управляемых систем можно представить в форме уравнений Лагранжа [2,3].
Динамика сосредоточенных систем описывается системой дифференциально-алгебраических уравнений, составленных непосредственно из уравнений динамики и уравнений связей, наложенных на обобщённые координаты и скорости [4]. Поскольку число уравнений связей, как правило, оказывается меньше числа переменных, определяющих состояние системы, то построение уравнений динамики сводится к решению систем линейных алгебраических уравнений с прямоугольной матрицей коэффициентов и последующему построению дифференциальных уравнений, решения которых обладали бы заданными свойствами. Неоднозначность решения этой задачи позволяет учитывать возможные дополнительные условия, определяющие свойства исследуемой системы.
2. Моделирование динамики системы
Уравнения динамики механической системы могут быть построены, если известны её кинетическая энергия Т = Т* {дг,ср, Ь), потенциальная энергия Р = Р*(дг,Ь), диссипативная функция И = В*(дг,д>,£), обобщённые непотенциальные внешние силы Qs = Qs (дг ,д°, и управляющие силы из = из {дг ,д°, , где 4 = ж, 1,3,8 = !,...,п.
Статья поступила в редакцию 9 июня 2008 г.
Работа выполнена при финансовой поддержке РФФИ грант №06-01-00664.
Управляющие силы из, действующие на систему, призваны обеспечить выполнение уравнений связей
Г (дг, *) = 0, г = г:дг + Я = 0, Г (дг, г) = о, (1)
ги д/и : 9/: . .
Л = Л = = 1,...,щ ^ = 1,...,m, р = т + 1,...,г.
В (1) и в дальнейшем по одинаковым индексам предполагается суммирование.
Рассмотрим механическую систему, кинематическое состояние которой определяется координатами дг, дп+и и скоростями д^, дп+К, где ( к = 1,..., г). В данном случае дп+и, дп+К — избыточные координаты, которые вводятся с целью обеспечения устойчивости связей. Определим динамические показатели системы: кинетическую энергию, потенциальную энергию и диссипативную функцию соответственно как функции переменных: дг, ^, дп+и, дп+К, ¿:
Т = Т (д\ дп+и, дп+к, ^ , Р = Р(д\ дп+и, В = Б(д\ дп+и, д>, дп+к, г).
Представим связи, наложенные на систему, уравнениями:
дп+и - (д\ *) = о, дп+К - /к (д\ д>,г) = о. (2)
Предполагается, что функции Т, Р, В по крайней мере дважды дифференцируемы по всем переменным, и при дп+и = 0, дп+К = 0 расширенная система совпадает с исходной заданной системой, и выполняются уравнения (1).
Будем предполагать, что отклонения от уравнений связей (1), которые измеряются переменными дп+и, дп+К достаточно малы, то есть: существует положительное число ео = 0, такое, что для любого положительного е ^ £о выполняется |дп+и| ^ £, |дК| ^ £. Будем предполагать, что функции Т, Р, В можно представить разложением в ряд по степеням дп+и, дп+К и соответствующие ряды не содержат членов первой степени относительно дп+и, дп+К:
т = т* (д,д\ ^ + (д\ д\ ^ дп+Кдп+ +..., (3)
Р = Р* (д\ I) + 1Ь^ (д\ г) дп+идп+" + ..., (4)
в = Б*(д\д>,г) +1 (д\?,ь) дп+Кдп+ +..., (5)
к=1,..., г; = 1,...,т.
Из (2) следует система линейных алгебраических уравнений относительно 5дг,
§дп+К
¡гК8дг = 8 дп+К, (6)
где
ги = Ж р = Э^
Ь ддí, и Э .
Согласно принципу Даламбера-Лагранжа с учётом решения системы уравнений (6) и условиями идеальности связей (2) можно построить уравнения динамики расширенной системы. Они состоят из двух частей. Первая группа уравнений соответствует уравнениям динамики управляемой системы
<ид ^ дф = а? д с?+Уг +кХ- (7)
(8)
вторая группа представляет собой систему уравнений возмущений связей
ё дТ дТ _ дР дБ
(Ид дп+- дд п+>1 = дд п+>1 д дп+- '
Последнее слагаемое правой части уравнения (7) /-Хк представляет собой управляющие силы и3.
3. Исследование устойчивости системы
Решение системы уравнений (7)—(8) сводится к определению множителей Хк из (7) и выражениями избыточных переменных дп+дп+- из (2), (8) через обобщённые координаты дг и скорости д3.
Для получения выражений для Хк представим уравнения (7)-(8) в виде, разрешённом относительно старших производных дг, дп+-:
= ?, = тi (дк, д8, А + Г (дк, д8, А Х- (дк, д8, дп+1\ дп+-, А , (9)
^пт = 4п+1Л, = К , у8, + Щ (V, ч8, *) чп+ ■ (10)
Продифференцируем второе уравнение (2) по времени :
(дп+К д_г^.г _ о/-=0 (11)
( дд1 (1 д^ ( а . ()
А 'п + к А
В последнем уравнении необходимо заменить — и -—т правыми частями уравнений (9)-(10):
к-дп+1 + Ъ-дп+ - д-<? - ^ (т + ГХо) - = 0. (12)
В результате приходим к системе линейных алгебраических уравнений относительно множителей ХК:
а-еХе = а-. (13)
а-е = ^Г не а- = к-¿п+у + ь-ап+1 _ ^И^ _ ЁИт _ а = дд*1 , а =кУ4 + ^ Зд- 4 дд1 дЪ '
Система (13) имеет решение:
Хе = ае-а-, (14)
где аеК — матрица, обратная а-е : аеК = , а а- определяется суммой
слагаемых а- = а-(0) + а-1 +... , распределённых по степеням переменных дп+^, дп+- : а-(0) = — ^щгдг — ^цгтг — . При множителях ХК, определённых таким образом, уравнения связей (1) исходной системы оказываются справедливыми для решений дг = дг ({) дифференциальных уравнений динамики управляемой системы при дп+1Л = 0, дп+- = 0:
ж = * • ж = т +'а'-а-■ (15)
Отклонения решений системы (15) от уравнений связей (1) описываются уравнениями возмущений связей
^п+и ыдп+к
4 = дп+и 4 = кК дп+- + ЬКг,п+^ (16)
т =д , Ы =кк'д + Ькд . (16)
Необходимым условием стабилизации связей (1) является асимптотическая устойчивость соответствующего интегрального многообразия системы (15). Пусть уравнения (16) являются линейными относительно переменных дп+и, дп+К. В силу равенств (2) устойчивость интегрального многообразия в новых переменных (с учётом избыточных) можно рассматривать как устойчивость по части переменных [5] уравнений (15)-(16). Если в пространстве переменных дг, д определить расстояние до интегрального многообразия (1) величиной
1И1, 2= (г1,..., гт+г) , ги = дп+и, гт+К = дп+К,
то об устойчивости многообразия (1) можно судить по свойствам устойчивости тривиального решения дп+и = 0, дп+К = 0 системы (16).
Если в выражениях (3)-(5) коэффициенты ар,, Ър,, ср, — постоянные, то система (16) состоит из линейных дифференциальных уравнений с постоянными коэффициентами. В этом случае судить об устойчивости её тривиального решения можно по корням характеристического уравнения (ц25" — рЬ" — к') = 0. Если все корни имеют отрицательные действительные части, то интегральное многообразие устойчиво асимптотически.
Для исследования устойчивости тривиального решения в общем случае можно воспользоваться методом функций Ляпунова [6]. Если функция V = V ^дв, дк, дп+и, дп+к, ^^ является положительно определённой по переменным дп+и, дп+к, а её производная, вычисленная в соответствии с уравнениями (15), (16),
• = ду дУт + д^рва аК ду -п+и дУ кКдп+у , ду ЬК п+у, +ду
у = дф д + ддэт + дV } авка + ддп+ид + ддп+'к,д + эдп+' Ьу>д + т
отрицательно определённая функция, и функции V, дп+и, дп+К допускают бесконечно малый высший предел [6], то интегральное многообразие (1) системы (15) устойчиво асимптотически. Учитывая тождество т^ + адКак(0) = 0, последнее равенство принимает вид:
у = ^У дг + у "У + дп+и + Г+У< + ^УгЬиг+ + д~У,
ддг ^дд^ ддп+ич ддп+К ддп+К дЬ'
О 1.
В качестве функции Ляпунова можно также использовать положительно определённую квадратичную форму относительно переменных дп+и, дп+К:
2 V = 8 и „дп+ид" + 2ри Кдп+и д п+К + 1КГ, дп+' д1 ■ (17)
Коэффициенты 8ии, рик, предполагаются непрерывными, дифференцируемыми, ограниченными функциями переменных , • , . Производная функции V (17) приводится к виду
V = ^дп+ид»+р'иКдп+и дп+К + и дп+К д%,
где 8'и1,, р'ик, I' определяются через ар,,, Ьр,,, ср,,, к', Ь' уравнений (3), (4), (5), (16). Тогда для исследования устойчивости тривиального решения можно использовать критерий Сильвестра, в соответствии с которым можно судить о знакоопределённости V и V.
4. Заключение
Учёт возможных отклонений от уравнений связей управляемых систем на этапе составления уравнений позволяет решить задачу стабилизации связей. Для этого уравнения возмущений связей должны быть построены в соответствии с требованием асимптотической устойчивости по отношению к уравнениям связей.
Приведённые модели решения реализованы автором в системе символьных вычислений Maple и могут быть использованы для решения задач управления программным движением.
Литература
1. Layton R. A. Principles of Analytical System Dynamics. — N.Y.: Springer, 1998. — 158 p.
2. Мухарлямов Р. Г. О построении систем дифференциальных уравнений движения механических систем // Дифф. уравнения. — № 3, вып. 39. — 2003. — С. 343-353.
3. Мухарлямов Р. Г. О численном решении дифференциально-алгебраических уравнений // Вестник РУДН. Сер. Прикладная математика и информатика. — № 1. — 1999. — С. 33-37.
4. Горшков Е. А. Построение уравнений динамики управляемой системы // Вестник РУДН, сер. Прикладная математика и информатика. — № 3-4. — 2007. — С. 20-24.
5. Румянцев В. В., Озиранер А. С. Устойчивость и стабилизация движения по отношению к части переменных. — М.: Наука, 1987. — 256 с.
6. Мухарлямов Р. Г. О построении множества систем дифференциальных уравнений устойчивого движения по интегральному многообразию // Дифференциальные уравнения. — Т. 5, № 4. — 1969. — С. 688-699.
UDC 531.36, 590.86
Research of Stability of Integral Variety for Perturbation of Constraints with Constant Coefficients
E. A. Gorshkov
Department of Theoretical Mechanics Poeples' Friendship University of Russia 6, Miklukho-Maklaya str., Moscow, Russia, 117198
The method of constructing of mathematical model of dynamics of operated mechanical system is offered. The algorithm of updating of the equations of dynamics which allows to solve a problem of stabilization of communications is under construction and to provide demanded accuracy of the numerical decision of corresponding system of the differential-algebraic equations describing imposed on communication system, its kinematics and dynamics. The offered method can be used for research of dynamics of systems of the various physical nature.