Оценка погрешности методов Рунге-Кутты в задаче уменьшения кинетического момента микроспутника Текст научной статьи по специальности «Физика»

УДК 629.7

ОЦЕНКА ПОГРЕШНОСТИ МЕТОДОВ РУНГЕ-КУТТЫ В ЗАДАЧЕ УМЕНЬШЕНИЯ КИНЕТИЧЕСКОГО МОМЕНТА МИКРОСПУТНИКА

© 2015 В.В. Любимов, В.С. Лашин

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

Поступила в редакцию 30.07.2015

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

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

В современной космонавтике при управлении вращательным движением микроспутника наиболее часто используются магнитные системы управления, которые основаны на взаимодействии магнитного поля спутника с геомагнитным полем [1, 2]. Такие системы в основном используются на микроспутниках, поскольку для них характерны небольшими управляющими моментами. Электромагнитные катушки применятся в современных системах управления ориентацией спутников также для уменьшения кинетических моментов двигателей-маховиков [3].

При движении по орбите на спутник оказывают влияние различные возмущающие силовые факторы [4-6]. Кроме того, в процессе отделения от базового космического аппарата микроспутник приобретает остаточные угловые скорости [7], которые требуется минимизировать для обеспечения требуемой его ориентации. Этап управления, при котором происходит минимизация угловых скоростей, называется предварительным успокоением микроспутника. Предварительное успокоение вращательного движения может производится за счет взаимодействия магнитного поля электромагнитов с геомагнитным полем.

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

Математическая модель системы управления успокоением микроспутника включает динамические и кинематические уравнения Эйлера. Ди-

Любимов Владислав Васильевич, доктор технических наук, профессор, заведующий кафедрой высшей математики. E-mail: math@ssau.ru Лашин Вячеслав Сергеевич, аспирант.

намические уравнения вращательного движения спутника имеют вид:

СК

-+ ЮхК = М + М > С1)

С с р

где К = (К , К , К ) - ве ктор кинетического

Х У 2 —»■ ^ ч

момента микроспутника; Ю = (Юх, Ю , Ю2 ) - ректор угловой скорости микроспутника; М = (М , М , М ) - вектор управляющего

С V СХ Су ' С2 ^

момента; М = (М , М , М ) - вектор

' р \ рх' ру ' рк ' *

возмущающих моментов.

Механический управляющий момент, действующий на микроспутник спутник МС, в проекциях на оси X,Y,Z главной связанной системы координат OXYZ имеет вид:

M„

LB - LB,

y z z y'

M = LB

cy z x

L B ,

xz

M

LB - LB,

x y y x'

(2)

где Ь = (Ьх, Ьу, Ьг ) - вектор магнитного момента микроспутника.

Предположим, что закон управления вектором магнитного момента электромагнитных катушек при предварительном успокоении имеет вид [8]:

Ь 3

С Б2

Здесь кс - коэффициент управления, B - модуль вектора магнитной индукции геомагнитного поля.

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

(4)

dв / = az cos p + ay sin <р + а>ш sin у -«ZM cos у; dy/ dt = (< sin p-<y cos p) / sin в -

-<<+<<vctge cos у+<<с&в sin у;

dp/ dt = <x- dy / di cose — <5zs in в sin у.

Здесь ccM, Шум, &>ZM - угловые скорости вращения магнитной системы координат относительно неподвижной геоцентрической системы координат [5].

Численное решение системы из уравнений (1), (4) рассматривается совместно с решением уравнений, характеризующих движение микроспутника по орбите. Решение осуществляется посредством использования двух явных методов: четырехэтапного и семиэтапного методов Рун-ге-Кутты. В семиэтапном методе используется формула Дормана-Принса [9]. Четырехэтапный метод является классическим, и в дальнейшем будем использовать сокращенное название РК4.

Конкретный метод определяется числом 5 и коэффициентами Ъ1, а^ и С1. Эти коэффициенты упорядочивают в таблицу [9], называемую таблицей Бутчера. Таблица Бутчера при использовании формулы Дормана-Принса содержится в [9]. Используя коэффициенты из указанной таблицы, составим формулы для семиэтапного метода Рунге-Кутты:

к = / (хп, Уп);

