Системный анализ и управление
УДК 681.58
В.Н. Круглов, Л.Г. Доросинский, Ю.В. Чирышев
разработка модели следящей системы контроля местоположения изображения струи в пределах телевизионного кадра
V.N. Krouglov, L.G. Dorosinsky, Yu.V. Chiryshev
DEVELOPMENT OF SERVO SYSTEM MODEL TO CONTROL STREAM POSITION wITHIN VIDEO Frame
Рассмотрена структурная схема модели следящей системы, приведены расчетные методики для определения параметров входных воздействий и следящего привода на базе электромеханического шагового двигателя. Изучена математическая модель следящей системы, даны результаты динамического имитационного моделирования с применением программного пакета Simulink Response Optimization.
СЛЕДЯЩАЯ СИСТЕМА; МОДЕЛЬ; ПЛАВИЛЬНАЯ ПЕЧЬ; ШАГОВЫЙ ДВИГАТЕЛЬ; ПИД-РЕГУЛЯТОР; ТЕЛЕВИЗИОННАЯ КАМЕРА; СИСТЕМА ТЕХНИЧЕСКОГО ЗРЕНИЯ; РОБАСТ-НОСТЬ.
The block diagram and mathematical model of servo-system is considered and calculation methods for determination of parameters of input actions and servo-system on the base of step-by-step motor are given. The dynamic simulation method is used for the system description. Optimal controller parameters were found using the Simulink Response Optimization toolbox.
SERVO-SYSTEM; MODEL; MELTING FURNACE; STEPPER MOTOR; PID CONTROLLER; TV CAMERA; COMPUTER VISION SYSTEM; ROBUSTNESS.
Анализ проблемы энергосбережения показывает, что один из наиболее эффективных путей ее решения заключается в сокращении потерь тепла через ограждающие конструкции зданий, сооружений и тепловых сетей. Этим объясняется интенсивное развитие производства теплоизоляционных материалов, в частности, минеральной ваты. При ее производстве одним из основных контролируемых параметров является дебит расплава, вытекающего из плавильного агрегата.
В настоящее время наиболее перспективными средствами измерения дебита расплава металлургических печей, работающих в непрерывном режиме технологического цик-
ла, являются приборы на основе систем технического зрения [1]. Особенность процесса истечения высокотемпературной струи расплава заключается в постоянном изменении ее местоположения. В зарубежных аналогах таких систем, например системе фирмы Gamma Meccanica SpA (Италия), проблема определения местоположения струи в пределах видеокадра решается за счет увеличения угла поля зрения объектива видеокамеры. Однако это приводит к ухудшению разрешающей способности оптической системы и, как следствие, к увеличению погрешности оценки диаметра струи и ее дебита [2]. Вот почему для сохранения необходимой разрешающей способности оптической си-
стемы необходимо включение в состав рассматриваемых приборов следящей системы, контролирующей местоположение струи, а оценка параметров такой системы является актуальной задачей.
В статье предлагается модель следящей системы, предназначенной для поддержания углового значения линии визирования телевизионной камеры, направленной на струю расплава, вытекающую из плавильной печи. Управлению подлежит угол поворота линии визирования видеокамеры, другими словами, система должна управлять углом поворота вала двигателя, связанного с видеокамерой.
В процессе проектирования следящей системы необходимо учитывать широкий круг технических требований, влияющих в конечном счете на единый критерий качества системы. К названным требованиям следует отнести:
динамические характеристики, определяющие ошибки в системе управления при различных режимах работы;
надежность работы системы, устойчивость к влиянию внешних воздействий; характер эксплуатации системы; допустимое потребление энергии и т. д. Для исследования различных характеристик системы должна использоваться компьютерная модель, т. к. в процессе проектирования возникнет потребность в проведении числовых экспериментов в условиях, приближенных к реальным. В результате того что модель окажется в тех же условиях и будет подвергнута тем же внешним воздействиям, при которых будет работать реальная система, мы сможем получить количественные оценки поведения системы.
Адекватным подходом к исследованию следящей системы, находящейся под влиянием случайных воздействий, является применение статистических методов, которые дают возможность учесть случайный характер воздействия.
Для проектируемой следящей системы критерием оптимального управления будем считать такой режим, когда установившееся значение дисперсии ошибки минимально.
Характер входного воздействия проектируемой системы со временем меняется
и в общем случае представляет собой нестационарный случайный процесс. Есть интервалы времени, когда струя течет «спокойно», при этом можно достоверно определить статистические характеристики этого случайного процесса. Когда прочищается выпускное отверстие плавильной печи (летка), входное воздействие системы представляет собой случайный процесс с параметрами, изменяющимися во времени.
«Спокойное» поведение струи не представляет практического интереса, поскольку дисперсия случайного процесса является очень малой величиной. струя расплава постоянно находится в кадре. Следить за струей нужно в те периоды времени, когда струя отклоняется от центра кадра или движется за пределы видимости кадра, например, при расширении летки. Поведение струи в такие моменты представляет собой нестационарный процесс.
При использовании для целей динамического синтеза методов оптимальной фильтрации необходимо априорное знание статистических характеристик полезного сигнала и помехи. Построение простой, достаточно достоверной модели входного воздействия в данном случае сделать не удается, т. к. анализ свойств нестационарных случайных процессов представляет собой сложную задачу. Переход к структуре, содержащей устройства адаптации, лишь усложнит систему. Вот почему представляется целесообразным построение робастной системы управления, которая обеспечит требуемое качество слежения при изменении статистических характеристик входных воздействий в широких пределах.
В качестве исполнительного элемента в следящих системах широкое распространение получили шаговые двигатели (ШД). Связано это с тем, что они обладают некоторыми свойствами, которые делают их исключительно удобными для применения:
угол поворота ротора определяется числом импульсов, поданных на двигатель;
точность составляет 3—5 % от величины шага, причем ошибка не накапливается от шага к шагу;
старт/остановка/реверсирование осуществляются практически безынерционно;
надежность определяется только сроком службы подшипников.
Применение шаговых приводов обусловлено еще и тем обстоятельством, что они естественным образом сочетаются с цифровыми управляющими машинами, программными устройствами и микропроцессорами. Именно по этим соображениям использование шагового двигателя в качестве исполнительного устройства для рассматриваемой системы является наиболее предпочтительным.
В общем случае поведение т-фазного шагового двигателя описывается системой из т + 1 уравнений, причем т уравнений являются уравнениями электрического равновесия и одно — уравнением движения ротора [3]:
я • т ц) + +Ет = ит а)
йг
й 2а(г) й а(г) и--+ К-
йг
йг
й а(г) (1)
где ит (г) — мгновенное значение напряжения на фазе статора; Я — сопротивление статора; 1т (г) — мгновенное значение тока в фазе; Ь — индуктивность фазы; Ет — наведенная ЭДС в фазной обмотке; а — угол поворота ротора; / — приведенный момент инерции привода; К — коэффициент демпфирования силами вязкого трения; Р — момент сухого трения; Мн(г) — общий момент нагрузки; М (г) — мгновенное значение синхронизирующего момента шагового двигателя.
Дальнейшее упрощение заключается в линеаризации уравнений, переходу к операторной форме и записи передаточной функции шагового двигателя при отработке одного шага в виде Жоу (р) = , где
р — оператор Лапласа; ©0 — преобразование лапласа действительного углового положения; — преобразование Лапласа требуемого значения угла [3, 4].
переходный процесс ротора шагового двигателя при каждом шаге учитывается введением непрерывного звена. передаточ-
ная функция такого звена при единичном ступенчатом воздействии может быть определена следующим образом:
1
^ (р) = =_
н т2 • р2 + 2 т • р +1
, (2)
т 1
где Т = — =
и
Р • М„
е = 9,81 10-Мт щ0 • к ; щ0 - собствен" р 12 я 0 ш ? 0
ная частота малых колебаний ротора около положения равновесия; Р — число пар полюсов двигателя; Мт - максимальный статический синхронизирующий момент двигателя; / — момент инерции привода; кш — коэффициент, учитывающий демпфирование обмоток ШД (изменяется в пределах от 1 до 2); Я — сопротивление фазы статора.
Электромагнитные переходные процессы могут быть учтены звеном чистого запаздывания е~ХДВр. Передаточная функция при этом принимает вид:
(р) =
—ТДВ р
Т2 • р2 + 2 • е • т • р + 1
(3)
где тДВ — постоянная времени, определяемая отношением индуктивности и активного сопротивления цепи обмоток управления.
Перед непрерывным звеном должен быть включен формирующий элемент. В эквивалентных схемах используется формирующий элемент типа фиксатора уровня, что соответствует особенности двигателя затормаживаться после каждого шага. Передаточная функция фиксатора имеет вид:
Ф( р) =
1 - е
Т) р
(4)
Таким образом, эквивалентную схему можно представить в виде последовательного соединения звеньев Ф(р) и Т¥н (р).
Очевидно, реакция модели на любой входной сигнал не должна отличаться от реакции шагового двигателя на этот сигнал, т. к. только при этом условии возможна замена шагового двигателя эквивалентной схемой в общей структурной схеме следящей системы. Реакция 90 на шаговую функцию 9(. приведена на рис. 1.
е
15
§
£ «
05
-1-1-1-1- -
/............I.............~ 1
1 1$ 3
Время
29
Рис. 1. Перемещения ротора
перемещение ротора является колебательным затухающим процессом, амплитуда и время затухания которого зависят от параметров конкретного привода и нагрузки.
Цель синтеза робастной системы заключается в том, чтобы выходная переменная системы за минимальное время и с требуемой точностью воспроизводила все изменения входного сигнала. Любая практическая система содержит элементы емкостного и индуктивного типа, накапливающие энергию. Именно эти элементы в сочетании с элементами, рассеивающими энергию, обусловливают динамику системы слежения. Поэтому цель управления сформулируем следующим образом: амплитудную характеристику системы надо стремиться сделать как можно ближе к единице в как можно большем интервале частот.
В качестве устройства управления используется ПИД-регулятор с передаточной функцией [5]:
^пид( Р) = Кп + ^ + КД Р.
(5)
Его структурная схема приведена на рис. 2.
Во временной области выходная переменная регулятора и(1) и его входная переменная е(/) связаны уравнением:
п(1) = Кпва) + Ки | е(1)А + КД ^. (6)
Выбор трех коэффициентов ПИД-регулятора по сути и есть задача отыскания такого устройства управления, при котором обеспечивались бы желаемые показатели качества системы.
Рис. 2. Структурная схема ПИД-регулятора
Рис. 3. Структурная схема приближенной модели системы управления:
Ф — текущее значение угла поворота исполнительного вала
Система управления цифровым следящим приводом на основе электромеханического шагового двигателя может быть описана структурной схемой, приведенной на рис. 3. Данная модель представляет собой описание каскадной одноконтурной следящей системы автоматического управления позиционированием углового положения видеокамеры. Соответствующая ей модель, построенная из блоков пакета 81ши1тк, показана на рис. 4. В фиксированные моменты времени, следующие с периодом дискретизации Т, управляющая программа
сравнивает текущее положение струи расплава, приведенное к угловому положению вала двигателя, с угловым положением видеокамеры.
Значение Тд определяется аналого-цифровыми преобразователями Quantizer и Quantizeri.
Если струя находится в центре кадра или ее отклонение от середины не превышает определенного порога (струя находится в зоне нечувствительности), блок управления ошибками (Control of errors) игнорирует такое незначительное отклоне-
Рис. 4. Модель системы управления: Step — источник ступенчатого воздействия; Quantizer, Quantizeri — аналогово-цифровые преобразователи; Control of errors — подсистема управления ошибками; PID Controller — дискретный ПИД-регулятор; Stepper Motor — подсистема шагового привода; Scope — осциллограф
Derivative
Рис. 5. Структура цифрового ПИД-регулятора: Integrator — блок интегрирования; Derivative — блок дифференцирования; Kp, Ki, Kd — коэффициенты регулятора
ние и подает на регулятор (PID Controller) нулевую ошибку. При этом на шаговый двигатель не подается управляющего воздействия, он не меняет своего положения. Если программа определяет, что струя находится вне зоны нечувствительности, текущее значение ошибки попадет на вход регулятора, на выходе которого формируется сигнал управления шаговым двигателем для минимизации ошибки отклонения.
Структура цифрового ПИД-регулятора, показанная на рис. 5, отличается от структуры аналогового (см. рис. 2) лишь тем, что операторы дифференцирования и интегрирования в классическом уравнении ПИД-регулятора (6) заменены конечными разностями и конечными суммами. В результате получаем модель управляемого объекта и структуру системы управления. Далее необходимо получить линеаризован-
ную модель, которая будет использоваться для настройки системы управления. Для этого используем пакет Simulink Control Design. В нем выбираем те объекты, которые необходимо настроить. Так как устройством управления данной системы является ПИД-регулятор, он и будет выбран. Затем нужно задать частотные и временные характеристики — требования к системе.
Основным показателем качества работы системы является обеспечение отработки входного воздействия с минимальной ошибкой. Для того чтобы заданное воздействие воспроизводилось с ошибкой не более ^ ЛАХ проектируемой системы не должна заходить в запретную область, ограниченную двумя прямыми с наклонами —20 дБ и —40 дБ [6, 8]. Для определения требований к низкочастотной части лАх необходимо знать максимальное по модулю
Исходные данные
Максимальная скорость входного воздействия, град/с 1,6
Максимальное ускорение входного воздействия град/с2 0,4
Максимально допустимая ошибка слежения е^, град 3
Период дискретизации Т, с 0,3
М
чд
ей £
н н
30
20
10
'«Г
К
10
Частота, рад/с
Рис. 6. Контрольные точки для ЛАХ из условия точности воспроизведения входного сигнала
значение скорости gmax, максимальное по модулю значение ускорения gmax входного сигнала g(t) и максимально допустимую ошибку слежения emax.
Детальный анализ динамических характеристик струи расплава проводился с помощью базы видеоизображений, содержащей более 3500 кадров. В результате проделанной работы были определены скорость и ускорение перемещения струи в горизонтальном направлении. Исходя из условий технологического процесса, дальности расположения видеокамеры от струи расплава, угловых размеров поля зрения объектива задается максимально допустимая ошибка слежения, равная трем градусам. Полученные данные сведены в таблицу.
В соответствии с [6, 7] построена запретная область ЛАХ системы, удовлетворяющей требованию, чтобы максимальная ошибка не превышала заданного значения ^ = 3 град, приведенная на рис. 6.
Кроме этого введем дополнительные ограничения:
время переходного процесса при отработке ступенчатого воздействия не более 1 с; перерегулирование — не более 15 %.
Для нахождения оптимальных параметров ПИД-регулятора использовался пакет численной оптимизации Simulink Response Optimization с предварительно заданными ограничениями на ступенчатый входной сигнал и низкочастотную часть лАх. Достоинства пакета заключаются в том, что оптимизация параметров осуществляется в автоматическом режиме. В результате проведенного моделирования определены оптимальные параметры ПИД-регулятора, которые имеют следующие численные значения:
KP = 0,87, KI = 4,48, KD = 0,001.
В результате проделанной работы получены следующие результаты:
разработана математическая модель системы слежения на основе электромеханического шагового двигателя;
выполнен анализ динамических характеристик процесса истечения струи расплава на базе предприятия ОАО «АКСИ» (г. Челябинск), в результате которого определена максимально допустимая ошибка слежения, не превышающая трех градусов углового отклонения линии ви-
зирования видеокамеры от местоположения струи расплава в горизонтальном направлении;
для выбранного критерия оптимизации
следящей системы путем математического моделирования определены оптимальные коэффициенты ПИД-регулятора следящей системы.
СПИСОК ЛИТЕРАТУРЫ
1. Sirok B., Blagojevich B. Mineral wool. WoodHead Publishing limited, Cambridge, England, 2008. P. 184.
2. круглов В.Н. Способ определения производительности плавильного агрегата // Патент РФ № 2469962. Приоритет 06.07.2011.
3. Гумен В.ф., калининская Т.В. Следящий шаговый электропривод. Л.: ЛО «Энергия», 1980. 168 с.
4. кенио Т. Шаговые двигатели и их микропроцессорные системы управления: Пер. с англ. М.: Энергоатомиздат, 1987. 200 с.
5. Олссон Г., Пиани Д. Цифровые системы автоматизации и управления. СПб.: Невский Диалект, 2001. 557 с.
6. Бесекерский В.А., Изранцев В.В. Системы автоматического управления и микроЭВМ. М.: Наука. Гл. ред. Физматлит, 1987. 320 с.
7. Бесекерский В.А. цифровые автоматические системы. М.: Наука. Гл. ред. Физматлит, 1976. 576 с.
8. коновалов Г.ф. Радиоавтоматика: Учеб. для вузов по спец. «Радиотехника». М.: Высш. шк., 1990. 335 с.
REFERENCES
1. Sirok B., Blagojevich B. Mineral wool. WoodHead Publishing limited, Cambridge, England, 2008, P. 184.
2. Kruglov V.N. Sposob opredeleniya proizvoditel-nosti plavilnogo agregata. Patent RF, No. 2469962 [.Patent for the method of determining the performance of the melting unit], Prioritet 06.07.2011. (rus)
3. Gumen V.F., kalininskaya T.V. Sledyashchiy shagovyy elektroprivod [Tracking stepper actuator], Leningrad: Energiya Publ., 1980, 168 p. (rus)
4. kenio T. Shagovyye dvigateli i ikh mikroprotsessornyye sistemy upravleniya. Moscow:
Energoatomizdat Publ., 1987, 200 p. (rus)
5. Olsson G., Piani D. Tsifrovyye sistemy avtomatizatsii i upravleniya. St. Petersburg: Nevskiy Dialekt Publ., 2001, 557 p. (rus)
6. Besekerskiy V.A., Izrantsev V.V. Sistemy avtomaticheskogo upravleniya i mikroEVM. Moscow: Nauka Publ. Gl. red. Fizmatlit, 1987, 320 p. (rus)
7. Besekerskiy V.A. Tsifrovyye avtomaticheskiye sistemy. Moscow: Nauka Publ. Gl. red. Fizmatlit, 1976, 576 p. (rus)
8. Konovalov G.F. Radioavtomatika. Moscow: Vysshaya shkola Publ., 1990, 335 p. (rus)
кРУГлОВ Василий Николаевич — доцент кафедры информационных технологий Института радиоэлектроники и информационных технологий Уральского федерального университета имени первого Президента России Б.Н. Ельцина, кандидат технических наук.
620002, Россия, г. Екатеринбург, ул. Мира, д. 19.
E-mail: [email protected]
KROUGLOV, Vasilii N. Ural Federal University named after first president of Russia B.N. Yeltsin.
620002, Mira Str. 19, Ekaterinburg, Russia.
E-mail: [email protected]
ДОРОСИНСкИй леонид Григорьевич — профессор кафедры информационных технологий Института радиоэлектроники и информационных технологий Уральского федерального университета имени первого Президента России Б.Н. Ельцина, доктор технических наук.
620002, Россия, г. Екатеринбург, ул. Мира, д. 19.
E-mail: [email protected]
DOROSINSKY, leonid G. Ural Federal University named after first president of Russia B.N. Yeltsin.
620002, Mira Str. 19, Ekaterinburg, Russia.
E-mail: [email protected]
ЧИРыШЕВ Юрий Владимирович — аспирант кафедры информационных технологий Института радиоэлектроники и информационных технологий Уральского федерального университета имени первого Президента России Б.Н. Ельцина.
620002, Россия, г. Екатеринбург, ул. Мира, д. 19. E-mail: [email protected]
CHIRYSHEV, Yurii V. Ural Federal University named after first president of Russia B.N. Yeltsin. 620002, Mira Str. 19, Ekaterinburg, Russia. E-mail: [email protected]
© Санкт-Петербургский государственный политехнический университет, 2014