9138

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

Контрольная

Математика и математический анализ

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

Русский

2013-02-24

1.49 MB

68 чел.

Устойчивость САУ

Нули и полюсы передаточной функции

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

Если корни числителя и знаменателя передаточной функции лежат в левой полуплоскости (при этом корни числителя и знаменателя лежат в верхней полуплоскости), то звено называется минимально-фазовым.

Соответствие левой полуплоскости корней р верхней полуплоскости корней   (рис.2.2.1) объясняется тем, что   , или , т.е. вектор получается из вектора поворотом на угол по часовой стрелке. В результате все векторы из левой полуплоскости приходят в векторы в верхней полуплоскости.

Неминимально-фазовые  и  неустойчивые  звенья

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

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

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

Существуют звенья, у которых этот процесс выражен еще заметнее. Это объясняется наличием положительных вещественных или комплексных корней с положительной вещественной частью в характеристическом уравнении (знаменателе передаточной функции, приравненом нулю), в результате чего звено будет относиться к категории неустойчивых звеньев.

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

   ,

т.е. большее по абсолютной величине значение.

В связи с этим неустойчивые звенья относятся к группе не минимально-фазовых звеньев.

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

Например, звено с передаточной функцией    относится к группе не минимально–фазовых звеньев. Модуль частотной передаточной функции совпадает с модулем частотной передаточной функции звена, имеющего передаточную функцию   . Но фазовый сдвиг первого звена по абсолютной величине  больше:

                              .

Минимально-фазовые звенья имеют меньшие фазовые сдвиги по сравнению с соответствующими звеньями, имеющими такие же амплитудные частотные характеристики.

Говорят, что система устойчива или обладает самовыравниванием, если после снятия внешнего возмущения она возвращается в исходное состояние.

Так как движение системы в свободном состоянии описывается однородным дифференциальным уравнением, то математическое определение устойчивой системы можно cфоpмулировать следующим образом:

Система называется асимптотически устойчивой, если выполняется условие                                                                 (2.9.1)

Из анализа общего решения (1.2.10) вытекает необходимое и достаточное условие устойчивости:

Для устойчивости системы необходимо и достаточно, чтобы все корни характеристического уравнения имели строго отрицательные вещественные части, т.е.                                              Repi < 0,     I = 1…n.                                    (2.9.2)

Для наглядности корни характеристического уравнения принято изображать на комплексной плоскости рис.2.9.1а. При выполнении необходимого и достаточного условия (2.9.2) все корни лежат слева от мнимой оси, т.е. в области устойчивости.

Поэтому условие (2.9.2) можно сформулировать следующим образом.

Для устойчивости необходимо и достаточно, чтобы все корни характеристического уравнения располагались в левой полуплоскости.

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

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

Условие  Стодолы

Условие получено словацким математиком Стодолой в конце 19-го столетия. Оно интересно в методическом плане для понимания условий устойчивости системы.

Запишем характеристическое уравнение системы в виде

                                              D(p)  = a0pn + a1pn-1  +…an= 0.                   (2.9.3)  

По Стодоле для устойчивости необходимо, но недостаточно, чтобы пpи    a0 > 0 все остальные коэффициенты были строго положительны, т.е.

 a1 > 0,..., an > 0.

Необходимость можно сформировать так:

Если система устойчива, то все корни характеристического уравнения имеют , т.е. являются левыми.

Доказательство необходимости элементарное. По теореме Безу характеристический полином можно представить в виде

                             .

Пусть , т.е действительное число, а  – комплексно-сопряженные корни. Тогда

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

Недостаточность условия Стодолы заключается в том, что условие  не гарантирует, что все . В этом можно убедиться на конкретном примере, рассмотрев полином степени .

Заметим, что в случае условие Стодолы одновременно необходимо и достаточно. Из вытекает . Если , то и , чтобы .

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

Из условия Стодолы вытекает два важных следствия.

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

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

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

Это подтверждает и следующий пример.

Пример 2.9.1. Применить условие Стодолы к схеме рис. 2.9.2.

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

   D(p ) =  p2 + k1k2  =  0.

