Построение стабилизирующего управления для стационарных движений гиростата с полостью с вязкой жидкостью Текст научной статьи по специальности «Физика»

УДК 62.534 (031)

ПОСТРОЕНИЕ СТАБИЛИЗИРУЮЩЕГО УПРАВЛЕНИЯ ДЛЯ СТАЦИОНАРНЫХ ДВИЖЕНИЙ ГИРОСТАТА С ПОЛОСТЬЮ С ВЯЗКОЙ ЖИДКОСТЬЮ

© 2013 А.В. Алексеев, С.П. Безгласный, В.С. Красников

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

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

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

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

ВВЕДЕНИЕ

Задачи о пространственной ориентации спутников и летательных аппаратов на орбите имеют важное прикладное значение и широко рассматриваются исследователями во многих работах. Сферические движения летательных аппаратов относительно центра масс моделируются движениями твердых тел или систем тел, в частности, гиростатов. Например, движения под действием гравитационного момента на круговой орбите орбитальной станции Мир в 19992001 гг. и спутника Фотон М-3 в 2007 г. относительно своих центров масс в работе [1] моделируется с помощью осесимметричного гиростата с гиростатическим моментом, направленным по его оси симметрии. Основные способы и принципы управления вращательным движением тел и их систем были обозначены уже давно, например, в работах [2-4]. Современными отечественными и зарубежными авторами активно исследуются задачи о резонансных режимах и бифуркациях стационарных движений спутников [5, 6], о хаотических движениях и методах их устранения [7, 8], о стабилизации заданных программных движений гиростатов различной структуры [9, 10] .

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

Алексеев Алексей Владимирович, кандидат технических наук, доцент. E-mail: [email protected] Безгласный Сергей Павлович, кандидат физико-математических наук, доцент, докторант. E-mail: [email protected] Красников Виктор Сергеевич, аспирант. E-mail: [email protected]

кой: развитие ракетной техники, содержащей большое количество жидкого топлива, проблемы сейсмостойкости резервуаров для хранения жидкости и другие. Основополагающие результаты теории твердых тел с полостями, содержащими жидкость, изложены, например, в работах [11, 12]. На движение космического аппарата (КА) вокруг центра масс существенно влияют вращающиеся роторы и жидкое топливо, обуславливая необходимость изучения динамики и свойств устойчивости множества движений механической системы, представляющей собой однороторный гиростат с полостью, заполненной жидкостью. Этой проблеме посвящена предлагаемая работа.

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

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

Рассмотрим пространственное движение гиростата, представленного на рисунке 1, как системы двух динамически симметричных связанных тел с общей осью вращения, совпадающей с осями динамической симметрии обоих тел: тело-носитель, содержащее полость, целиком заполненную жидкостью большой вязкости, и ротор. Пусть А1 = В, С1 - главные моменты инерции тела-носителя вместе с жидкостью, А2 = В2, С2 - главные моменты инерции ротора. Координата а представляет собой угол поворота ротора относительно носителя. Неподвижная точка О гиростата совпадает с общим центром масс системы и лежит на оси динамической симметрии обоих тел.

7, 7

У

* У

Рис. 1. Гиростат с полостью с жидкостью

Уравнения движения однороторного гиростата с полостью, заполненной жидкостью, в проекциях на оси связанной с носителем системы координат Оху2 , следуя [13], запишем в следующем виде:

Ар + (С - В)дг + С2да = тх. Вд + (А - С) рг - С2ра = ту Сг + (В - А )рд + С2а = т2,

с 2 (г + а )= м2,

(1)

(2)

где А = А+ А, В = В1 + В2, С = С1 + С2 главные моменты инерции гиростата, вычисленные в системе координат Охуг. Выражение (2) представляет собой уравнение относительного вращения ротора, М2 - момент, действующий со стороны несущего тела на ротор. В дальнейшем будем рассматривать движение при предположении о равенстве нулю указанного момента. Тогда уравнение (2) приводится к следующему виду:

