Например, Бобцов

МЕТОДИКА РАСЧЕТА ЗВЕЗДНОЙ ВЕЛИЧИНЫ МЕЖДУНАРОДНОЙ КОСМИЧЕСКОЙ СТАНЦИИ

ОПТИЧЕСКИЕ И ОПТИКО-ЭЛЕКТРОННЫЕ ПРИБОРЫ И СИСТЕМЫ
УДК 524.3
В. М. ТЫМКУЛ, Л. В. ТЫМКУЛ, Ю. А. ФЕСЬКО, К. В. КУДРЯШОВ, С. В. ЛУГОВСКИХ
МЕТОДИКА РАСЧЕТА ЗВЕЗДНОЙ ВЕЛИЧИНЫ МЕЖДУНАРОДНОЙ КОСМИЧЕСКОЙ СТАНЦИИ
Рассматривается методика расчета звездной величины Международной космической станции по отраженному солнечному излучению. Методика основана на использовании коэффициента габаритной яркости как оптической характеристики отраженного излучения станции и представлении ее поверхности в виде суперпозиции цилиндрических, конических, сферических и плоских поверхностей. Составлены алгоритм и программа расчета звездной величины станции в зависимости от углов ее ориентации на Солнце и наблюдателя. Приведены результаты расчета и их сравнительный анализ с данными из сети Интернет и результатами эксперимента.
Ключевые слова: звездная величина, Международная космическая станция, отражение, коэффициент габаритной яркости.
В Международном космическом эксперименте уже несколько лет принимает участие Международная космическая станция (МКС) — самый крупный из имеющихся искусственных спутников Земли (ИСЗ), звездная величина m которого достигает значений, сравнимых с звездной величиной, например, Сириуса (m = –1,47) и Канопуса (m = –0,72). Однако звездная величина МКС не является постоянной за счет изменения ориентации станции относительно точки наблюдения и Солнца.
Знание звездной величины как МКС, так и других крупногабаритных ИСЗ необходимо, в первую очередь, при расчетах и проектировании оптико-электронных систем ориентации и навигации [1, 2]. Кроме того, звездная величина искусственных космических объектов представляет интерес для исследований в астрономии и астрофизике, а также ряде разделов прикладной физики и оптики. Так, например, в работе [3] рассматривается проблема использования в космической геодезии свойств отражения солнечного света поверхностями крупногабаритных космических аппаратов и ИСЗ, в том числе спутников системы ГЛОНАСС, для анализа влияния светового давления на их орбиты.
Настоящая статья посвящена разработке методики определения звездной величины МКС в зависимости от ее ориентации в пространстве. Для решения поставленной задачи геометрическая модель поверхности МКС представляется как суперпозиция локальных элементов цилиндрической, конической, сферической и плоской поверхностей. При этом отражение солнечного излучения элементами корпуса станции носит диффузный или направленнорассеивающий характер, а отражение от элементов солнечных батарей является зеркальным.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 5

6 В. М. Тымкул, Л. В. Тымкул, Ю. А. Фесько и др.

Поток солнечных лучей отражается от освещенных элементов поверхности МКС и дос-

тигает наблюдателя, находящегося на поверхности Земли. Благодаря этому наблюдатель мо-

жет видеть солнечный свет, отраженный станцией. Используя такую оптическую характери-

стику отражения от объемных тел, как коэффициент габаритной яркости (КГЯ) [4], можно

определить силу отраженного станцией света через освещенность Солнца вне атмосферы (ЕС) по формуле

∑IΣ (ΘC , ϕC, Θн,ϕн )

=

EС π

k
ρiβi (ΘC, ϕC , Θн , ϕн ) Аi (Θн , ϕн ) ,
i=1

(1)

где ρi — коэффициент отражения i-го элемента МКС; Аi — площадь проекции i-го элемента

МКС в направлении наблюдателя; βi — нормированный КГЯ i-го элемента МКС; Θн, ΘС, ϕн,

ϕС — сферические координаты направлений на Солнце (ΘС, ϕС) и на наблюдателя (Θн, ϕн). Зная расстояние l от МКС до поверхности Земли, можно перейти от силы света к осве-

щенности Ea, формируемой станцией в точке наблюдения и определяемой по формуле

Ea = IΣτa / l2 ,

(2)

где τа — коэффициент пропускания атмосферы Земли. На основании закона Погсона [5], который устанавливает логарифмическую зависи-

мость ощущения глаза от раздражения его светом, звездная величина МКС, регистрируемая

на уровне моря, выражается через освещенность Ea:

m = −14, 01 − 2,51⋅ lg Еа .

(3)

В формуле (1) присутствуют нормированные значения коэффициента β локальных по-

верхностей, определяемые для элементов корпуса с диффузным покрытием следующим обра-

зом. Так, индикатриса нормированного КГЯ для цилиндра с коэффициентом сжатия tц=rц/hц,

где rц и hц — радиус основания и высота цилиндра, определяется как [6]

βц (ΘС, ϕС , Θн , ϕн )

=

