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

книги / Теоретические основы процессов получения и переработки полимерных материалов

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

Рис. 3.1. Схема течения потока между параллельными пластинами

Обозначим касательное напряжение на пластине как τw. Тогда касательное напряжение на поверхности слоя потока толщиной h

τ

 

= τ

 

dy

.

(3.2)

 

 

 

 

xy

 

w 0,5h

 

Из реологического закона течения определяется градиент скорости потока как функция касательного напряжения:

γ =

dvx

= − f (τxy ).

(3.3)

 

 

dy

 

Поэтому для степенного закона течения Освальда де Виля, которому подчиняются большинство полимерных материалов (см. кривые течения тиксотропных или псевдопластичных сред), можно записать

dv

 

 

τxy 1/ n

 

 

x

=

 

.

(3.4)

 

K

dy

 

 

 

Подставляя значение τxy из уравнения (3.4) в формулу (3.1), получим

31

P

= K

 

v

n

(3.5)

x

 

 

 

x .

 

 

 

y

y

 

Если пластины неподвижны, то после интегрирования выражения (3.5) с учетом (3.2) и (3.3) имеем закон распределения скоростей по координате y:

 

n

 

1

 

P 1/ n h

(n+1) / n

 

(n 1) / n

 

 

vx =

 

 

 

 

 

 

 

 

 

 

y

+

.

(3.6)

n +1

 

 

2

 

 

KL x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На рис. 3.2 показаны эпюры скоростей параллельных слоев потока, построенные на основе уравнения (3.6) при различных значениях показателя n. Видна смена фронтов течения после точки их пересечения по мере увеличения или уменьшения динамического коэффициента вязкости в соответствии с реологическим законом Освальда де Виля.

Рассматриваемый вид течения имеет место, например, в зазорах технологической оснастки при изготовлении полимерных композиционных материалов строительного назначения методом «глухого» прессования. Сходный характер течения наблюдается и при ремонтной «заделке» трещин.

Рис. 3.2. Распределение скоростей при течении между параллельными неподвижными пластинами: 1 – n = 0,5; 2 – n = 1,0; 3 – n = 1,5

Если верхняя пластина движется вдоль оси x с некоторой скоро-

стью v0 и Px ≠ 0 , то из уравнения (3.5) имеем

32

vx

=

 

1

P 1/ n

y1/ n + C

(3.7)

 

 

 

 

y

 

K

x

 

 

или после интегрирования

 

1 P 1/ n

1/ n

 

 

vx =

 

 

y

+ C1 y + C2 ,

(3.8)

 

 

K x

 

 

 

где С1 и С2 – постоянные интегрирования, определяемые из следующих краевых условий: при y = h vx = v0; при y = 0 → vx = 0:

C = v0

 

1

P 1/ n

 

 

 

1

h

 

 

 

 

 

K x

C2 = 0.

 

 

 

 

Подставляя указанные выражения С1 чим окончательное выражение для vx:

vx =

1 P

1/ n

y

 

v

 

 

 

1/ n

+ y

0

 

 

 

 

 

 

 

 

K

x

 

 

 

h

 

 

 

 

 

 

 

 

 

 

h1/ n ;

и С2 в формулу (3.8), полу-

1 ∂P 1/ n

1/ n

 

 

 

 

 

 

h

.

(3.9)

 

 

 

K x

 

 

 

 

 

 

 

 

 

 

На рис. 3.3 представлены эпюры скоростей полимерного потока для различных значений n при vx = v0, когда y = h/2 и ∂P / ∂x ≠ 0.

Видно, что при n > 1 (кривая 1) эпюра скоростей более вытянута по сравнению со случаем n = 1 (кривая 2). При безнапорном течении

( ∂P / ∂x = 0 )

Рис. 3.3. Эпюры скоростей потока при течении Куэтта: 1 – n = < 1; 2 – n = 1; 3 – n = 0

33

полимерного материала (кривая 3) константы С1 и С2 соответственно равны 2v0 / h и v0, а уравнение (3.8) принимает вид vx = v0 y / h; при

этом эпюра скоростей не зависит от реологических свойств исследуемого потока.

Стабилизированное течение полимерных материалов в цилинд-

рической трубе. Как и для случая течения между параллельными пластинами, при движении неньютоновского потока в трубе круглого сечения справедливо уравнение (2.8), которое для ламинарного стабилизированногоиизотермического теченияможнозаписатьтак:

∂τxy

+

∂τ

xz

=

P

,

(3.10)

y

 

