CC BY
55
УДК 629.783
АНАЛИЗ ДИНАМИЧЕСКИХ РЕЖИМОВ ПРОСТРАНСТВЕННОГО ДВИЖЕНИЯ КОСМИЧЕСКОГО АППАРАТА С ТВЕРДОТОПЛИВНЫМ РАКЕТНЫМ
ДВИГАТЕЛЕМ
© 2012 В. Л. Балакин, А. В. Дорошин, М. М. Крикунов
Самарский государственный аэрокосмический университет имена академика С.П. Королёва (национальный исследовательский университет)
Рассматривается движение космического аппарата (КА) переменного состава вокруг центра масс. Проводится анализ движения и синтез условий реализации требуемых режимов нутационно-прецессионного движения КА переменного состава с помощью качественного метода, основанного на вычислении кривизны фазовой траектории. Определяются возможные эволюции КА и причины этих эволюций. Проводится численное моделирование движения КА для нелинейных случаев изменения инерционно-массовых параметров.
Космический аппарат, динамическая система, фазовое пространство, переменный состав, синтез режимов движения, численное моделирование.
Введение
Исследование динамики космического аппарата (КА) переменного состава является одной из задач механики космического полёта. Классические случаи движения КА переменного состава были рассмотрены ранее [1,2]. Некоторые аспекты динамика движения КА переменного состава исследовались для случаев нелинейного изменения инерционно-массовых параметров при наличии возмущений [3-6].
В ракетно-космической технике находят применение космические аппараты, на которых устанавливаются
твёрдотовливные ракетные двигатели [7]. Заряды для таких двигателей могут иметь сложную форму [8, 9]. От формы и массы заряда, то есть от его инерционно-массовых характеристик, будет зависеть вектортяги КА и динамика его движения вокруг текущего положения центра масс. Для исследования динамики движения КА переменного составас зарядами сложной
формыцелесообразно исследовать динамику КА, имеющего заряд простой формы. В статье рассматривается заряд, имеющий форму цилиндра.
1. Постановка задачи
Рассматривается КА переменного сотава, который удовлетворяет следующим четырём допущениям.
Допущение 1. Масса и геометрия КА переменного состава изменяется так, что при этом он все равно остается динамически симметричным в каждый момент времени t:
A(t) = B(t) Vte[t0;tk]. (1)
Здесь A(t) и B(t) - зависимости моментов инерции с учётом движения центра масс вдоль оси Oz связанной системы координат, относительно его места нахождения в начальный момент времени.
Допущение 2. Движение
рассматривается в отсутствии внешних возмущающих сил и моментов Ме:
M (p, q, r )=м; (p, q, r ) =
= Mez (p, q, r) = 0.
Допущение 3.
Допущение 3. Рассматривается движение с КА переменного состава с гироскопической стабилизацией:
s =
<<1.
(3)
Допущение 4. Считается, что центр масс корпуса КА в начальный момент времени совпадает с центром масс твёрдотопливного заряда.
2. Математическая модель
Для исследования динамики движения КА переменного состава построим математическую модель. Рассмотрим
прецессионное движение космического аппарата переменного состава. Уравнения невозмущённого движения записываются в виде [4]:
A (t) p + ( C (t)-A (t)) qr = 0,
A (t) q -(C (t)- A (t)) pr = 0, (4)
C (t) r = 0.
Здесь C(t) - момент инерции КА относительно оси Oz связанной системы координат.
Сделаем замену переменных [6]:
p = G(t) sm F(t), q = G(t) cos F(t). (5)
Продифференцировав (5) по времени и подставив в (4), получим уравнения движения в новых переменных:
F =
1 с(0
л (t).
r, r = 0, G = 0.
(6)
У = p srnp + q cosp,
W =
1
p = r
cosy sin y
cosy
( p cosp- q sin p), ( p cosp- q sin p).
(7)
интерпритировать как годограф вектора тяги ракетного двигателя твёрдого топлива (РДТТ). Следовательно, для повышения точности импульса, создаваемого РДТТ КА, необходтмо добиться, чтобы продольная ось КА монотонно приближалась к оси процессии и, следовательно, годограф апекса продольной оси (и вектора тяги) КА переменного состава, как фазовой траектории в плоскости (у, у), представлял собой скручивающуюся спираль. Для обеспечения такого движения необходимо монотонное увеличение квадрата (модуля) кривизны этой фазовой траектории (годографа).
Запишем выражение для квадрата кривизны параметризованной кривой [6]:
k2 =
_ (w -w )2
Далее будем рассматривать движение КА переменного состава в углах Крылова [6]:
(у2 W (8)
)3
С учетом перепишется в виде:
к2 = Ф / G2.
Найдем производную функции кривизны, заданной неявно:
(9)
выражение (9) (10)
кк = ФФй-дФ2
(11)
Допущение 3 позволяет считать углы y и w малыми (рис. 1).
Для того чтобы модуль кривизны был монотонно возрастающей функцией, необходимо, чтобы функция Р = кк была больше нуля:
к ^ кк > 0 ^ ... .. 2 (11) ^ФФО - ОФ2 = Р (t )> 0.
С учётом (6) и (8), условие (12) перепишется в виде:
Р (t) = (СА - АС)> 0. (12) А
Так как
CGn
>
0 Vt e[í0; tk ] , то
Рис. 1. Углы Крылова Принимая во внимание (5) и допущение 3, запишем уравнения (7) в виде: у = G cos Ф( t), W = G sin Ф( t) ,p = r, (8) где Ф0 ) = F (t )-p(t). Фазовый портрет системы (8) в плоскости углов Крылова (у, ^) можно
необходимо добиваться выполнения условия:
(14)
)= C (t )A(t)- A(t )C (t) =
= -1С ]> 0.
сЧ \ А )
Из этого можно сделать следующий важный вывод: годограф вектора тяги РДТТ с цилиндрическим пороховым зарядом может изменить свой вид со
3
скручивающейся спирали на
раскручивающуюся (и наоборот) тогода и только тогда, когда отношение моментов
инерции С достигает экстремального
значения.
Задача синтеза динамики движения КА переменного состава сводится к отысканию таких параметров инерционно-массовых характеристик заряда РДТТ, при которых функцией P(t) была бы положительно определённой на
рассматриваемом промежутке времени [7-9].
3. Анализ динамики движения
Применим построенную модель для анализа движения КА переменного состава при горении заряда цилиндрической формы. Моменты инерции относительно точки, в которой находился центр масс в начальный момент времени, записываются следующим образом:
R2 H2
mR2
mzc =
A (t ) = 4 + 4ЧЖ + = A0 +aH + ßH3, где a = nR4p/4, ß = яЯ2р/3 ;
C (t ) = Co + [H; где л = nR 4р/ 2;
A (t ) = aH + 3ßH2 Hi; C (t ) = [Hi.
Подставим (15)-(18) в (14):
(15)
(16)
P (t ) =
= ^ ( Ca«" H 2ß( 2ßH + 3C0)).
Если моменты инерции корпуса не учитывать и рассматривать только горение заряда (A0 = 0, C0 = 0), то функция
P (t)
будет иметь вид:
2ßyH2 Hl
P (t ) = -
a + ßH7
(20)
Пусть цилиндрический заряд поджигается с основания, и его высота изменяется по линейному закону:
H = H (t) = H0 - Ät,
(21)
Л = ^ = А + т| — +-|, с = С0 +
^ I 4 12 ) 0 2 '
где А0, С0 - моменты инерции корпуса К; т = pV = ж Я2Нр. - масса заряда в текущий момент времени; Я -радиус заряда; Н -высота заряда в текущий момент времени ; р - плотность материала, из которого изготовлен заряд.
В силу изменения геометрии заряда центр масс будет перемещаться вдоль оси Оz связанной системы координат:
Zс = 2 (Н 0 - Н (г )) = 1 Л/.
Выпишем в явном виде моменты инерции и их производные:
(17)
(18)
где Н0 - высота цилиндра в начальный
момент времени; Ы - скорость выгорания заряда.
Рассмотрим два предельных случая:
1) моменты инерции корпуса пренебрежимо малы по сравнению с моментами инерции заряда (слу-чай когда рассматривается только горение заряда);
2) моменты инерции заряда пренебрежимо малы по сравнению с моментами инерции корпуса КА (случай КА постоянного состава).
Рассмотрим сначала функцию Р(г) без учёта моментов инерции корпуса КА. С учётом (21) выражение (20) примет следующий вид:
р (/ )= 2ВДЧН „- Ы)3 (22)
а + р(Н0 - Ыг)
Так как на промежутке го-рения заряда, то из (22) следует, что при линейном сгорании твёрдотопливного за-ряда его параметры не могут повлиять на эволюции движения годографа вектора тяги, который будет представлять из себя скручивающуюся спираль.
В качестве примера на рис. 2 представлены графики зависимостей моментов инерции от времени для значений параметров: А0 = 8 кг • м2 ; С0 = 8 кг • м2 ; Я = 0,3 м; Н0 = 0,5 м; Ы = 0,01 м / с.
п
я
Й S 10 15 20
£ с
0 5 10 15 21 ^ с
Рис. 2. Моменты инерции для заряда цилиндрической формы
На рис. 3 представлен график функции
Р(t) для рассматриваемого случая.
О 5 101 15 2й 1 г
Рис. 3. Функция Р (t) для заряда цилиндрической формы
Функция Р (t) в рассматриваемом случае является положительно определённой на всём рассматриваемом интервале (от 0 до 20 с).
На рис. 4 представлен годограф вектора тяги РДТТ для рассматриваемого случая.
у/,рад! с
Рис. 4. Годограф вектора тяги РДТТ с зарядом цилиндрической формы в плоскости углов Крылова у - у , у0 = 0,05 рад , у0 = 0,05 рад
Из следует, что функция P(t) может иметь нули на рассматриваемом промежутке времени:
при
H = const . следует, что P(t) = 0.
у/, рад! с
-0.4 -0.2 0:2 0 4
0 1',рад!с
-0.2 \
-0* ]
-ОБ 1
-0.8
А = 8 кг • м2
С„ = 8 кг • м
Л = 0,3 м
Н0 = 0,5 м; Л = 0,01 м /с.
Рис. 7. Функция Р (t) для заряда цилиндрической формы
Функция Р (t) в рассматриваемом случае является положительно определённой на всём рассматриваемом интервале (от 0 до 20 с).
На рис. 8 представлен годограф вектора тяги РДТТ для рассматриваемого случая.
Рис. 5. Годограф вектора тяги РДТТ с зарядом цилиндрической формы в плоскости углов Крылова у - у , у 0 = 0,05 рад , у0 = 0,05 рад
Вывод: в отсутствии заряда эволюций движения КА происходить не может. Таким образом, эволюции движения КА могут происходить только в тех случаях когда инерционно-массовые характеристики КА и заряда одного порядка.
Приведём два примера, когда характер движения может поменяться и когда он остаётся постоянным.
Пример 1. На рис. 6 представлены графики зависимостей моментов инерции от времени для значений параметров:
Рис. 8. Годограф вектора тяги РДТТ с зарядом цилиндрической формы в плоскости углов Крылова у - у , у0 = 0,05 рад , у0 = 0,05 рад Пример 2. На рис. 9 представлены графики зависимостей моментов инерции от времени для значений параметров: А0 = 8 кг• м2 ; С0 = 8 кг• м2 ; Я = 0,3 м ; н0 = 0,5 м; п = 0,01 м / с.
На рис. 10 представлен график функции Р(t) для рассматриваемого случая.
Рис. 6. Моменты инерции для заряда цилиндрической формы
На рис. 7 представлен график функции
Р(t) для рассматриваемого случая.
Рис. 10. Функция Р (t) для заряда цилиндрической формы
Функция р (г) в рассматриваемом случае пересекает ось абсцисс г = 17,5 с.
На рис. 11 представлен годограф вектора тяги РДТТ для рассматриваемого случая.
уг,рад!с
-0r6 -04 -02 04 0,6 0,8
7',pt¿Mc
-0.2
-0.4
-0.6
-0.B
-1
-1.2
14
Рис. 11. Годограф вектора тяги РДТТ с зарядом цилиндрической формы в плоскости углов Крылова у -у , у0 = 0,05 рад , у0 = 0,05 рад
Заключение и выводы
В работе применён качественный метод анализа фазового пространства
неавтономных динамических систем, основанный на определении кривизны фазовой траектории. С помощью этого метода проведён анализ движения и синтез условий реализации требуемых режимов нутационно-прецессионного движения КА переменного состава. Определены возможные эволюции КА и причины этих эволюций.
Проведено численное моделирование движения КА для линейного случая изменения инерционно-массовых
параметров. Построены фазовые портреты для рассматриваемых случаев движения КА с твердотопливным ракетным двигателем.
№ 5 (36), 2012
Работа выполнена при поддержке Российского фонда фундаментальных исследований (РФФИ № 11-08-00794-a).
Библиографический список
1. Охоцимский Д.Е., Сихарулидзе Ю.Г. Основы механики космического полёта. М.: Наука, 1990.
2. Космодемьянский А.А. Курс теоретической механики. Часть 2. М.: Просвещение, 1966. 398 с.
3. Аншаков Г.П., Асланов В.С., Балакин
B.Л., Дорошин А.В. и др. Динамические процессы в ракетно-космических системах // Вестник Самарского государственного аэрокосмического университета. Самара: СГАУ, №1, 2003 г.
4. Дорошин А.В. Эволюции прецессионного движения неуравновешенных гиростатов переменного состава // Прикладная математика и механика, Т. 72. Вып. 3, 2008.
5. В.С. Асланов, А.В. Дорошин Влияние возмущений на угловое движение космического аппарата на активном участке спуска // Космич. исслед. 2008, том 46, №2,
C. 168-173.
6. A.V. Doroshin // International Journal of Non-Linear Mechanics # 45 (2010) p. 193-205.
7. Фахрутдинов, И. Х. Ракетные двигатели твердого топлива [Текст] / И. Х. Фахрутдинов. — М.: Машиностроение, 1981. — 223 с.
8. Виницкий, А. М. Конструкция и отработка РДТТ [Текст] / А. М. Виницкий, В. Т. Волков, И. Г. Волковицкий, С. В. Хо-лодилов; Под ред. А. М. Виницкого. —М., Машиностроение, 1980.—230 с.
9. Петренко, В. И. РДТТ с регулируемым модулем тяги [Текст] / В. И. Петренко, В. Л. Попов, А. М. Русак, В. И. Феофилактов. — Миасс: издательство ГРЦ «КБ им. академика В.П. Макеева». — 1994. — 245 с.
SYNTHESIS OF DYNAMIC MODES OF ATTITUDE MOTION OF SPACECRAFT WITH
SOLID PROPELLANT ROCKET ENGINE
© 2012 V. L. Balakin, A. V. Doroshin, M. M. Krikunov
Samara State Aerospace University 18
named after academician S. P. Korolyov (National Research University)
An attitude motion of a spacecraft (SC) with variable structure around mass-center is considered. The qualitative method for the phase space analysis is applied. The analysis of motion and synthesis of modes of SC motion with variable structure is carried out. Possible evolutions of the SC motion are defined. Numerical modeling of the SC motion for nonlinear cases of inertia-mass parameters dependences is carried out.
Spacecraft, dynamic system, phase space, variable structure, synthesis of motion modes, numerical modeling.
Информация об авторах
Балакин Виктор Леонидович, доктор технических наук, профессор кафедры космического машиностроения, Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет).
E-mail: [email protected]. Область научных интересов: динамика и управление движением летательных аппаратов.
Дорошин Антон Владимирович, кандидат технических наук, доцент кафедры космического машиностроения, Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет).
E-mail: [email protected]. Область научных интересов: динамика пространственного движения космических аппаратов. постоянного и переменного состава.
Крикунов Михаил Михайлович, аспирант кафедры космического машиностроения, Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет). E-mail: [email protected]. Область научных интересов: динамика пространственного движения космических аппаратов.
Victor Leonidovich Balakin, doctor of science (engineering), professor of space mechanical engineering department, Samara State Aerospace University named after academician S. P. Korolyov (National Research University). E-mail: [email protected]. Areas of research: dynamics and aircraft motion control.
Anton Vladimirovich Doroshin, candidate of science (engineering), associated professor of space mechanical engineering department, Samara State Aerospace University named after academician S. P. Korolyov (National Research University). E-mail: [email protected]. Areas of research: dynamics of spacecraft attitude motion of constant and variable structure.
Mikhail Mikhailovich Krikunov, post-graduate student of space mechanical engineering department, Samara State Aerospace University named after academician S. P. Korolyov (National Research University). E-mail: [email protected]. Areas of research: dynamics of spacecraft attitude motion.
