Устойчивость регулярной прецессии гиростата-магнетика в магнитном поле

Бесплатный доступ

Рассматриваются условия существования и устойчивости регулярной прецессии гиростата, движущегося относительно центра масс и обладающего магнитными свойствами, в стационарном поле магнитного диполя.

Гиростат, магнетик, магнитное поле, устойчивость, регулярная прецессия

Короткий адрес: https://sciup.org/14729944

IDR: 14729944

Текст научной статьи Устойчивость регулярной прецессии гиростата-магнетика в магнитном поле

Под магнетиками понимаются любые материальные объекты, обладающие магнитными свойствами. Иначе, магнетики - это материальные объекты, наличие которых способно видоизменить или возбудить магнитное поле (МП) [1]. В определённом смысле магнетики являются аналогами диэлектриков.

Твёрдые тела-магнетики, движущиеся в околоземном космическом пространстве, динамически взаимодействуют с магнитным полем Земли ( геомагнитным полем, ГМП ). Основная часть этого поля может быть аппроксимирована магнитным диполем , ось которого образует с осью вращения Земли некоторый фиксированный угол [2, 3]. На основе данной дипольной модели построена теория движения орбитального космического аппарата-магнетика, моделируемого гиростатом, в ГМП [3].

Одной из задач этой теории является видоизменённая ограниченная задача о движении относительно центра масс гиростата-магнетика в стационарном ГМП, вектор напряжённости которого H сохраняет неизменное направление в инерциальном пространстве. В данной задаче ГМП моделируется прямым магнитным диполем, ось которого совмещена с осью вращения Земли [3].

Предполагается, что с носителем гиростата-магнетика неизменно связан постоянный магнит с магнитным моментом I 0 . Этот магнит может быть применён в качестве демпфера собственных колебаний гиростата. При этом демпфирующий момент, действующий со стороны демпфера-магнита на гиростат, равен

M = I 0 x H . (1)

Соотношение (1) имеет место при условии G1 /G <<  1, где G , G i - модули кинетических моментов гиростата-магнетика и магнитного демпфера относительно их центров масс. Это равенство может нарушаться лишь на очень малом интервале времени при импульсивном режиме изменения состояния механической системы [4].

1.    Основные предпосылки

Рассматривается движение относительно центра масс С в дипольном магнитном поле (МП) свободного от связей гиростата с заданным постоянным результирующим гиро-статическим моментом. Воздействие всех внешних сил, кроме сил МП, не учитывается.

Введём правые координатные ортобазисы с общим началом в полюсе С: неподвижный базис Z (Cz1z2z3), неизменно связанный с инерциальным пространством, и подвижный базис X (Cx1x2x3), оси которого направлены по главным в полюсе С направлениям тензора инерции гиростата.

Пусть s ( s 1 , s 2 , s 3 ) – направляющий орт, для которого

( s 1 , s 2 , s 3 ) = (sin 0 sin ф , sin 0 cos ф , cos 0 ). (2)

Здесь и всюду далее θ , ϕ , ψ – углы Эйлера, определяющие ориентацию базиса X относительно базиса Z.

Обозначим: A = diag ( A 1 , A 2 , A 3 ) - матрица тензора инерции гиростата в полюсе С ; ω ( ω 1 , ω 2 , ω 3 ) – абсолютная угловая скорость носителя гиростата; k ( k 1, k 2 , k 3) – постоянный гиростатический момент относительно полюса С , заданный в базисе X .

Постоянный магнит, жёстко закреплённый на носителе гиростата ( магнитный демпфер ) имеет магнитный момент I 0 (0, 0, I 0 υ 0 ), где I 0 = | I 0 |, υ 0 = const ≠ 0 (| υ 0 | ≤ 1). Носитель гиростата может намагничиваться, порождая магнитный момент I p . При этом влиянием вихревых токов (токов Фуко) и образованием гистерезиса пренебрегаем.

