Скачать .docx | Скачать .pdf |
Курсовая работа: Расчёт одноконтурной системы автоматического регулирования
ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ОБРАЗОВАНИЮ РФ
Государственное образовательное учреждение высшего профессионального образования
«ТОМСКИЙ ПОЛИТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ»
Факультет: Теплоэнергетический
Кафедра: Автоматизации теплоэнергетических процессов
Специальность: 220301 «Автоматизация технологических процессов и производств (в теплоэнергетике)»
Курсовая работа по ТАУ
Расчёт одноконтурной системы автоматического регулирования
Вариант №7
Исполнитель
студент гр.6241: Коростелев А.А.
Руководитель
преподаватель: Татарников А.А.
Томск 2007
Аннотация
В данной курсовой работе представлены расчёт и построение границы заданного запаса устойчивости, одноконтурной АСР с ПИ-регулятором, корневым методом с использованием РАФЧХ. Рассмотрен процесс определения оптимальных параметров настройки регулятора, произведены расчёт и построение переходных процессов в замкнутой АСР при возмущении f, идущем по каналу регулирующего воздействия, и при сигнале задания S. После каждого из графиков данных переходных процессов произведена оценка качества этих процессов.
Содержание
Введение……………………………………………………………………………………….……….4
1. Расчёт оптимальных параметров настройки(ОПН) ………………………………………….…..5
1.1 Расчёт и построение границы заданного запаса устойчивости АСР …………………..……...5
1.2 Обоснование и выбор ОПН регулятора……………………………….…………………….….10
2. Расчёт, построение и оценка качества переходного процесса по каналу S -Y ……………..…11
3. Расчёт, построение и оценка качества переходного процесса по каналу f -Y ……………...…15
Заключение…………………………………………………………………………………………....20
Введение
Данная курсовая работа посвящена расчёту одноконтурной системы автоматического регулирования. Для оценки систем регулирования с точки зрения их практической пригодности необходимо определить, в каких условиях эти системы можно использовать, какие настроечные параметры регулятора требуется установить, чтобы процесс регулирования, осуществляемый при помощи различных регуляторов систем, был оптимальным.
В настоящее время системы регулирования получили широкое применение в различных отраслях промышленности. В связи с этим проблема определения оптимальных параметров настройки регуляторов систем остаётся актуальной, даже несмотря на то, что разработано большое количество приёмов и методов, позволяющих решать эти проблемы. В частности, существует два инженерных метода расчёта систем регулирования: корневой (с использованием РАФЧХ) и частотный по максимуму АЧХ замкнутой системы (метод В.Я. Ротача).
В данной курсовой работе приводятся расчёта заданной АСР, исходные данные и структурная схема которой представлены в задании на выполнение курсовой работы. Первый пункт посвящен расчёту и построению границы заданного запаса устойчивости АСР с ПИ-регулятором и объектом регулирования, корневым методом. А также обоснование и выбор оптимальных параметров настройки. Второй пункт посвящён расчёту переходного процесса по каналу регулирующего воздействия S-Y, и прямой оценки качества этого процесса. Третий пункт содержит расчёт переходного процесса при возмущении f, идущему по каналу воздействия. А также произведены оценки прямых критериев качества.
1. Расчёт оптимальных параметров настройки (ОПН).
1.1 Расчёт и построение границы заданного запаса устойчивости АСР.
Для расчёта и построения границы заданного запаса устойчивости АСР с ПИ-регулятором, представленной на рисунке 1, воспользуемся корневым методом параметрического синтеза систем автоматического регулирования с применением расширенных амплитудно-фазовых частотных характеристик (РАФЧХ).
Используя исходные данные, приведенные в таблице 1, можем записать, что для заданной системы регулирования установлены следующие требования к запасу устойчивости системы: степень затухания переходного процесса в системе.
Исходя из этого можно определить, зависимость между степенью затухания переходных процессов в заданной системе регулирования ψ и степенью колебательности переходных процессов в заданной системе регулирования m, по таблице соответствия оценок запаса устойчивости приведённой ниже.
0 |
0.75 |
0.80 |
0.265 |
0.90 |
0.95 |
0.998 |
1.0 |
|
m |
0 |
0.221 |
0.265 |
0.305 |
0.366 |
0.478 |
1.0 |
∞ |
Эта таблица была получена на основе следующего соотношения:
(1)
где ψ - степенью затухания;
m – степень колебательности;
Передаточная функция объекта регулирования согласно исходных данных определяется по формуле:
(2)
где Р – оператор Лапласа;
К – коэффициент передачи;
При n=2 выражение для примет вид:
(3)
Используя данные таблицы 1 подставляем значения параметров в выражение (3). После подстановки значений параметров получаем окончательное выражение для передаточной функции объекта регулирования:
(4)
Определим расширенные частотные характеристики объекта регулирования. Расширенные частотные характеристики какого-либо звена можно получить подстановкой в передаточную функцию этого звена W(P), оператора или , в выражениях для оператора Лапласа ω – частота, с-1 . В первом случае расчётные формулы метода обеспечивают получение границы заданной степени колебательности системы m, а во втором - получение границы заданной степени устойчивости системы в пространстве параметров настройки регулятора.
Так как заданно значение колебательности, заменяем в формуле (4) оператор , в результате получаем выражение для РАФЧХ объекта регулирования:
(5)
Используя математический пакет MAthCad, предварительно задав начальное значение частоты =0 с-1 и шаг по частоте с-1 , рассчитываем расширенные частотные характеристики объекта при изменении частоты до ω=0,20 с-1 .
Расширенная вещественная частотная характеристика (РВЧХ):
Reоб (m,ω)=Re(Wоб (m,iω)) (6)
Расширенная мнимая частотная характеристика (РМЧХ):
Imоб (m,ω)=Im(Wоб (m,iω)) (7)
Расширенная амплитудно-частотная характеристика (РАЧХ)
(8)
Расширенная фазо-частотная характеристика (РФЧХ):
(9)
Результаты расчётов сведём в таблицу 2, приведенную ниже.
Таблица 2 – Расширенные частотные характеристики объекта регулирования
частота ω, с-1 |
Reоб (m,ω) |
Imоб (m,ω) |
Аоб (m,ω) |
φоб (m,ω), рад |
0,001 |
1.548 |
-0.178 |
1.558 |
-0.114 |
0,003 |
1.562 |
-0.47 |
1.631 |
-0.292 |
0,004 |
1.493 |
-0.772 |
1.681 |
-0.477 |
0,006 |
1.341 |
-1.049 |
1.702 |
-0.664 |
0,007 |
1.118 |
-1.273 |
1.695 |
-0.85 |
0,008 |
0.852 |
-1.425 |
1.648 |
-1.032 |
0.01 |
0.571 |
-1.499 |
1.604 |
-1.207 |
0.011 |
0.301 |
-1.501 |
1.531 |
-1.373 |
Окончание таблицы 2
частота ω, с-1 |
Reоб (m,ω) |
Imоб (m,ω) |
Аоб (m,ω) |
φоб (m,ω), рад |
0.013 |
0.06 |
-1.446 |
1.448 |
-1.529 |
0.014 |
-0.142 |
-1.352 |
1.359 |
-1.675 |
0.016 |
-0.303 |
-1.233 |
1.269 |
-1.812 |
0.017 |
-0.425 |
-1.102 |
1.181 |
-1.938 |
0.019 |
-0.512 |
-0.97 |
1.097 |
-2.057 |
0.021 |
-0.57 |
-0.841 |
1.017 |
-2.166 |
0.022 |
-0.605 |
-0.721 |
0.942 |
-2.269 |
0.024 |
-0.622 |
-0.612 |
0.872 |
-2.364 |
0.025 |
-0.624 |
-0.513 |
0.808 |
-2.454 |
0.027 |
-0.616 |
-0.426 |
0.749 |
-2.537 |
0.028 |
-0.601 |
-0.349 |
0.695 |
-2.615 |
0.03 |
-0.58 |
-0.283 |
0.645 |
-2.688 |
0.031 |
-0.556 |
-0.225 |
0.6 |
-2.757 |
0.033 |
-0.531 |
-0.176 |
0.559 |
-2.822 |
0.034 |
-0.504 |
-0.134 |
0.521 |
-2.883 |
0.036 |
-0.477 |
-0.097 |
0.487 |
-2.94 |
0.037 |
-0.451 |
-0.067 |
0.455 |
-2.995 |
0.039 |
-0.425 |
-0.041 |
0.427 |
-3.046 |
0.04 |
-0.4 |
-0.019 |
0.4 |
-3.095 |
0.042 |
-0.376 |
1.212e-4 |
0.376 |
3.141 |
Расчётные формулы корневого метода для ПИ- регулятора имеют следующий вид:
(10)
(11)
В вышеприведенных формулах (10) и (11) - коэффициент передачи ПИ- регулятора, - постоянная интегрирования ПИ- регулятора или время изодрома.
Зададим диапазон изменения частоты с-1 с шагом c-1 , определим настройки регулятора и Кр в заданном диапазоне частот. Результаты расчётов сведём в таблицу 3.
частота ω, с-1
Таблица 3 –Результаты расчёта настройки ПИ- регулятора в заданном диапазоне частот
w |
Kp |
Kp/Tu |
0,0010 |
-0.611 |
0,0001 |
2.5e-3 |
-0.522 |
0,0005 |
0,0040 |
-0.429 |
0,0013 |
5.5e-3 |
-0.33 |
0,0024 |
0,0070 |
-0.227 |
0,0035 |
8.5e-3 |
-0.12 |
0,0049 |
0.01 |
-8.723e-3 |
0,0066 |
0.011 |
0.106 |
0,0084 |
0.013 |
0.224 |
0.01 |
0.014 |
0.345 |
0.012 |
0.016 |
0.468 |
0.014 |
0.017 |
0.593 |
0.016 |
Окончание таблицы 3
частота ω, с-1 |
Kp |
Kp/Tu |
0.019 |
0.721 |
0.017 |
0.021 |
0.85 |
0.019 |
0.022 |
0.98 |
0.02 |
0.024 |
1.112 |
0.021 |
0.025 |
1.244 |
0.022 |
0.027 |
1.376 |
0.023 |
0.028 |
1.509 |
0.023 |
0.03 |
1.641 |
0.023 |
0.031 |
1.773 |
0.022 |
0.033 |
1.904 |
0.021 |
0.034 |
2.034 |
0.019 |
0.036 |
2.163 |
0.017 |
0.037 |
2.301 |
0.013 |
0.039 |
2.415 |
9.737e-3 |
0.04 |
2.537 |
5.243e-3 |
0.042 |
2.658 |
-4.031e-5 |
По данным таблицы 3 построим график зависимости =f(Kp ) ,т.е укажем границу заданного запаса устойчивости системы регулирования на рисунке 3.
Рисунок 3 - Область параметров настройки ПИ- регулятора
Полученная кривая является линией заданной степени затухания Ψ= Ψзад =0,9 процесса регулирования, что соответствует степени колебательности m=0.366. Таким образом, все значения и Kp , лежащие на этой кривой, обеспечивают определенную степень затухания (в данном случае Ψ= Ψзад =0,9).
Значения и Kp , лежащие внутри области, ограниченной данной кривой и осями координат, обеспечат процесс регулирования со степенью затухания больше заданного (Ψ1 > Ψзад ), а лежащие вне этой области – со степенью затухания меньше заданной (Ψ1 <Ψзад ).
3.2 Обоснование и выбор ОПН регулятора.
Поиск оптимальных параметров настройки регулятора осуществляется вдоль границы заданного запаса устойчивости системы регулирования, представленной на рисунке 3, до достижения экстремума заданного критерия качества. В задании на курсовую работу в качестве принятого критерия качества указан второй интегральный критерий.
Минимуму первого интегрального критерия на графике (рисунок 3) соответствует точка, в которой принимает значение равное 0,95 от максимального в сторону увеличения частоты. Эта точка и определит оптимальные параметры настройки ПИ- регулятора. Используя данные таблицы 3 и рисунка 3, находим, что точке максимума соответствуют значения:
, Kp = 1.509 при ω = 0.028 с-1 .
Поэтому оптимальные параметры настройки ПИ- регулятора имеют значения:
, Kp ∙ 0,95= 1.433 , с.
2. Расчёт, построение, и оценка качества переходного процесса по каналу регулирующего воздействия S - Y
Для одноконтурной системы регулирования, приведенной на рисунке 1, определим передаточную функцию замкнутой АСР по каналу S-Y из соотношения:
(12)
где передаточная функция объекта регулирования ,
передаточная функция ПИ- регулятора .
После подстановки значения в формулу (12), получаем окончательное выражение для передаточной функции замкнутой АСР по каналу S-Y:
(13)
Получим выражение для АФЧХ замкнутой системы путём замены оператора p в формуле (13) на , в результате получаем:
(14)
Используя математический пакет MathCad, предварительно задав диапазон изменения частоты с-1 с шагом c-1 , рассчитываем вещественную частотную характеристику замкнутой АСР при регулирующем воздействии: ReЗ.С..1 (ω). Результаты расчёта сведём в таблицу 4.
Таблица 4 – Результаты расчёта ВЧХ замкнутой АСР при регулирующем воздействии
частота ω, с-1 |
Reоб (m,ω) |
0,001 |
1.001 |
0,009 |
1.041 |
0.017 |
1.043 |
0.025 |
0.369 |
0.033 |
-0.674 |
0.041 |
-0.641 |
0.049 |
-0.452 |
0.057 |
-0.32 |
Продолжение таблицы 4
частота ω, с-1 |
Reоб (m,ω) |
0.065 |
-0.235 |
0.073 |
-0.178 |
0.081 |
-0.138 |
0.089 |
-0.109 |
0.097 |
-0.088 |
0.105 |
-0.071 |
0.113 |
-0.058 |
0.121 |
-0.048 |
0.129 |
-0.04 |
0.137 |
-0.033 |
0.145 |
-0.027 |
0.153 |
-0.023 |
0.161 |
-0.019 |
0.169 |
-0.015 |
0.177 |
-0.012 |
0.185 |
-0,0098 |
0.193 |
-0,0078 |
По данным таблицы 4 строим график ВЧХ замкнутой АСР, который приведен на рисунке 4.
Рисунок 4 - График ВЧХ замкнутой АСР при регулирующем воздействии
Переходный процесс в замкнутой АСР по каналу S-Y можно рассчитать по методу трапеций, используя график ВЧХ замкнутой АСР, приведенный на рисунке 4.
Установлено, что переходная характеристика какой- либо системы y(t) связана с ВЧХ этой системы Re(ω) выражением:
(15)
где t – время переходного процесса в замкнутой АСР.
Для более точного расчёта в качестве верхнего предела интеграла для y(t) принимают не , а значение частоты, при которой график Re(ω) стремится к 0, т.е частоту среза ωСР . По графику, приведенному на рисунке 4, определяем, ωСР =0,2 с-1 . Поэтому переходный процесс в замкнутой АСР по каналу S-Y можно рассчитать по формуле:
(16)
Задав диапазон изменения времени переходного процесса с и шаг с, рассчитываем переходный процесс в замкнутой АСР по каналу S-Y. Результаты расчета сведём в таблицу 5.
Таблица 5 – Результаты расчёта переходного процесса в замкнутой АСР по каналу S-Y
t, c |
Ys-y(t) |
0 |
0,000 |
30 |
0.189 |
60 |
0.722 |
90 |
1.169 |
120 |
1,330 |
150 |
1.245 |
180 |
1,070 |
210 |
0.937 |
240 |
0.898 |
270 |
0.931 |
300 |
0.986 |
330 |
1.024 |
360 |
1.032 |
390 |
1,020 |
420 |
1.002 |
450 |
0.992 |
480 |
0.99 |
510 |
0.994 |
540 |
1,000 |
570 |
1.003 |
600 |
1.003 |
630 |
1.001 |
660 |
1,000 |
690 |
0.999 |
По данным таблицы 5 строим график переходного процесса в замкнутой АСР по
каналу S-Y, который приведён на рисунке 5.
Рисунок 5 - График переходного процесса в замкнутой АСР по каналу S-Y
Используя данные таблицы 5 и рисунка 5, произведём оценку качества переходного процесса в замкнутой АСР по каналу S-Y.
Прямые критерии качества:
1.Максимальная динамическая ошибка: А1 =0,34;
2.Перерегулирование: (17)
где - уровень установившегося значения регулируемой величины при времени переходного процесса , равного ;
3.Динамический коэффициент регулирования Rд не определяется для такого типа процессов;
4.Степень затухания переходного процесса: (18)
где - второй максимальный выброс регулируемой величины;
5.Статическая ошибка: (19)
где S – сигнал регулирующего воздействия 1(t);
6.Время регулирования: при величине , значение которой задают для контроля переходного процесса с заданной степенью точности.
Все приведенные выше критерии качества указаны на рисунке 5.
3. Расчёт, построение и оценка качества переходного процесса по каналу f - Y
Для одноконтурной системы регулирования, приведенной на рисунке 1, определим передаточную функцию замкнутой АСР по каналу f -Y по формуле:
(20)
После подстановки выражения для в формулу (7), получаем окончательное выражение для передаточной функции замкнутой АСР по каналу f -Y:
(21)
Получим выражение для АФЧХ замкнутой системы путём замены оператора p в формуле (18) на , в результате получаем:
(22)
Используя математический пакет MAthCad, предварительно задав диапазон изменения частоты с-1 с шагом c- , рассчитываем вещественную частотную характеристику замкнутой АСР: ReЗ.С.2 (ω). Результаты расчёта сведём в таблицу 6.
Таблица 6 - Результаты расчёта ВЧХ замкнутой АСР при возмущении f
частота ω, с-1 |
Reоб (m,ω) |
0.01 |
0.315 |
0.02 |
0.772 |
0.03 |
0.158 |
0.04 |
-0.306 |
0.05 |
-0.256 |
0.06 |
-0.185 |
0.07 |
-0.135 |
0.08 |
-0.101 |
0.09 |
-0.077 |
0.1 |
-0.06 |
0.11 |
-0.047 |
0.12 |
-0.037 |
0.13 |
-0.03 |
0.14 |
-0.024 |
0.15 |
-0.019 |
0.16 |
-0.015 |
0.17 |
-0.012 |
0.18 |
-0,0093 |
0.19 |
-0,0071 |
0.2 |
-0,054 |
По данным таблицы 6 строим график ВЧХ замкнутой АСР при возмущении f, который приведен на рисунке 6
Рисунок 6 – График ВЧХ замкнутой АСР при возмущении f
Переходный процесс в замкнутой АСР по каналу f-Y можно рассчитать по методу трапеций, используя график ВЧХ замкнутой АСР при возмущении f (рисунок 6).
Поэтому переходный процесс в замкнутой АСР по каналу F-Y можно рассчитать по формуле:
(23)
Как уже было сказано выше, для более точного расчёта в качестве верхнего предела интеграла для YF - Y (t) принимают значение частоты среза ωСР. По графику, приведенному на рисунке 6, определяем, что ωСР =0,2 с-1 .
Задав диапазон изменения времени переходного процесса с и шаг с, рассчитываем переходный процесс в замкнутой АСР по каналу f-Y. Результаты расчета сведём в таблицу 7, приведенную ниже.
Таблица 7 - Результаты расчёта переходного процесса в замкнутой АСР по каналу f-Y
t, c |
Ys-y(t) |
0 |
0 |
30 |
0.116 |
60 |
0.37 |
90 |
0.472 |
120 |
0.374 |
150 |
0.181 |
180 |
0.014 |
210 |
-0.063 |
240 |
-0.059 |
270 |
-0.017 |
300 |
0.02 |
330 |
0.033 |
360 |
0.024 |
390 |
0,0034 |
420 |
-0,0047 |
450 |
-0,0087 |
480 |
-0,0067 |
510 |
-0,013 |
540 |
0,021 |
570 |
0,034 |
600 |
0,018 |
630 |
0,0024 |
60 |
-0,0065 |
690 |
-0,0093 |
По данным таблицы 7 строим график переходного процесса в замкнутой АСР по каналу f-Y, представленный на рисунке 7.
Рисунок 7 - График переходного процесса в замкнутой АСР по каналу f-Y
Используя данные таблицы 7 и рисунка 7, произведём оценку качества переходного процесса в замкнутой АСР по каналу F-Y.
Прямые критерии качества:
1.Максимальная динамическая ошибка: А1=0,47;
2.Перерегулирование: (24)
где - первое минимальное отклонение регулируемой величины;
3.Динамический коэффициент регулирования RД : (25)
где - коэффициент передачи объекта;
4.Степень затухания переходного процесса: ; (26)
5.Статическая ошибка: ;
6.Время регулирования: при величине .
Так как в заданной АСР, представленной на рисунке 2, имеется звено чистого транспортного запаздывания с передаточной функцией , то переходные процессы в этой системе имеет запаздывание на величину 4 с относительно их начала. Для наглядности указанного факта изобразим начальные части графиков переходных процессов по каналам S-Y и f-Y соответственно на рисунке 8 и 9.
Рисунок 8 – Начальный участок графика переходного процесса в замкнутой АСР по каналу S-Y
Рисунок 9 – Начальный участок графика переходного процесса в замкнутой АСР по каналу f-Y
Заключение
Определение оптимальных параметров настройки регуляторов, расчёт различных систем автоматического регулирования, без сомнения, являются одной из главных задач любого инженера. Использование современных систем регулирования требует знания различных методов и приёмов расчёта этих систем, определения и установки требуемых параметров настройки регулятора, основных недостатков и преимуществ разного рода регуляторов по сравнению друг с другом.
В процессе написания курсовой работы был изучен один из двух инженерных методов расчёта одноконтурных систем регулирования: корневой метод (с использованием РАФЧХ). Было выяснено, что оптимальными параметрами настройки какого-либо регулятора считают те параметры, при которых обеспечивается близкий к оптимальному процесс регулирования. Под оптимальным процессом регулирования обычно понимают процесс, удовлетворяющий требованиям к запасу устойчивости системы. Поиск оптимальных параметров настройки осуществляется вдоль границы заданного запаса устойчивости системы регулирования до достижения экстремума принятого критерия качества. В данной курсовой работе, согласно заданию, был принят второй интегральный критерий.
В результате проделанной работы, были получены переходные процессы по каналам S-Y и f-Y. Оценка качества этих процессов показала, что они удовлетворяют требованиям к запасу устойчивости системы, приведенных в исходных данных.
Можно заметить, что переходный процесс по каналу f-Y имеет прямые критерии качества лучше, чем переходный процесс по каналу S-Y:
Таблица 8 – Прямые критерии качества переходных процессов по каналам S-Y и f-Y
S-Y |
f-Y |
|
Максимальная динамическая ошибка |
0,34 |
0.47 |
Перерегулирование (%) |
34 |
14 |
Степень затухания переходного процесса |
0,88 |
0,91 |
Время регулирования tp , с |
270 |
175 |
Статическая ошибка для этих процессов |
0 |
0 |
Следовательно регулятор установленный в канале обратной связи способствует лучшей работе системы нежели он будет установлен в основном канале.