Переориентация космического аппарата дистанционного зондирования Земли с помощью роторов

Полный текст

Введение

С каждым днем появляется все больше потребностей в дистанционном зондировании Земли (ДЗЗ). Дистанционное зондирование – это процесс наблюдения и изучения нашей планеты, при котором используются специальные инструменты и технологии, используемые на спутниках или самолетах. Существуют и другие задачи, в которых требуется нацеливание космического аппарата (КА) на определенный объект.

Чаще всего при переориентации КА для нацеливания антенны или камеры используются двигатели малой тяги, которые подразделяются на несколько категорий: на сжатом газе, на монотопливе, на двухкомпонентном топливе; специальные установки и двигатели на закиси азота, аммиаке, газообразном водороде и кислороде; электроракетные системы на базе импульсных плазменных, ионных или стационарных плазменных двигателей [1].

Все вышеперечисленные инструменты для переориентации КА подразумевают расход рабочего тела, запас которого на борту КА ограничен. Предлагается использовать двигатели-маховики на основе управляемого моментного бесконтактного двигателя постоянного тока [2]. Данный способ является более эффективным из-за отсутствия необходимости в расходуемом рабочем теле, и, как следствие, уменьшения рабочей массы КА. Подобные КА называют спутниками-гиростатами.

Исследованию движения гиростатов посвящено достаточно много работ отечественных и зарубежных ученых, например, Румянцева В.В. [3; 4], Кошлякова В.Н. [5], Нейштадта А.И. и Пивоварова М.Л. [6], Виттенбурга Й. [7], Асланова В.С. и Дорошина А.В. [8; 9] и многих других. В отмеченных работах, в основном, исследуется динамика и устойчивость движения гиростатов, но мало внимания уделяется определению законов движения и управляющих воздействий, обеспечивающих данные движения.

В настоящей работе предпринимается попытка решения конкретной прикладной задачи: определения закона ориентации КА ДЗЗ с частичной закруткой для наблюдения за объектом, а также управляющего внутреннего момента, обеспечивающего полученный закон.

Математическая модель

Предположим, что спутник-гиростат представляет собой динамически несимметричное несущее тело (тело 4) с тремя одинаковыми динамически симметричными роторами (тела 1, 2, 3), оси вращения которых совпадают с главными осями инерции всей системы (рисунок 1). На рисунке 1 роторы изображены схематично (вне несущего тела) для наглядности и удобства обозначения координатных осей и других параметров: на самом деле роторы находятся внутри КА.

 

Рисунок 1. Схема спутника-гиростата и используемые системы координат

 

Динамические уравнения движения спутника-гиростата строятся на основании теоремы об изменении кинетического момента механической системы [3]:

$\begin{array}{c}Ap+(C-B)qr+I{\sigma }_1+I(q{\sigma }_3-r{\sigma }_2)=M_x; \\ Bq+(A-C)pr+I{\sigma }_2+I(r{\sigma }_1-p{\sigma }_3)=M_y; \\ Cr+(B-A)pq+I{\sigma }_3+I(p{\sigma }_2-q{\sigma }_1)=M_z.\end{array}$ (1)

где $A=\sum _{i=1}^4{A}_i$, $B=\sum _{i=1}^4{B}_i$, $C=\sum _{i=1}^4{C}_i$; A_i, B_i, C_i $(i=\overline{1,3})$ – главные моменты инерции роторов, вычисленные в своих связанных системах координат $O{x}_i{y}_i{z}_i(i=\overline{1,3})$, а моменты инерции тела-носителя в системе координат $Oxyz-A_4, B_4, C_4$; $p, q, r$ – проекции угловой скорости несущего тела на оси системы координат Oxyz; $I=A_1=B_2=C_3$ – осевые моменты инерции роторов. Правые части уравнений (1) представляют собой проекции момента внешних сил.

Кинематические уравнения для углов Эйлера и углов относительного закручивания:

$\begin{array}{c}\phi =r-ctg\theta (p \sin \phi +q \cos \phi );\\ \theta =p \cos \phi -q \sin \phi ;\\ \psi =\frac{1}{\sin \theta }(p \sin \phi +q \cos \phi );\\ \alpha ={\sigma }_1; \beta ={\sigma }_2; \gamma ={\sigma }_3.\end{array}$ (2)

