Файл: Балахонцев Б.Г. Сближение в космосе.pdf

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

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

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

Добавлен: 01.07.2024

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

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

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

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

Дифференцированием уравнения (2.25) по времени найдем значения составляющих углового ускорения ор­

битальной системы

координат:

 

 

*х = Ъ

 

 

, _

гРаг + Риг'

РагГР .

(2.26)

 

 

2р]Лг0 р

 

 

*оР

 

г

2г* ]/\0р

>

 

Воспользовавшись выражением

 

 

г

р

. еЬ sin

е

cos

 

 

 

г

р

1 +

е cos О

 

и уравнениями

(2.22),

исключим

из уравнения

(2.26)

р и г. Тогда

будем

иметь:

 

 

 

*, =

0;

 

 

 

 

 

 

 

 

 

 

 

 

(2.27)

 

 

 

 

е sin 8 +

ра

 

Видно,

что

рЦу

не

оказывает

влияния на

угловое

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

щая углового ускорения ez зависит только

 

от р ц ж ,

в

то

время как гѵ

зависит от рц*, рш

и р Ц г .

 

 

 

 

 

Перейдем к рассмотрению методики определения со­

ставляющих

угловой

скорости

со и углового

ускорения е

при управлении движением центра масс

цели в

абсо­

лютной системе координат.

В

этом

случае

выбранный

закон управления

позволяет

установить

значение

со­

ставляющих

вектора

управляющего

ускорения

р ц

по

осям абсолютной

системы

координат

Axayaza

(рис. 2.1).

Запишем уравнения движения цели в этой системе ко­ ординат:

*"

^о-^а

 

х» — Pax,

 

рау

гз > ;

(2.28)


где г — 1-r х\ + у\ + г*; ха, уа, г а —текущие значения координат центра масс цели в абсолютной системе коор­

динат; P u , , Рцу , Р и г —составляющие

вектора р п по

осям абсолютной системы координат.

_

Угловая скорость орбитальной системы координат ш может быть представлена в виде векторной суммы угло­ вой скорости вращения вектора г и угловой скорости вращения орбитальной системы координат относитель­ но оси у:

= «ог + "V

(2.29)

Запишем выражения для составляющих вектора угловой скорости по осям абсолютной системы коор­ динат:

Ѵа^а

ту

 

'хя^й

(2.30)

 

г'1

 

 

О)

хаУа

У&хъ

 

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

гх

 

 

 

со3

 

(2.31)

 

 

 

ГУ

 

 

- m « _

 

шгга

_

 

где

 

 

 

 

 

 

 

«и а.12

a\3

 

 

 

[Л]

22

ß 2 3

 

 

 

 

32 « 3 3 -

 

On = —

 

va

 

 

Кг

« 1 2 =

ya

« 1 3 -

Vf,

 

v n

 

 

 

31


Уага

?аУз ;. Ö 8 2

Z&xa. ' xaza

 

•*а Va Vaxa ^

ГѴ.

 

 

 

 

rVn

ѵ пх ~ л а

~рГ (xaXa

-fr y a y a

-f~

ZaZa),

УЯпу =

Уа -

- ^ -

+ УаУа

+

Z a Z a ) ;

V \ z

=

Za~-pr

+ УаУа +

ZaZa)\

В результате решения матричного уравнения (2.31) найдем значения составляющих:

= 0;

(2.32)

Второй элемент правой части уравнения (2.29) опре­ деляет угловую скорость вращения плоскости орбиты цели вокруг радиуса-вектора г:

Риг

(2.33)

уа. '

ГДе Рцг = « 3 1 / Ѵ а + ö 32/>«y a + a»»Pmt.

С учетом приведенных выше зависимостей (2.29), (2.32), (2.33) запишем выражения для составляющих угловой скорости ш по осям орбитальной системы коор­ динат:

ш, = 0;

(2.34)

п

32

 

Составляющие

углового

ускорения

s получим диф­

ференцированием

соотношений

(2.34):

 

 

 

 

 

= 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

е у = -^Г

«га (Уа2а— **Уш) + Р^ш*—

 

 

 

* Ä )

+

 

 

 

 

 

 

 

 

 

хаха

+ УаУа

+

 

ZaZ3

+

 

ѵа

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.35)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

/>ц ,а

(*аУа -

Уа**)] + Р ^

(Уа*а — Z a y a

)

+

 

+

Л , У а

(^а^а -

- * Ä ) +

 

ауа

-

у а Л а ) } J

 

 

«, = -f" [ ^ * +

 

 

+ У а ^ а+

2 а ^ а ) ] »

 

 

 

 

 

где.»

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^ л , = *а +

"^рг (*а*а +

УаУа +

 

 

 

 

-

T T [•*. (•*•*. +

УаУа +

« À ) +

 

 

 

+

X, [х* + УІ + il +

 

хаха

+

у а

у а

+

* , z , ) ] ;

 

 

 

 

= Уа +

 

С * * * . + У а у а

+ « , Z , ) -

 

 

 

-

