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

Теоретические основы переработки полимеров (механика процессов)

..pdf
Скачиваний:
45
Добавлен:
15.11.2022
Размер:
21.7 Mб
Скачать
\ЧЧ\\\\ч\\ч\\^
Рис. VUI. 4. Схема движения расплава в дозирующей зоне червяка:
а Поступательный поток; б — циркуля­
ционное течение; в — итоговое винтовое движение.

добно течению жидкости в винто­ вых насосах — по винтовой тра­ ектории (рисVIII.4). \Принято представлять это течение как сумму двух независимых движе­ ний [2, 3]: поступательного движе­ ния расплава вдоль оси винтового канала (рис. VIII.4,б) и цирку­ ляционного (кругового) движе­ ния ,в плоскости xmfr нормаль­ ной ' к оси винтового канала (рис. VIII.4,б).

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

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

т^^.МШ1^.л^с-и--аа--с.чет тедл^выдеМющегося вследствие интен­ сивной деформации сдвига. Одновременно идет процесс гомогенизш^ш1^атпТлав?ГТТроисходит окончательное расплавление мелких включений и выравнивание температурного поля. Для нормальной работы экструдера необходимо, чтобы расплав, поступающий к ра­ бочему инструменту (к головке), имел заданную, однородную по сечению температуру. Поэтому время пребывания расплава в зоне

дозирования должно быть достаточно для его прогрева и гомоге­ низации.

Для правильного анализа процесса экструзии надо рассматВивать совокупность всех его стадий, соблюдая условие постоянс£ва материального расхода полимера для любого сечения червяка.

VIII. 3. ОСНОВНЫЕ УРАВНЕНИЯ ДВИЖЕНИЯ РАСПЛАВА В ЗОНЕ ДОЗИРОВАНИЯ ЧЕРВЯКА ЭКСТРУДЕРА. ПОСТАНОВКА ЗАДАЧИ

Математическая модель зоны дозирования является результатом совместного решения системы уравнений, которые выражают законы сохранения массы, энергии и количества движения при да_ минарном течении, с уравнениями, описывающими физическое Со_ стояние перекачиваемой жидкости [40—48].

Наиболее последовательно модель винтового изотермическ0го течения ньютоновской жидкости была рассмотрена в 1953 г. в се_ рии работ американских авторов Карлея, Маллука и Мак-Ке^ви [44—47], положивших начало гидродинамическому подходу к оци_ санию процесса экструзии.

Использование упрощенной модели позволило разобраться в механизме движения расплава. Экспериментальные исследовация Маддока [33—35], Эккера [48] и других показали, что качеств^, ная картина хорошо согласуется с данными опыта. Однако попц1тки применения этой модели для целей инженерного расчета оца_ зались безуспешными, поскольку при ее построении игнорироца_ лись основные особенности расплава, такие как аномалия вязц0. сти, температурная зависимость вязкости, наличие процессов тегь лообмена.

Впоследствии задача о винтовом движении аномально-вязкой жидкости неоднократно рассматривалась в ряде отечественных и зарубежных работ. Так, изотермическое течение «степенной» Жид­ кости рассматривалось в работах Балашова [49], Торнера [25, gg 51], Бостанджияна [52], Мори [53], Глайда [54], Якоби [55], Крю­ гера [56], Кроссера [57], Мак-Келви [58] и других авторов. Винто­ вое течение бингамовой жидкости рассмотрено в работах Столица и Бостанджияна [59], Тябина [60], Де-Ховена [61] и других. Дела­

лись также попытки построения моделей течения эластичных жид­ костей [62]. Однако и эти работы не позволяли моделировать ре_ альный процесс, поскольку в них

не учитывалось существование теплообмена.

Чтобы удовлетворительно опи­ сать процесс экструзии, матема­ тическая модель течения распла­

ва в пределах зоны дозирования должна учитывать основные осо­

