Проектирование системы терморегулирования с двухфазным теплоносителем для космического аппарата производительностью 7 кВт

Обложка

Цитировать

Полный текст

Аннотация

В работе представлен методологический подход к проектированию системы терморегулирования космического аппарата (КА) с прокачкой теплоносителя хладопроизводителъностью до 7,0 кВт. Рассмотрены два варианта конструкций. В компоновочной схеме построения КА отсутствует герметичный приборный контейнер, а все тепловыделяющее оборудование располагается непосредственно на панелях силовой конструкции, поэтому отвод избыточного тепла от КА осуществляется непосредственно с наружной стороны приборных панелей.

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

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

Второй вариант отличается от первого тем, что для выравнивания температур панелей модуля полезной нагрузки каждая тепловая труба панели «Север» соединена с тепловыми трубами панели «Юг».

Из результатов сравнительного анализа бюджета масс и энергетической эффективности следует вывод: для систем терморегулирования с ДФК с капиллярной прокачкой наиболее предпочтительным является второй вариант, который при остальных одинаковых характеристиках имеет меньшую массу. Удельная массо-энергетическая характеристика такой системы составляет ~ 22,9 кг/кВт.

Полный текст

Введение

Система терморегулирования космического аппарата (СТР КА) предназначена для поддержания в требуемых диапазонах рабочих температур для всего установленного оборудования, на участке выведения и орбитальном функционировании, а также для обеспечения теплового режима оборудования КА совместно с технологическими средствами термостатирования при наземных испытаниях в диапазоне как рабочих температур электронных систем полезной нагрузки [1–4], так и специальных криогенных температур для работы инфракрасных рефлекторов [5].

Наиболее распространенные жидкостные СТР имеют существенный недостаток в части удельных массо-энергетических характеристик из-за большей массы заправки теплоносителя, использующего только теплоемкостную аккумуляцию тепла, и, как следствие, недопустимости паровой фазы на контурном центробежном насосе, хотя модели и тепловые балансы таких систем достаточно проработаны [6–9].

Методологический подход к проектированию элементов СТР, основанных на тепловых трубах, и определению коэффициентов теплоотдачи от жидкостно-паровых фаз в достаточном объеме представлен в источниках [10–18].

В настоящей работе рассматриваются пассивные СТР на основе тепловых труб (ТТ) и капиллярных насосов (КН) с использованием теплоты фазового перехода, т. е. двухфазные контуры (ДФК) без активных насосных агрегатов.

Постановка задачи исследования

Обычно КА, работающий на геостационарной орбите, функционально структурирован двумя модулями (рис. 1):

– модулем служебных систем (МСС);
– модулем полезной нагрузки (МПН).

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

Для эффективного отбора тепла от наиболее теплонагруженных приборов и переноса его на излучающую поверхность применяются либо ТТ, либо ДФК.

Трехосная ориентация КА на стационарной орбите предопределяет тот факт, что каждая из панелей со стороны осей ±Z поочередно в течение примерно полугода подвержена воздействию солнечного потока, снижая тем самым эффективность одной из радиационных поверхностей.

Кроме того, модульный принцип будущих КА приводит к тому, что большая часть тепловой нагрузки (принято до 7,0 кВт) приходится именно на МПН. Отвод такого количества тепловой энергии от МПН требует значительных радиационных поверхностей, которые в составе конструкции МПН, как правило, создать невозможно.

Перечисленные задачи – отбор тепла от оборудования, перенос тепла с наиболее нагретых поверхностей на менее нагретые и перенос тепла от оборудования МПН на радиационные поверхности МСС – можно решить с помощью как ДФК, так и ТТ.

Дополнительные радиационные поверхности могут быть выполнены в виде откидных панелей-радиаторов, при этом угол наклона раскрываемых радиаторов предпочтительно должен составлять ~20° к панелям ±Z.

