МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОЦЕССА МЕМБРАННОГО КОНЦЕНТРИРОВАНИЯ КРАХМАЛЬНОГО МОЛОКА НА ОСНОВЕ ПЕРЕДАТОЧНЫХ ФУНКЦИЙ
1
Хачатрян Л.Р.
1
1 ФГБОУ ВО «Кемеровский технологический институт пищевой промышленности (университет)», Кемерово
Настоящая статья посвящена разработке математической модели на основе передаточных функций, описывающей динамику процесса мембранного концентрирования крахмального молока в ультрафильтрационном аппарате. Описаны конструкция и принцип работы ультрафильтрационного аппарата с гидродинамическими элементами. Выявлены основные входные воздействия и выходные параметры объекта моделирования. Составлена структурная схема процесса мембранного концентрирования. Выбран вид входных воздействий. Установлен вид передаточных функций по каналам воздействия исходной концентрации, рабочего давления и температуры среды на удельную производительность. В процессе мембранного концентрирования крахмального молока экспериментально определены диапазоны входных воздействий. Посредством обработки экспериментальных зависимостей получены параметры передаточных функций. Математическая модель процесса мембранного концентрирования реализована программными средствами. Приведены результаты оценки адекватности модели экспериментальным данным.
математическое моделирование
передаточная функция
мембранное концентрирование
ультрафильтрация
крахмальное молоко
1. Антонов, А.В. Системный анализ. – М.: Высшая школа, 2004. – 454 с.
2. Волкова, В.Н. Основы теории систем и системного анализа / В.Н. Волкова, А.А. Денисов. – СПб.: СПбГПУ, 1999. – 512 с.
3. Ерофеев, А.А. Теория автоматического управления. – СПб.: Политехника, 2003. – 302 с.
4. Лурье, Б.Я. Классические методы автоматического управления / Б.Я. Лурье, П.Дж. Энрайт. – СПб.: БХВ-Петербург, 2004. – 640 с.
5. Пат. 152744 Российская Федерация. МПК7 B01D 63/00 (2006.01). Мембранный аппарат / Хачатрян Л.Р., Котляров Р.В.; заявитель и патентообладатель Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Кемеровский технологический институт пищевой промышленности». – № 2014105407/05; заявл. 13.02.14; опубл. 20.06.15, Бюл. № 17(Пч.). – 1 с.
6. Хачатрян, Л.Р. Исследование технологических режимов работы мембранного аппарата при концентрировании крахмального молока / Л.Р. Хачатрян, Р.В. Котляров, Б.А. Лобасенко // Техника и технология пищевых производств. – Кемерово. – 2015. – № 2(37). – С. 61–66.
Математическое моделирование является неотъемлемой частью современного научного исследования. Математическая модель представляет собой эффективный инструмент познания реального объекта и позволяет определить его характеристики, получить оценку показателей эффективности качества, осуществить поиск оптимальной структуры и параметров объекта.
Универсальным методом математического моделирования является использование передаточных функций, которые позволяют описывать динамические свойства системы практически любой структуры и сложности. Передаточная функция представляет собой динамический оператор, характеризующий прохождение сигналов через систему. Совокупность передаточных функций образует математическую модель системы, позволяющую оценивать и прогнозировать процессы, происходящие в ней.
Целью данной статьи является описание этапов моделирования процесса мембранного концентрирования крахмального молока в ультрафильтрационном аппарате с гидродинамическими элементами: анализ объекта моделирования, структурная и параметрическая идентификация математической модели, проверка адекватности модели.
Материалы и методы исследования
Объектом моделирования является процесс мембранного концентрирования крахмального молока в ультрафильтрационном аппарате с гидродинамическими элементами (рис. 1) [5]. Аппарат содержит корпус 1, выполненный в виде цилиндра. С одной стороны корпуса располагается патрубок 2 для подачи исходного потока, с другой – патрубок 3 для отвода конечного продукта в виде концентрата, а также патрубок 4 для отвода фильтрата. Внутри корпуса располагается коаксиально полупроницаемая мембрана 5. В мембране располагается вставка, состоящая из несущего стержня 6, гидродинамических элементов 7 и фиксирующих колец 8.
Рис. 1. Ультрафильтрационный аппарат: 1 – корпус; 2 – патрубок для подачи исходного потока; 3 – патрубок для отвода концентрата; 4 – патрубок для отвода фильтрата; 5 – фильтрующий элемент (мембрана); 6 – несущий стержень вставки; 7 – гидродинамические элементы вставки; 8 – фиксирующие кольца
Исходный раствор под давлением подается через патрубок 2 во внутренний канал мембраны 5. В процессе мембранного концентрирования на внутренней поверхности фильтрующего элемента формируется слой осадка, который снижает производительность аппарата. По ходу движения потока внутри мембраны, за счет уменьшающегося зазора между гидродинамичеким элементом 7 и внутренней стенкой мембраны 5, происходит увеличение его скорости. Благодаря этому осуществляется очистка поверхности фильтрующего элемента от образующихся отложений. Прошедший через мембрану фильтрат отводится через патрубок 4. Оставшаяся часть исходного раствора отводится через патрубок 3.
Результаты исследования и их обсуждение
На начальном этапе проведен анализ объекта моделирования [1, 2]. Структурная схема процесса мембранного концентрирования представлена на рис. 2. Входным воздействием является концентрация задерживаемых веществ в исходном растворе С(S), % мас. Выходным параметром является удельная производительность по фильтрату G(S), м3/(м2?с). Основным каналом типа «вход – выход» является канал «исходная концентрация С(S) – удельная производительность G(S)», передаточные свойства которого описывает передаточная функция WС(S).
На процесс также оказывают влияние основные параметры технологического режима: рабочее давление p(S), температура концентрируемого раствора T(S). Передаточные свойства каждого канала системы определяются соответствующей передаточной функцией: канал «температура – удельная производительность» характеризуется передаточной функцией WТ(S), канал «давление – удельная производительность» передаточной функцией WР(S).
Рис. 2. Структурная схема процесса мембранного концентрирования
На этапе структурной идентификации установлено, что изменение удельной производительности под воздействием исходной концентрации наиболее точно описывает зависимость (1) [3, 4].
(1)
где k – коэффициент усиления объекта, м3/(м2·с); Т – постоянная времени объекта, с; k0 – удельная производительность по фильтрату в начальный момент времени, м3/(м2·с).
Изменение удельной производительности под воздействием температуры и давления наиболее точно описывают зависимости (2) и (3) соответственно [3, 4].
WT(S) = kT; (2)
Wp(S) = kp, (3)
где kТ – коэффициент влияния температуры концентрируемого раствора на удельную производительность по фильтрату, м3/(м2?с??С); kp – коэффициент влияния рабочего давления на удельную производительность по фильтрату, м3/(м2?с?МПа).
Вектор входных воздействий системы определен в виде
X(S) = {C(S), T(S), p(S)}. (4)
Характер изменения входного воздействия имеет ступенчатый вид, при этом величина в момент времени t0 мгновенно возрастает от начального Xнач до конечного Xкон постоянного значения и далее остается неизменной [3]. Такому воздействию соответствует функция
(5)
На этапе параметрической идентификации по результатам экспериментальных исследований процесса мембранного концентрирования крахмального молока в ультрафильтрационном аппарате [6] установлены значения Xнач и Xкон входных воздействий
(6)
(7)
(8)
(9)
Для определения реакции системы на входные воздействия p1(t), p2(t), T(t) необходимо их совместное внесение в систему с воздействием С(t).
Экспериментально установлено, что коэффициент усиления k выражения (1) зависит от концентрации исходного сырья
(10)
Параметрическая идентификация модели процесса мембранного концентрирования в ультрафильтрационном аппарате проводна на основе экспериментальных зависимостей (разгонных характеристик), полученных при внесении в систему выбранных воздействий.
Передаточные свойства канала «исходная концентрация С(S) – удельная производительность G(S)» описывает передаточ ная функция
(11)
В диапазоне изменения давления от 0,15 до 0,25 МПа коэффициент влияния давления на удельную производительность составил
kp1 = –0,0009 % мас./(м2?с?МПа).
В диапазоне изменения давления от 0,25 до 0,3 МПа коэффициент влияния давления на удельную производительность составил
kp2 = 0,0012 % мас./(м2?с?МПа).
Влияние температуры на удельную производительность характеризуется коэффициентом усиления
kТ = 0,000007 % мас./(м2?с?°С).
Итоговая модель процесса мембранного концентрирования имеет вид
(12)
На заключительном этапе реализация математической модели программными средствами MatLAB позволила оценить ее адекватность экспериментальным данным. Расхождение экспериментальных и теоретических данных во всем диапазоне изменения технологических параметров процесса не превышает 4,8 %.
Таким образом, математическая модель процесса мембранного концентрирования в ультрафильтрационном аппарате с гидродинамическими элементами, разработанная на основе передаточных функций, с достаточной степенью адекватна экспериментальным данным и может быть использована для моделирования, расчета и прогнозирования работы мембранного оборудования.
Рецензенты:
Карташов В.Я., д.т.н., профессор кафедры прикладной математики, ФГБОУ ВПО «Кемеровский государственный университет», г. Кемерово;
Полтавцев В.И., д.т.н., профессор кафедры общенаучных дисциплин, ФГБОУ ВПО «Кемеровский государственный сельскохозяйственный институт», г. Кемерово.
Хачатрян Л.Р. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОЦЕССА МЕМБРАННОГО КОНЦЕНТРИРОВАНИЯ КРАХМАЛЬНОГО МОЛОКА НА ОСНОВЕ ПЕРЕДАТОЧНЫХ ФУНКЦИЙ // Фундаментальные исследования. – 2015. – № 11-5.
– С. 928-931;
URL: https://fundamental-research.ru/ru/article/view?id=39535 (дата обращения: 23.11.2024).
