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

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

..pdf
Скачиваний:
49
Добавлен:
15.11.2022
Размер:
21.7 Mб
Скачать

где Ux = nDN sin q>; с — отношение расхода про­ тивотока к расходу вынужденного потока в цир­ куляционном течении, равное

с = Т Т Г - ¥ -

(VIII. 24)

ШхЦ дх

 

Рис. VIII. 10. Распределение скоростей циркуляционного те­ чения.

Если пренебречь утечками, суммар­ ный расход в направлении оси х будет равен нулю, так как с обеих сторон поток

ограничен стенками канала.

Поэтому

с = 1, и уравнение (VIII. 22)

принимает

вид:

 

f* = t / * [ 3 ( f - ) 2 - 2 | - ] (VIII. 25)

Из уравнения (VIII. 25) видно, что эпюра скоростей циркуляци­ онного течения (рис. VIII. 10) не зависит от давления в головке,

а полностью определяется размерами канала и скоростью враще­ ния червяка.

Результирующее распределение скоростей. Вектор истинной скорости в трехмерном потоке определится как векторная сумма

компонент вектора

скорости в данной точке

потока [I]

(рис. VIII. 11). (На

рис. VIII. 11,6 эти же диаграммы

распределе­

ния скоростей изображены в перспективе [48]).

 

Значения поступательной (параллельной оси z ) и циркуляцион­ ной (параллельной оси х) компонент вектора скорости рассчитаны для различных значений отношения yjh при режиме свободного

выхода

(а = 0), закрытого выхода, или нулевого расхода (а =

1)

и для

промежуточного

режима

(а = 0,5). Эти компоненты пока­

заны пунктиром на рис. VIII. 11,а. Векторы, изображающие фак­

тическую скорость в

каждой

точке

потока, показаны

на

рис. VIII. 11,6 сплошными стрелками. Из

этих диаграмм видно,

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

Отметим, что расположенная в плоскости 1оу компонента ско­ рости Vi всегда положительна при любых значениях а. Следова­ тельно, весь находящийся в червяке материал всегда движется вдоль его оси в направлении к головке. Исключением является режим нулевого расхода, при котором элементы находящейся в канале среды будут совершать круговое движение в плоскости, перпендикулярной оси /, совершенно не перемещаясь. При этом слой материала, расположенный на расстоянии 2/з высоты канала от его дна, остается неподвижным.

Если а = 0 (полностью

открытый

выход), находящийся

в ка­

нале расплав

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

(см.

рис. VIII. 4).

Увеличение а

приводит

к уменьшению шага

этой

Рис. VIII. 11. Распределение скоростей в винтовом канале

червяка:

а —векторная сумма скоростей; материал движется вдоль

оси/; б —пространственная диа­

грамма распределения скоростей в канале червяка (Л—поверхность корпуса; В —сердечник

червяка).

 

спирали и продлению пребывания расплава в червяке. В предель­

ном случае, когда а == 1, шаг спирали становится равным нулю, и

поступательное движение расплава вдоль оси I прекращается.

ПОЛБ СКОРОСТЕЙ И ДИФФЕРЕНЦИАЛЬНОЕ УРАВНЕНИЕ ОБЪЕМ­ НОЙ ПРОИЗВОДИТЕЛЬНОСТИ ЧЕРВЯКА ПРИ ЭКСТРУЗИИ АНО­ МАЛЬНО-ВЯЗКОЙ ЖИДКОСТИ

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

(d v x / d y )г I / d v z \ ( d v j d y У J V d y )

(VIII. 26)

Поскольку напряжения сдвига определяются выражением вида (VIII. 4), можно заменить отношение скоростей деформаций, вхо­ дящее в уравнение (VIII.20), отношением напряжений. Действи­ тельно, с одной стороны, из уравнений (VIII. 16) и (VIII. 16а) по­ лучается:

Рху

(дР!дх) (и — Пои)

Рги = (дР/dz) (и — По)

С дру(юй стороны, из уравнений (VIII. 3) и (VIII. 4) по опреде­ лению следует:

1 л

Рхи = Ц [у Ь (Л)]

dvx/dy

л и . 28)

 

 

1-л

 

 

Ргу = Ц [у Ь (Л)]

dv,/dy

(VIII. 28а)

Таким образом

 

 

 

Рху

dvx/dy

 

 

(VIII. 29)

Pzy

dVg/dy

 

 

 

 

 