2(2 sin ϕн

1 + πtц

cos ϕн

) {sin ϕн

sin ϕС

cos ΘС

cos Θн

×

×[π − (ΘС + Θн ) + (sin ΘС cos ΘС ) + sin Θн cos Θн ] + sin ϕС sin ϕн sin(ΘС + Θн )×

×(− cos2 Θн + cos2 ΘС ) + sin ϕн sin ϕС sin ΘС sin Θн ×

×[π − (ΘС − Θн ) − (sin ΘС cos ΘС − sin Θн cos Θн )] + 2πtцχ(ϕн , ϕС ) cos ϕн cos ϕС},

(4)

где χ(ϕн, ϕС)=1, когда ϕС и ϕн одновременно заданы в области [π/2, π]; в противном случае
χ(ϕн, ϕС)=0; ϕС, ϕн — угол между осью цилиндра, совпадающей с осью OZ декартовой системы координат станции и направлениями соответственно на источник света и на наблюдателя;

ΘС, Θн — угол между осью ОХ декартовой системы координат и проекцией данных направлений на координатную плоскость XOY.
Индикатриса нормированного КГЯ для усеченного конуса с коэффициентом сжатия

tк = (r2–r1)/hк (r2 > r1), где r1, r2 — радиусы большого и малого оснований конуса, определяется выражением [6]

βк (ΘС , ϕС , Θн , ϕн ) =

tк2

tк + 1 cos ϕн

⎡ ⎢⎣

tк2 2

sin

ϕн

sin

ϕС

cos(ΘС



Θн

)

+

cos

ϕн

cos

ϕС

⎤ ⎥⎦

.

(5)

При этом выполняются условия π–arctg t < ϕС и ϕн ≤ π, т.е. освещается и визируется только боковая поверхность усеченного конуса.

Индикатриса нормированного КГЯ для сферы имеет вид [7]

βсф (ΘС , ϕС, Θн , ϕн ) = (2 3π) ((π − δ) cos δ + sin δ) ,

(6)

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 5

Методика расчета звездной величины Международной космической станции
где cos δ = sin ϕн sin ϕС cos(Θн − ΘС ) + cos ϕн cos ϕС .
Индикатрисы нормированного КГЯ солнечных батарей определяются по формуле

7 (7)

βб (ΘС , ϕС , Θн , ϕн ) = β(α0, αн )δ(α0 − αн ) ,

(8)

где β(α0, αн) — индикатриса коэффициента яркости материала поверхности солнечных батарей; δ(α0 – αн) — дельта-функция Дирака, численно равная

δ(α0



αн

)

=

⎧1, ⎨⎩0,

α0 α0

= αн; ≠ αн,

(9)

причем α0, αн — это углы между нормалью к поверхности батареи и направлением соответственно на Солнце и наблюдателя.

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

при условии α0=αн и нахождении этих углов в одной плоскости. Методика расчета КГЯ сложных поверхностей с направленно-рассеивающим покрытием

и, в частности, цилиндрических, сферических и конических поверхностей подробно рассмот-

рена в работе [8].

Зная геометрические размеры всех элементов МКС [9], площади их проекций в зависи-

мости от углов Θн, ϕн можно определить следующим образом: — площадь проекции цилиндра — по формуле:

Ац (Θн , ϕн ) = 2rцhц sin ϕн + πrц2 cos ϕн ;

(10)

— площадь проекции усеченного конуса — по формуле

Ак (Θн , ϕн ) =

=

⎡ ⎢a ⎢⎣

2

⎛ ⎜ ⎝

π 2

+

arcsin



a tg ϕн

⎞ ⎟





cos ϕн

+ ahк

1−

a2 hк2tg2ϕн

sin

ϕн

⎤ ⎥ ⎥⎦

+

δ(ϕн

)πa2

cos ϕн

,

(11)

где

δ(ϕн

)

=

⎧1 ⎨⎩0

при при

ϕн > π 2, ϕн ≤ π 2;

а = r2 + r1;

(12)

— площадь проекции сферы

Acф = π rc2ф ,

(13)

где rcф — радиус сферы.

Таким образом, методика расчета звездной величины МКС при различных углах ориентации станции на Солнце и наблюдателя основывается на соотношениях (3)—(13).

Рассмотрим алгоритм расчета звездной величины m.

1. Выбираются значения углов ΘС, ϕС и Θн, ϕн. 2. Задаются значения коэффициентов отражения ρi локальных элементов МКС. 3. На основе формул (4)—(6), (8) и (10), (11), (13) применительно к геометрической мо-
дели станции рассчитывается значение суммы в соотношении (1).

4. Вычисляется солнечная постоянная EC вне атмосферы:

EC

=

W

(TC π

)

2π(1 −

cos

α)

,

где W(TC) — светимость Солнца при абсолютной температуре TC, вычисляемая по формуле Стефана — Больцмана: W (TC ) = σTC4 (σ — постоянная Стефана — Больцмана); 2α =32′ — угловой размер солнечного диска.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 5

