Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

VVS-LabRabota-06-TrassyPoletaKA

.pdf
Скачиваний:
71
Добавлен:
16.03.2015
Размер:
793.77 Кб
Скачать

Курс «Введение с специальную технику

ЛАБОРАТОРНАЯ РАБОТА №6

Орбиты КА в трехмерном пространстве. Трассы полета КА.

ТЕОРЕТИЧЕСКАЯ ЧАСТЬ

1)Параметры орбиты КА в трехмерном пространстве:

Впредыдущей лабораторной работе рассматривались орбиты КА в двухмерном пространстве. Если рассмотреть движение КА по орбите относительно Земли в трехмерном пространстве, то к ранее указанным параметрам орбиты добавятся еще несколько (см. рис. 1).

Рисунок 1 – Параметры орбиты КА в трехмерном пространстве

Из рисунка 1 видно,

что в плоскости орбиты КА ее форма также описывается с помощью

 

 

 

радиус векторов перигея

( r ) и

апогея ( r ). Угловое положение КА в плоскости орбиты

определяется по углу истинной аномалии .

Кроме того для описания параметров орбиты КА в трехмерном пространстве вводятся следующие дополнительные параметры:

X г.о.Yг.о. Zг.о. - геоцентрическая прямоугольная система координат;

- угол долготы восходящего узла;

i - угол наклонения плоскости орбиты КА к плоскости экватора (наклонение орбиты);

- аргумент перигея (угол между линией пересечения плоскостей экватора и орбиты и радиус-вектором перигея орбиты, см. рис. 1).

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

 

u .

(1)

 

2) Расчет приращения скорости для изменения плоскости орбиты КА:

 

Во многих случаях требуется не только изменение плоскостных параметров орбиты КА

 

 

( r , r ), но также

изменения угла наклонения плоскости орбиты к экватору Земли. К примеру,

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

космодрома Байконур наклонение орбит составляет 51,6° ( i 51,60 ), что связано с районами падений отработавших ступеней ракет-носителей.

1

Курс «Введение с специальную технику

Изменение угла наклона плоскости орбиты проводится путем изменения направления вектора скорости разгонного блока КА. Для этого необходимо сообщить КА приращение импульса скорости. Схема для расчета приращения характеристической скорости, необходимой для изменения плоскости орбиты на угол , представлена на рисунке 2.

Рисунок 2 – Схема для расчета приращения скорости КА, необходимой для изменения плоскости орбиты

 

На рисунке 2 указаны следующие векторы:

 

-

первоначальный вектор скорости КА,

 

V1

- вектор скорости КА после изменения плоскости орбиты на угол

i . Считаем,

что модуль

V2

скорости КА не изменяется (V1 V2 ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i (чтобы

 

Из схемы видно,

что для изменения направления вектора скорости V

на угол

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

получить вектор скорости

V2

) необходимо добавить вектор скорости

V3 .

Из геометрических

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

соотношений нетрудно получить модуль вектора V3 :

 

 

 

 

 

 

 

 

 

V3

 

2 V1

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для спутников, запускаемых с космодрома Байконур с углом наклона плоскости орбиты

i 51,60 , приращение скорости для перевода КА в плоскость экватора составит:

 

 

 

 

 

 

 

 

 

i

 

 

 

51,60

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V 2 V sin

 

 

2 V sin

 

 

 

 

 

0,87 V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

1

 

2

 

1

 

 

2

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим пример перелета КА с низкой опорной круговой орбиты с наклонением i1 на

высокую круговую орбиту с наклонением i2 .

 

 

 

 

 

 

 

 

 

 

 

 

Пример 1. Требуется определить оптимальное суммарное приращение скорости КА V ,

требуемое для перевода его с низкой опорной орбиты с высотой

Hкр1 200 км и наклонением

i

51,60

на геостационарную орбиту (

H

кр 2

36000

км, i

2

00 ).

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим три варианта схем выведения:

 

 

 

 

 

 

 

 

 

 

1) Первый импульс

V1

поворачивает плоскость низкой опорной круговой орбиты высотой

 

H

кр1

200

км на угол

i i

i

51,60

00 51,60 ,

второй импульс V

уже в плоскости

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

экватора переводит КА на эллиптическую орбиту с высотой перигея H

200 км и высотой

 

