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

книги из ГПНТБ / Баринов К.Н. Теория полета космических аппаратов

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

 

 

120

 

COS

(760 +

= COS {l80+

-&г).

Можно показать, что в энергетическом отношении этот маневр

эквивалентен первому. В самом деле,

 

A V p =

e 2 s l n й[-l/^e,sLn=

]f~e,[(sin(7flb£f)-

Управляющий импульс при этом направлен по радиусу. Это не противоречит, однако, сделанному ранее выводу р том, что вели­ чина А\!г равна удвоенному значению радиальной скорости и на­ правлена в обратную сторону, поскольку на нисходящей ветви орбиты радиальная скорость отрицательна.

Изменение периода обращения. Период обращения Т зависит только от величины большой полуоси а и связан с ней формулой

Т = 2 ас „

'аз'

к

Отсюда заключаем, что для

изменения периода достаточно изме­

нить соответственно величину большой полуоси й .

Установим

зависимость

между АГи Аа . Имеем

откуда

 

г

Аа

= а,

(4.17)

Большая полуось, в свою очередь, связана с перигейным и апогейным расстояниями соотношением

+ г„

а-

2

Поэтому изменение периода обращения сводится в рассмотренным ранее маневрам, связанным с изменением высоты перигея или апогея на некоторую величину А г , равную 2 Д а •

Остается выяснить вопрос о тон, являются ли эти маневры равноценными с точки зрения энергозатрат, для чего исследуем данные маневров в энергетическом отношении. Больная полуось а связана с величиной механической энергии формулой

К

а-

 

121

 

 

 

где 3T= ~2~ + K(J^~ JT)- механическая энергия единицы

массы. Чем

больше механическая

энергия Э }

, тем меньше разность

— - э, ,

тем больше величина

а . При

V — V„

Э — — ,

а

оо.

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

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

сти AV направлен либо по вектору

скорости

,

либо против

него. На основании указанных предпосылок можно записать

(V0 + AV)2

Vo"

 

 

ду2

 

— = V ДV +

 

г

0

 

г

Из формулы (4.18) следует, что приращение энергии АЭ,,обу­ словленное импульсом скорости AV, зависит от скорости V 0 в точке приложения импульса на орбите. Поскольку скорость в пе­ ригее больше скорости в апогее

V

n

7 + е

 

 

 

е

 

то отсюда приходим к заключению, что максимальное изменение

энергии вызывает импульс, приложенный7-

в перигее орбиты.

Таким образом, с точки зре­

 

 

ния минимума энергозатрат управ­

 

 

ляющий импульс при изменении пе­

 

 

риода обращения должен.быть при­

 

 

ложен в перигее орбиты. При этом

 

 

для увеличения периода обращения

 

 

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

4

 

приложен в направлении скорости,

 

 

а для уменьшения периода - про­

 

 

тив вектора

скорости.

Рис.4.5

Количественная зависимость ве­

 

личины AV от

требуемого изменения периода обращения АГ приведе­

на для исходных круговых орбит на рис.4.5.В относительных вели- чинах • — , & эта зависимость универсальна в том смысле,что она

122

применима к любым круговым орбитам. Для эллиптических орбит с тем же периодом обращения Т0 величина импульса скорости ДУ всегда меньше, чем для круговых, поэтому для околокруговых ор­ бит приведенная на рис.4.5 зависимость может использоваться в качестве верхней оценки требуемых энергозатрат.

Поворот плоскости орбиты относительно радиуса орбиты в точке приложения импульса. Поворот плоскости орбиты относится к некомпланарному маневру. Необходимость изменения плоскости

орбиты возникает при межорбитальном переходе, когда

требуется

перейти на орбиту, лежащую в новой плоскости.

 

 

 

 

Рассмотрим простейший некомпланар­

ный маневр, когда

требуется

повернуть

,Vt плоскость орбиты

на некоторый угол

у

без изменения геометрии орбиты. Для вы­

полнения этого маневра управляющий им­

пульс должен быть приложен так, чтобы

скорость V и угол наклона

вектора

ско­

рости к местному

горизонту

б

,

от

кото­

рых зависят параметры р ,

е

,

были бы

одинаковы до и после приложения управля­

 

ющего импульса (рис.4.б). В соответствии

 

с основным свойством импульсного манев­

 

ра радиус-вектор точки приложения импуль­

 

са является общин для исходной и новой

 

орбиты. В данном случае относительно

Рис.4.6

этого радиуса и происходит поворот плос-

кости орбиты на угол 3- .

 

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

По условию выполнения маневра имеем

 

V = v r , + A

v ;

 

И8 равнобедренного треугольника непосредственно находим

искомую зависимость

 

 

A V =

2 V,., s i n ~ •

(4.19)

123

 

 

При малых углах гг , когда допустимо

считать s i n — —

,

будем иметь

 

 

4 V = V r i j - t

i

(4.20)

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

трансверсальной составляющей скорости в точке приложения импуль­ са на орбите. 5 свою очередь эта составляющая скорости зависит от истинной аномалии точки приложения импульса

V l c i = Vj" ( f + e c o s £ , ) .

Функция VT l (i5",) имеет минимум в апогее и максимум в перигее. Минимальные энергозатраты для совершения рассматриваемого ма­ невра будут в том случае, когда в точке приложения импульса трансверсальная составляющая скорости имеет минимум[см.(4.20_ . Поэтому энергетически оптимальный маневр характеризуется пово­ ротом плоскости орбиты в точке апогея. Для орбит с большим экс­ центриситетом разница скорости в точках апогея и перигея весь­ ма существенна. Например, для орбиты типа " Молния":

Уп _

I + е_,

I

+

0.75

_ -

V A "

1-е ~

I

-

0,75

 

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

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

рядка

300 - 500 км на 1° требуется сообщить величину скорости

AV =

135 м/сек. Для поворота плоскости орбиты на угол j - = 90°

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

Бели исходная

орбита круговая, то при этом

 

 

 

AV = v T V K p

= V n a p .

 

 

Определим

теперь направление управляющего импульса относи­

тельно трансверсали. Из треугольника скоростей

(рис.4.б)

имеем:

 

о/, = 90 +

.

(4

.21)

124

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

§ 4 . 3 . ПЕРЕХОД С ОРШТЫ В ЗАДАННУЮ ТОЧКУ

Предположим, что на некоторой плоскости

(рис.4.7)

опреде­

лена исходная орбита своими параметрами

 

р,

, е , . На этой ор­

 

бите

точка I

приложения управ­

 

ляющего импульса

имеет

истинную

 

аномалию

&t ,

а точка 2

опреде­

 

лена

полярными координатами

гг ,

 

• f t j .

Требуется найти орбиту

пе­

 

рехода из точки I в точку 2,ве­

 

личину и направление

управляю­

 

щего

импульса

AV .

 

 

 

 

 

 

 

Обозначим параметры переход­

 

ной орбиты через

р

 

и

е .

За­

 

пишем в принятых

обозначениях

 

уравнения связей, налагаемых

 

на маневр:

 

 

 

 

 

 

 

 

 

 

 

 

-ecostf*Р

 

 

 

 

 

 

 

 

Р,

7+<

 

'

 

 

(4.22)

 

 

 

 

 

 

 

— у

 

 

 

 

 

 

 

7+e,cosiX

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис.4.7

 

 

+ ecos(3+<£>)

 

2?

 

 

 

2

 

 

 

 

 

 

где

 

приложения импульса на пере­

-истинная аномалия точки '

7

 

 

 

 

 

 

 

 

 

ходной

орбите; Ф = д'г- ^ - дальность

 

перехода.

 

 

 

 

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

AV_

V h

\TYP^~

\ (4.?3)

 

AV,

 

 

AV,

COSoJ, = AV,

 

r,

AV,

 

 

 

 

 

 

 

125

 

 

 

 

 

 

 

где

 

of, - угол

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

AV, и трансверсальным направ­

лением.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Функциональную

зависимость

AV, от параметров,

подлежащих

определению,

будем

называть функцией импульса. В данном

случае

требуется определить параметры переходной орбиты р

,

е

и истин­

ную аномалию точки приложения импульса на переходной орбите.

 

Поскольку три неизвестных

величины связаны между

собой дву­

мя условиями

( 4 . 2 2 ) , то представляется возможным свести

функ­

цию импульса (4.23)

к функции одного переменного. Возьмем в ка­

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

р

и преобразуем

выражение (4.23) с учетом

( 4 . 2 2 ) . Исключим

сперва величину

т9"* . Из первого уравнения

(4.22) имеем

 

 

 

 

 

 

 

 

 

Р

1

 

 

 

 

 

 

 

 

 

 

 

 

cos

Л

 

 

 

 

 

 

(4.24)

 

 

 

 

 

 

 

 

 

 

 

 

 

Второе

уравнение (4.22)

представим в виде

 

 

 

 

 

 

 

c o s ^ c o s # ~ s L n

Л

sLn<$ -

 

 

 

 

Решая их совместно

относительно

si-nt?*,

находим

 

 

 

 

 

sin

0 ; = (

-

/

 

 

 

 

7

 

 

 

(4.25)

 

 

е

cos Ф

е

J

sin Ф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Зависимости (4.24) и (4.25)

при заданных

р а е

однозначно

определяют третий параметр

. Возводя в квадрат

левые и

правые части равенств

(4.24)

и (4.25) и складывая

их,получим

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

(4.26)

 

 

 

 

 

 

 

 

 

 

 

 

s i n

Ф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Наконец, для радиальной составляющей скорости на переходной

орбите в точке I

имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

е s i n #* = ••/—

 

 

. .

 

• .

(4.27)

 

 

 

 

 

 

VP

 

 

s m Ф

 

 

 

 

 

Подставив найденную

зависимость

для радиальной

скорости

(4.27) F исходное

выражение

( 4 . 2 3 ) ,

получим

 

 

 

 

Щр)=

J ^

V p -

VpSf+к

(f-7)cos$-^->)

e,S Ln^

 

(4.28)

 

 

 

 

 

 

 

 

Vp

sin <£

 

V/T

 

 

 

126

Функция (4.28) определена на переходной орбите, удовле­

творяющей заданный уравнениям связей

( 4 . 2 2 ) .

Поэтому

она каж­

дому значению фокального параметра р

ставит

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

некоторое

значение AV, , при котором

обеспечивается

переход

из точки I

в заданную точку 2.

 

 

 

Таким образом поставленная задача перехода с орбиты в за­ данную точку не имеет единственного решения. Для получения единственного решения деобходимо задать еще одно какое-нибудь

условие, характеризующее

маневр. В качестве такого

условия

примем

 

 

A V = m l n A V ( p ) .

(4.29)

Найденная из условия

(4.29) переходная орбита

обеспечи­

вает выполнение заданного маневра при наименьших энергозатра­ тах и в этом смысле решение поставленной задачи является опти­ мальным.

Для отыскания минимума функции (4.28) необходимо проварь-

ировать значение

фокального параметра

р

,вычислить последова­

тельность значений AV,, и выбрать из нее

наименьшее значение

AV7 .Тогда

р о п т

при ДУ7 / п : п полностыо

определит собой перехЪд»-

ную орбиту

в соответствии с формулами

( 4

. 2 4 ) , (4.25) и ( 4 . 2 6 ) .

Остается лишь найти область возможных изменений фокального пар; метра р . Для этого обратимся к формуле

я

'пар

I

Если предположить, что запас характеристической

скоро­

сти АУ^о п расходуется полностью, то область определения допу­ стимых значений р найдется следующим образом:

где в левой части неравенства стоит минимально возможное зна­ чение фокального параметра а в правой г максимально возмож­ ное.

Функция импульса

A V

( р )

(4.28) выпуклая. Поэтому в

области p m i n - р

- ртах

о

ы а и а е

е г единственный минимум

(рис.4.8), и для

отыскания его

можно использовать стандарт­

ные процедуры градиентного метода при вычислении на ЭЦВМ.

127

До сих пор предполагалось, что начальная точка перехода

фиксирована (задана истинной аномалией

я?",

) . В случае,

ког­

да требуется перейти с ис­

 

 

 

 

 

 

ходной орбиты в заданную

 

 

 

 

 

 

точку, возникает

необходи­

 

 

 

 

 

 

мость онределеню

наряду

с

 

 

 

 

 

 

управляющим импульсом

и

 

 

 

 

 

 

положения точки.старта

с

 

 

 

 

 

 

исходной орбиты. При этом

 

 

 

 

 

 

функция импульса

имеет

тот

 

 

 

 

 

 

же вид ( 4 . 2 8 ) , но является

 

 

Pmin

Ропт

 

 

уже функцией двух перемен­

 

 

 

 

 

 

Рис.4.8

 

 

ных:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Vp

sin Ф

^

J

 

 

 

 

 

 

 

где

 

 

Л

 

 

 

 

 

 

Л

=

 

 

 

 

 

 

7 + e7

COST?,.

 

 

 

 

7

 

 

 

 

Результатом

решения задачи являются величины р и v}

,при

которых функция импульса достигает наименьшего значения. Есте­ ственно, что минимизация функции по двум параметрам представ­ ляет собой сложную задачу с вычислительной точки зрения. Для

упрощения решения этой

задачи исследуем

частный случай,

ког­

да дальность перехода

Ф

равна

К

,

тем более,

что при

Ф =

= 5t - в выражении (4.28) имеет

место

неопределенность

типа ~ .

При Ф = ЗС условия

(4.22) принимают вид

 

 

 

 

7+есозтЗ^ =

П

7

+ е с о з ^ + З Е )

