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

     СОДЕРЖАНИЕ

 

      1 АНАЛИЗ  И СИНТЕЗ САУ

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

     1.1.1 Преобразование  структурной схемы САУ к эквивалентной.

     Исходная  схема имеет вид:

       

       
 
 

     Буквенным обозначениям параметров соответствуют  следующие числовые значения:

     K2 = 5;

     K3 = 2;

     K4 = 0.1;

     T1 = 0.02;

     T2 = 0.04;

     T3 = 0.12.

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

       

     

     Раскрывая скобки, получаем следующий вид этой передаточной функции:

     

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

     

     

 

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

 

       

     Производим замену переменной s на j·w: 

       

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

       

      ,

     откуда

     

     .

     

 

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

 

       

     Производим замену переменной s на j·w. 

       

     Для удобства работы по разбиению передаточной функции на вещественную и мнимую составляющие, заменим слагаемые знаменателя буквами и произведём с ними необходимые действия: 

      ; 

      ,

     откуда 

     

     

 

     

     1.2 Построение частотных  характеристик исходной САУ

 

     Построение  частотных характеристик выполняется в среде MathCAD.

     1.2.1 Частотные характеристики  разомкнутой исходной системы.

 

     Амплитудно-фазо-частотная  характеристика. 

     ω (0 ; 1000)

     Рисунок 1.2.1.1 — АФЧХ разомкнутой системы 

Таблица 1.2.1.1 — Данные для построения АФЧХ разомкнутой системы

ω 0 2 6 10 20 40 60 80 100
Up(ω) 1 0.908 0,437 0.068 -0.171 -0.084 -0.032 -0.013 -0.006
Vp(ω) 0 -0.337 -0.650 -0.579 -0.220 -0.008 0.011 0.009 0.006

 

      Амплитудно-частотная характеристика. 

       

     ω (0 ; 100)

     Рисунок 1.2.1.2 — АЧХ разомкнутой системы 

Таблица 1.2.1.2 — Данные для построения АЧХ разомкнутой системы

ω 0 10 20 30 40 50 60 80 100
Up(ω) 1 0.583 0.279 0.147 0.084 0.052 0.034 0.016 0.009

 

      Фазо-частотная характеристика. 

       

     ω (0 ; 100)

 

     Рисунок 1.2.1.3 — ФЧХ разомкнутой системы 

Таблица 1.2.1.3 — Данные для построения ФЧХ разомкнутой системы

ω 0 10 20 30 40 50 60 80 100
φp(ω),рад. 0 -1.454 -2.231 -2.716 -3.052 -3.298 -3.485 -3.747 -3.921
 

 

      Логарифмическая амплитудно-частотная  характеристика. 

       

     ω (0,1 ; 1000) 

 

     Рисунок 1.2.1.4 — ЛАЧХ разомкнутой системы 

Таблица 1.2.1.4 — Данные для построения ЛАЧХ разомкнутой системы

ω 0.1 1 4 10 20 60 100 400 1000
Lр(ω), дБ -0.00 -0.07 -1.04 -4.69 -11.09 -29.40 -40.91 -151.71 -199.32
 

 

      Логарифмическая фазо-частотная характеристика. 

     ω (0,1 ; 1000)

 

     Рисунок 1.2.1.5 — ЛФЧХ разомкнутой системы 

Таблица 1.2.1.5 — Данные для построения ЛФЧХ разомкнутой системы

ω 0.1 1 4 10 20 40 60 100 1000
Lp(ω), рад. -0.018 -0.179 -0.686 -1.454 -2.231 -3.05 -3.485 -3.92 -4.63

 

      1.2.2 Частотные характеристики замкнутой исходной системы.

     Амплитудно-фазо-частотная  характеристика.

     ω (0 ; 1000)

     Рисунок 1.2.2.1 — АФЧХ замкнутой системы 

Таблица 1.2.2.1 — Данные для построения АФЧХ замкнутой системы

ω 0 2 6 10 20 40 60 100 200
Uз(ω) 5 4.917 4.224 2.763 -1.268 -0.917 -0.328 -0.064 -0.005
Vз(ω) 0 -0.898 -2.614 -3.923 -2.994 -0.090 -0.122 0.064 0.011

 

      Амплитудно-частотная характеристика. 

       

     ω (0 ; 100)

     Рисунок 1.2.2.2 — АЧХ замкнутой системы 

Таблица 1.2.2.2 — Данные для построения АЧХ замкнутой системы

ω 0 10 20 30 40 50 60 80 100
Uз(ω) 5 4.798 3.251 1.692 0.922 0.548 0.350 0.166 0.091
 

 

      Фазо-частотная характеристика. 

       

     ω (0 ; 100)

 

     Рисунок 1.2.2.3 — ФЧХ замкнутой системы 

Таблица 1.2.2.3 — Данные для построения ФЧХ замкнутой системы

ω 0 10 20 30 40 50 60 80 100
φз(ω), рад. 0 -0.957 -1.972 -2.646 -3.044 -3.307 -3.497 -3.757 -3.927
 

 

      Логарифмическая амплитудно-частотная  характеристика. 

       

     ω (0.1 ; 1000)

 

     Рисунок 1.2.2.4 — ЛАЧХ замкнутой системы 

