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

книги из ГПНТБ / Балахонцев Б.Г. Сближение в космосе

.pdf
Скачиваний:
183
Добавлен:
23.10.2023
Размер:
7.74 Mб
Скачать

чивается переход из точки M в точку N за данное

вре­

мя т, можно принять значения

Гі = гм, r2

= rN, ф и т.

При

использовании этих величин

в качестве

исходных

дан­

ных для решения задачи сближения мы тем самым сво­ дим ее к решению задачи определения орбиты КА, про­ ходящей через две заданные точки пространства M и N, причем время движения от точки M до точки N должно быть равно т.

§ 3.6. О П Р Е Д Е Л Е Н И Е О Р Б И Т Ы

КА ПО

Д В У М

З А Д А Н Н Ы М

П О Л О Ж Е Н И Я М

 

 

В общем случае переход

из точки M в точку

/V мо­

жет быть произведен по гиперболической,

параболиче-

Рис. 3.11.

В о з м о ж н ы е орбиты

пере­

хода

м е ж д у точками M и

N

ской или эллиптической орбите. Тип орбиты опреде­ ляется взаимным положением точек M и N, а также временем перехода т. Однако при рассмотрении задачи сближения с орбиты ожидания в условиях ближнего космоса наиболее вероятными орбитами перехода будут эллиптические вследствие того, что переход по парабо­ лическим или гиперболическим орбитам потребует боль­ шого импульса скорости АѴ\, реализация которого бу­ дет практически невозможна. Поэтому все дальнейшее изложение мы будем проводить для эллиптических ор­ бит сближения и наша задача будет заключаться в оп-

90

ределении параметров эллипса, проходящего через точ­ ки M и N, положение которых относительно одного из фокусов F известно. Если известно и положение второ­ го фокуса F*, то задача определения траектории, про­ ходящей через точки M и N, решается однозначно.

Очевидно, что положение второго фокуса нельзя вы­ бирать произвольно, поэтому предварительно найдем геометрическое место вторых фокусов всех эллипсов, по которым возможен переход из точки M в точку N без ограничения на время перехода.

Если через гх и г*2 (рис. 3.11) обозначить расстоя­ ния от точек M и ІѴ до второго фокуса, то, используя свойство эллипса, можем записать следующее:

 

r\ =

2a~rl;

r; =

2a-r2.

 

(3.71)

Из формул (3.71) следует, что

если

большая

полу­

ось эллипса

равна

а, то

второй

 

фокус

орбиты

F* яв­

ляется точкой пересечения двух окружностей с

центра­

ми в точках

M и N

и радиусами

2

а г{

и 2 а г2.

Введем в рассмотрение величину S как расстояние между точками M и N:

S =

+

2 / у г cos <!>.

(3.72)

В зависимости от значения 2 а возможны следующие случаи [38].

1) г\ - fr* < S. Из выражения (3.71) следует, что это равносильно условию

2a<±-(r1

+ r, + S).

(3.73)

В этом случае окружности не пересекаются и пере­ ход из точки M в точку N по эллиптической орбите не­ возможен.

2) г* + r'2 — S. Это равносильно условию

2a=\(rl

+r2

+ S).

(3.74)

В этом случае окружности

касаются

и второй фо­

кус F* орбиты находится

на прямой MN.

Очевидно, что

91

это наименьшее

значение большой

полуоси а = ат,

при

котором эллиптическая

орбита еще

возможна.

 

3) r\ + r'2>

S. Это

равносильно условию

 

 

2a>T(rl

+ ri +

S).

(3.75)

В этом случае окружности пересекаются в двух точ­ ках F* и F*. Следовательно, существуют две эллипти­ ческие орбиты, связывающие точки M и N. Эти орбиты имеют одну и ту же большую полуось, а их вторые фо­ кусы располагаются симметрично относительно пря­ мой MN.

Как

видно из рис. 3.11, переход из точки

M

в точ­

ку N по эллиптической орбите с заданным

значением

большой полуоси

может быть

произведен по

 

четырем

различным траекториям: MNXN,

MN2N, MMXN

и

MM2N.

Если же направление обхода фокуса F задано, то оста­

ются

только две

возможные

траектории,

например

MNXN и MMiN.

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

Из рис. 3.11 также видно, что эллиптическая орбита MM\N характеризуется тем, что ее второй фокус F* ле­ жит вне эллиптического сегмента, образуемого прямой

MN и

рассматриваемой орбитой. Подобные орбиты бу­

дем

называть

э л л и п т и ч е с к и м и

о р б и т а м и

(траекториями)

п е р в о г о

р о д а

[38].

Если второй

фокус

F*, отвечающий орбите

MN\N,

лежит внутри со­

ответствующего эллиптического сегмента, то такую ор­

