Принципы построения компьютерной модели

Использование ступенчатого метода расчёта [8] позволило разработать модель сушки плотного слоя зерна на базе компьютерной модели элементарного слоя [9].

При построении компьютерной модели СВЧ — конвективной сушки зернового слоя учитывали результаты исследований распределения микроволнового поля в СВЧ активной зоне. Поскольку удельная энергия, поглощаемая зерном, изменяется в зависимости от влажности зерна и расстояния до магнетрона, то был проведён двухфакторный эксперимент по определению такой зависимости [10]. Результаты эксперимента частично приведены в таблице 1.

Таблица 1 Экспериментальные данные по влиянию влажности зерна и расстояния до магнетрона на удельную энергию поглощаемую зерном.

С использованием полученных экспериментальных данных выполнен регрессионный анализ и получена следующая регрессионная модель:

. (1).

Регрессионный анализ проводили с использованием пакета прикладных программ MATLAB.

Коэффициенты модели значимы по критерию Стьюдента. Адекватность модели проверялась по критерию Фишера. Точность модели оценивается следующими показателями: SSE: 1.409e+06; R-square: 0.9731; Adjusted R-square: 0.9698; RMSE: 185.4, что характеризует хороший уровень достоверности.

Уравнение регрессии (1) использовано в компьютерной модели, что позволяет моделировать изменение величины удельной энергии микроволнового поля в зерновом слое в процессе его сушки.

Было принято, что толщина зернового слоя, который подвергается СВЧ — конвективному воздействию, составляет 15 см. Это сопоставимо с глубиной проникновения поля СВЧ в зерновой слой пшеницы кондиционной влажности (14%). Приняли допущение, что зерновой слой может быть разделён на участки по 5 см толщиной, в пределах которых микроволновое поле распределено равномерно. Компьютерная Simulink — модель зернового слоя толщиной 15 см, разделённого на три слоя по 5 см, приведена на рис. 1.

Рис. 1 — Компьютерная модель тепло и — влагообмена в зерновом слое при СВЧ — конвективном воздействии

Зерновой слой представлен тремя блоками-слоями «Слой зерна», в каждом из которых поле СВЧ имеет разную напряжённость, по аналогии с реальным зерновым слоем. Агент сушки последовательно проходит через все слои. Его входные параметры задают при помощи блоков Tvh и Fvh.

С помощью блока Relay управляют включением/выключением магнетронов. Магнетроны отключают, когда температура зерна, в точке наибольшей мощности поля СВЧ, достигает 55^оС. На осциллографы выводятся графики изменения температуры зерна в каждом слое (блок Q), влажности зерна в слое (блок W), относительной влажности воздуха на выходе из каждого слоя зерна (Fvih), температуры агента сушки на выходе из каждого зернового слоя (Tvih) [11].

При разработке конструкции СВЧ — конвективных зон большое значение имеет направление движения агента сушки, его температура, равномерность распределения в зерновом слое. Исходя из этого, проводили моделирование процесса СВЧ — конвективной сушки зерна. На рис. 2 показана условная схема расположения слоёв зерна, распространения микроволнового поля и движения агента сушки.

а) первый вариант воздухораспределения.

б) второй вариант воздухораспределения Рис. 2 — Схема расположения зерновых слоёв и распределения в них поля СВЧ и агента сушки

При первом варианте движения агента сушки сначала он проходит через Слой 1 зерна. Поскольку данный слой ближе других расположен к источнику микроволнового поля, то он должен иметь большую температуру, чем Слой 2 и Слой 3. Поэтому, при прохождении воздуха последовательно Слой 1> Слой 2 > Слой 3 он должен нагреваться от зерна в Слое 1. Затем его температура может, как уменьшаться, так и поддерживаться на постоянном уровне. Всё зависит от конкретных режимов сушки зерна.

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

Результаты моделирования. На рис. 3 показаны результаты моделирования изменения температуры зерна по слоям. Моделирование проводили при следующих исходных параметрах:;; ^oC;. При относительной влажности воздуха равновесная влажность зерна составит 14%. Такое значение относительной влажности воздуха принято для того, чтобы оценить, как будет изменяться сушка зерна в слоях дальних от магнетрона.

а) первый вариант воздухораспределения.

б) второй вариант воздухораспределения Рис. 3 — Изменение температуры зерна по слоям в процессе СВЧ — конвективной обработки

Из графиков видно, что схема воздухораспределения существенно влияет на изменение температуры зерна в СВЧ — конвективной зоне. При первом варианте воздухораспределения в начале процесса сушки наблюдается существенная неравномерность нагрева зерна по слоям. Затем неравномерность нагрева зерна уменьшается, но в пределах 3 ^оС остаётся до высыхания зерна. Максимальная температура не превышает 55 ^оС. Это позволяет говорить, что процесс СВЧ — конвективной сушки достаточно управляем для того, чтобы обеспечить сохранность технологических качеств зерна.

При втором варианте воздухораспределения неравномерность нагрева зерна по зонам ещё выше, чем при первом варианте. Причём неравномерность нагрева зерна между первым Слоем и Слоями 2 и 3 изменяется в процессе сушки незначительно. Слой 3 в меньшей степени подвергается воздействию микроволнового поля, от которого температура зерна значительно не повышается. Плюс на него постоянно действует агент сушки при температуре 20^oC. Поэтому температура Слоя 3 не превышает 35^oC. Подогретый зерном в Слое 3 воздух попадает в Слой 2. Удельная мощность микроволнового поля в Слое 2 также выше, чем Слое 3. Поэтому температура нагрева зерна в нём уже достигает 37…39^oC. В Слое 1 зерна выделяется максимальная мощность поля СВЧ, поэтому магнетрон отключают, когда температура зерна достигает 55 ^oC. Да и воздух в Слой 1 подаётся уже подогретым из Слоя 2, поэтому колебания температуры нагрева зерна в Слое 1 минимальны.

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

а) первый вариант воздухораспределения.

б) второй вариант воздухораспределения Рис. 4 — Кривые сушки зерна по слоям

Скорость сушки существенно отличается по слоям зерна. Причём максимальная удельная мощность поля СВЧ, поглощаемая зерновым слоем, совершенно не обеспечивает максимальной скорости сушки. Если первый слой не высыхает до кондиционной влажности 14% и за 4 часа, то второй высыхает за 1,5 часа, а третий — за 1,3 часа. При этом средняя влажность всех слоёв достигает 14% за 1,7 часа. Неравномерность влажности по слоям составляет 3%. компьютерный зерновой магнетрон При втором варианте воздухораспределения Слой 1 высыхает до влажности 14% за 1,4 часа. Ещё 0,3 часа он продолжает сушиться, а затем начинает увлажняться. Слой 3 высыхает до кондиционной влажности за 2,5 часа. В Слое 2 на первом этапе происходит увлажнение зерна, а затем сушка до кондиционной влажности за 4 часа. Такая система воздухораспределения даёт достаточно большую неравномерность высыхания зерна по слоям и увеличивает время сушки зерна. Так средняя влажность зерна по слою достигнет 14% только за 4 часа, что более чем в два раза больше, чем при первом варианте воздухораспределения.