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

Содержание

Слайд 2

САУ устойчива если после кратковременного возмущения она возвращается в прежнее или занимает новое устойчивое положение.

САУ устойчива если после кратковременного возмущения она возвращается в прежнее или

занимает новое устойчивое положение.
Слайд 3

ПРИМЕРЫ Устойчива ∆х р

ПРИМЕРЫ Устойчива

∆х

р

Слайд 4

Неустойчива ∆х р ∆х

Неустойчива

∆х

р

∆х

Слайд 5

Устойчива “в малом”

Устойчива “в малом”

Слайд 6

Линейные САУ описываются линейными дифференциальными уравнениями (ДУ). Для решения ДУ следует найти корни его характеристического уравнения:

Линейные САУ описываются линейными дифференциальными уравнениями (ДУ). Для решения ДУ следует

найти корни его характеристического уравнения:
Слайд 7

a) Вещественные корни

a) Вещественные корни

Слайд 8

б) Комплексно-сопряженные корни

б) Комплексно-сопряженные корни

Слайд 9

Слайд 10

Итак, чтобы САУ была устойчивой необходимо и достаточно, чтобы все корни

Итак, чтобы САУ была устойчивой необходимо и достаточно, чтобы все корни

ее характеристического уравнения находились в левой полуплоскости комплексной плоскости.
Слайд 11

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

Критерии устойчивости САУ – это некоторые признаки позволяющие не решая характеристического

уравнения оценить устойчивость САУ.
Слайд 12

ВНИМАНИЕ Характеристическое уравнение замкнутой САУ – это знаменатель ее передаточной функции

ВНИМАНИЕ

Характеристическое уравнение замкнутой САУ – это знаменатель ее передаточной функции

Ф(s) приравненный “0”.
Характеристическое уравнение разомкнутой САУ - это знаменатель ее передаточной функции Wp(s) приравненный “0”.
Слайд 13

А. Алгебраические критерии устойчивости САУ Критерий Гурвица (1895г.). Пусть дано ХУ замкнутой САУ anpn+an-1pn-1+…+a0=0 (1)

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

Критерий Гурвица (1895г.).
Пусть дано ХУ замкнутой САУ


anpn+an-1pn-1+…+a0=0 (1)
Слайд 14

Составим определитель Гурвица из коэффициентов ХУ: ∆n= (2)

Составим определитель Гурвица из коэффициентов ХУ:
∆n=

(2)

Слайд 15

Как видно из (2): На главной диагонали определителя Гурвица располагаются сверху

Как видно из (2):

На главной диагонали определителя Гурвица располагаются сверху вниз

коэффициенты ХУ начиная со второго.
Выше элементов главной диагонали выписываются коэффициенты при младших степенях “р” по мере их убывания.
Ниже элементов главной диагонали выписываются коэффициенты при старших степенях “р” по мере их возрастания.
Остальные элементы определителя Гурвица равны “0”.
Слайд 16

Составим главные диагональные миноры ∆1= an-1 ∆2 = ∆3=

Составим главные диагональные миноры

∆1= an-1
∆2 =
∆3=

Слайд 17

1. Критерий Гурвица: Для устойчивости линейной САУ необходимо и достаточно, чтобы

1. Критерий Гурвица:

Для устойчивости линейной САУ необходимо и достаточно, чтобы при

аn>0 все главные диагональные миноры определителя Гурвица были бы положительны.
Слайд 18

Примеры 1. n=1 a1p+a0=0 Условия устойчивости a1>0 ∆1=a0>0 2. n=2 a2p2+a1p+a0=0

Примеры

1. n=1 a1p+a0=0
Условия устойчивости
a1>0 ∆1=a0>0
2. n=2 a2p2+a1p+a0=0
Условия устойчивости
a2>0 ∆1=a1>0
∆2= =

a1a0>0
Слайд 19

3. n=3 a3p3+a2p2+a1p=0 Условия устойчивости a3>0 ∆1=a2>0 ∆2 = =a2a1-a3a0>0 ∆3= =a0*∆2>0