Угол наклона в ~20° раскрываемых радиаторов к панели ±Z выбран исходя из того, чтобы, с одной стороны, максимально уменьшить влияние струи двигателей коррекции на терморегулирующее покрытие панелей-радиаторов и, с другой стороны, минимизировать радиационный теплообмен между раскрываемыми радиаторами и панелью ±Z. Раскрываемые радиаторы представляют собой сотопанели со встроенными трактами ДФК. С обеих сторон панели покрыты терморегулирующим оптическим покрытием (ОСО-С). В сложенном состоянии раскрываемые радиаторы располагаются вдоль панелей МСС, расположенных со стороны осей ±Y. Радиаторные панели «Север – Юг» МПН и МСС также представляют собой сотопанели со встроенными ТТ или трактами ДФК.

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

Для облегчения раскрытия панелей с ДФК после вывода КА на орбиту, а также обеспечения возможности стыковки гидравлических разъемов, входящих в разные сборочные единицы конструкции КА, в состав СТР вводятся гибкие секции.

Сравнительная оценка различных схем СТР для автоматических КА, работающей на стационарной орбите, проводилась при следующих исходных данных (ИД):

– максимальная тепловая нагрузка МПН 7,0 кВт;
– тепловая нагрузка МСС постоянна и равна 700 Вт;
– рабочая температура МПН (расчетная) 40 °С (диапазон работоспособности от –30 до +60 °С);
– теплоноситель – аммиак.

Принятые допущения при расчетах:

– внешний тепловой поток, поглощаемый одним квадратным метром панели «Север» МПН, МСС и раскрываемых панелей-радиаторов составляет 160 Вт;
– температура окружающей среды для панелей МПН и МСС составляет 3 °К, для раскрываемых панелей –125 °С;
– коэффициент полезного действия при излучении с панелей МПН и МСС равен 0,88, а при излучении с раскрываемых панелей 0,85;
– площадь каждой панели МПН составляет 5 м2; площадь панелей МСС – 3,5 м2;
– шаг между ТТ на панелях ~75 мм, на раскрываемом радиаторе ~150 мм.

Описание исследуемых компоновочных решений СТР

Общая компоновка КА приведена на рис. 1.

 

Рис. 1. Общая схема размещения радиационных поверхностей на автоматических КА

Fig. 1. General layout of radiation surfaces on automatic spacecraft

 

В варианте 1, представленном на рис. 2, ТТ МПН соединены с МСС через КТТ раскрываемого радиатора (или через капиллярный испаритель) и тепловая нагрузка МПН распределяется между панелью МСС и раскрываемыми радиаторами.

Для удобства компоновки теплообменного узла корпус капиллярного испарителя имеет форму профиля с двумя полочками. Ширина полочек корпуса испарителя принята равной удвоенной по ширине полочки ТТ. К одной полочке корпуса испарителя подводятся ТТ от МПН, к другой – от МСС. Это обстоятельство вводит ограничение по компоновке панелей – расстояние между панелями.

Панели «Север – Юг» МПН должны быть на две высоты капиллярного насоса больше, чем расстояние между панелями «Север – Юг» МСС.

 

Рис. 2. Схема пассивной ДФК СТР на 7 кВт. Вариант 1 - схема соединения ТТ МПН с ТТ МСС через КТТ

Fig. 2. The scheme of the passive DF PAGE for 7 kW. Option 1 – scheme of connection of TT MPN with TT MSS using CTT

 

Каждый раскрываемый радиатор состоит из двух КТТ, конденсаторы которых сформированы алюминиевыми трубами (внутренним диаметром 3 мм), соединенными с испарительными зонами ТТ раскрываемых радиаторов. Конденсаторные трубки одного контура переплетаются
с трубками другого контура, тем самым обеспечивая эффективность панели, позволяющую удовлетворять температурным требованиям при выходе из строя одного контура. Конденсаторные трубки и ТТ одного контура связаны с одной обшивкой раскрываемого радиатора, конденсаторные трубки и ТТ другого контура – с другой обшивкой. Таким образом получается два ДФК на раскрываемом радиаторе. Конденсаторные трубки и ТТ одного контура сдвинуты относительно трубок другого контура примерно на 75 мм.

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

 

