Система поддержки принятия решений по управлению экструдерами с перенастраиваемой конфигурацией

А.Н. Полосин,

доц. каф. САПРиУ, к.т.н., доц., polosin-1976@mail.ru

Т.Б. Чистякова,

зав. каф. САПРиУ, д.т.н., проф., chistb@mail.ru

СПбГТИ(ТУ), г. Санкт-Петербург

Разработана система поддержки принятия решений (СППР) по управлению аппаратно гибкими экструдерами различных типов в многоассортиментных производствах полимерных материалов. СППР помогает управленческому производственному персоналу решать задачи структурного и параметрического синтеза экструдеров для обеспечения требуемого качества продукции при перенастройке линий на новые типы материалов и производительность. Она включает подсистемы синтеза модульных шнеков и поиска режимных параметров экструдеров, обеспечивающих заданные показатели качества экструдата. СППР позволяет сформировать 3D модель шнека путем компоновки моделей типовых элементов на основе правил проектирования, рассчитать геометрические параметры секций шнека, осуществить структурный синтез и параметрическую настройку математической модели процесса экструзии, на базе которой решается задача управления качеством экструдата.

 

Decision support system (DSS) for hardware flexible extruders control in multi-assortment productions of polymeric materials has been developed. DSS helps operative factory personnel to solve structural and parametric synthesis tasks of various types extruders for required product quality assurance under lines change-over on new types of materials and throughput. It includes subsystems of modular screws synthesis and determination of extruder’s regime parameters ensuring required extrudate quality. DSS allows to create screw 3D model by configuration of standard elements models based on design rules, to calculate geometric parameters of screw sections, to perform structural synthesis and parametrical setting of mathematical model for extrusion based on which extrudate quality control task is solved.

 

Современные производства полимерных материалов (фармацевтических и пищевых упаковочных пленок, теплоизоляционных пеноплит) характеризуются широким ассортиментом продукции, жесткими требованиями к качеству полуфабрикатов, изделий, неполнотой информации о прямых показателях качества и их зависимостью от множества технологических параметров, сложностью протекающих физических процессов, высокой стоимостью брака из-за дорогостоящего, дефицитного сырья, больших энергетических затрат на производство. Ключевой стадией производства, определяющей качество целевой продукции, является экструзия, заключающаяся в подготовке расплава полимера и его предварительном формовании с целью получения экструдата – однородной пластичной массы, подлежащей окончательному формованию в готовое изделие (пленку, лист, плиту). Эта стадия характеризуется множеством вариантов аппаратурного оформления, что обусловлено использованием экструдеров различных типов и их аппаратной гибкостью, связанной с переменной модульной конфигурацией шнеков и формующих головок. Для выполнения широкого спектра потребительских заказов в современных производствах параллельно функционируют, как правило, от 2 до 15 экструдеров. Это, с одной стороны, позволяет перенастраивать производство на различные типы полимерных материалов и производительность, а с другой стороны, усложняет управление качеством продукции, к которому предъявляются жесткие требования. Вследствие неполноты информации об объекте (из всего массива производственных данных анализируются только 5% технологических параметров) операторы вынуждены принимать решения по управлению на основе субъективной визуальной оценки качества экструдата, исходя из собственного производственного опыта и экспериментально подобранного технологического регламента. Это приводит к увеличению доли брака и возвратных отходов, энергетических затрат, времени перенастройки линии на новое задание, и, как следствие, к снижению рентабельности и конкурентоспособности производства. По данным экспертов, вследствие ошибочных решений операторов до 30% дефектов изготовленных полимерных материалов, наименее приемлемыми из которых для заказчиков являются включения нерасплавленного полимера (белые точки), деструкционные полосы и черные точки, возникает именно на стадии подготовки экструдата. Время перенастройки линии при изменении толщины полимерной пленки на 0,1 мм составляет 10 минут, при переходе на другой тип полимера – 30 минут, при переходе на новый цвет – до 3 часов. С ростом времени перенастройки уменьшается производительность. Поэтому актуальна разработка компьютерной СППР, которая помогает управленческому производственному персоналу решать задачи структурного и параметрического синтеза экструдеров с перенастраиваемой конфигурацией для обеспечения качества экструдата (по требуемой степени однородности и допустимой глубине протекания термической деструкции), гарантирующего отсутствие дефектов целевой продукции, при перенастройке линии на новый тип полимерного материала и производительность.

