Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Календерьян В.А. Теплоотдача плотного движущегося слоя и методы ее интенсификации

.pdf
Скачиваний:
2
Добавлен:
23.10.2023
Размер:
7.9 Mб
Скачать

Обобщение данных по теплоотдаче материалов с хорошими сы­ пучими свойствами приведено на рис. IV.4. Все данные, полученные для смесей, лежат в первой области, данные для фракций охватыва­ ют обе области.

Ни

 

 

 

 

 

 

 

400

600 6001000 1500 Ре

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

о - )

4 - 8

<-- 15

 

 

 

 

 

 

 

 

 

 

° - 2

» - 5

*•

/6

 

 

 

 

 

 

 

 

 

 

ѵ-з

»-m

t

17

 

 

 

 

 

 

 

 

 

4-4

« - / /

 

 

 

 

 

 

 

 

 

 

 

 

"-5

»-12

 

 

 

 

 

 

+

 

i

 

 

 

+ - «

- / 3

 

 

 

 

ï

 

 

r 1 и 1

i

u

 

- 7

 

 

 

 

 

S i

 

 

 

 

 

 

r

 

T

? 0 à о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

—t—' DO *

 

 

 

 

 

 

 

 

 

 

 

 

0,8

 

 

 

 

 

 

 

 

 

 

 

 

Ou

0,6

 

 

 

 

 

 

 

 

 

 

 

 

fo

го

 

40

SO 80 100

§

200

400

600 8001,

 

2000

Ре

Р и с . I V . 4. Обобщенны е зависимост и по теплообмен у слоя

с горизонталь ­

 

 

 

 

 

 

 

ным и

цилиндрами :

 

 

 

 

а — с м е с и :

сГ=0,52 мм;

1 — D-21,5;

2 — D=33,5;

3 — D=42,5; 4 — D=48;

5 — 0 - 57 ;

6 — 0 - 76 . 5

л я ; rf=l,2D

« л ; 7 — .0=21,5 мм; б — фракционированные

слои

(усл. обозн,

 

 

см. рис. IV.1 и ІѴ.З): / — Р е < Р е п ;

/ / — Р е > Р е п .

 

 

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

а) для смесей (песок, агломерат, полукокс, электрокорунд) при

170 <

Ре <

Ре п , 47 <

 

< 238

 

 

 

 

 

№Ï

=

0 , 8 P e ° ' 2 8 g ) ° ' 3 3 ;

(IV. 2)

б) для фракционированных слоев (песок,

карбид кремния) при

10 <

Ре <; Ре п ,

13 <

 

< 147

 

 

 

 

 

Ш

=

0 , 3 9 Р е ° . 2 в ( £ ) М 8 ;

(IV. 3)

пр и Ре„ <

Ре

< 3200, 13 <

<

147

 

 

 

 

Nu

=

0,84

Р е 0 ' 0 6 8

( f f 6 8 .

(IV. За)

6*

 

 

 

 

 

 

 

83

Определяющими параметрами являются диаметр цилиндра, скорость слоя в минимальном сечении и его средняя температура. Средний размер частиц смесей определен по уравнению (П.7), что

обеспечило наилучшее обобщение опытных данных.

 

Уравнения (IV.2) и (IV.3) с вероятной ошибкой ± 6 %

справедли­

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

 

Для получения единой зависимости по материалам

с плохими

сыпучими свойствами имеющихся данных недостаточно. В связи с

этим для них

рекомендуем

предварительно

частную

формулу

 

Nu

=

с Ре",

 

 

 

(IV. 2а)

где для концентрата с =

1, 2,

я

= 0,26; для

шлиф- и

микропорош­

ков карбида кремния с =

1,

5,

п =

0,26.

 

 

 

Уравнения

(IV.2) и (IV.3) могут быть использованы

при

нагреве

и охлаждении слоя и уровне температур до 500—600° С.

 

Учитывая

проиллюстрированное

выше влияние времени

контак­

та на теплоотдачу, можно преобразовать уравнение (VI.2), введя в

него в качестве определяющего критерий Фурье,

 

 

N u = l , 3 4

F o - ° - 2 8 g ) 0

> 3

3 ,

(IV. 4)

где 3,4 -10

<

Fo = — - — = — <

5,

7-10

. К аналогичному

 

 

D2

 

Ре

 

 

 

виду могут

быть приведены

уравнения (IV.3).

 

Представляет

интерес

сравнение

экспериментальных зависи­

