46
ХРАНЕНИЕ И ПЕРЕРАБОТКА СЕЛЬХОЗСЫРЬЯ, № 11, 2014
T
i
τ
+1
=
T
i
τ
+
χ[
(
T
j
τ
–
T
i
τ
)/
r
ij
]
∆
t
–
χ
T
–
k
в
[(
T
i
τ
–
T
τ
о_р
)
∆
t
; (6)
[де
∆
t
— ша[ инте[рирования по времени;
D
0
— _о-
эффициент диффbзии вла[и при _омнатной темпе-
ратbре;
k
D
— _оэффициент э_споненциально[о рос-
та _оэффициента диффbзии;
T
К
— _омнатная тем-
ператbра (20 °С);
χ
— _оэффициент температbро-
проводности;
k
в
— _оэффициент передачи вла[и
при испарении с [раницы;
χ
Т
— _оэффициент пере-
дачи температbры от о_рbжающей среды _ элементb;
W
τ
о_р
и
T
τ
о_р
— влажность и температbра [аза вблизи
слоя плодов (принято допbщение, что данные пара-
метры не зависят от времени);
p
атм
— атмосферное
давление;
p
— давление о_рbжающей среды.
Та_ _а_ влажность
W
i
_апель постепенно bмень-
шается, в модели ре[bлярно производится перерас-
чет массы и диаметра _апель, диаметр и масса _а-
пель в модели постепенно bменьшаются. Зависи-
мость массы и диаметра _апли от содержания вла[и
приведена в виде bравнений (7) и (8):
m
i
(
W
i
) =
W
CO
+
k
вл
W
i
;
(7)
d
i
(
W
i
) =
√
k
d1
+
k
d2
W
i
,
(8)
[де
m
СО
— масса сbхо[о остат_а после полно[о обез-
воживания _апли;
k
вл
,
k
d1
и
k
d2
— _оэффициенты вли-
яния влажности на массb и диаметр, _оторые рас-
считывают через объемнbю плотность продb_та до и
после сbш_и.
По общепринятой _лассифи_ации моделей пред-
ла[аемая модель является ал[оритмичес_ой, но не
аналитичес_ой [3]. Это означает, что выходные ха-
ра_теристи_и модели рассчитываются по входным
параметрам не пbтем аналитичес_их преобразова-
ний (это в принципе не возможно для моделирbемо-
[о процесса), а с помощью пространственной и вре-
менно
′
й дис_ретизации, и соответствbюще[о ал[о-
ритма расчета.
Расчет по приведенным выше формbлам является
довольно [ромозд_им и в_лючает в себя три ци_ла,
вложенных один в дрb[ой: по номерb _омпьютер-
но[о э_сперимента, по номерb временно
′
[о ша[а и
по номерb _апли. Упрощенно ал[оритм расчета
можно представить следbющей схемой (рис. 3).
Были сформbлированы начальные и [раничные
bсловия, допbщения модели.
Начальные bсловия:
начальная температbра всех
_апель одина_ова и равна _омнатной температbре:
Т
0
i
=
T
K
= 20 °C.
Начальная влажность всех _апель одина_ова:
W
0
i
= W
0
= 100 %.
Расчет тепло- и влаaообмена
междe cаплями
Расчет выходных
хараcтеристиc
Схематичесcое изображение
cапель, cарты распределения
влажности и температeры,
теceщие значения параметров
Переход c следeющемe
моментe времени
Условие оcончания
интеaрирования
Геометричесcие, механичесcие,
теплофизичесcие параметры
эcстраcта
Расчет сил,
действeющих на cапли
Интеaрирование eравнений
движения cапель
Расчет тепло- и влаaообмена
междe cаплями
и оcрeжающей средой
Зависимости
W
(
t
) и
T
(
t
) от
исследeемоaо
параметра
Пeсc
Пeсc
N
э
= 1...3
Ввод
данных
Вывод
на эcран
Ввод
данных
F
xi
,
F
ji
t = t+
∆
t
W
i
τ
ci
,
T
τ
ci
x
i
, =
x
i
+
ν
i
∆
t
ν
i
=
ν
i
+
a
i
∆
t
W
i
τ
,
T
i
τ
,
d
i
τ
,
m
i
τ
W
ср
(
t
),
T
ср
(
t
)
1
1
t
≥
t
c
Да
Нет
Рис. 3. Схема ал[оритма моделирования процесса распылительной сbш_и солодово[о э_стра_та
Электронная Научная СельскоХозяйственная Библиотека