Курсовая работа: Расчет кондиционирования промышленного здания

МИНИСТЕРСТВО ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ

Пермский государственный технический университет

Строительный факультет

Кафедра отопление и вентиляция

Курсовой проект

на тему “Расчет кондиционирования промышленного здания”

Выполнил гр. ПГС_д -07_уск

Волков А.Н.

Руководитель

Матушкина Е.Н.

г. Пермь 2009

1. ИСХОДНЫЕ ДАННЫЕ

Район строительства г.Винница, расчетное количество людей в помещении n = 25. Поступление технологических вредностей:

• тепла Qт=80000кДж/ч;

количество пара прорвавшегося в помещение G_пара =6кг/ч при Р=5 атм. Выполняемая работа – тяжелая.Параметры воздуха в помещении:

• температура воздуха в помещении t_в =22⁰С;

• относительная влажность φ_в =55%;

• рабочий перепад температур в летний период Δt_y =3⁰С;

L_мо =20000м³ /ч Высота помещения Н_р =8м воздух подаваемый в помещение необходимо очищать от пыли. В помещениях смежных с кондиционируемым в теплый период года температура воздуха на 3⁰С выше расчетной наружной температуры, зимой в этих помещениях температура воздуха t_сн =16⁰С.

2. УСТАНОВЛЕНИЕ РАСЧЕТНЫХ ПАРАМЕТРОВ НАРУЖНОГО И ВНУТРЕННЕГО ВОЗДУХА

Расчетные параметры наружного воздуха устанавливаются в зависимости от географического расположения объекта и от вида и назначения систем кондиционирования воздуха. При расчете данной курсовой работы, принимается, что рассматриваемые СКВ первого класса. Вследствие этого, значения наружного воздуха принимаем по параметрам Б для теплого и холодного периодов года.

Значения параметров воздуха в рабочей зоне промышленного помещения равны для летнего и зимнего периодов года и принимаются согласно заданию.

3. БАЛАНС ТЕПЛА И ВЛАГИ В ЛЕТНИЙ И ЗИМНИЙ ПЕРИОДЫ ГОДА

При составлении теплового баланса по полным и явным теплопо ступлениям в летний период года пользуются следующими выражениями:

Здесь Q_т - теплопоступления от технологического оборудования в помещении;

- количество тепла, поступающего в помещение от солнечной радиации:

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

∑F - суммарная площадь остекления, м² ;

- теплопоступления, Вт/м, от солнечной радиации для вертикального заполнения световых проемов.

Здесь количество теплоты соответственно прямой и рассеянной солнечной радиации, Вт/м² .

К_обл - коэффициент облучения, зависящий от углов :

(6)

К_обл.г. и К _обл.в определяется_. в зависимости от углов :

(7)

К_отн - коэффициент относительного проникания солнечной радиации через заполнение светового проема;

τ₂ - коэффициент, учитывающий затенение светового проема переплетами.

Коэффициент инсоляции вертикального заполнения световых проемов:

(8)

Здесь

L_г , L_в – соответственно размер горизонтальных и вертикальных элементов затенения (откосы, стационарные элементы фасада), м;

L_г =0,3, L_в =0,3.

с, а - расстояние от горизонтального и вертикального элементов затенения, до откоса светового проема, м;

с =0, а=0. А_с.о. - солнечный азимут остекления, град, определяется для проемов, ориентированных на С и Ю, а также на В для первой и на З для второй половины дня, по формуле:

Для световых проемов, ориентированных на В для второй половины дня ( после 12 ч) и на З для первой половины дня ( до 12 ч), солнечный азимут определяется по формуле:

Здесь А_о , А_с - соответственно азимут остекления световых проемов и азимут солнца.

β - угол между вертикальной плоскостью остекления и проекцией солнечного луча на вертикальную плоскость, перпендикулярную рассматриваемой плоскости остекления:

где h- высота стояния солнца, град. Теплопоступления, обусловленные теплопередачей, Вт/м² :

Условная температура наружной среды при вертикальном заполнении световых проемов

где t_н.ср. - средняя температура наружного воздуха наиболее жарких суток;- суточная амплитуда температуры наружного воздуха;

β_{2 -} коэффициент, учитывающий гармоническое изменение температуры наружного воздуха;

