Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
проектирование и конструирование / Raketno-kosmicheskaya_tekhnika_Mashinostroenie_En.pdf
Скачиваний:
1301
Добавлен:
09.03.2016
Размер:
14 Mб
Скачать

Аджян А.П., Аким Э.Л., Алифанов О.М., Андреев А.Н. Ракетно-космическая техника. Машиностроение. Энциклопедия. T. IV-22 В двух книгах. Книга первая

КОСМИЧЕСКИЙ МУСОР

441

 

 

При пиковой активности величины лив невых метеороидных потоков могут превысить фоновый поток в пять и более раз. Изредка очень значительные потоки (метеороидные штормы превышают фоновую активность в 10 000 и более раз) могут неожиданно возник нуть в короткие периоды времени (1…2 ч) в те годы, когда кометы источники проходят свои перигелии. Например, такой шторм был пред сказан для потока «Леонид» в ноябре 1998 и

гг. (17–20 ноября).

Считается, что метеороидные ливневые потоки состоят из относительно больших час тиц с массами более чем 1 10 8…1 10 6 г с низ кими плотностями (0,5…1,0 г/см3).

3.13.3.КОСМИЧЕСКИЙ МУСОР

Вотличие от микрометеороидов космиче4 ский мусор (КМ), имеющий искусственное про исхождение (отработавшие нефункционирую щие КА, отработавшие ступени ракет, целевые объекты, продукты, образовавшиеся в результате износа поверхности КА, а также фрагменты взрывов КА и ракетных ступеней), вращается вокруг Земли и будет оставаться на орбите до тех пор, пока из за атмосферного торможения и других возмущающих сил со временем их орби ты не опустятся к атмосфере. Так как атмосфер ное торможение ослабевает с увеличением вы соты орбиты, большие объекты КМ на орбитах 600 км и более могут оставаться на орбите в те чение десятков, тысяч и даже миллионов лет [2].

Понятие космического мусора

Международная академия астронавтики IAA (The International Academy of Astronautics) дает следующее определение техногенного «космического мусора» (в нотации ЕКА) или «орбитального мусора» (в нотации НАСА):

КМ — любые объекты, образовавшиеся в результате человеческой деятельности и нахо дящиеся на орбите вокруг Земли, которые яв ляются нефункционирующими с неприемле мым ожиданием начала выполнения или во зобновления предназначенных им функций или любых других функций, для которых они могут быть предназначены, включая фрагмен ты и части таких объектов [3].

Модель техногенной обстановки SDPA 2000

Принятая в РКА модель анализа и про гнозирования загрязнения ОКП техногенными объектами SDPA 2000, разработанная А.И. На заренко, является наиболее известной из отече

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

Модель SDPA 2000 — полуаналитическая стахостическая и служит для среднесрочного и долгосрочного прогнозирования техногенных частиц размером более 1 мм, построения про странственных распределений концентраций и характеристик скорости, а также оценки риска столкновений. На основе модели SDPA 2000 создана инженерная модель SDPA Е для бы строго и визуального представления характе ристик КМ [15, 16].

Вклад каждого объекта в концентрацию техногенных частиц в точке с геоцентрически ми координатами (r, Ι) определяется по фор муле [17]:

7(r,Ι)

1

 

 

,

 

 

 

 

 

2 3 ra (sin2 i sin2 Ι)(r r

p

)(r r)

 

 

 

a

 

 

 

 

(3.13.11)

где a (rp ra ) / 2 — большая полуось орбиты. Концентрация принимается не зависящей

от долготы.

Далее результаты расчетов для разных то чек околоземного космического пространства суммируются по объектам КМ.

В модели SDPA в качестве исходных дан ных рассматривают статистические нормиро ванные распределения трех элементов орбит: высоты перигея p(hp), эксцентриситета р(е) и наклонения р(i). Для функции концентрации выведена следующая формула:

7(r,Ι) N8F (Ι) 2 2r 2 h

9(hp ,e,r)Μ(hp ,e,r)p(hp)p(e) hp e, (3.13.12)

h p e

где N8 — общее число объектов;

