WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!


На правах рукописи

Коверга Александр Юрьевич

ПРИМЕНЕНИЕ ЛОКАЛЬНЫХ МЕТОДОВ В ИССЛЕДОВАНИИ КОЛЕБАТЕЛЬНЫХ РЕШЕНИЙ НЕКОТОРЫХ УРАВНЕНИЙ С ЗАПАЗДЫВАЮЩИМ АРГУМЕНТОМ

Специальность 01.01.02 – дифференциальные уравнения, динамические системы и оптимальное управление

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата физико-математических наук

Ярославль – 2012

Работа выполнена на кафедре математического моделирования федерального государственного бюджетного образовательного учреждения высшего профессионального образования Ярославский государственный университет имени П.Г. Демидова Научный доктор физико-математических наук, руководитель профессор Кубышкин Евгений Павлович Официальные оппоненты: Рудый Александр Степанович – доктор физико-математических наук, профессор, Ярославский государственный университет им. П.Г. Демидова, заведующий кафедрой микроэлектроники Старков Сергей Олегович – доктор физикоматематических наук, профессор, Обнинский институт атомной энергетики, заведующий кафедрой компьютерных систем, сетей и технологий Ведущая ФГБОУ ВПО Самарский государственный организация университет

Защита состоится мая 2012 г. в часов на заседании диссертационного совета Д 212.002.05 в Ярославском государственном университете им. П.Г. Демидова по адресу: 150000, г. Ярославль, ул. Советская, д. 14.

С диссертацией можно ознакомиться в библиотеке Ярославского государственного университета им. П.Г. Демидова по адресу: 150000, г. Ярославль, ул. Полушкина роща, д. 1.

Автореферат разослан 2012 г.

Ученый секретарь Глызин С.Д.

диссертационного совета

Общая характеристика работы



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

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

Принцип сведения использовал А.М. Ляпунов2 при изучении устойчивости решений в критических случаях, хотя понятие инвариантного многообразия он не использовал. Различные вопросы теории инвариантных многообразий и принципа сведения для систем обыкновенных дифференциальных уравнений рассматривали Д.В. Аносов3, В.А. Плисс4, S. Sternberg5, A. Kelley6, Ю.Н.

Бибиков, Дж. Хейл7.

Пуанкаре, А. О кривых, определяемых дифференциальными уравнениями / А. Пуанкаре М.-Л.:

Гостехиздат, 1947. 392 с.

Ляпунов, А.М. Общая задача об устойчивости движения / А.М. Ляпунов М.-Л.: Гостехиздат, 1950.

473 с.

Аносов, Д.В. Многомерный аналог одной теоремы Адамара / Д.В. Аносов // Науч. докл. высшей школы (физ.-мат. н.). 1959. №. 1. С. 3–12.

Плисс, В.А. Принцип сведения в теории устойчивости движения / В.А. Плисс // Изв. АН СССР. Сер.

матем. 1964. Т. 28, Вып. 6. С. 1297–1324.

Sternberg, S. Local constructions and a theorem of Poincare / S. Sternberg // Amer. J. Math. 1957.

V. 79 P. 175–187.

Kelley, A. The stable, center-stable, center-instable, unstable manifolds / A. Kelley // J. Diff. Equat.

1967. V. 3 P. 546–570.

Хейл, Дж. Колебания в нелинейных системах / Дж. Хейл. М.: Мир, 1969. 232 с.

Эти результаты систематизированы в монографиях Ф. Хартмана8, Ю.А.

Митропольского и О.Б. Лыковой9, а также А.М. Самойленко10 во введении которой имеется достаточно подробный обзор по указанной тематике.

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

Куликова, M. Hirch, C. Pugh11. Систематизированное изложение данных вопросов можно найти в монографиях Дж. Марсдена, М. Мак-Кракена12, Д.

Хенри13, Б. Хэссарда, Н. Казаринова, И. Вэна14. Там же можно найти многочисленные приложения указанной теории.

Метод построение системы обыкновенных дифференциальных уравнений, описывающих поведение траекторий на критическом инвариантном многообразии (нормальной формы) для уравнений с запаздывающим аргументом был впервые предложен Ю.С. Колесовым15. Построение ведется в амплитудной форме (полярных координатах). В работе Е.П. Кубышкина16 предложен более удобный способ построения нормальных форм уравнений с запаздывающим аргументом. Этот метод также использовался в работе С.Д. Глызина, Е.П. Кубышкина17. С различных позиций в квазилинейной постановке коХартман, Ф. Обыкновенные дифференциальные уравнения / Ф. Хартман М.: Мир, 1970. 720 с.

Митропольский, Ю. А. Интегральные многообразия в нелинейной механике / Ю.А. Митропольский, О.Б. Лыкова М.: Наука, 1973. 512 с.