Рис. 3. Схема раскрываемого радиатора на КТТ (вариант 1): КН – капиллярный насос; ГШ – гибкий шланг; ТТ – тепловая труба_____ контур 1; -------- контур 2

Fig. 3. The diagram of the radiator being opened on the CTT.  Option 1 ______ contour 1; ------- contour 2

 

Вариант 2 отличается от варианта 1 тем, что для выравнивания температур панелей МПН каждая ТТ панели «Север» МПН соединена с ТТ панели «Юг» МПН. Выравнивание температур панелей МСС производится по половине ТТ из-за малости передаваемой мощности каждой тепловой трубой (менее 35 Вт).

Принципиальная схема варианта 2 приведена на рис. 4.

 

Рис. 4. Схема пассивной ДФ СТР на 7 кВт. Вариант 2 – схема с тепловыми трубами между панелями «Север» и «Юг»

Fig. 4. The scheme of the passive DF PAGE for 7 kW. Option 2 – а scheme with heat pipes between the North and South panels

 

Алгоритм и результаты расчетов СТР

Расчеты проводятся в следующей последовательности:

  • определяется форма панелей МПН, МСС и раскрываемых панелей. Из конструктивных соображений принято, что каждая теплоотводящая панель МПН и МСС содержит в среднем 16 ТТ, расстояние между которыми составляет 160 мм. Тогда панели МПН («Север» и «Юг») имеют размеры 2560´1950 мм, панели МСС – 2560´1370 мм при площади панелей 5 и 3,5 м2 соответственно;
  • определяется мощность, отводимая от панели МПН к КТТ или панели МСС;
  • ;
  • определяется мощность, отводимая одной ТТ;
  • определяется перепад температур на ТТ МПН. По отводимой мощности и перепаду температур уточняется профиль, из которого изготавливается ТТ;
  • определяется тепловая нагрузка, распределяемая между МСС и раскрываемыми панелями;
  • вычисляются размеры раскрываемых панелей;
  • определяются геометрические параметры капиллярного испарителя КТТ;
  • определяются характеристики ДФК.

Результаты расчетов приведены табл. 1 и 2. Основное отличие в массе связано с разницей площадей раскрываемых панелей и масс ТТ.

 

Таблица 1. Бюджет масс для варианта 1 

Наименование

Кол-во в СТР,
площадь

МАССА

одной единицы*

суммарная

1. Сотопанели:

1.1. Раскрываемая панель Р1

2,1 м2

5,78 кг/м2

12,14 кг

1.2. Раскрываемая панель Р2

2,1 м2

5,78 кг/м2

12,14 кг

1.3. Раскрываемая панель Р3

2,1 м2

5,78 кг/м2

12,14 кг

1.4. Раскрываемая панель Р4

2,1 м2

5,78 кг/м2

12,14 кг

2. Тепловые трубы:

 

 

 

2.1. «Север» МПН

16 шт.

0,572 кг

9,15 кг

2.2. «Юг» МПН

16 шт.

0,572 кг

9,15 кг

2.3. «Север» МСС

16 шт.

0,376 кг

6,02 кг

2.4. «Юг» МСС

16 шт.

0,376 кг

6,02 кг

2.5. Раскрываемая панель Р1

7 шт.

0,443 кг

3,10 кг

2.6. Раскрываемая панель Р2

7 шт.

0,443 кг

3,10 кг

2.7. Раскрываемая панель Р3

7 шт.

0,443 кг

3,10 кг

2.8. Раскрываемая панель Р4

7 шт.

0,443 кг

3,10 кг

3. КТТ:

3.1. «Север»

8 шт.

2,9 кг

23,2 кг

3.2. «Юг»

8 шт.

2,9 кг

23,2 кг

4. Механизм раскрытия

4 шт.

5,0 кг

20,00 кг

5. Электрообогреватели (ЭО) трубопроводов двигателей ориентации (ДУОС)

6 шт.

