Синергетика_Глава_1_03_Бифуркация_Хопфа
.pdf1.3. Бифуркация Хопфа. Рождение предельного цикла. Мягкая и жесткая потеря устойчивости. Аналогия с фазовыми переходами I и II рода 1.3.1. Бифуркация Хопфа
Бифуркация Хопфа, при которой изменение параметра управляющего параметра приводит к тому, что оба характеристических значения p пересекает мнимую ось, встречается во многих физических системах. При бифуркации Хопфа происходит переход от Re p 0 к Re p 0 . Например, часто она описывает порог генерации в лазерах. Поэтому рассмотрим ее ниже более подробно. Причем не будем ограничиваться какой-то конкретной системой, а разберем сразу достаточно общий случай. Итак, предположим,
что в результате бифуркации Хопфа |
комплексная амплитуда |
X t exp pt , характеризующая состояние |
системы (ее отклонение от |
положения равновесия, ставшего неустойчивым) стала экспоненциально
нарастать со временем t : |
|
p Re p i Im p i ; |
(BHE.1) |
X t Aexp t exp i t ; |
(BHE.2) |
с инкрементом нарастания γ > 0 и амплитудой A. Этот результат, однако, справедлив лишь в линейном приближении, когда отклонения от состояния равновесия малы. Важно отметить, что нелинейные эффекты приводят к тому, что модуль амплитуды нестационарного движения не растет
неограниченно по закону |
|
X t |
|
2 |
|
|
A |
|
2 exp 2 t , а стремится к некоторому |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
конечному пределу. При значениях управляющего параметра λ, близких к критическому значению c , этот предел может оказаться все еще мал, и для его определения можно воспользоваться следующим приемом (Л. Д. Ландау,
1944). Заметим, что в случае, |
рассмотренном в предыдущем |
параграфе, |
||||
величина критического значения |
c 2 . |
|
||||
|
Найдем производную |
по |
времени от квадрата модуля |
амплитуды |
||
|
X t |
|
2 . Для малых временах |
t |
, когда еще применима теория возмущений |
|
|
|
|||||
(линейное приближение), справедливо соотношение
|
|
|
|
|
2 |
|
|
|
|
|
|
|
d |
X |
t |
|
2 |
X t |
2 . |
(BHE.3) |
|||||
|
|
|
|
|||||||||
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|||||||
|
dt |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
35 |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
Это выражение справа следует рассматривать как первый член разложения в
ряд по степеням X и X * . При учете нелинейных эффектов появляются более высокие степени в этом разложении. Ближайшие следующие члены − третьего порядка по X . Нас, однако, будет интересовать не точное значение производной, а ее среднее по времени значение. Причем усреднение подразумевается по промежуткам времени T , с одной стороны, большими по
сравнению с периодом T 2 , а с другой стороны, малыми по сравнению с характерным временем нарастания амплитуды 1/γ. Промежуток времени
усреднения T удовлетворяет соотношению: T T |
1 |
. Поскольку члены |
||||||||||||
|
||||||||||||||
|
|
|
|
|
|
|||||||||
|
|
|
|
|
||||||||||
третьего порядка, имеющие малость O X 3 , такие как X 3 , X * 3 , X |
|
|
|
X |
|
|
2 ; |
|||||||
|
|
|
|
|
||||||||||
X * |
|
X |
|
2 содержат периодические множители exp i t и exp 3i t они |
||||||||||
|
|
|||||||||||||
при усреднении по времени T пропадают. Среди членов четвертого порядка |
||||||||||||||
единственный, не выпадающий при усреднении, это, очевидно, член |
|
X |
|
4 . |
||||||||||
|
|
|||||||||||||
Таким образом, с точностью до членов четвертого порядка малости,
дифференциальное уравнение для |
|
X |
|
2 (BHE.3) изменится: |
||||||||||||||||
|
|
|||||||||||||||||||
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d |
X |
t |
|
2 |
X t |
2 |
X |
4 , |
|
|
|
(BHE.4) |
|||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
dt |
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
где постоянная величина , называемая постоянной Ландау, может быть как положительной, так и отрицательной. Нас будет интересовать ситуация,
когда при λ>λc решение для X t 2 становится неустойчивым уже при сколь
угодно малом возмущении. Для того, чтобы реализовалась такая ситуация должно выполняться условие α>0.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
1. |
Случай |
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
Рассмотрим |
случай 0 . Поскольку |
при λ>λc инкремент |
γ>0, |
то |
|||||||||||||||||
уравнение (BHE.4) имеет устойчивое стационарное решение. Это означает, |
||||||||||||||||||||||
что |
при |
0 |
происходит |
ограничение роста. В уравнении |
(BHE.4) |
|||||||||||||||||
существует слагаемое «прихода» 2 |
|
X t |
|
2 , имеющее |
положительный знак, |
|||||||||||||||||
|
|
|||||||||||||||||||||
и слагаемое ухода |
|
|
|
|
X t |
|
4 , |
имеющее при |
0 |
отрицательный знак. |
||||||||||||
|
|
|
|
|||||||||||||||||||
|
|
|||||||||||||||||||||
Если |
в |
уравнении (BHE.4) |
ввести переменную |
y t |
|
X t |
|
2 , |
то |
|||||||||||||
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
36 |
|
|
|
|
|
|
|
|
|
|
||
дифференциальное уравнение для |
y t имеет вид |
|
dy t |
2 y t |
|
|
|
, так |
|||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
dt |
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
как правая часть уравнения (BHE.4) |
|
d |
|
X |
|
|
d |
|
|
1 |
|
|
|
1 |
|
dy |
, |
|
|
а левая |
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
y2 t dt |
|
|
|||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
dt |
|
y t |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
часть |
уравнения |
|
|
|
(BHE.4) |
|
принимает |
|
|
вид |
|
|
2 |
|
X t |
|
2 |
|
2 |
|
, |
||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
y t |
|||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
X t |
|
4 |
|
|
|
|
|
. |
|
|
Решение |
такого |
дифференциального |
|
уравнения |
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
y2 t |
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
dy t |
|
|
|
|
|
|
|
|
|
|
|
вид y t |
|
|
|
B exp 2 t , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
2 y |
t |
|
|
|
|
имеет |
|
|
|
|
где |
|
|
|
B |
− |
||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||
|
|
|
dt |
|
|
|
|
|
2 |
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
постоянная величина. |
На больших временах t 1 стационарное решение |
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
для y t |
|
|
|
|
|
1 |
|
|
, |
следовательно положительное стационарное |
|||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||
2 |
|
X max |
|
2 |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
значение для |
|
X |
|
2max 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
X 2max
|
|
|
2 |
|
|
. |
(BHE.5) |
|
|
|
|
|
|||
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
При этом значении X max 2 в уравнении (BHE.4) производная обращается в
ноль |
|
d |
|
X max |
|
|
2 |
0 . |
Если |
|
|
|
0 , |
|
то |
при |
|
|
X |
|
2 |
|
|
X max |
|
2 |
|
|
производная |
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||
|
d |
|
X |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
становится |
|
отрицательной, |
|
так |
как |
|
при |
таких |
|
|
значениях |
|||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
dt |
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
X |
|
4 |
|
2 |
|
X |
|
2 . Это |
|
означает, что |
в последующие |
моменты времени |
||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
величина |
|
X |
|
2 будет уменьшаться, |
|
приходя из состояния |
|
X |
|
2 |
|
X max |
|
2 в |
||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
равновесное состояние |
|
X max |
|
2 . При |
|
X |
|
2 |
|
X max |
|
2 производная |
d |
|
X |
|
2 |
будет |
||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
положительной, что означает, что в последующие моменты времени величина X 2 будет возрастать, также приходя в равновесное состояние
X max |
|
2 . Следовательно данное положение равновесия |
|
X max |
|
2 действительно |
|
|
|
является устойчивым.
37
Величина инкремента γ зависит от управляющего параметра λ. Вблизи
своего критического значения c функция γ(λ) может быть разложена по |
|||||
степеням разности c . Поскольку для критического значения |
с 0 , |
||||
то, используя ряд Тейлора, можно записать |
|
||||
|
|
D |
|
c , |
(BHE.7) |
|
|
||||
где величина D играет роль некоторой постоянной. Заметим что функцияявляется отрицательной величиной 0 при значениях c .
Если значения удовлетворяют условию c , то величина 0 . Из
выражения (BHE.5) получается зависимость установившейся амплитуды колебаний от надкритичности при c
X max 
c . (BHE.8)
Такую зависимость можно изобразить графически так рис. 5
Рис. BHE_1. Мягкая потеря устойчивости
2. Случай 0
Рассмотрим случай 0 . Для определения предельной амплитуды в
уравнении (BHE.4) надо учесть отрицательный член более высокого порядка. Для определенности пусть это будет отрицательное слагаемое,
пропорциональное |
|
|
|
X |
|
6 , где коэффициент |
0 : |
|
|||||||||||||||||||
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d |
X |
t |
|
2 |
|
X t |
|
2 |
|
|
|
X t |
|
4 |
|
X t |
|
6 . |
(BHE.9) |
|||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
dt |
|
|
|
|
|||||||||||||||||||||
|
|
|
|
38 |
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
Равновесное значение |
|
X max |
|
2 получается из условия нулевого |
значения |
||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||
правой части уравнения (BHE.9). Разделив правую часть (BHE.9) на |
|
X max |
|
2 , |
|||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||
получаем уравнение для определения равновесного значения |
|
|
|
|
|||||||||||||||||||||||||||
|
|
4 |
|
|
|
|
|
X max |
|
2 |
|
|
|
|
|
|
|
|
|||||||||||||
|
X max |
|
|
|
|
2 0. |
(BHE.10) |
||||||||||||||||||||||||
Корни этого уравнения равны: |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
X max |
|
2 |
|
|
|
|
|
|
|
|
|
2 8 |
(BHE.11) |
||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Эту зависимость можно изобразить следующим образом рис. BHE_2.
Рис. ВНЕ_2. Жесткая потеря устойчивости.
Больший корень (со знаком + перед квадратным корнем) описывает верхнюю часть этой кривой, а меньший (со знаком −) соответственно нижнюю.
Значение |
'c |
|
|
|
/ 2 |
соответствует |
равенству |
|
нулю |
|
|
подкоренного |
|||||||||||||||||||||||
|
|
|
|
|
|||||||||||||||||||||||||||||||
выражения |
в (BHE.11) |
2 8 0 , |
где |
|
|
D |
|
c |
то |
|
есть |
||||||||||||||||||||||||
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
2 8 |
|
D |
|
c |
(BHE.7). При с |
|
|
|
два |
|
значения |
для |
|
X max |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
2 |
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
c |
|
|
|
|
|
|
|
|
|
|
||||||||
(верхняя и |
нижняя |
ветвь) совпадают. |
Величина |
|
|
|
по-прежнему |
||||||||||||||||||||||||||||
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
39
соответствует значению γ=0. При c |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
верхняя ветвь кривой |
X max |
|
|||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
X max |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
c |
|
|
|
|
|
||||||||||||||||||
соответствует |
|
|
величине |
|
|
|
|
|
|
|
|
|
, а |
нижняя |
ветвь |
при |
|
|
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
соответствует |
|
|
X max |
|
|
0 . Легко убедиться в том, |
что верхняя часть кривой |
||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
является устойчивой |
|
|
|
|
X max |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
2 8 |
, а |
нижняя |
часть |
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
неустойчивой |
|
|
X max |
|
2 |
|
|
|
|
|
|
|
|
2 8 |
. |
Она поэтому изображена на рис. |
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X max |
|
|
|
|
|
|
|
||||||||||||||||||
BHE_2 пунктиром. Кроме того, можно показать, что решение |
|
|
= 0, |
||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||
которое тоже обращает в ноль правую часть уравнения (BHE.11), устойчиво |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||
при значениях λ < λc. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Отсюда следует, что при пересечении λ значения c |
(слева направо) в |
||||||||||||||||||||||||||||||||||||||||||||||||||||||||
системе ничего не происходит и она остается в состоянии равновесия c X max . При прохождении же значения c амплитуда колебаний скачком возрастает от нуля до некой конечной величины, определяемой верхней веткой кривой. В интервале c <λ< c , где в наличии имеются два устойчивых решения, состояние системы метастабильно. Одно из решений устойчиво по отношению к бесконечно малым, но неустойчиво по отношению к возмущениям конечной амплитуды.
При уменьшении параметра λ система возвращается в состояние равновесия (скачком) не при λ= c , а только, при λ = c . Другими словами,в
поведении системы имеется так называемый гистерезис: эволюция системы (при изменении управляющего параметра) в одну сторону не совпадает с ее эволюцией в обратную сторону.
Таким образом, в обоих рассмотренных случаях в результате потери устойчивости стационарное состояние переходит в колебательный периодический режим. В первом случае после потери устойчивости амплитуда колебаний пропорциональна квадратному корню из надкритичности. Этот вид потери устойчивости называется мягкой потерей устойчивости, так как устанавливающийся колебательный режим при малой надкритичности мало отличается от состояния равновесия.
40
Во втором случае перед тем, как установившийся режим теряет устойчивость, область притяжения этого режима становится очень малой, и всегда присутствующие случайные возмущения выбрасывают систему из этой области еще до того, как область притяжения полностью исчезнет.Этот вид потери устойчивости называется жесткой потерей устойчивости. При этом система уходит со стационарного движения скачком и перескакивает на иной режим движения. Другой характерной особенностью поведения системы в этом случае является уже упомянутый выше гистерезис.
1.3.2. Аналогия с фазовыми переходами
Между двумя перечисленными выше механизмами потери устойчивости и фазовыми переходами можно провести определенную аналогию. Для этого
обозначим z |
|
X t |
|
2 |
и запишем уравнения (ВHE.4) и (BHE.9) в несколько |
|
|
ином виде, введя потенциальную энергию V z .
1.Для «мягкого» режима» 0 уравнения (BHE.4) можно переписать в виде
dz dV z ; dt dz
V z z2 |
z3 |
; |
|
|||
3 |
|
|
||||
|
|
|
|
|
|
|
|
dz |
2 z z2 |
. |
(BHE.12) |
||
|
dt |
|||||
|
|
|
|
|
|
|
2.Для «жесткого» режима 0 ; 0 уравнения (BHE.9) можно переписать в виде
|
|
dz |
|
dV z |
; |
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
dt |
|
|
dz |
|
|
|
|
|
|
|
|
||||||||||||
V z z2 |
|
|
|
|
z3 |
|
|
|
|
|
z4 |
; |
|||||||||||||
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
2 |
|
|
|
4 |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
dz |
2 z |
|
|
|
z2 |
|
|
|
|
z3 . |
|
|
|
|
(BHE.13 |
||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||
dt |
|
|
|
|
|||||||||||||||||||||
41 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
В выражениях (BHE.12), (BHE.13) функция , зависит от управляющего параметра по закону (BHE.7): D c . Отметим, что уравнения (BHE.12) и (BHE.13) описывают передемпфированное движение частицы в
поле с потенциалом V z . Величина |
F z |
dV z |
играет роль силы |
|
dz |
||||
|
|
|
трения. При такой формулировке становится ясным, что устойчивые состояния равновесия − это точки минимума функции V z рис. (BHE_3). При этом возбуждение в мягком режиме аналогично фазовому переходу II
рода, где |
параметр |
порядка M |
Tc T |
не равен нулю ниже точки |
переxода, |
где T - |
текущая температура, |
а Tc - температура фазового |
|
перехода. |
рис. BHE_3. |
|
|
|
Рис. BHE_3: V z для мягкого режима.
А жесткий режим напоминает фазовый переход I рода с присущим ему гистерезисом рис. BHE_4.
42
Рис. BHE_4. V z для жесткого режима.
43
