CC BY
83
УДК 531.36+ 629.78
ДИНАМИКА БУКСИРОВКИ КОСМИЧЕСКОГО АППАРАТА НА ТРОСЕ В ЦЕНТРАЛЬНОМ ПОЛЕ
© 2012 В. С. Асланов, В. В. Юдинцев
Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет)
Рассматривается движение космического аппарата как твёрдого тела, соединённого вязкоупругим тросом с космическим буксиром в центральном гравитационном поле. На основе формализма Лагранжа получены уравнения пространственного движения космического аппарата и орбитального движения центра масс системы в оскулирующих элементах. Изучено влияние на движение космического аппарата параметров системы и начальных условий движения. Показана возможность существования критических режимов движения, приводящих к запутыванию троса.
Космический буксир, упругий трос, твёрдое тело, центральное гравитационное поле, оскулирую-щие элементы.
1. Введение
Статья посвящена изучению пространственного движения космического аппарата (КА) как твёрдого тела на тросовой связи с КА-буксиром. Динамике твёрдого тела в космическом пространстве посвящено большое число работ. Отметим, в частности, монографию Белецкого В.В. [1] об орбитальном движении космических аппаратов относительно центра масс, монографии Ярошевского В. А. [2] и Асланова В. С. [3] о движении твёрдого тела в атмосфере. Движение КА относительно центра масс на упругой связи в составе космической тросовой системы рассматривалось в работах Асланова В. С. [4,5,6], Асланова В. С. и Ледкова А. С. [7], Асланова В. С. и Стратилатова Н. Р. [8], Забо-лотнова Ю. М. и Наумова О. Н. [9]. В этих работах описана и проанализирована ситуация, когда одно из концевых тел тросовой системы совершает заданное орбитальное движение.
В представленной статье рассматривается движение тросовой системы, когда одно из тел системы - КА-буксир является активным, оказывающим влияние как на движение всей системы, так и на движение вокруг центра масс твёрдого тела -пассивного КА. Это необходимо при детальном исследовании динамики системы
в целом и отдельных её элементов. На основе формализма Лагранжа формируются уравнения пространственного движения твёрдого тела на вязкоупругом тросе вокруг центра масс тела и относительно центра масс всей системы, а оскулирую-щие элементы используются для изучения эволюции орбиты системы. Предложенный в работе выбор обобщённых координат позволяет упростить анализ относительного движения пассивного КА и троса.
2. Уравнения движения
Будем рассматривать активный КА в виде материальной точки, пассивный КА как осесимметричное твёрдое тело, а трос упруговязким и невесомым. Предполагаем, что центры масс тел в процессе движения остаются в плоскости орбиты. Для построения математической модели используются четыре правые системы координат:
1. Инерциальная система координат О'хрУр^р, расположенная в точке перицентра с плоскостью О'хр>Р, совпадающей с плоскостью орбиты, осью О'ур, совпадающей с прямой, проведённой через центр притяжения и перицентр и осью
3. Подвижная система координат СхсУс^с с началом координат в центре масс пассивного КА: ось Схс направлена по прямой, соединяющей центры масс пассивного и активного КА, ось Сус параллельна плоскости Ох0уо.
4. Четвёртая система координат Сху2 является главной центральной системой координат пассивного КА (рис. 1).
Рис. 1. Схема тросовой системы
О 'хр, направленной по вектору орбитальной скорости в перицентре.
2. Подвижная орбитальная система координат Ох0у0^0 с началом в центре масс системы, осью Оу0, направленной от притягивающего центра и осью Ох0, лежащей в плоскости орбиты и направленной в сторону направления орбитальной скорости.
Для построения уравнений движения системы вокруг центра масс будем использовать формализм Лагранжа, а для уравнений движения центра масс системы - оскулирующие элементы.
2.1. Уравнения движения вокруг центра масс
Для определения положения системы относительно общего центра масс и положения пассивного КА относительно-
го собственного центра масс выберем следующие обобщённые координаты:
д1 = а - угол поворота системы координат Схсус2с относительно Ох0у0^0 вокруг оси Oz0 (рис. 1); д2 ° 5 - расстояние между центрами КА; д3 = у, д4 = в, q5 ° р -
углы, задающие положение базиса Cxyz, связанного с твёрдым телом, относительно базиса Схсу^с. Уравнения движения
запишем в форме уравнений Лагранжа 11го рода [10]:
й дТ дТ _ . 1 5 (1)
------------= о,. = 1,...5. (1)
Ж д^, дql 1 W
Кинетическая энергия системы
имеет вид:
Т = Т + Тю , (2)
где Ту, Тт - кинетические энергии поступательного и вращательного движений, соответственно.
Кинетическую энергию поступа-
тельного движения будем определять в виде:
Т
m1v1
+
mv
2
2 2
(3)
2 2 '
где индекс 1 отвечает активному КА, индекс 2 - пассивному КА.
Абсолютная скорость центров масс КА определяется как сумма скорости центра масс системы и скорости движения центров масс КА относительно центра масс системы:
У; = у0 + і = 1,2, (4)
где &Кі / & - локальная производная радиус-вектора положения і-го КА относительно центра масс системы в орбитальном подвижном базисе.
Координатный столбец вектора скорости центра масс системы г0 в проекциях
на оси орбитальной подвижной системы координат определяется следующим образом:
1 + e cos J
' (5)
v O =
e sin J 0
где i9 - угол истинной аномалии, p - фокальный параметр, e - эксцентриситет орбиты, m - гравитационный параметр. Радиус-векторы пассивного и активного КА относительно центра масс системы в проекциях на оси OxoyoZo имеют вид:
(6)
s, cos a -s, cos a
R 2 = s, sin a , Ri = s, sin a
O O
m2 m,
s1 = s--------2—, s2 = s------------------1
m1 + m2
m1 + m2
(7)
где s1, &2 - расстояния от центра масс системы до центров масс КА 1 и 2, соответственно; ^ - расстояние между центрами масс КА. Дифференцируя (6) с учётом формул (3), (4) и (7), получим кинетическую энергию поступательного движения:
m2m,
і f 2e 2m+4em cos j+2 m+p (s2 (a - j )2+s
+
+
2 (m + m2) p (В)
mm1 (e2 + 2ecos( J) +1) + mm22 (e2 + 2e cos( J) +1)
2 (mi + m2) p
Для определения кинетической энергии вращательного движения Тю в
выражении (2) найдём абсолютную угловую скорость твёрдого тела как сумму угловой скорости базиса Схсус2с и угловой скорости тела относительно этого базиса:
«2 = «с + «2с . (9)
Угловая скорость базиса Схсусхс определяется изменением угла а и угла истинной аномалии ®:
«с = к 0(а - (10)
где к0 - орт оси Сг0 базиса Охсусгс. Угловая скорость твёрдого тела относительно собственного базиса Схсус2с записывается следующим образом [11]:
j+y cosО О cos j+y sin О sin j -О sin j+y sin О cos j
2Tw = «2c • J • «2c .
(12)
(11)
Для осесимметричного твёрдого тела с двухосным эллипсоидом инерции 1=ё1а§(А, В, В) с учётом (8)-(10) и (11)
Кинетическая энергия вращательно- выражение Та примет вид: го движения твёрдого тела относительно
центра масс определяется как
Ta = — \a2 (- (со + J)sin О cosy + y cos О + j)2 +
+ B2 ((со +J)(cos О sin f cosy + cos f siny) + y sin О sin f + О cos + B2 ((со + J)(cos О cos f cos y - sin f sin y) + y sin О cos f - О sin
+
Обобщённые силы включают в себя силу тяги реактивного двигателя, приложенную к активному КА, силу натяжения троса и центральную гравитационную силу. Силу натяжения троса определим как
^ = И, [с, (I - /о) + к,/],
где I - длина троса, 10 - свободная длина троса, с, - жёсткость троса, к, - коэффициент демпфирования троса; функция Хевисайда:
Hv =
0, l < l
O ’
1, l > Io.
Радиус-вектор Я.12, соединяющий точку крепления троса на твёрдом теле с активным КА, определяется следующим образом:
К12 = + Р - ^1>
где р - радиус-вектор точки крепления троса на твёрдом теле относительно его центра масс.
l = R
12
Длина троса в общем случае равна = ^s2 - 2s[px cos 0 + sin 0(py sin j + pz cos j)] + p2.
Будем считать, что длина троса l существенно больше модуля вектора p:
l » s - px cos О - sin О (py sin j + pя cos j), скорость деформации троса равна l = s + О •\px sinО -(py sin j+pz cosj)cos О]+
+ j sinО (pz sin j - py cos j)
Вектор силы натяжения троса, действующий на тела 1 и 2, определим следующим образом:
F =-F = i P
2 lc1 ’
где ic - орт оси Cxc. Радиус-вектор точки приложения силы натяжения троса на пассивном КА равен:
Г 2 = Г2 + P .
Точка приложения силы натяжения троса на первом теле определяется радиус-вектором гп = г1.
Обобщённые силы записываются с помощью выражения [10]:
0г = ^ ^ + Я2 ^ + ^•-дГ1 +
Эqг
+¥ •* + р 2 А 1 = 1...7,
я1 дЧг *2 дч: ’
где Я1 - активная сила; Гё1, Гё2 - силы притяжения, действующие на тела 1 и 2:
р =-г 1 = 12
Г3 Г1, 1 1,2.
Найденное значение кинетической энергии системы и обобщённые силы по-
«2c
2
2
зволяют сформировать уравнения движения рассматриваемой системы.
2.2. Уравнения движения центра масс системы
Для записи уравнений движения центра масс системы воспользуемся оску-лирующими элементами [12]. Для возмущающей орбиту силы Я1=1сР, действую -щей вдоль оси Охс орбитальной подвижной системы координат, уравнения
движения примут вид:
d w 1 Ip [ r0
-----= -J—ax I 1 + —
dt e
m у p ,
W = W0 = const, i = i0 = const, d- = 2«xroV p / m,
de
dt
dJ
dt
= ax.\—
m
e rO
- +
p
;os J
-a
1 p
e\ ml
і+-0-p
ЗШ },
твёрдого тела и прямой СС1 в0 = р/3, угол ас = 0 и в = у& = 0.
3.1. Движение твёрдого тела при постоянно натянутом тросе и начальной угловой скоростьи (р0 = 0,05 с-, момент
инерции А = 3500 кг-м2
Пусть система совершает движение при постоянно натянутом тросе: 1>1с. Г рафик изменения угла нутации в на протяжении половины орбитального периода приведен на рис. 2. Начальная угловая скорость пассивного КА ф0 приводит к возникновению гироскопического эффекта, что при наличии возмущающей силы натяжения троса вызывает увеличение угла в. На рис. 3 представлен график изменения высоты орбиты под действием силы, приложенной к КА-буксиру.
где о - аргумент перицентра, О - долгота восходящего узла, 1 - наклонение орбиты, р - аргумент перицентра, е - эксцентриситет орбиты, } - угол истинной аномалии, ах=Р /(т1! т2) - проекция на ось Охс возмущающего ускорения.
3. Численное моделирование
С целью изучения влияния начальных условий и параметров системы на её поведение рассмотрим четыре случая, отличающиеся начальными значениями расстояния между КА, начальной угловой скоростью пассивного КА и его моментами инерции.
Пусть на КА-буксир действует постоянная сила, направленная в сторону, противоположную оси Ох0: Р = - 1 Н.
Система имеет следующие неизменные параметры: В = С = 15000 кг-м2, т1 = 500 кг, т2 = 3000 кг, рх = 1 м; с = 83 ГПа, к = 0,1 Н-с/м, 1о = 200 м, т = 3.987 1014 и элементы орбиты: ес = 0.01, рс = 6871 км, 10 = р/3, О0 = 0, ю0 = р/2, $0=0. Во всех рассматриваемых случаях в начальный момент угол между осью симметрии
Рис. 2. Изменение угла в по времени для первого случая
с
Рис. 3. Изменение высоты орбиты по времени для первого случая в течение 4 орбитальных периодов
3.2. Движение твёрдого тела при периодических ослаблениях троса
2
r
о
ІЗ
Рассмотрим случай, когда начальное расстояние между телами меньше свободной длины троса я0=190 м. Трос воздействует на пассивный КА в течение коротких промежутков времени. Импульсные воздействия троса приводят существенным изменениям угла в, и этот угол может превысить некоторое критическое значение в* (рис. 4), которое определяется формой пассивного КА.
Рис. 5. Изменение угла в по времени для третьего случая 3.4. Движение твёрдого тела при постоянно натянутом тросе и начальной угловой скорости <р0 = 0,1 с
Результаты численного моделирования, представленные на рис. 6, показывают, что при увеличении начальной угловой скорости собственного вращения КА в два раза по сравнению с первым случаем наблюдаются эффекты аналогичные тем, что представлены на рис. 5.
Таким образом, при некоторых значениях начальной угловой скорости собственного вращения пассивного КА, величины продольного момента инерции и начального натяжения троса могут возникать режимы движения, когда угол между линией троса и продольной осью КА превышает некоторое критическое значение в*, что может привести к зацеплению троса за элементы конструкции КА.
а.
Рис. 4. Изменение угла в по времени для второго случая
3.3. Движение твёрдого тела при постоянно натянутом тросе и большем моменте инерции пассивного КА
Рассмотрим случай, отличающийся от первого большим значением продольного момента инерции КА: А = 12000 кгм . Как показано на рис. 5, увеличение кинетического момента пассивного КА приводит к большему отклонению угла в от начального значения. Как и во втором случае, значение угла в в процессе движения превышает некоторое критическое значение в*.
7Г/3
\. конус Предельных Й >о* положений \ "Л£ троса трос 1 / /
ч
V У V
V
500
1000
1500
2000
2500 I, С
Рис. 6. Изменение угла в по времени для четвёртого случая
4. Заключение
На основе формализма Лагранжа построена математическая модель движения тросовой системы, состоящей из пассивного КА, вязкоупругого троса и космического буксира в центральном гравитационном поле. Анализ результатов численного моделирования показывает, что при буксировке КА необходимо учитывать его движение вокруг центра масс. В процессе движения пассивного КА возможны режимы движения, приводящие к большим колебаниям продольной оси КА, что может привести к запутыванию троса. Эти режимы определяются с начальными условиями движения и параметрами буксируемого КА.
Работа выполнена при поддержке Российского фонда фундаментальных исследований (проект 12-01-00317-а).
Библиографический список
1. Белецкий, В. В. Движение спутника относительно центра масс в гравитационном поле [Текст] / В. В. Белецкий -М.: Наука, 1975.
2. Ярошевский, В. А. Движение неуправляемого тела в атмосфере [Текст] / В. А. Ярошевский - М.: Машиностроение, 1978.
3. Асланов, В. С. Пространственное движение тела при спуске в атмосфере [Текст] / В. С. Асланов - М.: Физматлит, 2004.
4. Асланов, В. С. Колебания тела с орбитальной тросовой системой [Текст] /
B. С. Асланов // Прикладная математика и механика. - 2007. - Т. 71, Вып. 6. -
C. 1027-1033.
5. Асланов, В. С. Влияние упругости орбитальной тросовой системы на колебания спутника [Текст] / В. С. Асланов // Прикладная математика и механика. -2010. - Т. 74. Вып. 4. - С. 582-593.
6. Асланов, В. С. Колебания спутника с вертикальным упругим тросом на орбите [Текст] / В. С. Асланов // Известия РАН «Механика твердого тела». - 2011. -№5. - С. 3-15.
7. Асланов, В. С. Хаотические колебания КА с упругим радиально ориентированным тросом [Текст] / В. С. Асланов, А. С. Ледков // Космические исследования. - 2012. - Т. 50, № 2. - С. 194-204.
8. Асланов, В . С. Малые колебания
осесимметричного космического аппарата с тросовой системой [Текст] /
В. С. Асланов, Н. Р. Стратилатов // Вестник СамГУ - Естественнонаучная серия. Механика. - 2008. - 6 (65). - С. 202-208.
9. Заболотнов, Ю. М. Движение спускаемой капсулы относительно центра масс при развертывании орбитальной тросовой системы [Текст] / Ю. М. Заболотнов, О. Н. Наумов // Космические исследования. - 2012. - Т. 50 № 2. - С. 183-193.
10. Бухгольц, Н. Н. Основной курс теоретической механики [Текст]: учеб. для гос. ун-тов / Н. Н. Бухгольц. - 4-е изд., перераб. и доп. С. М. Таргом. - М.: Наука, 1966. Ч. 2: Динамика системы материальных точек.
11. Виттенбург, Й. Динамика систем твёрдых тел [Текст]/ Й. Виттенбург - М.: Мир, 1980.
12. Охоцимский, Д. Е. Основы механики космического полёта [Текст]: учеб. пособие / Д. Е. Охоцимский, Ю. Г. Сихарулидзе - М.: Наука. 1990.
DYNAMICS OF TOWING TETHERED SPACECRAFT IN THE CENTRAL
GRAVITATIONAL FIELD
© 2012 V. S. Aslanov, V. V. Yudintsev
Samara State Aerospace University named after academician S.P. Korolyov (National Research University)
The motion of spacecraft as a rigid body connected to a space tow with a viscoelastic tether in the central gravitational field is discussed in the paper. The equations of spacecraft spatial motion and the orbital motion of the system’s center of mass in osculating elements are obtained on the basis of Lagrange formalism. The influence of the system’s parameters and the initial conditions of motion on the motion of the spacecraft is analyzed. The possibility of critical modes of motion leading to the entanglement of the tether is shown.
Space tow, elastic tether, rigid body, central gravitational field, osculating elements.
Информация об авторах
Асланов Владимир Степанович, доктор технических наук, профессор, заведующий кафедрой теоретической механики, Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет). E-mail: aslanov [email protected]. Область научных интересов: динамика твёрдого тела, динамика систем тел, хаотическая динамика, динамика космических тросовых систем.
Юдинцев Вадим Вячеславович, кандидат технических наук, доцент кафедры теоретической механики, Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет). E-mail: [email protected]. Область научных интересов: динамика твёрдого тела, динамика систем тел, хаотическая динамика, динамика космических тросовых систем.
Aslanov Vladimir Stepanovich, doctor of technical science, professor, head of the department of theoretical mechanics, Samara State Aerospace University named after academician S.P. Korolyov (National Research University). E-mail: aslanov [email protected]. Area of research: rigid-body dynamics, dynamics of body systems, chaotic dynamics, dynamics of space tether systems.
Yudintsev Vadim Vyacheslavovich, candidate of technical science, associate professor, the department of theoretical mechanics, Samara State Aerospace University named after academician S.P. Korolyov (National Research University). E-mail: [email protected]. Area of research: rigid-body dynamics, dynamics of body systems, chaotic dynamics, dynamics of space tether systems.