F (Ι) :

 

p(i)di

 

при sini sinΙ; (3.13.13)

 

 

 

sin2 i sin2 Ι

i

 

 

9(hp , e, r) — нормированных в долях периода времени, в течение которого объект с элемента ми орбиты (hp , e) находится в высотном диапа зоне (r, r h);

 

(1 e)2

 

 

r

2

 

Μ(h , e, r)

 

 

 

 

&

 

)

. (3.13.14)

 

 

 

 

 

 

 

 

&

 

)

p

1

e2

 

 

 

 

% rp

(

 

При моделировании фрагментации в моде ли SDPA применяют зависимость среднего числа ежегодно образующихся частиц от их размеров,

Аджян А.П., Аким Э.Л., Алифанов О.М., Андреев А.Н. Ракетно-космическая техника. Машиностроение. Энциклопедия. T. IV-22 В двух книгах. Книга первая

442 Глава 3.13. ИССЛЕДОВАНИЯ РИСКА ПРОБОЯ КОСМИЧЕСКОГО АППАРАТА

которую принимают пропорциональной прирос ту числа каталогизированных объектов [17]:

n(d 6 dj ) k(dj )n(dcat ).

(3.13.15)

Значения коэффициента k(dj ) являются настраиваемыми параметрами модели.

Мгновенное значение плотности потока Q(t) равно произведению концентрации час тиц 7(t) на значение относительной скорости Vrel(t) в данной точке:

Q(t) 7(t)Vrel (t).

(3.13.16)

Значение относительной скорости опре деляется как

n

Vrel Vrel (Azj )[pVrel (Az j )], (3.13.17)

j 1

где pVrel (Azj ) — статистическое распределение направлений относительной скорости; Azj — угол между вектором относительной скорости и вектором скорости КА; Vrel (Azj ) — усредненная зависимость относительной скорости от ее на правления.

Модель техногенной среды 1991 г.

В модели НАСА [8] техногенная среда, как и метеороидная, описывается математиче ской формулой, определяющей плотность по тока частиц с размерами, большими рассмат риваемого, через единичную неориентирован ную площадку.

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

Средняя скорость перемещения осколков оценивается в 10 км/с. Структура распределе ния скоростей осколков различна в зависимо сти от высоты и наклонения. Средняя массо вая плотность космических осколков диамет ром до 1 см составляет 2,7 г/см3. Для более крупных объектов предлагается использовать меньшую плотность.

Основные положения. Внутри области 2 000 км, в пределах которой работает рассмат риваемая модель, по оценкам на основе данных на середину 1988 г. содержится 1,5…3 млн кг объектов человеческой деятельности на орби

те. Большинство из них находятся на орбитах с большим наклонением, где они движутся со средней скоростью 10 км/с.

Примерно 1 500 использованных ракет ных ступеней, пассивных полезных нагрузок и несколько активных полезных грузов состав ляют большую часть этой массы. Данные объ екты в настоящее время находятся под наблю дением Группы по космическому наблюдению США (US Space Command) наряду с 4 500 дру гими с суммарной массой 20 000 кг, большин ство из которых — фрагменты спутников или другие орбитальные металлические осколки. Недавние наблюдения показывают, что сум марная масса орбитальных с диамет ром 1 см и менее составляет около 1000 кг, а осколков менее 1 мм — около 300 кг. Такое распределение массы делает среду орбиталь ных осколков опаснее метеороидной по боль шинству вопросов, связанных с космическими кораблями на высотах ниже 2 000 км [8].

Поток орбитальных осколков на свободно вращающуюся поверхность

В осколочной модели 1991 г. задан в ана литическом виде поток осколков КМ за один год, имеющих значение диаметра более d через цилиндрическую поверхность, ось вращения, направленную по радиус вектору Земли, с пло щадью главного сечения 1 м2. Данный поток описывается выражением [18]:

F (d) 4H(d) (a, S) (i)[F1(d)g1(t) F2(d)g2(t)], (3.13.18)

где

 

 

 

 

 

lg( d ) 0,78

 

2

 

 

 

 

&

 

 

 

)

 

 

 

 

