Хаотические колебания космического аппарата с упругим вертикальным тросом Текст научной статьи по специальности «Физика»

УДК 629.78

ХАОТИЧЕСКИЕ КОЛЕБАНИЯ КОСМИЧЕСКОГО АППАРАТА С УПРУГИМ ВЕРТИКАЛЬНЫМ ТРОСОМ

© 2010 В.С. Асланов, Б.В. Иванов

Самарский государственный аэрокосмический университет

Поступила в редакцию 12.02.2010

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

Ключевые слова: космический аппарат, тросовая система, метод Мельникова, аналитическое решение, хаос.

ПОСТАНОВКА ЗАДАЧИ

Возможность применения тросовых систем для операций в космосе широко обсуждается [1,2], в том числе для доставки с орбиты полезных грузов на Землю [3,4] и для мониторинга Земли, с помощью опущенного на вертикальном тросе средства наблюдения. В этих задачах весьма важными являются вопросы, связанные с поведением космического аппарата (КА) относительно собственного центра масс [5]. В статье рассматривается движение вокруг центра масс КА с тросом, развернутым по местной вертикали, под действием гравитационного момента [6] и момента от силы натяжения троса. В зависимости от отношения моментов инерции КА, характеристик троса и массы концевого груза могут наблюдаться неустойчивые положения равновесия, а в силу действия периодических возмущений возникать хаотические явления, приводящие к значительным эвол юциям углового положения КА. Целью настоящей работы является изучение влияния механических характеристик троса и параметров КА на угловое поведение аппарата методами хаотической динамики [7].

УРАВНЕНИЕ ДВИЖЕНИЯ КА С ТРОСОВОЙ СИСТЕМОЙ

Рассмотрим механическую систему (рис. 1), состоящую из КА c центром масс в точке O, троса PP и привязной груз P2 . Будем считать, что

Асланов Владимир Степанович, доктор технических наук, профессор, заведующий кафедрой теоретической механики. E-mail: aslanov_vs@mail.ru Иванов Борис Всеволодович, аспирант. E-mail: boris063@mail.ru

сила натяжения троса мала по сравнению с гравитационной силой, поэтому ее влияние будем учитывать только при рассмотрении движение КА вокруг центра масс. Трос будем представлять как невесомый упругий стержень, длина которого изменяется согласно закону Гука. Масса КА существенно превышает массу концевого груза. Центр масс КА движется по круговой орбите. В плоскости орбиты введем систему координат Ох1 у1, в которой ось Ох1 совпадает с местной вертикалью, и связанную со спутником систему координат Оху2 , в которой плоскость Оху совпадает с плоскостью орбиты, ось Ох направлена по продольной оси КА. Запишем уравнения движения КА относительно центра масс в проекции на ось О2, перпендикулярную плоскости движения Оху (на рис. 1 орбита показана штриховой линией), в виде [5]:

Рис. 1. КА с тросовой системой

Ca = -3f p 3 (B -A) sin a cos a - ГА sin(a -ф), (1)

где a - угол между продольной осью тела и местной вертикалью, ф - угол между линией действия силы натяжения троса и местной вертикалью, A,B,C - главные компоненты тензора инерции тела в связанной системе координат Oxyz , Т - величина силы натяжения троса, А = OP1 (рис. 1), f - гравитационная постоянная, p - параметр орбиты, в = a - ф .

Движение концевого груза относительно КА описывается следующими уравнениями в полярных координатах [1]:

ф + а + 2r"V(ф + а) + 3с2 sin фcosф = 0, (2)

r -r [(ф + с)2 + С(3cos2 ф-1)] + Tm-1 = 0, (3)

где с = V = ffn p- 3/2, m - масса груза. Система уравнений (2-3) имеет две пары стационарных решений [1], одна из них отвечает устойчивому (вертикальному) расположению троса

ф = 0,п ; Т = 3a2mr , (4)

когда при ф = 0 груз движется ниже КА, а при ф = п - выше КА.

Трос будем представлять как невесомый упругий стержень, длина которого изменяется согласно закону Гука:

Т = E (rl -1 - 1), (5)

где E - модуль упругости, l - длина ненапряженного троса. Пусть груз движется строго по местной вертикали в окрестности положения равновесия (4). Тогда уравнение (3) с учетом (5) примет вид

- + Q 2r = Em -1. (6)

где Q2 = Em- 3с2 > 0 (для материалов, применяемых для космических тросовых систем [1]). Равновесное положение груза определяется формулой

r0 = E(E - 3me21)-11.

При начальных условиях: t0 = 0: r = r0,r = V0, решение уравнения (6) имеет вид

r = r0 + V0Q-1 sin Qt, (7)

Для того чтобы трос находился в растянутом положении, начальная скорость груза при прохождении положения равновесия не должно превышать следующего значения

V0 < 3с 2 m1/213/2 ( E - 3 ma21 )-1/2 .

Сила натяжения троса согласно (5) и (7) изменяется по гармоническому закону

Т = Т0 + Tm sin Qt, (8)

где Т0 = 3 ma2 IE (e - 3 ma21) , Т = V0 Em V2l-1/2 (E - 3a2 ml) Ш .

m 0 \ /

БИФУРКАЦИОННЫЙ АНАЛИЗ И СЕПАРАТРИСНЫЕ РЕШЕНИЯ

Будем рассматривать движение КА относительно центра масс с развернутым по местной вертикали тросом, тогда уравнение (1) с учетом (8) примет вид

a = - a sin a - c sin a cos a - s sin a sin Q t, (9)

где a =

АТ0

c =

B - A 3f

AT

с С р- <10>

