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

12. Демидович Б.П., Марон И.А., Шувалова Э.З. Численные методы анализа. Приближение функций, дифференциальные и интегральные уравнения. - СПб: Лань, 2010. - 400 с.

13. Мельников В.Г. Применение метода экономизации К. Ланцоша при исследовании нелинейных колебаний механических систем // Научно-технический вестник СПбГУ ИТМО. - 2004. - Т. 50. - № 15. -С. 16-18.

14. Матросов В.М., Румянцев В.В., Карапетян А.В. Нелинейная механика. - М.: Физматлит, 2001. - 432 с.

Мельников Виталий Геннадьевич - Санкт-Петербургский национальный исследовательский университет

информационных технологий, механики и оптики, кандидат технических наук, доцент, зав. кафедрой, melnikov@mail.ifmo.ru

УДК 531

АЛГОРИТМИЧЕСКАЯ РЕАЛИЗАЦИЯ МЕТОДА ИССЛЕДОВАНИЯ НЕЛИНЕЙНЫХ ДИНАМИЧЕСКИХ СИСТЕМ С.Е. Иванов

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

Введение

В работе предлагается алгоритм метода многочленных преобразований с последующими экспоненциальными преобразованиями и алгоритмические формулы, на основе которых автором разработан пакет программ для исследования нелинейных динамических систем с тремя степенями свободы [1]. Данный алгоритм метода реализован посредством универсальной системы МаШетайса с использованием математических библиотек. Рассматриваются нелинейные динамические системы, которые описываются системой дифференциальных уравнений шестого порядка с нелинейной частью в виде полинома четвертой степени относительно фазовых координат с постоянными и периодическими коэффициентами [2]. В результате применения модифицированного метода многочленных преобразований нелинейная математическая модель динамической системы приводится к автономному виду и находятся существенные константы, определяющие свойства динамической системы [3]. Метод применен к исследованию математической модели приборной системы с малыми нелинейными частями, которые существенно влияют на решение и полностью учитываются в предлагаемом методе [4]. Точность получаемых решений с помощью аналитического метода преобразований контролируется посредством сравнения с решением, получаемым численным методом Рунге-Кутта.

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

Нелинейная динамическая система

Рассматриваются динамические системы, которые описываются обыкновенными дифференциальными уравнениями с нелинейными членами в форме многочленов от фазовых переменных с постоянными и периодическими параметрами [6]. Для исследования таких систем применяется метод многочленных преобразований. Применение этого метода связано с большим объемом трудоемких выкладок, что приводит к необходимости создания программ, позволяющих выполнять расчеты средствами современной вычислительной техники. Автором составлен комплекс программ для реализации метода на языке программирования высокого уровня МаШетайса. Разработанные программы позволяют проводить исследование динамических нелинейных систем, получать аналитический вид решения и преобразовывать систему к автономному виду.

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

С.Е. Иванов

1д + Бд + Сд = I hv со8(ю/Г )v2 д/3 д2^ д3^ д/6 д2д3и

(1)

И=1

где д = [, д2, д3 ] - вектор обобщенных координат; I, Б, С - матрицы третьего порядка; V =(у1у2V3V4V5V6V7у8) - векторный индекс, |у| = У1 + V2 +... + у8 .

Для записи на языке программирования выражение суммы коэффициентов с векторным индексом представлено в следующем виде:

V7 ^8

4 >8 ¡7 >6 >5 '4 >3 '2

= ^ I I I I I I I \>2->1,

>8=1 >7 = о >6=о >5 = о >4=о >3=0 >2=о >1=о

'3 '2 >'4 3 > 5

Предполагается, что система содержит малые нелинейный части угБ^ < е .

Алгоритм метода

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

до = ехр(/ю/) и до = ехр(-/ю/), Х1 = >ю. Запишем в новых переменных исходную систему:

1 _ 1 _

соз(ю/) = - (до + до) и ) = — (до - до).

2 2/

(2)

Систему дифференциальных уравнений (1) можно представить в виде системы восьми дифференциальных уравнений первого порядка в нормальной форме Коши с фазовым вектором

Х = [ до дд дз дА дз ]

(3)

" Ш о о " " о "

X = о о I X + о

А-1Я -А-1С - Л-1 Б А-1в

где постоянная квадратная блочная матрица порядка имеет вид Ш =

. Выполним линейное

>ю о _ о ->ю

преобразование У = БХ , в результате которого линейная часть системы (3) приводится к диагональному виду:

