
Теория автоматического управления
..pdf100
ω2 = 24,9 с–1, ω1 = 24,254 с–1. Поскольку ω0 < ω1 < ω2 , САУ
снова устойчива.
Таким образом, область устойчивости расположена с внешней стороны границы, как показано штриховкой на рис. 4.8. Асимптоты a0 = 0 и a3 = 0 являются прямыми Kp = 0 , T1 = 0 , т.е. совпадают с осями координат.
Вариант 2 . Проведем D-разбиения в плоскости параметров x1 = τ1 , x2 = Kp . Тогда
a0 (x2 ) = x2 , a1(x1, x2 ) =1+ x2 (x1 + τ3 ) ,
a2 (x1, x2 ) =T1 +T2 + x2 x1τ3 , a3 (x1) =T1T2 = a3 .
Система уравнений (4.4) примет вид:
a |
(x ) |
−a |
2 |
(x |
, x |
)ω2 = x −(T +T + x x τ |
3 |
)ω2 |
= 0, |
||||||||
|
0 |
2 |
|
1 |
2 |
|
2 |
|
1 |
2 |
2 |
1 |
|
|
(4.7) |
||
|
|
|
|
|
|
|
|
|
(x |
|
) − a ω2 |
|
|
|
|||
ω[a (x ) − a ω2 ] = ω[1 |
+ x |
+ τ |
] = 0. |
|
|||||||||||||
|
|
1 |
2 |
|
3 |
|
|
2 |
1 |
3 |
|
3 |
|
|
|
|
|
Система (4.7) является нелинейной, поскольку в обоих ее уравнениях присутствует произведение варьируемых параметров, поэтому ее можно решить только путем подстановки.
Выразим из второго уравнения системы переменную x2 и подставим получившееся значение в первое уравнение:
|
|
|
|
|
|
|
|
|
|
a ω2 −1 |
|
|
|
|
|
|
|||||
|
|
|
|
|
x (x ,ω) = |
|
3 |
|
|
|
, |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
2 |
1 |
|
|
|
x1 +τ3 |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
a |
ω2 −1 |
|
|
|
|
a |
ω2 |
−1 |
|
|
|
|
|
|
||||||
|
3 |
|
|
|
− T |
+T |
+ |
3 |
|
|
|
x τ |
ω2 |
= 0 , |
|||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
x1 + τ3 |
|
1 |
2 |
|
|
|
|
1 |
3 |
|
|
|
|
|
||||||
|
|
|
|
|
x1 + τ3 |
|
|
|
|
|
|
||||||||||
отсюда |
|
|
|
|
[a − τ (T +T )]ω2 −1 |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
x |
(ω) = |
|
3 |
3 |
1 |
2 |
|
|
|
|
|
|
, |
|
|||||
|
|
ω2 (ω2a |
τ +T |
+T |
− τ |
|
|
|
|||||||||||||
|
|
1 |
|
|
|
3 |
) |
|
|
||||||||||||
|
|
|
|
|
|
|
|
3 |
3 |
|
1 |
2 |
|
|
|
|
|
||||
|
|
x |
(ω) = |
ω2 (ω2a |
τ +T |
+T |
−τ ) |
|
|||||||||||||
|
|
|
|
3 3 |
|
1 |
2 |
|
|
3 . |
|||||||||||
|
|
2 |
|
|
|
|
ω2τ 2 |
+1 |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|

