Исследование плоского движения относительно центра масс спускаемого аппарата с тригармонической моментной характеристикой при входе в атмосферу Текст научной статьи по специальности «Физика»

УДК 629.78.015

ИССЛЕДОВАНИЕ ПЛОСКОГО ДВИЖЕНИЯ ОТНОСИТЕЛЬНО ЦЕНТРА МАСС СПУСКАЕМОГО АППАРАТА С ТРИГАРМОНИЧЕСКОЙ МОМЕНТНОЙ ХАРАКТЕРИСТИКОЙ ПРИ ВХОДЕ В АТМОСФЕРУ

© 2010 Е. В. Баринова, И. А. Тимбай

Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет)

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

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

Рассматривается плоское движение неуправляемого осесимметричного спускаемого аппарата (СА) относительно центра масс на начальном участке траектории спуска в атмосферу, когда можно пренебречь изменением скорости центра масс и угла наклона траектории, а также гравитационным моментом и аэродинамическим демпфированием. Исследуются случаи, когда в процессе спуска происходит изменение характера движения: вращательное движение переходит в колебательное, колебательное движение «скачкообразно» переходит в колебательное движение с другими амплитудными характеристиками. В [1 - 3] рассмотрены переходные режимы движения СА с синусоидальной и бигармонической моментными характеристиками. В данной работе рассматривается движение СА с зависимостью восстанавливающего момента от угла атаки, имеющей вид тригармонического ряда. Такая зависимость характерна для аппаратов сложной геометрической формы, например, типа «Викинг» [4].

Наряду с тем, как наличие второй гармоники в моментной характеристике обуславливает возможность появления одного дополнительного положения равновесия СА по углу атаки по сравнению со случаем си-

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

Плоское движение спускаемого аппарата с тригармонической моментной характеристикой относительно центра масс при указанных выше допущениях описывается следующим уравнением [2]:

а + а(а + Ь(г^т2а + с(г^т3а = 0, (1)

где а - угол атаки; а(г), Ь(г), с(г) - коэффициенты моментной характеристики; г - медленно меняющийся параметр, переменность которого связана с медленным изменением плотности атмосферы в процессе спуска.

Коэффициенты уравнения движения

(1), если зависимость плотности атмосферы от высоты аппроксимировать экспонентой, могут быть представлены в виде [3]:

а ■■

■ а0 г, Ъ = Ь0 г, с

а0 =

-™аЫРо , Ъ0 =-ЩЯІРо ,

V2

с0 =-mcSlPo^A, г = eb(t-t0), b = AVo|sin0o|,

где ma, mb, mc - постоянные коэффициенты, S - характерная площадь, l - характерный размер, A - поперечный момент инерции СА, V0 - скорость, 0О - угол наклона траектории, р0 - плотность атмосферы в начальный момент времени t = t0, l - логарифмический градиент плотности атмосферы по высоте.

Для выяснения общих свойств движения СА, описываемого системой (1), воспользуемся методом фазовой плоскости. Интеграл энергии системы в случае невозмущённого движения, когда коэффициенты a, b, с постоянны ( г = const), имеет вид:

(X 2

Y - f (X) = h,

2 4 3

f (a) = (a - c)cos a + b cos a + 3 с cos a .

(2)

Экстремальные значения функции f (a) соответствуют состояниям равновесия уравнения (1), т.е. особым точкам на фазовой плоскости. Из (2) следует, что в зависимости от значения коэффициентов a, b, с могут существовать две, три или четыре особые точки на отрезке [0, p ]:

a1 = 0, а2 = p,

-b - m -b + m

cosa3 =--------, cos a4 =-------------------, (3)

3 4c 4 4c (3)

где m = Vb2 - 4ac + c2 .

На рис. 1 изображена номограмма для определения области, к которой относятся коэффициенты моментной характеристики СА, и таблица для определения числа и характера особых точек на фазовой плоскости системы. Номограмма разделена на пять об-

a = 2b

ластей следующими границами: — -----------3 ,

a = 2b a b2

— =--------3 , — = —2 +1. В таблице в зави-

c c c 4c

симости от знака коэффициента с и номера области показано, какие положения равно-

№ области a = 0 а3 а 4 а 2 = p

c <0

1А, 1В Седло Центр Седло Центр

2 Седло - - Центр

3 Центр - - Седло

4 Центр - Седло Центр

5 Седло Центр - Седло

c >0

1А, 1В Центр Центр Седло Седло

2 Центр - - Седло

3 Седло - - Центр

4 Седло Центр - Седло

5 Центр - Седло Центр

Рис. 1. Номограмма и таблица для определения числа и типа особых точек фазового портрета

весия являются устойчивыми («центр»), а какие неустойчивыми («седло»). В областях 2 и 3 имеются две особые точки на фазовом портрете системы: 0 и р - этот случай аналогичен случаю физического маятника. В областях 4 и 5 имеются три особые точки на отрезке [0, р] - случай аналогичен случаю бигармонической зависимости моментной характеристики. В областях 1А и 1В имеются четыре особые точки: две точки типа «центр» и две типа «седло».

В связи с изменением коэффициентов

а(г), Ь(г), с(г) в процессе движения происходит эволюция фазовых траекторий, в результате которой они могут пересекать сепаратрисы, попадая в различные области фазового портрета, что сопровождается качественным изменением характера движения.

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

max

I = I ada,

(4)

где а и а - соответственно минималь-

^ min max

ное и максимальное значения амплитуды угла

атаки (при вращении amin =-p и amax = p ),

величина а определяется из (2).

Для системы (1) равенство I = const справедливо для большинства начальных

условий с точностью O (e ln e) на временах

порядка 1/ e [5], где e - малый параметр, характеризующий скорость изменения параметра z. Исключительное множество начальных условий, для которых эта оценка не справедлива, имеют меру O (en), где n > l - любое наперёд заданное число. Режимы движения, соответствующие данным начальным условиям, называются режимами зависания СА в окрестности неустойчивого равновесия. Эти режимы подробно исследованы в [2].

В основу анализа движения системы положены аналитические выражения для интеграла действия (4). Моменты времени,

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

Величина угла атаки на границе перехода от одного типа движения к другому в общем случае зависит от начальных условий углового движения (от законов распределения начальных углов атаки и угловых скоростей), а также от скорости изменения коэффициентов а(z), Ь(z), с(2). Полагается, что

СА за время движения от / = 0 до границы перехода совершено несколько оборотов или колебаний.

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

внутренние области движения А1, А2 от внешней А3. Для выбора области продолжения движения А1 или А2 используется понятие вероятности Р1 или Р2 захвата в каждую из них. В соответствии с [6] эта вероятность определяется как доля фазового объёма малой окрестности начальной точки движения, «захватываемая» в рассматриваемую область в пределе, когда малый параметр е ® 0 и размер окрестности § ® 0 , е << 8 , причём Р + Р2 = 1. Отношение вероятностей вычисляется по формулам:

0 = jd[H (a’a’z)-H (°’a4,z)] Ф df

(5)

dz

i = 1,2,

(6)

где H = a— (a - с )cos a - b cos2 —-—- cos3a -2 V 7 3

гамильтониан; Фz = z = ßz; a = 0, a = a4 -

a

координаты седловой особой точки на фазовом портрете. Интегралы (6) вычисляются

вдоль сепаратрис /1 и /2 , параметризованных временем I невозмущённого движения по ним. Физический смысл функции 0. - ско-рость приближения площади, ометаемой фазовой траекторией, к площади, ограниченной сепаратрисой.

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

2 = Ь^ > 0, то будут иметь место только однократные проходы фазовой точки через сепаратрисы из внешних областей во внутренние.

Рассмотрим СА, для которых положение а1 = 0 является устойчивым. Поэтому в дальнейшем будем исследовать только те области, в которых а1 = 0 является центром. Это области 3 и 4 при с < 0 и области 1, 2, 5 при с >0 (рис. 1).

Исследуем области 1А, 1В. Они отличаются по количеству особых точек от случаев, когда на СА действует синусоидальный или бигармонический восстанавливающий момент. В области 1А внешняя сепаратриса

проходит через а2 = р, в области 1В - че-

рез а4 = агссоБ-

1А и 1В

а

с

-Ъ + т 4с

3Ъ2 Ъ_

16с2 2с

. На границе областей

эти сепаратрисы

сливаются. На фазовом портрете, соответствующем области 1А (рис. 2а), имеются пять областей движения: вращательная, внешняя колебательная, в которой содержатся три внутренних колебательных области. В области 1В (рис. 2б) также имеются пять областей движения: вращательная, две внешних колебательных, и в одной из внешних содержатся две внутренних. На границе областей на фазовом портрете имеются четыре области движения: вращательная и три колебательных.

Вычислим необходимые интегралы: 11 , 12 - интегралы действия, взятые соответственно вдоль сепаратрис /1 , /2 , проходящих через седловую точку а = а4; 13 - интеграл

действия, взятый вдоль сепаратрисы /3 , проходящей через седловую точку а = Р .

Введём дополнительные коэффициенты:

/ = л/ 16с2 + 48ас ,

-3-2-1 0 ¿к, а, а, ж - -3 -2 -1 0 ал а, а, ж

а а, рад О а, рад

Рис. 2. Фазовые портреты: а) область 1А; б) область 1В

.------------------- «4

dA,B = V• 9b2 • 24bc±48ac . (7) ^ = J «da = 2q + 2q2 (F( j,kj-K(k1)) +

-«4

+2q3 (E (j ki) - E (ki)),

Интегралы действия для области 1А. имеют вид:

(11)

4

I1 = J ada = 2q1 + 2q2 (F(j1, k1) - K(k1)) +

-a4

+2q3 (E (j k1) - E (k1)),

, ^4

I2 = J a da = 2q1 + 2q2 F (j k1) + 2q3 E(j1, k1),

a4

(8) (12)

(4c + b + 2да )(4c - b + да)

где j = arcsin,1

1

I2 = J ada = q1 + q2F(j, k1) + q3E(j, k1), (9) V 24дас

k = ' 8c

I 24mc

где j = arcsin

(4c + b + 2да )(4c - b + да)

i---------- q1 = /—

, 14c + b + 2да -b - 2да 9c^/да

k =J--------------, cos a* =---------, -----------

y 8c 4c V b + 2да + 4c (b + 2да - 4c)

q2 = 37ec •

= (4c+3b)Jlm{4c-m+b) q (4b-да)Vb + 2да + 4c ;

1 6c^ (4c - b + да) 3^c

16c2 + 3b2 - 6да2 - 16bc + 4даc + 3bда

q2" 1W5c • ,______

q 8b-2да I3 = J ada = ( 4c - 3b

3 3y/3c ’ a. 6c

V4c + 3b + d (4c + 3b)

+----------j=-----------1( E k) - E j k3)) +

6c V 3

r J a da 8 (7+8c - 6bW9c - 6b + 3a . -dj(4c + 3b + d) . .

13 =-Jada = 3J3 (f + 16c) + + V 6 я---------1 (K(k3)-F(j„k3)),

1 V 4c + b + 2да ’

(да + 2b )^ 2 (4c - да + b )(да + 4c - b)

a

+

Jf (3b + 4c) c

~~?= Ek3) + при a >—,

3V3c (10) 3

I7 (f - 4c -3b) F (j k)

- ^3c F(j3> ^ (13)

14c + 3b + d

где j3 = 2arctg

16c k= 7+3b + 4c где j = arcsl% 16c

_ 1T, , 14c + 3b - d -3b + 4c + d

Интегралы действия для области 1В k3 = J----------------—----;, cos a* =-------------

имеют вид:

4c + 3b + d’ 8c

a

a

p л/(12с - 3b - d) (4c - 9b +

I3 = I ada = ^----------. M

a* 3pc (12c - 3b + d)

V2d (3b + 4c). .

+------; r- \ E (k3) - E (j3, *3> )-

6cv 3

■V2d (3b + 4c + d)

d)

+

12cV3

c

(K (kä) - F (p„ ks)),

при a <-

3

(14)

где p3 = arcsin.,

2d

d - 3b + 12c

, id - 3b - 4c -3b + 4c - d

k3 = J----------------, cos a* =■

2d

8c

p (4c - 3b) 2

I ada =-¡=—, при a -

L 3c46

(15)

где соб а* = ■

-3b

4c

tA = Ain Ll tA = —ln-^ ß V 2 ß L + 212

= -l^ l0

ß

I1 + 12

= -ln l0 11

ß 2I3

(16)

где интегралы 11 , 12 , 13 вычисляются по формулам (8) - (15) при

a = ao, b = bo,

значение интеграла действия 10 определяется по начальным условиям и в общем случае вычисляется по формуле (4). В случае, когда

коэффициенты а0, Ь0, с0 малы по сравнению с угловой скоростью в начальный момент времени, интеграл действия 10 можно вычислить по формуле: 10 = 2ра0.

Пусть начальные условия соответствуют движению во внешней колебательной области. Тогда формула для вычисления времени перехода *2 из внешней колебательной области во внутреннюю имеет вид:

tA

I")

2

ln

I0

ß L + 212

-ln^°-ß 2I3

Пусть начальные условия соответствуют вращательному движению СА. Из (8) -(15), учитывая постоянство интеграла действия, можно определить значение коэффициента с в момент перехода из одной области движения в другую через начальные условия. Тогда, учитывая, что

с = с0ехР (р -*0))

и полагая ^ = 0, время *1 перехода из вращения в колебание и время *2 перехода из

внешней колебательной области во внутреннюю определяются по формулам (индексами А и В вверху обозначены области 1А и 1В соответственно):

I

Вероятности попадания в колебания относительно различных положений равновесия ( а = 0, а = ±а3, а = р) с учётом (5), (6) вычисляются следующим образом:

PA =

0

I

L1 + 212

PA =-

4-i-v

I

L1 + 212

PB =-

Г,

I1

I1 +12

PpB

I1 +12

PB.

(17)

Исследуем области 2 и 5 (рис. 3). На фазовом портрете, соответствующем области 2 (рис. 3 а), имеются две области движения: вращательная и колебательная. На фазовом портрете, соответствующем области 5 (рис. 3б), имеются три области движения: вращательная и две колебательные.

В области 2 имеется только одна сепаратриса, проходящая через а = р. Выражение для интеграла действия совпадает с выражением (10) для 13 в области 1А. Задача

B

B

t1

t

B

2

2

2

— о

С

-1

-2

2

á, 1 г^о

с

-1

-2

а

б

-3

1 «4 2

л

а, рад

Рис. 3. Фазовые портреты: а) область 2; б) область 5

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

