close

Вход

Забыли?

вход по аккаунту

?

Математическая модель энергосберегающей технологии сушки зерна..pdf

код для вставкиСкачать
УДК 51-74: 631.365.32
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ЭНЕРГОСБЕРЕГАЮЩЕЙ ТЕХНОЛОГИИ
СУШКИ ЗЕРНА
© 2015 г. А.Н. Васильев, О.В. Северинов
Снижение энергоѐмкости процесса сушки зерна возможно за счѐт совершенствования конструкций
оборудования, технологического процесса и за счѐт использования систем оптимального управления
процессом. Совершенствование управления процессом сушки требует наличия математической модели
объекта управления, на базе которой может быть построена модель оптимального управления. Поскольку
сушка – процесс тепло- и влагообмена между зерном и сушильным агентом, то представлены данные
закономерности и выполнены соответствующие преобразования для построения систем управления. Для
промежуточного решения дифференциальных уравнений в частных производных использовался метод
интегрального преобразования Лапласа. В результате получена такая запись дифференциального уравнения,
которая позволяет не выполнять сложные преобразования для записи результата в аналитическом виде, а
представить его в виде набора передаточных функций. Такое представление дифференциальных уравнений
позволяет разрабатывать компьютерные имитационные модели и получать численные решения при любых
начальных условиях. В зависимости от того, какая стратегическая задача стоит (необходимо сушить с
большей производительностью или с меньшими энергозатратами), определяют требуемую
производительность вентилятора и схему перемещения зерна по зонам бункера. Для решения этих задач
была разработана математическая модель. В основу модели положено уравнение теплового баланса для
элементарного слоя. Полная модель процесса сушки зерна в установках активного вентилирования с
системой оптимального управления построена с использованием модели сушки плотного слоя и модели
сушки элементарного слоя и позволяет снизить энергоѐмкость процесса. Предполагается, что в бункере
имеется возможность выпускать зерно из определѐнной зоны.
Ключевые слова: слой зерна; сушка зерна; тепло- и влагообмен; математическая модель;
дифференциальное уравнение; преобразование Лапласа; оптимальное управление.
It is possible to reduce the energy intensity of grain drying process by improving equipment design,
technological process and using systems of optimum process control. Improving the management of the drying
process requires a mathematical model of the controlled object basing on which there can be designed optimal
control model. As the drying process is the process of heat and moisture exchange between the grain and the drying
agent, there are presented these regularities and performed the appropriate conversion for constructing management
systems. For intermediate solutions of differential equations in partial derivatives of the integral method there was
used Laplace transformation. As a result, we obtain such a record of this differential equation that allows not to
conduct complex transformations to record result in an analytical form but to present it in the form of transfer
functions set. Such a notion of differential equations enables the development of computer simulation models and
numerical solutions for any initial conditions. Depending on what strategic task (it is necessary to dry with higher
performance or using less energy) there is determined the required fan capacity and the circuit for moving the grain
through the hopper zones. To solve these problems there was developed a mathematical model. The model is based
on the heat balance equation for the elementary layer. Complete grain drying model in active ventilation units with
optimal control system was built using drying dense layer model and the model of the drying unit layer and allows
to lower the energy intensity of the process. It is assumed that in the bunker it is possible to release grain from a
specific area.
Key words: grain layer; grain drying; heat- and moistureexchange; mathematical model; differential equation;
Laplace transformation; optimal control.
Введение и цель работы. Снижение
энергоѐмкости процесса сушки зерна
возможно за счѐт совершенствования
конструкций
оборудования,
технологического процесса и за счѐт
использования
систем
оптимального
управления
процессом
[1].
Совершенствование управления процессом
сушки требует наличия математической
модели объекта управления, на базе
которой может быть построена модель
оптимального
управления. Поскольку
сушка – процесс тепло- и влагообмена
между зерном и сушильным агентом, то
необходимо
рассмотреть
данные
закономерности и представить их в виде,
удобном
для
построения
систем
управления.
Описание математической модели
и ее обсуждение. Описание процесса
тепло- и влагообмена в элементарном
зерновом слое может быть представлено
следующей системой уравнений [2, 3, 4]:


ln( 1  F )
Wp  

6
 5,47  10  (T  273) 
0, 435
;
K  7,1  10 2  e0, 05T ;
dW
 KW  KW p ;
d
 с d
 з r
dT
dT
dW
;
 3600 V 
 з з

d
dR
 в св d 100   в св d
 с dT
 с dT
d
r  dW
;
 в в
 3600 V в в

d
 з с з d
 з с з dR c з 100 d
10 з dW
dD
dD
 3600  V

;
d
dR
 в d
F
745  D
(2)
(3)
(4)
(5)
(6)
,
(7)
7,5  T
T
238
(5).
После
этого
определяют
влагосодержание D (6) и относительную
влажность агента сушки на выходе из
зернового слоя F (7). Последовательное
выполнение
расчетов
для
каждого
элементарного слоя позволит иметь
полную картину динамики сушки зерна.
Необходимо учитывать, что величина
скорости V сушильного агента может быть
изменена на любом слое и на любом этапе
расчета.
Это позволит достаточно точно выполнять
расчет сушки зерна в установках с
радиальным распределением воздуха, при
уравнении расходом вентилятора.
Современное
прикладное
программное
обеспечение
позволяет
построить имитационную модель процесса
сушки зерна в платном слое на базе модели
элементарного слоя. Для построения
имитационной модели элементарного слоя
представим дифференциальные уравнения
(3), (4), (5), (6) в виде передаточных
функций.
Для уравнения (3) воспользуемся
(622  D) 100,622 
где T – температура агента сушки, °C;
D – влагосодержание сушильного
агента, г/кг;
W – текущая влажность зерна, %;
θ – температура зерна, °С;
V – скорость агента сушки, м/с;
св, сз – теплоемкость воздуха и зерна,
кДж/кг·°С;
ε – порозность зернового слоя;
 ′ – скрытая теплота парообразования
воды, кДж/кг;
γз – объемная масса зерна, кг/м3;
γв – удельный вес воздуха, кг/м3;
K – коэффициент сушки, 1/ч;
Wp – равновесная влажность зерна, %;
R – радиус зерновки, м;
τ – время, ч,
F – относительная влажность воздуха, %.
Уравнения
записаны
в
последовательности их применения. По
параметрам воздуха, подающегося в
зерновой
слой,
рассчитывают
его
равновесную влажность Wp (1). Затем
определяют коэффициент сушки К (2).
Уравнение
коэффициента
сушки
единичного зерна получено Г.С. Окунем
[5]. Далее рассчитывают влажность зерна
W (3). Потом температуру Т, которую
приобретает
агент
сушки
после
прохождения единичного слоя зерна (4).
Затем рассчитывают температуру зерна θ
(1)
d
 p и получим
d
pw( p )  Kw( p )  Kw p ( p ) ,
правилом замены [28]
откуда Ww w ( p ) 
p
K
, где
pK
Wwp w ( p) –

LT ( x,  )  T ( x, p)  e  p T ( x,  )d ,
0
передаточная
функция
по
каналу
можно интегрированием по частям
равновесная влажность – влажность зерна.
установить соотношение:
Уравнения (3), (4), (5) содержат
 dT 
частные
производные.
Для
L    pT ( x, p)  T ( x,0) ,
промежуточного
их
решения
 d 
воспользуемся методом интегрального
 dT  dT
L  
( x, p ) ,
преобразования Лапласа. В этом случае
 d  dx
преобразование частных производных
– начальное значение
T (x,0)
осуществляется по следующему правилу где
[6]: если T  T ( x, ) и преобразование температуры.
Выполним преобразование Лапласа
Лапласа производится по переменной
для (3) по τ:
 (  0) , то, обозначив
 c
 з  r
dT
pT ( x, p)  T ( x,0)  3600  V
  з з p ( x, p) 
pW ( x, p) ,
dR
 в cв
100   в cв
или
dT
pT ( x, p)  T ( x,0)  3600 V
  A1 p ( x, p)  A2 pW ( x, p) ,
