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

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

..pdf
Скачиваний:
125
Добавлен:
15.11.2022
Размер:
16.96 Mб
Скачать

мых

корней

pk =*= ±

/рл (колебательная

граница

устойчивости),

то функция

F (Jсо) при (о = 0 или со =

обратится в нуль.

В

практических

расчетах удобно применять

с л е д с т в и е

и з к р и т е р и я М и х а й л о в а :

 

 

система устойчива, если действительная и мнимая части ха­

рактеристической функции F (/со)

обращаются в нуль пооче­

редно (рис.

5.2, г),, т. е. если корни

уравнений

Р (со) = 0

и Q (со) = 0

(5.37)

перемежаются.

Это утверждение вытекает непосредственно из формулировки критерия Михайлова — из условия последовательного прохожде­ ния кривой F (/со) через п квадрантов.

Критерий Михайлова удобно применять для анализа устойчи­ вости систем высокого порядка (п > 5).

Пример. Определим с помощью критерия Михайлова устойчивость си­ стемы, характеристический полином которой имеет вид

 

F (р) =

а Лръ+ а гр4 +

а2р2 +

а3р2 + д4р + д5,

(5.38)

где

д0 = 0,16; д х =

1,7;

а 2 =

6,0; д3 =

9,6; д4 =

6,7; д 3 = 8.

ную

Подставляя в выражение (5.38) р =

j со, найдем мнимую и действитель­

составляющие

характеристической

функции:

 

 

Q (со) =

а0со5 — д2со3 +

а4со,

 

 

(5.39)

 

Р (со) =

Д 4СО4 — ДдСО2 +

05-

 

 

(5.40)

Подставляя в выражения (5.39) и (5.40) конкретные значения па­ раметра со, можно вычислить координаты ряда точек кривой F (/со) и по ее виду оценить устойчивость системы'. Однако этот путь требует сравнительно много вычислений (обычно 10— 15 точек в интервале со = 0 оо), поэтому воспользуемся следствием из критерия Михайлова.

Найдем положительные корни уравнений

д0соб — а2со3 +

д4со = 0,

(5.41)

Д] СО4 —

ДдСО2 +

д5 = 0.

(5.42)

Первое уравнение всегда имеет один нулевой корень соМ1 =

0. Вычис­

лим корни

уравнения

 

Д 0СО4 — Д 2С02 + д4 = 0,

для чего введем промежуточную переменную z = со2: д022 — a2z + д4 = 0.

Корни вспомогательного уравнения (5.44)

После подстановки числовых значений получим:

Z\ = 1,15;

= 35.

Положительные корни исходного уравнения (5.43):

соМ2 = 1» 1; ^мз = 5,9.

Аналогично найдем положительные корни уравнения (5.42):