Так как отсутствует член с р в первой степени (a1 = 0), то условие Стодолы не выполняется и система неустойчива. Данная система структурно неустойчива, так как ни при каких значениях параметров k1 и  k2 не может быть      устойчивой.

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

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

                     ,

                     .

Теперь условие Стодолы выполняется при любых . Так как в случае уравнения второй степени оно не только необходимо, но и достаточно, то система устойчива при любых положительных коэффициентах усиления .

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

                             .

Аналогично предыдущему система устойчива при любых положительных .

Критерий  устойчивости  Раусса-Гурвица

Математики  Раусс (Англия) и Гурвиц (Швейцария) разработали этот критерий приблизительно в одно время. Отличие заключалось в алгоритме вычислений. Мы познакомимся с критерием в формулировке Гурвица.

По Гурвицу для устойчивости необходимо и достаточно, чтобы при a0 > 0 определитель Гурвица  = n и все его главные миноры 1, 2,..., n-1 были строго положительны, т.е.

                                                                   (2.9.4)

                   

Cтруктура определителя Гурвица легко запоминается, если учесть, что по главной диагонали расположены коэффициенты а1,…n, в строчках расположены коэффициенты через один, если они исчерпаны, то свободные места заполняются нулями.

Пример 2.9.2. Исследовать на устойчивость по Гурвицу систему с единичной отрицательной обратной связью, в прямой цепи которой включены три инерционных звена и, следовательно, передаточная функция разомкнутой системы имеет вид                          (2.9.5)

Запишем характеристическое уравнение замкнутой системы как сумму числителя и знаменателя (2.9.5):

Следовательно,

Определитель Гурвица и его миноры имеют вид

     (2.9.6)

с учетом a0 > 0 из строгой положительности определителя Гурвица и миноров (2.9.6) вытекает условие Стодолы и, кроме того, условие a1 a2 - a0 a3  > 0, что после подстановки значений коэффициентов дает          

                              1Т2+ Т1Т32Т3)(Т123)>Т1Т2Т3(1+k).                (2.9.7)

Отсюда видно, что при увеличении k система из устойчивой может превратиться в неустойчивую, так как неравенство (2.9.7)  перестанет выполняться.

     Передаточная функция системы по ошибке равна

     (2.9.8)

Согласно теореме о конечном значении оригинала установившаяся ошибка отработки единичного ступенчатого сигнала будет равна 1/(1+k).  Следовательно, обнаруживается противоречие между устойчивостью и точностью. Для уменьшения ошибки надо увеличивать k, но это приводит к потере устойчивости.

Принцип  аргумента  и  критерий  устойчивости  Михайлова

Критерий Михайлова основан на так называемом принципе аргумента.

Рассмотрим характеристический полином замкнутой системы, который по теореме Безу можно представить в виде

                  D(p) = a0pn+ a1pn-1+…+ an = a0(p - p1)…(p - pn).

Сделаем подстановку p = j

D(j) = a0(j)n+ a1(j)n-1+…+ an = a0(j - p1)…(j - pn) = X()+jY().

Для конкретного значения имеет точку на комплексной плоскости, задаваемую параметрическими уравнениями  

Если изменять в диапазоне от - до , то будет прочерчена кривая Михайлова, т. е. годограф. Изучим поворот вектора D(j) при изменении от - до , т. е. найдем приращение аргумента вектора (аргумент равен сумме для произведения векторов): .

При  = - разностный вектор, начало которого в точке рi, а конец на мнимой оси, направлен вертикально вниз. По мере роста конец вектора скользит вдоль мнимой оси, а при  = вектор направлен вертикально вверх. Если корень левый (рис. 2.9.19а), то arg = +, а если корень правый, то arg = -.

Если характеристическое уравнение имеет m правых корней (соответственно n - m левых), то          .                               

Это и есть принцип аргумента. При выделении действительной части Х() и мнимой Y() мы отнесли к Х() все слагаемые, содержащие j в четной степени, а к Y() – в нечетной степени. Поэтому кривая Михайлова симметрична относительно действительной оси (Х() – четная, Y() – нечетная функция). В результате, если изменять от 0 до +, то приращение аргумента будет в два раза меньше. В связи с этим окончательно принцип аргумента формулируется следующим образом        .         (2.9.29)  

Если система устойчива, т.е. m = 0, то получаем критерий устойчивости Михайлова.

