книги из ГПНТБ / Балахонцев Б.Г. Сближение в космосе
.pdfоси 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
Для интегрирования системы (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