Исследование движения разгонного блока с тороидальным баком относительно центра масс в параметрах Кэли-Клейна Текст научной статьи по специальности «Физика»

УДК 629.78 + 531.36

ИССЛЕДОВАНИЕ ДВИЖЕНИЯ РАЗГОННОГО БЛОКА С ТОРОИДАЛЬНЫМ БАКОМ ОТНОСИТЕЛЬНО ЦЕНТРА МАСС В ПАРАМЕТРАХ КЭЛИ-КЛЕЙНА

© 2011 А. В. Алексеев, В. С. Красников

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

Предложена математическая модель движения вокруг центра масс разгонного блока (РБ) с тороидальным баком, заполненным жидким топливом. Определены зависимости параметров движения от времени численным и аналитическим методами. Исследованы предельные случаи ориентации в углах Эйлера.

Разгонный блок, параметры Кэли-Клейна, углы

Введение

В ряде задач динамики полёта космических аппаратов необходимо определить движение относительно центра масс. Традиционно для определения ориентации твёрдого тела в пространстве применяются углы Эйлера. Однако, определить аналитические зависимости указанных углов от времени в общем случае трудно. Поэтому используют в основном их численные зависимости, полученные с помощью численного интегрирования соответствующих уравнений движения. Целью данной работы является определение аналитических зависимостей от времени традиционных углов ориентации твёрдого тела с полостью, целиком заполненной жидкостью. Твёрдым телом с тороидальной полостью, целиком заполненной жидкостью, моделируется разгонный блок (РБ). Форма полости выбрана неслучайно: это один из наиболее распространённых типов баков наряду с цилиндрическими и сферическими, используемыми в РБ.

Под простейшими случаями движения твёрдого тела с полостями, наполненными жидкостью, будем подразумевать случаи, когда движение жидкости в полости можно полностью охарактеризовать конечным числом переменных [1]. Очевидно, это возможно лишь при полном заполнении жидкостью полости, когда отсутствует свободная поверхность. Известно, что если при этом движение жидкости является потенциальным или же однородным вихревым, то движение си-

Эйлера, кинематические уравнения.

стемы описывается обыкновенными дифференциальными уравнениями вместе с уравнением Лапласа [2, 3]. Исследование движения системы в этих случаях значительно упрощается.

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

— + ю X О = Ь . (1)

—

Суммарный момент количеств движения твёрдого тела и жидкости в его полости равен:

О = О1 + О* + я,

где вектор Я представляет собой линейную функцию главных циркуляций и является постоянным в связанной системе координат.

Уравнение движения системы относительно центра масс в векторном выражении имеет вид

— (О1 + О*) + юх(О1 + О* + Я) = Ь .

—

Такой же вид имеет уравнение движения системы около неподвижной точки. Если оси абсолютной системы координат направить по главным осям эллипсоида

r(0* + 0(1))r = 1 , то векторное уравнение (1) в проекциях на оси жёстко связанной с телом системы координат даёт следующие три скалярных обыкновенных дифференциальных уравнения:

Ap + (C - B)qr + R3q - R2r = L;

< Bq + (A - C)pr + R1r - R3p = M; (2)

Cr + (B - A)pq + R2p - R1q = N.

Здесь A, B, C - главные моменты преобразованного тела для точки О, равные суммам моментов инерции твёрдого тела и эквивалентного тела, относительно соответствующих осей координат, направленных по главным осям эллипсоида инерции для точки О; L, M, N - проекции внешнего момента на связанные оси координат; p, q, r - проекции вектора мгновенной угловой скорости w на

оси связанной системы координат; R1, R2, R3

- проекции вектора R момента количества циклического движения жидкости, которые определяются для всех кольцевидных полостей в форме тел вращения вокруг оси х3 через потенциал циклического движения жидкости [1]:

R = R = о, R3 = Mp. (3)

2p

где M2 - масса жидкости, c - циркуляция жидкости.

Таким образом, полная математическая модель, позволяющая получить зависимости параметров движения от времени, состоит из системы (2) и любых кинематических уравнений, например, кинематических уравнений Эйлера:

