CC BY
52
ИНФОРМАТИКА И ТЕХНИЧЕСКИЕ НАУКИ
МОДЕЛИРОВАНИЕ ПРОЦЕССОВ И УПРАВЛЕНИЕ В ТЕХНИЧЕСКИХ СИСТЕМАХ
УДК 681.53
Д. А. Гринюк, И. Г. Сухорукова, Н. М. Олиферович, И. О. Оробей
Белорусский государственный технологический университет
АНАЛИЗ ЭФФЕКТИВНОСТИ ПИД-РЕГУЛЯТОРОВ С ДВУМЯ СТЕПЕНЯМИ СВОБОДЫ С ПОМОЩЬЮ ИНТЕГРАЛЬНЫХ КРИТЕРИЕВ
В статье рассматриваются различные варианты настройки ПИД-алгоритма для одновременного решения задачи компенсации отклонений по сигналам задания и возмущения. Исследованы различные варианты конструкции регуляторов с двумя степенями свободы (of the two-degree-of-freedom structure, 2DoF). Для объекта управления с несколькими вариантами динамики произведен анализ как известной структуры регулятора с двумя степенями свободы, так и собственный вариант. Также были разработаны и рассмотрены варианты построения нелинейного ПИД-регулятора с двумя степенями свободы. В качестве основы настройки и последующего сравнения линейных и нелинейных вариантов 2DoF-регуляторов использованы два варианта интегральных критериев. Показано, что два варианта структур с двумя степенями свободы улучшают работу ПИД-регулятора, но не являются эквивалентными и не обеспечивают одновременной минимизации интегральных критериев возмущения и заданного значения. Сделан вывод, что при выборе структур для решения прикладных задач, следует обладать априорной информацией об амплитуде и интенсивности основных возмущений, которые будут воздействовать на объект управления. Среди использованных для комплексной настройки интегральных критериев модифицированный интегральный критерий приводит к более приемлемому для технологических объектов переходному процессу.
Обнаружено меньшее расхождение между работами двух нелинейных структур 2DoF. Для рассмотренных вариантов динамики объекта управления с точки зрения интегральных критериев нелинейная структура лучше справлялась с подавлением возможных отклонений.
Ключевые слова: ПИД-регулятор, регулятор с двумя степенями свободы, нелинейный регулятор, интегральные критерии настройки.
D. A. Hryniuk, I. G. Suhorukova, N. M. Oliferovich, I. O. Orobei
Belarusian State Technological University
ANALYSIS OF THE EFFICIENCY OF PID CONTROLLER WITH TWO-DEGREE-OF-FREEDOM STRUCTURE WITH INTEGRATED CRITERIA
The article considers various options for configuring the PID algorithm for simultaneous solution of the problem of compensation of deviations from the reference and disturbance signals. Different variants of the controller with two-degree-of-freedom structure (2DoF) are investigated. For a control object with several variants of dynamics, analysis of both the known structure of a controller with two degrees of freedom, so its own version, is performed. Also options for constructing a nonlinear PID controller with two degrees of freedom were developed and considered. Two variants of integral criteria are used as a basis for tuning and subsequent comparison of linear and nonlinear variants of 2DoF controllers. It is shown that two versions of structures with two degrees of freedom improve the operation of the PID controller, but are not equivalent and don't provide simultaneous minimization of integral disturbance criteria and a given value. It is concluded that when choosing the structures for solving applied problems, it is necessary to have a priori information about the amplitude and intensity of the main disturbances that will affect the control object. Among the integral criteria used for complex turning, the modified integral criterion leads to a more transient process that is more acceptable for technological processes.
A smaller discrepancy is found between the work of two nonlinear 2DoF structures. For the considered variants of the dynamics of the control object, from the point of view of integral criteria, the nonlinear structure copes better with suppression of possible deviations.
Key words: PID controller, controller with two-degree-of-freedom structure, non-linear PID, integral performance criteria.
Введение. Качественная настройка промышленных регуляторов способствует повышению стабильности технологических параметров и уменьшает энергозатраты на единицу продукции. Это обусловлено уменьшением переходных процессов и критических отклонений. Практика обучения специалистов теории управления, а также большинство учебной литературы используют наблюдение за переходным процессом вследствие изменения сигнала задания для проверки качества регулирования. В то же время большинство систем стабилизации технологических параметров предназначены для подавления возмущений в процессе функционирования. Иногда возникают задачи, при которых система подвергается частым возмущениям как со стороны канала возмущения, так и ввиду частого изменения сигнала задания.
Основная часть. Одним из вариантов решения проблемы является использование two-degree-of-freedom structure (2DoF) [1-7]. Существует классический вариант реализации структуры (рис. 1 [3]), но можно предложить и альтернативный вариант. Несмотря на введение дополнительных коэффициентов a и b, задачу поиска настроек регулятора можно легко разделить на два этапа: настройка основных коэффициентов KP, TI и TD, а затем поиск дополнительных коэффициентов a и b.
SP
iF®-
Td s
F
WAhS^
Wo
Ti S
Ws
Рис. 1. Классический вариант регулятора с двумя степенями свободы (КРДСС). БР - сигнал задания; Уоит - выход; КР - коэффициент усиления; Т1 - время интегрирования; Тв - время дифференцирования; - передаточная функция исполнительного механизма; - передаточная функция измерительного преобразователя; W0 - передаточная функция объекта; ¥ - сигнал возмущения; 5 - оператор передаточной функции
SP
Td S
Ws
Рис. 2. Альтернативный вариант регулятора с двумя степенями свободы (АРДСС)
В качестве критерия настройки, как отражено в работе [8], легко можно использовать интегральные критерии [9]. В работе были использованы два интегральных критерия настройки.
Первый критерий является классическим, в то же время второй обеспечивает компромиссный вариант между быстродействием, перерегулированием и запасом по устойчивости [9].
tf
J e(t)2 dt ^
min;
J t2 |e(t )dt
^ min.
(1)
(2)
При настройке по информационному каналу задания ошибка вычислялась как
е(() = БР - Уоит. (3)
При настройке по информационному каналу возмущения ошибка соответствовала выходному значению
e(t ) = Y
out■
(4)
Алгоритм настройки предполагает попеременное применение формул (3) и (4). Для структуры на рис. 1 сначала находятся настройки КР, Т1 и Тп по критерию с использованием формулы (4), а затем по тому же критерию уже с использованием выражения (3) рассчитываются дополнительные коэффициенты а и Ь. В случае использования структуры на рис. 2 КР, Т и Т0 определялись в соответствии с (3), а уже а и Ь по формуле (4).
а
b
out
1
Для исследования была выбрана система со следующими передаточными функциями:
ЖА () = () = •
1
55 +1
1 ; 105 + 1:
(5)
(6)
Жа (5 ) =
(155 + 1)255 + 1)
ех
р(-Т5). (7)
Поскольку транспортное запаздывание времени т существенно влияет на результат применения интегральных критериев [8], методика настройки была проверена для нескольких значений т = [0; 3; 7,5; 15].
Минимизация критериев производилась в Ма1;1аЬ. Время моделирования для каждого значения т было свое, оно определялось предварительным экспериментом. При использовании критерия (1) время (/ было в 1,5 раза больше, чем при использовании критерия (2), что являлось следствием большей длительности процесса затухания.
Кроме линейных структур использовался и ПИД-регулятор с нелинейной характеристикой [8, 10].
ц(е) = sign(e) 1п(1 + \в\)-
1
(8)
При этом нелинейное преобразование использовалось только для интегральной и пропорциональной составляющих регулятора (второй вариант нелинейного ПИД-регулятора на рис. 3).
БР
Рис. 3. Нелинейный альтернативный вариант регулятора с двумя степенями свободы
Соотношения между настройками линейного (ЬЯ) и нелинейного (КЯ) ПИД-регу-ляторов по каналу возмущения / и задания БР можно видеть в табл. 1-2, где КР,БР, Т1бБР, Т0,БР -настройки, которые были найденные при минимизации интегральных критериев по каналу сигнала задания; КР,р, Т1рр, Т0,р - настройки,
которые были найдены при минимизации интегральных критериев по каналу возмущения.
Таблица 1
Соотношения между параметрами регуляторов, которые найдены по критерию (1)
Тип ПИД Соотношение параметров т, с
0 3 7,5 15
ЬЯ КР,БР/КР,Р 0,8615 0,8686 0,8690 0,8847
Т1,БР/Т1,Р 16,4030 7,2977 4,4148 2,9489
ТДЯР/ТДР 0,5319 0,6145 0,6391 0,6593
КР,БР/КР,Р 0,9375 0,9188 0,9000 0,8647
Т1,БР/Т1,Р 11,9070 6,1048 3,8248 2,6693
ТДЯР/ТДР 0,5737 0,6377 0,6598 0,6812
Таблица 2 Соотношения между параметрами регуляторов, которые найдены по критерию (2)
Тип ПИД Соотношение параметров т, с
0 3 7,5 15
ЬЯ КР,БР/КР,Р 0,3689 0,5342 0,4882 0,6394
Т1,БР/Т1,Р 5,9179 3,4831 2,7151 1,9518
ТД5Р/ТДР 0,3486 0,5221 0,4102 0,5660
КЯ КР,БР/КР,Р 0,5797 0,7028 0,7083 0,7849
Т1,БР/Т1,Р 4,8036 2,8595 2,3341 1,6585
ТД5Р/ТДР 0,5107 0,6459 0,6205 0,7453
Табл. 3-7 получены путем деления значения интегрального критерия на минимальное значение интегрального критерия при данном запаздывании, где КЮБР - относительное значение интегрального критерия для настроек классического ПИД-регулятора, найденных при настройках по каналу задания; КЮР - относительное значение интегрального критерия для настроек классического ПИД-регулятора, найденных при настройках по каналу возмущения.
Таблица 3
Относительный интегральный критерий по возмущению по критерию (1)
Тип Соотношение т, с
ПИД параметров 0 3 7,5 15
АРДСС 3,9087 2,4091 1,8054 1,4594
ЬЯ КЮбр 5,2363 2,8780 2,0339 1,5660
КЮР 1 1,0114 1,0192 1,0246
АРДСС 3,1277 1,9829 1,5357 1,2920
КЯ КЮ8р 3,8849 2,3120 1,7251 1,4231
КЮр 1,0024 1 1 1
1
Таблица 4
Относительный интегральный критерий по заданию по критерию (1)
Тип ПИД Соотношение параметров т, с
0 3 7,5 15
ья КРДСС 2,4908 1,7292 1,4404 1,2781
КЮцр 1 1 1 1
КЮР 2,9579 2,1497 1,7578 1,4768
ж КРДСС 2,4406 1,8364 1,4698 1,2759
КЮхр 1,0576 1,0601 1,0504 1,0503
КЮр 2,4936 1,9342 1,5796 1,3642
Таблица 5
Относительный интегральный критерий
по возмущению по критерию (2)
Тип ПИД Соотношение параметров т, с
0 3 7,5 15
ья АРДСС 10,986 5,9113 4,4847 2,7973
КЮхр 28,975 14,531 7,0241 4,1731
К1БР 1,0260 1,0375 1,0036 1,0563
ж АРДСС 7,7601 4,8329 3,2140 2,3835
КЮ8р 27,575 12,215 5,8081 3,2652
К1БР 1 1 1 1
Таблица 6
Относительный интегральный критерий по заданию по критерию (2)
Тип Критерий т, с
ПИД 0 3 7,5 15
КРДСС 3,2499 2,6725 2,2413 2,0020
ья КЮхр 1,0964 1,0028 1,0798 1
К!БР 5,5760 4,6725 3,8762 3,2010
КРДСС 6,5172 4,1317 3,1874 1,5022
NN. КЮхр 1 1 1 1,0195
К!БР 6,5461 4,7148 4,0921 2,8739
отличаться в несколько раз. Наибольшее отличие наблюдается во времени интегрирования. Следует отметить, что использование нелинейного преобразования ошибки придает ПИД-регулятору большую универсальность при стабилизации технологических параметров.
Данные таблиц показывают, что дополнительные коэффициенты настроек позволяют уменьшить значение интегральных критериев, однако не всегда существенно. Наибольший эффект достигается для объектов без запаздывания или при небольшом его значении. Улучшения наблюдаются также и по перерегулированию, и по времени переходного процесса (рис. 4-9). Использование критерия (1), как обычно [9], приводит к высокой колебательности и длительности затухания.
Рис. 4. Переходной процесс по возмущению для запаздывания т = 0. Настройка по критерию (2)
Таблица 7
Относительное улучшение от использования структур КРДСС и АРДСС для обеспечения универсальности
Тип ПИД Критерий т, с
0 3 7,5 15
КЮр/КРДСС
ья (2) 1,7157 1,7483 1,7294 1,5989
(1) 1,1874 1,2430 1,2202 1,1554
ж (2) 1,0044 1,1411 1,2839 1,9131
(1) 1,0217 1,0533 1,0747 1,0692
КЮ8р!АРДСС
ья (2) 2,5328 2,3851 1,5598 1,4897
(1) 1,3565 1,1967 1,1271 1,0731
ж (2) 3,2352 2,4157 1,7805 1,3669
(1) 1,2483 1,1682 1,1252 1,1027
Анализ таблиц показывает, что отличия между оптимальными настройками регуляторов по каналу задания и возмущения могут
0,7 0,6 0,5 £ 0,4 0,3 0,2 0,1 0 -0,1
/^Л \ АРД (СС - ья
£Гч V \ #л.4 V ' \ 4 V - ж
\ 4 \ \ V V
КЮр V \\
ч \\ \ \\ N \\
"Ктп;р
V и ЧЧ
0
50
100
150
и С
200
250
300
350
Рис. 5. Переходной процесс по возмущению для запаздывания т = 15. Настройка по критерию (2)
В целом структура КРДСС за счет дополнительных коэффициентов способствует уменьшению перерегулирования и времени переходного процесса. Структура АРДСС, в первую очередь, способствует улучшению быстродействия и уменьшению переходного процесса.
Результаты работы нелинейного регулятора по анализу таблиц и переходных процессов свидетельствуют о его преимуществе в сравнении с линейным регулятором. С практической точки зрения он показывает выигрыш не только по интегральному критерию, но и по времени переходного процесса, перерегулированию.
Использование интегрального критерия (1) не способствует получению качественного переходного процесса для работы в прикладных условиях. Как было показано в публикации [9], использование данного критерия приводит к малому запасу по устойчивости, большому перерегулированию и длительному процессу затухания колебаний.
1,5
0,5
kidf я f 's* КРДСС
' 'я > ^ KIDsp
J»jf
50
100
t , s
150
200
250
Рис. 6. Переходной процесс по заданию для запаздывания т = 15. Настройка по критерию (2)
1,4
1,2 & 1
.о ^ 0,8
0,6
0,4
0,2
0
0
50
100
150
200
250
t, с
Рис. 7. Переходной процесс по заданию для запаздывания т = 3. Настройка по критерию (2)
1,6
1,4 &1,2
£ 1
0,8 0,6 0,4 0,2 0
KIDF
» y крДсс
ft/r^L
(/ ■*"
t \kid.}
0
50
100
150
t, с
200
250
300
Рис. 8. Переходной процесс по заданию для запаздывания т = 15. Настройка по критерию (1)
0,25 0,2
,0,15
S3 , о ^ 0,1
0,05
0
-0,05 -0,1
8гЛ KIDsp
1 uV ГАРДСС
V #v'\ _/cX
\KiqF
50
100 150
t, с
200
250
300
Рис. 9. Переходной процесс по возмущению для запаздывания т = 0. Настройка по критерию (1)
Следует отметить, что отдельно были проверены варианты использования 2БоР для ПИ-регулятора. Во многом они совпадают с результатами анализа работы ПИД-регулятора.
Методика последовательной настройки структуры 2БоБ ПИД-регулятора путем использования интегрального критерия дает устойчивый вариант комплексной настройки. Одновременная настройка всех пяти коэффициентов позволит добиться большего баланса. Но в этом случае потребуется сформировать обобщеный критерий для канала по возмущению и сигналу задания с весовыми коэффициентами. Вариация этих весовых коэфициентов позволит перераспределять эффективность функционирования исходя из требований прикладной технологии.
Заключение. В результате проведенных исследований сделаем следующее выводы:
1. При использовании ПИД-регулятора для стабилизации технологических параметров следует заранее определить основные источники возмущений и на их основании производить выбор структуры и настройку регулятора degree-of-freedom.
2. Предложенные нелинейные варианты two-degree-of-freedom чаще всего обеспечивают лучшее качество функционирования контуров с обратными связями.
Литература
1. Aström K. J., Hägglund T. Advanced PID control. ISA - the Instrumentation, Systems, and Automation Society. Durham: Research Triangle Park, 2006. 461 p.
2. Alfaro V. M., Vilanova R. Model-Reference Robust Tuning of PID Controllers. Cham: Springer International Publishing Switzerland, 2016. 202 p.
3. Денисенко В. В. Компьютерное управление технологическим процессом, экспериментом, оборудованием. М.: Горячая линия - Телеком, 2009. 608 с.
Y.
OUT
0
0
0
4. O'Dwyer A. Handbook of PI and PID controller tuning rules. London: Imperial College Press, 2009, 624 p.
5. Viteckova M., Vitecek A. 2DOF PI and PID Controllers Tuning // IFAC Proceedings. 2010. Vol. 43, issue 2. P. 343-348.
6. Latha K., Rajinikanth V. 2DOF PID controller tuning for unstable systems using bacterial foraging algorithm // SEMCCO 2012: Swarm, Evolutionary, and Memetic Computing, Lecture Notes in Computer Science. 2012. P. 519-527. DOI: 0.1007/978-3-642-35380-2_61.
7. Multi-objective optimization based tuning tool for industrial 2DoF PID controllers / C. Gamboa [et al.] // IFAC-PapersOnLine, July 2017. Vol. 50, issue 1. P. 7511-7516.
8. Hryniuk D., Suhorukova I., Orobei I. Non-linear PID controller and methods of its setting // Open Conference of Electrical, Electronic and Information Sciences (eStream), 2017, pp. 1-4. DOI: 10.1109/ eStream.2017.7950327.
9. Гринюк Д. А., Оробей И. О., Сухорукова И. Г. Модификация интегральных критериев для повышения запаса по устойчивости // Труды БГТУ. 2012. № 6: Физ.-мат. науки и информатика. C. 118-121.
10. Гринюк Д. А., Оробей И. О., Кузьмицкий И. Ф. Численное исследование алгоритмов уменьшения интегрального насыщения // Труды БГТУ. Сер. VI, Физ.-мат. науки и информатика. 2005. Вып. ХШ. С. 140-143.
References
1. Astrom K. J., Hagglund T. Advanced PID control. ISA - the Instrumentation, Systems, and Automation Society. Durham, Research Triangle Park, 2006. 461 p.
2. Alfaro V. M., Vilanova R. Model-Reference Robust Tuning of PID Controllers. Cham, Springer International Publishing Switzerland, 2016. 202 p.
3. Denisenko V. V. Komp'yuternoe upravlenie tekhnologicheskim protsessom, eksperimentom, oborudovaniem [Computer control of technological process, experiment, equipment]. Moscow, Goryachaya liniya - Telekom Publ., 2009. 608 p.
4. O'Dwyer A. Handbook of PI and PID controller tuning rules. London, Imperial College Press, 2009, 624 p.
5. Viteckova M., Vitecek A. 2DOF PI and PID Controllers Tuning. IFAC Proceedings, 2010, vol. 43, issue 2, pp. 343-348.
6. Latha K., Rajinikanth V. 2DOF PID controller tuning for unstable systems using bacterial foraging algorithm. SEMCCO 2012: Swarm, Evolutionary, and Memetic Computing, Lecture Notes in Computer Science, 2012, pp. 519-527. DOI: 0.1007/978-3-642-35380-2_61.
7. Gamboa C., Rojas J. D., Arrieta O., Vilanova R. Multi-objective optimization based tuning tool for industrial 2DoF PID controllers. IFAC-PapersOnLine, July 2017, vol. 50, issue 1, pp.7511-7516.
8. Hryniuk D., Suhorukova I., Orobei I. Non-linear PID controller and methods of its setting. Open Conference of Electrical, Electronic and Information Sciences (eStream), 2017, pp. 1-4. DOI: 10.1109/eStream.2017.7950327.
9. Hryniuk D. A., Orobey I. O., Sukhorukova I. G. Modification of integral criteria for increasing the margin of stability. Trudy BGTU [Proceedings of BSTU], 2012, no. 6, Physics and Mathematics. Informatics, pp. 118-121 (In Russian).
10. Hryniuk D. A., Orobey I. O., Kuz'mitskiy I. F. Numerical study of algorithms for decreasing integral saturation. Trudy BGTU [Proceedings of BSTU], series VI, Physics and Mathematics. Informatics, 2005, issue ХШ, pp. 140-143 (In Russian).
Информация об авторах
Гринюк Дмитрий Анатольевич - кандидат технических наук, доцент, доцент кафедры автоматизации производственных процессов и электротехники. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). Е-mail: [email protected]
Сухорукова Ирина Геннадьевна - ассистент кафедры информационных систем и технологий. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). Е-mail: [email protected]
Олиферович Надежда Михайловна - ассистент кафедры автоматизации производственных процессов и электротехники. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13 а, Республика Беларусь). Е-mail: [email protected]
Оробей Игорь Олегович - кандидат технических наук, доцент, доцент кафедры автоматизации производственных процессов и электротехники. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). Е-mail: [email protected]
Information about the authors
Hryniuk Dmitryy Anatol'yevich - PhD (Engineering), Associate Professor, Assistant Professor, the Department of Automation of Production Processes and Electrical Engineering. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: [email protected]
Suhorukova Irina Gennad'yevna - assistant lecturer, the Department of Information Systems and Technology. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: [email protected]
Oliferovich Nadezhda Mikhaylovna - assistant lecturer, the Department of Automation of Production Processes and Electrical Engineering. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: [email protected]
Orobei Igor' Olegovich - PhD (Engineering), Associate Professor, Assistant Professor, the Department of Automation of Production Processes and Electrical Engineering. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: [email protected]
Поступила 15.05.2018