3. n=3 a3p3+a2p2+a1p=0
Условия устойчивости
a3>0 ∆1=a2>0
∆2 = =a2a1-a3a0>0
∆3= =a0*∆2>0

Слайд 20

Недостаток критерия Гурвица С увеличением “n” раскрывать определители становится трудно.

Недостаток критерия Гурвица

С увеличением “n” раскрывать определители становится трудно.

Слайд 21

Пример для КСР Пусть дана структура замкнутой САУ

Пример для КСР

Пусть дана структура замкнутой САУ

Слайд 22

Необходимо с помощью критерия Гурвица определить устойчивость САУ. План исследований: 1.

Необходимо с помощью критерия Гурвица определить устойчивость САУ. План исследований:

1. Найти передаточную

функцию замкнутой САУ.
2. Определить ХУ замкнутой САУ и его коэффициенты.
3. Составить определитель Гурвица.
4. Определить все главные диагональные миноры и оценить устойчивость САУ по критерию Гурвица.
Слайд 23

2. Критерий Рауса Пусть дано ХУ замкнутой САУ “n”–го порядка: anpn

2. Критерий Рауса

Пусть дано ХУ замкнутой САУ “n”–го порядка:
anpn +

an-1pn-1 +… +a1p + a0 = 0
Слайд 24

Составим таблицу Рауса из коэффициентов ХУ:

Составим таблицу Рауса из коэффициентов ХУ:


Слайд 25

Критерий Рауса Для устойчивости линейной САУ необходимо и достаточно, чтобы коэффициенты

Критерий Рауса

Для устойчивости линейной САУ необходимо и достаточно, чтобы коэффициенты первого

столбца таблицы Рауса были положительны, т.е.:
С11> 0 c12 > 0 … c1,n+1 > 0
Слайд 26

Пример I для КСР Пусть ХУ замкнутой САУ: P6 + 6p5

Пример I для КСР

Пусть ХУ замкнутой САУ:
P6 + 6p5 + 21p4

+ 44p3 + 62p2 + 52 + 100 =0
Необходимо исследовать устойчивость этой системы используя критерий Рауса.
Слайд 27

План исследования Составим таблицу Рауса и заполним ее первые две строки.

План исследования

Составим таблицу Рауса и заполним ее первые две строки.
Вычислим последовательно

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

Итак, составим таблицу Рауса

Итак, составим таблицу Рауса

Слайд 29

Задание по КСР: Завершить заполнение таблицы Рауса и оценить устойчивость САУ.

Задание по КСР:

Завершить заполнение таблицы Рауса и оценить устойчивость САУ.

Слайд 30

Б. Частотные критерии устойчивости САУ Критерий Михайлова (1938) Дано ХУ замкнутой

Б. Частотные критерии устойчивости САУ

Критерий Михайлова (1938)
Дано ХУ замкнутой линейной САУ:
А(s)

= ansn + an-1sn-1 + … + a0 = 0 (1)
Слайд 31

Представим полином (1) в виде: A(s) = an (s – s1)

Представим полином (1) в виде:

A(s) = an (s – s1) (s

–s2) … (s - sn) (2)
Где si – корни ХУ
i = 1, 2 … n
Положим s = jω, тогда:
А(jω) = an (jω – s1)(jω – s2)… (jω - sn) (3)
Слайд 32

Каждая из скобок (3) представляет собой вектор, начало которого лежит в

Каждая из скобок (3) представляет собой вектор, начало которого лежит в

т. si, а конец находится на мнимой оси jω. При этом возможно два варианта:

При -∞ < jω < +∞
∆arg (jω – si) = +

При -∞ < jω < +∞
∆arg (jω – si) = -

1 – корень лежит в левой
полуплоскости

2-корень лежит в правой
полуплоскости

Слайд 33

Итак, если ХУ A(s) = 0 содержит L корней в правой

Итак, если ХУ A(s) = 0 содержит L корней в правой

полуплоскости, то
∆arg A(jω) = (n-2L) = (n-L) - L
Слайд 34

Для устойчивости линейной САУ необходимо, чтобы L = 0, т.е.: ∆arg A(jω) = n (4)