Самойленко, А. М. Элементы математической теории многочастотных колебаний. Инвариантные торы / А.М. Самойленко М.: Наука, 1987. 301 с.

Hirch, M. Stable manifolds and hyperbolic sets / M. Hirch, C. Pugh // Proc. Symp. Pure Math., XIV, Am.

Math. Soc. 1970. V. 14 P. 133–163.

Марсден, Дж. Бифуркация рождения цикла и ее приложения / Дж. Марсден, М. Мак-Кракен М.:

Мир, 1980. 368 с.

Хенри, Д. Геометрическая теория полулинейных параболических уравнений / Д. Хенри М.: Мир, 1985. 376 с.

Хэссард, Б. Теория и приложения бифуркации рождения цикла: Пер. с англ. / Б. Хэссард, Н. Казаринов, И. Вэн М.: Мир, 1985. 279 с.

Колесов, Ю.С. Метод нормальных форм для систем с запаздыванием / Ю.С. Колесов // Литовский математический сборник. 1980. Т. 20, №. 4. С. 73–78.

Кубышкин, Е.П. Некоторые вопросы динамики распределенных роторов / Е.П. Кубышкин // Математика в Ярославском университете: Сборник обзорных статей к 25-летию математического факультета.

Ярославль. 2001. С. 157–182.

Глызин, С.Д. Нелинейная динамика одного дифференциального уравнения второго порядка с периодически возмущенным запаздыванием / С.Д. Глызин, Е.П. Кубышкин // Моделирование и анализ лебательные решения уравнений с запаздывающим аргументом изучались в работах А.Д. Мышкиса, С.Н. Шиманова18, В.П. Рубаника19, В.Н. Фодчука20.

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

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

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

Научная новизна работы В диссертации выявлены условия возникновения колебательных решений математической модели Ланга-Кобаяши полупроводникового лазера с запаздывающей обратной связью; выявлены условия возникновения хаотических колебаний в зоне комбинационного параметрического резонанса в математической модели генератора электромагнитных колебаний с запаздывающей обратной связью; исследованы условия возникновения хаотических колебаний в математической модели распределенного ротора из материала с нелинейно наследственными свойствами, одна из опор которого испытывает вибрацию; изучена в нелинейной постановке задача двухчастотного параметрического воздействия, в случае основного резонанса. В этом случае выявлены условия генерации хаотических колебаний.





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

информационных систем. 2005. Т.12. № 1. С. 40–45.

Шиманов, С.Н. Некоторые вопросы теории колебаний систем с запаздыванием / С.Н. Шиманов // Пятая летняя математическая школа. Киев: Ин-т матем. АН УССР. 1968. С. 473–549.

Рубаник, В.П. Колебания квазилинейных систем с запаздыванием / В.П. Рубаник – М.: Наука, 1969.

287 с.

Фодчук, В.Н. О непрерывной зависимости решения дифференциальных уравнений с запаздывающим аргументом от параметра / В.Н. Фодчук // Укр. мат. журн. 1964. Т. 16, №. 2. С. 273–279.

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

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

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

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

Результаты диссертации могут быть использованы при получении научнообоснованных рекомендаций при проектировании генераторов хаотических электромагнитных колебаний.

Апробация работы Основные результаты работы были представлены на Второй Всероссийской научной конференции Математическое моделирование и краевые задачи (Самара, 2005), Международной научной конференции Математические методы в технике и технологиях (ММТТ-20) (Ярославль, 2007), Воронежской зимней математической школе С.Г. Крейна (Воронеж, 2008), Воронежской зимней математической школе С.Г. Крейна (Воронеж, 2012), Международной научной конференции, посвященной 35-летию математического факультета и 25-летию факультета информатики и вычислительной техники Ярославского гос. университета им. П.Г. Демидова (Ярославль, 2012).

Кроме того, результаты диссертации неоднократно докладывались на семинаре кафедры математического моделирования Ярославского государственного университета им. П.Г. Демидова, а также обсуждались на семинаре научно-образовательного центра Нелинейная динамика Ярославского государственного университета им. П.Г. Демидова.

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

Структура и объем диссертации Диссертационная работа состоит из введения, четырех глав, заключения и списка литературы, содержащего 66 наименований. Работа содержит рисунков. Общий объем диссертации составляет 80 страниц.

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

В первой главе работы исследуется математическая модель ЛангаКобаяши полупроводникового лазера с запаздывающей обратной связью. Она представляет собой систему обыкновенных дифференциальных уравнений с запаздывающим аргументом следующего вида = k(1 + i)(N - 1)E + exp(-i0)E(t - ), (1) = -||(N - J + |E|2N). (2) Здесь E(t) = Ex(t) + iEy(t), i = -1 - комплексная переменная, описывающая электромагнитное поле, N(t) - плотность носителей зарядов, J ток накачки, k - коэффициент затухания колебаний, 1/|| - время спонтанной эмиссии, - коэффициент, характеризующий лазер, - процент отраженного излучения, 0 - фазовый сдвиг излучения, - величина запаздывания, равная времени, которое необходимо излучению, чтобы достичь зеркала и вернуться обратно.

Изучаются автоколебательные решения системы уравнений (1)-(2), бифурцирующие из состояния равновесия E(t) = 0, N(t) = J (3) при изменении параметров системы уравнений.

Характеристическое уравнение линеаризованной на (3) системы уравнений имеет вид - A(1 + i) - exp(- - i0) = 0, A = k(J - 1). (4) Расположение корней характеристического уравнения исследуются методом D-разбиений21. На рис. 1 представлена характерная картина D-разбиений при указанных значениях параметров.

Неймарк, Ю.И. D-разбиение пространства квазиполиномов (к устойчивости лианеаризованных распределенных систем) / Ю.И. Неймарк // ПММ 1949. Т. 13, №. 4. С. 349–380.

Рис. 1. = 4, = 0.005, 0 = /2.

Как следует из рис. 1, потеря устойчивости состояния равновесия (3) может происходить с прохождением одного (i1), либо двух (i1, i2, |1| < |2|) корней характеристического уравнения (3) через мнимую ось комплексной плоскости. При этом оказывается, что резонансного соотношения |2|/|1| = 1 реализовано быть не может.

Изучается характер колебательных решений системы уравнений (1)-(2), бифурцирующих из (3) в случае потери устойчивости, связанной с прохождением двух корней характеристического уравнения (4) через мнимую ось комплексной плоскости. В этом случае поведение решений системы уравнений (1)-(2) в окрестности состояния равновесия (3) определяется следующей системой обыкновенных дифференциальных уравнений 1 = (i1 + (1))z1 + (d11|z1|2 + d12|z2|2)z1 +... Z1(z1, z2, z1, z2; ), (5) 2 = (i2 + (1))z2 + (d21|z1|2 + d22|z2|2)z2 +... Z2(z1, z2, z1, z2; ), (6) коэффициенты которой эффективно вычисляются, - малый параметр.

В диссертации выполнен анализ поведения решений системы уравнений (5)-(6), установлена взаимосвязь между решениями уравнений (5)-(6) и системы уравнений (1)-(2), построены асимптотические (по ) формулы для периодических и инвариантных торов системы уравнений (1)-(2).

