Исследование процесса терморадиационной сушки тканей

ИССЛЕДОВАНИЕ ПРОЦЕССА ТЕРМОРАДИАЦИОННОЙ СУШКИ ТКАНЕЙ

А.И. Ольшанский, В.И. Ольшанский

ВВЕДЕНИЕ

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

Сушка тканей термоизлучением является одним из самых перспективных методов термического обезвоживания тканей. Наряду с интенсификацией процесса сушки за счет возникающих больших удельных тепловых потоков получают значительный технологический и экономический эффект.

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

Механизм терморадиационной сушки тканей принципиально не отличается от кинетики конвективной сушки.

ПОСТАНОВКА ЗАДАЧИ

Сушка тканей термоизлучением проводилась на металлической поверхности стола экспериментальной установки при вынужденном и свободном теплообмене.

Ткань располагалась параллельно близко расположенным (200 мм) излучающим светлым источникам излучения. В качестве излучателей использовались лампы с нихромовой спиралью со степенью черноты окисленного нихрома е» 0,75- 0,8 и стекла е» 0,89 - 0,92. При принятом расположении ткани от излучающих поверхностей угловые коэффициенты излучения р, определяющие эффективность поверхности излучения, р»1. Степень черноты натуральных тканей е» 0,75 — 0,9 .

Экспериментальное исследование комбинированной конвективно-терморадиационной сушки тканей проводилось с пятикратной повторяемостью эксперимента для каждого режима сушки с целью получения достоверных результатов.

Сушка тканей проводилась при трех режимах:

1. = 150°С; tc = 80°C; скорость воздуха о» 0,9 —1,2 м/с;

2. tu3]i = 190 °С; tc = 100 0C; о » 0,9 -1,2 м/с;

3. t^ = 250 ос; tc = 120 ос; о = 0 м/с (свободный теплообмен).

На рис. 1 и 2 изображены кривые сушки W = f

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

Рисунок 1 - Кривые сушки

Ж = /

' ' и температурные кривые 4 ' при конвективно-

терморадиационной сушке ткани из вискозы для режимов: 1 — іс = 80 °С; іизл = 150 °С;

2 — і = 100 °С; ітл = 190 °С; у= 1,2 м/с; 3 — іс = 120 °С; ітл = 250 °С; у= 0 м/с

изл

О/

й/%

160

120

80

/Ієн

1  2 2 3

/ /

У / 1 / 1

/

/

/Л

80

60

1*0

20

О

О 0,5

1

1.5

2

3,5

1+ гмин

п о „ Ж = /

Рисунок 2 - Кривые сушки ' ' и температурные кривые 4 ' при конвективно-

терморадиационной сушке ткани изо льна. Режимы сушки указаны на рис. 1

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

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

Сушка тканей термоизлучением при принятом способе энергоподвода представляет крайне сложный комбинированный терморадиационно-конвективно-кондуктивный процесс сушки с переменной температурой на металлической поверхности экспериментального стола.

Исследованиями установлено [1], что максимальная глубина проникновения инфракрасных лучей зависит от свойств материала. Для натуральных тканей тепловые лучи проникают на толщину более четырех слоев. В большинстве случаев влажность материала снижает глубину проникновения инфракрасного излучения только в периоде постоянной скорости сушки и поглощение термоизлучения происходит поверхностью материала [1].

Для тонких тканей (S = 0,7 ^ 1,2 мм) в результате проникновения излучения через

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

Механизм комбинированной сушки тканей термоизлучением - сложный процесс тепломассообмена и строгая математическая постановка задачи на нахождение полей влагосодержаний и температуры на основе решения дифференциальных уравнений влаготеплопереноса представляет большие трудности. Поэтому задача исследования ограничивается только изучением внешнего механизма терморадиационной сушки тканей и кинетикой процесса на основе обработки экспериментальных данных.

