Файл: Баринов К.Н. Теория полета космических аппаратов.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 17.07.2024

Просмотров: 147

Скачиваний: 3

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

 

 

 

 

 

 

 

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 на новой орбите.

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