Во второй главе рассматривается задача параметрического возбуждения хаотических колебаний в дифференциальном уравнении второго порядка с запаздывающим аргументом следующего вида + A + x + [B + G(x, , x(t - h(t)), (t - h(t))] (t - h(t)) = 0, (7) где h(t) = h + a · sin(t); A, B, h, a, - положительные параметры (h > a), G(x1, x2, x3, x4) = g1x1 + g2x2 + g3x3 + g4x4 + g11x2 + g12x1x2 +...

достаточно гладкая нелинейная функция.

Уравнения вида (7) возникают при моделировании электронных устройств с активными нелинейными элементами и запаздывающей обратной связью.

Изучается возможность возбуждения за счет периодического изменения запаздывания сложных, в том числе хаотических колебаний. При этом предполагается, что при a = 0 нулевое решение уравнения (7) асимптотически устойчиво.

Положим a = 0 и рассмотрим характеристическое уравнение линейной части уравнения (7) P () 2 + A + 1 + B exp(-h) = 0. (8) Анализ расположения корней (8) проводится методом D-разбиений, из которого следует, что при определенных значениях параметров уравнение может иметь корни ±ij, (j > 0, j = 1, 2). При этом остальные корни имеют отрицательные вещественные части. При этом оказывается, что при A = A0 = 6/6, = B0 = 6/3, h = h0 = 4 2/3 уравнение (8) имеет корни B ±i1 = ±i 2/2, ±i2 = ±i 2, т.е. имеет место внутренний резонанс 1 : 2.

Указанный случай рассматривается во второй главе.

Положим A = A0 +A1, B = B0 +B1, h = h0 +h1, и a = a1, и выберем = 1 + 2 + , ( - расстройка резонанса).

Таким образом рассматривается случай комбинационного параметрического резонанса в нелинейной постановке при наличии внутреннего резонанса 1 : 2.

Уравнение (7) имеет в окрестности нуля фазового пространства C(-h(t, ), 0) C(-h(t, ), 0) четырехмерное 2/ периодическое локальное асимптотически устойчивое гладкое интегральное многообразие поведение решений на котором определяет поведение решений уравнения (7) из некоторого фиксированного шара с центром в нуле фазового пространства. Поведение решений на интегральном многообразии определяется поведением решений следующей системы обыкновенных дифференциальных уравнений 1 = (i1 + 1 + d11|z1|2 + d12|z2|2)z1 + b1z1z2 + c1z2 exp(i) +... Z1(t, z1, z2, z1, z2; ), ( = t), (9) 2 = (i2 + 1 + d21|z1|2 + d22|z2|2)z2 + b2z1 + c2z1 exp(i) +... Z2(t, z1, z2, z1, z2; ), (zj = zj(t) C, j = 1, 2). (10) В (9)-(10) точками обозначены слагаемые, имеющие по соответствующим переменным более высокий порядок малости. В явном виде приведены лишь главные“ слагаемые разложений.

” Комплексные коэффициенты системы эффективно вычисляются.

1 Рассмотрим линейную часть уравнений (9)-(10). Обозначим 1 = j +ij j и выберем A1, B1, h1 таким образом, чтобы j < 0 (j = 1, 2). Как следует из результатов работы22 в этом случае в плоскости параметров (, ) существует область параметрического (комбинационного) резонанса определяемая неравенством (1)(a1) + o() < 1 + 2 - < (2)(a1) + o(), (11) где 1 1 1 1 1 1 1 (1)(a1) = (2 - 1 ) Im(c1c2) + (1 + 2 ) B /(21 2 ) + 1 - 2, 1 1 1 1 1 1 1 (2)(a1) = (2 - 1 ) Im(c1c2) - (1 + 2 ) B /(21 2 ) + 1 - 2, а 1 1 1 B = Im2(c1c2) + 41 2 Re(c1c2) - 4(1 2 )2.

В условиях (11) исследуется поведение решений системы уравнений (9)(10).

1/2 -Положим в (9)-(10) zj = j exp(i), j 0, (j = 1, 2), t t счи1/тая при этом gj = g, (j = 1,..., 4), и выделим главную часть уравнений медленных“ переменных 1, 2, 1 = t - 1 - 2, 2 = 21 - 2.

” Сделаем это для конкретных значений параметров: g1 = 22, 0; g2 = 23, 0;

g3 = 55, 0; g4 = 6, 0; g33 = -1, 8. Остальные gjk = 0, (j, k = 1,..., 4). Нормируем j j(j /ajj)1/2, (ajj > 0). Выбрав теперь A1, B1, h1 таким образом, 1 чтобы 1 = 2 = 1, 1 = 2, получим систему уравнений 1 = (-1 - 2 + a122)1 + c1 cos(1 + 1)2 + b1 cos(-2 + 1)12, (12) 1 2 = (-1 + a212 - 2)2 + c2 cos(1 + 2)1 + b2 cos(2 + 2)2, (13) 1 2 1 = 1 - b112 + b122 - c1 sin(1 + 1)2/11 -c2 sin(1 + 2)1/2 - b1 sin(-2 + 1)2 - b2 sin(2 + 2)2/2, (14) 2 = 2 + b212 + b222 + 2c1 sin(1 + 1)2/11 -c2 sin(1 + 2)1/2 + 2b1 sin(-2 + 1)2 - b2 sin(2 + 2)2/2, (15) в которой a12 = -5.69; a21 = -0.705; b1 = 8.17; b2 = 4.39; c1 = a1 · 0.3;

c2 = a1 · 0.014; b11 = 0.138; b12 = 0.527; b21 = 0.454; b22 = 0.367; 1 = -0.222;

2 = 0.322.

Параметры 1 и 2 - характеризуют, соответственно, расстройку параметрического возбуждения и расстройку внутреннего резонанса. Это свободные параметры. Таким образом система зависит от трех параметров - 1, 2 и a1.

Кубышкин, Е.П. Параметрический резонанс в линейных системах с последействием / Е.П. Кубышкин // Исследования по устойчивости и теории колебаний: Межвузовский тематический сборник. Ярославль. 1976. С. 43–76.

Приведем некоторые результаты численного исследования системы (12)(15). Система численно анализировалась с использованием программы TracerВ рассматриваемом случае область параметрического резонанса в плоскости (1, a1) симметрична относительно оси a1. На рис. 2 приведена половина области неустойчивости (она заштрихована). Она определяется функциями (1)(a1) и (2)(a1), приведенными в (11). Положим b1 = b2 = 0, т.е исключим влияние внутреннего резонанса. Уравнения (12)-(13) в этом случае не зависят от 2. При малых a1 и любых 1 решения (12)-(13) стремятся к единственному состоянию равновесия 1 = 2 = 0, 1 = 0. При увеличении a1 и переходе границы области параметрического резонанса от указанного состояния равновесия ответвляется асимптотически устойчивое состояние равновесия вида 10 > 0, 20 > 0, 0. Дальнейшее увеличение a1 приводит к увеличению 10, 20. Отметим, что этому состоянию равновесия в уравнении (7) отвечает асимптотически устойчивый инвариантный тор.

Рис. 2. Область неустойчивости Рис. 3. Проекция аттрактора на (1, 2) Пусть теперь b1, b2 выбраны согласно (12)-(15). При a1, 1 принадлежащих области устойчивости, представленной на рис. 2, и произвольном 2 все решения (12)-(15) по 1 и 2 стремятся к нулю. При пересечении границы неустойчивости от нулевого решения ответвляется устойчивый цикл, размеры которого увеличиваются с ростом a1. При a1, 1 принадлежащих области параметрического резонанса возможно сложное поведение траекторий. Так, при a1 = 37.571, 1 = 5.0 существует хаотический аттрактор. Значения его ляпуновских показателей равны 1 = 0.02; 2 = 0; 3 = -0.01; 4 = -10.14, а ляпуновская размерность dL 3.00. На рис. 3 приведена проекция этого аттрактора на плоскость 1, 2.

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

Глызин, Д.С. Пакет программ для анализа динамических систем “Tracer”. Заявка №2008610548 от 14.02.2008г. Свидетельство о государственной регистрации программы для ЭВМ №2008611464. Зарегистрировано в Реестре программ для ЭВМ 24.03.2008г.

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

Предполагается, что одна из опор ротора испытывает периодическое воздействие (вибрацию). Материал ротора считается наследственно вязкоупругим и подчиненным следующей реологической модели вязкоупругого тела24:

f((t)) - R()f((t + ))d, (t) = E - где (t), (t) соответственно напряжение и относительная деформация, E модуль Юнга, R() - функция релаксации, f() = +f33+f55+..., (fj > 0) нелинейная функция деформации. Относительно функции R() (- < < 0) предполагается выполнение следующих условий:

dR() > 0, R() > 0, R()d < 1, d- R() M0 exp(0), (M0, 0 > 0), -. (16) Математической моделью рассматриваемой механической системы является следующая краевая задача2u u 2 u 2 2u 2 2u + a + b t2 t t s2 s2 s 2u(s, t + ) 2 2u(s, t + ) - R() exp(-i)b d = 0, (17) s2 s- u u|s=0 = = 0, (18) s s= b 2u 2 2u - R() exp(-i)b 2u(s, t + ) 2 2u(s, t + )d = s2 s2 s2 s- s== 1 exp(it + i1), (19) b 2u 2 2u - R() exp(-i)b 2u(s, t + ) 2 2u(s, t + )d = s s2 s2 s2 s- s=Работнов, Ю.Н. Элементы наследственной механики твердых тел / Ю.Н. Работнов М.: Наука, 1977. 384 с.

= 2 exp(it + i2), (20) где u(s, t) = ux(s, t) + iuy(s, t), i = -1, 0 < 1, 2 1, 0 1, 2 2, R.

Краевая задача приведена в безразмерных переменных s = z/l, u = u /l, t = t /t0, = t0, t0 = m1/2l2(EI0)-1/2, функции a(), b() являются аналитическими в окрестности точки = 0 и имеют вид a() = a0 + a1 +..., b() = 1 + b1 +....

Здесь u (z, t ) = u x(z, t ) + iu y(z, t ) смещения средней линии ротора в направлении осей OXи OY соответственно; ось OZ системы координат OXY Z, связанной с инерциальным пространством, направлена вдоль средней оси недеформированного ротора; t - время; m - погонная масса ротора;

2(j+1) aj = l4-2jm1-j(EI0)j-1a j, bj = (I0l2j)-1Ijf2j+1; Ij = x2(j+1)dx1 = y1 dyмоменты инерции поперечного сечения ротора относительно одной из осей соответствующих порядков; функция a () = a 0 + a 1 +... (a j > 0) характеризует внешнее нелинейное трение; j, j (j = 1, 2) и характеризуют амплитуду, фазу и частоту внешних изгибающего момента и периодической силы.

Изучается возможность и условия возникновения в краевой задаче (17)(20) хаотических колебаний (странных аттракторов).

Уравнение (17) является уравнением с бесконечным запаздыванием аргумента. Следуя16 дается определение фазового пространства для краевой задачи (17)-(20), пространства начальных условий и понятие решения.

Положим сначала 1 = 2 = 0 и исследуем устойчивость нулевого решения краевой задачи (17)-(20). Показано, что устойчивость нулевого решения (17)-(20) определяется характером расположения корней последовательности характеристических уравнений ln() 2 + a0 + n(1 - R() exp(( - i))d) = 0, (21) - n = 1, 2,..., n = n, где n - положительный корень уравнения ch n · cos n + 1 = 0. При этом потеря устойчивости уравнения (21) происходит по соответствующей функции en(s) = wn(s)/ wn(s), wn(s) = (sh n + Lsin n)(ch(ns) - cos(ns)) - (ch n + cos n)(sh(ns) - sin(ns)).

Расположение корней удобно исследовать методом D-разбиений. Положив = i и выделив вещественную и мнимую части, имеем nRS( - ) -2 + n(1 - RC( - )) = 0, a0 = -, где RC() = R(-) cos()d, RS() = R(-) sin()d, (22) 0 составляющие нормированного комплексного модуля упругости материала E() = (1 - RC() + iRS()), который определяется экспериментально.

Отметим, что согласно условию (16) при > 0, 0 < RC(), RS() < 1, RC(), RS() 0 при . Для R(), удовлетворяющих условиям (16), функции (22) качественно имеют вид, представленный на рис. 4. В соответствии с этим, кривые на плоскости (a0, ), соответствующие корням характеристического уравнения (21), расположенным на мнимой оси, качественно имеют вид, представленный на рис. 5. При этом каждая кривая является границей области устойчивости (неустойчивости) нулевого решения (17)-(20) по n-ой собственной функции en(s). Области неустойчивости заштрихованы, n0 = n(1 - RC(0))1/2. Как следует из рис. 5, потеря устойчивости решений может происходить по одному или по двум собственным функциям en(s). В последнем случае каждая функция en(s) имеет собственную частоту колебаний. Такой случай в дальнейшем и рассматривается.

Рис. 4. Вид функций RC и RS Рис. 5. Кривые, соответствующие корням Пусть теперь 1, 2 = 0. Точку пересечения кривых, исходящих из n и n+10 обозначим (an, n), а соответствующие им значения через n и n+1. Будем изучать поведение решений начально-краевой задачи (17)-(20) при изменении параметров в окрестности указанных точек. Введем для этого параметр 0 < 1 и положим a = an + an1, = n + n1, j = j01/2, (j = 1, 2), = 0 + = 2n - n+1 + 1. (23) Изучим характер установившихся колебательных решений краевой задачи (17)-(20) возникающих в окрестности нулевого решения при потере его устойчивости в предположениях (23). Краевая задача (17)-(20) имеет в окрестности нуля фазового пространства четырехмерное 2/-периодическое локальное устойчивое гладкое интегральное многообразие, поведение решений на котором определяет поведение решений начально-краевой задачи (17)(20). Система обыкновенных дифференциальных уравнений, описывающая поведение траекторий краевой задачи (17)-(20) на интегральном многообразии будет иметь вид 1 = (in + 1 + d11|z1|2 + d12|z2|2)z1 + 1/20A1z1z2 exp(it) +..., (24) n 2 = (in+1 + 1 + d21|z1|2 + d22|z2|2)z2 + 1/20A2z1 exp(-it) +..., (25) n+коэффициенты которой эффективно вычисляются через параметры краевой задачи (17)-(20).

В диссертации в качестве примера рассмотрено ядро вида ()- R() = exp(-) (0 < < 1) (1 - ) где (1 - ) - гамма функция Эйлера, (модифицированное ядро Абеля)24. В результате будем иметь RC() = cos((1 - ) arctg()), 1- (1 + 2) RS() = sin((1 - ) arctg()).

1- (1 + 2) Положим = 0.84 и обратимся к рис. 5. Соответствующие значения параметра равны 1 = 2.888 и 2 = 21.424.

Выберем n1 = 1.2, an1 = 1.2, a1 = 0.1, b1 = 0.21. Считая 10 = 20 = 0, вычислим коэффициенты системы (24)-(25). В (24)-(25) положим zj = 1/2jexp(ij), (j 0, j = 1, 2), t -1t и выделим главную часть уравнений медленных“ переменных 1, 2, = t - 21 + 2. Пронормировав теперь ” 1 0.051, 2 0.022 получим следующую систему дифференциальных уравнений:

1 = (3.94 - 2 - 1.252)1 + 1.8012 cos( - 0.46), 1 2 = (4.7 - 1.412 - 2)2 + 0.00302 cos(- + 0.55), 1 2 = + 2.282 + 1.932 - 3.602 sin( - 0.46) + 0.00301/2 sin(- + 0.55), 1 зависящую от двух параметров 0 и .

Система численно анализировалась при разных значениях параметров 0 и с использованием программы Tracer23. Система может иметь как устойчивые состояния равновесия, периодические решения, так и хаотические колебания. Так, при = 2.1, изменяя 0 имеем следующую динамику. При 0 = 2.8 имеем устойчивое состояние равновесия с координатами 10 = 1.624, 20 = 0.991, 0 = 8.329. Затем из состояния равновесия при 0 = 3.25 происходит рождение цикла и далее при 0 = 4.73, 0 = 4.961, 0 = 4.997 происходит серия бифуркаций удвоения периода. В результате чего при 0 = 5.образуется хаотический аттрактор. Его ляпуновские показатели равны 1 = 0.38; 2 = 0; 3 = -7.195, а ляпуновская размерность dL 2.053.

Рис. 6. Проекция аттрактора на (1, 2) Рис. 7. Сечение Пуанкаре = /Проекция аттрактора на плоскость (1, 2) приведена на рис. 6. При этом переменная неограниченно возрастает при t . На рис. 7 изображено сечение Пуанкаре плоскостью = /2. Просматривается фрактальная структура предельного множества.

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

Рассматривается дифференциальное уравнение с запаздывающим аргументом x(t) + x(t) + kf(x(t - h(t, ))) = 0, (26) где k > 0 некоторый параметр, f(x) = x + f2x2 + f3x3 +..., |x| < x0 гладкая функция, h(t, ) = h(1 + a1 cos(1t + 1) + a2 cos(2t + 2)) - величина запаздывания аргумента, в которой h, j, aj > 0; 0 j < 2, j = 1, 2; 0 < 1.

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

Положим в (26) aj = 0 (j = 1, 2) и рассмотрим характеристическое уравнение линейной части уравнения (26) + 1 + k exp(-h) = 0. (27) Выберем k = k0 таким образом, чтобы уравнение (27) имело корни = ±i0 (0 > 0), а остальные корни уравнения (27) при этом имеют отрицательные вещественные части. Положим k = k0 + k1 и обозначим (), (), (() = i0+1+...) соответствующие корни уравнения (26). Считаем, что k1 < 0. При этом Re 1 0.

Пусть теперь aj = 0 (j = 1, 2). Положим j = 20 + j (j 1, j = 1, 2).

Таким образом рассматривается случай двухчастотного параметрического резонанса.

В сформулированных предположениях уравнение (26) имеет в окрестности нуля фазового пространства C(-h(t, ), 0) локальное экспоненциально устойчивое интегральное многообразие, поведение решений на котором определяется поведением решений некоторой двумерной нелинейной периодической системы (нормальная форма уравнения (26)).

Приведем вид и численные результаты системы для конкретного значения параметров. Положим в (26) 1 = 2 = 0, h = 3/4, f2 = 0.1, f3 = -1, = 1, k1 = -1. При этом k0 = 2, 0 = 1. В результате с учетом некоторых нормировок нормальная форма уравнения (26) примет вид = (-1 + 1.6826a1 + 1.6826a2 cos(t))x+ +(-0.1906 + 1.6826a2 sin(t))y + (x2 + y2)(-x + 7.343y), = (1.6826a2 sin(t) + 0.1906)x+ +(-1 - 1.6826a1 - 1.6826a2 cos(t))y + (x2 + y2)(-y - 7.343x).

Рис. 8. Проекция аттрактора на плоскость (x, y) Положим сначала параметр a2 = 0 и будем изменять параметр a1 от нуля в сторону возрастания. Это соответствует периодическому воздействию на систему. При a1 0.6069 нулевое состояние равновесия теряет устойчивость, из которого рождаются два ненулевых устойчивых состояния равновесия. Зафиксируем теперь a1 0.8488 и будем изменять a2. Из этих ненулевых состояний равновесия одновременно бифурцируют при a2 4.244 · 10-3 два устойчивых цикла. Состояния равновесия при этом теряют устойчивость.

Дальнейшее увеличение параметра a2 приводит к увеличению амплитуды колебаний периодических решений. На рис. 8(а) представлены проекции периодических решений на плоскость (x, y) при a1 0.8488, a2 0.8488.

При a2 1.1972 оба цикла теряют одновременно устойчивость, и неустойчивое многообразие первого цикла пересекается с устойчивым многообразием второго, и, наоборот, неустойчивое многообразие второго цикла пересекается с устойчивым многообразием первого. Это приводит к образованию странного аттрактора (хаотического режима), проекция которого на плоскость (x, y) для случая a1 0.8488, a2 1.1973 представлена на рис. 8(б). Для этого случая с помощью программы Tracer23 были вычислены ляпуновские показатели и ляпуновская размерность; 1 0.1932, 2 -2.6211, dL 1.073.

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

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

Список публикаций по теме диссертации Статьи в ведущих рецензируемых журналах:

1) Кубышкин, Е.П. Параметрическое возбуждение хаотических колебаний в одном дифференциальном уравнении второго порядка с запаздывающим аргументом / Е.П. Кубышкин, А.Ю. Коверга // Моделирование и анализ информационных систем. 2008. Т.15. № 2. С. 67–71.

