УДК 681.51
А. В. Епифанцев1, В. П. Кривошеев2, М. А. Сачко3 ИССЛЕДОВАНИЕ КАСКАДНЫХ СИСТЕМ УПРАВЛЕНИЯ
Выводятся передаточные функции двухконтурной каскадной АСР для основной и вспомогательной регулируемых переменных. Рассматривается метод вещественных частотных характеристик для построения переходного процесса. Проводится сравнение переходных процессов одноконтурной и каскадной АСР. Проводится исследование зависимости регулируемых переменных от изменения времени запаздывания по каналам возмущения.
Ключевые слова: автоматическая система регулирования, каскадная система регулирования, передаточная функция, переходный процесс, формула Мейсона, метод вещественных частотных характеристик.
Введение
При автоматизации сложных технологических объектов управления (ТОУ) часто возникает задача применения многоконтурных автоматических систем регулирования (АСР), одними из представителей которых являются каскадные АСР. Такие АСР обеспечивают достаточно высокое качество регулирования ТОУ, имеющих помимо выходных регулируемых координат еще и дополнительные промежуточные координаты, которые могут быть использованы для организации вспомогательных контуров стабилизации. Вследствие этого каскадные АСР всегда имеют многоконтурную структуру и включают в свой состав два и более регулятора. При автоматизации химико-технологических объектов чаще всего используют двухконтурные каскадные системы [1, 2].
1 © Алексей Владимирович Епифанцев, аспирант кафедры информационных систем и компьютерных технологий Владивостокского государственного университета экономики и сервиса, ул. Гоголя, 41, г Владивосток, Приморский край, 690014, Россия, E-mail: aale-xey @rambler. ru.
2 © Владимир Петрович Кривошеев, д-р техн. наук, профессор кафедры информационных систем и компьютерных технологий Владивостокского государственного университета экономики и сервиса, ул. Гоголя, 41, г. Владивосток, Приморский край, 690014, Россия, E-mail: krivosheev@vvsu.ru.
3 © Максим Анатольевич Сачко, старший преподаватель кафедры информационных систем и компьютерных технологий Владивостокского государственного университета экономики и сервиса, ул. Гоголя, 41, г. Владивосток, Приморский край, 690014, Россия, E-mail: keeper@vvsu.ru.
На рисунке 1 изображена схема двухконтурной каскадной АСР в общем виде. На практике, как правило, рассматривается каскадная АСР с подачей возмущающего воздействия на вход ТОУ по каналу управления (базовый случай), схема которой изображена на рис. 2. Каскадная АСР общего вида может быть приведена к базовому виду, если взять передаточные функции по каналам возмущения равными:
Wв (^) = W0(s), (1)
ад=w0l(s) (2)
Рис. 1. Схема каскадной АСР в общем виде
Рис. 2. Схема каскадной АСР при подаче возмущающего воздействия на вход
ТОУ по каналу управления
В данной работе решаются следующие задачи:
- получение передаточной функции замкнутой каскадной АСР по каналу возмущения для расчета основной регулируемой переменной;
- получение передаточной функции замкнутой каскадной АСР по каналу возмущения для расчета вспомогательной регулируемой переменной;
- построение переходного процесса каскадной АСР с помощью метода вещественно-частотных характеристик;
- сравнение переходных процессов каскадной и одноконтурной АСР;
- исследование зависимости основной и вспомогательной регулируемых переменных двухконтурной каскадной АСР от времени запазды-
К (р)
вания каналов возмущения и инерционности канала вЧ ' при приведении каскадной АСР общего вида к базовому случаю.
Получение передаточных функций замкнутой каскадной АСР
Для получения передаточной функции замкнутой каскадной АСР по каналу возмущения используется формула Мейсона [3]:
1 т
ф- Е К,и,(4')д, (*) (3)
(,(р) ,
где ек (г) - передаточная функция 1-го прямого канала, связывающего
X у Д( р)
вход к с выходом Уе; т - число таких каналов; v ’ - специальный полином, который определенным образом характеризует совокупность всех замкнутых цепей системы, содержащих обратные связи, и вычисляется как сумма передаточных функций разомкнутых контуров этих цепей и произведений передаточных функций разомкнутых контуров пар, троек и т.д. не соприкасающихся друг с другом цепей с обратными связями:
Д(р)= 1 - Ек, (р) + (р)Ж}. (р)-Ек, (р)Ж}. (р)жг (р)+... (4)
г, ] г,3,г
Полином '(составляется по правилу, аналогичному (4), но только для цепей с обратными связями, не соприкасающихся с і-м прямым каналом. Знаки всех сигналов прямых каналов и обратных связей учитываются в формулах (3) и (4) перед соответствующими передаточными функциями.
Для получения передаточной функции для расчета основной регулируемой переменной был построен соответствующий ей сигнальный граф (рис. 3), на котором отображены основные узлы системы и переходы между ними.
Рис. 3. Сигнальный граф каскадной АСР
На построенном графе были выделены графы двух прямых каналов (рис. 4) и двух замкнутых цепей (рис. 5).
Поскольку все контуры соприкасаются (имеют общие ветви или вершины), то парные произведения и триады в формуле отсутствуют. Нет также в сигнальном графе контуров, не соприкасающихся с прямыми каналами, поэтому ^і (^)_ 1 и ^1 )_ 1.
Рис. 4. Сигнальные графы прямых каналов
Рис. 5. Сигнальные графы замкнутых цепей Таким образом:
Ж! (5) = (*), (5)
Ж (5) = -Жв1(5)ВДЖ0 (5), (6)
Д(5)=1+адж»+вдад, (7)
АI(5) = ДII(5) =1. (8)
Подставив (5), (6), (7) и (8) в (3), получим передаточную функцию замкнутой каскадной АСР:
ж(5) = Ж (5) - Ув1(5)Л;(5)У0 (5)
^ 1 + Л(5) Д (5)Ж0 (5) + Д (5)Жо1 (5) (9)
Аналогичным образом была получена передаточная функция для расчета вспомогательной регулируемой переменной. Был построен сигнальный граф (рис. 6) и выделены графы прямых каналов (рис. 7). Графы замкнутых контуров для этого случая те же, что и при расчете передаточной функции для основной переменной. В результате была получена передаточная функция каскадной АСР для расчета вспомогательной регулируемой переменной:
Ж (5):
Ж» - Ж (д)ОДадУо1(д)
1+ВДОДЖ (5)+ВДЖ»
Рис. 6. Сигнальный граф каскадной АСР
(10)
Рис. 7. Сигнальные графы прямых каналов
Построение переходного процесса
Эффективность работы АСР, как правило, определяют с помощью графика переходного процесса, показывающего изменение регулируемого параметра системы при ступенчатой подаче возмущающего воздействия. В данной работе для получения переходного процесса использовался метод вещественных частотных характеристик [4]:
а
-сію»-
(11)
Для построения переходного процесса от передаточных функций (9) и (10) переходят к амплитудно-фазовой характеристике (АФХ), заменяя
переменную s на jw, где w - круговая частота, а j - мнимая единица. В
(11) Re(w) есть вещественная частотная характеристика АФХ замкнутой каскадной АСР.
На практике в формуле (11) переменная 1 изменяется от нудя до
wmax / Dw wmax
max , где max - максимальная частота системы, а - шаг по частоте.
Для определения максимальной частоты системы строится график
вещественной частотной характеристики Re(w). За Wmax берется частота
затухания Re(w).
С целью автоматизации расчетов была разработана программа «Transient Cascade», предназначенная для построения графиков вещественночастотных характеристик и переходных процессов для основной и вспомогательной регулируемых переменных двухконтурных каскадных АСР. Алгоритмы расчета графиков вещественных частотных характеристик и переходных процессов каскадных АСР будут впоследствии интегрированы в учебно-методический комплекс «Studium» [5].
Сравнение переходных процессов каскадной и одноконтурной АСР
С использованием программы «Transient Cascade» был построен переходный процесс каскадной АСР для сравнения с переходным процессом одноконтурной АСР.
В качестве примера рассмотрена каскадная АСР со следующими ха-
1 < 1т
W (s) = 10 W (s) = 10——е -2 s
рактеристиками: 1 - 100s , 1 - 10s ,
11 W (s) = 4s (s) = 5-— e“1^
1 - 80s , 1 - 5s .
Для данной АСР были рассчитаны параметры регуляторов:
4 1,962 - 0,34s л . 0,289 - 0,676s
R(s) =------------ ^(У) =---------------
s , s .
Для того же ТОУ были рассчитаны параметры регулятора однокон-
Wo (s) = 10---1--e"5s W (s) = 10—1— e~ 4s
турной АСР: 0 1 - 100s , в 1 - 80s ,
4 0,131 - 1,67s
R(s) =------------
s.
После этого были построены переходные процессы для одноконтурной (рис. 8) и каскадной (рис. 9) АСР.
Рис. 8. Переходный процесс одноконтурной АСР
Рис. 9. Переходный процесс каскадной АСР Из рисунков 8 и 9 следует, что в данном примере время затухания графиков переходных процессов и максимальная амплитуда переходного процесса у каскадной АСР значительно ниже, чем у одноконтурной, что свидетельствует о преимуществе каскадной АСР перед одноконтурной.
Исследование зависимости регулируемых переменных двухконтурной каскадной АСР от времени запаздывания и инерционности каналов возмущения при приведении каскадной АСР общего вида к базовому случаю
Для исследования зависимости регулируемых переменных двухконтурной каскадной АСР от времени запаздывания каналов возмущения была взята АСР со следующими передаточными функциями:
К (*) = К (*) = 10^^^ е -5* к (*) = к» = 10- 1
-2s
1 +1005
1 +105
ч 1,962 + 0,345 . 0,289 + 0,6765
ад =---------------- ^(5) =----------------
5 5
Сначала изменялось время запаздывания передаточной функции Ж (5)
в14 , как в большую, так и в меньшую сторону. На рисунках 10, 11 и 12
изображены графики переходных процессов основной регулируемой пе-
“ У Жв1(5) .. 0 ,
ременной * при времени запаздывания канала «іч/, равном 1, 2 и 3 соответственно. На рисунках 13, 14 и 15 приводятся графики переходных
процессов вспомогательной регулируемой переменной у при времени
КМ 1 о о
запаздывания канала 61 , равном 1, 2 и 3 соответственно.
Рис. 10. Переходный процесс переменной У при времени запаздывания канала Кв1(^),
равном 1
Рис. 11. Переходный процесс переменной У при времени запаздывания канала Кв1(^),
равном 2
Рис. 12. Переходный процесс переменной У при времени запаздывания канала Ка1(^) ,
равном 3
Как видно из рис. 10, 11, 12, изменение времени запаздывания канала
КМ ф
814 ' не влияет на время затухания графика переходного процесса основной регулируемой переменной у и незначительно влияет на максимальную амплитуду графика. Однако при увеличении времени запаздывания существенно увеличивается значение переменной у в первой точке максимума.
Рис. 13. Переходный процесс переменной У1 при времени запаздывания канала К«1(я),
равном 1
V К (я)
Рис. 14. Переходный процесс переменной •'1 при времени запаздывания канала 81 ,
равном 2
Рис. 15. Переходный процесс переменной У1 при времени запаздывания канала К«1(я),
равном 3
На рисунках 13, 14, 15 видно, что при увеличении времени запазды-
Кв1(я) ф
вания канала в1 время затухания графика переходного процесса
вспомогательной регулируемой величины у также не изменяется. Однако заметно уменьшается значение переменной у в первой точке мини-
мума и второй точке максимума графика, таким образом, изменяется максимальная амплитуда графика переходного процесса.
После этого было построено несколько переходных процессов при
й ф Ж (¿)
изменении времени запаздывания передаточной функции канала в 4 '. На рисунках 16, 17 и 18 приведены графики переходных процессов основной регулируемой переменной У при времени запаздывания канала Ж (5)
вЧ , равном 3, 5 и 7 соответственно. На рисунках 19, 20 и 21 показаны графики переходных процессов вспомогательной регулируемой пе-
“ У Ж С?) -зет
ременной *1 при времени запаздывания канала вЧ , равном 3, 5 и 7
соответственно.
Рис. 16. Переходный процесс переменной У при времени запаздывания канала Ж (5) ,
равном 3
Рис. 17. Переходный процесс переменной У при времени запаздывания канала Ж (5) ,
равном 5
Рис. 18. Переходный процесс переменной У при времени запаздывания канала Ж (я),
равном 7
Из рисунков 16, 17, 18 видно, что при увеличении времени запазды-
Ж (я) ф
вания канала вЧ 7 время затухания графика не изменяется, но значительно уменьшается значение переменной У в первой точке максимума. Также при изменении времени запаздывания в большую или меньшую сторону относительно базового случая уменьшается значение величины У в первой и второй точках минимума и повышается во второй точке максимума, увеличивая тем самым максимальную амплитуду.
Рис. 19. Переходный процесс переменной У1 при времени запаздывания канала Ж (я).
равном 3
V Ж(я)
Рис. 20. Переходный процесс переменной 1 при времени запаздывания канала в
равном 5
Рис. 21. Переходный процесс переменной У1 при времени запаздывания канала Ж (я) ,
равном 7
Из рисунков 19, 20, 21 видно, что при увеличении времени запазды-
Ж (я) й У, й
вания канала вЧ 7 уменьшается значение переменной •У1 в первой точке
минимума. При изменении времени запаздывания в большую или меньшую сторону относительно базового случая увеличивается значение переменной У1 во второй и третьей точках максимума и уменьшается во второй точке минимума, следовательно, увеличивается максимальная амплитуда графика. На время затухания переходного процесса переменной
У Ж (я)
•У1 изменение времени запаздывания канала вЧ ' не влияет.
Затем были построены переходные процессы основной регулируемой
“ У “ АГП “ Ж (я)
переменной * каскадной АСР при разной инерционности канала в 7.
На рисунках 22, 23, 24 представлены графики переходных процессов пе-
ременной У при значениях переменной Т передаточной функции равных 80, 100 и 120 соответственно.
Ж (я)
Рис. 22. Переходный процесс переменной У при значении переменной Т канала Ж (я) .
равном 80
Рис. 23. Переходный процесс переменной У при значении переменной Т канала Ж (я) .
равном 100
Рис. 23. Переходный процесс переменной У при значении переменной Т канала Ж (я) ,
равном 120
Из рисунков 21, 22 и 23 видно, что при увеличении времени запаздыва-
Ж (я) й У й
ния канала в 4 7 значения переменной у в первой точке максимума и пер-
вой и второй точках минимума уменьшается, во второй и третьей точках максимума растет, таким образом, увеличивается максимальная амплитуда переходного процесса. При этом время затухания графика не изменяется.
Можно сделать вывод, что при увеличении времени запаздывания
й ф Жв1(я) й й
передаточной функции в1У ' двухконтурной каскадной системы регули-
рования незначительно возрастает максимальная амплитуда переходного процесса основной регулируемой переменной У и заметно увеличивается максимальная амплитуда переходного процесса вспомогательной переменной У1 . При этом время затухания переходных процессов обеих переменных не изменяется. При изменении времени запаздывания передаточной функции Ж (я) в большую или меньшую сторону относительно общего случая максимальная амплитуда переходных процессов регулируемых переменных У и У1 существенно возрастает, время затухания также
Ж (я)
не изменяется. При увеличении инерционности канала вЧ ’ значительно увеличивается максимальная амплитуда переходного процесса основной регулируемой переменной У , на время затухания изменение инерционно-Ж (я)
сти канала в ч / не влияет.
1. Стефани Е.П. Основы настройки регуляторов теплоэнергетических процессов / Е.П. Стефани. - М.: Энергия, 1972.
2. Автоматическое управление в химической промышленности / Е.Г. Дудников, А.В. Казаков, Ю.Н. Софиева, А.Э. Софиев, А.М. Цир-лин. - М.: Химия, 1987.
3. Лукас В. А. Теория управления техническими системами /
В. А. Лукас. - Екатеринбург, 2002.
4. Теория автоматического управления / С.Е. Душин, Н.С. Зотов, Д.Х. Имаев, Н.Н. Кузмин, В.Б. Яковлев. - М.: Высшая школа, 2003.
5. Епифанцев А.В., Сачко М. А., Кривошеев В.П. Разработка информационно-программного обеспечения для исследования комбинированных систем управления // Интеллектуальный потенциал вузов - на развитие Дальневосточного региона России и стран АТР: материалы XI международной конференции студентов, аспирантов и молодых ученых. -Владивосток: Изд-во ВГУЭС, 2009. - С. 75 - 79.