мостей типа (IV.2)—(IV.4), учитывающих реальные условия омывания цилиндра, с расчетными, полученными при рассмотрении слоя как неограниченного массива с цилиндрической полостью. Экспе­

риментальные зависимости Nu

= / (Fo) лежат ниже и выражены сла­

бее расчетной (для последней

Nu ^ F o ° ' 5 ) . Это объясняется тем,

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

но зависит от времени контакта. Отсюда вытекает

недопустимость

использования для поперечного омывания каких-либо

представлений

и расчетных зависимостей, игнорирующих действительный

характер

обтекания цилиндра (как это сделано, например, в [215]).

Сопоставим полученные результаты с данными

[57,

130]. На

рис. IV.5 приведены результаты расчетов по (IV.2), (ІѴ.З) и крите­

риальным зависимостям [57, 130]. Расчеты выполнены для

конкрет­

ного материала — смеси, фракционный состав которой приведен в

табл. I I . 1, причем

размер частиц усреднялся по

 

рекомендуемой

авторами

методике

(см. гл. I I ) . Наши

данные для

фракционирован­

ного слоя и смеси удовлетворительно

согласуются

с уравнением

С. В. Донскова

[57]: темп зависимости

от критерия

Пекле, влияние

 

D

 

 

 

 

 

симплекса

-^-

примерно одинаковы,

количественные

расхождения

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

84

нении с данными 11301 для сухих материалов: согласно ИЗО] влия­ ние на теплоотдачу критерия Пекле меньше, а симплекса-^- больше,

чем по уравнению (IV.2), влияние же диаметра цилиндра противо­ положно. В связи с этим количественные расхождения зависят от этих параметров и для принятых условий достигают 100% (при этом зависимости Ю. П. Курочкина систематически ниже наших).

 

 

 

 

fL

 

 

 

 

12

Р и с .

I V . 5. С р а в н е н и е с

л и т е р а т у р н ы м и дан -

8

 

нымп:

 

 

/ — по

уравнению (IV.2);

/ / — по

уравнению

^

(ІѴ.З); / / / — п о С. В. Донскову; IV— по

ІО. П. Ку-

 

 

рочкнну.

 

 

В целом эти выводы остаются справедливыми независимо от фракци­ онного состава материала. Указанные обстоятельства не позволяют рекомендовать приведенную в [1301 зависимость для расчетов, тем более, что она противоречит прямым опытным данным для мелко­ дисперсных материалов ^при ~ > 2о) о влиянии диаметра цилинд­ ра на теплоотдачу.

IV. 2. НАКЛОННЫЕ ЦИЛИНДРЫ

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

та теплообменников с различной компоновкой поверхности

нагрева.

Изучался средний теплообмен кварцевого песка (смесь d =

0,52jnjn)

с неоребренными цилиндрами диаметром 33,5; 48,2

и 57,5

мм,

длиной

280

мм при углах атаки г[) = 90° (горизонтальное

положение), 85,

75,

60, 30°, изменении скоростей слоя от 0,75 до 5 мм/сек,

симплекса

от 70 до 120. Температуру поверхности измеряли 12 медь-кон-

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

Анализ первичных экспериментальных данных показывает, что

85

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

'S

го

50

40

50

60

70

80

90І00В/І

 

 

 

 

 

D

 

 

 

Рис .

I V .

6. В л и я н и е

симплекса ~j

на

теплообмен

слоя

с н а к л о н н ы м и

ц и л и н д р а м и

при

ф,

р а в н о м :

 

 

/ — 90°; 2 — 83°;

3 — 75°;

4 — 60°;

5 — 30°.

 

Обобщение экспериментальных данных производилось в такой

же форме, как для горизонтальных цилиндров, но с учетом угла

ата­

ки.

 

На рис. IV.6 приведены частные зависимости Nu = f j-^j

при

Ре = idem, причем г|) принят в качестве параметра. Увеличение сим­ плекса интенсифицирует теплообмен, но с уменьшением угла ата­ ки это влияние сглаживается и при *ф = 30° практически исчезает. Аналогичные результаты получены при различных значениях кри­ терия Пекле.

Зависимость интенсивности теплообмена от критерия Пекле оп­ ределяется углом атаки и с уменьшением последнего становится бо­ лее существенной (рис. IV.7).

Все зависимости, приведенные на рис. ІѴ.7, описываются кри­ териальным уравнением вида

NÜ = cPen g)m , (IV. 5)

где с, п, m — переменные величины, зависящие от тр.