биту

будем

называть э л л и п т и ч е с к о й

о р б и т о й

(траекторией)

в т о р о г о

р о д а . Если

величина

боль­

шой

полуоси

определяется

из условия

формулы

(3.74),

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

92

Рассмотрим теперь характер изменения величины и направления скорости в точке М, при которой обеспечи­ вается переход из точки M в точку N по эллиптической орбите без учета времени. Связь между величиной ско­ рости Ѵм и ее направлением в точке M и заданными величинами г,, г2, ф дается уравнением семейства тра­ екторий одинаковой дальности [27]. Это уравнение после незначительных преобразований можно записать в виде

 

1/2 =

 

 

 

« 8 ( 1 - с о 8 ф )

 

 

(

3

7

6 )

 

 

м

Г] cos

Ѳ м

[sin

ф sin

ÖД І

+ (k — cos ф)

cos 9^]

'

v

"

'

где Ѳ— угол

между вектором

скорости

и

плоскостью

местного горизонта

в точке

M;

k — отношение

rjr2.

вве­

Для

проведения

дальнейшего

анализа

 

удобно

 

сти безразмерные

величины:

 

 

 

 

 

 

 

 

 

 

 

<• —V-.г

•—

V

'

 

 

 

(3.77)

 

 

 

 

 

Мк

 

ѵ

Мк

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Mr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

к

 

-И к

 

 

 

 

 

 

где

ѴМк

— круговая

скорость

на

орбите

радиуса

 

Г\\

ѴМп

ѴМгі—радиальная

 

и

трансверсальная

составляю­

щие вектора скорости в точке М.

 

 

 

 

 

 

 

Подставив значения выражений (3.77) в формулу

(3.76),

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

(A COS-TO ï2+

( и п ф ) б ч

— О -

c o s t ) ^ . = 0.

(3.78)

Уравнение

(3.78)

представляет

собой

уравнение кри­

вой

второго порядка с инвариантами:

 

 

 

 

 

 

 

 

 

Д =

-

L (1 _

cos ф) sin2 Ф >

О,

 

 

 

 

 

 

8 =

— \

sin2 ф <

0 при

Ф ф пк

(я =

0,1,...)-

 

(3-79)

Таким образом, для перехода из точки M в точку N по эллиптической орбите необходимо, чтобы конец век­ тора скорости в точке M лежал на гиперболе, задавае­ мой уравнением (3.78). Легко показать, что асимптота­ ми гиперболы являются ось ѵ\ и линия, соединяющая точки M и N. Отсюда следует, что гипербола «выпрям-

93

ляется»

при приближении

значения

угла

ф к

180° при

заданном k и с увеличением k

при фиксированном ф.

При

ф = «іс инварианты

Д

= 0 и

8 = 0,

а

уравнение

(3.78) дает две прямые. Так как случай ф = 2 п ъ (n = 0,ï,...) практического интереса не представляет, то определим

положение этих

прямых при ф = ( 2 / г + 1 ) і т

(я = 0, 1,

. . . ) .

Из формулы (3.79)

следует, что

при

ф=(2гс+1)тс

(п = 0,1,...)

 

 

 

 

 

Знак перед

корнем

определяется

направлением

об­

хода притягивающего центра. Поскольку мы условились рассматривать движение только по ходу часовой стрел­

ки, то из двух ветвей

гиперболы при ф Ф пъ и двух

пря­

мых при ф = ( 2 п + 1 ) т с

будем выбирать только те,

для

которых S>0.

Как известно, параболическая скорость Ѵп в точке,

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

притягивающего центра,

и круговая

скорость

Ѵ„ в этой

же точке связаны

соот­

ношением

Ѵп Ѵ^ 2

Ѵк

или

в безразмерных

величинах

в соответствии

с формулой (3.77)

 

 

 

 

 

 

 

Vn = VY.

 

 

 

(3.81)

Таким

образом, для

обеспечения

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

ор­

биты перехода

между

точками M и N величины

£ и к]

должны удовлетворять

условию

 

 

 

 

 

 

 

Ъ +

<

2.

 

 

(3.82)

Как видно из выражения (3.80), для перехода

из

точки M в точку

по

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

при

ф =

= (2п + 1)іг достаточно

обеспечить

значение

трансвер-

сальной скорости £ =

| / ] - q r ^

. тогда

как значение

ради­

альной скорости может быть любым, но удовлетворяю­ щим условию

которое следует из условия (3.82).

94

Введем новые безразмерные координаты £і, тц и при­ ведем уравнение (3.78) к каноническому виду. В ре­ зультате получим

«2

1.

(3.84)

 

В этом уравнении

-V

2

(

1 - cos 41)

2 ( 1

+

kS — k)

 

 

+