x

 

z

 

 

 

если считать, что течение происходит вдоль оси x, а плоскостью поперечного сечения является поверхность y–О–z. Здесь удобно использовать цилиндрическую систему координат. Тогда уравнение (3.10) примет вид

1

(rτ

rx

) =

P .

(3.11)

 

 

r r

 

 

x

 

Если радиус используемой трубы равен R, то, интегрируя зависимость(3.11) впределах 0 r R , получимформулу дляτrx настенке:

τrx =

r

P .

(3.12)

 

 

2 x

 

Подставляя выражение реологического закона Освальда де Виля, которому подчиняютсябольшинствополимерныхматериалов, ввиде

τ = K vx n

rx

r

в уравнение (3.12), можно получить формулу, описывающую течение степенной полимерной жидкости в цилиндрической трубе:

vx

 

r

P 1/ n

 

 

=

 

.

(3.13)

r

2K

 

x

 

 

 

34

 

 

Рис. 3.4. Распределение скоростей по сечению цилиндрической трубы при течении различных полимерных материалов (n < 1) в сравнении с потоком

ньютоновской жидкости (n = 1)

Интегрирование уравнения (3.13) по радиусу в пределах от r до R позволяет определить закон, описывающий распределение скоростей по сечению трубы (рис. 3.4):

 

R

 

r

vx =

r

 

 

2K

 

 

P 1/ n

dr =

n 1 1 ∂P 1/ n

 

 

 

 

 

 

 

 

 

 

 

 

2K

x

x

 

 

 

 

 

n

 

 

 

nR

 

R

 

 

P 1/ n

 

 

r (

=

 

 

 

 

 

 

 

 

1

 

 

 

n +1

 

 

 

 

 

2K x

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(R(n+1)/ n r(n+1) / n ) =

n+1)/ n

(3.14)

.

 

 

 

 

 

При этом максимальная скорость течения по оси трубы (r = 0)

vxmax

=

nR

 

R

P 1/ n

(3.15)

 

 

 

 

.

n +1

2K

 

 

 

x

 

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

R R

Q = 2πrvxdr = 2πrvx max

0 0

1−

r (R

n+1)/ n

 

n +1

2

 

 

dr =

 

 

πR

vxmax

. (3.16)

 

 

 

3n +1

 

 

 

 

 

 

 

 

 

 

35

Расход Q определяется как произведение средней скорости течения vxmid на площадь поперечного сечения πR2 :

 

 

Q = vx

πR2.

(3.17)

 

 

 

 

mid

 

 

 

 

 

Заметим также, что массовый расход потока (М) равен произ-

ведению ρQ .

 

 

 

 

 

 

 

 

Путем комбинации формул (3.16) и (3.17) легко определить за-

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

 

 

vx

 

=

n +1

 

vx .

(3.18)

 

 

3n +1

 

max

 

mid

 

 

В частности, при n

=

1

(ньютоновский поток)

имеем

vx

= 0,5vx .

 

 

 

 

 

 

 

mid

max

 

 

 

 

 

 

 

 

Учитывая зависимость (3.18), уравнение (3.14) можно записать

еще и так:

v

 

=

n +1

v

1

 

r

(n+1) / n

x

 

 

 

 

 

 

3n +1

x mid

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

. (3.19)

Подставляя в выражение (3.18) значение максимальной скорости из уравнения (3.15) и проводя алгебраические преобразования, можно получитьзависимостьмежду расходомиперепадомдавлениявтрубе:

Q =

πn(3n +1)

 

 

 

(n +1)

2

 

 

 

где P → ∆P = (P1 P2 ) , а x = ∆L

1 P 1/ n

 

3n 1 / n

 

 

 

 

R(

+ )

,

(3.20)

 

2K x

 

 

 

 

= (L1 L2 ).

Из формулы (3.19) следует, что при течении в цилиндрической трубе полимерного материала, реологическое поведение которого описывается законом Освальда де Виля, его объемный расход пропорционален перепаду давления P / x в степени 1/ n и радиусу трубы в степени (3n 1) / n .

Описанные зависимости используются в инженерных расчетах при выборе необходимого оборудования и разработке технологических процессов получения новых полимерных композиционных материалов (табл. 3.1).

36

Таблица 3 . 1

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

Освальда де Виля

Формула для Величина Расчетная формула ньютоновской

жидкости

Напряжение сдвига

 

 

 

 

 

 

 

 

 

 

 

τw =

 

R P

 

 

 

 

 

 

 

 

 

 

 

