
книги из ГПНТБ / Брандин, В. Н. Основы экспериментальной космической баллистики
.pdfJx 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