книги из ГПНТБ / Балахонцев Б.Г. Сближение в космосе
.pdfчивается переход из точки 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Ï |
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 |
