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

книги / Методы и приёмы качественного исследования динамических систем на плоскости

..pdf
Скачиваний:
16
Добавлен:
12.11.2023
Размер:
16.72 Mб
Скачать

 

АВТОПОДСТРОЙКА

ии ф'(й) = 0 , rti,mm(A) =

0 , т. е.

рт = 2 / Л 0(и +

Л0),

ft2

/ / _

i t

391

(6)

и Ф имеет указанное выше значение.

Условие ipi(0) = 0 дает бифуркационную кривую, соответст­ вующую возникновению петли сепаратрисы сшитого седла:

При убывании ао (или возрастании Ко) от значения, опреде­ ляемого (7), величина i|>i(0) становится положительной, и из петли сепаратрисы сшитого седла рождается неустойчивый пре­ дельный цикл.

Из выражений (6 ) и (7) видно, что при фиксированном [J обе бифуркационные кривые будут прямыми, проходящими че­ рез начало координат плоскости (ао, Ко)- Между ними располо­ жена область, для точек которой система ( 1 ) имеет два предель­ ных цикла, охватывающих фазовый цилиндр.

2. Рассмотрим систему

% - У , &(ф) = ^(ф + 2 л), (8)

отличающуюся от рассмотренной выше наличием члена —Ау в первом уравнении (А. > 0).

Сохраняя принятую аппроксимацию F(cp) кусочно-постоян­ ной функцией, принимающей значения ± 1 , и вводя малый пара­ метр ц: а = ра0, K = PTfo, А = рА0, для функции ф(А, А0), корни которой соответствуют предельным циклам, охватывающим ци­ линдр, получим выражение

Здесь i|>i(A) то же,

что

и в предыдущем примере, а (А, Ао)

дается выражением

2 1

 

ф2 (fc, А0) =

[(2h f * - (2л + 2h)3/2] < 0.

-J

Имеем также

 

 

— оо.

26*

392

МЕТОД ПОНТРЯГИНА В КУСОЧНО-СШИТЫХ СИСТЕМАХ

[ГЛ. 18

 

Кривая 1р2 = гМЛДо) для любых Яо, отличных от нуля,

ухо­

дит в бесконечность при h -*■°°, но будет лежать внутри полосы сколь угодно малой ширины е на любом заданном интервале

изменения h, если выбрать Хо достаточно малым.

\pi(/i) +

Проследим

за изменением

числа

нулей

функции

+ тЫ^До) = ф(ЛДо)

на интервале 0

< Л < ° °

при возрастании

параметра Хо от нуля.

е <

I 'ф1 т т (Л*) I. Если Ао

выбрано

так, что

Пусть е <

I ipi (0) I,

на интервале 0 < h < h*

кривая 'фг(ЛДо) лежит внутри полоски

шириной е (—е < ф 2

< 0 ), то

число

нулей функции ipi + тр2 по

сравнению с ipi не может измениться ни за счет изменения знака ipmin в точке минимума, ни за счет изменения знака \|3i(0 ). С дру­ гой стороны, как бы ни было мало Хо, функции г|ц + фг и г|ц име­

ют при достаточно большом h разные знаки Ш т

= 2 лу0»

\ h ~ * оо

 

lim (TJJJ + ф2) =

— оо). Очевидно поэтому, что при возрастании Хо

h - ю о

)

 

нуль функции ф] + фгТак

от нуля из бесконечности появляется

 

/

/

0 , то при этом у системы

как в нуле, очевидно, будет % +

ф2<

(8 ) из бесконечности появляется устойчивый предельный цикл.

Если ао

и

[} выбраны так,

что \|3i,mi n < 0 , i|Ji(0)>0 и Яо>0

и мало, то система (8 )

будет

иметь три предельных цикла (два

устойчивых п один неустойчивый).

 

 

§ 3. Автоколебания синхронного мотора. Рассмотрим систему

