В исследуемом случае критический коэффициент усиления определяется при касании L(wр) оси частот. Перенесем ЛАХ параллельно самой себе так, чтобы в точке w=wр она касалась оси частот и продлим первую асимптоту до пересечения с осью частот. В этой точке k=w=7.244, что соответствует значению (ku)кр=16.74.
Выделение областей устойчивости
Среди физических параметров, характеризующих САУ, всегда имеется несколько, легко поддающихся изменению и использующихся для определенной настройки системы. При конструировании системы весьма важно знать диапазоны значений изменяемых параметров, допустимые с точки зрения сохранения устойчивости САУ. Об этих диапазонах можно судить, если в пространстве изменяемых параметров построить область устойчивости, т.е. выделить область значений параметров, при которых система сохраняет устойчивость.
Область устойчивости в теории автоматического управления принято называть D – областью, а представление области параметров в виде областей устойчивости и неустойчивости называют D – разбиением.
Построение области устойчивости по алгебраическим критериям
Допустим, что коэффициенты характеристического уравнения
зависят от двух изменяемых параметров m и l. Для построения области устойчивости прежде всего нужно, в соответствии с необходимым условием устойчивости, выделить область изменяемых параметров при нахождении в которой, коэффициенты характеристического уравнения положительны. Это можно сделать, решив систему уравнений
(3.26)Для построения границы положительности коэффициентов аi необходимо из решений уравнений (3.26) выбрать те, которые обеспечивают положительность всех коэффициентов. Из всех границ положительности только две одновременно могут быть и границами устойчивости. Такими являются границы, уравнениями которых являются
(3.27)Доказано, что если d0 и dnприблизятся к нулю, то характеристическое уравнение будет иметь два действительных корня
(3.28)При дальнейшем уменьшении коэффициенты d0 и dnперейдут через ноль, станут отрицательными, а корни (3.28) окажутся положительными. Так как вещественные корни определяют апериодические составляющие решения дифференциального уравнения, то границы (3.27) называют апериодическими границами устойчивости. На самих границах устойчивости корни (3.28) равны соответственно ±¥ и 0. Стороны кривых, di(m,l)=0, примыкающие к области положительности соответствующих коэффициентов, штрихуются в сторону положительности. Может случиться так, что какой либо из коэффициентов, d0 или dn не зависит от изменяемых параметров. Это означает отсутствие соответствующей апериодической границы устойчивости.
Колебательной границей устойчивости называется кривая в плоскости изменяемых параметров, при переходе через которую пара комплексно – сопряженных корней изменяет знак своей вещественной части на обратный. Доказано, что колебательная граница устойчивости определяется выражением
(3.29)В этом выражении Dn-1 – (n-1) – й определитель Гурвица. Колебательная граница устойчивости штрихуется в сторону положительности Dn-1.
Пример. Построить область устойчивости в плоскости параметров kuи kwz системы стабилизации угла тангажа.
Характеристическое уравнение замкнутой системы имеет вид
Исследуем неравенства d2>0, d3>0, d4>0 . Из первого неравенства следует, что для положительности коэффициента d2 необходимо, чтобы выполнялось условие
Неравенство d4>0 определяет, что для положительности этого коэффициента необходимо, чтобы ku>0. Для выполнения неравенства d3>0 требуется, чтобы
При любых значениях передаточного числа по углу больших нуля, правая часть последнего выражения по модулю будет больше единицы. Таким образом, границами положительности коэффициентов будут
От изменяемых параметров зависит коэффициент dn=d4 и не зависит коэффициент d0. Поэтому уравнение ku=0 одновременно является и апериодической границей устойчивости.
Составив определитель Гурвица, для его Dn-1 минора получим
Подставим в это выражение значения коэффициентов d2, d3, d4, как функций параметров ku и kw , после преобразований получим квадратное уравнение, определяющее передаточное число по угловой скорости как функцию от передаточного числа по углу тангажа
По этому выражению строится колебательная граница устойчивости. График деления области исследуемых параметров на области устойчивости и неустойчивости показан на рис. 3.19.
Граница колебательной неустойчивости штрихуется в сторону положительности Dn-1- го определителя Гурвица, а прямая kwz=0 в сторону положительности этого коэффициента.
Для проверки полученных результатов выберем какие – либо значения параметров внутри заштрихованной области, например ku=5, kwz=0.6, вычислим значения коэффициентов характеристического уравнения и оценим устойчивость замкнутой системы по критерию Гурвица.
Получим, что при выбранных значениях передаточных чисел система устойчива. Это означает, что и вся область, внутрь которой обращены штрихи, является областью устойчивости.
D – разбиение в плоскости одного параметра
Пусть нас интересует влияние какого – либо одного параметра на устойчивость САУ и этот параметр входит в характеристическое уравнение линейно, так что это уравнение можно представить в виде
(3.30)Сделав замену s= jw , получим
(3.31)Задавая значения частоты от -¥ до +¥, можно построить кривую m(w), отображающую мнимую ось плоскости корней на плоскость m. Эта граница D – разбиения симметрична относительно вещественной оси. Поэтому вычисления можно вести в диапазоне частот от 0 до +¥, а затем дополнить полученную кривую ее зеркальным отображением на диапазон частот от -¥ до нуля. При движении по мнимой оси от -¥ до +¥ на плоскости корней область устойчивости остается слева. Поэтому при движении по кривой D – разбиения в сторону увеличения частоты ее штрихуют слева. Область, внутрь которой обращены штрихи, является предполагаемой областью устойчивости. Для окончательного решения, необходимо взять какое – либо вещественное значение параметра m в исследуемой области и воспользоваться каким – либо критерием устойчивости. Если при избранном значении параметра система устойчива, то рассматриваемая область является областью устойчивости.
Пример. Построить область устойчивости системы стабилизации угла тангажа в плоскости передаточного числа ku.
Характеристическое уравнение исследуемой системы можно записать в виде
где
В полученных выражения сделаем замену s=jw и получим
В этих выражениях
Построенная по этим выражениям кривая D – разбиения показана на рис. 3.20.
Так как необходимым условием устойчивости рассматриваемой системы является ku>0, то мнимая ось также является границей устойчивости и штрихуется в сторону положительности ku. Значение этого коэффициента, равное 5, находится внутри заштрихованной области и мы знаем, что при этом значении система устойчива. Значит и весь отрезок вещественной оси, расположенный внутри заштрихованной области, дает значения передаточного числа по углу, при которых система устойчива. Можно показать, что окончание этого отрезка находиться в точке, равной критическому значению коэффициента ku=16.56.
D – разбиение в плоскости двух параметров
Пусть коэффициенты характеристического уравнения линейно зависят от двух параметров m и l так, что его можно записать в виде