Главная страница  Приемы качественного исследования 

[0] [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] [66] [67] [68] [69] [70] [71] [72] [73] [74] [75] [76] [77] [78] [79] [80] [81] [82] [83] [84] [85] [86] [87] [88] [89] [90] [91] [92] [93] [94] [95] [96] [97] [98] [99] [100] [101] [102] [103] [104] [105] [106] [107] [108] [109] [110] [111] [112] [113] [114] [115] [116] [117] [118] [119] [120] [121] [122] [123] [124] [125] [126] [127] [128] [129] [130] [131] [132] [133] [134] [135] [136] [137] [138] [139] [140] [141] [142] [143] [144] [145] [146] [147] [148] [149] [ 150 ] [151] [152] [153] [154] [155] [156] [157] [158] [159] [160] [161] [162] [163]

2nv +

128 k* - k + i

15 fc5

64 ЗГ-А:«-2 15 fc5

2nv -

/ 128 fc* - fc + 1

15 fc5

3k-k*-2

Здесь F{n/2, k) и E{n/2, k)-полные эллиптические гнте-тралы первого и второго рода и введены обозначения

v = r/p, а = (а + Р)/р, А:2 = 2/(Л+1), К Л< оо.

Значения константы ho, выделяющие кривые С/, консервативной системы, охватывающей состояние равновесия, определяются как корни уравнения

аФ-) = .[ j {рч> + qy) dcp dy = - {а + Р cos 2ф - ysin ф) dj/ =

= - 2 /2р j (о - 2 sin2 ф /соз (р + h) dtp =

-Ро

= - 16Р {[о - (16/15) (X* - х2 + 1)] (л/2, X) -Ь [о (к2 - 1) -

- (8/15) (Зк2 - X* - 2)] F (1/2, лк)}, x2=i(/i+1)/2, -1<й<1, фо = arccos (-Л).

Корни уравнений i)i(/i) = 0, iJ)2(A) = 0, 1))з(Л) = 0 зависят от двух параметров с; и v. В плоскости (с;, v) можно получить разбиение на области, соответствующие различным возможным распределениям корней уравнений. Каждому распределению будет соответствовать определенная структура разбиения фазового пространства на траектории. Следующий набор условий (каждому условию соответствует некоторая кривая в плоскости (с;, v)) определяет все возможные в системе (9) бифуркации:

1) азз(-1) = 0, 2) a3i(oo) = 0 или а32(оо)=0,

3) з(1) = 0, 4) .pi(l) = 0, 5) al32(l) = 0,

6) i{h) = 0 и [{h) = 0,

7) i32W = 0 и я>;(/1) = 0,



Положим для определенности Р > О, v > О и будем рассматривать верхнюю полуплоскость (о, v) (при v<0 получаем разбиение пространства параметров, симметричное области v > О относительно оси о). В этом случае уравнение Tfii(ft) = 0 имеет не более одного корня, 2{h) = 0 и 1зз() = 0 -не более двух корней. Перечисленному набору условий соответствуют следую-гцие уравнения граничных кривых и бифуркаций:

1. 0 = 0 - при возрастании о из фокуса появляется неустойчивый предельный цикл.

2. 0 = 1 - при возрастании о устойчивый предельный цикл появляется из при убывании о неустойчивый предельный цикл появляется из - о°.

3. о = 16/15 = 1,066...- при возрастании о из сепаратрисы седла появляется устойчивый (так как в седле + Qy = S - р,ро <; о) предельный цикл, охватывающий состояние равновесия.

4. 2nv - 8о + 128/15 = 0 - при убывании о от петли сепаратрисы на верхнем полуцилиндре появляется устойчивый предельный цикл.

5. 2nv + 80 - 128/15 = о - при возрастании о от петли сепаратрисы на нижнем полуцилиндре появляется неустойчивый (если - Рро>0) или при убывании о от петли сепаратрисы появляется устойчивый (если -рРо<0) предельный цикл.

6. При V > о кривая не существует (в верхнем полуцилиндре не может быть двух предельных циклов при v >0).

7. Если обозначить 1з2(/г) = Р [2nv + 0Ф1 (ft) - Ф2(й)] = 0, то параметрические уравнения кривой будут

ф; (h) jh) Ф[ jh) - ф; (h) ф (h)

Ф[{к) 2яФ(Л)

Кривая проходит между точками А{0, (128/30)л~= 1,36) и 5(1, 0). При возрастании о двойной предельный цикл, возникший на нижнем полуцилиндре из сгущения траекторий, разделяется на два (нижний-неустойчивый, верхний - устойчивый) .

8. Из il:3(ft) = -16p[T,(ft)o-T2(ft)] = 0 и pз{h) = 0, исключая о, получаем уравнение для определения ft. Уравнение 21 - 21 = О имеет единственный корень ft = 0,86, соответствующий 0 = 1,09. При возрастании о исчезает двойной предельный цикл, охватывающий состояние равновесия.

Разбиение пространства параметров о, v на области, которым соответствуют различные качественные структуры разбиения фазового пространства, представлено на рис. 244. Штриховкой отмечены две тонкие области, для точек которых в фазовом пространстве есть два предельных цикла.



Исследование системы (9) при полигональных и разрывных характеристиках может быть аналогично проведено методам малого параметра (см. гл. 18) [66].

Система (9) с полигональными характеристиками (см. рис. 243, б) и малым параметром р имеет при р = О интеграл

(- (ф + пум - л < Ф < - л/2,

ф2/я-я/2 =ft, -(ф-я)7я

Н{ср,у) = \ +

- я/2<ф<я/2,

я/2<: ф<:я.

(12)

-\ 2

0,5 В

1,5 в

Рис. 244

Замкнутые кривые семейства (12) при -я/2</г<0 охватывают особую точку, при О < /i < о° - фазовый цилиндр. Состояния равновесия при р=?0 смещены с линии сшивания и будут 0,((1/2)рГя, 0)-фокус и 02(я -(1/2)рГя, 0) -седло.

Приведем (в том же порядке, что и для (10)) уравнения границ в пространстве параметров о и v:

1) 0 = 0, 2) 0=1,

3) о-(4У2/3)(1 + я/4)- = = 1,056...,

4) 2яу-(У2/2)яЗ/2(Ц-+ я/4)о + (4/3)яЗ/ = 0,

5) 2яу+>(У2/2)яЗ/2(1 + + я/4)о-(4/3)яз/2 о,

6) при V > о не существует,

7) кривая проходит между точками .4(0, (2/3) Уя = и 5(1, 0),

8) 0 = 1,07...

Разбиения фазового пространства и пространства параметров для полигональных характеристик останутся качественно тождественными разбиениям для характеристик рис. 243, а. Сохранятся и тонкие области, для точек которых в фазовом пространстве есть два предельных цикла. Их размеры лишь незначительно изменятся. На рис. 245, а, д представлены разбиения цилиндрического фазового пространства соответственно для областей i и 2 рис. 244.

Для полигональных характеристик (см. рис. 243, б) качественно эквивалентные рис. 245, а, д разбиения фазового пространства будут соответствовать областям пространства параметров, расположенных, как и на рис. 244, в полосе 0<а<1. Для обеих рассмотреппых аппроксимаций уравнения (9) области i и 2 в пространстве параметров будут разделены узкой полосой, для точек которой в фазовом пространстве есть два предельных

1,26...>




[0] [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] [66] [67] [68] [69] [70] [71] [72] [73] [74] [75] [76] [77] [78] [79] [80] [81] [82] [83] [84] [85] [86] [87] [88] [89] [90] [91] [92] [93] [94] [95] [96] [97] [98] [99] [100] [101] [102] [103] [104] [105] [106] [107] [108] [109] [110] [111] [112] [113] [114] [115] [116] [117] [118] [119] [120] [121] [122] [123] [124] [125] [126] [127] [128] [129] [130] [131] [132] [133] [134] [135] [136] [137] [138] [139] [140] [141] [142] [143] [144] [145] [146] [147] [148] [149] [ 150 ] [151] [152] [153] [154] [155] [156] [157] [158] [159] [160] [161] [162] [163]

0.0111