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

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

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

ходящего узла сГ^ и наклонением орбиты /, а положение перигея в плоскости орбиты — угловым расстоянием ш от восходящего узла. Текущее положение КА или цели

на орбите

может быть

определено

истинной

анома­

лией 9 и

радиусом

г.

В некоторых

случаях

вместо

угла 9 используется

угол

« = ш + 9.

 

 

 

 

Рис.

2.1.

Геоцентрическая

и

орбитальная

 

 

 

 

 

 

 

системы координат

 

 

 

 

 

 

2. Axryrzr

геоцентрическая

вращающаяся

система

с началом А в центре Земли

(рис. 2.1). Ось Ауг

системы

направлена по

текущему радиусу-вектору цели,

а

ось

Azv—по

вектору

ее

угловой

скорости

движения

на

орбите.

 

 

 

 

 

 

 

 

 

 

 

 

 

3. Oxyz— орбитальная система

координат

с

нача­

лом О в центре масс

цели (рис. 2.1). Оси этой системы

параллельны

соответствующим

осям

геоцентрической

вращающейся

системы

координат.

Орбита

цели

в об­

щем случае может быть оскулирующей.

 

 

 

 

 

4. OxnyHzH—орбитальная

 

невращающаяся

система.

Начало

этой

системы

совпадает

с началом

системы

Oxyz,

а ее оси перемещаются

поступательно. В началь­

ный

момент (£=0 или 9 = Э0)

система OxHynzB

 

совпадает

с орбитальной

системой.

 

D,

 

 

 

D,

 

 

5.

Сферические

координаты

<р,

Ѳ и

<pH, 0Н

(рис. 2.2). При этом индекс «н» означает,

что углы <рн

и Ѳн

определяются

относительно

невращающейся

си­

стемы

координат.

 

 

 

 

 

 

 

 

 

 

 

20

Заметим, что для описания относительного Движе­ ния космического аппарата и цели геоцентрическую вра­

щающуюся систему, а

также

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

иногда удобнее связывать

с КА.

При дальнейшем изло­

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

Рис. 2.2. Сферические координаты

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

димо

уметь

переходить

от

одной

системы координат

к другой.

 

 

 

 

 

 

 

D,

ф, Ѳ или D, фн ,

Переход от сферических

координат

Ѳн к прямоугольным X, у,

z

или ха,

уа,

zH

осуществляется

по формулам:

 

y = D cos Ѳ sin 9; z=

— D sin Ѳ; 1

X—D

cos Ѳ cos 9;

xH=Dcos8Hcos9H;

yH==Dcos 6Hsin 9H ; zH=

Dsin9H . j ^2 "1 ^

Формулы обратного перехода легко получить при

рассмотрении

рис. 2.2:

 

 

 

 

 

 

__.

 

 

 

D^Vx2+

 

 

y2

+

z*;

 

 

 

9 =

arctg •

=

— arctfif

 

 

 

 

 

 

D = Vxl

 

+

yl

+

z*;

 

(2.2)

 

 

 

 

 

 

 

9H =

arctg-^а-; Ѳн =

-

arctg

 

 

 

 

 

 

 

 

 

 

Ѵ4+Л

21

[Ajj]'-

Связь между двумя системами прямоугольных коор­ динат, имеющих общее начало, может быть осущест­ влена с помощью направляющих косинусов <x,j, под ко­ торыми понимаются косинусы углов между осями Оі 1-й системы (і=Хі, уи z*) и осями Oj /-й системы (/ =

Обычно направляющие косинусы записывают в виде таблицы

 

Х!

Уj

1 г,

Л ,

а и

а, 2

а 13

 

 

 

 

(2.3)

Уі

а 21

а 22

а ? з

 

а 81

а 32

а

з з

или матрицы

 

 

 

 

 

"12

"'13

 

 

*?1

а 22

а 23

(2.4)

 

і.м

а 3 2

а .

j

 

u . 3 3

Тогда координаты точки в системе Oxjt/jZj при пере­ ходе к системе OxiijiZi преобразуются с помощью

матрицы із]:

XJ

Уі = [A,j] • У) (2.5)

*1-

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

Хі

=

a n J C ; + а и у у

+ а 1 3 2 / , j

 

Уі

=

*21*у + а 2 2 У /

+ « 2 3 2 / ,

(2.6)

Если же необходимо перейти от системы 0%ji/jZä к си­ стеме Oxjt/jZj, то преобразование координат осущест­ вляется с помощью транспонированной матрицы

У}

Уі

(2.7)

Uzj J

Lzi

J

22

В координатной форме это преобразование будет иметь такой вид:

Xj

=

* n x t +

Ч1У1

+

S A !

 

у .

=

а 1 2 Х ; +

а 2 2 У «

+

a32z<"> {

( 2 - 8 )

В силу ортогональности преобразования на направ­ ляющие косинусы накладываются следующие условия:

сумма квадратов элементов любой строки или столбца равна единице;

сумма парных произведений соответствующих элементов, взятых в любых двух строках или столбцах, равна нулю;

любой из элементов равен соответствующему ему алгебраическому дополнению.

Запишем теперь матрицы направляющих косинусов между введенными выше прямоугольными системами координат.

1. Матрица 00 ] перехода от системы OxHyazH к системе Oxyz:

'c o s ( » - & o )

И о о н ] =

I

— sin (ft Ѳ0)

 

0

sin (&-&„)

О"

 

cos ( & - » „ )

0

(2.9)

O l .

 

Чтобы перейти от орбитальной системы координат Oxyz к геоцентрической вращающейся системе, доста­ точно воспользоваться формулами параллельного пере­ носа:

 

хг

= х; уг = у + г, zr = z,

 

где г текущий

радиус-вектор цели.

AxryrzT к си­

2. Матрица [Ла г ] перехода от системы

стеме

Axay.Az&:

 

 

 

cos fi, sin и +

И.г] =

 

cos S\j cos и

sin «TL sin /'

+

sin S\_, cos и cos i

— sin fi, sin и cos i

 

sin S\j sin и

 

sin <f^ cos и - f

-cos <f£, sin i

 