S_в и d_в - количество теплоты прямой и рассеянной радиации, поступающей в каждый час суток на вертикальную поверхность; ρ_n - приведенный коэффициент поглощения солнечной радиации заполнением световых проемов; α_н - коэффициент теплоотдачи наружной поверхности ограждения, Вт/(м²⁰ С), зависящий от скорости ветра, V, м/с:

Таблица 1. Результаты расчета теплопоступлений от солнечной радиации через вертикальное остекление, ориентированное на Юг.

Аналогично просчитывается теплопоступления через оконные проемы, ориентированные на Восток.

Таблица 2. Результаты расчета теплопоступлений от солнечной радиации через вертикальное остекление, ориентированное на Восток.

Далее определяется расчетный час максимума теплопоступлений в кондиционируемом помещении. Для этого суммируются последние строки табл.1 и 2 (см табл.3)

Таблица 3. Определение максимума теплопоступлений.

Теперь можем найти количество тепла, поступающее в помещение от солнечной радиации по формуле (3)

Q _пар - теплопоступление в помещение с паром:

где G_пар – количество прорвавшегося пара в кг/ч. і_пар – теплосодержание пара, принимается в зависимости от давления пара.

Из исходных данных, давление пара прорвавшегося в здание 5 атм., чему соответствует энтальпия і_пар = 2756 кДж/кг.

Qⁿ _люд ,Q ^я _люд – полные и явные теплопоступления от людей:

где - полные и явные теплопоступления от одного человека, зависит от тяжести работы; n - количество людей.

- теплопоступления из смежных помещений:

кДж / ч (21)

где К_с - коэффициент теплопередачи смежной стены, принимается равным 2,2 8 кДж/(м С); F_с - площадь ограждения, м² ; t_см , t_в - соответственно температуры воздуха в смежном помещении и кондиционируемом.

Для теплого периода года температура в смежных с кондиционируемым помещением

t_см = 27,3+3=30,3⁰С

Для холодного периода года температура в смежных с кондиционируемым помещением t_см =16 ⁰С. Определив все параметры, возвратимся к формуле (1) и (2):

Баланс по полным и явным теплопоступлениям в зимний период года производится по следующим формулам:

где Q_пот = Q_нар + Q_см - потери тепла через наружные и внутренние ограждения. Q_нар - потери через наружные ограждения и покрытие:

где t_B , t_н - соответственно температуры внутреннего и наружного воздуха, К_п =4,85, К_с = 3,2 кДж/(м⁰ С) - коэффициент теплопередачи покрытия и стены, Fп, F_с - площадь покрытия и наружной стены, м.

Q^пом _см в зимний период года учитываются лишь в том случае, когда перепад температур между расчетным помещением и смежным превышает 3⁰ С и рассчитывается по формуле(21), с последующим умножением на(-1).

Q_пот = Q_нар + Q_см =169815,6+7333,2=177148,8кДж/ч

Тогда баланс по полным и явным теплопоступлениям в зимний период года составит:

где W_люд - влаговыделения от людей, определяется как

Здесь n - количество людей в помещении,w - влаговыделения от одного человека.

W_T - поступающее в помещение количество влаги от технологических процессов. По заданию W_Т = кг/ч.

4. РАСЧЕТ ВОЗДУХООБМЕНА В ПОМЕЩЕНИИ И ПОЛНОЙ ПРОИЗВОДИТЕЛЬНОСТИ КОНДИЦИОНЕРА

Расход приточного воздуха L, м³ /ч, определяется расчетом только для летнего периода и принимается наибольший из величин, полученных по формулам (27)-(29):

а) по избыткам явной теплоты:

б) по избыткам полной теплоты:

в) по избыткам влаги ( водяного пара):

где L_мo - расход воздуха, удаляемого из обслуживаемой или рабочей зоны помещения системами местных отсосов, м³ /ч. ∑Q^полн _лет , ∑Q^явн _лет - суммарные полные и явные теплопоступления в летний период года, кДж/ч в помещение; tв, tп,tу - соответственно температуры воздуха в обслуживаемой или рабочей зоне помещения, подаваемого СКВ в помещение и удаляемого из помещения