Предполагается, что носитель гиростата обладает осью кинетической симметрии, совпадающей с главной центральной осью инерции гиростата Cx 3 . В силу этого имеем

I p = H - 1 n 2 s 3 е з ,               (3)

где H = | H |, n 2 ≥ 0 – постоянная, зависящая от магнитной проницаемости материала [1] носителя- парамагнетика и от характерной постоянной МП [4]; e 3 – орт оси Cx 3 .

Магнитный момент гиростата-магнетика представляется равенством

Aco + ( w x Ao ) + ( w x k ) = Q ( s 3) f ,   (6)

s + ( w x s ) = 0,               (7)

где f = [ - s 2, s 1 ,0] T .

Система уравнений (6)-(7) является аналитически замкнутой по переменным ω , s и обладает первыми независимыми алгебраическими интегралами:

J 1 ( o , s 3) - -( w T Aw ) - W ( s 3) = h , (8)

J2 (w, s) = (Ao + k) • s = h1,(9)

J3(s) - H2 = 1,(10)

где h , h 1 - постоянные интегрирования,

W (s 3) = | Q (s 3) ds 3 = n1 s 3 +1 n 2 s 32(11)

  • - характерная энергетическая (соленоидальная) функция [5].

Из многообразия возможных состояний гиростата, определяемого системой уравнений (6), (7), выделим множество движений, устанавливаемое системой условий:

0 ( t ) = 0 0 , ф ( t ) = a , V ( t ) = Q ,    (12)

(0 0 0 n ), ( 0 0, a, Q ) = const ^ 0.

Движение (12), если оно существует в МП, является регулярной прецессией (РП), совершаемой со скоростью Ω.

Поставим задачу: определить критерий существования РП, области её реализации, условия существования, а также свойства движения и области его устойчивости и неустойчивости. □

Модифицированная функция Лагранжа для данной задачи имеет вид

L = 1 ( w Ao ) + ( w k ) + W ( s 3),    (13)

где функция W ( s 3 ) (одномерный квазипотенциал) определяется равенством (11).

  • 2.    Существование и области реализации регулярной прецессии

Введём вектор углов Эйлера Φ ( θ , ψ , ϕ ) и структурно-динамические условия:

A 1 = A 2 = A , k 1 = k 2 = 0,        (14)

F ( Q , u 0) - Au 0 Q 2 - G Q- Q ( u 0) = 0,  (15)

где обозначено G = A 3 a 3 0 + k 3, u 0 = cos 0 0.

Соотношения (14) выражают структурно-динамическую симметрию гиростата относительно оси Сx 3 , при которой гиростатиче-ский вектор-момент k коллинеарен этой оси.

Величина G является значением при t = 0 проекции вектора кинетического момента гиростата на координатную ось Сx 3 .

Для РП (12) имеет место следующая

Теорема. Для того чтобы движение гиростата (12) существовало, необходимо и достаточно, чтобы выполнялись условия (14), (15).

Доказательство. Необходимость.

Пусть выполняются условия (12). Тогда, в силу скалярных уравнений системы (6) и зависимостей вида ω ( Φ , Φ ) (системы кинематических уравнений Эйлера в осях Сx j ), получаем равенства, являющиеся тождествами по переменным θ , φ , из которых непосредственно следуют ограничения (14), (15).

Достаточность. Пусть выполняются условия (14), (15). Тогда, в силу одного из скалярных уравнений системы (6), получаем первый интеграл ю3(t ) = ® ° ,а из остальных уравнений согласно зависимостям ω ( Φ , Φ ) следуют равенства, которые должны быть тождествами по переменным θ , φ . В силу тождественности этих равенств при 0 <  θ < π и уср # 0 получаем ( 0 , у , ср )= const , т.е. условия (12). □

Получим условия реализации РП в силу соотношения (15). Если выполняются ограничения u0 ^ 0, D(u0) = G2 + 4Au0Q(и0) > 0, (16)