В связи с этим целью работы является разработка проблемно-ориентированной СППР, позволяющей на базе многовариантных моделей описания экструзионных агрегатов (информационных моделей, математических моделей, моделей представления знаний) решить следующие задачи:

·       сформировать 3D модель шнека путем компоновки моделей типовых элементов на основе правил проектирования для выбранной из базы данных (БД) экструдеров марки экструдера и заданной конфигурации шнека;

·       рассчитать геометрические параметры секций сформированного шнека на основе параметров его элементов;

·       осуществить структурный синтез и параметрическую настройку функциональной математической модели (ММ) процесса экструзии для расчета производительности, энергопотребления экструдера, среднего времени пребывания в экструдере, показателей качества экструдата и поиска управляющих воздействий, обеспечивающих требуемое качество при заданном типе материала, производительности и энергопотреблении экструдера;

·       визуализировать результаты моделирования в виде таблиц и 3D графиков зависимостей показателей качества экструдата от управляющих воздействий с отображением значений управляющих воздействий, обеспечивающих требуемое качество, сохранить в БД результатов синтеза виртуальную модель синтезированного шнека и допустимые значения управляющих воздействий на экструдер.   

Материалы на основе высоковязких полимеров (например, полиэтилена низкой плотности, поливинилхлорида) изготавливаются методами плоскощелевой и раздувной экструзии, каландрования. Метод производства Mform определяет аппаратурное оформление стадий подготовки и формования экструдата. При изготовлении пленок, листов, пеноплит методом плоскощелевой экструзии для подготовки экструдата применяются одношнековые и двухшнековые (с встречным или односторонним вращением шнеков) экструдеры, а для формования экструдата в полотно используются плоскощелевые головки. В большинстве производств пленок методом раздувной экструзии применяются одношнековые экструдеры с кольцевыми головками, выполняющие функции подготовки и предварительного формования экструдата в трубчатую заготовку (рукав), подлежащую дальнейшему раздуву и растяжению. При каландровании экструдат подготавливается, как правило, в осциллирующих и двухшнековых (с односторонним вращением шнеков) экструдерах, оснащенных фильерами, а формуется в полотно на каландрах. Основным возмущением процесса являются колебания качества сырья, вызванные подачей в экструдер потока возвратных отходов (например, раздробленной кромки пленки или раздробленных и гранулированных некондиционных пеноплит) с расходом Gwast.

Характерной особенностью экструдеров, позволяющей перенастраивать их на различные типы полимерных материалов, является модульная конфигурация шнеков. Шнек экструдера представляет собой совокупность насаживаемых на общий сердечник взаимозаменяемых элементов нескольких стандартных типов (транспортные, смесительные, с непрерывной нарезкой и др.) [1], которые отличаются геометрическими параметрами (числом заходов нарезки, характером изменения глубины канала по его длине, числом осевых прорезей в нарезке) и, как следствие, обладают различной смесительной и транспортирующей способностью. Объемные геометрические модели типовых элементов шнеков одношнековых и осциллирующих экструдеров приведены на рис. 1.

 

1 – транспортный элемент с непрерывной однозаходной нарезкой и цилиндрическим сердечником; 2 – транспортный элемент с прерываемой одной осевой прорезью однозаходной нарезкой и цилиндрическим сердечником; 3 – смесительный элемент с прерываемой тремя осевыми прорезями двухзаходной нарезкой и цилиндрическим сердечником; 4 – транспортный элемент с непрерывной однозаходной нарезкой и коническим сердечником            

Рис. 1 Геометрические модели типовых элементов шнеков одношнековых и осциллирующих экструдеров

