Как найти запасы устойчивости

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

Запас
устойчивости по амплитуде – это минимальный
отрезок действительной оси h, характеризующий
расстояние между критической точкой
(-1, j0) и ближайшей точкой пересечения
годографом Найквиста вещественной оси.

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

Рисунок.
17. Запасы устойчивости САР по критерию
Найквиста

Определение
запасов устойчивости по логарифмическому
критерию.

Запас
устойчивости по амплитуде характеризует
отрезок
при том значении частоты, при котором
ЛФЧХ
(
)= –.
Запасу устойчивости по фазе
соответствует угол между ЛФЧХ и линией
–при частоте среза.

Рисунок18.
Запасы устойчивости САР по логарифмическому
критерию

6. Синтез системы автоматического управления методом логарифмических частотных характеристик.

Meтод
ЛЧХ является одним из наиболее
распространенных методов синтеза
автоматического управления, так как
построение ЛЧХ, как правило, может
выполняться практически без вычислительной
работы. Особенно удобно использовать
асимптотические «идеальные» ЛАЧХ.

Процесс
синтеза обычно включает в себя следующие
операции;

1.
Построение ЛАЧХ неизменяемой части
системы.

Неизменяемая
часть системы регулирования содержит
объект регулирования и исполнительный
элемент, а также основной элемент
обратной связи и элемент сравнения ЛАЧХ
неизменяемой части строят по передаточной
функции разомкнутой неизменяемой части
системы.

2.
Построение желаемой части ЛАЧХ.

График
желаемой ЛАЧХ делается на основе тех
требований, которые предъяв­ляются
к проектируемой системе управления.
Желаемую ЛАЧХ Lж
условно можно разделить на три части:
низкочастотную, среднечастотную и
высокочастотную.

2.1
Низкочастотную часть определяет
статическую точность системы, точность
в установившихся режимах. В статической
системе низкочастотная асимптота
параллельна оси абсцисс. В астатической
системе наклон этой
асим­птоты
составляет –20мдБ/дек, где
–
порядок астатизма (=1,2). Ордината низко­частотной частиLж
определяется значением передаточного
коэффициента К ра­зомкнутой системы.
Чем шире низкочастотная часть Lж,
тем больше высоких частот воспроизводится
системой без замкнутого ослабления.

2.2
Среднечастотная часть является наиболее
важной, так как она определяет
устойчивость, запас устойчивости и,
следовательно, качество переходных
процес­сов, оцениваемое обычно
показателями качества переходной
характеристики. Основные параметры
среднечастотной асимптоты -это её наклон
и частота среза
ср (частота при которойLж
пересекает ось абсцисс). Чем больше
наклон сред­нечастотной асимптоты,
тем труднее обеспечить хорошие
динамические свойст­ва системы.
Поэтому наиболее целесообразен наклон
-20дБ/ дек и крайне редко он превышает
-40дБ /дек Частота срезаср определяет быстродействие системы,
и значение величины перерегулирования.
Чем большеср,
тем выше быстро­действие, тем меньше
время регулирования Тпп переходной
характеристики, тем больше перерегулирование.

2.3
Высокочастотная часть ЛАЧХ незначительно
влияет на динамические свойства системы.
Лучше иметь возможно больше наклон ее
асимптоты, что уменьшает требуемую
мощность исполнительного opгана
и влияние высоко­частотных помех.
Иногда при расчете высокочастотную
ЛАЧХ не принимают во внимание.

При
построении желаемой ЛАЧХ рекомендуется
иметь наклон среднечастотной асимптоты
-20дб/дек, а частоту среза и частоты,
ограничиваю- щие среднечастотную
асимптоту,
иопределить по формулам (2;c
218)

где
–
коэффициент, зависящий от величины
перерегулирования,

,должен
быть выбран по графику приведенном на
рисунке 1.

Рисунок
18- График для определения по допустимому
перерегулированию коэффициента
.

Ордината
низкочастотной асимптоты определяется
соответственно коэффици­

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

3.
Определение параметров корректирующего
устройства.

3.1
График ЛАЧХ корректирующего устройства
получается вычитанием из зна­чения
графика желаемой ЛАЧХ значений графика
неизменяемой, после чего по ЛАЧХ
корректирующего устройства определяется
его передаточная функция [1;с415].

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

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

Пример:

6.Синтез системы
автоматического управления методом
логарифмических частотных характеристик.

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

• #
• #
• #
• #
• #
• #
• #
• #
• #
• #
• #

