Файл: Брандин, В. Н. Основы экспериментальной космической баллистики.pdf

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

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

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

Добавлен: 17.10.2024

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

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

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

r = V x l + x lj - x l, •и = К'и? + 'У2 + 'Уз;

}j (/=1, 2, 3) — составляющие вектора ускорения

/ = - ^ 0- ^

+ ^ * + — ( Р + Я ) .

гА

т

2.1.3. Модель движения в гринвичской экваториальной системе координат

Гринвичской экваториальной системой координат назовем правую прямоугольную систему координат Ох\*х2*х3*, начало ко­ торой помещено в центр Земли, ось Ох3* совпадает с осью враще­ ния Земли, ось Ох\* направлена в точку пересечения меридиана Гринвича с плоскостью экватора. Эта система координат враща­ ется с той же угловой скоростью 0 3, что и Земля. При записи модели движения центра масс объекта в гринвичской экватори­ альной системе координат в уравнениях (2.1.5) и (2.1.4) следует принять:

r 0= 0, w — 0, w = Q3=cons1, е= 0.

С учетом сказанного и принятых в п. 2.1.2 допущений о дей­ ствующих на объект силах модель (2.1.6) примет вид

V* -= - Н о

+ g * * + - ( P ± R ) - 2 Q 3 X

» * ;

( r f

т

 

 

r*= v \ .

(2.1.9)

где

gr** = gr* —Ц3 X (2 3 X г*).

Соответствующую модели (2.1.9) скалярную модель можно записать в виде

• *

А

* *

*

 

» 1

Х\

» 1

 

• *

= А У

*

*

( 2. 1. 10)

V2

» 2

*

А

• *

*

 

 

* 3

» 3

 

 

 

 

где Xj*, vj* (/= 1, 2, 3 ) — составляющие

радиуса-вектора г* и

вектора скорости»* объекта по осям гринвичской

экваториаль­

ной системы:

 

 

 

 

„•=Г(„;>ЧМ)ЧЩ)’. (Д)Ч(Д)Ч(^)!,

21


f f (/=1, 2, 3) — проекции ускорения

/ * = - ft,

+ Г * + — (Р + * ) - 2Йз X ®*

на оси той же системы.

2.1.4. Модель кеплерова движения

Модель кеплерова движения следует из модели (2.1.7) при допущении, что на космический объект действует лишь одна ос­ новная сила притяжения Земли. Тогда

v = — Ро—*— 5 r = v.

(2.1.11) ,

г 3

Результат проектирования левых и правых частей векторных уравнений модели (2.1.11) на оси основной экваториальной сис­ темы Oxj (/ = 1, 2, 3) представляется системой скалярных урав­ нений

Хл

Н ) — г

Г 3

.. х 2

f t) ~ ~ г

г 3

v a

-*3

Н*о

Г

 

Г3

( 2. 1. 12)

*1 Vi

х2 —

х3

г — \ ^ х \ - J - j cI ; " v = Y

Модель кеплерова движения может быть представлена в ко­ нечном виде. Непосредственное итегрирование уравнений (2.1.12) затруднительно. Интегралы исходных векторных уравне­ ний (2.1.11) проще получить, применяя теоремы о кинетическом моменте и кинетической энергии. Применение теоремы о кинети­ ческом моменте дает векторный интеграл уравнений движения вида

=

(2.1.13)

где с — постоянный вектор, показывающий, что

космический

объект движется в плоскости, проходящей через

центр Земли,

22

положение которой в пространстве не меняется. Абсолютная ве­ личина вектора с определяется по формуле

 

c = rv cos 6,

 

(2.1.14)

где

0 — угол между вектором v

и трансверсалью.

 

 

Применение теоремы о кинетической энергии дает скалярный

интеграл

 

 

 

 

2ft)

------- — ) ,

(2.1.15)

 

\

г

r0 J

 

где

Го, £>о — значения г я v в начальной точке.

рассматри­

 

В качестве интеграла уравнений

(2.1.11) можно

вать также вектор Лапласа

 

 

 

 

A = v X с — !\) — •

(2.1.16)

 

 

 

Г

 

Объединение интегралов (2.1.13) и (2.1.15) приводит к из­ вестному уравнению траектории космического объекта в плос­ кости, определяемой вектором с,

г =

Р

(2.1.17)

1

+ е cos ft

 

где обозначения уже оговорены при записи формулы (1.3.2). Для полного описания кеплерова движения в пространстве элементы р и е дополняют другими элементами. К ним можно отнести, на­ пример, долготу восходящего узла П, наклонение орбиты i, ар­ гумент перицентра со*, время прохождения через перицентр т*, Система кеплеровых элементов орбиты, включающая П, i, со*, р. е и тя, не является единственно возможной. Существует большое число других полных систем элементов. Часто вместо фокального параметра р используют большую полуось

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

3

7\,ск = —^ - 2я,

(2.1.18)

Wo

а вместо времени т* прохождения через перигей часто использу­ ют значение средней аномалии в заданный момент времени

^ —

—*«)= £о —e s in £ 0,

(2.1.19)

где to — некоторый заданный момент времени;

23


Ео — соответствующее этому времени значение эксцентриче­ ской аномалии, связь которой с истинной аномалией вы­ ражается формулой

 

lg^

= l / 4i f t g -

r ;