апогея H

36000 км, третий

импульс

 

V3

переводит КА

с эллиптической

орбиты на

 

геостационарную орбиту.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2) Первый импульс

V1

переводит

 

КА

 

на

эллиптическую

орбиту с

 

высотой перигея

 

H

 

200 км и

высотой

апогея

H

 

36000 км (наклонение

текущей

орбиты

i 51,60 ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

второй импульс V2 переводит КА с эллиптической орбиты на круговую орбиту с высотой

 

H

кр 2

36000 км (наклонение текущей орбиты

i 51,60 ), третий импульс

V

поворачивает

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

3

 

 

 

плоскость круговой орбиты высотой

 

H

кр 2

36000 км на угол i i

i 51,60

00 51,60 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Курс «Введение с специальную технику

3) Первый

импульс

 

V1

переводит

 

КА

 

на

эллиптическую

орбиту с

высотой перигея

H

 

200 км

и

 

высотой

апогея H

 

36000 км

(наклонение

текущей

орбиты i 51,60 ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

второй импульс V2

в точке апогея поворачивает плоскость эллиптической орбиты на угол

i i i

2

51,60 00

51,60 ,

третий импульс V

переводит КА с эллиптической орбиты с

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

наклонением i

2

 

00

на геостационарную орбиту.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Схема 1:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассчитываем скорость КА на опорной круговой орбите, воспользовавшись

выражениями из прошлой лабораторной работы:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

V

 

 

 

RЗ

 

7,91

6371

 

7,790 км/с.

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

кр1

 

 

 

 

rкр1

 

 

 

 

6371 200

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассчитываем первый импульс скорости V1 :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

51,60

 

 

 

 

 

 

V 2 V

кр1

sin

 

2 7,790 sin

 

 

 

6,777 км/с.

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассчитываем большую полуось эллиптической орбиты:

 

 

a

r r

 

6371 200 6371 36000

24471 км.

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассчитываем скорость в перигее эллиптической орбиты V :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

1

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

V VI

RЗ

 

 

 

 

7,910

6371

 

 

 

 

 

 

 

10,25

км/с.

 

 

 

 

r

 

a

 

6571

 

24471

 

 

 

Рассчитываем второй импульс скорости V2 :

V2 V Vкр1 10,25 7,79 2,46 км/с.

Рассчитываем скорость в апогее эллиптической орбиты V :

 

 

2

 

1

 

 

 

 

2

 

 

1

 

 

 

 

 

 

 

 

 

 

 

V VI

 

 

 

7,91

6371

 

 

 

 

 

1,589 км/с.

RЗ

 

 

 

 

 

 

 

r

 

a

 

 

 

42371

 

24471

 

 

Рассчитываем скорость Vкр 2

на круговой орбите высотой Hкр 2 36000 км:

V

V

 

 

 

RЗ

 

 

7,91

 

 

6371

 

3,067 км/с.

I

 

 

 

 

 

кр 2

 

 

 

r

 

42371

 

 

 

 

 

 

 

 

Рассчитываем третий импульс скорости V3 :

V3

Vкр 2

V

3,067 1,589 1,478 км/с.

Определяем суммарный импульс скорости V :

V V1 V2

V3

6,777 2,460 1,478 10,715 км/с.

Схема 2:

Рассчитываем скорость КА на опорной круговой орбите:

V

 

V

 

 

 

RЗ

 

 

7,91

 

6371

 

 

7,790 км/с.

 

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

кр1

 

 

 

 

rкр1

 

 

 

 

 

 

 

 

6371 200

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассчитываем большую полуось эллиптической орбиты:

 

 

a

r r

 

 

6371 200 6371 36000

24471 км.

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассчитываем скорость в перигее эллиптической орбиты V :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

V

VI

 

 

 

RЗ

 

 

 

 

 

7,910

6371

 

 

 

 

 

 

 

10,25

км/с.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

a

 

 

 

 

 

6571

 

24471

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

Курс «Введение с специальную технику

Рассчитываем первый импульс скорости V1 :

V1 V Vкр1 10,25 7,79 2,46 км/с.

Рассчитываем скорость в апогее эллиптической орбиты V :

 

 

 

 

 

2

 

 

1

 

 

 

 

 

 

 

2

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V VI

 

RЗ

 

 

 

 

 

 

7,91 6371

 

 

 

 

 

 

1,589

км/с.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

a

 

 

 

 

 

 

42371

 

24471

 

 

Рассчитываем скорость Vкр 2 на круговой орбите высотой Hкр 2 36000 км:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

V

 

 

 

RЗ

 

7,91

 

 

6371

 

3,067 км/с.

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

 

кр 2

 

 

 

r

 

 

 

 

42371

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассчитываем второй импульс скорости V2 :

 

 

 

 

 

V2

Vкр 2

V 3,067 1,589 1,478км/с.

 

 

 

 

 

 

Рассчитываем третий импульс скорости V3 :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

51,60

 

 

 

 

 

 

 

V 2 V

 

 

sin

 

2 3,067 sin

 

 

 

2,668 км/с.

 

 

 

 

 

 

3

 

 

кр 2

 

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определяем суммарный импульс скорости V :

V V1 V2 V3 2,668 2,460 1,478 6,606 км/с.

Схема 2:

Рассчитываем скорость КА на опорной круговой орбите:

V

 

V

 

 

 

RЗ

 

 

7,91

 

6371

 

 

7,790 км/с.

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

 

 

кр1

 

 

 

 

rкр1

 

 

 

 

 

 

 

 

6371 200

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассчитываем большую полуось эллиптической орбиты:

 

a

r r

 

 

6371 200 6371 36000

24471 км.

 

 

2

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассчитываем скорость в перигее эллиптической орбиты V :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

V

VI

 

 

 

RЗ

 

 

 

 

 

7,910

6371

 

 

 

 

 

 

 

10,25 км/с.

 

 

 

 

 

 

 

 

 

r

 

a

 

 

 

 

 

6571

 

24471

 

 

Рассчитываем первый импульс скорости V1 :

V1 V Vкр1 10,25 7,79 2,46 км/с.

Рассчитываем скорость в апогее эллиптической орбиты V :

 

 

2

 

1

 

 

 

2

 

 

1

 

 

 

 

 

 

 

 

 

 

V VI

 

 

 

7,91

6371

 

 

 

 

 

1,589 км/с.

RЗ

 

 

 

 

 

 

 

r

 

a

 

 

42371

 

24471

 

 

Рассчитываем второй импульс скорости V2 :

 

 

 

 

 

 

 

 

 

i

 

 

 

51,60

 

 

V

2 V

sin

2

1,589 sin

 

1,390

км/с.

 

2

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассчитываем скорость Vкр 2 на круговой орбите высотой Hкр 2 36000 км:

 

 

 

 

 

 

 

 

 

 

 

 

 

V

V

 

 

RЗ

 

7,91

6371

 

3,067 км/с.

 

I

 

 

 

кр 2

 

 

 

r

 

 

42371

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассчитываем третий импульс скорости V3 :

V3 Vкр 2 V 3,067 1,589 1,478 км/с.

Определяем суммарный импульс скорости V :

V V1 V2 V3 1,390 2,460 1,478 5,328 км/с.

4

Курс «Введение с специальную технику

Из расчетов видно, что наименее затратным вариантом выведения КА на геостационарную орбиту является схема 3 с V 5,328 км/с.

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

На рисунке 3 приведена оптимальная схема вывода КА на геостационарную орбиту с космодрома Байконур.

Рисунок 3 – Оптимальная схема выведения КА на геостационарную орбиту с космодрома Байконур

Пример 2. Рассчитайте начальную массу КА M

нач , который должен совершить перелет КА

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

КА

с низкой опорной орбиты с высотой

 

H

кр1

200 км и наклонением i 51,60 на геостационарную

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

орбиту ( H

кр 2

36000 км, i

2

00 ). Считаем, что на геостационарную орбиту требуется доставить

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

полезный

груз

 

 

массой

 

M

кон 1000 кг. Конструктивная характеристика разгонного блока КА

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

КА

 

 

 

 

 

 

 

 

 

s 12 . Скорость истечения двигательной установки КА 3000 м/с.

Решение:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В примере 1 мы определили, что суммарное приращение скорости КА V 5,328 км/с.

По формуле из предыдущих лабораторных работ определяем массу топлива M Т :

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

5328

 

 

 

 

 

 

 

 

e

 

1 s 1

e

 

 

1

12 1

 

 

 

 

 

 

 

 

 

 

 

3000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M T

 

 

 

 

 

 

 

 

 

 

 

 

 

M КАкон

 

 

 

 

 

 

 

 

1000 8751кг.

 

 

 

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

5328

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s e

 

 

 

 

 

 

 

 

 

12 e 3000

 

M

 

 

 

s MT

 

 

12 8751

9546 кг.

 

 

 

 

РБ

 

 

 

 

 

 

 

 

 

 

 

s

1

 

 

 

 

12 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M нач M

кон

M

РБ

1000 9546 10546 кг.

 

 

КА

 

 

 

 

KA

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3) Построение трассы полета КА:

 

 

 

 