p =y& sin0 sin j + 0 cos j;

<q = y& sin0cosj - 0 sinj; (4)

r = цг cos0 + j,

где y, 0, j - углы Эйлера. Одним из вариантов исследования движения системы является получение решения системы дифференциальных уравнений (2), (4) в углах Эй-

лера. Нелинейность уравнений (4) не позволяет получить зависимости углов от времени в аналитическом виде. На рис. 1, 3, 5 представлены зависимости углов Эйлера от времени, полученные в результате численного интегрирования системы дифференциальных уравнений (2), (4). На рис. 2, 4, 6 для анализа влияния жидкости на движение вокруг центра масс показаны те же углы, полученные аналогичным образом, но для тела без жидкости.

Решения, полученные в результате численного интегрирования, полезны для определения характера движения и его параметров для определённого набора начальных условий и характеристик системы. Проведение по ним качественного анализа движения и выявление влияния различных параметров системы затруднительно и требует больших затрат времени и ресурсов. Удобнее проводить анализ по аналитическим зависимостям.

2. Получение аналитических зависимостей параметров движения от времени

Для решения поставленной задачи введём следующие допущения, характерные для РБ и позволяющие упростить решение динамических уравнений.

1. Будем считать, что движение происходит под действием только тяги маршевых двигательных установок, которая не создаёт момента относительно центра масс. Тогда проекции внешнего момента на оси связной системы координат равны нулю:

N = М = Ь = 0.

2. В качестве топливных баков используются баки тороидальной формы, т.е. для них проекции момента циклического движения жидкости на оси связной системы координат имеют вид (3).

3. РБ является динамически симметричным:

А = В Ф С.

4. Топливный бак расположен перпендикулярно продольной оси симметрии РБ.

Перепишем систему (2) с учётом введённых допущений:

Рис. 1. Зависимость угла прецессии от времени для тела с жидкостью

Рис. 2. Зависимость угла прецессии от времени для тела без жидкости

Рис. 3. Зависимость угла нутации от времени для тела с жидкостью

Рис. 4. Зависимость угла нутации от времени для тела без жидкости

Рис. 5. Зависимость угла собственного вращения Рис. 6. Зависимость угла собственного вращения

от времени для тела с жидкостью от времени для тела без жидкости

Ap + (C - B)qr + R3q = G; Bq + (A - C)pr - R3p = G; Cr = G.

(З)

(7)

rI

р = - - pI+-(p+qI )a; 22 rI

y=- yI+-(p - qI )S; 22

rI S=- - dI + -(p + qI )y. 22

Подставляя в систему дифференциальных уравнений первого порядка (7) аналитические выражения для угловых скоростей (6) и решая задачу Коши для этой системы с

начальными условиями а0, Р0, 70, ^0, после введения замены

m = V- k 2 - rGG- 2К - qG- pG

получим следующие аналитические выражения для параметров Кэли-Клейна:

Из (5) получим аналитические зависимости от времени проекций вектора угловой скорости на оси связанной системы координат:

p(t) = p0 cos kt - q0 sin kt;

q(t) = p0 sin kt + q0 cos kt; (6)

r (t) = r^

где p0, q0, r0 - начальные условия,

k = — [r0(C - A) + R3]-циклическая частота. A

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

Для определения ориентации РБ будем использовать параметры Кэли-Клейна, в которых кинематические уравнения запишутся в виде

rI а =—aI +—(p - qI)P; 22

- ^( Ik ~№)t - ^( Ik+m)t

a(t) = Ag e 2 - BG e 2 ;

- ^( Ik +m)t - ^( Ik-m)t

CGe 2 + DGe 2

(8)

P(t) =

(-pg + Iq g )e

-Ikt

где константы т, А0, В0, С0, Б0 определяются из начальных условий системы.