0,17 кг

1,00 кг

6. ЭО замещающие

10 шт.

0,3 кг

3,00 кг

7. Экранно-вакуумная теплоизоляция (ЭВТИ)

Компл.

3,0 кг

3,00 кг

ИТОГО:

 

 

184,86 кг

Примечание. *Для раскрываемых паналей приведена удельная масса на 1 м2 поверхности панели.

 

Таблица 2. Бюджет масс для варианта 2

Наименование

Кол-во в СТР,
площадь

МАССА

одной единицы

суммарная

1. Сотопанели:

1.1. Раскрываемая панель Р1

1,25 м2

5,78 кг/м2

7,22 кг

1.2. Раскрываемая панель Р2

1,25 м2

5,78 кг/м2

7,22 кг

1.3. Раскрываемая панель Р3

1,25 м2

5,78 кг/м2

7,22 кг

1.4. Раскрываемая панель Р4

1,25 м2

5,78 кг/м2

7,22 кг

2. Тепловые трубы:

2.1. «Север» МПН

16 шт.

0,723 кг

11,57 кг

2.2. «Юг» МПН

16 шт.

0,723 кг

11,57 кг

2.3. «Север» МСС

16 шт.

0,600 кг

9,60 кг

2.4. «Юг» МСС

16 шт.

0,600 кг

9,60 кг

2.5. Раскрываемая панель Р1

7 шт.

0,253 кг

1,77 кг

2.6. Раскрываемая панель Р2

7 шт.

0,253 кг

1,77 кг

2.7. Раскрываемая панель Р3

7 шт.

0,253 кг

1,77 кг

2.8. Раскрываемая панель Р4

7 шт.

0,253 кг

1,77 кг

2.9. Соединительные панели МПН

16 шт.

0,490 кг

7,84 кг

2.10. Соединительные панели МПН

8 шт.

0,490 кг

3,92 кг

3. КТТ: 

3.1. «Север»

8 шт.

2,9 кг

21,68 кг

3.2. «Юг»

8 шт.

2,9 кг

21,68 кг

4. Механизм раскрытия

4 шт.

5,0 кг

20,00 кг

5. ЭО трубопроводов ДУОС

6 шт.

0,17 кг

1,00 кг

6. ЭО замещающие

10 шт.

0,3 кг

3,00 кг

7. ЭВТИ

Компл.

3,0 кг

3,00 кг

ИТОГО:

 

 

160,42 кг

 

В расчетах также учитывалось:

– наличие второй полочки в местах крепления ТТ как между собой, так и в местах крепления с трактами КТТ;
– отсутствие полочек в транспортных зонах, где нет соединения с обшивкой панели;
– соединительные ТТ принимались длиной 2,6 м, причем длины испарительной и конденсационной зон принимались равными 0,3 м, транспортная зона 2,0 м.

Оценка полученных результатов

По результатам сравнительного анализа массовых бюджетов по двум вариантам можно сделать следующий вывод: по массе из систем с капиллярной прокачкой преимуществом обладает вариант 2, который при остальных одинаковых характеристиках имеет меньшую массу. Удельная массо-энергетическая характеристика такой системы составляет ~22,9 кг/кВт. Это обусловлено выравниванием температур панелей, расположенных на сторонах «Север» и «Юг», что выравнивает нагрузку на ТТ. Благодаря этому уменьшается температурный перепад вдоль тепловой цепочки, что, в свою очередь, приводит к уменьшению площади раскрываемых радиаторов.

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

Сами по себе элементы ДФК с капиллярной прокачкой не потребляют энергию. Однако не исключено, что для запуска КН раскрываемых радиаторов потребуются электрообогреватели (ЭО) мощностью ~5 Вт, установленные на компенсационных камерах КН (гидроаккумуляторах), которые могут использоваться для регулирования сброса тепла раскрываемыми радиаторами.

Диапазон и точность регулирования обеспечиваются замещающими ЭО (при понижении нагрузки оборудования МПН), расположенными на панелях МПН и МСС. Точность регулирования температуры плат МПН и МСС зависит от системы управления включением – выключением замещающих ЭО.