τw

 

=

R

 

P

 

 

 

 

 

 

 

у стенки трубы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 ∂x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 ∂x

 

 

 

 

 

 

 

Напряжение сдвига

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

в потоке

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

τ = τw

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

τ = τw

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

Скорость потока

 

vxmax

 

=

 

 

 

 

n

 

 

 

 

 

1 P 1/ n

 

 

 

vxmax

=

 

R2

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 ∂x

 

 

 

 

 

в центре трубы

 

 

 

 

 

+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

2K x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R(n+1) / n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Скорость в

 

 

v

 

= v

 

 

 

 

1−

 

r (n1) / n

 

 

v

 

= v

 

 

 

 

1−

r

2

любой точке потока

 

 

 

x

 

 

 

 

 

xmax

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

xmax

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Средняя скорость

 

 

 

vxmid

 

= vxmax

 

 

n +1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

потока в трубе

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

xmid

= v

xmax

0,5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3n +1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vx

 

 

n +1

 

 

1/ n

 

 

 

 

dv

x

 

 

 

2vx

 

 

 

r

 

 

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

 

 

dv

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

max

 

 

 

 

 

сдвига потока

γ

=

 

 

x

 

=

 

 

 

max

 

 

 

 

 

 

 

 

 

 

 

 

γ =

 

dr

=

 

 

 

 

 

R

 

 

 

R

 

 

dr

 

 

 

 

 

 

 

 

 

n

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dvx

 

 

3n +1 Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

γw =

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dvx

 

 

 

 

 

4Q

 

 

 

 

dr

 

 

 

 

 

n

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

 

 

 

 

 

πR

 

 

 

 

 

 

 

 

 

 

 

 

сдвига потока

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γw

=

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dr

πR

3

 

у стенки трубы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

=

 

πR3n

R P 1/ n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Объемный расход

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

πR

4

P

 

 

 

 

 

 

 

 

 

3n

+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q =

 

 

 

 

 

 

 

полимерного потока

 

 

 

 

 

 

 

 

 

2K x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

37

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

Рис. 3.5. Двумерная область среды, выделенная как комплекс треугольных элементов

На рис. 3.5 показан пример разбиения на треугольные элементы некоторой двумерной области, например полимерного потока в канале произвольного сечения. Каждый элемент может иметь свои собственные размер и форму, которые выбирают так, чтобы они наилучшим образом соответствовали форме и размерам исследуемого сечения канала или элемента оборудования.

Метод МКЭ отличается от метода конечных разностей, при котором используется сетка с ячейками одинакового размера, описываемыми теми же координатами, что и геометрия сечения. Точки пересечения кривых, ограничивающих соседние элементы, называются узлами. Значения переменных, вычисленные в узлах, дают искомое решение.

38

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

Выбранные аппроксимирующие функции называются пробными функциями или пространственными изменяемыми моделями. Обычно в качестве таких функций используют усеченные полиномы. Число членов (коэффициентов) в таком полиноме должно быть по крайней мере равно числу степеней свободы, присущих каждому отдельному элементу.

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

39

Лекция IV

Турбулентное течение полимерных материалов

Определение критерия Рейнольдса. Разнообразие в реологиче-

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

Число Рейнольдса – это средняя скорость (vmid) потока полимерного материала в канале, умноженная на характерный размер канала (L) и деленная на динамический коэффициент вязкости (η):

Re =

vmid L

.

(4.1)

 

 

η

 

В качестве характерного размера можно использовать радиус (диаметр) трубы, если изучается течение в цилиндрической трубе. Если же поперечное сечение канала имеет произвольную геометрию, то за характерный размер может быть принят приведенный (эквивалентный) радиус Requ или гидравлический радиус Rr. Его величина определяется как радиус круга, имеющего такую же площадь исследуемого сечения канала. При этом гидравлический радиус представляет собой отношение площади поперечного сечения канала S к его периметру P: Rr = S/P. Заметим также, что гидравлический радиус круга в два раза меньше геометрического:

Rr = S = πR2 = R . P 2πR 2

Если известна функция реологического закона τ = f (τ) для ис-

следуемого полимерного материала, то в знаменателе формулы (4.1) можно записать значение так называемой эффективной средней вязкости (ηeffmid ) , определяемой как тангенс угла наклона (α) прямой,

соединяющей точку O с точкой А, соответствующей рассматриваемому градиенту скорости сдвига потока γA (рис. 4.1). На практике приходится пользоваться такими формами записи критерия Рей-

40