Федеральное агентство по образованию РФ
ГОУ ВПО
"Московский Государственный Текстильный Университет им. А.Н. Косыгина"
ФИТАЭ
Кафедра автоматики и промышленной электроники
Курсовая работа
Тема: "Синтез астатических систем"
Дисциплина: "Теория Автоматического Управления"
Москва 2009
Исходные данные:
ПИ-закон
1. По заданной системе уравнений получить передаточную функцию объекта управления и составить структурную схему замкнутой САУ, считая регулятор звеном с входной величиной Е, выходной U и передаточной функцией
2. Выбрать регулятор таким образом, чтобы система обладала свойством астатизма по отношению к постоянному задающему g(t) и возмущающему f(t) воздействию.
Для того чтобы система обладала астатизмом необходимо выполнение условия – величина установившейся ошибки по заданию и возмущению должна быть равна 0. Проверим это условие на регуляторах:
· Выбираем П – регулятор
Т.к. величина установившейся ошибки не равна нулю, то система не обладает астатизмом,, а следовательно считать ошибку по каналу f-E нет необходимости.
· Выбираем И – регулятор
Обе ошибки равны нулю следовательно система с И - регулятором является астатической. Данный регулятор нам подходит
· Проверим наш ПИ – регулятор
Система так же является астатической, что является тем, что выбранный нами ПИ – регулятор подходит к условию задания.
3. Найти область значений постоянной времени регулятора для И – закона управления, обеспечивающих устойчивость системы.
Найдем передаточную функцию замкнутой системы по каналу g-y и затем выделим ее характеристическое уравнение.
Характеристическое уравнение Q(p):
Уравнение 3 порядка, следовательно, для устойчивости системы необходимо чтобы произведение внутренних коэффициентов характеристического уравнения было больше произведения внешних коэффициентов.
Следовательно область значений для И – закона можно определить интервалом
.
4. Для И – закона регулирования найти минимальную возможную величину установившейся ошибки, если g(t)=2t и f(t)=-3t.
Для определения минимальной установившейся ошибки нам необходимо узнать при каком значении постоянной времени система находится на границе устойчивости. Для этого отыщем передаточную функцию и характеристическое уравнение системы с И – регулятором. И затем найдем значение
Находим значение постоянной времени на границе устойчивости:
Найдем величину установившейся ошибки при g(t)=2t и f(t)=-3t
Поскольку задающее воздействие у нас g(t)=2t, то используя преобразование Лапласа получаем:
тогда величина установившейся ошибки будет
Подставляем полученное значение и получаем
Найдем ошибку по каналу f-E
Подставляем и получаем
Тогда
5. Построить, с использованием ЭВМ, область устойчивости, на плоскости, параметров регулятора при использовании ПИ – закона, обосновать возможность и путь получения допустимой установившейся ошибки при и
.
Для построение области устойчивости необходимо найти характеристическое уравнение передаточной функции для данного регулятора.
Отсюда:
Запишем условие, при котором система находится на границе устойчивости:
Выразим зависимость
от
Строим область устойчивости по по лученной зависимости:
Ти | Кп |
0,01 | 4,00 |
0,02 | 1,50 |
0,03 | 0,67 |
0,04 | 0,25 |
0,05 | 0,00 |
0,06 | -0,17 |
0,07 | -0,29 |
0,08 | -0,38 |
0,09 | -0,44 |
Таким образом выбирая значения параметров регулятора над границей устойчивости – мы получаем устойчивую систему, и наоборот.
6. Вычислить значения параметров ПИ – регулятора, обеспечивающих устойчивость и установившуюся ошибку в системе = 0,06 при g(t)=2t и f(t)=0
Поскольку возмущение f(t)=0, то . Найдем
:
для этого найдем передаточную функцию замкнутой системы по каналу g-E
По условию , тогда подставим это значение в получившееся выражение:
Таким образом для получения в системе установившейся ошибке равной 0,06 необходимо задать параметру постоянной времени значение 0,03.
7. Для интегрального закона регулирования и начальных условий , выбрать оптимальное значение постоянной времени регулятора по критерию:
(Рассматривается движение в системе при g(t)=f(t)=0 и ненулевых начальных условиях).
Для решения мы будем использовать метод Мондельштама. Для этого нам необходимо найти передаточную функцию замкнутой системы и взять характеристическое уравнение:
Получаем уравнение:
Поочередно умножаем его на и на Е
Интегрируем полученное уравнение по частям
Получаем:
Интегрируем полученное уравнение по частям
Получаем:
Из полученных уравнений составим систему уравнений:
Выбираем =1,054
8. Для найденного в пункте 7 значения постоянной времени регулятора построить с помощью ЭВМ вещественную частотную характеристику P(ω), приняв входным воздействием g(t) и входной координатой E(t)/
Для нахождения вещественной характеристики нам понадобится передаточная функция замкнутой системы по каналу g-y.
Перейдем в частотную область p=jω:
Домножаем на сопряженное знаменателю число и получаем:
Отделяем действительную часть U(ω):
При =1,054
w | p(w) |
0 | 1 |
0,1 | 1,010058 |
0,2 | 1,039191 |
0,3 | 1,079471 |
0,4 | 1,099022 |
0,5 | 0,997092 |
0,6 | 0,58593 |
0,7 | -0,06976 |
0,8 | -0,48243 |
0,9 | -0,56794 |
1 | -0,5208 |
1,1 | -0,44696 |
1,2 | -0,3782 |
1,3 | -0,32081 |
1,4 | -0,27428 |
1,5 | -0,23666 |
1,6 | -0,20606 |
1,7 | -0,18095 |
1,8 | -0,16013 |
1,9 | -0,14269 |
2 | -0,12796 |
2,5 | -0,08003 |
3 | -0,05481 |
3,5 | -0,03991 |
4 | -0,03037 |
4,5 | -0,02389 |
5 | -0,01929 |
5,5 | -0,0159 |
6 | -0,01334 |
9. По вещественной характеристике P(ω) пункта 8 построить переходной процесс E(t) при единичном ступенчатом изменении g(t) и нулевых начальных условиях методом трапециидальных частотных характеристик.
Для построения переходного процесса нам необходимо разбить получившуюся вещественную характеристику на трапеции и построить переходный процесс для каждой из полученных трапеций.
R(0) | Wo | Wd | æ | ||
I | - | 0,099021688 | 0,38 | 0,1 | 0,263158 |
II | + | 1,666965285 | 0,88 | 0,43 | 0,488636 |
III | - | 0,567943597 | 6 | 0,95 | 0,158333 |
1-я трапеция
t табл | h(æ) | t=t табл/Wo | h=R(0)*h(æ) |
0 | 0 | 0 | 0 |
0,5 | 0,199 | 1,315789474 | -0,019705316 |
1 | 0,386 | 2,631578947 | -0,038222372 |
1,5 | 0,56 | 3,947368421 | -0,055452145 |
2 | 0,709 | 5,263157895 | -0,070206377 |
2,5 | 0,833 | 6,578947368 | -0,082485066 |
3 | 0,928 | 7,894736842 | -0,091892127 |
3,5 | 0,994 | 9,210526316 | -0,098427558 |
4 | 1,039 | 10,52631579 | -0,102883534 |
4,5 | 1,057 | 11,84210526 | -0,104665924 |
5 | 1,067 | 13,15789474 | -0,105656141 |
5,5 | 1,067 | 14,47368421 | -0,105656141 |
6 | 1,054 | 15,78947368 | -0,104368859 |
6,5 | 1,043 | 17,10526316 | -0,103279621 |
7 | 1,035 | 18,42105263 | -0,102487447 |
7,5 | 1,025 | 19,73684211 | -0,10149723 |
8 | 1,024 | 21,05263158 | -0,101398209 |
8,5 | 1,022 | 22,36842105 | -0,101200165 |
9 | 1,025 | 23,68421053 | -0,10149723 |
9,5 | 1,027 | 25 | -0,101695274 |
10 | 1,027 | 26,31578947 | -0,101695274 |
10,5 | 1,028 | 27,63157895 | -0,101794295 |
11 | 1,025 | 28,94736842 | -0,10149723 |
11,5 | 1,021 | 30,26315789 | -0,101101144 |
12 | 1,015 | 31,57894737 | -0,100507013 |
12,5 | 1,01 | 32,89473684 | -0,100011905 |
13 | 1,005 | 34,21052632 | -0,099516797 |
13,5 | 1 | 35,52631579 | -0,099021688 |
14 | 0,997 | 36,84210526 | -0,098724623 |
14,5 | 0,996 | 38,15789474 | -0,098625601 |
15 | 0,995 | 39,47368421 | -0,09852658 |
15,5 | 0,995 | 40,78947368 | -0,09852658 |
16 | 0,995 | 42,10526316 | -0,09852658 |
16,5 | 0,995 | 43,42105263 | -0,09852658 |
17 | 0,995 | 44,73684211 | -0,09852658 |
17,5 | 0,995 | 46,05263158 | -0,09852658 |
18 | 0,995 | 47,36842105 | -0,09852658 |
18,5 | 0,994 | 48,68421053 | -0,098427558 |
19 | 0,992 | 50 | -0,098229515 |
19,5 | 0,991 | 51,31578947 | -0,098130493 |
20 | 0,991 | 52,63157895 | -0,098130493 |
2-я трапеция
t табл | h(æ) | t=t табл/Wo | h=R(0)*h(æ) | |
0 | 0 | 0 | 0 | |
0,5 | 0,24 | 0,568181818 | 0,400071669 | |
1 | 0,461 | 1,136363636 | 0,768470997 | |
1,5 | 0,665 | 1,704545455 | 1,108531915 | |
2 | 0,833 | 2,272727273 | 1,388582083 | |
2,5 | 0,967 | 2,840909091 | 1,611955431 | |
3 | 1,061 | 3,409090909 | 1,768650168 | |
3,5 | 1,115 | 3,977272727 | 1,858666293 | |
4 | 1,142 | 4,545454545 | 1,903674356 | |
4,5 | 1,138 | 5,113636364 | 1,897006495 | |
5 | 1,118 | 5,681818182 | 1,863667189 | |
5,5 | 1,092 | 6,25 | 1,820326092 | |
6 | 1,051 | 6,818181818 | 1,751980515 | |
6,5 | 1,018 | 7,386363636 | 1,696970661 | |
7 | 0,993 | 7,954545455 | 1,655296528 | |
7,5 | 0,974 | 8,522727273 | 1,623624188 | |
8 | 0,966 | 9,090909091 | 1,610288466 | |
8,5 | 0,966 | 9,659090909 | 1,610288466 | |
9 | 0,97 | 10,22727273 | 1,616956327 | |
9,5 | 0,975 | 10,79545455 | 1,625291153 | |
10 | 0,982 | 11,36363636 | 1,63695991 | |
10,5 | 0,987 | 11,93181818 | 1,645294737 | |
11 | 0,993 | 12,5 | 1,655296528 | |
11,5 | 0,997 | 13,06818182 | 1,66196439 | |
12 | 0,997 | 13,63636364 | 1,66196439 | |
12,5 | 0,997 | 14,20454545 | 1,66196439 | |
13 | 0,997 | 14,77272727 | 1,66196439 | |
13,5 | 0,998 | 15,34090909 | 1,663631355 | |
14 | 1 | 15,90909091 | 1,666965285 | |
14,5 | 1,002 | 16,47727273 | 1,670299216 | |
15 | 1,005 | 17,04545455 | 1,675300112 | |
15,5 | 1,008 | 17,61363636 | 1,680301008 | |
16 | 1,011 | 18,18181818 | 1,685301904 | |
16,5 | 1,011 | 18,75 | 1,685301904 | |
17 | 1,012 | 19,31818182 | 1,686968869 | |
17,5 | 1,009 | 19,88636364 | 1,681967973 | |
18 | 1,008 | 20,45454545 | 1,680301008 | |
18,5 | 1,006 | 21,02272727 | 1,676967077 | |
19 | 1,001 | 21,59090909 | 1,668632251 | |
19,5 | 0,998 | 22,15909091 | 1,663631355 | |
20 | 0,996 | 22,72727273 | 1,660297424 | |
3-я трапеция
t табл | h(æ) | t=t табл/Wo | h=R(0)*h(æ) |
0 | 0 | 0,0000 | 0,0000 |
0,5 | 0,184 | 0,0833 | -0,1045 |
1 | 0,256 | 0,1667 | -0,1454 |
1,5 | 0,516 | 0,2500 | -0,2931 |
2 | 0,655 | 0,3333 | -0,3720 |
2,5 | 0,833 | 0,4167 | -0,4731 |
3 | 0,863 | 0,5000 | -0,4901 |
3,5 | 0,928 | 0,5833 | -0,5271 |
4 | 0,974 | 0,6667 | -0,5532 |
4,5 | 0,977 | 0,7500 | -0,5549 |
5 | 1,012 | 0,8333 | -0,5748 |
5,5 | 1,019 | 0,9167 | -0,5787 |
6 | 1,013 | 1,0000 | -0,5753 |
6,5 | 1,009 | 1,0833 | -0,5731 |
7 | 1,006 | 1,1667 | -0,5714 |
7,5 | 1,006 | 1,2500 | -0,5714 |
8 | 1,008 | 1,3333 | -0,5725 |
8,5 | 1,01 | 1,4167 | -0,5736 |
9 | 1,016 | 1,5000 | -0,5770 |
9,5 | 1,022 | 1,5833 | -0,5804 |
10 | 1,025 | 1,6667 | -0,5821 |
10,5 | 1,028 | 1,7500 | -0,5838 |
11 | 1,029 | 1,8333 | -0,5844 |
11,5 | 1,027 | 1,9167 | -0,5833 |
12 | 1,025 | 2,0000 | -0,5821 |
12,5 | 1,022 | 2,0833 | -0,5804 |
13 | 1,019 | 2,1667 | -0,5787 |
13,5 | 1,017 | 2,2500 | -0,5776 |
14 | 1,016 | 2,3333 | -0,5770 |
14,5 | 1,015 | 2,4167 | -0,5765 |
15 | 1,014 | 2,5000 | -0,5759 |
15,5 | 1,014 | 2,5833 | -0,5759 |
16 | 1,014 | 2,6667 | -0,5759 |
16,5 | 1,014 | 2,7500 | -0,5759 |
17 | 1,013 | 2,8333 | -0,5753 |
17,5 | 1,012 | 2,9167 | -0,5748 |
18 | 1,011 | 3,0000 | -0,5742 |
18,5 | 1,009 | 3,0833 | -0,5731 |
19 | 1,008 | 3,1667 | -0,5725 |
19,5 | 1,006 | 3,2500 | -0,5714 |
20 | 1,005 | 3,3333 | -0,5708 |
20,5 | 1,004 | 3,4167 | -0,5702 |
21 | 1,003 | 3,5000 | -0,5696 |
21,5 | 1,003 | 3,5833 | -0,5696 |
22 | 1,002 | 3,6667 | -0,5691 |
22,5 | 1,002 | 3,7500 | -0,5691 |
23 | 1,002 | 3,8333 | -0,5691 |
23,5 | 1,002 | 3,9167 | -0,5691 |
24 | 1,001 | 4,0000 | -0,5685 |
24,5 | 1 | 4,0833 | -0,5679 |
25 | 1 | 4,1667 | -0,5679 |
25,5 | 0,999 | 4,2500 | -0,5674 |
26 | 0,999 | 4,3333 | -0,5674 |
Суммируем графически полученные процесы и получаем
10. Определить время регулирования и максимальное перерегулирование в системе.
11. Разработать и начертить структурную схему комбинированной САУ, инвариантной по отношению к контролируемому возмущению .
Привести передаточную функцию устройства управления.
Проверить выполнение условия инвариантности.
Условия инвариантности:
, если
, если
При выборе корректирующих звеньев в виде обратных передаточных функций возникает проблема. Она обычно связана с тем, что порядок числителя корректирующего звена становится больше порядка знаменателя. Это означает, что частотные характеристики таких звеньев являются расходящимися, что говорит о том, что физически такие звенья нереализуемы. В тех случаях, когда корректирующие звенья применять необходимо, порядок числителя этих звеньев искусственно приравнивают к порядку знаменателя, отбрасывая в числители высшие порядки.
12. Предложить и обосновать методы компенсации действия неконтролируемых возмущений, если известен класс функций, которыми они описываются.
астатический автоматический управление моделирование программа
Решить проблему инвариантности можно, если известна предварительная информация о классе возмущающих воздействий. В частности, если известен математический аппарат, описывающий функцию f(t), заданный в виде решения некоторого дифференциального уравнения.
Процедура синтеза предусматривает:
1. восстановление вида дифференциального уравнения по заданному решению;
2. получение характеристического уравнения;
3. выбор передаточной функции регулятора, в которой знаменатель совпадает с видом полученного характеристического уравнения. Числитель передаточной функции регулятора выбирается того же порядка, что и знаменатель;
4. неизвестные коэффициенты числителя передаточной функции регулятора определяются из условий устойчивости замкнутой системы.
13. Провести моделирование в программе MatLab. Определить настройки регулятора, обеспечивающего минимизацию времени регулирования.
При нулевом задающем воздействии со значением регулятора, полученными в 7 пункте:
При единичном задающем воздействии:
Для снижения времени регулирования можно немного увеличить значения регулятора примерно до 1,085.
Так же значительно уменьшает время регулирование и введение пропорциональной составляющей, т.е. использование ПИ – регулятора. С его помощью легко можно уменьшить время регулирования примерно в 1,7 раза.
Список используемой литературы
· Лекции по курсу ТАУ, Румянцев Ю.Д.
· "Теория автоматического управления", Воронов А.А.
· "Теория систем автоматического управления", Бесекерский В.А.
Приложение