Критерий устойчивости Гурвица 1 страница

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

Вводная часть к критерию Гурвица.

Пусть замкнутая САУ имеет следующую передаточную функцию

Из коэффициентов характеристического многочлена

составляем следующую матрицу

(1.45)

Порядок заполнения матрицы следующий. Сначала по диагонали матрицы располагают коэффициенты от a₁ до a_n. Затем над диагональными элементами располагают коэффициенты с возрастающими индексами. Если коэффициенты в процессе заполнения все исчерпаны, то ставят 0. Далее под диагональными элементами располагают коэффициенты с убывающими индексами. Если коэффициенты в процессе заполнения все исчерпаны, то ставят 0.

Формулировка критерия Гурвица: САУ устойчива, если:

1) положительны все коэффициенты характеристического многочлена;

2) положительны все n главных (диагональных) определителей Гурвица матрицы (1.45):

(1.46)

Последний определитель не вычисляют, так как его знак совпадает со знаком Δ_n-1:

Δ_n=a_n×Δ_n-1

Если хотя бы один определитель Гурвица отрицателен, то САУ неустойчива. Если имеется хотя бы один определитель Гурвица равен нулю при остальных положительных, то САУ находится на границе устойчивости.

Числовой пример.

Определить устойчивость САУ с передаточной функцией

(1.47)

Составляем матрицу Гурвица и главные определители

Оба определителя положительны, поэтому САУ устойчива.

Определение допустимых настроек САУ.

Если передаточная функция САУ содержит хотя бы один буквенный коэффициент, значение которого может быть любым числом, то с помощью критерия Гурвица можно определить допустимые по условию устойчивости значения такого коэффициента. При двух буквенных коэффициентах возможно совместное определение допустимых значений таких коэффициентов и выделение областей устойчивости на плоскости этих коэффициентов. Покажем это на примере.

Пусть САУ управления курсом судна, представленная на рис.1.34, состоит из двух звеньев - авторулевого и судна. Передаточные функции авторулевого и судна имеют вид, соответственно,

Постоянная времени T судна зависит от загрузки судна и изменяется от 10 с при порожнем судне до 60 с при полностью загруженном. Параметром настройки авторулевого является коэффициент передачи K. Необходимо найти такие значения параметра K, при которых САУ устойчива при изменении загрузки судна.

Определяем передаточную функцию замкнутой САУ

(1.48)

Составляем матрицу Гурвица и вычисляем 2-й определитель:

(1.49)

С учётом положительности всех коэффициентов характеристического многочлена при условии (1.49) САУ будет устойчива при одновременном выполнении следующей системы неравенств

(1.50)

САУ будет находиться на границе устойчивости, если будет выполнено хотя бы одно из равенств

. (1.51)

Каждое из равенств (1.51) является на плоскости T-K границей области устойчивости (рис.1.35). Штриховкой обозначены области устойчивости по отношению к линиям границы устойчивости. Общая область a0d для всех заштрихованных областей является областью устойчивости САУ.

Пусть при порожнем судне с Т_порож установлен коэффициент передачи K₁ регулятора. Этот состояние САУ отмечено точкой 1, лежащей в области устойчивости. Если после загрузки судна значение T увеличится до T_груж, то при том же K₁ САУ в точке 2 будет неустойчива. Необходимо будет увеличить K до значения K₃, чтобы система оказалась в точке 3. Область, ограниченная ломаной линией abcd, будет областью устойчивости при любой загрузке судна.

Вопросы и задания

1. Дайте определение устойчивости. Назовите прямые методы устойчивости и дайте характеристику области их применимости.

2. Что представляют собой критерии устойчивости? В чем состоят их преимущества перед прямыми методами?

3. Приведите формулировку критерия устойчивости Гурвица и продемонстрируйте его использование для оценки САУ 3-го порядка с произвольной передаточной функцией.

4. Поясните процедуры определения допустимых настроек САУ по условию устойчивости.

1.10. Критерий устойчивости Михайлова.

Определение допустимых настроек САУ

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

Подготовительные операции.

Пусть замкнутая САУ имеет следующую передаточную функцию

Из коэффициентов характеристического многочлена

образуем годограф Михайлова

(1.52)

Вычисляем и строим годограф D(jω) при значениях ω = 0 … ∞ (рис.1.36).

Формулировка критерия Михайлова: САУ п -го порядка устойчива, если годограф Михайлова, начинаясь на положительной части действительной оси, проходит против часовой стрелки подряд ровно п квадрантов (рис.1.36а). Если годограф проходит п квадрантов не подряд или проходит меньше квадрантов, чем п, то САУ неустойчива (рис.1.36б). Если годограф при любой частоте проходит через начало координат, то САУ находится на грани устойчивости (рис.1.36в).

