УДК 52-323.8
М. Н. Бурдаев
Годографы и уравнение перелета в центральном поле тяготения
Аннотация. В статье приведен способ геометрического описания динамики космических полетов с помощью годографов. Представлен вывод уравнения годографа начальных скоростей перелета между двумя заданными точками центрального поля тяготения.
Ключевые слова и фразы: годограф, перелет, обучение, когнитивная графика.
Введение
Формы и методы, разрабатываемые в когнитивной графике, представляются перспективными для использования в процессах обучения различным специальностям и предметам. Основанием этого предположения служат такие свойства когнитивной графики как обобщенный характер представления информации, наглядность обобщенных форм ее отображения в графическом виде, возможность наглядно показывать на чертеже в обобщенном виде многомерные нелинейные зависимости, наличие у когнитивных изображений четкого физического смысла.
Одной из актуальных областей применения форм и методов когнитивной графики является динамика космического полета. По сложившейся традиции для описания движения космических объектов до настоящего времени используется аналитический математический
Работа выполнена при финансовой поддержке проекта РФФИ № 12—07— 00205—а «Разработка новых способов решения задач управления движениями космических аппаратов на всех этапах полетов и оперативного отображения получаемых результатов на основе методов годографов и когнитивной графики».
© М. Н. Бурдаев, 2012
© ФГБУ «Научно-исследовательский испытательный Центр подготовки космонавтов имени Ю. А. Гагарина», 2012
© Программные системы: теория и приложения, 2012
аппарат небесной механики. Его формулы не дают наглядного представления о взаимных связях параметров орбитальных движений космических объектов. Для успешной работы с этим аппаратом требуется развитое абстрактное мышление и отсутствие дефицита времени. Он эффективно применяется в научных исследованиях, в учебном процессе в высших учебных заведениях и в работе наземных служб навигационно-баллистического обеспечения космических полетов. Однако и в этих областях есть настоятельная потребность в повышении наглядности изложения теоретического материала, облегчении и повышении надежности восприятия информации, в разработке и применении оперативных методов решения прикладных задач. В учебном процессе подготовки космонавтов эти требования наиболее актуальны. Для их выполнения наиболее подходят формы и методы когнитивной графики. В динамике космического полета они нашли применение и развитие в виде специально разработанной отрасли теории годографов—системы обобщенных взаимно адекватных графических и аналитических форм представления предмета этой области знаний [1].
1. Годографы
Обобщенная форма уравнения конических сечений, отображающего траектории небесных тел, впервые появилась в теории годографов. Классическая форма этого уравнения
Р
г = -----------
1 + е cos V
имеет в правой части три независимых параметра: р (фокальный параметр орбиты), е (величина эксцентриситета траектории перелета) и $ (угол между местной вертикалью и направлением вектора скорости космического аппарата). Это создает значительные неудобства при попытках отображать на чертеже или экране дисплея более одной конкретной траектории космического объекта. Проблема наглядного отображения семейств траекторий была решена путем введения обобщенного уравнения конических сечений
Г_ = f = fx (е)
X 1 + cos $ 1 + cos $
В этом уравнении вместо х может стоять величина любой из характеристик эллипса, имеющих размерность длины: х = р, гп, га, а, Ь, с. Вид функций fx (е) = и ^ для различных х приведен в табл. 1.
Таблица 1. Вид функций fx (е) = ^ и ^ в обобщенном уравнении орбит
X р rn ra a b
fx (е) 1 1 + Є 1 — e 1 — e2 \/l — e2
Г l l+e 1-е l-e2 Vl-e2
X l+e cos 'д l+e cos'd l+e cos'd l+e cos'd l+e cos'd
Обобщенная форма уравнения конических сечений позволяет показать на чертеже семейства траекторий космических объектов, объединенных по какому-либо общему признаку. Совместное изображение семейства орбит, у которых радиус обобщен по какому-либо линейному параметру х, представляет собой сетку линий равных значений эксцентриситета е = Const в полярных координатах ^, $. Каждая линия в таком семействе отображает все орбиты с одинаковым эксцентриситетом, независимо от их размеров.
Например, для условия х = гп обобщенное уравнение орбит принимает вид
г 1 + е
гп 1 + е cos ё
Общий вид семейства орбит в этой обобщенной форме показан на рис. 1. На этом рисунке нанесены также оси прямоугольной горизонтальной системы обобщенных координат X и Y:
X = — cos $ ; Y = — sin $.
Линии равных значений параметров для краткости называют «изолиниями».
В тех же полярных координатах ^, $ можно отобразить изолинии любых параметров орбитального движения, которые контролируются и регулируются при пилотировании космического аппарата и решении навигационных задач.
Например, семейство изолиний истинных аномалий в полярных координатах ^, $ представлено пучком прямых линий, выходящих из начала координат.
Рис. 1. Семейство орбит, обобщенное по радиусу перигея
Совмещение в единой системе координат нескольких сеток изолиний позволяет наглядно представить количественные связи отображаемых параметров.
Графики, на которых совмещены два и более семейств изолиний, называют «годограммами».
Наиболее наглядно показать связи вектора орбитальной скорости с параметрами орбитального движения позволяет новая обобщенная форма годографов орбитальных скоростей. В этой форме рассматриваются все возможные векторы орбитальной скорости в произвольно выбранной точке поля тяготения.
Разделив величину вектора орбитальной скорости V в точке поля тяготения с геоцентрическим радиусом г
V = ^ — (е2 + 2е сов $ + 1)
на величину местной круговой скорости Укруг = \р^ на том же геоцентрическом радиусе г, получаем соотношение V:
V
1 + е
е + cos •&
^круг V 1 + Є cos $
Направление обобщенного вектора орбитальной скорости V совпадает с направлением вектора орбитальной скорости V и определяется углом ф, вычисляемым по формуле
е sin $
ctg Ф =—----------7 .
1 + е cos §
Два последних уравнения позволяют рассчитать и построить в орбитальной системе координат график, показанный на рис. 2.
250/1 | 280 290 300
Рис. 2. Годограмма обобщенных орбитальных скоростей в произвольной точке поля тяготения
На рис. 2 нанесены две сетки изолиний: е = Const и $ = Const. Эти сетки наглядно отображают связь между эксцентриситетом орбит е, истинной аномалией точки на орбите $, величиной V и направлением ф обобщенного вектора орбитальной скорости.
Каждой точке с истинной аномалией $ на орбите с эксцентриситетом е соответствует один обобщенный вектор орбитальной скорости с модулем V и углом к горизонту ф. И наоборот: каждой паре значений V и ф соответствует конкретная точка с истинной аномалией •&2 на орбите с эксцентриситетом Є2.
У этой обобщенной формы вектора орбитальной скорости есть полезное свойство: из нее очень просто определить величину орбитальной скорости в физической размерности. Для этого достаточно умножить величину обобщенной скорости на величину местной круговой скорости на том же геоцентрическом радиусе Укруг.
Эти свойства рис. 2 позволяют наглядно показать на нем процесс одноимпульсной коррекции орбиты и взаимные связи его параметров.
В тех же обобщенных координатах V и ф можно представить в обобщенных формах все геометрические параметры орбит.
На рис. 3 показана сетка изолиний обобщенного параметра радиуса перигея ^ = Const. В обобщенных координатах V и ф можно наглядно представить связь параметров орбитального движения в точке орбиты и времени полета от перицентра до рассматриваемой точки на обобщенных изображениях семейств орбит в виде сетки линий равных значений параметра времени Т = (рис. 4).
Рис. 3. Изолинии обобщенного радиуса перигея ^
Рис. 4. Изолинии обобщенного параметра времени Т
На основе показанных обобщенных векторов положения и скорости космических объектов разработаны два геометрических языка динамики космического полета. Для практики наиболее ценен язык обобщенных векторов скорости, непосредственно связывающий параметры орбитальных движений с вектором орбитальной скорости и его изменениями при выполняемых в космосе маневрах.
Далее представлен вывод уравнения годографа начальных скоростей перелета между двумя заданными точками центрального поля тяготения.
2. Уравнение перелета в центральном поле тяготения
Расчеты характеристик траекторий перелета между двумя точками центрального поля тяготения выполняются при решении многих задач практической космонавтики. Параметры траекторий перелета в поле тяготения с центром в точке О и гравитационной постоянной ц определяются пятью взаимно независимыми начальными условиями: величинами радиуса гм начальной точки М перелета,
радиуса гN конечной точки N перелета, угловой дальности Д$ перелета, величиной Ум скорости в начальной точке М перелета и углом фм, указывающим направление вектора этой скорости относительно местной вертикали (рис. 5).
Рис. 5. Начальные условия перелета между заданными точками центрального поля тяготения
При заданных величинах ^, гм, и Д$ параметры Ум и фм
начальных скоростей перехода могут непрерывно изменяться в некоторых пределах. Поэтому перелет между двумя заданными точками центрального поля тяготения возможен по бесконечному количеству траекторий, каждой из которых соответствует определенное сочетание величины Ум скорости и угла фм в начальной точке М перелета.
Семейство векторов начальной скорости Ум перелета для заданных начальных условий ц, гм, и Д$ представляет собой годограф
начальных скоростей перелета (рис. 6).
Выведем уравнение годографа начальных скоростей перелета. Решив совместно относительно величины эксцентриситета е уравнения конических сечений для граничных радиусов перелета гм и
Р
гм
(1)
получим
Гм =
1 + е 008 'вм р
(2)
1 + е 008 ('&М + ДІЇ)
1 _ 1М.
гы
( — 008 Д^ 008 $м + 8ІП $м 8ІП ДІЇ
е
Рис. б. Годограф начальных скоростей перелета
Подставим полученное выражение (2) в известное уравнение
^ + cos ём (1 - cos Д$) ctg +sinД^
и решим его относительно tg $м:
1 — cos Д-в
(3) tg'&M
(і - tgФм - sin Д$
г2 V2 . 2
Используя уравнение фокального параметра р = м sin фм и уравнение (І), получим равенство
(4) ^ ^ =1 + ecos^M.
^ (1 + ctg2 фм)
Получим другое представление. Подставим в правую часть ра-
венства (4) выражение (2) для эксцентриситета и выражение (3):
(1 - F7) cos^M____________
1 + е cos 'вм = 1 +
= 1 +
- cos Д^) cos 'дм + sin 'дм sin Д^
1 - ™ ) [(1 - ™) tg^M - sin Д^1
(5)
= 1 +
= 1 +
(tf - cos Д^) (1 - tg Фм + (1 - 7^) sin Дг& 1 - ™) К1 - ™) tg - sin н
(1 - - cos Д^) tg 'фм + sin Д$]
1 - tg-фм - sin Д-д^
^ - cos Д^ tg фм + sin Дг?]
(1 - cos/itf)tg^M _ УШУ
- cos^) tg 'Фм + sinAtf Л^Д^ +ctg фм
tg Дт
-cos Д$
-------------+ ctg фм
sin Д$
Знаменатель в окончательной форме выражения (5) имеет геометрический смысл.
В треугольнике OMN опустим перпендикуляр из вершины N на сторону OM (рис. Т). Из треугольника ONP определяем
NP = гм sin Дд,
в треугольнике MNP имеем
MP = гм - гм cos Д-в
и
NP = MP tg Дфм = (гм - гм cos Д$) tg (т - Дфм).
Приравнивая взаимно два полученных выражения для NP и решая это равенство относительно Д фм, находим
cos Д-в - —
Гы
(6) ctg ДФм = .до •
sm Др
Таким образом, обнаружено, что выражение (6) обозначает тригонометрическую функцию внешнего угла QMN треугольника OMN (рис. 6 и 7).
Рис. 7. К определению величины угла А'фм
Угол QMN = Дфм имеет физический смысл. Его сторона ММ указывает направление из начальной точки перелета М на его конечную точку N и соединяет эти граничные точки перелета. Продолжение стороны ММ является асимптотой годографа начальных скоростей перелета (рис. 6). Луч MQ направлен по местной вертикали в начальной точке М перелета и является исходным ориентиром для отсчета углов, указывающих направления векторов начальных скоростей перелета и второй асимптотой годографа начальных скоростей перелета. В пределах угла QMN лежит диапазон всех возможных значений начальных скоростей перелета по эллиптическим и гиперболическим траекториям. Границы этого диапазона будут определены ниже.
Подставляя (6) в (5) и (5) с этой подстановкой в (4), получаем уравнение годографа начальных скоростей перелета в форме
(7)
р (1+^2^м) ^фм - ^Дфм'
Представим это уравнение в обобщенном виде:
гмум „+„ Д^ _ (1 + ^2^м)
(8) — ___________==
р 2 ctg фм - Дфм '
В этой форме уравнения обобщенный вектор начальной скорости перелета гм^м ctg 4т является функцией двух взаимно независимых переменных: угла Д фм, указывающего направление из начальной на конечную точку перелета, и угла фм, определяющего направление начальной скорости перелета.
Уравнение (8) впервые позволило построить двумерное обобщенное изображение взаимных связей всех величин, характеризующих перелет между двумя заданными точками центрального поля тяготения (рис. 8).
Рис. 8. Семейство годографов обобщенного параметра ^ Ум начальной скорости перелета
Радиальные лучи на рис. 8 обозначают направления из начальной в конечную точку перелета и соответствуют постоянным значениям угла Афм. Изогнутые линии — годографы начальных скоростей перелетов для тех же значений углов Афм.
При необходимости, из уравнения (8) определяется величина начальной скорости перелета в физической размерности:
В пределах диапазона углов 0 < фм < Дфм существуют границы, разделяющие множества эллиптических и гиперболических орбит. Это углы фм, соответствующие параболическим траекториям.
Они определяются из условия гмУм = 2.
С использованием уравнения (8) это условие записывается в виде
2 _ гм Ум _ (1 + е^2 фм) Ад
М Фм - ^ Дфм 2 ’
откуда определяются граничные значения углов фм 1,2:
(9) фм 1,2 _ аге^
А$ / 2 А$ А'в
с^ё 2 ± у - 2^ — ^Афм - 1
или, после подстановки в (9) величины ctg Дфм из (6),
Дд
Фм 1,2 = агс^ё
Большее значение угла фм 1 определяет границу между эллиптическими и гиперболическими траекториями перелета. Меньшее значение угла фм2 ограничивает минимальное значение угла фм для эллиптических траекторий перелетов. Гиперболические траектории, соответствующие углам фм, меньшим угла фм2, уходят в бесконечность и не возвращаются в конечную точку перелета.
Максимальное значение угла фм для гиперболических траекторий теоретически равно величине угла Д фм. Этому случаю при конечных значениях величин начальных радиусов гм перелетов соответствуют бесконечные величины начальных скоростей перелетов, что противоречит законам физики. Поэтому максимальный возможный угол фм для перелетов по гиперболическим траекториям ограничен максимальной предельной скоростью перелета, равной скорости света.
Отметим, что скорости света в начальных точках перелетов с конечными значениями радиусов гм соответствуют по законам нерелятивистской механики гиперболические, следовательно, криволинейные траектории.
Заключение
Материал статьи является продолжением исследований, начатых автором в работе [2]. Полученные результаты используются в ФГБУ «Научно-исследовательский испытательный Центр подготовки космонавтов имени Ю. А. Гагарина» для обучения космонавтов теории перелетов.
Список литературы
[1] Бурдаев М. Н. Теория годографов в механике космического полета. М. : Машиностроение, 1975.— 151 c. f[]
[2] Бурдаев М. Н. Применение метода годографов к расчету времени перелета в центральном поле тяготения // Космические исследования, 2009. Т. 47, № 2, с. 204-208. f2
Рекомендовал к публикации д.т.н. В. М. Хачумов
Об авторе:
Михаил Николаевич Бурдаев
Доктор технических наук, профессор, академик Академии космонавтики, космонавт-испытатель, главный научный сотрудник 57 отдела 5-го Научного управления ФГБУ «Научно-исследовательский испытательный Центр подготовки космонавтов имени Ю. А. Гагарина».
e-mail: bmn@starcity.ru
Образец ссылки на эту публикацию:
М. Н. Бурдаев. Годографы и уравнение перелета в центральном поле тяготения // Программные системы: теория и приложения : электрон. научн. журн. 2012. T. 3, №3(12), с. 77-90.
URL: http://psta.psiras.ru/read/psta2012_3_77-90.pdf
M. N. Burdayev. Hodographs and the equation of flight in a central gravitational
field.
Abstract. The paper presents a method of geometrical description of space flight dynamics with the help of hodographs. Presented a derivation of the initial velocity hodograph flights between two given points of the central gravitational field. (In Russian).
Key Words and Phrases: hodograph, flight, teaching, cognitive graphics.