а3 =-г.

(3)

Правые части уравнений (1) представляют собой проекции момента сил, действующих на несущее тело со стороны полости с жидкостью, и согласно работе [11] определяются равенством:

т

Ж

ш х Ь

(4)

Гиростатический момент жидкости будем определять по формуле [11]:

Ь = - — Р • Ш = -—РЕ • Ш,

V V

где Ш = (р, д, г) - вектор углового ускорения несущего тела, Р - коэффициент, учитывающий форму полости, р - плотность и V - кинематическая вязкость жидкости. Учитывая вышесказанное, равенство (4) можно переписать в виде:

т = -—Р%,

V

где т = (тх, ту, тг) ,

% =

р + дг - гд д + гр - рг г + рд - др

Равенства (1), (3) представляют собой динамические уравнения движения рассматриваемой системы.

Согласно [13] введём в уравнениях (1), (3) новую переменную:

5 = [(С - А)г + С2а] А-1.

(5)

С учётом новой переменной и результатов интегрирования равенства (3) уравнения движения (1) примут вид:

р + эд д - ф--

С1

С - А

рэ [С15 - С2 (го +ао )] ,

— Р

V АС - А р 1 Р

V АС~^А [С15 - С2 ( +ао)]' э = -—Рэ (р 2 + д2).

(6)

V А

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

2. СТАЦИОНАРНЫЕ РЕШЕНИЯ И ИССЛЕДОВАНИЕ ИХ УСТОЙЧИВОСТИ

Определим для данной системы множества стационарных движений, под которыми понимаются движения гиростата, для которых проекции

вектора угловой скорости на оси связанной системы координат постоянны: p = q = r = 0. Решая полученную систему уравнений, получаем два стационарных движения:

1) p = q = 0, ^ = u = const;

2) p = a = const, q = b = const, s = 0. (7)

Исследуем на устойчивость найденные стационарные движения. Составим уравнения возмущённого движения однороторного гиростата и уравнения их первого приближения вида x = Ax или в скалярном виде для первого и второго движений соответственно:

1)

x = Mx1 - Fx 2, x2 = Fx1 + Mx2,

•¿3 = 0,

2)

X1 Tx3,

X = Vx

2 r л3> (8)

JC3 Kx3,

где обозначено

E = M =

T = -b-

p P

7 AQ - A

p 1 PC2 a

7 ~AC - a

[C1u2 - C2U ( +^0 )] , ( + ^0 ),

p 1 PC2b , ,

V = a(0 + a0),

K =

7 AC-A pP1

(9)

7 AC

(a2 +b2 )(A-C1),

Е = -# = и.

Составляя характеристические уравнения линейных систем (8) и решая их относительно X , получим следующие корни:

1) X = 0, Яе(Х2) = Е, Яе(Х3) = Е;

2) X = 0, Х2 = 0, Х3 = К. (10)

При Яе(Х) = Е < 0 и Яе(Х) = К < 0, ввиду присутствия нулевых корней, согласно теореме Ляпунова об асимптотической устойчивости по первому приближению выводов об асимптотической устойчивости стационарных движений сделать нельзя, и данные движения необходимо исследовать другими методами теории устойчивости. При выполнении неравенств Яе(Х) = Е > 0 и Яе(Х) = К > 0 имеем неустойчивость таких движений [15].

3. СТАБИЛИЗАЦИЯ СТАЦИОНАРНЫХ ДВИЖЕНИЙ ГИРОСТАТА

Решим задачу о стабилизации второго неустойчивого стационарного решения (7) при условии, что Яе(X) = К > 0. Согласно [14] добавим в систему линейное по отклонениям управление вида и = Вх, и определим множества значений