Для устойчивости линейной САУ необходимо, чтобы L = 0, т.е.:
∆arg

A(jω) = n (4)
Слайд 35

Критерий Михайлова является графической интерпретацией выражения (4). При этом рассматриваются лишь

Критерий Михайлова является графической интерпретацией выражения (4).
При этом рассматриваются лишь

положительные частоты, т.е.:
∆arg A(jω) = n * (5)
Слайд 36

Критерий Михайлова Для устойчивости линейной САУ необходимо и достаточно, чтобы годограф

Критерий Михайлова

Для устойчивости линейной САУ необходимо и достаточно, чтобы годограф Михайлова

(5) А(jω), начинаясь при ω = 0 на действительной оси с ростом “ω” от “0” до “∞” обходил последовательно “n” квадрантов против часовой стрелки (где n – порядок характеристического уравнения).
Слайд 37

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

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

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

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

Слайд 38

ПРИМЕР Определить предельный коэффициент Кпр при котором САУ теряет устойчивость, если ее структура имеет вид:

ПРИМЕР

Определить предельный коэффициент Кпр при котором САУ теряет устойчивость, если ее

структура имеет вид:
Слайд 39

Найдем передаточную функцию замкнутой САУ: ХY САУ – это знаменатель ее

Найдем передаточную функцию замкнутой САУ:
ХY САУ – это знаменатель ее передаточной

функции приравненный к 0 т.е.:
Слайд 40

3. Годограф Михайлова (при s = jω): А(jω) =D(jω) + К

3. Годограф Михайлова (при s = jω):
А(jω) =D(jω) + К

0 < ω < ∞
4. Построим, вначале, D(jω):
D(jω) = (Т1jω+1)(T2jω+1)(T3jω+1)=Re(ω) + jIm(ω)
Re(ω) = 1 – (T1T2 + T1T3 + T2T3)ω2
Im(ω) = (T1 + T2 + T3)ω –T1T2T3ω3
Слайд 41

Кпр определим из уравнений Кпр определим из уравнений

Кпр определим из уравнений

Кпр определим из уравнений

Слайд 42

НЕДОСТАТОК критерия Михайлова Годограф Михайлова не имеет физической сущности (его нельзя

НЕДОСТАТОК критерия Михайлова

Годограф Михайлова не имеет физической сущности (его нельзя получить

экспериментально). Между тем при исследовании сложных систем хотелось бы опираться на характеристики получаемые не только аналитически, но и экспериментально.
Слайд 43

2. Критерий Найквиста (1932) Основан на использовании wp(s), которую можно получить

2. Критерий Найквиста (1932)

Основан на использовании wp(s), которую можно получить экспериментально.
Пусть:

- ПФ разомкнутой САУ
Тогда: - ПФ замкнутой САУ
Слайд 44

Образуем функцию: - XY замкнутой САУ - XY разомкнутой САУ

Образуем функцию:
- XY замкнутой САУ
- XY разомкнутой САУ

Слайд 45

РАССМОТРИМ 1-й случай – разомкнутая САУ устойчива. Тогда, согласно критерию Михайлова: ∆arg N(jω) = n* 0

РАССМОТРИМ

1-й случай – разомкнутая САУ устойчива.
Тогда, согласно критерию Михайлова:
∆arg N(jω) =

n*
0 < ω < ∞
Слайд 46

Чтобы система и в замкнутом состоянии была устойчива необходимо чтобы: ∆arg

Чтобы система и в замкнутом состоянии была устойчива необходимо чтобы:
∆arg
0 <

ω < ∞
Это значит что: ∆arg F(jω)= 0
0 < ω < ∞
Слайд 47

Изобразим F(jω) на комплексной плоскости

Изобразим F(jω) на комплексной плоскости

Слайд 48

Сдвинем теперь F(jω) влево на “1” и получим, таким образом, wp(jω)

Сдвинем теперь F(jω) влево на “1” и получим, таким образом, wp(jω)

Слайд 49

Критерий устойчивости Найквиста: Если разомкнутая САУ устойчива, то для ее устойчивости

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

Если разомкнутая САУ устойчива, то для ее устойчивости в

замкнутом состоянии необходимо и достаточно, чтобы АФЧХ wp(jω) при 0<ω<∞ не охватывала точку с координатами (-1; j0).
Слайд 50

ПРИМЕР Найти, используя критерий Найквиста, предельный коэффициент Кпр при котором САУ потеряет устойчивость если:

ПРИМЕР

Найти, используя критерий Найквиста, предельный коэффициент Кпр при котором САУ потеряет

устойчивость если:
Слайд 51

По критерию Найквиста САУ находится на границе устойчивости если: Полагая Im(ω ) = 0 найдем:

По критерию Найквиста САУ находится на границе устойчивости если:
Полагая Im(ω )

= 0 найдем:
Слайд 52

Подставив в Re(ω ) найдем: Т.е. результат такой же, как и при использовании критерия Михайлова.

Подставив в Re(ω ) найдем:
Т.е. результат такой же, как и при

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

Запасы устойчивости САУ – это количественные оценки степени устойчивости систем. Удобнее

Запасы устойчивости САУ – это количественные оценки степени устойчивости систем. Удобнее

всего такие оценки (запасы) получить используя критерий Найквиста

Рассмотрим АФЧХ устойчивой САУ

Wp(jω) при 0<ω<∞ не охватывает т.(-1, j0).
Следовательно САУ устойчива.

Слайд 54

САУ может потерять устойчивость по двум причинам: а) увеличения К без

