D_2_S (Разработка и внедрение автоматизированных систем управления технологического оборудования минипекарень), страница 2

Nu = a ´ l / l - число Нуссельта.

Число Нуссельта характеризует собой условия теплопередачи между твердым телом и средой, оно содержит в себе искомую величину - коэффициент теплоотдачи a, коэффициент теплопроводности среды l и определяющий размер l, характеризующий собой геометрическое подобие.

Ре = u ´ l / a - число Пекле.

Число Пекле обычно преобразуется и представляется в виде двух критериев:

Число Рейнольдса Re содержит в себе скорость потока u и коэффициент кинематической вязкости n = m/r м²/с, где m - коэффициент динамической вязкости, характеризует собой ее внутреннее трение; r - плотность среды. Число Рейнольдса является критерием гидродинамического подобия, он характеризует собой условия вынужденного движения среды.

Множителями числа Прандтля Pr являются физические параметры - кинематическая вязкость и коэффициент температуропроводности - число Прандтля характеризует собой свойства среды. Оно практически не зависит ни от давления, ни от температуры. Так как коэффициент температуропроводности

a = l / (c ´ r ),

то Pr = c ´ r ´ n / l ,

где с - теплоемкость среды;

r - плотность среды;

l - коэффициент теплопроводности среды.

Так как мы имеем дело с теплоотдачей в потоке движущейся среды, то кроме теплового подобия, должны быть соблюдены условия гидромеханического подобия. Критерии гидромеханического подобия выделяются из дифференциального уравнения движения несжимаемой вязкой жидкости Навье-Стокса. Это то же число Рейнольдса, а также число Грасгофа:

Gr = b´g´l³´Dt/n²,

где g - ускорение свободного падения;

Dt - температурный перепад между средой и омываемой ею поверхностью;

b - функция, связывающая изменение плотности среды с температурой.

Число Грасгофа Gr характеризует свободное конвективное движение среды.

Критериальное уравнение теплопередачи конвекцией строится по типу:

Nu = f ( Re , Gr , Pr )

Здесь Nu содержит в себе искомую величину a и является неопределяющим критерием, тогда как критерии Re, Gr, Pr - определяющими.

Для газов одинаковой атомности и, в частности, для воздуха, когда Pr = const, будем иметь:

Nu = F ( Re , Gr ).

А при вынужденном турбулентном движении газа, что имеет место в расстойном шкафу при обтекании потоком воздуха нагревателей, когда естественной конвекцией можно пренебречь, выпадает число Грасгофа:

Nu = F ( Re ).

Значение критической скорости, при которой происходит переход от ламинарного режима течения воздуха к турбулентному, соответствующее числу Рейнольдса Re = 2200, равно:

u_кр = 2200 ´ n / d.

При работе расстойного шкафа в установившемся режиме в нем происходят постоянные колебания температуры в установленных пределах. Это объясняется работой системы управления. То есть не только при прогреве, но даже в установившемся режиме коэффициенты теплоотдачи поверхностей ТЭНов, тележек и стенок не являются постоянными и не подлежат однозначному точному математическому описанию.

Еще большую проблему представляет нахождение коэффициента теплоотдачи поверхностей тестовых заготовок. Это связано с тем, что при поступлении тестовых заготовок в расстойный шкаф они прогреваются значительно медленнее, чем циркулирующая в камере шкафа паровоздушная среда. Когда температура заготовок оказывается меньше температуры точки росы паровоздушной среды, на их поверхности конденсируется влага, многократно увеличивая коэффициент теплоотдачи и интенсифицируя процесс теплопередачи от паровоздушной среды к поверхности тестовых заготовок, в результате чего скорость прогрева их поверхности увеличивается. Влага, покрывающая тестовые заготовки, также предотвращает их от затвердевания и от образования трещин при увеличении тестовых заготовок. Конденсация влаги прекращается по достижении поверхностью тестовых заготовок температуры точки росы (в свою очередь зависящей от постоянно меняющихся температуры и влажности циркулирующего в камере расстойного шкафа воздуха). Коэффициент теплоотдачи поверхности тестовых заготовок при этом уменьшается, что влечет за собой уменьшение интенсивности их прогрева. Таким образом, строгое математическое описание коэффициента теплоотдачи поверхности тестовых заготовок не представляется возможным.

Модель поддержания заданной влажности воздуха

Относительная влажность воздуха в расстойном шкафу находится по уравнению:

j_возд = r_п / r_max ,

где r_max - максимально возможная абсолютная влажность воздуха при данной температуре;

r_п - действительная абсолютная влажность ненасыщенного воздуха, скорость изменения которой (dr_п/dt) может быть выражена как:

,

где V_возд - объем циркулирующего в расстойном шкафу влажного воздуха;

G_потерь - расход пара на конденсацию на стенках камеры расстойного шкафа и на поверхности тестовых заготовок;

G_пара - расход пара на увлажнение воздуха в камере расстойного шкафа:

G_пара = Р_{тен вл} / r ,

где r - теплота парообразования воды;

Р_{тен вл} - мощность ТЭНов, используемых для подогрева и испарения воды, с целью увлажнения воздуха в расстойном шкафу.

В связи с тем, что в процессе расстойки необходимо поддерживать заданную влажность, данные ТЭНы включены только пока относительная влажность воздуха в камере расстойного шкафа меньше заданной. Как только относительная влажность воздуха превышает заданный предел система управления подает сигнал на отключение ТЭНов увлажнения. При этом Р_{тэн вл} = 0. При падении относительной влажности ниже предельной, система управления подает сигнал на включение ТЭНов увлажнения. При этом Р_{тэн вл} = Р_{тэн вл зад}, где Р_{тэн вл зад} - номинальная мощность ТЭНов увлажнения.

Максимально возможная абсолютная влажность воздуха (r_max) зависит от температуры циркулирующего в камере расстойного шкафа воздуха, а теплота парообразования воды (r) зависит от температуры воды. И если последняя в установившемся режиме работы расстойного шкафа практически неизменна, то температура воздуха меняется в заданном диапазоне и в установившемся режиме работы расстойного шкафа. А r_max довольно существенно зависит от температуры воздуха. То есть даже в установившемся режиме работы расстойного шкафа r_max будет существенно меняться и эти изменения не описываются с большой точностью математически.

Потери пара на конденсацию (G_потерь) происходят не всегда, а только при условии, что внутренняя поверхность стенок камеры расстойного шкафа или поверхность тестовых заготовок имеют температуру меньшую, чем температура точки росы (t_р) при данных условиях.

Конденсацию пара на стенках можно практически предотвратить сделав достаточной теплоизоляцию стенок расстойного шкафа. Напротив, конденсация пара на поверхности тестовых заготовок является неотъемлемой частью технологического процесса расстойки тестовых заготовок, напрямую влияющей на качество готовой продукции, и происходит в первой половине процесса расстойки, до момента достижения поверхностью тестовых заготовок температуры точки росы. В свою очередь, температура точки росы зависит от влажности и температуры воздуха в камере расстойного шкафа. Таким образом, математическое описание потерь пара на конденсацию не представляется возможным.

Из всего вышеизложенного становится ясно, что полная математическая модель не пригодна для написания по ней алгоритма программы и самой программы для ЭВМ с целью моделирования процессов протекающих в расстойном шкафу и выбора параметров системы управления, удовлетворяющих заданным требованиям.