за пределами обслуживаемой, или рабочей зоны; ∑ W - суммарные поступления влаги в помещение, г/ч. d_B , d_n , d_y - соответственно влагосодержание воздуха в рабочей зоне помещения, подаваемого воздуха и удаляемого воздуха вне рабочей зоны, г/кг; c_B =1,2 кДж/(кГ °С) - теплоемкость воздуха; I_B , I_п , I_y - энтальпии воздуха соответственно в рабочей зоне помещения, подаваемого в помещение и удаляемого из него, кДж/кг.

Формула (27) преобразована к следующему виду:

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

Коэффициент К_воз = 1,45 для рекомендованных к расчету в курсовой работе конструкций воздухораспределителей. Найдем температуру в рабочей зоне помещения:

Так как получены все неизвестные решаем уравнение (30):

Из формулы (31)

найдем температуру удаляемого из помещения воздуха t_у :

Для определения значений d_в , d_п , d_у , I_B , I_п иI_y , используемыхв формулах (28) и (29), воспользуемся I-d диаграммой. Предварительно определяется угловой коэффициент луча процесса из соотношения:

На I-d диаграмме на пересечении изотермы t_B и кривой φ_B находим положение точки В, соответствующей заданному состоянию воздуха в помещении. Через точку В проводим луч процесса. Он проводится следующим образом: на поле I-d диаграммы нанесены риски с указанием численных значений луча процесса. Находим риску, численно равную величине, определенной по формуле (33), и соединяем ее с 0°С на оси температур. Параллельно полученной прямой проводим линию через точку В до пересечения с изотермами t_п и t_y . Точки пересечения П и У имеют координаты приточного и удаляемого из помещения воздуха. С I-d диаграммы считываются значения d и I в точках В, П и У.

После этого определяется воздухообмен по полным теплоизбыткам и влаговыделениям в помещении (формулы (28) и (29)).

За полезную производительность кондиционера L_пп принимается наибольшая из величин L_MO , L_Я , L_п , L_W .

L_пп = L_w =58583,33м³ /ч

При транспортировании воздуха, из-за неплотностей в воздуховодах, возможны утечки, которые должны восполняться за счет увеличения производительности кондиционера на 10-15%.

Для создания подпора в кондиционируемом помещении, необходимо увеличение L_пп . С учетом потерь на утечки и создание подпора в помещении, полная производительность кондиционера L_полн , составит:

где L - количество воздуха, вводимое в помещение на каждого человека, проходящего через дверь, для создания повышенного давления в помещении, м³ /ч; n - количество проходящих людей; V - объем кондиционируемого помещения, м³ ; К - кратность воздухообмена, для создания подпора в помещении; β - коэффициент, учитывающий утечки воздуха.

5. ВЫБОР СХЕМЫ ОБРАБОТКИ ВОЗДУХА И КОМПАНОВКИ КОНДИЦИОНЕРА

По L_полн определяется типоразмер кондиционера. Так как L_полн до 80 тыс. м³ /ч, то выбирается кондиционер КТЦ3-80. Последние две цифры в обозначении марки кондиционера показывают его производительность в тыс.м³ /час.

При расчете системы кондиционирования воздуха, пользуются массовыми единицами расхода воздуха, для чего необходимо использовать следующее выражение:

Здесь ρ _в =1,2кг/м³ – плотность воздуха при стандартных условиях. После перевода L_мо , L_полн иL_пп к массовым единицам, вводим для них соответствующие обозначения - G_мо , G_полн иG_пп :

Так как величина L_MO небольшая из L_Я , L_п , L_W , то применяется схема с рециркуляцией.

При рециркуляции необходимо определить количество наружного воздуха G_нар , кг/ч. Полная производительность кондиционера G_полн :

где G_рец - количество рециркуляционного воздуха, кг/ч.В помещении имеются местные отсосы, тогда

Окончательно скомпоновать кондиционер можно после построения процессов обработки воздуха в I-d диаграмме для летнего и зимнего периодов.

Кондиционер с первой рециркуляцией, построение процесса обработки воздуха на I-d диаграмме.

Летний период