САУ может потерять устойчивость по двум причинам:
а) увеличения К без изменения

фаз - все вектора wp(jω) увеличиваются и когда-нибудь САУ станет неустойчивой. Очевидно, что увеличивать К можно
в раз т.ч.
∆А= - запас устойчивости САУ по амплитуде.
Слайд 55

б) увеличения φ(ω) без изменения К – все вектора характеристики wp(jω)

б) увеличения φ(ω) без изменения К – все вектора характеристики wp(jω)

поворачиваются по часовой стрелке на некоторые углы ∆φ. На рисунке видно на какой угол ∆φ можно повернуть wp(jω) прежде чем САУ потеряет устойчивость.
Слайд 56

Проводя окружность радиусом “1” можно найти ту точку ω , которая

Проводя окружность радиусом “1” можно найти ту точку ω , которая

попадет в точку (-1; j0) если на частоте ω φ(ω ) увеличится на угол ∆φ.
Следовательно ∆φ – запас устойчивости САУ на фазе.
Слайд 57

Итак, существуют две количественные оценки степени устойчивости САУ Запас “по амплитуде”

Итак, существуют две количественные оценки степени устойчивости САУ

Запас “по амплитуде” -

∆А=
Запас “по фазе” - ∆φ
Недостаток частотных критериев устойчивости – сложно строить кривые А(jω) и wp(jω)
Слайд 58

Анализ устойчивости САУ по логарифмическим характеристикам АФЧХ можно построить в логарифмическом

Анализ устойчивости САУ по логарифмическим характеристикам

АФЧХ можно построить в логарифмическом масштабе

в виде двух характеристик:
L(ω) – логарифмической амплитудной частотной характеристики
φ(ω) – фазовой частотной характеристики.
Слайд 59

Тогда анализ устойчивости заметно упрощается. Особенно просто строятся т.н. асимптотические L(ω) – в виде кусочно-прямолинейных характеристик.

Тогда анализ устойчивости заметно упрощается. Особенно просто строятся т.н. асимптотические L(ω)

– в виде кусочно-прямолинейных характеристик.
Слайд 60

РАССМОТРИМ вначале два элемента таких кусочно-ломанных L(ω) Пусть: Тогда:

РАССМОТРИМ вначале два элемента таких кусочно-ломанных L(ω)

Пусть:
Тогда:

Слайд 61

Пусть: Тогда:

Пусть:
Тогда:

Слайд 62

Приближенно: при ω при ω > Итак L(ω) состоит из двух

Приближенно:
при ω <
при ω >
Итак L(ω) состоит