При малом положительном параметре е , уравнение (9) можно к классу квазиконсервативных с малым периодическим внешним возбуждением. Невозмущенное движение ( е = 0 > отвечает условию (4) и описывается уравнением:

а = - а 8та- с 8тасо8а . (11)

В силу (10) коэффициент а > 0 всегда больше нуля, а знак коэффициента с зависит от соотношения моментов инерции А и В. Положения равновесия невозмущенной системы (11) определяются как корни уравнения

8та(1 + /со8а)= 0 , (12)

где у = са-1. Из (12) следует, что два положения равновесия отвечают значениям: а* = 0, п, а третье положение равновесия: а* е (0, п) существует, при условии:

>1. (13)

Если это условие не выполняется, то наблюдается только два положения равновесия, причем центр всегда находится в точке а* = 0, а седло в точке - а* = п.

Очевидно, что для у < -1 точки а* = 0, п являются седлами, а промежуточная точка

(-Г-)

а =± arceos^ -

центром (левая ветвь диаграммы на рис.2). Если у > 1, то наблюдается обратная картина (правая ветвь диаграммы на рис. 2). Отметим, что бифуркационная диаграмма (рис. 2) для отрицательных углов а получается зеркальным отображением относительно оси абсцисс. Гиперболические точки существуют, когда выполняется условие (13). В этом случае периодическое внешнее возмущение £ sin а sin Q.t в возмущенной системе (9) может стать причиной хаоса [6].

Для исследования поведения динамической системы (9) в окрестности сепаратрис воспользуемся методом Мельникова [7], который позволяет найти толщину хаотического слоя и пересечение устойчивых и неустойчивых многообразий. Для применения метода Мельникова требуется получить аналитические выражения для сепаратисных орбит, другими словами, найти решения уравнения невозмущенного движе-

/ Седло ._. _ г \ { Центр ¡1 ... 1 х ' +Л i 4.

Стали \ \ __________i_i________ : """ ЦСЛТР 1

-i oí

Рис. 2. Бифуркационная диаграмма невозмущенной системы (11)

ния на сепаратрисах. Вид решений зависит от величины параметра у=есТ1 и начальных условий движения. Будем рассматривать два случая, определяющих величину параметра у. Случай 1. Пусть параметр

у = ес_ > 1, (14)

тогда седло

а8 =± агооо8 ) (15)

лежит между значениями 0 , п (рис. 2), а в точках = 0, ±п - (16)

при у ^да aS

л/2.

Уравнение (11) имеет интеграл энергии

W(а) = E , (17)

1 - 2 —а

потенциальная

2

где W(a) = -acosa-^-cos2a

энергия, E - полная энергия.

Если полная энергия E > WS, где WS = W(as ) , тогда система совершает вращение в одной из внешних областей. Если E < WS , то в зависимости от начальных условий наблюдаются колебания в одной из внутренних областей. При равенстве E = WS движение происходит по

гетероклиническим траекториям, которые соединяют гиперболические точки 51 и 5_1.

Сначала рассмотрим область А0 (рис. 3). Найдем решение уравнения (11) на сепаратрисах (гетероклинические траектории). Разделим переменные в интеграле энергии (17), тогда получим

t =

da

2[W (aS) + a cos a + 2 cos2 a]

(18)

где W(aS) = -a cosaS--cos2 aS =— . Замены

2 2c

переменных:

x=tana/2 (19)

в интеграле (18) и последующее интегрирование приводит к следующим гетероклиническим траекториям [8]

a± (t) = ±2 arctan | tan—tanh — )

I 2 2 ),

<r± (t) = (a )± =±

Asina

cosh At + cosa

(20)

аС

имеются центры. Отметим, что центр аС = _п совпадает с центром аС =п . При а ^ _ п и при а ^ п скорости а совпадают, поэтому можно утверждать, что фазовые траектории замкнуты на цилиндрическом фазовом пространстве. Будем рассматривать эволюцию фазовых траекторий на цилиндрическом пространстве для ае [_п,п]. Выделим две области А0 и А1, разделенные двумя седлами 51 и 5_1 (рис. 3). Необходимо отметить, что область А1 "разрезана" по вертикали а = п, _п . Из (15) следует, что седло 51 принадлежит интервалу а8 е (п / 2, п) и

где Я = (е2 - a2) c~1/2 - действительное число при выполнении условия (14).

Для области A1 гетероклинические траектории имеют аналогичный вид [8]

а±(t) = п±2arctan| cot—tanh—1

- I 2 2 )'

а± (t) = (á )± =

A sin aS

cosh At - cosa

(21)

Рис. 3. Потенциальная энергия и фазовый портрет при a = 1, c = 4

а

2

Случай 2. Пусть параметр

Y = ca—1 <—1, (22)

тогда центр

aC =± arceos ) (23)

лежит между значениями 0 , п (рис. 2), а в седло находится в точке

aS = 0 .

В данной работе гиперболические точки aS = ±п ( у = ca= —1) не рассматриваются.

В этом случае имеем две гомоклинические траектории, которые расположены справа и слева от гиперболической точки aS = 0 (рис. 4). Ввиду симметрии орбит будем рассматривать только правую гомоклиническую траекторию. Выполняя вычисления, аналогичные (17)-(19), получим следующие сепаратрисные решения

а± (t) = ±2 arctan | а± (t) = (а )± = +

d

ч cosh At J '

2 Ad sinh At (cosh At )2 + d2

(24)

)1/2 axn

тельные числа при выполнении условия (22)

где A = (-a - c)12, d = (-a - c) a"" - действи-

ФУНКЦИЯ МЕЛЬНИКОВА

Толщина хаотического слоя характеризуется функцией Мельникова [7], которая для возмущенной системы (9) имеет вид

M± (t0) = -s J ст± [sin а± sin(Qt + Qt0 )]dt =

J-ад

= -s cos Qt0 J c± sin а± sin Qtdt, (25)

Рис. 4. Потенциальная энергия и фазовый портрет при a = 1, c = -4

где а± (t),ст± (t) решения для гомо/гетероклини-ческих орбит (20), (21) или (24). Интеграл из (25)

/•ад

I = -1 ст± sin а± sin Qtdt.

J-ад

(26)

Определяет амплитуду изменения толщины хаотического слоя. Найдем его абсолютную величину для трех орбит:

- двух гетероклинических (20) и (21),

- одной гомоклинической (24).

Подставим решения (20), (21) и (24) в интеграл (26) и приведем полученные выражения к безразмерному времени, соответственно:

2 га

10 = sln as J_,

sinhr1

"sinQiTdT ,(27)

(coshr1 + cos aS)

2 Г ад - .

I1 = sin asJ 7—г-1-:rsinQTTdT,(28)

J-" (coshr1 -cosaS) 'v '

sinh t1

12 = 2 X 22 J-

sinh 2t2

(cosh2 т2 + x22)2

sin Q T n dr2

(29)

где T1 = Alt, т2 = Ajt - безразмерное время,

Q j = QA"1, Q 2 = QA2-1 - безразмерные частоты упругих колебаний троса.

Отметим, что как следует из (10) частоты

A=( c2 - a2 )"2 С-1/2 and A2 = (-a - c )1/2

зависят от параметров тросовой системы и величин моментов инерции КА. Зависимости максимальной толщины хаотического слоя (27), (28) и (29) для гомо/гетероклинических орбит (20), (21) и (24) от безразмерных частот возмущающей силы Qj = QA"1 и Q 2 = QA2-1 показаны на рис. 5.

Численные вычисления интегралов (27)-(29) показывают, что при Q1; Q2 > 6 толщина хаотического слоя стремится к нулю (рис.5) и возмущенная система (9) совершает регулярные движения без гомо/гетероклинических пересече-

1 / \ L

1 / "Ч. : ■' Р :/ i. 1

. ' У С. ^Hí

. S „Т——

ю Q

Рис. 5. Зависимости максимальной толщины

хаотического слоя (27), (28) и (29) от безразмерных частот возмущающей силы 01 = ОД"1 и О 2 = 0.1-1

ний. Это означает, что при О2 > 6 возмущающая сила гзтазт О, связанная с упругими колебания троса, не оказывает влияние по поведение КА относительно центра масс.

ЗАКЛЮЧЕНИЕ

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

Работа выполнена при поддержке Российского фонда фундаментальных исследований (09-01-00384).

СПИСОК ЛИТЕРАТУРЫ

1. Белецкий В.В., Левин Е.М. Динамика космических тросовых систем. М.: Наука, 1990. 329 с.

2. Сидоров И.М. Об использовании тросовых систем для создания постоянно действующего транспортного канала в космическом пространстве / / Полет. 2000. Вып. 8. С. 36-39.

3. MisraA.K. Dynamics and Control of Tethered Satellite Systems // Acta Astronáutica. 2008. V.63.P. 1169-1177.

4. Williams P., Hyslopb A, StelzerM, Kruijff M. YES2 optimal trajectories in presence of eccentricity and aerodynamic drag // Acta Astronautica. 2009. V. 64. P. 745-769.

5. Асланов В.С. Колебания тела с орбитальной тросовой системой // Прикладная математика и механика. 2007. Т. 71, Вып. 6. С.1042-1048.

6. Гукенхеймер Д., Холмс Ф. Нелинейные колебания, динамические системы и бифуркации векторных полей. Москва-Ижевск: Институт компьтерных исследований. 2002. 560 с.

7. Мельников В.К. Об устойчивости центра при периодических по времени возмущениях // Труды Московского математического общества. Вып.12. 1963. С.1-56.

8. Aslanov V.S. Chaotic behavior of the biharmonic dynamics system // International Journal of Mathematics and Mathematical Sciences (IJMMS). ID 319179. 2009. P. 1-18.

THE CHAOTIC OSCILLATIONS OF THE SPACECRAFT WITH A VERTICAL ELASTIC TETHER

© 2010 V.S. Aslanov, B.V. Ivanov

Samara State Aerospace University

The motion about a centre of mass of a spacecraft with a vertical elastic tethered system under the action of the gravitational moment and small periodic tethered force at a circular orbit is studied. The mathematical model is derived using Lagrange's equations including tether vibrations and oscillations of the spacecraft relatively of the proper mass center. The paper contains bifurcation analysis, phase space study, and analytic solutions for separatrixes. The considered mechanical system performs chaotic motion near separatrixes under the influence of small disturbances. The Melnikov method gives a criterion for homo/heteroclinic chaos in terms of system parameters. Results of the study can be useful for the analysis of gravitational stabilization systems with space tethers and for studying the behavior of a spacecraft with a deployed tether. Key words: spacecraft, tethered system, Melnikov method, the analytical solution, chaos.

Vladimir Aslanov, Doctor of Technics, Professor, Head at the Theoretical Mechanics Department. E-mail: aslanov_vs@mail.ru

Boris Ivanov, Graduate Student. E-mail: boris063@mail.ru