Построение процесса в I-d диаграмме начинаем с нанесения точки В, с параметрами внутреннего воздуха, через которую проводим луч процесса в помещении, до пересечения с изотермой заданной температуры t_п в точке П. Через точку П проводим луч подогрева d_п =const, до пересечения с кривой с φ_о = 90% в точке О. Её параметры соответствуют состоянию воздуха, покидающего КО. Наносим точки Н и В', с координатами наружного и рециркуляционного воздуха, перед входом его в камеру смешивания. Точки В' и Н соединяем прямой линией, которая является линией смеси наружного и рециркуляционного воздуха перед KO. Точка В' находится выше точки В на 1⁰ С по d_в =const.

Положение точки С, соответствующей состоянию смеси воздуха, можно найти из пропорции:

Откладывая от точки В' длину отрезка В' С, находим на прямой смеси В'Н положение точки С. Через точки О и С проводим луч процесса охлаждения и осушения воздуха, до пересечения с кривой φ= 100% в точке W. Координаты точки W- t_w и I_w используются при расчете КО. Все точки процесса обработки воздуха в кондиционере построены, выписываем параметры, характеризующие их:

• линия НВ' - смешение наружного и рециркуляционного воздуха;

• линия СО - обработка воздуха в камере орошения (КО), воздух охлаждается и осушается;

• линия ОП' - нагрев воздуха в секции второго подогрева при d=const, ( секция первого подогрева не работает);

• линия П' П - нагрев воздуха за счет трения о стенки воздуховодов и в результате работы рабочего колеса вентилятора на 1.0 ⁰ С (учитывается только в летний период года);

• линия ПВ - ассимиляция тепла и влаги в помещении;

• Линия ВВ' - нагрев воздуха за счет трения о стенки рециркуляционного воздуховода на 1⁰ С.

Зимний период

(с подмешиванием воздуха перед секцией первого подогрева)

На I-d диаграмму наносят точку В (параметры внутреннего воздуха) и проводят луч процесса

Так как ε<10000, то определяется величина ассимилирующей способности приточного воздуха по влаге:

Зная Δd_p , нетрудно определить величину влагосодержания приточного воздуха:

Пересечение луча процесса ε с линией d_П =const, определяет положение точки П. Через точку П проводят линию d_П =const до пересечения с кривой φ₀ = 90%. Пересечение определяет положение точки О, характеризующее состояние воздуха, покидающего КО. Наносится точка Н и проводится линия ВН. Положение точки смеси С на этой прямой может быть находится из соотношения:

Проведя через точку Случ нагрева в секции первого подогрева по линии d_с =const до пересечения с адиабатой I₀ =const, получим точку К, характеризующую состояние воздуха перед КО. По каждой точке выписываем параметры I, d, t, φ.

По ходу движения воздух в кондиционере с первой рециркуляциией воздух претерпевает следующие изменения:

• линия НВ – смешение наружного и рециркуляционного воздуха;

• линия СК – нагрев воздуха в секции первого подогрева;

• линия КО - обработка воздуха в камере орошения ( КО), воздух охлаждается и увлажняется;

• линия ОП - нагрев воздуха в секции второго подогрева при d=const;

• линия ПВ - процесс взаимодействия приточного воздуха со средой обслуживаемого помещения.

6. РАСЧЕТ ВОЗДУХОНАГРЕВАТЕЛЕЙ (ВН)

(1 подогрева)

Схема обвязки ВН по теплоносителю последовательная по фронту и параллельная по ходу движения воздуха.

1.Типоразмера кондиционера выбираются теплотехнические и гидравлические характеристики ВН.

2.Определяются относительные перепады температур по воздуху Θрв и воде Θрж для расчетных режимов по формулам:

Здесь tв.к.=16,7°С, tв.н.=15,3°С - соответственно конечная и начальная температура воздуха в процессе его обработки в ВН, tж.к.=70°С, tж.н.=130°С - конечная и начальная температура теплоносителя в ВН.

3. Вычисляется относительный расход воздуха G по формуле:

4. На номограмме находят точку с координатами Θ^р _ж и Θ^р _в . Положение построенной точки в поле номограммы определяет выбор поверхности ВН. Все ВН, теплотехнические характеристики которых при G увеличенном в 1,25 раза, расположены справа отрасчетных точек обеспечивают требуемую теплопроизводительность с некоторым запасом.

Минимальный запас по поверхности соответствует самой левой из всех кривых и определяет выбор ВН.

