щ/
КОСМОНАВТИКА
УДК 621:536.2:629.78
*М. Ж. Мурадимов, 1В. В. Двирный, Т. В. Двирный, 1С. Г. Кукушкин, 2В. В. Голованова, 3Е. С. Сидорова
1АО «Информационные спутниковые системы» им. акад. М. Ф. Решетнёва»,
г. Железногорск, Красноярский край, Россия 2Федеральное государственное унитарное предприятие Конструкторское бюро «Арсенал» им. М. В. Фрунзе, г. Санкт-Петербург, Россия 3 Сибирский федеральный университет, г. Красноярск, Россия
ТЕПЛОВАЯ СХЕМА МАЛОГО КОСМИЧЕСКОГО АППАРАТА ТИПА «ЮБИЛЕЙНЫЙ»
и определение параметров теплообменного устройства
Приведен расчет параметров теплообменного аппарата (радиатора) из состава системы терморегулирования (СТР) негерметичного телекоммуникационного спутника малого класса (масса спутника до 600 кг), находящегося на низкой круговой орбите -1500 км.
Ключевые слова: тепловая схема, расчет параметров радиатора, спутник «Юбилейный», «холодный» и «горячий» случаи, внутренние и внешние тепловые потоки, энергетический баланс.
1M. Z. Muradimov, 1V. V. Dvirny, 1 G. V. Dvirny, 1S. G. Kukushkin, 2V. V. Golovanova, 3E. S. Sidorova
JCS «Academician M. F Reshetnev» Information Satellite Systems»,
Zheleznogorsk, Russia
2The Arsenal Design Bureau named after M. V Frunze, St. Petersburg, Russia
3Siberian Federal University, Krasnoyarsk, Russia
THERMAL ciRcuiT oF «YuBILEYNY» sMALL
class spacecraft and heat exchanger parameters determination
Calculations of the heat exchanger (radiator) parameters of the group of temperature-control system in unpressurized small class telecommunication satellite (satellite mass < 600 kg), on low circular orbit -1500 km, are shown.
Key words: thermal circuit, calculation of radiator parameters, satellite «Yubileyny», «hot» and «cold» variant, inboard and off-board heat flow,
energy balance.
Спутниковые системы связи представля- находящих широкое применение в различных ют собой одну из наиболее представительных областях человеческой деятельности. и наиболее быстро развивающихся разновид- Космический аппарат (КА) «Юбилей-
ностей космических информационных систем, ный» запущен с космодрома «Плесецк».
«Юбилейный» создан в АО «Информа-
© МурадимовМ. Ж., Двирный В. В., Двирный Г. В., ционные спутниковые системы» им. акад. Кукушкин С. Г., Голованова В. В., Сидорова Е. С., М. Ф. Решетнёва» (АО «ИСС»). Спутник ОТ-2015 работал гарантированные 12 месяцев и про-
М. Ж. Мурадимов, В. В. Двирный, Г. В. Двирный, С. Г. Кукушкин, В. В. Голованова, Е. С. Сидорова
Тепловая схема малого космического аппарата типа «Юбилейный»
должил функционировать на орбите без отказов и замечаний. Бесперебойная работа КА обеспечена удачными техническими решениями, заложенными при его проектировании. В том числе тепловой схемы в негерметичном исполнении [1].
Одним из важнейших результатов работы КА «Юбилейный» стало обеспечение летной квалификации платформы, предназначенной для бесконтейнерных спутников микрокласса.
Успешная реализация проекта «Юбилейный» положила начало возрождению разработки и изготовления микроспутников в АО «ИСС».
Многофункциональная платформа КА «Юбилейный», сконструированная в АО «ИСС» для аппаратов массой от 30 до 100 кг, получив летную квалификацию, в настоящее время используется при создании всех малых спутников.
На проекте «Юбилейный» ИСС и партнерский вуз СибГАУ отработали новую образовательную технологию для подготовки специалистов-спутникостроителей с высоким уровнем профессиональной компетенции. Студенты СибГАУ принимали непосредственное участие в создании этого спутника. На его борту были размещены созданные в вузе приборы для испытаний в условиях космоса. Кроме того, на базе СибГАУ создан центр управления полетами, в котором студенты в режиме реального времени принимают телеметрическую информацию с «Юбилейного» и других подобных аппаратов, учат ее расшифровывать и управлять спутниками.
КА «Юбилейный» (рис. 1), первоначально задуманный как радиолюбительский спутник, в результате стал настоящей маленькой космической лабораторией и прообразом целой серии подобных аппаратов, с помощью которых ИСС испытывает космосом перспективные приборы [2].
19
В процессе годового движения спутника вместе с Землей вокруг Солнца угол между плоскостью орбиты и направлением на Солнце может изменяться в очень широком диапазоне - практически от 0° до 90°. Поэтому в тепловых анализах спутника выделяют два крайних случая: «горячий» и «холодный».
«Горячий» случай реализуется при максимальном значении солнечной постоянной и углом между плоскостью орбиты и направлением на Солнце около 90°, приведен на рис. 2.
«Горячий» случай характеризуется следующими основными параметрами:
- солнечная постоянная S0 = 1 449 Вт/м2;
- тепловыделение бортовой аппаратуры максимально;
- значения оптических коэффициентов на конец срока активного существования (САС);
- теневой участок орбиты отсутствует.
«Холодный» случай реализуется при минимальном значении солнечной постоянной (максимум имеет место в июне) и углом между плоскостью орбиты и направлением на Солнце около 0°, приведен на рис. 3.
«Холодный» случай характеризуется следующими основными параметрами:
- солнечная постоянная S0 = 1 351 Вт/м2;
- минимальное тепловыделение бортовой аппаратуры;
- значения оптических коэффициентов в начале САС;
- электрообогреватель отключён;
- продолжительность тени Земли максимальна.
Внешние тепловые потоки
Со стороны окружающего космического пространства на спутник действуют лучистые потоки, которые включают в себя прямое и отраженное от планет излучение Солнца, собственную тепловую радиацию
Плоскость орбиты спутника
Рис. 3. «Холодный» случай
планет и излучение открытого космического пространства, так называемое реликтовое излучение, температура которого составляет 3-4 К. Коэффициент черноты космического пространства равен единице.
Пониженное давление в космосе приводит к тому, что внешний теплообмен космических объектов осуществляется, как правило, только излучением. Исключение составляют низколетящие орбитальные объекты. Например, спутники, летящие ниже 200 км, для которых кроме лучистого теплообмена необходимо оценивать тепловое воздействие атмосферы планеты.
Различные виды притока тепла к элементу внешней поверхности спутника зависят от его ориентации относительно Солнца и планет, местоположения его на орбите, планетоцентрического расстояния, физических свойств атмосферы и поверхности планеты, вблизи которой находится спутник.
Отвод тепла от спутника определяется его тепловым излучением в космос. Таким образом, тепловое состояние спутника определяется балансом внутренних и внешних тепловых потоков.
Кроме перечисленных факторов, которые составляют основной энергетический баланс объекта, существуют такие факторы,
как микрометеоритное вещество и пыль планетного происхождения, ультрафиолетовое и рентгеновское излучение Солнца, протоны «солнечного ветра» и радиационные пояса планет, которые оказывают косвенное воздействие на тепловой режим объекта, изменяя свойства его внешних поверхностей (величину теплообмена с космосом).
Солнечный поток, падающий на поверхность спутника, рассчитывается по формуле
Qs=ЛЛ'^'С08Ф, (!)
где Qs - поглощенный солнечный поток, Вт; As - коэффициент поглощения радиационного покрытия внешней обшивки сотовой панели; S0 - плотность солнечного излучения, падающего на единичную поверхность в месте нахождения КА, Вт/м2; F - площадь поверхности, на которую падает солнечное излучение, м2; j - угол между направлением на Солнце и нормалью к поверхности F, град.
Собственное излучение с открытых поверхностей спутника рассчитывается по уравнению Стефана - Больцмана:
Q„ = s-<VF-n (2)
где Q - собственное излучение спутника, Вт; s - коэффициент черноты; s0 - постоянная Стефана - Больцмана, 5,6710-8 Вт/м2К4;
М. Ж. Мурадимов, В. В. Двирный, Г. В. Двирный, С. Г. Кукушкин, В. В. Голованова, Е. С. Сидорова
Тепловая схема малого космического аппарата типа «Юбилейный»
F - площадь излучающей поверхности, м2; T - температура излучающей поверхности, К.
В замкнутой автономной теплоэнергетической системе выделяется поток Q извне подводится Q , сам объект излучает в окружающее пространство Q
Уравнение теплового баланса в общем
виде
Ем тг ■ С ■d ■ T = (бнар - QBH)dt - QBJt> (3)
где m c.,-T. - масса, удельная теплоемкость и температура /-го элемента.
Величина внутреннего теплового потока оценивается по затратам мощности в узлах аппарата и зависит от программы работы оборудования, энергоустановок и системы управления объектом. Наибольшее количество тепла выделяется радиоэлектронной аппаратурой, элементами энергетических и двигательных установок. Отвод тепла осуществляется излучением от поверхности теплообменного устройства (радиатора-излучателя).
Моделирование теплового режима КА и принципы построения его тепловой модели рассмотрены в работе В. Е. Чеботарева и В. Е. Косенко [4].
Каждое тело, если его температура выше абсолютного нуля, непрерывно излучает в окружающее пространство тепло, которое, попадая на другие приборы, частично поглощается, частично переотражается и частично проходит сквозь тело. Поглощенная энергия переходит в тепло, которое, в свою очередь, частично излучается. Отраженная и прошедшая сквозь тело часть энергии также попадает на другие тела и частично ими поглощается, частично отражается и частично проходит сквозь, и так далее, до тех пор, пока излученная телом энергия полностью не распределится между окружающими телами.
Количество тепла, переданное телом 2 телу 1, вычисляется по формуле
Q21 = £1£2G0Flj21(T14 -T24X (4)
где e1 - степень черноты тела 1; е2 - степень черноты тела 2; а0 - постоянная Стефана -Больцмана, 5,6710-8 Вт/м2К4; F1 - площадь тела 1, м2; ф21 - средний угловой коэффициент, или коэффициент облученности; T1 - температура тела 1, К; T2 - температура тела 2, К.
Средний угловой коэффициент рассчитывается по формуле
cos(91) • соб(ф2)
Ф21
— f dF2 f -
F2 Jf2 ■>f 1
n^r
2
-dFi, (5)
где r - расстояние между элементами тел 1 и 2, м; ф1 и ф2 - углы нормалей с линией, соединяющей центры элементарных площадок
dF1 и dFT
Анализ построения тепловой схемы спутника
На внешнюю обшивку радиатора спутника падает солнечный поток. Внешняя обшивка является излучающей.
Тепловая нагрузка, которую должен отводить радиатор, должна быть не менее 300 Вт. Наиболее «комфортной» температурой посадочных мест приборов считаем температуру +30 °С. Радиатор имеет два режима работы: «горячий» и «холодный».
Требуется рассчитать оптимальную площадь радиатора и провести анализ полученных результатов.
Расчет параметров радиатора
В соответствии с исходными данными спутник имеет негерметичное построение, а все тепловыделяющее оборудование располагается непосредственно на панелях силовой сотовой конструкции, поэтому отвод избыточного тепла от спутника целесообразно осуществлять непосредственно с наружной стороны приборных панелей.
Для эффективного отбора тепла от наиболее теплонагруженных приборов и переноса его на излучающие поверхности необходимо применение тепловых труб и управляемого электрообогрева.
Тепловая схема спутника приведена на рис. 4.
Так как спутник имеет 3-осную ориентацию, то радиационные поверхности целесообразно разместить на сторонах ±Z. В этом случае на радиационные поверхности будут воздействовать следующие внешние тепловые потоки:
- собственное излучение Земли;
- отраженное от Земли солнечное излучение;
- прямое солнечное излучение под углом не более 5° (что обусловлено точностью ориентации).
Все поверхности спутника (кроме радиационных поверхностей с фотопреобразова-
21
№ 1 (и) январь-март 2015
22
телем и оптических окон приборов ориентации) закрыты электронно-вакуумной теплоизоляцией (ЭВТИ).
Расчёт площади радиатора и мощности электрообогревателя
Величина суммарного теплового потока на поверхность площадью F выражается как
Qy = Q + Q + Q + Q =
^сол ^.з ^отр ^пр
-пр
= As ■ q ■ F + As ■ q ■ F +
сол отр
+ 8 • q ■ F + Q . (6)
^зем ^пр v '
Значения q , q , q находим исходя
отр сол зем
из полученных путем интерполяции графических зависимостей (рис. 4-6):
- q0Tp ~ q0Tp ■ As = 24 • 0,27 = 6,48 Вт/м2;
- q™ = я™ ■ As =75 ■ °27 = 20,25 Вт/м2;
- qseM = £ ■ ^зем = 0,85 ■32 = 27,2 Вт/м2 •
q, Вт/м2 Средние за виток земные потоки для пластины
200
“а=0
-а=90
-а=45
200 1200 2200 3200 4200 5200 6200 Н, км
Рис. 5
q Вт/м2 Средние за виток солнечные потоки для пластины
450
400
350
300
250
200
150
100
50
200 1200 2200 3200 4200 5200 6200 7200 Н км
-а=45
-а=90
Рис. 6
0
Уравнение теплового баланса
бсол + Q3 + ботр + Qup =£’C0 • Fp 'Tp- (7)
Из выражения (6) находим площадь радиатора F:
F = —.---------&--------------, (8)
Т -г-а-(Js-4С0Л + As-q0Tp + e-q3eu)
________________300_________________
293,164 - 0,85 - а - (20,25 + 6,48 + 27,2)
= 1 м2. (9)
q Вт/м2 180
160
140
120
100
80
60
40
20
0
Средние за виток отраженные потоки
“а=0
-а=45
-а=90
2500 Н, км
0
500
1000
1500
2000
Рис. 7
М. Ж. Мурадимов, В. В. Двирный, Г. В. Двирный, С. Г. Кукушкин, В. В. Голованова, Е. С. Сидорова
Тепловая схема малого космического аппарата типа «Юбилейный»
Найдём температуру радиатора при минимальном тепловыделении приборов спутника Qmin:
rn min _ -
Tp =
(QC(M + Q3 + Qcxp ) +
Q
пр
r-’max
hv
г-a
4
1
100
(20,25 + 27,2 + 6,48) +-
=-35,43 °C. (10)
0,85-a
Найдём требуемое тепловыделение приборов QTpe6 при минимально допустимой температуре работы приборов tm :
ai6 =т4-s-o-F-Q, + Q, +QOTP) • F =
= 263,164 -0,85-o-l--(20,25 + 27,2 + 6,48) = 177,21 Вт. (11)
Найдём мощность электрообогревателя: бэо = Q%* - Qn7 = 177,21 -100 - 80 Вт. (12)
Библиогр афические ссылки
1. ОАО «Информационные спутниковые системы» имени академика М. Ф. Решетнёва» // Сибирский спутник. 2010. № 19 (213), 19 мая.
2. ОАО «Информационные спутниковые системы» имени академика М. Ф. Решетнёва» // Сибирский спутник. 2013. № 11 (340), 20 мая.
3. Двирный В. В., Краев М. В. Малорасходные автономные нагнетатели. Красноярск : Изд-во Краснояр. ун-та, 1985. 152 с.
4. Чеботарев В. Е., Косенко В. Е. Проектирование космических аппаратов информационного обеспечения : учеб. пособие / Сиб. гос. аэрокосм. ун-т. Красноярск, 2011. 488 с., [24] с ил.
Статья поступила в редакцию 25.03.2015 г.
23