|
|
|
|
|
|
|
101 |
|
|
|
|
|
|
|
|
В полученных соотношениях присутствуют две критиче- |
|||||||||||||||
ские |
частоты: |
ω |
= 0 |
и |
ω |
|
= |
τ3 −T1 −T2 |
(уравнение |
||||||
|
|
|
кр1 |
|
|
|
кр2 |
|
|
a3τ3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ω2τ 2 |
+1 = 0 дает мнимые решения), причем, для заданных зна- |
||||||||||||||
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
чений τ3 , T1, |
T2 частота ωкр2 |
также принимает только мнимые |
|||||||||||||
значения, поэтому при проведении дальнейших расчетов нужно |
|||||||||||||||
избегать только значения ω= 0 . |
|
|
|
|
|
x1(ω) |
|
||||||||
Рассчитанные значения |
варьируемых |
параметров |
и |
||||||||||||
x (ω) приведены втабл. 4.2 при изменении частотыот 7 до 21 с–1. |
|||||||||||||||
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 4.2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ω |
7 |
|
9 |
11 |
|
|
13 |
|
|
15 |
17 |
19 |
|
21 |
|
x1(ω) |
0,038 |
|
0,05 |
0,055 |
|
0,056 |
|
0,056 |
0,055 |
0,054 |
0,051 |
||||
x2 (ω) |
30 |
|
50 |
77,8 |
|
112 |
|
155 |
206 |
268 |
|
342 |
|||
На рис. 4.9 показана область устойчивости САУ для принятых |
|||||||||||||||
варьируемых параметров. Она также располагается с внешней сто- |
|||||||||||||||
роны границы устойчивости, |
поскольку при 50 < Kp < 400 |
нару- |
|||||||||||||
шается чередование частот ω2 |
и ω1 . |
|
|
|
|
|
|
|
|||||||
x2(ω) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
600 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
500 |
|
|
|
|
|
|
|
|
|
САУ устойчива |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
400 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
300 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
200 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
100 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
0,04 |
0,045 |
|
0,05 |
0,055 |
0,06 x1(ω) |
||||||||
|
|
|
|||||||||||||
Рис. 4.9 — Область устойчивости САУ при x1(ω) = τ1, |
x2(ω) = K р |
|
102
Следует отметить, что построение области устойчивости с помощью D-разбиений является более общим подходом к решению этой задачи. Так, например, если бы в примере 4.3 данная задача решалась посредством критерия Гурвица, то граница устойчивости представлялась бы двумя кривыми x21 = f1(x1) и
x22 = f2 (x1) , точку сопряжения которых определить весьма
сложно. При проведении D-разбиений этот вопрос разрешается автоматически, поскольку граница устойчивости — единая зависимость.
4.5 Критерий устойчивости Найквиста
Критерий устойчивости Найквиста позволяет оценить устойчивость замкнутой САУ по ее разомкнутой цепи. Для этого в передаточной функции Wрц( p) производят замену оператора p
на переменную jω и на комплексной плоскости при изменении
частоты от нуля (если это возможно) до бесконечности строят АФЧХ W ( jω) (годограф Найквиста).
Если разомкнутая цепь устойчива (а это всегда имеет место, если САУ не содержит неустойчивых неминимально-фазовых звеньев), то формулировка критерия Найквиста звучит следующим образом.
Для устойчивости замкнутой САУ необходимо и достаточно, чтобы годограф Найквиста при изменении ω от 0 до
∞ не охватывал точку с координатами (−1, j0) .
На рис. 4.10 изображены основные из возможных ситуаций прохождения годографа Найквиста на комплексной плоскости. Сплошная кривая 1 на рис. 4.10, а соответствует абсолютно устойчивой замкнутой САУ (системе, которая остается устойчивой при уменьшении коэффициента передачи разомкнутой цепи), а пунктирная кривая 2 — условно устойчивой САУ (системе, устойчивой только в некотором диапазоне изменения коэффициента передачи разомкнутой цепи, как это было в примере 4.4). Сплошная кривая 3 на рис. 4.3, б проходит через критиче-
скую точку с координатами (−1, j0) , и это означает, что замкнутая система находится на колебательной границе устойчивости.

