CC BY
43
DOI 10.23859/1994-0637-2017-6-81-6 УДК 681.5
Толпегин Олег Александрович
Доктор технических наук, профессор, Балтийский государственный технический университет «ВОЕНМЕХ» им. Д.Ф. Устинова (Санкт-Петербург, Россия) E-mail: [email protected]
Литвинова Полина Юрьевна
Ассистент, Балтийский государственный технический университет «ВОЕНМЕХ» им. Д. Ф. Устинова (Санкт-Петербург, Россия) E-mail: [email protected]
УПРАВЛЕНИЕ МАЛЫМ КОСМИЧЕСКИМ АППАРАТОМ С ИСПОЛЬЗОВАНИЕМ ДВИГАТЕЛЯ-МАХОВИКА НА ОСНОВЕ МЕТОДА УПРАВЛЕНИЯ С ПОВОДЫРЕМ
Аннотация. В статье рассматривается синтез системы стабилизации малого космического аппарата (МКА) с использованием двигателя-маховика при действии возмущений, статистические свойства которых неизвестны, но ограничены по абсолютной величине. Решение проводится с использованием математической модели системы стабилизации по углу тангажа. Управление выбирается с использованием эталонного движения - движения поводыря. Траектория движения поводыря выбирается с помощью метода фазовых плоскостей. Приведен алгоритм определения допустимого диапазона возмущений с заданными статистическими свойствами для решения исходной задачи с допустимой точностью, исходя из ограниченных возможностей системы управления движением объекта.
Ключевые слова: малый космический аппарат, двигатель-маховик, дифференциальные игры, метод управления с поводырем, случайное возмущение_
© Толпегин О.А., Литвинова П.Ю., 2017
Tolpegin Oleg Aleksandrovich
Doctor of Technical Sciences, Professor, Baltic State Technical University "VOENMEKH" named after D.F. Ustinov (St Petersburg, Russia) E-mail: [email protected]
Litvinova Polina Jur'evna
Assistant, Baltic State Technical University "VOENMEKH" named after D.F. Ustinov (St Petersburg, Russia) E-mail: [email protected]
A SMALL SPACECRAFT CONTROL WITH THE USE OF FLYWHEEL ON THE BASIS OF THE CONTROL METHOD WITH
A GUIDE
Abstract. The article deals with the synthesis of the stabilization system for a small spacecraft using a flywheel engine under the action of perturbations whose statistical properties are unknown but are limited in absolute magnitude. The solution is carried out using a mathematical model of the pitch angle stabilization system. Control is selected using the reference movement - the guide's movement. The path of the guide's movement is selected using the phase plane method. An algorithm is given to determine the permissible range of disturbances with given statistical properties for solving the initial problem with a permissible accuracy, taking into account limited capabilities of the object's motion control system.
Keywords: small spacecraft, flywheel engine, differential games, control method with guide, random perturbation
Введение
Рассматривается управление малым космическим аппаратом с использованием двигателя-маховика при действии возмущений, статистические свойства которых
неизвестны, но ограничены по абсолютной величине. Задача рассматривается как дифференциальная игра двух игроков: в качестве одного игрока, стратегия которого неизвестна, выступает возмущающий момент, а в качестве другого игрока - космический аппарат, управление которого нужно выбрать. Особенность решения состоит в том, что множество управлений малого космического аппарата разбивается на два подмножества: одно подмножество управлений используется для компенсации возмущений, а второе - для решения исходной задачи при отсутствии возмущений. Полученная траектория называется траекторией поводыря. Разделить множество управлений на два удается в редких случаях, поэтому в данной работе заранее выделяется подмножество управлений для решения исходной задачи при отсутствии возмущений, строится траектория поводыря, а затем в результате моделирования определяется допустимый диапазон возмущений с заданными статистическими свойствами, при которых возможно решение исходной задачи с допустимой точностью. Данный подход можно использовать как для линейных, так и нелинейных систем.
Основная часть
Рассмотрим систему стабилизации угла тангажа малого космического аппарата (МКА) с двигателем-маховиком. В полете на МКА действует возмущающий момент, статистические характеристики которого неизвестны, но известны пределы, в которых изменяется величина возмущающего момента.
Данную задачу будем рассматривать как дифференциальную игру двух игроков, так как в качестве одного игрока, стратегия которого неизвестна, выступает возмущающий момент, а в качестве другого игрока - космический аппарат, управление которого нужно выбрать.
Динамика МКА вокруг поперечной оси 02 при действии возмущений определяется системой уравнений:
=—к«-+*); (1)
а Уж = ю2 +«0,
где = Jz — I; и = Мупр2 = = I ■ ^^^ск ' - момент инерции космического
аппарата без учета момента инерции маховика; Jz - момент инерции космического аппарата по оси 02; I - момент инерции маховика; и - управляющий сигнал; Мупр2 -
управляющий момент двигателя-маховика по оси 02, который появляется при разгоне маховика и приложен к корпусу летательного аппарата; К - кинетический момент двигателя-маховика; - проекция абсолютной угловой скорости маховика по
оси 02; ю2 - угловая скорость космического аппарата по оси 02; ^ - угол тангажа;
«0 кт ки 0 - круговая скорость вращения космического аппарата; ® и и - передаточные
числа системы управления; £ - возмущающий момент [1].
Управляющий момент удовлетворяет ограничению:
\и\ < а. (2)
Возмущающий момент £ характеризует внешние неконтролируемые возмущения, действующие на космический аппарат и удовлетворяющие ограничению:
|3<Р = р-а,0<р<1.
(3)
Задание ограничения на возмущения в виде (3) дает возможность оценить влияние возмущений на управление МКА. Заданы граничные условия:
при г== о, юг(0)=ю3о = -«о> Э(0)=эо; (4)
при г = Т, юг(Т) = юЗТ = -ч>, 9(Т) =8З.
Требуется найти управление МКА, обеспечивающее минимум критерия
Т
3 = |М ^шти, (5)
о
но при этом необходимо учитывать действие возмущений.
Момент окончания переходного процесса Т не фиксирован.
Управление системой (1) при наличии возмущений будем рассматривать как антагонистическую дифференциальную игру двух игроков: первый игрок выбирает управление и(г), а второй игрок - возмущение .
Первый игрок стремится минимизировать (5), а второй - максимизировать этот критерий.
Для решения задачи используем игровой метод управления с поводырем [2].
Множество управлений и МКА, заданного ограничением (2), разобьем на два подмножества: и1 и и2.
Подмножество и1 используем для решения сформулированной задачи при отсутствии возмущений. Эту траекторию назовем траекторией поводыря. Для подмножества и1 ограничение на управление МКА зададим в виде:
N < а(1 -р), (6)
где коэффициент р изменяется в диапазоне от о до 1.
Подмножество и2 используем для компенсации действия возмущений. Для подмножества и2 ограничение будем задавать в виде:
\и(/)| < ар . (7)
Если р = 1 (см. формулу (6)), то все ресурсы управления используются для компенсации возмущений и и1=о. При р = о подмножество и2=о.
Для исследования действия возмущений предлагается следующий алгоритм.
1. Зададим коэффициент р, который изменяется в диапазоне от о до 1, и для множества и1 ограничение на управление МКА зададим в виде (6).
2. При ограничении (6) для системы (1) решим задачу оптимального быстродействия при отсутствии возмущений и построим оптимальную траекторию. Эту траекторию назовем траекторией поводыря и запишем в память БЦВМ.
3. Решаем задачу о максимальном сближении траектории движения МКА с траекторией движения поводыря при действии возмущений с ограничением (3) при заданном значении коэффициента р. При этом используем все возможности управления МКА, т.е. решаем задачу с ограничением (2).
Действие возмущений будем оценивать по отклонению реальной траектории движения МКА от траектории движения поводыря.
Рассмотрим расчет траектории поводыря при заданном значении коэффициента р. Для расчета траектории движения поводыря при отсутствии возмущения систему уравнений (1) с учетом (6) запишем в виде:
d ш
г/- =К,-юг - К • и-(1 -р)));
'dt'
d Vdt +«0,
(8)
где коэффициент р влияет на возможности управления поводыря.
Решим задачу оптимального быстродействия для системы (8) с критерием (5) при заданных начальных и граничных условиях (4). Для решения задачи используем принцип максимума Л.С. Понтрягина [3]. Составим функцию Гамильтона:
Н - •(-*„ -шг - ku • и-(1 -р))-J3-1 +Т 2 .(шг +«о)-1 где переменные ^1(t), удовлетворяют системе уравнений:
(9)
d^1/ - -5Н/ - W К/ -W ■ /dt~ /дюг " '/J
d Т2
--5Н/ - 0 /dt - /59 - 0'
Решение системы (10) имеет вид:
T1(t) - С1 • е J + С2,
К J3;
Т2 (t) - C2 - const. Функция Гамильтона достигает максимума при управлении вида:
У t с /
и --sign(C1 • е' 3 + / • J3)-а
(10)
(11)
(12)
Оптимальное управление й (?) является кусочно-постоянным и меняет знак не
более одного раза, так как функция (^ имеет не более одного нуля [3].
Для задачи синтеза используем метод фазовых траекторий, при этом решение сводится к определению линии переключения, разделяющей фазовое пространство координат ю2 и^ на две области: в одной области оптимальное управление
и (?) =+а, а в другой и (?) = -а . При попадании фазовой траектории на линию переключения знак оптимального управления изменяется на противоположный.
Для построения линии переключения первое уравнение системы (8) поделим на второе уравнение этой системы:
йъ/ =к-®г - к ■ и -(1 -р))/
/ё 9 /(+ Ц,))
и после разделения переменных
^ • (юг +«0)
-Г —ёЮг = Г ё 9; J кю -ю + к • и -(1 - р) 9
ш,0 ю 2 и \ г} 90
9
проинтегрируем от 0 до 9Г правую часть и от юг0 до левую часть:
»Г =9о -Ук • К -юго) + ^ • (ки •м•(1-р)-кю -о).
X 1П ((К •М • (1-Р) + кю-Ю* V
х 1П| /к • м • (1 -р)+кю-юго)^
(13)
где 90 и ю20 - точка, в которую должна прийти система на фазовой плоскости,
»Г - абсцисса на фазовой плоскости, находящаяся на линии переключения, ю^ -текущее значение угловой скорости.
Алгоритм выбора оптимального управления поводыря имеет вид:
и = а, если <; 2
V '
ю2 (О > 0 и 9(0 >»Г
ю2 < 0 и 9(0 >»Г ю2 (0 > 0 и 9(/) <9
и и =-а, если ^ Г. (14)
X (0 < 0 и 9(0 <9Г
Построение линии переключения определяет траекторию поводыря.
При управлении с поводырем управление МКА нужно выбирать с использованием алгоритма (14), так как в этом случае траектория МКА будет наискорейшим образом подходить к линии переключения, тогда после попадания на линию переключения, построенную для поводыря, траектория МКА будет колебаться относительно этой линии и отклонение МКА в момент Т от заданного значения будет не больше допустимой величины е, обусловленной дискретностью выбора управления.
После расчета траектории поводыря можно оценить величину возмущений, при которых возможно приведение МКА в заданную точку на фазовой плоскости за допустимое время при заданном значении коэффициента р. Для этого используем следующий алгоритм.
1. Задаем значение коэффициента р, с использованием рассмотренного выше алгоритма строим траекторию поводыря и находим время переходного процесса ТПсе.
2. Задаем возмущение £ с заданными статистическими характеристиками и подтягиваем траекторию движения МКА при действии этого возмущения к траектории поводыря с использованием алгоритма (14), при этом используем полные возможности управления МКА, т. е. используем управление с ограничением (2).
3. Находим при t = ТПОВ отклонение реальной траектории движения МКА от траектории движения поводыря, используя критерий:
0 = (Тпов ) - ^ )2 + (9 (Тпов ) - )2. (15)
4. В результате многократного решения этой задачи находим математическое ожидание 0 критерия (15).
5. Если 0 меньше допустимой величины £, то МКА устранит действия возмущений с данными статистическими характеристиками к моменту ТюВ. В этом случае статистические характеристики возмущений увеличиваем на некоторую дискретную величину А и находим новое значение 0 .
6. Если 0 больше допустимой величины е, то МКА не сможет компенсировать возмущения с данными статистическими характеристиками к моменту ТПОВ при данном р.
Изменяя значение р и вычисляя для каждого р максимально допустимые значения возмущений, найдем допустимый диапазон возмущений с заданными статистическими свойствами, при котором математическое ожидание 0 меньше допустимой величины.
Рассмотрим результаты моделирования с использованием данного алгоритма.
Исследование проводилось при следующих значениях параметров гипотетического МКА: тка = 31,6, кг - масса космического аппарата; 1м = 0,195 -10-2, 12 = 0,8953 кг • м2 - моменты инерции маховика и летательного аппарата с учетом момента инерции маховика по оси О2; и0 = 0,001 рад/с, кт = 0,083, кй = 1.
Моделирование проводилось при различных р, и ограничениях
й < а, а = Мупр _тах = 20 •Ю-3 Н • м, |^|<|3 = К-р-а.
Начальные условия имели следующие значения:
г0 = 0, юг (0) = -0,075 рад/с, 9(0) = 0,26 рад.
Фазовые координаты точки, в которую должна прийти система:
юг (Т) = ю3Т = -о0 = -0,001рад/ , 9(Т) = 93 = 0,022 рад.
Время движения поводыря определялось как первый момент времени, при котором значение критерия (15) становилось по абсолютной величине меньше
е = 10-3 рад/с.
Шаг численного интегрирования к = 0,01с.
На рис. 1 представлены графики линий переключения для р = 0; 0,2; 0,4; 0,6; 0,8 .
о.а 0.6 0-4 0.2 0 0.2 р=0. р-0 и-«, р=£>, и--.! Ч¥г<0 -1'-С- 2 и=а. ^>0
Р=0Д -р=0.4. и -м>а
-----р=СЛ. и=-«. Ч/^О -р=0.6. |.-<г wz>0
-----(.-О 0 1-- а г>2-'0 -р=0А и=л.
V -----р=0.3. и=- «.ч^чО
р-ов.
о-о, /
1 . .__1_I_I_I_I_I_1" -_I
-10 -й -6 -4 -2 0 г 4 6 в 10
ВД.гвй
Рис. 1. Линии переключения при различных значениях р
На рис. 2, 3 представлены графики изменения ю2 (?) и $(?), построенные при
различных р = 0,3; 0,5; 0,7.
Рис. 2. График переходного процесса Ю2 (?) при различных р
О 35
Рис. 3. График переходного процесса 9(0 при различных р
На основе рассмотренного выше алгоритма оценим величину возмущений, при которых возможно приведение МКА в заданную точку на фазовой плоскости за время движения поводыря при заданном значении коэффициента р. Возмущения задавались в виде случайных функций, распределенных по нормальному закону с нулевым математическим ожиданием и заданным среднеквадратическим отклонением с, причем величина 4 была ограничена < 0,006 Н • м.
В таблице представлены математические ожидания значения критерия оптимальности (15) при различных р и с. При вычислении 7 проводилось по 20 испытаний. Допустимое значение е = 103.
Таблица
Математические ожидания 7 значения критерия оптимальности при различных значениях среднеквадратического отклонения а и р
ст(Н/м ) 0,02 0,03 0,04 0,05 0,06 0,07 0,08 0,09 0,1
7 (р = 0,2) 0,00019 0,00021 0,00013 0,00033 0,00117 0,00142 0,00222 0,00351 0,00334
7 (р = 0,4) 0,00021 0,00019 0,00023 0,00024 0,00073 0,00102 0,00075 0,00149 0,00254
7 (р = 0,б) 0,00014 0,00018 0,00020 0,00025 0,00049 0,00114 0,00090 0,00156 0,00122
7 (р = 0,8) 0,00007 0,00013 0,00014 0,00038 0,00068 0,00090 0,00094 0,00123 0,00128
Из таблицы следует, что если допустимое значение 7 < е, то с случайного возмущения при р = 0,2 не должна превышать 0,05 Н/м, при р = 0,4 и р = 0,6 - 0,06 Н/м, а при р = 0,8 - 0,08 Н/м.
Выводы
Результаты моделирования показывают работоспособность предлагаемого игрового алгоритма для оценки возможностей управления МКА с двигателем-маховиком при действии разного вида возмущений, ограниченных по абсолютной величине.
Литература
1. Алексеев К.Б., Бебенин Г.Г. Управление космическим летательным аппаратом. М.: Машиностроение, 1964. 404 с.
2. Красовский Н.Н., Субботин А.И. Позиционные дифференциальные игры. М.: Наука, 1974. 456 с.
3. Понтрягин Л.С., Болтянский В.Г., Гамкрелидзе Р.В., Мищенко Е.Ф. Математическая теория оптимальных процессов. М.: Наука, 1969. 384 с.
References
1. Alekseev K.B., Bebenin G.G. Upravlenie kosmicheskim letatel'nym apparatom [Space flight control system]. Moscow: Mashinostroenie, 1964. 404 p.
2. Krasovskii N.N., Subbotin A.I. Pozicionnyie differencial'nyie igry [Positional differential games]. Moscow: Nauka, 1974. 456 p.
3. Pontriagin L.S., Boltjanskii V.G., Gamkrelidze R.V., Mishhenko E.F. Matematicheskaia teo-riia optimal'nyh processov [Mathematical theory of optimal processes]. Moscow: Nauka, 1969. 384 p.
Толпегин О.А., Литвинова П.Ю. Управление малым космическим аппаратом с использованием двигателя-маховика на основе метода управления с поводырем // Вестник Череповецкого государственного университета. 2017. №6 (81). С. 44-52.
For citation: Tolpegin O.A., Litvinova P.J. A small spacecraft control with the use of flywheel on the basis of the control method with a guide. Bulletin of the Cherepovets State University, 2017, no. 6 (81), pp. 44-52.
