Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Переработка полимеров

..pdf
Скачиваний:
10
Добавлен:
12.11.2023
Размер:
1.58 Mб
Скачать

Рис. 9.6. Сечение канала в зоне плавления

Рассчитывая зону плавления, прежде всего необходимо определить длину зоны плавления, что потребует расчета полей скоростей и температур в материале.

Обозначим ширину твёрдого слоя X. Необходимо получить зависимость X X Z . w X Z – ширина области расплава,

X Z H – локальный расход материала.

Рассмотрим элементарный объём в канале зоны плавления (рис. 9.7), расположенный перпендикулярно поверхности раздела пробка гранул–плёнка расплава. Материал пробки движется с локальной скоростью sz , направленной вдоль канала «червя-

ка». Скорость пробки относительно поверхности цилиндра Vj

находится как разность векторов – окружной скорости и локальной продольной скорости.

 

V

j

 

b

 

sz ,

(9.18)

111

модуль вектора искомой скорости определяется как

V

j

 

2

2

2

sz

.

 

 

bz

sz

 

bz

Рис. 9.7. Элементарный объем полимера в зоне плавления

Введем следующие обозначения:

sy – модельная скорость, вводимая в рассмотрение исхо-

дя из сделанных допущений (плавление происходит только на верхней стороне пробки);

b

– окружная скорость;

bz

– компонента окружной скорости в направлении оси z;

bx

– компонента окружной скорости в направлении оси x;

sz

– скорость пробки вдоль оси z;

sy

– скорость пробки вдоль оси y;

Tb

– температура корпуса;

Ts

– температура загружаемого материала.

Задачу определения длины зоны плавления разобьем на ряд задач:

1.Рассмотрим движение и теплообмен в расплаве полимера.

2.Рассмотрим теплообмен в твердой фазе.

112

3.Составим баланс потоков тепла на границе раздела фаз

(y = 0) (см. рис. 9.7).

4.Рассмотрим баланс масс в жидкой фазе.

5.Рассмотрим баланс масс в твердой фазе.

6.Определим длину зоны плавления.

1. Рассмотрим уравнение движения и уравнение энергии в расплаве. С учетом сделанных допущений они примут вид

 

 

 

2

j

0 .

 

(9.19)

 

 

y2

 

 

 

 

 

 

 

 

 

2T

2

0 ,

(9.20)

 

m y2

 

 

 

 

 

где j (y) – скорость жидкой фазы с учетом введения относитель-

нойскорости(9.18), позволившей упростить граничныеусловия. Граничныеусловия для скоростиитемпературы следующие:

j

 

y 0 0;

j

 

y Vj ;

T

 

y 0 Tm ;

Tj

 

y Tb .

 

 

 

 

 

 

 

После интегрирования уравнения (9.19) и подстановки граничных условий получим выражение для скорости (течение Куэтта)

 

 

 

Vj

y .

(9.21)

j

 

 

 

 

 

Продифференцируем уравнение (9.21), определив тем самым скорость сдвига, и подставим полученное выражение в уравнение энергии.

2T

 

 

V 2

y

2

 

 

i .

 

 

 

m

 

Дважды проинтегрируемполученное выражение, будем иметь

T

 

 

V 2

y C1 ,

y

 

 

i

 

 

m

 

 

113

 

 

 

 

 

V

2

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

i

 

y2 C y

C

2

.

(9.22)

 

 

 

 

 

 

 

 

2 m

 

 

 

 

 

 

 

 

 

 

 

 

 

Найдем С1 и С2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

T ; T

 

 

 

V

2

2

C T

2

 

 

 

 

 

i

 

 

 

 

 

 

m b

 

 

 

 

 

 

 

 

 

1

 

 

 

m

 

 

 

 

 

2 m

 

 

 

 

 

 

 

 

 

 

 

Ñ

T

T

 

 

 

 

 

V

2

.

 

 

 

 

 

b

 

 

m

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 m

 

 

 

 

 

Подставим С1 и С2 в (9.22), получим выражение для температуры в жидкой фазе

T

 

V

2

y2

T

T

 

 

i

 

 

b

m

 

 

 

 

 

 

 

 

 

2 m

 

 

 

 

 

Vi

2

 

 

 

 

y T

. (9.23)

 

 

 

 

 

 

 

m

 

2 m

 

 

 

 

 

2. Рассмотрим уравнение энергии в твердой фазе

c

 

 

 

dT

 

 

d 2T .

(9.24)

s

s

 

sy

dy

 

 

 

s

dy2

 

Граничные условия:

 

 

 

 

 

 

 

 

 

 

T

 

y 0

Tm ;

T

 

y Ts .

 

 

 

 

 

 

 