d(p/dt =

у,

dy/dt = Т simp — р (а + р cos 2<р — if sin Ц>)у,

 

 

 

 

а, р,

к, Т >

О,

 

где р — малый параметр

[59, 60, 67].

 

 

В качестве фазового пространства будем рассматривать по­

лосу, заключенную между прямыми <р = —я

и ср = я; точки этих

прямых, имеющие одинаковые координаты у, отождествляем.

Полагая

(рис. 203)

2 я -1<р — 2 ,

— я ^ ф < 1 — я/2 ,

 

 

 

sin ф «

s (ф)

 

2 я _1 ф,

— я/2

< ф < я / 2 ,

 

 

 

2 я _1ф -f- 2 ,

я / 2

ф я;

 

 

 

— 4я_1ф — 3,

— я ^ ф <: — я /2 ,

cos ф «

с (2 ф)

 

4я_1ф + 1,

- я/ 2 ^ ф ^ О,

— 4я_1ф + 1,

0

^ ф ^ 2 /л,

 

 

 

 

 

 

 

4я-1ф — 3,

 

«Ss Y •

получим систему, близкую к кусочно-линейной консервативной системе

dq>/dt = у = Р, dy/dt = Т - s(<p)- IL[<X + $с(2<р)~ 4s((p)]y = Q (1)

(р — положительный параметр).

§ 31

АВТОКОЛЕБАНИЯ СИНХРОННОГО МОТОРА

393

 

1. Консервативная система.

 

 

 

d<p/dt = у,

dy/dt = Т — s(q>).

(2)

Система (2) имеет замкнутые фазовые траектории при Т < 1, Замкнутые фазовые траектории системы (2) охватывают состоя­ ние равновесия. Другие траектории представляют собой спирали, накручивающиеся на фазовый цилиндр (рис. 204).

Сшивание траекторий консервативной системы внутри полосы

(—я, я)

(кусков эллипсов

и гипербол)

осуществляется

на пря­

мых ф = ±я/2.

Сшивание

траекторий

системы

(1)

происходит

также и на прямой ф = 0 .

 

 

 

 

 

 

 

 

Консервативная система (2) имеет интеграл

 

 

 

 

 

н ( ф, »)— £

+

 

 

 

 

 

 

 

 

 

 

 

 

 

F

(ф + я + ip)

+ 2лТ

 

 

 

.

.

 

Я

 

 

 

 

Я <

ф

<

— - у ,

+

 

 

 

 

 

 

 

 

 

Я _

,

 

 

 

 

 

 

 

=

й’

у <

ф

<

т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

--7 Г

Ф~ л + "г)

 

 

 

■f < Ф < Я .

 

 

 

 

 

 

 

 

 

 

Замкнутым

кривым

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

изменения

k

от

 

I

 

 

 

(состояние равновесия типа

центр)

до

h — 0

h — ----5 - я (1 — Г)2

(петля сепаратрисы седла). Точка пе­

 

 

 

 

 

 

ресечения

петли

сепаратрисы с

осью

 

 

 

 

 

 

ф лежит

в

полосе

(0 , я /2 ) , если

2

 

 

 

 

 

 

—У2 < Г < 1

; в полосе

(—я/2 , 0),

если

 

 

 

 

 

 

3 —2 У2

< Г <

2

—V2 ;

в

полосе

(—я,

 

 

 

 

 

 

—п/2), если 0 < Г < 3 —2У2. В зависи­ мости от Т изменяется число сшива­ ний для периодического движения. В дальнейшем ограничимся случаем, когда периодическое движение сшива­

ется из траекторий не более чем трех полос ( Т > 3 —2У2). При этом будут исчерпаны все возможные случаи разбиения фазового пространства на траектории.

394

МЕТОД ПОНТРЯГИНА В КУСОЧНО-СШИТЫХ СИСТЕМАХ

[ГЛ. 18

2. Состояния равновесия. Седловая величина. Система (1)

имеет при Т < 1 два состояния равновесия: Oi(nTI2,0) и 0г(я — яГ/2, 0). Состояние равновесия 0\ будет неустойчивым фокусом, если для точки 0 \ будет

(^ ф + = — р (а + Р — 2РТ — уТ) = — рог > 0, а < 0,

устойчивым фокусом, если а > 0, и центром, если а = 0. В послед­

нем легко убедиться:

при условии а = а + р —2$Т уТ = 0 си­

стема (1 ) на полосе 0

< ф < я / 2

может быть представлена в виде

а,

( г — И

“ - » 1 1 2 1 1 + *)■']

 

dq>

 

у

 

и проинтегрирована.

 

 

величина

Состояние равновесия Ог — всегда седло. Седловая

будет иметь значение

 

 

 

 

(^*ф + Qy)i = — F 7-

(3)

Существенно, что выражения £*Ф+ Qy для фокуса и для сед­ ла имеют одинаковое значение. Как будет показано, это опре­ деляет существенные особенности структуры разбиения фазового пространства и пространства параметров.

3. Предельные циклы. Если систему (1) записать в виде

dyldt = Ну, dy/dt = — Яф + р9 (ф, у),

(4)

где ф(ф, !/) = [—а — рс(2 ф)+фз(ф)]1/, то значения константы

ho,

выделяющие кривые Снаконсервативной системы, вблизи которых

при

малом р будут

предельные циклы системы (4), определя­

ются как корни уравнения ф(й) = 0 , где

Ф ( h )

= f J Я у d<f> d y =

—j j [a + Рс(2ф) — у«(ф)] cfy d y =

= — J fa + pc (2ф) — ys (ф)] у dq>. (5) ch

Чтобы выполнить интегрирование, нужно сделать определен­ ные предположения о расположении петли сепаратрисы консер­ вативной системы на плоскости и тем самым о числе кусков, из которых сшивается кривая

а) Пусть

2 —У 2 < Г < 1 .

При изменении h

в

интервале

—я(1 — Г)2^

< /г< —я(1 — Г)2/4

кривые

Ch будут

эллипсами,

целиком расположенными

в

полосе

(0 , я/2 ) ,

а

интервалу

(—л/4) (1 —Т)2< h < 0 будут соответствовать замкнутые кривые, сшитые из кусков эллипса и гиперболы, расположенных соответ­ ственно на полосах (0 , я /2 ) и (я/2 , я ) .

§ 3] АВТОКОЛЕБАНИЯ СИНХРОННОГО МОТОРА

Выполняя интегрирование, получим

ib(h) =

| — o / ? « ( fe) ’

— а 2 <

ft < — а 2/ 2 ,

'

i - o / ^ f t ) ,