из двух прямых (асимптот):
1 – совпадающей с осью ω при ω <
2 – имеющей наклон –20 дб/дек при ω >
Слайд 63

Частота ω = = ωс называется сопрягающей. На сопрягающей частота ωс

Частота ω = = ωс называется сопрягающей.
На сопрягающей частота ωс =
φ(ω)

= - arctg1 = -450
При ω→∞ φ(ω) = -arctg∞→ -900
ω→0 φ(ω) = -arctg0 →00
Слайд 64

Итак, чтобы построить L(ω) и φ(ω) для этого элемента – w(jω)

Итак, чтобы построить L(ω) и φ(ω) для этого элемента – w(jω)

= нужно:

Найти сопрягающую частоту:
Вдоль оси ω построить участок 1 для
ω < ωс
Построить участок 2 с наклоном -20дб/дек для
ω < ωс

Слайд 65

4. По формуле φ(ω)= -arctgωT задаваясь разными частотами 0

4. По формуле φ(ω)= -arctgωT
задаваясь разными частотами 0<ω<∞ построить фазовую

частотную характеристику L(ω)
Слайд 66

Пусть теперь: Тогда: Приближенно: при ω при ω >


Пусть теперь:
Тогда:
Приближенно: при ω <
при ω >

Слайд 67

Т.О. и здесь L(ω) состоит из двух участков: 1 – вдоль

Т.О. и здесь L(ω) состоит из двух участков:
1 – вдоль оси

ω до ω ≤ ωс =
2 – с наклоном +20дб/дек при ω >
Слайд 68

Слайд 69

Построение асимптотических L(ω) и φ(ω) для сложных САУ. Пусть например: Заменив S→jω получим амплитудно-фазовые частотные характеристики:

Построение асимптотических L(ω) и φ(ω) для сложных САУ.

Пусть например:
Заменив S→jω

получим амплитудно-фазовые частотные характеристики:
Слайд 70

Представим последнюю характеристику в виде произведения характеристик элементарных звеньев:

Представим последнюю характеристику в виде произведения характеристик элементарных звеньев:

Слайд 71

Определим сопрягающие частоты:

Определим сопрягающие частоты:

Слайд 72

ω1/с Построим участок 1: W(iω) = 100/jω A(ω)= 100/ω L(ω)= 20

ω1/с

Построим участок 1:
W(iω) = 100/jω
A(ω)= 100/ω
L(ω)= 20 lg 100 –

20 lg ω
При ω=1
L(ω)= 20 lg 100= 40дб/дск
Построив участок 1 до ω= ωс1, строим участок II. Изменив наклон на -20дб/дск (т.к. скобка (jω+1) - в знаменателе!!!).
Участок II продляем до ω= ωс2 с наклоном -40дб/дск.
На ω ≥ ωс2 снова изменяем наклон, но уже на +20дб/дск, т.к. скобка (0,1S+1) стоит в числителе Wp(S). Т.о. на участке III наклон снова становится -20дб/дск до частоты ω= ωс3. На частоте ωс3 наклон участка IV снова равен -40дб/дск, т.к.скобка (0,01S+1) стоит в знаменателе Wp(S)
Слайд 73

Фазовая характеристика φ(ω) САУ складывается из фазовых характеристик отдельных звеньев Wp(S): φ(ω) = –90°−(arctg1ω)+(arctg0,1ω) – (arctg0,01ω)

Фазовая характеристика φ(ω) САУ

складывается из фазовых характеристик отдельных звеньев Wp(S):
φ(ω)

= –90°−(arctg1ω)+(arctg0,1ω) – (arctg0,01ω)
Слайд 74

Для ее построения удобно построить таблицу Фазовая характеристика строится по точкам

Для ее построения удобно построить таблицу

Фазовая характеристика строится по точкам

под амплитудной,
причем масштаб по оси “ω” тот же.
Слайд 75

Об устойчивости САУ судят по расположению точек пересечения L(ω) оси частот

Об устойчивости САУ судят по расположению точек пересечения L(ω) оси частот

ωср –частота среза и φ(ω) оси -180°.