Сделаем замену переменной, обозначим dTdy u ,

тогда уравнение (9.24) примет вид

cs s sy u du .

s dy

Проинтегрируемпоследнееуравнение, разделив переменные.

cs s sy dy du ;s u

114

 

 

 

 

 

 

 

 

 

cs s sy

y

 

 

 

 

 

 

 

 

u C3

e

 

 

 

s

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тогда, сделав обратную подстановку,

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

cs s sy y

 

 

 

 

 

 

 

C3 e

 

s

 

 

 

 

 

 

.

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Проинтегрируем последнее выражение.

 

 

 

dT C3 e

 

cs s sy

y dy ,

 

 

 

 

s

 

 

 

 

 

 

 

s

 

 

 

 

 

 

cs s sy y

 

 

T C3

 

 

 

 

 

e

 

s

 

 

 

 

C4 .

c

s

 

sy

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найдем C3 и С4.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

4

T ;

С

 

 

Tm Ts

 

s

 

 

sy

c .

 

 

 

 

s

3

 

 

 

 

 

s

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Окончательно имеем выражение для температуры в твердой фазе

T Tm Ts e

cs s sy

y Ts .

 

s

(9.25)

3. Составим баланс потоков тепла на границе раздела фаз

(y = 0).

Разность между количеством тепла, подводимого к разделу фаз и отводимого в пробку гранул, есть тепло, расходуемое на плавление материала с единицы поверхности раздела фаз в единицу времени, что определяется скоростью sy.

qm y 0 qs y 0 sy s ,

где Джкг – скрытая теплота плавления, или количество

энергии, которое необходимо для перевода 1 кг вещества из твердого состояния в жидкое.

115

Раскроем выражения для потоков тепла, используя уравнения (9.23) и (9.25) на границе раздела фаз, то есть при y = 0.

q

 

 

 

 

 

 

T

 

 

 

 

T

T

 

 

Vj2

;

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

m

 

 

 

m

b

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

y 0

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

y

 

y 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

s

 

y 0

 

 

 

T

 

 

T

 

T

 

s

c

s

 

sy

.

 

 

 

 

 

 

s

 

 

 

 

m

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

y 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда уравнение баланса тепла примет вид

T T

 

 

Vj

2

 

T

T

 

c

 

 

 

 

 

 

 

b m

m

 

 

 

 

 

m

s

s

 

sy

 

sy

 

s

. (9.26)

 

 

 

 

 

2

 

 

s

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4. Для определения неизвестной, искусственно введенной

скорости sy рассмотрим баланс масс в жидкой фазе.

Количество материала, расплавляющегося c поверхности пробки на длине dz в единицу времени, определяется выра-

м кг

м

2

 

кг

жением, имеющим размерность расхода

с м

3

 

 

.

 

 

 

 

 

с

 

 

 

sy

s X dz ;

 

 

 

 

 

количество материала, уносимого в бассейн расплава в на-

правлении оси х в единицу времени,

 

 

 

 

 

 

 

 

bx m dz ;

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

количество

материала,

 

утекающего в направлении оси z

в единицу времени на участке длиной dz,

 

 

 

 

 

bz sz

m

X

 

 

bz sz

m

X

 

 

.

 

 

 

 

 

2

 

 

z dz

2

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда уравнение баланса масс выразится следующей алгеб-

раической суммой:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

bz sz m X

 

z dz X

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

(9.27)

 

bx m dz sy s X dz dz,

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

116

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где – удельная скорость плавления, размерность: кгм с –

количество

расплавившегося

 

материала в единицу времени

с единицы длины.

 

 

 

 

 

Разделим (9.27) на dz и перейдем к пределу. Пренебрегая

величинами малого порядка, получим

 

 

 

 

 

sz m

d X 0

 

 

sy s X . (9.28)

 

bz

 

 

 

 

bx m

 

2

dz

 

 

 

 

 

2

 

Первый член (9.28) равен нулю, так как величина мала и мало его изменение по длине канала, тогда определим выражение для искомой скорости.

sy

bx

m

.

(9.29)

 

 

2

s X

 

С другой стороны, из тождества (9.28) найдем удельную скорость плавления

 

 

 

 

 

bx m .

 

 

 

 

(9.30)

 

 

 

 

 

 

2

 

 

 

 

 

 

 

Подставим (9.29) в (9.26), получим

 

 

 

 

 

 

T

T

Vj 2

Tm Ts

cs s bx m

 

bx

m s

.

b

m m

 

 

 

 

 

 

 

 

 

2

 

2 s X

2

s

X

 

 

 

 

 

 

 

 

Приведем выражение к общему знаменателю и выразим .