— a 2/ 2 <

fe < 0 .

Здесь

395

( 6)

 

a2=

4- (! -

 

Л 2.

 

(A)V =2n 5/2(a2+'

A),

^ (ft) = / 2 ^ [ a

/ 2 f t +

a 3 +

 

 

?

 

+

(A +

 

+ arcsin

/ 2 (ft + a2)j

V 2 h + a2 + a . j

 

 

 

 

 

 

 

 

 

 

Имеем также

 

 

 

 

 

/ 2

й a2++ a

 

 

Л

,

Q,

. i

^(А ) = / 2 л -=- +

arcsin

w-___

=r + In

V —2 ft

 

 

2

 

V 2 ( h +

a2)

r," ,,.

 

lA 7 2 a /2 f t + a2,

 

 

 

 

 

------------ й г г ?)—

 

 

 

 

 

Рг(— a2/2) =

/

л/2а2л,

F[ (— a2/2) = / 2 я 3/г,

 

 

 

 

(— ®2 /2 ) =

0 .

 

 

Таким образом, функция ф(А)

сшивается из линейной функции

OFQ иа

интервале —а2< ft < —а2 / 2

и монотонной

на интервале

—а2 / 2 < ft < 0 , имеющей с функцией —о/о одинаковый знак про­

изводной.

Функция

Ч>(А)

буДет убывающей при a > 0 , возраста­

ющей при

о < 0 и

будет

тождественно обращаться в нуль при

о = 0 .

 

 

_

Предельных циклов в рассматриваемом случае 2 —У 2 ^ Г < 1 нет. При о = 0 существует внутри петли сепаратрисы континуум замкнутых кривых. При о ¥=0 характер устойчивости фокуса однозначно определяет структуру разбиения фазового цилиндра на траектории. На рис. 205, 1, 14, 4 изображены структуры разбиения фазового пространства соответственно для случаев о < 0 , о = 0 и о > 0 ._

б) Пусть 3 —2V2<7, *£2 —У2. Фазовые траектории консер­ вативной системы располагаются при этом не более чем на трех полосках (не выходят на полосу (—л, —л/2)). Положим также для определенности, что Т > 0,5 (случай Т < 0,5 аналогичен слу­ чаю Т > 0,5, но отвечает другой последовательности интервалов изменения ft, соответствующих определенному числу сшиваний); для Т < 0,5 будет

- (1 — т)2< — -£ (Г2 — 4Г + 2) < — (1 — т? < о;

396 МЕТОД ПОНТРЯГИНА В КУСОЧНО-СШИТЫХ СИСТЕМАХ [ГЛ. 18

для Г > 0 ,5 будет

_ JL (i _ Tf < - ■ £ ( 1 - T f < - - f ( Г 2 - AT + 2 ) < 0 .

Выполняя

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

формуле (5) для

случая

0,5 <

.< Т < 2 —V2, получим

 

 

— a2 < f t <

— а2 /2 ,

 

 

 

 

oF0(h),

 

 

 

ijf(ft) =

oFi(h),

 

— а2 / 2 <

ft <

— ft2 /2

,

/7\

 

 

 

F . W - o F ^ h ) ,

 

- f t 2 /2 < f c < 0 .

 

 

 

 

Здесь ft2 =

~ - л (Т2АТ + 2),

F0(h)

и Fi(h)

имеют

значение,

указанное в

(6 ), и Fi(h)

будет

 

 

 

 

 

 

 

 

 

Ft {К) =

2 (2 а2 -

ft2) /

2 Л +

ft2 + -i-(2A +

ft2 /

2

 

 

 

 

 

 

 

 

 

 

— 2(к + а2)У 2а2Ьг(лЦ 2 arcsin

 

——

Имеем также

 

 

 

 

 

 

 

 

 

 

 

 

V 2(h+a*)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F, (h) = 4 (2h +

ft2/

2 - 2 / 2 а 2

-

b2 л — 2

arcsin

/

2 а2 - ь 2 \

 

 

 

 

 

 

 

 

 

 

 

/

2

(fe+

а2)/

Fl (h)

2 /

2 fe +

ь2

 

 

 

 

 

 

 

 

 

 

 

М -а2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F2( -

ft2 /2) =

F't ( -

ft2 /2) = F; ( -

Ь212) = 0,

F%( -

ft2 /2) > 0,

Из последнего следует, что / ’г(Л) на

интервале

—b2/ 2 < h ^ 0

есть функция, принимающая положительные значения, монотон­ но возрастающая и обращающаяся в нуль в единственной точке h = - ft2/2 .

Выражения (7) определяют ф(Л) как непрерывную функцию, сшитую из трех кусков с различными аналитическими представ­ лениями на интервалах (—а2, —а2/2 ), (—а2/2 , —ft2 /2 ) и (—ft2/2 , 0 ). Функция г|з(Л) на интервале (—а2, —ft2 /2 ) имеет знак, противо­

положный знаку о, обращаясь при а = 0 в нуль

тождественно.

На интервале (—ft2 /2,0) функция ty(h) состоит из

двух слагае­

мых, одно из которых положительно, а другое имеет знак, про­ тивоположный знаку о. Выбором о значение ф(й) на всем интер­ вале (—ft2 /2 , 0 ) может быть сделано как положительным, так и отрицательным.

Проследим за корнями функции ф(й) и изменением качест­ венной структуры разбиения фазового пространства на траекто­ рии при возрастании о от отрицательных значений. При о < 0 функция ф(й) на всем интервале (—а2 , 0 ) будет положительной.

6 3] АВТОКОЛЕБАНИЯ СИНХРОННОГО МОТОРА 397

Предельных циклов нет. Фокус неустойчивый. Структура раз­ биения фазового пространства изображена на рис. 205,1.

При о = 0

функция O|J(A) тождественно обращается в нуль

на интервале

(—а2, —Ь2/2 ) и сохраняет положительное значение

на интервале (~Ь2/2,0). Существует континуум замкнутых кри­ вых, окружающих состояние равновесия типа центр. Петля се­ паратрисы не может существовать, так как о|з(0 )¥=0 ; (о-сепарат- риса седла скручивается с границы области, заполненной замк­

нутыми кривыми; «-сепаратриса седла

уходит в

бесконечность

по верхнему полуцилиндру (рис. 205,7—2).

 

о

достаточно

При

о > 0

фокус становится

устойчивым. Если

мало,

то

на

интервале

(—а2, —Ь2/2) будет

а|з(й)<0 ,

но сохра­

нится

тр (0) > 0.

На

интервале

(—Ь2/2,0)

будет

существовать

h = h\ — корень

функции т|з(h),

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

н е у с т о й ч и ­

во м у

(тр' (hi) > 0 )

предельному циклу, возникшему из границы

области, заполненной замкнутыми кривыми (рис. 205,2).

При возрастании

о кривая ар (А)

будет опускаться — при этом

корень

о|:(h)

будет

возрастать — и

для

достаточно

больших h

функция а|з(й) становится отрицательной на всем

интервале

(—а2,

0). Предельных

циклов

нет. Фокус

устойчив.

Структура

398

МЕТОД ПОНТРЯГИНА В КУСОЧНО-СШИТЫХ СИСТЕМАХ

[ГЛ. 18

разбиения фазового пространства

однозначно

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

(рис. 205, 4).

переходе

от структуры

 

Исчезновение корня h = hi при

рис. 205,2 к структуре рис. 205, 4 может происходить либо при возрастании ho до нуля (этому соответствовало бы влипание не­ у с т о й ч и в о г о предельного цикла в петлю сепаратрисы седла), либо при слиянии его с другим корнем функции i|i(h). Не про­ водя подробного исследования поведения функции г|>(/г) при раз­ ных о, можно по знаку седловой величины заключить, что при

возрастании о реализуется именно последний случай.

(3),

имеет

Седловая

величина, определяемая

выражением

при о > 0 отрицательное значение и,

следовательно,

к

петле

сепаратрисы

может стянуться или

из

нее появиться

лишь

у с т о й ч и в ы й предельный цикл. Отсюда

следует,

что

обраще­

ние в нуль величины г|>(0 ), соответствующее возникновению петли сепаратрисы, должно предшествовать исчезновению корня

h = h\. При

перемене знака ф(0)

(когда ф(0) становится

отри­

цательным)

появляется

второй

корень

h — Аг функции

ф(й),

убывающий

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

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

у с т о й ч и ­

в о м у

(г|/(Й2 ) < 0 )

предельному циклу. При дальнейшем

возра­

стании

о корни hi

и йг

сближаются,

сливаются

и исчезают.

Описанному процессу соответствует на фазовой плоскости смена

структур, представленных на рис. 205: рис. 205, 2—3

(ф(0) =

= 0)— возникла петля сепаратрисы; рис. 205,5 (AI </»2

< 0 ) —

из

петли сепаратрисы возник устойчивый предельный

цикл;

рис.

205,5—4 (hi = h2) — возник двойной полуустойчивый пре­

дельный цикл; рис. 205,4 (ф (/г)< 0 ) — циклов нет.

Приведенные простые рассуждения дают строгое доказатель­ ство существования области пространства параметров, для точек которой система (1 ) имеет по крайней мере два предельных цик­ ла. Непосредственное исследование функции г|>(/г) позволяет убедиться, что г|>(Л) не может иметь более двух корней. Это исследование, однако, довольно громоздко [67].

Г Л А В А 19

КАЧЕСТВЕННОЕ ИССЛЕДОВАНИЕ СШИТЫХ СИСТЕМ МЕТОДАМИ ТЕОРИИ БИФУРКАЦИИ

§ 1. Кусочно-линейная система с тремя параметрами. Рас­ смотрим дифференциальное уравнение [71]

х + а [1 — $F' (х) ] i + F(x) = 4 .

Полагая

F W =

"ЗГ х

при

я ____ л

 

2

^ ^ ^ 2

 

F =

~k х + 2

при

я

______ 3

:

2

2

 

и вводя новые переменные и параметры, приходим к системе вида

dx/dt = у,

 

при

я

_

 

dy/dt — — х — 2hxy +

а

2

 

^

' 2

dx/dt = у,

 

при

я .

.

3

( 1)

dy/dt = х 2h2y + а

— я

2

^

х ^

2

:

 

Будем рассматривать фазовую цилиндрическую поверхность склеенной системы, развернутую на часть плоскости, соответст­

вующей

неравенствам —я / 2 ^ х ^

< Зя/2 (рис. 206). Прямая

х = я/2

разбивает

рассматриваемую

часть

плоскости на области 1 и 2, в каж­ дой из которых фазовые траектории определяются соответственно линей­ ными системами. Прямые х = —я / 2 И 1 = Зя/2 отождествляются. Склеен­ ная система имеет два состояния равновесия: Oi(a, 0 ) и Ог(я — а, 0 ).

Точка 0{ — устойчивый фокус

при

0 < /г < 1, устойчивый узел при

1 . Точка 0 2 — всегда седло,

сепаратрисы которого определяются уравнениями

У ~ ( ^ j ± V + 1 ) [# — —- о )].

400 МЕТОДЫ ТЕОРИИ БИФУРКАЦИЙ В СШИТЫХ СИСТЕМАХ [ГЛ. 19

Обозначим через Si, 5г и S3 полупрямые х = —я /2 , х = я / 2 И 1 = Зя/2, соответствующие значениям у > 0.

Фазовые траектории системы осуществляют точечные преоб­

разования полупрямой

S i в S 2 (преобразование

L i )

и полупря­

мой S 2

в S 3 (преобразование L 2 ) . Пусть S i ,

S 2

и S 3

— ординаты

точек

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

полупрямых; т/ам,

О/002 — времена про­

бега изображающей точки (см. рис. 206) через области 1, 2, со­ ответствующие преобразованиям L \ и L 2 .

Величины S i , S 2 , т, 0 принимают положительные или нулевые значения. Интегрируя линейные уравнения (1) в областях 1 и 2, обычным образом получаем параметрические уравнения для функций соответствия.

Для области 1 получаем

Si (т, К

а) =

 

(

 

 

А

т

+

ctg( - тJ ++ hj),а ) К

S%(т, h

а) =

Si (т,

hu

d),

 

 

 

 

(Li)

Здесь cOj = V

 

i h\,

ki = hi/d)i

и 0 < h < i .

Выражения

для

производных могут быть представлены в виде

 

 

 

 

dS1

. P2hiT ll.

 

 

 

 

 

 

 

dSl

-e

5 i,

 

 

 

 

 

 

£ f i _

 

 

 

 

 

 

 

 

 

 

dS%

 

 

 

 

 

 

 

 

 

 

В случае

h >

1

выражения

для

Si (т, hi,

а)

и

Sz(T, hi,

а) и

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

жений заменим sin т и ctg т соответственно

через sh т и

cth т.

Параметр т меняется в пределах от т = 0

(5 2 (0 ) = 5 I (0) = 00)

до значения т = т*, при котором Si

или S2 обращаются в

нуль.

Тот или иной случай при 0 < h < 1

реализуется в зависимости

от знака выражения

 

 

 

/ 2 а) Л л — (я/ 2 + а).

 

(2 )

При h > 1 всегда реализуется случай 5г(т*) = 0, S i (т*)> 0 .

Кри­

вая 5 1 = 5 1 (5 2 ) имеет асимптоту 5 1 =

5 2 + 2 я/&1 . Некоторые

воз­

можные виды кривой, соответствующие преобразованию Li, изо­ бражены на рис. 207.

Рис. 207, а соответствует случаю (я/2 — а) ек*п — (я/2 + а) <С 0 й, в частности, всегда реализуется при малых hi. Промежуточ­ ный между изображенными случай, когда кривая проходит че­ рез начало координат, соответствует т* = я и обращению в нуль выражения (2 ).

Соседние файлы в папке книги