Оценка устойчивости с использованием графика годографа Михайлова наглядна, но требует громоздких вычислений по выражению (1.52). Если учесть, что для оценки устойчивости важен только факт прохождения годографом п квадрантов против часовой, а не сама форма годографа (рис.1.36а), то на этом можно существенно сократить вычисления. Достаточно найти частоты ω₁, ω₂, …, ω_п, при которых пересекаются оси координат, упорядочить их по возрастанию и для этих частот найти последовательно знаки (+ или -) выражений Р(ω₁), Q(ω₂), P(ω₃),…, где Р и Q - соответственно, действительные и мнимые части выражения W(jω). САУ будет устойчива, если знаки выражений Р и Q будут чередоваться согласно табл.1.6. Любое отступления от табл.1.6 свидетельствует о том, что САУ либо неустойчива, либо находится на грани устойчивости.

Числовой пример.

Определить устойчивость САУ с передаточной функцией (1.47). Характеристический многочлен от (1.47) и годограф Михайлова имеют вид

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

, откуда ω₁=0 и

Для этих частот рассчитываем значения Р (ω):

В точках пересечения годографа с мнимой осью выполняется равенство

, откуда

Для этой частоты рассчитываем значение Q (ω):

Составляем таблицу по типу табл.1.6.

Чередование знаков Р (ω) и Q (ω) правильное и, поэтому, САУ устойчива.

Определение допустимых настроек САУ.

Если передаточная функция САУ содержит хотя бы один буквенный коэффициент, значение которого может быть любым числом, то с помощью критерия Михайлова можно определить допустимые по условию устойчивости значения такого коэффициента. При двух буквенных коэффициентах возможно совместное определение допустимых значений таких коэффициентов и выделение областей устойчивости на плоскости этих коэффициентов. Покажем это на примере, рассмотрев САУ управления курсом судна, представленную на рис.1.34.

Характеристический многочлен согласно (1.48) и годограф Михайлова имеют вид

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

, откуда ω₁=0 и .

Для этих частот рассчитываем значения Р (ω):

Значение Р (ω₂) будет отрицательным числом при T<100K или K>0,01^.T, что совпадает с (1.50). При этом соотношении или .

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

, откуда

Для этой частоты рассчитываем значение Q (ω):

при K>0,01^.T.

Составляем таблицу по типу табл.1.6.

Чередование знаков Р(ω) и Q(ω) правильное и, поэтому, САУ при K>0,01^.T устойчива.

Вопросы и задания

1. Сформулируйте критерий устойчивости Михайлова и продемонстрируйте его использование для оценки САУ 3-го порядка с произвольной передаточной функцией.

2. Поясните процедуры определения допустимых настроек САУ по условию устойчивости.

1.11. Критерий устойчивости Найквиста.

Использование ЛАЧХ для оценки устойчивости САУ

Критерий Найквиста является частотным. САУ должна иметь структуру с единичной обратной связью (рис.1.37). Для оценки устойчивости используется частотная характеристика разомкнутой САУ.

Подготовительные операции.

Вычисляем и строим годограф W_раз (jω) частотной характеристики разомкнутой САУ при значениях ω = 0 … ∞ (рис.1.38).

Формулировка критерия Найквиста зависит от вида частотной характеристики W_раз (jω) разомкнутой САУ. Ниже рассмотрим только два вида W_раз (jω), которые чаще всего встречаются на практике.

1). Пусть разомкнутая САУ устойчива. Тогда замкнутая САУ будет устойчива, если годограф W_раз (jω) Найквиста, начинаясь на положительной части действительной оси, не охватывает точку –1 действительной оси (рис.1.38а). Если годограф охватывает точку –1 действительной оси (рис.1.38б), то САУ неустойчива (рис.1.38б), а если проходит через точку –1 действительной оси, то САУ находится на грани устойчивости (рис.1.38в).

2). Пусть разомкнутая САУ обладает астатизмом 1-го порядка, являясь САУ, находящейся на грани устойчивости. В этом случае начало годографа W_раз (jω) разомкнутой САУ при ω=0 находится на бесконечности (рис.1.39). Необходимо такой годограф дополнить дугой бесконечно большого радиуса, начинающейся на положительной части действительной оси и заканчивающийся на начале годографа W_раз (jω). Тогда к такому годографу применимы все правила, по которым оценивалась устойчивость как на рис.1.38.

Использование ЛАЧХ для оценки устойчивости САУ

Введем на построениях критерия устойчивости Найквиста косвенные показатели качества – запас по фазе γ и частоту среза ω_ср (рис.1.40). Запас по фазе γ характеризует удаленность годографа W_раз (jω) от точки с координатой –1 на вещественной оси. Для устойчивости САУ достаточно, чтобы запас по фазе был бы положительным. Чем больше запас по фазе, тем меньше перерегулирование на графике переходного процесса (см. рис.1.41). Считается для качественной САУ достаточно иметь запас по фазе в пределах 30^о…50^о.