Надежность СТР с ДФК с капиллярной прокачкой высокая и зависит, в основном, от надежности ЭО, которые можно резервировать без заметного увеличения массы.

Заключение

Из результатов сравнительного анализа бюджета масс и энергетической эффективности следует вывод: для СТР с ДФК с капиллярной прокачкой наиболее предпочтительным является вариант 2, который при оостальных динаковых характеристиках имеет меньшую массу. Удельная массо-энергетическая характеристика такой системы составляет ~ 22,9 кг/кВт.

Однако на этапе эскизного проектирования основные конструктивные параметры СТР с ДФК, характеристики основных базовых элементов выбранного варианта могут быть уточнены.

×

Об авторах

Олег Валентинович Шилкин

Акционерное общество «Информационные спутниковые системы» имени академика М.Ф. Решетнева»

Автор, ответственный за переписку.
Email: gift_23j@mail.ru

главный конструктор ОКР «СТР с ДФК», начальник сектора 3604

Россия, 662972, г. Железногорск Красноярского края, ул. Ленина, 52

Анатолий Петрович Колесников

Акционерное общество «Информационные спутниковые системы» имени академика М.Ф. Решетнева»

Email: shilkin@iss-reshetnev.ru

начальник отдела 360

Россия, 662972, Железногорск Красноярского края, ул. Ленина, 52

Александр Викторович Делков

Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева

Email: delkov-mx01@mail.ru

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

Россия, 660037, Красноярск, проси, им. газеты «Красноярский рабочий», 31

Александ Анатольевич Кишкин

Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева

Email: spsp99@mail.ru

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

Россия, 660037, Красноярск, проспект им. газеты «Красноярский рабочий», 31

Юлия Николаевна Шевченко

Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева

Email: gift_23j@mail.ru

аспирант кафедры холодильной, криогенной техники и кондиционирования

