Научная статья на тему 'Оптимальное управление космическим аппаратом с солнечным парусом и двигателем импульсной тяги'

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

CC BY
75
21
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Оптимальное управление космическим аппаратом с солнечным парусом и двигателем импульсной тяги»

появляется характерный безразмерный параметр N = umaxp3/c2. При численном решении полагалось, что N = 0.5.

Значение минимизируемого функционала в первом случае оказалось равным JHeocoб = 0.4897, а во втором - J0C06 = 0.4936.

Работа выполнена при финансовой поддержке РФФИ (проект Ks 12-01-00 165).

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1, Челноков Ю. Н. Оптимальная переориентация орбиты космического аппарата посредством реактивной тяги, ортогональной плоскости орбиты // Математика, Механика: сб. науч. тр. Саратов: Изд-во Сарат, ун-та, 2006, Вып. 8, С, 231-234,

УДК 629.78

Я. Г. Сапунков

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

С использованием KS-переменных (переменные Кустаанхеймо -Штифеля) с помощью принципа максимума Понтрягина решена пространственная задача оптимального управления о встрече двух космических аппаратов (КА), один из которых неуправляемый и движется только под действием силы притяжения к Солнцу, второй аппарат управляется с помощью солнечного паруса и импульсного двигателя. Функционал, определяющий качество процесса управления, представляет собой линейную комбинацию с весовыми множителями двух критериев: времени и суммы величин импульсов, затраченных на процесс управления. Приводятся результаты численного решения задачи.

1. KS-переменные u = (u0, ui, u2, u3), s = (s0, si, s2, s3) связаны с векторами положения центра масс КА и его скорости ну соотношениями [1, 2]

r = PT(u)u; у = (2/r)PT(u)s; r = |r| = u2, h = v2/2 - jM/r.

(1.1)

hM масса притягивающего центра, y - гравитационная постоянная.

Переменная т, функциями которой являются KS-переменные u, s, h, связана со временем t соотношением

dt/dT = u2. (1.2)

Тяга р,0/ солнечного паруса, отнесенная к единице массы КА, определяется по формуле, в которой п - единичный вектор нормали к плоскости паруса, обращенной от Солнца, в - угол между векторами г и п, Д -коэффициент, характеризующий площадь паруса, отнесенную к единице массы аппарата [3]:

,С082 в , -4 ч \ 2 (л

р,о1 = а г2 п = а (и) (Р (и)и, п) п. (1.3)

Если через Я обозначить характерный масштаб длины, например радиус орбиты Земли, на которой находится управляемый аппарат в начальный момент времени, то связь между размерными и безразмерными переменными, которые обозначаются верхним индексом «*», будет определяться соотношениями:

Н = (7М/Я)Н*; 8 = (7М)1/28*; и = Я1/2и*; т = (Я/(7М))1/2т*; г = Я(Я/(7М ))1/2г*; Д = 7МД*; г = Яг*; V = (7М/я)1/2у*.

Безразмерная величина Д* определяет отношение максимальной величины тяги солнечного паруса к силе притяжения аппарата к

Я.

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

т

С = (Со, С1, С2, Сэ) и Ю = (Л0, Р1, Р2, Д$) соотношениями

иа = C cos(kra) + D sin(kra); sa = k (D cos(kra) - C sin(kra));

1 т"а о

k = (-ha/2)1/2; ha = - (C2 + D2)-1 < 0; t = f (u)2dra

-1

(u

Та > Тан; C = const; D = const.

(1.4)

Если через К обозначить безразмерную величину импульса тяги, отнесенного к единице массы КА, то она будет связана с размерной величиной изменения вектора скорости |Av^| под действием импульса соотношением

|Дуг| = (7MR-1)1/2 К

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

мального управления, имеет вид

I = о^к + К (1.5)

¿=1

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

2. В декартовой системе координат Ох1 х2хз, начало которой находится в центре Солнца, а плоскость Ох1х2 совпадает с плоскостью орбиты Земли, в начальный момент времени в безразмерных переменных положение и скорость управляемого аппарата, который находится на орбите Земли, определяется координатами

Х1 = 1.0, Х2 = 0.0, хз = 0.0, VI = 0.0, = 1.0, Уз = 0.0. (2.1)

Неуправляемый аппарат находится вблизи орбиты Марса. Положение и скорость неуправляемого аппарата в начальный момент времени в безразмерных переменных, вектора С и О, значения вели чин к, Тан, входящих в закон движения в КБ-переменных для неуправляемого аппарата (1.4), определяются соотношениями

Х1 = 0.76, Х2 = 1.3164, хз = 0.05, VI = -0.7019, У2 = 0.4052, Уз

Со = 0.0, С1 = 1.0679, С2 = 0.6163, Сз = 0.0234, = -0.03196, А = -0.6171, Л2 = 1.0671, £з = -0.01484,

К = -0.3288, к = 0.4054, тан = 0.0.

Безразмерный множитель, характеризующий тягу солнечного паруса, ё, = 0.1, весовые множители в функционале (1.5) а0 = 1.0, а1 = 0.75. Расчеты показали, что для данной задачи в оптимальном управлении импульсный двигатель включался в начале и в конце процесса управления, обеспечивающего мягкую встречу аппаратов. Безразмерные координаты первого импульса, который осуществлялся при £ = 0, и его величина определяются соотношениями:

КН1 = -0.4247, Кн2 = 0.3965, Кн3 = 0.0021, К = 0.5810. (2.3)

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

= -0.0267. (2.2)

У1 = -0.4247, У2 = 1.3965, Уз = 0.0021.

(2.4)

Далее до момента гк = 1.4502 управляемый аппарат движется под действием солнечного паруса, в момент подлета к неуправляемому аппарату его положение и скорость в безразмерных координатах определяются соотношениями:

Х1 = 0.3756, Х2 = 1.4725, х = 0.0009, (2 5)

Vl = -1.1077, У2 = 0.5555, ^з = -0.0008. ( )

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

КК1 = 0.3224, Кк2 = -0.7564, Кк3 = -0.0369, К = 0.8231. (2.6)

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

V! = 0.7853, у2 = -0.2009, ^з = -0.0377. (2.7)

Функционал качества процесса управления (1.5), обеспечивающего мягкую встречу двух аппаратов, I = 2.5031.

Работа выполнена при финансовой поддержке РФФИ (проект Жд 1201-00 165).

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Штифель Е., Шейфеле Г. Линейная и регулярная небесная механика. М,: Наука, 1975. 304 е.

2. Сапунков Я. Г. Решение задач оптимального управления космическим аппаратом с ограниченной и импульсной тягой в КЯ переменных // Мехатроника, автоматизация, управление. 2010. 3. С. 73-78.

3. Поляхова Е. Н. Космический полет с солнечным парусом: проблемы и перспективы. М,: Наука, 1986. 304 с.

УДК 519.257

Г. П. Шиндяпин, A.A. Матутин

ВОЗНИКНОВЕНИЕ НЕВЫРОЖДЕННОЙ УДАРНОЙ ВОЛНЫ В ВОЗДУХЕ ПРИ РЕФРАКЦИИ УДАРНОЙ ВОЛНЫ В ОКЕАНЕ

Для различных режимов нелинейной рефракции ударных волн (УВ) на поверхности океана с образованием волны разрежения или отраженной УВ методами асимптотической теории коротких волн [1] установлены области существования режимов с невырожденным фронтом преломленной У В в воздухе.

i Надоели баннеры? Вы всегда можете отключить рекламу.