CC BY
57
УДК 621.315
А.Г. Щербинин, Е.В. Субботин, А.Е. Терлыч
Пермский национальный исследовательский политехнический университет
ИССЛЕДОВАНИЕ ВЛИЯНИЯ ПРОЦЕССОВ ТЕПЛОПРОВОДНОСТИ ШНЕКА НА ТЕПЛОВОЙ БАЛАНС ЭКСТРУДЕРА
Предлагается математическая модель по определению температуры шнека. Проведено исследование влияния процессов теплопроводности в шнеке на перераспределение энергетических и температурных характеристик одночервячного пластицирую-щего экструдера.
В современном кабельном производстве при изготовлении пластмассовой изоляции широкое применение находят одночервячные пластицирующие экструдеры. Качество изделий, получаемых на экструзионном оборудовании, во многом зависит от температурного режима переработки. Одним из регулируемых технологических параметров является температура корпуса. Температура шнека устанавливается естественным образом в процессе теплообмена и зависит от заданной температуры корпуса, числа оборотов шнека, расхода материала и его теплофизических и реологических свойств. При этом в процессе работы экструдера о распределении температуры в шнеке ничего неизвестно. Некоторые шнеки имеют продольное отверстие, которое предназначено для его охлаждения.
О распределении температуры в шнеке можно судить по экспериментальным данным или на основе соответствующих математических моделей. Поскольку шнек вращается, а длина отверстия в шнеке, предназначенного для его охлаждения, много больше диаметра (примерно 100:1), то использование промышленных устройств, контролирующих температуру, затруднительно. Поэтому наиболее перспективным способом определения температуры шнека и возможным способои ее регулирования является использование математического моделирования.
При построении модели для определения температуры шнека вводятся следующие допущения. Процесс считается стационарным. Винтовой канал (рис. 1) заменяется кольцевым зазором (рис. 2) с учетом соотношения объемов винтовой нарезки и винтового канала. Сложное движение полимера в винтовом канале червяка заменяется стержневым в коаксиальном зазоре. Величина скорости определяется из массового расхода. Теплофизические свойства стали и полимера являются величинами постоянными, не зависимыми от температуры.
Я,
я,
1 ї-'
■Ко.
Л
Рис. 2. Схема модели для расчета температуры в шнеке
На основании сделанных упрощающих предположений задача по определению температуры в шнеке в осесимметричной постановке сводится к решению уравнения энергии для полимера [ 1]
1 д ( дТл г— дг
г дг
V
б 2Т Л
"бС
(1)
и уравнения теплопроводности для шнека
1 д ( дТ г-
г дг
дг
д 2Т ~д^
= О,
(2)
где Т - температура; - средняя скорость движения полимера в кольцевом зазоре модели; рп - плотность полимера; Сп - коэффициент
теплоемкости полимера; X п - коэффициент теплопроводности; qv -
удельная мощность диссипативного источника тепла в результате действия сил вязкого трения, распределение которой по длине определяется по вышеприведенной математической модели процессов движения и теплообмена пластицирующего экструдера.
Для получения полной математической формулировки задачи уравнения (1), (2) дополняются граничными условиями по температуре. Температура корпуса определяется технологическими режимами переработки:
Т\г== °- <3>
Если шнек не имеет отверстия, то в силу осесимметричности
задачи на оси граничное условие имеет вид
?L = о- (4)
дг |г=о
Если рассматривается шнек с осевым продольным отверстием, которое может быть использовано для охлаждения водой, то применяется граничное условие третьего рода:
=__^(г R - Т0)> (5)
dry=я01 А,м |г 01
где а - коэффициент теплоотдачи, Т0 - температура подвижной среды в отверстии, Xм - коэффициент теплопроводности стального шнека.
Граничное условие на входе (Q = 0) задается следующим образом:
Т = Т
\г_^От ^0 ’
Т - т
\г =ЕШ п ’
где Гшо - температура шнека при Q = 0, Гп - температура полимера на входе.
В силу того, что длина шнека много больше диаметра, граничное условие на выходе экструдера определяется из условия установления теплового потока:
дТ |
где L - осевая длина шнека.
На границе раздела металл-полимер должно выполняться граничное условие четвертого рода:
X — = Х8-^ • (8)
ЛМ - п п
дг |г=яш дг |г=яш
Задача решается методом конечных элементов. Предлагается следующая методика расчета температуры шнека: полимер совершает стержневое движение в кольцевом зазоре между шнеком и внутренней поверхностью цилиндра при условии, что гребни отсутствуют, при этом тепловой поток от корпуса к шнеку через зазор между гребнями учитывается итерационно. На каждом шаге итерации по разности температур между корпусом и шнеком определяется плотность теплового потока от корпуса к гребню нарезки червяка через радиальный зазор. Вычисленный таким образом тепловой поток прикладывается к шнеку. Итерационный процесс продолжается до тех пор, пока максимальная разность по температуре шнека между текущей и предыдущей итерацией в узлах расчетной области по абсолютному значению не будет меньше некоторой наперед заданной малой величины.
Численные исследования проведены для экструдера с классической геометрией шнека, геометрические размеры которого представлены в табл. 1. В качестве полимера выбран полиэтилен (ПЭ), теплофизические и реологические свойства которого приведены в табл. 2 и на рис. 3.
Таблица 1
Базовая геометрия шнека
№ п/п Геометрические размеры Значение
1 Внутренний диаметр цилиндра (корпуса), мм 160,0
2 Наружный диаметр шнека, мм 159,4
3 Шаг винтовой нарезки, мм 160,0
4 Ширина канала Ж, мм 137,3
5 Ширина гребня винтовой нарезки 51, мм 15,3
6 Длины геометрических зон загрузки, сжатия и дозирования Д21 / Д22 / Д23, витки 10/10/7
7 Суммарная длина шнека витки 27
8 Глубина канала в зоне загрузки Н\, мм 16
9 Глубина канала в зоне дозирования Н2, мм 4
10 Угол подъема винтовой линии © 17°39'
11 Радиальный зазор между гребнем шнека и корпусом 5, мм 0,3
12 Диаметр отверстия в шнеке, мм 48
Зависимость вязкости от скорости сдвига определяется степенным законом [2]. Индекс я в табл. 2 соответствует твердому состоянию полимера, ат- расплавленному.
Таблица 2
Реологические и теплофизические свойства полимеров
Полимер п Цо То р т р, Р ш *, * и
- Па-сп ос 1/°С ос кг/м3 Вт/(м-°С)
ПЭ 0,44 10825 160 0,018 110 919,0 779,0 0,335 0,182
с,
2,0------------------------1---------------
и 50 ;оо 150 200 7:ч:
Рис. 3. Зависимость удельной теплоемкости от температуры
Начальная температура гранулята Т2о равна 20 °С. Число оборотов шнека Мш = 60 об/мин. Температура цилиндра экструдера Тц на первых семи витках линейно растет с 92 до 200 °С, далее до конца не изменяется. Массовая производительность экструдера Со = 0,06 кг/с.
Рабочие характеристики экструдера с учетом определяемой температуры шнека находятся путем совместной итерационной процедуры решения задачи тепломассопереноса полимера в канале экструдера [3, 4] и представленной здесь задачи по определению температуры шнека.
Энергетический баланс экструдера определяется соотношением
а=&+а, . (9)
где Qi - изменение внутренней энергии полимера в единицу времени; Qф - мощность диссипативного источника тепла; Qx - мощность теплового потока подводимого (отводимого) через стенки канала.
На рис. 4 приведено изменение Qi, Qф и Qx по длине канала. Из рисунка видно, что сразу после окончания зоны загрузки разогрев материала в канале происходит только за счет диссипативной энергии, а корпус работает в режиме отвода тепла.
Величина Qx определяется суммой двух составляющих:
<2\ = бщ + (Ю)
где Qxц - тепло, подводимое (отводимое) через стенки цилиндра, а Qxш - через стенки шнека.
На рис. 5 показано изменение подводимого (отводимого) тепла от цилиндра и шнека к полимеру в единицу времени по длине канала.
Рис. 4. Изменение внутренней энергии полимера, диссипативной и подводимой (отводимой) через стенки канала составляющих тепла в единицу времени по длине канала
Из рис. 5 видно, что по длине сначала возрастает. На этом участке шнек отдает тепло полимеру. Максимальное значение достигает на длине, где происходит пересечение зависимостей средней температуры ТСр и температуры шнека Гш (рис. 5). Величина
на этой длине составляет 30 % от величины . Далее до конца канала тепло переходит от более горячего полимера к шнеку. Поскольку
на внутренней поверхности шнека задана адиабата, то величина Qxш в конце шнека стремится к 0. Таким образом, массивный стальной шнек, имеющий большое значение коэффициента теплопроводности, вносит достаточно существенный вклад в перераспределение энергетических параметров экструдера.
По изменению Qxц на рис. 5 видно, что корпус практически на
всей длине работает в режиме отвода тепла, за исключением первых трех витков.
Рис. 5. Изменение энергетических параметров и температуры по длине канала
Температура шнека Тш на последних 12,5 витках превышает температуру цилиндра Тц. К выходу средняя температура расплава
полиэтилена достигает 226 °С, температура шнека - 209 °С при температуре корпуса 200 °С.
Представленная в данной работе модель совместно с моделью [3] позволяет получать наиболее полную картину процесса переработки полимеров на одночервячных экструзионных агрегатах с учетом геометрических и технологических параметров, а также свойств перерабатываемого материала. Кроме того, с помощью данного подхода можно подбирать режимы охлаждения шнека водой пли другими теплоносителями с целью обеспечения необходимого качества изделия и увеличения производительности экструдера.
Библиографический список
1. Щербинин А.Г., Труфанова Н.М., Янков В.И. Пространственная математическая модель одночервячного пластицирующего экструдера. Сообщение 2: Математическая модель по определению температуры шнека // Пластические массы. - 2004. - № 8. - С. 38-40.
2. Янков В.И., Труфанова Н.М., Щербинин А.Г. Изотермическое течение аномально-вязких жидкостей в винтовых уплотнениях с продольной циркуляцией // Химическое и нефтегазовое машиностроение. - 2006. - № 6. - С. 3-5.
3. Щербинин А.Г., Труфанова Н.М., Янков В.И. Пространственная математическая модель одночервячного пластицирующего экструдера. Сообщение 1: Математическая модель процесса тепло-массопереноса полимера в канале экструдера // Пластические массы. -2004,-№6.-С. 38-41.
4. Щербинин А.Г., Труфанова Н.М., Янков В.И. Пространственная математическая модель одночервячного пластицирующего экструдера. Сообщение 3: Проверка адекватности модели // Пластические массы. - 2005. - № 5. - С. 43-45.
Получено 05.09.2012