Последовательно соединенные элементы одного типа составляют секцию шнека. Комбинируя секции, можно управлять степенью деформационного воздействия на материал. Секции из элементов типа 1 используются в зонах питания и транспортировки расплава одношнековых экструдеров, секции из элементов типа 2 – в зонах питания и транспортировки расплава осциллирующих экструдеров, секции из элементов типа 3 – в зонах плавления и смешения осциллирующих экструдеров, секции из элементов типа 4 – в зонах плавления одношнековых экструдеров.  

Анализ характеристик экструдеров позволил сформировать информационное описание объекта управления, включающее входные X и варьируемые V параметры, возмущающие воздействия F, выходные параметры Y (рис. 2).

Рис. 2 Информационное описание экструдера как объекта управления

На рис. 2 использованы следующие обозначения: Y0 – задание; Tprod – тип продукции; G0 – производительность линии, кг/с; Emax, E – предельно допустимое и текущее удельное энергопотребление, Дж/кг; Sprod – вид продукции; Traw – вектор типов сырьевых материалов; Tpolym, Tba – тип полимера и вспенивающего агента; Gscr, Gdie – векторы геометрических параметров шнека и головки; Hpolym – вектор параметров свойств полимера; Qext – вектор показателей качества экструдата; g, g min – текущая и минимально допустимая степень смешения; Id, Idmax – текущий и максимально допустимый индекс термодеструкции, %; S – вектор структурных параметров экструдера; U – вектор управляющих воздействий; Textrud – тип экструдера; D, L – диаметр и длина шнека, м; Cscr – конфигурация шнека; Cdieконфигурация головки; Tej – тип j-го элемента шнека; Ne – число элементов; Nh, N – частоты вращения загрузочного шнека и шнека экструдера, об/с; Tbk – температура корпуса в k-й тепловой зоне, °С; nT – число тепловых зон корпуса; RS – множество допустимых для задания Y0 значений структурных параметров; Umin, Umax – регламентные ограничения на управляющие воздействия; G – производительность экструдера, кг/с; t – время пребывания в экструдере, с.

Показателями качества экструдата Qext, характеризующими интенсивность деформационного и теплового воздействия на материал, являются степень смешения g, отражающая материальную однородность экструдата и определяемая накопленной за время пребывания среднемассовой деформацией сдвига, а также индекс деструкции Id, отражающий степень необратимых изменений в материале под действием высокой температуры за время пребывания [2]. Поэтому время пребывания t должно быть достаточным для гомогенизации материала без деструкции.

При перенастройке на новое задание Y0 СППР позволяет решить задачу управления качеством экструдата:

выбрать из БД характеристик экструдеров марку экструдера, определяемую структурными параметрами {Textrud, D, L, Cdie}, и сформировать из БД технологических параметров диапазоны управляющих воздействий на экструдер [Umin; Umax] для получения из сырьевых материалов заданных типов Traw экструдата с производительностью G ³ G0 и энергопотреблением E £ Emax при реализуемом методе производства Mform продукции вида Sprod;

определить по ММ конфигурацию шнека Cscr Î RS и значения управляющих воздействий на экструдер U Î [Umin; Umax], которые обеспечивают получение экструдата требуемого качества g ³ gmin, Id £ Id max, гарантирующего соблюдение регламентных ограничений на число дефектов поверхности производимого полимерного материала.

Разработанная СППР включает подсистемы структурного и параметрического синтеза экструдера, информационное обеспечение (позволяющее настраивать СППР на новое задание Y0), интерфейс оператора экструдера (рис. 3).

Рис. 3 Функциональная структура СППР