Аналитические зависимости для 5(1), у(1) являются комплексным сопряжением а(Х), ((^ и, наоборот, с точностью до знака [4]:

У(!) = -((0;

5^) = Щ/).

Для определения ориентации РБ в традиционных углах Эйлера необходимо воспользоваться их связью с параметрами Кэли-Клейна [5]:

в = arccos(aS + Py);

Г Г ay ^ Re 7

PS.

(9)

y = — arccos

2

Re

yS

Графическое представление решения (9) полностью совпадает с графиками на рис. 1, 3, 5, что подтверждает адекватность полученных зависимостей.

Таким образом, получены аналитические зависимости углов Эйлера от времени, которые определить из кинематических соотношений Эйлера весьма сложно, а в некоторых случаях и невозможно. С помощью полученных соотношений (9) можно анализировать влияние параметров системы, в том числе и жидкости, на угловое движение РБ.

l

3. Предельные случаи ориентации в углах Эйлера

Как известно, при 9=0 или 9=п углы у, р не определены, а определена только их сумма у+р. Эта особенность углов Эйлера делает их малопригодными для исследования движения вблизи указанных значений угла нутации. В этом случае необходимо использовать другие углы. Однако, если использовать решение (9), полученное выше, то эта проблема устраняется. Например, численное решение для нулевого начального угла нутации получить невозможно из-за возникновения деления на ноль в кинематических уравнениях Эйлера. График зависимости угла нутации от времени, полученной при начальных условиях 9(0) = 0, р(0) = 0, у(0) = 0, приведён на рисунке 7.

Заключение

На основании работ Н. Е. Жуковского построена математическая модель движения вокруг центра масс осесимметричного РБ с тороидальным баком, заполненным жидким топливом. Для анализа влияния топлива на движение системы проведено численное интегрирование уравнений движения. Для исследования предельных случаев движения в углах Эйлера и качественной оценки влияния параметров системы на её движение с помощью параметров Кэли-Клейна построено аналитическое решение уравнений движения.

Данная статья написана по результатам проведения поисковой научно-исследовательской работы в рамках реализации ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009 - 2013 годы.

Библиографический список

1. Жуковский, Н. Е. О движении твердого тела, имеющего полости, наполненные однородной капельной жидкостью [Текст] / Н. Е. Жуковский. - М.: Гостехиздат, 1949.766 с.

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

3. Алексеев, А. В. Приведение спутни-

ка-гиростата с полостью с жидкостью к системам твердых тел с вязким трением [Текст] / А. В. Алексеев, А. В. Дорошин // Общероссийский научно-технический журнал «Полёт». - 2007. - № 9. - С. 26-33.

4. Лурье, А. И. Аналитическая механика [Текст] / А. И. Лурье. - М.:ГИФМЛ, 1961.825 с.

5. Голдстейн, Г. Классическая механика: пер. с англ., 2 изд.[Текст] / Г. Голдстейн. -М., 1975. - 413 с.

MOTION OF AN UPPER STAGE ROCKET WITH A TOROIDAL TANK AROUND THE MASS CENTER IN PARAMETERS OF CAYLEY-KLEIN

© 2011 A. V. Alekseev, V. S. Krasnikov

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

The paper presents a mathematical model of motion of an upper-stage rocket with a toroidal tank filled with liquid fuel around the centre of mass. Time dependencies of motion parameters are established using a numerical and an analytical method. Limiting cases of orientation in Eulerian angles are analyzed.

Eulerian angles, upper-stage rocket, kinematic equations, Cayley-Klein parameters.

Информация об авторах

Алексеев Алексей Владимирович, кандидат технических наук, доцент кафедры теоретической механики, Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет). E-mail: [email protected]. Область научных интересов: динамика систем твёрдых тел, динамика систем переменного состава, движение твёрдых тел с жидкостью.

Красников Виктор Сергеевич, магистрант, Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет). E-mail: [email protected]. Область научных интересов: движение твёрдых тел с жидкостью.

Alekseev Aleksey Vladimirovich, candidate of technical sciences, associate professor, Samara State Aerospace University, e-mail: [email protected]. Area of research: dynamics of rigid bodies, dynamics of variable composition systems, motion of tanks with liquids.

Krasnikov Victor Sergeevich, student of Samara State Aerospace University, e-mail: [email protected]. Area of research: motion of tanks with liquids.