Например, Бобцов

Течение неньютоновской жидкости в шнековом прессе

УДК 664.6.
Течение неньютоновской жидкости в шнековом прессе
Андреев А.Н.
andreevanatoly@yandex.ru
Санкт-Петербургский национальный университет информационных технологий, механики и оптики
Институт холода и биотехнологий
В работе для вывода дифференциальных уравнений течения и уравнения объемного расхода для канала шнекового пресса использована гидродинамическая теория неньютоновской жидкости которая позволяет установить количественные зависимости между геометрическими характеристиками рабочего пространства, свойствами материала и режимом обработки. Ключевые слова: неньютоновская жидкость, тесто соломки, шнековый пресс.
К группе хлебобулочных изделий пониженной влажности относится соломка, в виде округлых палочек диаметром 8 мм и длиной от 10 до 28 см., которая вырабатывается комплексно-механизированных линиях [1,2]. Специфическим процессом в производстве соломки является формование жгута на шнековом прессе. Вопросом течения вязкопластических масс в шнеке занимались многие ученые: Арет В.А., Мачихин Ю.А., Берман Г.К., Ворожцов Л.А., Азаров Б.М., Шенкель Г., Штруб Р., Сквайрс Р., Галт У. и др. [3,4,5].
В работе поставлена задача: рассмотреть течение неньютоновской жидкости к которой относится тесто соломки в условиях двумерного сжатия в шнековом прессе с целью разработки методики расчета и усовершенствования конструкции с учетом ресурсосбережения.
Для вывода дифференциальных уравнений течения и уравнения объемного расхода для канала шнека можно воспользоваться уравнениями движения вязкой несжимаемой жидкости (рис.1).
Рис. 1 Сечения канала шнека с выбранным направлением осей координат

Уравнение движения для стационарного ламинарного течения несжимаемой

изотропной жидкости вдоль оси винтового канала шнека имеет вид:

∂P ∂Z

=

µ

(

∂ 2ϑ ∂x

Z 2

+

∂ 2ϑ Z ∂ y2

)

+

∂ϑ Z ∂x



∂µ ∂x

+

∂ϑZ ∂y



∂µ ∂y

(1)

В этом уравнении вязкость может зависеть от положения элемента

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

скорости в этой точке канала. Если считать, что вязкость жидкости очень мало

изменяется вдоль оси в поперечном сечении канала (так как температура

жидкости в канале и величина градиента скорости при переходе от одной

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

можно приравнять нулю, тогда уравнение принимает вид:

∂P ∂Z

=

µ

(

∂ 2ϑ Z ∂ x2

+

∂ 2ϑ ∂y

Z 2

)

+

∂ϑ Z ∂y



∂ ∂

µ y

(2)

Уравнение (2) описывает неизотермическое течение неньютоновской

жидкости в канале шнека для всех практических случаев. Общее аналитическое

решение этого уравнения весьма сложно и до настоящего времени не получено.

Рядом авторов [3,4,5] получены частные решения уравнения, описывающие

изотермический и адиабатический режимы работы для шнеков с заданной

геометрией винтового канала при выпрессовывании различных материалов с

определенными реологическими свойствами.

Если предположить, что вязкость жидкости по высоте канала не меняется,

∂µ

то последний член у уравнения 2 исчезает, так как ∂ y равно нулю. При этом

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

учетом этого ограничения получается следующее уравнение:

∂P ∂Z

=

µ

(



2ϑ Z ∂ x2

+

∂ 2ϑ Z ∂ y2

)

(3)

Если в канале шнека существует только вынужденный поток, то

уравнение (3) должно удовлетворять следующим граничным уравнениям

ϑ Я = 0 при х=0 и х=W, y=0

ϑ Я = U Z при y=h Поскольку в данном случае градиент давления по оси канала равен нулю,

то уравнение (3) превращается в однородное дифференциальное уравнение:

∂P ∂Z

=

µ

∂ 2ϑ Z ∂ y2

(4)

Решение этого дифференциального уравнения, удовлетворяющее

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

независимо друг от друга в различных выражениях рядом авторов [3,4,5,].

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

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

введением в основное одномерное уравнение поправочных коэффициентов, на

которые соответственно умножились члены, определяющие величину

объемного расхода вынужденного потока Qd и величину расхода противотока Qp . Уравнение объемного расхода вынужденного потока может быть записано в виде:



=

n⋅UZ ⋅ ω 2

⋅ h ⋅ Fα

(5)

где n - число заходов шнека; UZ - скорость движения жидкости вдоль оси; ω - ширина винтового канала шнека; h - глубина канала шнека; Fα -

коэффициент форм для расхода вынужденного потока;

Величина Fα зависит только от геометрии канала шнека и полностью

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

h . Уравнение ω в результате

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

уравнения (5), удовлетворяющим следующим граничным условием:

ϑ Я = 0 при x=0; x=ω , y=0 и y=h

Решение, удовлетворяющее этим граничным условиям и описывающее

объемный расход противотока, выражается следующим уравнением:

Qp

=



n

⋅ω ⋅ 12µ

h3

(

∂ ∂

P Z

)



FP

(6)

где µ - вязкость жидкости; FP - коэффициент формы для расхода

противотока,

∂P ∂Z

- градиент давления вдоль оси.

Результирующий поток, возникающий вследствие наложения

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

уравнения 5 и 6.

Q = QB + Qn

(7)

Q=

n⋅UZ ⋅ ω 2

⋅ h ⋅ FB −

n

⋅ω ⋅ 12µ

h3



∂ ∂

P Z

(8)

Учитывая, что:

UZ = π ⋅ D ⋅ N ⋅ cosϕ

(9)

Z

=

e sinϕ

(10)

где D - диаметр шнека; N - число оборотов шнека в единицу времени; ϕ угол подъема винтовой нарезки шнека; e - ширина гребня канала шнека.

Уравнение (8) может быть записано в виде:

-

Q=

Fb ⋅ n ⋅ π

⋅ D⋅ N ⋅h⋅ω 2

⋅ cosϕ



F n⋅n ⋅ h ⋅ ω ⋅ sinϕ 12 ⋅ µ



∂P ∂e

(11)

или

Q=

Fb ⋅ α

⋅N−

Fn



β µ



∂P ∂e

(12)

где α , β , Fb и Fn зависят только от геометрических размеров шнека. В данном случае:

α

=

Fb ⋅ n ⋅ π ⋅ D ⋅ h ⋅ ω ⋅ cosϕ 2

(13)

β

=

Fn ⋅ n ⋅ h3 ⋅ ω 12

⋅ sinϕ

(14)

Следует отметить, что уравнение (12) широко используется для

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

пластических масс, резиновых смесей, пищевых продуктов [6]. Однако, при

выводе уравнения (12) , авторы делали допущение о том, что вязкость

материала в поперечном сечении шнека остается постоянной. В прессе с

глубоким каналом (как это имеет место при формовании соломки) по глубине

канала шнека могут иметь место градиенты температуры в перерабатываемом

материале и различные скорости сдвига на внешнем и внутреннем диаметре

шнека. Величину среднего градиента γ скорости в канале можно определить

по формуле:

γ

=

π

⋅ D⋅ N h

(15)

где h - глубина канала шнека.

Кроме того тесто соломки как вязко-пластичный материал, заполняющий

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

так называемую «эластичную пробку», которая равномерно продвигается вдоль

оси канала. Наличие «пробки» изменяет условие продвижения материала по

каналу шнека, ухудшает гомогенизацию, способствует росту градиента

температуры и скорости сдвига. Эти факторы, как показали многочисленные

исследования, проведенные для различных вязко-пластичных систем, в том

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

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

влияния их на изменение вязкости обрабатываемого материала.

До настоящего времени течение теста с переменной по глубине канала

вязкостью не было изучено. Впервые решение этой задачи дал Р. Штруб для

случаев линейного и экспотенциального изменения термопластичного

материала по глубине канала червячной шприц-машины. Однако Р Штруб,

интегрируя

уравнение

Навье-Стокса,

пренебрег

членом

(



ϑЯ ∂y



∂µ ∂y

),

который

появляется, если воспользоваться более общей формой уравнения движения.

Обоснованное решение задачи дано в работах Р. Сквайрса и У. Галта [5].

Авторы рассматривали случаи шприцевания материала с переменной по

сечению канала вязкостью, вызванной ненъютоновскими свойствами

термопластичного материала. Если предположить, что течение материала в

канале является одномерным, то уравнение (2) сводится к виду:

∂P ∂Z

=

µ



∂ 2ϑ Z ∂ y2

+

∂ϑZ ∂y



∂µ ∂y

(16)

В этом уравнении µ является переменной величиной, которая

изменяется по глубине канала, т.е. вдоль оси y. Интегрируя уравнение (16)

получено для распределения скоростей следующее выражение:

ϑZ

=

UZ ln µ /

⋅ ln(1 −

y h

(1



µ

/ )) −

y⋅ h µS − µb



∂ ∂

P Z

+

(µ S



h µ b ) ln µ

/



∂P ∂Z

ln(1 −

y h

(1



µ

/ ))

(17)

где µ S - вязкость материала на дне канала шнека; µ b - вязкость материала у стенок корпуса.

Обозначив

µ/

=

µb µS

,

проинтегрируем

уравнение

(17) по площади

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

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

Q=

Fµ α

⋅α

⋅N−



P



µ

β
CP



∂P ∂e

(18)

где µ CP - средняя вязкость материала, находящегося в канале шнека;

Fµ α , Fµ P - коэффициент, учитывающий влияние изменения вязкости материала по глубине канала на вынужденное течение и противоток. При этом:

Коэффициенты Fµ α и

µ CP =

µb+ µS 2

(19)

Fµ P зависят только от µ / и определяются

следующими уравнениями:

Fµ α

=

2( µ

µ
/

/


1



1 ln µ

/

)

(20)

Fµ P

=

3(1 + 1−

µ µ

/ /

)

(1



Fµ α

)

(21)

В работе дана графическая зависимость поправочных коэффициентов

Fµ α и Fµ P от величины µ / . Для изотермического режима величина градиента

давления в канале постоянного сечения определяется выражением: (2)

Фактическая производительность шнекового пресса определяется

взаимным влиянием характеристики шнека и матрицы. В общем случае

величина расхода через матрицу может быть описана уравнением вида:

g = K ⋅ ∆ Pµ µM

(22)

где K - постоянный коэффициент, зависящий от геометрических размеров

матрицы;

∆ PM - перепад давления в матрице; µ M - вязкость жидкости в матрице. Для матрицы, имеющей цилиндрический круглый канал, значение

может быть вычислено по формуле:

K

=

π ⋅ RM4 8 ⋅ LM

(23)

где RM - радиус канала; LM - длина канала.

Совместное решение уравнений (22) и (23) позволяет определить

давление в матрице и производительность рассчитываемого пресса, которая

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

Правильный подбор этих характеристик позволит повысить

производительность процесса при снижении количества сырья, что является

одним из важнейших факторов общего снижения ресурсоемкости производства

соломки.

Таким образом в настоящее время отсутствует общая математическая

теория течения материалов в шнековом прессе. Решение этой задачи во многом

связано с реологическими свойствами (упруго-пластично-вязкими) и

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

реологических свойств теста соломки и экспериментальные исследования

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

Заключение

Для вывода дифференциальных уравнений течения и уравнения

объемного расхода для канала шнекового пресса использована

гидродинамическая теория неньютоновской жидкости, которая позволяет

установить количественные зависимости между геометрическими

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

режимом обработки.

Список литературы 1. Андреев А.Н., Мачихин. С.А. Механизация производства соломки. М.: ЦНИИТЭИ хлебопродуктов. - 1992. - 48 с. 2. Андреев A.Н. Маслова Г.В и др. Моделирование узлов формования соломки из рыбного теста. – Л. Сб. науч. трудов ЛИСТ, вып. 64. 1977. С. 23-28. 3.Арет В.А., Мачихин Ю.А. Формование конфетных масс выдавливанием. - М.: МТИПП. - 1969. - 35 с. 4.Берман Г.К., Ворожцов Л.А., Мачихин Ю.А. Течение вязкопластических масс в шнеке. - Известия ВУЗов. Пищевая технология. 1970. - №3. - С. 160-161. 5. Шенкель Г. - Шнековые прессы для пластмасс. – М.: - Госхимиздат. 1966. – 465 с.

For non-Newtonian fluid in a screw press
Andreev A.N.
andreevanatoly@yandex.ru
Saint-Petersburg State University of information technologies, mechanics and optics Institute of refrigeration and biotechnologies
The output of the differential equations for the flow and volume flow equations for the feed screw press used in hydrodynamic theory of non-Newtonian fluid which allows you to set a quantitative relationship between the geometrical characteristics of the working space, material properties and processing mode. Keywords: non-Newtonian liquid dough straws, screw press.