бенности процесса: существова­

 

 

 

 

ние аномалии вязкости; взаимное

 

Ось чердяна

влияние

циркуляционного

и по­

 

 

 

 

ступательного

течений;

влияние

 

 

 

 

тепла,

выделяющегося в

резуль­

 

 

 

 

тате вязкого трения, и теплооб­

 

 

 

 

мена с окружающей средой на

 

 

 

 

температуру и эффективную вяз­

 

 

 

 

кость расплава.

 

 

 

 

 

 

 

 

Очевидно, что для работы экс­

 

 

 

 

трудера

безразлично, вращается

 

 

 

 

ли червяк внутри

неподвижного

 

 

 

 

корпуса

или

наоборот,

 

корпус

Рис. VI 11.6.

Развертка

винтового

канала

вращается относительно червяка.

По этой

причине для упрощения

червяка на

плоскости:

 

 

/ — стенка канала; 2—развертка

червяка;

будем

считать

корпус

вращаю­

3— поверхность корпуса;

4—канал.

щимся

 

[42—45].

Расположение

связанной с червяком неподвиж­ ной системы координатных осей х, у, z и вспомогательной оси / показано на рис. VIII. 5. Ось z ориентирована вдоль оси винтового канала червяка, ось I— вдоль оси червяка. Ввиду малости отно­ шения глубины винтового канала h к радиусу /?, составляющего [2] у большинства машин 0,005—0,050, можно пренебречь кривиз­ ной канала червяка и развернуть его на плоскость, как это пока­ зано на рис. VIII. 6. При этом корпус будет изображаться беско­ нечной плоскостью, движущейся над развернутыми каналами в

перпендикулярном к оси / направлении.

 

 

Течение считаем установившимся. Массовыми силами и силами

инерции пренебрегаем. С учетом этих допущений

получим, что

уравнения равновесия в напряжениях сведутся к виду:

др

дРух

 

дРгх

 

дх

ду

+

дг

 

дР

дрху

,

дргу

(VIII. 1)

ду

дх

+

дг

 

дР

друг

I

дРхг

 

Я-7

 

Л

 

Среду считаем несжимаемой:

 

■dvx

 

 

dv2

 

дх

 

 

dz = 0

(VIII. 2)

Связь между компонентами тензора напряжений рц и тензора скоростей деформаций гц определяется выражением:

Р ц ^ ^ а ^ и Ч . /' = *. У . Z)

(VIH.3)

или в развернутой форме

 