dR
Решим дифференциальное уравнение
 з cз
 з  r
где A1 
, A2 
.
относительно R:
 c
100   c
в в
в в
dT
3600 V
  A1 p ( x, p)  pT ( x, p)  T ( x,0)  A2 pW ( x, p) ,
dR
при этом р играет роль переменной. С
Решение уравнения с использованием
учетом, что T ( x,0)  T0 , T0 – температура пакета MATLAB [7] при начальных
T (0)  T0 дало
условиях
следующий
воздуха на входе в зерновой слой.
результат:
T ( p)  T0 (
1 1  p 1
 e )  A1 p ( p)  A2 pW ( p) ,
p p
R
.
3600  V
Уравнение (8) представляет собой
зависимость изображения температуры
теплоносителя T(р) от изображений
температуры агента сушки на входе в
зерновой слой Т0(р), от изображения
текущей температуры зерна θ(р), от
изображения текущей влажности зерна
W(p). Такая запись дифференциального
уравнения
позволяет
не
выполнять
сложные преобразования для записи
результата в аналитическом виде, а
представить
его
в
виде
набора
передаточных функций. Решение будет
представлено в виде трѐх блоков: входное
воздействие
Т0(р)
и
параллельное
где  1 
(8)
соединение двух ветвей – интегрирующего
1
звена
и интегрирующего звена с звеном
p
1
запаздывания e  p1 ; входное воздействие
p
θ(р) и последовательное соединение
усилительного
звена
A1
и
дифференцирующего звена р; входное
воздействие W(p) и последовательное
соединение усилительного звена A2 и
дифференцирующего звена p. Такое
представление
дифференциальных
уравнений
позволяет
разрабатывать
компьютерные имитационные модели и
получать численные решения при любых
начальных условиях.
Аналогичным
способом
осуществлялось
преобразование
дифференциальных уравнений (5) и (6). В
результате проведенных преобразований
получена полная система уравнений и
передаточных функций, позволяющих
выполнять расчет тепло- и влагообмена в
элементарном слое зерна.
В укрупненном виде структурная
схема расчета процесса тепло- и
влагообмена в единичном зерновом слое
приведена на рисунке 1.
Рисунок 1 – Структурная схема сушки элементарного слоя
Энергосберегающая
технология
сушки
предусматривает
определение
параметров зернового слоя в определѐнной
зоне бункера активного вентилирования,
определение энергозатрат на сушку слоя
зерна при известных параметрах агента
сушки
(относительная
влажность,
температура, скорость). В зависимости от
того, какая стратегическая задача стоит
(необходимо
сушить
с
большей
производительностью или с меньшими
энергозатратами), определяют требуемую
производительность вентилятора и схему
перемещения зерна по зонам бункера. Для
решения этих задач была разработана
математическая модель. В основу модели
положено уравнение теплового баланса для
элементарного слоя [8]:
L( I1  I 2 )  QН  QИ  QП ,
где L – массовый расход агента сушки,
кг/с;
L  SVp ,
S – площадь поперечного сечения участка
элементарного слоя (рисунок 2), м2;
V – скорость агента сушки, м/с;
p B – плотность агента сушки, кг/м3;
 – порозность зернового слоя, о.е.;
I1
–
энтальпия
агента
сушки,
поступающего в зерновой слой,
кДж/кг;
I 2 – энтальпия агента сушки на выходе из
зернового слоя, кДж/кг;
QН – количество теплоты, расходуемой на
d
нагрев зерна, QН   3с3
кДж;
d
QИ – количество теплоты, расходуемой на
 r  dW
