Текст произведения
(PDF):
Читать
Скачать
Введение Объективный контроль вязкоупругих свойств гелеобразных веществ является необходимым фактором в производстве многих пищевых продуктов. Одним из примеров служит контроль гелеобразования в молоке при производстве сыров. Существуют различные способы такого контроля [1], но большинство из них, прямо или косвенно, связаны с наблюдением изменения реологических параметров среды в процессе гелеобразования. В данной статье мы остановимся лишь на нескольких методах, основанных на непосредственном механическом взаимодействии измерительного устройства с коагулирующей средой. Одним из первых достаточно простых приборов, позволяющих непосредственно наблюдать и записывать изменение реологических характеристик при гелеобразовании в молоке стал Formagraph (Foss Electric, Hillerød, Denmark). Этот прибор широко использовался в конце XX века [2, 3] и благодаря своей простоте и надежности конструкции продолжает использоваться даже сегодня, когда появились более совершенные методы [4, 5]. Formagraph представляет собой небольшой маятник, обычно в форме петли, опущенный в кювету с коагулирующей жидкостью. Кювета совершает возвратно-поступательное движение, увлекая за собой маятник. Однако пока вязкость жидкости невелика, амплитуда колебаний маятника пренебрежимо мала. По мере образования структуры в жидкости воздействие на маятник усиливается и амплитуда его колебаний возрастает. Амплитуда колебаний записывается в течение всего процесса, позволяя сохранять и анализировать данные (в первых вариантах прибора запись производилась на фотопленку). Типичная зависимость размаха колебаний маятника от времени при коагуляции молока показана на рис. 1. Рис. 1. Диаграмма сычужного свертывания молока, полученная с помощью Formagraph [2] Как видно из рис. 1, анализ формы графика позволяет не только определить время коагуляции (r), но и ввести объективные параметры, характеризующие кинетику образования и прочность сгустка (k20, a30). Одним из самых распространенных реологических методов наблюдения структурообразования в пищевых продуктов на сегодняшний день является метод динамической реологии [6-8]. Он заключается в том, что образец подвергается сдвиговой синусоидальной деформации e небольшой амплитуды . Возникающее при этом в образце напряжение сдвига q совпадает с деформацией по частоте, но отличается по фазе колебаний . Разность фаз колебаний d определяется вязкоупругими свойствами образца. Так, для абсолютно упругого материала деформация и напряжение колеблются в одной фазе. Для абсолютно вязкого материала сдвиг фаз между напряжением и деформацией составляет p/2. Для вязкоупругих материалов сдвиг фаз принимает промежуточные значения. Так как сдвиг фаз между деформацией и напряжением обусловлен потерями на преодоление сил сопротивления, он характеризует потери механической энергии. Величина тангенса угла потерь определяется отношением двух динамических модулей: модуля потерь и модуля упругости. Физический смысл этих величин заключается в том, что модуль упругости определяет количество сохраненной энергии за каждый период колебаний. Модуль потерь определяет энергию, рассеянную за период в виде тепла. Изучение зависимости динамических модулей от времени в течение технологических процессов дает важную информацию о структурных изменениях сгустка и физико-химических факторах, лежащих в основе этих процессов. Тангенс угла потерь является важной физической величиной, характеризующей скорость релаксации связей в геле. Целью данного исследования является разработка способа исследования гелеобразования в пищевых продуктах, сочетающего в себе простоту устройства Formagraph и информативность динамической реологии. Прибор, в основе которого лежит этот способ, мы условно назвали «динамическим формографом», чтобы подчеркнуть связь с двумя рассмотренными выше методами. Объект и методы исследования Объектом исследования являлся лабораторный прототип динамического формографа, устройство которого описано ниже. В качестве «эталонов» для градуировки, настройки и проверки работоспособности прибора использовались: упругие элементы в виде стальных пружинок с известной жесткостью; вязкие элементы в виде растительного масла и глицерина; вязкоупругие элементы в виде водных растворов желатина. Лабораторный стенд Неподвижный зонд 3 (в виде шарика или пересекающихся пластинок) связан с датчиком механического напряжения 4. Датчик напряжения регистрирует силу F(t), действующую на зонд со стороны среды, перемещающейся относительно зонда по закону X(t). Сосуд 2 может совершать вертикальные синусоидальные колебания заданной частоты (0,1¸10,0 Гц) и амплитуды (0,5¸6,0 мм). Его положение определяется датчиком 5 и задается шаговым двигателем с редуктором и кривошипом 6. Управление шаговым двигателем осуществляется блоком 7. Сигналы с датчиков перемещения корзины и механического напряжения, испытываемого зондом, преобразуются и усиливаются блоком 8. Блоки усилителей-преобразователей 8 и управления шаговым двигателем 7 подключены к одному источнику питания 9. Низкоуровневые функции аналого-цифрового преобразования сигналов и программной реализации USB для связи с персональным компьютером 10 осуществляются аппаратным комплексом 10 (NATIONAL INSTRUMENTS ELVIS II+). Рис. 2. Блок-схема лабораторной установки (описание в тексте статьи) Датчик положения корзины представляет собой дифференциальный трансформатор (LVDT - linear variable differential transformer). Градуировочная кривая, снятая с датчика через усилитель при ходе корзины ±6 мм, показана на рис. 3. Как видно из рис. 3, характеристику датчика можно считать практически линейной при амплитуде смещения корзины до 3 мм. Рис. 3. Амплитудная характеристика датчика смещения корзины Датчик напряжения собран на основе тензометрических датчиков для настольных лабораторных весов и имеет чувствительность примерно 0,1 мН. На рис. 4 представлен градуировочный график зависимости силы, действующей на зонд со стороны прикрепленной к корзине тонкой стальной пружинки жесткостью 5 Н/м. Рис. 4. Зависимость силы, действующей на зонд, для упругой вынуждающей силы Шаговый двигатель с редуктором взят от принтера EPSON LX, кривошипная пара изготовлена самостоятельно с возможностью регулировки хода корзины. Блок управления шаговым двигателем представлен драйвером для униполярных шаговых двигателей на микросхеме КР561ТМ2 с регулируемым генератором П-образных импульсов на КР1006ВИ1. Усилитель-преобразователь служит для согласования уровней сигналов с датчиков с аппаратным комплексом NATIONAL INSTRUMENTS ELVIS II+. На персональный компьютер установлены драйверы аппаратного комплекса NATIONAL INSTRUMENTS ELVIS II+ и программа LabVIEW 8.5, с помощью которой непосредственно и производится обработка и запись оцифрованного аналогового сигнала. Математическая модель установки Первый опыт построения установки для измерения динамических характеристик сгустков на основе линейного возвратно-поступательного движения зонда в виде шарика был осуществлен несколько ранее [9]. Однако предыдущая версия установки требовала либо компенсации достаточно большой выталкивающей силы, либо учета сил инерции при движении шарика. Новая установка лишена этих недостатков, так как зонд в ней неподвижен. Рассмотрим упрощенную модель процесса, учитывая, что датчик напряжения, к которому прикреплен шарик не идеально твердый, а является упругим элементом с жесткостью k0. Пусть на подвешенный на пружине жесткостью k0 шарик массой m, уравновешенный в жидкости ( ), действует вынуждающая сила со стороны возвратно-поступательно движущейся относительно шарика вязкоупругой среды. Если x - координата шарика; X - координата, характеризующая смещение среды; k - константа, характеризующая упругость среды; a - константа, характеризующая вязкость среды, то уравнение движения шарика, при условии неразрушения шариком упругих связей среды (случай малых смещений), уравнение движения шарика может быть представлено следующим образом: , где - смещение шарика относительно среды, а - относительная скорость шарика . Пусть возвратно-поступательное движение среды представляет собой гармонические колебания с амплитудой A и циклической частотой W: . (1) Тогда уравнение движения шарика примет вид: , (2) где ; ; . Уравнение (2) можно решить, выбрав подходящие начальные условия, например, и . Однако очевидно, что общим решением этого уравнения является сумма затухающих собственных колебаний и вынужденных колебаний. Собственные колебания, являющиеся решением однородного уравнения типа (2), затухают за время порядка . Например, в случае ламинарного движения шарика в жидкости с вязкостью h, , где R - радиус шарика. Для шарика радиусом 2 см в глицерине при 20 °С . То есть при частоте колебаний, равной 0,1 Гц, затухание собственных колебаний происходит за время порядка одного периода колебаний. Таким образом, имеет смысл записать решение уравнения (2) для больших времен (t >> t) в виде вынужденных колебаний. Вынужденные колебания представляют собой частное решение уравнения (2), которое может быть представлено в форме суммы синуса и косинуса, изменяющихся с частотой вынуждающей силы: . (3) Подставляя это решение в дифференциальное уравнение (2), получим следующую систему алгебраических уравнений: Решением этой системы будет: (4) В соответствии с правилами тригонометрии решение (3) можно записать в виде синуса со сдвигом фаз: , (5) где - амплитуда смещения шарика, а - сдвиг фазы вынужденных колебаний шарика по отношению к колебаниям среды (1). Как видно из решения, в нашем случае сдвиг составляет: . (6) Очевидно, что в случае чисто упругой вынуждающей силы w ≠ 0, b = 0, а j = 0, то есть механическое напряжение на зонде совпадает по фазе с вынуждающей силой. Такой вывод полностью согласуется с данными, приведенными на рис. 4. В случае же чисто вязкой вынуждающей силы w = 0, b ≠ 0, а . Что не соответствует, однозначно, «классическому» значению j = p/2. Результаты и выводы Последнюю формулу для сдвига фаз между смещением среды и механическим напряжением, возникающем на зонде, можно представить в следующей форме: , где Fmax - максимальная сила механического напряжения, возникающего на зонде; mamax - максимальная сила инерции, действующая на зонд; Fc - часть вязкой силы, действующей на зонд из-за его собственного движения. Совершенно очевидно, что в случае стремления жесткости измерительной пружины к бесконечности, смещение зонда стремится к нулю, а сдвиг фаз к «классическому» значению j = p/2. Оценим параметры для нашей установки. Для глицерина, например, при 20 °С, при радиусе шарика зонда R = 2 см, его массе 30 г, частоте 0,3 Гц и амплитуде сдвига корзины Xmax = 5 мм, W = 1,9 с-1 и Fmax =3 мН. Оценка жесткости тензометрического подвеса дала величину k0≥ 105 Н/м. Откуда: xmax ≤ 3×10-5 мм. Тогда » 2×10-5 мН; »3×10-6 мН. Таким образом, в нашей установке для чисто вязких жидкостей реализуется «классический» сдвиг фаз между смещением среды и механическим напряжением, равным p/2. Рис. 5. Зависимость силы, действующей на зонд, для вязкой вынуждающей силы В качестве подтверждения полученного выше вывода на рис. 5 представлены графики зависимостей смещения среды X и силы механического напряжения F от времени для растительного масла. Так как вязкость масла примерно в 10 раз меньше, чем у глицерина, Fmax составляет около 0,3 мН, что находится вблизи предела чувствительности установки. Рис. 6. Зависимость амплитуды силы, действующей на зонд, в растворе желатина от времени Несмотря на малую амплитуду сигнала, соответствующего механическому напряжению на зонде, и тот факт, что при амплитуде смещения корзины, равной 5 мм, датчик смещения работает в нелинейном режиме (рис. 3), хорошо видно, что максимальной силе механического напряжения примерно соответствует нулевое значение сдвига корзины (то есть ее максимальная скорость). Таким образом, сдвиг фаз, действительно, соответствует примерно четверти периода. На рис. 6 показана работа установки в режиме формографа при гелеобразовании в растворе желатина. Видно, что зависимость амплитуды механического напряжения, приложенного к зонду, качественно соответствует сигналу, представленному на рис. 1. Анализируя рисунок, можно заметить, например, что в районе 460 с вблизи зонда произошло локальное разрушение (расслоение) сгустка. В заключение отметим, что проведенные исследования подтверждают возможность использования разработанного прибора как для определения момента начала структурообразования в коллоидных растворах, так и для исследования динамических характеристик.вязкоупругих.сред.