CC BY
13
www.sibac.info
УПРАВЛЕНИЕ И МОДЕЛИРОВАНИЕ СЛУЧАЙНОГО ДВИЖЕНИЯ КОСМИЧЕСКОГО АППАРАТА В ОКРЕСТНОСТИ РЕЗОНАНСА ПРИ СПУСКЕ В АТМОСФЕРЕ ВЕНЕРЫ
Любимов Владислав Васильевич
д-р техн. наук, доцент, заведующий кафедрой высшей математики, Самарский Государственный Аэрокосмический Университет,
РФ, г. Самара E-mail: vlubimov@mail.ru
Лашин Вячеслав Сергеевич
аспирант,
Самарский Государственный Аэрокосмический Университет,
РФ, г. Самара E-mail: glory665@mail.ru
CONTROL AND SIMULATION OF A RANDOM MOTION OF THE SPACECRAFT IN THE VICINITY OF THE RESONANCE AT THE DESCENT IN VENUS ATMOSPHERE
Vladislav Lyubimov
associate professor, doctor of technical sciences, Head of the Chair of Higher Mathematics, Samara State Aerospace University,
Russia, Samara
Vyacheslav Lashin
postgraduate student, Samara State Aerospace University,
Russia, Samara
АННОТАЦИЯ
Рассматривается спуск космического аппарата с малой инерционной и аэродинамической асимметрией в атмосфере Венеры. Исследуется случайный характер вращательного движения космического аппарата в окрестности главного резонанса при малых углах атаки. Показывается, что применение двух разных законов управления величиной инерционной асимметрии обеспечивает гарантированный
проход через резонанс. Приводится условие гарантированного прохода через резонанс.
ABSTRACT
We consider the descent of the spacecraft with a small inertial and aerodynamic asymmetry in the atmosphere of Venus. We study the random nature of the rotational motion of the spacecraft in the vicinity of the main resonance at small angles of attack. Recorded the condition of the guaranteed passage through resonance. It is shown that the application of the two different laws control the amount of inertial asymmetry provides a guaranteed passage through resonance.
Ключевые слова: твердое тело; космический аппарат; атмосфера; резонанс; вероятность захвата; асимметрия; угловая скорость; угол атаки.
Keywords: solid; spacecraft; atmosphere; resonance; probability of capture; asymmetry; angular velocity; angle of attack.
При отделении от базового космического аппарата спускаемому аппарату передается небольшое значение кинетического момента [1]. Это приводит к появлению малой начальной угловой скорости, которую имеет спускаемый аппарат при входе в плотные слои атмосферы. Известно, что наличие асимметрии у спускаемого аппарата может привести к реализации различных резонансных явлений, способных явиться причиной различных аварийных ситуаций. Кроме непосредственно длительных резонансов моменты от асимметрии могут быть причиной аварийных ситуаций, вызванных реализацией вторичных резонансных эффектов [2-4].
Начальные условия движения на сепаратрисе, отделяющей область вращения от области резонансных колебаний, расположены случайным образом. В связи с этим возникает задача о нахождении вероятности захвата системы в резонансную область [5].
Рассмотрим движение космического аппарата в атмосфере Венеры как твердого тела конической формы с малой инерционной и аэродинамической асимметриями при малых углах атаки. Квазилинейная система уравнений движения асимметричного космического аппарата относительно центра масс, полученная посредством известного подхода [6], имеет вид:
, А 2 2
dœ m œ 9а
^ = Е-J^sin(2e + 294), (1)
dt Ix
¿/а 2а
ш 21л Г
= +е-:—[2сот -со12рт(29 + 294) +
Л 2со
А л ч иа ^ш
ш
8-008(9 + 0!) -8
2ша 2ша йг
А
= шх-ш12 ±8^— 81и(9 + 91)± йг х 1,2 2аш а
т&ш1 2 8 ,
(2)
2ш
[2шх -ш12^ 008(29 + 294)
а
dш _ ш dq
аг 2q dt
(3)
(4)
Здесь 9=ф-я/2, ш - угловая скорость тела относительно
связанной с телом оси X, а - угол атаки, ф-угол аэродинамического крена, ш-известная функция переменных движения центра масс,
тА = ^(т{1 )+(тА ^ тА=^ (1у2 )2 + (А )2
2 2 А Ш ф — 2 А Ш ф — 2
тА =--тф0 - 1хгш2х т2 =--т|0 +1 шХ
т21 , ;
со52'4=-тА 8т2'4=-тА
тф0, тф, Сж1, тг1 - коэффициенты аэродинамической асимметрии
формы, лобового сопротивления и восстанавливающего момента.
При малых величинах момента a и малых углах атаки можно получить выражение для оценки вероятности захвата [7]:
2\Ф\
Рг =
2Аа + ф . (5)
42аа4Гх (2 - 1х)
где
„ 8(1 - 1х )ш2 Ра-А Здесь Ф = —)=- —-V тхт 0080 - 02) +
" " " ~х?П
а
+ 8
/2(1 - 1Х)(
(2 - 1х)
1-
1 - 1х Л
(2 - 1Х )2J 4ГХ
ж
а =
2(1 - 1х )1/2
2 - 1х Ж
Выражение (5) рассматривается при выполнении неравенства \ф\ < 2ж\а\
В итоге получаем, что величина параметра асимметрии тл при вероятности захвата равной нулю (гарантированный проход через резонанс) определяется следующим образом:
Л 21Х (1 - 1х Г2 И®/
т =-*-—---
2 - 31х
2 2 ш а
(6)
Здесь величины а,ш ,/ И вычисляются на сепаратрисе [7]. С целью уменьшения величины вероятности захвата в резонанс введем в рассмотрение экспоненциальный и синусоидальный законы уменьшения инерционной асимметрии:
тл = кехН, тл = к ът.(ЛН)
(7)
(8)
Здесь к и X - коэффициенты управления, Н- высота полета аппарата.
При атмосферном спуске наблюдается уменьшение высоты полёта Н. В результате осуществляется уменьшение параметра инерционной асимметрии тл
При численном моделировании предполагалось, что космический аппарат конической формы с массой 10 кг и малой асимметрией совершает спуск в атмосфере Венеры при малых углах атаки. Значения коэффициентов управления были следующими: к=0.029, X = 10-5 и к=0.3, X=10-6 (для экспоненциального и синусоидального законов управления, соответственно).
Часть результатов численного моделирования для экспоненциального закона представлена на рисунках 1-3. На рисунке 1 содержатся: изменение угловой скорости ш (1) и её резонансных
с
СибАК
Инновации в науке
www.sibac.info_№ 11 (48). 2015 г
значений шх (1) (кривая, пропорциональная изменению скоростного напора д(ф. Параметры управления выбраны таким образом, что
показанное на рис. 1 изменение угловой скорости шх (1)
г
и резонансного значения данной угловой скорости шх (1) характерно как для экспоненциального, так и для синусоидального закона изменения . На рисунке 1 показано, что при уменьшении инерционной асимметрии (в соответствии с законом (7)) на восходящей ветви скоростного напора происходит реализация прохода через резонанс. Проход через резонанс имеет место на 250 секунде полета в атмосфере. Дальнейшее уменьшение параметра асимметрии приводит к стабилизации угловой скорости и угла атаки. На рисунке 2 представлен результат изменения инерционной асимметрии в соответствии с экспоненциальным законом (7). На рисунке 3 показано, что уменьшение инерционной асиметрии способствует также уменьшению вероятности захвата. При достижении резонансных значений угловой скорости (250 с.) вероятность захвата равна малой величине (менее 0,15). На рисунках 4-5 представлены: уменьшение инерционной асимметрии и расчет вероятности захвата, соответственно. При этом, в качестве закона управления асимметрией использовался синусоидальный закон (8). Численное моделирование случайного движения космического аппарата в окрестности резонанса позволяет сделать вывод о том, что выбор коэффициентов в законах управления (7) и (8) позволяет обеспечить гарантированный проход через резонанс. Из сравнения результатов численного моделирования, представленных на рисунках 2 и 4 следует, что синусоидальный закон управления (по сравнению с экспоненциальным законом) позволяет достигнуть меньших величин инерционной асимметрии. Величины вероятности захвата в резонанс при достижении резонансных значений угловой скорости на рисунках 3 и 5 также несколько меньше при синусоидальном законе, чем при экспоненциальном законе. Следовательно, выбирая один из законов управления величиной инерционной асимметрии (7) или (8) можно уменьшить вероятность захвата системы в резонанс до величины, обеспечивающей гарантированный проход системы уравнений движения через резонанс (6).
Рисунок 1. Эволюция угловой скорости и резонансного
значения
Рисунок 2. изменение параметра по экспоненциальному закону
РГ 0.111
IV! 5
о.о;
\
ш ПГРППЧ 11^ 1 1 < д .
1 с
25(1 50(1 7$(! 1000
Рисунок 3. Вероятность захвата твердого тела в резонанс при экспоненциальном законе изменения
250 500 750 1 (9СН*
Рисунок 4. изменение параметра по синусоидальному закону
Рт
0,30
015
0.20
0 г50 500 750 ] ООО
Рисунок 5. Вероятность захвата твердого тела в резонанс при синусоидальном законе изменения
Результаты работы представляют определенный практический интерес, так как в рассмотренной системе малое уменьшение величины инерционной асимметрии позволяет обеспечить реализацию нерезонансного вращения космического аппарата. Представленная работа относиться к классу задач, в которых изучается управление в системах с динамическим хаосом [8].
Список литературы:
1. Калаев М.П., Любимов В.В., Сёмкин Н.Д. Полунатурное и имитационное моделирование процесса отделения микроспутника // Гироскопия и навигация. - 2014. - № 2 (85). - С. 52-60.
2. Любимов В.В. Асимптотический анализ вторичных резонансных эффектов при вращении космического аппарата с малой асимметрией в атмосфере // Известия вузов. Авиационная техника. - 2014. - № 3. -С. 23-28.
3. Любимов В.В., Лашин В.С., Численное моделирование резонансного эффекта с изменением направления вращения при атмосферном спуске асимметричного наноспутника // Информационные технологии и нанотехнологии (ИТНТ-2015): материалы международной конференции и молодежной школы. - Самара, 2015. - С. 250-253.
4. Заболотнов Ю.М., Любимов В.В. Вторичный резонансный эффект при движении КА в атмосфере // Космические исследования. - 1998. - Т. 36. -№ 2. - С. 206-214.
www.sibac.info
5. Нейштадт А.И. Захват в резонанс и рассеяние на резонансах в двухчастотных системах // Тр. МИАН. - 2005. - Т. 250. - C. 198-218.
6. Заболотнов Ю.М. Асимптотический анализ квазилинейных уравнений движений в атмосфере КА с малой асимметрией III. Космические исследования// - 1994. - Т. 32. - № 4-5. - C. 112-125.
7. Любимов В.В. Оценка вероятности захвата в резонанс при движении динамически несимметричного твердого тела // Вестник СамГТУ. Серия: Физико-математические науки. - 2007. - № 2. (15) - С. 110-115.
8. Магницкий Н.А., Сидоров С.В. Новые методы хаотической динамики. -М.: Едиториал УРСС, 2004. - 320 с.
