Файл: Нелинейные системы гидродинамического типа..pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 29.10.2024

Просмотров: 85

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

метра jx не нарушает общности рассмотрения, так как, делая замену переменных zJ^Xz^ z2= X 2z2 (X — вещественно), в новых переменных полупим ц' = Х2ц, и это значение может быть сделано равным любому заданному числу (кроме естественного значения единица может оказаться полезным значением fi= Ѵ 2, как это будет показано ниже).

Наряду с системой уравнений (6. 4) следует иметь в виду эквивалентную систему, получаемую операцией комплексного

сопряжения

1’

t =

(6.5)

Легко убедиться, что £'=(fxz1^1+ z 2z2)’является интегралом движе­ ния — фундаментальной формой.

Рассмотрим замену переменных

где фазы

Ѳх и

z' =

e<e*Zj,

z' =

efe»z2,

Ѳ2 — вещественны. При такой подстановке система

сохраняет вид,

при этом

к1 == Ä:e‘*(2ві+о2)

01=<р/3, Ѳ2= — 2 ср/3, можно

Если7с=|/кс'|е-<<р, то, положив, например,

сделать

вещественным.

 

системы является совокуп­

Таким образом, основной ггруппой

ность таких преобразований, при которых фазы удовлетворяют

условию

 

 

2Ѳ1+ Ѳ2к= 0.

(6.J6)

при этом условии

коэффициент

остается вещественным,

если

он был

таковым.

В дальнейшем,

не нарушая общности,

будем

считать

к =

1.

 

которые оставляют неизменным

Совокупность преобразований,

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

02 = О, 01== - А .

В отличие от вещественного случая комплексная система допускает еще один — кубический интеграл движения.

Умножим первое уравнение на 2zjZ2, а второе — на pzf и сложим. Тогда получим вещественную величину

Pgj(*!**) = —;212? I •12 ! I + V I z i I*

Таким образом

И = Im (zfz,)

(6.7)

91


является интегралом движения, который в развернутой форме имеет вид

II = 2 — i/f) у, + 1у1 • .г,.

(6. 8)

Уравнения движения (6. 4) также запишем в развернутой форме

dxi

I

•‘ ddt = —^1*2+ №

y,

х і2/з +

V- -^г

(6.9)

d'Ji _

dt — %С\Уі-

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

дхг

ду^

. дх2 . ді/.2

_ q

д хх "т- дуг

дх» г" ду2

 

Более того, при соответствующем выборе параметра р (р=1/г) соответствующие уравнения движения можно записать в форме Гамильтона. Чтобы убедиться в этом, вычислим производные от кубической формы Н

^ = 2

{хлу2

+

х 2ух),

^ - = 2 ( —

У і У . + х ^ ) ,

 

 

 

дН 0

 

 

 

дН_

 

 

dxz- ~ 2 x ilJ' ’

 

 

д уі— 0*2 .

у!-

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

Х і —

Рі7

УдНі =

/?дрі

1,.

*q22 =3дН2 .

У г = Р ѵ

 

Уі

/др2,

 

Pi

=

=

 

 

p2

=

dHjdq2,

(6. 10)

 

~ d H ld q v

 

 

 

 

 

 

 

