![](/user_photo/_userpic.png)
книги из ГПНТБ / Каргу Л.И. Системы угловой стабилизации космических аппаратов
.pdfВ соответствии с рис. 4.11 и принятой последователь ностью операций запишем проекции угловых скоростей КА на связанные с ним оси координат:
у + Ф sin &-j-«>0 sin ф cos &;
q= ф cos & cos Y-f-(& — % cost})) sin y—% sin ф sin & cos y;
r——<j>cos&sin y+(^—">0 cos ф) cos Y-(-a)0 sin фэш Ь siny.J
Соответствующие значения проекций угловых скоро стей гироскопов приведены в табл. 4.1.
В рассматриваемом примере за первый угол разворо та КА принят угол г|5. Считается также, что аппарат ори ентирован на Землю, поэтому в приведенных выражениях
фигурирует величина со0. |
(4.8) |
и табл. 4.1 уравнение для |
||||||||||||||||
Используя выражения |
||||||||||||||||||
кинетической энергии запишем следующим образом: |
|
|||||||||||||||||
2 Т |
Л |
, |
V + 1 г *'2 + 'Ы(Р |
cos р/ + |
q sin р/) 2 |
+ |
|
|||||||||||
+ |
/г (q cos |
р / -p |
sin P/)2 |
+ |
/ r |
(r + |
p i ) 2 + /„ (W |
+ |
|
|||||||||
+ |
/»cos JV - q sin px ')2 + |
Im{p |
|
cos fV + |
<7 sin Р Л 2 + |
|
||||||||||||
+ |
/r(<7cos p / ' - ^ s i n |
p1 ")2+/o (c p i "+j t 7 |
cos |
|
|
sin |
^ " ) 2 + |
|||||||||||
+ |
I k 0 (q co.s Р/ + |
r sin |
P2')2 + |
/ м (q cos p2" + г sin |
P2")2 |
+ |
||||||||||||
+ |
/ г (r cos |
P2' - ^ |
sin |
p2 ')2 + |
/ г |
(r cos |
p2" - |
q sin P2")2 + ' |
|
|||||||||
+ |
fT(P |
+ iV)2 + lo(ЪГ+Я |
|
cos К |
+ r sin |
P2 "f |
+ |
|
|
|||||||||
+ |
/ „ |
+ |
<7 cos p2 ' + |
г sin |
h'f |
+ |
Ik0(r |
cos p,' |
+ |
|
|
|||||||
+ |
p sin P3')2 + /,„ (r cos P3" + |
p sin |
P3")2 |
|
+ |
|
|
|
||||||||||
+ |
/Г (p cos p,' - |
r sin B3 ')2 + |
/ г (/7 |
cos |
p3 " - r sin |
p3 ")2 |
+ |
|||||||||||
+ |
/г (ft,' + |
+ h(?S |
+ Г C O S |
|
P,' + |
/7 |
Sin |
Рз')2 |
+ |
|
|
|||||||
+ |
/о (?з"+ r cos Рз" + |
p sin |
p/) 2 + /ЛР |
+ Kf |
+ |
|
|
|||||||||||
|
|
|
|
+ |
/r(<7 + |
№ |
|
|
|
|
|
|
|
(4-9) |
||||
где I r = ha + L ; |
їх*, Iy*, |
Iг*—моменты |
|
инерции |
КА отно |
|||||||||||||
сительно соответствующих осей. |
|
|
|
|
|
|
|
|
|
по
Элемент |
Оси |
конст |
коорди |
рукции |
нат |
|
X |
Гу |
|
Кожух |
У |
|
|
|
2 |
с |
|
Б" |
У |
н |
|
со |
|
|
г |
|
X |
С |
|
>> |
У |
о |
|
|
г |
|
X |
а. |
У |
о |
|
н |
|
|
z |
Спарка 1 |
Спарка 2 |
р cos Pi' + <7 sin р,'
? cos Pi'—р sin PJ'
tpi'+P cosp,' + + 9 sin Pi'
<7 cos Pi'—p sin PJ'
r+Pi
p cos Pi'+t? sin p,"
9 cos Pi"—p sin pt "
r + Pi"
<fl" + P COS Pi" +
+ q sin p,"
q cos Pi"—p sin p,"
Г+Рі"
P+P2 '
9 cos р г ' •—r sin p 2 '
r cos p z ' — ? s ; n Pa'
P+Pz'
4t' + 4 cos p 2 ' + + r sin P'2
г cos p 2 ' — q sin p 2 '
>+P*"
9 cos p 2 " + r sin p 2 °
r cos p2 "— ? sin p2 "
P+PV
<Pl" + <7 C ° S Рг" + +r sin p2 "
r cos p2 "—q sin p2 "
Таблица 4.1
Спарка 3
p cos pa'—г sin p 3 '
tf + Pa'
r cos p,'+p sin Рз'
p cos p a ' — r sin Pa'
<7 + PV
i„'+r COS Pa' —
—p sin Pa'
p cos Рз"—r sin p3 °
?+ Рз"
Гcos Рз"+р sin Pa"
p cos pa"—r sin p,"
•7 + Рз"
f,"+r cos p3 " + + P sin Pa"
Выражение (4.9) записано в предположении, что все
ГИРОСКОПЫ ИДеНТИЧНЫ И ИМеЮТ ОДИНаКОВЫе ОСеВы е / о , /feo
и экваториальные А,, /ь э моменты инерции роторов и ко жухов.
Если все гироскопы имеют .кинетические моменты, рав ные Н, а кинематические связи опарок являются абсолют но жесткими, то, выполнив необходимые вычисления, ли-
пеаризованные уравнения движения їв орбитальной си стеме координат можно привести к ,виду
/ Л + 2 / / ' & - ( / / - / , > 0 ф = Ж Л ;
/ ^ + 2 Я 3 1 ' - ( / / - / / ) % у = Му; |
(4_ ] 0 ) |
Irh~ |
h(> = |
Mr2; |
|
|
|
|
ІгЬ + 2Н% = |
Мг] |
|
|
|
|
|
где |
/ , = / / + |
2 / м + |
4/Г ; |
|
|
|
|
/„ = |
/ / + |
27*0 + |
4/,.; |
|
|
|
/ , = / / + 2 / м + 4/Г ; |
|
|
|||
Мх, Му, Mz— |
возмущающие моменты, |
действующие |
||||
|
но осям стабилизации; |
|
|
|||
Мн ( / = 1 , 2, 3) — моменты, |
действующие |
по |
осям пре |
|||
|
цессии. |
|
|
|
|
|
Каждая пара |
уравнений |
системы (4.10) |
описывает |
движение механической системы КА — ГИО по каналам вращения, рыскания и тангажа. Из этих уравнений вид но, что каналы вращения и рыскания взаимосвязаны, а
движение |
по каналу тангажа |
независимо. |
Если |
соо = 0, |
||
то система |
(4.10) |
распадается |
на три пары |
уравнений. |
||
При достаточно больших углах отклонения гироско |
||||||
пов необходимо пользоваться системой уравнений |
|
|||||
1ху + |
2 Я 8 3 |
cos 83 - |
2ЩХ sin 8Х - |
( / / - V ) <М>= Мх; |
|
|
/г%~Иу |
cos 83 + |
Я (& - ш 0 ) sin 83 = МГЗ; |
|
|
||
/ ^ + 2 / Y ? 1 c o s p 1 - 2 / y p a s i n f l , - ( / / - / / ) a ) 0 V = ^ 1 , ; |
ЕГ |
|||||
/ Д — Яф cos 8 Х + Я у sin З х = УИг 1 ; |
|
І |
||||
/ г & + 2 Я 8 2 с о з rfa—2H$3sm % = |
MZ\ |
|
|
/Г Єа — Я ( » - % ) cos ря + //ф sin В2 = Лїг 2 .
Из этих уравнений видно, что при больших углах от клонения гироскопов все шесть уравнений связаны между собой гироскопическими перекрестными связями.
Системы (4.10) и (4.11) записаны без учета сил вяз кого трения, которое неизбежно присутствует в опорах гироскопов. Обычно между корпусом КА и осями прецес сии гироскопов устанавливаются демпфирующие устрой ства, поэтому каждое четное уравнение систем (4.10) и
(4.11) необходимо дополнить членами (д.(3г, где ц — коэф фициент (вязкого трения. Демпфирующая связь по осям стабилизации в условиях; космического полета отсутст вует.
4.4. Анализ движения полупассивной
гироскопической системы
Исследуем движение канала тангажа полупассивной системы угловой стабилизации при условии малых углов отклонения объекта и спаренных двухстепенных гиро скопов. Решению подлежит последняя пара дифферен циальных уравнений системы (4.10), которая при (о0 = 0, Мг2 = 0 и \хф0 может быть представлена следующим об разом:
|
|
|
/r fj + |
jip — # 6 = |
0. |
|
(4.12) |
|||
|
|
|
|
|
||||||
При нулевых начальных условиях и Mz =oonst реше |
||||||||||
ние этой системы будет иметь вид |
|
|
||||||||
|
|
|
|
|
|
|
|
|
-Da>t |
|
|
9 W = |
^ ^ ( + |
|
i ( 1 _ 4 D |
» ) + |
- F ^ = r |
x |
|||
|
|
X sin |
( ш і / l - Z W |
- f - |
t,)]; |
|
со |
|||
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|||
P(0= |
|
2D |
|
|
-Dmt |
|
(ш]Л _ZJ"f + |
tp) |
||
2H |
|
• Vl—D? |
||||||||
|
|
|
|
|
|
|
sin |
|||
где |
|
|
|
|
|
|
|
|
|
|
|
|
м - |
V2H |
|
|
|
, |
( 4 D 2 — l ) ] A l _ D 2 |
||
|
2/ |
, |
|
; t s ^ a r c t g i |
— |
|
||||
|
г ш |
|
Vlzh |
|
|
|
D ( 4 D 2 — 3 ) |
|||
|
|
|
x9 |
= |
arctg 2 D V\ — D 2 |
|
||||
5 |
1981 |
|
|
|
|
|
2 D 2 _ 1 |
|
113 |
|
|
|
|
|
|
|
|
|
В этом решении под (о следует понимать частоту соб ственных (нутационных) .колебаний системы КА — ГИО.
Если демпфирующая связь между корпусом КА и ги роскопом отсутствует (fi = 0), то решение (4ЛЗ) транс формируется в выражения
»(0= Мг1г |
— cos |
|
|
|
|
2 Я 2 |
|
(4. |
14) |
|
М, |
1 |
||
Р ( ' ) = - |
sin ші |
|
||
2Н |
|
|
|
Как видно из решений уравнений (4.13) и (4.14), при действии по оси стабилизации постоянного возмущения
Mz |
угол |
|3 непрерывно |
возрастает. Систематическое |
воз |
|||||||||
|
|
|
|
|
|
растание угла |
•& при \іф0 |
в |
|||||
|
|
|
|
|
|
ц/2Я раз отличается от ско |
|||||||
|
|
|
|
|
|
рости |
возрастания |
угла |
(3. |
||||
|
|
|
|
|
|
Так, |
при |
Я = 10 |
кг-м2 /с |
и |
|||
|
|
|
|
|
|
ц. = 0,1 |
кг-м2 /с |
скорость |
си |
||||
|
|
|
|
|
|
стематического |
ухода по уг |
||||||
|
|
|
|
|
|
лу f) в 200 раз меньше, чем |
|||||||
|
— E D - |
|
|
по углу (3 и при сравнитель |
|||||||||
|
|
|
но слабом |
демпфировании |
|||||||||
Рис. |
4.12. |
Структурная |
схема |
||||||||||
может |
не |
приниматься |
|
во |
|||||||||
полупассивной системы |
со |
спа |
внимание. |
|
|
|
|
|
|||||
|
ренными |
гироскопами |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
Физически |
возрастание |
||||||
ясняется |
|
|
|
|
угла •& при действии Mz |
объ |
|||||||
наличием перекрестной связи между осями пре |
|||||||||||||
цессии |
и |
стабилизации, обусловленной |
вязким |
трением. |
Гироскоп, прецессируя, встречает сопротивление со сторо ны демпфера, а возникающий при этом момент вызывает отклонение аппарата по оси тангажа.
Периодические значения углов О и р являются зату хающими. Декремент затухания ц/21т определяется вели чиной коэффициента вязкости и инерционностью гиро скопов.
Системе уравнений (4.12) соответствует структурная схема полупассивной гироскопической системы, приве денная на рис. 4.12. Входными величинами в данном слу чае являются возмущающий МЛРЛ и управляющий мо менты Mz, а выходными — углы отклонения f} и (3.
Передаточные функции замкнутой системы для угла стабилизации найденные из структурной схемы, мож но записать в виде
|
0 ^ ( 5 ) = |
l i £ L = |
|
|
; |
|
( 4 . 1 5 ) |
||||
|
|
* |
w |
Mz(s) |
|
s(T2s2 + |
2iTs+l) |
|
|
v |
' |
|
ФА'Ш($) = - H |
S ) |
= |
k- |
|
; |
(4. |
16) |
|||
где |
|
|
|
|
|
|
|
|
|
|
|
k = ^r |
' |
k = |
^ |
f% |
— ^r^z |
' £ = |
^ |
і / |
JjL |
|
|
1 |
2 Я 2 |
' |
2 |
2 Я 2 |
' |
2 Я 2 |
2 |
J/2 |
| / |
Л ' |
|
|
|
|
|
|
|
|
|
|
|
||
а передаточные функции для угла прецессии |
|
|
|
||||||||
|
|
Ф*г(8) = Ш- |
= |
!£±h |
|
; |
( 4 . |
17) |
|||
|
<b*w(s) = - H |
s ) |
= |
|
|
. |
(4. |
18) |
|||
Если параметры гироскопа и аппарата выбраны та |
|||||||||||
ким образом, |
что 0 < £ < 1 , то эти передаточные |
функции |
соответствуют передаточной функции колебательного зве на. Очевидно, что при р = 0 система становится консер вативной ( | = 0). На практике стремятся проектировать системы с возможно большим значением коэффициента вязкого трения ц., используя для этого жидкостные и эле ктромагнитные демпферы.
При движении КА по круговой орбите с угловой ско
ростью ©о для случая, когда |
Mz |
= 0, р = 0, имеем |
|
—%[t |
— sin со/ I ; |
|
|
|
to |
• / |
19) |
|
|
(4. |
При ориентации аппарата на Землю необходимо ис ключить систематическое возрастание угла д. Это можно сделать при помощи датчика момента, установленного на оси прецессии гироскопа, прикладывая к гироскопу по стоянный момент такой величины, чтобы выполнялось ра венство
М м д =щН. |
(4.20) |
5* |
115 |
Если |
на |
КА |
действует |
периодическое возмущение |
||
Mz — Mz0s'm |
at, то |
|
|
|
"v |
|
|
|
|
Izo> (д2 _ <o2)(asin ші — ш sin at); |
(4.21) |
||
|
|
|
4- о»2 cos at—a2 |
cos |
||
|
|
' w |
|
|||
|
|
2Ha\ |
|
|
|
|
Из этого решения видно, что гироскопические систе |
||||||
мы работают без насыщения при периодических |
внеш |
|||||
них возмущениях. |
|
|
|
|||
Для каждой конкретной системы с ГИО должны быть |
||||||
заданы |
допустимые углы |
отклонения |
гироскопов. Так, |
|||
если M z |
—const, то в зависимости от величины Мг |
в пер |
вом приближении время насыщения может быть найдено из выражения
2 Я |
(4. 22) |
|
Mz |
||
|
В прецессионной теории индикаторных гироскопов не учитывают нутационные колебания, так как их круговая частота имеет порядок нескольких тысяч 1/с, а амплиту
ды |
не превышают |
10~5 рад. |
В рассматриваемом |
|
слу |
|||
чае, |
когда |
период |
нутационных |
колебаний |
равен |
не |
||
скольким |
секундам, и даже |
десяткам секунд, |
имеет |
|||||
смысл учитывать гармонические |
составляющие.' |
|
|
|||||
Действительно, |
из решения |
|
(4.14) видно, |
что |
с |
уче |
том амплитудного значения р насыщение может насту пить либо раньше, либо позднее, чем рассчитанное по
формуле (4.22). Значение |
поправки |
AtH можно опреде |
||
лить как |
|
|
|
|
|
д*„= |
+ — . |
|
(4.23) |
При |
М г = 1 н - м ; #=10 0 кг-м2 /с; |
р н = я/3; |
со = 0,4 1/с; |
|
^ = 209 |
с, А/н = 2,5 с. |
КА гармонического |
внешнего |
|
В случае действия на |
возмущения необходимо, чтобы амплитуда вынужденных колебаний гироскопа не превосходила бы значения рн - Используя решение (4.17) при условии со»а, получим
(4, 24)
Из этого неравенства следует, что максимальное зна чение возмущающего момента не должно превышать ве личины, определяемой неравенством,
М zO' :2На$н.
При а = 0,01 1/с и остальных значениях ранее задан ных параметров Мг о^'0,2 Н-м.
Заметим, что для полупассивных гироскопических си стем, работающих в режиме демпфирования колебаний, понятие насыщение гироскопов не имеет практического значения.
В режиме сброса кинетического момента работа си стемы аналогична случаю, когда на КА действует мо мент Mz—const. В этом режиме при помощи других ис полнительных органов, например реактивных сопел, соз дается момент Мр.с, направленный против момента Мг
ипревосходящий его по величине.
Втех случаях, когда аппарат необходимо развернуть на определенный угол, достаточно к гироскопу прило
жить |
момент М. Если М д м |
=oonst, то решение системы |
|||||||
|
|
|
/ г д + 2Яр = |
0; |
|
|
(4.25) |
||
|
|
|
/г £ + ^ - Я » = Л Г д м |
||||||
|
|
|
|
||||||
при нулевых начальных условиях имеет вид |
|||||||||
|
мдм |
|
2D |
-Don |
sin |
(uYl-DH+u) |
|||
|
|
|
|
|
|||||
|
|
|
|
|
|
- £ ) 2 |
|
|
|
|
мдм |
|
-Da>t |
sin |
|
|
|
||
|
1 —]Л — 1)2 |
|
|
|
|||||
|
/ Г 0 ) 2 |
|
|
|
|
||||
где |
|
|
|
|
|
|
|
arctg Vi — Pi |
|
|
|
|
, |
2DV\—D2 |
|
tp = |
|||
|
= |
arctg — |
|
|
|||||
|
|
|
|
|
D |
||||
|
|
|
& |
2D—1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
При |
p = 0 это решение преобразуется |
следующим обра |
|||||||
зом: |
|
|
|
мт |
|
|
|
|
|
|
|
|
|
Л, |
|
1 |
sin Ы |
|
|
|
|
|
|
н |
t~- |
|
|
||
|
|
|
|
|
|
|
|
14-27) |
|
|
|
|
|
|
(1 — cos Ы). |
||||
|
|
|
|
V>2 |
|
Из выражений (4.26) видно, что по окончании пере ходного процесса КА приобретает скорость
|
|
н |
к |
|
|
Интересно отметить, что в отличие |
от маховиков при |
||||
Ж д м =const угловая скорость f} = const и становится |
рав |
||||
ной нулю как только М д м |
=0 . |
|
|
||
Из |
первого уравнения |
системы (4.25) видно, что раз |
|||
ворот |
КА происходит под |
|
действием |
гироскопического |
|
момента |
|
|
|
|
|
|
MT = |
2H$(t). |
|
|
|
Определив из решения |
(4.27) $(t) |
и подставив его в |
|||
выражение для Мг, получим |
|
|
|
||
|
^ г = ^ д м j / ^ s i n W . |
(4.29) |
|||
Для начального момента |
разворота, когда sin at~ |
at, |
а вектор Н отклонился от исходного положения на не
значительный угол (cos (3^1), справедливо |
приближен |
ное равенство |
|
M^M^-^t. |
(4.30) |
' г |
|
Обычно Я^>/ г , поэтому, даже при малых t, у гиро скопических систем проявляется свойство усиления мо мента, которое заключается в том, что при помощи ма лых моментов, прикладываемых к гироскопу, можно соз давать весьма большие моменты, управляющие аппара том.
Маховики не обладают этим свойством, поскольку реактивный момент двигателя непосредственно приложен к корпусу аппарата. Действительно, из уравнения (3.3) следует, что при M A = const и Mz = 0,
откуда |
А = ' ^ М . |
< 4 ' 3 1 ) |
щ
В установившемся режиме угловая скорость разворо та этого же КА при помощи ГИО будет равна
ft — ^ Д м |
(4.32) |
|
г |
Н cos р |
|
|
|
Если f>M = 'Or, то получим
мл |
Ht0 cos |
(4.33) |
|
||
|
|
Очевидно, чем меньше время t0 действия момента М- тем это отношение больше. При
где п — 'число периодов нутационный, колебаний, |
будем |
||||||
иметь |
|
|
|
|
|
|
|
|
|
2лп |
cos ро |
|
|
|
|
Если |
/ z = 1 0 5 |
юг-м2 ; |
/ г = 1 |
иг - м 2 ; |
п=10 , |
(J0 = Jt/3, |
|
МД /АГД М |
= 10. |
|
|
|
|
|
|
Приведенный пример показывает, что при достаточно |
|||||||
больших отрезках времени действия |
на КА моментов М д |
||||||
и Мд м |
свойство |
усиления |
момента |
сохраняется. |
|
||
Определим, в каком случае отношение МЯЩдм |
равно |
||||||
единице. Из (4.33) получим |
|
|
|
|
|||
|
|
|
Н cos Ро |
|
|
|
|
При |
заданных |
параметрах |
/=33,4 |
мин (для # = |
= 100 кіг-м2 /с). Следовательно, только по истечении срав нительно большого времени разворота маховики и двух степенные гироскопы становятся равнозначными. Одна ко переходные процессы в СУС протекают значительно быстрее, поэтому всегда Мл/М д м ^ > 1.
Принцип усиления момента аналогичен принципу ры чага: выигрывая в силе (моменте), проигрываем в рас стоянии (времени разворота).