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

книги из ГПНТБ / Каргу Л.И. Системы угловой стабилизации космических аппаратов

.pdf
Скачиваний:
37
Добавлен:
24.10.2023
Размер:
9.13 Mб
Скачать

В соответствии с рис. 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 -

Х 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 ші

 

 

 

 

Как видно из решений уравнений (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 Н-м.

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

В режиме сброса кинетического момента работа си­ стемы аналогична случаю, когда на КА действует мо­ мент Mzconst. В этом режиме при помощи других ис­ полнительных органов, например реактивных сопел, соз­ дается момент Мр.с, направленный против момента Мг

ипревосходящий его по величине.

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

жить

момент М. Если М д м

=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.

Принцип усиления момента аналогичен принципу ры­ чага: выигрывая в силе (моменте), проигрываем в рас­ стоянии (времени разворота).

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