Abstract

Abstract

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

Full Text

Переработка полимерных материалов и создание необходимого для этого оборудования представляют собой одно из важнейших направлений современной промышленности. На сегодняшний день около 60 % мирового производства пластических масс перерабатываются методом экструзии, основным преимуществом которого являются непрерывность и возможность его совмещения с другими технологическими процессами. Наиболее широкое применение среди экструзионных машин при переработке полимеров нашли одночервячные пластифицирующие экструдеры. В зоне дозирования пластицирующего экструдера существует только жидкая фаза. При построении математической модели процессов течения и теплообмена расплава полимера в зоне дозирования винтового канала экструдера (рис. 1) вводятся следующие упрощающие предположения: процесс стационарный при постоянном массовом расходе; винтовой канал разворачивается на плоскость и используется принцип обращенного движения (рис. 2); массовыми силами по сравнению с силами вязкого трения пренебрегаем. Пренебрегаем утечками жидкости через зазор. Таким образом, процесс движения и теплообмена полимера в винтовом канале экструдера моделируется тепломассопереносом материала в длинном прямоугольном канале, верхняя стенка которого движется с постоянной скоростью , равной окружной скорости червяка, под углом нарезки винтовой линии к оси канала. Рис. 1. Схема винтового канала экструдера Рис. 2. Схема развернутого винтового канала С учетом указанных допущений система дифференциальных уравнений в декартовой системе координат, описывающая движение и теплообмен полимера в канале червяка и полученная на основании законов сохранения массы, количества движения и энергии, будет иметь вид [1–3]: (1) (2) (3) (4) (5) где , , – прямоугольные координаты; – компоненты вектора скорости; tij – компоненты тензора девиатора напряжений; F – диссипативный источник тепла; – градиенты давления. Связь тензора напряжений и тензора скоростей деформаций определяется реологическим уравнением состояния. Компоненты тензора напряжений для несжимаемой жидкости в декартовой системе координат имеют вид: (6) где – эффективная вязкость, являющаяся функцией скорости сдвига и температуры, (7) – температурный коэффициент вязкости; – показатель аномалии вязкости; – коэффициент консистенции расплава полимера при ; – квадратичный инвариант тензора скоростей деформаций, (8) Функция диссипации вычисляется по формуле (9) Искомое решение системы (1)–(9) в виде полей скорости и температуры определяется при заданной геометрии исследуемой области и совокупности краевых условий. На входе в канал задается известное распределение температуры. На внутренней поверхности корпуса и на шнеке задается распределение температуры, определяемое технологическими условиями переработки полимерного материала. Граничные условия для составляющих скоростей определяются из условия прилипания жидкости к твердым непроницаемым поверхностям (стенкам канала и поверхности раздела фаз) (см. рис. 2). Расчеты проводятся в режиме заданного расхода. Задача решалась численным методом конечных элементов с помощью программного комплекса ANSYS FLUENT. Разбивка на конечные элементы была произведена в программе ICEM CFD. Геометрические размеры канала приведены в табл. 1. Таблица 1 Геометрия шнека № п/п Геометрические размеры Значение 1 Диаметр шнека, мм 160,0 2 Шаг винтовой нарезки, мм 160,0 3 Ширина канала W, мм 137,3 4 Ширина гребня винтовой нарезки S, мм 15,3 5 Суммарная длина шнека, витки 4,7 6 Глубина канала, H, мм 15,25 7 Угол подъема винтовой линии Q 17°39¢ Реологические и теплофизические свойства исследуемого расплава полиэтилена приведены в табл. 2. Таблица 2 Реологические и теплофизические свойства n m0 T0 b r l С – Па×сn °С 1/°С кг/м3 Вт/(м·°С) Дж/(кг·°С) 0,44 10825 160 0,018 779,0 0,182 2400 Температура расплава полиэтилена на входе в канал равна 180 °С, на границах канала – 200 °С. Скорость вращения шнека равна 60 об/мин. Исследования проведены для переменной производительности экструдера , которая изменялась в диапазоне от 0,05 до 0,2 кг/с. На рис. 3 приведена напорно-расходная характеристика экструдера. Поле компоненты скорости складывается из напорной и вынужденной составляющих расходов. Увеличение производительности экструдера при неизменной вынужденной составляющей расхода приводит к росту напорной составляющей расхода и к падению градиента давления, в связи с чем давление на выходе снижается. Подпись: Давление на выходе, МПа Рис. 3. Зависимость изменения давления на выходе из канала при изменении расхода Изменение температуры расплава полиэтилена на выходе показано на рис. 4. Из рисунка видно, что температура материала на выходе превышает температуру, заданную на границе канала. Это связано с тем, что большая роль в разогреве материала связана с диссипацией тепла в результате работы сил вязкого трения. При увеличении производительности снижается среднее время пребывания материала в канале, в результате чего температуре на выходе падает. Подпись: Температура на выходе, С Рис. 4. Зависимость изменения температуры на выходе при изменении расхода Таким образом, представленная трехмерная математическая модель может быть использована для исследования процессов тепло- и массопереноса в винтовых каналах зоны дозирования пластицирующих экструдеров, шнековых насосов. Би