При больших углах атаки структура уравнения (ІѴ.5) такова же, как (ІѴ.2) для поперечного омывания, при малых углах — анало­ гична приведенной выше для продольного омывания при нестеснен­ ном движении. Таким образом, изменение угла атаки приводит не только к количественным, но и к качественным различиям.

Описанные явления объясняются изменением характера движе­ ния и, следовательно, теплообмена при переходе от поперечного омы-

86

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

играет

 

 

 

D

учитываю-

 

 

 

 

 

 

симплекс — ,

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

щий изменения

 

в

характере

 

 

 

 

 

 

омывания. При

вертикальном

 

 

 

 

 

 

расположении

цилиндр

омы­

 

 

 

 

 

 

вается

безотрывно,

интенси­

 

 

 

 

 

 

фицирующее

действие скорос­

 

 

 

 

 

 

ти одинаково сказывается на

 

 

 

 

 

 

теплоотдаче всей

поверхности,

 

 

 

 

 

 

темп зависимости

от

критерия

 

 

 

600 800 im im im

pe

Пекле возрастает. При нестес-

 

 

 

 

 

 

„ *

„ _ т

 

 

' , „ „ „ „ , , „

Рис . I V .

7.

О б о б щ е н н ы е з а в и с и м о с т и

по

ненном

продольном движении т е п л о о б м

е н у

с л о я с

н а к л о н н ы м и

ци-

соотношение между диаметром

 

 

л и н д р а м и :

 

 

цилиндра и

размером

частиц

/ - д=зз . 5 :

2 — £>=48,2:

з — 0=57.5

мм.

 

на теплообмен

не влияет.

 

 

 

 

 

 

Д л я

наклонных цилиндров закономерности движения и теплооб­

мена носят промежуточный характер, включая особенности и по­ перечного и продольного омывания. При больших углах атаки пре­ обладают первые, при малых — вторые, что и отражает формула (IV.5). Ухудшение теплоотдачи при уменьшении угла атаки можно объяснить увеличением времени контакта слоя с поверхностью из-за того, что слой не только омывает наклонный цилиндр по периметру,

но и «сползает» вдоль образующей. В связи с этим наблюдается

сни­

жение интенсивности

теплообмена по длине, тем более значитель­

ное, чем ближе расположение цилиндра к вертикальному.

Д л я

иллюстрации на рис. IV.8 приведены распределения относительных

локальных коэффициентов теплоотдачи по длине при разных

углах

атаки и скорости слоя

1,3 ммIсек (значения ах при -ф =^90° отнесе­

ны к соответствующей величине для той же точки и той же скорости

слоя

при т|з =

90°).

 

 

 

 

 

 

Ухудшение теплоотдачи по периметру в верхнем сечении (IV)

невелико и

при

я|э = 83° составляет

примерно

3—4%,

а

при

=

= 30°—20%. По

мере удаления от

верхнего

сечения

наблюдается

дальнейшее

падение интенсивности

теплообмена, особенно резкое

при малых углах атаки. Так, в нижнем сечении (I)

при

і|з =

83

коэффициент

теплоотдачи на 6—8%

ниже,

а

при

-ф =

30° —

на

45%

ниже,

чем для горизонтального

цилиндра.

Д л я всех

характер­

ных точек цилиндра (лобовая, кормовая и экватор) темп изменения теплоотдачи по длине практически одинаков и, следовательно, вы­ зван тем же эффектом «сползания», который, очевидно, имеет место по всему периметру. Аналогичные результаты получены и при дру­ гих скоростях слоя, диаметрах цилиндров.

87

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

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

MJ 1

 

 

 

 

 

,

-

 

 

 

 

 

t

 

 

 

 

 

п

п

 

Р и с . I V . 8. Р а с п р е д е л е н и е к о э ф ф и ц и ­

Р и с .

I V . 9.

З а в и с и ­

ентов т е п л о о т д а ч и

по д л и н е

цилинд ­

мость

т, п, с

от

угла

ра:

 

 

а т а к и .

 

 

/ — ф - 8 3 ° ; / / — ф=30°;

/— ф=0°:

2— ф~90°;

 

 

 

 

3 _ ф = 1 8 0 " .

 

 

 

 

 

слоя с наклонными цилиндрами единой формулой (IV.5), где коэф­ фициент и показатели степени зависят от угла атаки и определяются по рис. IV.9. В качестве нижней границы применения формулы для надежности принято значение 30°, так как при меньших углах атаки, возможно, потребуется учитывать относительную длину цилиндров