(k — COS ip)

v .

+

kS + k

} (3.85)

 

2

(

1 - COS <|>)

 

 

 

 

 

ftS

(ft — cos <|0

 

 

 

)

Перейдя в выражении (3.74) к безразмерным вели­ чинам, для граничной эллиптической орбиты можем за­ писать

 

 

г,

\

+k + kS

 

(3.86)

 

 

 

4 k

 

 

 

 

 

 

 

 

Затем, используя

зависимость

 

 

1 у 2

=

^

( 2 _

^ )

(

(3.87)

найдем

 

 

 

 

 

 

 

72

,

~2

_

2 ( l +

kS—k)

(3.88)

ïrp

Т

^гр

 

 

Z

'

 

 

 

 

1 +

kS + k

 

Из сравнения выражений (3.85) и (3.88) видно, что большая полуось годографа скорости равна той скоро­ сти, которую должен иметь КА в точке М, чтобы пе­ рейти из точки M в точку N по граничной эллиптиче­ ской орбите (рис. 3.12). Отсюда также следует, что гра­ ничная орбита обеспечивает переход при минимальной величине вектора скорости в точке М, и поэтому ее ино­ гда называют о р б и т о й м и н и м а л ь н о й э н е р г и и .

Используя зависимости для канонического преобра­ зования кривой второго порядка, найдем, что угол по-

95

ворота f системы O^yji относительно системы

опре­

деляется из выражения

 

 

 

 

kS — (k — соэФ) .

 

(3.89)

 

tg т —

к—г—. —

 

 

 

0 1

sm ф

 

 

Очевидно,

что угол у = ѲГр, т. е. углу

между

плоско­

стью местного

горизонта

и направлением

вектора ско­

рости для граничной эллиптической орбиты. При ф ф %п

 

 

Рис. 3.12. Годограф скорости и пре­

 

 

 

 

 

образование

координат:

 

 

 

 

 

 

ab — параболические

орбиты;

ас — эллип­

 

 

 

 

тические

орбиты

I

рода;

сЬ — орбиты

 

 

 

 

I I

рода; с — граничная

орбита

 

 

 

числитель в выражении

(3.89) всегда

больше

нуля и

знак

tg у будет определяться

знаком

sin ф. Отсюда

сле­

дует,

что при углах

перехода,

меньших

(2л+1)тс,

век­

тор скорости в точке M для

граничной

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

орбиты

направлен

выше,

а

при

ф > (2n + 1 )тс — ниже

плоскости местного

горизонта.

 

 

 

 

ф —іг«

Из

правила

Лопиталя

следует,

что при

tgY = 0

и граничная

орбита

для этого

случая представ­

ляет собой орбиту

Хомана.

 

 

 

 

 

 

 

Все геометрические соотношения и аналитические зависимости настоящего параграфа получены без учета

времени

перехода между заданными точками M и N.

Из этих

зависимостей следует, что при нефиксирован­

ном времени существует бесчисленное множество эллип­ тических орбит, обеспечивающих переход из точки M в точку N.

96

Дальнейшая задача состоит в том, ..чтобы из всего семейства возможных эллиптических орбит выбрать та­ кую (если этот выбор вообще возможен), которая обес­ печивала бы переход из M в N в течение заданного времени т. Для определения такой орбиты воспользуем­ ся уравнением Эйлера—Ламберта [32].

 

 

 

 

§ 3.7. Р А С Ч Е Т Т Р А Е К Т О Р И И

С Б Л И Ж Е Н И Я С

 

 

 

 

П О М О Щ Ь Ю У Р А В Н Е Н И Я Э Й Л Е Р А -

 

 

 

 

 

Л А М Б Е Р Т А

 

 

 

 

 

 

 

 

 

 

Уравнение Эйлера—Ламберта дает связь между

временем

движения от точки

M до точки

N,

величина­

ми r\, г2, ф и большой

полуосью

 

орбиты. В

зависимости

от типа орбиты

перелета

оно

записывается

в следую­

щем

виде [38]:

 

 

 

 

 

 

 

 

 

 

 

 

 

— для параболических

орбит

 

 

 

 

 

 

 

 

 

 

ія

• = — l

IVir.+r,

+ S)* +

 

 

 

 

 

 

 

 

6 У я 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ V{rx

+

r2

-

S)»] ;

 

 

 

 

(3.90)

— для эллиптических

орбит

первого

рода

 

 

/ э і =

j /

l

l [2ш +

-

sin е) + (8 -

sinS)], (3.91)

где

sm -g- =