2) Коверга, А.Ю. Хаотические колебания одной распределенной динамической системы с бесконечным запаздыванием / А.Ю. Коверга, Е.П. Кубышкин // Моделирование и анализ информационных систем. 2011.

Т.18. № 1. С. 46–55.

Другие публикации:

3) Коверга, А.Ю. Об одной математической модели полупроводникового лазера / А.Ю. Коверга, Е.П. Кубышкин // Математическое моделирование и краевые задачи: Труды Второй Всероссийской научной конференции / Самарск. гос. тех. ун-т. Самара, 2005. С. 127–130.

4) Коверга, А.Ю. Некоторые особенности поведения решений уравнений Ланга-Кобаяши / А.Ю. Коверга, Е.П. Кубышкин // Современные проблемы математики и информатики: Сборник научных трудов молодых ученых, аспирантов и студентов / Яросл. гос. ун-т. Ярославль, 2005.

Вып. 7. С. 146–150.

5) Коверга, А.Ю. Характер поведения решений нелинейной системы дифференциальных уравнений в одном критическом случае / А.Ю. Коверга, Е.П. Кубышкин // Математические методы в технике и технологиях - ММТТ-20: Сборник трудов ХХ Международной науч. конференции в 10 т. / Яросл. гос. техн. ун-т. Ярославль, 2007. Т. 1. Секция С. 87–88.