Трассы полета характеризуют

 

движение КА

относительно поверхности вращающейся

Земли. Трасса полета показывает движение подспутниковой точки О (см. рис. 4) по поверхности

Земли в течении

полета. Трасса полета строится в системе координат

 

, где

-

географическая

долгота положения подспутниковой точки (изменяется от

0° до

360°),

-

географическая широта положения подспутниковой точки (может изменяться в пределах от

минус 90° до +90°).

5

Курс «Введение с специальную технику

Рисунок 4 – Подспутниковая точка и линия трассы КА

Пример построения трассы полета КА приведен на рисунке 5.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рисунок 5 – Трасса полета КА

 

В общем виде система уравнений трассы записывается в виде:

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

arcsin sin i sin

T

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

arctg

 

 

 

tg

2

 

 

 

 

 

t

 

t

 

 

 

 

cos i

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

З

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

T

 

t t

н

, t

к

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

- географическая широта подспутниковой точки в момент времени t ;

- географическая долгота подспутниковой точки в момент времени t ;

i- наклонение орбиты спутника (может быть от 0° до 90°);

T - период обращения КА;

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Курс «Введение с специальную технику

З

- угловая скорость вращения Земли ( З

 

2

 

 

 

рад/сек.);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

24 * 3600

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- начальное значение долготы подспутниковой точки;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- изменение долготы восходящего узла орбиты за один виток;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t - текущее время, отсчитываемое от момента начала движения.

 

 

 

 

 

 

 

 