Таблица 1.2.2.4 — Данные для построения ЛАЧХ замкнутой системы

ω 0.1 1 4 10 20 40 100 400 1000
Lз(ω),дБ 13.98 13.98 13.96 13.62 10.24 -0.71 -20.85 -55.85 -79.66
 

 

      Логарифмическая фазо-частотная характеристика. 

       

     ω (0,1 ; 1000) 

 

     Рисунок 1.2.2.5 — ЛФЧХ замкнутой системы

Таблица 1.2.2.5 — Данные для построения ЛФЧХ замкнутой системы

ω 0.1 1 4 10 20 40 60 100 1000
φз(ω), рад. -0.009 -0.09 -0.364 -0.957 -1.97 -3.044 -3.497 -3.927 -4.63

 

      1.3 Анализ устойчивости САУ.

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

     Для построения годографа Михайлова, необходимо представить характеристическое уравнение передаточной функции замкнутой системы в комплексной форме, заменив переменную s на j·ω, и разбив получившееся представление на вещественную и мнимую части. Эта операция производилась на этапе разбиения передаточной функции замкнутой системы на вещественную и мнимую, поэтому, воспользуемся её результатами:

      — вещественная часть;

      — мнимая часть.

     Теперь, строим годограф Михайлова на комплексной плоскости:

     ω (0 ; 100)

Рисунок 1.3.1.1 — годограф Михайлова

Таблица 1.3.1.1 — Данные для построения годографа Михайлова

ω 0 2 10 20 40 60 100 400 1000
Cз(ω) 2 1.968 1.2 -1.2 -10.8 -26.8 -78 -1278 -7998
Dз(ω) 0 0.359 1.704 2.832 1.056 -9.936 -78 -6072 -95820
 

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

 

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

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

      . 

     Коэффициенты  характеристического уравнения  для определителя Гурвица нумеруем соответственно показателям степени переменной при них:

     a0=2; a1=0,18; a2=0,008; a3=0,000096;

     Определитель Гурвица:

     

     Подставляя  полученные значения, вычисляем его:

     

     Главный определитель Гурвица положителен. Аналогично исследуем все оставшиеся миноры.

        

        

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

 

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

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

      . 

     Коэффициенты  характеристического уравнения для таблицы Рауса нумеруем соответственно показателю степени переменной при них: 

     a0=2; a1=0,18; a2=0,008; a3=0,000096; 

     Таблица 1.3.1 — Таблица Рауса.

       

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

 

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

     Здесь используется АФЧХ разомкнутой системы:

   

     Рисунок 1.3.4.1 — годограф Найквиста

     При стремлении частоты в бесконечность, годограф приходит в начало координат, закручиваясь по часовой стрелке, и не охватывает точку с координатами (–1 ; j0), что свидетельствует об устойчивости как разомкнутой, так и замкнутой системы. 

     Все критерии оценки устойчивости показали, что система устойчива и в  замкнутом, и в разомкнутом состоянии.

 

      1.3.5 Построение области устойчивости САУ.

     Характеристическое  уравнение замкнутой системы  с общим коэффициентом усиления, принятым переменным (k), имеет вид:

     

     Выполним  преобразование :

     

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

     В данном случае, характеристическое уравнение  решается относительно коэффициента усиления k: 

       

     Действительная  часть:

     

     Мнимая  часть:

     

     На  графике D-разбиения наносится штриховка в сторону устойчивой области.

     Рисунок 1.3.5.1 — D-разбиение 

Таблица 1.3.5.1 — Данные для построения D-разбиения

ω 0 2 10 15 20 25 30 40 45 50 60
Ud(ω) -10 -8 -2 8 22 40 62 118 152 190 278
Vd(ω) 0 -8.88 -17.04 -23.76 -28.3 -30 -28.1 -10.56 6.48 30 99.4
 

     Как видно, вся плоскость по параметру  K разбивается на три зоны, разделяемые точками на оси Ud(ω) с абсциссами:

     

     

     Первая  область — (–∞ ; –10);

     вторая  область — (–10 ; 190,666667);

     третья  область — (190,666667 ; +∞).

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

 

        1.4 Построение переходного процесса системы методом трапеций

     Выполняем построение вещественной частотной характеристики замкнутой системы (рис. 1.4.1).

     

     Рисунок 1.4.1 — ВЧХ замкнутой САУ.

Таблица 1.4.1 — Данные для построения АФЧХ замкнутой системы

ω 0 2 10 21 25 30 50 70 100 120
Uз(ω) 5 4.92 2.76 -1.43 -1.67 -1.49 -0,54 -0,21 -0.64 -0.034

 

      Разбиваем ВЧХ на три трапеции.

     Рисунок 1.4.2 — Разбивка ВЧХ на трапеции.

     Синим цветом выделен контур первой трапеция, зелёным – второй, розовым – третьей.

     Определяем параметры трапеций: высоту и частоты начала и окончания наклонной стороны ( и соответственно).

     Для наглядности, совместим трапеции основаниями  с осью частот.

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