(0д1 = 1,0; (0д2 = 2,2.

(5.43)

(5.44)

(5.45)

(5.46)

(5.47)

(5.48)

171

Нетрудно заметить, что корни мнимой и действительной составляющих чередуются, следовательно, система устойчива.

Если хотя бы один из корней вспомогательного уравнения получится отрицательным или комплексным, то это будет означать, что не существует соответствующего положительного корня ш, при котором кривая F (/to) пересекает оси Р (со) или /Q (ш). Кривая F (/со) в этом случае не пройдет че­

рез л квадрантов, и система будет неустойчивой.

5.4. Критерий Найквиста

Критерий был сформулирован в 1932 г. американским физиком X. Найквистом, занимавшимся исследованием свойств электрон­ ных усилителей с обратной связью. Позднее А. В. Михайлов обос­ новал этот критерий и показал возможность применения его для анализа автоматических систем управления.

В отличие от критериев Гурвица, Рауса и Михайлова, которые основаны на анализе характеристического уравнения системы, критерий Найквиста позволяет судить об устойчивости системы по амплитудно-фазовой характеристике разомкнутого контура си­ стемы. В этом заключается существенное преимущество критерия, так как построение амплитудно-фазовой характеристики разомкну­ того контура для большинства реальных систем оказывается проще, чем построение годографа Михайлова. Особенно упрощается это построение для одноконтурных систем, состоящих из типовых звеньев. А в тех случаях, когда неизвестно математическое опи­ сание одного или нескольких конструктивных элементов системы и оценка их свойств возможна только путем экспериментального определения частотных характеристик, критерий Найквиста яв­ ляется единственно пригодным.

О с н о в н а я ф о р м у л и р о в к а к р и т е р и я Н а й ­ к в и с т а :

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

тудно-фазовая характеристика W (jco)

разомкнутого контура

не охватывает точку с координатами (— 1; /0).

Эта формулировка справедлива д л я

с и с т е м , к о т о р ы е

в р а з о м к н у т о м с о с т о я н и и

у с т о й ч и в ы . Тако­

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

На рис. 5.3, а изображены амплитудно-фазовые характеристики разомкнутого контура, соответствующие трем различным случаям: система устойчива (кривая /); система находится на колебательной границе устойчивости (кривая 2); система неустойчива (кривая 3).

Употребленное в формулировке критерия Найквиста понятие охвата точки имеет некоторую неопределенность, из-за чего в слу­ чаях сложной формы кривой W (/©) могут возникнуть затруднения в суждении об устойчивости системы. Поэтому для большей ясно­ сти рекомендуется следующий прием. Надо проследить мысленно

172

а

jQ(v)n

6

 

 

Рис. 5.3. Амплитудно-фазовые характеристики разомкнутого контура

(а)

и физическая трактовка (б) критерия

Найквиста

 

 

за движением вектора W 1 (/со) =

1 + W (/to),

вращающегося

во­

круг точки (— 1; /0) и скользящего по кривой

W (/со). Угол пово­

рота вектора W 1 {j<o), равный я,

означает охват точки (— 1; /0),

а угол, меньший я — неохват.

 

 

 

Для использования изложенного приема применительно к аста­ тическим системам, которые содержат интегрирующее звено, и ам­ плитудно-фазовые характеристики которых начинаются в — оо на мнимой оси, характеристику W (/to) предварительно дополняют в четвертом квадранте дугой окружности бесконечно большого радиуса. ^

Для суждения об устойчивости систем, имеющих а. ф. х. слож­ ной конфигурации, когда кривая а. ф. х. пересекает действительную ось левее точки (— 1; /0) несколько раз, можно также использовать п р а в и л о п е р е х о д о в , сформулированное советским уче­ ным Я. 3. Цыпкиным: а. ф. х. не охватывает точку (— 1; /0), т. е. система устойчива, если при возрастании to разность между числом положительных (сверху вниз) и отрицательных (снизу вверх) пере­ ходов а. ф. х. через ось абсцисс слева от точки (— 1; /0) равна нулю.

Если а. ф. х. начинается или заканчивается на отрезке (— оо;

1), то считают, что характеристика совершает полперехода. Частота, при которой амплитудная характеристика A (to) [мо­

дуль функции W (jin) ] принимает значение 1, называется частотой

среза и

обозначается toCp. Частоту, при которой фазовый сдвиг

Ф (to) =

— я, обозначают юя.

Пользуясь введенными обозначениями, можно записать условие

нахождения системы на границе устойчивости:

 

tocp = а>„.

(5.49)

Частота, с которой система колеблется на границе устойчивости, называется критической и обозначается toKp.

Дадим ф и з и ч е с к у ю т р а к т о в к у основной формули­ ровки к р и т е р и я Н а й к в и с т а . Предположим, что на входе

g (() =

системы (рис. 5.3, б) действует гармонический сигнал

= gmsin a t e малой амплитудой gm. Пусть частота © равна частоте ©я, при которой фазовый сдвиг <р (©), создаваемый звеном W (/со), равен — я. Тогда сигнал отрицательной обратной связи окажется в фазе с сигналом g (/), и мгновенные значения сигналов будут суммироваться.

Если на частоте со = соя модуль | W (/©я) | = 1, т. е. выполняется условие (5.49), то в контуре системы будут поддерживаться неза­ тухающие колебания даже после исчезновения внешнего воздейст­ вия g (t), т. е. система будет находиться на границе устойчивости. Характеристика W (/со) при этом проходит через точку (— 1; /0).

Если на частоте © = соя модуль | W (/©я) | < 1, то после исчез­ новения внешнего воздействия колебания в контуре затухнут, т. е. система устойчива, характеристика не охватывает точку (— 1; /0).

Если же модуль | W (/©„) | > 1, то амплитуда сигналов в кон­ туре будет неограниченно возрастать, т. е. система будет неустой­ чивой. Характеристика W (/©) в этом случае охватит точку (— 1;

/ 0).

Таким образом, особая роль точки (— 1; /0) заключается в том, что она, во-первых, соответствует превращению отрицательной обратной связи в положительную, и во-вторых, является гранич­ ной между режимами усиления и ослабления сигналов звеном W (/©).

Иногда на практике встречаются системы, в контуре которых имеется одно или несколько неустойчивых элементов. Такие с и - с т е м ы в р а з о м к н у т о м с о с т о я н и и н е у с т о й ­ ч и в ы . Для суждения об их устойчивости необходимо использо­ вать д р у г у ю ф о р м у л и р о в к у к р и т е р и я Н а й ­ к в и с т а :

автоматическая система управления устойчива, если ампли­ тудно-фазовая характеристика W (/©) разомкнутого контура охватывает 112-раз точку с координатами (— 1; /0), где I число правых корней характеристического уравнения разомкну­ того контура.

Количество охватов при этом можно определять по правилу Цыпкина как разность между числом положительных и отрицатель­ ных переходов.

Данная формулировка критерия Найквиста является более об­ щей, чем предыдущая. Действительно, если разомкнутая система устойчива (т. е. 1 = 0), то для устойчивости замкнутой системы а. ф. х. W (/©) должна точку (— 1; /0) охватывать нуль раз, т. е. не охватывать.

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

решается сравнительно легко, при помощи

критерия Гурвица:

для этого приравнивают к нулю знаменатель

передаточной функ-

174

ции W (р) разомкнутого контура и анализируют данное характе­ ристическое уравнение.

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

Пример 1. Определим с помощью критерия Найквиста максимально до­ пустимое значение общего передаточного коэффициента системы, состоящей из трех инерционных звеньев первого порядка с одинаковыми постоянными

времени Т х =

Т 2 =

Т 3 = Т.

 

 

Передаточная функция разомкнутого контура системы

 

W (р) = kl(Tp + 1)3.

 

 

 

(5.50)

Амплитудная

частотная функция

контура

 

А (со) = | W (/со) |

=

kl (V Т2со2 +

1 )з,

(5 .5 1 )

фазовая частотная

функция

 

 

Ф (со) = arg W (/со) =

— 3 arctg соТ.

(5.52)

Система будет находиться на границе устойчивости, если а. ф. х. разом­

кнутого контура пройдет через точку (— 1; /0 ), т. е. если при

некоторой ча­

стоте

со = сокр одновременно выполняются условия:

 

А

(^кр) —

 

ф (®кр) =

— я*

 

(5.53)

Для рассматриваемой системы условия (5.53) имеют вид

 

f t/( '\A ’4

p +

1

У =

I,

 

(5.54)

— 3 arctg сОкрТ =

— я.

 

(5.55)

Из условия (5.55) имеем

 

 

сокрГ = л/3 .

 

 

 

 

 

(5.56)

Подставляя выражение (5.56) в условие (5.54), получим

искомое зна­

чение передаточного коэффициента k = 8 .

значение пе­

Из приведенного

решения следует также, что предельное

редаточного коэффициента не зависит от абсолютного значения постоянных времени Т.

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

W (/со) = W' (/со)

(5.57)

где W ' (/со) — эквивалентная

амплитудно-фазовая функция ос­

тальных звеньев.

Характеристику W (/со) строят следующим образом. Вначале строят кривую W' (jсо), а затем каждый вектор, соответствующий частоте со* поворачивают на угол оуг (рис. 5.4, б).

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

CL

Рис. 5.4. Оценка устойчивости системы с запаздыванием:

а — структура системы; б — а. ф. х. разомкнутого контура

Пример 2 . Определим, будет ли устойчива статическая система, состоя­ щая из пропорционального управляющего устройства и статического инер­ ционного объекта первого порядка с запаздыванием (3.130), при следующих значениях параметров:

k = k yk 0 = 7;

Г0 = 100 с; то =

20

 

 

 

 

 

 

 

 

 

Амплитудно-фазовая характеристика разомкнутого контура

 

W (/(о) =

W y (/(D) W 0 (/(D) = k е - /шхо/(Т0/ш + 1).

 

 

 

 

(5.58)

Амплитудная

функция

контура

 

 

 

 

 

 

 

 

 

А (ш) =

 

A/д /

Т 2ш2 + 1 ,

 

 

 

 

 

 

 

 

 

 

 

(5.59)

фазовая функция

 

 

 

 

 

 

 

 

 

 

 

 

 

ср (©) =

 

— arctg ©Т0 — ©т0.

 

 

 

 

 

 

 

 

 

 

(5.60)

Найдем вначале частоту ©я , при которой ср (©я) =

л.

Решая методом

последовательных приближений

трансцендентное

уравнение

 

 

 

— arctg ©я 100 — ©л20 =

— я,

 

 

 

 

 

 

 

 

 

(5.61)

получим соя » 0,085

с-1.

 

А (со)

при

частоте сол =

0,085:

 

 

Теперь

 

вычислим

значение

 

 

А ((О) |ю='),085 = 7/УЮО2 (0.085)2 +

1

«

0,82 <

1.

 

 

 

 

(5.62)

Следовательно, а. ф. х. не охватит точку (— 1; /0).

Система устойчива.

Найдем

критическое запаздывание

т0. Кр — значение

т0,

при

котором

рассматриваемая

система

станет

неустойчивой.

Для

этого

из

условия

Л (©Ср) =

1

определим ©ср =

0,069

с-1.

При этой частоте

инерционная

часть разомкнутого контура создает фазовый сдвиг

фи (соСр) =

а гс tg

©срТ’о = — 82°. Соответственно

на

запаздывающую часть

приходится угол

фз (Юср) =

— 180° — (— 82°) =

— 98°.

Отсюда

критическое

запаздывание

то. кр =

 

фз (toCp)/57,3©cp « 25 с.

 

 

 

 

 

 

 

 

 

(5.63)

Если разомкнутый контур системы образован последовательным соединением типовых динамических звеньев, то целесообразно ча­

стотную

характеристику контура строить в

л о г а р и ф м и ч е ­

с к о й

с и с т е м е к о о р д и н а т

и об

устойчивости системы

судить по виду этой характеристики.

При этом используют р а з -

176

Рис. 5.5. Логарифмические частотные характерис­ тики статических систем:

1 — устойчивой; 2 — находящейся на границе устойчивос­ ти; 3 — неустойчивой

н о в и д н о с т ь о с н о в н о й ф о р м у л и р о в к и к р и ­ т е р и я Н а й к в и с т а :

система устойчива, если при достижении фазовой частотной характеристикой значения — 180° логарифмическая амплитуд­ ная характеристика будет отрицательной (рис. 5.5, кривые 1).

Действительно, если L (©) < 0, то А (ю) <

1. Поэтому отрица­

тельность L (со) при ф (©„) = — 180° свидетельствует о том,

что

а. ф. х. разомкнутого контура не охватывает точку (— 1; /0).

 

Логарифмические частотные характеристики

L (со) и ф (со)

ра­

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

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

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

Рис. 5.6. Пример расчета устойчивости по логари­ фмическим частотным характеристикам

первого порядка с

постоянными времени Т г = 0,2 с,

= 0 ,1 с, Tz =

= 0,05 с. Передаточный коэффициент разомкнутого контура k = 20.

Сопрягающие частоты звеньев

 

(дс1Т~\ = 5с- 1 ;

©с2 = Т ^ 1= 10с” 1; а с3 = Т3 = 20с- 1 .

 

В соответствии с этими значениями построены приближенная л. а. ч. х. L (о)) (рис. 5.6, а) и фазовая характеристика <р (©) (рис. 5.6, б) разомкнутого контура, которая получена согласно выражению

Ф (©) = — arctg

— arctg ®Т2— arctg ©Т3.

(5.64)

По графику видим, что при ф (©) = — 180° функция L (©) > 0, следо­ вательно, система неустойчива.

5.5. Построение областей устойчивости

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

а

к Область

неустойчивости

Плосность параметров Плоскость корней

Рис. 5.7. Области устойчивости в плоскостях коэффициентов (а) и корней

(б) характеристического уравнения

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

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

На рис. 5.7, а показана область устойчивости, построенная в про­ странстве двух коэффициентов характеристического уравнения ас и ak. Каждой точке, находящейся ниже заштрихованной кривой, например, точке А \ соответствуют только левые корни, например, точки р\ и р2 (рис. 5.7, б). Любой точке, находящейся выше кривой (точке А '") обязательно соответствуют хотя бы один действитель­ ный корень или пара комплексных корней, расположенных справа на плоскости корней (точки р{" и р2"). Границей области устойчи­ вости в данном примере является заштрихованная кривая. Каждой точке этой кривой, например, точке Л ", соответствует пара чисто мнимых корней (точки р\ и /?г).

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

кой путь связан с большим объемом вычислений. Эффективным спо­ собом отыскания границ областей устойчивости является метод D-разбиения, разработанный советскими учеными А. А. Соколовым и Ю . И. Неймарком.

Сущность метода D -разбиения заключается в следующем. Пусть анализируемая система описывается характеристическим уравне­ нием п-го порядка

а0рп+ ахр"-1+ <hPn~2+

-+ а п= 0.

(5.65)

Расположение всех п

корней характеристического

уравнения

на комплексной плоскости а—/р зависит от значений коэффициен­ тов а0, аг, , ап. В общем случае в пространстве варьируемых коэффициентов, например, а{ и ак (см. рис. 5.7), существуют такие значения коэффициентов, при которых / корней расположены справа, а (л—/) корней — слева от мнимой оси /р . Совокупность всех таких значений образует в пространстве коэффициентов об­ ласть, которую можно обозначить D (л—/; /). Очевидно, что могут

существовать

области и с

другим распределением корней:

D (л—/—1; / +

1), D (л—I—2;

/ + 2) и т. д. Область D (л; 0) яв­

ляется областью устойчивости.

Процесс построения в пространстве параметров (или коэффи­ циентов) областей с различным распределением корней называется D-разбиением. Линии, разграничивающие эти области, называются

кривыми D -разбиения.

Шереход из одной области в другую область пространства ко­ эффициентов соответствует переходу одного действительного или пары комплексных корней через мнимую ось /р . Следовательно, каждой точке, лежащей на границе между двумя областями, со­ ответствуют либо нулевой корень р = 0, либо пара чисто мнимых корней р = ± /р . Поэтому кривую D -разбиения можно рассматри­ вать как отображение мнимой оси плоскости корней. Указанная особенность кривых D -разбиения используется при отыскании их уравнений. Для этого в уравнение (5.65) подставляют р = / Р или р = /со и разрешают уравнение относительно варьируемых пара­ метров. Совокупность значений варьируемых параметров, соот­ ветствующих всем возможным значениям © (от — оо до + оо), дает все точки кривой D -разбиения. Выделение области D (п; 0) среди остальных областей производят при помощи специальной процедуры — штриховки кривых D -разбиения.

Построение области устойчивости по одному параметру■ При­

меним метод D -разбиения для построения области устойчивости в плоскости одного параметра. Предположим, что варьируемый параметр k входит в характеристическое уравнение линейно:

F(p) = kA(p) + B(p) = 0,

(5.66)

где А (р) и В (р) — степенные полиномы от р.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]