Обработка экспериментальных данных. При исследовании влаготеплопереноса при сушке тканей термоизлучением воспользуемся результатами аналитического решения задачи на нагревание влажной неограниченной пластины в среде с постоянной температурой tc = const А. В. Лыкова [2] и методом регулярного режима Г. М. Кондратьева [2, 3].

Основным признаком регулярного режима нагревания влажного тела при изменении температуры и влагосодержания являются соотношения

dt „

~^- = mt(tC —t);

dr

dW

_ _=m W _ WP),

где dt /dz и dW/dr - скорости изменения температуры и влагосодержания во времени; tc - температура среды; t и W - среднеинтегральные температура и влагосодержание материала; Wp - равновесное влагосодержание материала; m и m - темп нагревания влажного тела и темп убыли влагосодержания, (мин-1).

Темп нагревания влажного тела mt и темп убыли влагосодержания mu находятся из соотношений на основе эксперимента:

ln

tgpj = —L-C-----—----L-C--- = m = const;

r — r0

ln

T_To U

где ^ - среднеобъемная температура материала в первом периоде сушки; Wkp - критическое влагосодержание материала, определяемое по кривым сушки; t и W- значения температуры и влагосодержания, соответствующие текущему времени сушки r.

Отсчет времени сушки в периоде падающей скорости для всех режимов проводится от т = 0 и влагосодержания Wkp Среднеобъемная температура /0 в первом периоде вычисляется по соотношению [4]

^ = 10 + 0,29^.

На рис. 3 а и рис. 4 изображены зависимости 1п

для вискозной и льняной тканей при терморадиационно-конвективной сушке.

Видно, что в стадии регулярного режима эти графики имеют вид прямых.

Значения ш( и ши для некоторых влажных материалов даны в таблице 1.

Таблица 1 - Значения постоянных щ и ши для некоторых влажных материалов

8, мм W , кр ’ % Режим сушки т, -1 мин т, -1 мин

Материал г , °С с ’ и, м/с

90-150

Войлок 8-10 75-80 50-60 3-10 0,02 0,08

Фетр Подошвен-ная 4 35-40 40-60 0,5-0,8 0,02 0,08

кожа Вискозная и 4 57-60 80-120 3-5 0,04 0,055

льняная ткань 0,7-1,2 60-80 і им=150-250 0-1,2 0,11 0,45

Рисунок 3 а, б - Зависимость

1п

при комбинированной сушке

Л - г — 80оС г — 150°С вискозной и льняной ткани термоизлучением при режимах: а) с ; изли — 1,2 м/с 2 - г — 100°С г — 190°С и — 1,2 м/с 3 - г — 120°С г — 250°С

; ’ ; с ; изл. ; ’ ; с ; изл. ;

І I 15 5 9 9

и — 0 м/с. | - вискоза, II - лен;

б) зависимость / (г) / (^1п ) для комбинированной сушки вискозной и льняной

ткани термоизлучением при режимах, указанных на рис. 3 а

ш-ш

З

2

1

хо* II 1 - /

Д д  - //

АД •

• 4^9

тмин

Рисунок 4 - Зависимость

и льняной ткани термоизлучением при режимах: I - вискоза: о = 1,2 м/с 2 - г = 100оС. * = 190° С. о = 1,2 м/с. 3 - с = 120°С. *

С  иЗЛ. -7-е  изл

^ = о м/с .. 4 - ї = 80оС ї -1^°^ >'-п„/р * /-шло

II - лен: С - изл

ї = 190ос ^ = 1,2м/с 6-ї = 120ОС. ї = 250ОС. ^ = 0м/с

при комбинированной сушке вискозной 1 -^ = 80ОС. = 150ОС.

250О С

150оС V = 1,2м/с 5-ї = 100ОС.

Из решения уравнений (1), (2) с учетом выражений (3), (4) находится время сушки т и температура материала * в периоде убывающей скорости [4]:

1

ш.

ти =—Ы

гі - г л

1С 10

г — г

V гс г

1

тп =

1п

С учетом времени сушки в первом периоде т = РТв — №кр/М - продолжительность сушки тканей до заданного влагосодержания равна:

Уо — Укр 1 ,,

т =-----——- +-------1п

— У ^

гг кр р

N

ш..

V

У — У

р У

где N = dW / ит - максимальная скорость сушки в первом периоде;

W - текущее влагосодержание материала.

Из уравнения (6) находится среднеинтегральная температура в процессе сушки тканей:

г = г.

г — г

гс г0

ЄХр

где т - время сушки от начального Ж0 до текущего значения влагосодержания Ж При исследовании теплообмена в процессе сушки воспользуемся результатами аналитического решения дифференциального уравнения на нагревание влажной

неограниченной пластины в среде с постоянной температурой tc = const при граничных условиях третьего рода (сушка) А. В. Лыкова [2, 5].

Из решения задачи можно записать:

t — t 0 _ lC 1П

tC t о.с.

= f (exp Fo - Pd) = f (exp Pt),

где 1П, Ьдс - температуры поверхности влажного материала и окружающей среды; ¥о

- критерий Фурье; Рй - критерий Предводителева;

Произведение критериев Фурье ¥о и Предводителева Рй представляет соотношение [4, 5]

17 па ах а

Жо • Рй = —- •—----= Рт,

Я2 а

где Р - скорость изменения температуры на поверхности влажного тела, мин-1; 6 - относительная температура влажного тела, м; а - коэффициент

температуропроводности, м2/с; Я - характерный размер тела, м; т - время сушки, мин.

На рис. 5 представлена зависимость \%в = / (г) при сушке вискозы и льна.

о

0,2

ОА

0,6 -1р 0

1

2

3 4

7

тмин

ОІ N. д 1 • ОХ - /

3  II 1 <ча - //

і >

с

Рисунок 5 - Зависимость относительной температуры ^0 = / (г) для

комбинированной сушки вискозной и льняной ткани термоизлучением при режимах,

указанных на рис. 1 I - вискоза, II - лен

Обработкой зависимости ^0 = /(г) (рис. 5) получено значение коэффициента р при

сушке тканей р= 0,25 мин-1.

Из выражения (10) определяется температура тканей в процессе сушки:

* = *п = *с -

Изменение температуры тканей в процессе сушки можно также выразить и в виде степенной зависимости ^ = f (рт)’“.

На рис. 3 б дана зависимость = f(lgPт) в процессе комбинированной сушки

тканей термоизлучением при различных режимах.

Обработкой экспериментальных данных, представленных на рис. 3 (б) методом

средних, получено уравнение для вычисления температуры тканей в процессе

комбинированной сушки:

іп = В,

В = 0,7іі

С целью проверки достоверности полученных результатов обработки экспериментальных данных проведено сопоставление расчетных значений времени сушки и температуры тканей по уравнениям (8), (9), (11) и (12) с экспериментом, представленное в таблице 2. Видно, что расчетные значения параметров, определяемые по формулам, хорошо согласуются с экспериментом.

Таблица 2 - Сравнение расчетных значений времени сушки и температуры тканей по уравнениям (8), (9), (11) и (12) с экспериментом при комбинированной сушке___

Вискоза: ‘изл —150 °С; ‘с — 80 0С; о —1,2 м/с

W, % Тжс,мин т, мин

80 0,51 0,53 35 36,3 30,5 34,6

60 1,65 1,58 41,2 40,7 41 41,5

40 2,98 2,88 45,5 46,2 46,5 46,9

20 4,69 4,36 53,5 52,5 54,9 53,5

10 6 5,92 58,5 56,7 59,5 56,9

Вискоза: ‘изл —190 °С; ‘с —100 °С; о —1,2 м/с

80 0,75 0,70 42,5 42,5 39 41,8

60 1,35 1,25 51 49,8 48,7 52,5

40 2,25 2,1 58,5 57,8 57,5 60

20 3,85 3,78 67 66,6 68,9 67,5

10 5 5,1 70 68,5 73 71