По Михайлову для устойчивости необходимо и достаточно, чтобы

                                              ,                                 (2.9.30)

то есть кривая Михайлова должна последовательно проходить через n четвертей против часовой стрелки.

Очевидно, что для применения критерия Михайлова не требуется точного и детального построения кривой. Важно установить, каким образом она огибает начало координат и не нарушается ли последовательность прохождения n четвертей против часовой стрелки.

Пример 2.9.6. Применить критерий Михайлова для проверки устойчивости системы, показанной на рис.2.9.20.

Характеристический полином замкнутой системы при k1k2 > 0 соответствует устойчивой системе, так условие Стодолы   выполняется, а для n = 1 оно достаточно. Можно непосредственно найти корень р1 = - k1k2 и убедиться, что необходимое и достаточное условие устойчивости выполнено. Поэтому применение критерия Михайлова носит иллюстративный характер. Полагая p=j, получим

                                        D(j) = X()+jY(),

где                    Х() = ;   Y() = .       (2.9.31)

По параметрическим уравнениям (2.9.31) построен годограф Михайлова на рис.2.9.21, из которого видно, что при изменении от 0 до вектор D(j) поворачивается против часовой стрелки на +/2 , т.е. система устойчива.

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

Как уже было отмечено, критерий Найквиста занимает особое положение среди критериев устойчивости. Это частотный критерий, позволяющий определить устойчивость замкнутой системы по частотным характеристикам разомкнутой. При этом предполагается, что система разомкнута по цепи единичной отрицательной обратной связи (рис.2.9.22).          

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

Вывод критерия основан на использовании принципа аргумента. Передаточная функция разомкнутой системы (по цепи единичной отрицательной обратной связи на рис.2.9.22) равна   

Рассмотрим        .       (2.9.32)

В случае реальной системы с ограниченной полосой пропускания степень знаменателя передаточной функции разомкнутой системы п  больше степени числителя   , т.е. n > . Поэтому степени характеристических полиномов разомкнутой системы и замкнутой системы одинаковы и равны n. Переход от АФХ разомкнутой системы  к АФХ по (2.9.32)  означает увеличение вещественной части на 1, т.е. перенос начала координат в точку (-1, 0), как показано на рис.2.9.23.

Предположим теперь, что замкнутая система устойчива, а характеристическое уравнение разомкнутой системы А(р ) = 0 имеет m правых корней. Тогда в соответствии с принципом аргумента (2.9.29) получим необходимое и достаточное условие устойчивости замкнутой системы по Найквисту

   (2.9.33)

Т.е. для устойчивости замкнутой системы вектор W1(j) должен делать m/2 полных оборотов против часовой стрелки, что равносильно повороту вектора Wpaз(j) относительно критической точки (-1,0).

      На практике, как правило, разомкнутая система устойчива, т.е. m = 0. В этом случае приращение аргумента равно нулю, т.е. АФХ разомкнутой системы не должна охватывать критическую точку (-1,0).

Критерий  Найквиста  для  ЛАХ  и  ЛФХ

На практике чаще используются логарифмические характеристики разомкнутой системы. Поэтому целесообразно сформулировать критерий Найквиста для определения устойчивости замкнутой системы по ним. Количество оборотов АФХ относительно критической точки  (-1,0) и охват или не охват ее

зависят от количества положительных и отрицательных пересечений интервала (-,-1) действительной оси и соответственно пересечений фазовой характеристикой линии -180° в области L()  0 . На рис.2.9.24 изображены АФХ и показаны знаки пересечений отрезка (-,-1) действительной оси.

Справедливо правило

                                 

где                      -  число положительных и отрицательных пересечений.

По АФХ  рис.2.9.24в построены ЛАХ и ЛФХ, изображенные на рис.2.9.25, причем на ЛФХ отмечены положительные и отрицательные пересечения. На отрезке    (-,-1) модуль больше единицы, чему соответствует L() > 0. Поэтому Критерий Найквиста:

Для устойчивости замкнутой системы ЛФХ разомкнутой системы в области, где L() > 0, должна иметь положительных пересечений линии -180° на больше, чем отрицательных.