Г г

'

 

 

откуда получаем '

 

 

 

 

 

 

 

г, +

г 2

 

 

 

 

 

 

(4.30)

 

 

 

 

 

 

 

 

 

Таким образом, при дальности перехода, равной 5с

 

, фокаль­

ный параметр переходной орбиты

р

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

непосредствен­

но из граничных условий

(4 . 22) . Вместе

с тем оказывается опре­

деленной и трансверсальная

составляющая

импульса

 

 

 

 

128

Поэтому минимизация функции импульса сводится к минимизации радиальной составляющей импульса

4 у _

[ e s L n ^* -

е ?

&i

 

Если существует переходная орбита,

проходящая через

точки

I и 2, при отсутствии

радиальной составляющей импульса,

то

это решение и будет отвечать минимуму характеристической ско­ рости.

 

Исследуем радиальную составляющую импульса AVpHa предмет

обращения ее в

нуль:

 

 

 

 

 

 

 

 

 

 

 

 

esinfi*

_

er

sin

# i

_

 

 

 

 

 

 

~~W

 

 

v p ;

 

- °

-

 

 

Это равенство

будет

для переходной

орбиты выполняться,

ес­

ли

существует e s t n i ^ ,

удовлетворяющее

условию

 

 

 

 

 

e s l n i S ^ = е, s i n - д г у

 

 

 

 

Преобразуем

это равенство

с учетом первого условия

(4.22)

и

выражения ( 4 . 3 0 ) :

 

 

 

 

 

 

 

 

 

 

 

Р

 

 

 

_

 

2 г,

 

 

 

 

 

/ + e c o s # * ~ Р ?

^

 

г] + г 2