При перенастройке на задание Y0 оператор вводит его параметры: выбирает вид продукции Sprod (пленка, лист, пеноплита), типы сырья – полимера Tpolym и пеноагента Tba (если Sprod = «пеноплита»), метод производства Mform (в зависимости от вида продукции и типов сырья СППР выдает рекомендации по выбору подходящего метода производства). Также оператор задает производительность линии G0 и предельно допустимое энергопотребление экструдера Emax. Модуль выбора марки экструдера формирует из БД геометрических характеристик экструдеров на основе правил выбора множество марок экструдеров, каждая из которых определяется совокупностью структурных параметров {Textrud, D, L, Cdie} и является допустимой для переработки данных типов сырья при изготовлении заданного вида продукции выбранным методом с требуемой производительностью и энергопотреблением [3]. Правила выбора экструдеров формируются разработчиком (администратором) СППР путем извлечения и обработки знаний экспертов (технологов производств полимерных материалов) и с помощью соответствующего интерфейса заносятся в базу правил. Они строятся по схеме продукций типа «ЕСЛИ условие, ТО следствие». Условие состоит из утверждений о виде продукции, типах сырья, методе производства, производительности и энергопотреблении экструдера, соединенных связками И (конъюнкция – «Ù») и ИЛИ (дизъюнкция – «Ú»). Следствие состоит из одного или нескольких предложений, которые образуют выдаваемое правилом решение (тип экструдера, диаметр и длину шнека, конфигурацию головки). Следствие принимается при истинности условия. Примеры правил: ЕСЛИ Sprod = «плоская пленка» Ù Tpolym = «поливинилхлорид» Ù Mform = «каландрование» Ù G0 = 3000 кг/ч Î [2500 кг/ч; 4000 кг/ч] Ù E = 1000 кДж/кг < Emax = 1750 кДж/кг, ТО Textrud = «осциллирующий экструдер» Ù D = 0,2 м Ù L = 2,2 м Ù Cdie = «фильера»; ЕСЛИ Sprod = «пеноплита» Ù Tpolym = «полистирол» Ù Mform = «плоскощелевая экструзия» Ù G0 = 450 кг/ч Î [400 кг/ч; 550 кг/ч] Ù E = 700 кДж/кг < Emax = 1000 кДж/кг, ТО Textrud = «двухшнековый экструдер с односторонним вращением шнеков» Ù D = 0,18 м Ù L = 4,5 м Ù Cdie = «плоскощелевая головка».

Оператор задает конфигурацию шнека, формируя последовательность расположения элементов различных типов в шнеке при помощи правил компоновки. Эти правила разрешают присоединять новый (j+1)-й элемент к уже расположенному в шнеке j-му элементу, если его тип Tej+1 принадлежит множеству допустимых типов для сопряжения с jэлементом шнека данного типа экструдера Textrud и выполняется условие сопряжения элементов по глубине их каналов: Hinj+1 = Houtj, где Hinj+1 – глубина канала на входе в (j+1)-й элемент, Houtj – глубина канала на выходе из j-го элемента. Число элементов Ne определяется в зависимости от длины шнека L и длины каждого элемента Lej, зависящей от его типа Tej и диаметра шнека D. Для заданной конфигурации шнека формируется набор значений геометрических параметров его элементов Gelem. Вектор Gelem включает длину Lej, шаг Bj, угол наклона φj, число заходов zfj и толщину витков ej нарезки, глубину Hinj, Houtj, ширину Wj, длину Zj, объем Vej канала, зазоры различных типов (радиальный δ0j, боковой δsj, межвалковый δrj в зависимости от Textrud), число mj и ширину прорезей sj в нарезке элемента. Параметры {Lej, Bj, zfj, ej, Hinj, Houtj, δ0j, δsj, δrj, mj, sj} передаются из БД геометрических характеристик экструдеров и их элементов в зависимости от типа элемента Tej и диаметра шнека D. Параметры канала, в том числе его текущая глубина Hj, рассчитываются по уравнениям аналитической геометрической модели канала, построенной при допущении о малости кривизны канала (плоская модель канала):

, , , , , (1)

где zj – координата по длине канала j-го элемента, 0 £ zj £ Zj.    

Параметры элементов используются для расчета числа Ns и геометрических параметров секций шнека Gsect:

, ,                                                                        (2)

где Zsi, Vsi – длина (м) и объем (м3) канала i-й секции шнека; Nei – число элементов в i-й секции.  

Сформированная конфигурация шнека и геометрические параметры его элементов передаются в модуль сборки и визуализации 3D модели шнека, который осуществляет настройку параметров и компоновку 3D моделей типовых элементов, выбираемых из библиотеки моделей элементов шнеков в зависимости от типов элементов, входящих в шнек [4]. Таким образом, подсистема структурного синтеза формирует и передает в подсистему параметрического синтеза структурные параметры экструдера S, геометрические параметры шнека Gscr (в том числе характеризующие особенности экструдеров разных типов и зависящие от диаметра D: ход шнека осциллирующего экструдера S0, расстояние между осями шнеков двухшнекового экструдера Ascr), параметры формующей щели головки Gdie.