(' дох

,

до у-

 

Рху — Ча у. ду

"*■

дх .

 

 

'д о у

,

дог '

(VIII. За)

Руг =

Ча (ч дг

'

ду .

 

f dvz

,

дох -

 

Ргх =

Ча (, дх

1

дг .

 

Рхх — Руу — Pzz =

Р

 

В качестве реологического уравнения используем обобщенный степенной закон:

Ча =

Voe~b (r~Tl) ( 4- b ') 2"

(VIII. 4)

где

м-о, п — константы материала; b — температурный

коэффициент вязкости;

12 — квадратичный инвариант тензора скоростей деформаций, определяемый вы­ ражением (III. 23).

Уравнение энергетического баланса, составленное для устано­ вившегося режима в предположении, что все теплофизические ха­ рактеристики не зависят от температуры, имеет вид:

рсРI

дТ .

дТ ,

дТ Л _

—пь

\ д Ч

1

д Ч

д Ч

'■ ~dZ +

V y ~dj +

V z ~te

)

 

ктп

1. дх 2

ду 2

1 д г2

 

dvx

dvu

 

дог

1

п

f дох

 

д° у '

+ Рхх

 

 

 

1

 

дх + Р у у ду

+ ргг

дг

1

Рху \

1

дх

,) +

 

 

 

\ ду

 

( dvx

dvz \

f d v u

 

dvz \

 

 

 

(VIII. 5)

+ ^ { - i r + - d r ) + ^ y b r + w )

где p — плотность расплава; cp — теплоемкость расплава; Т — температура рас­ плава; km— коэффициент теплопроводности расплава.

Граничные условия для такой полной дифференциальной мо­ дели зоны дозирования имеют вид:

vx = vy = vz = 0

 

при

у —

О

 

 

 

Г = Г2

 

 

при

х = 0

и

х = w

(для всех у)

(VIII. 6)

Vx = Ux>

vz =

Uz

при 0 < х < w

и у = h

 

где Т2 и

Тw— температура

червяка истенки

корпуса, соответственно.

 

Уравнения (VIII. 1) — (VIII.5) образуют полную систему урав­ нений относительно восьми функций (Р, рц, Vu Т) и достаточно строго описывают винтовое движение «степенной» жидкости при наличии теплообмена. Однако их решение может быть получено только численным методом [24, 39, 64]. Для построения модели, допускающей аналитическое решение, сделаем следующие допу­ щения.

1. Течение в направлении оси у существует только в непосред­ ственной близости к стенкам канала. В остальной части сечения канала течение в направлении оси у отсутствует (а1/ = 0).

2.Размеры канала по всей длине постоянны, следовательно значения vx и vz не зависят от г.

3.Температурный градиент в поперечном направлении из-за циркуляционного течения пренебрежимо мал по сравнению с про­

дольным градиентом. Следовательно: dT/dx = dT/dy = 0; dT/dz =£ 0.

Вследствие большого отношения ширины винтового канала к его глубине, учитывая второе допущение, а также особенности циркуляционного течения, можно принять, что на некотором рас­ стоянии от стенок канала скорости vx и vz не зависят от х. Следо­ вательно

дах

О II

dvx

dvu

dvz

 

dvz

(VIII. 7)

И Г * * 0'’

- ^ = 0;

0:

дх " °

ду

 

Если vx — vx(y) И vz = vz(y), то поле скоростей зависит только от координаты у. Поэтому значение квадратичного инва­

рианта - - 12 оказывается равно:

т'-(ж )!+(-згУ

 

 

(vm'8)

Замена частных производных на обыкновенные есть следствие

условия

(VIII. 7). При

этом

из соотношений (VIII. За)

и (VIII. 4)

имеем:

 

 

 

 

р ху — Ла

dvx .

dvz .

 

(VIII. 9)

dy ' Р уг — Ла ‘dy

Ргх = 0

Итак, компоненты напряжений рху и pyz оказываются функ­ циями только у и г , причем последняя зависимость возможна в случае, если T = T ( z ) . Учитывая выражения (VIII. 9), получим вместо уравнений (VIII. 1):

дР

друх

дР

дргу

дР

 

дру2

(VIII. Ю)

дх

ду ’

ду

дг 9

дг

~~

ду

 

При этом интегрирование квазистагических уравнений равновесия выполняется независимо, поскольку число уравнений (VIII. 10) соответствует числу искомых функций.

При записи уравнений энергетического баланса будем считать, что теплопередача за счет теплопроводности вдоль оси канала пренебрежимо мала. В этом случае уравнение (VII 1.5) сведется к виду:

Pcpvx ^

Pxu^

+ Pzy^ .

(VIII. 11)

Представленное в такой форме уравнение (VIII.5), по суще­ ству, превращается в уравнение тепловыделения. При его инте­ грировании нужно удовлетворять температурным условиям на поверхности полимерной струи, учитывая, что температура по­ верхности зависит от теплообмена с окружающей средой.

Можно проинтегрировать уравнение (VIII. 11) подобно тому, как это сделано в работах [45, 46, 51]. В этом случае вместо диф­ ференциального выражения для элементарного объема получим дифференциальное уравнение для элементарного сечения:

dT

dW

 

(VIII. 11a)

Q9ер dz

dz

« i f

 

где dW — работа вязкого трения на участке длиной dz.

Уравнение (VIII. 11а) не является строгим, так как при его получении не учтены условия на границах. Поэтому для учета теп­ лообмена с окружающей средой введем дополнительную характе­ ристику процесса — коэффициент политропичности [26, 64—66] и представим уравнение (VIII. 11а) в виде:

QpcpdT = k d W + QdP

(VIII. 116)

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

Коэффициент политропичности k, входящий в уравнение (VIII. 116), характеризует направление теплового потока и количество тепла, отводимого в окружающую среду или, наоборот, подводимого извне.

VIII. 4. ОСНОВНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ, ОПИСЫВАЮЩИЕ ИЗОТЕРМИЧЕСКИЙ РЕЖИМ ЭКСТРУЗИИ

Распределение напряжений сдвига. В случае изотермического ре­ жима исходная система уравнений равновесия подвергается даль­ нейшим упрощениям. Вследствие неизменности температуры вели­ чина dpzyldz = 0 [см. уравнения (VIII. 4) и (VIII. 9)], и уравнение равновесия сводится к виду:

дР

друх

 

дх

ду

 

дР

0

(VIII. 12)

ду =

дР

дру2

 

dz

ду

 

Из второго уравнения следует, что P = P(xyz). Поскольку ле­ вые части двух других уравнений (VIII. 12) зависят только от х и г, а правые только от ууто

дР

д р ух

= const = С ]

дх

ду

дР

друг

(VIII. 12а)

 

dz ~~ ду

Интегрируя выражения (VIII. 12а) последовательно, получим

для компонентов тензора напряжений выражения:

 

дР

 

Р у * - - д Г У ~ с*

(VIII. 13)

дР

 

P « * = - d 7 y ~ Ci

 

Для функции P(x,z) получим:

дру2

р = - щ - x +

fi (z) = Cl* + /, (г)

(VIII. 13а)

 

друг

 

 

 

Р =

Z +

/ 2 ( X) =

C2Z + f 2 (*)

 

 

“ 3 7 “

 

 

 

Сопоставляя оба выражения для Р и учитывая уравнения

(VIII. 12), найдем:

 

 

 

. м

дР

 

 

 

 

fi {г)==Ж

г = С2г

 

(VIII.

14)

h(x) = ^ - x

= cxx

 

 

 

 

 

Следовательно

 

 

 

п ,

.

дР

, дР

 

(VIII.

15)

Р (х’ г)==- д 1 г + Ж

х = С22 + с'х

 

 

Полученные результаты показывают, что градиенты давления в рассматриваемом изотермическом случае плоского течения по­ стоянны по всему потоку.

Из уравнений (VIII. 13) следует, что напряжения сдвига, дей­ ствующие во взаимно перпендикулярных плоскостях хоу и гоул распределены линейно, т. е. зависят только от у.

Поле скоростей. Для определения поля скоростей для каждого из двух течений необходимо использовать уравнение состояния. С этой целью преобразуем уравнения (VIII. 13) к виду:

! f ( 4 - T

l o ) =

Pzylh

(V III. 16)

(Л -

Лоц) =

p Xy l h

(V III. 16а)

где г)о и г)0ц — безразмерные координаты сечений, в которых

