Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
Предоплата всего
Подписываем
В формулах конвекционной теплопередачи присутствуют коэффициенты теплоотдачи . Как было показано ранее, коэффициенты теплоотдачи зависят от многих факторов: от температур поверхности и омывающей ее среды, скорости движения последней, ее теплопроводности, вязкости, плотности и теплоемкости, от конфигурации и состояния поверхности и омывающей ее среды. В связи с невозможностью математического описания данных коэффициентов, для их нахождения пользуются экспериментальными данными, широко используя теорию подобия, позволяющую в известной степени обобщить полученные опытные результаты. Но используемые для нахождения коэффициентов теплоотдачи критериальные уравнения содержат критерии подобия (Nu, Pe, Re, Pr, Gr), которые зависят от многих параметров поверхностей и омывающей их среды, некоторые из которых зависят от температуры среды и от разности между ней и температур омываемых ею поверхностей. Данные зависимости не описаны математически. Конденсация влаги на поверхности тестовых заготовок в процессе их расстойки еще больше затрудняет точное нахождение коэффициента теплоотдачи их поверхности.
Конденсация влаги на поверхности тестовых заготовок, а также на внутренней поверхности стенок камеры расстойного шкафа приводит к уменьшению абсолютной и относительной влажности в камере расстойного шкафа. Для поддержания заданной относительной влажности воздуха применяется испарение воды, контролируемое проектируемой системой управления. Но вместе с паром в камеру расстойного шкафа попадает дополнительная энергия. Конденсацию влаги на внутренней поверхности стенок камеры расстойного шкафа можно свести к минимуму путем их лучшей теплоизоляции. Так как найти точное количество конденсируемой на поверхности тестовых заготовок влаги не представляется возможным, то точное количество испаряемой воды и зависящее от него количество вносимой с паром энергии не поддается математическому описанию. Следует учесть, что конденсация влаги на поверхности тестовых заготовок происходит только в период. пока температура поверхности тестовых заготовок не достигнет температуры точки росы для данных параметров среды, то есть в первой половине операции расстойки. Далее конденсация прекращается, и необходимость в испарении воды для увлажнения воздуха в расстойном шкафу отпадает.
Также не является математически описуемым и коэффициент теплоемкости влажного воздуха (свозд), зависящий от его температуры и влажности.
В связи с этим нами были приняты следующие упрощения и допущения:
Коэффициенты теплоотдачи расчитываются по экспериментальным критериальным уравнениям. Учитывая, что температура воздуха в расстойном шкафу в установившемся режиме работы поддерживается системой управления в установленных пределах относительно заданной температуры (Тзад), то параметры воздуха для нахождения критериев подобия берутся при неизменной температуре, равной заданной температуре (Тзад) в камере расстойного шкафа.
Коэффициент теплоемкости влажного воздуха расчитывается для заданных значений его температуры и относительной влажности.
Энергия, вносимая с паром, не учитывается. Это возможно благодаря допущению о полном отсутствии конденсации в установившемся режиме работы расстойного шкафа.
Камера расстойного шкафа считается абсолютно герметичной.
Давление воздуха в камере расстойного шкафа постоянное (p=const).
Рассматривается нагрев и охлаждение термически тонких тел ( ).
Система поддержания влажности не рассматривается.
Уравнение теплового баланса расстойного шкафа:
Qвозд = Qтэн - Qтеста - Qтел - Qст ,
где Qвозд - теплота затрачиваемая на прогрев воздуха;
Qтэн - тепловой поток с поверхности ТЭНов;
Qтеста - количество теплоты, идущее на прогрев теста;
Qтел - количество теплоты, идущее на прогрев тележек;
Qст - потеря тепла через стенки.
Распишем все составляющие этого уравнения.
может быть описана как:
Qвозд = cвозд mвозд (dTвозд / dt),
откуда:
,
где dTвозд/dt - скорость изменения температуры воздуха.
cвозд - теплоемкость воздуха:
cвозд = (св + cп dв/1000),
где св - теплоемкость сухого воздуха, при температуре 40С : св = 1005 Дж/(кггр);
сп - теплоемкость перегретого пара:
сп = 2000 Дж/(кггр);
dп - влагосодержание воздуха, при температуре 40С и относительной влажности 75% оно равно:
dп = 36,9 г/кг;
Таким образом:
cвозд= (1005+200036,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,9610-6 = 1769,
Следовательно:
Nu = 0,238 17690,6 = 21,15 ,
Откуда:
тэн = Nu / dтэн ,
где - коэффициент теплопроводности воздуха, при температуре 40С:
= 2,7610-2 Вт/(мгр),
Значит:
тэн = 21,15 2,7610-2 / 0,006 = 97 Вт/(м2 гр).
Таким образом:
Ктэн = 97 0,0377 = 3,6568 Вт/гр
и
Qтэн = 3,6568 (Ттэн - Твозд).
При этом, излишки энергии ТЭНов идут на изменение их температуры: ,
где 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теста = 2100,450,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,9610-6 = 5900,
Следовательно:
Nu = 0,216 59000,8 = 224,46 ,
Откуда:
теста = Nu / lтест ,
где - коэффициент теплопроводности воздуха, при температуре 40С:
= 2,7610-2 Вт/(мгр),
Значит:
теста = 224,46 2,7610-2 / 0,25 = 24,8 Вт/(м2 гр).
Таким образом:
Ктеста = 24,8 6 = 148,8 Вт/гр
и
Qтеста = 148,8 (Твозд - Ттеста),
где Ттеста - температура тестовых заготовок, скорость изменения которой, с учетом того, что при расстойке в тестовых заготовках выделяется энергия Qтеста выд, составляет:
,
где cтеста - теплоемкость тестовых заготовок:
cтеста = 3000 Дж/(кггр)
mтеста - масса тестовых заготовок:
mтеста = nтест заг mтест заг ,
где nтест заг =120 шт. - число тестовых заготовок;
mтест заг = 0,46 кг - масса тестовой заготовки;
Откуда: mтеста = 120 0,46 = 55,2 кг,
Аналогично, тепловой поток, получаемый тележками и используемый для их прогрева, также может быть описан формулой конвективного теплообмена:
Qтел = Ктел (Твозд - Ттел),
где Ктел = тел Sтел ,
где Sтел - площадь поверхности тележек:
Sтел = 22100,450,66 + 2440,021,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,9610-6 = 15566,
Следовательно:
Nu = 0,064 155660,8 = 144,52 ,
Откуда:
тел = Nu / lтел ,
где - коэффициент теплопроводности воздуха, при температуре 40С:
= 2,7610-2 Вт/(мгр),
Значит:
тел = 144,52 2,7610-2 / 0,66 = 6 Вт/(м2 гр).
Таким образом:
Ктел = 6 7 = 42 Вт/гр
и Qтел = 42 (Твозд - Ттел),
где Ттел - температура тележек, скорость изменения которой:
,
где cтел - теплоемкость тележек:
cтел = 500 Дж/(кггр);
mтел - масса тележек:
mтел = 50 кг.
рассчитываются по уравнению теплопередачи:
Qст = Кст (Твозд - Тос),
где Тос - температура окружающей среды.
Кст = kcт Sст ,
где Sст - площадь стенок камеры расстойного шкафа:
Sст = (1,85(1,4+0,7)+1,40,7)2 = 9,73 м2;
kст - коэффициент теплопередачи через стенки:
,
где ст - толщина стали стенок расстойного шкафа:
ст = 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 = ghст3t/2 - число Грасгофа,
где g - ускорение свободного падения;
t - температурный перепад между средой и омываемой ею поверхностью;
- функция, связывающая изменение плотности среды с температурой. Для газов можно принять:
= 1/T;
- коэффициент кинематической вязкости среды.
вторая составляющая общего коэффициента теплоотдачи:
,
где ст - степень черноты стенок:
ст = 0,9;
Тст - температура стенок, С;
ст - постоянная Стефана-Больцмана:
ст = 5,67 Вт/(м2К4).
Исходя из того, что температура на внутренней и внешней поверхности стенок расстойного шкафа является неизвестной величиной, принимаем в первом приближении:
1 = 2 = 10 Вт/(м2гр);
Тогда коэффициент теплопередачи через стенки расстойного шкафа составит:
откуда
Qст = 2 ( 40 - 20 ) 9,73 = 390 Вт.
При этих данных температура на внутренней поверхности стенок камеры расстойного шкафа составит
,
аналогично, на наружной поверхности
,
Во втором приближении:
Для внутренней поверхности стенок:
Pr1 = 0,699 (при T = 40 С)
Учитывая, что при T = 40С
возд = 16,9610-6 м2/c ,
получим:
Тогда:
Nu1=0,15(Gr1Pr1)1/3=0,15(2,75961090,699)1/3=186,724
Откуда, учитывая, что при T = 40С
возд = 2,75610-2 Вт/(мгр),
получим
кон1 = Nu1возд/hст=186,7242,7610-2/1,85 = 2,79 Вт/(м2гр)
Значение коэффициента теплоотдачи излучением:
Следовательно, общий коэффициент теплоотдачи к внутренней поверхности стенок расстойного шкафа составляет
1 = кон1 + изл1 = 2,79 + 6,14 = 8,93 Вт/(м2гр).
Аналогично, для внешней поверхности стенок расстойного шкафа:
Pr1 = 0,703 (при T = 20 С)
Учитывая, что при T = 20С
возд = 15,0610-6 м2/c ,
получим:
Тогда:
Nu2=0,15(Gr2Pr2)1/3=0,15(3,73881090,703)1/3=207
Откуда, учитывая, что при T = 20С
возд = 2,5910-2 Вт/(мгр),
получим
кон2 = Nu2возд/hст=2072,5910-2/1,85 = 2,898 Вт/(м2гр)
Значение коэффициента теплоотдачи излучением:
Следовательно, общий коэффициент теплоотдачи к внутренней поверхности стенок расстойного шкафа составляет
1 = кон1 + изл1 = 2,898 + 5,24 = 8,14 Вт/(м2гр).
Коэффициент теплопередачи через стенки расстойного шкафа во втором приближении составит:
откуда потери теплоты через стенки расстойного шкафа:
Qст = 1,87 ( 40 - 20 ) 9,73 = 363,8 Вт.
При этих данных температура на внутренней поверхности стенок камеры расстойного шкафа составит
,
аналогично, на наружной поверхности
,
Степень расхождения между первым и вторым приближениями для каждой из этих температур:
т’ = 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возд;
;
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возд /(10792,22).
Для расчета термодинамических процессов происходящих в камере расстойного шкафа при расстойке тестовых заготовок, а также для выбора параметров СУ обеспечивающих заданный режим, была разработана программа для ЭВМ, моделирующая работу системы управления расстойным шкафом. Блок-схема данной программы приведена на чертеже, а текст программы приведен в Приложении 1. По результатом работы программы были построены переходный процесс и фазовый портрет (см. графики). При этом мощность ТЭНов и допуск на отклонение температуры воздуха в камере расстойного шкафа от заданного значения были выбраны исходя из результатов исследований, изложенных в разделе 6. Из графика переходного процесса видно, что, после выхода в установившийся режим, температура циркулирующего в камере расстойного шкафа воздуха поддерживается на заданном уровне, не выходя за пределы заданного допуска, а температура поверхности тестовых заготовок достигает заданной к окончанию времени расстойки. Это говорит о правильности расчетов и верности выбора параметров СУ.
Также была проведена идентификация разработанной модели СУ расстойного шкафа с работающим образцом. Отклонения параметров работы модели от образца оказались небольшими, что указывает на правильный выбор допущений и упрощений, сделанных в процессе разработки данной модели.