Добавил:
tg: @Yr66gi4 Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Синергетика_Глава_2_04_Нелинейный_маятник_с_диссипацией_в_электрическом_поле

.pdf
Скачиваний:
4
Добавлен:
22.09.2025
Размер:
590.05 Кб
Скачать

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

частот,

Нетрудно показать,

что переходом от функции U x

к функции

Z x

(U (x) Z x ),

где

U x Z x exp x

линеаризованное

уравнение (EMC.17) можно записать для искомой функции

Z x

 

 

Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

a 2q cos

 

2x

Z

 

x

 

0

 

(EMC.17)

Константа a в

уравнении

(EMC.17)

равняется a b 2 .

Уравнение

(EMC.17) представляет собой классическое

 

дифференциальное уравнение

Матьё. Таким образом, для того чтобы решение уравнения (EMC.16) было неустойчивым, необходимо, чтобы соответствующее решение уравнения (EMC.17) нарастало бы по закону Z x exp px , где параметр p .

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

Исследуем возможность экспоненциального нарастания во времени решений классического уравнения Матьё y t exp t (EMC.17). Функция

частоты t в этом уравнении является периодической функцией времени

t T t .

Период этой функции

T

определяется

частотой

приложенного периодического поля

T 2 / .

Следовательно уравнение

(EMC.15) инвариантно

относительно

преобразования

t t T .

Отсюда

следует, что если

y t

- решение уравнения (EMC.15), то

и y t T также

есть решение этого уравнения. Дифференциальное

уравнение (EMC.15)

представляет собой обыкновенное дифференциальное

уравнение, у которого

есть два линейнонезависимых решения y1 t

и y2 t . Если y1 t

и y2 t

- два независимых решения дифференциального уравнения, то при замене

t t T

эти решения

преобразуются линейным образом

друг через друга.

Можно выбрать

эти решения y1 t и

y2 t

таким образом, чтобы их

изменение

при

замене

t t T сводилось

просто к

умножению на

постоянный множитель

80

y1 t T 1 y1 t ; y2 t T 2 y2 t . (EMC.18)

Наиболее общий вид этих функций, обладающих таким свойством, можно выбрать в виде

 

t

 

 

t

 

 

 

 

 

 

 

y1 t 1T 1

t ;

y2 t 2T 2 t ,

(EMC.19)

где функции 1 t и

2 t

- периодические

функции времени,

удовлетворяющие условиям

1 t 1 t T ;

2 t 2 t T .

Покажем, что постоянные 1 и 2 в этих функциях должны быть связаны друг с другом определенным соотношением. Мы имеем два уравнения y1 t 2 t y1 t 0 ; y2 t 2 t y2 t 0 . Умножим первое уравнение

на y2 t

и второе уравнение на

y1 t

и

вычтем эти

два

уравнения. В

результате

мы получим

соотношение

y1 y2 y2 y1

 

d

 

y1 y2

y1 y2 0.

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1 y2 y1 y2 Const .

 

 

 

 

 

 

(EMC.20)

Учтем свойство периодичности функций 1 t

и 2 t ,

и

запишем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

y1 t

 

 

 

 

 

t ;

значения функций в момент t и в момент времени

t T :

 

1T 1

 

 

t

 

 

t

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

y1 t T 1 1T 1 t ;

y2 t 2T 2 t ,

y2

t T 2 2 T 2

t .

Вычислим

производные

этих функций в момент времени

t ,

и в момент

времени

t T . Запишем равенство

(EMC.20)

в

момент

времени

t :

y1 t y2 t y1 t y2 t Const ,

и

в

момент

времени

 

 

 

 

t T :

y1 t T y2 t T y1 t T y2 t T Const .

Условие

 

сохранения

величины (EMC.20) требует чтобы между константами 1

и 2 выполнялось

бы соотношение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2 1.

 

 

 

 

 

 

 

 

 

(EMC.21)

Из условия вещественности решений y1 t и y2 t вытекает, что константы

1

и 2

являются вещественными числами.

Таким

образом, 1 ,

2

1 / ,

где 1. Следовательно одно решение

y1 t

растет со времени

 

 

81

 

 

 

t

t

 

 

 

 

 

 

 

 

 

 

 

 

 

по экспоненциальному закону y1 t T 1

t exp

 

ln

1

t , а

 

 

 

 

T

 

 

 

другое решение экспоненциально затухает со временем

 

 

 

t

 

 

 

 

y2

t exp

 

 

ln

2

t . Таким образом,

мы доказали, что уравнение

T

 

 

 

 

 

 

 

Матье проявляет

свойство неустойчивости, и

достаточно малое отклонение

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

которому решение y1 t

экспоненциально возрастает со временем.

Отметим основные

свойства решения уравнения Матье (EMC.17)

(Кузнецов А.П., Кузнецов С.П., Рыскин Н.М, Нелинейные колебания, издание 3, М., Ленард, 2020, Лекция 16; Мэтьюз Дж., Уокер Р. Математические метолы физики, Москва, Атомиздат, 1972 ) На некоторой плоскости параметров, для которых осью абсцисс является ось величин a , а осью ординат является величина q существуют зоны неустойчивостей, которые имеют вид характерных клювов, расположенных вблизи величины

a n2 , где n 1,2,3... Эти устойчивые области изменения параметров q и a уравнения Матье (EМC.17) закрашены зеленым цветом. Области неустойчивого решения уравнения Матье имеют белый цвет.

Рис. EMC_3. Границы зон неустойчивости на плоскости параметров a и q для классического уравнения Матьё. Цифры 1,2,3,4, 5 соответствуют номерам резонансов. Зеленым цветом закрашены устойчивые области. При учете затухания 0 границы этих областей сдвигаются наверх. Пунктирными линиями обозначены границы устойчивости при 0 .

Первый резонанс наблюдается при a 1. Область устойчивости при a 1 очень узкая и распространяется только для очень малых q . В наших

82

расчетах мы считали, что a b 2 =1.0009−(0.03)2=1. Эта величина в точности соответствует первому резонансу уравнения Матьё ( рис. 1_3_2). Для нашего численного расчета нелинейного уравнения Матье с затуханием мы получили, что границы устойчивости нелинейного уравнения Матье представляют величину q 0.5 ( случай 6).

Вывод

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

83