8 В. М. Тымкул, Л. В. Тымкул, Ю. А. Фесько и др. 5. На основе формул (2) и (3) вычисляются соответственно величина Ea и искомая звезд-
ная величина m. В качестве примера на рис. 1 приведены результаты компьютерного расчета распределе-
ния значений звездной величины МКС по предлагаемой методике при случайных значениях углов ΘС, Θн, ϕС, ϕн (по оси абсцисс показан номер Nк компьютерного эксперимента). При этом значения ρi для всех элементов корпуса станции принимались одинаковыми, а функция распределения случайных значений углов ΘС, Θн, ϕС, ϕн соответствовала равномерному закону.
m 4
3
2
1
0 10 20 30 40 Nк Рис. 1
Согласно расчетам значения звездной величины МКС в данном случае сосредоточены в интервале от +3 до +0,5.
При обработке запроса значений звездной величины МКС с сайта [10] (для наблюдений из Новосибирска в период с 30 апреля по 25 мая 2004 г.) были получены результаты распределения значений т, показанные на рис. 2 (здесь Nн — номер практических наблюдений).
m 2
1
0
–1 0 10 20 30 40 Nн Рис. 2
Согласно этому распределению значения звездной величины МКС сосредоточены в интервале от +1,8 до –0,5.
Для оценки корректности предлагаемой модели расчета звездной величины МКС в этот же период были проведены визуальные наблюдения пролета станции над Новосибирском с помощью телескопа с измерением звездной величины методом Аргеландера [11]. По результатам этих наблюдений значения m находились в интервале от +1 до 0. Имеющиеся различия результатов расчета с данными из сети Интернет и данными экспериментов можно объяснить, по-видимому, выбором случайных значений углов ориентации станции на Солнце
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 5

Методика расчета звездной величины Международной космической станции

9

и в направлении наблюдателя, а также использованием одинаковых значений коэффициентов отражения ρi для всех элементов корпуса станции.
В заключение следует отметить, что по предлагаемой методике можно проводить не только оценочные расчеты звездной величины МКС при случайном выборе углов ориентации станции, но и определять точные значения величины m, если будут заданы действительные
значения коэффициентов ρi и конкретные значения углов ΘС, ϕС и Θн, ϕн.

СПИСОК ЛИТЕРАТУРЫ

1. Изнар А. Н., Павлов А. В., Федоров Б. Ф. Оптико-электронные приборы космических аппаратов. М.: Машиностроение, 1972. 368 с.

2. Ивандиков Я. М. Оптико-электронные приборы для ориентации и навигации космических аппаратов. М.: Машиностроение, 1971. 200 с.

3. Сурнин Ю. В. История создания и развития межкафедральной научно-исследовательской лаборатории космической геодезии // Вестн. СГГА: науч.-техн. журн. Новосибирск: СГГА, 2010. Вып. 2(13). С. 128—145.

4. Холопов Г. К., Шуба Ю. А. Коэффициент габаритной яркости как количественный параметр для оценки отражательных свойств объемных тел // ОМП. 1974. № 1. С. 8 — 10.

5. Якушенков Ю. Г. Теория и расчет оптико-электронных приборов. М.: Логос, 1999. 480 с.

6. Тевяшов В. И., Тымкул В. М., Шуба Ю. А, Оптические характеристики отражения объемных тел в поляризованном свете // ОМП. 1979. № 10. С. 8—11.

7. Тымкул В. М., Шуба Ю. А. Рассеяние света диффузными квадратичными поверхностями в фотометрическом представлении // Там же. 1978. № 11. С. 11—13.

8. Тымкул В. М., Тымкул Л. В., Кудряшов К. В. Отражение оптического излучения телами с направленно– рассеивающим покрытием // Изв. вузов. Приборостроение. 2007. Т. 50, № 10. С. 58—63.

9. Схема МКС [Электронный ресурс]: .

10. [Электронный ресурс]: .

11. Методика наблюдений переменных звезд [Электронный ресурс]: .

Василий Михайлович Тымкул Любовь Васильевна Тымкул Юрий Александрович Фесько Кирилл Васильевич Кудряшов Светлана Валерьевна Луговских

Сведения об авторах — канд. техн. наук, профессор; Сибирская государственная геодезиче-
ская академия, кафедра наносистем и оптотехники, Новосибирск; E-mail: kaf.oep@ssga.ru — канд. техн. наук, профессор; Сибирская государственная геодезическая академия, кафедра наносистем и оптотехники, Новосибирск; E-mail: kaf.oep@ssga.ru — аспирант; Сибирская государственная геодезическая академия, кафедра наносистем и оптотехники, Новосибирск; E-mail: y.a.fesko@gmail.com — магистр; Сибирская государственная геодезическая академия, кафедра наносистем и оптотехники, Новосибирск; E-mail: kir-kudryashov@yandex.ru — бакалавр; Сибирская государственная геодезическая академия, кафедра наносистем и оптотехники, Новосибирск; E-mail: mumak2004@yandex.ru

Рекомендована кафедрой наносистем и оптотехники

Поступила в редакцию 08.02.13 г.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 5