ИССЛЕДОВАНИЕ ДИНАМИКИ ПРОСТРАНСТВЕННОГО ДВИЖЕНИЯ ТЕЛ ПЕРЕМЕННОГО СОСТАВА Текст научной статьи по специальности «Физика»

Электронный журнал «Труды МАИ». Выпуск № 41

УДК 531.36

www.mai.ru/science/trudv/

Исследование динамики пространственного движения тел

переменного состава

М.М. Крикунов

Аннотация

Исследование пространственного движения твердых тел постоянного и переменного состава остается одной из центральных проблем динамики твердого тела и их систем, а также имеет важное значение для прикладных задач механики космического полета.

В работе проводится исследование твердого тела переменного состава с заданными законами изменения моментов инерции. Выводятся уравнения движения твердого тела переменного состава на основе формализма Гамильтона. На основе выведенных уравнений строится математическая модель движения твердого тела переменного состава в канонических переменных Андуайе-Депри. Для анализа динамики проводится аналитическое и численное моделирование движения тела. Показывается существование хаотических режимов в системе на основе метода Мельникова-Виггинса.

Результаты работы могут быть использованы для анализа и синтеза режимов углового движения космического аппарата на активных участках с изменением массы.

Ключевые слова

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

Введение

Задача моделирования движения твердого тела переменного состава была и остается одной из центральных задач динамики твердого тела [1-7]. В зависимости от того как представляется твердое тело существует множество моделей и методов описания его движения. Одним из самых эффективных методов исследования является формализм Гамильтона [1], позволяющий записать 2п (где п — количество степеней свободы механической системы) дифференциальных уравнений первого порядка. В данной работе

строится модель тела переменного состава на основе формализма Гамильтона в канонических переменных Андуайе-Депри. Так как гамильтониан тела переменного состава, записанный в этих переменных, содержит только одну позиционную координату, это позволяет существенно упростить исследование динамики движения твердого тела.

В работе рассматривается движение твердого тела для различных случаев изменения состава тела. Переменность состава моделируется заданием законов изменения инерционно-массовых характеристик.

Построение математической модели. Общий вид канонических уравнений движения

тела переменной массы

Для построения математической модели используются уравнения Гамильтона, записанные для твердого тела переменного состава [1]:

дИп дИп

дра , (1.1)

где , — переменные действие-угол; р — обобщенная реактивная сила,

отнесенная к ° -ой обобщенной координате; — внешняя обобщенная сила, отнесенная к

координате .; — функция Гамильтона, которая может быть записана в переменных действие-угол. В качестве таких переменных могут быть выбраны переменные Эйлера, Крылова и др. В данной работе формализм Гамильтона применен для переменных Андуайе-Депри [2].

Запишем гамильтониан свободного тела переменного состава в общем виде:

Нп = Т = — ¿51(5

2 , (1.2)

1 =

где 5 = (Р % г) — вектор угловой скорости; у х у "у — компоненты тензора инерции в связанной системе координат.

В данной работе рассматриваются случаи, когда внешние силы, действующие на тело переменной массы, отсутствуют, но абсолютные скорости отбрасываемых частиц не равны

П 14 О = ° Р * 0 нулю, следовательно, в уравнениях (1.1) , ° .

Запишем в общем виде обобщенную реактивную силу [1]:

(I I I Л

■■ ху Х2

I I I

ху УУ У?

I I I

у Х? У?

ра=±^йdvV

у=\ <31 (1 3)

где т - масса у -ой отбрасываемой частицы; - абсолютная скорость отбрасывания

V -ой частицы.

Случай 1. Рассмотрим частный случай, когда скорость отбрасываемых частиц является линейной функцией скорости:

й = К( t) v

Подставим (1.4) в (1.3)

V -WJ' V (1.4)

^=É *V (') dH =£-± (t) drf

v=i dt dqa v=i dt dqa dqv v=i dt 2 (1 5) Введем в рассмотрение функцию:

n = ¿ V t)

"' ' ,(1.6)

V=! dt 2

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

Р = —

Таким образом, если закон изменения массы известен, то задача сводится к

V ~ г^т

отысканию вектора г, который имеет вид [2]:

Я =ахру

ё = ( p q г) _ Mimn 11гт„пй V .

(1.7)

)T

где - / — вектор угловой скорости; ' V — радиус-вектор к -ой

отбрасываемой частицы.

В переменных Депри компоненты вектора угловой скорости записываются следующим образом [3]:

p = — VG2 - — sin l, q = — VG2 - — cos l, r = —

ABC (1.8)

dVv

Для нахождения производных по обобщенным скоростям в (1.5), выразим

вектор v через производные углов Депри:

Л y a\ (/, t ) a2 (/, t ) B^ aj2 (/, t ) a2 (/, t ) с aL(l, t )

(1.9) где

1 sin2 / cos2 / /, ч sin2 / cos2 /

a

/ , \ 1 sin l cos l i 1 \ sin l

(/'t ) = CÇ) " Л(7) " ~ВЩ a(/'t ) =+ B (t )

(1.10)

Таким образом, для частного случая, когда справедливо (1.6), задача о записи канонических уравнений Гамильтона в переменных Андуайе-Депри сводится к отысканию

v

вектора _.

Случай 2. Рассмотрим теперь более общий случай, когда скорость отбрасываемых частиц U_ в формуле (1.3) не является линейной функцией скорости, а имеет вид:

U = v + ^ (1.11)

V V TV

V V

где " вектор относительной скорости движения частиц; ™ - относительная

скорость отбрасываемой частицы.

Подставим выражение (1.11) в (1.3), получим:

Р у у -д^

„=— Л " д%а „=! Л ГУ д4а (1 12)

где ^ будем считать заданными, так как нам известно, как именно отбрасываются частицы.

Введем обозначения:

„=— Л дс[а (1 13)

"dm_ - дч

г°2 2 ^ VT

1 dt r_ dqa (114)

Тогда (1.12) примет вид: P = P + P

' a a1 a2

(1.15)

Подставим теперь в (1.13) выражение для вектора скорости (1.7):

dm,,

i

dt

Так как

да - l

f ^ \

да

dq t

J

(116)

V

(а хД)

да

Л

= а •

У

то (1.13) принимает вид:

да

КдЧо

• Ру ~Ру

Ру

да

дд,

йт , 3 3 \

Ро'=Т ИТ(а хр')

да

[ддо

'Р у

т

йту I йг

а

о У У

да

(1.17)

дд,

— •Ру ~Ру

з да

Ру

Х\

о УУ

(118)

Расписывая (1.17) и приводя подобные слагаемые, перегруппировывая члены получим:

(

а

да я 2 я

— Ру ~Ру

дЯ,

\ 1о

дС1аУУ

Р(Л2+)+(х2+х2)+г^ (х2 + )■

дЧа дЧ* дЧ*

дд

дг

^ др дд^ д—+ р

ддо ддо

X V

УУУ

^ др дг ^ г —+ р

[

о У

ддо дд

Ху

дд дг

г—т~ + д

V 1о

оУ

ддо ^

Уу*у

V

оУ

(119)

Подставим (1.19) в (1.18):

Рг1 =Т

у=1

йту йг

■а •

да з 2 з

— •Ру ~Ру

удЧо 2 , 2 '

Ру

да

ддо

Х\

V У У

др -А йту / 2 2 \ дд -А йту ( 2 2 \ дг А йту ( 2 2 \

= рТ~г( У1 + *2) + д^т- Ч X2 + *2)+г— Т~г( + V2 )

дда 7=1 йг у ' ддо у=1 йг 4 7 ддо йг 4 7

^ др дд ^

д—+ р —

V ^о

т

оУ

йту 1 йг

X V

у У у

др дг + р-

ддо ддо

л

т

йту 1 йг

дд дг —- + д-

ддо дд<

т

о

йту 1 йг

-Уу *у

т

т.

V ^Чо ^ЧоУ

Так как тело является твердым, то относительные смещения его точек отсутствуют

й Ру

йг

= о

, следовательно

Т Р £1= й (м Рс)

, (1.20)

м = м (г)

Заметим, что в силу (1.20) имеет место следующее равенство:

м = м Рп = Р (г)

где 4 у — масса тела; ^с — радиус-вектор центра масс.

I йт„

-X,

й

й

^ ^ уУу=Т "Г (туХуУу) = Т т ХуУу = 1 ХТ =1 йг у=1 йг аг у=1 (121)

С учетом (1.21):

п

п

V

р = + + -

о\ Р 1 XX ъ д 1 и ъ ' 1 ZZ

дЧа дЧа дЧа

(

\

др да

д + р —

V дЧа дЧеУ

(

I -

1ХУ

др дг г+ р —

V ддо ддоу

л ( дд дгл

I -

1 Х2

г--ъ д—

V ддо ддоу

Запишем последнее выражение в матричной форме:

Ра— =¿4

Т} д(

ддо (1.22)

Заметим, что если главные оси инерции остаются неподвижными в связанной системе координат, то

дР л , „ д<

дг

Ра! = Р~А ъ д—в + г —с

дС

дЧа ддо

. (123)

Обозначим компоненты вектора относительной скорости отбрасывания частиц

V.

V =(v1 V2 V3 )Т

V V V V '

Вычислим производные обобщенным скоростям:

дVV

дЧо

¿V/ д , 5 у Ч д( у у

-^ =-(5 х Ру) =-х Ру+ 5

ддо ддо ддо ддо

в (1.14), для этого продифференцируем (1.7) по (124)

дА

Так как радиус-вектор точки явно не зависит от обобщенной скорости, то Следовательно

0

ддо

ду„ д / у у ч д( у

(¿ХРу)=~ХРу

ддо ддо дд

Подставим (1.25) в (1.14):

(125)

Р =ут V

Ро2 £ л Уп

д(

Л

ддо

'ХРу

У^Чо У

(1.26)

и выпишем соответствующие частные производные по обобщенным скоростям Андуайе-Депри:

д( _(дР дд дг"

~д1 "1"д7 И ¥

у / \Т

д( дР дд дг

дф2 ^дф2 дф2 дф:

у [ V

д( дР дд дг

2

дф3 \дфъ дфъ дф-_

з У

(1.27)

Допущение 1. Предположим, что все частицы отбрасываются в одном и том же направлении. Тогда вектор скорости отбрасывания частиц будет одинаковым для всех точек тела переменного состава:

V = V, уъ)т .

Выражение (1.26) принимает вид:

Ра2 Е & п

да

д¿1.

~ХРу

V е

V Ю

1 &

да

д¿Ч.

~ХРу

(1.28)

Проводя преобразование:

р.,=V Е

&ш„

1

да

д<1.

~ХРу

У

? Е

у=1

г да &ш„ у ^ 'Ру

^¿а &

= V ■

да

а V у

1 ~йТРу

(1.29)

и подставляя (1.20) в (1.29), получим:

Рет2 = V ■дах &- (М Рс ) = а2 д1 &г с)

у да & , ч

К • ^ ■ 1, (М рс )

(1.30)

С учетом (1.22) и (1.30) для обобщенной силы запишем следующее выражение:

р=т1да

дЧа

у да & , ч

V' ^' й (М рс )

(131)

Подставляя (1.31) в (1.1) получим УРАВНЕНИЯ ДВИЖЕНИЯ ДЛЯ ТЕЛА ПЕРЕМЕННОГО СОСТАВА:

¿¿а =

дЩ

дРа

Ра =

н

дЧа

+ ат I

1 да

дЧа

у да & , ч

к ■ ыТ • &(М рс)

(1.32)

Х

У

V

Основные допущения модели

Сформулируем допущения, которые будут использованы для исследования движения твердого тела переменного состава.

Допущение 1. Предположим, что все частицы отбрасываются в одном и том же направлении. Тогда вектор скорости отбрасывания частиц будет одинаковым для всех точек тела переменного состава:

V = ^1 V2 V))

Допущение 2. Будем считать, что центр масс остается на том же месте.

й

— (м рс )= М рс

йг

. Тогда (1.31) примет вид:

„ зт тТ -тт да ■ -г да - -т т да Р = ат 1Ь = ат I — + М •V •——х рс =ат I-

ддо

ддо

ддо

+М

- да -Vрс

ддо

(2.1)

где

- да -V' Рс

ддо

смешанное произведение векторов.

Уравнения движения с учетом этого допущения запишутся:

Чо =

Н

дро

ро =

дН

ддо

о Iда

ддо

+м

— да -VРс

ддо

(2.2)

Допущение 3. Пусть главные оси тела остаются неизменными в процессе изменения массы. Физически это может означать, что масса должна изменяться равномерно по всему объему тела, либо изменяться симметричным образом, в том числе с различными коэффициентами расширения (сжатия) по направлениям осей связанной системы координат. Данное допущение несколько сужает круг рассматриваемых тел. Например, оно не подходит для случая отделения частей тела произвольной формы.

Канонические уравнения движения тела переменной массы в переменных Андуайе-

Депри

Для исследования движения твердых тел постоянного и переменного состава целесообразно использование переменных Андуайе-Депри [3-4], характеризуемых связью с основными динамическими параметрами.

Рисунок 1 - Переменные Андуайе-Депри

Построим уравнения движения (2.2) для переменных Андуайе-Депри. Находя соответствующие производные в формуле (1.31), запишем получившиеся векторы:

Г б1п / Ь 1

Ау/О2 - Ь аЬ Г б1п / О 1

да б1п / Ь А^О2 - Ь2 ао

д1 А^О2 - Ь2 аЬ да соб / О да дфз Г 01

1 дф2 Ву/О2 - Ь2 ао 0

V СаЬ V 0 (0)

(3.1)

где аЬ , а° вычисляются по формуле (1.10).

Используя компоненты вектора угловой скорости (1.8), вычислим необходимые скалярные и смешанные произведения в формуле (1.31). Пусть связанная система является

главной центральной, тогда = (0 0 0). При условии неизменности главных осей обобщенные реактивные силы запишутся:

*=Ь

а

р = О

Ф2

Ф а

1 ■ б1П2 / ■

с2 - - в2

соб2/

В

Аб1П2 / ■ СОБ2 / ^

А2

-А + ■

В2

■В

. (3.2)

(3.3)

Гамильтониан тела в переменных Андуайе-Депри при условии неизменности главных осей инерции записывается следующим образом [2]:

Ь2

Н =

О2 - Ь2 Г 81П2 / 0082 Л

А

В

2 С

(3.4)

где А, В, с — главные моменты инерции; 1, Ь, О — переменные Андуайе-Депри. Вычисляя частные производные от функции Гамильтона и подставляя (3.2), (3.3) и (3.4) в (2.2) получим УРАВНЕНИЯ ДВИЖЕНИЯ ДЛЯ ТЕЛА ПЕРЕМЕННОГО СОСТАВА В ПЕРЕМЕННЫХ АНДУАЙЕ-ДЕПРИ:

/ = Ь

Ь = -

1 б1П / СОБ /

^ Л

С

В

ф = 12

2 /Л

б1П / СОБ /

В

Фз = 0, 1з = 0,

1 (О2 - Ь)(

1 1 1 • о, ---Б1П2/ +

А В)

Ь

а

V ^ ~ У

( 1 л б1П2 / СОБ2 О —тС--— А--— В

С2 А2 В2 У

о=0.

а

б1П / СОБ /

2/ >

А+

V

В2

В

Из (3.5) следует, что

ф = еотг /3 = еотг

Заметим, что в случае, когда моменты инерции являются постоянными, то уравнения для тела переменного состава вырождаются в в уравнения для тела постоянного состава.

Исследование движения твердого тела переменного состава

Рассмотрим ряд частных случаев движения твердого тела переменного состава.

Случай линейного изменения моментов инерции

Пусть моменты инерции тела изменяются по линейным законам:

А ф в ф С А = А - а,, В = В0 - Ь,, С = С0 - с,. Тогда уравнения движения (3.5) принимают вид:

1 = - Ь

^ вт21 сов21 ^

А - а, В - Ь

+ -

ь

С - С,

ь=-2 (122 -ь)[ Х-

1

Л

ь

а.

А - а, В - & у 8т21 -

а

ф2 = О

21 -Ь

^ 8Ш2 I СОВ2 I ^

А - а, В - Ь

V о

О = -О

( Со - С, ) ( Ао - а, ) 8Ш2 I + -

(Во - Ь,)

-СО82I

а

Л

а

-СО82I

У

( А - ах) (В - Ы )

Если при этом тело обладает динамической симметрией:

А = В = А - аг, С = С - с, Тогда уравнения (4.1) записываются следующим образом:

(4.1)

! = —

ь

ь

Ь =

А - а, С0 - С

ь

, 9 =

О

О = -

аО

А - а,' А - а,'

^ С - с, А - а^ а—--с

(А - С ) - (а - с), ^ А - а, С0 - с,

У

(4.2)

Второе и третье уравнение имеют аналитические решения:

О = ( А - а,) ■ сот,, 9 = сот, ■, + сот,2 Случай квадратического закона изменения моментов инерции

Пусть тело имеет форму цилиндра и ориентировано таким образом, что ось 2 связанной системы координат совпадает с осью цилиндра.

т-т ~ Н = Н - к,

Пусть высота цилиндра меняется линейно: о и пусть при этом центр масс

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

<

цилиндра равномерно распределена по всему объему и остается постоянной, тогда имеет

М = рБН (,) М = РН = -Р^к = сот,

место следующий закон изменения массы: , тогда .

Осевые моменты инерции относительно главных осей инерции, проходящих через

центр масс цилиндра, и их производные имеют вид:

M f

R

4 V

pSh(

4

H

2Л

C =

MR2

. (4.3)

(R2 +(H0 - ht)2) C

pSHR2

-= const

2 (4.4)

В этом случае уравнения движения запишутся:

; L L . G ■ L

l =---1—, < =—, L =-

A C A A - C

C

G

— A--C I, G = —A C A ) A (4.5)

Подставляя (4.3) и (4.4) в (4.5) получим окончательных вид уравнений:

l =

2L

г

1

pS(Hо -htR2 3R2 +(Щ-ht)2

L =

V

Lh

(Р2 =

12G

pS (H0 - ht)(3R2 + (H0 -ht)2)'

G = -

(H0 - ht)((H0 - ht)2-3R2) 3h ( R2 +( H--ht )2)

9R2 -(Щ-ht)2-

18R 4

3R2 +(H0 - ht)

( H0-ht)( 3R2 +( H0-ht )2)

G

(4.6)

Решение данных уравнений возможно только численно.

Динамика движения твердого тела при изменении одного момента инерции по гармоническому закону

Пусть законы изменения моментов инерции имеют следующий вид, характерный для описания малых упругих колебаний [4-6]

A = A (1 + £ sin Qt) B = B0 = const C = C0 = const ^ ^

Подставляя (4.7) в (3.5), получим следующие уравнения движения:

L sin21

l = LaL, <p2 = GaG, G =

L=-2 (12L)

aG A2 (1 + £ sin Qí):

-sA Q c°s Qt,

Л

A (1 + s sin Qt) B(

sin 2l

L

sin2 l

о У

aL Al (1 + s sin Qí)'

-sA Q c°s Qt

<

1 sin21 eos21 _ sin21 cos21

C0 A (1 + ^sinQt) B0 12 = A (1 + 6sinQt) B0

(4.9)

Выделим в системе (4.8) невозмущенную и возмущенную части, раскладывая правые части уравнений в ряд по степеням 6 и отбрасывая члены второго и более порядков малости:

l = f (l, L) + 6g, (l, L, t), ф2 = Д (l, G) + 6gñ (l, G, t), <L = fL (l, L, G ) + 6gL (l, L, G, t), G = 6gr (l, G, t).

L , (4.10)

(1 т 4 sin21 T . _ , ч g, (l,L, t) =-L sin Qt (i r\~r

Jl (l' L) = LaLnopoxd Ao fV2 G) = О(Опорожд

где , , ,

1 P

gjl, G, t)= G sin Qt fL (l, L, G ) = —1 (i2 — L )1

A 2 v A Bo j

sin 2l

gL (l, L, G, t) =1 (G2 - L2 ) sin 2l—sin Qt--1-Q L sin21 cos Qt

2 Л ^Lnopoxd. A

gG (l, G, t) =-1-—G sin21 cos Qt

XОпорожд. Л0

1 sin21 cos21 sin21 cos21

X =-----X =--1--

Lпорожд. ^ i г, Gпорожд. i

Co Ao Bo Ao Bo

Анализ хаотизации возмущенного движения на основе метода Мельникова-Виггинса

Для анализа хаотизации динамики систем [4-7] в последнее время активно применяется метод Мельникова-Виггинса [8-11], который заключается в определении расстояния между возмущенным устойчивым и неустойчивым многообразием седловых точек гомоклинических траекторий. Если функция Мельникова имеет простые нули, то возмущенная устойчивая и неустойчивая ветви имеют точки пересечения, что свидетельствует о наличии в системе локальных хаотических режимов движения в окрестности сепаратрисы.

В общем виде функция Мельникова строится следующим образом [5]:

M(t), b) = J [(DxЯ, gx) + (DgЯ, gG)](x* (t), G*, t +1)0dt

—ад

Dg я (г(О), G), J" gG (x* (t), G*, t +1) )dt

Г(О) - О н = о ^

где г 4О ' — решение системы уравнений * . Для нашей системы это точки

(о;0) и (^;о)

Следуя формализму Виггинса, запишем функцию Мельникова (5.1) для системы (4.10):

(о*2 -Ь2)|—

М (,о) = Т{ 1 (О*2 - Г2)(

Ао В,

2*

б1п2/ -Ь Б1П0(, + ,0) +

о У

Ао

Т*

+ь а

Ьпор.

1 (О *2-Г2) 8т2Г —в1п о(, + ,0) 1 °

Ао

а* пор. Ао

Г б1П2 /* СОБ 0(, + ^ )

+

+

а*

а

1 о

Опор. „* лО* ^П2 1 * СО8 , + О

1 о

аОпор. Ао

Оаопор. (I = о; ь = о)|

-

у

-ш аОпор. Ао

О* б1П2 /* СОБ о (, + ,0 ) Ж =

Ао

1 1 171

ч А В V 2

V о о у-ш

11 (о*2 - Г2) ь МП 21* МП2 Г о (, + ,0) Ж +

+

1 1

| 1 ( О*2 - Г2) ь 8т2/а-вт 0(, +1,)-Ъа^ *

Ао аьпор. Ао

Ьпор. ° ь ^ 1" СОВ 0 (, + ,о )

Ж +

+— I О*2 б1п2 /* СОБ о (, + ,0) Ж, - 0аОпор (/ = о; ь = о) — | —I— О* б1п2 /* соб О (, + ,0) Ж,

Ао -ш Ао —ш аОпор.

о+ш 1

где

Ох = VО2 - ь б1П /, О = VО2 - ь2 СОБ /, О = ь,

следовательно,

. 2 7* Ок

Б1П2 / = Х

*2 *2

2 7* О7

СОБ / =

*2 *2

(5.2)

1 -2 7* 2 7*

1 Б1П / СОБ /

аьпор ~

Со Ао

• 2 7* 2 7*

* Б1П / СОБ /

а, =--1--

/2пор. А Вп

Перепишем (5.3) с учетом (5.2):

1

а

ьпор.

Л

V Оо Д Д> у

(5.3)

а

12пор. С}*2

( г<*2- Г<*2 Л

Ох , Ог

V

Ао Во

о У

(5.4)

* /-к 7* /-к У-Г* • 7* 7* 2^7^707

■ Б1П2/ = 20* Б1П / СОБ / ^ х 7 2

Т*2 = О*2

ао (о,о) = ^

Во

Окончательно функция Мельникова будет иметь вид:

1

М ^ ) = л

1 1

1

п т по

*2

л+—л2—4+■

4) Во у А> А А)

-14 -по*215

(5.5)

где

л. Ч

*3 * *

охогог . *2 81П

п( г + Го ) Л 12 = -1 121 --1122 --1 Л23

С) Д #0

Л21 = | 0*х0*01 вШ п (г + го ) Ж Л22 = |

ОхОо^в1п п( г + Го) л = л

> Z

Л 23 = /

п(г+(о) л

Л3 = 1

/—г *2 *2

у 07

о*

а*2°*2 собп(г + (о)Ж 14 = 1 0.*оХс*2 собп(г + (о)л

Л =

1 +ш 1

— Г — 1Й ^ *

°х

А Во -а«.р. о*2 - о?

_ ^*2 (^7*2

соб п(г+г0) Жг

А Во

а

'12

Во°х +

(5.6)

С учетом (5.6) перепишем интеграл

о*2

^5 = |

-"х

соб п(г + г0) Жг

-» Во°х + ^оС1!-

Известны классические решения для гетероклинических сепаратрисных полодий на

в*

сфере кинетического момента твердого тела

О* - (_;[)Кк-1)/2] * /«3 - ^2

Бес

а3 - «

■Ь (п2г)

о* =(-1) о^Ь (п2г)

о* ={-1)[к/ 2] о

а2 - а,

а3 - а1

- БесЬ (п2г)

(5.7)

к = 1,2,3,4; П2 = ^(а2-а1)(а3-а2); а1 = -1; а2 = -1; а3 = 77-;

А Во Со

Подставим в (5.7) значение к =1, соответствующее одной из ветвей гетероклинической сепаратрисы:

1

*

X у a3 - a

sech (n2t), G* = G*tgh (n2t), G* = G*

«2 - a

a3 -a

1 sech (n2t).

(5.8)

L - L

Используя (5.8), преобразуем интегралы 11 15

sech4 (n2t) tgh ( n2t) .

*3 a - a sech4 (n

L = G

2 n2 ,

(a3 - ai) i - ^iZ^L sech2 (n2t) G*3n

J

121 =

sin Q(t +10) dt

a3 - a

(5.9)

a - a

J sech2 (n2t)tgh(n2t)sin Q(t +10)dt

(5.10)

L23

_a3 - a2

1 Л

G*3n2 +r sech2 (n21) tgh3 (n21) a3 - ai -J1 - sech2 (n2t)

a - a

sech2 ( n2t)

sin Q(t +10) dt

(5.11)

J

a - ax

a3 - ai -1 - sech2 (n2t)

sech2( n2t)

cos Q(t + t0 ) dt

a - a9 L =—-2

a3 - ai -L B a3 - a2 B0

a - a

sech2 (n2t) + Л tgh2 (n2t)

(5.12) ■cos Q(t + t0 ) dt

(5.13)

Вычисляя интегралы (5.9)-(5.13), получим:

cosh

( ( q(

Ii = G

*з a^ Ж

a - a • 1 ( жО

2 1 sinh —

I 2n2,

V 2 v

ж-2arccos. i«2——

(a - a

3 1 yy

Q

a - a 2n2 (a - a)

2 (a2- ai) cos (Qto),

(5.14)

121 =

G Отж

2 (a - a) nl sinh

V 2n2 У

cos

(Qto)

G*3n„ ( 1

123

—г23 — I cos Qt,

a - a v 4 4

(5.15)

2

где

¥ =--L

11 а2

2л

лО

п2 sinh

лО

V 2П2 ,

1 л

cosh

Я«

V П 2

1 л

cosh

О«

V П 2

1

а 2n\ sinh

л

Ол | 2а п.,

V 2 П2

j 2 sinh

1 ( 2лО

лО

V n2

23 а2 • , (лО | а n2 sinh

V 2П2

2а п.,

Ол

sinh

лО

V n2

Ол

П cosh

V

'Ол

V 2П2 у

2n2 sinh

(Ол| „ 3 . , (Ол 2Пз sinh

V 2 n2

V 2П2 JJ

+ -

л

2а n2 sinh

( лО

V n2

(О | 4 sin « cosh —«

V n2 J

+ 4sin2 « cosh

V n2 J

+ 3 cosh

V n2 J

+ 3 cosh

(О || — «1 •

V n2 J J

j =(;*2 (a3 - a2 )(a2 - a1 )

(a3 - a1)

1

Ол

а 2 ■ u(Ол n2 sinh -

V 2n2,

■ +

+

2л

2а^ n7^/га sinh —

2 sinh

(О (л /—I--arccosv а I

V n 2 V 2 J

^ (Ол|

cosh -

2n

J V 2 n 2 J

n

(5.16)

1

1 =

а3 - а2

а3 - а

2л

nj1

- а sinh

Ол

2 sinh

(О(л У—, (Ол — I--arccosv а

V n2 V 2 JJ

cosh

V 2 n2 J

n

соб

(5.17)

а - а 2

1 =— 2

л sinh

(_О_

V 2n2

Л

arccos у

J

а - аг

A n2^J 1 -у2 sinh

Ол Cos О 2n

(5.18)

С учетом (5.9)-(5.18) из общего вида функции Мельникова (5.8) следует, что она является линейной комбинацией гармонических функций с нулевым средним, что свидетельствует о наличии у функции бесконечного числа простых нулей. Это означает, что функция расстояния меняет свой знак бесконечное количество раз, и, следовательно, в системе возникает хаос в окрестности сепаратрисы.

Численный анализ хаотизации возмущенного движения

Для качественного анализа динамических систем одним из эффективных методов является метод Пуанкаре построения фазовых сечений [3]. Метод Пуанкаре дает сечение

фазового портрета системы, что позволяет судить о его структуре. На рисунках 1-2 изображены сечения Пуанкаре для случая с гармоническим изменением одного момента инерции. В качестве условия сечения для отображения Пуанкаре выбрано сечение по

sin (Oí/2) = 0

моментам времени, соответствующим целому периоду возмущения: v ' ' .

т, л п т-г A =5, Я =9, С =12, O = 1 s = 0 00001

Рисунок 1 - Сечение Пуанкаре при ^ ' 0 ' 0 s 0,00001

г> о п тг А =5, Вп =9, С =12, П = 1 е = 0 001

Рисунок 2 - Сечение Пуанкаре при ^ ' 0 ' 0 е 0,001

Как видно из полученных изображений в рассматриваемом случае с увеличением

малого параметра е возникают хаотические слои в окрестности сепаратрисы, что

подтверждает расщепление сепаратрисы.

Изучение движения твердого тела переменного состава с гармоническими моментами

инерции на основе форализма Гамильтона метода Мельникова-Виггинса осуществлялось

ранее в работах [4-7]. Принципиальное отличие результатов настоящей работы заключается в том, что построена принципиально новая математическая модель для анализа динамики твердого тела с учетом переменности его состава.

ЗАКЛЮЧЕНИЕ

В работе построена математическая модель пространственного движения твердого тела переменного состава на основе гамильтонова формализма. Записаны уравнения движения в канонических переменных Андуайе-Депри. Проведено моделирование движения твердого тела переменного состава с различными законами изменения инерционно-массовых параметров. На основе построенной математической модели изучен ряд частных случаев движения:

Движение динамически симметричного тела с линейными законами изменения моментов инерции;

Движение динамически симметричного тела с квадратичными законами изменения моментов инерции;

Движение тела с гармоническими законами изменения моментов инерции.

Для аналитического доказательства наличия хаотических режимов построена функция Мельникова и показано наличие нулей этой функции.

Для иллюстрации возможности возникновения хаотических режимов в динамике системы построены сечения Пуанкаре для возмущенного движения твердого тела с гармоническими законами изменения моментов инерции.

Результаты работы могут быть использованы для анализа и синтеза динамики углового движения космических аппаратов на активных участках орбитального движения.

Работа поддержана Грантом Президента РФ МК-1497.2010.8.

Библиографический список

1. Космодемьянский А.А. Курс теоретической механики. — М.: Просв., 1955. — 656 с.

2. Архангельский Ю.А. Аналитич. динамика твердого тела. — М.: Наука, 1977. — 328 с.

3. Борисов А.В., Мамаев И.С. Динамика твердого тела. Гамильтоновы методы,

интегрируемость, хаос. — Москва-Ижевск: Институт компьютерных исследований,

2005. —576 с.

4. M. Iñarrea, V. Lanchares, V.M. Rothos, J.P. Salas, Chaotic rotations of an asymmetric body with time-dependent moments of inertia and viscous drag, Int. J. Bifurcation Chaos 13 (2003)393-409.

5. Дорошин А.В. Эволюции прецессионного движения неуравновешенных гиростатов переменного состава // Прикладная матем. и мех.. — 2008. — № 3. — С. 385-398.

6. Асланов В.С., Дорошин А.В. Хаотическая динамика неуравновешенного гиростата // Прикладная матем. и мех.. — 2010. —№ 5, Т. 74. — С. 735-750.

7. Anton V. Doroshin Analysis of attitude motion evolutions of variable mass gyrostats and coaxial rigid bodies system, Int. J. Non-Linear Mech. 45 (2010) 193-205.

8. V.K. Melnikov, On the stability of the center for time periodic perturbations, Trans. Moscow Math. Soc. 12 (1963) 1-56.

9. S. Wiggins [1988] Global Bifurcations and Chaos: Analytical Methods (SpringerVerlag, NY).

10. P.J. Holmes, J.E. Marsden, Horseshoes and Arnold diffusion for Hamiltonian systems on Lie groups, Indiana Univ. Math. J. 32 (1983) 273-309.

11. Лоскутов А.Ю. Динамический хаос. Системы классической механики // Успехи физических наук. — 2007. — №9, Т. 177. — С. 989-1015.

Сведения об авторах

Крикунов Михаил Михайлович, Самарский государственный аэрокосмический университет имени академика С.П. Королева (национальный исследовательский университет), инженер научно-исследовательской части (17г-Р101-090), телефон +7-927-72965-28, e-mail:: Krikunov_MM@mail.ru.