103
Пунктирная кривая 4 охватывает критическую точку, поэтому замкнутая САУ неустойчива.
Физический смысл критерия Найквиста заключается в том, что при увеличении частоты входного воздействия сигнал, проходящий по цепи обратной связи, оказывается в противофазе с входным. Это равносильно замене отрицательной обратной связи на положительную. Если же при этой частоте разомкнутый
контур обладает усилением (т.е. Kp >1), то замкнутая САУ ста- |
||||
новится неустойчивой (любое увеличение сигнала на выходе |
||||
приводит к увеличению сигнала на входе по цепи обратной свя- |
||||
зи, что вызывает дальнейший рост выходного сигнала и т.д.). |
||||
Im[W(jω)] |
|
|
|
|
|
|
Im[W(jω)] |
|
|
|
|
Re[W(jω)] |
|
|
–1 |
0 |
ω = 0 |
|
Re[W(jω)] |
|
|
–1 |
0 |
ω = 0 |
|
|
ω |
|
ω |
|
|
1 |
|
|
|
|
2 |
|
4 |
|
|
а |
|
3 |
|
|
|
б |
|
Рис. 4.10 — Варианты годографа Найквиста для устойчивой (а), |
||||
неустойчивой САУ и САУ на границе устойчивости (б) |
Для аналитических расчетов с помощью критерия Найквиста условия нахождения системы на границе устойчивости можно записать, используя вещественную и мнимую частотные функции разомкнутой цепи:
Re[W ( jω )] = −1, |
|
|
|
π |
|
Im[W ( jω )] = 0, |
(4.8) |
|
|
π |
|
|
|
|
где ωπ — частота, соответствующая повороту радиус-вектора
АФЧХ разомкнутой цепи на угол −π, т.е. до совпадения с отрицательной вещественной полуосью; ее называют частотой пе-
реворота фазы.

104
При решении практических задач для оценки устойчивости САУ не обязательно строить годограф Найквиста, достаточно в частотной передаточной функции разомкнутой цепи W ( jω)
приравнять к нулю мнимую часть и определить из получившегося уравнения частоту переворота фазы ωπ (или ее квадрат). Затем подставить получившееся значение в вещественную часть W ( jω) и вычислить ее модуль. Если Re[W ( jωπ)] <1, то систе-
ма устойчива, в противном случае — неустойчива.
Несмотря на наглядность критерия Найквиста и его физическую прозрачность, он имеет один существенный технический недостаток — вычислительные трудности, возникающие при разделении вещественной и мнимой частей W ( jω) . Осо-
бенно это проявляется, если САУ содержит форсирующие звенья. Если же таких звеньев в структуре нет, то задача анализа устойчивости по критерию Найквиста решается просто. Продемонстрируем это на примере.
Пример 4.5
По критерию Найквиста получить выражение для оценки устойчивости и определить граничный коэффициент передачи для САУ, изображенной на рис. 4.3.
В примере 4.1 была получена передаточная функция разомкнутой цепи
Wрц( p) = W ( p) = |
|
|
k1k2koc |
= |
|||
|
(T1 p +1)(T2 p + |
1)(Toc p +1) |
|||||
|
|
|
|
||||
= |
|
|
Kp |
|
. |
|
|
(T p + |
1)(T p |
+1)(T p +1) |
|
||||
1 |
2 |
oc |
|
|
|
Произведем в Wрц( p) замену оператора p на переменную jω и выделим в знаменателе получившегося выражения мнимую и вещественную части:

105
|
|
W ( jω) = |
|
|
|
|
Kp |
|
|
|
|
= |
|
|
||
|
|
(1+ jωT1 )((1+ jωT2 ))(1+ jωToc ) |
|
|
||||||||||||
|
|
= |
|
|
|
|
|
Kp |
|
|
|
= |
|
|
|
|
|
|
[1 |
−ω2T T |
+ jω(T |
+T )](1+ jωT |
) |
|
|
|
|||||||
|
|
|
|
|
|
1 2 |
|
1 |
2 |
|
oc |
|
|
|
|
|
= |
|
|
|
|
|
|
|
Kp |
|
|
|
|
|
|
= |
|
1 |
−ω2[T T |
−(T |
+T |
)T |
] + jω[T |
+T |
+T |
(1 |
−ω2T T |
)] |
||||||
|
|
1 2 |
1 |
2 |
oc |
|
1 |
2 |
oc |
|
1 2 |
|
|
= Kp . u(ω) + jy(ω)
По правилам деления комплексных чисел [1] числитель и знаменатель полученного выражения нужно умножить на комплексную сопряженную функцию u(ω) − jy(ω) . Однако по-
скольку при определении частоты переворота фазы мнимая часть приравнивается к нулю и знак перед ней не играет роли, то операцию умножения на комплексную сопряженную функцию можно не проводить, приняв
Im[W ( jωπ)] =T1 +T2 +Toc (1−ω2πT1T2 ) = 0 ,
Re[W ( jωπ)] = |
|
|
|
Kp |
|
. |
|
−ω2 |
[T T |
−(T |
|
||
1 |
+T )T ] |
|||||
|
|
π |
1 2 |
1 |
2 oc |
Выразив из первого уравнения квадрат частоты ωπ , полу-
чим ω2 |
= |
T1 +T2 +Toc |
, тогда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
π |
|
T1T2Toc |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
Re[W ( jωπ)] |
|
= |
|
|
|
|
|
|
|
Kp |
|
|
|
|
|
|
|
= |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
1− |
T1 +T2 +Toc |
[T T + (T +T )T |
] |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
T1T2Toc |
|
|
1 2 |
1 |
|
|
2 oc |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
= |
|
|
|
|
|
|
|
|
Kp |
|
|
|
|
|
|
|
|
. |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
1 |
−(T |
+T |
+T |
) |
1 |
+ |
1 |
+ |
1 |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
1 |
2 |
oc |
|
|
|
|
T2 |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
T1 |
|
|
Toc |
|
|
|
|

106
При Re[W ( jωπ)] <1 система будет устойчивой. На границе устойчивости Re[W ( jωπ)] = −1 , т.е.
|
|
|
Kгp |
|
|
|
|
|
|
|
|
|
|
= −1 , |
|||
|
|
|
|
1 |
|
1 |
|
1 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ + |
|
|
|
|
|
|
|||||
1−(T1 +T2 +Toc ) |
|
|
|
|
|
|
|
|
|
||||||||
отсюда |
|
T1 |
|
T2 |
Toc |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
1 |
|
1 |
|
|
1 |
|
|
|||||
Kгp |
= (T1 +T2 |
+Toc ) |
+ |
|
+ |
|
|
−1 . |
|||||||||
|
|
|
|
|
|||||||||||||
|
|
|
|
|
T1 T2 |
|
|
Toc |
|
Зададимся конкретными значениями коэффициентов передачи и постоянных времени. Пусть k1 = 5, k2 =10, koc = 0,5 ,
T1 = 0,5 с, |
T2 = 0,1 с, |
Toc = 0,05 с. Тогда по полученным соот- |
||
ношениям |
при K |
p |
= 5 10 0,5 = 25 получим ω2 |
= 260 c−2 , |
|
|
π |
|
Re[W ( jωπ)] =1,263 >1, следовательно, система с такими пара-
метрами будет неустойчивой. Ее граничный коэффициент передачи
|
|
1 |
|
1 |
|
1 |
|
|
Kгp |
= (0,5 +0,1+0,05) |
|
+ |
|
+ |
|
|
−1 =19,8 < Kp . |
|
0,1 |
0,05 |
||||||
|
0,5 |
|
|
|
|
Критерий Найквиста позволяет оценить устойчивость САУ, содержащих звенья чистого запаздывания.
Пусть звено чистого запаздывания с передаточной функци-
ей e−pτ (при единичном коэффициенте передачи) включено последовательно с системой без запаздывания с передаточной
функцией W0 ( p) .
Результирующие передаточная и частотная передаточная функции разомкнутой цепи будут иметь вид
W ( p) =W0 ( p)e−pτ , W ( jω) =W0 ( jω)e− jωτ .
Поскольку W0 ( jω) = A0 (ω)e jϕ0 (ω) , то
W ( jω) = A0 (ω)e j[ϕ0 (ω)−ωτ] .
Таким образом, звено чистого запаздывания вносит лишь дополнительный фазовый сдвиг. При этом изменяется АФЧХ,