то, в силу условия (15), РП может быть в общем случае реализована со скоростями

Q=Q , + B о 7 D ( и о ) , (17)

где Q * = B о G , B о = (2 Au о ) - 1 . Здесь Q=Q , - первая критическая (по отношению к изменению знака параметра Ω) скорость РП, значение которой определяется действительным корнем уравнения д F /dQ = 0.

Согласно равенству (17) при D ( u 0 ) 0 могут быть реализованы РП с двумя скоростями: Ω = Ω 1 и Ω = Ω 2 > Ω 1 , существующими при

Q p = Q 2 + 2 B о Q ( u o ) 0.

Эти скорости соответствуют двум режимам – "медленной" и "быстрой" прецессиям. Если u 0 Q ( u 0) 0, то Q 1 Q 2 0 и РП совершаются со скоростями Ω 1 , Ω 2 в противоположных направлениях. При D ( u 0) = 0, когда Q p = 0, скорость РП является единственной (Q = Q , ).

Если u 0 = 0, то, согласно равенству (15), при G # 0 имеем Q = - n 1 G 1 и скорость РП однозначно определяется значениями параметров n 1 , k 3 , ω 30 .

Определим область L реализации РП. Обозначим p2 = (4A)-1G2, 31 = n2 - 4n2p2   (17*)

и пусть u 1 u 2 - действительные корни уравнения D ( u 0) = 0, существующие при 3 1 0. При 3 1 = 0 имеем u 1 = u 2 = u * , где обозначено u * = - (2 n 2) - 1 n 1 ^ 0, | n 1| < 2 n 2, причём для парамагнетика имеем n 2 > 0 [6].

Пусть δ1 > 0. На координатной плоскости (u0, v) зависимость v (u 0) = n 2 u 0 + n1 u 0 + p2          (18)

определяет параболическую границу области L с осью симметрии u 0 = u * . В силу условия (16) эта область определяется ограничением v ( u 0 ) ≥ 0 и на данной плоскости для значений u 0 ей соответствуют полуинтервалы (- 1, u 1 ], [ u 2 , 1), где значения u 1 , u 2 есть

( u 1 , u 2) = u * (1 + n 1 - 1 -3 1 ).

Если 3 1 0, то области L на данной плоскости соответствует полоса - 1 <  u 0 < 1.

Определим области знакопостоянности значений параметра Ω. Из интеграла (9) и дополнительного интеграла системы (6), (7) в силу условий (14) для РП (12) имеем

Q= [A(1 - uo2)]-1 f (uo),(19)

где обозначено f ( u ) = h 1 - Gu , u = cos 0 .

Условие существования (15) РП (12) может быть представлено в виде m1 + m 2 u 0 = 0,(20)

где обозначено m1 = n1 + GQ, m2 = n2 - AQ2,(21)

причём m J m 2| .

Можно показать, что уравнение (20) с параметрами (21), эквивалентное соотношению (15), определяет множество стационарных значений в = в0 е (0, п ), реализующих движение вида (12). Для этого в силу интегралов (9), дополнительного и лагранжиана (13), выраженных в углах Эйлера, составляется выражение для потенциальной энергии U ( θ ) приведённой по Раусу системы [7, c. 87; 8]. Уравнение для стационарной точки U '(0 0 ) = 0 определяет двупараметрическое множество решений, зависящее от параметров h 1 , G и соответствующее уравнению (20).

В силу равенства (19) изменение знака величины Ω при G ≠ 0 происходит в точке u 0 = h 1 G -1 = u p . При G = 0, u 0 ≠ 0 знак величины Ω сохраняется и совпадает со знаком h 1 . Интервалы 1 1 = ( - 1, u p ), 1 2 = ( u p , 1) зна-копостоянности величины Ω представлены в табл. 1.

Таблица 1. Области знакопостоянности величины скорости прецессии

G

Ω > 0

Ω < 0

G > 0

и 0 е 1 1

и 0 е 1 2