6) Коверга, А.Ю. Поведения решений нелинейной системы дифференциальных уравнений в критическом случае кратной пары чисто мнимых корней с жордановой клеткой / А.Ю. Коверга, Е.П. Кубышкин // О работе семинара Нелинейная динамика. Моделирование и анализ информационных систем. 2007. Т.14. № 2. С. 86.

7) Коверга, А.Ю. Параметрическое возбуждение хаотических колебаний в одном линейном дифференциальном уравнении с запаздывающим аргументом / А.Ю. Коверга, Е.П. Кубышкин // Воронежская зимняя математическая школа С.Г. Крейна – 2008. Тезисы докладов / Воронеж:

ВорГУ, 2008. С. 90.

8) Параметрическое возбуждение хаотических колебаний в RC-генераторе с запаздывающей обратной связью / Л.Н. Казаков, [и др.] // Вестник Ярославского государственного университета им. П.Г. Демидова. Серия Естественные и технические науки. 2009. № 1 С. 59–63.

9) Казаков, Л.Н. Хаотические колебания генератора, обусловленные периодическим изменением запаздывания в обратной связи / Л.Н. Казаков, А.Ю. Коверга, Е.П. Кубышкин // Вестник Ярославского государственного университета им. П.Г. Демидова. Серия Естественные и технические науки. 2011. № 1 С. 70–74.

