Энциклопедия АСУ ТП Спонсор проекта: Skip Navigation LinksЭнциклопедия АСУ ТП : 5 ПИД-регуляторы : 5.4 Особенности реальных регуляторов : 5.4.3 Запас устойчивости и робастность Соспонсор:




Робот BotEyes




Промышленные контроллеры RealLab!

5.4.3. Запас устойчивости и робастность

Возможность потери устойчивости является основным недостатком систем с обратной связью. Поэтому обеспечение необходимого запаса устойчивости являются самым важным этапом при разработке и настройке ПИД-регулятора.

Устойчивость системы с ПИД-регулятором - это способность системы возвращаться к слежению за уставкой после прекращения действия внешних воздействий. В контексте данного определения под внешними воздействиями понимаются не только внешние возмущения, действующие на объект, но любые возмущения, действующие на любую часть замкнутой системы, в том числе шумы измерений, временная нестабильность уставки, шумы дискретизации и квантования, шумы и погрешность вычислений. Все эти возмущения вызывают отклонения системы от положения равновесия. Если после прекращения воздействия система возвращается в положение равновесия, то она считается устойчивой. При анализе устойчивости ПИД-регуляторов обычно ограничиваются исследованием реакции системы на ступенчатое изменение уставки , шум измерений и внешние возмущения . Потеря устойчивости проявляется как неограниченное возрастание управляемой переменной объекта, или как ее колебание с нарастающей амплитудой.

Рис. 5.77. Структура разомкнутой системы управления с ПИД-регулятором для анализа устойчивости

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

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

Наиболее полную информацию о запасе устойчивости системы можно получить, решив дифференциальное уравнение, описывающее замкнутую систему при внешних возмущениях. Однако этот процесс слишком трудоемок, поэтому для линейных систем используют упрощенные методы, позволяющие дать оценку запаса устойчивости без решения уравнений [Воронов]. Мы рассмотрим два метода: оценку с помощью годографа комплексной частотной характеристики разомкнутого контура (критерий Найквиста) и с помощью логарифмических АЧХ и ФЧХ (диаграмм Боде).

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

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

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

Рассмотрим систему, состоящую из контроллера и объекта управления (рис. 5.77), которая получена путем исключения из классической системы с ПИД-регулятором (рис. 5.34) сигнала уставки . Будем считать, что обратная связь разомкнута, а для ее замыкания достаточно соединить точки и . Предположим теперь, что на вход подан сигнал

.

(5.86)

Тогда, пройдя через регулятор и объект управления, этот сигнал появится на выходе с измененной амплитудой и фазой, в виде

,

(5.87)

где - комплексная частотная характеристика (КЧХ) системы, - аргумент КЧХ, - модуль КЧХ. Таким образом, при прохождении через регулятор и объект амплитуда сигнала изменится пропорционально модулю, а фаза - на величину аргумента КЧХ.

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

, , т. е. .

(5. 88)

Рис. 5.78. Три годографа КЧХ разомкнутых систем для объекта второго порядка при , и пропорциональном коэффициенте регулятора

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

Комплексную частотную характеристику графически изображают в виде годографа (диаграммы Найквиста) - графика в координатах и (рис. 5.78). Стрелка на линии годографа указывает направление движения "карандаша" при возрастании частоты. Точка , которая соответствует условию существования незатухающих колебаний в системе, на этом графике имеет координаты и . Поэтому критерий устойчивости Найквиста формулируется следующим образом [Ротач]: контур, устойчивый в разомкнутом состоянии, сохранит устойчивость и после его замыкания, если его КЧХ в разомкнутом состоянии не охватывает точку с координатами [-1, j0]. Более строго, при движении вдоль траектории годографа в направлении увеличения частоты точка [-1, j0] должна оставаться слева [Astrom], чтобы замкнутый контур был устойчив.

На рис. 5.79 показаны реакции замкнутых систем с тремя различными годографами (рис. 5.78) на единичный скачок уставки. Все три системы устойчивы, однако скорость затухания колебаний и форма переходного процесса у них различная. Интуитивно понятно, что система с параметрами наиболее близка к тому, чтобы перейти в состояние незатухающих колебаний при небольшом изменении ее параметров. Поэтому при проектировании ПИД-регулятора важно обеспечить не столько устойчивость, сколько ее запас, необходимый для нормального функционирования системы в реальных условиях.

Запас устойчивости оценивают как степень удаленности КЧХ от критической точки [-1, j0]. Если , то можно найти, во сколько раз осталось увеличить передаточную функцию, чтобы результирующее усиление вывело систему в колебательный режим: , откуда

.

(5. 89)

Запасом по усилению называется величина, на которую нужно умножить передаточную функцию разомкнутой системы , чтобы ее модуль на частоте сдвига фаз 180˚ стал равен 1.

Если на частоте сдвига фаз 180˚ коэффициент усиления разомкнутого контура равен (рис. 5.78), то дополнительное усиление величиной переведет систему в точку [-1, j0], поскольку .

Аналогично вводится понятие запаса по фазе: это минимальная величина , на которую нужно увеличить фазовый сдвиг в разомкнутой системе , чтобы суммарный фазовый сдвиг достиг 180˚, т.е.

.

(5.90)

Рис. 5.79. Переходная характеристика замкнутой системы, которая имеет годограф, показанный на рис. 5.78

Знак "+" перед стоит потому, что .