Так как точка построенная по найденным координатам не попадает в диапазон с кривыми, мы прекращаем расчет ВН секции первого подогрева.

РАСЧЕТ ВОЗДУХОНАГРЕВАТЕЛЕЙ (ВН)

(2 подогрева)

Схема обвязки ВН по теплоносителю последовательная по фронту и параллельная по ходу движения воздуха.

1.Типоразмера кондиционера выбираются теплотехнические и гидравлические характеристики ВН.

2.Определяются относительные перепады температур по воздуху Θ^р _в и воде Θ^р _ж для расчетных режимов по формулам:

Здесь t_в.к. =22,25°С, t_в.н. =13,25°С - соответственно конечная и начальная температура воздуха в процессе его обработки в ВН, t_ж.к. =40°С, t_ж.н. =70°С - конечная и начальная температура теплоносителя в ВН.

3. Вычисляется относительный расход воздуха G по формуле:

4. На номограмме находят точку с координатами Θ^р _ж и Θ^р _в . Положение построенной точки в поле номограммы определяет выбор поверхности ВН. Все ВН, теплотехнические характеристики которых при G увеличенном в 1,25 раза, расположены справа от расчетных точек обеспечивают требуемую теплопроизводительность с некоторым запасом.

Минимальный запас по поверхности соответствует самой левой из всех кривых и определяет выбор ВН.

5. Для определения величины запаса по теплообменной поверхности через расчетную точку проводят луч из начала координат до пересечения с выбранной кривой Θ_ж (Θ_в ). Абсцисса и ордината точки пересечения определяют значения Θ^у _ж и Θ^у _в , соответствующих реальному процессу теплопередачи при расчетном расходе воды, определяемом значениями Θ^р _в и Θ^р _ж .

Θ^у _ж =0,535; Θ^у _в =0,165

Запас по поверхности вычисляется по формуле:

где 0,1 - коэффициент, учитывающий допустимое по техническим условиям на кондиционеры КТЦ3 отклонение коэффициента теплопередачи ВН.

6. Находится фактический расход теплоносителя G ^ф _ж , обеспечивающий заданную конечную температуру воздуха:

Где Θ^ф _ж – значение Θ_ж , полученное на точки пересечения вертикальной прямой Θ^р _в с выбранной кривой Θ_ж (Θ_в ); с_в , с_ж - теплоемкость воздуха и теплоносителя.

7.Определяется расчетный расход воды:

8. Вычисляется массовая скорость воздуха в ВН:

значения Fфр=6,63м2 принимаются по номограмме.

9.По номограмме определяется величина потерь давления ΔРа.

ΔРа=26Па

10. Определяются гидравлические потери теплоносителя в ВН, кПа:

где Бп- коэффициент гидравлического сопротивления ВН, принимаемый по номограмме.

7. РАСЧЕТ КАМЕРЫ ОРОШЕНИЯ КОНДИЦИОНЕРОВ КТЦ3

Расчет КО ОКФ-3, двухрядной, исполнения 2, для политропного процесса в летний период года.

1. Вычисляем коэффициент адиабатной эффективности

Ib.k., Ib.h.- энтальпия воздуха перед КО и после нее кДж/кг, Iw - энтальпия воздуха при пересечении кривой ф= 100% с линией процесса изменения состояния воздуха в КО.

2. Находится коэффициент орошения μ и коэффициент энтальпийной эффективности Еп для принятого типоразмера КО по графику.

μ=1,65, Еп=0,55

3. Вычисляется относительный перепад температур воздуха:

4. Вычисляется начальная температура воды:

5. Определяется конечная температура воды:

6. Определяется расход разбрызгиваемой воды:

7.

По графику определяются потери давления по воде в КО - ΔРж=32кПа

8. РАСЧЕТ ВОЗДУХОРАСПРЕДЕЛИТЕЛЯ (ВР)

В качестве воздухораспределителя используем плафон регулируемый, многодиффузорный (ПРМ)

Методика расчета ПРМ:

1.Эффективное воздухораспределение неизотермическими струями обеспечивается для вертикальной подачи при условии:

где Н - геометрическая характеристика струи:

2.Выбирается соотношение сторон (а1,b1) ячейки помещения и ее площадь

