Text (PDF):
Read
Download
Введение При производстве сладкого йогурта в Кемеровской области на ООО «Деревенский молочный завод» широко применяется резервуарный метод. В отличие от классической технологии стабилизатор находится в сухом виде (модифицированный крахмал) и вносится в молоко совместно с сахаром и сухим молоком. Эти ингредиенты в необходимых пропорциях смешиваются в отдельном бачке ручным способом. Качество такого процесса не удовлетворяет современным требованиям, в результате чего при внесении сухой йогуртной основы в молоко стабилизатор и сухое молоко слипаются в прочные и нерастворимые конгломераты. В этих условиях готовый продукт не соответствует заданной рецептуре и желаемому конечному качеству [3]. Для решения данной проблемы предлагается включить в технологическую линию по производству йогурта на стадию смешивания его сухой основы смесительный агрегат, включающий в себя дозаторы объемного типа и барабанный смеситель непрерывного действия (Патент РФ № 2508937) [1, 4, 6]. Последний обладает низкой энергозатратностью, эффективностью смешивания, сохранением структуры смешиваемых компонентов при их соотношении в диапазоне от 1:10 до 1:40, а также возможностью совмещенного продольного и поперечного смешивания сыпучих компонентов за счет наличия различных контуров рециркуляции смеси по всей длине барабана. Модернизированная технологическая схема производства йогурта представлена на рис. 1. Рис. 1. Модернизированная технологическая схема производства йогурта: 1 - емкость для сырого молока; 2 - насосы; 3 - балансировочный бачок; 4 - пластинчатая пастеризационно-охладительная установка; 5 - пульт управления; 6 - оборотный клапан; 7 - сепаратор-нормализатор; 8 - гомогенизатор; 9 - емкость для выдерживания молока; 10 - емкость для йогурта; 11 - барабанный смеситель; 12 - блок дозаторов; 13 - пульт управления; 14 - емкость для предварительного смешивания молока с сухой смесью Цель работы - определение рациональных конструктивных и технологических параметров работы нового барабанного смесителя на основе регрессионного анализа результатов теоретических и экспериментальных исследований при получении сухой йогуртной основы. В соответствии с поставленной целью решались следующие задачи: - исследовать усовершенствованную конструкцию смесителя барабанного типа с целью выявления рациональных конструктивных и технологических параметров работы, обеспечивающих производство качественной смеси; - провести регрессионный анализ полученных данных для выявления степени влияния варьируемых параметров на качество получаемой сухой йогуртной основы; - получить регрессионную модель, позволяющую прогнозировать качество конечного продукта. Объекты и методы исследований В качестве объекта исследований была взята новая конструкция барабанного смесителя с установленными внутри Г-образными лопастями для создания внутренней рециркуляции получаемой смеси при соотношении компонентов 1:40. Предметом исследования являлось выявление рациональных конструктивных и технологических параметров барабанного смесителя, влияющих на качество смешения при получении сухой йогуртной основы. Для реализации поставленных задач использованы разработанные ранее на кафедре «Технологическое проектирование пищевых производств» методики по определению концентрации ключевого компонента в смеси с применением статистических методов обработки экспериментальных данных [1]. Результаты и их обсуждение Для достижения поставленной цели были проведены исследования по выявлению степени влияния частоты вращения барабана n, коэффициента заполнения барабана К и расположения Г-образных лопастей на качество многокомпонентных смесей (сухая йогуртная основа), оцениваемое коэффициентом вариации Vc. В ходе экспериментальных исследований изменялась частота вращения барабана n в диапазоне 10-40 (об/мин), коэффициент заполнения барабана К в диапазоне 10-30 (%), шахматное и спиралевидное расположение Г-образных лопастей. Базовые (нулевые) точки и шаги варьирования приведены в табл. 1. Таблица 1 Базовые точки и шаги варьирования Фактор и его обозначение Верхний уровень Нижний уровень Центр плана Интервал варьирования Частота вращения барабана Х1 40 10 25 15 Коэффициент заполнения барабана Х2 30 10 20 10 Рецептура сухой йогуртной основы и физико-механические характеристики смешиваемых компонентов, при которых осуществлялось смешивание, представлены в табл. 2. Таблица 2 Рецептура сухой йогуртной основы Наименование сырья Сахар Сухое молоко Стабилизатор Количество сырья в 100 кг молока 5,78 кг 1,66 кг 0,25 кг Плотность, кг/м3 900 450 560 Влажность, % 0,05...0,1 3...4 18...20 Диаметр частиц, мкм 940...960 80...100 150...250 Варьируемые параметры и полученное качество смешивания, оцениваемое при помощи коэффициента неоднородности Vc, сведены в табл. 3. Таблица 3 Варьируемые параметры и значения коэффициента неоднородности Vc В кодированном виде В натуральном виде Коэффициент неоднородности Vc, % Х1 Х2 n, об/мин К, % Шахматное расположение лопас-тей Спиралевидное расположение лопастей - - 10 10 9,61 9,11 + - 40 10 10,25 8,27 0 - 25 10 11 7,38 - + 10 30 7,25 7,40 + + 40 30 8,29 6,5 0 + 25 30 6,89 5,4 - 0 10 20 6,08 6,7 + 0 40 20 9,32 6,12 0 0 25 20 7,54 5,05 Из табл. 3 видно, что наилучшее качество получается при спиралевидном расположении Г-образных лопастей. Это объясняется тем, что при данном расположении лопастей осуществляется большая рециркуляция смеси в рабочей зоне аппарата, приводящая к значительному усреднению материальных потоков и увеличению времени пребывания частиц в аппарате. Также видно, что лучшие результаты получаются при К, равном 20 %, это объясняется тем, что при малой степени заполнения не происходит достаточного послойного перемешивания, характерного для барабанных смесителей, а при большой степени заполнения начинает создавать помехи центральный вал, на котором закреплены лопасти. Для лучшего восприятия полученных данных представили их в виде поверхностей отклика, полученных в программе Statistica 8. Для шахматного расположения Г-образных лопастей поверхность отклика представлена на рис. 2. Рис. 2. Зависимость коэффициента неоднородности от частоты вращения и коэффициента заполнения барабана Рис. 3. Зависимость коэффициента неоднородности от частоты вращения и коэффициента заполнения барабана Из рис. 2 видно, что при работе барабанного смесителя с расположением Г-образных лопастей в шахматном порядке достигается наилучшее качество сухой йогуртной основы при частоте вращения барабана 10 об/мин и коэффициенте заполнения аппарата 20 %. Рис. 3 показывает, что при смешивании сухой йогуртной основы на барабанном смесителе со спиралевидным расположением Г-образных лопастей следует выбирать частоту вращения барабана 25 об/мин и коэффициент заполнения аппарата 20 %. Дальнейший регрессионный анализ [2, 5] был направлен на определение степени влияния варьируемых параметров на качество сухой йогуртной основы, смешанной на барабанном смесителе с поочередно установленными Г-образными лопастями в шахматном или спиралевидном порядке. Для этого применялась статистическая программа Statis-tica 8 с использованием модуля «Общие регрессионные модели». На первом этапе моделирования провели анализ работы смесителя с шахматным расположением Г-образных лопастей. В табл. 4 приведены оценки данной модели. Таблица 4 Оценка модели Множественный коэффициент корреляции, R 0,846 Множественный коэффициент детерминации, R2 0,715 Скорректированный коэффициент корреляции, R 0,621 F-критерий 7,552 p-уровень 0,0229 Коэффициент корреляции (R), равный 0,846, близок к единице, что говорит о сильной зависимости выходной переменной от входных переменных. Коэффициент детерминации (R2) данной модели, равный 0,715, также высок. Из этого следует, что доля дисперсии зависимой переменной, объясняемая рассматриваемой моделью зависимости, равна 71,5 %. F-критерий Фишера имеет достаточное значение (F = 7,552), чтобы утверждать, что модель является адекватной и может быть применима для прогнозирования коэффициента Vc. Рассматриваемая модель является статистически значимой, так как p-уровень составляет 2,29 %. Это показывает, что модель с вероятностью 0,0229 будет являться лишь случайным совпадением для данной выборки. В табл. 5 приведены коэффициенты регрессии модели. Таблица 5 Коэффициенты модели Vc, % Param Vc, % t Vc, % p Vc, % Beta Св. член 9,913 8,680 0,00013 n, об/мин 0,054 1,959 0,0978 0,426 K, % -0,141 -3,357 0,0153 -0,731 Статистическая значимость (p-уровень) у коэффициентов низкая, в пределах 0,01-1,5 %, лишь у одного коэффициента чуть выше - 9,78 %. Это показывает, что каждый найденный коэффициент с вероятностью, равной соответствующему ему p-уровню, будет говорить, что найденная зависимость является лишь случайной особенностью данной выборки. Аналогичные результаты отображает t-критерий Стьюдента. Он достаточно высок, что также подтверждает статистическую значимость коэффициентов. В соответствии с этим оценены коэффициенты Beta, показывающие меру чувствительности одной переменной к другой. Из табл. 5 видно, что наибольшее влияние на качество сухой йогуртной основы оказывает коэффициент заполнения барабана K, причем это влияние обратно пропорционально, т.е. с возрастанием К коэффициент неоднородности Vc снижается. Влияние второй переменной - частоты вращения в исследуемом диапазоне на Vc менее заметно по отношению к коэффициенту заполнения, однако это не дает оснований для ее иcключения из регрессионного анализа. Сложив коэффициенты (параметры модели Param) из табл. 5, получаем итоговую регрессионную модель (уравнение регрессии), позволяющую прогнозировать значения зависимой переменной Vc: . (1) По данному уравнению регрессии построили поверхность отклика, представленную на рис. 4. Рис. 4. Зависимость коэффициента неоднородности от частоты вращения и коэффициента заполнения барабана С помощью полученной модели сравнили наблюдаемые значения (полученные в ходе эксперимента, табл. 3) зависимой переменной с предсказанными (полученные с помощью математической модели). Адекватность предсказанных значений коэффициента неоднородности оценили с помощью относительной погрешности, рассчитанной по формуле . (2) Полученные результаты сведены в табл. 6. Таблица 6 Сравнение наблюдаемых и предсказанных значений № Vc, % Наблюдаемые Vc, % Модельные Относительная погрешность, % 1 9,61 9,05 5,77 2 10,25 10,69 4,34 3 11,0 9,87 10,22 4 7,25 6,24 13,86 5 8,29 7,88 4,88 6 6,89 7,06 2,53 7 6,08 7,65 25,82 8 9,32 9,29 0,32 9 7,54 8,47 12,33 Из табл. 6 видно, что среднее значение относительной погрешности составляет 8,9 %, а как известно, инженерные расчеты допускают не более 10 %, следовательно, модель можно использовать для предсказания качества сухой йогуртной основы, полученной на барабанном смесителе с шахматным расположением Г-образных лопастей. На втором этапе моделирования провели анализ работы смесителя со спиралевидным расположением Г-образных лопастей. В табл. 7 приведены оценки модели. Таблица 7 Оценка модели Множественный коэффициент корреляции, R 0,998 Множественный коэффициент детерминации, R2 0,996 Скорректированный коэффициент корреляции, R 0,992 F-критерий 271,68 p-уровень 0,00004 Коэффициент корреляции (R), равный 0,998, практически приблизился к единице, что говорит о сильной зависимости выходной переменной от входных переменных. Помимо этого, p-уровень составляет 0,004 %, что подтверждает статистическую значимость модели. В табл. 8 приведены коэффициенты регрессии модели. Таблица 8 Коэффициенты модели Vc, % Param Vc, % t Vc, % p Vc, % Beta Св. член 16,937 50,272 0,000001 n, об/мин -0,338 -18,920 0,000046 -3,357 n2, об/мин 0,00625 17,738 0,000059 3,147 K, % -0,646 -20,146 0,000036 -4,270 K2, % 0,0139 17,486 0,000063 3,706 Статистическая значимость (p-уровень) у всех коэффициентов находится в пределах 0,001-0,006, что говорит о высокой статистической значимости рассматриваемых коэффициентов. В соответствии с этим оценены коэффициенты Beta, показывающие, что все параметры практически в равной степени влияют на Vc. Сложив коэффициенты (параметры модели Param) из табл. 8, получаем итоговую регрессионную модель (уравнение регрессии), позволяющую прогнозировать значения зависимой переменной Vc: . (3) По данному уравнению получили поверхность отклика, представленную на рис. 5. Рис. 5. График зависимости коэффициента неоднородности от частоты вращения и коэффициента заполнения барабана Аналогичным образом с помощью полученной модели сравнили экспериментальные значения (табл. 3) коэффициента неоднородности с полученными с помощью математической модели. Адекват-ность предсказанных значений Vc оценили по формуле (2). Результат сравнения приведен в табл. 9. Таблица 9 Сравнение наблюдаемых и предсказанных значений № Vc, % Наблюдаемые Vc, % Модельные Относительная погрешность, % 1 9,11 9,108 0,01 2 8,27 8,34 0,80 3 7,38 7,32 0,87 4 7,40 7,29 1,50 5 6,50 6,52 0,24 6 5,40 5,49 1,77 7 6,70 6,81 1,67 8 6,12 6,04 1,33 9 5,05 5,02 0,62 Из табл. 9 видно, что среднее значение относительной погрешности составляет всего 0,98 %. Следовательно, модель может считаться адекватной и применяться с высокой степенью точности для предсказания качества сухой йогуртной основы, полученной на барабанном смесителе со спиралевидным расположением Г-образных лопастей. Ввиду приведенных выше рассуждений полученные модели адекватны и могут быть применимы для дальнейшего предсказания коэффициента неоднородности при получении сухой йогуртной основы в зависимости от исследуемых технологических и конструктивных параметров работы барабанного смесителя непрерывного действия. Выводы 1. Получены рациональные параметры работы смесителя непрерывного действия барабанного типа при получении сухой йогуртной основы: n = 25 об/мин, K = 20 %, спиралевидное расположение Г-образных лопастей, при которых коэффициент неоднородности имеет свое наименьшее значение 5,5 %. 2. При работе барабанного смесителя с шахматным расположением Г-образных лопастей наилучшие показатели качества (Vc = 6,08 %) достигаются при частоте вращения барабана 10 об/мин и коэффициенте заполнения 20 %. 3. На основе инструмента программы Statistica получили математические регрессионные модели, позволяющие с высокой степенью точности предсказывать качество получаемой сухой йогуртной основы. Значения относительной погрешности между экспериментальными значениями зависимой переменной Vc и предсказанными составили 8,9 и 0,98 % для шахматного и спиралевидного расположения Г-образных лопастей соответственно.