Для оценки запаса устойчивости используют также минимальное расстояние от кривой годографа до точки [-1, j0] (рис. 5.78).

На практике считаются приемлемыми значения =2...5, =30˚...60˚, =0,5...0,8 [Astrom].

Для графика на рис. 5.78 эти критерии имеют следующие значения:

  • для случая , =12,1; =15˚; =0,303.
  • для случая , =11,8; =47,6; =0,663.
  • для случая , =1,5; =35,2; =0,251.

Если кривая годографа пересекает действительную ось в нескольких точках, то для оценки запаса устойчивости берут ту из них, которая наиболее близка к точке [-1, j0]. При более сложном годографе может быть использована оценка запаса устойчивости как запас по задержке [Astrom]. Запас по задержке- это минимальная задержка, при добавлении которой в контур он теряет устойчивость. Наиболее часто этот критерий используется для оценки запаса устойчивости систем с предиктором Смита.

Частотный критерий устойчивости

Рис. 5.80. Оценка запаса по фазе и усилению для системы с годографом, показанным на рис. 5.78

Для графического представления передаточной функции разомкнутой системы и оценки запаса устойчивости могут быть использованы логарифмические АЧХ и ФЧХ (рис. 5.80). Для оценки запаса по фазе сначала с помощью АЧХ находят частоту ("частота среза" или "частота единичного усиления"), при которой , затем по ФЧХ находят соответствующий запас по фазе. Для оценки запаса по усилению сначала с помощью ФЧХ находят частоту , на которой фазовый сдвиг равен 180˚, затем по АЧХ находят запас по усилению. На (рис. 5.80) приведены примеры графических построений для оценки запаса по фазе и усилению для системы, годограф которой показан на рис. 5.78.

Если запас по фазе разомкнутого контура равен 0˚ или запас по усилению равен 1, после замыкания контура обратной связи система окажется неустойчивой.

Функции чувствительности

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

.

(5. 91)

Поделив обе части этого равенства на и подставив в правую часть , получим

.

(5.92)

Из последнего соотношения виден смысл коэффициента - он характеризует степень влияния относительного изменения передаточной функции объекта на относительное изменение передаточной функции замкнутого контура, то есть является коэффициентом чувствительности замкнутого контура к вариации передаточной функции объекта. Поскольку коэффициент является частотозависимым, его называют функцией чувствительности [Astrom].

Как следует из (5.92),

.

(5. 93)

Введем обозначение .

(5.94)

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

Величина называется комплементарной (дополнительной) функцией чувствительности [Astrom], поскольку .

Заметим, что функция чувствительности входит во все три слагаемые уравнения замкнутой системы с ПИД-регулятором (5.42).

Функция чувствительности позволяет оценить изменение свойств системы после замыкания обратной связи. Поскольку передаточная функция разомкнутой системы равна , а замкнутой , то их отношение . Аналогично, передаточная функция от входа возмущений на выход замкнутой системы равна (см. (5.42)), а разомкнутой - , следовательно, их отношение также равно . Для передаточной функции от входа шума измерений на выход системы можно получить то же отношение .

Таким образом, зная вид функции (см, например, рис. 5.81), можно сказать, как изменится подавление внешних воздействий на систему для разных частот после замыкания цепи обратной связи. Очевидно, шумы, лежащие в диапазоне частот, в котором , после замыкания обратной связи будут усиливаться, а шумы с частотами, на которых , после замыкания обратной связи будут ослаблены.

Наихудший случай (наибольшее усиление внешних воздействий) будет наблюдаться на частоте максимума модуля функции чувствительности (рис. 5.81):

.

(5.95)

Максимум функции чувствительности можно связать с запасом устойчивости (рис. 5.78). Для этого обратим внимание, что представляет собой расстояние от точки (-1, j0) до текущей точки на годографе функции . Следовательно, минимальное расстояние от точки (-1, j0) до функции равно

.

(5.96)

Сопоставляя (5.95) и (5.96), можно заключить, что .

Если с ростом частоты модуль уменьшается, то, как видно из (рис. 5.78), . Подставляя сюда соотношение , получим оценку запаса по усилению, выраженную через максимум функции чувствительности:

.

(5.97)

Аналогично, но с более грубыми допущениями можно записать оценку запаса по фазе через максимум функции чувствительности [Astrom]:

.

(5.98)

Например, при получим и .

Робастность

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

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

Для оценки влияния больших изменения параметров объекта представим передаточную функцию объекта в виде двух слагаемых

,

(5.99)

где - расчетная передаточная функция, - величина отклонения от , которая должна быть устойчивой передаточной функцией. Тогда петлевое усиление разомкнутой системы можно представить в виде . Поскольку расстояние от точки (-1, j0) до текущей точки на годографе невозмущенной системы (для которой ) равно (см. рис. 5.82), условие устойчивости системы с отклонением петлевого усиления можно представить в виде

,

откуда , или ,

где - дополнительная функция чувствительности (5.94). Окончательно можно записать соотношение

,

(5.100)

Рис. 5.82. Пояснение к выводу соотношения (5.100)

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

5.4.4. Сокращение нулей и полюсов

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

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

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

5.4.5. Безударное переключение режимов регулирования

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

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

или .

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

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


© RLDA Ltd. info@rlda.ru  Рейтинг@Mail.ru Спонсоры проекта: , а также