Пример 3. Составить алгоритм построения трассы КА на одном витке с периодом

обращения T с помощью программы MS Exscel.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как в уравнениях трассы КА присутствуют функции

tg и

arctg , то при угловых

положениях КА на орбите

и

3

возникают участки разрыва математической функции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t . Поэтому уравнения трассы удобно представить в виде системы:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t 0, T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

Если

0

 

 

 

 

t

 

 

 

 

, то arctg cos i tg

 

 

t

 

 

 

t t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

2

 

 

 

 

 

 

 

 

T

 

 

 

 

З

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

3

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если

 

 

 

 

 

 

 

 

t

 

 

 

 

, то

arctg cos i tg

 

 

 

t

 

З t t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

T

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

T

 

 

 

 

 

 

3

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

t

Если

 

 

 

 

 

 

 

 

 

 

 

t

2

, то 2 arctg cos i tg

 

 

 

 

t З

t t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

arcsin sin i sin

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Считаем, что

 

 

0

,

0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Алгоритм построения графика трассы КА на одном витке будет следующим:

1)Задаем исходные значения параметров орбиты КА ( r , r , i )

2)Рассчитываем угловую скорость вращения Земли:

 

 

 

2

7,269 10 5

рад/сек.

 

 

 

З

 

 

 

 

 

 

 

 

 

 

24 * 3600

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3) Рассчитываем период обращения КА:

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

T 2

a 2

, где a

r

 

 

3,986

105

км32

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

З

 

З

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4) Задаем числовой ряд из N точек, то есть:

 

 

 

 

 

 

 

N 101, где

n 0..100 ( n 1,2,3,...,98,99,100 ).

 

5)Задаем зависимость текущего времени на витке от номера текущей точки графика: t 100T n ,

6)После выполнения пунктов 1-5 в MS Exscel должна быть создана форма таблицы как на рисунке 6.

7

Курс «Введение с специальную технику

Рисунок 6 – Форма таблицы для построения трассы полета КА

7) Программируем столбик значений долготы при 0 t 2 :

T 2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

arctg cos i tg

 

 