А_р =a₁ b₁ =6x6=36м² ,

обслуживаемая одним ВР. Определяется количество ВР в помещении

N=A_пл /A_р = 648/36=18шт

и производительность одного плафона

L_ВР = L_полн /N=72467,67/18=4026м³ /ч

Скорость истечения воздуха из ВР определяют по формуле:

V₀ = L_вр /(3600хA₀ ), м/с

• для ПРМ1 - V₀ =22,37м/с;

• для ПРМ2 - V₀ =13,98м/с;

• для ПРМ3 - V₀ =8,6м/с;

• для ПРМ4 - V₀ =5,59м/с;

Исключаем из дальнейшего расчета плафоны ПРМ1,ПРМ2 и ПРМ3, так как для них скорость истечения воздуха из плафона значительно выше рекомендуемой в 6 м/с.

Определяем геометрическую характеристику струи для ПРМ3 по минимальным и максимальным значениям m и n:

Принимаем среднее Н=16,56м

Определяем высоту расположения плафона относительно рабочей зоны:

Тогда Н=9,45, что не удовлетворяет требования задания.

Принимаем h=5м, (h-h_рз )=3м, Н=6,7м

3.Найдем относительную площадь приточной струи ПРМ4, которая должна находится в диапазоне 0,2-0,5, по формуле:

где для С1=0, С2=-10

Уточним значение коэффициентов m и n

m=1,245, n=1,08

4.Проверим равномерность распределения скоростей температур по площади рабочей зоны:

5. Определяют параметры воздуха при входе в рабочую зону:

где Н_р , h_pз , h - высоты соответственно здания, рабочей зоны и расположения ВР от уровня пола, соответственно, м; А₀ - расчетная площадь ВР, м; m и n - скоростной и температурный коэффициенты ВР, определяемые интерполированием; Δt₀ - избыточная температура приточного воздуха, ^⁰ С; К_соп – коэффициент стеснения струи:

где L_соп - расход воздуха, удаляемого из РЗ, приходящегося на один ВР (при равномерном распределении L_соп = L_мо /N=20000/18=1111,11м³ /ч), м3/ч;

Кп - коэффициент неизотермичности струи:

Тогда параметры воздуха при входе в рабочую зону:

9. ХОЛОДОСНАБЖЕНИЕ ЦЕНТРАЛЬНЫХ СКВ

При использовании КО в центральных СКВ необходимо расположить сборный бак (с отсеками теплой и холодной воды БХВ и БТВ) ниже бака поддона КО. Это дает возможность отепленной воде через переливное устройство самотеком поступать в отсек БТВ.

Далее отепленная вода подается в кожухообразный испаритель холодильной машины. Охлаждаемая вода проходит по трубкам, а в межтрубном пространстве испарителя кипит хладагент. При испарении температура воды понижается на 4-8 С. Охлажденная вода с температурой t = 6-7 °С по соединительному трубопроводу поступает бак холодной воды (БХВ) сборного бака.

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

Объем баков БТВ и БХВ определяются из условия 5-10 часового расхода воды, подаваемого форсунками КО и принимаются равными 0,1 от величины расхода воды. Напор насосов, подающих воду к форсункам, складывается из давления воды перед форсунками, давления, необходимого для преодоления сопротивлений на всасывающей и нагнетательной стороне, а также на преодоление высоты подъема воды.

ЛИТЕРАТУРА

1. СНиП 2.04.05-91* Отопление, вентиляция и кондиционирование. Госстрой СССР - М.:ЦИТП Госстрой СССР, 1991 - 64 с

2. Руководящий материал по центральным кондиционерам ч.ы, Альбом II. Методика расчета ВН., Харьков 1989

3. Руководящий материал по центральным кондиционерам ч.ы, Альбом I. Методика расчета камеры орошения., М.: 198 9

4. Внутренние санитарно-технические устройства. В 3ч. Ч.3 Вентиляция и кондиционирование воздуха Кн.2/ Б. В. Баркалов, Н.Н.Павлов, С. С. Амирджанов и др.; Под ред. Н.Н.Павлова и Ю.И. Шиллера.-4-е изд., перераб. и доп. - М.: Стройиздат, 1992. - 416 с.: ил.-( Справочник проектировщика).