УПРАВЛЕНИЕ РОБОТОМ-МАНИПУЛЯТОРОМ ПРИ СМЕНЕ ТОПЛИВНЫХ БАКОВ ДВИГАТЕЛЬНОЙ УСТАНОВКИ ГЕОСТАЦИОНАРНОГО СПУТНИКА Текст научной статьи по специальности «Механика и машиностроение»

УДК 629.78 : 681.51

УПРАВЛЕНИЕ РОБОТОМ-МАНИПУЛЯТОРОМ ПРИ СМЕНЕ ТОПЛИВНЫХ БАКОВ ДВИГАТЕЛЬНОЙ УСТАНОВКИ ГЕОСТАЦИОНАРНОГО СПУТНИКА

© 2022 С .А. Бутырин12, Е.И. Сомов1-2, С.Е. Сомов1-2, Т.Е. Сомова2

1 Самарский федеральный исследовательский центр Российской академии наук, г. Самара, Россия 2 Самарский государственный технический университет, г. Самара, Россия

Статья поступила в редакцию 15.01.2022

Разработаны методы управления космическим роботом-манипулятором при смене топливных баков двигательной установки при одновременной стабилизации углового движения механической связки робота с геостационарным спутником. Определяются требования к электромеханическим приводам робота-манипулятора и представляются результаты компьютерной имитации динамических процессов.

Ключевые слова: геостационарный спутник связи, смена топливных баков, управление космическим роботом.

Б01: 10.37313/1990-5378-2022-24-1-96-104

Работа поддержана РФФИ, грант 20-08-00779.

ВВЕДЕНИЕ

В связи с возросшими требованиями к глобальности, скорости и объемам информационного обмена в настоящее время многие страны приступили к реализации низкоорбитальных спутниковых систем связи, состоящих из тысяч мини-спутников на орбитах высотой от 500 до 1600 км. Достойный ответ традиционной космиче -ской связи на основе геостационарных платформ противостоит этим вызовам. Такие перспективные платформы собираются на геостационарной орбите (ГСО) космическими роботами-манипуляторами (КРМ) из сменных и пополняемых компонентов, а затем регулярно обслуживаются КРМ в течение нескольких десятилетий.

Традиционные геостационарные спутники связи (ГСС) имеют гарантированный срок службы до 15 лет. При продлении срока службы

Бутырин Сергей Анфимович, старший научный сотрудник отдела «Динамики и управления движением» СамНЦ РАН; начальник лаборатории «Моделирования систем управления» НИИ Проблем надежности механических систем СамГТУ. E-mail butyrinsa@mail.ru Сомов Евгений Иванович, ведущий научный сотрудник отдела «Динамики и управления движением» СамНЦ РАН; начальник отдела «Навигации, наведения и управления движением» НИИ Проблем надежности механических систем СамГТУ. E-mail e_somov@mail.ru Сомов Сергей Евгеньевич, научный сотрудник отдела «Динамики и управления движением» СамНЦ РАН; научный сотрудник отдела «Навигации, наведения и управления движением» НИИ Проблем надежности механических систем СамГТУ. E-mail s_somov@mail.ru Сомова Татьяна Евгеньевна, научный сотрудник отдела «Навигации, наведения и управления движением» НИИ Проблем надежности механических систем СамГТУ. E-mail te_somova@mail.ru

ГСС важнейшей является проблема дозаправки топливом их электрореактивных двигательных установок (ЭДУ). В системе управления движением (СУД) КРМ применяются исполнительные органы в виде кластера гиродинов (ГД) - двухстепенных силовых гироскопов и ЭДУ на основе как плазменных, так и каталитических электрореактивных двигателей (ЭРД) [1,2]. В штатном режиме СУД измерение координат движения космического аппарата (КА) выполняется бесплатформенной инерциальной навигационной системой (БИНС) с коррекцией сигналами от навигационных спутников GPS/ ГЛОНАСС и звездных датчиков.

В статье рассматриваются вопросы управления КРМ при дозаправке ЭДУ на основе смены топливных баков. Предлагается стратегия и разрабатываются методы управления КРМ в этом режиме при одновременной стабилизации углового движения механической связки КРМ с ГСС. Определяются требования к электромеханическим приводам робота-манипулятора и представляются результаты компьютерной имитации динамических процессов.

