Министерство образования и науки Российской Федерации
Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования
"Национальный исследовательский
Томский политехнический университет"
Наименование института ЭНИН
Наименование специальности Теплоэнергетика и теплотехника
Наименование
выпускающей кафедры Автоматизация теплоэнергетических процессов
Курсовая работа
Расчёт
одноконтурной системы автоматического регулирования
Исполнитель
студент группы 5Б1В
Броцман А.А.
Руководитель
Андык В.С.
Томск
- 2014
Содержание
Введение
. Структурная схема одноконтурной АСР
. Расчёт и построение границы заданного запаса устойчивости АСР
. Определение оптимальных параметров настройки ПИ-регулятора
. Расчёт, построение и оценка качества переходных процессов
.1 Переходный процесс по каналу регулирующего воздействия S-Y
.2. Переходный процесс при возмущении f, идущем по каналу регулирующего воздействия F-Y
Заключение
Список использованных источников
регулирование автоматизированный
регулятор устойчивость
Введение
Данная курсовая работа посвящена расчёту одноконтурной системы автоматического регулирования. Для оценки систем регулирования необходимо определить, в каких условиях эти системы можно использовать, какие настроечные параметры регулятора требуется установить, чтобы процесс регулирования, осуществляемый при помощи различных регуляторов систем, был оптимальным.
В настоящее время системы регулирования получили широкое применение в различных отраслях промышленности. В связи с этим проблема определения оптимальных параметров настройки регуляторов систем остаётся актуальной, даже несмотря на то, что разработано большое количество приёмов и методов, позволяющих решать эти проблемы. В частности, существует два инженерных метода расчёта систем регулирования: корневой (с использованием РАФЧХ) и частотный по максимуму АЧХ замкнутой системы (метод В.Я. Ротача).
В данной курсовой работе в первом пункте
приводятся исходные данные для расчёта заданной АСР, структурная схема которой показана
во втором пункте. Третий пункт посвящен расчёту и построению границы заданного
запаса устойчивости АСР с ПИ-регулятором и объектом регулирования корневым
методом. В четвёртом пункте определяются оптимальные параметры настройки
ПИ-регулятора. Последний пункт посвящён расчётам переходных процессов по двум
каналам: по каналу регулирующего воздействия S-Y и при возмущении f, идущем по
каналу регулирующего воздействия, здесь же представлены графики этих процессов
и произведены оценки их качества.
. Структурная схема одноконтурной АСР
Структурная схема системы регулирования,
приведенная в задании, имеет вид:
Рисунок 1 - Структурная схема
заданной системы регулирования
Учитывая исходные данные, приведенные ранее структурную
схему системы регулирования можно преобразовать к виду:
Рисунок 2 - Преобразованная структурная схема
заданной системы регулирования
. Расчёт и построение границы заданного запаса
устойчивости АСР
Для расчёта и построения границы заданного запаса устойчивости АСР с ПИ-регулятором, представленной на рисунке 1, воспользуемся корневым методом параметрического синтеза систем автоматического регулирования с применением расширенных амплитудно-фазовых частотных характеристик (РАФЧХ).
Используя исходные данные, приведенные в таблице
1, можем записать, что для заданной системы регулирования установлены следующие
требования к запасу устойчивости системы: степень затухания переходного
процесса в системе
.
Исходя из этого, зная зависимость между степенью
затухания переходных процессов в заданной системе регулирования
ψ и
степенью колебательности переходных процессов в заданной системе регулирования
m, можно определить значение заданной степени колебательности m системы по
формуле:
(1)
где ψ - степень затухания переходных процессов в заданной системе регулирования.
Передаточная функция объекта регулирования
согласно исходных данных определяется по формуле:
(2)
где Р - оператор Лапласа.
При n=1 выражение
примет
вид:
(3)
По данным таблицы 1 определяем значения
неизвестных параметров: К=1,
, T1=40. Тогда
после подстановки значений выше приведенных параметров получаем окончательное
выражение для передаточной функции объекта регулирования:
(4)
Определим расширенные частотные характеристики
объекта регулирования. Расширенные частотные характеристики какого-либо звена
можно получить подстановкой в передаточную функцию этого звена W(P) оператора
или
,
в выражениях для оператора Лапласа ω - частота,
с-1. В первом случае расчётные формулы метода обеспечивают получение границы
заданной степени колебательности системы m, а во втором - получение границы
заданной степени устойчивости системы
в
пространстве параметров настройки регулятора.
Заменим в формуле (4) оператор
,
в результате получаем выражение для РАФЧХ объекта регулирования:
(5)
Используя математический пакет MathCad,
предварительно задав начальное значение частоты
=0
с-1 и шаг по частоте
с-1, рассчитываем
расширенные частотные характеристики объекта при изменении частоты до ω=0,055
с-1.
Расширенная вещественная частотная характеристика (РВЧХ):
об(m,ω)=Re(Wоб(m,iω))
(6)
Расширенная мнимая частотная характеристика (РМЧХ):
об(m,ω)=Im(Wоб(m,iω))
(7)
Расширенная амплитудно-частотная характеристика
(РАЧХ)
(8)
Результаты расчётов сведём в таблицу 2,
приведенную ниже.
Таблица 2 - Расширенные частотные характеристики объекта регулирования
|
частота ω, с-1 |
Reоб(m,ω) |
Imоб(m,ω) |
Аоб(m,ω) |
|
1,00E-09 |
1 |
-4.8e-8 |
1 |
|
0.01 |
0.9 |
-0.484 |
1.022 |
|
0.02 |
0.539 |
-0.724 |
0.903 |
|
0.03 |
0.228 |
-0.714 |
0.749 |
|
0.04 |
0.04 |
-0.621 |
0.622 |
|
0.05 |
-0.067 |
-0.522 |
0.526 |
|
0.06 |
-0.129 |
-0.436 |
0.455 |
|
0.07 |
-0.165 |
-0.365 |
0.4 |
|
0.08 |
-0.187 |
-0.306 |
0.358 |
|
0.09 |
-0.2 |
-0.256 |
0.325 |
|
0.1 |
-0.208 |
-0.214 |
0.298 |
|
0.11 |
-0.212 |
-0.177 |
0.276 |
|
0.12 |
-0.213 |
-0.145 |
0.258 |
|
0.13 |
-0.212 |
-0.116 |
0.242 |
|
0.14 |
-0.21 |
-0.091 |
0.229 |
|
0.15 |
-0.206 |
-0.068 |
0.217 |
|
0.16 |
-0.202 |
-0.046 |
0.207 |
|
0.17 |
-0.196 |
-0.027 |
0.198 |
|
0.18 |
-0.19 |
-9.296e-3 |
0.19 |
|
0.19 |
-0.183 |
7.038e-3 |
0.183 |
|
0.2 |
-0.176 |
0.022 |
0.177 |
|
0.21 |
-0.168 |
0.036 |
0.172 |
|
0.22 |
-0.159 |
0.049 |
0.167 |
|
0.23 |
-0.151 |
0.06 |
0.162 |
|
0.24 |
-0.141 |
0.071 |
0.158 |
|
0.25 |
-0.132 |
0.081 |
0.154 |
|
0.26 |
-0.122 |
0.089 |
0.151 |
|
0.27 |
-0.112 |
0.097 |
0.148 |
|
0.28 |
-0.101 |
0.104 |
0.145 |
|
0.29 |
-0.091 |
0.11 |
0.143 |
|
0.3 |
-0.08 |
0.115 |
0.14 |
|
0.31 |
-0.069 |
0.12 |
0.138 |
|
0.32 |
-0.058 |
0.123 |
0.136 |
|
0.33 |
-0.047 |
0.126 |
0.134 |
|
0.34 |
-0.036 |
0.128 |
0.133 |
|
0.35 |
-0.025 |
0.129 |
0.131 |
Расчётные формулы корневого метода для
ПИ-регулятора имеют следующий вид:
(10)
(11)
В вышеприведенных формулах (10) и (11)
-
коэффициент передачи ПИ-регулятора,
- постоянная интегрирования
ПИ-регулятора.
Зададим диапазон изменения частоты
с-1
с шагом
c-1,
определим настройки регулятора
и Кр в заданном
диапазоне частот. Результаты расчётов сведём в таблицу 3.
Таблица 3 -Результаты расчёта настройки ПИ- регулятора в заданном диапазоне частот
|
частота ω, с-1 |
Кр/Ти |
Кр |
|
1,00E-09 |
0 |
-1 |
|
0.01 |
4.859e-3 |
-0.759 |
|
0.02 |
0.019 |
-0.465 |
|
0.03 |
0.04 |
-0.125 |
|
0.04 |
0.067 |
0.252 |
|
0.05 |
0.099 |
0.659 |
|
0.06 |
0.133 |
1.088 |
|
0.07 |
0.167 |
1.532 |
|
0.08 |
0.2 |
1.981 |
|
0.09 |
0.229 |
2.429 |
|
0.1 |
0.252 |
2.868 |
|
0.11 |
0.268 |
3.291 |
|
0.12 |
0.275 |
3.692 |
|
0.13 |
0.271 |
4.063 |
|
0.14 |
0.255 |
4.4 |
|
0.15 |
0.226 |
4.697 |
|
0.16 |
0.182 |
4.95 |
|
0.17 |
0.123 |
5.154 |
|
0.18 |
0.048 |
5.306 |
|
0.19 |
-0.042 |
5.403 |
|
0.2 |
-0.147 |
5.444 |
|
0.21 |
-0.268 |
5.426 |
|
0.22 |
-0.404 |
5.351 |
|
0.23 |
-0.553 |
5.216 |
|
0.24 |
-0.714 |
5.023 |
|
0.25 |
-0.887 |
4.773 |
|
0.26 |
-1.068 |
4.469 |
|
0.27 |
-1.257 |
4.112 |
|
0.28 |
-1.451 |
3.705 |
|
0.29 |
-1.648 |
3.252 |
|
0.3 |
-1.845 |
2.757 |
|
0.31 |
-2.04 |
2.224 |
|
0.32 |
-2.23 |
1.657 |
|
0.33 |
-2.413 |
1.063 |
|
0.34 |
-2.586 |
0.445 |
|
0.35 |
-2.747 |
-0.19 |
По данным таблицы 3 построим график зависимости
,
т.е. укажем границу заданного запаса устойчивости системы регулирования на
рисунке 3.
Рисунок 3 - Область параметров настройки
ПИ-регулятора
Полученная кривая является линией заданной
степени затухания Ψ=Ψзад=0,75
процесса регулирования, что соответствует степени колебательности m=0.221.
Таким образом, все значения
и Kp, лежащие на
этой кривой, обеспечивают определенную степень затухания (в данном случае Ψ=Ψзад=0,75).
Значения
и
Kp, лежащие внутри области, ограниченной данной кривой и осями координат,
обеспечат процесс регулирования со степенью затухания больше заданного (Ψ1>Ψзад),
а лежащие вне этой области - со степенью затухания меньше заданной (Ψ1<Ψзад).
3. Определение оптимальных параметров настройки
ПИ-регулятора
Поиск оптимальных параметров настройки регулятора осуществляется вдоль границы заданного запаса устойчивости системы регулирования, представленной на рисунке 3, до достижения экстремума принятого критерия качества. В задании на курсовую работу в качестве принятого критерия качества указан первый интегральный критерий.
Минимуму второго интегрального критерия
на
графике (рисунок 3) соответствует точка, в которой
принимает
значение равное 0,95 от максимального в сторону увеличения частоты. Эта точка и
определит оптимальные параметры настройки ПИ- регулятора. Используя данные
таблицы 3 и рисунка 3, находим, что точке максимума соответствуют значения:
, Kp = 3.692
при ω = 0.12 с-1.
Поэтому оптимальные параметры настройки ПИ-
регулятора имеют значения:
, Kp∙
0,95=3.5074,
с.
. Расчёт, построение и оценка качества
переходных процессов
.1 Переходный процесс по каналу регулирующего
воздействия S-Y
Для одноконтурной системы регулирования,
приведенной на рисунке 1, определим передаточную функцию замкнутой АСР по
каналу S-Y по формуле:
(12)
где передаточная функция объекта регулирования:
,
передаточная функция ПИ- регулятора:
.
После подстановки значения
в
формулу (12), получаем окончательное выражение для передаточной функции
замкнутой АСР по каналу S-Y: