CC BY
24РАКЕТНО-КОСМИЧЕСКАЯ ТЕХНИКА
УДК 536.2; 539.3
DOI 10.26732/2618-7957-2018-3-147-156
ТЕМПЕРАТУРНОЕ СОСТОЯНИЕ И ОЦЕНКА ОТКЛОНЕНИЯ ФОРМЫ СФЕРИЧЕСКОЙ ОБОЛОЧКИ КАЛИБРОВОЧНОГО КОСМИЧЕСКОГО АППАРАТА НА ТЕНЕВОМ УЧАСТКЕ ОКОЛОЗЕМНОЙ ОРБИТЫ*
В. С. Зарубин, В. Н. Зимин, Г. Н. Кувыркин
Московский государственный технический университет имени Н. Э. Баумана,
г. Москва, Российская Федерация
Представлена комплексная математическая модель, описывающая изменение формы сферической оболочки калибровочного космического аппарата, вызванное неравномерным распределением температуры по ее поверхности на теневом участке околоземной орбиты. Тепловая часть модели учитывает воздействие на оболочку собственного излучения поверхности Земли и радиационный теплообмен в полости оболочки. Неравномерность облучения поверхности оболочки приводит к неравномерности распределения ее температуры и, как следствие, к отклонению ее формы от сферической, что может повлиять на функциональные характеристики аппарата. Для определения искажения сферической формы оболочки использованы основные положения безмоментной теории оболочек вращения. Количественный анализ представленной модели проведен для установившегося неравномерного распределения температуры по поверхности оболочки, возникающего на различной высоте расположения аппарата над поверхностью Земли. Полученные оценки искажения формы оболочки калибровочного космического аппарата могут быть использованы для сопоставления с допустимыми отклонениями, не нарушающими его функциональные характеристики на теневом участке околоземной орбиты.
Ключевые слова: калибровочный космический аппарат, сферическая оболочка, излучение поверхности Земли, радиационный теплообмен.
Введение
Калибровочные космические аппараты (ККА) применяют для оценки и контроля энергетического потенциала радиолокационного канала наземного комплекса мониторинга движения космических объектов [1; 2]. Распространенным вариантом геометрической формы ККА, называемого также калибровочно-юстировочным, является выполненная с высокой точностью сферическая оболочка. Аналогичную форму имеют пассивные ретрансляторы сигналов и используемые в целях калибровки и юстировки радиолокационной аппаратуры некоторые типы малоразмерных эталонных отражателей [1; 3], а также крупноразмерные надувные конструкции [4-6].
Зарубин В. С., Зимин В. Н., Кувыркин Г. Н., 2018
Работа выполнена в рамках государственных
заданий Минобрнауки РФ
(проекты 9.7784.2017/БЧ и 9.24.22.2017/ПЧ)
Орбиты ККА могут быть как круговыми с высотой около 1000 км, в том числе близкими к полярным, так и эллиптическими с апогеем до 2200 км [2; 3]. Основным фактором, влияющим на форму оболочки ККА, является распределение температуры по ее поверхности, зависящее от солнечного излучения и собственного излучения Земли [7; 8]. На теневом участке околоземной орбиты температурное состояние оболочки ККА наряду с собственным излучением Земли в некоторой степени может зависеть от энергии соударения с молекулами и атомами атмосферных газов и рекомбинации этих молекул [9]. На высоте 300 км над поверхностью Земли оценки плотности локальных тепловых потоков, вызванных этими процессами, не превышают 10 % от плотности потока излучения Земли, а на высоте более 350 км составляют менее 1 % [10].
Оболочку ККА выполняют из полимерной пленки толщиной несколько десятков микрометров, покрытой тонким слоем (толщиной несколько нанометров) напыленного алюминия [11].
*
тнношга
ii'T i','!'
После вывода ККА на околоземную орбиту оболочка принимает сферическую форму путем наполнения ее газом при сравнительно невысоком давлении. При этом возможно получить сферическую оболочку достаточно большого диаметра, что характерно для современных тенденций развертывания на орбите крупногабаритных трансформируемых конструкций [12].
Для орбиты ККА с теневым участком в силу периодической смены интенсивности облучения оболочки происходит существенное изменение ее температурного состояния. Это приводит не 148 только к периодическому изменению формы оболочки, но может повлиять на ресурс работоспособности ее материала в силу циклического термомеханического воздействия. Количественный анализ термомеханического состояния оболочки ККА может быть проведен методами математического моделирования [13], используя ее тепловую модель [9] и определяющие соотношения для безмоментной сферической оболочки [14].
Условия облучения оболочки
Интенсивность собственного излучения Земли, определяющего температурное состояние оболочки ККА на теневом участке орбиты, зависит от температуры участков поверхности Земли, влажности и температуры атмосферы и от состояния облачного покрова [15]. Плотность потока собственного излучения различных участков поверхности Земли, измеренная оптико-электронными приборами, установленными на орбитальной обсерватории, имеет существенный разброс и достигает значения 357 Вт/м2 [16]. Для экваториальной области Земли в летний период плотность этого потока составляет от 175 (при сплошной облачности с верхней границей на высоте 9 км) до 247 Вт/м2 при отсутствии облаков, а для участков в зоне северного полярного круга при тех же условиях от 145 до 225 Вт/м2 [17].
Среднее значение плотности потока собственного излучения Земли можно оценить из условия ее радиационного баланса по формуле [10; 15] д0 = д. (1 - а.) / 4, где д. - солнечная постоянная для среднего расстояния Земли от Солнца (по рекомендации Всемирного центра радиационного мониторинга (Давос, Швейцария) =1368 Вт/м2), а,, - альбедо Земли (среднее значение коэффициента отражения солнечного излучения поверхностью Земли). Однако значение д. (несмотря на принятое название этой величины) не является строго постоянным. Интенсивность солнечного излучения зависит от текущего состояния поверхности Солнца и от положения Земли на ее эллиптической орбите при движении вокруг Солнца (эксцентриситет этой орбиты ~ 0,0167). Многолетние наблюдения показали, что изме-
Том 2
нения солнечной активности слабо влияют на значение д., вызывая отклонения в пределах 1 %, тогда как влияние положения Земли на ее орбите более существенно.
С изменением расстояния Ь от Солнца плотность потока солнечного излучения изменяется пропорционально квадрату отношения Ь0/Ь, где Ь0 ~ 1,49 108 км - среднее расстояние Земли от Солнца. Поэтому даже при достаточно малом эксцентриситете орбиты Земли в перигелии д.=1414 Вт/м2, а в афелии - 1322 Вт/м2 [10; 18]. В инженерных расчетах в случаях, когда увеличение может ухудшить характеристики КА или повлиять на его работоспособность, обычно принимают завышенное значение д. (например, 1420 Вт/м2 [19] и 1400 Вт/м2 [20]). Для последующих числовых расчетов температурного состояния оболочки ККА примем значение д. = 1420 Вт/м2. Тогда при среднем значении а. = 0,37 [18] из условия радиационного баланса Земли получим д0 ~ 215 Вт/м2.
Рис. 1. Расчетная схема
При наблюдении с Земли угловой размер ККА даже при радиусе г0 его сферической оболочки порядка 10 м начиная с высоты H = 200 км аппарата над поверхностью Земли не превышает ё = 10-4 рад. Поэтому при заданном значении H (расстояние между точками O и O" на рис. 1) можно считать, что на оболочку ККА падает собственное излучение Земли с участка ее поверхности, площадь которого зависит от угла
ут = arccos(1/(1+^)), где ц = H/R0, R0 ~ 6371 км - средний радиус поверхности Земли. Выделим на этом участке поверхности кольцевой слой бесконечно малой ширины R0dy (у е [0; yj) площадью
dS' = 2%R sin y dу (на рис. 1 этот слой заштрихован).
Пренебрежимо малое значение ё позволяет положить угол ф между нормалью к выделенному кольцевому слою и направлением на различные участки поверхности оболочки равным в + у,
Температурное состояние и оценка отклонения формы сферической оболочки
Я(0) = Я'(0)/ д0, дх = д[ / д0 и < = 2 /(7ст02^0)
а квадрат расстояния между этим слоем и ККА равным.
12 = + (1 + л)2-2(1 + у) .
Если принять, что распределение по отдельным направлениям излучения с поверхности Земли определяется законом Ламберта [21], то суммарный поток этого излучения, падающий на поверхность оболочки ККА, составит [9]:
Q = 2(1 -^(2 + Л) / (1 + Л))
Но распределение плотности потока
этого излучения по поверхности оболочки будет неравномерным, зависящим от угла у, отсчитываемого от ближайшей к поверхности Земли точки на ее поверхности (рис. 1) вдоль любой дуги большого круга радиусом проходящей через эту точку, в которой плотность потока будет равна [9] д'0 = д'(0) = д0/(1 + л)2 .
С возрастанием угла у функция убывает и в соответствии с принятым выше допущением о возможности пренебречь угловым размером оболочки примет нулевое значение при ут = п - ут, т. е. д'(щт) = 0 и, кроме того, dq,(у)/dу\у=у = 0. Из симметрии распределения д'(у) относительно прямой О'О следует dq' (у)^у\^0 = 0. Этим условиям удовлетворяет аппроксимирующее соотношение
q\ш) = q'0 cos
2 кш . . 2 кш +q sin2 r
(1)
где
Tm '.—i
2 Лу . , — | cos -sin у dу —
0 2У m
1 - cos уm - 2(у m / л)2
2 - 2(уm / л)2
rm
=| sin
2 • J о 1" cos у m 2 Sin ydy = 2- "
0
В итоге получим
4" (Vm /ЯГ
4i= q
4 - (Wm / Л)'
2(1 - cos Wm )
fx Ул(2 + л) 1 - cos wm - 2(Wm / Л)2 1 + л 2(1 + л)2(1 -(Wm / Л)2
Так как ут = л-ут =%-arctg^л(2 + л) и cos =-co s ym = -1/ (1 + л) , отношения
и угол ут можно представить функциями лишь одного аргумента п (рис. 2). На рис. 3 построены рассчитанные с использованием формулы (1) графики зависимости от отношения у = У / ут относительной плотности 1 (V) = 1 '(V) /10 падающего на поверхность оболочки собственного излучения поверхности Земли для различных значений параметра п.
2¥т ¥ш
в котором коэффициент найдем из баланса потоков падающего излучения в виде равенства
Q = 2^г02(д'01о + д[/,) ?
Рис. 2. Графики зависимостейот параметра 1 - д(0); 2 - ; 3 - Q ; 4 - Ут
1.0 0.8 0.6 0.4 0.2
0
■2
ч ^
\ Оч \
—-6 Ч.Ч\
5 V s
\\\ *» V4
6 . 'ч \ iv
0.2
0.4
0.6
0.8
¥
_Рис. 3. Графики зависимостей д (у) от отношения у = у / у т при различных значениях параметра 1 - 0,03; 2 - 0,05; 3 - 0,1; 4 - 0,2; 5 - 0,5; 6 - 1
Температурное состояние сферической оболочки
Оценки термического сопротивления оболочки из тонкой алюминированной полимерной пленки в направлении нормали к поверхности по-
149
тнношга
казывают, что температуру оболочки можно принять однородной по ее толщине. Действительно, оболочка из каптона [11], аналогом которого является полиимидная пленка ПМФ-352 с двусторонним покрытием политетрафторэтиленом (тефлоном) [22; 23] с коэффициентом теплопроводности X = 0,17 Вт/(мК), при толщине Н = 30 мкм имеет весьма малое термическое сопротивление Н / X ~ 1,8-10-4 м2 К/Вт. Слой напыленного алюминия толщиной Н' =10 нм при коэффициенте теплопроводности X' = 206 Вт/(мК) [24] обладает термическим сопротивлением Н'/Х'~ 5• 10-11 м2 К/Вт. Даже 150 при завышенной оценке в 200 Вт/м2 возможного значения плотности теплового потока, проходящего через оболочку, получим также завышенную оценку перепада температуры по толщине оболочки, не превышающую 0,04 К.
Эффективный коэффициент теплопроводности оболочки в тангенциальных направлениях, оцениваемый по формуле X* = (ХН + Х'Н) / (Н + Н1) составит не более 0,18 Вт/(м-К). Это позволяет при формировании тепловой модели оболочки не учитывать перенос тепловой энергии путем теплопроводности вдоль поверхности оболочки. Малое значение полной теплоемкости оболочки позволяет считать ее температурное состояние установившимся при движении ККА на теневом участке орбиты.
Собственное излучение Земли относится к инфракрасной части спектра и для средней температуры излучающей поверхности Земли, оцениваемой значением 288 К, наибольшая интенсивность этого излучения соответствует, согласно закону Вина [21], длине волны около 10 мкм. При близких к указанному значению температурах допустимо считать коэффициент Л' поглощения излучения поверхностью равным ее коэффициенту е' излучения, часто называемому «степенью черноты».
Материал оболочки примем непрозрачным для излучения. Свойства внутренней поверхности оболочки примем соответствующими модели диффузно-серой поверхности [21], что позволяет считать равными коэффициенты Л поглощения и е излучения этой поверхности. Все упомянутые оптические характеристики внешней и внутренней поверхности оболочки полагаем неизменными на каждой из поверхностей. При наличии в оболочке оборудования ККА считаем, что занимаемый им объем достаточно мал. Это дает возможность при анализе переноса излучения внутри оболочки не учитывать возможного взаимного затенения участков ее внутренней поверхности. Отсюда следует, что при осесимме-тричном относительно прямой О'О (рис. 1) распределении плотности падающего собственного излучения Земли искомое установившееся распределение температуры в оболочке также долж-
Том 2
но быть осесимметричным, описываемым функцией Т(у) лишь одного аргумента у. Поскольку давление газа в полости оболочки ККА диаметром даже порядка метра не превышает значения в 10 Па [7], а для более крупных оболочек имеет еще меньшее значение, переносом тепловой энергии в этой полости за счет теплопроводности газа (даже в случае использования высокотеплопроводного гелия) можно пренебречь и учитывать теплоперенос лишь путем излучения.
Принятые допущения позволяют перейти к рассмотрению модели непрозрачной для излучения сферической оболочки нулевой толщины радиусом г0, равным радиусу средней поверхности реальной оболочки ККА. Двумя близкими сечениями, перпендикулярными прямой О'О (рис. 1), выделим из такой условной оболочки кольцевую полосу, положение которой характеризует угол Ее внешняя поверхность поглощает часть Л'д'(щ) потока, падающего на оболочку излучения Земли плотностью д'(^), и испускает поток собственного излучения плотностью е'а0Т4(р), где а0 ~ 5,6693 • 10-8 Вт/(м2 • К) - постоянная Стефана-Больцмана. Разность плотностей этих потоков составляет плотность д*(р) результирующего излучения с внутренней поверхности этой полосы, поступающего в полость оболочки. Тогда уравнение баланса потоков для единичной площадки внешней поверхности кольцевой полосы можно представить в виде
(2)
A'q'(W) -s'g0T4(v) = q(w).
В силу вогнутости внутренней поверхности оболочки на рассматриваемую полосу от соседних кольцевых полос падает поток излучения плотностью q*. Часть Aq* этого потока поглощается, а часть (1 - A) q*, отражаясь, снова поступает в полость оболочки (коэффициенты поглощения и отражения для непрозрачной оболочки в сумме составляют единицу). Особенность полости сферической оболочки состоит в том, что для нее q* = const [7-9] т.е. не зависит от угла В данном случае при е' = A' имеем
q* = q= Q / (4w02)
- среднее по внешней поверхности оболочки значение плотности потока падающего на нее излучения Земли. В итоге уравнение баланса тепловых потоков для единичной площадки внутренней поверхности кольцевой полосы примет вид
Aq' + q\y) = so0T4 (у)
(3)
Исключив из соотношений (2) и (3) величину д*, при е = Л получим формулу
= (Ад'+А'д'(у))/((А + А')а0), (4)
Температурное состояние и оценка отклонения формы сферической оболочки
определяющую установившееся распределение температуры по поверхности оболочки. Из этой формулы, в частности, следует, что температура участка оболочки при у > ут, на которую не падает излучение Земли, не зависит от угловых координат. Эту температуру определяют параметр П и значения А и А'.
Оценка отклонения формы оболочки от сферической
Оболочку ККА при неизменной по ее толщине температуре и нагружении только внутренним давлением р допустимо считать безмомент-ной [14], а ее полимерный материал изотропным. Поскольку температурное состояние оболочки осесимметрично относительно прямой О'О (рис. 1), она будет симметричной относительно оси, совпадающей с этой прямой, т. е. безмо-ментной оболочкой вращения относительно этой оси. Для такой оболочки уравнение равновесия в проекции на нормаль к ее средней поверхности (уравнение Лапласа) имеет вид
NJ R + N2 / R2 = p
(5)
где N и Ы2 - погонные усилия соответственно в меридиональном и окружном направлениях; а и Я2 - радиусы кривизны средней поверхности оболочки в этих направлениях. Для полных деформаций безмоментной оболочки вращения в этих направлениях справедливы соотношения [14]:
Sj = (du / d0 + w) / R , s2 (w + u ctg 0) / R = ur / r,
(6)
поверхности. При нагружении оболочки вращения только внутренним избыточным давлением p в ее поперечном сечении радиусом r суммарная сила в направлении оси будет равна P = 2nrN1sinO = nr2p. Отсюда, с учетом равенства r = R2 sinO и уравнения (5), следует N1 = pr / (2 • sinO) = pR2 / 2и N2 = (pR2/2)(2 - R2/Ri). Для сферической оболочки в этом случае R1 = R2 и N1 = N2 = pR/2. Если при осесимметричном распределении температуры по поверхности оболочки ККА в первом приближении принять R1 ~ R2 ~ r0, то получим N2 ~ pR2 / 2 = Nb что, согласно соотношениям (7), приводит к приближенному равенству е1 ~ г2. Тогда из формул (6) следует du/dO ~ u • ctg O, т. е. u ~ C • sinO, где C = const. В итоге принятого приближения из первых равенств (6) и (7) находим
w - pr¡ (1 -v) / (2Eh) + ar0 AT(у) - C cos 0
Перемещение точек средней поверхности оболочки в направлении ее оси равно [14] u1 = u • sinO - w • cosO, что с учетом формул для перемещений u и w дает
щ» C — (pr (1—v) / (2Eh)+ aAT(^))rocos0
Если за нуль отсчета перемещения u1 принять ближайшую к Земле точку при у = O = 0 (рис. 1), то получим
С=p (r0 + w0)2(l - v) / (2 Eh)+a(r0 + w0) AT (щ)
U » (p (r0 + w0) (1 - v) / (2Eh) + aAT(у) ) •
■( ro + wo)(1 - cos 0)
Из вторых равенств (6) и (7) следует
(8)
где и и w - перемещения точки средней поверхности в меридиональном направлении и в направлении нормали к этой поверхности; 9 - угол между осью оболочки и этой нормалью; иг = и • cos9 + w • sin9 - перемещение перпендикулярно оси оболочки; г - радиус контура поперечного сечения средней поверхности плоскостью, перпендикулярной оси оболочки.
В предположении упругой работы материала оболочки, согласно обобщенному закону Гука, для полных деформаций в меридиональном и окружном направлениях запишем
= -^Жй+аДГ(¥)), е2 = (Щ - уЩ) / (Eh + а ДТ(у)),
где V, Е и а - соответственно коэффициент Пуассона, модуль Юнга и температурный коэффициент линейного расширения материала оболочки; ДТ(у) = Т(у) - Т0 - приращение температуры оболочки по сравнению с однородной по ее поверхности температурой Т0, при которой оболочка имела сферическую форму с
и - (рг0 (1 - V) / (2ЕИ) + аАТ(у)) г0 sm 0 (9)
Из формул (8) и (9) видно, что изменение давления в рамках принятого приближения не искажает форму сферической оболочки, а лишь приводит к одинаковому во всех точках средней поверхности некоторому приращению ее первоначального радиуса г0. Отклонение формы оболочки от сферической непосредственно связано с неравномерностью распределения температуры по ее поверхности, причем эта неравномерность возникает, согласно формуле (4), только на облучаемой Землей части этой поверхности. Из соотношений (8) и (9) следует, что перемещения, вызванные лишь неравномерностью распределения температуры, при отождествлении в первом приближении угла 9 с углом у можно оценить по формулам
u* = aAT (у) r0 sin у, u* = aAT(у) г0 (1 - cos у)
151
и
■жив
АППАРАТЫ И
ПХШПЁ^^Н
№ 3 (25) 2018
152
Результаты расчетов
Степень влияния излучения Земли на температурное состояние оболочки ККА зависит от значения параметра п. Для сравнения проведем вычисления при четырех значениях этого параметра, равных 0,03; 0,05; 0,1 и 0,2, что соответствует значениям высоты Н над ее поверхностью Земли, равным 191,3; 318,6; 637,1 и 1274,3 км.
Пусть оболочка ККА изготовлена из упомянутой выше полиимидной пленки ПМФ-352 с напылением алюминием на внешней поверхности. Коэффициент поглощения А внутренней поверхности такой оболочки будет не менее 0,5 [25]. Увеличение значения А приводит к выравниванию распределения температуры по поверхности оболочки и уменьшению наибольшего значения Т(0) температуры при у = 0 в точке, ближайшей к поверхности Земли. Для сравнения при проведении расчетов примем два крайних значения А = 0,5 и А = 1 (второе значение для полости оболочки соответствует модели абсолютно черного тела [21]).
Для коэффициента поглощения слоя напыленного алюминия в его исходном состоянии после напыления можно использовать линейную зависимость А'(Т) = 0,025 - 0,00007 • (300 - Т) от температуры Т [26], установленную в промежутке от 100 К до 300 К для охлажденного расплава алюминия с неокисленной поверхностью. Однако при контакте слоя алюминиевого покрытия с атмосферным воздухом происходит окисление поверхности этого слоя и увеличение этого коэффициента. Кроме того, в процессе пребывания в космосе поверхность оболочки подвержена микрометеоритному воздействию, что приводит к возникновению шероховатости поверхности покрытия и возрастанию его коэффициента излучения. В силу отсутствия определенности величины коэффициента А' расчеты проведены для трех его дискретных значений, равных 0,015; 0,05 и 0,15.
Результаты расчетов распределения температуры по поверхности сферической оболочки ККА в зависимости от отношения у = у / ут для указанных значений параметра п приведены на рис. 4-7. На всех этих рисунках сплошные кривые построены при значении А = 0,5, а штрих-пунктирные - при значении А = 1. Кривые 1, 2 и 3 отвечают значениям соответственно 0,015; 0,05 и 0,15.
При фиксированном сочетании коэффициентов поглощения А и А' наибольшая температура сферической оболочки, соответствующая значению у = 0, и температура изотермического участка, отвечающая значению у = 1, убывают по мере возрастания параметра п. Это вызвано увеличением высоты ККА над поверхностью Земли
Том 2
и связано с уменьшением потока излучения, падающего на поверхность оболочки. Но даже при П = 0,2 общий уровень температур (рис. 7) остается допустимым для такого полимерного материала, как пленка ПМФ-352 [27; 28]. Из представленных результатов видно, что при фиксированных значениях параметра п и коэффициента А неравномерность распределения температуры увеличивается с увеличением коэффициента А', а при неизменных значениях п и А' переход от А = 0,5 к А = 1, как и следовало ожидать, в связи с выравниванием распределения температуры по поверхности оболочки температура изотермического участка оболочки возрастает.
Найденные распределения температуры по поверхности оболочки ККА могут быть использованы для оценки по формулам (10) отклонения формы этой оболочки от исходной сферической.
Рис. 4. Графики зависимости Т(у) от у при значении параметра п = 0,03
Г К
Рис. 5. Графики зависимости Т(у) от у при значении параметра п = 0,05
Температурное состояние и оценка отклонения формы сферической оболочки
На рис. 8 и рис. 9 штрихпунктирной дугой полуокружности представлен контур продольного сечения средней поверхности оболочки при однородном распределении температуры Т0, равной рассчитанному значению Т(п/2) (отсчет угла у происходит от положительного направления оси Охх).
198г
194 190 186 182 178 174 170
------- л.
>2
\ -------
\
>3
0
0.2
0.4 0.6 0.8 \|/
Рис. 6. Графики зависимости Т(у) от у при значении параметра п = 0,1
Рис. 7. Графики зависимости Т ( у) от значении параметра п = 0,2
У
У
при
I
X,
\
"о.
-г0 О
Рис. 8. Отклонение формы контура средней поверхности оболочки при п = 0,05
го
Светлыми кружками на этой полуокружности отмечена часть контура, на которую падает излучение Земли, а светлыми квадратами - его часть с постоянной температурой. Смещение темных кружков и квадратов на сплошной кривой по отношению к светлым позволяет сравнить степень локального искажения исходного контура средней поверхности при рассчитанном неоднородном распределении температуры Т(у). На участке контура с темными квадратами температура постоянна. Поэтому эта часть оболочки сохраняет сферическую форму, но с меньшим по сравнению с г0 радиусом г*, определяемым отношением (г0 - г) /г0 =а(Т(п / 2) - Т(п)). Для наглядности относительное изменение этого радиуса и смещение темных кружков по отношению к светлым увеличено на рис. 8 в 0,02 • К-1/а раз, а на рис. 9 - в 0,05 • К-1/а раз.
Рис. 9. Отклонение формы контура средней поверхности оболочки при п = 0,2
Из сравнения штрихпунктирной полуокружности со сплошной кривой видно, что в окрестности наиболее близкой к Земле точки оболочки ее радиусы кривизны увеличиваются по сравнению с исходным значением, равным г0. Затем по мере возрастания угла у при продолжающемся росте значения Я2 радиус кривизны Ях начинает убывать. При у > п/2 после достижения радиусом Я2 максимальной величины его значение начинает убывать и сближается с возрастающим значением Ях до их совпадения при у ~ 2,83 ~ 162° на рис. 8 и при у ~ 2,56 ~ 146° на рис. 9.
Заключение
Тепловая часть представленной комплексной математической модели, описывающая установившееся температурное состояние алюми-нированной полимерной сферической оболочки ККА, находящегося на затененном участке околоземной орбиты, позволяет установить зависимость распределения температуры по поверхности этой оболочки от высоты космического аппарата над поверхностью Земли и от коэффициентов излучения внешней и внутренней по-
153
тнношга
ii'T
Том 2
верхностей оболочки. Эта зависимость исполь- зультаты могут быть полезны для прогноза рабо-
зована для оценки отклонения формы оболочки тоспособности ККА на затененном участке око-
от сферической, вызванного неравномерностью лоземной орбиты. распределения температуры. Полученные ре-
Список литературы
1. Малые космические аппараты информационного обеспечения / Под ред. В. Ф. Фатеева. М. : Радиотехника, 2010. 320 с.
2. Ракеты и космические аппараты конструкторского бюро «Южное» / А. Н. Мащенко, В. Н. Паппо-Корыстин,
B. А. Пащенко, В. Г. Васильев; под общей редакцией Генерального конструктора С. Н. Конюхова.
154 Днепропетровск : ГКБ «Южное» им. М.К. Янгеля, 2000. 240 с.
3. Тарасенко М. В. Военные аспекты советской космонавтики. М. : Агентство Российской печати, ТОО «Николь», 1992. 164 с.
4. Freeland R. E., Bilyeu G. D., Veal G. R. Development of flight hardwarefor a large, inflatable-deployable antenna experiment // Acta Astronautica, 1996, vol. 38, no. 4-8, pp. 251-260. doi: 10.1016/0094-5765(96)00030-6.
5. Chodimella S. P., Moore J., Otto J., Fang H. Design evaluation of alarge aperture deployable antenna // AIAA Papers, 2006, no. 1603, pp. 1-20. doi: 10.2514/6.2006-1603.
6. Зимин В. Н., Крылов А. В., Мешковский В. Е., Сдобников А. Н., Файзуллин Ф. Р., Чурилин С. А. Особенности расчета раскрытия крупногабаритных трансформируемых конструкций различных конфигураций // Наука и образование. 2014. №10. С. 179-191. doi: 10.7463/1014.0728802.
7. Зарубин В. С., Зимин В. Н., Кувыркин Г. Н. Температурное состояние и отклонение формы сферической оболочки космического калибровочноюстировочного аппарата // Аэрокосмический научный журнал. 2016. № 1. С. 27-45. doi: 10.7463/aersp.0116.0831867.
8. Zarubin V. S., Zimin V. N., Kuvyrkin G. N. Temperature distribution in the spherical shell of a gauge-alignment spacecraft // Journal of Applied Mechanics and Technical Physics, 2017, vol. 58, no. 6, pp. 1083-1090. doi: 10.1134/S0021894417060141.
9. Зарубин В. С., Зимин В. Н., Кувыркин Г. Н. Температурное состояние оболочки калибровочного космического аппарата на теневом участке околоземной орбиты // Аэрокосмический научный журнал. 2016. № 4.
C. 22-37. doi: 10.7463/aersp.0416.0846458.
10. Основы теплопередачи в авиационной и ракетно-космической технике / Под ред. В. К. Кошкина. М. : Машиностроение, 1975. 624 с.
11. Optical Calibration Sphere Experiment [Электронный ресурс]. URL: https://directory.eoportal.org/web/eoportal/ satellite-missions/j/jawsat# footback10 (дата обращения 25.08.2018).
12. Зимин В. Н. К вопросу моделирования и расчета динамики раскрытия трансформируемых космических конструкций // Оборонная техника. 2006. № 1. С. 123-127.
13. Zarubin V. S., Kuvyrkin G. N. Mathematical Modeling of Thermomechanical Processes under Intense Thermal Effect // High Temperature, 2003, vol. 41, no. 2, pp. 257-265. doi: 10.1023/A:1023390021091.
14. Балабух Л. И., Колесников К. С., Зарубин В. С., Алфутов Н. А., Усюкин В. И., Чижов В. Ф. Основы строительной механики ракет. М. : Высшая школа, 1969. 496 с.
15. Марчук Г. И., Кондратьев К. Я., Козодеров В. В. Радиационный баланс Земли: ключевые аспекты. М. : Наука, 1988. 224 с.
16. Павлов А. В. Оптико-электронные приборы. М. : Энергия, 1974. 360 с.
17. Кондратьев К. Я., Дьяченко Л. Н., Козодеров В. В. Радиационный баланс Земли. Л. : Гидрометеоиздат, 1988. 350 с.
18. Gilmor D. G. Spacecraft thermal control handbook. vol. 1. Fundamental Technologies / El Segundo, California, Aerospace Press, 2002. 836 p.
19. Комарова М. А. Температурные условия на корпусе узлового модуля на этапе автономного полета к международной космической станции // Известия РАН. Энергетика. 2012. № 2. С. 23-30.
20. Гукало А. А., Грибков А. С. Оптимизация температуры плоского и крестообразного холодильника-излучателя космической ядерной энергетической установки с учетом внешнего теплового излучения // Известия РАН. Энергетика. 2012. № 2. С. 103-110.
21. Зигель Р., Хауэлл Дж. Теплообмен излучением / Пер. с англ. М. : Мир, 1975. 934 с.
22. Бюллер К. У Тепло- и термостойкие полимеры / Пер. с нем. М. : Химия, 1984. 1056 с.
23. Аналитический портал химической промышленности [Электронный ресурс]. URL: http://www.newchemistry. ru/production.php?cat-id=52&catparent=7&level=3 (дата обращения 28.08.2018).
Температурное состояние и оценка отклонения формы сферической оболочки
24. Чиркин В. С. Теплофизические свойства материалов ядерной техники : справочник. М. : Атомиздат, 1968. 484 с.
25. Многослойные и комбинированные пленочные материалы и изделия из них: интернет-ресурс [Электронный ресурс]. URL: rccgroup.ru/plast/wp-content/uploads/3.9.-Ivanenko.ppt (дата обращения 28.08.2018).
26. Излучательные свойства твердых материалов : справочник / Под общ. ред. А. Е. Шейндлина. М. : Энергия, 1974. 472 с.
27. Обзор базовых свойств и сфер применения политетрафторэтилена (PTFE) [Электронный ресурс]. URL: http://www.engplast.ru/entry.php?198&id=2 (дата обращения 28.08.2018).
28. Комков М. А., Сабельников В. В., Баслык К. П. Конструктивно-технологический анализ трубопроводов, намотанных из полиимидофторопластовой пленки // Вестник МГТУ им. Н. Э. Баумана. Сер. Машиностроение. 2012. С. 78-86.
155
TEMPERATURE STATE AND ESTIMATION OF DEVIATION IN THE SHAPE OF THE SPHERICAL SHELL OF THE CALIBRATION SPACECRAFT IN THE SHADOW SECTION OF THE NEAR-
EARTH ORBIT
V. S. Zarubin, V. N. Zimin, G. N. Kuvyrkin
Bauman Moscow State Technical University, Moscow, Russian Federation
A complex mathematical model describing alteration in a shape of a spherical shell of a calibration spacecraft is presented. The alteration is caused by the uneven temperature distribution over the shell surface when it is in the shadow section of the near-earth orbit. The thermal part of the model takes into account influence upon the shell of the own Earth radiation and the radiation heat transfer in the shell cavity. Unevenness of the shell surface radiation leads to the unevenness of the temperature distribution over the shell surface and thus to the alteration in shape of the spherical shell, that can affect the device Junctional characteristics. To determine the distortion in shape of the spherical shell, the fundamentals of the membrane theory of rotational shells were used. The quantitative analysis of the discussed model was carried outfor cases of the steady non-uniform temperature distribution along the shell surface arising at different heights of the spacecraft above the Earth's surface. Obtained estimations of the distortion in shape of the shell of the calibration spacecraft can be compared with the permissible deviations that do not violate the spacecraft functional characteristics in the shadow section of the near-earth orbit.
Keywords: calibration spacecraft, spherical shell, Earth surface radiation, radiation heat transfer.
References
1. Fateev V F. Malye kosmicheskie apparaty informatsionnogo obespecheniya [Small spacecrafts for information support]. Moscow, Radiotekhnika Publ., 2010, 320 p. (In Russian)
2. Mashchenko A. N., Pappo-Korystin V. N., Pashchenko V. A., Vasil'ev V. G. Rakety i kosmicheskie apparaty konstruk-torskogo byuro «Yuzhnoe» [Rockets and spacecraft from the engineering department «Yuzhnoe»]. Dnepropetrovsk, The Yangel' state engineering department «Yuzhnoe» Publ., 2000, 240 p. (In Russian)
3. Tarasenko M. V. Voennye aspekty sovetskoy kosmonavtiki [Military aspects of the soviet astronautics]. Moscow, Russian press agency special partnership «Nikol'» Publ., 1992, 164 p. (In Russian)
4. Freeland R. E., Bilyeu G. D., Veal G. R. Development of flight hardwarefor a large, inflatable-deployable antenna experiment // Acta Astronautica, 1996, vol. 38, no. 4-8, pp. 251-260. doi: 10.1016/0094-5765(96)00030-6.
5. Chodimella S. P., Moore J., Otto J., Fang H. Design evaluation of alarge aperture deployable antenna // AIAA Papers, 2006, no. 1603, pp.1-20. doi: 10.2514/6.2006-1603.
6. Zimin V. N., Krylov A. V, Meshkovskiy V E., Sdobnikov A. N., Fayzullin F. R., Churilin S. A. Osobennosti rascheta raskrytiya krupnogabaritnyh transformiruemyh konstrukcij razlichnyh konfiguracij [Special features of deploy calculation for large transformable structures of different configuration] // Nauka i obrazovanie, MGTU im. N. E. Baumana, 2014, no. 10, pp. 179-191. doi: 10.7463/1014.0728802. (In Russian)
7. Zarubin V. S., Zimin V N., Kuvyrkin G. N. Temperaturnoe sostoyanie i otklonenie formy sfericheskoj obolochki kos-micheskogo kalibrovochnoyustirovochnogo apparata [Temperature state and deviation in a shape of a spherical shell of a space calibrator] // Aerokosmicheskiy nauchny zhurnal, 2016, no. 1, pp. 27-45. doi: 10.7463/aersp.0116.0831867. (In Russian)
55ft
8. Zarubin V S., Zimin V. N., Kuvyrkin G. N. Temperature distribution in the spherical shell of a gauge-alignment
№ 3 (25) 2018
TOM 2
spacecraft // Journal of Applied Mechanics and Technical Physics, 2017, vol. 58, no. 6, pp. 1083-1090. doi: 10.1134/ S0021894417060141. (In Russian)
9. Zarubin V. S., Zimin V. N., Kuvyrkin G. N. Temperaturnoe sostoyanie obolochki kalibrovochnogo kosmicheskogo apparata na tenevom uchastke okolozemnoj orbity [Temperature state of a shell of a gauge spacecraft on the shadow section of the near-earth orbit]. Aerokosmicheskiy nauchny zhurnal, 2016, no. 4, pp. 22-37. doi: 10.7463/aer-sp.0416.0846458. (In Russian)
10. Koshkin V. K. Osnovy teploperedachi v aviatsionnoy i raketno-kosmicheckoy tekhnike [Fundamentals of heat-transfer in aviation and space-rocket engineering]. Moscow, Mashinostroenie Publ., 1975, 624 p. (In Russian)
11. Optical Calibration Sphere Experiment. Available at: https://directory.eoportal.org/web/eoportal/satellite-missions/j/ jawsat# footback10 (accessed 25.08.2018).
12. Zimin V N. K voprosu modelirovaniya i rascheta dinamiki raskrytiya transformiruemyh kosmicheskih konstrukcij 156 [On modeling and calculation of deployment dynamics of transformable space structures] // Oboronnaya tekhnica,
2006, no. 1, pp. 123-127. (In Russian)
13. Zarubin V. S., Kuvyrkin G. N. Mathematical Modeling of Thermomechanical Processes under Intense Thermal Effect // High Temperature, 2003, vol. 41, no. 2, pp. 257-265. doi: 10.1023/A:1023390021091.
14. Balabukh L. I., Kolesnikov K. S., Zarubin V S., Alfutov N. A., Usyukin V. B., Chizhov V F. Osnovy stroitel'noy mekhan-iki raket [Fundamentals of the rocket structural mechanics]. Moscow, Visshaya shkola Publ., 1969, 496 p. (In Russian)
15. Marchuk G. I., Kondrat'ev K. Ya., Kozoderov V. V. Radiatsionny balans Zemli: klyuchevye aspekty [Earth-radiation budget: main aspects]. Moscow, Nauka Pabl., 1988, 224 p. (In Russian)
16. Pavlov A. B. Optiko-elertronnyepribory [Optoelectronic devices]. Moscow, Energiya Pabl., 1974, 380 p. (In Russian)
17. Kondrat'ev K. Ya., D'yachenko L. N., Kozoderov V. V. Radiatsionny balans Zemli [Earth-radiation budget]. Leningrad, Gidrometeoizdat Pabl., 1988, 350 p. (In Russian)
18. Gilmor D. G. Spacecraft thermal control handbook. vol. 1. Fundamental Technologies / El Segundo, California, Aerospace Press, 2002. 836 p.
19. Komarova M. A. Temperaturnye usloviya na korpuse uzlovogo modulya na ehtape avtonomnogo poleta k mezhdun-arodnoj kosmicheskoj stancii [Temperature conditions at the main module airframe on the stage of autonomy flight to the international space station] // Izvestiya RAN. Energetika, 2012, no. 2, pp. 23-30. (In Russian)
20. Gukalo A. A., Gribkov A. S. Optimizaciya temperatury ploskogo i krestoobraznogo holodil'nika-izluchatelya kosmicheskoj yadernoj ehnergeticheskoj ustanovki s uchetom vneshnego teplovogo izlucheniya [Temperature optimization for a flat cross-formed cooler-emitter of a space nuclear power system that takes into account the outer heat radiation] // Izvestiya RAN. Energetika, 2012, no. 2, pp. 103-110. (In Russian)
21. Zigel' R., Khauell G. Teploobmen izlucheniem [Radiative heat exchange]. Moscow, Mir Publ., 1975, 934 p. (In Russian)
22. Byuller K. U. Teplo- i termostoykie polimery [Heat-resistant polymers]. Moscow, Khimiya Publ., 1984, 1056 p. (In Russian)
23.Analiticheskiy portal khimicheskoy promyshlennosti [Analytical portal of the chemical industry]. Available at: http:// www.newchemistry.ru/production.php?cat-id=52&catparent=7&level=3 (accessed 28.08.2018). (In Russian)
24. Chirkin V. S. Teplofizicheskie svoystva materialov yadernoy tekhniki. Spravochnik [Heat-transfer properties of the nuclear engineering materials. Reference book]. Moscow, Atomizdat Publ., 1968, 484 p. (In Russian)
25. Mnogosloynye i kombinirovannye plenochnye materialy i izdeliya iz nikh [Multilayer and combined film materials and product made from them]. Available at: http://rccgroup.ru/plast/wp-content/uploads/3.9.-Ivanenko.ppt (accessed 28.08.2018). (In Russian)
26. Sheydlin A. E. Izluchatel'nye svoystva tverdykh materialov [Emitting properties of the hard materials]. Moscow, Energiya Publ., 1974, 472 p. (In Russian)
27. Obzor bazovykh svoystv i sfer primeneniya politetraftoretilena (PTFE) [Review of fundamental properties and applications of polytetrafluoroethylene (PTFE)]. Available at: http://www.engplast.ru/entry.php?198&id=2 (accessed 28.08.2018). (In Russian)
28. Komkov M. A., Sabel'nikov V. V, Baslyk K. P. Konstruktivno-tekhnologicheskij analiz truboprovodov, namotannyh iz poliimidoftoroplastovoj plenki [Constructively-technological analysis of pipe conduits wound with polyimidofluoro-plastic film] // Vestnik MGTU im. N. E. Baumana, Mashinostroenie, 2012, pp. 78-86. (In Russian).