10.1. Понятие структурной устойчивости. АФЧХ астатических САУ

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

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

Например, если САУ состоит из любого количества инерционных и колебательных
звеньев, она имеет вид, показанный на рис.72

.

При увеличении коэффициента усиления
САУ K каждая точка ее АФЧХ удаляется от начала координат, пока при некотором
значении K_крит АФЧХ не пересечет точку
(-1, j0). При дальнейшем увеличении K, САУ будет неустойчива.
И наоборот, при уменьшении K такую САУ в принципе возможно сделать устойчивой,
поэтому ее называют структурно устойчивой.

Если САУ астатическая, то при ее размыкании характеристическое уравнение можно
представить в виде: pD₁p(p)
= 0, где n – порядок астатизма, равный количеству последовательно
включенных интеграторов. Это уравнение имеет нулевые корни, поэтому при
0, АФЧХ стремится
к (рис.71в и 71г). Например,
пусть W_р(p) = , здесь
= 1, тогда АФЧХ разомкнутой САУ:

W(j)
= = P()
+ jQ().

Так как порядок знаменателя больше порядка числителя, то при
0 имеем P()
–,
Q()
-j. Подобная
АФЧХ представлена на рис.73.

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

Если САУ имеет два интегрирующих звена (порядок астатизма
= 2), ее АФЧХ уходит в бесконечность во втором квадранте (рис.74).

Например, пусть W_р(p)
= , тогда АФЧХ САУ:

W(j)
= = P()
+ jQ().

При
0 имеем P()
–,
Q()
+ j. Такая САУ не
будет устойчива ни при каких значениях параметров, то есть она структурно неустойчива.

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

10.2.
Понятие запаса устойчивости

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

Согласно критерия Найквиста, чем дальше
АФЧХ от критической точки (-1, j0), тем больше запас устойчивости. Различают
запасы устойчивости по модулю и по фазе.

Запас устойчивости по модулю
характеризует удаление годографа АФЧХ разомкнутой САУ от критической точки в направлении
вещественной оси и определяется расстоянием h от критической точки до точки
пересечения годографом оси абсцисс (рис.75).

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

Как уже отмечалось, с ростом коэффициента
передачи разомкнутой САУ растет модуль каждой точки АФЧХ и при некотором значении
K = K_кр АФЧХ пройдет через критическую
точку (рис.76) и попадет на границу устойчивости, а при K > K_кр
замкнутая САУ станет неустойчива. Однако в случае “клювообразных” АФЧХ (получаются
из-за наличия внутренних обратных связей) не только увеличение, но и уменьшение
K может привести к потере устойчивости замкнутых САУ (рис.77). В этом
случае запас устойчивости определяется двумя отрезками h₁
и h₂, заключенными между критической
точкой и АФЧХ.

Обычно при создании САУ задаются требуемыми запасами устойчивости h
и , за пределы
которых она выходить не должна. Эти пределы выставляются в виде сектора, вычерчиваемого
вокруг критической точки, в который АФЧХ разомкнутой САУ входить не должна (рис.78).

10.3.
Анализ устойчивости по ЛЧХ

Оценку устойчивости
по критерию Найквиста удобнее производить по ЛЧХ разомкнутой САУ. Очевидно, что
каждой точке АФЧХ будут соответствовать определенные точки ЛАЧХ и ЛФЧХ.

Пусть
известны частотные характеристики двух разомкнутых САУ (1 и 2), отличающихся друг
от друга только коэффициентом передачи K₁
< K₂. Пусть первая САУ устойчива в замкнутом
состоянии, вторая нет.(рис.79).

Если W₁(p) – передаточная
функция первой САУ, то передаточная функция второй САУ W₂(p)
= KW₁(p),
где K = K₂/K₁.
Вторую САУ можно представить последовательной цепочкой из двух звеньев с передаточными
функциями K (безынерционное звено) и W₁(p),
поэтому результирующие ЛЧХ строятся как сумма ЛЧХ каждого из звеньев.

Поэтому ЛАЧХ второй САУ: L₂()
= 20lgK + L₁(),

а ЛФЧХ: ₂()
=₁().

Пересечениям АФЧХ вещественной оси соответствует значение фазы
= –. Это
соответствует точке пересечения ЛФЧХ
= – линии
координатной сетки. При этом, как видно на АФЧХ, амплитуды A₁()
< 1, A₂()
> 1, что соответствует на САЧХ значениям L₁()
= 20lgA₁()
< 0 и L₂()
> 0.

