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

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

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

оси 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

Для интегрирования системы (2.45) принимаем ве­ личины (л и е в качестве малых параметров и предста­ вим ее решение следующим образом:

JC — JCf£

"~|~" ^*^2*

(2.46)

2 = zK + pzx + ez2.

Подставив уравнения (2.46) в систему (2.45) и при­ равнивая члены, содержащие р. и е в одинаковых сте­ пенях, получим системы дифференциальных уравнений для определения интересующих нас величин:

 

 

 

хк-2пукх;

 

 

 

 

 

ук

+ 2пхк

— 3п2у=ру;

(2.47)

 

 

 

zK + n2zv

= рг\

 

*і — 2/zy\ =

Зп2хкук;

 

у\

+ 2пхх

-

Зп?Уі =

4 - " 2 С*2 - 2 Л 2 +

(2.48)

Zj +

« 2 z t

=

3/г2 ^к гк ;

 

 

Jc2 2пу2

= /г (4^/к -f- лл:к ) cos /г (^ — / п ) —

 

— 2п2ук

sin n(t

— tn);

 

 

y2 + 2nx2

— 3n2y2

= n(\0nyK — 4xK) cos n(t— }

(2.49)

— tR) - f 2n2xK sin n(t — t„);

 

z2 + n2z2

= — 3n2zK cos n(t — / п ) .

 

Системы (2.47) — (2.49) должны решаться при сле­ дующих начальных условиях:

хк(0)=х0;

ук(0)=у0;

 

z K ( 0 ) = z 0 ;

x,(0) = Xal

Ук(0)=у0,

zK(0) = z0,

Xt ( 0 ) = ^ ( 0 ) = г < ( 0 ) = 0 ;

(2.50)

 

х.(0)=уі(0)=2і(0)

 

= 0;

2. j

37

Задание условий в виде выражений (2.50) предпола­ гает, что в момент начала сближения ^ = 0, а в момент его окончания t = i. Поэтому в системе (2.49) t n <!0 и по модулю равно времени движения цели от перигея своей орбиты до точки, в которой она находится в момент на­ чала сближения.

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

новки,

создающей ускорения рх, pv,

Pz-

Тогда

измене­

нием

координат

за

время работы

двигательной

уста­

новки

 

можно пренебречь

и, положив рх = рѵ — рг = 0,

счи­

тать,

что изменение

скорости КА происходит мгновенно

(импульсно). При начальных условиях

(2.50)

решение

системы

(2.47)

можно записать в таком

виде:

 

 

 

 

 

^

— Зуо) sin nt -

- ^ - cos nt

+

 

 

 

+

6y,

 

3 -v ° W - K o +

2 ^ ;

 

 

(2.51)

 

 

 

 

 

 

 

 

 

 

yK

=

 

sin nt +

_ 3y0 ) cos nt +

4 y 0 -

 

 

 

zK = -^- sin nt -f- z0 cos nt.

Уравнения (2.51) представляют собой решение ли­ неаризованной системы уравнений относительного дви­

жения для

круговой орбиты

цели.

В этом случае п =

= ш= ]^оІг3

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

собой

угловую скорость

движения цели по орбите. В дальнейшем при

рассмо­

трении задач сближения с целью, движущейся

по кру­

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

только

линеаризованное решение в виде

уравнений

(2.51).

В этом случае индекс «к» мы будем

опускать.

 

Подставляя полученные выражения для хк, ук, zK в уравнения (2.48) и (2.49), найдем поправки, учитываю­ щие влияние квадратичных по относительному расстоя­ нию членов выражения (2.48) и эксцентриситета орбиты цели как величины первого порядка равенств (2.49).

38

/

Выражения для этих решений удобно записать в сле­ дующем виде:

xt

=

А'0

•+ А\

sin nt +

Л'2

cos

nt

- f Л 3

sin 2л/ -f-

 

 

cos 2л/ +

A'5nt + A'6 nt

sin nt

- f

л/ cos nt;

 

yt

=

ßj, -+- ß'j sin л/ +

£ 2

cos

л/ + JS3

sin 2nt

- f

 

- f

ߣ cos 2nt

+

Bl5 nt

+ B'6 nt

sin л/

+

 

 

|

(2.52)

+

B'7nt

cos nt

- f ß< (л/)2 ;

 

 

 

 

 

 

 

2 ;

=

Cl0

- f C'j sin л/ - f C'2

cos

л/ + C3 sin 2л/

+

 

+

C\ cos 2л/ +

Q nt sin л/ -f-

nt

cos л/,

 

 

где индекс i

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

1 и 2.

 

 

При

і—\

 

получим

решение

для

круговой

орбиты

цели,

а

при

і = 2—поправку

 

на

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

Коэф­

фициенты в выражениях для Х\, уи Z\ вычисляются по формулам:

А\

=

- 3 - ^ 0 ^ +

36 А- у 0 -

Ю[%)

-

2

 

 

 

- 3 0 у > - 3 ^ - г 2 - 2

 

 

 

 

Л 2

= 2 - ^ - - 6 у 0 ^ - - 3 . * 0 у 0 - 2 г 0 ^ - ;

 

4

=

 

+

Т

Г7Г +

^Уо-Т

 

Г Уо +

(2.53)

 

 

 

 

 

+ - г Н5-

-

'0»

 

 

 

 

А\

= *оУо

_1 и JL л. J _ , iL -

 

 

 

А\

=

3[xl +

11

 

 

 

 

 

2

 

 

 

 

 

 

 

 

2

- 2 1 у 0

^ +

2 * о ^ -

 

 

 

 

 

 

 

39

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