G < 0

и 0 е 1 2

u 0 е 1 1

3.    Достаточное условие устойчивости

Первые интегралы (8)-(10) при ограничениях (14) принимают вид:

2 J 1 = A ( ю2 + ю 22) - 2 W ( s 3 ) = 2 h ,

J 2 = A (ю 1 s 1 + m 2 s 2) + Gs 3 = h 1 ,     (22)

J3 = s2 + s 2 + s 32 = 1, а дополнительный интеграл, существующий при условиях (14), представляется в виде

J ^ A 3 a3 + k 3 = G .            (23)

Применяя прямой метод А.М.Ляпунова [9], составим характерную функцию задачи, согласно равенствам (22), в виде

V ( to , s ) = V 5 + 2 J 1 - 2AJ 2 + AA2 J 3 , (24)

где V 0 - произвольная постоянная, а множитель λ имеет размерность угловой скорости.

Обозначим p j Asj (j = 1,2).    (25)

В силу выражений (22)-(25) имеем

V ( P 1 , P 2 , s 3 ) = V 0 + A ( P 2 + P 22 ) +

+ n 3 s 2 - 2 m 3 s 3,

где обозначено n3 = AA2 - n2, m3 = n1 + AG, а множитель λ подлежит определению.

Из условий стационарности функции V (26) по всем переменным для движения гиростата в режиме (12) имеем

Pj = 0 ( j = 1,2), u о = n 3 - 1 m 3 .      (27)

Последнее равенство (27) имеет место при 2 2 # О к , где

О k = V A -1 n 2               (28)

  • - вторая критическая (по отношению к изменению характера устойчивости РП) скорость прецессии гиростата.

Полагая λ = Ω, из последнего равенства (27) получаем соотношение (15); при этом n 3 = - m 2 , m 3 = m 1 .

Найдём условие устойчивости состояния (12). Принимая это состояние за невозмущённое, положим в возмущённом движении s 3 = u 0 - q , сохраняя в силу равенств (27) для остальных переменных прежние обозначения. Полагая V 0 = m 1 и 0, представим в возмущениях соотношение (26) с учётом условия (20) и последнего равенства (27)

V (P1, p 2, q) = A (P12 + p 22) - m 2 q2, где величина q принята достаточно малой.

Функция (29) положительно определённа при условии m2 < 0.(30)

Это ограничение в силу равенства (28) эквивалентно следующему

IQ > О к.(31)

Поскольку существование состояния (12) в МП как стационарного движения установлено, то, на основании метода Четаева построения функции Ляпунова [10, c. 21] заключаем, что стационарное движение (12) устойчиво по отношению к переменным p1, p2, s3. Отсюда, согласно интегралам (22), (23) и соотношениям (25), находим, что РП гиростата при данных предпосылках устойчива и по отношению к величинам ю12 + ю2, s12 + s2, ю3, s3.

Таким образом, ограничение (30) выражает достаточное условие устойчивости РП (12) по указанным величинам.

4.    Области устойчивости прецессии

Установим области реализации условия устойчивости (30). Введём характерную угловую скорость

Q r = - n 1 G - 1 ^ - A 3 1 k 3 ), (32) соответствующую значению u 0 = 0, а также интервалы 1 3 = (0, П 2), 1 4 = ( П 2, n ).

Определяющее соотношение (15) представим в виде

Q 2 Q k = 2 Q ( Q-Q r ), (33)

где Q , - величина, содержащаяся в равенстве (17). Тогда области выполнения условия (30) в силу равенства (33) могут быть представлены табл. 2.

Таблица 2. Области устойчивости прецессии

G

Q 0 в 1 3

Q 0 в 1 4

G >  0

Q>Q r

Ω < Ω r

G <  0

Ω < Ω r

Ω > Ω r

5.    Необходимое условие устойчивости