Сравнивая АФЧХ и ЛФЧХ можно заключить, что система в замкнутом состоянии будет
устойчива, если значению ЛФЧХ
= – будут
соответствовать отрицательные значения ЛАЧХ и наоборот. Запасам устойчивости
по модулю h₁ и h₂,
определенным по АФЧХ соответствуют расстояния от оси абсцисс до ЛАЧХ в точках,
где =
–, но
в логарифмическом масштабе.

Особыми точками являются точки пересечения АФЧХ с единичной окружностью. Частоты
_c1
и _c2,
при которых это происходит называют частотами среза.

В точках пересечения A()
= 1 = > L()
= 0 – ЛАЧХ пересекает горизонтальную ось. Если при частоте среза фаза АФЧХ
_c1
> –
(рис.79а кривая 1), то замкнутая САУ устойчива. На рис.79б это выглядит так,
что пересечению ЛАЧХ горизонтальной оси соответствует точка ЛФЧХ, расположенная
выше линии
= –. И
наоборот для неустойчивой замкнутой САУ (рис.79а кривая 2) _c2
< –,
поэтому при
= _c2
ЛФЧХ проходит ниже линии
= –. Угол
₁
= _c1-(-)
является запасом устойчивости по фазе. Этот угол соответствует расстоянию от
линии
= – до
ЛФЧХ.

Исходя из сказанного, критерий устойчивости Наквиста по логарифмическим
ЧХ, в случаях, когда АФЧХ только один раз пересекает отрезок вещественной оси
[-;-1], можно
сформулировать так: для того, чтобы замкнутая САУ была устойчива необходимо
и достаточно, чтобы частота, при которой ЛФЧХ пересекает линию
= –, была
больше частоты среза.

Если АФЧХ разомкнутой САУ имеет сложный вид (рис.80), то ЛФЧХ может несколько
раз пересекать линию
= –. В
этом случае применение критерия Найквиста несколько усложняется. Однако во многих
случаях данной формулировки критерия Найквиста оказывается достаточно.

Вопросы

1. Какие
   САУ считаются структурно устойчивыми и структурно неустойчивыми?
2. В каком
   квадранте уходит в бесконечность АФЧХ разомкнутой САУ если порядок астатизма равен
   трем? Является ли такая САУ структурно устойчивой в замкнутом состоянии:
3. Как
   сделать устойчивой структурно неустойчивую САУ?
4. Что называется запасом
   устойчивости по модулю?
5. Что называется запасом устойчивости по фазе?
6. В
   чем особенность определения запасов устойчивости для клювообразных САУ?
7. Как
   влияет коэффициент усиления САУ на запасы устойчивости?
8. Чему соответствуют
   на АФЧХ пересечение ЛАЧХ оси w?
9. Чему соответствуют на АФЧХ пересечение
   ЛФЧХ значения j = -p?
10. Что называется частотой среза?
11. Сформулируйте
    критерий Найквиста для логарифмических характеристик.
12. В чем особенность
    логарифмических характеристик, если АФЧХ имеет клювообразный характер?

Далее…

Но сначала краткое содержание предыдущих серий:

1. Введение в теорию автоматического управления.2. Математическое описание систем автоматического управления 2.1 — 2.3, 2.3 — 2.8, 2.9 — 2.13. 

3. ЧАСТОТНЫЕ ХАРАКТЕРИСТИКИ ЗВЕНЬЕВ И СИСТЕМ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ РЕГУЛИРОВАНИЯ. 3.1. Амплитудно-фазовая частотная характеристика: годограф, АФЧХ, ЛАХ, ФЧХ. 3.2. Типовые звенья систем автоматического управления регулирования. Классификация типовых звеньев. Простейшие типовые звенья. 3.3. Апериодическое звено 1–го порядка инерционное звено. На примере входной камеры ядерного реактора. 3.4. Апериодическое звено 2-го порядка. 3.5. Колебательное звено. 3.6. Инерционно-дифференцирующее звено. 3.7. Форсирующее звено.  3.8. Инерционно-интегрирующее звено (интегрирующее звено с замедлением). 3.9. Изодромное звено (изодром). 3.10 Минимально-фазовые и не минимально-фазовые звенья. 3.11 Математическая модель кинетики нейтронов в «точечном» реакторе «нулевой» мощности. 

4. Структурные преобразования систем автоматического регулирования.

5. Передаточные функции и уравнения динамики замкнутых систем автоматического регулирования (САР).