значения напряже­

ний ргУ и рХу равны нулю, соответственно в поступательном

и циркуляционном

течениях

(точки А и В на рис. VIII. 7).

 

Для большей наглядности изобразим пространственную эпюру напряжений сдвига, действующих на верхнюю грань элементар­ ного объема. Такую эпюру можно построить, векторно суммируя компоненты напряжений сдвига ргу и рху, существующие в посту­ пательном и циркуляционном течениях. Типичная пространствен­ ная диаграмма распределения тангенциальных напряжений (рис. VIII. 7) показывает, что в каждом из составляющих течений

 

 

 

существует

сечение нулевых

на­

 

 

 

пряжений. При этом поскольку в

 

 

 

циркуляционном

течении

всегда

 

 

 

существуют области с различным

 

 

 

направлением движения и, следо­

 

 

 

вательно,

градиент скорости

ме­

 

 

 

няет

знак,

сечение,

в

котором

 

 

 

рху — 0, располагается внутри ка­

 

 

 

нала

(0 <с; т)оц <

1)-

 

 

 

 

 

 

Циркуляционное течение в ка­

 

 

 

нале

возникает вследствие того,

 

 

 

что

направление

относительного

 

 

 

движения червяка и корпуса не

 

 

 

совпадает

с осью

винтового