У = ЛУ + я\ , , где Л = Шая^Д^.-АА].

(4)

Применим многочленное преобразование к системе (4), содержащее многочлены четвертой степени:

Уе = ^ +1 а^V, (5 = 3,...,8), 7V = г;1 ¿V2...¿V,

М=2

где а5 - неизвестные коэффициенты преобразования.

В результате многочленного преобразования система (4) приводится к автономному виду:

¿е = + V, (е = 3,...,8), (5)

М=2

где д* - искомые коэффициенты преобразованной системы.

Автором получены алгоритмические формулы для расчета коэффициентов преобразования и преобразованной системы, которые применены в составленном пакете программ:

I V +1 (а* гV (IV, - X*)) +1 (а? гVII Vkzk-1 II -) = (г), (5 = 3,...,8).

2 ^1= 2 к=1 2 к =3 |-|=2

Для нахождения особых значений векторного индекса при фиксированном 5 необходимо найти решение двух уравнений [7]:

IКVк -* о, IVк = 2,3,4.

(6)

Постоянные д* приравниваем к нулю при неособых значениях индексов; при таких значениях вычисляют постоянные а*. И наоборот, при особых значениях индексов полагают коэффициенты а*

равными нулю и вычисляют д^, .

В нерезонансном случае из уравнений (6) находим следующие особые индексы:

V=1

при : V = (00100011), V = (00101100), V = (00210000), V = (11100000), при д4: V = (00010011), V = (00011100), V = (00120000), V = (11010000), при : V = (00001011), V = (00002100), V = (00111000), V = (11001000), при д^: V = (00000111), V = (00001200), V = (00110100), V = (11000100), при д7: V = (00000021), V = (00001110), V = (00110010), V = (11000010), при д8 : V = (00000012), V = (00001101), V = (00110001), V = (11000001).

Для получения аналитического решения автономной системы (5) выполним замену переменных: zs = р*ехр(1т+ 0*)), = р*ехр(/(г 1тХх+1 -0*)), * = 3,5,7; ^ = ехр(+/гю). В результате получаем дифференциальное уравнение

рх = )рх + Яе(V), = Рх'1т(V), £ = 3,5,7

4 (7)

= Е<рГ4РГ6рГ8ехр(/(03^ -V4)+05(V5 -V6)+07(V7 -V8)-е8)).

М =2

В нерезонансном случае [8] экспонента не входит в систему (7), так как ее степень равна нулю. Приравнивая к нулю правые части системы (7), находим стационарные решения. Для получения решения системы (1) в исходных переменных выразим вектор У по формулам замены, обратной линейной (2):

X = Б-1У .

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

Заключение

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

Приведенный алгоритм реализован в созданных программах на языке программирования высокого уровня МаШетайса. Для приборных динамических систем с тремя степенями свободы построены графики, которые показали точность получаемых результатов и применимость метода к широкому кругу задач общего вида (1).

Работа выполнена при поддержке гранта РФФИ 10-08-01046-а.

Литература

1. Бабаков И.М. Теория колебаний. - М.: Дрофа. - 2004. - 591 с.

2. Зиновьев Н.М., Мяснянкин Ю.М. Введение в теорию колебаний конструкций: Учебное пособие. -Воронеж: Изд-во ВГУ, 2005. - 35 с.

3. Иванов С.Е. Определение установившихся режимов работы виброзащитной системы с двумя степенями свободы // Научно-технический вестник СПбГУ ИТМО. - 2010. - № 4 (68). - С. 44-46.

4. Иванов С.Е. Исследование нелинейных динамических систем с тремя степенями свободы // Научно-технический вестник СПбГУ ИТМО. - 2011. - № 4 (74). - С. 62-64.

5. Иванов С.Е. О реализации численно-аналитического метода многочленных преобразований на компьютере // Современные технологии: Труды молодых ученых ИТМО. - СПб: СПб ГИТМО (ТУ), 2001. - С. 138-141.

6. Кузнецов А.П. Нелинейные колебания: Учебное пособие. - М.: Издательство физико-математической литературы, 2002. - 292 с.

7. Мельников Г.И. Динамика нелинейных механических и электромеханических систем. - Л.: Машиностроение, 1975. - 198 с.

8. Мельников В.Г., Мельников Г.И., Иванов С.Е. Компьютерные технологии в механике приборных систем: Учебное пособие. - СПб: СПбГУ ИТМО, 2006. - 127 с.

Иванов Сергей Евгеньевич - Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кандидат физ.-мат. наук, доцент,