Лен: ‘изл —150 ОС; ‘с — 80 °С; о —1,2 м/с

120 0,8 0,8 38 37,5 36,1 38,6

80 1,65 1,75 40 41,5 40 43,5

40 2,65 2,58 42,5 43,5 44,2 45,6

10 4,15 4,6 53 51,5 55,1 54,4

Лен: ‘изл —190 °С; ‘С —100 °С; о —1,2 м/с

100 0,65 0,6 44 43,5 40 44,4

80 0,75 0,76 47 45,5 46,5 47,6

60 1,35 1,31 50 47,5 48,5 52,5

40 1,55 1,95 52 49,7 53 49,8

10 4,25 4,7 71 70 73,5 72,5

Из решения этой же задачи А. В. Лыковым [2] получен также критерий, который представляет отношение психрометрической разности (‘с-‘м) к избыточной

температуре окружающей среды (‘с — ‘ос) :

‘ —ь

/з с ‘ м .

°М --------—»

‘С ‘о. с.

Ом - ,

акр

где ‘м - температура мокрого термометра, равная температуре поверхности материала в периоде постоянной скорости сушки; }0 - интенсивность испарения влаги в первом периоде, кг/м2ч; а - коэффициент теплообмена в первом периоде, Вт/м2-°С; г -

теплота парообразования, кДж/кг; 0М - относительная температура в первом периоде сушки.

Теплота парообразования вычисляется по формуле

Я = 2490 - 2,3 ‘

Из уравнения (14) и (15) имеем:

і = і — -п±-

1М 1Є —

а

кр

Таким образом, относительная температура 0М приобретает значение

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

Уравнения (15) и (17), установленные А. В. Лыковым [2] для периода постоянной скорости сушки, используем в расчетах периода падающей скорости.

Запишем уравнения (15) и (17) в виде

л _ ІС І П &М ~

По

ІС - 1о.е. а

VI

і = і ——

1П 1С —

а

где $ и а - интенсивность испарения влаги кг/м2ч и коэффициент теплообмена в периоде падающей скорости, Вт/м2 °С.

Коэффициенты теплообмена а при комбинированной сушке тканей термоизлучением определялись по формуле [6, 7] для вынужденного и свободного теплообмена:

N4 = 0,87 Яв

0,5

( т \ Тс 2 г т ^ изм 0,45 Ґ \ Ж 0,5

Т V ТМ ) Т Ж V кр

N4 = 0,75^г • Рг)

0,25

^ т > ТС 2 (Т Л изм 0,45 Г Ж Ї

Т V ТМ Т V 1 С V Жкр)

0,5

где Ми - критерий Нуссельта; Яв - критерий Рейнольдса; О - критерий Грасгофа; Рг -критерий Прандтля; Тс , Тм, Ттл - абсолютные температуры среды мокрого термометра, излучателя.

Поправка (Ж/Ж )0,5учитывает уменьшение коэффициента теплообмена а с

кр

уменьшением влагосодержания. Интенсивность испарения влаги ] определялась в экспериментах весовым способом и по кривым сушки

] =

,кг/м2ч,

где dW/dт - скорость сушки во втором периоде; О0 - абсолютно сухой вес материала, кг; F - поверхность испарения ткани, м2.

Анализ экспериментальных данных по комбинированной конвективно-радиационно-кондуктивной сушке ткани показал, что температура металлической поверхности стола оказывалась выше температуры ткани на 5 - 12 °С (рис. 1) за счет максимальной глубины проникновения излучения через ткань и возникновения кондуктивного дополнительного подвода тепла от металлической поверхности к ткани. При этом значительно возрастает интенсивность испарения влаги за счет совпадения градиентов влагосодержания и температуры по направлению к поверхности ткани.

Обработка опытных данных показала, что определяющей температурой в уравнениях (18) и (19) при сушке термоизлучением является не температура среды *с, а

температура излучателя /изл . Уравнение (19) примет вид:

ґ п — ґ — ґ —

П изл

