Вестник СПбГУ. Математика. Механика. Астрономия. 2020. Т. 7 (65). Вып. 3 УДК 517.93:531.391.3 МБС 70Л25, 70К40
Влияние неоднородности магнитного поля Земли на динамику электродинамической тросовой системы*
А. П. Дериглазов, А. А. Тихонов
Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
Для цитирования: Дериглазов А. П., Тихонов А. А. Влияние неоднородности магнитного поля Земли на динамику электродинамической тросовой системы // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2020. Т. 7(65). Вып. 3. С. 539-551. https://doi.org/10.21638/spbu01.2020.316
Рассматривается электродинамическая тросовая система, функционирующая в режиме натянутой связки и движущаяся по круговой экваториальной околоземной орбите в условиях неоднородности магнитного поля Земли. Найден возможный режим движения связки, представляющий собой близкое к вертикальному равновесное положение в плоскости орбиты. Исследована устойчивость найденного положения относительного равновесия по линейному приближению. Представлены результаты численного моделирования некоторых соотношений для тросовой системы при варьировании ее параметров.
Ключевые слова: электродинамическая тросовая система, устойчивость, геомагнитное поле, космический мусор.
1. Введение. Космические тросы рассматриваются в настоящее время не только как способ организации механической связи между орбитальными объектами, но и как источники дополнительных сил, действующих со стороны геофизических полей [1, 2]. Так, проводящий трос, движущийся в геомагнитном поле, испытывает воздействие амперовых сил, которые могут оказывать тормозящее действие на трос и присоединенные к нему орбитальные объекты [1]. Это обстоятельство породило целое направление исследований, в которых изучаются динамические свойства орбитальных связок с проводящими тросами, называемых электродинамическими тросовыми системами (ЭДТС). Особую актуальность этим исследованиям придает возможность использования ЭДТС для удаления с околоземных орбит нефункциониру-ющих спутников, последних ступеней ракет-носителей и других крупных объектов, называемых космическим мусором и представляющих опасность для действующих спутников [3—11]. При этом наибольшей эффективностью отличается проводящий трос, функционирующий в режиме проводника с током, ориентированного в околоземном пространстве по местной вертикали [1, 3, 12]. Такой режим равновесной ориентации для обычного весомого троса обеспечивается силами гравитационного притяжения Земли и силами инерции. Однако для ЭДТС устойчивая вертикальная ориентация троса разрушается под действием момента сил Ампера [1, 8, 12-14].
* Работа выполнена при финансовой поддержке РФФИ (грант № 19-01-00146-а). (¡5 Санкт-Петербургский государственный университет, 2020
Рис. 1. Схема ЭДТС.
Для борьбы с этим нежелательным явлением в работах [15-17] был предложен способ, основанный на разделении зарядов, протекающих в ЭДТС, и концентрации их на концевых коллекторах зарядов. Суть этого способа поясняется с помощью рис. 1. К концам проводящего изолированного троса 5 присоединены коллекторы — устройства для сбора электрических зарядов. Коллектор 7, находящийся на верхнем конце троса (на том, который дальше от Земли), с помощью электроизолирующих креплений 2 соединен с концевым телом 6 проводящего троса 5. Коллектор 7 получает положительный заряд, поддерживаемый одним или несколькими электронными эмиттерами (например, полевыми электронными эмиттерами, холодными электронными эмиттерами на основе нанопористого углерода, или холловскими ионными источниками) 3, передающими отрицательный заряд на концевое тело 6. На противоположном конце троса (на том, который ближе к Земле) концевое тело 4 аналогичным образом соединено с коллектором 1. Коллектор 1 получает отрицательный заряд, поддерживаемый одним или несколькими электронными эмиттерами 3, передающими отрицательный заряд с концевого тела 4. Избыток отрицательных зарядов стекает с коллектора 1 в окружающую околоземную плазму, что обеспечивает в тросе наличие тока, протекающего по псевдоцепи, замыкающейся через ионосферные токи, текущие вдоль силовых линий геомагнитного поля [1].
В результате разделения зарядов и концентрации их на коллекторах возникает дополнительный восстанавливающий момент сил Лоренца, который, как показано, может быть использован для стабилизации проводящего троса в околоземном пространстве вдоль местной вертикали [15-17]. Исследования, выполненные в этих работах, а также в других работах, процитированных выше, проводились, конечно, с учетом градиентности гравитационного поля Земли и возникающего вследствие этого гравитационного момента. Вместе с тем, градиентность магнитного поля Земли (МПЗ) не учитывалась ни в одной из указанных работ. Поэтому представляется естественным поставить вопрос о влиянии градиентности МПЗ на динамику ЭДТС, особенно с учетом большой протяженности троса — порядка нескольких километров. В данной работе рассматривается ЭДТС, взаимодействующая с МПЗ за счет протекающего по тросу тока и наличия разноименных зарядов на концевых коллекторах. Анализируется влияние градиентности МПЗ на существование положений равновесия ЭДТС, в том числе вертикального, их устойчивость и колебания.
2. Постановка задачи. Рассматривается ЭДТС, центр масс которой движется по круговой кеплеровой экваториальной орбите в гравитационном и магнитном полях Земли. В качестве гравитационного поля рассматривается центральное ньютоновское поле, МПЗ моделируется прямым магнитным диполем. Такая простейшая постановка задачи позволяет исключить нестационарности почти периодического типа [18, 19] и сосредоточиться именно на выявлении влияния градиентности МПЗ [20] на динамику ЭДТС. Сама ЭДТС представляет собой натянутый проводящий трос с массами, прикрепленными к его концам. В дальнейшем для краткости данная
Рис. 2. ЭДТС в орбитальной системе координат.
конструкция будет называться связкой. Концевые массы могут представлять собой два искусственных спутника Земли. Ввиду малых размеров спутников по сравнению с длиной троса, спутники рассматриваются как материальные точки, каждая из которых обладает массой и электростатическим зарядом. Будем предполагать, что трос находится в режиме, близком к натянутому состоянию вдоль местной вертикали благодаря градиенту гравитационного поля Земли. Данный режим движения троса будем называть номинальным. В этом режиме движения ось г, направленная вдоль натянутого троса, коллинеарна оси С£, направленной вдоль радиуса-вектора Е = О—С = Е(о центра масс троса относительно центра Земли Ое (рис. 2).
Оси СС (орт Со) и Сп (орт г/о), направленные соответственно по касательной к орбите в сторону движения точки С и по нормали к плоскости орбиты, образуют с осью СС (орт /о) орбитальную систему координат ССпС В инерциальном пространстве орбитальная система координат поворачивается с угловой скоростью Шо = ¿оСо. Упругие деформации троса не рассматриваются [1]. Таким образом, рассмотрению подлежит ЭДТС, моделируемая тонким проводящим однородным стержнем с погонной массой ро и точечными грузами на концах. Масса нити то = ро(г2 — гх), где г1 и г2 — координаты концевых масс т\ и т-2 относительно центра масс системы (положения концевых масс отмечены точками Сх и С2 на рис. 2). Координаты центров зарядов и ^2 будем считать совпадающими с координатами точечных масс. В системе главных центральных осей инерции Схуг с ортами ¿1, ¿2, ¿з тензор инер-
т 1- / л л л А г2 + гхг2 + г2 2 2
ции связки имеет вид .У = ша£(А, А, Со), где А = то---- + т^^ + т^гн, а
Со пренебрежимо мало по сравнению с А. Взаимную ориентацию осей систем координат ССпС и Схуг зададим с помощью матрицы направляющих косинусов
(«1 «2 «з\
в1 в2 вз I (1)
71 72 7з /
так, что Со = «1«1 + «2«2 + «з«з, По = в1С1 + в2г2 + взгз, Со = 71«1 + 72«2 + 7з«з. Поскольку рассматривается симметричная относительно продольной оси конструкция троса, то для устранения неопределенности в выборе осей х, у представляется целесообразным ввести в рассмотрение сопутствующие оси (оси Резаля) Схуг с ортами '1,3,к так, что ось С г направлена вдоль натянутого троса, а трехгранник Схуг не
участвует в повороте троса вокруг оси Cz на угол ф. Кинематическими характеристиками вращательного движения связки являются: абсолютная угловая скорость ш, угловая скорость сопутствующего трехгранника относительно орбитальной системы координат jji = pi + qj + rk, угловая скорость связки относительно сопутствующего трехгранника фк, угловая скорость связки относительно орбитальной системы координат ш = uj\ + фк. Эти величины связаны соотношением j = ш + jо, которое в проекциях на оси Cxyz имеет вид
Шх = Р + шовь jy = q + шов2, jz = Г + ф + шовз, (2)
где p, q,r — проекции относительной угловой скорости ш .
3. Дифференциальные уравнения вращательного движения связки.
Оси Резаля остаются главными центральными осями инерции ЭДТС во все время движения. Согласно теореме об изменении кинетического момента K = Ajxi + Ajyj под действием главного момента M внешних сил относительно центра масс имеем
+ (ji + ш0) х K = M.
xyz
Отсюда в проекциях на оси x, y, z получим следующую систему дифференциальных уравнений:
' Ашх - Ajy jz = MGx + MLx + Max
Aj y + Ajz Шх = MGy + MLy + MAy, (3)
jz = jz0 = const.
Совместно с кинематическими уравнениями Пуассона
a i + jy аз - jz «2 = - joYi,
$1 + jyвз - jzв2 =0, (J ^ У ^ ¡) (4)
Yi + jyY3 - jzY2 = joai
уравнения (3) образуют замкнутую дифференциальную систему, описывающую вращательное движение связки. Чтобы записать уравнения (3) в развернутом виде, нужно определить их правые части.
В центральном ньютоновском гравитационном поле на связку действует гравитационный момент Mg, который в проекциях на оси Cxyz имеет вид [21]
Mgx = -3j2AY2 Y3, MGy = 3j2AYi Y3, MGZ = 0.
Для вычисления моментов сил Ампера и Лоренца с учетом градиентности МПЗ определим сначала вектор магнитной индукции B в виде [22]
B = BC + ABB,
где Bс — магнитная индукция МПЗ в центре масс связки, а AB — дополнительное слагаемое, вычисляемое для данной точки пространства. Для МПЗ, моделируемого
( RE\ 3
прямым магнитным диполем [21-23], Вс = — д® ( J ffo, где дJ — гауссова постоянная, имеющая на эпоху 2015 значение gi = -29442.0 нТл. В точке троса, определяемой координатой z, составляющая AB, порождаемая градиентностью МПЗ, в
осях Схуг имеет вид [22]
ез (виз + взц\ АВ (г) = Зд^г /327з + /З372 • \ 2вз7з )
Как показано в [24], для принятой модели геомагнитного поля при вычислении момента сил Лоренца достаточно рассмотреть следующее его представление:
Мь = М ю + ДМ 1 + М2о + ДМ2.
Здесь первая составляющая лоренцева момента обусловлена орбитальным движением связки, вторая составляющая учитывает градиентность МПЗ, третья составляющая обусловлена движением ЭДТС относительно ее центра масс, четвертая составляющая учитывает градиентность МПЗ при движении ЭДТС относительно ее центра масс. При дипольном приближении МПЗ (прямой диполь) аналитические выражения для указанных составляющих таковы [25]:
Мю = Р х [На х Во) , М2о = ^2 + [(ш — ¿е) х к] (Вс • к) ,
ДМ = QlCtJl х (¿а х ДВ1) + Q2CC2 х (¿а х ДЁ^ ,
ДВ2 • к .
ДМ2 = д1г2 (ш — сСе) хк (АВ1 • к) + ^2г| (сС — ССе) хк
Здесь Р = (^1г1 + Q2г2)k, г>а = Е(шо — шЕ)Со, шЕ — угловая скорость суточного вращения Земли, ДВ^ = ДВ(г^) (г = 1, 2), векторы С—С^ направлены от центра масс ЭДТС к концевым массам (рис. 2), р = гк.
Момент сил Ампера вычислим по формуле [26]
МА = МА1о + ДМА1 = J р х (1к х Ва )йг + ^ р х (1к х ДВ )йг, где I — сила тока, протекающего по тросу. Вводя обозначения
С = 3Аш2
о, а = ~91:
о 1 { Ее\ / 2 2) от ЕЕ („з
2Е
Ез
Е2 6до ( Ее
Ее
Ее
Е
Ь
о=
691 ( Ее
Е \ Е
(Qlг3 + Q2z3) ,
запишем систему (3) в виде
Ашх —Ашушг = —С^21з — Ь^2 — Ь (2в2вз — 3727з) + Ьо^вз—
—Ьв2вз7з — Ьо^вз7з — ав1 + а (в17з + вз71), Аш у + Ашг шх = Су17з + Ь71 + Ь (2в1вз — 3717з) — Ьорвз+
+Ьв1вз7з + Ьорвз!з — ав2 + а (в27з + вз72),
Ш2 = ш2о = СОПЭ^
(5)
з
1
з
з
з
з
z
С
í
X
Рис. 3. Углы ф и в, определяющие ориентацию связки.
п
4. Положения относительного равновесия. Для решения вопроса о существовании возможных положений равновесия связки в орбитальной системе координат, положим проекции относительной угловой скорости p, q равными нулю в уравнениях (2), (4), (5) и перейдем от направляющих косинусов к углам ф, в (рис. 3) по формулам
= sin ф, в2 = cos ф cos в, вз = — cos ф sin в, 7i =0, 72 = sin в, 7з = cos в. Тогда уравнения (5) примут вид
{—Aj2 cos2 ф cos в sin в = G cos в sin в + L sin в—
—L sin в cos в (3 + 2 cos2 ф) — L cos2 ф sin в cos в + a sin ф — a sin ф cos в, (g)
Aw2 sin ф cos ф sin в = 2L sin ф cos ф sin в + L sin ф cos ф sin в cos в+ +a cos в cos ф — a cos ф cos 2в.
Положению троса, вытянутого вдоль местной вертикали, соответствует значение угла в = 0. Подставляя это значение в систему (6), получаем уравнение
откуда следует, что данный режим движения является практически нереализуемым, поскольку уравнение (7) накладывает неприемлемо жесткие ограничения на значения z\ и х2. Действительно, поскольку R на несколько порядков больше величин |zi| и х2, то при самых малых отклонениях от требуемых значений х\, Х2 левая часть уравнения (7) сильно отличается от нуля.
Рассмотрим возможность существования «наклонных» положений равновесия связки в орбитальной системе координат. В силу сферической симметрии гравитационного поля и осевой симметрии МПЗ, следует проанализировать возможность существования положений равновесия связки в двух плоскостях орбитальной системы координат: в плоскости орбиты (£, Z) ив плоскости (ц, Z), ортогональной плоскости орбиты.
Полагая sin ф = 0 в системе (6), получим следующую систему уравнений для нахождения положений равновесия связки в плоскости (r¡, Z), перпендикулярной к
(7)
плоскости орбиты:
'coseh2=(R(z1 + z2)±v/(z1 + z2)2R2 + 32(z2 + z1z2 + z2)2^/8(z21+z1z2 + z¡), cos в = —L/ - 5L - L
Несложно заметить, что во втором уравнении системы (8) правая часть является отрицательной в силу того, что L — положительная величина при любых параметрах связки, соблюдающих физически обоснованное расположение положительного заряда выше отрицательного, а знаменатель — положительная величина в силу того,
что Aw2 на несколько порядков больше, чем L и L. Таким образом в плоскости, перпендикулярной плоскости орбиты, отсутствуют равновесные положения связки, при которых положительный заряд выше отрицательного. «Перевернутые» же положения связки либо не представляют практического интереса в рамках поставленной задачи, либо вообще не реализуемы вследствие того, что такая связка будет смята кулоновыми и лоренцевыми силами, и поэтому далее не рассматриваются.
Положив cos ф = 0 в системе (6), перейдем к рассмотрению возможных положений равновесия связки в плоскости орбиты (£, Z). В данном случае одно из уравнений (6) удовлетворяется тождественно и система преобразуется к уравнению
G cos в sin в + L sin в — 3L sin в cos в ± а + S cos в = 0, (9)
где верхний и нижний знаки соответствуют ф = ±п/2. Перейдем в (9) от sin в к cos в, для краткости обозначим cos в = x. Для переменной x G [0; 1] получим следующее алгебраическое уравнение четвертого порядка:
p(x) = x4 (ü — 3L)2 + 2x3L (g — 3L) + x2 L2 — (g — 3L)2 + S2
— 2x l(g — 3:Lj — aS + a2 — L2 = 0. (10)
Для определения максимума возможного числа вещественных корней уравнения (10) на заданном интервале, воспользуемся теоремой Бюдана — Фурье [27]. Составим последовательность Фурье для полинома p (x) и его производных:
fseq(x) = |p(x),p'(x),p"(x),pm(x),p(IV^(x)| .
Число вещественных корней уравнения p (x) = 0, размещенных между 0 и 1, не может быть больше, чем число перемен знаков, потерянных в fseq (x) при переходе от подстановки x := 0 к подстановке x := 1. Для связок, симметричных (z2 = —zi) относительно центра масс, а = L = 0 и, как несложно проверить, получаем следующую таблицу знаков:
х = 0 х = 1
р(х) - +
р'(х) - +
р"( х) знак не важен +
р"'(х) + +
p(IV){x) + +
По представленным в таблице знакам можно утверждать, что симметричная относительно центра масс связка имеет единственное положение равновесия в плоскости орбиты.
Поскольку космическая тросовая система является сильно протяженным космическим объектом и условие 22 = —21 может выполняться не совсем точно, представляет интерес вопрос, существуют ли положения равновесия в плоскости орбиты для несимметричных относительно центра масс связок. Для этого проварьируем часть параметров ЭДТС при фиксированных значениях остальных. Параметрами ЭДТС являются: концевые массы и заряды, длина троса, погонная масса троса и сила тока. Будем считать фиксированными следующие из них: погонная масса троса р0 = 0.002 кг/м, сила тока I =2 А, заряды = —= 10-4 Кл. Поскольку за начало системы координат берется центр масс, можем определить связь
между координатами концевых масс через выражение 2-координаты центра масс:
г° = ~——I р°г<1г = ~——гпхгх + т2г2 + у- -
то + Ш1 + т2 у то + Ш1 + т2
ь
0.
Пусть величины и г2 связаны соотношением г-~т = М- Будем иметь следующее
ы
уравнение связи между массами и координатами 21 и 22:
+ ^ + (п)
р22 22 2 \ р2
Численные эксперименты с вариацией параметров выполнялись в двух вариантах. В первом из них производилась вариация параметра р на отрезке [0.9; 1.1] с шагом = 0.002 и вариация 22 на отрезке [500; 1500] м с шагом Л-2 =5 м. Строилась прямоугольная сетка (размера 100 х 100) параметров, для которых были вычислены корни уравнения (10) при фиксированном значении одной из концевых масс: т-2 = 100 кг. Величина т-1 вычислялась для заданных р и 22 по формуле (11). Результат вычислений показан в виде поверхности 0(р, 22) на рис.4, а.
Во втором варианте производилась вариация зарядов = —= Q на отрезке [10-5; 10-4] Кл с шагом = (10-4 — 10-5)/100 и масс т1 = т2 = т на отрезке [100; 600] кг с шагом Л-2 = 5 кг. При этом считались фиксированными следующие величины: сила тока I = 2 А, погонная масса троса 0.002 кг/м, 22 = —21 = 500 м. Результат вычислений показан в виде поверхности 0^, т) на рис.4, б. Выбор симметричной относительно центра масс связки обусловлен тем, что случай 22 = —21 обеспечивает интересующий нас минимум отклонения связки от вертикали (см. рис. 4, а).
По полученным значениям в пакете MATLAB были составлены кубические сплайны, для которых построены соответствующие поверхности. Из рис. 4, а, б видно, что для широкой области значений параметров существует положение равновесия ЭДТС в плоскости орбиты, которое соответствует значению угла 0, близкому к нулю. На графиках изменения угла при вариации 22 или Q являются неразличимыми, но стоит отметить, что при увеличении Q наблюдается уменьшение угла 0, ровно так же как и при увеличении 22 Следует также отметить, что в процессе расчетов проверялось условие натяжения троса. При всех рассмотренных значениях параметров связка находится в натянутом состоянии, несмотря на силы кулонова притяжения зарядов.
2
г
1
5. Устойчивость равновесного положения связки. Для определения условий устойчивости положения равновесия связки в плоскости орбиты вернемся к
уравнениям (5) и перейдем в них к углам 0,0 с учетом равенств p = 0, q = 0sin 0, • п
г = фсоъв. Будем полагать, что при значениях 0 = 0*, ф = ± — связка находится в
положении равновесия. Введем малые отклонения от этого положения Д0 и Д0 так, п
что 0 = 0*+ А0, ф = ±— + А0. С учетом малости углов А0 и А0 линеаризованная система (5) примет вид
АД0 = - (g - 3L) (cos 0* sin 0* + Д0 cos 20*) --L (sin 0* + Д0 cos 0*) T a ± a (cos 0* - Д0 sin 0*),
(АД0 + A^gД0) sin0* = = ( 2L sin 0* + L cos 0* sin 0* ± a cos 0* T a cos 20* ) Д0.
(12)
С учетом уравнения (9), которому удовлетворяют координаты положения равновесия, система (12) допускает упрощение до следующего вида:
( АД0 + ЬД0 = 0, \АД0 + = 0.
(13)
Здесь b = (Ъ - 3L) cos 20* + L cos 0* ± a sin 0*, d = A^g - 2L - L cos 0* sin 0* t a cot 0* ±
cos 20*
a--—. Условием устойчивости колебаний по углу АО будет положительность ко-
sin 0*
эффициента b, поскольку A > 0 всегда. Так как G на порядок больше, чем L, коэффициент L всегда положительный, а a умножается на малую величину, то неравенство b > 0 выполняется.
Для колебаний по углу Д0 условием устойчивости будет положительность коэффициента d. В силу малости 0* слагаемое Ta cot 0* будет иметь наибольший по-
рядок в выражении d. Поэтому следует предположить, что при ф* = п/2 и выполнении условия Z2 < |zi|, а также при ф* = —п/2 и выполнении условия Z2 > |zi| коэффициент d будет положительным в силу положительности слагаемого cot в* (знаки + соответствуют знакам ф* = ±п/2). В ходе численного моделирования с аналогично построенной сеткой, как для уравнения (10), и с теми же значениями фиксированных параметров ЭДТС было установлено, что при ц =1 происходит разрыв значения d, и этот коэффициент будет положительным при ^ ^ 1 для ф* = п/2 и при ^ ^ 1 для ф* = —п/2.
Были также численно проанализированы частоты колебаний связки по углам в и ф. Установлено, что оптимальной является симметричная относительно центра масс связка, поскольку такой конфигурации отвечают наибольшие частоты одновременно нутационных (по углу 9) и прецессионных (по углу ф) колебаний. Частота нутационных колебаний к\ = \/Ь/А увеличивается при увеличении зарядов коллекторов. Если |Qi| > Q2 в симметричном случае, то также наблюдается увеличение частот нутационных и прецессионных колебаний при увеличении зарядов. Полученные зависимости подтверждают тот факт, что момент сил Лоренца оказывает стабилизирующее воздействие на ЭДТС, а влияние градиентности МПЗ, приводящее к исчезновению вертикального положения равновесия связки, тем не менее не препятствует реализации устойчивого рабочего режима функционирования ЭДТС в положении, близком к вертикальному.
6. Заключение. В работе рассмотрено влияние градиентности магнитного поля Земли на положения равновесия, устойчивость и колебания электродинамической тросовой системы, содержащей заряженные коллекторы на концах троса. В ходе анализа динамических уравнений движения установлено, что градиентность геомагнитного поля приводит к исчезновению вертикального положения равновесия связки даже при условии симметричности связки относительно центра масс. Обнаруженное наклонное положение равновесия связки в плоскости орбиты для широкой области параметров соответствует малому отклонению троса от местной вертикали, что позволяет также эффективно использовать данную конфигурацию ЭДТС для выполнения задачи по уводу космического мусора с орбиты Земли. Установлена устойчивость по линейному приближению найденного положения равновесия при выполнении соответствующих соотношений для координат концевых масс. Подтверждено, что момент сил Лоренца оказывает стабилизирующее воздействие на ЭДТС и расширяет область устойчивости рабочего режима ЭДТС.
Литература
1. Белецкий В. В., Левин Е. М. Динамика космических тросовых систем. М.: Наука, 1990.
2. Родников А. В., Красильников П. С. О пространственных движениях орбитальной леерной связки // Нелинейная динам. 2017. Т. 13, №4. С. 505-518. https://doi.org/10.20537/nd1704004
3. Forward R. L. Electrodynamic drag terminator tether, Appendix K of high strength-to-weight tapered Hoytether for LEO to GEO payload transport. Final Report on NASA SBIR Phase I Contract NAS8-40690, 10 July 1996.
4. Forward R.L., Hoyt R.P., Uphoff C. Application of the Terminator Tether™ electrodynamic drag technology to the deorbit of constellation spacecraft // 34th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibition, Cleveland, OH, USA, July 13-15, 1998. Art. no. 98-3491. https://doi.org/10.2514/6.1998-3491
5. Forward R.L., Hoyt R.P. Terminator Tether™ : a spacecraft deorbit device // Journal of Spacecraft and Rockets. 2000. Vol.37. P. 187-196. https://doi.org/10.2514/2.3565
6. Tethers in Space Handbook / Eds. by M.L.Cosmo, E. C. Lorenzini. 3rd ed. Cambridge, MA, USA: Smithsonian Astrophysical Observatory, 1997.
7. Vannaroni G., Dobrowolny M., De Venuto F. Deorbiting with electrodynamic tethers: comparison between different tether configurations // Space Debris. 1999. Vol. 1. P. 159—172. https: / / doi.org/ 10.1023/A:1012570808078
8. less L., Bruno C., Ulivieri C., Ponzi U., Parisse M., Laneve G., Vannaroni G., Dobrowolny M., De Venuto F., Bertotti B., Anselmo L. Satellite de-orbiting by means of electrodynamic tethers. Part I: general concepts and requirements // Acta Astronautica. 2002. Vol. 50. Iss. 7. P. 399—406. https: //doi.org/10.1016/S0094-5765 (01)00180-1
9. less L., Bruno C., Ulivieri C., Vannaroni G. Satellite de-orbiting by means of electrodynamic tethers. Part II: system configuration and performance // Acta Astronautica. 2002. Vol. 50. Iss. 7. P. 407— 416. https://doi.org/10.1016/S0094-5765(01)00181-3
10. Ishige Y., Kawamoto S., Kibe S. Study on electrodynamic tether system for space debris removal // Acta Astronautica. 2004. Vol.55. Iss. 11. P. 917-929. https://doi.org/10.1016/j.actaastro.2004.04.015
11. Yamaigiwa Y., Hiragi E., Kishimoto T. Dynamic behavior of electrodynamic tether deorbit system on elliptical orbit and its control by Lorentz force // Aerospace Science and Technology. 2005. Vol. 9. Iss. 4. P. 366-373. https://doi.org/10.1016/j.ast.2004.09.005
12. Zhong R., Zhu Z. H., Libration dynamics and stability of electrodynamic tethers in satellite deorbit // Celestial Mechanics and Dynamical Astronomy. 2013. Vol. 116. P. 279-298. https://doi.org/10.1007/s10569-013-9489-4
13. Levin E. M. Dynamic Analysis of Space Tether Missions. San Diego, California: American Astronautical Society, 2007. (Vol. 126 of Advances in the Astronautical Sciences.)
14. Pelaez J., Lorenzini E. C., Lopez-Rebollal O., Ruiz M. A new kind of dynamic instability in electrodynamic tethers // Advances in the Astronautical Sciences. San Diego, CA: Spaceflight Mechanics, AAS Publications, 2000. Vol. 105. P. 1367-1386.
15. Тихонов А. А. Об одной конструктивной схеме электродинамического троса для расширения возможностей и повышения эффективности решения задачи спуска ИСЗ с орбиты // Сборник трудов X международной конференции «Современные методы прикладной математики, теории управления и компьютерных технологий» (ПМТУКТ-2017), 18-24 сентября 2017, Воронеж. 2017. С. 347-350.
16. Тихонов А. А. Патент RU — №2666610 на изобретение «Устройство стабилизации электродинамической тросовой системы для удаления космического мусора» по заявке №2017129790. Приоритет 22.08.2017. Дата гос. Регистрации в Гос. Реестре изобретений РФ 11.09.2018.
17. Tikhonov A. A., Shcherbakova L. F. On equilibrium positions and stabilization of electrody-namic tether system in the orbital frame // The International Scientific Conference on Mechanics "The Eighth Polyakhov's Reading", 29 January — 2 February 2018, Saint Petersburg, Russia. 2018. Vol. 1959, art. no. 020001. (AIP Conference Proceedings.) https://doi.org/doi:10.1063/1.5034626
18. Тихонов А. А., Петров К. Г. Мультипольные модели магнитного поля Земли // Космические исследования. 2002. Т. 40, №3. С. 203-212.
19. Александров А.Ю., Тихонов А. А. Одноосная стабилизация вращательного движения твердого тела при наличии возмущений с нулевыми средними значениями // Вестник С.-Петерб. ун-та. Математика. Механика. Астрономия. 2019. Т. 6(64). Вып. 2. С. 270-280. https://doi.org/10.21638/11701/spbu01.2019.209
20. Антипов К. А., Тихонов А. А. Мультипольные модели геомагнитного поля: построение N-го приближения // Геомагнетизм и аэрономия. 2013. Т. 53, №2. С. 257-267.
21. Белецкий В. В. Движение искусственного спутника относительно центра масс. М.: Наука,
1965.
22. Петров К. Г., Тихонов А. А. Момент сил Лоренца, действующих на заряженный спутник в магнитном поле Земли. Ч. 1: Напряженность магнитного поля Земли в орбитальной системе координат // Вестник С.-Петерб. ун-та. Сер. 1. Математика. Механика. Астрономия. 1999. Вып. 1. С. 92-100.
23. Ovchinnikov M. Yu., Penkov V.I., Roldugin D.S., Pichuzhkina A. V. Geomagnetic field models for satellite angular motion studies // Acta Astronautica. 2018. Vol.144. P. 171-180. https://doi.org/10.1016/j.actaastro.2017.12.026
24. Тихонов А. А. Уточнение модели «наклонный диполь» в задаче об эволюции вращательного движения заряженного тела в геомагнитном поле // Космические исследования. 2002. Т. 40, №2. С. 157-162.
25. Петров К. Г., Тихонов А. А. Момент сил Лоренца, действующих на заряженный спутник в магнитном поле Земли. Ч.2: Вычисление момента и оценки его составляющих // Вестник С.-Петерб. ун-та. Сер. 1. Математика. Механика. Астрономия. 1999. Вып. 15. С. 81-91.
26. Corsi J., less L. Stability and control of electrodynamic tether for de-orbiting applications // Acta Astronautica. 2001. Vol.48. Iss.5-12. P.491-501. https://doi.org/10.1016/S0094-5765(01)00049-2
27. Akritas A. Elements of Computer Algebra with Applications. New York: Wiley-Interscience,
1989.
Статья поступила в редакцию 4 декабря 2019 г.;
после доработки 25 декабря 2019 г.; рекомендована в печать 19 марта 2020 г.
Контактная информация:
Дериглазов Алексей Павлович — студент; alexeideriglazov@gmail.com
Тихонов Алексей Александрович — д-р физ.-мат. наук, проф.; a.tikhonov@spbu.ru
The effect of nonuniformity of the Earth's magnetic field on dynamics of electrodynamic space tether system
A. P. Deriglazov, A. A. Tikhonov
St. Petersburg State University, 7-9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation
For citation: Deriglazov A. P., Tikhonov A. A. The effect of nonuniformity of the Earth's magnetic field on dynamics of electrodynamic space tether system. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy, 2020, vol. 7(65), issue 3, pp. 539-551. https://doi.org/10.21638/spbu01.2020.316 (In Russian)
An electrodynamic space tether system in a circular equatorial orbit under conditions of nonuniformity of the Earth's magnetic field is considered. The dynamic equations of motion are analyzed to search for equilibrium modes of motion. The oblique equilibrium position of the tether in the orbit plane is found. The stability of the found equilibrium is investigated in linear approximation. The results of numerical modeling are presented. Keywords: electrodynamic tether system, stability, geomagnetic field, space debris.
References
1. Beletsky V.V., Levin E. M., Dynamics of Space Tether Systems (Nauka Publ., Moscow, 1990). (In Russian)
2. Rodnikov A. V., Krasilnikov P. S., "On spacial motions of an orbital tethered system", Nelineinaya Dinamika 13 (4), 505-518 (2017). https://doi.org/10.20537/nd1704004 (In Russian)
3. Forward R. L., Electrodynamic drag terminator tether, Appendix K of high .strength-to-weight tapered Hoytether for LEO to GEO payload transport (Final Report on NASA SBIR Phase I Contract NAS8-40690, 10 July 1996).
4. Forward R. L., Hoyt R. P., Uphoff C., "Application of the Terminator Tether™ electrodynamic drag technology to the deorbit of constellation spacecraft", 34th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibition, Cleveland, OH, USA, July 13-15, 1998, art. no. 98-3491. https://doi.org/10.2514/6.1998-3491
5. Forward R. L., Hoyt R. P., "Terminator Tether™: a spacecraft deorbit device", Journal of Spacecraft and Rockets 37, 187-196 (2000). https://doi.org/10.2514/2.3565
6. Tethers in Space Handbook (3rd ed., M.L.Cosmo, E. C. Lorenzini (eds.), Smithsonian Astro-physical Observatory, Cambridge, MA, USA, 1997).
7. Vannaroni G., Dobrowolny M., De Venuto F., "Deorbiting with electrodynamic tethers: comparison between different tether configurations", Space Debris 1, 159-172 (1999). https://doi.org/10.1023/A:1012570808078
8. Iess L., Bruno C., Ulivieri C., Ponzi U., Parisse M., Laneve G., Vannaroni G., Do-browolny M., De Venuto F., Bertotti B., Anselmo L., "Satellite de-orbiting by means of electrody-namic tethers. Part I: general concepts and requirements", Acta Astronautica 50, iss. 7, 399-406 (2002). https://doi.org/10.1016/S0094-5765(01)00180-1
9. Iess L., Bruno C., Ulivieri C., Vannaroni G., "Satellite de-orbiting by means of electrodynamic tethers. Part II: system configuration and performance", Acta Astronautica 50, iss. 7, 407-416 (2002). https://doi.org/10.1016/S0094-5765(01)00181-3
10. Ishige Y., Kawamoto S., Kibe S., "Study on electrodynamic tether system for space debris removal", Acta Astronautica 55, iss. 11, 917-929 (2004). https://doi.org/10.1016/j.actaastro.2004.04.015
11. Yamaigiwa Y., Hiragi E., Kishimoto T., "Dynamic behavior of electrodynamic tether deorbit system on elliptical orbit and its control by Lorentz force", Aerospace ¡Science and Technology 9, iss. 4, 366-373 (2005). https://doi.org/10.1016/j.ast.2004.09.005
12. Zhong R., Zhu Z.H., "Libration dynamics and stability of electrodynamic tethers in satellite deorbit", Celestial Mechanics and Dynamical Astronomy 116, 279-298 (2013). https://doi.org/10.1007/s10569-013-9489-4
13. Levin E. M., Dynamic Analysis of Space Tether Missions (American Astronautical Society, San Diego, California, 2007, vol. 126 of Advances in the Astronautical Sciences).
14. Pelaez J., Lorenzini E. C., Lopez-Rebollal O., Ruiz M., "A new kind of dynamic instability in electrodynamic tethers", Advances in the Astronautical Sciences 105, 1367-1386 (Spaceflight Mechanics, AAS Publications, San Diego, CA, 2000).
15. Tikhonov A. A., "Ob odnoy konstruktivnoy skheme elektrodinamicheskogo trosa dlya rasshireniya vozmozhnostey i povysheniya effektivnosti resheniya zadachi spuska ISZ s orbity", Sbornik trudov X mezhdunarodhoy nauchoy konferencii "Sovremennyye metody prikladnoy matematiki, teorii upravleniya i komp'yuternykh tekhnologiy", September 18-24, 2017, Vorohezh, 347-350 (2017). (In Russian)
16. Tikhonov A. A., Device of stabilization of electrodynamic cable system for removing space waste. RU-2666610-C1, B64G-001/32; B64G-001/34 (Web of Science Derwent Collection, No2018-81465E).
17. Tikhonov A. A., Shcherbakova L. F., "On equilibrium positions and stabilization of electrodynamic tether system in the orbital frame", The International Scientific Conference on Mechanics "The Eighth Polyakhov's Reading", 29 January — 2 February 2018, Saint Petersburg, Russia 1959, 020001 (AIP Conference Proceedings, 2018). https://doi.org/doi:10.1063/1.5034626
18. Tikhonov A. A., Petrov K. G., "Multipole models of the Earth's magnetic field", Cosmic Research 40 (3), 203-212 (2002). https://doi.org/10.1023/Aa015916718570
19. Aleksandrov A. Yu., Tikhonov A. A., "Uniaxial Attitude Stabilization of a Rigid Body under Conditions of Nonstationary Perturbations with Zero Mean Values", Vestnik St. Petersb. Univ. Math. 52, iss. 2, 187-193 (2019). https://doi.org/10.1134/S106345411902002X
20. Antipov K. A., Tikhonov A. A., "Multipole Models of the Geomagnetic Field: Construction of the N-th Approximation", Geomagnetism and Aeronomy 53 (2), 257-267 (2013). https://doi.org/10.1134/S0016793213020023
21. Beletsky V. V., Motion of an Artificial Satellite about its Center of Mass (Israel Program for Scientific Translation, Jerusalem, 1966).
22. Petrov K. G., Tikhonov A. A., "The moment of Lorentz forces, acting upon the charged satellite in the geomagnetic field. Part 1. The strength of the Earth's magnetic field in the orbital coordinate system", Vestnik of St. Petersburg University. Ser. 1. Mathematics. Mechanics. Astronomy, iss. 1, 92-100 (1999). (In Russian)
23. Ovchinnikov M.Yu., Penkov V. I., Roldugin D.S., Pichuzhkina A.V., "Geomagnetic field models for satellite angular motion studies", Acta Astronautica 144, 171-180 (2018). https://doi.org/10.1016/j.actaastro.2017.12.026
24. Tikhonov A. A., "Refinement of the Oblique Dipole Model in the Evolution of Rotary Motion of a Charged Body in the Geomagnetic Field", Cosmic Research 40 (2), 157-162 (2002). https://doi.org/10.1023/A:1015149420500
25. Petrov K. G., Tikhonov A. A., "The moment of Lorentz forces, acting upon the charged satellite in the geomagnetic field. Part 2. The determination of the moment and estimations of its components", Vestnik of St. Petersburg University. Ser. 1. Mathematics. Mechanics. Astronomy, iss. 15, 81-91 (1999). (In Russian)
26. Corsi J., Iess L., "Stability and control of electrodynamic tether for de-orbiting applications", Acta Astronautica 48 (5-12), 491-501 (2001). https://doi.org/10.1016/S0094-5765(01)00049-2
27. Akritas A., Elements of Computer Algebra with Applications (Wiley-Interscience, New York, 1989).
Received: December 4, 2019 Revised: December 25, 2019 Accepted: March 19, 2020
A u t h o r s' i n fo r m a t i o n:
Alexey P. Deriglazov — alexeideriglazov@gmail.com
Alexey A. Tikhonov —a.tikhonov@spbu.ru