№•

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

IV. 3. СФЕРИЧЕСКИЕ ПОВЕРХНОСТИ

Исследование локального и среднего теплообмена слоя со сферой было проведено в Институте технической теплофизики А Н УССР при участии одного из авторов [164]. Сфера диаметром 25 мм располага­ лась в вертикальном цилиндрическом канале диаметром 100 мм. Локальные коэффициенты теплоотдачи определяли с помощью дат-

88

чика удельного теплового потока 141) и рассчитывали по уравнению (II.4). Датчик размером 8 X 8 X 1,2 мм, представлявший собой ба­ тарею дифференциальных микротермопар с плотностью укладки око­ ло 2000 шт. /см2, располагался на поверхности латунной сферы, в которой размещался электронагреватель (рис. IV . 10). Для измере­ ния температурного напора между слоем и поверхностью служила медь-константановая дифтермопара, горячий спай которой зачеканен

Рис. IV .

10. Сферический

калориметр :

/ — термопара; 2

— пластина; 3 — датчик теплового потока;

4—изоляция;

5 — латунная сфера;

6 — нагреватель.

в латунной пластинке, покрывавшей датчик. Плотность теплового потока датчика определялась по тарировочной кривой q = f (Е). Ошибка измерения локальных коэффициентов теплоотдачи по дан­ ной методике составляет около 6%. Средние по поверхности сферы коэффициенты теплообмена подсчитывали как среднеинтегральные (уравнение (П.2)).

Сыпучим материалом служили частицы силикагеля размером

3—Ъмм(а — 4мм), характеристики которого приведены в табл. I I . 1 - Температура сферы в опытах не превышала 100° С, скорость слоя в минимальном сечении изменялась в пределах 3,2—43,2 мм/сек, определяющие геометрические симплексы оставались неизменными и

составляли -^- = 6,45, 9,35. При каждой скорости опреде­ ляли локальные коэффициенты теплообмена в 13 точках при поворо­

те сферы от 0 до 180° (через 15°).

 

 

 

 

На рис. IV . 11 приведены распределения а ф по поверхности,

полу­

ченные

при V = idem. Интенсивность

теплообмена

максимальна в

лобовой точке, а по мере удаления от нее монотонно

убывает. При

этом ухудшение

теплоотдачи

до ср =

70 90° незначительно (на­

пример, при V =

12 мм/сек

с 270 до 240 втІм2град),

а затем ста­

новится

очень существенным

(до 60

вт/м2град

при ф =

180°).

С увеличением скорости слоя коэффициенты теплоотдачи возрастают, и картина и распределения остается практически неизменной. Темп изменения коэффициента теплоотдачи на верхней половине сферы

89

практически одинаков (п « 0,23), в области 90° < ф < 165° он становится меньше (п » 0,13), а в кормовой точке ф = 180° снова возрастает (я » 0,30). Описанные выше закономерности позволяют судить о характере омывания сферы. Можно полагать, что ее верхняя половина омывается безотрывно, застойная зона в лобовой части от­ сутствует; в нижней части наблюдается отрыв слоя (судя по

\90°

Р и с . IV . 11. Р а с п р е д е л е н и е к о э ф ф и ц и е н т о в теп­

л о о т д а ч и

по

поверхности

с ф е р ы :

/ _ t,-0,3;

/ / — 0-1,2;

/ / / — » = 3,1; IV— о - 4 , 3 мм/сек.

рис. I V . 11, при ф =

85 ~

95°), образуется

воздушный мешок. Эта

часть поверхности (около 50%) весьма неэффективна в тепловом от­ ношении.

Отношение средних коэффициентов теплоотдачи верхней и ниж­ ней половин цилиндра составляет ~ 3 , 3 . Степень неравномерности распределения коэффициентов теплообмена по поверхности сферы

несколько уменьшается с ростом скорости — в исследован­

'мин

ном диапазоне от 5,5 до 4,15. Благодаря отсутствию застойной зоны распределение теплоотдачи по сфере заметно отличается от описан­ ного выше для цилиндра. Эти выводы справедливы при значениях

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

Д л я мелких частиц / высоких і возможно появление застойных зон

в лобовой части.

90

Зависимость, описывающая средний теплообмен слоя со сферой при 173 < Ре < 2350 и указанных выше геометрических характе­ ристиках, имеет вид

Nu

= 3,75

Р е 0 ,

2 5 .

(IV. 6)

Уравнение (IV.6) может быть представлено в виде