107
т.е. меняются условия устойчивости (характеристика «закручивается» по часовой стрелке). При некотором τ система станет неустойчивой.
По АФЧХ системы без запаздывания можно определить граничное (предельное) значение запаздывания τгр , что поясня-
ется построением на рис. 4.11.
Пусть АФЧХ устойчивой САУ без запаздывания W0 ( jω) пересекает окружность единичного радиуса на частоте ω= ωcp
(как будет показано ниже, это частота среза) при повороте ради- ус-вектора АФЧХ на угол ϕcp . При введении в САУ звена чис-
того запаздывания на границе устойчивости конец этого радиусвектора совпадет с точкой (−1, j0) и будет справедливым соотношение ϕcp −ωcpτкp = −π, отсюда
τгp = π+ ϕcp .
ωcp
Im[W0(jω)]
W0(jω)
|
|
|
–1 |
0 ϕ Re[W0(jω)] |
ω = ωср
Рис. 4.11 — Определение граничного запаздывания
108
4.6Оценка устойчивостиСАУпо логарифмическим частотным характеристикам. Запасы устойчивости
В соответствии с критерием Найквиста об устойчивости можно судить не только по АФЧХ, но и совместно по амплитудной и фазовой частотным характеристикам разомкнутой цепи. Обычно при этом пользуются логарифмическими характеристиками, что представляет большое удобство в силу простоты их построения. Но если ЛАЧХ используется асимптотическая, то расчеты будут достаточно грубыми.
Если АФЧХ не охватывает точку (−1, j0) , то при частоте, на которой А(ω) =1, абсолютное значение фазы больше −π. Но значение A(ω) =1 соответствует G(ω) = 20lg A(ω) = 0 . Поэтому
для устойчивости замкнутой САУ необходимо, чтобы ЛАЧХ разомкнутой цепи пересекла ось абсцисс раньше, чем фаза, спадая, окончательно перейдет за значение –π.
Однако G(ω) = 0 на частоте среза ωcp , а ϕ(ω) = −π на частоте переворота фазы ωπ . Следовательно, система будет абсолютно устойчива, если ωcp < ωπ .
На рис. 4.12 для ЛАЧХ, показанной сплошной линией, это неравенство соблюдается (ωcp1 < ωπ ) , значит, САУ устойчива, Kp1 < Kгp . При увеличении коэффициента передачи Kp ЛАЧХ
смещается вверх и на границе устойчивости (пунктирная кривая) при Kp2 = Kгp частоты ωcp2 и ωπ равны друг другу. Дальнейшее увеличение величины Kp до Kp3 приводит к неустой-
чивости системы, тогда ωcp3 > ωπ (штрихпунктирная кривая).
Возможен и более сложный случай. Как было показано в примере 4.4, САУ устойчива при Kp1 = 40 с–1. При этом ее
АФЧХ будет дважды пересекать отрицательную вещественную полуось, как показано на рис. 4.10, а, при частотах ωπ1 и ωπ2 . ЛАЧХ такой системы изображена сплошной линией на рис. 4.13, и

109
ее частота среза ωcp1 меньше любой из частот ωπ1 и ωπ2 . При Kp2 =100 с–1 частота среза ωπ1 < ωcp2 < ωπ2 (на рис. 4.13 ЛАЧХ показана пунктирной линией) и САУ становится неустойчивой. При Kp3 = 500 с–1 САУ снова станет устойчивой, хотя частота
среза ωcp3 (ЛАЧХ изображена штрихпунктирной линией на рис. 4.13) больше частот ωπ1 и ωπ2 .
G(ω) 20lgKp3
20lgKp2 20lgKp1
0 |
ωcp1 ωcp2 ωcp3 lgω |
ϕ(ω)
0 |
ωπ |
lgω |
−π
Рис. 4.12 — Оценка устойчивости САУ по ЛАЧХ и ЛФЧХ
G(ω) 20lgKp3 20lgKp2 20lgKp1
0
ϕ(ω 0
−π
cp1 |
cp2 |
cp3 |
π1 π2
lg
lg
Рис. 4.13 — Варианты ЛАЧХ и ЛФЧХ для условно устойчивой САУ