Из построений, приведенных на рис.1.40, следует, что:

1) запас по фазе можно рассчитать по выражению частотной характеристики W_раз (jω), подставив в него значение частоты среза ω_ср.

2) на частоте среза ω_ср значение модуля частотной характеристики W_раз (jω) равно единице - , логарифм от такого модуля равен нулю - и, следовательно, значение ЛАЧХ от частотной характеристики W_раз (jω) на частоте среза также равно нулю - ;

3) запас по фазе γ и значение фазочастотной характеристики, вычисленное на частоте среза удовлетворяют равенству

, (1.53)

и, следовательно, запас по фазе определяется по выражению фазочастотной характеристики φ(ω), вычисленное при одной частоте – частоте среза ω_ср.

Описанные свойства построений, приведенных на рис.1.40, позволяют перейти к определению частоты среза из графика ЛАЧХ разомкнутой САУ (рис.1.41).

Частота среза определяется графически как точка пересечения ЛАЧХ L(ω) с осью частоты ω. Далее рассчитывается только одно значение фазы φ(ω_сз) и по выражению (1.53) определяется запас по фазе γ.

Простота построения ЛАЧХ (см. тему 1.4), а также простота вычисления запаса по фазе, обусловливает простое использование критерия устойчивости Найквиста в следующей формулировке:

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

Вопросы и задания

1. Сформулируйте критерий устойчивости Найквиста для статической САУ.

2. Сформулируйте критерий устойчивости Найквиста для астатической САУ.

3. Какие косвенные показатели качества вводятся на построениях Найквиста?

4. Сформулируйте критерий устойчивости Найквиста, основанный на построениях ЛАЧХ.

1.12. Прямые показатели качества САУ. Расчет ошибок

регулирования. Статические и астатические САУ

Прямые показатели качества подразделяются на показатели качества динамического и установившегося режимов.

Показателями качества динамических режимов определяются из графика переходного процесса и основными из них являются (рис.1.42):

- перерегулирование или заброс σ, равный максимуму отклонения значения переходного процесса относительно установившегося значения процесса h_ycm;

- время первой установки t₁, определяемое моментом первого пересечения графиком переходного процесса установившегося значения h_ycm;

- время переходного процесса t_ПП, определяемое момент окончательного входа графика переходного процесса в зону допуска, равную ±5% от установившегося значения процесса h_ycm.

Для всех названных динамических показателей качества невозможно в общем случае получить формулы для их расчета. Это является существенным препятствием для решения задач анализа и синтеза САУ.

Показателями качества установившихся режимов являются ошибки регулирования, равные абсолютной величине разности между заданным и фактическим значениями сигналов САУ и которые в зависимости от вида входного сигнала САУ подразделяются на статические (ε_СТ) и скоростные ошибки (ε_СК) и ошибки (ε_m) при отработке гармонического входного сигнала.

Для всех названных ошибок регулирования можно в общем случае получить формулы их расчета.

Из структурной схемы замкнутой САУ (рис.1.43) следуют выражения передаточной функции САУ W_ε(p) по ошибке и изображения ε(р) ошибки регулирования:

Расчет ошибки ε_m отработки гармонического входного сигнала x=X_msinωt производится по формуле

где - модуль комплексного числа .

Статическая (ε_СТ) и скоростная (ε_СК) ошибки равны установившимся значениям оригинала и , или в общем виде, по формуле . Значение вычисляют через изображение ε (р) по доказываемой в теории операционного исчисления формуле предельного перехода,

(1.54)

Выражение передаточной функции разомкнутой САУ в общем случае может быть приведено к виду:

(1.55)

где К – общий коэффициент усиления разомкнутой САУ:

ν - порядок астатизма САУ, причем ν является целым неотрицательным числом.

Для удобства вычислений по формуле (1.54) подставим в нее выражение W_РАЗ(р) из (1.55) и выполним предельный переход:

(1.56)

Статическая ошибка регулирования ε_СТ рассчитывается при постоянном входном сигнале x (t) =X=const, а скоростная ε_СК - при входном сигнале x=Vt, изменяющемуся во времени с постоянной скоростью V=const. Далее расчеты статической (ε_СТ) и скоростной (ε_СК) ошибок выполним раздельно.

Расчеты статической ошибки ε_СТ регулирования

Входной сигнал x(t)=X=const и изображением его является . В соответствии с (1.56) статическую ошибку ε_СТ следует вычислять по формуле

(1.57)

1). Пусть в (1.57) значение порядка ν астатизма САУ равно нулю: ν=0. Такая САУ называется статической. Тогда статическая ошибка ε_СТ будет равна