Россия, 660037, Красноярск, проси, им. газеты «Красноярский рабочий», 31

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

  1. Meseguer J., Perez-Grande I., Sanz-Andres A. Spacecraft thermal control. Cambridge. UK: Woodhead Publishing Limited, 2012. 413 p.
  2. Gilmore D. G. Spacecraft thennal control handbook. Tire Aerospace Corporation Press. 2002. 413 p.
  3. Крушенко Г. Г.. Голованова В. В. Совершенствование системы терморегулирования космических аппаратов // Вестник СибГАУ. 2014. № 3 (55). С. 185-189.
  4. Chebotarev V. Е., Zimin I. I. Procedure for evaluating the effective use range of the unified space platforms // Сибирский журнал науки и технологий. 2018. Т. 19, № 3. С. 532-537. doi: 10.31772/2587-6066-2018- 19-3-532-537.
  5. Проектирование системы пассивного охлаждения бортового комплекса космического аппарата / О. В. Шилкин. А. А. Кишкин, А. А. Зуев и др. // Вестник Московского авиац. ин-та. 2021. Т. 28. № 2. С. 96-106. doi: 10.34759/vst-2021-2-96-106.
  6. Two-dimensional thermal model of the thermal control system for nonhermetic formation spacecraft / F. V. Tanasienko. Y. N. Shevchenko, A. V. Delkov и др. // Сибирский журнал науки и технологий. 2018. Т. 19. № 3. С. 445-451. doi: 10.31772/2587-6066-2018- 19-3-445-451.
  7. Определяющие тепловые сопротивления в модели жидкостного контура системы терморегулирования космического аппарата / Ю. Н. Шевченко, А. А. Кишкин, Ф. В. Танасиенко и др. // Сибирский журнал науки и технологий. 2019. Т. 20. № 3. С. 366-374.
  8. Delcov А. V., Hodenkov A. A.. Zhuikov D. A. Mathematical modeling of single-phase thermal control system of the spacecraft // Proceedings of 12th Intern. Conf, on Actual Problems of Electronic Instrument Engineering. APEIE 2014. 2014. P. 591-593.
  9. Delcov A. V., Hodenkov A. A.. Zhuikov D. A. Numerical modeling and analyzing of conjugate radiation-convective heat transfer of fin-tube radiator of spacecraft // IOP Conference Series: Materials Science and Engineering. 2015. Vol. 93. No. 012007. P. 1-6.
  10. Зуев А. А., Назаров В. П.. Арнгольд А. А. Определение локального коэффициента теплоотдачи с использованием модели температурного пограничного слоя в полостях вращения газовых турбин // Вестник Московского авиац. ин-та. 2019. Т. 26, № 2. С. 99-115.
  11. Определение локального коэффициента теплоотдачи с использованием модели температурного пограничного слоя с конвективной составляющей в полостях вращения ТНА ЖРД / А. А. Зуев, В. Ю. Пиунов. В. П. Назаров. А. А. Арнгольд // Вестник Юж.-Урал. гос. ун-та. Серия: Машиностроение. 2019. Т. 19. № 2. С. 30-44. doi: 10.14529/enginl90203.
  12. Weybume D. W. Approximate heat transfer coefficients based on variable thermophysical properties for laminar flow over a uniformly heated flat plate // Heat andMass Transfer. 2008. Vol. 44, Iss. 7. P. 805-813. doi: 10.1007/s00231-007-0306-z.
  13. Weybume D. W. New thickness and shape parameters for the boundary layer velocity profile // Experimental Thermal and Fluid Science. 2014. Vol. 54. P. 22-28. doi: 10.1016/j.expthermflusci. 2014.01.008.
  14. Triple diffusive mixed convection from an exponentially decreasing mainstream velocity / P. M. Patil, M. Roy, A. Shashikant et al. // International Journal of Heat and Mass Transfer. 2018. Vol. 124. P. 298-306. doi: 10.1016/j.ijheatmasstransfer.2018.03.052.
  15. Improved velocity and temperature profiles for integral solution in the laminar boundary layer flow on a semi-infinite flat plate / S. M. Seyyedi, A. S. Dogonchi, M. Hashemi-Tilehnoee et al. // Heat Transfer - Asian Research. 2019. Vol. 48, Iss. l.P. 182-215. doi: 10.1002/htj.21378.
  16. Denarie A., Aprile M.. Motta M. Heat transmission over long pipes: New model for fast and accurate district heating simulations // Energy. 2019. Vol. 166. P. 267-276. doi: 10.1016/j.energy. 2018.09.186.
  17. Прямолинейное равномерное течение газов с теплоотдачей в энергетических установках летательных аппаратов / М. И. Толстопятов, А. А. Зуев, А. А. Кишкин и др. // Вестник СибГАУ. 2012. №4(44). С. 134-139.
  18. Analysis of efficiency of systems for control of the thermal regime of spacecraft / A. V. Delkov. A. A. Kishkin, N. A. Lavrov et al. // Chemical and Petroleum Engineering. 2016. No. 9. P. 714-719.

Дополнительные файлы

Доп. файлы
Действие
1. JATS XML
2. Рис. 1. Общая схема размещения радиационных поверхностей на автоматических КА

Скачать (256KB)
3. Рис. 2. Схема пассивной ДФК СТР на 7 кВт. Вариант 1 - схема соединения ТТ МПН с ТТ МСС через КТТ

Скачать (948KB)
4. Рис. 3. Схема раскрываемого радиатора на КТТ (вариант 1):

Скачать (526KB)
5. Рис. 4. Схема пассивной ДФ СТР на 7 кВт. Вариант 2 – схема с тепловыми трубами между панелями «Север» и «Юг»

Скачать (976KB)

© Шилкин О.В., Колесников А.П., Кишкин А.А., Шевченко Ю.Н., 2023

Creative Commons License
Эта статья доступна по лицензии Creative Commons Attribution 4.0 International License.

Данный сайт использует cookie-файлы

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

О куки-файлах