Если разомкнутая система устойчива, то число положительных и отрицательных пересечений фазовой характеристикой линии -180° в области L() > 0 для устойчивости замкнутой системы должно быть одинаковым или пересечений не должно быть.

Критерий  Найквиста  для  астатической  системы

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

                           .

В этом случае  при 0, т. е. амплитудно-фазовая характеристика (АФХ) разомкнутой системы уходит в бесконечность. Раньше мы строили АФХ при изменении от - до и это была непрерывная кривая, замкнутая при  = 0. Теперь она также замыкается при  = 0, но на бесконечности и при этом не ясно, с какой стороны действительной оси (на бесконечности слева или справа?).

Рис.2.9.19в иллюстрирует, что в этом случае возникает неопределенность в подсчете приращения аргумента разностного вектора. Он теперь все время расположен вдоль мнимой оси  (совпадает с  j). Только при переходе через нуль изменяется направление (при этом поворот вектора против часовой стрелки на  или по часовой стрелке на - ?), Для определенности считаем условно, что корень левый и огибание начала координат происходит по дуге бесконечно малого радиуса против часовой стрелки (поворот на +). Соответственно  в окрестности  = 0  представим в виде

                  ,

где = + при изменении  от – 0 до + 0. Последнее выражение показывает, что при таком раскрытии неопределенности АФХ поворачивается при изменении  от – 0 до + 0 на угол - по часовой стрелке. Соответственно построенную АФХ надо при  = 0 дополнить дугой бесконечности радиуса на угол , т. е. против часовой стрелки до положительной действительной полуоси.

Запасы  устойчивости  по  модулю  и  фазе

Чтобы гарантировать устойчивость при изменениях параметров системы вводятся запасы устойчивости по модулю и фазе, определяемые следующим образом.

Запас устойчивости по модулю показывает во сколько раз или на сколько децибел допустимо увеличивать или уменьшать коэффициент усиления, чтобы система оставалась устойчивой (оказывалась на границе устойчивости). Он определяется как min(L3,L4) на рис.2.9.25. Действительно, если не менять ЛФХ, то при подъеме ЛАХ на L4 частота среза ср переместится в точку 4 и система окажется на границе устойчивости. Если опустить ЛАХ на L3, то частота среза сместится влево в точку 3 и система также окажется на границе устойчивости. Если опустить ЛАХ еще ниже, то в области L() > 0 останется только отрицательное пересечение ЛФХ линии -180°, т.е. по критерию Найквиста система станет неустойчивой.

Запас устойчивости по фазе показывает, на сколько допустимо увеличить фазовый сдвиг при неизменном коэффициенте усиления, чтобы система оставалась устойчивой (оказалась на границе устойчивости). Он определяется как дополнение (ср) до -180°.

На практике  L  12-20 дБ,      20-30°.

Рис.2.2.1. Расположение корней минимально-фазового звена

Рис.8.12. Плоскость корней

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

уравнения A(p) = 0

ОУ- область устйчивости

Рис.2.9.1б. Плоскость корней

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

= 0. ОУ- область

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

Рис.2.9.1а. Плоскость корней

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

A(p) = 0. ОУ- область

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

Рис.2.9.2.  Структурная

схема  системы

Рис.2.9.3. Схема с дополнительной связью

Рис.2.9.4. Схема с корректирующим

звеном

ис.2.9.19. Определение приращения аргумента вектора:

а) случай левого корня; б) случай правого корня;

в) случай чисто мнимого корня

Рис.2.9.20.  Структурная

схема  системы

Рис.2.9.21. Годограф Михайлова для

системы рис.2.9.20

Рис.2.9.22. Место размыкания системы при проверке устойчивости

по критерию Найквиста

Рис.2.9.23. АФХ для определения

устойчивости по Найквисту

Рис.2.9.24. Примеры АФХ

Рис.2.9.25. ЛАХ и ЛФХ для

АФХ на рис.2.9.24в


 

А также другие работы, которые могут Вас заинтересовать

39953. Течение газа в соплах 182.5 KB
  В рамках этой модели течения невязкий газ и пограничный слой при отсутствии отрыва потока представляется возможным с достаточной точностью определить оптимальное сопло для заданных конструктивных условий габариты масса тяга. Основные недостатки сопел Лаваля связанные с их большой длинной массой и низкой эффективностью при перерасширении потока становятся особенно ощутимыми при больших степенях расширения сопла в этом случае размеры и масса сопла могут быть на порядок больше размеров и массы камеры сгорания а потери тяги...