t

З t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8) Программируем столбик значений долготы при

 

t

2

 

 

3

:

2

 

T

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

arctg cos i tg

 

 

 

 

t

З t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9) Программируем столбик значений долготы при

 

3

 

t

 

2

 

2 :

2

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 arctg cos i tg

 

 

t

З t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10) С помощью логической функции «ЕСЛИ» программируем суммарное значение текущей долготы КА.

11) Программируем столбик текущего значения широты КА в течение полета на витке:

 

 

2

 

arcsin sin i sin

 

t .

 

 

T

 

11)Переводим значения долготы и широты КА из радиан в градусы, и оформляем в MS Exscel в виде крайних правых двух столбцов.

12)Строим точечную диаграмму, где значения аргумента – это долгота КА, а значения функции – это широта КА (см. рис. 7)

Рисунок 7 – Трасса полета КА на одном витке

8

Курс «Введение с специальную технику

КОНТРОЛЬНЫЕ ВОПРОСЫ

1.Перечислите параметры орбиты КА в трехмерном пространстве.

2.Объясните, как можно поменять наклонение орбиты КА.

3.Объясните, какая наиболее оптимальная схема перелета между двумя орбитами с разными наклонениями.

4.С какой широты выгоднее всего запускать геостационарные спутники? Поясните ответ.

5.Что такое географическая система координат?

6.Что такое трасса полета КА?

 

 

 

 

 

 

 

ЗАДАНИЕ

 

 

 

 

 

 

 

1) КА расположен на стартовой орбите с высотой перигея

H 1 , высотой апогея

H 1 и

наклонением орбиты i1 .

Требуется перевести КА на новую орбиту с высотой перигея

 

H 2 ,

высотой апогея H 2

и наклонением орбиты i2 .

 

 

 

 

 

 

 

 

Рассчитаете

оптимальное суммарное приращение скорости V , требуемое КА

для

перелета с начальной орбиты на конечную орбиту.

 

 

 

 

 

 

 

 

2) Постройте в MS Excel трассу полета КА на одном витке на орбите с высотой перигея

H 2 ,

высотой апогея H 2

и наклонением орбиты i2 .

 

 

 

 

 

 

 

 

Исходные данные для расчетов приведены в таблице 1.

 

 

 

 

 

 

 

Таблица 1 – Исходные данные для расчетов

 

 

 

 

 

 

 

 

H 1 , км

 

H 1 , км

i1 , град

H 2 , км

H 2 , км

i2 , град

 

 

 

варианта

 

 

 

 

 

 

 

 

 

 

 

 

1

 

200

 

 

200

30

3000

3500

35

 

 

 

 

2

 

250

 

 

250

60

5000

5000

10

 

 

 

 

3

 

200

 

 

250

56

400

10000

60

 

 

 

 

4

 

210

 

 

220

45

36000

36000

0

 

 

 

 

5

 

205

 

 

210

56

10000

20000

74

 

 

 

 

6

 

200

 

 

300

70

400

400

60

 

 

 

 

7

 

190

 

 

200

70

500

700

86

 

 

 

 

8

 

200

 

 

200

70

600

700

87

 

 

 

 

9

 

200

 

 

200

71

360

500

60

 

 

 

 

10

 

200

 

 

200

71

405

410

87

 

 

 

 

11

 

300

 

 

300

2

36000

36000

0

 

 

 

 

12

 

250

 

 

300

10

1000

1100

71

 

 

 

 

13

 

250

 

 

300

71

670

670

74

 

 

 

 

14

 

200

 

 

250

56

720

720

60

 

 

 

 

15

 

200

 

 

250

10

300

20000

15

 

 

 

 

16

 

210

 

 

220

56

500

30000

71

 

 

 

 

17

 

220

 

 

220

56

400

36000

0

 

 

 

 

18

 

230

 

 

250

56

300

600

74

 

 

 

 

19

 

230

 

 

240

71

400

2000

72

 

 

 

 

20

 

300

 

 

300

56

300

2000

74

 

 

 

 

21

 

300

 

 

310

71

1500

1500

74

 

 

 

 

22

 

310

 

 

315

56

1500

2000

40

 

 

 

 

23

 

200

 

 

400

56

310

600

64

 

 

 

 

 

 

 

 

 

 

9