испарение влаги, QИ  3
кДж;
100 d
QП – количество теплоты, теряемой в
процессе сушки, кДж.
Рисунок 2 – Схематичное представление части элементарного слоя зерна площадью 1 м2,
используемого для расчѐтов
Энтальпия агента сушки определится
как [9]:
I  (0,579  D  cвT )  4,19 ,
cв  (0,24  0,00047  D)  4,19
где
теплоѐмкость воздуха.
–
Тогда I1  I 2  0,597  ( D1  D2 )  1,0056 (T1  T2 )  0,00197 ( D1T1  D2T2 ) ,
При
допущении,
что
где T1 , T2 – температура агента сушки
соответственно на входе и производительность подачи агента сушки
регулируют таким образом, что вся энергия
выходе зернового слоя;
D1 , D2 – влагосодержание агента сушки агента сушки используется полностью,
получено базовое уравнение сушки с
соответственно на входе и учѐтом энергетического потенциала агента
выходе зернового слоя.
сушки:
SV p в   0,597  ( D1  D2 )  1,0056 (T1  T2 )  0,00197 ( D1T1  D2T2 ) 
 r  dW
d
 n  d пр  S   з c з
 n  d пр  S  з
,
d
100 d
(9)
где dпр – приведѐнный диаметр зерновки, м;
n – количество элементарных слоѐв.
Основным показателем процесса
сушки является влажность зернового слоя
W. Запишем уравнение (9) относительно W:
vpв  100
100 c з d
dW

 0,597  ( D1  D2 )  1,0056 (T1  T2 )  0,00197 ( D1T1  D2T2 ) 
.
d  з r   d пр  n
r  d
После преобразований уравнение примет следующий вид:
dW V
d
 K1  D1 (0,00197  T1  0,597)  1,0056  T1  D2 (0,597  0,00197  T2 )  1,0056  T2   K 2
,
d
n
d
свойства влажности зерна по каналам
100  c з
v  pв  100 V
 K2 .
 K1 ,
где
«параметры агента сушки» и температуру
 з r   d пр  n n
