Анализ и синтез систем автоматического управления и исследование нелинейной системы СОДЕРЖАНИЕ 1 АНАЛИЗ И СИНТЕЗ САУ.. 4 1.1 Определение передаточной функции разомкнутой и замкнутой САУ. 4 1.1.1 Преобразование структурной схемы САУ к эквивалентной. 4 1.1.2 Разбиение на вещественную и мнимую составляющие передаточной функции разомкнутой системы. 6 1.1.3 Разбиение на мнимую и вещественную составляющие передаточной функции замкнутой системы. 7 1.2 Построение частотных характеристик исходной САУ.. 8 1.2.1 Частотные характеристики разомкнутой исходной системы. 8 1.2.2 Частотные характеристики замкнутой исходной системы. 13 1.3 Анализ устойчивости САУ. 18 1.3.1 Критерий Михайлова. 18 1.3.2 Критерий Гурвица. 20 1.3.3 Критерий Рауса. 21 1.3.4 Критерий Найквиста. 22 1.3.5 Построение области устойчивости САУ. 23 1.4 Построение переходного процесса системы методом трапеций. 26 1.5 Определение параметров и построение желаемой ЛАЧХ.. 31 1.5.1 Параметры для построения исходной ЛАЧХ. 31 1.5.2 Параметры для построения желаемой ЛАЧХ. 32 1.6 Определение передаточной функции корректирующего устройства. 35 1.7 Построение переходного процесса с использованием ПЭВМ.. 36 2 Исследование нелинейной системы. 39 2.1 Построение фазового портрета нелинейной САУ.. 40 2.2 Оценка устойчивости нелинейной САУ по критерию В.М. Попова. 42 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 — Разбивка ВЧХ на трапеции. Синим цветом выделен контур первой трапеция, зелёным – второй, розовым – третьей. Определяем параметры трапеций: высоту и частоты начала и окончания наклонной стороны ( и соответственно). Для наглядности, совместим трапеции основаниями с осью частот. Рисунок 1.4.3 — Трапеции, совмещённые по оси частот. На основании полученных результатов, строим таблицу. Таблица 1.4.2 — Параметры трапеций 1 трапеция 2 трапеция 3 трапеция Wd = 1.25 Wп = 25 Wd = 25 Wп = 50 Wd = 50 Wп = 100 Х = 0,05 R0 = 6.67 Х = 0,5 R0 = -1.13 Х = 0,50 R0 = -0,54 В таблице h -функций находим соответствующую каждому функцию . Искомую составляющую получаем из этой функции путём умножения ординат на величину . Время получаем как частное от деления величины на . ;. Таблица 1.4.3 — Значения для построения переходного процесса. t табл, h табл, t h1(t) t табл, h табл, t h2(t) t табл, h табл, t h2(t) 0 0,0000 0,0000 0,0000 0 0,0000 0,0000 0,0000 0 0,0000 0,0000 0,0000 0,5 0,1760 0,0200 1,1733 0,5 0,2400 0,0100 -0,2704 0,5 0,2400 0,0050 -0,1296 1 0,3400 0,0400 2,2667 1 0,4610 0,0200 -0,5194 1 0,4610 0,0100 -0,2489 1,5 0,4940 0,0600 3,2933 1,5 0,6650 0,0300 -0,7492 1,5 0,6650 0,0150 -0,3591 2 0,6280 0,0800 4,1867 2 0,8330 0,0400 -0,9385 2 0,8330 0,0200 -0,4498 2,5 0,7390 0,1000 4,9267 2,5 0,9670 0,0500 -1,0895 2,5 0,9670 0,0250 -0,5222 3 0,8280 0,1200 5,5200 3 1,0610 0,0600 -1,1954 3 1,0610 0,0300 -0,5729 3,5 0,8920 0,1400 5,9467 3,5 1,1150 0,0700 -1,2562 3,5 1,1150 0,0350 -0,6021 4 0,9380 0,1600 6,2533 4 1,1420 0,0800 -1,2867 4 1,1420 0,0400 -0,6167 4,5 0,9600 0,1800 6,4000 4,5 1,1380 0,0900 -1,2821 4,5 1,1380 0,0450 -0,6145 5 0,9770 0,2000 6,5133 5 1,1170 0,1000 -1,2585 5 1,1170 0,0500 -0,6032 5,5 0,9860 0,2200 6,5733 5,5 1,0920 0,1100 -1,2303 5,5 1,0920 0,0550 -0,5897 6 0,9820 0,2400 6,5467 6 1,0510 0,1200 -1,1841 6 1,0510 0,0600 -0,5675 6,5 0,9800 0,2600 6,5333 6,5 1,0180 0,1300 -1,1469 6,5 1,0180 0,0650 -0,5497 7 0,9790 0,2800 6,5267 7 0,9930 0,1400 -1,1188 7 0,9930 0,0700 -0,5362 7,5 0,9800 0,3000 6,5333 7,5 0,9740 0,1500 -1,0974 7,5 0,9740 0,0750 -0,5260 8 0,9850 0,3200 6,5667 8 0,9660 0,1600 -1,0884 8 0,9660 0,0800 -0,5216 8,5 0,9890 0,3400 6,5933 8,5 0,9660 0,1700 -1,0884 8,5 0,9660 0,0850 -0,5216 9 0,9970 0,3600 6,6467 9 0,9700 0,1800 -1,0929 9 0,9700 0,0900 -0,5238 9,5 1,0040 0,3800 6,6933 9,5 0,9750 0,1900 -1,0985 9,5 0,9750 0,0950 -0,5265 10 1,0090 0,4000 6,7267 10 0,9820 0,2000 -1,1064 10 0,9820 0,1000 -0,5303 10,5 1,0130 0,4200 6,7533 10,5 0,9870 0,2100 -1,1120 10,5 0,9870 0,1050 -0,5330 11 1,0150 0,4400 6,7667 11 0,9970 0,2200 -1,1233 11 0,9970 0,1100 -0,5384 11,5 1,0160 0,4600 6,7733 11,5 0,9970 0,2300 -1,1233 11,5 0,9970 0,1150 -0,5384 12 1,0150 0,4800 6,7667 12 0,9970 0,2400 -1,1233 12 0,9970 0,1200 -0,5384 12,5 1,0130 0,5000 6,7533 12,5 0,9970 0,2500 -1,1233 12,5 0,9970 0,1250 -0,5384 13 1,0120 0,5200 6,7467 13 0,9970 0,2600 -1,1233 13 0,9970 0,1300 -0,5384 13,5 1,0110 0,5400 6,7400 13,5 0,9980 0,2700 -1,1244 13,5 0,9980 0,1350 -0,5389 14 1,0110 0,5600 6,7400 14 1,0000 0,2800 -1,1267 14 1,0000 0,1400 -0,5400 14,5 1,0120 0,5800 6,7467 14,5 1,0020 0,2900 -1,1289 14,5 1,0020 0,1450 -0,5411 15 1,0120 0,6000 6,7467 15 1,0050 0,3000 -1,1323 15 1,0050 0,1500 -0,5427 15,5 1,0140 0,6200 6,7600 15,5 1,0080 0,3100 -1,1357 15,5 1,0080 0,1550 -0,5443 16 1,0150 0,6400 6,7667 16 1,0110 0,3200 -1,1391 16 1,0110 0,1600 -0,5459 16,5 1,0160 0,6600 6,7733 16,5 1,0110 0,3300 -1,1391 16,5 1,0110 0,1650 -0,5459 17 1,0160 0,6800 6,7733 17 1,0120 0,3400 -1,1402 17 1,0120 0,1700 -0,5465 17,5 1,0150 0,7000 6,7667 17,5 1,0090 0,3500 -1,1368 17,5 1,0090 0,1750 -0,5449 18 1,0150 0,7200 6,7667 18 1,0080 0,3600 -1,1357 18 1,0080 0,1800 -0,5443 18,5 1,0150 0,7400 6,7667 18,5 1,0060 0,3700 -1,1334 18,5 1,0060 0,1850 -0,5432 19 1,0150 0,7600 6,7667 19 1,0010 0,3800 -1,1278 19 1,0010 0,1900 -0,5405 19,5 1,0140 0,7800 6,7600 19,5 0,9980 0,3900 -1,1244 19,5 0,9980 0,1950 -0,5389 20 1,0130 0,8000 6,7533 20 0,9960 0,4000 -1,1222 20 0,9960 0,2000 -0,5378 20,5 1,0120 0,8200 6,7467 20,5 0,9950 0,4100 -1,1210 20,5 0,9950 0,2050 -0,5373 21 1,0110 0,8400 6,7400 21 0,9950 0,4200 -1,1210 21 0,9950 0,2100 -0,5373 21,5 1,0110 0,8600 6,7400 21,5 0,9960 0,4300 -1,1222 21,5 0,9960 0,2150 -0,5378 22 1,0110 0,8800 6,7400 22 0,9960 0,4400 -1,1222 22 0,9960 0,2200 -0,5378 22,5 1,0110 0,9000 6,7400 22,5 0,9970 0,4500 -1,1233 22,5 0,9970 0,2250 -0,5384 23 1,0110 0,9200 6,7400 23 0,9980 0,4600 -1,1244 23 0,9980 0,2300 -0,5389 23,5 1,0100 0,9400 6,7333 23,5 0,9990 0,4700 -1,1255 23,5 0,9990 0,2350 -0,5395 24 1,0100 0,9600 6,7333 24 1,0000 0,4800 -1,1267 24 1,0000 0,2400 -0,5400 24,5 1,0090 0,9800 6,7267 24,5 1,0000 0,4900 -1,1267 24,5 1,0000 0,2450 -0,5400 25 1,0080 1,0000 6,7200 25 1,0000 0,5000 -1,1267 25 1,0000 0,2500 -0,5400 25,5 1,0080 1,0200 6,7200 25,5 1,0000 0,5100 -1,1267 25,5 1,0000 0,2550 -0,5400 26 1,0070 1,0400 6,7133 26 1,0000 0,5200 -1,1267 26 1,0000 0,2600 -0,5400 При помощи функции линейной интерполяции linterp, встроенной в математический пакет MathCAD, находим функции переходных характеристик, соответствующих каждой из трапеций, и производим их сложение. Рисунок 1.4.4 — Графики переходных процессов замкнутой САУ: h1(t) — переходная характеристика первой трапеции; h2(t) — переходная характеристика второй трапеции; h3(t) — переходная характеристика третьей трапеции; H(t) — суммарная переходная характеристика. 1.5 Определение параметров и построение желаемой ЛАЧХ 1.5.1 Параметры для построения исходной ЛАЧХ. T1 = 0.02с T2 = 0.04с T3 = 0.12с ω1 = 50 с-1 ω2 = 25 с-1 ω3 = 8.333 с-1 20·log K = 20·log 1 = 0 дБ 1.5.2 Параметры для построения желаемой ЛАЧХ. В задании выданы следующие показатели качества регулирования САУ, которые необходимо соблюсти: σmax = 25% . Условно, весь частотный диапазон, в котором будет построена желаемая ЛАЧХ, разбивается на три поддиапазона: поддиапазон низких, средних и высоких частот (рисунок 1.5.1). ЛАЧХ поддиапазона средних частот строится с наклоном –20 дБ/дек, проходя через частоту среза до пересечения с линиями запас по амплитуде. Низко- и высокочастотная асимптота ЛАЧХ проводится с наклоном, кратным –20 дБ/дек, до соединения с ЛАЧХ исходной системы в этом диапазоне, или проводится параллельно ей (с одинаковым наклоном). Среднечастотная асимптота проводится через точку с наклоном –20 дБ/дек. Определим границы этой асимптоты. Для этого, по номограммам (рисунок 1.5.2) найдем запасы устойчивости по амплитуде. Рисунок 1.5.2 Номограммы для определения запасов устойчивости и частоты среза При σmax = 30%: Lmin = 14.5 дБ, ωс = (0.6 ÷ 0.9)·ωп ; ωс = (0.6 ÷ 0.9)·73,3 = (43,98 ÷ 65,97) с–1 . Низко- и высокочастотную асимптоты проводим параллельно исходной ЛАЧХ Построение исходной асимптотической ЛАЧХ производится исходя из значений частот нулей и полюсов разомкнутой передаточной функции. Определение ЛАЧХ корректирующего устройства производится по разности ЛАЧХ исходной и желаемой системы. 1.6 Определение передаточной функции корректирующего устройства. Производится по параметрам её асимптотической ЛАЧХ. Параметры корректирующего устройства: T1 = 0.04c; T2 = 0.005c; K = 5.3. Таким образом, передаточная функция корректирующего устройства принимает следующий вид: Передаточной функции WКУ (s) соответствует следующая схема корректирующего устройства: Рисунок 1.6.1 Схема электрическая принципиальная корректирующего устройства 1.7 Построение переходного процесса с использованием ПЭВМ При построении переходного процесса и оценки качества регулирования итоговой САУ используется математический пакет MathCAD. Переходный процесс должен удовлетворять следующим показателям качества: ≤0,15с, ≤30%. Рисунок 1.7.1 Переходный процесс скорректированной САУ Анализируя переходной процесс системы управления (рисунок 1.7.1), можем сказать, что время регулирования и перерегулирование, не выходит за пределы значений, заданных „коробочкой Солодовникова“. Следовательно, переходный процесс удовлетворяет предъявленным условиям качества регулирования САУ. 1.8 Анализ устойчивости скорректированной САУ Производится по критерию устойчивости Михайлова. Передаточная функция разомкнутой скорректированной САУ имеет вид: Передаточная функция замкнутой скорректированной САУ определяется следующим образом: Раскроем скобки в знаменателе передаточной функции: Заменяем переменную s на jω: Разобьем это выражение на действительную и мнимую составляющие. — вещественная часть; — мнимая часть. По этим данным строится годограф Михайлова. Для устойчивости САУ, необходимо и достаточно, чтобы вектор годографа Михайлова последовательно обошёл вокруг начало координат и в 3 квадранте ушёл в бесконечность. ω ∈ (0 ÷ 300) Рисунок 1.8.1 — Годограф Михайлова для скорректированной системы Таблица 1.8.1 — Данные для построения годографа Михайлова ω 0 10 30 50 60 70 80 100 200 500 Cз (ω) 6.3 6.00 3.52 -1.44 -4.85 -8.88 -13.5 -24.7 -117.7 -768.7 Dз (ω) 0 1.438 4.03 5.75 6.108 6.03 5.46 2.5 -67 -1427 Вектор Михайлова обошел вокруг начала координат и в 3 квадранте ушел в бесконечность. Отсюда следует, что скорректированная САУ устойчива. 2 Исследование нелинейной системы. Согласно заданию, структурная схема нелинейной САУ выглядит следующим образом: Рисунок 2.1 – структурная схема нелинейной системы. Рисунок 2.2 – передаточная характеристика нелинейного звена: с = 2; b = c/K4 = 2/0.1 = 5. 2.1 Построение фазового портрета нелинейной САУ Выполняется вручную методом изоклин при помощи математического пакета MathCAD, в котором производится построение изоклин.. Уравнения изоклин получаем исходя сначала из передаточной функции линейной части системы: Здесь, N – коэффициент угла наклона фазовой траектории при прохождении через изоклину; X – отклонение выходной величины от её заданного значения; X1 (t) – функция, зависящая от свойств нелинейного звена. X1 (t) принимает следующие значения: b при X < –c –b при X > c K4∙X при –c≤X≤c. Рисунок 2.1.1 — Изоклины фазового портрета. Талица 2.1.1 — Данные для построения изоклин. C -100 -50 -37.5 -25 X -40 -20 20 40 -40 -20 20 40 -40 -20 20 40 -40 -20 20 40 Y -137.5 -75 75 137.5 -550 -300 300 550 -2200 -1200 1200 2200 1100 600 -600 -1100 C -20 -15 -10 -5 X -40 -20 20 40 -40 -20 20 40 -40 -20 20 40 -40 -20 20 40 Y 687.5 375 -375 -687.5 500 272.7 -272.7 -500 392.9 214.3 -214.3 -392.9 323.5 176.5 -176.5 -323.5 C 0 5 10 15 X -40 -20 20 40 -40 -20 20 40 -40 -20 20 40 -40 -20 20 40 Y 275 150 -150 -275 239.1 130.4 -130.4 -239.1 211.5 115.4 -115.4 -211.5 189.7 103.4 -103.4 -189.7 C 20 25 50 100 X -40 -20 20 40 -40 -20 20 40 -40 -20 20 40 -40 -20 20 40 Y 171.9 93.8 -93.8 -171.9 157.1 85.7 -85.7 -157.1 110 60 -60 -110 68.8 37.5 -37.5 -68.8 2.2 Оценка устойчивости нелинейной САУ по критерию В.М. Попова. Для определения устойчивости по критерию Попова, по передаточной характеристике линейной части системы строится частотная характеристика вида и определяется возможность проведения хотя бы одной прямой, проходящей через точку с координатой (–1/ k ; j0) , и не пересекающей график частотной характеристики. Коэффициент k здесь — произведение коэффициента усиления линейной части системы и тангенса угла наклона прямой, определяющей класс нелинейности системы. Если такая прямая существует, то система абсолютно устойчива. Передаточная функция линейной части нелинейной САУ: ; Произведём замену переменной s на j∙ω : Разложим получившееся выражение на вещественную и мнимую части: Значение k определяется следующим образом: . Рисунок 2.2.1 — частотная характеристика критерия Попова Таблица 2.2.1 — Данные для построения частотной характеристики к критерию абсолютной устойчивости Попова ω 0 1 1 6 9 20 30 80 200 Up (ω) 2 1.96 1.479 1.03 0.5 -0.166 -0.195 -0.057 -0.01 Vp (ω) 0 -0.315 -0.193 -7.17 -10.59 -11.55 -8.46 -1.956 -0.341 Судя по расположению характеристики и точки с координатой (-1/k ; j0), возможно провести через эту точку прямую, которая не пересечёт частотную характеристику. Следовательно, система абсолютно устойчива.