ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 01.07.2024
Просмотров: 145
Скачиваний: 0
оси z. Составляющие рпх и рцу не оказывают влияния на угловую скорость орбитальной системы координат.
Дифференцированием уравнения (2.25) по времени найдем значения составляющих углового ускорения ор
битальной системы |
координат: |
|
|
*х = Ъ |
|
|
|
, _ |
гРаг + Риг' |
РагГР . |
(2.26) |
|
|
2р]Лг0 р |
|
|
*оР |
2г |
|
г |
2г* ]/\0р |
> |
|
Воспользовавшись выражением
|
|
г |
р |
. еЬ sin |
— е |
cos |
|
|
|
г |
р |
1 + |
е cos О |
|
|
и уравнениями |
(2.22), |
исключим |
из уравнения |
(2.26) |
|||
р и г. Тогда |
будем |
иметь: |
|
|
|
||
*, = |
0; |
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.27) |
|
|
|
|
е sin 8 + |
ра |
|
|
Видно, |
что |
рЦу |
не |
оказывает |
влияния на |
угловое |
ускорение орбитальной системы координат. Составляю
щая углового ускорения ez зависит только |
|
от р ц ж , |
в |
то |
||||||
время как гѵ |
зависит от рц*, рш |
и р Ц г . |
|
|
|
|
|
|||
Перейдем к рассмотрению методики определения со |
||||||||||
ставляющих |
угловой |
скорости |
со и углового |
ускорения е |
||||||
при управлении движением центра масс |
цели в |
абсо |
||||||||
лютной системе координат. |
В |
этом |
случае |
выбранный |
||||||
закон управления |
позволяет |
установить |
значение |
со |
||||||
ставляющих |
вектора |
управляющего |
ускорения |
р ц |
по |
|||||
осям абсолютной |
системы |
координат |
Axayaza |
(рис. 2.1). |
Запишем уравнения движения цели в этой системе ко ординат:
*" |
^о-^а |
|
х» — Pax, |
7Г |
|
рау |
гз > ; |
(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 |
|
- 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