T T [Уа (*а*а +

УаУа +

 

*,«',)

+

 

 

 

+

Уа (*а +

Уа +

^а +

*а*« +

УаУа +

« . « . ) ] î

 

 

 

 

= г* + Цг

(***• + УаУа +

 

 

-

 

 

 

 

— 7 г [ г , ( а д + у а у 8

+ z a z a ) +

 

 

 

+

2 а {х* +

у \ +

Z\ +

ЗД

+

у а

у а

+

 

* , « , ) ] .

 

Для

определения

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

 

угловой

скорости ш

и углового ускорения

s как при управлении

движением

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

33


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

Дополнив систему уравнений (2.16) зависимостями

(2.22), (2.25), (2.26) или (2.27), (2.34), (2.35) при из­

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

§ 2.4. У Р А В Н Е Н И Я О Т Н О С И Т Е Л Ь Н О Г О Д В И Ж Е Н И Я П Р И С Б Л И Ж Е Н И И С Н Е М А Н Е В Р И Р У Ю Щ Е Й Ц Е Л Ь Ю .

Р А З Л О Ж Е Н И Е В Р Я Д Ы Р Е Ш Е Н И Й У Р А В Н Е Н И Й Д В И Ж Е Н И Я

Дифференциальные уравнения (2.16), описывающие относительное движение цели и КА, были получены в предположении, что цель, с центром масс которой свя­ зана орбитальная система координат, совершает управ­ ляемый полет и ее траектория представляет собой кри­ вую двоякой кривизны. Однако на практике очень ча­ сто приходится рассматривать решение задачи сближе­ ния с целью, совершающей свободный орбитальный по­ лет в центральном поле. В этом случае на нее дейст­ вует только_сила притяжения, а сила тяги реактивного двигателя Я ц =0 . С учетом этого условия уравнение (2.15) перепишем в виде

D = р + g - g l l - Wa - Wi - Wc.

(2.36)

Чтобы записать уравнения движения в проекциях на оси орбитальной системы координат, связанной с целью, найдем выражения для проекций векторов Wm, Ws, Wc.

При движении по кеплеровой орбите проекции век­ тора угловой скорости со и углового ускорения s на оси орбитальной системы координат на основании (2.25) и (2.26) будут иметь такой вид:

:0; с«у = 0; шг = 9 _ г ,

(2.37)

= 0; еу = 0; e, = » = - - ^ - e s i n d .

34

В выражении (2.37) величина 0 представляет собой истинную аномалию, соответствующую текущему поло­ жению цели на орбите.

Подставляя

выражения

(2.37)

в равенства (2.19) —

(2.21), найдем

следующие

соотношения:

 

 

 

(2.38)

( Wc)x = -

2 by, ( Wc)y

= 2bx;

( Wc)z = 0.

Используя выражения (2.38), а также выражения (2.17) и (2.18), можем записать уравнения относитель­ ного движения в проекциях на оси орбитальной системы координат:

Х = Рх

+ Рх + by + 2Ьу;

 

 

У = Ру-

•*о(г +у) +

-^т +

Ѵу-І>х-20х;

(2.39)

Z = P z

•к0г

 

 

 

 

 

 

 

В этих уравнениях

 

 

 

 

rx=Vx2+(y

+ ry +

z\

(2.40)

Выражения (2.39) представляют собой систему не­ линейных дифференциальных уравнений с переменными коэффициентами. Методы получения аналитического ре­ шения таких систем в общем случае неизвестны. Однако если предположить, что расстояние между двумя кос­ мическими аппаратами D мало по сравнению с расстоя­ нием от центра притяжения до цели (/)<С),то нелиней­ ные члены можно разложить в степенные ряды, сохра­ нив только линейные и квадратичные члены [3]. Тогда можем записать

1 — З ц у •

•H2 2 - 4j/2 + г 2

(2.41)

где

(2.42)

2*

35


Подставляя выражение (2.41) в уравнения (2.39) и

(2.40) и удерживая члены, линейные относительно у., получим систему:

— Oy 29j> +

а -

82 ) X -

3fl п?ху

= рх

у + Ьх + 2Ьх -

2 +

2 ) у

[І. п (х•2

- 2 y 2 + z 2 ) = / 7 y ;

 

 

(2.43)

 

 

 

z + nïz Spzy = рг.

 

 

 

Если эксцентриситет

орбиты

цели достаточно мал,

то зависимость угловой

скорости

9 и ее

производной 9,

а также радиуса орбиты г от времени можно предста­

вить в виде рядов по степеням

эксцентриситета

е [14].

Тогда с точностью до первой степени е получим

выра­

жения:

 

 

 

 

 

 

 

 

 

è =

ft[l

+ 2ecosft (* — *„)];

 

 

 

»=

— 2п2е

sin ft (t~

t„);

 

 

 

л ' =

я2[1 +

3* cos я (/ — /„)];

(2.44)

 

 

 

 

 

ц =

ц[1 -f- е cos п (/ — /„)],

 

где

n2 =ito/a3 ; |а=1/а;

а — большая

полуось орбиты цели;

ta — время прохождения

перигея.

 

 

 

Используя разложения (2.44), систему (2.43) можно

теперь записать в таком виде:

 

 

 

л: — 2пу =

3(і/г2л*у +

е [п7х cos n{t~

t„) —

 

— 2п7у sin п (t — tp)

+ Any cos n (t — („)] + px;

 

у + 2nx -

Ъп?у = -J3- ця 2 (л;2 -

2y2

+ z2 ) +

(2.45)

+

<? [2п?х sin « (/ — ^n ) + I0n2y

cos я (* — <n) —

 

— 4nx cos я (i! — «"„)] - f

py;

 

 

 

z+

n2z = 3pn2yz

— ЪепЧ cos n {t — *n ) -f- /?r

 

36