21 7о

t = — ln-^

b V

(18)

Pp =■

12

I + 12

(19)

L =

я

I a da ■-

W3c - 2b-

В области 5 имеется одна сепаратриса, проходящая через а = а4. Выражения для интегралов действия совпадают с выражениями (11), (12) для 1 и 12 в области 1В. Время перехода из вращения в колебание и вероятности попадания в колебания относительно а = 0 или а = р можно определить по формулам

г =—1п—^, Р0 = —^,

Р 11 +1/ 0 11 +12 ’

■4ld (3b + 4c) 3*J3c ■\Í2d (4c + 3b + d) 6л/3с

(E (£3) - E j £3) )-

(F (£3) - F (^3, £3)),

(20)

где j3 = arcsin „

12c - 3b + d -4c - 3b + d

£3 =

Id - 3b - 4c

2d

Исследуем области 3 и 4 (рис. 4). На фазовом портрете, соответствующем области 3 (рис. 4а), имеются две области движения: вращательная и колебательная. На фазовом портрете, соответствующем области 4 (рис. 4б), имеются три области движения: вращательная и две колебательные.

В области 3 имеется одна сепаратриса, проходящая через а = р. Выражение для интеграла действия имеет следующий вид:

Задача выбора области продолжения движения не возникает, а время перехода из вращения в колебание можно определить аналогично области 2 по формуле (18).