10) Коверга, А.Ю. Хаотические колебания одной нелинейной распределенной динамической системы / А.Ю. Коверга // Воронежская зимняя математическая школа С.Г. Крейна - 2012. Материалы международной конференции / Воронеж: Издательско-полиграфический центр ВГУ, 2012.

С. 102–105.

11) Коверга, А.Ю. Хаотическое поведение решений одного нелинейного дифференциально-разностного уравнения в случае двухчастотного параметрического резонанса / А.Ю. Коверга, Е.П. Кубышкин // Воронежская зимняя математическая школа С.Г. Крейна - 2012. Материалы международной конференции / Воронеж: Издательско-полиграфический центр ВГУ, 2012. С. 105–107.

12) Коверга, А.Ю. Специфика возникновения хаотических колебаний в одной нелинейной распределенной динамической системе / А.Ю. Коверга // Моделирование и анализ информационных систем. Труды международной научной конференции, посвященной 35-летию математического факультета и 25-летию факультета информатики и вычислительной техники Ярославского государственного университета им.П.Г. Демидова. Ярославль, 2012. С. 116–118.

13) Коверга, А.Ю. Некоторые особенности двухчастотного параметрического возбуждения колебаний в нелинейных динамических системах / А.Ю. Коверга, Е.П. Кубышкин // Моделирование и анализ информационных систем. Труды международной научной конференции, посвященной 35-летию математического факультета и 25-летию факультета информатики и вычислительной техники Ярославского государственного университета им.П.Г. Демидова. Ярославль, 2012. С. 119–121.

Подписано в печать 10.04.12. Формат 60х84/16.

Бумага оф. Отпечатано на ризографе.

Тираж 100 экз. Заказ 19/12.

Отдел оперативной полиграфии ЯрГУ 150000, Ярославль, ул. Советская, 14.






© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.