ка­

 

 

 

нала. Поэтому уровень действую­

 

 

 

щих в плоскости хоу напряжений

 

 

 

сдвига зависит от угла подъема

 

 

 

винтового канала, возрастая сего

 

 

 

увеличением. В

случае

ньюто­

 

 

 

новской жидкости взаимное влия­

 

 

 

ние

поступательного

и циркуля­

Рис. VIII. 7.

Типичная

пространственная

ционного течений ограничивается

только этой зависимостью, и

до­

диаграмма

распределения напряжений

сдвига при двумерном

течении.

статочно

определить

граничные

 

 

 

значения

 

компоненты

Ux

=

= f/ sin qp, чтобы рассчитать все параметры циркуляционного те­ чения.

В случае псевдопластичной жидкости эффективная вязкость, определяющая значения скоростей деформации в каждом из тече­ ний [уравнение (VIII. 3)], зависит от квадратичного инварианта тензора скоростей деформаций, в который входят компоненты обоих течений. Поэтому независимое интегрирование выражений (VIII. 16) и (VIII. 16а) для определения функций vx(y) и vz(y), строго говоря, оказывается невозможно.

РАСПРЕДЕЛЕНИЕ СКОРОСТЕЙ В СЛУЧАЕ ЭКСТРУЗИИ НЬЮТО­ НОВСКОЙ ЖИДКОСТИ. ПОСТУПАТЕЛЬНЫЙ ПОТОК

Известно несколько различных способов представления результа­ тов интегрирования уравнения (VIII. 16) для случая ньютонов­ ской жидкости [41—45, 67, 68]. Наибольшее распространение полу­ чила следующая форма решения [1]:

(X, у) = vzd (ху у) + vzP (Ху у)

(VIII. 17)

Выражение (VIII. 17)

показывает, что скорость поступатель­

ного потока в любой точке канала представляет собой алгебраи­ ческую сумму скоростей вынужденного потока vzd и противотока

vzP. Аналитические выражения, описывающие распределение ско­ ростей в плоскости уог, имеют вид:

 

 

оо

 

 

1

sh (gny/w) .

gnx

Vzd =

4Uz

у

 

 

 

 

g

sh (gnh/w) ~

w

 

Я

Z-I

 

 

 

g-U з, 5

 

 

 

 

 

V zp =

1

дР

J

У2

hy

,

4h2

V

ц

dz

\

 

2

2

+

n*

Z ,

 

 

 

 

 

[

 

 

 

 

 

g - 1, 3, 5

 

 

(VIII. 18)

 

 

gny

(

gnw \

h

I

2h )

(VIII. 18a)

 

 

ОБЪЕМНЫЙ РАСХОД ПРИ ПОСТУПАТЕЛЬНОМ ТЕЧЕНИИ НЬЮТОНОВСКОЙ ЖИДКОСТИ

Объемная производительность экструдера (расход поступатель­ ного потока) определяется интегрированием уравнения (VIII. 17) по площади поперечного сечения канала червяка:

w

h

 

Q = $

\ о г ( х , 1 / ) д х д у

(VIII. 19)

о о

Результаты интегрирования принято представлять в форме [I, 44, 45]:

Uzwh

wh3

дР_

(VIII. 20)

2

Fd~ 12т) t p