В области 4 имеется одна сепаратриса,

проходящая через a = a4. Выражения для интегралов действия I1, I2, взятых вдоль сепаратрис l1 и l2, вычисляются по следующим формулам:

а4

I1 = I ada = 2q1 + 2q2 (F(j1, k1) - K(k1)) +

-a4

+2q3 (E(j ki) - E(ki)),

а а, рад V а, рад

Рис. 4. Фазовые портреты: а) область 3; б) область 4

I2 = J ada = 2q1 + 2q2F(j1,k1) + 2q3E(j1,k^,

a4

(22)

где j1 = arcsin.,

\m - b + 4c

8c

k = 8c

1 V 4c - b - 2m ’

qj

q2 =

qs =

О ¥ <N - m + b) (4с - b + m)

1 о 40 6c -5b2 - 4bm)

6су1 (4b- Ь (ь- m и + 2m - 4c)

^2 (b + 2m - 4c)

6^c

Время перехода из вращения в колебание и вероятности попадания в колебания относительно а = 0 или а = ж можно определить аналогично области 5 по (19).

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

а0 = 0,05 с

b0 = - 0,1 с-

Со = 0,1 с"

а0 = 1,5 рад, a&0 = 0,8 рад / с, ß = 0,03 с j

СА, совершив три полных оборота ( бж), в некоторый момент времени начинает совершать колебательные движения относительно устойчивого положения равновесия а = 0, охватывая при этом два неустойчивых положения равновесия

±а4 =±0,819 и два устойчивых положения равновесия ±а3 =±1,755. Затем СА в следующий момент времени 12 начинает совершать колебания относительно устойчивого положения равновесия -а3 =-1,755 . Аэродинамические характеристики СА относят-

а Ь

ся к области 1А ( —- = 0,5 — = -1, рис. 1),

С0 С0

фазовый портрет аналогичен приведённому на рис. 2а. Следовательно, вычисления интегралов действия необходимо производить по формулам (8) - (10), а времён перехода и вероятности попадания в колебания относительно различных положений равновесия

( а = 0, а = ±а3) - по формулам (16), (17) с индексом А. В результате расчётов:

^ = 26,576 с, ^ = 118,661 с, Р0 = 0,282,

Ра = 0,359, Ра = 0,359.

(а3 ^ Сл3 ^

Как видно из рис. 5, результаты численного интегрирования уравнения (1) для вре-

'О 'О

Рис. 5. Характер изменения углового движения СА а) фазовая траектория; б) зависимость угла атаки от времени

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

