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

книги из ГПНТБ / Брандин, В. Н. Основы экспериментальной космической баллистики

.pdf
Скачиваний:
50
Добавлен:
22.10.2023
Размер:
12.04 Mб
Скачать

Jx i направим в центр Р визируемой планеты,

а ось Jx i в

сторону визируемой звезды. Обозначим через

половину

углового диаметра визируемой планеты и введем единичные

щекторы е{, е{, е , ев ,

направленные

соответственно

по

осям

Jx i и Jx i и линиям визирования на звезду * и край

В

осве­

щенной части планеты.

Введем тэ,кже

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

r(t) и

rp{t), определяющие положение объекта и планеты относитель­ но центра Земли, и вектор pp(t), определяющий - положение центра планеты относительно объекта. Положение объекта

Рис. 3.6.4. К измерению угла между ли-

Рис. 3.6.5. К визированию на

ниями визирования на звезду и кромку

противоположную освещенную

освещенной части планеты

часть планеты

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

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

(7 = 1 ,

2, 3) (рис.

3.6.4).

Обозначим на­

правляющие косинусы единичных векторов е{, е{,

е* в основной

экваториальной системе

через

е2}, е)

(у'=

1,

2, 3). Пред­

ставим измеряемую функцию угла j\t) через тригонометриче­ ские функции соsx(t) и sinx7):

со,sx{t)=eBe*= e x/B Xl