Оценим устойчивость РП по углу θ . Представим скалярные динамические уравнения (6) при условиях (14) в углах Эйлера и выделим из множества возможных состояний, определяемых этой системой, движение по переменной θ . Принимая движение (12) за невозмущённое и полагая в возмущённом движении θ = θ 0 p , где p – малое возмущение, произведём линеаризацию уравнений в малой окрестности значения θ = θ 0 .

В линейном по p приближении примем cos Q ~ s + rp, sin Q « r - sp, где (s, r) = (cos в0, sin θ0). Тогда в возмущённом движении согласно условию (15), представленному в форме (27), в линейном приближении получаем p - A 1 m2 r2 p = 0,            (34)

где m 2 определяется равенством (21).

Согласно уравнению (34) для устойчивости по углу θ движения гиростата в режиме РП необходимо, чтобы выполнялось ограничение (30), полученное ранее как достаточное условие устойчивости.

6.    Условие неустойчивости

Неустойчивость РП гиростата по углу θ определяется ограничением m 2 > 0, из которого непосредственно следует соотношение

I Q|к, (35)

противоположное по смыслу условию (31).

При условии (35) на плоскости ( p , p ) уравнения (34) траектории фазовой точки являются неустойчивыми, за исключением двух асимптотических траекторий - сонаправлен-ных полупрямых, проходящих через седловую точку (0, 0) [11, c. 46]. Эти траектории соответствуют особому режиму РП, который здесь не рассматривается, и являются сепаратрисами, разделяющими фазовую плоскость на четыре области, содержащие множества гиперболических траекторий.

Согласно соотношениям (33), (35) области неустойчивости РП гиростата представляются табл. 3.

Таблица 3. Области неустойчивости прецессии

G

Q 0 в 1 3

Q 0 в 1 4

G >  0

Qr

Ω > Ω r

G <  0

Ω > Ω r

Ω < Ω r

Здесь Ω r определяется равенством (32).

Таким образом, при переходе значений величины Ω через критическое значение Ω k , определяемое равенством (28), характер устойчивости РП гиростата изменяется на противоположный.

7.    Устойчивость вырожденной регулярной прецессии

Невырожденный режим РП гиростата в МП имеет место при 0 <  θ 0 π . При значениях θ 1 = 0, θ 2 = π РП вырождается в перманентное вращение вида

( ш , s ) = [(0,0, ю );(0,0, s )],         (36)

где ю = ю 0, s = s 0 = ± 1.

При движении в режиме (36) в граничных точках отрезка 0 ≤ θ π происходит дерегуляризация РП как стационарного состояния. Это приводит к вырождению циклических координат исходной системы уравнений и к связанному с ним вырождению соотношений (19) и U ( θ ) – потенциальной энергии приведённой по Раусу системы, которые имеют полюсы при s = ± 1. Данное вырождение в этом движении обусловлено совпадением ψ - и φ - осей, соответствующих этим координатам. В результате этого движение (36) по отношению к РП, для которой 0 <  θ 0 π , является вырожденным особым движением (термин работы [12]). В связи с этим следует отметить, что определяющее условие (15) имеет место только для невырожденной РП.

Получим достаточное условие устойчивости состояния (36). Полагая G ≠ 0, составим функцию Ляпунова в виде линейной связки интегралов

V = J 1 - V 2 - 1 Л 2 J 3 ,       (37)

где λ 1, λ 2 – множители Лагранжа.

В силу соотношений (36), (37) и условия стационарности функции V по s имеем

2 2 = - ( n 2 + ms ),           (38)

где m = n 1 + A 1 G .

Принимая движение (36) за невозмущённое, в возмущённом движении положим

( o * , s * ) = [( r i , r 2 , a), (q i , q 2 , s + q 3 )]

и представим характерную функцию – интеграл V (37) в новых переменных – возмущениях r ( r 1 , r 2 , 0), q ( q j ) ( j = 1, 2, 3) (с учётом определяющего соотношения для множителя Лагранжа (38)) в виде