В [2] рассмотрено плоское движение неуправляемого тела со следующей момент-ной характеристикой:

m(a) = 0,694sin а + 0,342sin 2а - 0,126sin3a.

Коэффициенты, стоящие при гармониках разложения, соответствуют случаю области 4 (— = -5,508 — = -2,714, рис. 1), С0 С0

фазовый портрет аналогичен приведённому на рис. 4б. Вычисляется вероятность попадания в колебания относительно устойчивого положения равновесия а = ж путём расчёта по интегральным формулам, а также многократным численным интегрированием дифференциального уравнения (1). Получено значение вероятности Рж = 0,05. Расчёт, проведённый по формулам (19), (21), (22), дал такой же результат.

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

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

1. Кузмак Г. Е. Динамика неуправляемого движения летательных аппаратов при входе в атмосферу [Текст]/Г. Е. Кузмак - М.: Наука, 1970.

2. Ярошевский В. А. Движение неуправляемого тела в атмосфере [Текст]/ В. А Ярошевский - М.: Машиностроение, 1978.

3. Асланов В. С., Тимбай И. А. Переходные режимы углового движения КА на верхнем участке траектории спуска [Текст]/ В. С. Асланов, И. А. Тимбай// Космические исследования. - 2001. - Т. 35. - № 3. - С. 279 -286.

4. Асланов В. С. Пространственное движение тела при спуске в атмосфере. [Текст]/ В. С. Асланов - М.: ФИЗМАТЛИТ, 2004.