Заменяя в выражении (VIII. 26) отношение скоростей дефор­ мации отношением напряжений сдвига, получим две новые формы записи квадратичного инварианта:

- х 2

(Л - Лоц)21

( d vz V

г

1 + v 2

(Л — Ло)2

1 ( d v x \ 2

(VIII. 30)

(Л — Ло)г J

\ dy )

L

(Л — Лоц)2

J \ dy )

где X и v — отношения градиентов давлений, действующих в циркуляционном и продольном течениях, соответственно;

_

дР/дх

я

щ_ дР/дг

(VIII. 30а)

г

дР/дг

V ~~ дР/дх

 

Выше было показано, что градиенты давлений постоянны, а скорости vz и vx могут считаться функциями только одной пере­ менной у. Поэтому решение задачи приводится к квадратурам.

Подставляя в уравнение (VIII. 16а) значения напряжений из соотношения ( (VIII. 3) и учитывая, что V2I2 в уравнении (VIII. 4) вьфажается из равенства (VIII. 30), получим после интегрирова­ ния два уравнения, описывающих распределение скоростей в по­ ступательном vz (т]) и циркуляционном их(т|) течениях:

Г)

л-1

М л ) = $ [(Л-Ло)2 + х 2(л-Лоц)21 2

(л-Ло)^Л

(V III. 31)

 

о

 

 

 

 

 

 

 

 

л-_1

 

(л) = X ^[(л — Ло)2 + X2 (л — Лоц)2]

2

(VIII. 31а)

2 (л — Лоц) <*л

где Vz и Vx — безразмерные скорости;

 

 

 

б*(Л) =

VzV"

Vx (л) =

 

Vx\l

 

п hn+1

 

П г.Л+1 ‘

 

 

(dP/dz)n к

 

(дР/дг)п К

 

Введем дополнительные обозначения:

 

 

 

л -1

 

 

 

{ (fl) =

[(л — Ло)2 + X2 — Лоц)2] 2

 

 

(VIII. 32)

 

 

л+1

 

 

Ф (л) =

[(л Ло)2 + X2 Лоц)2] 2

(VIII. 32а)

ац = ^

[(л — Ло) + X2(л — Лоц)] (Я + 1)

(VIII. 326)

Тогда, интегрируя выражения (VIII. 31) и (VIII. 32а), получим:

vz (л) =

0 ( T J ) — Ф ( 0 )

%vx (л)

(VIII. 33)

п + 1

 

б* (Л) ( ^

+ * ) =

Ф (Пя + -f -(°-) ~ (Лоц ~ Ло) \ I (Л) *1

(VIII. 33а)

 

 

 

0

 

 

 

 

Л

 

вг(л) (х г + 1 ) =

*хЧ«+Ф1)0) + (Т,01Г_^ 5 f (n)dr]

(VIIL33б)

 

 

 

О

 

По аналогии с выражением (III. 125), введем безразмерный градиент давлений в поступательном течении (Вг) и в циркуля­ ционном (В*):

, __ ( д Р/ д г ) я

 

 

_ д Р / д х

 

 

 

 

 

 

(VIII. ЗЗв)

2 ( д Р / д г ) * ’

 

*

( д Р / д х )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

fdP\* _ „ Г (п + 1) UzVtn.

/дР\* _

м Г (” +1)0* I 1'"

(VIII. ЗЗг)

VdzJ

Ч

hn+[

J

 

 

Ч

hn+l

J

 

 

 

 

 

Из

уравнения

(VIII. 33) следует,

что между

 

величиной % и

параметрами г)оц и ло существует однозначная зависимость:

В\

[0 (1 ) - 0

(0)]-1

 

 

 

 

 

 

 

(VIII. 34)

 

tg<p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

 

 

п+\

 

 

п+1

} \(п

В г { [(1 -

Ч э ) У

+ (1 - Чоц)2]

2 ~

( Г\У + л2ц)

2

}

(VIII. 34а)

 

 

 

 

 

(V +

t g ф ) 1^

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

Ф = a r c t g (Ux/Uz)-

 

 

 

 

 

 

 

 

 

 

Параметры

т|0

и т)оц могут быть определены из уравнений

(VIII. 33а) и (VIII.33)

при г| =

1 из условия

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

^

(л)

= О

 

 

 

 

 

 

 

 

(VIII. 35)

о

 

 

 

 

 

 

 

 

 

 

 

 

