Цель лабораторной работы - изучение процессов тепломассопереноса полимера в зоне плавления червячного пресса. Задачей лабораторной работы является исследование влияния геометрических, технологических и физических факторов на изучаемый процесс с использованием метода математического моделирования
Краткие теоретические сведения
Экструзия - это технологический процесс, сущность которого состоит в придании материалу нужной формы путем продавливания материала через профилирующий инструмент (матрицу).
Процесс экструзии включает в себя четыре элементарных стадии, разделенные по зонам: загрузки, задержки плавления, плавления и дозирования.
В зоне плавления материал переходит из твердого состояния в жидкое под воздействием тепла, подводимого от стенок корпуса, и тепла, выделяющегося в результате вязкого трения. Многочисленные исследования показывают, что процессы, происходящие в зоне плавления, для всех пластифицирующих экструдеров примерно одинаковы.
Полимер в твердом состоянии, называемый пробкой гранул, двигается по винтовому каналу червяка и попадает в участок корпуса, где температура внутренней поверхности корпуса выше температуры плавления материала. Контактирующая со стенкой поверхность полимера начинает плавиться, и на внутренней поверхности корпуса экструдера образуется тонкая пленка расплава. В тот момент, когда толщина этой пленки d (см.рис.6) оказывается больше, чем величина радиального зазора e между гребнем нарезки червяка и стенкой корпуса, толкающая стенка канала червяка начинает соскребать слой расплава с внутренней стенки корпуса и собирать его в области, непосредственно примыкающей к передней грани толкающей стенки. По мере движения пробки ширина области, заполненной расплавом (W‑X), постепенно увеличивается, а ширина твердой пробки X, соответственно, уменьшается. Этот механизм плавления представлен на рис. 1. Процесс плавления заканчивается в тот момент, когда пробка исчезает.
Рис.1. Схематическое представление механизма плавления полимера: 1 - червяк ; 2 - бассейн расплава ; 3- пленка расплава; 4 - корпус экструдера; 5 - твердая пробка полимера; е - зазор между гребнем червяка и корпусом экструдера; d - толщина слоя расплава; w - ширина канала; H - высота канала; x- ширина твердой фазы.
Очевидным следствием такого механизма плавления является зависимость длины зоны плавления от основных параметров технологического режима: производительности и заданного распределения температур.
Естественно, что интенсивность процесса плавления должна быть достаточно высока для того, чтобы процесс плавления успел закончиться и канал червяка оказался заполненным расплавом прежде, чем материал подойдет к формующей головке. Поэтому основным назначением математической модели процесса плавления в экструдере является предсказание количества переплавленного полимера в любом сечении канала экструдера; длины экструдера, необходимой для полного плавления полимера и зависимость этой величины от геометрии экструдера, условий переработки и свойств материала.
Для описания процессов тепломассопереноса в зоне плавления будем использовать принцип обращенного движения. Кроме того, т. к. высота канала много меньше, чем внутренний диаметр шнека, то кривизной канала можно пренебречь. В этом случае, развернув канал, как ленту серпантина, можно рассматривать процесс в длинном прямоугольном канале, дно которого и боковые стенки, образованные поверхностью спиральной нарезки червяка, неподвижны, а верхняя стенка, образованная внутренней поверхностью корпуса, движется с постоянной скоростью Vk, равной окружной скорости вращения червяка и направленной под углом j к продольной оси канала z. Этот угол соответствует углу подъема винтовой линии шнека (рис. 2).
Рис.2. Геометрическая модель канала.
Для математического описания процессов движения и плавления в канале пластицирующего экструдера возможно использовать различные подходы, основанные на тех или иных упрощающих предположениях. Наиболее простой математической моделью процесса плавления является модель З.Тадмора, позволяющая в одномерной постановке определить зависимость длины зоны плавления от различных технологических, геометрических характеристик и свойств перерабатываемого материала.
Допущения модели Тадмора:
1) кривизной канала пренебрегаем;
2) процесс стационарный;
3) расплав является ньютоновской жидкостью;
4) в целом задача одномерная;
5) характеристики материала постоянны;
6) пробка гранул имеет прямоугольную форму;
7) плавление происходит только у внутренней поверхности корпуса;
8) температура пробки гранул изменяется только по высоте канала;
Представление процессов движения и теплообмена полимеров основывается на законах сохранения массы, количества движения и энергии.
Главной задачей исследования является нахождение длины зоны плавления, которая определяется длиной канала, где ширина твердой фазы обращается в нуль.
Обозначим ширину пробки гранул через X, получим зависимость X=X(z), т. е. изменение ширины пробки по длине канала. Выделим из пробки гранул элементарный объем. На рис. 3 представлен элементарный фрагмент пробки и распределение температуры в поперечном сечении канала
Рис. 3. Элементарный объём пробки гранул и температурный профиль пробки гранул
(1)
(2)
Где Vb – окружная скорость; Vbz – компонента окружной скорости в направлении оси z; Vbx – компонента окружной скорости в направлении оси x; Vsz – скорость пробки вдоль оси z; Vsy – скорость пробки вдоль оси y; Tb – температура корпуса; Ts – температура загружаемого материала, n – частота вращения шнека.
Пробка гранул движется с локальной скоростью Vsz, направленной вдоль канала червяка. Скорость твердой пробки относительно поверхности (относительная скорость Vj) цилиндра находится как:
(3)
Введение относительной скорости Vj позволяет упростить задачу (в частности, это касается граничных условий).
Гидродинамика расплава, с учетом указанных ранее допущений, описывается следующими уравнениями движения и граничными условиями:
(4)
(5)
Интегрируя уравнение (4) с учетом (5), получим изменение скорости расплава в пленке в зависимости от координаты y по высоте:
(6)
Рассмотрим уравнение энергии для жидкой фазы в пленке толщиной d:
(7)
Граничные условия:
(8)
Отметим, что все теплофизические характеристики для расплава будем обозначать индексом m (от слова melt – расплав) , а для твердой фазы – s (stone – камень).
Проинтегрируем (7) с учетом (6) и (8):
(9)
Уравнения (6) и (9) описывают процессы тепломассообмена в пленке расплава. Для твердой фазы уравнения энергии имеют следующий вид:
, (10)
. (11)
Решение уравнения (10) с учетом (11) позволяет определить изменение температуры пробки по высоте:
(12)
где
Отсюда видно, что распределение температуры в пробке не зависит от координаты z, следовательно, не изменяется по длине, сохраняя один и тот же вид зависимости.
Используя уравнения (12) и (9), определим величины потоков тепла на границе раздела твердой и жидкой фаз.
(13)
где l – теплота фазового перехода, Дж/кг.
Уравнение (13) связывает толщину слоя расплава d со скоростью Vsy. Для определения последних составим уравнение баланса масс в жидкой фазе.
Пренебрегая малыми величинами, будем иметь:
где w – скорость плавления материала на единичной длине канала, кг/м/с.
Откуда:
(14)
Подставляя уравнение (14) в уравнение (13) и выражая величину d, получим:
(15)
Для скорости плавления имеем выражение:
(16)
Приближаясь к решению задачи, т. е. к определению зависимости X=X(z), запишем уравнение баланса фаз, теперь для твердой фазы на участке dz и, переходя к пределу, получим уравнение:
(17)
Где выражение для w(X), полученное на основании вышеизложенных выкладок, определяется уравнением (16).
Проинтегрируем уравнение (17) с учетом условия X(0)=w, для канала постоянной глубины будем иметь:
(18)
Из выражения (18), положив X=0, определяем длину зоны плавления:
(19)
Для канала переменной глубины, высота которого изменяется по закону : H=Hl-A×z , где H(0)=Hl, A- угол конусности червяка, решение выглядит так:
(20)
Длина зоны плавления в этом случае определится:
(21)
Теплофизические и реологические характеристики некоторых полимеров смотри в приложении.
Задание
В соответствии со своим вариантом задания выполнить следующие расчеты.
1) Для заданного номинального технологического режима:
– разработать алгоритм и расчетную программу;
– рассчитать компоненты скоростей и относительную скорость;
– определить длину зоны плавления, определить Ф и Y; определить длину зоны плавления для канала постоянной глубины, для канала переменной глубины;
– рассчитать ширину твердой пробки в зависимости от длины зоны плавления или для канала постоянной или переменной глубины соответственно;
– определить распределение скорости плавления по длине;
– рассчитать распределение температуры по высоте канала в пяти различных сечениях в зоне плавления для жидкой и твердой фаз соответственно.
2) Исследовать влияние на процесс плавления полимера следующих факторов:
– температуры корпуса;
– начальной температуры материала;
– частоты вращения шнека;
– расхода материала;
– физико-механических свойств материала;
– угла конусности винтовой нарезки червяка.
3) Провести анализ полученных закономерностей процессов тепломассопереноса полимера в зоне плавления экструдера.
4) Построить графики полученных зависимостей.
5) Оформить отчет.
Исходные данные
Свойства материалов
№
λm/ λsДж/м/с/°С
ρm / ρsкг/м3
cm / csДж/кг/°С
μПа с
Ts/Tm°С
0.24/0.28
800 / 884
2000 / 2030
15/110
0.18/0.22
870 / 921
1940 / 2000
20/110
0.27/0.28
930 / 1019
2010/ 2060
25/110
0.25/0.30
760 / 835
2000/ 2070
16/110
0.20/0.25
810 / 870
1990/2040
17/110
0.24/0.28
750 / 850
2020/2080
18/110
0.30/0.33
850 / 950
2040/2120
19/110
0.24/0.28
790 / 890
2100/2200
21/110
0.28/0.32
750 / 850
2150/2300
22/110
0.24/0.30
800 / 900
1900/2000
23/110
0.18/0.24
850 / 960
1950/2030
24/110
0.22/0.27
720 / 820
1800/1900
25/110
0.27/0.31
730 / 830
2030/2100
10/110
0.25/0.31
880 / 980
1940/2010
11/110
0.20/0.27
710 / 807
2010/2080
12/110
0.24/0.29
820 / 920
2000/2090
13/110
0.30/0.35
740 / 830
1990/2060
14/110
0.28/0.32
810 / 855
2020/2100
15/110
0.34/0.36
830 / 930
2040/2110
16/110
0.20/0.27
700 / 800
2100/2200
17/110
0.25/0.29
880 / 980
2150/2210
18/110
0.27/0.34
810 / 910
1900/1990
19/110
0.22/0.31
820 / 920
1950/2050
20/110
0.28/0.35
780 / 880
1800/1950
21/110
0.30/0.37
810 / 900
2100/2200
22/110
λ=130000 Дж/кг
Геометрические и технологические характеристики экструдеров
Номер варианта
Диаметр шнека, Db м
Угол нарезки, гр.
Шаг нарезки, м
Ширина гребня, м
Высота канала в з. п., м
Скорость вращения, об/мин
Расход мате-риала, кг/с
Темпе-ратура в з. п., °С
0.09
17.67
0.09
0.008
0.0055
0.0185
220-250
0.12
17.67
0.12
0.011
0.005
0.0457
220-250
0.15
17.95
0.15
0.015
0.0047
0.07041
220-250
0.16
17.67
0.16
0.017
0.0179-0.0033
0.0786
220-250
0.09
17.67
0.09
0.009
0.015-0.007
0.01944
220-250
0.0381
17.67
0.04
0.00635
0.0061-0.002
0.002707
184-204
0.1
17.67
0.1
0.01
0.01325-0.007
0.02142
220-250
На рис.4-рис.6 приведены характерные зависимости теплофизических характеристик полимеров от температуры.
Рис. 4. Зависимость теплопроводности полимеров от температуры
Рис. 5. Зависимость удельной теплоемкости полимеров от температуры
Рис. 6. Зависимость плотности полимеров от температуры
3.3 Контрольные вопросы
Что такое экструзия?
Какие допущения для расчета по модели Тадмора принимаются?
Для чего вводят скорость Vj и каков ее смысл?
Каков механизм плавления полимера в червячном прессе?