Курсовая работа
«Расчет выпарнойустановки»
Задание
Спроектироватьтрехкорпусную выпарную установку для упаривания 18 т/ч раствора NH4NO3 от начальнойконцентрации 18% до концентрации 48%. Давление греющего пара 6,2 ат.Температура охлаждающей воды 100С, давление в конденсаторе 0,2 ат.Раствор поступает в первый корпус подогретым до температуры кипения. Выпарныеаппараты принять с выносной камерой. Конденсатор смешивающего типа сбарометрической трубой.
Выпарные установки
Выпариваниепредставляет собой термический процесс кипения раствора с выделением пароврастворителей в практически чистом виде (если не считаться с уносом жидкости);при этом растворимое нелетучее вещество (твердое тело, например, соль, иливязкая жидкость, например вазелин) остается в концентрированном виде ваппарате. Получаемые при выпаривании пары удаляются в атмосферу или вконденсирующее устройство.
Превращатьсяв пар растворитель может при кипении жидкости или при поверхностном испаренииее. В выпарных аппаратах применяется более интенсивный из этих способовпревращения растворителя в пар, а именно кипение. впервые выпаривание получилопромышленное применение в производстве сахара, а в дальнейшем – в химическойпромышленности. При концентрировании растворов вода иногда удаляется до 90%первоначального веса.
Предположим,что имеются два открытых сосуда, обогреваемых паром через паровые рубашки. Впервом сосуде находится чистая вода, а во втором – 70%-ный раствор селитры NH4NO3.Пусть давление греющего пара составляет 3,92.105 Па (4кгс/см2). Вода закипит при температуре 100 оС; растворпри том же давлении закипит только при температуре 120 оС. Однако образующиесяиз этого раствора водяные пары будут иметь температур ту же, что и в случаекипения чистой воды, т.е. около 100 оС.
Понижениетемпературы образующихся из раствора водяных паров по сравнению с температуройкипения раствора называют физико-химической температурной депрессией. Обозначивее через ∆t, можем написать:
/>,
где tр– температура кипения раствора, оС; υ – температураобразующихся паров воды, оС.
Физико-химическаятемпературная депрессия различна для разных растворов. Она больше у раствороввеществ с малым молекулярным весом. Для раствора одного и того же веществафизико-химическая температурная депрессия увеличивается с повышением его концентрации.
Подконцентрацией раствора понимают отношение массы сухого вещества в растворе кобщей массе раствора в процентах
/>,
где b– массовая концентрация раствора, %; W – количество растворителя иливоды в растворе, кг; Gсух – количество растворенного илисухого вещества в растворе, кг.
Классификациявыпарных аппаратов и установок
По принципуработы выпарные установки разделяются на периодически и непрерывно действующие.В периодически действующих установках жидкость подается в аппарат, выпариваетсядо необходимой более высокой концентрации, затем упаренный раствор удаляется изаппарата. Опорожненный аппарат вновь наполняется неконцентрированным раствором.Периодическое выпаривание применяется при небольшой производительностиустановки или когда сгущенная жидкость не поддается откачке насосом или в техслучаях, когда требуется выпарить весь растворитель. В аппаратах непрерывногодействия неконцентрированный (слабый) раствор непрерывно отводится из него. Посравнению с аппаратами периодически действующими аппараты непрерывного действияболее экономичны в тепловом отношении, так как в них отсутствуют потери,связанные с расходом тепла на периодический разогрев аппарата.
В большинствеслучаев аппараты непрерывного действия компонуются в так называемыемногокорпусные выпарные установки, в которых упариваемый растворпоследовательно проходит через ряд отдельных аппаратов. В каждом последующемаппарате устанавливается большая концентрация раствора, чем в предыдущем.
По давлениювнутри аппарата различают выпарные аппараты, работающие при избыточном иатмосферном давлениях и вакууме.
Вакуумприменяется в следующих случаях: а) когда раствор под влиянием высокойтемпературы разлагается, изменяется цвет, запах (например, сахар, молоко); б)когда раствор при атмосферном давлении имеет высокую температуру кипения, т.е.обладает большой физико-химической температурной депрессией, и требует высокихпараметров греющего пара (например, раствор аммиачной селитры, едкого калия ит.п.); в) когда греющий теплоноситель имеет низкую температуру и,следовательно, нужно снижать температуру кипения раствора; г) для увеличениярасполагаемого температурного перепада в многокорпусной установке.
Конструкциинаиболее распространенных выпарных аппаратов
1. Аппаратыс паровым обогревом. Наибольшее распространение получили вертикальныевыпарные аппараты с трубчатой поверхностью нагрева, которые хорошо компонуютсяи занимают меньшую площадь.
Во всехконструкциях для облегчения очистки поверхности нагрева от накипеобразованийпар поступает в межтрубное пространство, а раствор подогревается и кипит втрубках.
Выпарныеаппараты с паровым обогревом можно разбить на три группы: с естественнойциркуляцией раствора, с принудительной циркуляцией раствора и пленочныеаппараты.
Аппараты сестественной циркуляцией раствора. Движущей силой естественной циркуляции раствораявляется разность весов столба жидкости в опускных трубах и парожидкостнойэмульсии в подъемных за счет разницы плотностей ρж и ρэ.
Тепловойрасчет
1. Производимтепловой расчет в первом приближении.
Количествоводы, выпариваемой всей установкой:
/>/>
где Gо– количество исходного раствора, кг/ч; bо – начальная концентрация,%; bк – конечная концентрация раствора, %.
На 1 кгначального раствора выпарено />, кг/кгр-ра.
/>/>
Количествоводы, выпариваемое по корпусам, предварительно принимаем одинаковым; W=idem
/>/>/>
Концентрацияраствора на выходе i-ой ступени:
/> или />
/> /> />
2. Примем впервом приближении одинаковые перепады давлений по корпусам и найдем давления вкорпусах
/>
/>
Давлениевторичного пара по корпусам:
/> />
/> />
/> />
3. Полнаяразность температур для всей установки
/> />
где tIГ– температура греющего пара в 1-ом корпусе при заданном давлении. Находится порр1 на линии насыщения; />= f(рп3)– температура вторичного пара на выходе 3-й ступени (температура пара вконденсаторе при давлении в конденсаторе). Температуры взяты по таблицамнасыщенных водяных паров М.П. Вукаловича.
Потери общейразности температур определяем как сумму депрессионных физико-химическихпотерь, потерь от гидростатического эффекта и гидравлических потерь втрубопроводах.
Потери температурв установке:
/>
/>13,81+28,07+3=44,88
где ∆ji– потери температурного напора по ступеням.
а) Σ∆1 – потери общей разноститемператур за счет физико-химической депрессии. Для i-ой ступени:
/>
/> /> />
где Т –температура кипения воды при данном давлении, К; r – теплота парообразованияпри данном давлении, кДж/кг; ∆i1Н – нормальнаятемпературная депрессия (при нормальном давлении 760 мм рт. ст.);определяется по таблицам 2.22, 2.24 │7│, определенаэкспериментально и затабулирована для различных веществ.
Или ∆1iопределяется по упрощенной формуле Тищенко
/>,
где К =f(t) – поправочный коэффициент, принимается по табл. 2.2.
Общиедепрессионные физико-химические потери определяются по формуле:
/>
/>
б) Σ∆2– потери от гидростатического эффекта; зависят от высоты уровня раствора,плотности ρ парожидкостной эмульсии и скорости циркуляции.
Для i-ойступени:
/>
/> 147,18 />144,76
/> 126,79 />=122,7
/> 81,27 />59,8
/>147,18–144,76=2,42
/>126,79–122,7=4,09
/>81,27–59,8=21,47
где tк.с.в-температуракипения воды при давлении />; tк.в- температура кипения воды при давлении рвт; рвт– давление вторичного пара над раствором; ∆рг –гидростатическое давление раствора у середины греющих труб.
Давлениераствора в середине греющих труб больше, чем давление пара на величину ρэgh.
Определяемгидростатическое давление раствора у середины греющих труб ∆рг:
/>,
/>0,5 м />,
/> 4 м />,
/> 0,5 + 4/2 =2,5 м />,
где />; ρi– плотность раствора в зависимости от его концентрации, находится по графику2.9 или таблицам 2.21, 2.23 │7│; h – расстояние от верхнегоуровня раствора до середины греющих труб; hизб–расстояние от уровня раствора до трубной доски, принимаем 0,25–0,5 м; hтр– высота греющих труб, принимается 3–5 м.
Давлениевторичного пара по корпусам было определено выше.
Давлениерастворов у середины греющих труб:
/>
/>
/>
/>
Общие потериза счет гидростатического эффекта:
/>
/>2,42+4,09+21,47=27,98
в) Гидравлическаятемпературная депрессия связана с потерями давления при движении пара потрубопроводам. В выпарных установках гидравлические потери при прохождении параиз парового пространства предыдущего корпуса в греющую камеру последующегосоставляют 1,0–1,5 оС.
Принимаем:
λ=0,03; l=15 м; W=20 м/с; ρ=2,2кг/м3; d=0,3 м. Из этого Δр будет равно примерно 1 кПа, чтосоответствует потерям в 0,1–0,2 0С, но по опыту предыдущих расчетовпринимаем:
/> оС
Общиегидравлические потери: />= 3 оС.
4. Полезнаяразность температур для всей установки:
/>
/>99,81–43,56=56,25
Будемпроектировать установку исходя из равенства поверхностей нагрева по ступеням установки,тогда суммарная полезная разность температур должна быть распределена поступеням пропорционально отношениям тепловых нагрузок к коэффициентамтеплопередачи.
Изпрактических данных эксплуатации установок известно:
КI:КII: КIII = 1: 0,7: 0,4
Примем приэтом, что количество тепла, передаваемое через греющую поверхность будет равнымдля всех корпусов: QI = QII = QIII.
Тогдаполезная разность температур i-ой ступени:
/>
При нашихдопущениях имеем:
/>; />; />
/>; />;
/>
Проверить(округлить), чтобы />.
5.Температура кипения раствора у середины греющих труб и у верхнего уровня в 1-ойступени:
/>при 6,2 ата по таблицеВукаловича
/>
/>159,61–11,413=148,197
/>
/>148,197 – 2,42=145,777
Температуравторичного пара в 1-ой ступени:
/>/>
/>/>145,777 – 2,66=143,117
Температурагреющего пара во 2-ой ступени:
/>
/>143,117 – 1=142,117
Температуракипения раствора у середины греющих труб и у верхнего уровня во 2-ой ступени:
/>
/>142,117 – 16,304=125,813
/>
/>125,813 – 4,09=121,723
Температуравторичного пара во 2-ой ступени:
/>
/>121,723 – 4,63=117,093
Температурагреющего пара в 3-ей ступени:
/>/>
/>117,093 – 1=116,093
Температуракипения раствора у середины греющих труб и у верхнего уровня в 3-ей ступени: />
/>116,093 – 28,533=87,56
/>
/>87,56 –21,47=66,09
Температуравторичного пара в 3-ей ступени:
/>
/>66,09 – 5,29=60,8
Температурапара в конденсаторе:
/>
/> 60,8 – 1=59,8
Полезныеперепады температур /> должны быть неменее
10–15 оСпри /> ≤ 2.10-6Па.с,
18–24 оСпри 2.10-6
где /> – динамический коэффициентвязкости раствора при средней концентрации.
Потемпературам паров находим в таблицах М.П. Вукаловича давления и энтальпиипаров и конденсата, а по концентрациям растворов определяем теплоемкости иинтегральные теплоты растворения (берем из графиков), полученные данные сводимв таблицу.Наименование параметров Обозначение Корпус (ступень) I II III Концентрация, вх/вых, % b 18/22,7 22,7/30,9 30,9/48
Полезная разность температур, оС
∆tп 11,413 16,304 28,533
Температура греющего пара, оС
tн, />/> 159,61 142,117 116,093 Температура кипения раствора у середины греющих труб
tкс 148,197 125,813 87,56
Температура кипения раствора у верхнего края труб, оС
tк 145,777 121,723 66,09
Гидростатические потери, оС
∆2 2,42 4,09 21,47
Физико-химическая дисперсия, оС
∆1 2,66 4,63 5,29
Гидравлические потери в трубопроводах, оС
∆3 1 1 1
Температура вторичного пара, оС
/> 143,117 117,093 60,8 Давление греющего пара, атм
рГ 6,2 3,89 1,77 Энтальпия греющего пара, ккал/кг
hГ 658,59 653,67 645,13 Энтальпия конденсата, ккал/кг
hк 160,93 142,89 116,38 Давление вторичного пара, атм
рвт 4 1,83 0,21 Энтальпия вторичного пара, ккал/кг
hвт 653,97 645,48 623,62 Теплоемкость раствора, вх/вых, кДж/кг*град
сi 3,6/ 3,4 3,4/ 3,1 3,1/ 2 Интегральная теплота растворения, кДж/кг
qRн/qRк -100/-120 -120/-180 -180/-215
∆qR= qRн-qRк 20 60 35
6. Удельныйрасход пара на выпаривание 1 кг раствора без учета теплоты дегидратации иравенстве = 1 коэффициента самоиспарения во всех корпусах
/> кг/кг р-ра
/>
/>
β1 = 0, если раствор вводятв 1-ую ступень с температурой кипения.
/>; /> – это коэффициентысамоиспарений.
/>; />
Расход парана 1, 2 и 3 ступени, кг/с:
/>, кг/с
/>
/>
Расходы параможно подсчитать и следующим образом. Расчет начинаем с 3-ей ступени.
/>
/>
где kзап= 1,03 – коэффициент, учитывающий потери тепла аппаратом в окружающую среду; снi,скi – теплоемкости раствора при начальной и конечной концентрациираствора в ступени аппарата (по составленной нами таблице); ∆qRi– разность интегральных теплот растворения вещества между существующей ипредыдущей концентрациями растворов.
/> – количество раствора,перетекающего из одной ступени в другую. Таким образом:
/> – для 1 ступени; />18000 – 3750=14250
/> – для 2 ступени; />14250 – 3750=10500
/> – для 3 ступени. />10500 – 3750=6750
7. Количествотеплоты, передаваемой через поверхность нагрева i-го корпуса
/>, кВт
/> 3422,61 (645,13 – 116,38)=1809705ккал/ч =2104867 Вт
8.Коэффициент теплопередачи от конденсирующегося пара к кипящему раствору
/>
/>
/>
/>
Проверка: q3=∆tп3*k3=28,533*1041,2=29708
/>
/> м²
/>3422,61 />
/> ккал/ч =2043364Вт
/>
/>
Проверка: q2=∆tп2*k2=16,304*1408,1=22958
/>м²
/>3439,79
/>
/> ккал/ч=1922588 Вт
/>
/>
Проверка: q1=∆tп1*k1=11,413*1743,5=19898
/>м²
9.Ориентировочные значения поверхности нагрева i-го корпуса:
/>
Если Fiдалеки друг от друга, или W1≠D2, а W2≠D3,то сделать перерасчет.
Произведемперерасчет количества воды, выпариваемой по ступеням:
в 1-ойступени:
/>
/>
/>
во 2-ойступени:
/>
/>
в 3-ойступени:
/>
/>
Для полученияболее точного значения поверхностей нагрева произведем расчет во второмприближении.
Концентрациирастворов:
в 1-м корпусе
/>
/>
/>
/>
Общаяразность температур находится по той же формуле, что и в первом приближении.
/>=99,81
Гидростатическиедавления растворов у середины греющих труб:
/>
/>,
/>,
/>
Плотности ρiвзяты при новых концентрациях растворов по корпусам.
Давлениявторичных паров по корпусам будут те же, что в первом приближении.
Давлениярастворов у середины греющих труб пересчитываются по известной формуле:
/>
Далее всепересчитываем по уже известным формулам, но подставляя новые значения,полученные при пересчете.
Всеполученные данные сводим в таблицу, как и при расчете в первом приближенииНаименование параметров Обозначение Корпус (ступень) I II III Концентрация, вх/вых, % b 18/22,6 22,6/30,2 30,2/48
Полезная разность температур, оС
∆tп 13,306 17,515 24,399
Температура греющего пара, оС
tн, />/> 159,61 140,134 112,489 Температура кипения раствора у середины греющих труб
tкс 146,304 122,619 88,09
Температура кипения раствора у верхнего края труб, оС
tк 143,764 117,949 66,13
Гидростатические потери, оС
∆2 2,54 4,67 21,96
Физико-химическая дисперсия, оС
∆1 2,63 4,46 5,33
Гидравлические потери в трубопроводах, оС
∆3 1 1 1
Температура вторичного пара, оС
/> 141,134 113,489 60,8 Давление греющего пара, атм
рГ 6,2 3,68 1,58 Энтальпия греющего пара, ккал/кг
hГ 658,59 653,07 643,85 Энтальпия конденсата, ккал/кг
hк 160,93 140,86 112,73 Давление вторичного пара, атм
рвт 3,78 1,63 0,21 Энтальпия вторичного пара, ккал/кг
hвт 653,38 644,21 623,62 Теплоемкость раствора, вх/вых, кДж/кг*град
сi 3,6/ 3,4 3,4/ 3,1 3,1/ 2 Интегральная теплота растворения, кДж/кг
qRн/qRк -100/-120 -120/-180 -180/-215
∆qR= qRн-qRк 20 60 35
/>
/> – для 1 ступени; />18000 – 3646,46=14353,54
/> – для 2 ступени; />14353,54 – 3628,24=10725,3
/> – для 3 ступени. />10725,3 – 3975,3=6750
Количествотеплоты, передаваемой через поверхность нагрева i-го корпуса
/>, кВт
/> 3691,47 (643,85–112,73)=1960612ккал/ч=2280192 Вт
Коэффициенттеплопередачи от конденсирующегося пара к кипящему раствору
/>
/>
/>
/>
Проверка: q3=∆tп3*k3=24,399*1040,64=25222
/>
/> м²
/>3691,47/>
/> ккал/ч =2195706 Вт
/>
/>
Проверка: q2=∆tп2*k2=17,515*1408,36=24667
/>м²
/>3685,92
/>
/> ккал/ч =2037887 Вт
/>
/>
Проверка: q1=∆tп1*k1=13,306*1730,81=22490
/>м²
Средняяповерхность нагрева:
/>м²
Проектированиеаппарата
По табл. 2.16|7| принимаем поверхность нагрева F=80 м2; диаметр корпусааппарата Dвн=0,8 м. Число труб в греющей камере:
/>,
/>
где Нтр– длина (высота) трубки, м. Нтр = 3 м – подвесная камера; Нтр= 4 м – с выносным сепаратором; dср – средний диаметр трубок, d= 38÷50 мм.
Произведемрасчет штуцеров выпарного аппарата. Диаметр штуцера определим по формуле
/>,
где /> – объемный расходтеплоносителя, м3/сек; G – массовый расход теплоносителя,кг/ч; γ – плотность пара, кг/м3; w – скоростьпара, м/сек.
Скорость парапринять 20 м/сек.
Расчетысводим в табл.
Таблица расчетовштуцеров выпарной установкиНаименование штуцера Расход пара, кг/ч Давление пара, ат
Плотность, кг/м3
Секундный расход, м3/с Скорость пара, м/с Диаметр, мм расчетный принятый Вход греющего пара 3521 6,2 3,24 0,302 20 138 150 Выход вторичного пара…………. 3685,9 4,47 2,21 0,46 20 169 170 Вход раствора 18000 1076 0,005 1 76 80 Выход упаренного раствора………. 6750 1218 0,0015 0,5 62 70
Расчетбарометрического конденсатора
Определяемконечную температуру охлаждающей воды при давлении в конденсаторе/> ат, tп =59,8º C, удельный объем пара uп = 7,749 м3 /кг.Температура охлаждающей воды t´2 =10 0С.
Температуравыходящей охлаждающей воды меньше tп на δ = 1 – 30 С; вследствиенесовершенства теплопередачи принимаем δ = 30 С, тогда
t»2 = tп – δ = 59,8 – 3=56,8º C.
Кратностьохлаждения составляет
m =W/D = (i – t»2)/(t»2 – t´2)= (623,62 – 56,8)/(56,8 –10)=12,112 кг/кг
Часовойрасход охлаждающей воды при количестве конденсируемого пара после 3-го корпусасоставляет:
W=Dm=3691,47 ∙12,112 =44709,4кг/ч
Диаметрконденсатора при скорости в конденсаторе ω, равной 15 м/сек:
dк =0,0188∙√(D∙uп/ω)=0,0188∙√(3691,47∙7,749/15)=0,79 м
Принимаем dк =800 мм.
Согласно табл.2–20 барометрический конденсатор имеет следующие размеры: высота конденсатора H= 5088 мм, ширинаполки b=500 мм и высота борта равна 40 мм. Число полок – 6.
Диаметрбарометрической трубы определяем из расчета на пропуск смеси воды и конденсата.Из уравнения
W + D = (πd2/4)∙ω,
полагая ω= 1 м/сек, получаем:
/>м
Принимаем d=150 м.
Высотаводяного столба, соответствующая заданному вакууму,
H1=10,33∙B/760 = 10,33∙560/760=7,6 м
Принимаемпредварительно полную высоту трубок H=9 м
ЧислоРейнольдса для трубок при коэффициенте кинематической вязкости воды притемпературе 59,8 º C, равном ν = 0,517 м2/с
Re = ω∙d/ν = 1∙0,15∙106/0,517=232108,3
Коэффициент трениядля гладких труб при значениях Re = 105 — 103 определяетсяпо формуле Никурадзе
λ =0,0032 + 0,221/(Re0.237)=0.0032+0,221/(232108,3)0,237=0,015/>
Потеря напорана трение и местные сопротивления в барометрической трубе
H2=/> м.вод. ст.
где d и l – диаметр и длинабарометрической трубы; 2,5 – коэффициент, учитывающий потери на местныесопротивления.
Полная высотатрубы
H=H1+H2+H3=7,6+0,18+0,5=8,28 м
где H3= 0,5 м – поправка,учитывающая возможные колебания вакуума в конденсаторе или уровня воды в водоприемнике.
Принимаемвысоту трубы Н=9 м
Определениепроизводительности вакуум-насоса:
GВ=(0,25*(D+W)+100D)/10000=34,4 кг/ч
tВ=10+4+0,1 (59,8–10)=18,980С
рК=0,2*100000=2000 ммвд. ст.
рП=200 ммвд. ст.
рВ=2000–200=1800 ммвд. ст.
VВ=(29,27*GВ*(273+tВ))/рВ=163,4 м3/ч=2,7 м3/мин
Принимаемротационный водокольцевой вакуум-насос РМК-3 производительностью 5 м3/мин.
Проведемрасчет выпарного аппарата на прочность.
Толщинастенок цилиндрической обечайки греющей камеры:
/> см
Принимаем S=10 мм
σдоп=1340*0,9=1206кгс/см2==118 МПа
Толщинастенок цилиндрической обечайки сепаратора:
/>см
Принимаем S=10 мм
Толщинаверхней крышки сепаратора:
/>см
Принимаем S=10 мм
Толщинакрышки корпуса:
/>см
Принимаем S=10 мм
Толщина днищасепаратора:
Hэкв=3,79/0,001071=3725,5 см=3,7 м
h=1000 мм=100 см
Нобщ=Нэкв+h=3826 см ст.жидкости
/>см
Принимаем S=10 мм
Проверканеобходимости крепления вырезов под патрубки:
Максимальнодопустимый диаметр неукрепленного отверстия в греющей камере:
/>
/>
Вырез вгреющей камере диаметром d=400 мм надо укрепить.
Максимально допустимыйдиаметр неукрепленного отверстия в сепараторе:
/>
/>
Площадьукрепления выреза для патрубка d=400 мм кольцом толщиной δ=12 мм:
Fукр= δ*d+ S*d+2*a2=12*400+10*400+2*702=18600 мм2
Площадьотверстия выреза:
Fотв=S/0,9*(2*d-50)=10/0,9*(2*400–50)=8333,3 мм2
Fукр >Fотв
Литература
1. П.Д. Лебедев.Теплообменные, сушильные и холодильные установки. – М.: Энергия, 1972. – 320 с.
2. А.М. Бакластов и др.Промышленные тепломассообменные процессы и установки. – М.: Энергоатомиздат,1986. – 327 с.
3. А.М. Бакластов.Проектирование, монтаж и эксплуатация теплоиспользующих установок. – М.:Энергия, 1970. – 568 с.
4. Б.Н. Голубков и др.Теплотехническое оборудование и теплоснабжение промышленных предприятий. – М.:Энергия, 1979. – 541 с.
5. Теплотехническийсправочник. Под ред. В.Н. Юренева и П.Д. Лебедева. Т. 1 и 2. – М.:Энергия, 1975 и 1976. – Стр. 743 и 896.
6. Н.Б. Варгафтик.Справочник по теплофизическим свойствам газов и жидкостей.
7. П.Д. Лебедев, А.А. Щукин.Теплоиспользующие установки промышленных предприятий. Курсовое проектирование.– М.: Энергия. 1970. – 408 с