(

V (r, q) = ^ I -p-^1 Qj I + n2sq3 + j =1 k 2          7              (39)

+ 1Ф ( s ) q 3 2 + 0 ( q 3 3).

Здесь обозначено ( j = 1, 2,)

P j ( r , q j ) = Ar j + Q ( s ) sq j

Q j ( r j , q j ) = Ar i- q j - “Gsq j ,

Ф ( s ) = A 1 Gs + Q ( s ) s - n 2.

Для квадратичной формы (39) применим теорему Финслера [13]. Искомое достаточное условие устойчивости сводится к требованию положительной определённости квадратичных форм P j на множестве значений

Q j = 0   ( j = 1,2).           (40)

Согласно (40) при n2 G ^ 0 имеем qj = 2AG -1 srj (j = 1,2).        (41)

Полагая 21 = Qr, получим Ф (s) = 0. Из интеграла J3 (10), выраженного в возмущениях rj, qj (j = 1, 2, 3), имеем q3 = - s ± 71 - p2   (p2 = q2 + q2). (42)

Из равенств (42) следует, что для возмущений таких, что | q31 < 1, при s = 1 имеем q3 < 0, а при s = – 1 получаем q3 > 0. В силу этого, требуя выполнения условия n2sq3 > 0 для данного слагаемого, содержащегося в равенстве (39), положим n2 < 0.                    (43)

Внося выражения (41) в квадратичные формы P j , получим

P j = Arj [1 + 4 AG - 2 Q ( s ) s ] ( j = 1,2). (44)

Итак, согласно теореме Финслера и в силу соотношений (44), для устойчивости состояния (36) достаточно, чтобы выполнялось условие (43) и, кроме того, чтобы

D ( s ) = G 2 + 4 AQ ( s ) s >  0.        (45)

При этом имеет место равенство

^ 2 = - n 2 , cуществующее согласно соотношениям (38) и λ 1 = Ω r .

Ограничение (45) заведомо выполняется при

Q ( s ) s 0,                 (46)

и тогда соотношение (46) приводится к виду n1 s + n2 > 0.                (47)

Неравенство (47) является усилением условия устойчивости (45) и выполняется лишь в случаях, при которых величины n1, s имеют одинаковые знаки. В частности, если n1 > 0, то устойчивость состояния (36) достигается лишь при s = 1, а в случае, при котором n1 < 0 – лишь в положении s = – 1. При этом знак параметра n1 совпадает со знаком косинуса угла между вектором магнитного момента и осью Cx3.

Следует отметить, что условие (43) соответствует случаю, при котором материал носителя гиростата является диамагнетиком в отличие от парамагнетика , для которого, как указано выше, n 2 > 0 [1, c. 295].

Область устойчивости состояния (36) на плоскости параметров G , s определяется условием (45), а её граница R - уравнением

D ( s ) = 0     ( s = - 1,1).         (48)

Получим достаточное условие устойчивости данного состояния на границе R , являющейся в силу уравнения (48) гиперболой. Этому уравнению соответствует вырожденное особое состояние гиростата, принадлежащее множеству (36).

Интегральная связка (39) для данного вырожденного состояния в силу условия (48) принимает вид

V ( r , q ) = (2 A ) - 1 ^ [ A 2( r j - . q ) - j = 1

I G - .-.A ] 2 q j ] + n,sq 3 + (49)

+ - ф ( s ) q з 2 + 0 ( q з 3).

Требуя положительной определённости квадратичной формы (49) с учётом возмущений до второго порядка включительно, выделим множество значений, для которого

.= (2 A) -1 Gs.(50)

Тогда, согласно равенствам (42), (50), квадратичная форма V (49) примет вид

  • V(r, q) = (8A)-1 ^ (2Arj- Gsqj) — j =1(51)

  • - 1 n 2 P 2 + 1 ( P 2 - n 2 ) P 4 + 0 ( P 6),

22       8

где параметр p определяется формулой (17 * ).