39954. Одномерные течения несжимаемой жидкости. Ламинарное и турбулентное течения 344.5 KB
  При увеличении скорости воды картина изменялась струйка красителя сначала приобретала синусоидальную форму а дальнейшее увеличение скорости приводило к ее размыву что свидетельствовало о беспорядочном движении. Рейнольдс предположил что увеличение скорости потока приводит к возникновению какихто возмущений дестабилизирующих его структуру. Ускорение есть изменение скорости в единицу времени = u t. Одномерными называются течения в которых основные параметры потока зависят лишь от одной координаты направление которой совпадает с...
39955. Основы теории пограничного слоя 73.5 KB
  Основы теории пограничного слоя. Понятие пограничного слоя 8. Толщина пограничного слоя 8. Отрыв пограничного слоя.
39956. Основы теории подобия 362.5 KB
  Основы теории подобия План. На эти вопросы и отвечает теория подобия являющаяся основой современного физического эксперимента. В общем случае различают три вида подобия: геометрическое кинематическое и динамическое. Для площадей S и объемов V ; Применительно к физическим явлениям элементарные представления геометрического подобия расширяются и распространяются на все величины характеризующие данный процесс.
39957. Газодинамика как раздел механики сплошных сред 907.5 KB
  Краткий очерк развития механики жидкости и газа. Математический аппарат используемый в механике жидкости и газа [1. Газодинамика как раздел механики сплошных сред Многие машины и аппараты созданные к настоящему времени характеризуются перемещением газа или жидкости внутри их или перемещением самого аппарата в среде газа или жидкости. Целью курса Газодинамика является изучение явлений протекающих в газе и жидкости и закономерностей которым эти явления подчиняются.
39958. УРАВНЕНИЯ ГАЗОВОЙ ДИНАМИКИ ДЛЯ ЕДИНИЧНОЙ СТРУЙКИ 401.5 KB
  Предельная скорость движения газа. Уравнение неразрывности Выведем основные уравнения газовой динамики для элементарной струйки газа поперечные размеры которой настолько малы что в каждом ее сечении можно считать постоянными все основные параметры потока: скорость давление температуру и плотность газа. Чтобы получить уравнение неразрывности рассмотрим стационарное установившееся движение элементарной струйки газа рис. Элементарная струйка Рассмотрим некоторый участок струйки между двумя нормальными к поверхности тока сечениями 1 и...
39959. Элементы гидродинамики 441 KB
  Cилы действующие в жидкости 3.1 – Элементарный параллелепипед в потоке жидкости Грани бесконечно малой частицы жидкости имеющей в начале движения форму прямого параллелепипеда с ребрами dx dy dz с течением времени могут скашиваться и растягиваться рис.8 представляет собой уравнение неразрывности жидкости.9 Здесь под плотностью жидкости понимается предел отношения массы частицы к ее объему 3.
39960. ГАЗОДИНАМИЧЕСКИЕ ФУНКЦИИ 81 KB
  ГАЗОДИНАМИЧЕСКИЕ ФУНКЦИИ План лекции. Зависимость параметров потока в функции числа M. Зависимость параметров потока в функции скоростного коэффициента. Зависимость параметров потока в функции числа M.
39961. ДЕТАЛИ МАШИН И ОСНОВЫ КОНСТРУИРОВАНИЯ 10.06 MB
  1 а е: Ft Н окружная сила на барабане ленточного или на звездочке цепного конвейера; V м с скорость движения ленты или цепи; Dб мм диаметр барабана; Zзв число зубьев тяговой звездочки; Рзв мм шаг тяговой цепи.2 Вид передачи Твердость зубьев Передаточное число Uрек Uпред Зубчатая цилиндрическая: тихоходная ступень во всех редукторах uт 350 НВ 40. Термообработка зубчатых колес редуктора улучшение твердость зубьев 350НВ. Первая группа колеса с твердостью поверхностей зубьев Н  350 НВ Применяются в слабо и...