где $\phi , \theta , \psi$ – эйлеровы углы, которые характеризуют положение несущего тела относительно кениговой системы OXYZ; $\alpha , \beta , \gamma$ – углы относительного закручивания роторов; ${\sigma }_1, {\sigma }_2 , {\sigma }_3$ – относительные угловые скорости роторов.

Уравнения относительного движения роторов:

$\begin{array}{c}{A}_1(p+{\sigma }_1)=M_{1x}+M_{\alpha };\\ B_2(q+{\sigma }_2)=M_{2y}+M_{\beta };\\ C_3(r+{\sigma }_3)=M_{3z}+M_y.\end{array}$ (3)

где $M_{\alpha }, M_{\beta }, M_{\gamma }$ – моменты внутреннего взаимодействия тел вдоль оси вращения; $M_{1x}, M_{2y}, M_{3z}$ – моменты внешних сил, воздействующие на роторы.

Компоненты векторов угловых скоростей роторов $\overrightarrow{{\varpi }_1}(p_1, q_1, r_1)$, $\overrightarrow{{\varpi }_2}(p_2, q_2, r_2)$, $\overrightarrow{{\varpi }_3}(p_3, q_3, r_3)$, выраженные через компоненты угловой скорости несущего тела, имеют вид:

$\left\{\begin{array}{l}{p}_1=p+{\sigma }_1,\\ q_1=q \cos \alpha +r \sin \alpha ,\\ r_1=r \cos \alpha -q \sin \alpha .\end{array}\right.$ (4)

$\left\{\begin{array}{l}{p}_2=p \cos \beta +r \sin \beta ,\\ q_2=q+{\sigma }_2,\\ r_2=r \cos \beta +p \sin \beta .\end{array}\right.$ (5)

$\left\{\begin{array}{l}{p}_3=p \cos y+q \sin y,\\ q_3=q \cos y-p \sin y,\\ r_3=r +{\sigma }_3.\end{array}\right.$ (6)

Таким образом, уравнения (1) – представляют собой замкнутую математическую модель движения спутника-гиростата относительно центра масс. Приведенные уравнения при добавлении начальных условий позволяют определить все параметры движения системы в зависимости от времени.

Постановка задачи и преобразование математической модели

На данном этапе рассматривается так называемая плоская задача, для которой принимается ряд основных допущений:

1. Спутник движется прямолинейно с постоянной скоростью центра масс ν, находясь в одной вертикальной плоскости с наблюдаемым объектом, на высоте H, в начальный момент на расстоянии L (рисунок 2).
2. Для нацеливания спутника-гиростата используется только один маховик, ось вращения которого совпадает с осью Ox.
3. Значение угла $\theta$ в начальным момент принято равным 30°, что соответствует нацеливанию на объект наблюдения (остальные геометрические размеры соответствуют данному начальному углу).
4. Так как расстояние, на котором происходит съемка объекта, значительно меньше длины окружности земли ($\approx \frac{1}{46}$ длины экватора), то сферичностью земной поверхности пренебрегаем.
5. Внешние силы и моменты на систему не действуют.

 

Рисунок 2. Схема движения спутника-гиростата

 

С учетом перечисленных допущений, движение спутника-гиростата описывается следующими уравнениями, полученными из системы (1)–(6):

$\left\{\begin{array}{l}Ap+I{\sigma }_1=0;\\ \theta =p;\\ I(p+{\sigma }_1)=M_{\alpha };\\ \alpha ={\sigma }_1.\end{array}\right.$ (7)

Таким образом, ставится задача определения внутреннего управляющего момента , который обеспечивает нацеливание КА ДЗЗ на объект наблюдения в течение некоторого интервала времени.

Определение управляющего момента

В течении всего времени t активного участка полета (съемки) угол нутации может быть определен по формуле (рисунок 2):

$\theta =arccot\frac{L-vt}{H}$. (8)

Определим необходимую угловую скорость КА согласно второму уравнению системы (7):

$p=\theta =\frac{Hv}{H^2+{(L-vt)}^2}$. (9)

Найдем производную угловой скорости (9):

$p=\frac{2H{v}^2(L-vt)}{(H^2+{{(L-vt)}^2)}^2}$. (10)

Подставляя производную (10) в первое уравнение системы (7), выразим относительное угловое ускорение ротора:

${\sigma }_1=\frac{2H{v}^2 A(L-vt)}{I(H^2+{{(L-vt)}^2)}^2}$. (11)

Зная угловые ускорения (10) и (11), найдем из третьего уравнения системы (7) зависимость управляющего момента от времени:

$M_{\alpha }=\frac{2H{v}^2 (A-1)(vt-L)}{(H^2+{{(vt-L)}^2)}^2}$. (12)

 

Рисунок 3. График зависимости управляющего момента от времени

 

На рисунке 3 показана зависимость управляющего момента Mα от времени.

 

Рисунок 4. График зависимости проекции угловой скорости p от времени

 

Рисунок 5. График зависимости угла нутации от времени

 

Для проверки результатов, полученных выше, проведено численное интегрирование системы дифференциальных уравнений (7) при соответствующих исходных условиях и внутреннем моменте (12). На рисунках 4 и 5 представлены зависимости угловой скорости несущего тела p и угла нутации Θ от времени. При этом на рисунках представлены результаты аналитического расчета по формулам (8) и (9) и численного интегрирования системы (7). Так как результаты полностью совпали, для наглядности численное решение несколько смещено вверх.

Выбор двигателя-маховика

Для обеспечения закона (12) изменения внутреннего управляющего момента, действующего на ротор, необходимо на КА установить двигатель с соответствующими параметрами. Основным параметром решаемой задачи является величина управляющего момента. Для выбранных параметров системы значение максимального управляющего момента лежит в диапазоне 0,15- 0,20 Нм.

В таблице 1 представлены основные характеристики применяемых в настоящее время двигателей-маховиков для управления ориентацией КА, масса которых составляет менее 500 кг [10]. Двигатель-маховик представляет собой электрический двигатель, на валу которого установлен массивный маховик.

 

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

Параметры	ДМ1-20	ДМ5-20	ДМ-10-25	ДМ20-250	ДМБ
Диапазон изменения кинетического момента, Н	1,0	5,0	10,0	20,0	29,4
Управляющий момент Mx, $I·i$	0,02	0,05	0,025	0,25	0,35
Диапазон изменения частоты вращения n, об/мин.	65000	63500	66000	61470	61200
Число импульсов на оборот, N	36	36	36	56	960
Масса, кг	1,4	3,8	4	11,5	17,5
Максимальная потребляемая мощность, не более, Вт	15	31	31	70	100
Потребляемая мощность в установившемся режиме, не более, Вт	3	4	5	6	7

 

Таким образом, из представленных моделей двигателей для решения поставленной задачи подходит, например, модель ДМ20-250, так как его значения управляющего момента соизмеримы с теми, которые были получены в ходе вычисления.

Заключение

В качестве основных результатов работы следует отметить преобразование имеющейся математической модели под плоскую задачу, определение закона изменения угла нутации, обеспечивающего нацеленность КА на объект наблюдения, получение зависимости внутреннего управляющего момента от времени, оценку величины управляющего момента, которая позволила выбрать требуемый двигатель-маховик. Кроме того, проведено численное моделирование движения для проверки адекватности аналитических зависимостей.

Результаты работы могут быть использованы при проектировании КА ДЗЗ, с частичной закруткой для переориентации (маховики), а также для расчетов и моделирования движения существующих спутников-гиростатов.

Работа выполнена при финансовой поддержке Российского научного фонда (Проект № 19-19-00085).

Об авторах

Алексей Владимирович Алексеев

Самарский национальный исследовательский университет имени академика С.П. Королева

Автор, ответственный за переписку.
Email: alekseev.av@ssau.ru

к.т.н., доцент кафедры теоретической механики 

Россия, Самара

Алена Олеговна Голушкова

Самарский национальный исследовательский университет имени академика С.П. Королева

Email: alenagol15935745@gmail.com

студент кафедры теоретической механики 

Россия, Самара

Список литературы

Дополнительные файлы