sin l{t) = {eB X О , : в

В

ХчеУХгч.

*

в

*

(3.6.16)

"^x^xt,

 

где eBt, ев2, е*х,,

—направляющие

косинусы векторов ев и

е* в измерительной системе.

 

 

 

 

 

 

Найдем eBt и ев2. Из рис.

3.6.4

следует,

что

 

в

cos

iji ( О

в

.

sin

Х 2

( О

(3.6.17)

eXl =

 

; еХг= ± :

 

2

 

 

 

 

 

 

 

 

90

t

Знак плюс в последней формуле соответствует схеме визи­ рования, приведенной на указанном рисунке, а минус — схеме, представленной на рис. 3(6.5. Необходимо отметить также, что

Х2<0 R P

 

 

 

(3.6.18)

где ЯР—-известный радиус визируемой планеты.

 

Определим теперь

$

*

 

eXl, ех%. Очевидно,

 

eXl=

е*е{=

е*еи+ е\еп-f е\е13;

(3.6.19)

e Jfi =

e * e 2 = ^ 1 ^ 2 1 _ t- e 2 ^22~f_ e 3^23>

(3.6.20)

Для определения косинусов е3* (у = 1, 2, 3) можно восполь­ зоваться астрономическим ежегодником, где приводятся прямое восхождение а* и склонение 6* звезды, выбранной для визиро­ вания (см. рис. 3.3.7). Зная а* и 6*, найдем

е* = cos 8* cos а*; е \= cos Ь* sin а*; е*ъ—sin 5*.

(3.6.21)

Направляющие косинусы ец (/=1, 2, 3) определяются соот­ ношениями

'О :

хр! № ~

x i W

( j = 1, 2, 3),

(3.6.22)

Рр (0

 

 

 

где X p j ( t ) , X j ( t ) (/= 1, 2, 3 ) — координаты планеты ческого объекта в основной экваториальной системе:

РР(*)=‘У

У=1

и косми­

(3.6.23)

При этом координаты планеты считаются известными функция­ ми времени. Наконец, для определения е2] (у’=1, 2, 3) запишем

 

(е{ Х е * ) Хе{

(3.6.24)

 

<?2 =

 

 

| (< г( X е*) X е{\

 

Но

 

 

 

 

( e i x е*) X

е'-е{(е{е*).

(3.6.25)

Кроме того,

 

| sin (еГе*)\ =V1 —cos2(e( е*),

I

(е{ X е*) X e i 1=

где c6s(e{e*)=^eie*=eXl.

91

Поэтому

 

 

 

( е{ Х е * )

X е{ ■ V

\ - ( е ху .

(3.6.26)

С учетом полученного

запишем

формулу для

определения

e2j в виде

 

 

 

e j ел

и = 1, 2, 3) .

(3.6.27)

-2;

¥ = =

Таким образом получены необходимые соотношения для свя­ зи измеряемого угла %(t) с ' параметрами, характеризующими положение объекта в основной экваториальной системе коор­ динат. При выводе этих соотно­ шений принято допущение о сферической модели визируе­ мой планеты. Если сжатие пла­ неты существенно, что имеет, например, место при визирова­ нии Земли из ближнего космо­ са, то в зависимости от точно­ сти измерения угла %(t) иногда необходимо учитывать поправ­ ку Ах(i), обусловленную сжа­ тием (рис. 3.6.6). К рассмот­

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

ния, приведенные в § 3.3. Запишем теперь измеряемые функции углов р(£) и y(t) для

случая, когда определяется направление на объект из базисной точки, расположенной вне объекта. Если с базисной точкой свя­

зана измерительная система координат

J x j

(/=1,

2, 3) с из­

вестным

началом и направлением осей

в опорной системе

Ох$

(/= 1>

2,

3), то для записи измеряемых функций

р(/) и у(/)

можно

воспользоваться

очевидными

соотношениями

(см.

рис.

3.2.1)

 

 

 

 

 

 

 

 

 

 

 

хъ W

 

 

 

(3.6.28)

 

 

 

(0 < ? < 3 6 0 °);

 

 

 

x{{t )

 

 

 

 

 

 

 

sin у (t)~

x{(t)

(0 < y < 9 0 ° ),

(3.6.29)

 

 

 

 

 

/ j2=i (4 ) 2

 

 

 

 

 

где

 

1, 2, 3)

координаты

объекта

в измерительной

системе.

 

 

 

 

 

 

 

92

Для связи координат xj (t) с координатами Xj(t) в опорной системе можно воспользоваться формулами линейного неодно­ родного преобразования координат

rj(t) = Aj'0(t)r(t) + rJQ(t),

(3.6.30)

где

 

 

 

 

4 [t)

Xi it)

 

М * )=

x{[t)

; r{t)= x 2it)

 

 

х{ it)

*3 it)

 

г 0 it)— вектор-столбец

с компонентами, равными координа­

там начала опорной системы в измерительной;

A JM(t) — матри­

ца перехода от опорной системы координат к измерительной.

§3.7. и з м е р я е м ы е ф у н к ц и и с к о р о с т и

Кизмеряемым функциям скорости yci(t) =Uci{x, t) относятся измеряемые функции радиальной скорости, линейных комбина­ ций радиальной скорости и измеряемые функции угловых скоро­ стей '.

3.7.1.Измеряемые функции радиальной скорости

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

 

Продифференцируем выражение (3.6.1):

 

r(t) = rj(t) + P(t)-

(3.7-1)

Но

 

p(t) = QJ X?(t) + P*(t),

(3.7.2)

где QJ — угловая скорость вращения измерительной

системы

координат в опорной системе, которая считается здесь инерци­

альной; p*(t)

— локальная производная вектора p(t).

 

Обозначим

r[t) через v(t), p*{t) через

v* (t)

и гД /)

через

Vj(t). Тогда вместо (3.7.1) запишем

 

 

 

 

®(0 = ®/W + 2 / Xp( 0 + ®*W

(3-7.3)

или

®(*)= ®пЮ + ®*(*).

 

 

(3.7.4)

 

 

 

где ®п(/) = ®,(*) + Й'/ X р (Д— переносная

скорость.1

 

1 Здесь имеются в виду производные по времени от направляющих коси­ нусов или углов, характеризующих изменение направления линии визиро­ вания.

93

Вектор®*(t)

можно представить в виде

 

 

 

 

 

®* (t)-.

d{ p(Qgp]*

:P(0«p + PW е?>

(3.7.5)

 

 

 

 

aIt

 

 

 

 

 

 

 

где вр—

р (t)

 

 

о

**

 

 

 

единичным вектор;

ер — локальная производная

единичного

вектора.

 

I

 

Умножим левую и правую части выражения (3.7.5)

скалярно

на p ( t ) = p ( t ) e f :

 

 

 

 

 

 

 

 

р(0 ( t ) = P ( O P ( t ) + P 2(t)efeP

 

или

 

 

 

 

 

 

 

 

 

 

 

 

р(П®* ( 0 = Р(0 Р(0.

(3.7.6)

поскольку

 

 

 

 

 

 

йрер = 0.

Учитывая результат (3.7.4), получим

(3.7.7)

Но

Р (0 [t) = Р (0 К (0 + & X р (г1)]= Р (0 Я/ (0.

так как

р (0 [а / Хр(Ю ]=о.

Поэтому

р W р ( 0 = р ( 0 И 0 — ® / ( 0 1 -

Отсюда измеряемая функция радиальной скорости

(t) г} (0) (0 — Vj (Q]

р (0 =

р(0

(3.7.8)

(3.7.9)

При измерении радиальной скорости из начала опорной си­ стемы координат

?{t) = r{ty

r(t)v(t)

(3.7.10)

r(t)

Используя полученный результат и предполагая наличие двух базисных точек /i и / 2, легко записать измеряемые функ­ ции суммы и разности радиальной скорости

Pi+a(*)=PiW ± Рг(*),

(3.7.11)

 

94

где

[г (<)—гл

(0][г>(0 —

(^)]

Pi(*) =

P i V)

(3.7.12)

 

 

[r(t) — rji

( / ) ] [ * > (t)— Vj

( 0 ]

Ра(*) =

P2 (0

(3.7.13)

 

 

Векторной записи (3.7.9) —(3.7.11) соответствует следующая запись в координатной форме:

3

■V

[-*,■ (0 — xjj (0)

[Xj (0 — xjj у)]

2 j

р ( 0 = - ^ -------------------

-

(3.7.14)

 

/ )

Ь

(О — XJj (О]2

 

 

 

 

 

 

 

2 ]

X At) х At)

 

р(^):

7=1

 

 

(3.7.15)

 

 

l / i

ч (0

 

 

 

r

7~i

 

 

2 j

IX- W — ■*/,; (01 [xj (t)— Xj j

(0]

----------------------------------------

 

1 /

2

 

(0)2

 

 

r

7=1

 

 

 

3

(0 — XJJ (01 [X. (t) Xj j (01

 

2

 

7=1

 

 

 

 

(3.7.16)

 

 

 

 

 

j/ i i*,( 0 - х Л/ (OJ2

Г7=1

3.7.2.Измеряемые функции угловых скоростей

Запишем вначале измеряемые

функции

вида

[ c o s 6^. (г)3

(/=1, 2, 3). Из выражения (3.7. 5)

имеем

 

 

?(t)e; =v'{ t) - p (t) et.

 

(3.7.17)

Обозначим е*9 через ve . Тогда (3.7.17)

примет вид

 

Р(0®«р = ®*(*)—Р (*)«(>•

 

(3.7.18)

95

Но из соотношения (3.6.12) получаем

 

 

 

 

=S

^

[C° S W l * ' -

 

(3.7.19)

 

 

 

 

 

 

7 = 1

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

4

[cos 0 ,( 7 ) ] = * ^

( y = l ,

2,

3).

(3.7.20)

at

 

 

 

 

 

 

 

 

Учитывая результат

(3.7.18), запишем вместо (3.7.20)

 

[cos 0,(7)] =

 

(7)

 

p(Q

 

(7 = 1 , 2, 3)

(3.7.21)

P(0

 

p(7)

 

dt

.

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

[cos 0,(7)]:

г «* (7)

p (О

p(*)

"7>

(3.7.22)

 

p

( 0

P2 (t)

 

 

 

 

 

где

 

p(7) =

r ( 7 ) - r y (7);

 

 

 

 

 

(7) =

®.(7) —®n(7);

 

 

 

»nW = fyW + S ; X p(7);

p(0 (*)];

Р(*)= И 0 —Г, (*)]eP-

Это — выражение для измеряемой функции производной на­ правляющих косинусов. Если из базисной точки /, расположен­ ной вне объекта, измеряются углы |3(7) и y(t), то выражения из­

меряемых функций угловых скоростей (3(7) и y(t) можно полу­

чить путем дифференцирования соотношений (3.6.28) и (3.6.29). Поскольку

P(7) =

x{(t)

arctg

 

 

 

х [ (t)

у (7)=

arcsin

•*г(7)

то

 

Р(7)

 

 

 

х{ (t) x j (t)

х{

(t) х{ (t) _

[x[(tj\2 +

 

(3.7.23)

[xi (7)]2

Р (7 )

х { (t) — x{{t) р (t)

 

 

 

(3.7.24)

Р2 ( 7 ) V

[x{(t) ] 2 +

[ т ? з ( 0 ] 2

96

где

 

 

 

Р (0 =

 

 

 

2

х\

 

Р (0 = ' н

р(0

 

 

 

 

 

.*/(/) — составляющие вектора

относительной скорости

объекта

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

 

Для связи x Jj(t) (/= 1, 2, 3) с составляющими ij(t)

(/= 1 ,2 ,

3) скорости в опорной системе координат воспользуемся выраже­ нием (3.6.30). Дифференцируя его, получим

 

О V) =

^у,о it) г ( t ) + A Ji0 (0 г (*)+ r i (*),

(3.7.25)

где

производная матрицы

 

 

 

х {

(t) 1

■*i(0

xio {t)

 

 

r j ( t ) = х {

(t) \ ; r { t ) =

x 2'(t) ; r i { t ) =

x { 0 (t)

 

 

x i

(t) |

x 3 {t)

xio.(t)

 

§3.8. ИЗМЕРЯЕМЫЕ ФУНКЦИИ ОРИЕНТАЦИИ

ИВРАЩЕНИЯ

Кизмеряемым функциям ориентации yoi(t) = u0i(x, t) и вра­ щения yBi(t) = uhi(x, t) относятся измеряемые функции углов и угловых скоростей, характеризующих ориентацию и вращение объекта.

3.8.1.Измеряемые функции углов

Ориентация осей связанной с объектом системы координат Sxy1 (/ = 1, 2, 3) по отношению к опорной системе осей Ох, (/ = 1, 2, 3) определяется тремя углами Эйлера, получающимися как результат трех независимых поворотов триедра осей Sxy1 из на­ чального положения в конечное. Такими углами могут быть уг­ лы прецессии v, нутации б, чистого вращения ф или углы рыска­ ния фь тангажа крена yt [42]. Если на борту объекта установ­ лено два трехстепенных гироскопа, то, как уже было отмечено, можно измерить перечисленные выше или другие углы ориента­ ции объекта. Именно эту информацию использует система угло­ вой стабилизации объекта на участке выведения, орбитальном участке полета и участке спуска с орбиты. Следовательно, сво­ бодный трехстепенной гироскоп во многих случаях позволяет непосредственно измерить параметры ориентации объекта, и для

4—356

97

этого случая измеряемые функции углов запишутся в виде

УоА*)=*№ ( / , у = 1 ,2 ,3 ),

(3.8.1)

где x j(0 — углы ориентации объекта в опорной системе коор­ динат.

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

\ хз

Рис. 3.8.1. Углы, измеряемые

Рис. 3.8.2. К определению

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

направления известного век­

 

тора в опорной системе

 

координат

При применении солнечных датчиков известным является на­ правление на Солнце (см. рис. 3.3.8), при применении магнито­ метров * — вектор Н магнитной напряженности Земли

(рис. 3. 8. 1).

Будем считать, что независимо ot типа используемых,датчи­ ков на борту космического объекта измеряются углы Ai(0 и Д2(0, определяющие направление известного вектора Н в свя­ занной системе осей. Выразим измеряемые углы Ai(() и Аг(0 через параметры движения космического объекта относительно центра масс. Предположим, что движение относительно центра масс описывается невозмущенной моделью Эйлера — Пуансо и Полностью определяется неизвестными постоянными углами ць и о*, характеризующими направление вектора К кинетического

* Применяются магнитометры различной конструкции. На рис. 3.8.1 при­ ведена схема магнитометра, состоящего из двух рамок — внешней и внутрен­ ней. Внутренняя рамка такого магнитометра устанавливается перпендику­ лярно направлению вектора магнитной напряжённости Земли.

98

момента в опорной системе осей (см. рис. 2.2.5), и параметрами

бйо, Vfco, фы>, vfe, фл, характеризующими ориентацию и вращение объекта относительно вектора К (см. п. 2.2.5 § 2.2). Зададим на­ правление вектора Я в опорной инерциальной системе осей Oxj (/=1, 2, 3) известными углами г)Н(0 и aH(t) (рис. 3.8.2). Вве­

дем также единичные векторы ен , е\, el,

направленные соот­

ветственно по вектору Я и связанным осям

и S x i

Кроме

того, введем единичный вектор ен ', направленный по проекции

Н' вектора Я на плоскость

Sx\xl

Имея в виду,

что

О

< Д! (О <

180°;

0 < д2 {t) < 360°,

запишем

 

 

 

 

 

 

 

cos A1(t)=eHeh

cos A2(t)=e\eH';

(3.8.2)

sin д2( 0 = | e\ X ен '\={е\ X

eH') ,.

 

 

 

 

 

 

■*3

Остановимся более подробно

на получении

выражения для

A] (t). Очевидно,

 

 

 

 

 

 

 

 

ене \~ е \

 

е\ъ-\-вг

'

(3.8.3)

где e f, е\] (/ =

1, 2,

3) —направляющие

косинусы ен и е\

в опорной системе координат.

 

 

 

 

Легко видеть (см. рис. 3.8.2), что

 

 

 

в\

= sin

(^sin вн (t);

 

 

 

 

e f = cos Ря (0;

 

(3.8.4)

 

e” =sm4„{t)cos aH(t).

 

Для получения

е \ ]

(/=1,

2, 3) необходимо иметь матрицу

Лo,i перехода от связанной системы координат к опорной систе­ ме Oxj (/=1, 2, 3). Эту матрицу можно представить в виде

^o,i =

(3.8.5)

где Лол — матрица перехода от системы координат, связанной с вектором К, к опорной системе; Ah,i — матрица перехода от свя­

занной системы

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

связанной с

вектором К.

Е сли систему SK1K2K, связанную с

вектором К,

ввести так,

как это показано на рис. 3.8.3,

то элементы a ff

матрицы A0,k

соответствуют табл. 3.8.1.

 

 

 

 

Элементы

af}1 матрицы

Ah,1 можно найти,

совершив три

последовательных поворота на углы Vh, ди, фл с целью перехода

от системы осей SK1K2K к связанной системе 5 а 1л'2Хз ( с м . рис. 2.2.3). В результате получим табл. 3.8.2 направляющих ко­ синусов.

4* 99

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