3066
.pdfМинистерство науки и высшего образования РФ Федеральное государственное бюджетное образовательное учреждение высшего образования
«Воронежский государственный лесотехнический университет имени Г.Ф. Морозова»
Теория автоматического управления
Методические указания к выполнению курсового проекта для студентов по направлению подготовки
15.03.04 Автоматизация технологических процессов и производств
Воронеж 2018
УДК 681.5
Грибанов, А. А. Теория автоматического управления [Электронный ресурс] : метод. указания к выполнению курсового проекта для студентов по направлению подготовки 15.03.04 – Автоматизация технологических процессов и производств / А. А. Грибанов; ВГЛТУ. – Воронеж, 2018. – 22 с.
Печатается по решению редакционно-издательского совета ВГЛТУ
Рецензент: д.т.н., профессор, зав. кафедрой электротехники и автоматики ФГБОУ ВО «Воронежский государственный аграрный университет имени императора Петра I» Афоничев Д.Н.
2
Введение
Целью курсового проекта является закрепление теоретических знаний в области анализа объектов управления, выбора средств управления, анализа и синтеза системы комбинированного автоматического управления. Каждый студент, в соответствии с вариантом задания, получает технологический объект автоматического управления с входным параметром по каналу управления, регулируемым параметром и возмущающим воздействием, а также передаточные функции объекта по каналу управления и каналу основного возмущающего воздействия. При выполнении курсового проекта студент проводит анализ и синтез одноконтурных, многоконтурных систем автоматического регулирования, используя пакет прикладных программ MathCAD, а также знания и навыки исследования звеньев и устойчивости автоматических систем регулирования, полученные на лабораторных и практических занятиях по курсу «Теория автоматического управления» (ТАУ).
При выполнении курсового проекта каждый его раздел иллюстрируется схемами, графиками в формате А4 и по каждому разделу студент дает необходимое заключение о показателях качества средств и систем регулирования.
Графики статических, динамических и частотных характеристик необходимо строить в математическом пакете прикладных программ MathCAD.
В конце курсового проекта нужно представить развернутое заключение на 1-2 страницы, которое должно содержать техническую характеристику предложенной системы комбинированного управления, ее технический уровень, показатели качества, надежности, практической значимости.
Курсовой проект заканчивается списком используемой литературы.
3
1 Порядок выполнения курсового проекта
1.1 Введение
Студенту с помощью литературного источника необходимо отметить актуальность предложенной ему комбинированной системы автоматического управления для повышения эффективности соответствующих технологий лесного комплекса. Необходимо отметить также цель курсового проекта и перечислить основные научно-теоретические и прикладные результаты курсового проектирования.
1.2 Описание работы технологического объекта автоматического управления
Студент для своего объекта управления в учебной и научно-технической литературе находит описание технологического процесса, работы оборудования и дает соответствующую технологическую схему. При этом представляет достаточно подробное описание и отмечает важность, актуальность разработки рассматриваемой системы управления. Желательно сделать критический анализ существующих систем автоматического регулирования по отклонению, в том числе в вопросах быстродействия, точности этих систем и их возможностей эффективной компенсации возмущающих воздействий.
1.3 Составление структурной схемы одноконтурной системы автоматического регулирования по отклонению
|
|
|
|
|
|
|
|
f(t) |
|
|
|
||
g(t) |
|
|
x(t) |
|
xp1(t) |
|
xp(t) |
|
y1(t) |
|
y(t) |
||
W(P)рег |
W(P)им |
W(P)об |
W(P)дат |
||||||||||
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 1. Структурная схема одноконтурной АСР
где g(t) – задающее воздействие задатчика регулятора; x(t) – ошибка регулирования;
xр1(t) – регулирующее воздействие регулятора на исполнительный меха-
низм;
xр(t) – регулирующее воздействие исполнительного механизма на объект регулирования, изменяющее подачу вещества или энергии на объект регулирования;
y1(t) – регулируемая величина на выходе объекта регулирования;
y(t) – измеряемая датчиком выходная величина объекта, которая сравнивается с задающим воздействием g(t);
x(t)=g(t)-y(t) – ошибка регулирования;
f(t) – возмущающее воздействие на объект регулирования.
4
Определяются передаточные функции W(P)рег, W(P)им, W(P)дат. Передаточная функция объекта регулирования W(P)об показана в задании Вашего варианта.
1.4 Выбор технических средств системы автоматического регулирования
Выбор датчиков, регуляторов, исполнительных механизмов должен осуществляться исходя из требований надежности, быстродействия и экономичности. При этом необходимо учитывать диапазон возможного изменения регулируемого параметра, согласованность технических характеристик и физических сигналов на входе и выходе датчика, регулятора, исполнительного механизма и объекта. Из каталогов средств автоматизации, справочной литературы, электронных ресурсов сети Интернет выбрать отечественные современные средства автоматизации, определить коэффициенты усиления и постоянные времени. Выбранные технические средства определяют конкретную по каждому варианту структуру одноконтурной системы регулирования. После этого необходимо также выбрать датчик измерения основного возмущающего воздействия.
В результате строится принципиальная электрическая схема АСР и дается описание работы автоматической системы с указанием технических характеристик элементов автоматики и объекта регулирования.
После выбора типов, конструкций регуляторов, исполнительных механизмов, датчиков необходимо проставить размерности величин g(t), x(t), xр1(t), xр(t), y1(t), y(t) (размерности параметров объекта xр(t) и y1(t) известны по Вашему варианту задания). Размерности:
k рег |
xр1 (t) |
, kим |
xр (t) |
|
|
|
y (t) |
|
|
y(t) |
|
||||
|
|
|
, |
kоб |
1 |
|
, kдат |
|
|
. |
|||||
x(t) |
x |
р1 (t) |
xp |
(t) |
y1 |
(t) |
|||||||||
|
|
|
|
|
|
|
|
Перемножить размерности коэффициентов передачи
|
x |
р1 (t) |
|
x р (t) |
|
y (t) |
|
y(t) |
|
|||||
kрег∙ kим∙ kоб∙ kдат= |
|
|
∙ |
|
|
|
∙ |
|
1 |
|
∙ |
|
. |
|
x(t) |
x |
р1 |
(t) |
x |
p |
(t) |
y (t) |
|||||||
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
1 |
|
После сокращения получим безразмерный коэффициент усиления системы
в разомкнутом состоянии k |
|
|
y(t) |
|
y(t) |
. |
раз |
|
|
||||
|
|
x(t) |
|
g(t) y(t) |
||
|
|
|
|
1.5 Исследование статических, переходных и частотных характеристик объектов управления
Передаточная функция объекта по каналу управления имеет следующий
вид
W |
(P) |
Y (P) |
, |
|
|||
об g |
|
X p (P) |
|
|
|
где Y(P) – изображение по Лапласу регулируемого параметра объекта;
Xp(P) – изображение по Лапласу входного параметра объекта по каналу управления.
Рассмотрим передаточные функции объектов 1-го и 2-го порядка
5
|
Wоб g (P) |
|
kоб |
|
|
Y (P) |
, |
|
|||
|
Tоб P 1 |
X p (P) |
|
||||||||
|
|
|
|
|
|
|
|||||
для 1, 8, 12 варианта будет 2-го порядка |
|
||||||||||
Wоб g |
(P) |
|
|
kоб |
|
|
|
|
Y (P) |
, |
|
(Tоб1P |
|
|
|
|
|
|
X p (P) |
||||
|
|
1) (Tоб2 P 1) |
|
|
где kоб – коэффициент усиления;
Тоб, Тоб1, Тоб2 – постоянные времени объекта по каналу управления.
Применив обратное преобразование Лапласа, получим дифференциальные уравнения
Tоб y (t) y(t) kоб x p (t) ,
Tоб1Tоб2 y (t) (Tоб1 Tоб2 ) y (t) y(t) kоб x p (t) ,
где y(t) – оригинал по Лапласу регулируемого параметра объекта;
xp(t) – оригинал по Лапласу входного параметра объекта по каналу управления. Для анализа переходного и установившегося процесса в объекте регулирова-
ния нужно показать сумму общего и частного решения дифференциальных уравнений первого и второго порядка
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y(t) k |
об |
x |
p |
(t) (1 e Tоб ) , |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
для 1, 8, 12 варианта будет 2-го порядка |
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
t |
|
|
|
|
Tоб |
|
|
|
|
|
|
|
Tоб |
|
|
|
||||
|
|
|
|
T |
|
|
|
|
2 |
T |
||||||||
y(t) kоб x p (t) |
1 |
|
|
|
1 |
|
e |
об1 |
|
|
|
|
e об 2 |
. |
||||
|
|
|
|
Tоб1 Tоб2 |
||||||||||||||
|
|
|
Tоб1 Tоб2 |
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Вматематическом пакете прикладных программ MathCAD провести расчет
ипостроить переходную характеристику y(t) объекта для единичного входного
воздействия по каналу управления xр(t) = 1(t). Для установившихся режимов построить статическую характеристику объекта у(хр). Сделать вывод об инерционных и статических свойствах объекта.
Вматематическом пакете прикладных программ MathCAD выполнить анализ свойств объекта в частотной области. Для этого в передаточных функциях объекта по каналу управления перейти от преобразования Лапласа к преобразованиям Фурье заменой P = jω, в результате получим частотную передаточную функцию объекта
Wоб g ( j ) Y ( j ) U ( ) jV ( ) ,
X p ( j )
где ω – круговая частота входных воздействий;
U(ω) – вещественная часть частотной передаточной функции; jV(ω) – мнимая часть частотной передаточной функции;
j 1 .
Амплитудно-частотная характеристика (АЧХ), логарифмическая ампли- тудно-частотная характеристика (ЛАЧХ) и фазово-частотная характеристика (ФЧХ) объекта имеют следующий вид
6
A( ) U 2 ( ) V 2 ( ) ,
Ai ( ) 20lg A( ) ;
( ) arctg V ( ) .
U ( )
Таким образом для объекта регулирования первого порядка получим следующие выражения частотных характеристик
АЧХ |
A( ) |
|
|
kоб |
|
|
; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|||||
|
|
T |
2 |
2 1 |
|
|
|
||||
|
|
|
об |
|
|
|
|
|
|
|
|
ЛАЧХ |
A ( ) 20lg k |
об |
20lg |
T 2 2 |
1 ; |
||||||
|
i |
|
|
|
|
об |
|
|
|||
ФЧХ |
( ) arctg (Tоб ) . |
|
|
|
Для объекта регулирования второго порядка частотные характеристики имеют следующий вид
АЧХ |
A( ) |
|
|
|
|
kоб |
|
|
|
|
; |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
T 2 |
2 |
1 T 2 2 1 |
|
|
|
|||||||||
|
|
|
об1 |
|
|
|
|
об2 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|||||||
ЛАЧХ |
A ( ) 20lg k |
об |
20lg |
|
T 2 2 |
1 20lg |
T 2 2 |
1 ; |
|||||||||
|
i |
|
|
|
|
|
об1 |
|
|
|
|
об2 |
|
|
|||
ФЧХ |
( ) arctg (Tоб1 ) arctg (Tоб2 ) . |
|
|
|
Графики частотных характеристик объекта строить в частотном диапазоне ω = 0…1000 Гц, с шагом 0,01 Гц. Сделать вывод о приемлемом диапазоне пропускания частот входных воздействий по каналу управления.
1.6 Исследование статических, переходных и частотных характеристик датчика регулируемого параметра
Для определения передаточной функции датчика необходимо в справочной литературе, каталогах и сети Интернет найти диапазон изменения входных и выходных сигналов, и на этой основе определить его коэффициент усиления. Также необходимо определить постоянную времени датчика, если его передаточная функция описывается апериодическим звеном первого порядка.
После того как определены передаточная функция датчика и дифференциальные уравнения, в математическом пакете прикладных программ MathCAD строятся графики статической и переходной характеристик, а также графики частотных характеристик: АЧХ, ЛАЧХ и ФЧХ. В конце пункта необходимо сделать вывод о приемлемом диапазоне пропускания частот, инерционности и быстродействии датчика регулируемого параметра.
1.7 Исследование статических, переходных и частотных характеристик датчика возмущающего воздействия
Определив передаточную функцию датчика, в математическом пакете прикладных программ MathCAD необходимо построить графики статической и переходной характеристик, а также графики частотных характеристик: АЧХ, ЛАЧХ и ФЧХ. В конце сделать вывод о приемлемом диапазоне пропускания частот,
7
инерционности, быстродействии датчика возмущающего параметра.
1.8 Исследование статических, переходных и частотных характеристик автоматического регулятора
Для создания высококачественной системы автоматического регулирования необходимо из каталогов средств автоматизации, справочной литературы, электронных ресурсов сети Интернет выбрать современный ПИД-регулятор. Для анализа начального состояния системы регулирования нужно по технической характеристике регулятора взять средние параметры настройки – коэффициент усиления (предел пропорциональности), время изодрома, время предварения. По этим начальным параметрам настройки регулятора определить дифференциальное уравнение и передаточную функцию регулятора. В математическом пакете прикладных программ MathCAD построить графики статической, переходной характеристики регулятора, а также графики частотных характеристик: АЧХ, ЛАЧХ и ФЧХ.
Сделать заключения о динамических и частотных свойствах регулятора при его средних настройках: коэффициента усиления, времени изодрома, времени предварения.
1.9 Исследование статических, переходных и частотных характеристик исполнительного механизма
Студентом из справочной литературы, электронных ресурсов сети Интернет выбирается исполнительный механизм, для которого необходимо определить коэффициент усиления и постоянную времени. Для этого следует учитывать размерность входных и выходных сигналов механизма. Чаще всего исполнительный механизм представляет собой интегрирующее инерционное звено.
По этим данным составляются передаточная функция и дифференциальное уравнение. В математическом пакете прикладных программ MathCAD строятся графики статических, переходных и частотных характеристик: АЧХ, ЛАЧХ и ФЧХ. По результатам этих расчетов делается вывод о быстродействии механизма, о работоспособном в диапазоне частот входных воздействий регулятора на исполнительный механизм.
1.10 Анализ работоспособности, устойчивости системы регулирования методом Найквиста-Михайлова
Для получения передаточной функции системы в разомкнутом состоянии Wраз(P) необходимо выполнить произведение передаточных функций объекта по каналу управления Wобg(P), датчика регулируемого параметра Wдат(P), регулятора Wрег(P) и исполнительного механизма Wим(P)
Wраз (P) Wоб g (P) Wдат (P) Wрег (P) Wим (P) .
Осуществив в этой передаточной функции переход от преобразования Лапласа к преобразованию Фурье, получим частотную передаточную функцию системы в разомкнутом состоянии
Wраз ( j ) Wоб g ( j ) Wдат ( j ) Wрег ( j ) Wим ( j ) .
8
Используя инструменты для преобразования мнимых чисел математического пакета прикладных программ MathCAD, получим
Wраз( j ) U ( ) jV ( ) .
Впакете прикладных программ MathCAD на комплексной плоскости нужно построить амплитудно-фазовую частотную характеристику (АФЧХ) системы
вразомкнутом состоянии в частотном диапазоне ω = 0…1000 Гц, с шагом 0,01 Гц. На этом графике, проведя окружность единого радиуса, определяется запас устойчивости системы по амплитуде а и по фазе γ.
Если система регулирования окажется неустойчивой (ωс > ωπ), то необходимо изменить параметры настройки автоматического регулятора (коэффициент усиления, время изодрома, время предварения) в направлении вывода системы в устойчивое состояние и привести новые расчеты запасов устойчивости по фазе и амплитуде. Вывод должен отражать состояние системы с позиций устойчивости при окончательно принятых параметрах настройки регулятора.
1.11Анализ показателей качества системы
в замкнутом состоянии
Целью анализа является определение быстродействия (времени переходного процесса), динамического заброса (перерегулирования), колебательности, статической ошибки регулирования системы в замкнутом состоянии.
Используя полученную передаточную функцию системы с новыми параметрами настройки регулятора в разомкнутом состоянии по каналу управления Wраз(P), определим передаточную функцию системы в замкнутом состоянии
Wзам(P) |
|
Wраз(P) |
|
Y (P) |
, |
|
Wраз(P) |
|
|||
1 |
|
G(P) |
|
где Y(P) – изображение по Лапласу регулируемого параметра;
G(P) – изображение по Лапласу управляющего (задающего) воздействия.
Из этого выражения (при P d / dt ) получается дифференциальное уравнение системы в замкнутом состоянии
a |
|
d n y |
a |
d n 1 y |
... a |
|
dy |
a |
|
y b |
d m g |
b |
d m 1g |
... b |
dg |
b g , |
0 dt n |
|
n 1 dt |
|
|
|
|
||||||||||
|
1 dt n 1 |
|
|
n |
0 dt m |
1 dt m 1 |
m 1 dt |
m |
где a0, a1, …, an-1, an, b0, b1, …, bm-1, bm – коэффициенты дифференциального уравнения;
y – оригинал по Лапласу регулируемого параметра; g – оригинал по Лапласу управляющего воздействия.
В математическом пакете прикладных программ MathCAD найти корни характеристического уравнения, приравняв нулю правую часть и определить показатели качества системы. Построить график изменения регулируемого параметра y(t) при скачкообразном единичном изменении управляющего воздействия g(t) = 1(t). Проанализировать показатели качества, определить необходимость некоторого повышения быстродействия и снижения динамического заброса.
Корни характеристического уравнения
Pi i j i i=1, 2, …, n,
9
отобразить на комплексной плоскости в пакете прикладных программ MathCAD (α, jβ), определить величину степени устойчивости
min |
|
i |
|
||
i |
|
|
и колебательность системы
I max |
i |
, |
i |
i |
|
где αi и βi – положительные коэффициенты. |
|
|
1.12 Анализ устойчивости системы регулирования методом годографа Михайлова
Устойчивость системы в замкнутом состоянии определяется методом годографа Михайлова, что позволяет определить близость системы к колебательной и апериодической границе устойчивости. Для этого следует воспользоваться характеристическим уравнением систем в замкнутом состоянии
D(P) 1 Wраз(P) a0Pn a1Pn 1 ... an 1P an 0 .
Из состояния системы во временной области нужно перейти в частотную область, заменой P = jω.
После этого преобразования характеристический полином будет иметь вид
D( ) Rе( ) jIт( ) ,
где действительная часть
Rе( ) an an 2 2 ...,
мнимая часть
Iт( ) an 1 an 3 3 ....
Произвести разделение характеристического уравнения на действительную и мнимую часть в математическом пакете прикладных программ MathCAD и построить на комплексной плоскости (Re(ω); jIm(ω)) при изменении частоты от 0 до ∞ годограф Михайлова и определить насколько близка система при принятых настройках регулятора к апериодической и колебательной границе устойчивости.
1.13 Определение необходимых настроек регулятора при синтезе систем комбинированного управления
Целью синтеза является определение необходимых настроек регулятора, параметров корректирующего устройства охвата обратной связью инерционных объектов регулирования и исполнительных механизмов, а также параметров блока коррекции, который делает систему инвариантной к самому сильному заданному возмущающему воздействию.
Необходимо внимательно проанализировать график амплитудно-фазовой частотной характеристики системы в разомкнутом состоянии (при анализе системы на устойчивость методом Найквиста-Михайлова), годограф Михайлова и принять решение об изменении коэффициента усиления регулятора kpег, величина которого входит в значение коэффициента аn характеристического уравнения.
10