РАДИОТЕХНИКА И СВЯЗЬ
УДК 621.3 (075)
ОТОБРАЖЕНИЕ НА ФАЗОВОЙ ПЛОСКОСТИ МЕДЛЕННЫХ ПРОЦЕССОВ В КОНСЕРВАТИВНЫХ ЦЕПЯХ С ОДНОЙ СТЕПЕНЬЮ СВОБОДЫ ПРИ НЕЛИНЕЙНОМ РЕЗОНАНСЕ
В.Ф. Судаков
МГТУ им. Н.Э. Баумана, Москва, Российская Федерация e-mail: [email protected]
Рассмотрены слабо нелинейные электрические цепи без резисторов с одним источником гармонического возбуждения. Дифференциальные уравнения цепей идентичны типовому нелинейному уравнению второго порядка. Амплитуда возбуждения и его частотная расстройка считаются малыми. Методом усреднения получены уравнения для амплитуды и фазы резонансного колебания. Для анализа плоскости амплитуда-фаза предложен простой и эффективный метод, применимый только в случае консервативных цепей. С его использованием построена амплитудно-частотная характеристика и выделены режимы одно-и бистабильности. Каждому режиму соответствует фазовый портрет: в режиме одностабильности он имеет одну особую точку, в режиме бистабильности — три. Определены координаты особых точек, характер их устойчивости, фазовые траектории в их окрестности. В режиме бистабильности качественно построен полный фазовый портрет, позволяющий представить переходные процессы при различных начальных условиях.
Ключевые слова: нелинейный резонанс, бистабильность, фазовая плоскость, фазовый портрет, особые точки на фазовой плоскости.
DISPLAY ON THE PHASE PLANE OF THE SLOW PROCESSES IN CONSERVATIVE CHAINS WITH ONE DEGREE OF FREEDOM AT NONLINEAR RESONANCE
V.F. Soudakov
Bauman Moscow State Technical University, Moscow, Russian Federation. e-mail: [email protected]
The weakly nonlinear electric chains without resistors with a single source of harmonic excitation are discussed. Differential equations of the chains are identical to the typical nonlinear equation of second order. The amplitude of the excitation and its frequency detuning are small. The equations for the amplitude and phase of resonant oscillation are obtained by the method of averaging. A simple and effective method (applicable only in the case of conservative chains) is proposed for the analysis of amplitude-phase plane. Using it amplitude-frequency characteristic is presented and modes of one- and bistability are selected. The phase portrait corresponds to each mode: in the mode of one-stability it has one particular point, in the mode of bistability - three points. The coordinates of special points, the nature of their stability and phase trajectory in their surroundings are determined. In the mode of bistability the full phase portrait is built qualitatively and it enables to present transient processes at various initial conditions.
Keywords: nonlinear resonance, bistability, phase plane, phase portrait, special point on the phase plane.
Консервативные системы специфичны в том отношении, что не предполагают наличия диссипативных сил (не имеют энергетических потерь) и могут испытывать только воздействия потенциального типа (к тому же не зависящие явно от времени). Область приложения теории точно консервативных (гамильтоновых) систем невелика: движение заряженных частиц в постоянных электрическом и магнитном полях, движение небесных тел (в том числе космических аппаратов). В большинстве технических приложений модель консервативных систем в чистом виде не работает хотя бы потому, что диссипация не только всегда есть, но зачастую она играет принципиальную роль. Тем не менее, консервативные системы изучаются давно и очень активно. Это объясняется рядом обстоятельств, в частности тем, что потери обычно малы и могут быть учтены в первом порядке теории возмущений. В этом случае нулевое приближение - консервативная система. Другая причина интереса к консервативным системам следующая: именно в них отчетливо проявляется ряд нелинейных эффектов, например нелинейный резонанс. При учете малых потерь резонанс этого типа тоже имеет место, но его основные особенности все-таки наследованы от нелинейного резонанса в приближенной (модельной) консервативной системе.
Нелинейный резонанс — наиболее устойчивое состояние движения в природе. Резонансные механизмы, устройства и резонансные технологии имеют КПД ~ 100 % и являются наиболее оптимальными и перспективными в физике, химии, биологии и медицине, для техники и производства. Нелинейный резонанс позволяет достаточно просто объяснить многочисленные аномальные явления: ускорение биохимических реакций, трансмутацию элементов, холодный термоядерный синтез, бесконтактную активацию жидкостей, запись информации на воде, БАД, вихревые генераторы теплоты и энергии, сверхновые торсионные передатчики информации, уникальные слабоминерализованные стерилизующие, дезинфицирующие, моющие растворы, ан-тиоксидантные свойства и резонансную микрокластерную структуру водных растворов с временем релаксации от секунд до нескольких суток и более, и т.д.) и наметить дальнейшие пути решения назревших проблем методами нелинейной классической физики.
В электронике нелинейный резонанс может быть получен в очень простой цепи с нелинейной индуктивностью (феррорезонас) или с нелинейной емкостью, в качестве которой часто используется барьерная емкость запертого р-п-перехода полупроводникового диода. При контакте р- и п- областей часть электронов из п-области переходит в область р, а часть дырок из р-области переходит в п-область. Вблизи границы п-области образуется некомпенсированный положительный заряд неподвижных ионизированных доноров. Следовательно, в р-
области из-за ухода дырок образуется отрицательный объемный заряд неподвижных ионизированных акцепторов. Между этими зарядами возникает электрическое поле, которое препятствует дальнейшему переходу электронов в область р и дырок в область п. Разность потенциалов между областями р и п, образующаяся при контакте, равна разности работ выхода электронов. При подаче запирающего напряжения высота потенциального барьера между областями р и п возрастает. Это приводит к расширению р-п-перехода, вследствие чего увеличивается положительный объемный заряд в области п и отрицательный в области р. При подаче на р-п-переход прямого напряжения высота потенциального барьера уменьшается, переход сужается, заряды уменьшаются. Таким образом, изменения напряжения, приложенного к переходу, приводят к изменению объемного заряда в переходе, т.е. переход действует как емкость. Значение этой емкости зависит от напряжения на р-п-переходе, т.е. она нелинейная.
В машиностроении используются упругие системы с нелинейным элементом (обобщенно говоря, "пружиной"), для которого связь между деформацией и упругой силой нелинейная, т.е. нарушается закон Гука. В таких системах также возможен нелинейный резонанс при внешнем воздействии (как полезный, так и вредный).
В технике передачи информации альтернативой цифровым методам обработки сигналов является функциональная электроника. Ее операции обычно включают свертку, корреляцию, преобразование Фурье, согласованную фильтрацию и задержку сигналов. Эти операции выполняют твердотельные сигнальные процессоры, основанные на нелинейном резонансном взаимодействии электромагнитного поля с веществом. Выходной сигнал этих устройств представляет собой спиновое, фононное или фотонное эхо. Такие процессоры могут использоваться в качестве управляемых линий задержки и управляемых согласованных фильтров, осуществлять корреляционный и спектральный анализ сигналов.
Приведенный краткий и далеко не полный перечень применений нелинейного резонанса в консервативных системах позволяет оценить важность изучения этого явления. Более того даже простейшие нелинейные системы в условиях резонанса приобретают принципиально новое качество — в них возможна хаотическая динамика. Это, например, путь получения шумоподобных (т.е. сверхзащищенных) сигналов или систем с непредсказуемым (так называемым хаотическим) поведением.
В настоящей работе акцент сделан на консервативной электрической цепи. Ясно, что это только вопрос выбора технической модели. Существо же дела в используемой математической модели, общей для любых систем упомянутого типа.
Явления, возникающие в нелинейных электрических цепях при гармоническом воздействии, имеют большое прикладное и теоретическое значение. Поэтому им уделялось пристальное внимание [1, 2]. В таких цепях не всегда возникает резонанс, но если он есть, то его свойства отличаются от свойств резонанса в линейной цепи. При нелинейном резонансе колебания на всех частотах внутри резонансной области имеют конечные амплитуды. Нелинейный резонанс может возникать и в тех случаях, когда частота внешнего воздействия близка к кратному собственной частоты или ее дробной части. При нелинейном резонансе всегда имеет место гистерезис как в зависимости от расстройки, так и в зависимости от амплитуды воздействия. Исследованию нелинейного резонанса способствует удачный выбор применяемых для этого математических методов. Как правило, это методы усреднения [3-5]. Применительно к теории нелинейного резонанса методом усреднения всегда рассматривается наиболее общий случай нелинейной цепи с потерями. При такой общности постановки задачи и сложности исследования особенности полного фазового портрета цепи в условиях резонанса обычно не рассматриваются. Но если ограничиться только консервативными цепями, то можно построить их фазовый портрет на плоскости медленных параметров (амплитуды и фазы). Фазовый портрет позволяет представить переходные процессы для этих параметров при различных начальных условиях. В известных нам публикациях переходные процессы исследовались исключительно на аналоговых моделях или с использованием ЭВМ. Цель настоящей работы — приближенное изучение переходных процессов на фазовой плоскости аналитическими средствами.
Чтобы максимально упростить определение особых элементов фазовой плоскости, предложен метод исследования укороченных уравнений, полученных после применения метода усреднения. Этот метод основан исключительно на рассмотрении консервативных цепей и отличается простотой и эффективностью.
Цепи и их модельные уравнения. Рассмотрим одноконтурную цепь с гармоническим источником напряжения без резисторов.
dф
По законам Кирхгофа получим два уравнения: e (t) = — + uc,
C^UC = i где e (t) — ЭДС источника; ф — потокосцепление катушки; uc — напряжение на конденсаторе; i (ф) = L-1 (ф + аф3) — вебер-амперная характеристика катушки. Эти уравнения эквивалентны одному уравнению второго порядка:
d2ф
-ф = -C-1L-1 (ф + аф3) + uE0 cos ut.
Рассмотрим двухузловую цепь без резисторов с гармоническим источником тока. По законам Кирхгофа получим два уравнения
J (t) = C^^ + i (p), = , где J (t) — ЭДС источника. Эти dt dt
уравнения эквивалентны другому уравнению второго порядка:
= -C-1L-1 (Р3 + ap) + C-1J (t).
Оба уравнения второго порядка идентичны типовому уравнению
виДа d2x
— + = f (x,wt). (l)
К выражению (l) приводятся уравнения и всех других консервативных слабо нелинейных цепей с гармоническим воздействием. Назовем поэтому (l) уравнением цепи.
Укороченные уравнения. Специфика уравнения (1) в том, что ц — малый параметр, а нелинейная зависимость от х и зависимость от времени (внешнее воздействие) разделяются на f (x, wt) = k (x) + + H cos wt. Порождающее уравнение (при ^ = 0) имеет собственную частоту w1 и описывает автономную цепь без потерь. Рассмотрим глав-
Wi 1
ный резонанс, когда — = (1 + ^п)2. Немалая величина п называется
w
относительной расстройкой. Главный резонанс изучался многими авторами и описан во многих публикациях (например [6] и др.). Известно следующее: ни при какой расстройке п амплитуда вынужденных колебаний не стремится к бесконечности, что объясняется анизохронизмом этих колебаний; вынужденные колебания ангармоничны; в диапазоне расстроек пгр < П < ^ возможен бистабильный режим вынужденных колебаний, т.е. существуют два типа устойчивых колебаний, а возможное колебание третьего типа неустойчиво. Далее вернемся к этому вопросу, изменив известный подход и получив тем самым возможность продвинуться в изучении главного резонанса.
Перейдем к безразмерному времени ^ = w1t и параметру v =
_ i
= (1 + ^п) 2 в исходном уравнении (1). Тогда
d2 x
— + x = ^w-2f (x . (2)
Поскольку автономная цепь (с одной степенью свободы) находится под гармоническим воздействием, результирующая цепь имеет 3/2 степени свободы и должна рассматриваться в трехмерном фазовом пространстве. На фазовой плоскости можно изучать только медленно меняющиеся переменные (амплитуду и фазу), уравнения для которых могут быть получены методом усреднения [4].
Представим решение уравнения (2) в виде x = y cos (z + <^), где амплитуда y (t) и фаза z (t) есть функции времени со скоростями изменения порядка Из (2) известным образом [4] получим систему уравнений первого порядка (стандартного типа, по терминологии
Н.Н. Боголюбова):
dy . / ч dz g
— = ^g sin (z + , — = cos (z + , (3)
a^ a^ y
где g (y, z, = ny cos (z + — ш-2/ (y cos (z + ,
После усреднения по ^ на интервале [0, 2п] система (3) переходит в укороченную автономную систему вида
I=(y,z), I=(y,z)=. (2+f), (4)
где 2n
Y (y,z) = 2П/ g (У, ^ ^ sin (z + ^ =
0 2n
=--22 / (у cos (z + , sin (z + d<£,
2n
1
0
2n
Z (y, z) =--[k (y cos (z + + H cos cos (z + d<£.
2n J 0
Учтем, что к (ж) = ж3, где 7 < 0. Интегралы У (у, и . (у, могут быть вычислены явно, в результате чего получим
dy H . dz / n 3 2 H —— = — u—„ sin z, —— = u —|— Yy--
Эти уравнения известны. Они позволяют найти как свойства стационарных колебаний при резонансе (особые точки системы (5)), так и переходные процессы при различных начальных условиях на фазовой плоскости (y,z).
Вместо параметров H, ^1,7, n удобнее использовать два обобщенных параметра
Р = 37^ П < 0 при n < 0, q = 37^ H > 0 (6)
и приведенную "быструю" фазу s = 37^ <£.
С учетом этих обозначений укороченные уравнения (5) записываются следующим образом:
dy q . dz 1 2 q
= 2 Те = 2^Р + у - У с°8 (7)
Параметры стационарных колебаний на фазовой плоскости
{у, г}. Они определяются как координаты особых точек на фазовой
плоскости (у, г). Предлагаемый далее метод исследования, пригодный только в случае консервативных цепей, отличается от общего традиционного подхода (используемого во всех известных нам литературных источниках и изложенного, например, в [6, 7]). В указанном случае предлагаемый метод достаточно эффективен.
Заметим, что правая часть первого уравнения системы (7) обращается в нуль при г = 0, п. Тогда второе уравнение этой системы обращается в нуль при значениях у таких, что
Р + У2 - - = 0, г = 0, (8)
у
р + у2 + - = 0, г = п. (9)
у
Изобразим графически (рис. 1) зависимости р (у), которые следуют из (8), (9):
Я. я.
Р0 (У) = -У2 + -, Рп (У) = -У2 --. (10)
У У
В совокупности зависимости (8), (9) неявно определяют амплитудно-частотную характеристику (АЧХ) цепи. Из приведенного графика следует, что при одних значениях расстройки существует только один корень (по у) системы (8), (9), а при других — три корня. Другими словами, в зависимости от значения частотной расстройки могут существовать одно или три стационарных колебания. Покажем аналитически, что (9) может иметь два или ни одного корня, а уравнение (8) всегда имеет один корень.
Представим (8) в виде неполного кубического уравнения
у3 + РУ - — = 0, (11)
корни которого определяются аналитически [8]. Их вид зависит от
|р|43
с = ^ = *_Ы!_> о
Z /ТЧ 2 Q3 /о \ / тт\ 2 >
■П2 3V 3 WH
i'7U
—— _
Пусть 0 < С < 1. Рассмотрим = — \/1 — С < 0. Тогда уравне-
2
ние (11) имеет только один положительный корень [8]:
у., = (2 + 3 + (2 — 1 =
1
(i - уг-с < °)3 + (1 + /Л-С < о)
(12)
3
Пусть Z > 1. Рассмотрим угол а в качестве более удобного параметра:
cos а =- 4 = Z- 2 < 1. (13)
Тогда уравнение (11) имеет три вещественных корня [8]:
^ /Ы а /Ы (а п
yi,0 = ^ycos3' У2,0 = + 3
^ /Ы (а п У3'0 = -2V У COV3 - 3
(14)
Проанализируем эти корни, используя (14). Из (13) следует, что
п а п
0 < а < — и поэтому 0 < — < —, т.е. у1 о > 0. С другой сторо-
2 3 6'
папп папп
ны, 3 < з + з < 2 и поэтому у2 0 < 0, - 3 < 3 - 3 < - 6 и поэтому у3 0 < 0. Таким образом, только один корень у^о > 0, два
других отрицательны и физически недопустимы.
Итак, при всех значениях параметра 0 < С < на фазовой плоскости {у, г} существует единственная (физически допустимая) стационарная точка {у10, 0}, но значения стационарной амплитуды у10 в зависимости от значения £ следует вычислять либо по (12), либо по (14).
Представим (9) в виде неполного кубического уравнения,
у3 + ру + д = 0, (15)
корни которого также определяются аналитически. Их вид зависит от величины £ (при р < 0). Как следует из графика (см. рис. 1), вещественные корни (15) могут существовать только при достаточно больших по модулю расстройках |р| = |п|. Найдем граничный параметр £гр, соответствующий граничному значению расстройки.
По (10) (у) = 2у - у2 и равна нулю при уф = 3. Тогда
2
Рп (угр) = — угр--— = —3 ( 3. Корни (15) могут существовать
угр 2
только при |р| > |рп (угр)|. Параметр ^гр, который соответствует точке на графике [угр,рп (угр)], находится по формуле
>п (угр)Г 3
Z = у 3 / =1
Zrp / q \ 2 1.
2
Соответственно, корни (15) существуют только при £ > £гр = 1. Это обстоятельство позволяет найти указанные корни аналитически.
При С > Сгр = 1, Р < 0 корни уравнения (15) находятся по формулам [8]
„ /Ы а /Ы (а п
yi,o = ^yCOS3, У2,0 = + 3
^ /Ы (а п
y3,0 = у co^3 - 3
где
— 1 < cos а =
q
d 3
= -Z-2 < о.
(16)
(17)
тг .... пап
Из (17) следует, что п < а < —. Поэтому — < — < — и, как следует
2 3 3 2
2п а п 5п а п п
из (16), > 0. Кроме того, — < - + - < — и 0 < - - - < -, что
3 3 3 6 3 3 6
позволяет легко установить из (16): у2,п > 0, у3 п < 0. Отрицательная
амплитуда у3 п < 0 физически недопустима.
Итак, при допустимых значениях параметра £ > 1 на фазовой плоскости {у, г} всегда существуют только две стационарные точки {У\п, п} , {у2,п, п}. Значения стационарных амплитуд у1п,у2 п следует вычислять по (16).
Полученные расчетные формулы позволяют без труда выполнить расчет (например, с использованием пакета МайСАБ). Сведем в таблицу значения стационарных амплитуд для различных £ (зависит от расстройки) при д = 2.
Значения стационарных амплитуд
Z 0,1 0,5 0,7 0,8 0,9 1 1,5 2 10 25 150 300
У1,0 1,621 1,859 1,924 1,952 1,977 2 2,095 2,126 2,685 3,07 4,054 4,53
У2,п 1 1,425 1,587 2,371 2,84 3,928 4,431
У1 ,п 1 0,67 0,581 0,314 0,229 0,126 0,1
При г = 0 с ростом модуля £ значение резонансной амплитуды у10 возрастает. При г = ж с ростом £ (при £ > 1) возрастает стационарная амплитуда у1п и убывает стационарная амплитуда у2 п. Нетрудно убедиться в этих закономерностях и по графику рис. 1, если
учесть, что параметр £ связан с параметром р следующим образом: 1 2
|р| = 3£3 3. Роль найденных корней (стационарных амплитуд резонансных колебаний) в динамике цепей различна: различен характер
устойчивости соответствующих особых точек на фазовой плоскости {y, z}.
Тип устойчивости стационарных колебаний на фазовой плоскости
{y, z}. Обратимся к уравнениям фазовых траекторий (4). Какова бы ни была нелинейность реактивного элемента консервативной цепи (не только кубической, как в рассматриваемой задаче), уравнения (4) имеют специфическую зависимость правых частей от y, z (Y0, Z — постоянные):
Y (y, z) = Yo sin z, Z (y, z) =
= n + Zo (y) + — cos z. 2 y
Отсюда следует, что
Рис. 1. Зависимость частотной расстройки от амплитуды колебаний в форме р (у) по формулам (10)
dY dy
= 0,
q
dZ _ 1/
dy 2 I У2
cos z
dY 1
= - 2 q COS Z
dZ 1 q . — = --sin z. dz 2 y
Стационарным точкам соответствует sin z = 0 (z = 0,п) и, сле-dY dZ dZ 1 / a\ dY 1 !
довательно, sy = 0, A* = 0, ly = 2 (2y ± V4, "dZ = T 2q.
Характер устойчивости стационарных точек определяется коэффици-
dY dZ A dZ dY
ентами а = —--h -г— = 0 и Д = — характеристического
dy dz dy dz
уравнения A2 + ¿A + Д = 0. В случае кубической нелинейности, как следует из (7),
а=±1 ч 2y ± £
(18)
Корни зависят от знака А: при А < 0 корни вещественны, при А > 0 корни чисто мнимые. Следовательно, стационарная точка, для которой А < 0, есть седло, а стационарная точка, для которой А > 0, есть центр. Исследуем стационарные точки {у10, 0} , {у1>п, п} , {у1>п, п}, найденные ранее, на устойчивость, определив для них знак параметра А.
Рассмотрим стационарную точку {у1?0, 0} , 0 < С < го. Тогда
9 3
А = °2,0 = 7-V (2y3,0 + q) . 4y1,0
(19)
1 Верхний знак соответствует z = 0, нижний — z = п.
д
Ранее было получено, что —р = у20 +--> 0. Следовательно,
' У1,о
> 0, т.е. точка {у1>0,0} во всем диапазоне рас-
q
Д = ^ = — Р 4У2,0 строек |п| = С — это центр.
Рассмотрим стационарную точку {у1>ппри £ > Сгр = 1 .
Тогда
д (п..3
А =
4y f
(2y3 - q) .
(20)
a
Найдем 2y3,n — q для = 2^/ — cos —. Поскольку в силу (17) имеем
3 3
'И
—q = 2\1 ( — | cos a, то
2y3,n — q = 2 ^2
После некоторых преобразований получим
—3cos I1 +2
cos a.
2»L — q = ЦЦ 12cos 1 (cos 1 + 5) (cos 1 — 2
(21)
Для этой стационарной точки, как было показано при обсуждении
ч п a п _ a a 1 a 1
(17), — < — < —. Отсюда cos — > 0, cos —I— > 0, cos---< 0.
3 3 2 3 32 32
nan
Таким образом, из условия — < — < — существования стационарной
332 точки {y1,n,п} следует, что в силу (20), (21)
А = =
q
4yL
(2y3,n — q) > 0,
(22)
т.е. эта точка есть центр.
Рассмотрим стационарную точку {у2,ппри С > Сгр = 1.
I р (а
Здесь у2п = — 2л1 — - + — у. Поскольку в силу (17) имеем
Ъ|
—q = 2\1 ( — | cos a, то
p
2y3,n — q = 2 —2 .'-cos - + -
a
п
+ 2i Иг
3 V3 3
Отсюда после некоторых преобразований получим
cos a.
2y3,n — q = 2
cos
a
a п
+ 77
3 3
1"
з) — 2
X
a п\ 1 — co^ 3 + 3) — 2
(23)
3
3
3
Для этой стационарной точки, как было показано при обсуждении
пап 2п а п 5п (17), можно взять — < — < — или —- < — + — < —-. Следовательно,
3 3 2 3 3 3 6 /а п\ /а п\ 1 /а п\ 1 „ ^
с°э - + - < 0, — с°э — + - — > 0, с°э — + - — < 0. Таким \3 3) \3 3) 2 \3 3) 2
пап
образом, из условия — < — < — существования стационарной точки
332 {у2,п, п} следует в силу (21), (23), что
А = — А2,п = — ту^ (2у3,п — д) < 0, (24)
т.е. эта точка — седло.
Отметим важную особенность исследования устойчивости стационарных режимов в случае резонанса для нелинейных консервативных цепей с одной степенью свободы: какой бы ни была нелинейность реактивного элемента характер устойчивости определяется зна-2 (у)
ком у2—--± ¿1 при стационарном значении амплитуды у. В част-
«у
ности, если нелинейность кубическая, для исследования устойчивости следует оценить только знак величины 2у3 ± д (при стационарном значении амплитуды). Фактически ранее это и было проделано.
Фазовый портрет цепи на плоскости {у, г}. Фазовый портрет рассматриваемых цепей на плоскости {у, г} можно в общих чертах
описать, исходя из общих соображений [9]. Пусть параметры цепи
2
(|р| > |ргр| = 3 ^3 ^ таковы, что на плоскости {у, г}, могут существовать три стационарные точки. Если начальные значения амплитуды и фазы колебаний цепи совпадают с координатами стационарных точек, то эти амплитуды и фазы могли бы оставаться такими же и при £ ^ то. Реально "точные" начальные условия не существуют, неизбежны хотя бы незначительные отклонения. В результате фазовая траектория покинет сколь угодно малую окрестность неустойчивой стационарной точки, однако останется в окрестности устойчивых стационарных точек. Из общей теории следует, что любые начальные условия, отличные от координат устойчивой стационарной точки типа центр, вызовут периодические изменения как амплитуды, так и фазы стационарных колебаний в цепи. Изменения такого типа следует рассматривать как незатухающий (в цепи нет резисторов) переходный процесс. Через неустойчивую особую точку типа седла {у2,п,п} проходят только две замкнутые фазовые траектории, сепаратрисы. Они не реализуемы, так как период движения по ним фазовой точки бесконечен. Однако две замкнутые сепаратрисы играют важную роль как границы двух областей притяжения в окрестности двух устойчивых стационарных точек. Любую область притяжения фазовая траектория не может покинуть при £ ^ то. Области притяжения каждого из цен-
тров плотно покрыты замкнутыми фазовыми траекториями. Параметры колебаний существенно различны в зависимости от того, в какую из областей притяжения попадут начальные условия.
В данной цепи возможен и другой случай, когда есть только одна стационарная точка типа центр. Ее областью притяжения является вся фазовая плоскость, т.е. при любых начальных условиях возникает периодический переходной процесс. Множество замкнутых траекторий плотно покрывает всю фазовую плоскость.
Рассмотрим малую окрестность центра {у^0,0}: у = у1>0 + £у1>0, г = Уравнения малых колебаний нетрудно получить как уравнения в вариациях:
d£yi,o qr m-1„2 г
= --OZi,o, j- = ö ^io^yi,o.
^ 2 ^ ^ \2
На фазовой плоскости (в малой окрестности центра) эти уравнения представляют семейство (при различных С20) фазовых траекторий в форме эллипсов
с^Й С20 (Т^-1
Отношение полуосей этих эллипсов равно
o + -1=1. (25)
\
q^-2 / 4q2 / 16qy2,o
-i v^?,o V2y?,o + q
4
> 1,
где использовано (19). Следовательно, эллипсы растянуты в направлении оси у. Временные зависимости £у10 (£), (£) представляют
собой гармонические колебания с частотой П10 = ( - -7 ) П1>0.
4
Рассмотрим малую окрестность центра {у1,п,п}: у = у1п+,
г = п + . Уравнения малых колебаний также получим как уравнения в вариациях:
= 2 ^, = 2) ^^•
В малой окрестности центра на фазовой плоскости эти уравнения представляют семейство (при различных С2П) фазовых траекторий также в форме эллипсов:
i
,2 ^ ^
-i = 1. (26)
С2,п (2)-i §)п-:
Отношение полуосей этих эллипсов составляет
\
^ П-2 2)
У*,*
-1~\1 qo„.3 ' „ < 1
2yi,n- q
2 ' '
где использовано (22), следовательно, эллипсы растянуты в направлении оси г. Временные зависимости (£), (£) представляют
собой гармонические колебания с частотой П 1п = 3.
Рассмотрим малую окрестность седла {у2,п,п}: у = у2,п + ^у2,п,
г = п + £г2,п. Уравнения малых колебаний запишем как уравнения в вариациях
^У2,п Ях / д \-1л2 г
=2 ^, =Ы Д2>п ^.
На фазовой плоскости (в малой окрестности седла) эти уравнения представляют семейство гипербол при разных :
¿yL ¿4
2'П
CL q C22 л -1Л2 .
= 1. (27)
2 С2 ,п У2) Л2,п
Асимптоты этих гипербол (сепаратриссы седла) представляют собой в малой окрестности седла прямые. Угловой коэффициент такой прямой
$у2,п = ,п в квадранте £у2 ,п > 0,£г2 ,п > 0 равен
Y =
\
q/2
d-1Л2;
Если использовать (24), то получим
7 < 1- (28)
Следовательно, гиперболы с осью в направлении оси у заключены в угле меньшем п/2. Изменение во времени вдоль фазовой траектории определяется апериодическими зависимостями
¿У2,п (*) = А2,пеА2'п(-37Ь4 + В2,пе-А2'п(-37Ь4,
( \-1 (*) (29)
^ (') = и —¿г- .
Постоянные А2 п, В2 п определяются начальными условиями. В частности, можно получить закон движения по сепаратриссе и убедиться, что скорость движения по устойчивым усам сепаратрисы стре-
0,06
-0,015 -0,010 -0,005 0 0,005 0,010 0,015
-0,06
Рис.2. Сепаратрисы седла {у2,п} для £ = 300, д = 2; полуось 8г2,п = 0, 8у2,п < 0 направлена на центр {у\,ж}, полуось 8г2,п < 0, 8у2,п = 0 направлена на центр {ух,о, 0}
миться к нулю по мере приближения к седлу. На рис. 2 показаны фазовые траектории в малой окрестности седла {у2,п, п}.
Использованы параметрические зависимости (9), хотя тот же результат был получен и с помощью формулы (27) для фазовой траектории. Вершины гипербол, приведенных на рис. 2, настолько близки к нулевой точке (к седлу), что дают представление о направлении сепа-ратрисс.
Анализ положения и устойчивости стационарных точек позволяет качественно построить фазовый портрет цепи на плоскости медленных параметров {у, г} (рис. 3).
В силу периодической зависимости от фазы г фазовая плоскость представляет собой плоскую развертку цилиндра —п < г < п,у > 0. Кроме того, все фазовые траектории симметричны относительно оси г = 0, что следует из исходных дифференциальных уравнений. Эти
z
&2>
У
{У1,7С> -Л)
{У2,я> -Л}
Рис. 3. Фазовый портрет цепи на плоскости {у, г}
особенности фазовой плоскости отражены в фазовом портрете, в частности, в изображении замкнутых траекторий: они показаны разомкнутыми с краями на разрезе фазового цилиндра.
В малых окрестностях точек типа центр амплитуда и фаза совершают гармонические колебания с частотами, зависящими от положения центра и независящими от начальных условий. С удалением начальных значений от центров колебания остаются периодическими (фазовые траектории замкнуты), но их частоты уменьшаются с увеличением площади, ограниченной замкнутыми траекториями. Форма замкнутых траекторий с удалением начальных значений от центров все сильнее отличается от эллипсовидной, что говорит об ангармоничности колебаний. Из седловой точки выходят две фазовые траектории (устойчивая и неустойчивая ветви замкнутой сепаратрисы1), разделяющие фазовую плоскость на две области притяжения. Траектории, начинающиеся во внутренней области, будут всегда в ней оставаться. То же относится и к внешней области, в которой все траектории колеблются вокруг центра {у1>п, п}. Только изменением начальных условий можно объяснить переход из области притяжения одного центра в область притяжения другого.
Заключение. Показано, что динамические уравнения для двух видов (а также других этого типа) консервативных цепей при одинаковом виде нелинейности идентичны, что позволило изучать одно нелинейное дифференциальное уравнение общего вида.
Применение предложенного метода исследования укороченных уравнений для амплитуды и фазы позволило на фазовой плоскости этих параметров указать координаты особых точек. Установлено, что приведенный параметр £ цепи имеет бифуркационное значение: когда параметры цепи таковы, что £ < 1, на фазовой плоскости только одно устойчивое состояние — центр; если же £ > 1 (за счет увеличения частотной расстройки или/и амплитуды внешнего воздействия), то на фазовой плоскости три особых состояния (два центра и седло). В последнем случае амплитудные координаты особых точек чередуются: между двумя центрами находится седло. Такое расположение особых точек равносильно появлению гистерезиса.
Фазовая плоскость эквивалентна полосе —п < г < п на цилиндрической поверхности. В указанной полосе качественно построен фазовый портрет для наиболее интересного случая £ > 1, когда имеются три особых точки. Сепаратрисы седла отделяют области притяжения каждого из центров: переходные процессы, начавшиеся в одной из них, останутся в ней навсегда. Найдены формы фазовых траекторий (полуоси эллипсов) в окрестности центров и угловые скорости (частоты) на них; определены угловые коэффициенты сепаратрис в
1 Сепаратрисы показаны жирными линиями.
окрестности седла. Для всех величин указана связь с исходными параметрами цепи. Характер фазового портрета позволяет сделать вывод: резонансное колебание всегда имеет амплитудно-фазовую модуляцию с непредсказуемыми (но измеримыми) параметрами.
ЛИТЕРАТУРА
1. Бессонов Л.А. Теоретические основы электротехники. Электрические цепи. М.: Наука, 2007. 701 с.
2. Данилов Л.В., Матханов П.Н., Филиппов Е.С. Теория нелинейных электрических цепей. М.: Энергоиздат, 1990. 256 с.
3. Розо М. Нелинейные колебания и теория устойчивости. М.: Наука, 1971. 228 с.
4. Боголюбов Н.Н., Митропольский Ю.А. Асимптотические методы в теории нелинейных колебаний. М.: Наука, 1974. 504 с.
5. Моисеев Н.Н.Асимптотические методы нелинейной механики. М.: Наука, 1969. 381 с.
6. Бутенин Н.В., Неймарк Ю.И., Фуфаев Н.Л. Введение в теорию нелинейных колебаний. М.: Наука, 1987. 382 с.
7. Ландау Л.Д., Лифшиц Е.М. Механика. М.: Наука, 1988. 214 с.
8. Корн Г., Корн Т. Справочник по математике для научных работников и инженеров. М.: Наука, 1984. 831 с.
9. Андронов А.А., Витт А.А., Хайкин С.Э. Теория колебаний. М.: Наука, 1981. 918 с.
REFERENCES
[1] Bessonov L.A. Teoreticheskie osnovy elektrotekhniki. Elektricheskie tsepi [Theoretical basis of Electrical Engineering]. Moscow, Nauka Publ., 2007. 701 p.
[2] Danilov L.V., Matkhanov P.N. Filippov E.S. Teoriya nelineynykh elektricheskikh tsepey [The theory of nonlinear electrical circuits]. Moscow, Energoizdat Publ., 1990. 256 p.
[3] Rozo M. Nelineynye kolebaniya i teoriya ustoychivosti [Nonlinear vibration and stability theory]. Moscow, Nauka Publ., 1971. 228 p.
[4] Bogolyubov N.N., Mitropol'skiy Yu.A. Asimptoticheskie metody v teorii nelineynykh kolebaniy [Asymptotic methods in the theory of nonlinear vibration]. Moscow, Nauka Publ., 1974. 504 p.
[5] Moiseev N.N. Asimptoticheskie metody nelineynoy mekhaniki [Asymptotic methods of nonlinear mechanics]. Moscow, Nauka Publ., 1969. 381 p.
[6] Butenin N.V., Neymark Yu.I., Fufaev N.L. Vvedenie v teoriyu nelineynykh kolebaniy [Introduction to the theory of nonlinear vibration]. Moscow, Nauka Publ., 1987. 382 p.
[7] Landau L.D., Lifshits E.M. Teoreticheskaya fizika. V 10 t. T. 1. Mekhanika [Theoretical physics. Ten-volume set. Vol. 1. Mechanics]. Moscow, Nauka Publ., 1988. 214 p. (Eng. Ed.: Landau L.D., Lifshitz E.M. Mechanics. Volume 1 (Course of Theoretical Physics). Second Ed. Oxford, New York, Pergamon Press, 1969.)
[8] Korn G.A., Korn T.M. Mathematical handbook for scientists and engineers: definitions, theorems, and formulas for reference and review. USA, McGraw-Hill, 1968. 1130 p. (Russ. ed.: Korn G., Korn T. Spravochnik po matematike dlya nauchnykh rabotnikov i inzhenerov. Moscow, Nauka Publ., 1984. 831 p.)
[9] Andronov A.A., Vitt A.A., Khaykin S.E. Teoriya kolebaniy [Theory of oscillations]. Moscow, Nauka Publ., 1981. 918 p. (Eng. ed.: Andronov A.A., Vitt A.A. and Khaykin S.E. Theory of Oscillations. Translated from the Russian by F. Immirzi. Oxford-N.Y.-Toronto, Ont., Pergamon Press, 1966).
Статья поступила в редакцию 5.11.2013
Судаков Владимир Федорович — д-р техн. наук, профессор кафедры "Теоретические основы электротехники" МГТУ им. Н.Э. Баумана. Автор более 140 научных работ в области теоретической электротехники, квантовой электроники, радиолокации и радионавигации.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
Soudakov V.F. — Dr. Sci. (Eng.), professor of "Theoretical Fundamentals of Electrical Engineering" department of the Bauman Moscow State Technical University. Author of more than 140 publications in the field of theoretical electrical engineering, quantum electronics, radar and radio navigation.
Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.