Параметрический синтез экструдера осуществляется на основе функциональной ММ. Для учета аппаратной гибкости экструдеров, основных особенностей протекающих в них процессов и перерабатываемых материалов реализован комплексный подход к моделированию. Он заключается в синтезе статической ММ, используемой для расчета производительности, энергопотребления экструдера и показателей качества экструдата [2, 5], и динамической ММ для оценки времени пребывания в экструдере, характеризующемся сложной структурой потоков [6]. Структура SM = {Ψ, fG, fτ} и коэффициенты K = {Ks, Kd} ММ процесса экструзии настраиваются на конфигурацию {Textrud, Cscr, Cdie} и параметры {X, U} объекта управления. Ядром статической ММ является базовая модель, описывающая плавление и течение материала в канале шнека и представляющая собой систему уравнений в частных производных, построенных на основе законов сохранения физических субстанций и реологии полимерных материалов. Она предназначена для расчета распределений параметров состояния Ps – скоростей циркуляционного и поступательного потоков, давления, вязкости и температуры материала – по ширине W (0 x W), глубине H (0 y H), длине Z (0 z Z) канала. Синтез структуры статической ММ заключается в формировании системы, включающей уравнения базовой модели Ψ(X, U, Ks, Ps, G, x, y, z) = 0 и уравнение G = fG(D, Ne, Γelem, Γdie, Hpolym, N, Ps, G) для расчета производительности экструдера. Параметрическая настройка ММ осуществляется путем формирования значений коэффициентов ММ Ks с использованием библиотеки ММ в зависимости от типа Textrud и конфигурации экструдера {Cscr, Cdie } и типа материала Traw. Метод синтеза динамической ММ является обобщением описанного в литературе подхода [7] и заключается в компоновке уравнений гидродинамических моделей, описывающих течение в типовых секциях шнека, с включением между зонами рециклов, учитывающих осевое перемешивание, вызванное утечками и осцилляцией шнека (в экструдере соответствующего типа). В результате формируется система уравнений τ = fτ(D, Ns, Γsect, Hpolym, Nh, N, Kd, G) для расчета среднего времени пребывания τ и значения коэффициентов ММ Kd (в зависимости от типов экструдера Textrud и материала Traw). Сложность ММ процесса экструзии обусловлена сильной взаимной зависимостью параметров состояния, нелинейностью, наличием рециклов. Кроме того, ММ включает уравнения различных классов. Поэтому в СППР интегрирована библиотека методов решения (явные и неявные конечно-разностные методы, методы теории плоских несимметричных потоков [2]), обеспечивающих сходимость вычислительного процесса при соблюдении требований к точности и экономичности расчета параметров состояния для различных типов и конфигураций экструдеров и типов материалов.

В зависимости от задания Y0 и марки экструдера из БД технологических параметров формируются регламентные диапазоны управляющих воздействий [Umin; Umax], передаваемые в блок вычислительных экспериментов СППР. Решение системы уравнений статической ММ Ψ–fG позволяет рассчитать параметры состояния Ps и производительность G. При найденной производительности по динамической модели определяется время пребывания τ материала. При рассчитанных параметрах состояния и времени пребывания вычисляются показатели качества экструдата – степень смешения γ = fγ(Ne, Γelem, Ps, τ), индекс деструкции Id = fId(Hpolym, Ps, τ) – в зависимости от управляющих воздействий, варьируемых в регламентных диапазонах. Обработка этих зависимостей позволяет определить значения управляющих воздействий, обеспечивающие требуемое качество экструдата g ³ gmin, Id £ Id max при выполнении ограничений на производительность G ³ G0 и энергопотребление E £ Emax экструдера. Если критериальные ограничения не выполняются, оператор формирует новую конфигурацию, и вычислительный эксперимент на вновь синтезированной ММ повторяется. Если исследованы все допустимые для данной марки экструдера конфигурации, а критериальные ограничения не выполнены, выбирается другая марка (из множества допустимых для задания Y0), для которой осуществляется поиск конфигурации шнека и управляющих воздействий, обеспечивающих получение экструдата требуемого качества. Результаты синтеза экструдера отображаются на интерфейсе оператора в виде 3D модели и геометрических параметров синтезированного шнека, таблиц, 3D графиков зависимостей показателей качества экструдата от управляющих воздействий (рис. 4), совета оператору, содержащего значения управляющих воздействий на экструдер синтезированной конфигурации, обеспечивающих требуемое качество экструдата при условии выполнения ограничений на производительность и энергопотребление. Результаты синтеза и параметры задания Y0 сохраняются в БД результатов.

