Размещено на http://www.allbest.ru/

Отопления и вентиляции свинарника

1. Теплотехническая оценка наружных ограждений

1.1 Проверка теплозащитных свойств наружных ограждений

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

Чтобы удовлетворить всем этим требованиям, необходимо правильно подобрать и толщину самой стены, и толщины слоев утеплителей, которые приходится применять, чтобы обеспечить выполнение названных требований.

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

где R_в=; , - толщина и теплопроводность каждого слоя изоляции стены; R_н=; _в и _н - нормированные значения коэффициентов теплоотдачи от внутреннего воздуха к стене и от стены к наружному воздуху, [4].

R₀=1/8,7+0,64/0,7+0,043 = 1,07 м²К / Вт

Чтобы удовлетворить условиям теплоснабжения, величину требуемого термического сопротивления R_тр, меньше которого R₀ не может быть принята, определяется по нормам в зависимости от величины расчетного значения градусо-суток отопительного периода , рассчитываемой по формуле:

=(, (1.2)

где t_в - температура воздуха внутри помещения, определяемая в зависимости от назначения помещения по нормативным документам, °С; t_оп - средняя температура наружного воздуха за время отопительного периода в заданной климатической зоне, °С; Z_on - продолжительность отопительного периода в сутках [2].

= (16 + 4,4)253= 5161,2^оС•сут.

Найдем величину R_тр по приложению П6 [3]. R_тр=2,8 Сопоставив R₀ и R_тр (1,07<2, 8), вычислим невязку

Рассчитаем толщину слоя утеплителя по формуле:

(1.3)

В качестве утеплителя примем энергофлекс (?=0,038

1.2 Проверка на отсутствие конденсации на внутренней поверхности наружных стен

Наиболее вероятно выпадение конденсата в угловых помещениях, а также в помещениях с низкой температурой внутреннего воздуха. Температура внутренней поверхности ограждения определяется из выражения:

где t_в - расчетная температура внутреннего воздуха, принимаемая в зависимости от назначения помещения, °С; t_н - расчетная зимняя температура, равная -31 ^оС; R₀ - общее термическое сопротивление, ; R_в - сопротивление теплоотдаче внутренней поверхности, .

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

где _в-относительная влажность внутреннего воздуха.

Найдем общее термическое сопротивление наружной стены по формуле (1.1):

R₀=1/8,7+0,64/0,7+0,07/0,038+0,043 =2,9 м²К / Вт

Определим температуру внутренней поверхности наружных стен помещений по формуле (1.2), в которых температура внутреннего воздуха 20^оС:

Вычислим температуру точки росы для этих помещений по формулам (1.5) - (1.7):

,

.

Т.к. температура внутренней поверхности больше температуры точки росы (14,1>5,5), то конденсации на внутренней поверхности наружных стен не будет.

2. Тепловой баланс помещения

Тепловой баланс помещения сводится к определению тепловой мощности системы отопления, затраченную на компенсацию теплонедостатка в помещении, и находится по следующей формуле [8]:

Q_со=?Q_пот - ?Q_пост

где Q_со - расчетная мощность системы отопления, ?Q_пот - суммарные тепловые потери помещения, ?Q_пост - суммарные теплопоступления в помещении.

Так как по заданию дано жилое здание, то расчетная формула упрощается:

Q_со= Q_огр+ Q_и - Q_быт - Q_ч (4)

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

2.1 Теплопотери через ограждающие конструкции

Теплопотери через ограждающие конструкции условно разделяют на основные Q_осн и добавочные Q_доб [6]:

Q_огр= Q_осн+ Q_доб (5)

Основные теплопотери определяются теплопередачей через конструкцию ограждения и рассчитываются по следующей формуле [6] с учетом добавочных потерь на ориентацию наружной стенки:

Q_осн=F (t_в - t_н) (1+ ??)nk_о (6)

где F - площадь ограждения; n - коэффициент, зависящий от положения наружной поверхности ограждения по отношению к наружному воздуху; ? - коэффициент, зависящий от ориентации поверхности по сторонам света.

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

Q_осн= Q_нс+ Q_пп + Q_оп + Q_ок (7)

Расчет основных теплопотерь через наружные стены

Определим теплопотери через наружные стенки комнаты №1. Для этого рассчитаем тепловые потери каждой наружной стены комнаты в отдельности, то есть:

Q_нс = Q_осн1 + Q_осн2 (8)

По заданию высота комнаты составляет h=3,5 м; длина - a= 2,14 м; ширина - b₁ = 3,7 м.

Рассчитаем тепловые потери каждой наружной стенки отдельно. Размер первой стенки 2,14х3,7 м. Площадь стены будет равна:

F =3,73,5= 12,95 м².

Так как стенка угловая, то температуру наружного воздуха будем принимать равной нормируемой температуре внутри здания плюс 2 градуса [2]:

t_в = 16 + 2 = 18 °С (9)

Стена ориентирована на запад, значит коэффициент ? будет равен ?=0,05 [2]. Так как стена наружная, то коэффициент n=1 [2].

Значение коэффициента теплопроводности стены k_о из выше найденных данных будет иметь вид:

k_ст = 1/R_о = 1/2,9 = 0,34 Вт/м²К

Подставим все полученные значения в формулу (6) и получим значение тепловых потерь через первую стенку Q_осн1:

Q_осн1 = 12,95 (18 - (-31)) (1+0,05)10,34 = 226,53 Вт.

Размер второй стенки 2,14х3,5. Площадь стены будет равна:

F = 2, 143, 5= 7,5 м².

Стена ориентирована на юг, поэтому ?=0 [2]. Стена наружная и угловая, значит n=1 [2] и t_в=18 °С. Коэффициент теплопроводности равен k_ст = 0,34 Вт/м²К.

Подставим в формулу (6) вышеуказанные данные и найдем тепловые потери через вторую стенку Q_осн2:

Q_осн2 = 7,5 (18 - (-31)) (1+0)10,34 = 124,95Вт.

По формуле (8) определим потери тепловой энергии через наружные стены:

Q_нс = 226,53+ 124,95 = 351,9 Вт.

Расчет основных теплопотерь через окна

Определим потери тепловой энергии через окна. Так как, окна не предусмотрены конструкцией потерь через окна нет.

Расчет основных теплопотерь через пол

Теперь определим тепловые потери через половое покрытие дома, бетон покрытый слоем асфальта [2]. Разобьем весь пол на две полосы, параллельные наружным стенам и равные по ширине а=2 м. Пол состоит из материалов, коэффициент теплопроводности которых ?>1,16, поэтому будем считать, что полы не утепленные.

Общие потери через пол буду определяться суммой потерь через первую, вторую и третью полосы:

(10)

где F₁ - площадь первой полосы, F₂ - площадь второй полосы, F₃ - площадь третей полосы.

Определим площадь первой полосы F₁:

F₁ = 2,142 + 3,72 = 11,7 м²

Определим площадь второй полосы F₂:

F₂ = 0,141,7 = 0,23 м²

Определим площадь третьей полосы F₃:

F₃ = 0 м²

По формуле (10):

= 276,5 Вт.

Расчет основных теплопотерь через потолочное перекрытие

Посчитаем потери тепловой энергии через потолочное перекрытие и крышу. Размер потолочного перекрытия 3,7х2,14 [7]. Площадь потолка F_п =7,92 м². Коэффициент n для чердачных перекрытий будет равен n=0,9 [2]. Коэффициент ? для потолка не учитывается.

По формуле (1) определим термическое сопротивление потолочного перекрытия R_пп:

R_пп = 1/8,7+0,043+0,04/0,58+0,14/0,18+0,005/0,17+0,06/1,69= 1,07 м²К / Вт

Коэффициент теплопроводности перекрытия k_пп будет равен:

k_пп = 1/R_пп = 1/1,07= 0,93 Вт/м²К.

По формуле (6) определим тепловые потери через потолочное перекрытие Q_пп:

Q_пп =7,92 (18 - (-15))0,90,7= 182,95 Вт.

Теперь по формуле (7) найдем все основные тепловые потери комнаты №1 Q_осн:

Q_осн = 351,9 + 276,5+ 327= 955,4 Вт.

2.2 Расчет добавочных теплопотерь через ограждающую конструкцию

Рассчитаем добавочные тепловые потери Q_доб, они будут складываться из потерь на врывание холодного воздуха через входные двери Q_вхв; количество тепла, необходимое для нагревания воздуха в объеме, обеспечивающем минимальный нормативный воздухообмен здании Q_в [6].

Q_доб = Q_вхв + Q_в (13)

Потери на врывание холодного воздуха через входную дверь Q_вхв для комнаты №16 рассчитаем как [8]:

Q_вхв = 0,8Q_дв (14)

где Q_дв - теплопотери через дверное ограждение.

Размер двери 2,6х2,4х0,06 [7]. Площадь двери будет равна F_дв = 6,24 м². Коэффициент n для дверного перекрытия будет равен n=1 [2]. Коэффициент ?=0,1 [2], так как дверь ориентирована на северо-запад.

Рассчитаем теплопотери через дверные ограждения комнаты №9. Примем входную дверь за дверь, сделанную из сосны, коэффициент теплопроводности которой равен ?=0,14 Вт/м²К [10]. Тогда термическое сопротивление R_дв двери найдем по формуле (1):

R_дв = 1/8,7+0,06/0,14+1/12 = 0,63 м²К / Вт

Коэффициент теплопроводности перекрытия k_пп будет равен:

k_дв = 1/R_дв = 1/0,63 = 1,59 Вт/м²К.

Теплопотери через дверные ограждения Q_дв будут равны:

Q_дв= 6,24 (16 - (-31))11,59 (1+0,1) = 512,94 Вт.

По формуле (14) определим тепловые потери на врывание воздуха через входную дверь Q_вхв:

Q_вхв = 0,8512,94 = 410,36 Вт.

Количество тепла, необходимое для нагревания воздуха в объеме, обеспечивающем минимальный нормативный воздухообмен в комнате №1 Q_в [8]:

Q_в = (t_вн - t_н)F_п (18)

где F_п - площадь пола.

Размер комнаты №1 3,7х2,14 [7]. Площадь всего пола равна F_п = 7,92 м².

По формуле (18) определим количество тепла, необходимое для нагревания воздуха в объеме, обеспечивающем минимальный нормативный воздухообмен, Q_в:

Q_в = (18 - (-31))7,92 = 388,1 Вт.

Для комнаты №1 добавочные потери тепловой энергии будут равны:

Q_доб = 388,1 Вт.

Теперь по формуле (5) определим полные тепловые потери комнаты №1 Q_огр: Q_огр = 955,4 + 388,1 = 1343,5 Вт.

2.3 Расчет тепла на нагрев воздуха, поступающего инфильтрацией

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

Q_и = 0,278 (А_оG_o (F_ок + F_дв)) (t_вн - t_н) (15)

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

Расход воздуха поступающего путем инфильтрации через световые проемы окон и балконных дверей, G_o определяется по формуле [6]:

G_o = (16)

где ?р - разность давлений воздуха на наружной наветренной и внутренней поверхностях ограждающих конструкций; R_вп - сопротивление воздухопроницанию световых проемов.

Разность давлений воздуха ?р определяется по следующему выражению [6]:

?р = 9,8 ((H - h_o) (?_н - 1,27) + 0,05?_н? ² (С_н - С_з)k) (17)

где H - высота от поверхности земли до верха вытяжной шахты; h_o - высота от поверхности земли до центра окон или дверей; ?_н - плотность воздуха при наружной температуре; ? - скорость ветра; С_н и С_з - аэродинамические коэффициенты для наветренной и заветренной поверхностей соответственно; k - коэффициент, учитывающий изменение скорости напора в зависимости от высоты здания и типа местности.

Определим основные размеры по чертежу: H=7 м; h_o=1,2 м [7]. Плотность воздуха ?_н=1,4 кг/м³ [6]; скорость ветра для г. Архангельск ? = 3,7 м/с [9]; аэродинамический коэффициент для наветренной С_н=0,8 [6] и заветренной С_з=0,6 поверхностей; коэффициент, учитывающий изменение скорости напора, примем k =0,75 [6].

Подставим все значения в формулу (17) и рассчитаем разность давлений воздуха ?р:

?р = 9,8 ((7-1,2) (1,4-1,27)+0,051,43,7² (0,8-0,6)0,75) = 8,8 кг/м³.

Сопротивление воздухопроницанию проемов дверей R_вп= 1,2 м²ч Па / Вт [10].

Из формулы (16) найдем расход воздуха поступающего путем инфильтрации G_o:

G_o= 3,55 кг/ч.

Определим общую площадь дверей. Количество дверей в комнате №16 - Z = 1, [7]. Площадь дверей F_дв будет равна 6,24.

Коэффициент А_о примем для дверей с раздельными переплетами А_о=0,9 [6].

Подставим все найденные значения в формулу (15) и определим потери теплоты на подогрев воздуха, поступающего инфильтрацией, Q_и:

Q_и = 0,278 (0,93,55 (6,24)) (16 - (-31)) = 260,5 Вт.

2.4 Расчет бытовых теплопоступлений

Бытовые тепловыделения в комнате №1 рассчитаем по следующей формуле [8]:

Q_быт = 21F_п = 217,9 = 165,9 Вт.

Зададим количество свиней в комплексе равное n = 50. Тепловыделения от одного животного определим, предварительно приняв коэффициент, учитывающий интенсивность движения животного, ?₁ равным ?₁=1 [6], коэффициент, подвижность воздуха в помещении равным = 0,1 [6], по формуле [6]:

Q_ж1 = ?₁ (2,5 + 10,3) (35 - t_в) = 1 (2,5 + 10,30,1) (35 - 16) = 67,1 Вт

Q_ж = 67,1 50=3355 Вт

3. Расчет воздухообмена в помещении

Воздухообмен в помещении призван поддерживать в заданных пределах содержание влаги, СО₂ и других вредных примесей в воздухе помещения. Расчетный воздухообмен L находят суммированием соответствующих составляющих: воздухообмена при условии удаления из помещения избыточной влаги L_W и воздухообмена для понижения концентрации углекислоты L_СО2.

L = L_W + L_СО2 (18)

Воздухообмена при условии удаления из помещения избыточной влаги L_W определим по следующей формуле [2]:

L_W = (19)

где W_ч - масса водяных паров, выделяющихся в помещении от свиньи; d_в - влагосодержание внутреннего воздуха; d_н - влагосодержание наружного воздуха; ? - плотность воздуха при температуре помещения; n - колличество свиней в помещении.

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

W_ч = 1250 г./ч

Количество свиней в комплексе n = 50.

Плотность воздуха при температуре внутри помещения равной t_вн = 16 ?С будет равна:

? = 1,369 кг/м³

Влагосодержание внутреннего воздуха определим по i - d диаграмме воздуха [2]:

d_в = 3,5 кг/м³

Влагосодержание наружного воздуха определим по i - d диаграмме воздуха [2]:

d_н = 0,1 кг/м³

Подставим значения в формулу (19) и определим воздухообмен при условии удаления из помещения избыточной влаги:

L_W = = 13428 м³/ч

Воздухообмена для понижения концентрации углекислоты L_СО2 найдем из следующей формулы [2]:

L _СО2 = (20)

где с - колличество СО₂, выделяемое одним животным; с₁ - предельно допустимая концентрация СО₂ в воздухе помещения; с₂ - концентрация СО₂ в наружнем воздухе.

Колличество СО₂, выделяемое одним животным равно:

с = 150 л/ч

Предельно допустимая концентрация СО₂ в воздухе помещения по [10] равна: с₁ = 2,5 л/ м³

Концентрацию СО₂ в наружнем воздухе примем равной по [2]:

с₂ = 0,3 л/ м³

Подставим значения в формулу (20):

L _СО2 = = 3409 м³/ч

По формуле (18) определим расчетный воздухообмен здания:

L = 3409 + 13428 = 16837 м³/ч.

4. Определение мощности системы теплоснабжения

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

N = 1,1?Q = 1,117344 = 19079 Вт

Расчет теплопотерь зданием по укрупненным показателям

Для ориентировочных расчетов потребности в тепле на отопление здания в целом используют удельную тепловую характеристику здания q_о. Зная q_о, можно определить количество тепла, необходимое для компенсации теплопотерь зданием, , не делая подробного расчета, по формуле [8]:

= q_оV (t_вн - t_н)? (21)

Где V - объем всего здания; ? - поправочный коэффициент к средней разности температур.

Средняя высота здания h = 3, 5 м, длина - а= 18 и 10 м, ширина - b = 10 и 10 м [7], значит объем здания V найдем по следующей формуле:

V = 3, 5 (1810+1010) = 980 м³.

Удельную тепловую характеристику здания q_о определим по табличным данным [1]:

q_о = 0,43 Вт/м³К.

Поправочный коэффициент ? определим также по табличным данным [9] в зависимости от средней температуры наиболее холодной пятидневки методом экстраполирования:

? = 1,17 + (31 - 29)/(34 - 31) (1,17 - 1,08) = 1,23.

Подставим полученные данные в формулу (21) и найдем количество тепла, необходимое для компенсации теплопотерь зданием, :

= 0,43980 (18 - (-31))1,23 = 16086 Вт.

Ориентировочные расчеты потребности в тепле на отопление помещений дома в целом оказались ниже вышеприведенного полного расчета здания на потребности отопления. Это говорит о том, что характеристика не учитывает такие теплоты, как потери через наружные стенки, через потолочное покрытие, через полы. Рассчитаем процент ошибки ориентировочных расчетов необходимой теплоты :

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

Задача гидравлического расчета состоит в подборе диаметров труб на отдельных участках циркуляционного кольца таким образом, чтобы по ним протекали расчетные расходы теплоносителя со скоростями, близкими к оптимальным (0,3 - 1,5 м/c). За главное циркуляционное кольцо примем всю систему отопления.

Расход воды на участке:

где: с - теплоемкость воды, равная 4200 Дж/(кг·^оС);

Скорость воды в системе отопления зададим равной v = 0,6 м/с. Определяем диаметр трубопровода:

где ? - плотность воды, кг/м³, при t = 72,5 ^oC.

Ближайшее стандартное условное проходное сечение трубы d = 0,032 м. Пересчитаем скорость воды в принятом диаметре трубы:



Определяем потери давления на прямолинейном участке:

где: k_ш - абсолютная шероховатость труб теплопровода, равная для водяных тепловых сетей k_ш = 5*10^-4 м [3]; l - дина подающего трубопровода расчетного участка, м.

Определяем падение давления в местных сопротивлениях:

где ? - коэффициент местных сопротивлений.

На данном участке около узла ввода поставим двадцать шаровых кранов с коэффициентом местного сопротивления ?₁ = 0,2 и 8 углов-поворота на 90? с коэффициентом местного сопротивления ?₁ = 0,21 [1].

По формуле Дарси - Вейсбаха [3] определяем циркуляционное давление:

?р_ц = ?р_л + ?р_м; (39)

?р_ц = 10068 + 400,3 = 10468,3 Па.

Определяем давление от естественной конвекции:

?р_ск = h · g · (?_х - ?_г) (40)

?р_ск = 2 · 9,8 · (981,1 - 971,8) = 175 Па,

где: h - высота по вертикали от середины нагревательного прибора до сборной магистрали; g - ускорение свободного падения, 9,8 м/с²; ?_х и ?_г - плотность воды на выходе и входе из циркуляционного кольца соответственно, кг/м³.

Так как ?р_ск < ?р_ц, то потребуется насос, который будет обеспечивать перепад давления не менее чем:

?р_н = ?р_ц - ?р_ск; (41)

?р_н = 10468,3 - 175 = 10293,3 Па.

Принимаем циркуляционный насос для систем отопления Grundfos серия UPS 32-80 180, со следующими характеристиками:

- максимальный напор (Н), м: 7,5;

- максимальная величина протока (Q), м³/ч: 1,4;

- установочная длина, мм: 180;

- мощность min, Вт: 135;

- мощность max, Вт: 220;

- максимальное давление в системе, бар: 10.

Таблица 4. Результаты гидравлического расчета теплопроводов системы водяного отопления

Тепловая нагрузка Nс.о., Вт	Длина трубопровода l, м	Расход теплоносителя G, кг/ч	Диаметр d, мм	Скорость движения теплоносителя ?, м/с	Потеря давления на трение ?рл, Па	Суммарные потери давления ?рц, Па
19077,3	72,27	0,3	32	0,38		10468,3

Заключение

теплозащитный стена трубопровод конденсация

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

Расчетное термическое сопротивление заданных наружных ограждений здания для данных климатических условий были ниже требуемых термических сопротивлений. Вследствие чего, была увеличена толщина стенки за счет утеплителя. За счет утепления наружных ограждений удалось сократить потери через наружные ограждения здания до минимальных значений.

Список литературы

1. Тихомиров К.В. Теплотехника, теплогазоснабжение и вентиляция/ К.В. Тихомиров. - М.: Стройиздат, 1981. - 272 с.

2. Захаров А.А. Применение тепла в сельском хозяйстве/ А.А. Захаров. - М.: Колос, 1974. - 311 с.

3. Теплоэнергетическое установки и системы сельского хозяйства/ Р.А. Амерханов, А.С. Бессараб, Б.Х. Дрананов, С.П. Рудобашта, Г.Г. - М.: Колос-Пресс, 2002. - 424 с.

4. Теплоснабжение, отопление и вентиляция: Метод. указ. к курсовой работе по дисциплине «Основы теплотехники и применения тепла в сельском хозяйстве»/ В.И. Ляшков, С.Н. Кузьмин. - Тамбов.: Тамбовск. ин-т хим. машиностр., 1993. - 48 с.

5. Солнечная водонагревательная установка: Метод. разр./ Сост.: В.И. Ляшков, С.Н. Кузьмин. - Тамбов.: Тамб. гос. тех. ун-т, 2004. - 20 с.

6. Отопление и вентиляция промышленных и гражданских зданий: Метод. указ. / Сост.: В.И. Быченок. - Тамбов.: ТИХМ, 1987. - 28 с.

7. Отопление и вентиляция промышленных и гражданских зданий: Метод. указ. / Сост.: В.И. Быченок. - Тамбов.: Тамбовск. гос. техн. ун-т., 1995. - 50 с.

8. Расчеты теплоснабжения, отопления и вентиляции: Метод. разраб. / Сост.: В.И. Ляшков, В.И. Быченок. - Тамбов.: ТГТУ, 1998. - 44 с.

9. СНиП 23-01-99 Строительная климотология/ М.: Госстрой России, 2000. - 67 с.

10. МДК 4-05.2004 «Методика определения потребности в топливе, электрической энергии и воде при производстве и передаче тепловой энергии и теплоносителей в системах коммунального теплоснабжения».

11. Краснощеков Е.А. Задачник по теплопередаче: Е.А. Краснощеков, А.С. Сукомел. - М.: Энергия, 1980. - 288 с.

Размещено на Allbest.ru