которое означает, что объемный расход циркуляционного течения равен нулю (утечками через зазор пренебрегаем).

Интегрируя по г| уравнения (VIII. 33) и (VIII. 336) и учитывая условие (VIII. 35), получим два выражения для определения нор­ мированного расхода поступательного течения:

I

Q«— [Ф (n) - ♦ (0)1 *1 (Viii. 36)

0

1 Ч

Qz =

(Лоц — Ло) ^

(VIII. 36а)

где

о о

 

 

 

Qz =

QzfiPzhB* (п + 1)]

(VIII. 366)

Построение математической модели зоны дозирования состоит в использовании результатов интегрирования этих двух уравнений и определении из граничных условий (VIII. 16) значений парамет­ ров т]о, г]оц и v Основное возникающее при этом затруднение свя­ зано с тем, что квадратуры в выражениях vz и vx не удается найти в общем случае; это связано с тем, что при вычислении выражений (VIII. 31) и (VIII. 31а) или (VIII. 33а) и (VIII. 336) необходимо разложить в ряд подынтегральное выражение и после этого вы­ полнить почленное интегрирование.

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

муле

(VIII. 36)

для

п = 2—8

при

изменении

г)0

в диапазоне

- 2 0 ^

т)о^

0,3

с

результатами,

полученными

по формуле

одно­

мерного течения.

 

 

 

 

 

 

 

 

Из таблицы видно, что максимальные отклонения наблюдаются

при т]0 = 0,2—0,3

и

достигают 15—20%. В области значений

— 20 ^

т]0 ^

— 0,3,

 

представляющей

наибольший

интерес

для

практических расчетов, ошибка не превышает 8%. Учитывая, что многие применяемые в расчете данные определены с аналогичной точностью, можно согласиться с таким приближением, так как его использование существенно упрощает процедуру расчета внешних характеристик червяка.

В этой же таблице приведены значения цоц и х> соответствую­ щие различным г]0. Сравнивая их со значениями т]оц для одномер­ ного течения (рис. VIII. 12), видим, что в двумерном течении они несколько больше и слабо зависят от т\0. Очевидно, что безразмер­ ный градиент давления в двумерном течении будет несколько меньше, чем в одномерном, при одинаковом расходе. Аналогичным образом будут различаться и напряжения сдвига, так как в дву-

Сопоставление результ ат ов точного и приближенного расчет а

 

 

п= 2

 

 

 

 

71=4

Но

X

v ° 3 QJ K

 

v ° 3 X

 

—20,0

17,40

374

0,99

1,95

375

15,80

0,92

—5,0

8,74

388

1,00

1,88

388

11,80

0,94

—3,0

3,57

388

1,01

1,82

399

3,80

0,96

-1 ,0

2,87

386

1,03

1,63

427

2,86

0,97

-0,5

2,00

386

1,06

1,46

443

2,48

1,07

-0 ,3

1,81

385

1,07

1,35

450

2,29

1,08

-0,2

1,59

385

1,08

1,28

454

2,09

1,09

-0,1

1,28

385

1,09

1,20

457

1,80

1,11

0,05

0,98

385

1,09

1,04

462

1,40

1.13

0,2

0,76

385

1,09

0,82

459

1,12

1,14

0,3

0,65

385

1,16

0,52

447

1,08

1,16

 

 

 

71= 6

 

 

 

71= 8

 

v z

оЯ'О' X

Ч г К

v z

v ° 3

X

* г К

 

 

 

 

 

 

 

 

1

2,67

392

30,80

0,86

3,27

406

17,40

0,99

1,96

2,48

410

7,80

0,90

3,00

427

8,70

1,00

1,88

2,34

423

4,90

0,92

2,71

442

3,60

1,01

1.83

1,88

457

2,70

1,02

1,98

390

2,90

1,04

1,63

1,57

385

2,00

1,05

1,60

390

2,00

1,06

1,46

1,41

385

1,83

1,09

1,42

390

1,81

1,07

1,35

1,31

385

1,66

1,12

1,34

390

1,59

1,08

1,28

1,22

385

1,40

1,14

1,25

390

1,28

1,09

1,20

1,07

385

1,14

1,19

1,13

390

0,98

1.10

1,04

0,91

385

0,98

1,25

0,99

390

0,76

1.11

0.82

0,73

385

0,91

1,29

0,88

390

0,65

1,18

0,52