r
зерна, которые используем в полной
Воспользовавшись
принципом модели сушки зерна (рисунок 3).
суперпозиции, получим передаточные
функции, описывающие динамические
Рисунок 3 – Структурная схема полной модели процесса сушки зерна
в секционных установках активного вентилирования
Полная модель процесса сушки зерна
в установках активного вентилирования
строится с использованием модели сушки
плотного слоя с использованием модели
сушки элементарного слоя (рисунок 1).
Предполагается, что в бункере имеется
возможность
выпускать
зерно
из
определѐнной зоны.
Схема работает следующим образом.
Начальные параметры зернового слоя W0,
θ0 задают для каждой зоны бункера
активного
вентилирования.
Задают
начальные
параметры
атмосферного
воздуха T0, D0, F0 и начальную скорость
агента сушки V0. Параметры атмосферного
воздух
изменяют
в
процессе
моделирования сушки в соответствии со
среднесуточными графиками изменения
температуры и влажности. Параметры
агента сушки на выходе из зоны T2, D2, F2 и
T3, D3, F3 являются входными для
следующей зоны. Скорость агента сушки в
каждой
зоне
задаѐтся
отдельно
(уменьшается) в зависимости от размеров
зоны и типа зерна. Для задания скорости
агента сушки в каждой зоне входную
скорость
Vo
умножают
на
соответствующий коэффициент K3, K4.
В
процессе
моделирования
параметры воздуха на входе в зону T0, D0,
Vo и параметры агента сушки на выходе из
зоны T2, D2 подают в блок энергетического
баланса. Сюда же подают данные о
средней величине температуры зерна по
слою в зоне θср. С использованием
полученных
данных
в
блоке
энергетического баланса выполняется
расчѐт энергозатрат Эн на сушку зерна при
данной скорости агента сушки в первой
зоне бункера активного вентилирования и
время τс, в течение которого зерно
достигнет
требуемой
влажности. В
зависимости от того, по какому критерию
ведут управление процессом сушки
(минимум энергозатрат или максимум
производительности)
вносится
корректировка в скорость агента сушки.
Расчѐт
скорости
агента
сушки,
обеспечивающей
оптимальные
энергозатраты, или скорость сушки
осуществляет блок оптимизации (БОП).
Блок энергетического баланса с
блоком
оптимизации
получают
информацию и с других двух зон бункера
активного вентилирования. В зависимости
от стратегических задач, стоящих перед
зернопунктом, может быть выбрана схема
сушки зерна в зонах бункера. Это может
быть перемещение зерна из зоны в зону по
мере его подсушки или заполнение зон
свежим зерном. С учѐтом данных задач
блок оптимального управления будет
определять
требуемую
подачу
вентилятора.
Выводы
1.
Представленная
система
алгебраических
и
дифференциальных
уравнений позволяет разработать модель
сушки элементарного зернового слоя.
2.
Полученные
с
помощью
преобразования
Лапласа
дифференциальные
уравнения
в
операторном виде позволяют использовать
для построения модели зернового слоя
передаточные функции.
3. Разработанная модель зернового
слоя в виде передаточных функций
позволяет моделировать процесс тепло – и
влагообмена
при
изменяющихся
параметрах агента сушки.
4. Полная модель сушки зерна
позволяет
организовать
управление
процессом с учѐтом стратегических задач:
максимальная производительность или
минимальная энергоѐмкость.
Литература
1. Васильев, А.Н. Электротехнология
и управление при интенсификации сушки
зерна
активным
вентилированием
/ А.Н. Васильев. – Ростов-на-Дону: ТерраПринт, 2008. – 240 с.
2.
Васильев,
А.Н.
Задачи
совершенствования управления процессом
активного вентилирования зерна / А.Н.
Васильев, О.В. Северинов, А.С. Галов //
Инновации в сельском хозяйстве. – 2014. –
№ 1(6) http://ej.viesh.ru/wp-content/uploads
/2014/ 04/insel6.pdf (электронный журнал).
3. Васильев, А.Н. Моделирование
управления активным вентилированием
зерна / А.Н. Васильев, О.В. Северинов,
А.С. Галов // Инновации в сельском
хозяйстве.
–
2014.
–
№
2(7)
http://ej.viesh.ru/wpcontent/uploads/2014/04/insel6.pdf
электронный журнал).
4.
Васильев,
А.Н.
Вариант
определения
динамических
свойств
зернового
слоя
/ А.Н. Васильев, О.В. Северинов //
Энергообеспечение и энергосбережение в
сельском
хозяйстве:
труды
9-й
Международной
научно-технической
конференции (21–22 мая 2014 г.) в 5 ч. Ч. 5:
Инфокоммуникационные технологии и
нанотехнологии. – Москва: ГНУ ВИЭСХ,
2014. – С. 224–228.
5. Окунь, Г.С. Методы расчѐта
продолжительности сушки отдельного
зерна пшеницы и зернового слоя / Г.С.
Окунь // Машины для послеуборочной
поточной обработки семян. Теория и
расчѐт машин, технология и автоматизация
процессов; под ред. З.Л. Тица. – Москва:
Машиностроение, 1967. – С. 290–308.
6.
Мартинсон,
Л.К.
Дифференциальные
уравнения
математической физики: учебник для вузов
/
Л.К.
Мартинсон,
Ю.И. Малов; под ред. В.С. Зарубина,
А.П. Крищенко. – Изд. 4-е, стер. – Москва:
Изд-во МГТУ им. Н.Э. Баумана, 2011. –
367 с.
7. Ануфриев, И.Е. MATLAB 7
/ И.Е. Ануфриев, А.Б. Смирнов, Е.Н. Смирнова. – Санкт-Петербург: БХВ-Петербург,
2005. – 1104 с.
8. Гинзбург, А.С. Сушка пищевых
продуктов в кипящем слое / А.С. Гинзбург,
В.А. Резчиков. – Москва: Пищевая
промышленность, 1966. – 196 с.
9. Гинзбург, А.С. Основы теории и
техники сушки пищевых продуктов
/ А.С. Гинзбург. – Москва: Пищевая
промышленность, 1973. – 527 с.
References
1. Vasil'ev A.N. Jelektrotehnologija i
upravlenie pri intensifikacii sushki zerna
aktivnym ventilirovaniem [Electrotechnology
and managing at an intensification of grain
drying by active ventilation], Rostov-naDonu, Terra-Print, 2008, 240 p.
2. Vasil'ev A.N., Severinov O.V., Galov A.S. Zadachi sovershenstvovanija upravlenija processom aktivnogo ventilirovanija
zerna [The objectives of improving the
process control of grain active ventilation],
Teoreticheskij
i
nauchno-prakticheskij
zhurnal: Innovacii v sel'skom hozjajstve, 2014,
No. 1(6), available at: http://ej.viesh.ru/wpcontent/uploads/2014/04/insel6.pdf.
3. Vasil'ev A.N., Severinov O.V., Galov A.S. Modelirovanie upravlenija aktivnym ventilirovaniem zerna [Simulation
of managing grain active ventilation],
Innovacii v sel'skom hozjajstve, 2014, No.
2(7), available at: http://ej.viesh.ru/wp-content
/uploads/2014/04/insel6.pdf.
4. Vasil'ev A.N., Severinov O.V.
Variant opredelenija dinamicheskih svojstv
zernovogo sloja [Option of determining the
dynamic properties of the grain layer]
Jenergoobespechenie i jenergosberezhenie v
sel'skom
hozjajstve.
Trudy
9-j
Mezhdunarodnoj
nauchno-tehnicheskoj
konferencii (21-22 maja 2014 g.) v 5-ti
chastjah. Chast' 5. Infokommuni-kacionnye
tehnologii i nano-tehnologii, Moscow, GNU
VIJeSH, 2014, pp. 224–228.
5. Okun' G.S. Metody raschjota
prodolzhitel'nosti sushki otdel'nogo zerna
pshenicy i zernovogo sloja. Mashiny dlja
posleuborochnoj potochnoj obrabotki semjan.
Teorija i raschjot mashin, tehnologija i
avtomatizacija processov [Methods for
calculating the duration of a separate drying of
wheat grain and grain layer. Machinery for
postharvest flow seeds treatment. Theory and
calculation of machinery, technology and
process
automation],
Moscow,
Mashinostroenie, 1967, pp. 290–308.
6. Martinson L.K., Malov Ju.I.,
Differencial'nye uravnenija matematicheskoj
fiziki: uchebnik dlja vuzov [Differential
equations of mathematical physics: a textbook
for high schools], Izd. 4-e, ster., Moscow, Izdvo MGTU im. N.Je. Baumana, 2011, 367 p.
7. Anufriev I.E., Smirnov A.B., Smirnova
E.N. MATLAB 7 [MATLAB 7], SanktPeterburg, BHV-Peterburg, 2005, 1104 p.
8. Ginzburg A.S., Rezchikov V.A.
Sushka pishhevyh produktov v kipjashhem
sloe [Drying foodstuffs in boiling layer],
Moscow, Pishhevaja promyshlennost', 1966,
196 p.
9. Ginzburg A.S. Osnovy teorii i tehniki
sushki pishhevyh produktov [Fundamentals of
the theory and technique of drying foods],
Moscow, Pishhevaja promyshlennost', 1973,
527 p.
Сведения об авторах
Васильев Алексей Николаевич – доктор техн. наук, профессор, заместитель директора
по научной работе, ФГБНУ «Всероссийский научно-исследовательский институт
электрификации сельского хозяйства» (г. Москва, Россия). Тел: 8-495-709-34-46.
E-mail: vasilev-viesh@inbox.ru.
Северинов Олег Викторович – инженер, ФГБНУ «Всероссийский научноисследовательский институт электрификации сельского хозяйства» (г. Москва, Россия).
Е-mail: viesh@dol.ru.
Information about the authors
Vasiliev Alexey Nikolaevich – Doctor of Technical Sciences, professor, deputy director,
FSBSI «All-Russian Research Institute for Electrification of Agriculture» (Moscow, Russia).
Phone: 8-495-709-34-46. E-mail: vasilev-viesh@inbox.ru.
Severinov Oleg Viktorovich – engineer, FSBSI «All-Russian Research Institute for
Electrification of Agriculture» (Moscow, Russia). E-mail: viesh@dol.ru.
Документ
Категория
Без категории
Просмотров
20
Размер файла
859 Кб
Теги
энергосберегающие, зерна, технология, математические, pdf, модель, сушка
1/--страниц
Пожаловаться на содержимое документа