элементов матрицы В , стабилизирующих устойчивые и неустойчивые стационарные движения гиростата до асимптотической устойчивости. Уравнения первого приближения управляемой системы, соответствующие уравнениям (10), примут вид X = Ах + Вх .

Выбор структуры матрицы В обусловлен двумя соображениями. Во-первых, возьмем наибольшее возможное число нулевых элементов для упрощения искомого управления. Во-вторых, матрица В должна удовлетворять критерию полной управляемости линейных автономных систем [15]. Выберем матрицу В следующим образом (легко убедиться, что критерий полной управляемости выполнен, то есть ранг соответствующей матрицы управляемости равняется размерности фазового пространства - трём):

(13)

Для того чтобы управляемое движение было асимптотически устойчивым, достаточно, чтобы все корни соответствующего характеристического уравнения

(¿11 -Х)[Х2 -Х(К + Ъ22) + [КЬ22 -УЬЪ2}] = 0,(15)

имели отрицательные вещественные части. Используя критерий Гурвица, определим области допустимых значений коэффициентов Ь11, Ь21, Ь33 управления, при которых достигается асимптотическая стабилизация. А именно, для устойчивости второй линейной системы (10) необходимо и достаточно, чтобы все главные миноры матрицы Гурвица были положительны. Таким образом, мы получаем систему трех неравенств:

Г ьи 0 01

B = 0 Ь22 0

V 0 Ь32 0,

bn < 0;

(-K - b22) > 0; (Kb22 - Vb32) > 0.

(18)

4. АНАЛИЗ МНОЖЕСТВА КОЭФФИЦИЕНТОВ

СТАБИЛИЗИРУЮЩЕГО УПРАВЛЕНИЯ

Найдем решения последних двух неравенств системы (18), учитывая следующие замечания:

- будем рассматривать любые значения величины K;

- величины p = a = const, q = b = const могут принимать значения любого знака, определяемого различными наборами значений параметров системы и начальных условий движения.

Решая систему второго и третьего неравенств из (18) относительно b11 и b21, получим совокупность трех множеств решений (19) - (21), удовлетворяющих критерию Гурвица:

1)

V > О; Ъ22 < -K; Kb22

Ъ32 <—22

32 V

(19)

Система (19) совместно с неравенством Ь11 < 0 задает целое множество значений коэффициентов Ъ11, Ъ21, Ъ33 для синтеза стабилизирующего управления. Несложный анализ, проведенный ниже, этого решения дает следующие возможности для наиболее простого выбора элементов матрицы управления В, упрощающие структуру полученного стабилизирующего управления. А именно,

A) При К > 0 можно выбрать произвольные значения

КЪ

Ъп < 0, Ъ22 < -К <0 и Ъз2 < КЪ2 < 0.

Б) При К < 0 можно при Ъ11 < 0 взять Ъ22 = 0, тогда Ъ32 < 0. И наоборот, если взять Ъ32 = 0, то Ъ22 < 0.

B) При К = 0 имеем Ъ22 < 0 и Ъ32 < 0.

Проведем анализ второй системы неравенств,

следуемой из критерия Гурвица

2)

V = О;

Ъ22 < -K ;

О < KK, VK,.

(20)

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

A) При К = 0 решения системы (20) не существует;

Б) При К > 0 последнее неравенство (20) несовместно;

B) При К < 0 можно взять Ъ22 < 0 и, в частности, Ъ32 = 0 .

Для третьего случая, при Ъ11 < 0, имеем систему неравенств:

3)

V < О; Ъ22 < -K; Kb22

к, >•—22

(21)

V

В частности, из нее вытекают простейшие решения:

А) При К = 0 Ъ22 < 0 , Ъ32 > 0 ;

Б) При К > 0 управлять необходимо обо-

ими коэффициентами, выбираемыми согласно (21);

В) При К < 0 можно взять коэффициент Ъ22 = 0 , тогда Ъ32 > 0 . И наоборот, если взять Ъ32 = 0 , то ему будет соответствовать любой