V.

а

Для тонких тканей среднеобъемная температура ? и поверхности ^ практически одинаковы.

Рассмотрим влияние критерия Фурье Го на распределение температуры во влажном теле О . А. В. Лыков [8] отмечает, что если распределение температуры в телах подобны между собой, то для них величина ат/Я2 должна иметь одинаковое значение. Поэтому при заданном начальном распределении температуры последующие распределения зависят от сочетания величин ат/Я2, которые образуют обобщенную переменную. Критерий Фурье Го приобретает смысл обобщенного времени. Следовательно, распределение температуры для одного вида материала с одинаковым характерным размером Я = 5 является подобным и величина критерия Фурье Го должна иметь одинаковое значение для различных режимов сушки при одинаковых значениях времени т.

На рис. 6 представлена зависимость относительной температуры О = /(Го) в процессе комбинированной сушки тканей термоизлучением при различных режимах.

0

А • • ?_ А П ' /

А 9 • • - • ♦ А  //

.4  А Ж а _ X # • *

А А

0,8

0,6

ОА

О

5

10

15 20 25

30 35 /?;

о

Рисунок 6 - Зависимость относительной температуры 0 от критерия Фурье ¥о при

комбинированной сушке вискозной и льняной ткани термоизлучением. Режимы сушки

указаны на рис. 4

Обработка опытных данных показала, что для материала с одним характерным размером Я все опытные точки с допустимым разбросом укладываются на одну кривую (рис. 6). Зависимость 0 — /(Го) приближенно можно описать уравнением

0 — /

—0,18

Для вычисления критерия Фурье Го необходимо знать закономерности изменения коэффициента температуропроводности а от температуры и влажности. Коэффициент температуропроводности тканей

Я

а —

ср

где Я, с, р - соответственно теплопроводность, Вт/м2 -°С; теплоемкость, Дж/кг -°С; плотность влажной ткани, кг/м3.

Зависимость коэффициента теплопроводности Я тканей от температуры в границах от t«20—70°С незначительна и при расчетах не учитывается [9].

Влияние влажности на теплопроводность текстильных тканей хорошо исследовано [9] и выражается формулой

Я = Я0+ 0,00144W, Вт/мград,

где Я0 - коэффициент теплопроводности тканей в сухом состоянии;

W - влажность тканей на сухой вес.

Теплоемкость влажных тканей определяется выражением [6]

W

с + с ----

0 в 100 с = 100 > Дж/кг -°С,

1 +--

100

где со - теплоемкость сухой ткани; се = 4186 кДж/ кг°С - теплоемкость воды.

Плотность р влажной ткани в процессе сушки измерялась в экспериментах весовым способом.

Вычисленные значения коэффициента температуропроводности а приведены в таблице 3.

Таблица 3 - Зависимость коэффициента температуропроводности для тканей от влажности W

Вискоза Я = 0,058 Вт ; со = 1600Дж / кг • °С 0 м °С Лен Я„ = 0,046 Вт ; со = 1500Дж / кг • °С 0 м °С

W, % Вт кг Р> — м с> Дж / кг '°С а • 104 м /ч W, % , Вт кг Р — м с> Дж / кг -°С а • 104 2 м 2 /ч

Я, м °С м °С

90 0,187 759 2825 3,14 160 0,276 675 3150 4,67

80 0,173 713 2750 3,17 140 0,248 637 3060 4,58

60 0,144 640 2569 3,15 120 0,219 587 2960 4,537

50 0,13 600 2462 3,16 100 0,19 537 2840 4,49

40 0,116 556 2339 3,24 80 0,161 482 2690 4,48

20 0,087 488 2031 3,16 60 0,132 430 2500 4,47

10 0,072 453 1835 3,11 40 0,104 374 2260 4,41

20 0,075 321 1946 4,32

10 0,06 294 1744 4,27

При исследовании теплообмена в процессе комбинированной сушки тканей термоизлучением использовались критериальные уравнения для вынужденного и свободного теплообмена (20), (21) с поправками на излучение.

