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

Т.М. Зубкова, А.Ш. Насыров

УЧЕТ ХАРАКТЕРА ДВИЖЕНИЯ МАТЕРИАЛА В КАНАЛЕ ШНЕКА ПРИ МАТЕМАТИЧЕСКОМ МОДЕЛИРОВАНИИ ЭКСТРУДИРОВАНИЯ

РАСТИТЕЛЬНОГО СЫРЬЯ

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

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

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

Шнековый канал каждой секции прессующего механизма, как показано на рисунке 1, можно рассматривать как часть пространства между бесконечными параллельными плоскостями [2]. Предположим, что верхняя плоскость 1 движется относительно нижней плоскости 2 со скоростью ус .

Если пренебречь силами инерции, уравнение движения для модели шнекового канала, изображенного на рисунке 1, имеет вид [2]

da , ч

= ~dx(у-Уо

где Ту - напряжение сдвига в прессуемом материале;

da

--градиент нормальных напряжений в прес-

dx

суемом материале (сжимающим напряжениям приписаны положительные значения); Уо - координата плоскости, на которой касательные напряжения тХу = 0.

Анализ решений уравнения (1), проведенный С.А. Бостанджияном и А.М. Столиным [2], позволяет выделить два возможных случая распределения касательных напряжений в канале шнека:

- касательные напряжения не изменяют направления, то есть уо < 0;

- касательные напряжения меняют направление на противоположное, то есть 0 < уо < Ъш .

Рассмотрим расход материала через секцию прессующего механизма.

Будем считать, что зависимость напряжения сдвига тХу от скорости сдвига jx (градиента ско-

dvx ,

рости —у) удовлетворительно описывается уравнением Оставальда - -де Виля

t = mrx = m

/ JL

dvx

dy

Для случая уо < 0 с учетом направления касательного напряжения т> 0 уравнение (1) преобразуется в уравнение

dv

x 2

dy

= -аш

где

\Ш

da'

dx

= m

Интегрируя уравнение (3), приняв на нижней пластине условие прилипания материала по поверхности [3], получим

vx2 =-

m +1

/ \m+1 / \m+1

1].

Объемный расход прессуемого материала на единичной ширине пространства между пластинами для этого случая имеет вид

б = =

(-Уо Г+1 К -

аш =

n

/, \ш+2 / \ш+2

йу

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

йо

йх

о21 -1 -о21

где

Аш1

= Аш1

(-Уо» )щ +1 К»

йу

х1

йу

= аш

секция 1111) секция

1 2

Рисунок 1. Схема модели шнекового канала: 1 - плоскость, замещающая шнековый цилиндр; 2 - плоскость, замещающая дно шнекового канала.

х 2

йу

ш

= -аш

где о21 -1 и о21 - нормальные напряжения соответственно в начале и в конце i -й секции прессующего механизма;

хш» - протяженность шнекового канала i -й секции.

Производительность i -й секции шнекового прессующего механизма численно равна модулю расхода прессуемого материала в канале шнека. Тогда выражение производительности для i -й секции в случае Уо < о можно переписать в виде

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

Принимая граничные условия

Ух1 = о пРи У=о, Ух2 = ус пРи У=К, (13) получим из уравнений (11) и (12) интегрированием уравнения распределения скоростей

Ух1 = ^ [Уош+1 - (Уо - У)ш+1 ] при у < Уо, (14) ш + 1

Ух2 = Ус + №ш - Уо)ш+1 - (У- Уш)ш+1] ш + 1

при У > Уо.

Из формул (14) и (15) объемный расход О прессуемого материала на единичной ширине пространства между пластинами для случая о < уо < Иш можно определить интегрированием элементарных потоков. Получим уравнение

Уо К

О = \Ух1йУ + |Ух2^У = ус ('ш - Уо) +

о Уо

\ш+1п

[Уош+2 + (К - Уо)ш+2]

(¡хш1) (ш + 1)' Для случая о < уо < Ъш уравнение (1) имеет два решения

+ 1 .

ш + 2 у '

Выражение производительности для i -й секции в случае Уо > о можно переписать в виде

О = Сшг

где Сш»=-

ш + 2 . п ш + 2

Уо! +

Вш1 = ус»

Обеспечивая условие непрерывности скоростей ух1 = Ух2 при у = уо, из уравнений (14) и (15), можно получить соотношение

Уош+1

ш +1 ш + 1

которое имеет действительные корни при о < Уо < К.

Анализ уравнения (20) показывает, что это уравнение имеет действительные корни при скорости верхней пластины ус , заключенной в пределах

аш г.ш+1 < У < аш г,ш+1

ш + 1

уе — , 1 ш ш +1

Нас в этом интервале интересует абсолютная

величина ус .

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

1

к >

1 + п

К

1+п

т

do

dx

где п - индекс течения экструдируемого материала;

¡1 - коэффициент консистенции экструдируе-мого материала.

Касательные напряжения изменяют направление на противоположное, если выполняется условие

1

Ы £ т

+ п

¡1 dx

°21 = °2(г+1)-1

Аш1 {021 -1- аъ )щ = Аш

Если в соседних секциях выполняются условия (22) и (23), уравнение неразрывности записывается следующим образом

Сш1 (02г-1 -02г ) + Ашг =

= А,.

ут+1

Скорость ус является внешней величиной процесса, которая не изменяется по длине шнека. Обычно внешними величинами являются и другие параметры, входящие в выражения (22) и (23), кроме нормальных напряжений. Анализ показывает, что с увеличением температуры материала его свойства приближаются к свойствам ньютоновской жидкости. Это выглядит как уменьшение коэффициента консистенции и увеличение (приближение к единице) индекса течения [3], что приводит к уменьшению правой части в выражениях (22) и (23). Как было показано раньше [4], реологические свойства холодного корма обеспечивают выполнение условия (22) в области шнека, примыкающей к загрузочному отверстию, однако в области, примыкающей к матрице, возможен переход к условию (23). На рисунке 2 показано поведение границы между условиями (22) и (23) при различных значениях градиента напряжений и глубины шнекового канала. С ростом градиента напряжений граничная скорость резко увеличивается, причем с увеличением глубины канала данное значение граничной скорости достигается при большем градиенте напряжений. Граничная скорость сильно зависит от реологических свойств эк-струдируемого материала.

Если предположить постоянство плотности экструдируемого материала между секциями, уравнение неразрывности между соседними секциями при выполнении в секциях условия (22) и условия непрерывности нормальных напряжений на границе между соседними секциями

Рисунок 2. Диаграммы зависимостей, соответствующих

значению У0 = 0, при движении продукта без проскальзывания по нижней пластине при реологических параметрах прессуемого материала ¡1 = 0,0047 Па-с11; п = 0,2 и значениях глубины канала шнека: 1 -Кш1 = 0,012 м; 2 - Нш2 = 0,009 м; 3 - Кш3 = 0,006 м, а также ¡1 = 0,062 Па-с1; п = 0,2 и значениях глубины канала шнека: 4 - Кш1 = 0,012 м; Нш2 = 0,009 м; Кш3 = 0,006 м.

Наконец, если в соседних секциях выполняется условие (23), уравнение неразрывности имеет вид

Сшг (°2г-1 - °2г ) + Ашг = \Щ+1

= С,,

+ А.

.

имеет вид

Полученные зависимости (24), (25), (26) позволяют приближенно вычислять распределение напряжений по длине шнека, а на основании этого производить расчет комплекса параметров эффекта процесса экструдирования.

Ранее было показано [3], что при экструдиро-вании кормов реологические свойства прессуемого материала могут вызывать движение прессуемого материала с проскальзыванием по шнеку. Математически это выражается граничным условием

УХ = У0 Ф 0 при у = 0.

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

ВЕСТНИК ОГУ 1'2003 149

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

Получим интегрированием уравнений (11) и (12) уравнения распределения скоростей

Ух1 = уо + [

ш + 1

Ух2 = Ус + От[('ш - Уо)ш+1 - (У-Уо)ш+1], У > Уо. (29) ш +1

Можно выразить неизвестную скорость скольжения, приравняв скорости по формулам (28) и (29) при у = уо. Получим

уо = Ус +-

Ä - У0Г+1 - y0m+1].

уо ku

Q = J yxldy + J Ух2dy = УоК + <

o уо

(hu - yo)m+1 Уо

m +1

п„ ,, \m+2 м+2 (hu - У0) - У0

m + 2

Выражение производительности для ' -й секции в случае Уо > о можно переписать в виде

О = °ш1 (О2»-1 - О2» )) + Рш»'

где Gu„- =

1

mi

м+2 /1 „ \m+2

У0 - (hu - У0) m + 2

(hu - yo)m+1 Уо m +1

F  = -y h

1 ui vci'lui

hu I

Q = J yx2 dy = ychu +0

(hu - yo)m+2 - (-yo)m+2]

au hu (hu - Уо>

m +1

m+1

Тогда выражение производительности для i -й секции в случае Уо < о можно переписать в виде

О = Мшг

('ш» - Уо»Л +1

где Mui =

ш +1

Объемный расход О прессуемого материала при наличии проскальзывания по нижней пластине на единичной ширине пространства между пластинами для случая о < уо < Нш можно определить интегрированием элементарных потоков, используя формулы (28) и (29). Получим уравнение

/, + 2 / \П

(тхш ) (ш + 1)(ш» + 2) Анализируя зависимость (31), можно сделать заключение, аналогичное условию (21), на основании которого напряжение сдвига изменяет свое направление в канале шнека при выполнении условия для расхода материала

hu

/ ; ; \т hu da

m dx

m + 2

£ Q £ hu

/ ; 1 \m hu a G

m' dx

hu

.

Это условие можно использовать в случае проскальзывания по нижней пластине для выбора расчетных уравнений определения производительности данной секции. В канале шнека производительность не может превысить величины -Ус'ш. Тогда неравенство (38) примет вид

-h2

m

hu da

m' dx

hu m + 2

Пусть Уо < о , тогда слой прессуемого материала движется между пластинами со скоростью, описываемой только уравнением (29), причем при у = о ух ф о. Объемный расход О прессуемого материала на единичной ширине пространства между пластинами для случая проскальзывания по нижней пластине имеет вид

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

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

Если предположить постоянство плотности экструдируемого материала между секциями, уравнение неразрывности между соседними сек-

Рисунок 3. Диаграммы зависимостей, соответствующих значению Уо = 0, для движения продукта с проскальзыванием по нижней пластине при реологических параметрах прессуемого материала / = 0,0047 Па-с11; п = 0,2 и значениях глубины канала шнека: 1 -Нш1 = 0,012 м; 2 - Нш2 = 0,009 м; 3 - НшЪ = 0,006 м, а также // = 0,062 Па-с1; п = 0,2 и значениях глубины канала шнека: 4 - Нш1 = 0,012 м; 5 - Нш2 = 0,009 м; 6 - НшЪ = 0,006 м.

циями, если в секциях не выполняется условие (39), при непрерывных нормальных напряжениях на границе между соседними секциями

°21 =^2(г+1)-1

имеет вид Мш1 (-1 -^2г ) + =

! \Щ+1

= Мш

Если в I -й секции выполняется условие (38), а в (г + 1)-й не выполняется, уравнение неразрывности записывается следующим образом

= M„

Grni (s2i-1 - s2i)) + Frni = \Щ+1

(i+1) [s2(i+1)-1

Г

S2

+ F

.

Наконец, если в соседних секциях выполняется условие (39), уравнение неразрывности имеет вид

(2г-1 -^2г Л' + Рш1 = ^ш(г +1)'

/ \Щ+1

'

Полученные зависимости (40), (41), (42) позволяют приближенно вычислять распределение напряжений по длине шнека, а на основании этого производить расчет комплекса параметров эффекта процесса экструдирования.

Список использованной литературы:

1. Полищук В.Ю. Особенности шнекового прессующего механизма экструдера. Тракторы и сельскохозяйственные машины, 1993, № 5. - С. 19-21.

2. Бостанджиян С.А., Столин A.M. Течение неньютоновской жидкости между двумя параллельными плоскостями. Изв. АН СССР, Механика, № 1, 1965. - С. 185-188.

3. Карташов Л.П., Полищук В.Ю., Зубкова Т.М. Моделирование процесса экструдирования в одношнековых прессующих механизмах. Техника в сельском хозяйстве, 1998, № 6. - С. 12-14.

4. Карташов Л.П., Полищук В.Ю., Зубкова Т.М., Ханин В.П. Учет изменяющейся температуры в математической модели экструдера. Техника в сельском хозяйстве, 2000, № 1. - С. 12-14.

5. Карташов Л.П., Полищук В.Ю., Зубкова Т.М. Формирование математической модели движения материала в одношнековом прессующем механизме. Техника в сельском хозяйстве, 1996, № 5. - С. 7-9.