Вестник СамГУ — Естественнонаучная серия. 2013. № 3(104)
МЕХАНИКА
УДК 531.36
ДИНАМИКА БУКСИРОВКИ ТВЕРДОГО ТЕЛА НА УПРУГОМ ТРОСЕ В БЕЗГРАВИТАЦИОННОМ
ПРОСТРАНСТВЕ1
© 2013 В.С. Асланов, В.В. Юдинцев2
Рассматривается движение твердого тела, соединенного упругим тросом с космическим аппаратом - буксиром, в безгравитационном пространстве. Получены уравнения пространственного движения твердого тела и изучено влияние на его поведение начальных условий движения. Показана возможность существования критических режимов движения, которые могут приводить к запутыванию троса.
Ключевые слова: твердое тело, космический буксир, упругий трос.
Введение
Статья посвящена изучению пространственного движения твердого тела на тросовой связи с космическим буксиром. Динамике твердого тела в космическом пространстве посвящено большое число работ: отметим монографию В.В. Белецкого [1], посвященную орбитальному движению космических аппаратов (КА), монографии В.А. Ярошевского [2] и В.С. Асланова [3], в которых изучается движение твердого тела в атмосфере. Движение КА относительно центра масс на упругой связи в составе космической тросовой системы рассматривались в работах В.С. Асланова [4-6], В.С. Асланова и А.С. Ледкова [7], В.С. Асланова и Н.Р. Стратилатова [8], Ю.М. Заболотнова и О.Н. Наумова [9]. В работах [4-9] одно из концевых тел системы совершало заданное орбитальное движение. В представленной статье рассматривается движение тросовой системы, когда одно из тел системы - КА-буксир — является активным, оказывающим влияние как на движение всей системы, так и на движение вокруг центра масс твердого тела -пассивного осесимметричного КА. Уравнения пространственного движения твердого тела на вязкоупругом тросе вокруг центра масс и относительно центра масс всей системы строятся на основе формализма Лагранжа. Обобщенные координаты выбраны так, что упрощаются анализ относительного движения твердого тела и
1 Работа выполнена при поддержке Российского фонда фундаментальных исследований (проект 12-01-00317-а).
2Асланов Владимир Степанович (aslanov_vs@mail.ru), Юдинцев Вадим Вячеславович (yudintsev@classmech.ru), кафедра теоретической механики Самарского государственного аэрокосмического университета (Национального исследовательского университета), 443086, Российская Федерация, г. Самара, Московское шоссе, 34.
решение прикладных задач, связанных с оценкой безударности транспортных операций, совершаемой активным КА. Численное моделирование позволяет выявить режимы движения твердого тела, характеризующиеся большой амплитудой колебаний, приводящей к запутыванию троса.
1. Уравнения движения
При построении модели используются четыре системы координат. Инерциаль-ная система координат О' х'у'г': центры масс тел в процессе движения остаются в плоскости О'х'0у'0. Подвижная система координат Охоуо^о, расположенная в центре масс системы и движущаяся поступательно с центром масс, относительно О' х'0у'0г'0 (оси Ох, Оу, Ох параллельны соответствующим осям О'х', О 'у', О 'г'). Подвижная система координат, связанная с прямой, соединяющей центры масс тел Схсусгс: ось хс направлена от центра масс пассивного тела к активному телу, Сус параллельна плоскости Охоуо, а ось Сгс дополняет систему координат до правой. Четвертая система координат Схуг связана с твердым телом.
Для записи уравнений движения системы воспользуемся формализмом Лагран-жа [10]:
а дТ дТ ^ . 1
аЩ - дц = ^ = 7 (1)
В качестве обобщенных координат выберем (см. рис. 1): две координаты центра масс системы в О'х^'у'г' (х = ц\,у = ц^); угол поворота системы координат Схсусгс относительно Ох'у'г' вокруг оси г' (а = цз); расстояние между центром масс твердого тела и активным КА (в = Ц4) и угловое положение базиса Схуг, связанного с твердым телом, относительно базиса Схсусгс (ф = Ц5, в = Цб, ф = 9г).
1.1. Кинетическая энергия системы
Рассматривается космическая тросовая система, состоящая из двух тел: активного КА — материальной точки 1, пассивного КА — твердого тела 2 и вязкоупру-гого невесомого троса. Кинетическая энергия системы имеет вид:
Т = ТУ + Тш, (2)
где Тш - кинетическая энергия вращательного движения твердого тела вокруг своего центра масс; Ту - кинетическая энергия поступательного движения системы:
_ т1у21 | (3)
= 2 + 2 • (3) Абсолютные скорости центров масс КА VI, V2 определяются как суммы скоростей центра масс системы и скоростей движения центров масс КА относительно центра масс системы:
Vi = V' + ,1 = 1, 2. (4)
аг
Радиус-вектор скорости центра масс системы V' в проекциях на оси неподвижной системы координат О х' у' г' определяется следующим образом:
У 0
Здесь и далее при записи координатных столбцов векторов верхним индексом в скобках будем обозначать систему координат, в которой записывается вектор. Координатный столбец радиуса-вектора положения пассивного КА относительно центра масс системы в проекциях на оси Ох^уо^о имеет вид:
r:
(0)
S2 cos а s2 sin а 0
активного КА:
R
(0)
— si cos а si sin а 0
где si,s2 — расстояния от центра масс системы до центров масс КА 1 и ответственно равные:
m2 mi
si = s-
S2 = s-
(6)
(7)
со-
(8)
mi + m-2 mi + m-2
где s — расстояние между центрами масс КА. Подставляя (6), (7) и (8) в (4) получим вектор скорости активного КА:
m-2 (ás sin а + s( — cos а) + X) + mix m-2 (—s sin а — аs cos а + y) + miy , (9)
0
V'
(0)
1
mi + m-2
и вектор скорости пассивного КА:
mi (—as sin а + S cos а + x) + m^xc mi (S sin a + as cos a + y) + m2y 0
V.
(0)
1
TO1 + Ш2
(10)
Подставляя (9) и (10) в (3), получим кинетическую энергию поступательного движения системы:
T
T v
то2TOi (á2s2 + s2 + 2^2) + v2(to2 + to2)
(11)
2 (т1 + т2)
Для определения составляющей Тш в выражении (2) найдем угловую скорость твердого тела, которая складывается из угловой скорости базиса Схсусгс и угловой скорости этого тела относительно базиса Схсусгс:
ш2
22) = + = aMc) +
ш
(2)
(12)
где верхний индекс в 2 скобках указывает на систему координат Схуг, связанную с твердым телом; А2с - матрица преобразования координат из базиса Схсусгс в базис Схуг [11]:
A2
ce se sv se cv
ses-ф сфсф ce s^ ce
se сф cesvc^ + сф s-ф ce сфсф sф s^
где ce = cos в, se = sin в, c^ = cos ф, s^ = sin ф, сф = cos ф, вф = sin ф. Вектор угловой скорости базиса Cxcyczc определяется следующим образом:
,(c)
Ac
0 0
а
(13)
где Ac
матрица преобразования координат из базиса Ox0y0z0 в базис Cxcyczc:
— cos а — sin а 0 Aco = sin а cos а 0 . (14)
0 0 1
Вектор угловой скорости пассивного КА относительно системы координат Cxcyczc запишем, используя углы Эйлера и их производные [11]:
ш
(2)
ф + ф cos в в cos ф + ф sin в sin ф — в sin ф + ф sin в cos ф
(15)
Кинетическая энергия движения твердого тела относительно собственного центра масс будет определяться формулой:
Тш = 42)т л(2Ц2), (16)
где Л(2) — тензор инерции твердого тела в базисе Схуг: Л(2) = а%ад(Л,В,В). Подставляя в (16) (12), с учетом (13), (15) и (11) окончательно получим кинетическую энергию всей системы в виде:
2T =
i
mi+m2
[то2то1(а2в2 + s2 + 2v°) + vO (то2 + TO2 )] +
+А2(ф cos в — а sin в cos ф + ф)2 + В2(ф2 sin2 в + в2)+ +В2 а2 (cos2 в cos2 ф + sin2 ф) + а(2в sin ф + ф sin2в cos ф) .
2
c
c
c
2
c
1.2. Обобщенные силы
На активное тело действуют сила тяги Fi и сила натяжения троса Fti, на пассивное тело только сила натяжения троса Ft2. Сила натяжения вязкоупругого троса определяется формулой:
Ft = Hv [ct(l — lo)+ M, (18)
где l - длина троса; l0 - длина троса в недеформированном состоянии; ct - жесткость троса; kt - коэффициент демпфирования троса; l = dl/dt; Hv - функция Хевисайда
0, l < l0, 1, l>lo.
Радиус-вектор Ri2, соединяющий центры масс КА, запишем в виде:
RÍO) = r20) + AocAc2р(2) - Ri0),
где р(2) - координатный столбец радиуса-вектора точки крепления троса на твердом теле относительно системы координат Cxyz. Длина троса равна
|R(2)| = l = \Js2 — 2s [px cos в + sin в (py sin ф + pz cos ф)] + p2.
Будем предполагать, что длина троса существенно больше модуля вектора р: s ^ |р|, тогда приближенно будем полагать, что
l « s — px cos в — sin в (py sin ф + pz cos ф). (19)
Дифференцируя (19), получим скорость деформации троса
l = s + в [px sin в — cos в (py sin ф + pz cos ф)] + ф sin e(pz sin ф — py cos ф). (20)
Тогда вектор силы натяжения имеет вид:
F(0) = —F(0) = AoceXc)Ft,
где eXc) - орт оси x системы координат Cxcyczc. Координатный столбец точки приложения силы на твердом теле определяется в виде:
rt0) = r20) + AocAc2p(2). (21)
По определению обобщенной силы [10] имеем:
, , ЯГ(0) „ „ dr(0) „ dr(0)
Q = F(0) drt1 + F(0) t2 + F(0) dr1 ■ = . 7
Qi = Fti • + Ft2 ■ d" + Fi • ,г = 7.
Будем рассматривать случай закрепления троса на оси симметрии твердого тела
P(2) = (px, 0,0)T, тогда окончательно выражения для обобщенных сил примут вид:
Q1 = px sin в sin ф ct (L0 — s + px cos в) — k(s + вpx sin в) Hv +--2—--—, (22)
L J mi + m2
ct (s — L0 — cos^)px) + ks + вк sin(в)px Hv, (23)
Q3 = Q5 = Qr = 0, (24)
Q4 = sin вpx ct (L0 — s + cos вpx) — ks — вк sin^p Hv, (25)
Qe = Fix. (26)
Найденное значение кинетической энергии системы (17) и обобщенные силы (22)-(26) позволяют сформировать уравнения движения рассматриваемой системы (1).
m-2 cos(a)Fxi
q2 = -¡-
mi + m2
2. Численный анализ
Для проведения численного моделирования будем использовать следующие параметры системы: A = 3400 кг-м2, B = C = 15000 кг-м2, mi = 500 кг, m2 = 3000 кг, px = 1 м; ct = 10 Н/м, kt = 1 Н-с/м, 1о = 200 м. На КА-буксир действует постоянная сила, направленная в сторону противоположную оси Oxo: Fix = — 1Н, Fiy = 0. Рассмотрим четыре случая, отличающиеся начальными значениями расстояния между КА и углом а между линией действия активной силы и прямой, соединяющей центры масс КА. В начальный момент угол между осью симметрии твердого тела и прямой CCi во = 30°, и начальная угловая скорость тела равна: фо = 0,1 рад/с, фо = во = 0.
2.1. Движение твердого тела при периодическом ослаблении троса, а0 = 0
Рассмотрим случай, когда начальное расстояние между телами меньше свободной длины троса so = 190 м, и начальный угол между линией, соединяющей центры масс тел, и линией действия активной силы Fi равен: ао = 30°. На рис. 2, а приведена зависимость изменения длины троса от времени. Из графика следует, что при выбранных параметрах системы и начальных условиях большую часть времени трос находится в свободном состоянии, растягиваясь и воздействуя на твердое тело в течение коротких временных интервалов. Эти импульсные воздействия приводят и изменению угла между продольной осью твердого тела и линией, соединяющей центры масс тел, что иллюстрирует график на рис. 2, б. Отметим, что малая начальная угловая скорость вращения тела фо = 0,1 рад/с приводит к возникновению гироскопического эффекта и значительному увеличению угла в: после 1700 секунды этот угол становится больше 90°, что может привести к запутыванию троса для некоторых конструктивных схем пассивного КА.
I, с ^ с
а б
Рис. 2. Изменение длины троса и угла в для первого случая
2.2. Движение тела при периодическом ослаблении троса,
а0 = 0
Второй случай отличается от первого начальной ориентацией троса относительно линии действия активной силы: в начальный момент времени ао = 0. Ре-
зультаты расчетов приведены на рис. 3, где показан характер изменения длины троса и угла в.
200 198 196 194 192 190
0 500 1000 1500 2000 2500 3000 ^ с
0 500 1000 1500 2000 2500 3000 ^ с
а б
Рис. 3. Изменение длины троса и угла в для второго случая
Как следует из рис. 3, б, амплитуда колебаний угла между продольной осью твердого тела и линий СС1, как и в первом случае, достигает больших значений.
2.3. Движение тела с постоянно натянутым тросом, а = 0
Рассмотрим движение системы в случае, когда трос всегда натянут: / > /о. Начальные условия движения, за исключением величины во, которую здесь примем равной 200,9 м, соответствуют первому случаю, т. е. ао = 30°. Как следует из рис. 4, а, трос на всем рассматриваемом интервале времени остается растянутым, затухание высокочастотной составляющей колебания троса обусловлено наличием демпфирования ^ > 0. Величина угла в также достигает больших значений (см. рис. 4, б).
200.12 200.10 200.08 -- 200.06 200.04 200.02 200.00
..[ЩИ!. I .....
Ш11111||. т.__^«ч :
100
80
60
40
20
0 500 1000 1500 2000 2500 3000 ^ с
0 500 1000 1500 2000 2500 3000 ^ с
б
Рис. 4. Изменение длины троса и угла в для третьего случая
а
2.4. Движение тела постоянно натянутым тросом, а = 0
Рассмотрим движение системы, когда трос всегда натянут: / > /о и начальное значение угла ад = 0. Результаты, приведенные на рис. 5, а, показывают,
что в этом случае колебания троса не возмущаются дополнительными колебаниями угла а, что приводит к существенно меньшей амплитуде колебаний угла в (см. рис. 5, б) по сравнению с предыдущими случаями.
0 500 1000 1500 2000 2500 3000 t, c
0 500 1000 1500 2000 2500 3000 t, c
б
Рис. 5. Изменение длины троса и угла в для четвертого случая
а
Заключение
На основе формализма Лагранжа построена математическая модель движения твердого тела на вязкоупругой тросовой связи с КА-буксиром в безгравитационном пространстве. Анализ результатов численного моделирования показал, что при разработке транспортных тросовых космических систем необходимо учитывать движение вокруг центра масс неуправляемого твердого тела, перемещаемого при помощи космического буксира. В процессе движения пассивного твердого тела возможны режимы движения, приводящие к большим отклонениям угла в между продольной осью твердого тела и прямой, соединяющей центр масс твердого тела и буксир, что может приводить к запутыванию троса. Возникновение таких нежелательных режимов движения определяется начальными условиями движения твердого тела и взаимным положением тел системы в момент начала действия активной силы.
Литература
[1] Белецкий В.В. Движение спутника относительно центра масс в гравитационном поле. М.: Наука, 1975.
[2] Ярошевский В.А. Движение неуправляемого тела в атмосфере. М.: Машиностроение, 1978. 167 с.
[3] Асланов В.С. Пространственное движение тела при спуске в атмосфере. М.: Физматлит, 2004. 160 с.
[4] Асланов В.С. Колебания тела с орбитальной тросовой системой // Прикладная математика и механика. 2007. № 71(6). С. 1027-1033.
[5] Асланов В.С. Влияние упругости орбитальной тросовой системы на колебания спутника // Прикладная математика и механика. 2010. № 74(4). С. 582-593.
[6] Асланов В.С. Колебания спутника с вертикальным упругим тросом на орбите // Известия РАН. Сер.: Механика твердого тела. 2011. № 5. С. 3-15.
[7] Асланов В.С., Ледков А.С. Хаотические колебания КА c упругим радиаль-но ориентированным тросом // Космические исследования. 2012. 50. (2). С. 194-204.
[8] Асланов В.С., Стратилатов Н.Р. Малые колебания осесимметричного космического аппарата с тросовой системой // Вестник СамГУ. Естественнонаучная серия. 2008. № 6(65). С. 202-208.
[9] Заболотнов Ю.М., Наумов О.Н. Движение спускаемой капсулы относительно центра масс при развертывании орбитальной тросовой системы // Космические исследования. 2012. 50 (2). С. 183-193.
[10] Бухгольц Н.Н. Основной курс теоретической механики. Т. 2. Динамика системы материальных точек. 4-е изд. М.: Наука, 1966.
[11] Виттенбург Й. Динамика систем твердых тел. М.: Мир, 1980.
Поступила в редакцию 22/X7/2012; в окончательном варианте — 22/X7/2012.
THE DYNAMICS OF TOWAGE BY THE ELASTIC ROPE IN A GRAVITATIONLESS SPACE
© 2013 V.S. Aslanov, V.V. Yudintsev3
The dynamics of the rigid body connected with the space tug by the elastic rope in a gravitationless space is considered. The equations of motion are derived from a Lagrange formalism. The influence of initial conditions on the rigid body behaviour is studied. The posibility of existence of critical behaviour of movement which can cause the entanglement of rope is shown.
Key words: rigid body, space tug, elastic rope.
Paper received 22/XI/2012. Paper accepted 22/XI/2012.
3Aslanov Vladimir Stepanovich (aslanov_vs@mail.ru), Yudintsev Vadim Vyacheslavovich (yudintsev@classmech.ru), the Dept. of Theoretical Mechanics, Samara State Aerospace University, Samara, 443086, Russian Federation.