Проверка адекватности ММ экструзии (сравнением рассчитанных и измеренных значений температуры материала, производительности и времени пребывания) и тестирование СППР (сравнением результатов синтеза с данными экспертов-технологов) выполнены на примерах одношнековых экструдеров в производствах рукавных полиэтиленовых пленок, осциллирующих и двухшнековых экструдеров в каландровых производствах поливинилхлоридных пленок, двухшнековых экструдеров в производствах полистирольных пеноплит. Результаты подтвердили адекватность ММ и работоспособность СППР для данного класса объектов управления, а также показали повышение эффективности производства за счет улучшения качества продукции, ресурсо-, энергосбережения, снижения времени перенастройки на новые задания.

Рис. 4 Примеры результатов работы СППР (результатов структурного и функционального моделирования экструдеров)  

Таким образом, разработана гибкая проблемно-ориентированная СППР, которая позволяет на базе синтезируемых ММ процессов экструзии в аппаратно гибких агрегатах решать задачи поиска конфигураций экструдеров и управляющих воздействий на процесс, обеспечивающих получение экструдата требуемого качества при перенастройке производственных линий, реализующих различные методы изготовления продукции из полимеров, на новые типы материалов и производительность. СППР может быть использована в качестве советчика оператора для ресурсо- и энергосберегающего управления экструдерами с перенастраиваемой конфигурацией в многоассортиментных производствах полимерных пленок, листов, пеноплит, а также (при соответствующем расширении БД, баз правил и библиотек ММ) в аналогичных высокотехнологичных производствах, включающих стадию экструзии (например, производства линолеума, сорбентов и катализаторов, керамической плитки).

Литература

1.  Раувендааль К. Экструзия полимеров. СПб. : Профессия, 2008. 768 с.

2.  Полосин А. Н. Математическое моделирование и оптимизация процессов смешения термопластичных материалов в одношнековых экструдерах // Информ. технологии моделирования и упр. 2009. № 2. С. 297–305.

3.  Подсистема поддержки принятия решений по выбору экструдеров при автоматизированном проектировании производств полимерных пленок / А. Н. Полосин [и др.] // Соврем. информ. системы и технологии : материалы I междунар. науч.-практ. конф. Сумы, 2012. С. 194–195.

4.  Защиринский С. В., Иванов А. Б. Полосин А. Н. Подсистема геометрического моделирования осциллирующих экструдеров для многоассортиментных производств полимерных пленок // Мат. методы в технике и технологиях : сб. тр. XXVIII междунар. науч. конф. Саратов, 2015. Т. 10. С. 367–371.

5.  Чистякова Т. Б., Полосин А. Н. Методы и технологии синтеза математических моделей процессов экструзии в гибких производствах полимерных материалов // Вестн. Сарат. гос. техн. ун-та. 2011. № 4. С. 170–180.

6.  Чистякова Т. Б., Полосин А. Н. Система имитационного моделирования для управления экструдерами в гибких производствах полимерных материалов // Автоматизация в пром-сти. 2010. № 7. С. 3540.

7.  Hoppe S., Detrez C., Pla F. Modeling of a cokneader for the manufacturing of composite materials having absorbent properties at ultra-high-frequency waves. Part 1 // Polym. Eng. a. Sci. 2002. Vol. 42, № 4. P. 771–780.