6. Устойчивость систем автоматического регулирования. 6.1 Понятие об устойчивости САР. Теорема Ляпунова. 6.2 Необходимые условия устойчивости линейных и линеаризованных САР. 6.3 Алгебраический критерий устойчивости Гурвица. 6.4 Частотный критерий устойчивости Михайлова. 6.5 Критерий Найквиста.

6.6 Понятие об областях устойчивости

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

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

Предположим, что разомкнутая САР имеет передаточную функцию вида:

Если САР замкнута, то ее передаточная будет иметь вид:

Где– характеристически полином замкнутой САР.

Характеристическое уравнение такой системы имеет вид:

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

Предположим, что коэффициент изменяется от до . Будем давать фиксированные значения и определять значения всех “n” корней уравнения. Тогда при некоторых значениях все корни либо отрицательные, либо имеют отрицательную вещественную часть (расположены в левой полуплоскости). И наоборот: появляются корни с положительной вещественной частью (т.е. корни, расположенные в правой полуплоскости). Тогда на можно построить полуось изменения , определены области устойчивости, т.е. отрезки значений , при которых САР устойчива.

Если изменять 2 коэффициента в характеристическом уравения, например и (или коэффициенты и в передаточной функции САР), то мы получим области устойчивости на плоскости параметров:

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

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

Например, для системы 2-го порядка c характеристическим полиномом согласно критерию Гурвица необходимым и достаточным условием устойчивости является положительность коэффициентов и .

Вспоминая лекцию 6, “Если корни комплексные, то процесс колебательный, если корни реальные, то процесс апериодический (затухающий)”, можно сказать что:

ось ординат cоответсвует апериодической границе устойчивости, т.к. если, то корни характеристического уравления: ось абсцисс соответсвует колебательной границе устойчивости, т.к. если, то корни характеристического уравления: ;

Разбиение пространства коэффициентов (параметров) характеристического уравнения на области устойчивости и неустойчивости называется Д-разбиением. Данный подход, впервые предложил Ю.И. Неймарк (СССР) в 1948 г.

6.7 Д-разбиение плоскости по одному (комплексному параметру)

Рассмотрим замкнутую САР. Характеристическое уравнение имеет вид:

Где – корень характеристического полинома, при котором .

Предположим, что существует некоторый параметр САР, который входит линейно в один, несколько, или даже все все коэффициенты уравнения 6.7.1 (например, коэффициент усиления некоторого звена, или постоянная времени Т некоторого звена). Тогда, те коэффициенты, которые линейно зависят от можно представить как:

Выражение 6.7.2 может бы применено даже для одного слогаемого в полиноме. Подставляя выражение 6.7.2 в уравнение 6.7.1 получаем уравнение для определения корней характеристического полинома в виде:

где:
– часть полинома, которая не содержит множетель .
– часть полинома, которая содержит множитель .

Тогда можно записать выражение для :

В плоскости корнейсистема явялется устойчивой для корней в левой полуплоскость, а границей устойчивости – ось ординат (см. лекцию 6.1). Рассмотрим границу устойчивости – ось ординат, где у корнейнет действительной части.

Подставим в уравнение 6.7.3 вместо , значение , и рассмотрим изменение от до . При этом каждому значению можно найти такое значение множителя , что бы выражение 6.7.3 было равно нулю. В этом случае используя выражение 6.7.4 можно получить выражение для комплексного числа k, соотвествующего числу расположенному на граинице устойчивости:

где и – некторые функции от . Таким образом выражение 6.7.5 можно изобразить на комплексной полоскости.

Подставляя в выражение (6.7.5) значение   получаем точку 1 на плоскости ; подставляя значение получим точку 2 на плоскости ; подставляя   – третью и т.д. Соединив изображающие точки на плоскости k получаем линию Д-разбиения, т.е. отображение мнимой оси на плоскости  на плоскость k.

Линия Д-разбиения на плоскости  должна разделять области устойчивости и неустойчивости.

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

Для оценки областей устойивости штирхуется область слева от линии Д-разбиения.

Рассмотри точку А2 (см. рис. 6.7.1). Поскольку линия Д- разбиения формируют петлю, то непонятно точная принадлежность соотвесвующего корня к устойчивым или неустойчивым. Однако можно сказать, что переход из точки А2 в точку А3 (см. рисунок 6.7.1) не изменяет количества корней, лежащих в правой полуплоскости, поскольку нет перехода черех линию Д- различения;

