Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
Предоплата всего
Подписываем
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ
им. Н.Э.БАУМАНА
Факультет МТ
Кафедра МТ8
КУРСОВОЙ ПРОЕКТ
РАСЧЁТ И КОНСТРУИРОВАНИЕ ОБОРУДОВАНИЯ ДЛЯ ТЕПЛОВОЙ ОБРАБОТКИ
ТЕМА:
Лабораторная высокотемпературная вакуумная печь сопротивления
Студент: ______________(Каримова К.Л.) Группа_МТ8-82_
Руководитель проекта: ___________________(Горничев А.А.)
2008 г.
Введение
1.2.1 Определение S
1.2.2 Определение α
1.3 Определение критерия Био
4.Определение производительности печи
6.4 Расчёт потерь через теплоизоляционный слой
6.5 Расчёт Qкл
6.11 Определение КПД печи
6.12 Определение удельного расхода энергии
8. Определение удельного расхода энергии
10.Расчет механизма подъёма
11.Определение ориентировочной стоимости нагревательной камеры
12. Составление таблицы ТЭП печи
13. Список литературы
Техническое задание
В данном курсовом проекте требуется спроектировать лабораторную высокотемпературную вакуумную печь сопротивления для отжига.
Материал садки – ниобий, Тпл=2468°С.
Максимальная масса садки 5 кг.
Коэффициент заполнения садки К=0,5.
Максимальная температура печи 1250 °С.
Давление в нагревательной камере Рраб=2*10-5Па.
Вакуум создаётся при помощи диффузионного насоса.
Огнеупорный материал - шамот ультралегковес.
Теплоизолирующий материал – вакуум и экраны.
Нагреватель - свернутая в спираль проволока из Х20Н80, защищённая от возможного окисления кварцевым стеклом.
Форма камеры – цилиндрическая.
На боковой поверхности камеры расположен глазок из кварцевого стекла диаметром 2см.
Для загрузки и выгрузки садки предусмотрено отверстие, которое закрывается крышкой.
Крышка открывается с помощью механизма подъёма.
Для контроля рабочей температуры используются термопары.
Управление технологическими параметрами ведется при помощи цифровой системы контроля посредством клавиатуры и дисплея.
1.Определение теплотехнических характеристик и тепловой массивности садки
- коэффициент теплопроводности λме , для реальной садки λс = λмеК,
- коэффициент теплоёмкости сме, для реальной садки сс= сме,
- плотность ρме, для реальной садки ρс = ρмеК,
- коэффициент температуропроводности аме= λме /смеρме, для реальной ас= аме,
где К - коэффициент заполнения садки
Теплотехнические характеристики определяем по справочным материалам.
Определение исходных данных для заданной температуры нагрева произведено методом линейной интерполяции:
Температуры печи, начала и конца нагрева садки:
Тп=1250+273=1523 К,
T1=20+273=293 К,
T2=1200+273=1473 К,
1=52 Вт/мК, с1=268 Дж/кг*К
2=76 Вт/мК, с2=324 Дж/кг*К
Средняя температура периода нагрева:
, Тср=610+273=883К.
Òåïëîïðîâîäíîñòü ìàòåðèàëà ïðè ñðåäíåé òåìïåðàòóðå:
;
Òåïëîåìêîñòü ìàòåðèàëà ïðè ñðåäíåé òåìïåðàòóðå:
ρме= 8570 кг/м3.
Пересчёт параметров с учётом коэффициента заполнения садки:
λс= 640,5 = 32 (Вт/мК)
ρс= 85700,5= 4285 (кг/м3)
сс=сср=296 (Дж/кгК)
Определение коэффициента температуропроводности ас:
ас=32/2964285 =2510-6 м2/сК.
1.2 Определение тепловой массивности садки
При определении ТМС используем критерий Био:
Bi = αS /λс
Для этого требуется рассчитать коэффициент теплоотдачи α и уточнить величину S (наименьшее расстояние между самой горячей и самой холодной точками).
1.2.1 Определение характеристического размера садки
Отжигу подвергается садка цилиндрической формы, у которой диаметр равен высоте. Зная максимальную массу садки и плотность материала, можем рассчитать габариты детали.
m=ρ*π*D3/4, => D=H=3√(4m/πρ)=0,09м.
Из анализа габаритов садки следует, что характеристический размер будет равен половине толщины садки S=0,045м.
1.2.2Определение коэффициента теплоотдачи
В данном проекте требуется разработать вакуумную печь, то есть при расчёте α будем учитывать только составляющую излучения (αизл).
Коэффициент теплоотдачи αизл определяем по закону Стефана – Больцмана:
где спр -приведенный коэффициент (излучательная способность материала)
спр=εпрс0
с0=5,67 Вт/м2К4 - излучательная способность АЧТ
ε1– степень черноты материала садки (ε=0,8)
ε2– степень черноты материала печи (ε=0,9)
- тепловоспринимающая поверхность садки; - печи; принимаем = 0,5,
εпр=0,77
спр=0,775,67=4,37 Вт/м2*К4
Находим αизл для температурначала и конца нагрева садки:
α1=
α2=
αср.= (α1 + α2) / 2 =(191 +588)/2=390 (Вт/м2*К)
1.3 Определение критерия Био:
Bi =3900,045/32=0,548 >0,5 => садка массивная.
2. Расчёт времени нагрева садки
Рис.1 Диаграмма нагрева ТМТ
На этапе нагрева печи при постоянной мощности (τ1)рекомендуется принять:
Tп =0,85Тз=0,85 1523 =1295K
T пов = 0,5Тз=0,51523=761,5K
Тепловой поток q= αср*( Tп - T пов )=390*(1295-761,5)=208065 Вт/м2
Температура, до которой нагревается поверхность садки за время :
, t’’пов=990-273=717°С,
τ '=0,3S2/ ас =0,30,0452/ 2510-6= 24,3 с
Для ∞ цилиндра Δt= qR/2λс = 1082080,045/232 = 76 К,
t’пов=1,3*Δt=1,3*76=98,8°С,
t’ц=0,5*Δt=0,5*76=38°С;
Т’’пов> Δt => полученные результаты можно использовать.
Относительная температура нагрева
=q/ сс*ρс*S=208065/296*4285*0,045=3,645,
τ ''=( t’’пов- t’пов)/ =(717-98,8)/3,645=169,6 с,
τ 1= τ '+ τ ''=24,3+169,6=193,9с=3,23 мин.
=
– температурный критерий
tпеч – температура печи при установившемся нагреве (tпеч=1250 °С)
tнач – начальная температура центра пластины (tнач= t’’пов=717°С)
tц – требуемая температура центра садки (tц=1200°С)
=(1250-1200)/(1250-717) = 0,094
Определяем критерий Био для нужного интервала температур:
αср.= (α717 + α1200) / 2 =(362 +588)/2=475 (Вт/м2К)
Bi = αср S/λс= 4750,045/32 =0,67
Рис.2 График Будрина.
Из графика получаем значение критерия Фурье F0=4,5 ,
с=6,075мин.
Время нагрева садки:
с=9,3мин=0,155ч.
3.Определение продолжительности цикла работы печи
τцикла = τнагрева + τвыдержки + τохлаждения + τвспомогательное
τнагрева=0,155 ч,
τвыдержки=2 ч,
τохлаждения=(1200-20)/7=169мин=2,81 ч,
τвспомогательное=0,25 ч.
τцикла=0,155+2+2,81+0,255,2 ч
4. Определение производительности печи
Зная массу садки и общее время цикла, можем найти производительность печи:
g =mc/ τцикла = 5/5,20,96 кг/ч.
5. Определение основных размеров нагревательной камеры
Учитывая коэффициент заполнения К=0,5 и габариты садки назначаем размеры рабочего пространства:
диаметр рабочего пространства D=dc+2*15=90+30=120мм,
высота рабочего пространства H=hc+20=90+20=110мм.
6. Расчёт теплового баланса
Уравнение теплового баланса:
Qрасх=Qпол+Qвсп+Qкл+Qакк+Qотв+Qктз+ Qнеуч
Qрасх –всё тепло, потребляемое печью.
Qпол - тепло, идущее на нагрев садки.
Qвсп- потери тепла, идущие на нагрев тары и атмосферы печи.
Qкл - тепло, пропускаемое кладкой.
Qакк- тепло, аккумулируемое всеми составляющими печи.
Qотв- потери тепла через отверстия.
Qткз- тепло, отводимое через источники теплового короткого замыкания.
Qнеуч – неучтённые потери.
Qподв- энергия, потребляемая из сети.
6.1 Рассчитываем потери тепла на нагрев садки Qпол
Qпол =mсc(tкон-tн)/ τнагрева= 5296(1200-20)/0,155=1,03107Дж/ч=2864,5 Вт
6.2 Потери на нагрев тары и атмосферы печи Qвсп
Qвсп= Qтары + Qатм=0,1 Qпол+vг cг (tкон-tнач г).
В нагревательной камере вакуум, поэтому вторая составляющая очень мала;
тара отсутствует.
Qвсп=0 Вт.
6.3.1 Расчет теплового потока через боковые стенки
Назначим минимальную ширину огнеупорного слоя =0,027м, по условию материал этого слоя – шамот-ультралегковес. Стенки печи цилиндрические.
Определяем площади:
Из-за кварцевого стекла, защищающего нагреватели, Dвн=0,12+20,009=0,138 м,
Dнар= Dвн+2=0,138+20,027=0,192 м
Fвн=Dвн hвн =0,1380,11=0,048 м2
Fнар= Dнар( hвн+2)=0,192 (0,11+20,027)=0,099 м2
Площадь смотрового отверстия Fотв = D2отв/4=0,022/4=0,0003 м2 мала, поэтому её не учитываем.
Fнар / Fвн =0,099 /0,048=2,06<2, поэтому Fср=√(FвнFнар )=0,069 м2.
Назначаем температуры:
на внутренней поверхности шамота t1=1200°C,
на внешней поверхности шамота t2=1020°C,
tср=1110°C.
Определяем коэффициент теплоотдачи шамота-ультралегковеса :
λшул=1,363(0,08+0,000141110)=0,321 .
αвн.= α1200 =588 (Вт/м2*К),
αнар= α1020=.
Расчет термических сопротивлений:
R1=1/ αвн Fвн= 1/5880,048=0,035 К/Вт
R2=/ λшул Fср=0,027/0,3210,069=1,216 К/Вт
R3=1/ αнар Fнар=1/4640,099=0,022 К/Вт
Тепловой поток:
Qклбок=( t1- t2)/ΣR=(1200-1020)/1,273=141,4 Вт.
Поскольку при выборе граничных температур очень велика вероятность ошибки, необходима проверка правильности назначенных температур.
Проверка слева
t1’= t1- Qклбок R1=1200-141,40,035=1194,99°C
t2’= t1- Qклбок (R1+ R2)=1200-141,4(0,035+1,216)=1023,08°C
Проверка справа
t1’’= t2+ Qклбок (R3+ R2)=1020+141,4(0, 022+1,216)= 1194,99°C
t2’’= t2+ Qклбок R3=1020+141,40, 022=1023,08°C
Разница между минимальным и максимальным значениями должна быть:
Δ t1=5,01°C, Δ t2=3,08°C, => назначенные температуры подходят.
Назначим ширину огнеупорного слоя =0,027м, материал слоя – шамот-ультралегковес. Форма крышки круглая.
Определяем площади:
Fвн=D2вн/4=0,1382/4 =0,015 м2
Fнар= (Dвн+2)2/4= (0,138+20,027) 2/4=0,029м2
Fнар / Fвн =0,029 /0,015=1,93<2, поэтому Fср=(Fвн+Fнар )/2=0,022м2.
Назначаем температуры:
на внутренней поверхности шамота t1=1200°C,
на внешней поверхности шамота t2=1020°C,
tср=1110°C.
Определяем коэффициент теплоотдачи шамота-ультралегковеса :
λшул=1,363(0,08+0,000141110)=0,321 .
αвн.= α1200 =588 (Вт/м2*К), αнар= α1020=464 (Вт/м2*К).
Расчет термических сопротивлений:
R1=1/ αвн Fвн= 1/5880,015=0,113 К/Вт
R2=/ λшул Fср=0,027/0,3210,022=3,824 К/Вт
R3=1/ αнар Fнар=1/4640,029=0,074 К/Вт
Тепловой поток:
Qклсвод=( t1- t2)/ΣR=(1200-1020)/4,012=44,9 Вт.
Проверка слева
t1’= t1- Qклсвод R1=1200-83,70,055=1194,97°C
t2’= t1- Qклсвод (R1+ R2)=1200-83,7(0,055+2,054)=1023,33°C
Проверка справа
t1’’= t2+ Qклсвод (R3+ R2)=1020+83,7(0,042+2,054)= 1194,97°C
t2’’= t2+ Qклсвод R3=1020+83,70,042 =1023,33°C
Разница между минимальным и максимальным значениями должна быть:
Δ t1=5,03°C, Δ t2=3,33°C, => назначенные температуры подходят.
QшΣ= Qклбок+ Qклсвод+ Qклпод=141,4+44,9+44,9=231,2 Вт.
6.4 Расчёт потери тепла через теплоизоляционный слой
6.4.1 Расчёт теплового потока через боковые стенки
Òîëùèíà ýêðàíà :
Òîëùèíà êåðàìè÷åñêèõ èçîëÿòîðîâ:
Площадь смотрового отверстия:
Габариты и площадь боковой наружной поверхности шамота:
Ïëîùàäè ïîâåðõíîñòей ýêðàíîâ:
Ïðèâåäåííûé êîýôôèöèåíò ëó÷åèñïóñêàíèÿ
Потери тепла через боковые экраны:
6.4.2 Расчёт теплового потока через крышку
Ïëîùàäü ïîâåðõíîñòè шамота:
Ïëîùàäè ïîâåðõíîñòè ýêðàíîâ:
Ïðèâåäåííûé êîýôôèöèåíò ëó÷åèñïóñêàíèÿ
Ïîòåðè òåïëà ÷åðåç верхние экраны:
6.4.3 Расчёт теплового потока через под
6.4.4 Суммарные потери тепла через экранную теплоизоляцию
6.5 Расчёт Qкл
Qкл= QшΣ+ QэΣ =321,2+3773,93=4095,13 Вт.
6.6 Расчёт Qакк
Рассчитываем тепло, затрачиваемое на нагрев теплоизоляции Qтеплоиз.
ρш-у=400 кг/м3 - плотность шамота-ультралегковеса.
Объём шамота Vш=δ(Fсрб+2Fсркр)=0,027(0,069+20,022)=3,05110-3 м3.
Масса шамота Мш= ρш-у Vш=4003,05110-3=1,22 кг.
ρэ=1800 кг/м3 - плотность материала экранной теплоизоляции, его объём
масса Мэ= ρэ Vэ= 18000,5810-3= 1,044 кг.
C := 0.21+ 0.000055⋅(tвн − tнаруж)/2 – теплоёмкость слоя из шамота-ультралегковеса.
с=0,24 Дж/кгК, Qшул=сМш(tвн-tнар)=0,241,22 (1200-1020)=48,3 Дж.
Удельная теплоёмкость экранного материала – углеграфиового композита
Qэ=сМэ(tвн-tнар)=1674,51,044 (1020-60)=1,68106Дж=466,2 Вт
Qтеплоиз= Qшул+ Qэ=48,3+466,2=514,5 Вт.
Рассчитываем тепло, расходуемое на нагрев кожуха Qкожух.
Кожух изготовлен из стали 12Х18Н10Т.
ρ=7145 кг/м3 - плотность кожуха
δк=8·10-3 м – толщина кожуха
dк=dвн+2 δк=0,24+2·8·10-3 =0,256 м2
Объём кожуха
V=π· dк·δк·hк+2·π·dк2/4·δк = π· dк·δк(hк+ dк/2)=π·0,256·0,008·(0,2+0.256/2)=2,11·10-3 м3.
Масса кожуха М=V·ρ=2,11·10-3·7145=15,08 кг.
Теплоемкость кожуха ск=483 Дж/кг°С.
Qкожух.=М·ск·( tвн-tнар)=15,08·483·(60-20)=291346 Дж=80,93 Вт.
Рассчитываем Qакк
Qакк= Qтеплоиз+ Qкожух=514,5+80,93=595,43 Вт
6.7 Расчёт Qткз
Средний диаметр тепловых замыканий
d := 0.03
F =πd2/4
λ := 0.64 -для муллито-коррундовой трубки (90%Al2O3)
Длины трубочек:
L1 := 0.027
L2 := 0.039
Температуры слоев, где происходит ТКЗ:
tокр1 := 105°С
tокр2 := 300°С
Потери на ТКЗ:
Qткз1отв =Fλ(tвн − tокр1)/L1
Qткз1отв = 19.1847
Qткз2отв =Fλ(tвн − tокр2)/L2
Qткз2отв = 11.0197
Qткз := Qткз1отв + 2Qткз2отв
Qткз = 41.22Вт
6.8 Расчёт потерь тепла через глазок Qотв
Ñòåïåíü ÷åðíîòû èçëó÷àþùåãî òåëà (ñòåêëà)
Êîýôôèöèåíò äèàôðàãìèðîâàíèÿ
Ïëîùàäü ñìîòðîâого îêîíа
Òåìïïåðàòóðà ïå÷è
6.9 Расчёт тепла на неучтенные потери Qнеуч
Для печей периодического действия
Qнеуч=0,10…0,15(Qпол +Qкл +Qакк +Qткз+Qотв)
Qнеуч =0,12(2864,5+4095,13+595,43+41,22+1,65)=0,127597,9=911,75 Вт
6.10 Окончательный расчёт теплового баланса
Qрасх=Qпол+Qвсп +Qкл+Qакк+ Qткз+ Qотв+Qнеуч
Qрасх=2864,5+0+4095,13+595,43+41,22+1,65+911,75= 8509,65 Вт
Уравнение теплового баланса
Qподв=1,016Qрасх=1,0168509,65= 8645,8 Вт
6.11 Определение КПД печи
η=Qпол/Qрасх
η= 2864,5/8645,8=0,33
7. Определение мощности печи
Руст=к1(Qпол+Qвсп)+к2(Qрасх-Qпол-Qвсп)
к1- учитывает колебания напряжения в сети (1,1-1,3), к1=1,2
к2- учитывает потери на «старение» печи (1,2-1,4), к2=1,3
Pхх=Pуст-Qпол=10776-2864,5=7911,5 Вт
8. Определение удельного расхода энергии
А= Руст (τнагрева+τвыд)/mс=107762,155/5=4644 Втч/кг=4,6 кВтч/кг
9. Расчёт и размещение нагревателей
Нагреватель спиральный проволочный, его проще всего разместить в печи цилиндрической формы. Расположен в пазах футеровки.
Руст<25 кВт, =>питание однофазное.
Площадь поверхности садки в зоне, воспринимающей излучение,
d := 0.09 м
h := 0.09 м
Fизд=πd2/4+πhd , Fизд = 0.0318 м2
Tнmax := 1390°С
Tн := 1350+ 273=1623 К
Tc := 1240+ 273=1513 К
εнагр := 0.8 -для нихромовых нагревателей
εизд := 0.8 -для ниобия
γ := 1.24 -удельное сопротивление, Оммм2/м
Площадь поверхности, занятой нагревателем:
Fст= πDН=π0,1380,11 =0,0477 м2
Спр
Спр=4 Вт/м2К4
Wид = 7 Вт/см2 (по рис.4 из методических указаний)
Определение W:
α эф := 0.32 – коэффициент эффективности спиралевидной системы нагревателей,
α г := 1,8 – коэффициент шага нагревателя (t/d=4,5)
(по рис.6 из методических указаний),
α с := 1,2 – коэффициент, зависящий от приведённого коэффициента лучеиспускания (по рис.9 из методических указаний),
α р := 0,92 – коэффициент, зависящий от площадей, испускающих и принимающих тепло, (по рис.10 из методических указаний), Fизд/Fст=0,67
W := Wидα эфαгαсαр, W = 4,5 Вт/см2=4,5*104 Вт/м2
d=
d = 0,216 мм
Принимаем в соответствии со стандартом:
d = 0,8 мм.
Сопротивление нагревателя, Ом:
R=U2/Pуст, R = 13,4
Длина нагревателя, м:
l=π d2R/4γ, l =5,43
Проверка по удельной мощности
W=Pуст/ πdl=0,79 Вт/см2
Проверка по температуре:
Pуст= Pнагр
Площадь поверхности одного нагревателя, м2:
Fнагр =πdl/1000, Fнагр = 0.0136
Tнагр = 1514 К, Tнагр<Tн , следовательно, условие выполняется.
Назначаем шаг витка нагревателя: t=4,5*d, t=4мм.
L0=lt/πD=1,73м.
Расчет выводов нагревателей.
По правилам техники безопасности длина части нагревателя, входящего во
внутреннее пространство печи, должна быть не менее 50 мм. Аналогично для
вывода в окружающее пространство. С учетом этого: Lвыв=lt/πD=0,03м.
L=L0- Lвыв=1,7м, n=L/πD=1,7/π0,146=3,7.
10. Расчёт механизма подъёма
Масса крышки
Мкр=ρшVш+ρэVэ+ρкVк=4005,9410-4+18001,110-4+71454,1210-4=3,38 кг.
Смета стоимости основных материалов, готовых узлов и деталей
|
№ |
Наименование |
Количество |
Цена, |
Сумма, руб |
|
1 |
Шамот-ультралегковес |
1,5 кг |
200 руб/кг |
300 |
|
2 |
Углеграфитовый композит |
1 кг |
14000 руб/кг |
14000 |
|
3 |
Керамические изделия |
0,5 кг |
100 руб/кг |
50 |
|
4 |
Сталь 12Х18Н10Т |
15 кг |
150 руб/кг |
2250 |
|
5 |
Резина вакуумная |
0,3 кг |
333 руб/кг |
100 |
|
6 |
Термопара (кабель) |
2 м |
700 руб/м |
1400 |
|
7 |
Нагреватель Х20Н80 |
6м |
100 руб/м |
600 |
|
8 |
Кварцевое стекло |
1шт. |
100 руб/шт |
100 |
|
9 |
Электроника |
1 комплект |
3000 руб/компл. |
3000 |
|
10 |
Насос |
1шт. |
2000 руб/шт |
2000 |
|
Всего: |
23700 |
Смета расходов на основную зарплату основным рабочим
|
№ по пор |
Специальность рабочего |
Разряд |
Объем работы, ч |
Расценка, руб/час |
Зарплата, руб |
|
1 |
Сварщик |
5 |
8 |
200 |
1600 |
|
2 |
Слесари по приспособлениям |
5 |
6 |
200 |
1200 |
|
3 |
Слесарь-сборщик |
6 |
18 |
200 |
3600 |
|
4 |
Вакуумщик |
5 |
18 |
200 |
3600 |
|
5 |
Инженер по наладке |
4 |
6 |
200 |
1200 |
|
Всего: |
11200 |
Смета заводской себестоимости печи
|
№ |
Статья расходов |
Сумма, руб |
|
1 |
Стоимость основных материалов и готовых узлов |
23700 |
|
2 |
Зарплата основных рабочих |
11200 |
|
3 |
Цеховые расходы |
35840 |
|
4 |
Заводские расходы |
8960 |
|
Заводская себестоимость: |
79700 |
Плановая себестоимость: Спл=1.03·Сзавод=82091 руб
Плановая цена: Цплан=1.03·Спл=84554 руб =3382 $
12. Составление таблицы ТЭП печи
|
№ |
Название параметра |
Единицы |
Значение |
|
1 |
Назначение печи |
– |
Отжиг |
|
2 |
Размеры рабочего пространства |
мм |
110×120 |
|
3 |
Габаритные размеры печи |
мм |
200×256 |
|
4 |
Стоимость печи |
рублей |
79750 |
|
5 |
Масса печи |
кг |
20 |
|
6 |
Максимальная температура |
°С |
1250 |
|
7 |
Температура выдачи металла |
°С |
20 |
|
8 |
Время нагрева, выдержки и охлаждения металла |
ч |
5,2 |
|
9 |
Производительность печи |
кг/ч |
0,96 |
|
10 |
Установленная мощность |
кВт |
10,8 |
|
11 |
Напряжение |
В |
380 |
|
12 |
Число фаз |
– |
1 |
|
13 |
КПД печи |
– |
33% |
|
14 |
Удельный расход энергии |
кВтч/кг |
4,6 |
1. В. Г.Сорокин, М.А. Гервасьев. Стали и сплавы. Марочник.
2. С. Л. Рустем. Оборудование и проектирование термических цехов. Государственное научно-техническое издательство машиностроительной литературы. Москва, 1962.
3. Б.Н. Арзамасов Методика расчёта печей Учебное пособие МГТУ, 1973.
4. А.П. Альтгаузен Электротермическое оборудование: Справочник, М.:Энергия, 1980.