close

Вход

Забыли?

вход по аккаунту

?

Расчет распределения температуры расплава полимера по высоте канала шнека экструдера..pdf

код для вставкиСкачать
Вестник СГТУ. 2011. № 4 (62). Выпуск 4
УДК 678.057
М.Г. Хаметова
РАСЧЕТ РАСПРЕДЕЛЕНИЯ ТЕМПЕРАТУРЫ РАСПЛАВА ПОЛИМЕРА
ПО ВЫСОТЕ КАНАЛА ШНЕКА ЭКСТРУДЕРА
Решена задача распределения температуры расплава полимера по
высоте канала шнека одношнекового экструдера с учетом противотока
расплава по длине шнека.
Одношнекового экструдер, экструзионные машины, распределение
температуры расплава полимера.
M.G. Khametova
CALCULATION OF TEMPERATURE DISTRIBUTION POLYMER
MELTS ALONG EXTRUDER SCREW CHANNAL HEIGHT
Solved the problem of temperature distribution in the polymer melt channel
height screw single screw extruder with a melt backflow into account the length of
screw.
Single-screw extruder, extrusion machine, the temperature distribution of
the polymermelt.
Введение
Современные экструзионные машины, предназначенные для переработки полимерных материалов и изготовления пленки, листа, профилей, труб и т.д., должны подавать в
формующий инструмент экструдера термически однородный расплав с постоянной высокой
производительностью. Для того, чтобы обеспечить выполнение этого требования, управление процессом экструзии должно быть основано на полном понимании основных закономерностей работы этих машин.
Статья посвящена нахождению стационарного распределения температуры течения
неньютоновского расплава полимера по высоте прямоугольного канала шнека одношнекового экструдера. Предполагается, что расплав полимера подчиняется реологическому степенному закону.
Расчет распределения температуры проводится в рамках упрощенной теории потока [1]. В работах [1-4] исследовалась такая температурная зависимость, в том числе с помощью ЭВМ, но в них не был предложен метод расчета распределения температуры расплава
по высоте канала шнека y . Необходимость такого исследования связана с выбором технологических режимов переработки полимеров. В данной статье приводится аналитическое решение этой задачи, а именно, выводится уравнение распределения температуры расплава по высоте канала шнека и дается его решение. Из него следует, что имеет место эффект саморазогрева расплава, который наиболее выражен на нижней стенке канала и вызван вращением
шнека экструдера. Важно отметить, что саморазогрев расплава приводит к его деструкции.
122
Новые материалы и технологии
Постановка задачи
Обозначим z − положительное направление оси шнека, y ∈ [0, H ], H − высота канала шнека экструденра, η − вязкость расплава полимера, v z − скорость течения расплава полимера вдоль оси z , ∂v z / ∂y − скорость сдвига, λ коэффициент теплопроводности расплава,
T ( y ) − температура расплава в точке y .
Предположим, что вязкость расплава описывается выражением
η ( y ) = η 0 b(∂v z / ∂y ) n−1 exp{− α (T ( y ) − T ( H ))},
где α , b − эмпирические константы; n − индекс течения неньютоновской жидкости; η 0 −
начальная вязкость.
Ранее в работе [5] было дано теоретическое обоснование эффекта саморазогрева расплава полимера в одношнековом экструдере в результате неизотермического неньютоновского стационарного течения расплава, вызванного вращением шнека экструдера. Однако
решение этой задачи не учитывало наличие противотока расплава по длине шнека.
Известно [1], что в рамках упрощенной теории потока скорость сдвига ∂v z / ∂y и T ( y )
удовлетворяют системе уравнений:
∂p / ∂z = ∂ (η∂v z / ∂y ) / ∂y ,
(1)
λ d 2T dy 2 + η (dv z dy )2 = 0 ,
(2)
где p − давление.
Константу b можно определить, если известны значения вязкости расплава η ( H ) и
скорости сдвига dv z ( H ) / dy при y = H :
b = [η ( H ) / η o ][dv z ( H ) / dy ] .
Приведем граничные условия [1] для уравнения (1):
∂v z ( y ) / ∂y y =0 = 0 , v z ( y ) y =H = Vz ,
(3)
а для уравнения (2) они имеют вид [1]:
∂T ( y ) / ∂y y =0 = 0 ,
(4)
1− n
T ( y)
y=H
= T (H ) ,
где Vz − скорость движения верхней стенки канала шнека экструдера; T (H ) − температура
термостатирования расплава на верхней стенке канала экструдера.
Построение решения краевой задачи.
Пусть ξ ∈ [0,1] − новая независимая переменная, причем:
y (ξ ) = Hξ 2 n /( 3n +1) .
(5)
Положим:
θ (ξ ) = {α (T ( y ) − T ( H ))}/ n y = y (ξ ) .
(6)
Отметим, что θ (ξ ) − имеет смысл безразмерной температуры. Тогда из уравнений (1),
(2) следует, что θ (ξ ) удовлетворяет уравнению:
d 2θ / dξ 2 + β exp(θ (ξ )) = 0 ,
(7)
где β = (α / nλ ){p
/(η 0 b) } (2n /(3n + 1)) H
− константа. При этом θ (ξ ) удовлетворяет следующим граничным условиям:
dθ (ξ ) / dξ ξ =0 = 0, θ (ξ ) ξ =1 = 0 .
(8)
( n +1) / n
z
1/ n
2
( 3 n +1) / n
Сначала понизим порядок уравнения (7). Для этого умножим левую и правую части
(7) на dθ (ξ ) / dξ , имеем
0,5 d [dθ (ξ ) / dξ ] / dξ = − β exp θ (ξ )dθ (ξ ) / dξ .
2
123
Вестник СГТУ. 2011. № 4 (62). Выпуск 4
Интегрируя последнее равенство по ξ в пределах от нуля до ξ , имеем:
(dθ (ξ ) / dξ ) 2 − (dθ (0) / dξ ) 2 = −2 β (exp θ (0) − exp θ (ξ )) .
Отсюда следует, что θ (ξ ) удовлетворяет уравнению
dθ (ξ ) / dξ = ± (dθ (0) / dξ ) 2 + 2 β (exp θ (0) − exp θ (ξ )) ,
которое с учетом первого граничного условия примет вид:
dθ (ξ ) / dξ = ± 2 β (expθ (0) − expθ (ξ ) .
(9)
Отметим, что
− если θ (0) ≥ θ (ξ ) для любого ξ ∈[0,1] , то правая часть является действительной функцией;
− если существует ξ * ∈ [0,1] такое, что θ (0) < θ (ξ *) , то в этом случае dθ (ξ ) / dξ будет
комплексно значимым. Следовательно, θ (ξ ) будет осциллирующей функцией ξ и практически значимого интереса не представляет.
Поэтому будем рассматривать случай, когда θ (0) ≥ θ (ξ ) для любого ξ ∈ [0,1] , т.е. точка ξ = 0 является точкой максимума функции θ (ξ ) .
Требуется выбрать знак в правой части выражения (9). С этой целью рассмотрим (7) и
выполним интегрирование его левой и правой частей в пределах от нуля до ξ . Имеем с учетом граничного условия (8):
dθ (ξ ) / dξ = − (2β (expθ (0) − expθ (ξ )) .
(10)
Из (10) следует, что
dθ (ξ ) / dξ ξ = 0 = 0.
Найдем теперь общее решение нелинейного уравнения (10). Сначала заметим, что для
любого С ≠ 0 и x ≥ C :
d (arch( x / C ) / dx = 1 / x 2 − C 2 .
Поэтому неопределенный интеграл будет иметь вид:
2
2 −1 / 2
∫ ( x − C ) dx = archx / h + C ,
(11)
где С – постоянная интегрирования.
Обозначим 1 / C = exp θ (0) / 2 . Поскольку переменные в уравнении (10) делятся, то его
можно переписать в виде
-0,5 [exp(−θ / 2)dθ ] / (exp(−θ / 2) 2 − C 2 = β / 2C −1 dξ .
Отсюда в силу (11) имеем
exp[(θ (0) − θ (ξ )) / 2] = ch( β / 2 exp(θ (0) / 2ξ ) + C1 ,
(12)
где C1 − постоянная интегрирования. Очевидно, что C1 =0.
Таким образом, из (12) следует, что
θ (ξ ) = θ (0) − 2 ln ch( β / 2ξ exp(θ (0) / 2)) ,
(13)
Осталось найти значение θ (0) . Для этого воспользуемся вторым граничным условием
(8), из которого следует, что
expθ (0) / 2 = ch( β / 2 exp(θ (0) / 2)) .
(14)
Из выше приведенных выкладок следует, параметр θ (0) − это значение обезразмеренной температуры расплава на нижней стенке канала шнека экструдера.
Отметим, что нелинейное уравнение (14) имеет решение при β / 2 < 1. Если это условие выполняется, то уравнение (14) имеет два положительных корня, которые обозначим как
θ1 (0) и θ 2 (0) . В действительности на практике реализуется наименьший из них по величине
124
Новые материалы и технологии
в силу следующего соображения, а именно: в стационарном режиме любая динамическая система занимает такое состояние, которому соответствуют минимальные затраты энергии
(температуры).
Приведем оценки сверху и снизу для корней θ1 (0) и θ 2 (0) . Без ограничения общности можно считать, что θ1 (0) < θ 2 (0) <0.
Заметим, что из (14) следует неравенство:
expθ (0) / 2 = ch( β / 2 exp θ (0) / 2) ≥ 1 + ( β / 2) expθ (0) .
(15)
Из (15) следует, что имеет место неравенство:
2 ln (2 − 2 1 − β ) / β ≤ θ1 (0) < θ 2 (0) ≤ 2 ln (2 + 2 1 − β ) / β .
Отсюда
θ1 (0) ≈ 2 ln 2 (1 − 1 − β ) / β .
(16)
Таким образом, обезразмеренная температура θ (ξ ) принимает максимальное значение в точке ξ = 0 , т.е. θ (0) ≥ θ (ξ ) . Последнее означает, что происходит саморазогрев расплава полимера, и температура расплава принимает наибольшее значение у нижней стенки
канала шнека экструдера.
Следует отметить, что полученные здесь результаты полностью подтверждаются исследованиями, проведенными Раувендаалем К. [6]. Сначала экспериментально, затем с помощью программы с использованием метода конечных элементов Раувендааль К. установил профиль температуры расплава в канале шнека одношнекового экструдера. В качестве
объекта изучения был выбран полиэтилен высокого давления, который подвергали переработке в экструдере с диаметром рабочего цилиндра 38 мм и скоростью вращения шнека 100
об/мин. Эта программа позволила проводить трехмерный расчет скоростей, давлений и температур для любой точки в канале шнека. Им было показано, что температура расплава полимера в любой точке канала шнека значительно превышает температуру рабочего цилиндра
(температуру термостатирования). Самые высокие значения, примерно на 310 С выше температуры рабочего цилиндра, наблюдались в канале вблизи сердечника шнека (на нижней
стенке канала). Время пребывания расплава больше всего в нижней части канала. Сочетание
высоких температур расплава и большого времени пребывания в этой области канала делают
деструкцию расплава весьма вероятной.
Пример расчета величины саморазогрева расплава поликарбоната. Была рассчитана
относительная величина саморазогрева расплава поликарбоната при экструзии в одношнековом экструдере по данным, приведенным в [7-9]. Приведем данные о расплаве поликарбоната и его течении:
− энергия активации вязкого течения E = 108 Кдж/моль ;
− коэффициент теплопроводности λ =0,1884 Вт/мК 0 ;
− начальная вязкость η 0 =1,641 10 3 н сек/м 2 ;
[
]
[
[
]
]
− константа b =57,07 сек n −1 ;
− температура термостатирования 553 0 К ;
− производная от давления вдоль оси шнека p z =6 10 6 н/м 3 ;
− высота канала шнека экструдера 0,0063 м ;
− индекс течения n =0,4.
Поскольку
β = α ( p z ) ( n+1) / n /(η 0 b)1 / n [2n /(3n + 1)]2 ( H ) (3n+1) / n / λn ,
где значение константы α можно получить с помощью преобразования ФранкаКаменецкого формулы Аррениуса [10]:
α = E / R[T ( H )]2 = 0,04245 рад −1 .
[
]
125
Вестник СГТУ. 2011. № 4 (62). Выпуск 4
Здесь R − универсальная газовая постоянная, R =8,31 Кдж/кмоль град .
Отсюда следует, что
следует, что
β / 2 ≅ 0,1556 . Заметим теперь, что с одной стороны, из (17)
θ (0) ≅ 2 ln 2(1 − 1 − β / β ) = 8,803 ,
с другой стороны,
θ (0) = αT ( H )[(T (0) − T ( H )] / nT ( H ) .
Отсюда следует, что
[T (0) − T ( H )] / T ( H ) = n[θ (0) / αT ( H )] .
Поскольку
αT ( H ) / n = E / nRT ( H ) = 58,69 .
Следовательно,
(T (0) − T ( H )) / T ( H ) ≅ 0,15,
т.е. величина саморазогрева составляет 15% от температуры термостатирования и составляет
для расплавов поликарбонатов 636 0 К . Хорошо известно [9], что при такой температуре поликарбонат деструктирует.
Выводы
Показано, что при экструзионной переработке полимеров возникает саморазогрев
расплава. Причем наибольшее значение температуры расплава достигается на нижней стенке
канала шнека экструдера.
Предложена методика расчета величины саморазогрева при экструзионной переработке полимера.
ЛИТЕРАТУРА
1. Mckelvey J.M. Polymer Processing / J.M. Mckelvey. M.: Chemistry, 1965. 444 p.
2. Bernhardt E.C. Processing of Thermopastic Materials / E.C. Bernhardt. M.: 1962, 747s.
3. Middleman S. The Flow high Polymers / S. Middleman. M.: 1971. 256 s.
4. Truphanova N.M. Development of calculation methods of plastification process and
working parts extruder equipment for plastics. Jhesis… doctor of science. M.: MSACM, 1994.
5. Khametova M.G. Polymers melts self-heating in estrusion process / M.G. Khametova,
E.A. Zukov, S.N. Rukavizin // ММТТ-22 : sp. MNK. V.9. Pskov: PPI, 2009. 73 p.
6. Rauvendaal K. Influence and removal extrusion problems / K. Rauvendaal, Pilar Noryega,
H. Haris. SP.: Profession. 2008, 328 p.
7. Каlinchev E.L. Properties and processing of thermoplastics / E.L. Каlinchev,
M.B. Sacovcheva. М.: Chemistry, 1983. 286 p.
8. Thermophisical and reological polymers characteristics / Under reduction of Lipatov J.S.
Кiev.: Science idea, 1977. 244 p.
9. Khametova M.G. Calculation of extruder productivity in the light of destruction processing polymer / M.G. Khametova // Тheoretical basic of chemical technology. М.: V.24, 1990. №6. P.775-783.
10. Frank-Kamenesky D.А. Diffusion and heat transfer in chemical kinetics / D.А. FrankKamenesky. М.: Science, 1967. 491 p.
Хаметова Маргарита Григорьевна –
кандидат технических наук, доцент, декан машиностроительного факультета Московского
государственного университета инженерной экологии
Статья поступила в редакцию 5.07.11, принята к опубликованию 24.10.11
126
Документ
Категория
Без категории
Просмотров
19
Размер файла
152 Кб
Теги
температура, канал, расплавах, шнека, pdf, экструдера, расчет, распределение, высота, полимеры
1/--страниц
Пожаловаться на содержимое документа