'

 

 

откуда

 

е cos Л

 

г,-

п

 

 

 

 

 

 

 

 

 

'?

 

 

 

 

 

 

 

 

 

 

 

рг +

г,

 

 

 

 

 

Исключая истинную аномалию

тЭ"*

,

получим

 

 

 

е ~ ( e , s „ n t f j

 

 

+

[-jtt^-J.

(4.31)

Правая часть равенства (4.31)

при любых значениях ej%

,

rj

г2

положительна. Отсюда следует,

что всегда

существует

пере­

ходная орбита, удовлетворяющая граничным условиям (4.22) и

 

условию

( 4 . 3 1 ) .

 

 

 

 

 

 

 

 

 

 

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

дальности перехода ф = ЗС осуществляется

с помощью трансвер-

сального

импульса

 

 

 

 

 

 

 

 

 

129

Варьирование истинной аномалии тУ, точки старта с исходной орбиты эквивалентно варьированию дальности перехода Ф = г9"г-тЭ' Покажем, что для круговой исходной орбиты переход в задан­

ную точку строго оптимален именно при дальности § = 5t . Если

исходная орбита круговая, то начальные уело-!

 

вия не

зависят от

положения точки

старта на

2

орбите

( г = const,

V = c o n s t ) .

Поэтому

зада­

 

ча оптимального перехода с орбиты в

задан­

 

ную точку трансформируется в задачу измене­

 

ния орбиты с минимальным изменением механи­

 

ческой энергии. Минимуму же изменения энер­

 

гии Э,

отвечает минимум изменения

большой

 

полуоси

а . Последнее условие

выполняется

 

лишь в

том случае,

когда заданная

точка

2

 

лежит на самой оси апсид переходной

орбиты,

 

аэто соответствует дальности перехода § =

=ЭГ (рис.4.9). Поэтому, если истинная ано­

малия точки старта по условию задачи не задана, то достаточно ограничиться решением при $ = 5t , и это решение при старте с круговых орбит будет строго оптимальным, а при старте с эллип­ тических, орбит с небольшим эксцентриситетом - приближенно оп­ тимальным.

. § 4 . 4 . МЕЖОРБИТАЛЬНЫЙ ПЕРЕХОД

Предположим, что заданы исходная орбита и новая орбита, на которую требуется перевести КА. Как и в § 4 . 3 , рассмотрим переход в плоскости исходной орбиты.

Пусть исходная орбита

задана

параметрами р}

, е 7 ,

со,, а

новая орбита - параметрами

р 2 ,

е2 , с о 2 . Кроме

того,

опреде­

лены точки старта с исходной орбиты и прибытия на новую орби­ ту аргументами широты uJt иг соответственно. Требуется опре­ делить переходную орбиту и управляющие импульсы скорости.

Для совершения заданного маневра необходимо приложить ми­ нимум два управляющих импульса. Первый импульс прикладывается в точке старта I с исходной орбиты и служит для достижения точки 2 на новой орбите, а второй импульс прикладывается в мо­ мент достижения точки 2 и служит для уравнивания скорости при­ бытия с орбитальной скоростью точки 2 на новой орбите.

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