В статической САУ имеется статическая ошибка ε_СТ, которую можно только уменьшить путем увеличения общего коэффициента усиления К разомкнутой САУ, но обратить в ноль ее нельзя.

2). Пусть в (1.57) значение порядка ν астатизма САУ равно 1: ν=1. Такая САУ называется астатической 1-го порядка. Тогда статическая ошибка ε_СТ будет равна

В астатической САУ 1-го порядка статическая ошибка ε_СТ равна нулю, т.е САУ является абсолютно точной. Можно проверить, что при астатизме САУ выше 1, статическая ошибка регулирования всегда будет нулевой.

Расчеты скоростной ошибки ε_СТ регулирования

Входной сигнал x(t)=Vt и изображением его является . В соответствии с (1.56) скоростную ошибку ε_СК следует вычислять по формуле

(1.58)

1). Пусть в (1.58) значение порядка ν астатизма САУ равно нулю: ν=0. Такая САУ называется статической. Тогда скоростная ошибка ε_СК будет равна

В статической САУ скоростная ошибка ε_СК бесконечно большая и, поэтому, такая САУ неработоспособна.

2). Пусть в (1.58) значение порядка ν астатизма САУ равно 1: ν=1. Такая САУ называется астатической 1-го порядка. Тогда скоростная ошибка ε_СК будет равна

В астатической САУ 1-го порядка имеется скоростная ошибка ε_СК, которую можно только уменьшить путем увеличения общего коэффициента усиления К разомкнутой САУ, но обратить в ноль ее нельзя.

3). Пусть в (1.58) значение порядка ν астатизма САУ равно 2: ν=2. Такая САУ называется астатической 2-го порядка. Тогда скоростная ошибка ε_СК будет равна

В астатической САУ 2-го порядка скоростная ошибка ε_СК равна нулю, т.е САУ является абсолютно точной.

Выводы по расчетам статической и скоростной ошибок регулирования:

1. Ошибки регулирования могут быть уменьшены путем увеличения общего коэффициента усиления К и порядка астатизма ν разомкнутой САУ.

2. При увеличении К ошибки регулирования только уменьшаются. но не обращаются в ноль.

3. При увеличении ν САУ становится абсолютно точной - ошибка регулирования становится нулевой.

Вопросы и задания

1. Перечислите основные прямые показатели, которыми оценивается качество работы САУ.

2. Приведите вывод выражения установившейся ошибки регулирования.

3. Приведите расчеты статической ошибки регулирования для статической и астатической САУ.

4. Приведите расчеты скоростной ошибки регулирования для статической и астатической САУ.

5. Какие существуют универсальные способы уменьшения ошибок регулирования?

1.13. Косвенные показатели качества САУ и их связь

с прямыми показателями качества. Использование

ЛАЧХ для оценки качества САУ

Невозможность получения формул для расчета динамических показателей качества (рис.1.42), а также требования задач синтеза САУ, обусловило разработки комплексных показателей качества. Косвенные показатели качества, в большинстве своем, являются частотными, которые определяются из ЧХ, АЧХ, ФЧХ и ЛАЧХ. Косвенные показатели качества должны удовлетворять следующим требованиям:

1. Косвенные показатели должны просто вычисляться или определяться из частотных характеристик разомкнутой САУ.

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

3. Косвенные показатели должны быть приспособлены для эффективного решения задач синтеза САУ.

4. Косвенные показатели должны давать возможность просто анализировать влияние параметров настроек регуляторов САУ и характеристик любых других звеньев САУ на прямые показатели качества.

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

Прежде всего рассмотрим набор косвенных показателей качества, полученных из построений Найквиста (см. тему 1.12): частоту среза ω_СР и запас по фазе γ. Частота среза ω_СР просто определяется из ЛАЧХ (рис.1.41). Запас по фазе γ рассчитывается по выражению ФЧХ φ (ω) только при одном значении частоты ω_СР: γ=φ (ω_СР).

Основой применения косвенных показателей качества - частоты среза ω_СР и запаса по фазе γ - являются графические зависимости (рис.14.1) между косвенными и прямыми показателями качества - перерегулированием σ, временем первой установки t₁ и временем переходного процесса t_ПП.

По правой оси ординат отложены значения перерегулирования σ, в процентах от установившегося значения h_ycm (рис.1.42). По левой оси ординат записаны формулы, по которым рассчитываются t₁ и t_ПП в зависимости от частоты среза ω_СР. Если из частотных характеристик определены значения запаса по фазе γ и частоты среза ω_СР, то по графикам можно определить значения перерегулирования σ, времени первой установки t₁ и времени переходного процесса t_ПП. Например, пусть заданы значения γ=30^о и ω_СР=1,5 с^-1. Тогда, согласно приведенным на рис.1.44 построениям, получим:

Поиск по сайту: