Математическая модель процесса коэкструдирования пищевых масс
Автор: Фатыхов Ю.А., Шуманов В.А., Зарудный В.А.
Журнал: Вестник Воронежского государственного университета инженерных технологий @vestnik-vsuet
Рубрика: Процессы и аппараты пищевых производств
Статья в выпуске: 3 (57), 2013 года.
Бесплатный доступ
Рассмотрен процесс коэкструдирования пищевых масс. Получены уравнения движения в кольцевой и цилиндрической трубе двух соприкасающихся вязкопластичных жидкостей. Определены постоянные коэффициенты уравнений, зависящие от параметров процесса и конструктивных размеров коэкструзионной насадки.
Коэкструзия, пищевые массы, вязкость, модель, уравнения движения
Короткий адрес: https://sciup.org/14040088
IDR: 14040088
Текст научной статьи Математическая модель процесса коэкструдирования пищевых масс
В работах [1, 2] рассмотрен механизм и физическая картина получения комбинированных фаршевых изделий методом коэкс-трузии. В данной работе будет представлена математическая модель данного процесса.
Процесс формирования готовой продукции происходит в камере коэкструзион-ной насадки машинной системы Vemag893.
Камера формирования (рисунок 1) представляет собой полость цилиндрической формы, куда подается фаршевая оболочка. Коаксиально в данную полость подается начинка через кольцевой цилиндрический канал с клапаном, регулирующим подачу. Внизу камеры формирования установлена диафрагма с ножом в виде лепестков.
Рисунок 1 – Нижняя часть коэкструзионной насадки
1 – насадка; 2 – начинка; 3 – клапан; 4 – оболочка; 5 – камера формирования; 6 - диафрагма
Процесс формирования готового изделия связан с началом раскрытия ножей (лепестков) диафрагмы. По мере раскрытия диафрагмы поднимается клапан, связанный
Зарудный В.А., 2013
с пневматическим исполнительным механизмом, и по цилиндрическому кольцевому каналу внутрь оболочки выдавливается начинка. Формирование готового продукта связано с полным раскрытием и закрытием диафрагмы и клапана.
На рисунке 2 схематично представлено движение начинки и оболочки в камере форми
рования. С индексом 1 связаны величины, характеризующие начинку, с индексом 2 - оболочку.

Рисунок 2 - Схема движения начинки и оболочки в камере формирования
За основу математической модели возьмем уравнение Навье-Стокса в цилиндрических координатах для движения очень вязких жидкостей [3] в трубе:
равную нулю.
Окончательно получаем обыкновенное дифференциальное уравнение:
Эр
Эt
i ЭР + ^ / Э2 v + 1 Эг \ р Эг Эг 2 г Эг
v ^ T U + ”Т + а = 0 (5)
аг2 г аг уравнение для удобства интегрирования, представим в виде:
где v(r, t)- скорость, р - плотность, — — градиент давления, v - кинематическая вязкость.
Стационарная задача (^ = 0) коаксиального движения в трубе двухслойной жидкости при открытых клапане и диафрагме ранее была решена [4]. При решении нестационарного уравнения (1) используем зависимость градиента давления, меняющегося со временем по закону [5]:
- 1Эр = а |sin(^t)| , (2)
р Эг а скорость представим как:
v(r,t) = |sin(wt)| •U(r), (3)
Подставляя (2) и (3) в (1) и выделяя временную часть, получаем:
где U =
— (r—) + — = 0, аг аг v
( U
i
(r), 0
(U2(r), R
i
а граничные условия:
f U i (R i ) = U 2 (R i ),
U 2 (R 2 ) = 0,
R i V
dU t | аг Rt
au
= ^^ ,
a2u vau
a+v—7+^— dr2 г dr
MU
= |ct^(wt)|.
Уравнение (4) может иметь решение, если правая его часть равна постоянной величине. В нашем случае примем постоянную
где Ц 1 , ц2 - динамические вязкости.
Решение уравнения (6) с условиями (7) будет иметь следующий вид:
— т2 + СА ln r + С2,
4vt i 2'
— r2 + C3 In r + C4.
4V2 3 4
где Ci = 0 - из условия конечности скорости на оси трубы;
С3 определяется, исход из третьего условия (7), С4 - из второго, С2 - из первого.
Г
С 2
1 dp i
4^i dz
Rl +
1 ЭР 2
4^ 2 dz
Rl +
R2 / dp i
2 к dz
—
dpr )'" R 1
1 9P 2 R l
—
R l /ЭУ1
2 V dz
—
^ Ь R2-
V
C3
R
. 2
l
/ d p l к dz
—
®-
С 4 = — ^t^
+
Rj / dp i
2 к dz
—
at)'" 4
dp, dp2 .
В (9) "d^1 и "dz " - амплитуды градиентов давлений.
Окончательный вид нахождения скорости истечения жидкостей определяется подстановкой (8) и (9) в (3).
Объем жидкости (фарша), протекающий в единицу времени через поперечное сечение трубы, определим по формуле:
Q = 2 /J^1 2nw1dr + jRR2 2Hrv2dr). (10)
Коэффициент 2 получаем путем усреднения за время одного цикла функции |sin(^t)|:
т
-J02|sin[^t)|dt = -. (11)
Среднюю скорость течения можно опре- делить по формуле:
Q
-R2"
Коэффициенты С 2 , С 3 , С 4 , обозначенные в уравнениях (9), находят, исходя из известных параметров процесса и конструктивных размеров коэкструзионной насадки.
В заключение, следует сделать некоторые выводы:
-
- получены уравнения, описывающие движение в цилиндрической трубе двух соприкасающихся вязкопластичных жидкостей;
-
- определены выражения для постоянных коэффициентов уравнений движения, зависящих от параметров процесса и конструктивных размеров коэкструзионной насадки.