cos €{_, cos и cos i

-\- cos (f^, sin и cos i

 

 

cos и sin г

 

sin и sin г

cos г

 

 

 

 

(2.10)

23

 

§ 2.2. У Р А В Н Е Н И Я

О Т Н О С И Т Е Л Ь Н О Г О

 

Д В И Ж Е Н И Я В О Р Б И Т А Л Ь Н О Й С И С Т Е М Е

 

К О О Р Д И Н А Т

 

 

 

При изучении

относительного

движения

КА

и цели

в наиболее

общей

постановке

необходимо

учитывать,

что оба

аппарата

могут совершать управляемый

полет.

 

 

 

 

 

 

Относительное движение цели и КА будем

рассма­

тривать в орбитальной

системе координат Oxyz (рис. 2.1),

начало которой связано с центром масс цели.

 

 

Для записи уравнений относительного движения вос­ пользуемся теоремой Кориолиса, устанавливающей за­ висимость между ускорениями материальной точки в

абсолютном и относительном

движениях:

 

Wt=We+W,+

We,

(2.11)

где Wgj— вектор абсолютного ускорения центра масс КА; We — вектор переносного ускорения; Wr — вектор относительного ускорения КА; Wc — вектор ускорения Кориолиса.

Абсолютное движение рассмотрим в неподвижной системе координат Axay.dza (рис. 2.1). Векторы, входя­ щие в уравнение (2.11), равны:

We =

Г ц +

е X D + ш X (о) X D)\ '

( 2 1 2 )

 

WC

= 2^X

Ѣ,

 

 

где Fe—главный

вектор внешних сил,

действующих

на КА; m — текущее значение

массы КА;

 

Ц—вектор

ускорения начала орбитальной системы координат, т. е. вектор абсолютного ускорения, действующего на цель; ,D — вектор относительной дальности, определяющий положение центра масс КА относительно цели; D,

24

первая и вторая

производные по времени от

вектора

относительной дальности;

£—вектор

углового

ускоре­

ния орбитальной

системы

координат;

швектор угло­

вой скорости вращения орбитальной системы координат. Главный вектор внешних сил Fe в общем случае представляет собой сумму следующих основных элемен­ тов: тяги реактивного двигателя /\_ полной аэродинами­ ческой силы R и силы тяготения G. В дальнейшем при изучении относительного движения будем полагать, что движение происходит в центральном ньютоновском поле тяготения и на таких высотах, где действием аэродина­

мических сил можно пренебречь. Поэтому

 

 

m

 

 

 

 

или

 

 

 

 

 

 

Wa

= P + g,

 

 

(2.13)

где

р— отношение Р/т;

g—вектор

ускорения

от сил

тяготения, действующих на КА.

 

 

 

 

 

Аналогично может быть записано

выражение для №ц :

 

^

= р л + і ю

 

 

(2.14)

где

рц — отнѳшение Рц/тц*;

Рц

тяга

реактивного

двигателя цели; тц—текущее

значение

массы цели;

gii — вектор от сил тяготения, действующих

на цель.

 

С учетом приведенных выше

соотношений

(2.11) —

(2.14) векторное уравнение относительного движения

центра

масс

КА может быть записано в таком

виде:

 

b^p-Pb+g-g^-W^-W,

(2.15)

1 В

случае

необходимости в вектор р, так ж е как и

в рц,

мо­

гут быть включены ускорения от других возмущающих

сил,

дейст­

вующих

на КА.

 

 

 

25

где

І , „ = ш х (ш х23) ;

(2.15а)

U7 ^ s X D. У

Запишем это векторное уравнение в проекциях на оси орбитальной системы координат:

y=Py-p«y+g-g.y-(WJ-(W£)-(Wc)y.

| (2.16)

z =p-Pm+g~g,~

( W J , - (Ю -

(Wdl

где Px,_py, Pz, Рях, Род, Рчг — проекции управляющих уско­ рений р и рц на оси орбитальной системы координат.

Проектируя векторы £ ц и g, получим следующее:

 

 

gax

= 0;

 

 

 

 

 

 

(2.17)

 

 

gn

•0;

 

 

 

l * 2

+ (г + у)* + г*

 

 

gy =

 

 

•гср {г + у)

(2.18)

 

[*2

+ (г + ><)2 + г2 ,

 

 

 

 

 

 

 

1Го£

 

 

 

1*" + (Г + У)" + « " ] ' ' " '

>

где

тсо—произведение

 

гравитационной

постоянной на

массу

Земли, равное 0,39896« 101 5 м3 2 .

 

26

Используя выражения (2.15а), а также выражение для Wc в равенстве (2.12), найдем:

 

 

 

 

 

 

 

(2.19)

 

 

 

 

 

 

<u2z;

 

 

 

 

 

 

 

 

(2.20)

 

 

 

( ^ b

=

2((o,z-a,,y);

 

 

 

 

( ^ ) у = 2 ( ш ^ - ш > ) ;

 

(2.21)

где

я,

г/, z — текущие

координаты

центра

масс КА

в орбитальной

системе

 

координат; шх, ш„, шг проекции

вектора

угловой скорости со на оси

орбитальной систе­

мы

координат;

гх, еѵ,

ez — проекции

вектора

углового

ускорения s на

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

Чтобы система дифференциальных уравнений (2.16) была замкнутой, необходимо иметь_ соотношения, опре­ деляющие проекции векторов рц, р, ш, е на оси _орбитальной системы координат. Проекции векторов рц и р определяются принятыми методами управления объек­ тами. Для каждого конкретного метода или программы управления могут быть получены свои соотношения, определяющие проекции этих векторов на оси орбиталь­ ной системы координат. Более подробно этот вопрос будет рассмотрен в последующих главах. Для определе­ ния проекций векторов ш и s на оси орбитальной си­ стемы координат необходимо получить выражения, уста­ навливающие зависимость характера вращения орби­ тальной системы координат от управляющего ускоре­ ния р ц .

27

§ 2.3. О П Р Е Д Е Л Е Н И Е С О С Т А В Л Я Ю Щ И Х У Г Л О В О Й С К О Р О С Т И И У Г Л О В О Г О

У С К О Р Е Н И Я О Р Б И Т А Л Ь Н О Й С И С Т Е М Ы К О О Р Д И Н А Т

При записи соотношений, устанавливающих зависи­ мость характера вращения орбитальной системы коор­ динат от управляющего ускорения рц, необходимо учи­ тывать, в какой системе координат осуществляется уп­ равление движением цели. В том случае, когда на борту имеется аппаратура построения орбитальной системы координат и управление движением центра масс аппа­ рата производится в этой системе координат, для опре­ деления составляющих векторов со и s целесообразно ис­ пользовать метод оскулирующих элементов. Если же управление целью осуществляется в абсолютной си­ стеме координат, тогда для определения угловой скоро­ сти и углового ускорения орбитальной системы коорди­ нат необходимо использовать соотношения абсолютной системы координат. Рассмотрим последовательно оба указанных метода.

Запишем систему дифференциальных уравнений дви­ жения цели в оскулирующих элементах [38]:

 

 

dt

Ра-Х

У

ъ0 zr>

É1

У%

{р.у sinЪ

 

+ - L . ) cos 8 +

dt

 

 

 

 

 

d(ù

 

— Pixy

cosd

~

Pax

 

V - k t

e

(2.22)

dt

28

где г — отношение р/(1+<? cos ft) ; " — сумма

со + Э; р, е,

Ш , <Л, і, 3 — текущие значения фокального

параметра,

эксцентриситета, аргумента перигея, долготы восходя­ щего узла, наклонения плоскости орбиты и истинной аномалии соответственно.

Решая систему уравнений (2.22) при заданном за­ коне изменения вектора управляющего ускорения рц, определим движение цели, а следовательно, и характер

вращения

орбитальной

системы

координат,

связанной

с ее центром масс.

 

 

 

 

 

Угловая_скорость

вращения

орбитальной

системы

координат ео_будет определяться

значениями

угловых

 

dS\,

dl

du

 

 

 

скоростей

' If

' ~~dt т. е.

 

 

 

 

 

dSl

+

di

da

(2.23)

 

 

 

dt

 

 

Переход к проекциям вектора ш на оси орбитальной системы координат может быть осуществлен с помощью матрицы [Ла г ] (2.10), элементы которой записываются при условии, что £> =0. Тогда

 

 

di .

 

 

d£\,

sin i cos щ

 

 

 

w*

— ~dt S

l n u

 

dt

 

 

 

 

di

cos и +

 

sin / sin«;

\

(2.24)

 

wy

=

dt

 

 

 

 

 

 

 

 

 

 

dfi

 

. .

 

du

 

)

 

 

 

 

 

 

 

 

 

 

 

 

Подставим

в уравнение

(2.24)

значения

-^-,

dfl

du

du> ,

 

 

 

 

 

 

 

 

dt

 

,

из уравнения

(2.22):

 

St > Р а в н о е St

+

 

 

 

 

 

 

: =

0;

 

 

)

 

 

 

 

 

 

 

 

V~*oP

 

 

 

(2.25)

 

 

 

 

 

 

 

 

)

 

 

 

Следовательно, вращение орбитальной системы ко­

ординат

при управляемом

полете

цели

происходит

только

вокруг

осей

у

и z, а при р ц г = 0 только

вокруг

29

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