2 m Tb Tm Vj2 X Tm Ts cs bx m 2 0 ,

 

 

 

 

 

 

 

2

 

 

 

 

 

2 m Tb

Tm Vj

 

X

А X ,

(9.31)

 

T

T

c

s

 

bx

 

m

 

 

 

 

 

 

m

s

 

 

 

 

 

 

 

где А – константа, для определения которой все величины известны.

117

Рис. 9.8. Площади поперечного сечения

Определим выражение для скорости плавления через ширину твердой фазы, используя зависимости (9.30) и (9.31).

 

1

1

 

bx

m A X 2

Ф X 2 .

2

 

 

 

5. Рассмотрим баланс масс в твердой фазе.

На рис. 9.8. схематично показаны две площади поперечного сечения твердой фазы, отстоящие друг от друга на величину dz. Очевидно, что различие между ними обусловлено процессом плавления полимера, интенсивность которого обусловлена реализуемой скоростью плавления.

Закон сохранения масс для твердой фазы может быть записан так:

sz s X H z sz s X H z dz dz ,

где первый член уравнения – это массовый расход материала в сечении z, второй член уравнения в сечении z + dz. Для неизменной глубины канала (H = const), разделив каждый член уравнения на dz и перейдя к пределу, получим дифференциальное уравнение

 

 

 

 

1

 

 

 

 

dX

 

ФX 2

 

.

 

dz

sz s H

 

 

 

 

 

 

Разделим переменные.

 

 

 

 

 

 

 

 

dX

 

Фdz

.

 

1

 

 

 

 

sz s H

 

 

 

X 2

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегрируем.

 

 

 

 

 

 

 

 

1

 

 

Фz

 

 

 

 

 

2X 2

 

 

C ;

 

 

 

 

 

sz s

H

 

 

 

1

 

 

 

 

 

 

118

X sz Фs zH 2 C1 2 .

Найдем константу интегрирования C1 из условия: в начале зоны плавления ширина твердой фазы равна ширине канала w

(при z = 0 и X = w).

 

 

Фz

2

 

X W 1

 

 

.

(9.32)

1 2

 

 

sz s H 2w

 

 

Выражение (9.32) определяет зависимость ширины пробки гранул от продольной координаты и различных факторов.

6. Для определения длины зоны плавления, которая характеризуется тем, что ширина твердой фазы становится равной нулю (полное расплавление материала), приравняем нулю выражение, стоящее в левой части (9.32).

X 0

z zпл .

 

0 1

 

 

Фzпл

,

 

 

1 2

 

 

 

2 sz s H w

 

тогда

 

 

 

 

 

zпл

2 sz s H w1 2 .

(9.33)

 

 

 

Ф

 

Из полученных уравнений видно, что профиль пробки в канале постоянной глубины имеет параболическую форму.

Длина зоны плавления пропорциональна массовому расходу и обратно пропорциональна интенсивности плавления Ф. Влияние технологических параметров можно оценить через параметр Ф. Так, увеличение частоты вращения шнека при постоянном расходе материала приводит к возрастанию скорости плавления, поскольку, во-первых, увеличиваются тепловыделения за счет вязкого трения, во-вторых, улучшаются условия отвода расплава полимера от горячей поверхности корпуса (возрастает ско-

119

рость bx ). При увеличении температуры корпуса сначала про-

исходит возрастание интенсивности плавления, так как количество тепла, поступающего за счет теплопроводности, пропор-

ционально выражению Tb Tm m возрастает. Однако в связи

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

иведет к сокращению длины зоны плавления.

Кподобным выводам можно прийти, рассмотрев процесс плавления в канале «червяка» с коническим сердечником, с начальной высотой канала Н и углом конусности.

Канал с переменной глубиной (Н = var)

В зоне сжатия глубина канала изменяется линейно относительно продольной координаты z.

 

 

 

 

 

H H1 Az ,

 

 

 

(9.34)

при z 0,

H H1 .

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение массового баланса в твердой фазезапишется в виде

 

 

 

 

 

d HX

 

 

 

ФX

1

 

 

 

 

 

 

 

 

 

 

2

,

 

 

 

 

 

 

 

 

dz

sz s H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

H1

Az dX

AX

ФX 2

.

(9.35)

 

 

 

 

 

 

 

 

dz

 

 

 

 

sz s

 

 

Разделим

каждый член

уравнения (9.35)

на величину

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 X 2 H1 Az , получим

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

1

 

dX

 

AX 2

 

 

 

 

 

Ф

 

 

 

 

 

dz

 

 

 

.

 

 

1

2 H1 Az

sz s 2 H1 Az

 

2X 2

 

 

 

 

 

 

 

 

 

 

 

 

120