мерном течении эффективная вяз­ кость меньше, чем в аналогичных условиях в одномерном течении. В результате расчета необходимо определить такие характеристики двумерного течения, которые представляют интерес для анали­ за процесса экструзии. К их чис­ лу относятся: объемный расход поступательного течения, опреде­ ляющий производительность про­ цесса; напряжения сдвига в по­ ступательном и циркуляционном течениях, действующие на внут­ ренней поверхности корпуса и оп­ ределяющие величину диссипируемой в расплаве энергии; гра­ диенты давлений в поступатель­

ном и циркуляционном течениях, определяющие давление и осевое усилие, развивающиеся в пределах зоны дозирования.

Рассмотрим все эти задачи последовательно. В качестве основ­ ного параметра расчета воспользуемся величиной т]0. Теоретиче­ ски значение т]о может лежать внутри одной из двух областей:

QZ < 2 /(/Z + 2);

— оо<110< 0>5

(VIII. 37)

Q z> 2/(n + 2);

+ <x>>r\0 ^ O t5

(VIII. 37а)

При этом условие (VIII. 37) соответствует положительному гра­ диенту давлений, а условие (VIII. 37а) —отрицательному.

Безразмерный градиент давлений в поступательном течении можно было бы определить из выражения (VIII. 34):

Д г = (1 +

х t g ф ) 1/л/ [ Ф ( 1 ) — Ф

( 0 ) ] ,/Л

 

 

(VIII. 3 8 )

Однако

предварительно

необходимо

рассчитать

величину

Для

этого

приходится применять

метод

последовательных приближе­

ний.

Вначале рассчитывают

величины В'х и 5'

как функции т]'

и ^оц по формулам одномерного течения (III. 130) и (III. 131). Значение %в первом приближении определяется выражением:

Xl = 1иЧ>ипВ'х/в'г

(VIII. 39)

Полученное из выражения

(VIII. 39) значение xi подставляется

в соотношение (VIII. 34), которое интегрируется для вычисления

У\штеНН0Г0 значения Х> которое затем подставляют в выражение (VIII. 38) для определения из последнего В2. Величина Вх рассчи­

тывается из очевидного соотношения

= B2x/tg ф1/л

(VIII. 40)

У Заказ 483

257

Напряжение сдвига на стенке корпуса экструдера определяется из выражения (VIII. 16) с уч _ом соотношений (VIII.ЗЗв) и (VIII. ЗЗг):

Ргу =

р Д « [ ^ -

й- )-^ £] '/И (1 -

По)

(VIII. 41)

Рху =

рВх [ {п

(1 _

%ц)

(VIII. 41а)

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

Циркуляционное течение. Основным критерием, используемым для определения параметров циркуляционного течения, является условие нулевого расхода (VIII. 35).

Выше отмечалось, что используя это условие можно определить численным методом значения х\0ц и % как функцию т)0. Пример расчетов такого рода приведен в таблице VIII. 1.

Более простой способ — это определение величины г]0ц как функции индекса течения из уравнения (III. 133). Результирующее выражение имеет вид:

(1 - <S+2+ (Чэц)"+2 - (я +

2) ( i u) n+1 = О

(VIII. 42)

где т|0Ц— координата сечения

нулевых напряжений сдвига, определенная для

модельного течения в плоскости хоу.

 

Уравнение (VIII. 42) позволяет определить величину

т)'ц, кото­

рая оказывается однозначной функцией индекса течения:

 

Поц={1 + [(« + 2)/т1оц-1]1/П+2}

(VIII. 43)

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

п—1

этого разложим величину [1 + — По)2 v /(л — Лоц)2] 2 в ряд Тэйлора относительно т) = 1:

(Л — Ло)д УП

п-\

(1 — Ло>2 У2

2

(Л — Лш' J

 

U Tlcu)J

Ч- <n2,--F"(1) +

тУ~J + h - l) F ' (l ) +

(VIII. 44)

F'(l) = (« —1)(1 —Лои)п 3 [l + (1—Ло)°v?/(l—Лоц)2] 2 i(l + v?)ri + т1„ц-т|о]<

258

Пренебрегая членами порядка (т) — I)2, преобразуем выражение (VIII. 31а) к виду.:

ч

 

 

 

 

 

 

 

 

 

 

(Л ~ 1ц)п+1 -

Л?ц+1 +

$ - Пои)" [F (1) +

(Л -

1) F' (1)1 *1 -