к, = /(Хп + 5 А, Уп + 5 к1к);

3 3 9

к3 = /■ (х + — А, у + А(— к +— к,)); п 10 40 40

, г, 4, , 44, 56, 32

к4 = г (х + — А, у + А(— к--к, +—к3));

4 7 п 5 7п 45 1 1 5 2 9 3

, rr 8, ,,19372 25360,

к = f (x + — h, y + h(-k--к

5 n 9 6561 1 2187 2

■64448 кз-312к4)); (5)

6561 3 729 4

, fl , 9017 355 46732,

к = f (x + h, y + h(-к--к H--к

n n 3168 33 2 5247 3

49 , 5103 ,

--k4--к,));

176 4 18656

35 500 125 2187 11

k7 = f(x + h, y + h(-к +-k3 +-к--к +— k6));

7 n n 384 1 1113 3 192 4 6784 5 84 6

35 500 125 2187 11

y . = y + h(—к +-к +—к--к +—к)■

Л.+1 Л 384 1 1113 3 192 4 6784 5 84 6

ЯВНЫЙ одношаговыи семиэтапныш

МЕТОД РУНГЕ-КУТТЫ (РК7)

Семейство явных одношаговых методов Рун-ге-Кутты задается формулами

д

Уп +1 = Уп Ьгкг ,

I=1

где к - величина шага интегрирования. Здесь вычисление нового значения проходит в 5 этапов:

к1 = ¥(х, Уп X

к, = ¥(хп + с,А, Уп + а2к),

к,д = ¥(Хп + СА Уп + 0д 1к1 + 0д2к2 + - + 0д,д-А-1 )-

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

Четырехэтапный и семиэтапный (5) методы Рунге-Кутты для решения системы дифференциальных уравнений, включающей уравнения (1), (4) реализуются с помощью программы, интерфейс которой показан на рис. 1.

В интерфейсе программы вводятся значения моментов инерции (1х, I , 12 ) микроспутника, начальные условия (юх (0), (Оу (0), (й2 (0), 0(0), ¡//(0), (р(0)) для интегрирования дифференциальных уравнений (1) и (4), шаг (к) и конечное время интегрирования ^Кон ) .

Результаты расчетов с помощью программы имеют табличный вид в файле формата «1x1.».

Приведем пример расчета в программе с

Введете значения: Исходные данные:

15

Начальные условия: GmegaíK GrnegaY GrnegaZ tetta psi phi

0.3

0.5

0.4

[if

0.524

1.047

0.785

Конечное время: 1400

Uar интегрирования

P-

Посчитать PK4

Посчитать PK7

Закрыть

Рис. 1. Интерфейс программы

параметрами по умолчанию, которые показаны на рис. 1, методами РК4 и РК7, но шаг интегрирования зададим равным h = 0,4. Результаты численного интегрирования приведены на рис. 2, 3.

Как следует из рис. 2, 3, результаты расчётов с

использование двух методов РК4 и РК7 совпадают, что говорит о малой погрешности численного интегрирования.

Рассмотрим, как зависит погрешность результатов, которые находятся с помощью методов РК4

200

250

300

350

400

I, С

I, С

-го*4

-(0x7

-<¡0^4 -Шу7

в3 рад/с 0,41 0,405 0,4 0,395

Д/\А-

■ Из4 ■гаЛ

50

100

150

200

250

300

350

400

I, С

Рис. 2. Графики Юх (^), Ю (?), Ю2 ^)

Рис. 3. Графики в(г), у/(г), ср(г)

Таблица 1. Расчет погрешностей интегрирования при использовании РК4 и РК7

h ^кон С ^>% 4 ^>%

0,01 400 0 0 0 0 0 0

0,05 400 9,35E-06 3,95E-05 1,99E-07 1,71E-05 2,00E-06 2,70E-06

0,1 400 1,92E-03 2,77E-03 2,04E-06 1,55E-05 2,47E-05 1,28E-05

0,2 400 3,08E-03 1,44E-01 3,78E-05 1,54E-03 3,60E-04 4,70E-04

0,3 400 1,53E-02 1,27E-01 1,55E-04 3,46E-03 3,39E-03 3,12E-03

0,4 400 6,92E-01 3,86E-01 4,45E-04 8,27E-03 3,62E-03 9,87E-03

0,5 400 2,90E+00 3,60E-01 9,58E-04 1,61E-02 3,79E-03 1,68E-02

0,6 400 1,65E+01 1,04E+00 1,72E-03 2,39E-02 5,69E-03 2,55E-02

0,7 400 2,08E+01 6,16E+00 2,68E-03 2,95E-02 1,27E-02 3,60E-02

и РК7, от величины шага интегрирования. Для это-то составим табл. ^где £Фж ,£а ,етг ,£в,£ч/, - расхождение результатов, поЛучаемых методами РК4 и РК7 в процентах.

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

Из табл. 1 следует, что чем больше шаг интегрирования, тем больше величина расхождения решений между «РК4» и «РК7». При шаге 0,01 погрешность сводится к величинам порядка 10-9. При шаге большем 0,6 погрешность превышает 5%.

ВЫВОДЫ

Явные одношаговые методы Рунге-Кутты позволяют эффективно решать задачи по моделированию процесса успокоения вращательного движения микроспутника. Результаты численного моделирования показывают, что при шаге h < 0,4 интегрирования системы дифференциальных уравнений расхождение результатов при использовании семиэтапного метода Рунге-Кут-ты и классического метода РК4 не превышают и 1%. Однако, при моделировании движения микроспутника на значительных интервалах времени следует использовать методы более точные, чем РК4.

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

1. Белецкий В.В., ХентовА.А. Вращательное движение намагниченного спутника. М.: Наука, 1985. 288 с.

2. Коваленко А.П. Магнитные системы управления космическими летательными аппаратами. М.: Машиностроение, 1975. 248 с.

3. Lyubimov V.V., Malyshev V.I., Semkin N.D. Attitude control of small spacecraft with allowance made for failures in the system for decreasing angular momentums of flywheels // Gyroscopy and Navigation, 2013. Volume 4. Issue 4. Pages 216-222.

4. Любимов В.В. Об особенностях в возмущенном вращательном движении спутника с сильным магнитом на борту // Изв. вузов. Авиационная техника. 2009. №2. С.29-31.

5. Любимов В.В. Внешняя устойчивость резонансов при движении асимметричного твердого тела с сильным магнитом в геомагнитном поле// Известия РАН. Механика твердого тела. 2010. № 1. С.13-27.

6. Лебедев А.С., Любимов В.В., Сёмкин Н.Д. Моделирование управляемого движения микроспутника с магнитными и гравитационными исполнительными органами // Полёт. 2012. №7. С.39-44.

7. Калаев М.П., Любимов В.В., Сёмкин Н.Д. Полунатурное и имитационное моделирование процесса отделения микроспутника // Гироскопия и навигация. 2014. № 2 (85). С.52-60.

8. Семкин Н.Д., Любимов В.В., Малышев В.И. Моделирование законов функционирования магнитных исполнительных органов при ориентации микроспутника по местной вертикали// Физика волновых процессов и радиотехнические системы. 2012.

Том.15. № 1. С.103-108. венных дифференциальных уравнений. Нежесткие

9. Хайрер Э., Нёрсетт С., ВаннерГ. Решение обыкно- задачи. М.: Мир, 1990. 512 с.

THE ESTIMATION ERROR OF THE RUNGE-KUTTA METHOD IN THE TASK OF REDUCING AN ANGULAR MOMENTUM OF A MICROSATELLITE

© 2015 V. V. Lyubimov, V. S. Lashin

Samara State Aerospace University named after Academician S. P. Korolev (National Research University)

We consider the problem of numerical modeling of damping of a microsatellite with electromagnetic executive bodies. At the numerical integration of equations of motion of microsatellite used a seven-step Runge-Kutta method. The estimation error in the numerical integration of a equations of a motion through the application a seven-step and four-stage explicit Runge-Kutta methods is performed. Keywords: error, angular momentum, microsatellite, numerical integration, Runge-Kutta methods.

Vladislav Lyubimov, Doctor of Technics, Professor, Head at the Mathematics Department. E-mail: math@ssau.ru Vyacheslav Lashin, Graduate Student.