(/

— 4 |

 

;

sin — =

\

 

 

 

 

 

 

 

- ^ - 1 > - | - ^ 0 J ;

n — число

полных

оборо­

тов,

совершенных

при переходе

 

из точки

M в точку N

(при 0<ф<2тсга=0);

 

 

 

 

 

 

 

 

 

 

 

— для

граничных

эллиптических

орбит

 

 

 

 

' г Р = | / Г - £ І [ ( 2 Л

+ 1 ) т

е +

( 8 г р - з і п 8 г р

) ] >

(3.92)

 

8П> ' 1 / ri + r2

— S / i t

 

8гр .

\

 

 

 

где

an-g-

-

У

r ;

+ r ; +

s

 

\ ~ Т > - 2 > 0 ) '

 

 

 

 

— для

эллиптических

орбит

второго

рода

 

*вИ

= У^-

[2(л +

- ( е -

sin в) +

(й - sin 8)].

(3.93)

4

Сближение в космосе

97

В формулах (3.90) — (3.93) величина S определяется

спомощью равенства (3.72), знак «—» соответствует

случаю, когда ф2тт*Стс, а знак « + » — случаю, когда ф2 тс«>тс.

Уравнение Эйлера—Ламберта является трансцен­ дентным, и для его решения относительно искомой ве­ личины а могут быть использованы различные числен­ ные методы. Однако при

Рис. 3.13.

полета

/ 2

3

C0D,

C,D2Cz

0,5

;,о % г р

Зависимость времени

і от величины х = 1 / а

использовании

 

числен­

ных

методов

необходимо

знать,

имеет

ли

 

данное

уравнение

 

решение,

а

если имеет, то одно или

несколько.

 

 

 

 

Исследование

 

числа

возможных решений урав­

нения

Эйлера — Ламбер­

та

для

эллиптических

орбит

подробно

проведе­

но в работе [38], для

чего

используется

 

разложе­

ние

в

ряд Тейлора

вы-

ражений

(3.91)

и

(3.93)

по

переменной

 

х=1/а,

при

возрастании

которой

происходит

переход

от

параболической

 

орбиты

(а = о о ,

х = 0)

к

эллип­

тическим

(х>0) .

 

 

 

 

Зависимость

при ф<2тс имеет вид кривой /, изо­

браженной на рис. 3.13, где по оси

ординат

отложена

пропорциональная

времени

безразмерная

величина

х =

j / ^ T T ' : ^

П Р И ^=

0,7 и ф = 45°. На этой

кривой

точ­

ка

А0 соответствует

 

параболической

орбите,

интервал

А0В0— эллиптическим

орбитам первого

рода,

точка

Во — граничной

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

орбите

и

интервал

В0Со — эллиптическим орбитам второго рода.

На рис. 3.14 для ф = 30° и различных k изображены

гиперболы годографа скорости в координатах £, ц, а также изохроны для различных значений безразмерно-

98

го времени

т. При

этом

 

дугам

гипербол,

для

которых

t<trp,

соответствуют

эллиптические орбиты

первого

рода, а

при

т > т Г р — орбиты

второго

рода.

 

 

 

1

 

 

 

 

ф = 3 0°

 

 

 

 

0,8

Ч=о,б

 

 

 

 

 

 

 

 

 

 

 

 

 

.

0,t7

/

 

 

 

 

 

 

 

 

 

 

 

0,6

 

\ V л

 

 

 

 

 

 

 

т=о,зѴД

 

 

 

 

Орбиты

 

 

 

 

. 4 °

 

 

 

 

 

 

 

 

первого

 

 

 

Орбиты

\ V

 

 

 

рода

 

 

 

0,4

ѳтораго

\

\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

рода

Yin

 

 

 

 

 

 

 

02

T=0,2Y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т*0,І5

 

 

 

 

 

 

0,2

0,4

0,6

0,8

 

 

 

Рис. 3.14. Годограф скорости в точке M

с и з о *

 

 

 

 

 

 

хронами

 

 

 

 

Из рассмотрения кривой / на рис. 3.13 и кривых

рис. 3.14 следует, что при ф<2тг каждому значению t соответствует одно вполне определенное значение х или одно значение скорости в точке М. Так как величины х или Ѵм однозначно определяют величину а, то в этом случае уравнение Эйлера—Ламберта всегда имеет ре­ шение и 'притом единственное.

В том случае, когда до момента встречи с целью КА

должен совершить более

одного

оборота (n^s 1),

воз­

можны только эллиптические орбиты. Зависимость

t(x)

для ф > 2 т т =1, 2,

...)

имеет вид кривых 2 (п~\)

и

3 (п = 2)

на рис. 3.13.

На

этих

кривых интервалы

DnBn

( t t = l ,

2)

соответствуют

орбитам первого

рода,

точ­

ки Вп

— граничным орбитам,

а

интервалы

ВПС„ — ор-

А*

99

Соседние файлы в папке книги из ГПНТБ