в новых переменных H имеет вид: Н = р 2 2q\) ~{-2p1q2q1. Таким образом, подстановка z1 = p 1 -\-iq1, z2 = q 2 -\-ip2 приводит нас к си­ стеме, в которой функция Н является функцией Гамильтона.

Иногда бывает удобным рассматривать в качестве основной переменной сопряженное значение второй координаты — w = z %.

Каноническая форма уравнений движения имеет вид

dz/dt

=

—2

wZ,

dw/dt

(6. И )

 

=

z2.

 

Дифференцируя второе уравнение, получаем для w уравнение второго порядка с переменным, но вещественным коэффициентом

92


cPw/dt*-\-4\dw/dt\w—Q. Решив это уравнение при заданных на­ чальных условиях w\t=Q=w0, (dwldt)\i^—zl, находим, что

z = +\Jdw/dt-

Квадратичный интеграл, выраженный через функцию w, имеет вид

 

 

і? =

И

2 +

 

dw

 

 

т dt

а гамильтонова функция

записывается в форме

 

 

Н —

 

WZ

 

Im

[w ijp j.

Полагая

w=pei0,

 

Im ( 2) =

 

 

получим

 

H ~

Im

 

d w ___dp jO + ZpeiOdO

 

 

 

dt

d t 6

 

d t »

 

 

Наличие двух интегралов движения позволяет получить реше­

ние в квадратурах *.

 

 

вид

Уравнения движения в новых переменных =имеюто,

d2p

 

0.

 

dt'2 Р

 

(6. 12)

dO

d20

 

Z dt dt +

Pdt2

 

Второе уравнение (6. 12) непосредственно интегрируется и при­ водит к интегралу /f=const. Пользуясь этим интегралом, по­

лучаем

dt2

^ я * + 4 | / я м - р * ( і ) ! = 0 .

(6. IS)

Легко

d2p

проверить,

что «энергия»

 

E=f'+T,YH'+ K £ f

действительно является интегралом уравнения (6. 13). Рассмотрим частный случай Н = 0, тогда

S = 0 , О= 0о = const.

Не нарушая общности, благодаря возможности фазовых преоб­ разований, будем считать Ѳ0= 0 . Таким образом, система ока­ зывается вещественной, а для этого случая мы уже имеем общее решение

w

,у

111 at,

=

а 1

2

2 ch at ’

где а2 = 4Е, причем ршіп= 0 .

 

(В. 14)

* Эти квадратуры выражаются через эллиптические интегралы.

93


Фаза 0 меняется скачком

от —тс до тс при

переходе

через точку

і=0.

H=jt=

 

 

 

 

 

При

начало

отсчета

времени

в

момент,

 

0 выберем

когда р=р*, определяемому из кубического

уравнения на

основе

интеграла

энергии при условии

= О

 

 

(6.15)

 

 

 

 

 

 

 

Положим

P.= f

r,

ы /г==^-^>

0< г*< 1

 

Тогда

гі —

+

2h— 0,

2іг = ?-у(1 — rf).

 

Отсюда следует, что

\ІІ\

не может

превзойти

IInmx — h*a3l (6.16)

 

 

 

2, где

 

 

 

 

 

II

тал

 

 

 

=

1/\/3, т.

е.

/г* = х/3\/3 — значение, соответствующее

 

 

 

 

 

 

 

 

аЗ

 

 

 

 

 

В этом случае

dr/dt=

 

 

 

 

бТІ'

и система совершает гармо­

 

 

0, й2ѲДЙ2= 0 ,

нические

 

 

 

 

 

Н ш а х ____ I

2 д

 

 

 

колебания с частотой

~ - ^ з '

 

 

 

 

 

 

со =

 

II,р2

 

 

 

 

w = —^=ae,a>l,

 

dw

 

 

 

 

 

 

 

Знак m определяется знаком

 

~~ 3

 

'

si gn//,

 

 

 

' 2 ѴЗ “ “

 

 

 

 

 

 

 

 

 

 

 

 

dt

=

-тс-e

 

 

 

 

 

 

-%[dw

,

а

•(лг+г)

 

 

при 77 )> 0,

 

z = K B - = ± aj-f(,!2 .

 

 

 

 

 

 

 

 

 

 

4

\ 3

"T" 3

 

4 '

 

Энергия

* = м « + т | г г

 

 

 

2 : 1 между

w- и z - k o m -

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

 

в

отношении

понентами.

 

ситуация

показана на рис. 26.

При 0 < | Я |< Я пих

Амплитуда р* меняется за один полупериод в пределах от гх (/г)

до г2 (Л), являющихся корнями уравнения

(6. 16). При этом

 

 

 

о <

Гі <

 

1 /ѴЗ <

Ра < 1.

 

 

 

Полупериод обмена энергией между компонентами выражается

интегралом

 

 

і_

(

___________dy___________

 

 

 

~ 2 а

</2= 2(/*)

 

 

 

Г _

 

г

Іу (1 — г/)2— 4*2

(6.17)

 

2

 

 

)

 

 

94


После того как определено р (t) обращением интеграла

 

а

У=(Р/«)’

'SУ(1 - U)~ - 4Аа ’

 

1

г

dij

t

_________

 

2

]

 

S'j-'-jtA)

изменение фазы находится дополнительной квадратурой

0 = // \

(6.18)

J Р"

О

При уменьшении h период системы возрастает, и при Л—>0 оно переходит в апериодический режим, рассмотренный выше.

Рис. 26. Графическое реше­ ние уравнения (6. 16) (г* как функция h)

Описанная модель может быть применена к задаче о колеба­ ниях двух параметрически связанных вибраторах, рассмотренной в работах Цельмана (1970, 1971).

Рассмотрим теперь комплексное расширение триплета. Пара­ метры его будем считать вещественными и равными по абсолют­ ной величине единице. Это всегда можно сделать, выполнив масш­ табное преобразование.

Соответствующие

уравнения движения запишутся в форме

 

 

 

 

dz^jdt

-

——

 

 

 

 

 

 

 

dzjdt

 

 

 

 

 

 

 

dz3/dt

=

z3zv

 

 

(6.19)

 

 

 

 

z3z3

 

Эта

система

имеет

 

=

 

движения

два независимых

интеграла

 

 

 

Е = \ { \ * А І +

 

Ц * , ?

+ \*3\1),

 

(6.20)

Система не меняет своего вида при замене переменных

z[ = eiQ'Zlt

z' =

e*'° Z2,

 

= e‘'°3z3

5 = К І а- К І а-

 

z'3

при условии

 

 

 

 

 

 

2

 

®i +

+

^3 =

0.

 

 

 

 

 

 

 

 

Проводя фазовое преобразование, можно добиться, чтобы в задан­

ный момент времени при

t =

0 z1=|z1|=p(°> и z3=p(°). Фаза «сред-

 

95