5. Нейштадт А. И. Об изменении адиабатического инварианта при переходе через сепаратрису [Текст]/А. И. Нейштадт//Физи-ка плазмы. - 1986. - Т. 12. - Вып. 8. - С. 992.

6. Лифщиц И. М., Слуцкин А. А., На-бутовский В. М. Об особенностях движения заряженных частиц в переменном и неоднородном электромагнитном поле [Текст]/ И. М. Лифшиц, А. А. Слуцкин, В. М. Набу-товский// Журн. эксперим. и теорет. физики.

- 1961. - Т. 41. - Вып. 3. - С. 939.

References

1. Kuzmak, G. Ye. Dynamics of spacecraft uncontrolled motion during reentry into the atmosphere / G. Ye. Kuzmak - Moscow: Nauka, 1970.

2. Yaroshevsky, V. A. Motion of uncontrolled bodies in the atmosphere / V. A. Yaroshevsky - Moscow: Machinostro-yeniye, 1978.

3. Aslanov, V. S., Timbay I. A., Transient modes of spacecraft angular motion at the upper section of the reentry trajectory / V. S. Aslanov, I. A. Timbay // Kosmitcheskiye Issledovanya (Space Research) 2001. Vol. 35 No. 3, pp. 279 -286.

4. Aslanov, V. S. Spatial motion of a body during its reentry into the atmosphere / V. S. Aslanov - Moscow: Physmatlit, 2004.

5. Neuschtadt, A. I. Variation of Adiabatic Invariant during the crossing of a separatrix / A. I. Neuschtadt // Physica Plasmy (Plasma Physics), 1986. Vol.12. Issue 8. p. 992.

6. Lifshits, I. M., Slutskin A. A., Nabutovsky V. M. Peculiarities of charged particle motion in an alternative inhomogene-ous electromagnetic field / I. M. Lifshits, A. A. Slutskin, V. M. Nabutovsky // Journal of experimental and theoretical physics. 1961. Vol. 41, Issue 3. p. 939.

ANALYSIS OF PLANE MOTION RELATIVE TO THE CENTRE OF MASS OF A DESCENT CAPSULE WITH THE THREE-HARMONIC MOMENT CHARACTERISTICS DURING ITS REENTRY

© 2010 Ye. V. Barinova, I. A. Timbay

Samara State Aerospace University named after academician S. P. Korolyov (National Research University)

The paper analyses plane motion of an uncontrolled descent capsule about its centre of mass. The aerodynamic stabilizing moment of the spacecraft by an odd Fourier series by the angle of attack with the first three harmonics. A nomogram for determining the type and number of singular points of the system phase portrait has been constructed depending on the relation of coefficients at decomposition harmonics. Analytical formulae for the action integral taken along the separatrix, expressed through elementary functions and elliptical integrals of the first and second kind have been established. The moments of transition between various areas of the phase plane have been determined. Formulae for assessing the possibility of capture into one region or another have been established for the cases of motion when the phase point may fall into various oscillation regions when crossing the separatrix.

Descent capsule, aerodynamic stabilizing moment, phase portrait, separatrix, action integral, probability.

Информация об авторах Баринова Елена Витальевна, аспират кафедры высшей математики Самарского государственного аэрокосмического университета имени академика С. П. Королёва (национального исследовательского университета). Область научных интересов: динамика летательных аппаратов. E-mail: l5545@yandex.ru.

Тимбай Иван Александрович, доктор технических наук, профессор, заведующий кафедрой высшей математики Самарского государственного аэрокосмического университета имени академика С. П. Королёва (национального исследовательского университета). Область научных интересов: динамика летательных аппаратов. E-mail: timbai@ssau.ru.

Barinova Yelena Vitalyevna, post-graduate student of the department of higher mathematics, Samara State Aerospace University named after academician S. P. Korolyov (National Research University), l5545@yandex.ru. Area of research: aircraft dynamics.

Timbay Ivan Alexandrovitch, doctor of technical sciences, professor, head of the department of higher mathematics, Samara State Aerospace University named after academician S. P Korolyov (National Research University), timbai@ssau.ru. Area of research: aircraft dynamics.