CC BY
20
DOI: 10.24143/2072-9502-2019-4-46-58 УДК 681.5:517.443
МОДАЛЬНОЕ УПРАВЛЕНИЕ ПРОЦЕССАМИ ДОЗИРОВАНИЯ В СРЕДЕ ПРОСТРАНСТВА СОСТОЯНИЙ И ВЕЙВЛЕТ-ПРЕОБРАЗОВАНИЙ
Д. Б. Федосенков1, А. А. Симикова2, Б. А. Федосенков2
1 Сибирская генерирующая компания, Москва, Российская Федерация
2 Кузбасский государственный технический университет имени Т. Ф. Горбачева, Кемерово, Российская Федерация
На примере системы управления технологическими процессами дозирования в смесепри-готовигельных агрегатах представлены и описаны необходимые положения для проектирования и рабочего функционирования автоматизированных систем модального управления. Объекты исследования - структура блока дозирующего устройства и переходные процессы при дозировании. Предметом исследования являются вопросы проектирования и создания автоматизированной системы модального управления (АСМУ), функционирующей в режиме перестраиваемых параметрически нестационарных по частоте переходных процессов (ПП). При этом скалярные сигналы материалопотоков расхода дозирующих устройств (ДУ) формируются в терминах пространства состояний (в формате внутренних моделей «вход -состояние - выход»), а мониторирование этих процессов производится во время-частотной вейвлет-среде. Последнее позволяет проводить анализ протекающих в АСМУ процессов на семантически прозрачном и информационно насыщенном уровне. Представлены параметры (скважности) Фурье-модели сигналов дискретного дозатора, от которых зависит характер процесса дозирования. Сформирована модель блока дозаторов в терминах пространства состояний. Рассмотрены Ш-сигналы расхода и их отображения в виде время-частотных распределений (^-карт). Разработана новая концепция модального управления расходом блока дозирующих устройств смесеприготовительного агрегата путем введения понятия параметрической нестационарности внутренней среды дозирующих устройств, учитывающей варьирование многомерного коэффициента обратной связи по состоянию в режиме реального времени в ходе переходных процессов в системе, причем с целью экспериментальной оценки параметров (здесь: постоянных времени) переходных процессов последние определяются по отображениям сигналов расхода на ^-картах. Используемая в системе вейвлет-среда оперирует время-частотными распределениями класса Коэна (преобразованиями Вигнера - Вилля и Чуи - Уилльямса). Сочетание моделей ПП в терминах пространства состояний и их представление в формате распределений Вигнера позволяет оценивать, задавать и переназначать параметры внутренней структуры объекта управления (отдельного дозатора или блока дозирующих устройств) в режиме реального времени с целью достижения эффекта стабилизации расхода ингредиентов в номинальном режиме на предсмесительной стадии. Возмущенные режимы представлены в виде процессов, восстанавливающих номинальные режимы (так называемые обратные ПП). Разработана и описана процедура реализации алгоритма модального управления процессом дозирования, в том числе технология принудительной локализации и релокализации полюсов объекта в условиях нестационарных процессов в системе «исполнительный механизм - ДУ», на примере дозатора непрерывного действия в возмущенном режиме при сбросе нагрузки (при нисходящем обратном переходном процессе). Для данного режима представлен расчет нестационарной матрицы модального регулятора.
Ключевые слова: модальное управление, переходные процессы при набросе и сбросе нагрузки, параметрически нестационарные чирп-сигналы, время-частотная вейвлет-среда, распределение Вигнера, ^-карта, матрица модального регулятора.
Для цитирования: Федосенков Д. Б., Симикова А. А., Федосенков Б. А. Модальное управление процессами дозирования в среде пространства состояний и вейвлет-преобразований // Вестник Астраханского государственного технического университета. Серия: Управление, вычислительная техника и информатика. 2019. № 4. С. 46-58. DOI: 10.24143/2072-9502-2019-4-46-58.
Введение
Применение в промышленности непрерывно действующих смесительных аппаратов не получило широкого распространения в силу нерешенности ряда вопросов. В частности, недостаточно полно изучено влияние режимных параметров блока дозирующих устройств на характер и структуру выходных материальных потоков, формируемых на предсмесительной стадии, от которых в конечном счете зависит качество готовой смеси. Кроме того, не установлены и не сформулированы технологические причины возникновения нестабильности процесса дозирования вследствие такого влияния и, как следствие, не разработаны условия и способы введения принудительной параметрической нестационарности режимов работы блока дозаторов с целью динамической стабилизации материалопотоков при смене входных воздействий. Ликвидировать такого рода недостатки можно с помощью модального управления, дающего возможность (путем перестройки параметров замкнутой системы с обратной связью по состоянию) поддерживать выходные потоки блока дозаторов на требуемом номинальном уровне.
С учетом реальных условий возникает необходимость в разработке и адаптации альтернативного метода для моделирования, исследования и усовершенствования разнохарактерных нестационарных по частоте процессов дозирования, который, будучи реализован на компьютерной основе, позволил бы автоматизировать процесс управления режимами работы в целях интенсификации смесеприготовления. Таким методом, дающим возможность проводить высокоточный анализ разнохарактерных процессов, является метод всплескового (вейвлет) преобразования на основе адаптивной аппроксимации исследуемых сигналов базовыми локальными волновыми составляющими (вейвлет-функциями).
Объектом исследования является структура блока дозирующих устройств и процессы дозирования в смесеприготовительном агрегате непрерывного действия.
Предметная область исследования охватывает вопросы, связанные с проектированием и созданием автоматизированной системы модального управления процессами дозирования, функционирующей в режиме перестраиваемой параметрической нестационарности и в условиях мониторирования этих процессов в вейвлет-среде.
Цель работы - разработка и исследование аппаратной и программной основы комплекса для моделирования и текущего мониторинга процессов дозирования в составе смесеприготови-тельного агрегата на основе вейвлет-преобразований, разработка системы автоматизации блока дозирования на базе методов модального управления переходными процессами дозирования.
Научная новизна работы заключается:
- в формировании математических моделей объекта управления (блока дозирующих устройств - БДУ) высокого порядка, включающего дозирующие устройства (ДУ) непрерывного и дискретного действия - в терминах метода пространства состояний;
- создании и реализации концепции модального регулирования расхода БДУ на базе теории пространства состояний в сочетании с аппаратом вейвлет-преобразований;
- разработке процедуры многомерно-точечного отображения одномерных материалопо-токовых сигналов расхода время-частотными вейвлет-распределениями класса Коэна;
- разработке алгоритма динамического пересчета параметров модального регулятора с учетом идентификации режимов дозирования по вейвлет-картам при спорадическом проявлении параметрической нестационарности.
Математические модели сигналов дозирующих устройств непрерывного и дискретного действия
В работе представлены математические модели сигналов расхода дозирующих устройств непрерывного и дискретного действия. Дозаторы непрерывного действия (НДУ) - спиральные и шнековые - описываются моделями третьего порядка.
На основе опытных данных, полученных в ходе исследований порционных дозаторов разных типов (шлюзового, с возвратно-поступательным механизмом, тарельчатого, карусельно-ста-канчатого, а также микродозаторов), в качестве общего подхода был выбран способ точной формализации [1] сигнала порционного расхода. Сигнал порционного, или дискретного, дозирующего устройства (ДДУ), сформированный с учетом такого подхода, представлен на рис. 1.
Xd(t) Td(k(l -1) + 1]
Рис. 1. Параметризация сигнала порционного дозирования для общего режима
При таких условиях поток на выходе из порционного дозатора на протяжении произвольного количества циклов описывается функциональной зависимостью
X (t)=!■
1=1
Xmd(t - (1 - ) при (1 - < t < T W-1) + 1 ,
T
Xm,
Хцу
при Td (Хцу(1 -1) +1) < t < Td (Хц(1 -1) +1)
f Tr (t - Td (1 -1))-1Л
V Td
0
Хцу
Хц
1 - ц
при Td (Хц(1 -1) + 1) < t < Td (Х(1 -1) +1)
при
Хц
Td (1(1 -1) + 1)
Х
< t < 1T.
(1)
где Т - период дозирования; Хта - весовой расход материала через дозатор; k - количество формируемых доз (циклов); I - произвольный номер цикла дозирования; X, ц, V - соответственно значения скважностей: порционного дозирования, интервала формирования дозы до начала отсечки и интервала достижения режима номинального дозирования.
Разложение дозирующего сигнала (рис. 1) в тригонометрический ряд Фурье с произвольным периодом имеет вид:
A
2kn
. 2kn
— + V A cos-1+Bk sin-1
/"y ^^ k T k rf'
T
d J
2 Т"
A = ^ J X
t., 2 ч
d 0 Td
A = - J Xd (t)
T J
cos
d0 Td
f Л V Td J
dt,
D 2 ( v .. . 2kn Bk = —J Xd (t)si
T
)sin
d0
V Td J
dt,
t
где А0, Лк, - коэффициенты Фурье-модели; t - текущее время; Х<() - сигнал расхода порционного дозатора.
Определены коэффициенты Фурье-разложения как функции режимных параметров дозирования (X, ц, V, Td, Xmd) для наиболее типичных режимов работы порционных дозаторов.
Итоговая расчетная модель блока дозаторов (рис. 2), состоящая из спирального и шнеко-вого ДУ и трех ДУ порционного типа (ПДУ), рассматривается как векторная динамическая система в пространстве состояний.
Рис. 2. Блок ДУ в виде структурного фрагмента схемы моделирования смесеприготовительного агрегата в технологическом пространстве состояний
При моделировании процессов в блоке ДУ в пространстве состояний [2, 3] на основе операторных функций звеньев составляются система дифференциальных уравнений, которая затем преобразуется в систему дифференциальных уравнений состояния, записанных в нормальной форме Коши, и линейное алгебраическое уравнение выхода, выражающее выходной сигнал системы через комбинацию переменных состояния и сигнала управления. Представляя набор переменных в виде векторов, а набор параметров в виде соответствующих матриц, получаем векторно-матричную модель (ВММ) агрегата в пространстве состояний.
Блок дозирования на структурном уровне представлен системой SISO (англ. single input -single output). Для численной обработки сигналов и реализации функций управления использовались персональные компьютеры ППЭВМ IBM-совместимого типа. В качестве базового программного обеспечения были приняты универсальные интегрированные пакеты: Mattaad 7.0 Professional, Maple V Release 4.0, MathLab 5.2, а также инструментальный пакет Classic 3.0.
Итоговая модель дозирующего сигнала спирального дозатора (рис. 2) в терминах пространства состояний [4] имеет следующий вид:
X 1( t) = X 2( t),
x 2 (t) = - га d 2 x 1( t) + X
x 3 (t) = X d o1 u (tX u 1 (t) = x 1 (t) + x 3 ( t).
dm i ^ d i
u ( t ),
Векторно-матричная модель дозирующего сигнала шнекового дозатора соответственно выглядит так:
X 4 ( t ) = X 5( t ),
х5(t) =- ш л2 х4(t) + х т 2ш Л2и (tX X б( t) = X т 2 и (t), и 2( t ) = X 4( t ) + X 6( t ).
Дозирующий сигнал 1-го порционного ДУ, представленный Фурье-моделью, соответствующей общему режиму загрузки смесителя (см. рис. 1), характеризуется системой дифференциальных уравнений в форме Коши, которая в пространстве состояний примет вид
x 7( t) = А u (t), x 21 (t) = x 21 + 1( t)
+ i(t) = -((q - 1)юj)2 • x21(t) + u(t)
- (1 - ( -1) q )
q = 2,11 , l = -i-(-^ + 2 q,
ub (t)
7(t) +X x2i(t) , i = 4,12
где I - нормированный индекс переменной состояния; q - переопределенный номер полной гармоники.
x
По ВММ сформирована полная расчетная структурная схема БДУ в соответствии с рис. 2.
Получаемые первичные осциллограммы сигналов зашумлены и нестационарны по частоте, следовательно, в любой момент времени обладают случайно меняющимся мгновенным динамическим спектром. Отобразить такие спектры на визуальном уровне существующими традиционными способами не представляется возможным.
В связи с этим в работе реализован новый подход, в основе которого лежат математические модели и алгоритмы аппроксимации, идентификации и коррекции контролируемых переменных в определенном пространстве. Последнее осуществляется на платформе вейвлет-функций и сформированных на их основе время-частотных вейвлет-тезаурусов Габора [5]. При этом используемая информация о состоянии объекта отображается в вейвлет-среде в формате двумерных/трехмерных квадратичных распределений класса Коэна [6, 7]. В качестве рабочих использовались распределения Вигнера и Чуи - Уилльямса [5, 7, 8].
Суть процедуры управляющего мониторинга сводится к следующему.
При изменении структуры спектра временного вектора сигнала изменяется время-частотная локализация соответствующих элементов (атомов) на карте Вигнера [8], что должно быть визуализировано и отслежено мониторинговым комплексом, после чего формируется управляющее воздействие на электропривод дозатора.
Учитывая, что в БДУ отдельные дозаторы спорадически меняют свои выходные сигналы (вследствие варьирования входных воздействий за счет спонтанного изменения физико-механических характеристик поступающих в ДУ ингредиентов), представляется целесообразным наблюдать за общим выходным сигналом расхода БДУ, на вейвлет-карте которого запечатлены все сигналы отдельных ДУ. При этом на мониторе формируются «номинальные» фреймы (специфические окна-рамки - рис. 3), на которые проецируются время-частотные атомы вейвлет-распределения (например, распределения Вигнера), соответствующие номинальному режиму дозирования.
г\с 20
15
10
5
0 Гц
30
20 10
f = 1 TD 0
г\с 40
30 20 10 0
Гц 0
30
fj+u = 2 Гц
0 0,5 1 1,5 2 2,5 3 3,5 4 с 0 0,5 1 1,5 2 2,5 3 с
Atj = 1 c Atj+u = 0,5 c '
б
а
Рис. 3. Осциллограммы и их время-частотные карты, соответствующие номинальному (а) и текущему (возмущенному) (б) режимам работы дозатора порционного типа
Для НДУ спирального и шнекового типа такие распределения со стабильной частотой представляют собой тонкие по частоте (по оси ординат) и вытянутые по времени (по оси абсцисс) время-частотные атомы (ВЧА) карты Вигнера. Поэтому в реальных условиях работы дозаторов на вейвлет-карте (Ж-карте) происходит спорадический выход (девиация) того или иного ВЧА (для НДУ) или группы атомов (для ДДУ) за границы фрейма/фреймов. При этом автоматизированной системой фиксируется очередность нарушения номинальных режимов тех или иных ДУ. Поэтому время обработки Ж-карты соответствующего ДУ должно быть меньше промежутка времени между актами возникновения двух смежных по времени девиаций ВЧА за границы фреймов двух каких-либо разных ДУ либо того же дозатора.
Экспериментально установлено, что такие временные промежутки (/-гэпы) обычно составляют не менее 20 периодов для НДУ (при T = 0,5-1,0 с это составит 10-20 с) и 15 периодов для ДДУ (при Tmax = 1,0 с, т. е. fmm = 1 Гц, это равно 15 с; для 10-й гармоники имеем 10 Гц, при этом г-гэп составит 1,5 с). Подобная оценка распределений случайной величины «/-гэп» для обоих типов ДУ произведена в терминах отображения квартальных распределений в формате «Box&Whiskers». Следовательно, тактовая частота CPU современного PC, используемого в автоматизированной системе для обработки сигналов, может быть сколь угодно малой.
При фиксации недопустимой девиации ВЧА за границы фреймов начинает действовать технология модального управления.
Управление процессами дозирования в автоматизированной системе производится в вейвлет-среде на основе алгоритма модального управления. Цель алгоритма - стабилизировать процесс дозирования на предсмесительной стадии на номинальном уровне при флуктуаци-ях выходных потоков вследствие спорадических ступенчатых или частотно-время-зависимых изменений входных воздействий на дозаторы.
Суть проектирования системы модального управления заключается в синтезе замкнутой системы автоматического управления (ЗСАУ) с обратной связью (ОС) по полному вектору состояния путем задаваемого размещения полюсов объекта с последующим определением матрицы ОС (модального регулятора), формирующей задаваемую локализацию созвездия полюсов.
Задача синтеза заключается в определении желаемого положения корней характеристического уравнения системы и нахождении коэффициентов Ki (параметров модального регулятора), обеспечивающих заданное размещение корней (см. раздел «Алгоритм пересчета матрицы (параметров) модального регулятора К»).
Модель блока ДУ, состоящего из двух дозаторов - спирального и порционного дискретного (последний формирует материалопоток, задаваемый функцией расхода П-образного вида со скважностью X = T / т = 2, где T - период дозирования; т - длительность формирования дозы), -задается соответствующими импульсными переходными функциями. В редуцированном М-графе (сигнальном графе Мэйсона) [1, 3] отражены все составляющие изображений сигналов, что позволяет промоделировать полученную структуру в любой инструментальной системе (здесь моделирование проведено в пакете Classic).
Результаты моделирования и исследования реальных режимов работы автоматизированной системы модального управления процессом дозирования
Так как бункерные емкости каждого ДУ заполнены материалами исходных ингредиентов, обладающими в каждой емкости различными объемными физико-механическими свойствами (например, влажностными характеристиками и/или степенью однородности по своему составу), в процессе дозирования происходит спорадическое варьирование воздействия материала дозируемого ингредиента на входную (заборную) часть активного органа соответствующего дозатора. Вследствие этого происходит смена частоты дозирующего сигнала (т. е. частоты вращения исполнительного механизма (ИМ) дозатора ДУ) в переходном режиме, что в итоге приводит процесс дозирования к новому установившемуся частотному режиму, который не соответствует номинальному, требуемому на предсмесительной стадии.
При поступлении более однородного и сухого (менее увлажненного) вещества ИМ ДУ раскручивается благодаря сбросу нагрузки на рабочий орган дозатора, а частота сигнала материа-лопотока при этом возрастает до тех пор, пока не стабилизируется на новом установившемся уровне с новой, неноминальной, производительностью. При повышении же нагрузки на рабочий орган ДУ картина противоположна. В том и другом случаях требуется воздействовать на замкнутую САУ (ЗСАУ) с модальным управлением с целью стабилизации материального потока на номинальном уровне. Но поскольку воздействие на входе и реакция на выходе ДУ уже стабилизировались на новом уровне, требуется воздействовать на внутреннюю структуру ЗСАУ таким образом, чтобы изменить параметрические свойства системы. Таким образом, требуется создать эффект принудительной параметрической нестационарности за счет смены локализации «созвездия» полюсов ЗСАУ.
Следовательно, при раскрутке ИМ ДУ (т. е. при повышении частоты сигнала в результате ¿-сброса на рабочий орган ДУ) необходимо - с целью стабилизации режима расхода на номинальном уровне - снизить общий коэффициент передачи ЗСАУ, а также перманентно уменьшать значение мнимой сопряженной пары полюсов (для НДУ) в ходе переходного процесса (1111). Что касается обратного 1111 дозирования в ПДУ, то снижать значение мнимых полюсов следует только для наиболее мощных субгармоник первой гармоники Фурье-модели. Модель при этом имеет вид В^т ю^ + ^соб ю^. Режекция высших гармоник выполняется средствами низкочастотной фильтрации или с помощью высокочастотного режекторного фильтра.
При реализации обратного ПП после ¿-наброса средствами модального управления восстановление номинального режима дозирования производится по двухэкспоненциальной зависимости вида
х(0ь/И = х(()оь/<п + хфтЬ/а - Ть) - фо],
где и =Л [1^ - Ть)] + Ле -Л) ([1^ - Те)] - ехр[- ^ - Ть) /Тз] Тз / (Тз - Г4) + ехр[ - ^ - Ть) Т Г4 / (Г3 - Т) -изменение частоты чирп-сигнала [5] восходящего обратного ПП; индексы Ь и е обозначают начало и конец переходного процесса; ■ - номер дозатора; х(()тье = хтЫ - t / Т; ть < t > Те -изменение амплитуды переменной составляющей ВОПП; x(t)0Ь/eг■ = хоь [1(t - ты)] + + (хое, - хоы) да - Те)] - ехр[- (t - Ть) / Тз] Тз / (Тз - Т4) + ехр[- (t - Ть) / Т4] х Т4 / (Т3 - Т4)) -изменение апериодической составляющей восходящего обратного ПП.
Таким образом, полный график переходных процессов при двух циклах варьирования нагрузки (¿-наброса и ¿-сброса), включающий нисходящий прямой/восходящий обратный переходные процессы (последний восстанавливает номинальный режим дозирования), а также восходящий прямой ПП (ВППП) и нисходящий обратный ПП (НОПП), имеет вид (рис. 4):
а
б в
Рис. 4. Циклы модального управления процессом дозирования при ¿-набросе и ¿-сбросе: два цикла ПП (а); ВОПП - сигнал расхода в ходе восходящего обратного ПП (б); НОПП - сигнал расхода при нисходящем обратном ПП (в); НППП - нисходящий прямой ПП
На вейвлет-карте восходящий обратный переходный процесс (ВОПП) для ДУ непрерывного действия после Х-наброса отображается в виде распределения Вигнера (рис. 5).
ш
6,13
3,24 .
2,50' 2,03 0,99 0,82_
0 10 20 30 с
Рис. 5. Время-частотное отображение восходящего обратного переходного процесса, восстанавливающего требуемый режим дозирования (частота дозирования ^)шг = 3,24^6,13 Гц)
Для демонстрации разрешающей способности Ж-карт на рис. 6 представлены осциллограмма регистрируемого материалопотокового сигнала и ее Ж-карта (карта Вигнера) для БДУ, состоящего из одного порционного, двух спиральных, работающих непрерывно, с рабочими частотами 0,20; 3,23 и 4,02 Гц, и спирального дозатора, функционирующего дискретно во времени с заполняющей в пределах интервала дозы частотой 6,89 Гц и частотой дозирования 0,077 Гц.
г/с
Гц 0 5 10 15 20 25 30 с
10
7,5
5
2,5
0
0 5 10 15 20 25 30 с
Рис. 6. Материалопотоковый сигнал БДУ и его отображение на Ж-карте Вигнера
На рис. 7, а показана осциллограмма сигнала расхода ПДУ в виде чирп-сигнала с время-зависимой мгновенной падающей частотой в режиме наброса (увеличения) нагрузки, который в аналитическом виде может быть записан следующим образом:
где а - 2 Хатах
X
• Л _ * тах
- 5 а —
k ■ п
sm
х(0 = ДЛ, Ак, Вк, ю^
X
^ 2 X,
• о __а тах
• В -
ю k -
2п k Т
2
X
Рис. 7. Динамика сигналов порционного дозирующего устройства в режиме наброса нагрузки в одномерном и многомерном виде: исходный сигнал расхода (а); сигнал, восстановленный по спектральным составляющим в виде вейвлет-ряда (б); время-частотное распределение сигнала (в)
После разложения сигнала по алгоритму вейвлет-поиска соответствия [1, 5] получаем его время-частотное распределение (рис. 7, в), адекватность которого установлена путем сравнения исходного сигнала (рис. 7, а) с восстановленным по спектральным составляющим (рис. 7, б). При этом ошибка аппроксимации сигнала расхода в максимуме и на паузе составляет менее 4,3 %, а в точках Гиббса (в точках разрыва первого рода) - 8,9-9,2 %.
Рис. 8 демонстрирует материалопотоковый сигнал двух дозаторов непрерывного действия, работающих в режимах наброса и сброса нагрузки, реконструированный по ВПС-алгоритму сигнал, ошибку аппроксимации и карту Вигнера данного процесса.
Рис. 8. Переходные режимы работы двух дозаторов непрерывного действия при набросе и сбросе нагрузки: одномерный суммарный сигнал расхода (а); его аппроксимация по ВПС-алгоритму (б); ошибка аппроксимации (в); многомерные распределения на Ж-карте (г)
На карте Вигнера этого режима ниспадающая цепочка время-частотных атомов соответствует прямому переходному процессу (НППП) при набросе нагрузки в одном дозаторе, а восходящая цепочка - прямому переходному процессу (ВППП) в другом.
Следует отметить, что все переходные процессы (1111) в объектах (ДУ) являются параметрически нестационарными, т. е. протекают с время-зависимыми мгновенными частотами в виде так называемых чирп-сигналов. При этом нестационарность определяется частотой /шг ПП и меняющейся амплитудой Хтшг переменной составляющей процесса. Первый параметр меняется под действием непрерывно варьируемых полюсов Sp2 и Sp3 (мнимых сопряженных), а второй параметр Хта - в зависимости от соотношения/аг и постоянных времени в произвольный момент.
Алгоритм пересчета матрицы (параметров) модального регулятора К
Процедура реализации алгоритма модального управления процессом дозирования, в том числе технология принудительной локализации и релокализации полюсов разомкнутой САУ в условиях нестационарных процессов в системе «исполнительный механизм - ДУ» (ИМ - ДУ), рассмотрена на примере НДУ в возмущенном режиме при сбросе нагрузки (при переходном процессе типа НОПП). Подобную процедуру целесообразно применять и в многомерно-точечных системах, использующих вейвлет-среду как средство эффективного отображения одномерных сигналов. Расчет параметров модального регулятора, т. е. определение матрицы обратной связи К по полному вектору переменных состояния, сводится к ее расчету по формуле Аккермана [9]
К -[0 0 0 1]Q-Cí ■ а(А),
где QC-1 - обратная матрица управляемости: а(А) - матричный полином, сформированный согласно теореме Кэйли-Хэмилтона [2, 3].
Желаемый характеристический полином определяется по формуле а(з) = - а^3 + а^2 -
- + *4, где = (¿р1 + sp2 + sp3 + sp4), где spj - полюсы объекта; *2 = [¿р ¿р2 sp3 х (sp1 + sp2) + sp4 х
x(sp1 + ¿р2 + sp3 + ¿р4)]; *3 = [sp1 sp2 sp3 + sp1 ¿р2 sp3 sp4 ^р1 + sp2)]; *4 = sp1 sp2 sp3 ¿р4-Г'------------- -----------^---------------------„./• — А 4 л л3 I Л л.
Следовательно, матричный полином имеет вид а(А) I - единичная матрица.
*1 А + *2 А - *3 А + *4 I, где
а(А) -
0 -Т1ю2 +1 -ю2 +1 0 -ю2(Г2 +Г1) -ю2 +1 0 юТ1 -ю2(Т2 +Г1) ю2(ю2 -1)
-Т-1 +1 (-ю2 +1)+(Г2 -Г1) ю2(Г! -1)-Г3 +(Г2 -Г1)+1
0 ю2(ю2 -1)(Г2 +Г1) ю2(ю2 -1) ю2(ю2 -Г2)-ю^Т4 -1)+Г4 -Т-3 +(Г2 -Г1)
По виду а(А) нетрудно заметить, что в системе РСАУ влияние на переходный процесс имеют полюсы ¿р2, ¿р3 и ¿р4 (два мнимых сопряженных и вещественный отрицательный), соответствующие параметрам ю и Т, которые создают частотную составляющую и определяют быстродействие переходного процесса.
В результате получаем матрицу К модального регулятора вида
К = [0 0 0 1] Qc~1 а(А) = [0 ((Г1®2 - 1) |(ю2 - 1) КГ1 - 1)] = [К! \К2 |К3 К].
Данная матрица непрерывно пересчитывается в реальном масштабе времени с произвольным шагом миллисекундного диапазона в ходе переходного процесса, чтобы обеспечить его время-частотный (чирп-) характер.
По аналогичным процедурам автоматизированная система рассчитывает параметры модального регулятора для любого ДУ или блока дозаторов.
Полученные параметры позволяют задавать требуемые нестационарные полюсы замкнутой САУ с полной обратной связью по состоянию, которые определяет автоматизированная система модального управления по разработанному алгоритму.
Таким образом, разработанная автоматизированная система модального управления процессами дозирования выполняет функции стабилизации номинального режима работы БДУ на предсмесительной стадии, что способствует получению смесей высокого качества.
Заключение
1. Разработаны математические модели процессов дозирования в терминах технологического пространства состояний блока дозирующих устройств, с помощью которых можно эффективно отрабатывать широкий спектр текущих режимов мультиингредиентного дозирования устройствами непрерывного и дискретного действия.
2. Сформирована и отработана концепция способа модального управления расходом блока дозирующих устройств (объекта управления) на базе теории пространства состояний в сочетании с аппаратом вейвлет-преобразований, реализованная в замкнутой системе регулирования на основе принудительной параметрической нестационарности обратной связи по состоянию. В условиях спорадического нарушения номинальных потоков дозирования, возникающего вследствие изменения входных воздействий, это позволяет варьировать вектор внутренних параметров объекта, приводя выходные потоки последнего к требуемому состоянию.
3. Представлен алгоритм динамического пересчета параметров многомерной обратной связи по состоянию объекта управления в ходе процесса дозирования (в режиме реального времени), обеспечивающий стабильность материалопотока на предсмесительной стадии.
4. В качестве базовых - при проведении моделирования, расчетов и при обработке сигналов и их отображений - использованы вейвлет-функции Габора, адекватно соответствующие характеру технологических сигналов в подсистеме дозирования и в системе смесеприготовления в целом.
5. Испытания разработанной автоматизированной системы модального управления процессом дозирования показали, что она позволяет интенсифицировать процесс смесеприготовления за счет гармонизации материалопотоков на предсмесительной стадии на базе модального управления в пространстве состояний и вейвлет-преобразований, а также повысить точность дозирования по сравнению с существующими производственными системами. Последний фактор ведет к снижению коэффициента неоднородности конечных смесей, т. е. к повышению их качества.
1. Федосенков Б. А. Научно-технические основы создания и моделирования автоматизированных систем управления непрерывными смесеприготовительными процессами: автореф/ дис. ... д-ра техн. наук. М.: Издат. комплекс МГУПП, 2005. 55 с.
2. Ким Д. П. Теория автоматического управления. М.: Физматлит, 2004. Т. 2. Многомерные, нелинейные, оптимальные и адаптивные системы. 464 с.
3. Филлипс Ч. Л., Харбор Р. Д. Системы управления с обратной связью. М.: Лаборатория базовых знаний, 2001. 616 с.
4. Федосенков Д. Б., Симикова А. А., Федосенков Б. А. Синтез системы модального управления блоком дозирующих устройств // Машиностроение: инновационные аспекты развития: материалы Междунар. науч.-практ. конф. СПб.: СПбФ НИЦ МС, 2019. № 2. С. 162-163.
5. Малла С. Вейвлеты в обработке сигналов. М.: Мир, 2005. 673 с.
6. Cohen L. Time-frequency analysis. Englewood Cliffs: Prentice Hall, 1995. 299 p.
7. Auger F., Chassande-Mottin E. Quadratic time-frequency analysis I: Cohen's class // Time-frequency analysis: concepts and methods. ISTE. 2008 (January). P. 131-163.
8. Debnath L. Recent development in the Wigner-Ville distribution and time-frequency signal analysis // PINSA. 2002, (January). 68, A. N. 1. P. 35-56.
9. Ackermann J. E. Der Entwurflinearer Regelungssysteme // Zustandsraum, Regelungstechnik und Prozessdaten Verarbeitung. 1972. N. 7. Р. 297-300.
Федосенков Денис Борисович — Россия, 115054, Москва; Сибирская генерирующая компания; канд. техн. наук, доцент; начальник отдела управления активами; [email protected].
Симикова Анна Алексеевна — Россия, 650000, Кемерово; Кузбасский государственный технический университет им. Т. Ф. Горбачева; старший преподаватель кафедры информационных и автоматизированных производственных систем; [email protected].
Федосенков Борис Андреевич — Россия, 650000, Кемерово; Кузбасский государственный технический университет им. Т. Ф. Горбачева; д-р техн. наук, профессор; профессор кафедры информационных и автоматизированных производственных систем; [email protected].
СПИСОК ЛИТЕРА ТУРЫ
Статья поступила в редакцию 22.08.2019
ИНФОРМАЦИЯ ОБ АВТОРАХ
MODAL CONTROL OF FEEDING PROCESSES IN THE MEDIUM OF STATE SPACE AND WAVELET TRANSFORMS
D. B. Fedosenkov1, A. A. Simikova2, B. A. Fedosenkov2
1 Siberian Generating Company, Moscow, Russian Federation
2 T. F. Gorbachev Kuzbass State Technical University, Kemerovo, Russian Federation
Abstract. The article presents the essential aspects of designing and operational functioning of computer-aided modal control systems, using the example of a control system for feeding processes in a mixture-producing set. The object of study is the structure of feeding devices unit and transient processes during feeding. The subjects of research are the issues of designing and creating a computer-aided modal control system (CMCS) which operates in the mode of tunable parametri-cally non-stationary in frequency transient processes (TP). In this case, the scalar signals of the feeding devices (FD) material flows are formed in terms of the state space method (in the format of the internal models "input-state-output"), and monitoring of such processes is carried out in the time-frequency wavelet medium. The latter makes it possible to analyze the processes occurring in the CMCS at a semantically transparent and information-saturated level. The parameters (pulse signal off-duty ratios) of the Fourier model for the signals of a discrete-type feeder are presented, on which the nature of the dosing process depends. A model of the feeders unit is formed in terms of the state space. The ID-flow rate signals and their representations in the form of time-frequency distributions (W-maps) are considered. A new concept was developed for flow rate modal control of the feeding devices unit in the frame of the mixture-producing set by introducing the definition of parametric non-stationary state of the feeding devices internal environment that takes into account the real-time variation of the multidimensional state feedback coefficient during transients in the system and for the purpose of experimental estimation of transients parameters ( here: time constants), they are determined by flow signal distributions on W-maps. The wavelet medium used in the system deals with Cohen's class time-frequency distributions (Wigner-Ville and Choi-Williams transformations). The combination of TP-models in terms of the state space and their representation in the Wigner distribution format allows one to evaluate, set and reassign the parameters of the internal structure of the object controlled (a separate feeder or feeding devices unit) in real time in order to achieve the effect of stabilizing the ingredients flow rate in a rated mode at pre-mixing stage. The transient modes are presented in the form of processes restoring the rated modes (the so-called reverse TP). A procedure has been developed and described for the implementation of a modal control algorithm for the feeding process, including the technology of forced localization and relocalization of the object poles under non-stationary processes in the "Executive mecha-nism-FD" system, using the example of a continuous-type feeder in a transient mode on load decrease (with a downward reverse transitional process). For this mode, the calculation of the non-stationary matrix of the modal controller is presented.
Key words: modal control, transient processes on load increase and decrease, parametrically non-stationary chirp-signals, time-frequency wavelet medium, Wigner distribution, W-map, modal controller matrix.
For citation: Fedosenkov D. B., Simikova A. A., Fedosenkov B. A. Modal control of feeding processes in the medium of state space and wavelet transforms. Vestnik of Astrakhan State Technical University. Series: Management, Computer Science and Informatics. 2019;4:46-58. (In Russ.) DOI: 10.24143/2072-9502-2019-4-46-58.
REFERENCES
1. Fedosenkov B. A. Nauchno-tekhnicheskie osnovy sozdaniya i modelirovaniya avtomatizirovannykh sistem upravleniya nepreryvnymi smeseprigotovitel'nymi protsessami. Avtoreferat dis. ... d-ra tekhn. nauk [Scientific and technical fundamentals for creating and modeling computer-aided control systems for continuous mixing processes. Diss.Abstr... Dr.Tech.Sci.]. Moscow, Izdat. kompleks MGUPP, 2005. 55 p.
2. Kim D. P. Teoriya avtomaticheskogo upravleniya [Theory of automation control]. Moscow, Fizmatlit Publ., 2004. Vol. 2. Mnogomernye, nelineinye, optimal'nye i adaptivnye sistemy. 464 p.
3. Phillips Ch. L., Harbor R. D. Sistemy upravleniya s obratnoi sviaz'yu [Feedback control systems]. Moscow, Laboratoriya bazovykh znanii Publ., 2001. 616 p.
4. Fedosenkov D. B., Simikova A. A., Fedosenkov B. A. Sintez sistemy modal'nogo upravleniya blokom doziruyushchikh ustroistv [Synthesis of modal control system for metering unit]. Mashinostroyenie: inno-vatsionnye aspekty razvitiya: materialy Mezhdunarodnoi nauchno-prakticheskoi konferentsii. Saint-Petersburg, SPBF NITs MS, 2019. No. 2. Pp. 162-163.
5. Mallat S. Veivlety v obrabotke signalov [A wavelet tour of signal processing]. Moscow, Mir Publ., 2005. 673 p.
6. Cohen L. Time-frequency analysis. Englewood Cliffs, Prentice Hall, 1995. 299 p.
7. Auger F., Chassande-Mottin E. Quadratic time-frequency analysis i: Cohen's class. Time-frequency analysis: concepts and methods. ISTE, 2008 (January). Pp. 131-163.
8. Debnath L. Recent development in the Wigner-Ville distribution and time-frequency signal analysis. PINSA, 2002, January, 68, A, no. 1, pp. 35-56.
9. Ackermann J. E. Der Entwurflinearer Regelungssysteme. Zustandsraum, Regelungstechnik und Prozessdaten Verarbeitung, 1972, no.7, pp. 297-300.
Fedosenkov Denis Borisovitch - Russia, 115054, Moscow; Siberian Generating Company; Candidate of Technical Sciences, Assosiate Professor; Chief of the Department of Assets Administering; [email protected].
Simikova Anna Alekseyevna - Russia, 650000, Kemerovo; T. F. Gorbachev Kuzbass State Technical University; Senior Lecturer of the Department of Information and Computer-Aided Manufacturing Systems; [email protected].
Fedosenkov Boris Andreyevitch - Russia, 650000, Kemerovo; T. F. Gorbachev Kuzbass State Technical University; Dr. Sc; Professor; Professor of the Department of Information and Computer-Aided Manufacturing Systems; [email protected].
The article submitted to the editors 22.08.2019
INFORMA TION ABOUT THE A UTHORS