Ъ22 < О .

5. ЗАКЛЮЧЕНИЕ

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

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

1. Сазонов В.В. Периодические движения спутника-гиростата относительно центра масс под действием гравитационного момента // Космические исследования. 2013. Т. 51. № 2. С. 145-158.

2. Раушенбах В.В., Токарь В.И. Управление ориентацией космических аппаратов. М.: Наука, 1974. 589 с.

3. Летов А.М. Динамика полета и управление. М.: Наука, 1969. 359 с.

4. Румянцев В.В. Об устойчивости движения гиростатов // ПММ. 1961. Т. 25. Вып. 1. С. 9-16.

5. Маркеев А.П. Нелинейные колебания спутника при резонансе 1:1:1 // ПММ. 2012. Т. 76. Вып. 1. С. 52-68.

6. Харламов М.П., Рябов П.Е. Бифуркация первых интегралов в случае Ковалевской-Яхьи // Регулярная и хаотическая динамика. 1997. Т. 2. № 2. С. 25-40.

7. Hall C., Rand R. Spinap Dynamics of Axial Dual-Spin Spacecraft // J. Guidance Control Dyn. 1994. V.17. Issue 1. pp. 30-37.

8. El-Gohary A. Chaos and optimal control of steady-state rotation of a satellite-gyrostat on a circular orbit // Chaos, solitons & fractals. 2009. V. 42. Issue 5. Pp. 28422851.

9. Безгласный С.П., Мысина О.А. Об одноосной и трехосной ориентации составного тела // Вестник Сам-ГУ. Естественно-научная серия. 2011. № 2 (83). С. 80-90.

10. Безгласный С.П., Худякова М.А. Об ориентации гиростата с переменными моментами инерции // Автоматизация процессов управления. 2013. № 2 (32). С. 22-28.

11. Черноусько Ф.Л. Движение твердого тела с полостями, содержащими вязкую жидкость. М.: Изд. ВЦ АН СССР, 1968.

12. Моисеев Н.Н., Румянцев В.В. Динамика тела с полостями, содержащими жидкость. М.: Наука, 1965. 439 с.

13. Алексеев А.В. Движение спутника-гиростата, содержащего полость с жидкостью большой вязкости // Известия Самарского научного центра РАН. 2007. № 9.. С.671-676.

14. Алексеев А.В., Безгласный С.П., Красников В.С. Стабилизация стационарных движений однороторного гиростата с полостью, заполненной жидкостью большой вязкости // Вестник Самарского государственного аэрокосмического университета имени акаде-

мика С.П. Королева (национального исследовательского университета). 2012. № 6. С. 13-18.

15. Александров В.В., Болтянский В.Г., Лемак С.С., Парусников Н.А., Тихомиров В.М. Оптимальное управление движением. М.: ФИЗМАТЛИТ, 2005. 376 с.

BUILDING STABILIZING CONTROL FOR STATIONARY MOTIONS OF A GYROSTAT WITH A CAVITY VISCOUS FLUID

© 2013 A.V. Alekseev, S.P. Bezglasnyi, V.S. Krasnikov

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

Spherical motion relative to the center of mass of the free single rotor gyrostat with a spherical cavity is entirely filled with a fluid of high viscosity are researched. Equations of motion are found and investigated the stability of stationary motion. The set of controls on the basis of feedback which stabilize unstable and stable stationary motion of gyrostat to asymptotically stable, are constructed. Keywords: gyrostat, asymptotic stability, active control, feedback, stabilization, cavity with fluid.

Aleksey Alekseev, Candidate of Technical Sciences, Associate

Professor. E-mail: [email protected]

Sergey Bezglasnyi, Candidate of Physical and Mathematical

Sciences, Associate Professor, Doctoral Candidate.

E-mail: [email protected]

Viktor Krasnikov, Post Graduate Student.

E-mail: [email protected]