Nu

= 5 , 8 3 F o - 0

' 2 5

(IV. 6а)

при 2,67 • Ю - 3 < F o = - ^ — =

— <

3,63

• 10— .

^

D

Р е

 

 

Экспериментальные данные были сопоставлены с результатами расчета процесса нестационарной теплопроводности неограничен­ ного массива с шаровой полостью. Расчет выполняли по прибли­ женной методике [34]. Сравнение экспериментальных и расчетных данных позволяет судить о влиянии на теплоотдачу фактора движе­ ния, характера омывания шара слоем. Это влияние приводит к более медленному падению интенсивности теплообмена с ростом критерия Фурье вследствие перемешивания частиц у поверхности. Значитель­ ные количественные и качественные расхождения с расчетным урав­ нением говорят о том, что процесс конвективного теплообмена слоя со сферой нельзя рассматривать как процесс прогрева (охлаждения) неограниченного массива с шаровой полостью.

IV. 4. ВЛИЯНИЕ ФОРМЫ ПОВЕРХНОСТИ

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

омывания

поверхности

на интенсивность теплообмена. На рис. I V . 12

показаны

зависимости

для

продольного

нестесненного

омывания

 

 

 

 

L

 

пластин (уравнение (III.8))

и цилиндров

(III.5) п р и — =

1, попереч­

но омываемых цилиндров (IV.2), (IV.3) и сферы (IV.6). При продоль­

ном омывании

теплоотдача

заметно улучшается с

ростом

скорости

и не зависит от симплекса - ^ - При поперечном омывании

влияние

скорости

меньше, существенную роль

играет

симплекс

Увели­

чение определяющего

размера

во всех

случаях приводит к ухуд­

шению теплоотдачи,

однако

степень

этого

влияния различна.

В

табл.

I V . 1

приведены

показатели

степени

в зависимостях

 

.

_

_

I D " 1

 

 

 

 

 

a

c=*D

, а »

ѵ", а ^

^ - ^ - j

 

для различных

поверхностей.

 

 

Все

указанные

обстоятельства

обусловлены

особенностями

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

91

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

 

 

 

Nu.-п =

ЛРе'0,22

 

 

 

 

(IV.

7)

где

для

смесей

А =

1,95;

m = 0,33,

для

фракций

А =

4,0; m

=

=

0,48.

Связь

между

определяющими

размерами пластины L и ци­

 

 

 

 

 

линдра

D,

при

которых

коэффици­

 

 

 

 

 

енты

теплоотдачи

для

определен­

 

 

 

 

 

ного материала

и при ѵ = idem бу­

 

 

 

 

 

дут одинаковы, можно легко полу­

 

 

 

 

 

чить из уравнения

(IV.7)

 

 

І = Л 1 о 1 ' 4 4 Ы

( Н

- (IV. 8)

 

эф'

 

 

500

700 /ООО

гооо

зооо Рв

Рис .

I V .

12.

С р а в н е н и е д а н н ы х

д л я

тел различной

ф о р м ы :

/ — пластина;

/ / — продольно

омы­

ваемый цилиндр, о"*3 ';

Ш-

IV —

поперечно

омываемый

цилиндр,

Ü

 

(фракции);

V — то

же,

10,70

D

(смесь);

 

D

-ѵ=70

V/ — сфера

-6,5.

'пред

где для смесей Ах 3,8; т1 =0,66, для фракций Ах 15; т, = 0,96.

Теплоотдача вертикальных ци­ линдров заметно ухудшается с рос-

L

том симплекса — и при преьыше-

D

нии определенного значения і — і

\ D I пред

становится ниже, чем при попереч­ ном омывании. Выражение для пре­ дельной относительной длины име­ ет вид

(IV. 9)

Т а б л и ц а I V .

Показатели степени,

характеризующие влияние

 

на

теплообмен определяющих

ф а к т о р о в

 

Омывание

Поверхность

m

л

к

Нестесненное

Пластина

0

0,5

0,5

продольное

Цилиндр

0

0,38

0,4

 

Lг

 

 

 

 

( n p H £ j = i d e m )

0,33(смесь)

0,28

О.ЗЭ(смесь)

Поперечное

Цилиндр

 

0,48( фрак ­

»

0,24( фрак­

 

 

ции)

ции)

Смешанное

Сфера

0,25

0,75

Н а к л о н н ы е

0,12

0,35

0,53

 

цилиндры

(Ч>=зо°)

92

Соседние файлы в папке книги из ГПНТБ