exp&

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

 

0,637

 

 

 

H(d) 10 %

 

 

 

 

 

 

( ; (3.13.19)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

S

1,5

 

 

 

 

 

(a, S)

10200

140

 

 

 

 

;

(3.13.20)

 

 

a

 

 

 

S

1,5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10200

140

1

 

 

 

 

 

 

 

 

 

 

(i) — коэффициент, заданный в таблич ном виде, графически он представлен на рис. 3.13.3, для наклонения i 51,6 (наклоне ние орбиты МКС) 1,028;

 

 

F (d) 122, 10 5 d 2,5 ;

(3.13.21)

 

 

 

1

 

 

 

 

F

(d) 81, 1010

(d 700) 6 ;

(3.13.22)

 

2

 

 

 

 

 

 

!

t 1988

при t + 2011;

 

102,

 

 

 

 

g1(t)

23

 

 

t 2011

 

(3.13.23)

 

104,

при t 2011;

 

#102,

 

 

 

 

 

g2(t) 1 0,05(t 1988),

(3.13.24)

Аджян А.П., Аким Э.Л., Алифанов О.М., Андреев А.Н. Ракетно-космическая техника. Машиностроение. Энциклопедия. T. IV-22 В двух книгах. Книга первая

КОСМИЧЕСКИЙ МУСОР

443

 

 

Рис. 3.13.3. Зависимость коэффициента y от наклонения орбиты

где t — календарный год; a — высота в км; S — индекс солнечной активности F10,7.

Коэффициент g1(t) показывает влияние технической политики на рост числа частиц КМ. До 2011 г. их ежегодный прирост по мо дели составляет 2 %, а 4 % — после 2011 г.

На рис. 3.13.2 проводится сравнение потока орбитальных осколков и потока ме теороидов по моделям НАСА для парамет ров h 500 км, t 1995 г., i 28,5 , S(t1) 97,0.

Усредненная форма и массовая плотность частиц КМ

Информация о форме обломков очень скудна. Реальная их форма является непра вильной (стержни, пластины и т.д.). При раз работке принята сферическая форма объектов Плотность орбитальных осколков принимается 1993 г. в виде функции, завися щей от диаметра частицы d в см:

!7 2,8d 0,74

г /см3 при d 0,62;

 

 

 

при d+ 0,62.

#

7 4 г /см3

Распределение по скорости и направлению.

Усредненное по всем высотам ненормирован ное распределение (функция плотности веро

ятности) относительной скорости КА и части цы КМ описывается [18] следующим уравне нием:

 

 

 

!

 

V AV 0

 

2

 

V DV 0

 

2

 

 

 

 

 

 

 

&

 

)

 

 

&

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

)

 

 

&

 

)

 

 