IF (О -

F' (1)1

 

 

 

 

 

 

 

 

 

 

 

я + 1

 

 

(Л - Л 0ц)'1+2-Л ?ц 2

. Лоц[(Л-Лоц)п+1-Л?ц+1]

(VIII. 45)

+ F' (1)

 

п + 2

 

 

 

 

п

+ 1

+

 

 

 

 

 

 

 

 

 

 

Для определения т]0ц проинтегрируем уравнение (VIII.45) в

пределах от г\ =

0 до r\ =

1 и приравняем результат нулю:

 

F (1) — F' (1)

(1 -

Лоц)л+2 +

Л?ц+2 ~ ( я

+ 2)

Я+1

 

 

 

п + 1

 

 

 

 

п + 2

 

 

 

 

 

+

ЛзцР'О)

( 1 - Л оц)"+2 +

Л?ц+ 2 - ( «

+

2)ЛоПц+1

 

 

п + 1

 

 

 

п + 2

 

 

^

 

 

 

F ' (1)

(1 -

Лоц)п+3 +

Лоц+3 - (« +

3) Л о Г

п

 

(VIII. 46)

 

---:—тг"*------------------- ;—к------------------ —U

 

 

п + 2

 

 

 

ti

3

 

 

 

 

 

Первый и второй члены уравнения (VIII. 46) совпадают с выра­ жением (VIII.42), полученным в результате интегрирования урав­ нения, описывающего поле скоростей одномерного моделирующего течения. Следовательно, если T]O4 определяется из уравнения (VIII. 43), эти члены обращаются в нуль. Третий член уравнения (VIII. 46) представляет собой, по существу, выражение (VIII. 42), вычисленное при значении индекса течения, увеличенного на еди­ ницу. Следовательно, если пренебречь разницей между соседними значениями г)оц, которая с увеличением п стремится к нулю (см. рис. VIII. 12), то можно принять, что и этот последний член равен нулю. Поэтому для определения координаты т]оц реального течения будут использованы результаты, полученные при замене циркуля­ ционной составляющей двумерного течения одномерным течением с одинаковыми кинематическими характеристиками (h и U sin ф) и нулевым значением расхода. Величина градиента давлений модели­ рующего циркуляционного течения при этом равна:

 

 

(VIII. 47)

где

 

 

( п +

l)1/rt

(VIII. 48)

Е888

 

[ ( I - V T ' - C ' r

Для удобства расчетов результаты численного решения уравне­ ний (VIII. 43) и (VIII. 48) представлены в виде графиков (рис. VIII. 12).

Полученное приближенное решение позволяет рассчитать объемный расход и градиенты давлений в продольном и поперечном направ­ лениях для рассматриваемого случая плоского моделирующего те­ чения. Объемный расход определится выражением:

Q = h w F d $ о 2 ( т 1 )Л |« • ^ р ^ - Ч г (Ло)

(V III. 49)

О

где ф(т]0) находится из соотношений (III.'133) —(III. 134).

Если объемный расход известен и необходимо рассчитать все остальные параметры процесса, то выражение (VIII. 38) рассматри­ вается как трансцендентное уравнение, решая которое находят соответствующее значение rjo. Коэффициент Fd вводится для учета тормозящего влияния стенок. Принимаем, что численное значение для аномально-вязких жидкостей не отличается от значения, полу­ ченного при интегрировании уравнений ньютоновских жидкостей Г66].

Прежде чем перейти к определению параметров процесса экс­ трузии, выразим скорости Ux и Uz через характеристики червяка и частотой его вращения:

Uz = n D M cos qp;

UX = я D N sin ф

(V III. 50)

а> =

atD (1 — j - )

sin cp/i

(V III. 51)

f =

nD tgq>

 

(V III. 52)

Уравнение производительности при экструзии псевдопластичной жидкости приобретает вид:

Q

F d N

* ы

 

(VIII. 53)

п + 2

 

где

а — коэффициент подачи, численно равный половине объема

одного витка

канала;

 

 

 

а я® 0,5n 2D 2h (1 — ie/t) sin ф cos ф

(V III. 54)

 

Если расплав обладает свойствами ньютоновской жидкости,

уравнение производительности принимает вид:

 

Q =

a F d N

 

dP

(V III. 55)

ц

dl

где

 

 

 

 

 

 

Р =

n D h 3 (1 — ie/t) sin2 ф

(V III. 56)

 

 

12

 

 

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

2G0

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