ГЛАВА 17. ТОЧНЫЕ МЕТОДЫ ИССЛЕДОВАНИЯ УСТОЙЧИВОСТИ И АВТОКОЛЕБАНИЙ

§ 17.1. Фазовые траектории и метод точечных преобразований

Понятие о фазовом пространстве, о фазовых траекториях и их типах было уже дано выше. В данном параграфе на примерах построения фазовых траекторий для простейших систем второго порядка будут проиллюстрированы некоторые важные особенности процессов в нелинейных системах автоматического регулирования.

Пример 1. Возьмем систему автоматического регулирования с объектом без самовыравнивания и с приводом регулирующего органа, имеющим постоянную скорость. Уравнение регулируемого объекта без самовыравнивания будет

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

где — относительные изменения регулируемой величины, смещений чувствительного элемента, регулирующего органа, элемента обратной связи и управляющего золотника (рис. 10.11, а), — коэффициент. Привод регулирующего органа пусть имеет постоянную скорость в двух вариантах; 1) с мгновенным переключением (рис. 16.22, ж) при переходе управляющего элемента (золотника, струйной трубки) через нейтральное положение с зоной нечувствительности (рис. 16.22, з) вследствие наличия «перекрытия» золотника или струйной трубки. В первом случае уравнение привода регулирующего органа будет

а во втором

Возьмем фазовую плоскость приняв

Из уравнений (17.1), (17.2) и (17.5) имеем

Следовательно, переключения привода в первом варианте будут

иметь место при

что соответствует прямой (рис. 17.1, а) на фазовой плоскости, причем согласно (17.16) значениям соответствует часть плоскости слева от прямой , а — справа.

На основании первого из соотношений (17.6) с учетом (17.3) при получаем

а из (17.5)

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

или, после интегрирования,