(2 Л '20)

(оор —.

У м

угловая скорость

движения космическо-

w----- средняя

аIT

го объекта по орбите.

Указанные полные системы кеплеровых элементов орбиты не всегда оказываются удобными. Некоторые из элементов этих систем (оь при е->0, П при г->0) теряют физический смысл, что указывает на возможность описания движения объекта в упомя­ нутых случаях меньшим числом параметров. В связи с этим Лагранж для почти круговых орбит предложил использовать систему элементов

Q, i, k, a, h, М х,

где &= ecos<i)x; h = e sin иц; М х = тк~\-М0,

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

Q, cos г,

k,

a, h, М х,

qu

рх,

hx,

k x, 'Мъ а ,

где ^ ^ s i n i cos 2 ;

sin г sin £2;

hx — е sin (ш„-(-2);

k x = e cos((b„-)-2);

J\42=<лжЧ- 2 -j- M0.

Представляет также интерес для описания кеплерова движе­ ния использовать составляющие с\, с%, с3 и Ль Лг, Лз векторов с, Л по осям основной экваториальной системы координат. Помимо рассмотренных полных систем элементов, в небесной механике используются канонические системы элементов орбиты, напри­ мер, система Якоби, Делоне, первая и вторая системы элементбв Пуанкаре [25, 57, 58].

2.1.5. Модель движения в оскулирующих элементах

Для модели кеплерова движения элементы орбиты остаются постоянными. При действии возмущений реальная орбита отли­ чается от кеплеровой. Для удобства исследований возмущенного движения космического объекта каждой точке реальной орбиты можно поставить в соответствие кеплерову орбиту, которая име­

24


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

ее элементы Q(0 >i(0 >“ 4 0 >/4 0 >£(0 >** (0 , теперь уже зави­ сящие от времени, — оскулирующими элементами орбиты. Для

установления характера изменения оскулирующих элементов при изменении времени t (или какой-либо другой независимой переменной) составляют систему дифференциальных уравнений в оскулирующих элементах. Для указанных выше элементов мо­ дель движения при независимой переменной t имеет вид [57]

Р = 2г / Х;

e — f r smb -\- cos &-]-(e i~ cos ft)— / х;

i = — ■f u cos u\

P

2 = — / sin« cosec i;

<ол =

где

 

 

 

P

 

 

 

 

Cos'S* •> ,

sin ft

(

. , г

\ у

г у .

. .

( 2. 1. 21)

-------/ , +

------

\

1+ * —) A -------/« Sin UClgi;

e

e

p

)

p

 

 

V

P_

 

(eQ sin ft — cos ft) f r

— Q /Xj

r 2

- Ho

L

 

 

 

 

P2

/

, = l

/

~

 

,

 

/ , = \ f - t - u .

7 . = ] /

 

v

HoНо

 

УV

Но

уУ Ho

r = p{ 1 -\-e cos &)_1;

m=

юл -|-H — аргумент

широты;

 

 

 

 

 

COS erfe

 

 

 

 

 

 

(1 + e cos e)3

* — переменная интегрирования.

,

Истинная аномалия связана со временем t уравнением

 

 

 

 

JL

 

 

 

 

t - x — р2

г __ — а* ___ •

 

 

 

*

/и 7 } * + * « » '*

/ г,/т>/и — соответственно проекции возмущающего ускорения на продолжение радиуса-вектора г, на трансверсаль в плоскости

25


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

Удобной формой модели в спекулирующих элементах являет­ ся модель вида [2]

С — X^f3

xbfV

 

сч —xzf \

' x \ f z'i

 

сзxi f ч

x i f \

( 2. 1. 22)

A-l—(^2%

2Сг

/ з С2)>

A2 = (‘D3Ci —^CaJ + f/aCi —/iC3);

A 3— (® 1C2 'V2Cl) ~\~{f\ C1

f 2C\)’>

где fj (/= 1, 2, 3) — проекции возмущающего ускорения на оси основной экваториальной системы координат; Cj, Aj (j = 1, 2, 3) — составляющие векторов с и Л по осям той же системы.

В число уравнений модели (2.1.22) должно быть включено также уравнение для кинематического элемента; в качестве ко­ торого для эллиптических орбит может быть использовано вре­ мя т* прохождения через перигей, для круговых — время т * прохождения заданного аргумента широты и0.

2.1.6. Другие нелинейные модели движения

Помимо рассмотренных выше моделей движения, распростра­ нены модели в прямоугольных объектоцентрических системах координат (начало этих систем в центре масс объекта) и топоцентрических (начало систем на поверхности Земли). Кроме прямоугольных, используются криволинейные координатные, системы — цилиндрическая и сферическая. Модель (2.1.11) кеплерова движения в цилиндрических кординатах запишется в виде

p - p X * = - J ^ - p ;

(Р2^) = 0; -х3= ---- ^ - * 3,

(2.1.23)

/•з

dt

г3

 

 

где р в данном случае — проекция радиуса-вектора г

на

плос­

кость Ох\Х2 опорной системы координат Ох 1X2X3; X — угол,

обра­

зуемый этой проекцией с положительным направлением оси Охь

а координата х3— расстояние объекта от плоскости

0х\х2

(рис. 2.1.3). Та же модель в сферических координатах

 

Г — ГФ2— Д 2COS2ф — ----- ^ 2-

;

 

/-2

 

 

——- (/"2<р)-|-/'2Хsin срcos ср=

0;

(2.1.24)

26