ОПИСАНИЕ ЗАДАЧИ И СТРАТЕГИЯ УПРАВЛЕНИЯ

Рассматривается задача управления КРМ при смене топливных баков двигательной установки ГСС. Предполагается, что КРМ и ГСС жестко состыкованы, а механические манипуляции с баками происходят в автоматическом режиме. Концевое звено манипулятора содержит специальный механизм - ключ, который при замыкании жестко удерживает очередной топливный бак.

Целями исследования являются: анализ стратегии управления - последовательности действий манипулятора при захвате/освобождении баков и их перемещениях в заданные положения; выбор кинематической схемы манипулятора, числа его звеньев, их геометрических и инерционных характеристик, а также параметров электромеханических приводов; расчет рабочей зоны манипулятора; оценка динамических характеристик СУД при угловой стабилизации связки КРМ и ГСС в процессе смены баков. Компоновка жесткой сцепки КРМ с ГСС, для компактности без отображения панелей солнечных батарей (СБ), представлена на рис. 1, где все размеры указаны в метрах.

На рисунке 1 приведены используемые формы и размеры компонентов механической связки: корпус КРМ - цилиндр диаметром 1 м и длиной 2.5 м с установленными на нем манипулятором и контейнерами с баками, частично вы-

ступающими за его обвод; корпус ГСС - цилиндр диаметром 1.1 м и длиной 2.7 м. Было принято, что каждый сменный топливный бак имеет размеры 0.4 х 0.4 х 0.5 м, объем 80 литров, масса полного бака равна 90 кг, а пустого бака - 10 кг. При выбранных геометрических параметрах на каждой из двух продольных сторон корпуса КРМ допускается размещение до пяти контейнеров. Для начала последовательности смены топливных баков один из контейнеров в корпусе КРМ должен быть пустым, см. например, контейнер а на рис. 1.

На рисунке 2 представлена кинематическая схема бортового манипулятора антропоморфной структуры с 4 звеньями (/ = 1,2,...4 = 1 + 4) и 6 степенями свободы, которым соответствуют угловые координаты qs, 5 = 1 + 6. Здесь указаны концы звеньев манипулятора в точках А, В, С и концевая точка D ключа замкового механиз-

Рис. 1. Схема управления манипулятором и его рабочая зона при смене топливных баков

Рис 2. Кинематическая схема манипулятора

+

ма. В точках А и С применяются двухстепенные шарниры, которые однозначно формализуются двумя стандартными одноосными шарнирами при значениях длины фиктивных звеньев /01 = !03 = 0, см. рис. 2. Последний шарнир обеспечивает вращение ключа открытия/закрытия замка вокруг продольной оси бака, а остальные - совмещение точки С манипулятора с заданной в рабочей зоне при произвольной ориентации последнего звена.

Манипулятор закреплен на корпусе КРМ

стойкой (master) в точке Om системы координат Omxmymzm , которая фиксирована в системе координат Or xryrzr, связанной с корпусом КРМ

(robot) в его центре масс Or, см. рис. 1. Захват бака достигается перемещением ключа по нормали к плоскости крышки бака (см. рис. 1) с последующим его поворотом на угол к /2.

Применяемая стратегия смены топливных баков ЭДУ использует такие положения:

1) пустые баки поочередно перемещаются манипулятором из контейнеров ГСС в контейнеры КРМ, а заполненные топливом баки - наоборот, из контейнеров КРМ в контейнеры ГСС;

2) установка баков в контейнеры завершается их автоматической механической фиксацией;

3) включение /отключение системы подачи топлива от каждого бака к ЭДУ осуществляется роботизированным мехатронным модулем, встроенным в каждый контейнер ГСС.

МАТЕМАТИЧЕСКИЕ МОДЕЛИ И ПОСТАНОВКА ЗАДАЧИ

При описании движения механической связки КРМ с ГСС применяются (i) экваториальная инерциальная системы координат (ИСК)

I ф с началом в центре Земли Оф , (ii) связанная система координат (ССК) B (Oxyz) с началом в полюсе O, которая совпадает с системой координат Orxryrzr, КРМ (iii) орбитальная система координат (ОСК) O (Ox0 y° z0) с началом в полюсе O и ортами ro, т°, no, которая имеет следующие направления осей и ортов: ось Ox0 направлена по радиали r°, ось Oz0 - по нормали n° к плоскости орбиты, а ось Oy0 - по транс-

