Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Теория автоматического управления

..pdf
Скачиваний:
11
Добавлен:
05.02.2023
Размер:
1.36 Mб
Скачать

100

ω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 c2 ,

 

 

π

 

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

 

 

 

 

Критерий Найквиста позволяет оценить устойчивость САУ, содержащих звенья чистого запаздывания.

Пусть звено чистого запаздывания с передаточной функци-

ей epτ (при единичном коэффициенте передачи) включено последовательно с системой без запаздывания с передаточной

функцией W0 ( p) .

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

W ( p) =W0 ( p)epτ , W ( jω) =W0 ( jω)ejωτ .

Поскольку 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 — Варианты ЛАЧХ и ЛФЧХ для условно устойчивой САУ