g(V ) (2VV

 

 

 

% BV0

(

% EV 0

(

 

0

V 2) Ge

 

 

 

Fe

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

#

 

 

 

 

 

 

 

 

 

 

 

 

 

 

HC(4VV0 V 2),

 

(3.13.25)

где V — относительная скорость, км/с;

 

 

 

 

!7,25 0,015(i 30) при

i + 60;

 

 

 

 

 

 

 

V0

 

 

i 6 60;

 

 

 

 

 

 

 

 

 

 

 

#7,7 при

 

 

 

 

 

 

 

 

 

 

 

!18,7 при

i + 60;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

60+ i + 80;

 

 

 

 

G 18,7 0,0289(i 60)3 при

 

 

 

 

 

250,0 при

i 6 80;

 

 

 

 

 

 

 

 

 

#

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

!0,5 при

 

i + 60;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A 2,5; B 0,5 0,01(i 60) при 60+ i + 80;

 

 

 

 

 

0,3 при

 

i 6 80;

 

 

 

 

 

 

 

 

 

 

 

#

 

 

 

 

 

 

 

 

 

 

 

 

 

!0,3 0,0008(i 50)2 при

i + 50;

 

 

 

 

 

 

 

 

 

 

50) при 50+ i + 80;

 

 

 

 

 

 

F 0,3 0,01(i

 

 

 

 

 

 

0,0 при

i 6 80;

 

 

 

 

 

 

 

 

 

 

 

 

#

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D 1,3 0,01(i 30); E 0,55 0,005(i 30); H 1,0 0,0000757(i 60)2;

Аджян А.П., Аким Э.Л., Алифанов О.М., Андреев А.Н. Ракетно-космическая техника. Машиностроение. Энциклопедия. T. IV-22 В двух книгах. Книга первая

444 Глава 3.13. ИССЛЕДОВАНИЯ РИСКА ПРОБОЯ КОСМИЧЕСКОГО АППАРАТА

!0,0125

при i + 100;

 

C

0,00125(i

100) при i 6100;

#0,0125

i — наклонение орбиты КА, град. Вероятность того, что с КА частица КМ

будет сталкиваться со скоростями, значения которых лежат в диапазоне [V1; V2], описыва ется нормировочным выражением

V 2

1

P : g(V )dV

: g(V )dV . (3.13.26)

V 1

0

Число столкновений КА с орбитальными осколками, летящими из заданного направле ния, можно оценить, используя распределение по скорости g(V). За направление столкнове ния принимается пересечение вектора скоро сти КА и выбранной круговой орбиты орби тального осколка, т.е. векторы относительной скорости могут быть получены векторным сложением в тангенциальной по отношению к Земле плоскости.

Если скорость КА равна 7,7 км/с, то на правление вектора относительной скорости за дается соотношением

cos (/ ) V/V V , (3.13.27)

max

15,4

 

где — угол между вектором скорости столкно вения и вектором скорости КА в системе коор динат КА; Vmax — максимально возможная ско рость столкновения КА с КМ.

Модель техногенной обстановки ORDEM 2000

Разработанная НАСА модель техноген ной среды ORDEM 2000 основана на точной эмпирической оценке популяций частиц КМ с орбитами до 2 000 км. Основными источ никами для создания модели служили лог наблюдаемых орбитальных объектов

Space Command Catalog и данные, получен

ные телескопом Хайстэк (Haystack Radar) и полученные после изучения поверхности возвращенного с орбиты модуля длительного экспонирования LDEF. Данные использова лись для оценки распределений по разме рам, по наклонениям, по высоте и по экс центриситету. Минимальный размер частиц, рассматриваемый моделью 1 10 3 мм. В даль

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

(Haystack Auxiliary Radar) и Goldstone, а так же с возвращенной солнечной панели теле скопа Хаббл (Hubble).

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

Модель обстановки MASTER ЕКА

Модель MASTER описывает техноген ную и природную обстановку вокруг Земли в виде потока частиц с диаметрами более 1 мкм, падающего на определенный пользо вателем орбитальный объект. В модели учте ны все возможные источники появления КМ: отработавшие полезные нагрузки и верхние ступени, фрагменты столкновений и взрывов на орбите, пыль и шлак от продук тов сгорания твердотопливных двигателей, частицы теплоносителя из калий натриевой смеси спутников «Космос», на которых ис пользовался ядерный реактор, кусочки крас ки деградирующих поверхностей, а также выбросы.

При описании техногенной обстановки в модели MASTER орбиты по высоте разбива

3.13.2. Число осколков на всех орбитах в зависимости от размера по модели 1996 г. EKA MASTER (для GEO i 15 )

Размер

Общее число

Число объектов

Число объектов

Число объектов

осколка

объектов

на LEO

на MEO

на

GEO

 

 

 

 

 

 

6 0,1 мм

1,2 1011

1,0 1010 (8,3 %)

7,8 1010 (65,0 %)

9,8 108

(0,8 %)

6 1 мм

1,8 108

3,3 107 (18,3 %)

1,0 108 (55,6 %)

1,3 106

(0,7 %)

6 1 см

363 000

160 000 (44,0 %)

153 000 (42,1 %)

4 600 (1,3 %)

 

 

 

 

 

6 10 см

18 900

15 100 (80,0 %)

2 218 (11,7 %)

950 (5,0 %)