версали т° и направлена в сторону орбитального движения, см. рис 1, где представлена также

система координат Otxtytzt, связанная с корпусом ГСС (цель, target) в его центре масс Ot.

Далее используются общепринятые обозначения col(-) = {•}, line(-) = [■], (•,•), [*], (•)*, [ax] и °для векторов, матриц и кватернионов, матриц [a] i стандартного элементарного поворота вокруг i -ой оси на угол a, i = 1 ^ 3, а

также Ca = cos a, Sa = sin a .

Рис. 3. Схема СГК и область вариации его КМ

Для управления ориентацией связки КРМ с ГСС применяется силовой гироскопический кластер (СГК) на основе четырех ГД по схеме 2-БРЕ с ортами векторов кинетических моментов (КМ) h (Р ), р = 1 ^ 4, см. рис. 3, где приведена также область вариации нормированного вектора КМ ^Р) = ^ (Р ) кластера со столбцом Р = (Рр} и ее проекции на плоскости

базиса Oxgуг^. Применяемый явный аналитический закон настройки СГК (распределения трехмерного вектора его управляющего момента М8 = {М между четырьмя ГД) позволяет исключить избыточность данного кластера с вектором кинетического момента Н = ^И(Р), где Н - одинаковое для всех четырех ГД постоянное значение модуля собственного КМ.

Ориентация ССК В в ИСК I ф определяется кватернионом Л = (Х0,Х), где X = (X}, вектором модифицированных параметров Родрига (МПР) а = (а} = е £ (Ф /4) с ортом е оси Эйлера и углом Ф собственного поворота.

В ИСК I ф кинематические уравнения для вектора го расположения КРМ, кватерниона Л и вектора МПР а имеют вид

г0 = г0 + ЮХ го;

1 ^ 2ч 1 1

су = -j(l -о2)ю+т охю+-у (о, ю)о,

Л = Л о и/2;

1

2'

где вектор Ю представляет абсолютную угловую скорость корпуса КРМ и используется обозначение (•)' локальной производной по времени. Кватернион Л0 ориентации базиса О относительно базиса I9 определяется

уравнением Л° —Л° о ю°/2, углами рыскания

ф1 = у, крена ф2 = ф и тангажа ф3 =0 в последовательности 132, а также матрицей направляющих косинусов координатного перехода от ОСК к ССК С0 = [ф2]2 [ф3]3 [фД- В свою очередь, погрешность ориентации базиса В в орбитальном базисе О определяется кватернионом Е = Л° ° Л = (е0, е), где е0 = Сфе/2, е = 5фе/2ее с ортом оси Эйлера ее и углом ошибки Фе, матрицей Се = 13 - 2[ex]Qe, где Qe = 13е0 + [ех], вектором угловой погрешности 8ф = {8фг.} — {2в0вЛ} и вектором МПР

0 — {а®} — ее tg(Фe /4). При этом вектор 5ш погрешности угловой скорости определяется

как 5ш = ш - Сеш° (¿) .

При моделировании движения связки КРМ с ГСС применяется векторная форма классических уравнений Эйлера-Лагранжа [3].

В ССК Охуг с полюсом О векторы р.,

1 = 1 ^ 4 определяют положения центров масс

ег звеньев манипулятора с массами т. и собственными тензорами инерции Jгс, а векторы рг = 0 и р 1 - положения центров масс Ог и О( робота (индекс г) и цели (индекс 1) с массами и собственными тензорами инерции тг, IС

и т4,1С соответственно. Положение центра

масс С связки КРМ и ГСС (робот, манипулятор с 4 звеньями и цель, см. рис. 1) суммарной массы т = тг + Ет. + т4 определяется вектором Рс = К, ус, 2С} по соотношению

L = т Рс = тг Рг +Ет, рг + т1 р1,

где введен вектор статического момента L. Тензор инерции I механической системы в полюсе О вычисляется по соотношению

Iл = 10 + Е10 +10,

где при единичном тензоре Е

10 = I г; 1° = I с + т.(Е р1 рг-ргр1);

з; — зч+т( (Б р; Рг-Рг р|).

При векторе V 0 скорости полюса О поступательное движение связки КРМ с ГСС описывается векторным уравнением

тV0 - Ь ХЮ — ЮХ(Ь Хй) -Х, (тг (2юхр' + р*)) + Ре + V®, (2)

где V0 — у'0 + ЮХ\0; р' — Е,((Эр,./Эд,)д,); рГ—Х, ((Эр, / д д,) д, (д2 р, / д дк д д,)),

векторы Ре и Р8Г представляют силы тяги ЭДУ, центрированной в полюсе О , и гравитации.

Через ш. обозначим вектор угловой скорости -го звена манипулятора в ССК. Производная этого вектора по времени имеет вид ш . = ш ■ + ш х ш.. В ССК Охуг с полюсом О вращательное движение связки указанных твердых тел описывается векторным уравнением

Ь Х V 0 + 1Ю —-Ь Х (ЮХ уо)

— ■

ЮХ (1Ю+ Н) - И + М8 + Ме + М8Г, (3)

где

R = Е. (I.со. + ш х Iш.. + ш.. х I(ш + ш..) +т р.х (ш х (ш х р.)+2ш х р;+рг"),

векторы М8 = -Н; и Ме = М представляют управляющие моменты СГК с вектором КМ Н и ЭДУ на основе 8 каталитических ЭРД, а вектор М8Г - гравитационный момент.

Векторные уравнения Эйлера (2), (3) дополняются стандартными уравнениями Лагранжа по степеням подвижности манипулятора, где в правых частях наряду с обобщенными силами Qs, представляющими моменты электромеханических приводов, имеются механические моменты [3], которые обусловлены пространственным движением связки КРМ с ГСС.

При цифровом управлении и8(^) — {и8^ )} СГК с периодом Ти, где для к е N = [0,1,2,...)

компоненты (/)— V/ е [1к, , /к+1 = 4 + ти формируют вектор управляющего момента СГК

М£(0 = -\Аь(0(0 и£(0 ; р(0 = и8(0, (4)

где матрица Якоби Ай (Р) — ЭИ(Р)/ЭР.

В процессе угловой стабилизации связки КРМ с ГСС в ОСК при законе наведения

Л°(0,Ю°(0, (00(0 = 80(0 после дискретной фильтрации измеренных значений вектора

углового рассогласования 8^ = -5ф;, I е N с периодом Тд = Ти /4 формируются значения вектора , к е N, которые применяются в алгоритме управления СГК с периодом Ти в виде

£+1 — Ва§: + СаеГ; Ч — Ка(§; + Рае;Т);

М£ — Ю; Х С; + 3(С;е: + [С;Ю:Х]Ю; + Ш; ). (5)

Здесь вектор Gк = Jшk + Нк и используются постоянные диагональные матрицы Ка, Ва , Са и Ра. Далее вектор М£ с помощью явного

закона распределения команд между четырьмя ГД [4] «пересчитывается» в столбец и | — {и^ } командных сигналов управления гиродинами, которые фиксируются на полуинтервалах цифрового управления СГК с периодом Ти при формировании его управляющего момента М£ (^) по соотношениям (4).

В статье решаются следующие задачи:

(1) обоснование кинематической схемы, параметров манипулятора и его рабочей зоны;

(и) синтез законов наведения схвата манипулятора;

(III) оценка потребных характеристик электромеханических приводов манипулятора;

(IV) компьютерный анализ динамики СУД при смене топливных баков ЭДУ.

ПАРАМЕТРЫ МАНИПУЛЯТОРА И ЕГО РАБОЧАЯ ЗОНА

Будем считать, что начало От системы координат (СК) манипулятора Отхтутгт определяется вектором гт = {-0.2,-1,-0.6} м в СК Ог хг уг гг корпуса КРМ, см. рис. 1, где также показаны контейнеры а (пустой) и Ь с центрами крышек {0.55, -1.75, - 0.1} м и {0.55, -1.35, - 0.1}м, соответственно. Здесь центры крышек баков лежат в одной плоскости, т.е. рабочая зона манипулятора должна располагается над плоскостью крышек (zm > 0) и в ней не должно быть мешающих элементов конструкции КРМ и ГСС.

Положение центра масс ГСС относительно полюса О определяется вектором р( = {0, 2.8, 0} м. Пустой (заменяемый) бак f в системе координат ГСС О4 х1 у1 zt имеет координаты центра крышки {0,1,0} м. Расстояние от контейнера а КРМ до бака в контейнере f ГСС составляет 5.6 м. Из простых геометрических расчетов следует, что для суммы длин первых двух звеньев манипулятора должно выполняться условие + 12 > 3.3 м. Этим условиям удовлетворяет манипулятор с кинематической схемой на рис. 2 при длинах звеньев

1Х = 2 м, 12 = 1.5 м, 13 = 0.5 м, когда 14 = 0 (без бака) и /4 — 0.25 м (с баком), и с рабочей зоной в виде полусферы с радиусом 3 м, выделенной на рис. 1 голубым цветом, при следующих ограничениях на угловые координаты:

- п < qx < п, 0 < qг < п;

-п< <п, 1 — 3^6. (6)

С каждым .-ым (. = 1 ^ 4) звеном связывается СК с началом в точке сопряжения -ой пары звеньев, см. рис. 2. Кинематическая модель манипулятора обязательно включает формулы решения прямой и обратной задачи кинематики [5].

Прямая задача состоит в расчете векторов декартовых координат r(q), линейных скоростей v(q, q) и ускорений w(q, q, q) заданных точек в кинематической цепи манипулятора, а также кватернионов либо матриц ориентации, в системе координат манипулятора Omxmymzm при заданных значениях векторов угловых координат q = {qs}, s = 1 ^ 6, скоростей q и ускорений q звеньев. Для решения прямой задачи кинематики используется обратный рекуррентный переход из i -ой в (i -1) -ую СК, начиная с последнего звена до стойки при i = 1 [6].

Обратная задача кинематики манипулятора состоит в определении вектора угловых координат q = {qs} по заданным декартовым координатам вектора rD = (xD,yD,zD} точки D схвата и ориентации СК последнего звена манипулятора относительно СК Omxmymzm стойки (master) при ограничениях (6). В общем случае решение обратной задачи неоднозначно. Для схемы на рис. 2 единственное решение этой задачи получается однозначно при наличии возможности перемещения звена 3 с точкой D в отрицательном направлении оси Om Zm, см. рис. 1.

ЗАКОНЫ НАВЕДЕНИЯ СХВАТА И УПРАВЛЕНИЯ ПРИВОДАМИ

Синтез законов программного наведения схвата манипулятора при переходах между целевыми состояниями с краевым условиям по каждой из угловых координат qs включает два аспекта: (i) расчет значений времени начала tx (initial), завершения tf (final) и длительности Tm = tf - t; выполнения каждого этапа и (ii) формирование временных зависимостей углов qs (t), скоростей qs (t) и ускорений qs (t) при ограничениях на углы (6), а также на модули скоростей и ускорений

| qs (t) |< q*, | qs(t) |< q Vt e Tm = [t;, tf ] (7)

при заданных значениях параметров q и q . Здесь каждый из последовательных этапов перемещения включает стандартные участки разгона, движения с постоянной скоростью и торможения.

При цифровом управлении редукторными электроприводами по степеням подвижности манипулятора применятся адаптивно-робастные алгоритмы с эталонной моделью наведения [7].

ОЦЕНКА МОМЕНТНЫХ СВОЙСТВ ПРИВОДОВ МАНИПУЛЯТОРА

Для простейшей оценки потребных момент-ных характеристик электромеханических приводов вычисляются обобщенные силы - кру-

тящие моменты на выходных валах приводов манипулятора, когда их нагрузка представлена только обобщенными силами инерции при движении звеньев, а разнообразные моменты сил сопротивления такому движению не учитываются. В этом случае обобщенные силы Q = Q(/) вычисляются по соотношению

Q(t) = А (Ч (г)) ) + Ъ^ (г), <1 (г)), (8)

где А^) = ) - симметричная матрица инерции и ), 4(г) - вектор центробеж-

ных и кориолиссовых сил, приведенных к выходным валам приводов в шарнирах манипулятора.

АНАЛИЗ ДИНАМИКИ СИСТЕМЫ УПРАВЛЕНИЯ ДВИЖЕНИЕМ

нической системы твердых тел:

тг = 2800 кг; I° = 103 ^{5,3,4}кгм2;

т( = 3500 кг; I ° = 103 diag{20,18,13}кгм2; массы первых трёх звеньев манипулятора равны 50, 35 и 15 кг; масса 4-го звена составляет 10 кг (с пустым баком) или 90 кг (с заполненным баком), а сцепка КРМ и ГСС со всеми топливными баками имеет массу т = 6800 кг.

Было принято, что период цифрового управления СГК Ти = 4 с, каждый ГД имеет собственный КМ = 30 Нмс, а ограничения для приводов манипулятора таковы: ¿д*= 0.025 с"1; д* = 0.0125с-2.

При смене баков манипулятор выполняет типовые действия, поэтому далее кратко приводятся результаты компьютерной имитации замены только первого топливного бака, когда

1) схват переводится из транспортного по-

ложения г

= (хв, ув, гв} = (0,0,0}

сначала к

Рис. 4. Траектория схвата манипулятора при замене первого топливного бака ЭДУ

Анализ динамики СУД при смене топливных баков ЭДУ выполнен на основе методов и средств компьютерной имитации для следующих значений инерционных параметров меха-

контейнеру ГСС с пустым баком f (см. рис. 4), который затем перемещается в пустой контейнер а КРМ;

2) схват переходит к контейнеру Ь КРМ и далее перемещает полный бак Ь в контейнер f ГСС.

Декартовые координаты целевых точек рабочей зоны манипулятора в его системе координат От хт ут гт и последовательность действий во времени представлены в табл. 1.

На рисунках 5 и 6 представлены изменения угловых координат звеньев манипулятора и оценки соответствующих им потребных моментов на выходных валах приводов при замене первого топливного бака. При использовании редукторов с передаточным отношением

Таблица 1. Состояния манипулятора, координаты схвата, время выполнения этапов и их длительность

№ Последовательность состояний манипулятора и положения хвата Координаты схвата при завершении этапа, м Отметки времени и длительность, с

х0 Уо ^ к к Т ^ т

1 Развертывание, переход к баку £ и зависание 0.25 2.8 0 0 102 102

2 Вставка ключа в пустой бак £ ГСС 0.25 2.8 - 0.1 102 108 6

3 Закрытие замка и захват бака £ 0.25 2.8 - 0.1 108 175 67

4 Вытягивание бака £ из контейнера £ ГСС 0.25 2.8 0.25 175 186 11

5 Бак £ выведен из контейнера £ ГСС 0.25 2.8 0.5 186 195 9

6 Переход к баку а КРМ и зависание 0.25 -1.75 0.5 195 309 114

7 Вставка бака £в контейнер а 0.55 -1.75 - 0.1 309 333 24

8 Открытие замка бака £ в контейнере а КРМ 0.55 -1.75 - 0.1 333 384 51

9 Зависание над контейнером а КРМ 0.55 -1.75 0 384 391 7

10 Переход к баку Ь КРМ и зависание 0.55 -1.35 0 391 405 14

11 Вставка ключа и захват бака Ь КРМ 0.55 -1.35 - 0.1 405 476 71

12 Вытягивание бака Ь из контейнера Ь КРМ 0.55 -1.35 0.5 476 503 27

13 Переход к контейнеру £ ГСС и зависание 0.25 2.8 0 503 614 111

14 Вставка бака Ь в контейнер £ ГСС 0.25 2.8 - 0.1 614 725 111

Рис. 5. Изменение угловых координат манипулятора при замене первого топливного бака

а

100

200

300

400

500

600

700

1,5

Рис. 6. Оценки потребных моментов на выходных валах приводов манипулятора

хЮ

Рис. 7. Угловое отклонение ССК относительно ОСК при смене топливного бака

Рис. 8. Изменение вектора угловой скорости связки КРМ и ГСС при смене бака

Рис. 9. Управляющий момент силового гироскопического кластера при смене бака

600

1, Б

Рис. 10. Изменения угловых скоростей гиродинов при смене топливного бака

Рис. 11. Изменение вектора положения центра масс связки КРМ и ГСС при смене бака

Рис. 12. Изменение вектора скорости центра масс связки КРМ и ГСС при смене бака

спутника // Известия Самарского научного цен-

1:400 потребные моменты с учетом возможных потерь не превысят 10 Нм, что позволяет использовать шаговые двигатели с номинальным моментом 0.025 Нм и максимальной скоростью 100 об/мин.

Изменения кинематических параметров ориентации связки КРМ и ГСС относительно ОСК при смене бака приведены на рис. 7 и 8, а соответствующие изменения вектора управляющего момента СГК и угловых скоростей гиро-динов - на рис. 9 и 10.

Наконец, на рисунках 11 и 12 представлены изменения векторов положения и скорости поступательного движения центра масс связки КРМ и ГСС при смене бака.

ЗАКЛЮЧЕНИЕ

Кратко представлены стратегия и методы управления космическим роботом-манипулятором при смене топливных баков электрореактивной двигательной установки геостационарного спутника связи. Для конкретных исходных данных установлены требования к электромеханическим приводам манипулятора и представлены результаты компьютерной имитации динамических процессов при смене топливных баков с одновременной угловой стабилизацией механической связки робота и геостационарного спутника в орбитальной системе координат.

СПИСОК ЛИТЕРАТУРЫ

1. Сомов Е.И., Бутырин С.А., Сомов С.Е. Управление космическим роботом-манипулятором при встрече и механическом захвате пассивного

тра РАН. 2018. Том 20, № 6. С. 202-209.

2. Somov Ye., Butyrin S., Somov S., Somova T. Control of robot-manipulator during its preparation and capture of a passive satellite // Mathematics in Engineering, Science and Aerospace. 2019. Vol. 10, no. 3, pp. 421-432.

3. Лурье А.И. Аналитическая механика. М.: Физмат-лит. 1961. 824 с.

4. Сомов Е.И. Анализ сингулярных состояний и синтез явных законов настройки гирокомплексов кратных схем // Гироскопия и навигация. 2013. № 1(80). С. 134-148.

5. Фу К., Гонсалес Р., Ли К. Робототехника. М.: Мир, 1989. 624 с.

6. Елисеев С.В., Бутырин С.А., Засядко А.А. Функциональные модули пакета прикладных программ для задач проектирования манипуляционных роботов // Пакеты прикладных программ. Методы и разработки. Новосибирск: Наука, 1981. С. 206-216.

7. Крутько П.Д. Управление исполнительными системами роботов. М.: Наука. 1991. 336 с.

CONTROL OF THE ROBOT-MANIPULATOR WHEN CHANGING FUEL TANKS OF THE GEOSTATIONARY SATELLITE PROPULSION UNIT

© 2022 S.A. Butyrin12, Ye.I. Somov1-2, S.Ye. Somov12, T.Ye. Somova2

1 Samara Federal Research Scientific Center, Russian Academy of Sciences, Samara, Russia 2 Samara State Technical University, Samara, Russia

Abstract. Methods have been developed for control of a space robot-manipulator when changing fuel tanks of the propulsion unit while simultaneously stabilizing the angular motion of the robot mechanical mated with a geostationary satellite. The requirements for the electromechanical drives of the robot manipulator are determined and the results on computer simulation of dynamic processes are presented.

Key words: geostationary communications satellite, fuel tank change, space robot control DOI: 10.37313/1990-5378-2022-24-1-96-104

Sergey Butyrin, Senior Researcher of Department "Dynamics and Motion Control", Samara Federal Research Centre, Russian Academy of Sciences; Head of Laboratory for "Modeling of Control Systems", Research Institute for Problems of Mechanical Systems Reliability, Samara State Technical University. E-mail butyrinsa@mail.ru Yevgeny Somov, Leading Researcher of Department "Dynamics and Motion Control", Samara Federal Research Centre, Russian Academy of Sciences; Head of Department for "Navigation, Guidance, and Motion Control", Research Institute for Problems of Mechanical Systems Reliability, Samara State Technical University. E-mail e_somov@mail.ru Sergey Somov, Researcher of Department "Dynamics and Motion Control", Samara Federal Research Centre, Russian Academy of Sciences; Researcher of Department "Navigation, Guidance, and Motion control", Research Institute for Problems of Mechanical Systems Reliability, Samara State Technical University. E-mail s_somov@mail.ru Tatyana Somova, researcher of department "Navigation, Guidance, and Motion Control", Research Institute for Problems of Mechanical Systems Reliability, Samara State Technical University. E-mail te_somova@mail.ru