Плотность потока тепла в периоде убывающей скорости сушки находится из уравнения теплообмена

Ц = &

Для проверки достоверности полученных при исследовании результатов проведен расчет теплообмена по формуле излучения Стефана-Больцмана [6, 7]

а = е С

ч. изл пр о

^ ^ N 4 4 " ( \ ж

— ф

1100 J 110(0 J ж К кР

0,5

где Со - коэффициент излучения абсолютного черного тела, Со = 5,67 Вт/м2К4; ф -угловой коэффициент излучения ф = 1; ^/ Wкp)0’5 - поправка на уменьшение тепловых потоков с уменьшением влагосодержания.

Приведенная степень черноты системы определяется выражением [6, 7]

1

епр =

еГ

+ф

-1

\е2

где е1 - степень черноты излучателя (нихром) е = 0,80,81; е2 - степень черноты натуральных тканей е = 0,75 — 0,9.

Результаты сопоставления расчетов по уравнениям (20), (22), (23) (25) и (26) с экспериментом представлены в таблице 4.

Таблица 4 - Сопоставление расчетов по уравнениям (20), (22), (23) (25) и (26) с экспериментом

W, % а , Вт/м2 °С

Вискоза: = 150 0С,ЬС = 80 0 С,и = 1,2 м/с , ен/, = 0,845

80 31,6 1,42 42 40 1264 1228

60 29,7 1,33 43 44 1069 1050

40 24,5 1,07 46 46 833 839

20 21 0,85 55 54 546 572

10 17 0,65 59 58 374 383

Лен: Ьшл = 190 0 С, Ьс = 100 0 С, и = 1,2 м/с , ещ = 0,86

120 38,4 2,36 43 41,5 2166 2086

100 35,1 2,15 45 46 1895 1785

80 31,9 1,9 47 47 1691 1562

60 28,9 1,72 48,5 49 1474 1369

40 25,1 1,44 53 51,5 1217 1105

80 22,2 1,21 61 60 888 759

10 17,6 0,92 68,5 67 581 539

1

1

Сравнение расчетных значений температуры # в процессе сушки тканей с экспериментом и тепловых потоков по двум разным методам является подтверждением достоверности результатов, полученных при исследовании кинетики процесса сушки ткани термоизлучением.

ВЫВОДЫ

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

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

Список использованных источников

1. Лыков, А. В. Теория сушки / А. В. Лыков. - Москва : Энергия, 1968. - 470 с.

2. Лыков, А. В. Теория теплопроводности / А. В. Лыков. - Москва : Высшая школа, 1967. - 600 с.

3. Лыков, А. В. Тепломассообмен (справочник) / А. В. Лыков. - Москва : Энергия, 1971.

- 651 с.

4. Ольшанский, А. И. Исследование процесса сушки тонких материалов методом регулярного режима / А. И. Ольшанский, В. И. Ольшанский. // Известия НАН Беларуси. Серия физико-технических наук. - 2011. - № 2. - С. 75-81.

5. Пехович, А. И. Расчеты теплового режима твердых тел / А. И. Пехович, В. М. Жидких. - Ленинград : Энергия, 1968. - 303 с.

6. Теплоиспользующие установки в текстильной промышленности : учебное издание для вузов / под общ. ред. Е. А. Ганина. - Москва : Легпромиздат, 1989. - 390 с.

7. Исаченко, В. П. Теплопередача / В. П. Исаченко, В. А. Осипова, А. С. Сукомел. -Москва : Энергия, 1969. - 439 с.

8. Лыков, А. В. Теория тепло- и массопереноса / А. В. Лыков, Ю. А. Михайлов. -Москва. - Ленинград : Госэнергоиздат, 1963. - 535 с.

9. Колесников, П. А. Теплозащитные свойства одежды / П. А. Колесников. - Москва : Легкая индустрия, 1965. - 345 с.

Статья поступила в редакцию 12.12.2011