Это есть семейство парабол, показанное на рис. 17.1, а справа от линии (они симметричны относительно оси Так как (17.8) и (17.9) являются проекциями скорости изображающей точки М на оси х и у, то имеем а знак совпадает со знаком у.

Рис. 17.1.

В соответствии с этим на рис. 17.1, а укажем стрелочками направление движения изображающей точки М по фазовым траекториям. Аналогичным путем легко строятся параболы слева от прямой

В результате, как видно из общего расположения фазовых траекторий (рис. 17.1, а), получается устойчивая система с затухающим колебательным переходным процессом. Но число колебаний будет конечным. В самом деле, здесь имеется особый отрезок в который вливаются все фазовые траектории. Чтобы выявить поведение системы на этом отрезке, вспомним, что для него согласно (17.7) и (17.5)

Следовательно, попав на отрезок изображающая точка не может с него уйти, и система будет апериодически приближаться к установившемуся

состоянию, т. е. изображающая точка будет сползать по отрезку к началу координат О. Таким образом, имевший место вначале колебательный переходный процесс после конечного числа колебаний вырождается в этот так называемый скользящий процесс.

Крайние точки особого отрезка определяются, очевидно, как точки, в которых прямая касается одной из парабол соответственно правого и левого семейств. Поэтому, подставив значения — из (17.7) в выражение (17.10), найдем точку С:

По найденной картине расположения фазовых траекторий можно качественно представить себе кривую переходного процесса при любых начальных условиях. Начальными условиями определяется начальное положение изображающей точки М и тем самым — определенная фазовая траектория, иллюстрирующая протекание процесса. Она показывает (рис. 17.1, а) максимальное отклонение регулируемой величины максимальную скорость а также все последующие отклонения, число колебаний и т.

Рассмотрим теперь ту же систему, но с учетом зоны нечувствительности.

В этом случае переключениям привода (при на фазовой плоскости соответствуют согласно (17.6) две наклонные прямые (рис. 17.1, б):

Между этими прямыми правее их левее их а (причем

При из (17.4), (17.6) и (17.5) получаем

откуда (при )

(прямые, параллельные оси х в полосе на рис, 17.1, б).

При получим прежние параболы. В результате снова система оказывается устойчивой и имеет колебательный переходный процесс, но вместо особой точки О получаем особый отрезок , т. е. установившееся состояние определяется неоднозначно. Это соответствует тому, что регулятор может находиться в равновесии в любом месте внутри зоны нечувствительности. Здесь точно так же возможен скользящий процесс, как и в случае рис. 17.1, а.

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

Пример 2. Допустим, что требуется стабилизировать угловое положение некоторого тела, когда сопротивлением среды его вращению можно пренебречь. Уравнение объекта будет

где — момент инерции тела — угол поворота тела, — его угловая скорость, М — управляющий момент со стороны исполнительного органа системы стабилизации.

Уравнение регулятора (системы стабилизации) запишем в виде

где — постоянная положительная величина, — нелинейный закон регулирования, осуществляемый при помощи логического устройства по тому же простейшему принципу, что и на рис. 16.26, с той лишь разницей, что по углу фазовая плоскость ограничена значениями так как. это составляет один полный оборот тела (рис. 17.2).

Рис. 17.2.

Изобразим процесс регулирования на фазовой плоскости. Уравнение всей системы согласно (17.11) и (17.12) будет

где обозначено

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

Умножив почленно уравнение (17.13) на выражение

получим дифференциальное уравнение фазовой траектории

Это уравнение легко интегрируется внутри участков, на которых . В результате для каждого отдельно взятого участка уравнение фазовой траектории будет

где — значения в начальной точке данного участка.

Зададим начальные условия процесса:

Для данной начальной точки процесса (см. рис. 17.2) имеем Поэтому на первом участке процесса согласно (17.15) уравнение фазовой траектории

будет

Этот участок движения с постоянной скоростью заканчивается в точке 1 (рис. 17.2), где происходит включение исполнительного органа . Следовательно, для второго участка процесса (после точки 1) из (17.15) получим уравнение фазовой траектории

так как в начальной точке 1 этого участка Фазовая траектория (17.16) — парабола, ось которой совпадает с координатной осью

Это соответствует равнозамедленному движению Изображая параболу графически, доводим ее до границы (участок 2—2 на рис. 17.2), причем в точке 2 согласно (17.16)

Это значение переносим в точку 2 (для вращающегося тела — это одна и та же точка). Здесь происходит выключение исполнительного органа Поэтому дальнейшее движение согласно (17.15) пойдет с постоянной скоростью

до точки 3 (рис. 17.2). Таким образом, в рассмотренной начальной части процесса регулирования тело совершило один полный оборот, но в конце этого оборота скорость вращения его стала меньше начальной.

В точке 3 снова включается исполнительный орган в результате чего фазовая траектория будет

так как в точке Допустим, что соответствующая уравнению (17.18) парабола 3—4 не доходит до границы Это означает, что тело больше не совершит полного оборота, а начнет (с точки А) возвращаться в сторону нулевого положения.

В точке 4 (рис. 17.2) имеем скорость Следовательно, из (17.18) угловая координата ее будет

где определяется по формуле (17.17). Дальше (4—5) процесс пойдет с постоянной скоростью (так как после чего тело войдет в установившийся автоколебательный режим, определяемый предельным циклом (5—6—7—8). Уравнение параболы 7—8 согласно (17.15) будет

Отсюда амплитуда угловых автоколебаний акак значение при будет

а амплитуда колебаний скорости

Она равна зоне нечувствительности измерителя угловой скорости в то время как амплитуда угловых колебаний (17.19) несколько больше зоны нечувствительности измерителя угла

Период автоколебаний можно вычислить как сумму времен:

где — времена участков соответственно. По законам равномерного и равнозамедленного движений соответственно получаем

Итак, установившийся режим стабилизации в данной системе является автоколебательным. Однако уравнение системы (17.13) справедливо только для идеальной системы стабилизации. Всякое реально имеющееся запаздывание в работе усилительно-преобразовательного и исполнительного устройств приведет к увеличению амплитуд автоколебаний по сравнению с полученными здесь значениями. Решение задачи с учетом постоянных времени системы управления будет дано в следующей главе.

Пример 3. Уравнения системы автоматического регулирования курса водяной торпеды в упрощенном варианте имеют вид: линейная часть (16.40) и (16.41), т. е.

и нелинейное звено (возьмем сначала один случай — рис. 16.18, в)

Покажем, что здесь равновесное установившееся состояние системы с постоянным значением неустойчиво, но будет иметь место устойчивый автоколебательный процесс.

Возьмем фазовую плоскость с координатами (угол отклонения и угловая скорость отклонения оси торпеды от заданного курса). Уравнения (17.20) и (17.21) перепишутся в виде

Из сравнения этих уравнений с упрощенными уравнениями системы регулирования температуры в конце § 16.1 видна их полная аналогия. Поэтому здесь, так же как и в случае рис. 16.15, установившийся процесс движения торпеды будет автоколебательным, причем картина фазовых траекторий будет иметь вид, показанный на рис. 17.3, а.

При этом кривая предельного цикла, соответствующая автоколебательному процессу, определяется из уравнения (16.31) с таким значением произвольной постоянной чтобы выполнялось условие

так как только в этом случае и получится замкнутая кривая предельного цикла (рис. 17.3, а). Определив таким образом найдем амплитуду автоколебаний а как значение х при т. е. согласно (16.31)

Значения же (17.23) дают амплитуду q колебаний скорости у. Можно все это определять и графически прямо по чертежу (рис. 17.3, а). Период автоколебаний остается неизвестным.

Введем теперь в характеристику нелинейного звена (рулевой машинки) зону нечувствительности, как показано на рис. 17.3, 6, в. Так, на том участке характеристики (рис. 17.3, б), где из (17.22) следует, что

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

Рис. 17.3.

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

При дальнейшем увеличении зоны нечувствительности характеристика нелинейного звена и картина фазовых траекторий принимают вид, показанный на рис. 17.3, в. Здесь автоколебания отсутствуют и становится устойчивым установившийся процесс с постоянным значением Ранее неустойчивый особый отрезок теперь стал устойчивым. Дальнейшее увеличение зоны нечувствительности приводит к расширению отрезка т. е. к увеличению установившейся ошибки системы из-за слишком широкого участка равновесия.

Пример 4. Рассмотрим вибрационный регулятор напряжения, уравнения которого были составлены в § 16.2, а именно;

причем уравнение нелинейного звена (регулирующего органа)

В качестве ординаты фазовой плоскости здесь удобнее взять не скорость отклонения регулируемой величины как делалось раньше, а вторую

переменную Итак, примем для этой задачи

Тогда уравнения (17.24) преобразуются к виду

где согласно (17.25), (17.26) и (17.28) имеем

следовательно, первое из этих условий имеет место ниже прямой (рис. 17.4), а второе — выше нее.

Рис. 17.4.

В первом случае переключение реле происходит при , т. е. на прямой (рис. 17.4), а во втором случае — при т. е. на прямой Чертеж сделан в предположении, что . В результате получаем, что выше линии будет

а ниже линии

Рассмотрим сначала верхнюю область. Для нее, деля (17.28) на (17.27), с учетом (17.30) получим уравнение фазовых траекторий

которое можно представить в виде

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

где z — новая переменная вместо у. В результате найдем следующее уравнение фазовых траекторий

где — произвольная постоянная,

при решение будет иметь другой вид, а при будет эти решения не будут исследоваться).

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

и ко всем ординатам этой прямой добавлять

придавая произвольные значения (каждому значению будет соответствовать определенная фазовая траектория). Это будут параболы степени с осью

и с единым началом в точке Н (рис. 17.4), имеющей координаты

На рис. 17.4 показаны все ветви этих парабол, лежащие выше линии (так как только там справедливы данные выкладки). Направления стрелок на полученных фазовых траекториях определяются тем, что проекция скорости изображающей точки — справа от прямой (17.37) согласно (17.27) будет отрицательна, а слева — положительна; проекция же согласно (17.28) выше прямой будет отрицательна, а ниже — положительна (во всех точках прямой касательные к фазовым траекториям горизонтальны).

Аналогично строятся и все фазовые траектории ниже линии так как их дифференциальное уравнение отличается от (17.32) только заменой на согласно (17.31).

В результате на рис. 17.4 видим, что все фазовые траектории, исходящие из особого отрезка расходятся, а все траектории, идущие от краев чертежа, сходятся. Как те, так и другие асимптотически приближаются к установившемуся предельному циклу, обозначенному на чертеже жирной замкнутой кривой (чичевицеобразной). Это соответствует тому, что установившийся процесс в системе является автоколебательным, причем размеры предельного цикла представляют собой амплитуды автоколебаний соответственно регулируемого напряжения и тока в обмотке электромагнита реле

Определить фазовую траекторию, образующую этот предельный цикл, можно как такую кривую (17.33), у которой

чем определяется значение произвольной постоянной Значение для этой кривой и дает искомую амплитуду Амплитуда же определяется как ордината пересечения кривой предельного цикла с прямой (ибо, как было показано ранее, в точках этой прямой касательные к фазовым траекториям горизонтальны).

Из чертежа (рис. 17.4) видно, что предельный цикл лежит левее точки и охватывает точку С. Поэтому , т. е. амплитуда автоколебаний регулируемого напряжения заключена в интервале

где определяются формулами (17.34). Амплитуда же будет немного больше .

Пример 5. Рассмотрим следящую систему с сухим трением в управляемом объекте, для которой уравнения были написаны в § 16.3.

Рис. 17.5.

Уравнение регулируемого объекта (16.52) как нелинейного звена при отсутствии линейного трения имеет вид

При написании уравнения линейной части системы (16.53) пренебрежем постоянными времени (чтобы иметь возможность рассматривать уравнение всей системы как уравнение второго порядка), а именно;

Подставив это в уравнения объекта (17.39) и обозначив

получим уравнение всей следящей системы в целом:

За координаты фазовой плоскости примем, как обычно, Условие при котором согласно (17.42) будет т. е. система будет в равновесии, изображается на фазовой плоскости отрезком (рис. 17.5).

Вне этого отрезка согласно (17.41) необходимо отдельно рассмотреть два случая; т. е. верхнюю и нижнюю половины фазовой плоскости. При из (17.41) имеем

Это уравнение совпадает с уравнением (16.23), но со сдвигом на величину Следовательно, ниже оси х надо нанести такие же кривые, как на рис. 16.9, б (если или как на рис. (если но со сдвигом начала координат в точку А, что и сделано на рис. 17.5, а и соответственно.

Аналогичные кривые наносятся и выше оси х, но только со сдвигом начала координат в точку В (рис. 17.5), так как согласно (17.41) при имеем уравнение

В обоих случаях (рис. 17.5, а и б) система устойчива, причем в первом случае переходный процесс состоит из конечного числа затухающих колебаний управляемого объекта, а во втором случае имеем апериодическое движение. Положение равновесия объекта определяется неоднозначно, объект может остановиться в любой точке особого отрезка (рие. 17.5), как это было уже ранее при наличии зоны нечувствительности (см. пример 1). Особый отрезок определяется соотношением где с — абсолютное значение момента сухого трения при движении управляемого объекта.

Заметим, что произведенное здесь упрощение уравнений системы хотя и позволило решить их точно, но это решение, дающее в результате устойчивость системы при любых числовых значениях параметров системы, неполно отражает действительную картину явлений в данной нелинейной системе.

Метод точечного преобразования.

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

Пусть (рис. 17.6, а) граничными линиями между кусками фазовых траекторий являются ось х, линия и линия

Рис. 17.6.

Возьмем начальное положение изображающей точки где-нибудь на полуоси Один этап движения системы состоит в переходе изображающей точки на линию ограничивающую этот этап, в некоторое положение (рис. 17.6, а). Следующий этап переводит изображающую точку в положение на полуоси затем в положение на кривой и, наконец, в положение на исходной полуоси

Каждому положению на полуоси соответствует определенное положение точки на кривой Это называется точечным преобразованием полупрямой в кривую . Для краткости ему присваивают название, например: преобразование Дальше (рис. 17.6, а) идет

точечное преобразование кривой в полупрямую названное затем — точечное преобразование полупрямой в кривую и преобразование кривой в исходную полуось

Все это в целом (или, как говорят, преобразование называется точечным преобразованием полупрямой самой в себя. Это преобразование в данном случае записывается в виде определенной зависимости:

где через обозначены абсциссы точек (рис. 17.6, а). Если при любом оказывается то в системе будет затухающий процесс, а если — расходящийся процесс. Если же возможно равенство то на фазовой плоскости получится предельный цикл, который, как известно, может изображать либо устойчивый автоколебательный процесс, либо границу устойчивости системы в малом, либо может соответствовать особому случаю бифуркации (см. ниже).

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

В рассматриваемом случае верхняя полуплоскость симметрична нижней относительно начала координат. Поэтому достаточно рассмотреть только первую половину преобразования т. е. точечное преобразование в полупрямую и выразить его в виде зависимости

причем условие наличия предельного цикла на фазовой плоскости будет при

Пусть, например, зависимость (17.43) имеет вид кривой, показанной на рис. Проведем на этом графике еще прямую из начала координат под углом 45° к координатным осям. Если она пересечет кривую, то в точке пересечения получим о Чтобы определить, какому типу предельного цикла это соответствует, надо взять на оси абсцисс начальную точку сначала слева, а затем справа от точки пересечения и проследить ход точечного преобразования, как показано стрелками на рис. 17.6, 6. В данном случае процесс сходится с обеих сторон к точке пересечения. Следовательно, здесь будет устойчивый предельный цикл, соответствующий автоколебательному процессу в системе. При этом абсцисса точки пересечения (рис. 17.6, б) дает амплитуду автоколебаний.

Можно поступить и иначе. Допустим, преобразование от точки к точке выполняется достаточно просто, но оказывается, что по полученной определенной точке находить соответствующую точку труднее, чем по заданному положению определять Тогда будем подходить к кривой с двух сторон, задавая одновременно точку на полуоси и точку на полуоси и находя соответствующие точки (рис. 17.6, в). В результате получим точечные преобразования полупрямых и в кривую выраженные некоторыми зависимостями:

Изобразив это в виде двух кривых (рис. 17.6, г) анализируем их тем же способом, как и кривую с прямой на рис. 17.6, б.

Такие графики называют диаграммами точечного преобразования. Они соответствуют в данном случае устойчивому предельному циклу, т. е. наличию установившегося автоколебательного процесса в системе. Другие возможные типы диаграмм точечного преобразования показаны на рис. 17.7. При этом рис. 17.7, а соответствует неустойчивому предельному

циклу; он ограничивает область начальных условий при которых система оказывается устойчивой относительно установившегося состояния с постоянным значением регулируемой величины При начальных же условиях выходящих за контур этого предельного цикла, система неустойчива (система устойчива в малом и неустойчива в большом).

Рис. 17.7, в соответствует двум предельным циклам, из которых меньший неустойчив, а другой (больший) устойчив. Следовательно, при начальных условиях расположенных внутри первого предельного цикла, система устойчива, как и в предыдущем случае, а при всяких других начальных условиях она стремится к установившемуся автоколебательному процессу, который определяется вторым предельным циклом.

Рис. 17.7.

Этот случай может выродиться в случай, изображенный на рис. 17.7, б, когда оба предельных цикла сливаются в один полу устойчивый. Подобные особые случаи называются бифуркационными.

Наконец, на рис. 17.7, г и д изображены случаи, когда на диаграмме точечного преобразования кривая не пересекается с прямой, проведенной под углом 45° к осям. Это означает соответственно устойчивость и неустойчивость системы при любых начальных условиях (до которых справедливы исследуемые уравнения системы).

Заметим, что изложенное выше является лишь качественным рассмотрением, так как в нем отсутствует время Остается неизвестным течение процессов во времени, период (частота) автоколебаний. Чтобы получить полное решение задачи, нужно к данному рассмотрению добавить еще решение дифференциальных уравнений на отдельных участках во времени (как в методе припасовывания, см. § 16.1). Поэтому в методе точечного преобразования вводится соответствующий временной параметр (здесь это не рассматривается).

Для систем выше второго порядка вместо фазовой плоскости придется иметь дело с фазовым пространством и с точечным преобразованием не линий, а поверхностей. Там появляются новые особенности процессов. Однако ввиду чрезвычайной сложности таких построений рассматривать их не будем.

Метод изоклин.

Выше были рассмотрены такие примеры нелинейных систем второго порядка, для которых фазовые траектории легко находятся интегрированием уравнений по участкам. В тех случаях, когда интегрирование затруднено, ход фазовых траекторий, хотя бы качественно, можно проследить с помощью так называемого метода изоклин (без интегрирования уравнений). Количественно этот способ имеет сравнительно низкую точность. Применение его пока ограничено системами второго порядка.

Изоклиной называется такая линия, во всех точках пересечения которой с фазовыми траекториями последние наклонены под одним и тем же определенным углом к оси абсцисс х. Так, если известно дифференциальное уравнение фазовых траекторий

то для получения изоклины нужно положить

Рис. 17.8.

Уравнение изоклины, следовательно, будет

где с обозначает определенный тангенс гугла наклона фазовых траекторий. Каждому заданному значению с соответствует своя изоклина.

Например, часто встречается нелинейное уравнение

которое можно записать в виде

тогда дифференциальное уравнение фазовых траекторий будет выглядеть так:

а уравнение изоклин

Задавая различные значения с (при заданном к), для каждого из них строим по этому уравнению кривую на фазовой плоскости — изоклину (сплошные кривые на рис. 17.8). Затем на каждой кривой наносим стрелочки под углами с к оси абсцисс (на рис. 17.8 указаны значения с для каждой кривой). Из этих стрелочек и составляются искомые фазовые траектории; некоторые из них изображены на рис. 17.8 пунктиром. В данном случае получается устойчивый предельный цикл, что соответствует автоколебаниям в системе.