Переход из точки А2 в точку А1 уменьшает количество неустойчивых корней уравнения (6.7.1) на 1, а переход из точки А2 в точку А4 – наоборот увеличивает на 1.

Для многих САР, при линия Д-разбиения проходит через начало координат (но не для всех бывают исключения);

Замкнутая область со штриховкой внутри (см. рис. 6.7.1), явялется кандидатом на область устойчивости. Чтобы удостоверится, что данная область – область устойчивости, необходимо какую-нибудь точку на вещественной оси внутри заштрихованной области подставить в исходное характеристическое уравнение  (6.7.1) и либо решить его, либо используя какой-либо критерий сделать вывод об устойчивости или неустойчивости САР.

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

Линия Д-разбиения обладает свойством зеркальной симметрии относительно оси абсцисс:

Поэтому линию D-разбиения можно строить только для одной половины, например для и зеракльно потом зеркальн отобразить.

Линию D-разбиения можено строить так же с учетом запаса устойчивости. В этом случае используется не ось ординат – граница устойчивости, а сдвинутая в устойчивую область линия. Вместо в уравнение 6.7.5 подставляется значение . Тогда при изменение от до мы получим область значении , при которых корни уравнения 6.7.1 будут расположены левее линии , см. рисунок 6.7.2

Пример 1.

Определить область устойчивости замкнутой системы при варьировании параметра  если САР имеет вид:

где: ;

Решение:

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

Характеристическое уравнение замкнутой САР имеет вид:

Выражение для коэффициента из характеристического уравнения

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

Определяем вещественную и минимую части выражения зависимости

Строим линию Д-разбиения по формулам (П.3) и (П.4), наносим штриховку с левой сторный от линии, если идти от до см. рис. 6.7.4

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

Подставляя в характеристический многочлен  

Воспользумеся критерием Гурвица :

лители:

Все определители матрицы Гурвица больше нуля, следовательно САР устойчива. Следовательно весь отрезок – область устойчивости.

Решение в среде структурного моделирования

Для постороения линии Д-разбиения можно воспользоваться расчетом в среде динамического моделирование SimInTech. На рисунке 6.7.6 Представлен скриншот проекта построения линии Д-разбиения, для примера 1.

С помощью главного скрипта проекта записываем выражения для и , (формулы П.3 и П.4). В скрипте переменные ReK и ImK. Используя переменную w в качестве переменной . Осуществляем присвоение времени моделирования переменной w (выражение w = time;). Это позволяет нам построить график от 0 до T_end – заданного значение вермени моделирования. Используя симметричность линии Д-разбиения относительно оси реальных значений, записываем выражения для части линии, соотвествующе отрицательным заначениям от – T_end, до 0. (см. рис. 6.7.5)

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

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

Соберем схему, для проверки устойчивости САР используя три значения для: из разных областей полученных при построеи Д-разбиения. Для этого введем новые переменные в галваном скрипте программы и присвоим им значения в секции инициализации:

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

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

Пример 2

Найти область устойчивости по коэффиценту для передаточной функции характеристичекий полином которой определяется выражением:

Характеристическое уравнение:

Выражения для коэффициента заменяем

Выражение для реальной и мнимой части комплексного параметра :

Строим линию Д-разбиения:

Область возможной устойчивость . Проверяем устойчивось путем моделирования при . Видим, что САР устойчива при .

Дальше будет метод расчета САР с циркулем штангелем и рабиновичем на милиметровой бумаги. Я предупредил!

6.8  Метод Д-разбиений на плоскости 2-х действительных параметров

Предположим, что структурная схема некоторой САР имеет вид:

Необходимо определить параметры регулятора (например К_рег и Т_рег) обеспечивающих устойчивость системы. В этом случае задача анализа сводится к поиску области устойчивости САР при варьировании 2-х параметров (в данном случае К и Т):

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

Поскольку рассматриваются линейные (или линеаризованные) САР, то характеристическое уравнение имеет вид:

Предположим, что один или несколько коэффициентов   линейно зависит от и :

Подставляя значения коэффициентов   в уравнение (6.8.1), получаем:

Где – полиномы по степеням

Где – постоянные коэффициенты.

В плоскости корнейсистема явялется устойчивой для корней в левой полуплоскость, а границей устойчивости – ось ординат (см. лекцию 6.1). Рассмотрим границу устойчивости – ось ординат, где у корнейнет действительной части. Подставляя в (6.8.3)  получаем:

g«Отображаем» точки мнимой оси на плоскости  на плоскость (K,T) необходимо решить уравнение:

Совершенно очевидно, что подстановка  в уравнение (6.8.4) обращается в тождество при некоторых значениях К и Т. (например, точке при   соответствует точка 1,  – точка 2,  – точка 3 и так далее.) Чтобы найти значение К и Т, соответствующие границе Д-разбиения, необходимо решить уравнение (6.8.4). Приравнивания по отдельности чисто вещественную и чисто мнимую части, получаем:

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

где: – главный определитель, – вспомогательные определители системы:

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

Свойства симметрии при изменении заключаются в том, что линия при измененеи от 0 до линия Д-разбиения совпадает с линией при измени от до 0.

Например, если  изменяется от  (точка С на рис) до нуля (точка А), то линия Д-разбиения – СВА, причем если построить линию Д-разбиения при  изменяющейся от нуля до  , то эта линия совпадает с линией АВС, т.е. путь от  до 0 и путь от 0 до  проходит по одной и той же линии.

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

Правила штриховки

Направление штриховки зависит от знака главного определителя системы .

1) Если при движении по кривой Д-разбиения от меньших значений   (от ) к большим (до ) главный определитель больше 0 (), то штриховка слева.

2) Если при движении по кривой Д-разбиения от меньших значений   (от ) к большим (до ) главный определитель меньше 0 () то штриховка справа.

Учитывая, что точки на плоскости (К,Т) при  и совпадают (что позволяет строить только одну ветвь Д-кривой, например от 0 до ), легко показать, что главный определитель  при  и имеет разные знаки. Это означает, что на плоскость (К,Т) с какой-то из сторон линии Д-разбиение будет нанесена двойная штриховка.

Важной особенностью Д-разбиения определяемой на основании решения системы (6.8.5) является наличие случаев, когда все определители ( ) обращаются одновременно в нуль. То есть при определенных значениях

Это означает, что вместо системы (6.8.5) – системы 2-х уравнений с 2-мя неизвестными получается одно уравнение с двумя неизвестными, т.к. уравнения становятся линейно-зависимыми (одно получается из другого линейной операцией, например умножением на какое-то число).

Такая же ситуация возникает если все определители равны бесконечности.

Наиболее частый случай – особая прямая проходит через точку плоскости (К,Т) при или при   (Прямая 2 на рисунке).

К разряду особых прямых необходимо отнести и прямые 1 и 3 на упомянутом рисунке (проходящие через точки Bи E, соответственно).

Правила штриховки в случае особых прямых учитывают, как ведет себя главный определитель  при «прохождении» соответствующей точки на Д-кривой, где он равен нулю:

-если главный определитель  меняет знак с + на – или наоборот, то особая прямая штрихуется одинарной штриховкой в соответствии с зонами соприкосновения, т. е. учитывается штриховка на кривой Д-разбиения вблизи этой точки. (Прямые 1 и 2 соответственно).

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

Определение возможных зон устойчивости

Далее выявляются замкнутые (или полузамкнутые области), образованные штриховкой вовнутрь. Такая область – возможная область устойчивости см. рис. 6.8.4

Далее необходимо взять какую-нибудь точку М из области (К=К_м, T=T_м) и подставив значение К и Т в характеристический полином системы   определить устойчива или нет САР (например по Гурвицу или непосредственным вычислением всех полюсов или корней). Если при К=К_м, T=T_м САР устойчива, то она будет устойчива и при любых значениях К и Т из этой области, если нет, то возможная область не стала областью устойчивости.

Дополнение

1) Выше предполагалось, что параметры К и Т линейно входят в коэффициенты характеристического уравнения. Если это не так, то система, подобная (6.8.5), принимает вид:

где  – нелинейная функция содержащая реальную часть уравнения 6.8.4;

 – нелинейная функция содержащая мнимую часть уравнения 6.8.4;

Решение этой системы (например по Ньютону-Рафсону) может дать неединственность решений. Тем не менее можно построить линия Д-разбиения. Правило штриховки – аналогичное вышеприведенному, т.е. если >0, то слева и наоборот. В качестве главного определителя  выступает якобиан:

2) Метод Д-разбиения в полной мере применим и для систем с постоянным запаздыванием, когда характеристическое уравнение является не полиномом, а квазиполиномом:

Использовали вы когда нибудть миллиметровку для графических расчетов САР?

Проголосовали 11 пользователей.

Воздержался 1 пользователь.