Text (PDF):
Read
Download
Применение в промышленности непрерывно действующих смесительных аппаратов не получило широкого распространения в силу нерешенности ряда вопросов. В частности, недостаточно полно изучено влияние режимных параметров блока дозирующих устройств на характер и структуру выходных материальных потоков, формируемых на предсмесительной стадии, от которых в конечном счете зависит качество готовой смеси. Кроме того, не установлены и не сформулированы технологические причины возникновения нестабильности процесса дозирования вследствие такого влияния и, как следствие, не разработаны условия и способы введения принудительной параметрической нестационарности режимов работы блока дозаторов с целью динамической стабилизации материалопотоков при смене входных воздействий. Ликвидировать такого рода недостатки можно с помощью модального управления, дающего возможность (путем перестройки параметров замкнутой системы с обратной связью по состоянию) поддерживать выходные потоки блока дозаторов на требуемом номинальном уровне. С учетом реальных условий возникает необходимость в разработке и адаптации альтернативного метода для моделирования, исследования и усовершенствования разнохарактерных нестационарных по частоте процессов дозирования, который, будучи реализован на компьютерной основе, позволил бы автоматизировать процесс управления режимами работы в целях интенсификации смесеприготовления. Таким методом, дающим возможность проводить высокоточный анализ разнохарактерных процессов, является метод всплескового (вейвлет) преобразования на основе адаптивной аппроксимации исследуемых сигналов базовыми локальными волновыми составляющими (вейвлет-функциями). Объектом исследования является структура блока дозирующих устройств и процессы дозирования в смесеприготовительном агрегате непрерывного действия. Предметная область исследования охватывает вопросы, связанные с проектированием и созданием автоматизированной системы модального управления процессами дозирования, функционирующей в режиме перестраиваемой параметрической нестационарности и в условиях мониторирования этих процессов в вейвлет-среде. Цель работы - разработка и исследование аппаратной и программной основы комплекса для моделирования и текущего мониторинга процессов дозирования в составе смесеприготовительного агрегата на основе вейвлет-преобразований, разработка системы автоматизации блока дозирования на базе методов модального управления переходными процессами дозирования. Научная новизна работы заключается: - в формировании математических моделей объекта управления (блока дозирующих устройств - БДУ) высокого порядка, включающего дозирующие устройства (ДУ) непрерывного и дискретного действия - в терминах метода пространства состояний; - создании и реализации концепции модального регулирования расхода БДУ на базе теории пространства состояний в сочетании с аппаратом вейвлет-преобразований; - разработке процедуры многомерно-точечного отображения одномерных материалопотоковых сигналов расхода время-частотными вейвлет-распределениями класса Коэна; - разработке алгоритма динамического пересчета параметров модального регулятора с учетом идентификации режимов дозирования по вейвлет-картам при спорадическом проявлении параметрической нестационарности. Математические модели сигналов дозирующих устройств непрерывного и дискретного действия В работе представлены математические модели сигналов расхода дозирующих устройств непрерывного и дискретного действия. Дозаторы непрерывного действия (НДУ) - спиральные и шнековые - описываются моделями третьего порядка. На основе опытных данных, полученных в ходе исследований порционных дозаторов разных типов (шлюзового, с возвратно-поступательным механизмом, тарельчатого, карусельно-ста-канчатого, а также микродозаторов), в качестве общего подхода был выбран способ точной формализации [1] сигнала порционного расхода. Сигнал порционного, или дискретного, дозирующего устройства (ДДУ), сформированный с учетом такого подхода, представлен на рис. 1. Рис. 1. Параметризация сигнала порционного дозирования для общего режима При таких условиях поток на выходе из порционного дозатора на протяжении произвольного количества циклов описывается функциональной зависимостью (1) где Td - период дозирования; Xmd - весовой расход материала через дозатор; k - количество формируемых доз (циклов); l - произвольный номер цикла дозирования; l, m, n - соответственно значения скважностей: порционного дозирования, интервала формирования дозы до начала отсечки и интервала достижения режима номинального дозирования. Разложение дозирующего сигнала (рис. 1) в тригонометрический ряд Фурье с произвольным периодом имеет вид: где А0, Ak, Bk - коэффициенты Фурье-модели; t - текущее время; Xd(t) - сигнал расхода порционного дозатора. Определены коэффициенты Фурье-разложения как функции режимных параметров дозирования (l, m, n, Td, Xmd) для наиболее типичных режимов работы порционных дозаторов. Итоговая расчетная модель блока дозаторов (рис. 2), состоящая из спирального и шнекового ДУ и трех ДУ порционного типа (ПДУ), рассматривается как векторная динамическая система в пространстве состояний. Рис. 2. Блок ДУ в виде структурного фрагмента схемы моделирования смесеприготовительного агрегата в технологическом пространстве состояний При моделировании процессов в блоке ДУ в пространстве состояний [2, 3] на основе операторных функций звеньев составляются система дифференциальных уравнений, которая затем преобразуется в систему дифференциальных уравнений состояния, записанных в нормальной форме Коши, и линейное алгебраическое уравнение выхода, выражающее выходной сигнал системы через комбинацию переменных состояния и сигнала управления. Представляя набор переменных в виде векторов, а набор параметров в виде соответствующих матриц, получаем векторно-матричную модель (ВММ) агрегата в пространстве состояний. Блок дозирования на структурном уровне представлен системой SISO (англ. single input - single output). Для численной обработки сигналов и реализации функций управления использовались персональные компьютеры ППЭВМ IBM-совместимого типа. В качестве базового программного обеспечения были приняты универсальные интегрированные пакеты: Mathсad 7.0 Professional, Maple V Release 4.0, MathLab 5.2, а также инструментальный пакет Classic 3.0. Итоговая модель дозирующего сигнала спирального дозатора (рис. 2) в терминах пространства состояний [4] имеет следующий вид: Векторно-матричная модель дозирующего сигнала шнекового дозатора соответственно выглядит так: Дозирующий сигнал i-го порционного ДУ, представленный Фурье-моделью, соответствующей общему режиму загрузки смесителя (см. рис. 1), характеризуется системой дифференциальных уравнений в форме Коши, которая в пространстве состояний примет вид где l - нормированный индекс переменной состояния; q - переопределенный номер полной гармоники. По ВММ сформирована полная расчетная структурная схема БДУ в соответствии с рис. 2. Получаемые первичные осциллограммы сигналов зашумлены и нестационарны по частоте, следовательно, в любой момент времени обладают случайно меняющимся мгновенным дина-мическим спектром. Отобразить такие спектры на визуальном уровне существующими тради-ционными способами не представляется возможным. В связи с этим в работе реализован новый подход, в основе которого лежат математические модели и алгоритмы аппроксимации, идентификации и коррекции контролируемых переменных в определенном пространстве. Последнее осуществляется на платформе вейвлет-функций и сформированных на их основе время-частотных вейвлет-тезаурусов Габора [5]. При этом используемая информация о состоянии объекта отображается в вейвлет-среде в формате двумерных/трехмерных квадратичных распределений класса Коэна [6, 7]. В качестве рабочих использовались распределения Вигнера и Чуи - Уилльямса [5, 7, 8]. Суть процедуры управляющего мониторинга сводится к следующему. При изменении структуры спектра временного вектора сигнала изменяется время-частотная локализация соответствующих элементов (атомов) на карте Вигнера [8], что должно быть визуализировано и отслежено мониторинговым комплексом, после чего формируется управляющее воздействие на электропривод дозатора. Учитывая, что в БДУ отдельные дозаторы спорадически меняют свои выходные сигналы (вследствие варьирования входных воздействий за счет спонтанного изменения физико-меха-нических характеристик поступающих в ДУ ингредиентов), представляется целесообразным наблюдать за общим выходным сигналом расхода БДУ, на вейвлет-карте которого запечатлены все сигналы отдельных ДУ. При этом на мониторе формируются «номинальные» фреймы (специфические окна-рамки - рис. 3), на которые проецируются время-частотные атомы вейвлет-распределения (например, распределения Вигнера), соответствующие номинальному режиму дозирования. Рис. 3. Осциллограммы и их время-частотные карты, соответствующие номинальному (а) и текущему (возмущенному) (б) режимам работы дозатора порционного типа Для НДУ спирального и шнекового типа такие распределения со стабильной частотой представляют собой тонкие по частоте (по оси ординат) и вытянутые по времени (по оси абсцисс) время-частотные атомы (ВЧА) карты Вигнера. Поэтому в реальных условиях работы дозаторов на вейвлет-карте (W-карте) происходит спорадический выход (девиация) того или иного ВЧА (для НДУ) или группы атомов (для ДДУ) за границы фрейма/фреймов. При этом автоматизированной системой фиксируется очередность нарушения номинальных режимов тех или иных ДУ. Поэтому время обработки W-карты соответствующего ДУ должно быть меньше промежутка времени между актами возникновения двух смежных по времени девиаций ВЧА за границы фреймов двух каких-либо разных ДУ либо того же дозатора. Экспериментально установлено, что такие временные промежутки (t-гэпы) обычно составляют не менее 20 периодов для НДУ (при T = 0,5-1,0 с это составит 10-20 с) и 15 периодов для ДДУ (при Tmax = 1,0 с, т. е. fmin = 1 Гц, это равно 15 с; для 10-й гармоники имеем 10 Гц, при этом t-гэп составит 1,5 с). Подобная оценка распределений случайной величины «t-гэп» для обоих типов ДУ произведена в терминах отображения квартильных распределений в формате «Box&Whiskers». Следовательно, тактовая частота CPU современного PC, используемого в автоматизированной системе для обработки сигналов, может быть сколь угодно малой. При фиксации недопустимой девиации ВЧА за границы фреймов начинает действовать технология модального управления. Управление процессами дозирования в автоматизированной системе производится в вейвлет-среде на основе алгоритма модального управления. Цель алгоритма - стабилизировать процесс дозирования на предсмесительной стадии на номинальном уровне при флуктуациях выходных потоков вследствие спорадических ступенчатых или частотно-время-зависимых изменений входных воздействий на дозаторы. Суть проектирования системы модального управления заключается в синтезе замкнутой системы автоматического управления (ЗСАУ) с обратной связью (ОС) по полному вектору состояния путем задаваемого размещения полюсов объекта с последующим определением матрицы ОС (модального регулятора), формирующей задаваемую локализацию созвездия полюсов. Задача синтеза заключается в определении желаемого положения корней характеристического уравнения системы и нахождении коэффициентов Ki (параметров модального регулятора), обеспечивающих заданное размещение корней (см. раздел «Алгоритм пересчета матрицы (параметров) модального регулятора К»). Модель блока ДУ, состоящего из двух дозаторов - спирального и порционного дискретного (последний формирует материалопоток, задаваемый функцией расхода П-образного вида со скважностью λ = T / τ = 2, где T - период дозирования; τ - длительность формирования дозы), - задается соответствующими импульсными переходными функциями. В редуцированном M-графе (сигнальном графе Мэйсона) [1, 3] отражены все составляющие изображений сигналов, что позволяет промоделировать полученную структуру в любой инструментальной системе (здесь моделирование проведено в пакете Classic). Результаты моделирования и исследования реальных режимов работы автоматизированной системы модального управления процессом дозирования Так как бункерные емкости каждого ДУ заполнены материалами исходных ингредиентов, обладающими в каждой емкости различными объемными физико-механическими свойствами (например, влажностными характеристиками и/или степенью однородности по своему составу), в процессе дозирования происходит спорадическое варьирование воздействия материала дозируемого ингредиента на входную (заборную) часть активного органа соответствующего дозатора. Вследствие этого происходит смена частоты дозирующего сигнала (т. е. частоты вращения исполнительного механизма (ИМ) дозатора ДУ) в переходном режиме, что в итоге приводит процесс дозирования к новому установившемуся частотному режиму, который не соответствует номинальному, требуемому на предсмесительной стадии. При поступлении более однородного и сухого (менее увлажненного) вещества ИМ ДУ раскручивается благодаря сбросу нагрузки на рабочий орган дозатора, а частота сигнала материалопотока при этом возрастает до тех пор, пока не стабилизируется на новом установившемся уровне с новой, неноминальной, производительностью. При повышении же нагрузки на рабочий орган ДУ картина противоположна. В том и другом случаях требуется воздействовать на замкнутую САУ (ЗСАУ) с модальным управлением с целью стабилизации материального потока на номинальном уровне. Но поскольку воздействие на входе и реакция на выходе ДУ уже стабилизировались на новом уровне, требуется воздействовать на внутреннюю структуру ЗСАУ таким образом, чтобы изменить параметрические свойства системы. Таким образом, требуется создать эффект принудительной параметрической нестационарности за счет смены локализации «созвездия» полюсов ЗСАУ. Следовательно, при раскрутке ИМ ДУ (т. е. при повышении частоты сигнала в результате L-сброса на рабочий орган ДУ) необходимо - с целью стабилизации режима расхода на номинальном уровне - снизить общий коэффициент передачи ЗСАУ, а также перманентно уменьшать значение мнимой сопряженной пары полюсов (для НДУ) в ходе переходного процесса (ПП). Что касается обратного ПП дозирования в ПДУ, то снижать значение мнимых полюсов следует только для наиболее мощных субгармоник первой гармоники Фурье-модели. Модель при этом имеет вид B1sin ω1t + A1cos ω1t. Режекция высших гармоник выполняется средствами низкочастотной фильтрации или с помощью высокочастотного режекторного фильтра. При реализации обратного ПП после L-наброса средствами модального управления восстановление номинального режима дозирования производится по двухэкспоненциальной зависимости вида x(t)b/ei = x(t)0b/ei + x(t)mb/ei sin[2πfvar(t - τb) - φ0], где fvar = fb [1(t - τb)] + (fe - fb) ([1(t - τe)] - exp[- (t - τb) / T3] T3 / (T3 - T4) + exp[ - (t - τb) / T4] T4 / (T3 - T4)) - изменение частоты чирп-сигнала [5] восходящего обратного ПП; индексы b и e обозначают начало и конец переходного процесса; i - номер дозатора; x(t)mb/ei = xmbi - t / T; τb ≤ t ≥ τe - изменение амплитуды переменной составляющей ВОПП; x(t)0b/ei = x0bi [1(t - τb)] + + (x0ei - x0bi) ([1(t - τe)] - exp[- (t - τb) / T3] T3 / (T3 - T4) + exp[- (t - τb) / T4] × T4 / (T3 - T4)) - изменение апериодической составляющей восходящего обратного ПП. Таким образом, полный график переходных процессов при двух циклах варьирования нагрузки (L-наброса и L-сброса), включающий нисходящий прямой/восходящий обратный переходные процессы (последний восстанавливает номинальный режим дозирования), а также восходящий прямой ПП (ВППП) и нисходящий обратный ПП (НОПП), имеет вид (рис. 4): Рис. 4. Циклы модального управления процессом дозирования при L-набросе и L-сбросе: два цикла ПП (а); ВОПП - сигнал расхода в ходе восходящего обратного ПП (б); НОПП - сигнал расхода при нисходящем обратном ПП (в); НППП - нисходящий прямой ПП На вейвлет-карте восходящий обратный переходный процесс (ВОПП) для ДУ непрерывного действия после L-наброса отображается в виде распределения Вигнера (рис. 5). Рис. 5. Время-частотное отображение восходящего обратного переходного процесса, восстанавливающего требуемый режим дозирования (частота дозирования f(t)var = 3,24÷6,13 Гц) Для демонстрации разрешающей способности W-карт на рис. 6 представлены осциллограмма регистрируемого материалопотокового сигнала и ее W-карта (карта Вигнера) для БДУ, состоящего из одного порционного, двух спиральных, работающих непрерывно, с рабочими частотами 0,20; 3,23 и 4,02 Гц, и спирального дозатора, функционирующего дискретно во времени с заполняющей в пределах интервала дозы частотой 6,89 Гц и частотой дозирования 0,077 Гц. Рис. 6. Материалопотоковый сигнал БДУ и его отображение на W-карте Вигнера На рис. 7, а показана осциллограмма сигнала расхода ПДУ в виде чирп-сигнала с время-за-висимой мгновенной падающей частотой в режиме наброса (увеличения) нагрузки, который в аналитическом виде может быть записан следующим образом: x(t) = F(A0, Ak, Bk, ωk), где ; ; ; . Рис. 7. Динамика сигналов порционного дозирующего устройства в режиме наброса нагрузки в одномерном и многомерном виде: исходный сигнал расхода (а); сигнал, восстановленный по спектральным составляющим в виде вейвлет-ряда (б); время-частотное распределение сигнала (в) После разложения сигнала по алгоритму вейвлет-поиска соответствия [1, 5] получаем его время-частотное распределение (рис. 7, в), адекватность которого установлена путем сравнения исходного сигнала (рис. 7, а) с восстановленным по спектральным составляющим (рис. 7, б). При этом ошибка аппроксимации сигнала расхода в максимуме и на паузе составляет менее 4,3 %, а в точках Гиббса (в точках разрыва первого рода) - 8,9-9,2 %. Рис. 8 демонстрирует материалопотоковый сигнал двух дозаторов непрерывного действия, работающих в режимах наброса и сброса нагрузки, реконструированный по ВПС-алгоритму сигнал, ошибку аппроксимации и карту Вигнера данного процесса. Рис. 8. Переходные режимы работы двух дозаторов непрерывного действия при набросе и сбросе нагрузки: одномерный суммарный сигнал расхода (а); его аппроксимация по ВПС-алгоритму (б); ошибка аппроксимации (в); многомерные распределения на W-карте (г) На карте Вигнера этого режима ниспадающая цепочка время-частотных атомов соответствует прямому переходному процессу (НППП) при набросе нагрузки в одном дозаторе, а восходящая цепочка - прямому переходному процессу (ВППП) в другом. Следует отметить, что все переходные процессы (ПП) в объектах (ДУ) являются параметрически нестационарными, т. е. протекают с время-зависимыми мгновенными частотами в виде так называемых чирп-сигналов. При этом нестационарность определяется частотой fvar ПП и меняющейся амплитудой Xmvar переменной составляющей процесса. Первый параметр меняется под действием непрерывно варьируемых полюсов Sp2 и Sp3 (мнимых сопряженных), а второй параметр Xmvar - в зависимости от соотношения fvar и постоянных времени в произвольный момент. Алгоритм пересчета матрицы (параметров) модального регулятора К Процедура реализации алгоритма модального управления процессом дозирования, в том числе технология принудительной локализации и релокализации полюсов разомкнутой САУ в условиях нестационарных процессов в системе «исполнительный механизм - ДУ» (ИМ - ДУ), рассмотрена на примере НДУ в возмущенном режиме при сбросе нагрузки (при переходном процессе типа НОПП). Подобную процедуру целесообразно применять и в многомерно-точечных системах, использующих вейвлет-среду как средство эффективного отображения одномерных сигналов. Расчет параметров модального регулятора, т. е. определение матрицы обратной связи K по полному вектору переменных состояния, сводится к ее расчету по формуле Аккермана [9] где Qc-1 - обратная матрица управляемости: α(A) - матричный полином, сформированный согласно теореме Кэйли-Хэмилтона [2, 3]. Желаемый характеристический полином определяется по формуле α(s) = s4 - d1s3 + d2s2 - - d3s + d4, где d1 = (sp1 + sp2 + sp3 + sp4), где spj - полюсы объекта; d2 = [sp1 sp2 sp3 x (sp1 + sp2) + sp4 × ×(sp1 + sp2 + sp3 + sp4)]; d3 = [sp1 sp2 sp3 + sp1 sp2 sp3 sp4 (sp1 + sp2)]; d4 = sp1 sp2 sp3 sp4. Следовательно, матричный полином имеет вид α(A) = A4 - d1 A3 + d2 A2 - d3 A + d4 I, где I - единичная матрица. По виду α(A) нетрудно заметить, что в системе РСАУ влияние на переходный процесс имеют полюсы sp2, sp3 и sp4 (два мнимых сопряженных и вещественный отрицательный), соответствующие параметрам ω и T, которые создают частотную составляющую и определяют быстродействие переходного процесса. В результате получаем матрицу K модального регулятора вида K = [0 0 0 1] Qc-1 α(A) = [0 |(T-1ω2 - 1) |(ω2 - 1) |(T-1 - 1)] = [K1 |K2 |K3 |K4]. Данная матрица непрерывно пересчитывается в реальном масштабе времени с произвольным шагом миллисекундного диапазона в ходе переходного процесса, чтобы обеспечить его время-частотный (чирп-) характер. По аналогичным процедурам автоматизированная система рассчитывает параметры модального регулятора для любого ДУ или блока дозаторов. Полученные параметры позволяют задавать требуемые нестационарные полюсы замкнутой САУ с полной обратной связью по состоянию, которые определяет автоматизированная система модального управления по разработанному алгоритму. Таким образом, разработанная автоматизированная система модального управления процессами дозирования выполняет функции стабилизации номинального режима работы БДУ на предсмесительной стадии, что способствует получению смесей высокого качества. Заключение 1. Разработаны математические модели процессов дозирования в терминах технологического пространства состояний блока дозирующих устройств, с помощью которых можно эффективно отрабатывать широкий спектр текущих режимов мультиингредиентного дозирования устройствами непрерывного и дискретного действия. 2. Сформирована и отработана концепция способа модального управления расходом блока дозирующих устройств (объекта управления) на базе теории пространства состояний в сочетании с аппаратом вейвлет-преобразований, реализованная в замкнутой системе регулирования на основе принудительной параметрической нестационарности обратной связи по состоянию. В условиях спорадического нарушения номинальных потоков дозирования, возникающего вследствие изменения входных воздействий, это позволяет варьировать вектор внутренних параметров объекта, приводя выходные потоки последнего к требуемому состоянию. 3. Представлен алгоритм динамического пересчета параметров многомерной обратной связи по состоянию объекта управления в ходе процесса дозирования (в режиме реального времени), обеспечивающий стабильность материалопотока на предсмесительной стадии. 4. В качестве базовых - при проведении моделирования, расчетов и при обработке сигналов и их отображений - использованы вейвлет-функции Габора, адекватно соответствующие характеру технологических сигналов в подсистеме дозирования и в системе смесеприготовления в целом. 5. Испытания разработанной автоматизированной системы модального управления процессом дозирования показали, что она позволяет интенсифицировать процесс смесеприготовления за счет гармонизации материалопотоков на предсмесительной стадии на базе модального управления в пространстве состояний и вейвлет-преобразований, а также повысить точность дозирования по сравнению с существующими производственными системами. Последний фактор ведет к снижению коэффициента неоднородности конечных смесей, т. е. к повышению их качества.