Разработка и внедрение автоматизированных систем управления технологического оборудования минипекарень
Разработка и внедрение автоматизированных систем управления технологического оборудования минипекарень
Разработка и идентификация упрощенной математической модели процессов
в расстойном шкафу
Принятые упрощения и допущения
В формулах конвекционной теплопередачи присутствуют коэффициенты
теплоотдачи (. Как было показано ранее, коэффициенты теплоотдачи зависят от
многих факторов: от температур поверхности и омывающей ее среды, скорости
движения последней, ее теплопроводности, вязкости, плотности и
теплоемкости, от конфигурации и состояния поверхности и омывающей ее среды.
В связи с невозможностью математического описания данных коэффициентов, для
их нахождения пользуются экспериментальными данными, широко используя
теорию подобия, позволяющую в известной степени обобщить полученные опытные
результаты. Но используемые для нахождения коэффициентов теплоотдачи
критериальные уравнения содержат критерии подобия (Nu, Pe, Re, Pr, Gr),
которые зависят от многих параметров поверхностей и омывающей их среды,
некоторые из которых зависят от температуры среды и от разности между ней и
температур омываемых ею поверхностей. Данные зависимости не описаны
математически. Конденсация влаги на поверхности тестовых заготовок в
процессе их расстойки еще больше затрудняет точное нахождение коэффициента
теплоотдачи их поверхности.
Конденсация влаги на поверхности тестовых заготовок, а также на
внутренней поверхности стенок камеры расстойного шкафа приводит к
уменьшению абсолютной и относительной влажности в камере расстойного шкафа.
Для поддержания заданной относительной влажности воздуха применяется
испарение воды, контролируемое проектируемой системой управления. Но вместе
с паром в камеру расстойного шкафа попадает дополнительная энергия.
Конденсацию влаги на внутренней поверхности стенок камеры расстойного шкафа
можно свести к минимуму путем их лучшей теплоизоляции. Так как найти точное
количество конденсируемой на поверхности тестовых заготовок влаги не
представляется возможным, то точное количество испаряемой воды и зависящее
от него количество вносимой с паром энергии не поддается математическому
описанию. Следует учесть, что конденсация влаги на поверхности тестовых
заготовок происходит только в период. пока температура поверхности тестовых
заготовок не достигнет температуры точки росы для данных параметров среды,
то есть в первой половине операции расстойки. Далее конденсация
прекращается, и необходимость в испарении воды для увлажнения воздуха в
расстойном шкафу отпадает.
Также не является математически описуемым и коэффициент теплоемкости
влажного воздуха (свозд), зависящий от его температуры и влажности.
В связи с этим нами были приняты следующие упрощения и допущения:
. Коэффициенты теплоотдачи расчитываются по экспериментальным
критериальным уравнениям. Учитывая, что температура воздуха в
расстойном шкафу в установившемся режиме работы поддерживается
системой управления в установленных пределах относительно заданной
температуры (Тзад), то параметры воздуха для нахождения критериев
подобия берутся при неизменной температуре, равной заданной
температуре (Тзад) в камере расстойного шкафа.
. Коэффициент теплоемкости влажного воздуха расчитывается для заданных
значений его температуры и относительной влажности.
. Энергия, вносимая с паром, не учитывается. Это возможно благодаря
допущению о полном отсутствии конденсации в установившемся режиме
работы расстойного шкафа.
. Камера расстойного шкафа считается абсолютно герметичной.
. Давление воздуха в камере расстойного шкафа постоянное (p=const).
. Рассматривается нагрев и охлаждение термически тонких тел ( ( (( ( (
( ).
. Система поддержания влажности не рассматривается.
Уравнение теплового баланса расстойного шкафа
Уравнение теплового баланса расстойного шкафа:
Qвозд = Qтэн - Qтеста - Qтел - Qст ,
где Qвозд - теплота затрачиваемая на прогрев воздуха;
Qтэн - тепловой поток с поверхности ТЭНов;
Qтеста - количество теплоты, идущее на прогрев теста;
Qтел - количество теплоты, идущее на прогрев тележек;
Qст - потеря тепла через стенки.
Распишем все составляющие этого уравнения.
Теплота, затрачиваемая на прогрев воздуха
может быть описана как:
Qвозд = cвозд ( mвозд ( (dTвозд / dt),
откуда:
[pic],
где dTвозд/dt - скорость изменения температуры воздуха.
cвозд - теплоемкость воздуха:
cвозд = (св + cп ( dв/1000),
где св - теплоемкость сухого воздуха, при температуре
40(С : св = 1005
Дж/(кг(гр);
сп - теплоемкость перегретого пара:
сп = 2000 Дж/(кг(гр);
dп - влагосодержание воздуха, при температуре 40(С и
относительной влажности 75% оно равно:
dп = 36,9 г/кг;
Таким образом:
cвозд= (1005+2000(36,9/1000) =1079 Дж/(кг(гр); mвозд -
масса воздуха в расстойном шкафу;
mвозд = (возд ( Vвозд ,
где (возд - плотность влажного воздуха в камере
расстойного шкафа, при температуре 40(С и
относительной влажности 75%:
(возд = 1,11 кг/м3;
Vвозд - объем воздуха в камере расстойного шкафа:
Vвозд = 2 м3; (Лен)
Таким образом:
mвозд = 1,11 ( 2 = 2,22 кг; (Лен)
Тепловой поток с поверхности ТЭНов
описывается с помощью уравнения конвективной теплопередачи:
Qтэн = Ктэн ( (Ттэн - Твозд),
где Ттэн - температура ТЭНов;
Твозд - температура циркулирующего воздуха.
Ктэн - коэффициент, расчитываемый по формуле:
Ктэн = (тэн ( Sтэн ,
где Sтэн - площадь поверхности ТЭНов:
Sтэн = lтэн ( ( ( dтэн ,
где lтэн - длина ТЭНов;
dтэн - диаметр ТЭНов,
Откуда:
Sтэн = 2 ( ( ( 0,006 = 0,0377 м2;
(тэн - коэффициент теплоотдачи ТЭНов. Данный
коэффициент расчитывается по критериальному уравнению:
Nu = 0,238 ( Ref0,6 ,
где Ref - число Рейнольдса, вычисляемое:
Ref = ( ( dтэн / (,
где ( - скорость потока воздуха:
( = 5 м/c
dтэн - диаметр ТЭНов - их определяющий размер:
dтэн = 0,006 м;
( - коэффициент кинематической вязкости, для воздуха
при температуре 40(С :
( = 16,96 ( 10-6 м2/с.
Таким образом:
Ref = 5 ( 0,006 / 16,96(10-6 = 1769,
Следовательно:
Nu = 0,238 ( 17690,6 = 21,15 ,
Откуда:
(тэн = Nu ( ( / dтэн ,
где ( - коэффициент теплопроводности воздуха, при
температуре 40(С:
( = 2,76(10-2 Вт/(м(гр),
Значит:
(тэн = 21,15 ( 2,76(10-2 / 0,006 = 97 Вт/(м2 ( гр).
Таким образом:
Ктэн = 97 ( 0,0377 = 3,6568 Вт/гр
и
Qтэн = 3,6568 ( (Ттэн - Твозд).
При этом, излишки энергии ТЭНов идут на изменение их температуры:
[pic],
где dT/dt - скорость изменения температуры ТЭНов;
Ртэн - мощность ТЭНов.
Pтэн = 2000 Вт
Обоснование выбора такой мощности ТЭНов приведено в
разделе 6
Qтэн - тепловой поток с поверхности ТЭНов;
cтэн - теплоемкость материала ТЭНов, для ТЭНов изготовленных из
кантала А-1:
cтэн = 470 Дж/(кг(гр);
mтэн - масса ТЭНов:
mтэн = (тэн ( lтэн ( ( ( dтэн2 / 4 ,
где lтэн = 2 м - длина ТЭНов;
dтэн = 0,006 м - диаметр ТЭНов,
(тэн = 7100 кг/м3;
Откуда:
mтэн = 7100 ( 2 ( ( ( (0,006)2 / 4 = 0,4 кг,
В связи с тем, что в процессе расстойки необходимо поддерживать
заданную температуру, ТЭНы включены только пока температура воздуха в
камере расстойного шкафа меньше заданной. Как только температура воздуха
превышает заданный предел на величину допустимого отклонения, система
управления подает сигнал на отключение ТЭНов. При этом Ртэн = 0. При
падении температуры за нижний предел система управления подает сигнал на
включение ТЭНов. При этом Ртэн = Ртэн зад , где Ртэн зад - номинальная
мощность ТЭНов.
Тепловой поток, получаемый тестовыми заготовками
Тепловой поток, получаемый тестовыми заготовками и используемый для их
прогрева, может быть описан формулой конвективного теплообмена:
Qтеста = Ктеста ( (Твозд - Ттеста),
где Ктеста = (теста ( Sтеста ,
где Sтеста - площадь поверхности тестовых заготовок:
Sтеста = 2(10(0,45(0,66 = 6 м2;
(теста - коэффициент теплоотдачи поверхности тестовых
заготовок, расчитывается по экспериментальной критериальной
формуле:
Nu = 0,216 ( Re0,8,
где Re - число Рейнольдса, вычисляемое по
формуле:
Re = ( ( lтест / (,
где ( - скорость потока воздуха:
( = 0,4 м/c
lтест - определяющий размер тестовых заготовок:
lтест = 0,25 м;
( - коэффициент кинематической вязкости, для воздуха
при температуре 40(С :
( = 16,96 ( 10-6 м2/с.
Таким образом:
Re = 0,4 ( 0,25 / 16,96(10-6 = 5900,
Следовательно:
Nu = 0,216 ( 59000,8 = 224,46 ,
Откуда:
(теста = Nu ( ( / lтест ,
где ( - коэффициент теплопроводности воздуха, при
температуре 40(С:
( = 2,76(10-2 Вт/(м(гр),
Значит:
(теста = 224,46 ( 2,76(10-2 / 0,25 = 24,8 Вт/(м2 ( гр).
Таким образом:
Ктеста = 24,8 ( 6 = 148,8 Вт/гр
и
Qтеста = 148,8 ( (Твозд - Ттеста),
где Ттеста - температура тестовых заготовок, скорость изменения
которой, с учетом того, что при расстойке в тестовых
заготовках выделяется энергия Qтеста выд, составляет:
[pic],
где cтеста - теплоемкость тестовых заготовок:
cтеста = 3000 Дж/(кг(гр)
mтеста - масса тестовых заготовок:
mтеста = nтест заг ( mтест заг ,
где nтест заг =120 шт. - число тестовых заготовок;
mтест заг = 0,46 кг - масса тестовой заготовки;
Откуда: mтеста = 120 ( 0,46 = 55,2 кг,
Тепловой поток, получаемый тележками
Аналогично, тепловой поток, получаемый тележками и используемый для их
прогрева, также может быть описан формулой конвективного теплообмена:
Qтел = Ктел ( (Твозд - Ттел),
где Ктел = (тел ( Sтел ,
где Sтел - площадь поверхности тележек:
Sтел = 2(2(10(0,45(0,66 + 2(4(4(0,02(1,8 = 7 м2;
(тел - коэффициент теплоотдачи поверхности тележек,
расчитывается по экспериментальной критериальной формуле:
Nu = 0,064 ( Re0,8,
где Re - число Рейнольдса, вычисляемое по
формуле: Re = ( ( lтел / (,
где ( - скорость потока воздуха:
( = 0,4 м/c
lтел - определяющий размер тележек:
lтел = 0,66 м;
( - коэффициент кинематической вязкости, для воздуха
при температуре 40(С :
( = 16,96 ( 10-6 м2/с.
Таким образом:
Re = 0,4 ( 0,66 / 16,96(10-6 = 15566,
Следовательно:
Nu = 0,064 ( 155660,8 = 144,52 ,
Откуда:
(тел = Nu ( ( / lтел ,
где ( - коэффициент теплопроводности воздуха, при
температуре 40(С:
( = 2,76(10-2 Вт/(м(гр),
Значит:
(тел = 144,52 ( 2,76(10-2 / 0,66 = 6 Вт/(м2 ( гр).
Таким образом:
Ктел = 6 ( 7 = 42 Вт/гр
и Qтел = 42 ( (Твозд - Ттел),
где Ттел - температура тележек, скорость изменения которой:
[pic],
где cтел - теплоемкость тележек:
cтел = 500 Дж/(кг(гр);
mтел - масса тележек:
mтел = 50 кг.
Потери теплоты через стенки расстойного шкафа
рассчитываются по уравнению теплопередачи:
Qст = Кст ( (Твозд - Тос),
где Тос - температура окружающей среды.
Кст = kcт ( Sст ,
где Sст - площадь стенок камеры расстойного шкафа:
Sст = (1,85((1,4+0,7)+1,4(0,7)(2 = 9,73 м2;
kст - коэффициент теплопередачи через стенки:
[pic],
где (ст - толщина стали стенок расстойного шкафа:
(ст = 0,001 м;
(утепл - толщина утеплителя:
(утепл = 0,03 м;
(ст - коэффициент теплопроводности стальных стенок
расстойного шкафа:
(ст = 45 Вт/(м(гр);
(утепл - коэффициент теплопроводности утеплителя:
(утепл = 0,1 Вт/(м(гр);
(1 - общий коэффициент теплоотдачи к внутренней поверхности
стенок расстойного шкафа;
(2 - общий коэффициент теплоотдачи от наружной поверхности
стенок расстойного шкафа.
Общие коэффициенты теплоотдачи методически оцениваются
одинаково - как сумма коэффициентов теплоотдачи конвекцией ((кон) и
излучением ((изл),
(общ = (кон + (изл ,
где первая составляющая:
(кон = Nu ( ( / hст ,
где ( - коэффициент теплопроводности воздуха;
hст - определяющий размер стенок камеры расстойного
шкафа - их высота:
hст = 1,85 м;
Nu - коэффициент подобия Нуссельта:
Для омывания газами вертикальных поверхностей:
Nu = 0,15((Grвозд(Prвозд)1/3,
где Pr - число Прандтля характеризует собой свойства
среды;
Gr = ((g(hст3((t/(2 - число Грасгофа,
где g - ускорение свободного падения;
(t - температурный перепад между средой и
омываемой ею поверхностью;
( - функция, связывающая изменение плотности среды с
температурой. Для газов можно принять:
( = 1/T;
( - коэффициент кинематической вязкости среды.
вторая составляющая общего коэффициента теплоотдачи:
[pic],
где (ст - степень черноты стенок:
(ст = 0,9;
Тст - температура стенок, ( С;
(ст - постоянная Стефана-Больцмана:
(ст = 5,67 Вт/(м2(К4).
Исходя из того, что температура на внутренней и внешней поверхности
стенок расстойного шкафа является неизвестной величиной, принимаем в первом
приближении:
(1 = (2 = 10 Вт/(м2(гр);
Тогда коэффициент теплопередачи через стенки расстойного шкафа
составит:
[pic]
откуда
Qст = 2 ( ( 40 - 20 ) ( 9,73 = 390 Вт.
При этих данных температура на внутренней поверхности стенок камеры
расстойного шкафа составит
[pic],
аналогично, на наружной поверхности
[pic],
Во втором приближении:
Для внутренней поверхности стенок:
Pr1 = 0,699 (при T = 40 ( С)
Учитывая, что при T = 40(С
(возд = 16,96(10-6 м2/c ,
получим:[pic]
Тогда:
Nu1=0,15((Gr1(Pr1)1/3=0,15((2,7596(109(0,699)1/3=186,724
Откуда, учитывая, что при T = 40(С
(возд = 2,756(10-2 Вт/(м(гр),
получим
(кон1 = Nu1((возд/hст=186,724(2,76(10-2/1,85 = 2,79 Вт/(м2(гр)
Значение коэффициента теплоотдачи излучением:
[pic]
Следовательно, общий коэффициент теплоотдачи к внутренней поверхности
стенок расстойного шкафа составляет
(1 = (кон1 + (изл1 = 2,79 + 6,14 = 8,93 Вт/(м2(гр).
Аналогично, для внешней поверхности стенок расстойного шкафа:
Pr1 = 0,703 (при T = 20 ( С)
Учитывая, что при T = 20(С
(возд = 15,06(10-6 м2/c ,
получим:
[pic]Тогда:
Nu2=0,15((Gr2(Pr2)1/3=0,15((3,7388(109(0,703)1/3=207
Откуда, учитывая, что при T = 20(С
(возд = 2,59(10-2 Вт/(м(гр),
получим
(кон2 = Nu2((возд/hст=207(2,59(10-2/1,85 = 2,898 Вт/(м2(гр)
Значение коэффициента теплоотдачи излучением:
[pic]
Следовательно, общий коэффициент теплоотдачи к внутренней поверхности
стенок расстойного шкафа составляет
(1 = (кон1 + (изл1 = 2,898 + 5,24 = 8,14 Вт/(м2(гр).
Коэффициент теплопередачи через стенки расстойного шкафа во втором
приближении составит:
[pic]
откуда потери теплоты через стенки расстойного шкафа:
Qст = 1,87 ( ( 40 - 20 ) ( 9,73 = 363,8 Вт.
При этих данных температура на внутренней поверхности стенок камеры
расстойного шкафа составит
[pic],
аналогично, на наружной поверхности
[pic],
Степень расхождения между первым и вторым приближениями для каждой из этих
температур:
(т’ = 100 ( ( 36 - 35,8 )/ 36 = 0,6%;
(т’ = 100 ( ( 24,6 - 24 )/ 24 = 2,5%.
Это допустимо. В этой связи результаты второго приближения принимаем
за окончательные.
Для них выполним проверку на наличие или отсутствие конденсации пара
из парогазовой среды на внутренней поверхности стенок камеры расстойного
шкафа. Во избежание нежелательной конденсации пара необходимо, чтобы
температура на внутренней поверхности стенок Т’cт превышала температуру
точки росы Тр:
Т’cт ( Тр.
Для оптимальных (расчетных параметров) расстойки - температуры
парогазовой среды 40 (С и относительной влажности 75%, согласно данным
таблиц, температура точки росы
Тр = 34,5(С.
Отсюда следует, что в нашем случае конденсация пара на внутренней
поверхности стенок в установившемся режиме работы расстойного шкафа
отсутствует.
Окончательная формула потери теплоты через стенки расстойного шкафа, с
учетом того что
Кст = kcт ( Sст = 1,87 ( 9,73 = 18,2 Вт/гр,
запишется как
Qст = 18,2 ( (Твозд - Тос),
где Тос - температура окружающей среды.
Система дифференциальных уравнений
Таким образом, для моделирования работы системы управления расстойным
шкафом необходимо решить систему дифференциальных уравнений:
(T = Tзад - Tвозд;
[pic];
Qтэн = 3,6568 ( (Ттэн - Твозд);
dTтэн/dt = (2000 - Qтэн)/(470 ( 0,4);
Qтеста = 148,8 ( (Твозд - Ттеста);
dTтеста/dt = (Qтеста + 100)/( 3000 ( 120);
Qтел = 42 ( (Твозд - Ттел);
dTтел/dt = Qтел / (500 ( 50);
Qст = 18,2 ( (Твозд - Тос);
Qвозд = Qтэн - Qтеста - Qтел - Qст ;
dTвозд/dt = Qвозд /(1079(2,22).
Расчет и идентификация процессов протекающих в расстойном шкафу
Для расчета термодинамических процессов происходящих в камере
расстойного шкафа при расстойке тестовых заготовок, а также для выбора
параметров СУ обеспечивающих заданный режим, была разработана программа для
ЭВМ, моделирующая работу системы управления расстойным шкафом. Блок-схема
данной программы приведена на чертеже, а текст программы приведен в
Приложении 1. По результатом работы программы были построены переходный
процесс и фазовый портрет (см. графики). При этом мощность ТЭНов и допуск
на отклонение температуры воздуха в камере расстойного шкафа от заданного
значения были выбраны исходя из результатов исследований, изложенных в
разделе 6. Из графика переходного процесса видно, что, после выхода в
установившийся режим, температура циркулирующего в камере расстойного шкафа
воздуха поддерживается на заданном уровне, не выходя за пределы заданного
допуска, а температура поверхности тестовых заготовок достигает заданной к
окончанию времени расстойки. Это говорит о правильности расчетов и верности
выбора параметров СУ.
Также была проведена идентификация разработанной модели СУ расстойного
шкафа с работающим образцом. Отклонения параметров работы модели от образца
оказались небольшими, что указывает на правильный выбор допущений и
упрощений, сделанных в процессе разработки данной модели. |
