Синергетика_Глава_2_04_Нелинейный_маятник_с_диссипацией_в_электрическом_поле
.pdfисследовано математиками. Решения такого уравнения составляют особый класс специальных функций - функции Матьё. Для наших целей важно отметить свойства неустойчивости такого уравнения. На плоскости параметров - амплитуда частота существуют зоны неустойчивости, которые имеют вид характерных клювов, расположенных в окрестности резонансных
частот, |
Нетрудно показать, |
что переходом от функции 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