dz

 

Выражение (VIII. 20) по форме идентично результату, полу­ чающемуся при интегрировании уравнений плоско-параллельного течения ньютоновской жидкости [см. гл. III, уравнение (III. 136)]. Единственное отличие заключается в присутствии безразмерных коэффициентов Fd и Fp, учитывающих тормозящее влияние стенок, приводящее к уменьшению эффективной ширины канала. Оба этих коэффициента определяются относительной шириной канала, характеризуемой отношением h/w.

Коэффициент формы поступательного потока Fd определяется из выражения:

Fd

16

£-1,z3, 5

( gn

h \

(VIII. 21)

n 3h/w

\ 2

w )

 

 

Коэффициент формы противотока Fv определяется выраже­ нием, которое для прямоугольного канала h принимает вид:

FP- I

192h/w

'Г'

_1_ ,,

gn

(VIII. 21а)

я 2

2-1

g 5

2h/w

g -i. з, 5

Значения Fd и Fp, рассчитанные при изменении h/w в интер­ вале от 0 до 2, приведены на рис. VIII. 8.

 

 

Профиль скоростей

поступа­

 

тельного

течения

определяется

 

уравнением (VIII. 17). Если пре­

 

небречь влиянием стенок, то урав­

 

нение (VIII. 17)

сводится к выра­

 

жению:

 

 

 

 

 

Vz =

Uz [(1 - З а ) - | +

За ( ! - ) ']

 

где

 

 

 

 

(VIII. 22)

 

h2 tg ф

дР

 

 

 

Рис. VIII. 8. Коэффициенты формы вынуж­

__

 

 

 

а

QnDNr\

ol

 

 

 

денного потока (кривая 1) и противото-

 

 

 

к\ Fp (кривая 2).

 

ФиЗИЧеСКИЙ

СМЫСЛ

КОЭффи-

 

циента а состоит в том,

что а

это отношение, расхода противотока qp к расходу вынужденного потока qd. В общем случае, если к обоим концам червяка приложе­ ны произвольные внешние давления, а может принимать любые — как положительные, так и отрицательные значения. Однако если внешние давления отсутствуют, значение а изменяется от нуля до единицы. При этом для червяка, работающего в режиме свобод­ ного выхода, а = 0; для червяка, работающего с полностью за­ крытым выходом, но при отсутствии утечек, а = 1; в этих же условиях при наличии утечек а < 1.

Типичные эпюры скоростей, рассчитанные по уравнению (VI11.22), представлены на рис. VIII.9.

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

вотока.

Эпюры

результирующего профиля

скоростей,

соответ­

 

 

 

 

ствующие различным

значениям

 

 

 

 

отношения противотока к расходу

 

 

 

 

вынужденного потока, приведены

 

Поверхность корпуса

в нижней части

рис. VIII. 9.

,,6JB- вынуж-t

/ ^ . П р о т и в о т о к

Профиль скорости циркуляци­

y£-денный p

онного течения. Поле

скоростей в

уЫ

^zp

n бл - роток E

0,2- V za

S

 

 

сечении

канала,

отстоящем на

U

У\ ^

 

 

 

Сердечник червяка

расстоянии

до/2

от

его

стенок,

 

 

 

 

 

Р е зу л ь ти р ую щ и й п о т о к

вследствие

условий

симметрии

1,0 ^

 

 

_а§Ii=1

 

 

определится

выражением, полу­

г » 1

 

 

 

f

ченным

интегрированием

диффе­

 

 

u a=0F_F п-2

ренциального уравнения одномер­

(Свободный

* 4

3

(З акр ы ты й

ного течения [69]:

 

 

Выход)

 

 

вы ход)

 

 

 

 

 

 

Рис. VIII. 9. Типичные

эпюры скоростей

 

 

 

 

 

 

п п г т у п Я Т Р n k i m m т р и о и и п

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]