Квадратичная часть функции (51), находящаяся под знаком суммы, является знакопостоянной величиной вследствие вырожденности движения. Однако положительная определённость функции V достигается путём наложения ограничения (43).

Таким образом, на границе R области устойчивости, определяемой уравнением (48), устойчивость состояния (36) также обеспечивается условием (43).

Заключение

Приведённая частная задача о существовании и устойчивости РП гиростата в МП логически связана с общей фундаментальной проблемой движения относительно центра масс системы твёрдых тел в стационарном силовом поле.

В данной задаче установлена принципиальная возможность существования устойчивой РП гиростата в стационарном МП, моделируемом полем прямого магнитного диполя. При этом показано, что интервал изменения значений скорости устойчивой прецессии ограничен значением некоторой её критической скорости, зависящей от величины магнитной постоянной n 2 .

Некоторыми характерными свойствами обладает вырожденное особое движение гиростата. В частности, каждое из состояний множества (36) порождает два вырожденных расслоения инвариантных многообразий, являющихся вырожденными расслоениями с параметром ω .

Вырожденному состоянию (36) соответствует множество перманентных вращений гиростата вокруг его главной центральной оси инерции. Максимальным порождающим интегралом для этого множества является полная линейная связка интегралов (22), (23) системы уравнений вида

V = J + £ J

J = 1

где λ j – неопределённые множители Лагранжа. Эту интегральную связку можно трактовать как некоторое линейное пространство первых интегралов над собственно евклидовом пространством параметров, выбирая данные интегралы за базисные.

Следует отметить, что состояние (12) при выполнении условия (31) является не только устойчивым, но и равномерно устойчивым. Действительно, для системы уравнений (6), (7), отнесённой к области пространства ( ω , s ), являющейся открытым связным множеством, имеет место теорема о равномерной устойчивости [14, c. 20]. Эта устойчивость непосредственно следует из установленной устойчивости РП как стационарного движения.

Список литературы Устойчивость регулярной прецессии гиростата-магнетика в магнитном поле

  • Тамм И.Е. Основы теории электричества. М.: Наука, 1966. 624 с.
  • Яновский Б.М. Земной магнетизм. Л.: Изд-во Ленинградского ун-та, 1978. 592 с.
  • Белецкий В.В., Хентов А.А. Вращательное движение намагниченного спутника. М.: Наука, 1985. 288 с.
  • Садов Ю.А. Быстрое вращение спутника с магнитным демпфером. 1//Космические исследования. 1970. Т. 8, № 4. С. 547-556.
  • Макеев Н.Н. Интегрируемость гиростатических систем в магнитном поле//Проблемы механики и управления. Нелинейные динамические системы/Пермский гос. ун-т. Пермь, 2003. Вып. 35. С. 49-70.
  • Вонсовский С.В. Магнетизм. М.: Наука, 1971. 1032 с.
  • Меркин Д.Р. Введение в теорию устойчивости движения. М.: Наука, 1976. 320 с.
  • Парс Л.А. Аналитическая динамика/пер. с англ. М.: Наука, 1971. 636 с.
  • Ляпунов А.М. Общая задача об устойчивости движения. М.; Л.: Гостехиздат, 1950. 471 с.
  • Рубановский В.Н., Самсонов В.А. Устойчивость стационарных движений. М.: Наука, 1988. 304 с.
  • Обморшев А.Н. Введение в теорию колебаний. М.: Наука, 1965. 276 с.
  • Иртегов В.Д. Особенности многообразия стационарных движений. Пробл. аналитич. механики, теорий устойчивости и управления. М.: Наука, 1975. С. 154-160.
  • Кузьмин П.А. Малые колебания и устойчивость движения. М.:Наука, 1973. 207 с.
  • Руш Н., Абетс П., Лалуа М. Прямой метод Ляпунова в теории устойчивости/пер. с англ. М.: Мир, 1980. 302 с.
Еще
Статья научная