В небесной механике , А орбита Кеплер (или кеплеровская орбиты , названный в честь немецкого астронома Иоганна Кеплера ) является движение одного тела относительно другого, как эллипса , параболы или гиперболы , который образует двухмерную плоскость орбиты в трех- пространственное пространство. Орбита Кеплера также может образовывать прямую линию . Он учитывает только точечное гравитационное притяжение двух тел, пренебрегая возмущениями, вызванными гравитационным взаимодействием с другими объектами, атмосферным сопротивлением , давлением солнечного излучения и т. Д.сферическое центральное тело и т. д. Таким образом, говорят, что это решение частного случая проблемы двух тел , известной как проблема Кеплера . Как теория классической механики , она также не принимает во внимание эффекты общей теории относительности . Кеплеровские орбиты можно параметризовать на шесть орбитальных элементов различными способами.
В большинстве случаев используется большое центральное тело, центр масс которого считается центром масс всей системы. Путем разложения орбиты двух объектов одинаковой массы можно описать как орбиты Кеплера вокруг их общего центра масс, их барицентра .
Вступление
С древних времен до XVI и XVII веков считалось, что движения планет следуют идеально круговым геоцентрическим траекториям, как учили древнегреческие философы Аристотель и Птолемей . Вариации в движении планет объяснялись меньшими круговыми траекториями, наложенными на больший путь (см. Эпицикл ). По мере того, как измерения планет становились все более точными, были предложены поправки к теории. В 1543 году Николай Коперник опубликовал гелиоцентрическую модель Солнечной системы , хотя он все еще считал, что планеты движутся по идеально круговым траекториям с центром на Солнце. [1]
История Кеплера и телескопа
Кеплер переехал в Прагу и начал работать с Тихо Браге . Тихо поручил ему просмотреть всю информацию, которую Тихо имел о Марсе. Кеплер отметил, что положение Марса было подвержено большим ошибкам и создавало проблемы для многих моделей. Это побудило Кеплера сконфигурировать 3 закона движения планет.
Первый закон: планеты движутся по эллипсу с Солнцем в одном фокусе.
Закон изменил бы эксцентриситет на 0,0. и фокус больше с эксцентриситетом 0,8. которые показывают, что круговая и эллиптическая орбиты имеют одинаковый период и фокус, но разные размеры области, определяемой Солнцем.
Это приводит ко второму закону: радиус-вектор описывает равные площади в равное время.
Эти два закона были опубликованы в книге Кеплера Astronomia Nova в 1609 году.
Для кругового движения равномерно, однако для эллиптического, чтобы охватить область с равномерной скоростью, объект перемещается быстро, когда радиус-вектор короткий, и медленнее, когда радиус-вектор длинный.
Кеплер опубликовал свой Третий закон движения планет в 1619 году в своей книге « Гармоники мира» . Ньютон использовал третий закон для определения своих законов тяготения.
Третий закон: квадраты периодических времен относятся друг к другу как кубы средних расстояний. [2]
Развитие законов
В 1601 году Иоганн Кеплер получил обширные и тщательные наблюдения за планетами, сделанные Тихо Браге . Кеплер потратит следующие пять лет, пытаясь подогнать наблюдения планеты Марс к различным кривым. В 1609 году Кеплер опубликовал первые два из трех своих законов движения планет . Первый закон гласит:
В более общем смысле, путь объекта, претерпевающего кеплеровское движение, может также следовать параболе или гиперболе , которые, наряду с эллипсами, принадлежат к группе кривых, известных как конические секции . Математически расстояние между центральным телом и вращающимся телом можно выразить как:
где:
- это расстояние
- является большой полуосью , которая определяет размер орбиты
- это эксцентриситет , который определяет форму орбиты
- - истинная аномалия , представляющая собой угол между текущим положением орбитального объекта и местом на орбите, в котором он находится ближе всего к центральному телу (называемому перицентром ).
В качестве альтернативы уравнение может быть выражено как:
Где называется полу-латусной прямой кишкой дуги. Эта форма уравнения особенно полезна при работе с параболическими траекториями, у которых большая полуось бесконечна.
Несмотря на то, что эти законы были разработаны на основе наблюдений, Кеплер так и не смог разработать теорию, объясняющую эти движения. [3]
Исаак Ньютон
Между 1665 и 1666 годами Исаак Ньютон разработал несколько концепций, связанных с движением, гравитацией и дифференциальным исчислением. Однако эти концепции не были опубликованы до 1687 г. в « Началах» , в которых он изложил свои законы движения и закон всемирного тяготения . Его второй из трех законов движения гласит:
Ускорение тела параллельно и прямо пропорциональна чистой силы , действующей на тело, находится в направлении результирующей силы, и обратно пропорциональна массе тела:
Где:
- вектор силы
- масса тела, на которое действует сила
- - вектор ускорения, вторая производная по времени вектора положения
Строго говоря, эта форма уравнения применима только к объекту постоянной массы, что верно на основе упрощающих предположений, сделанных ниже.
Закон всемирного тяготения Ньютона гласит:
Каждая точечная масса притягивает любую другую точечную массу силой, направленной вдоль линии, пересекающей обе точки. Сила пропорциональна произведению двух масс и обратно пропорциональна квадрату расстояния между точечными массами:
где:
- величина гравитационной силы между двумя точечными массами
- является гравитационным постоянным
- масса первой точечной массы
- масса второй точечной массы
- расстояние между двумя точечными массами
Из законов движения и закона всемирного тяготения Ньютон смог вывести законы Кеплера, которые характерны для орбитального движения в астрономии. Поскольку законы Кеплера были хорошо подтверждены данными наблюдений, эта последовательность обеспечила сильную поддержку справедливости обобщенной теории Ньютона и единой небесной и обычной механики. Эти законы движения легли в основу современной небесной механики, пока Альберт Эйнштейн не представил концепции специальной и общей теории относительности в начале 20 века. Для большинства приложений кеплеровское движение приближает движения планет и спутников с относительно высокой степенью точности и широко используется в астрономии и астродинамике .
Упрощенная проблема двух тел
- См. Также анализ орбиты.
Чтобы решить вопрос о движении объекта в системе двух тел , можно сделать два упрощающих предположения:
- 1. Тела сферически симметричны и могут рассматриваться как точечные массы.
- 2. На тела не действуют никакие внешние или внутренние силы, кроме их взаимного тяготения.
По форме крупные небесные тела близки к сферам. По симметрии чистая гравитационная сила, притягивающая материальную точку к однородной сфере, должна быть направлена к ее центру. Теорема оболочек (также доказанная Исааком Ньютоном) утверждает, что величина этой силы такая же, как если бы вся масса была сосредоточена в середине сферы, даже если плотность сферы меняется с глубиной (как это происходит для большинства небесных тел). тела). Из этого немедленно следует, что притяжение между двумя однородными сферами происходит так, как если бы масса обоих была сосредоточена в их центре.
Более мелкие объекты, такие как астероиды или космические корабли, часто имеют форму, сильно отличающуюся от сферы. Но гравитационные силы, создаваемые этими неоднородностями, обычно невелики по сравнению с силой тяжести центрального тела. Разница между неправильной формой и идеальной сферой также уменьшается с увеличением расстояния, и большинство орбитальных расстояний очень велико по сравнению с диаметром небольшого орбитального тела. Таким образом, для некоторых приложений неравномерностью формы можно пренебречь без значительного влияния на точность. Этот эффект весьма заметен для искусственных спутников Земли, особенно на низких орбитах.
Планеты вращаются с разной скоростью и поэтому могут принимать слегка сжатую форму из-за центробежной силы. При такой сжатой форме гравитационное притяжение будет несколько отличаться от притяжения однородной сферы. На больших расстояниях эффект сжатия становится незначительным. Движение планет в Солнечной системе можно вычислить с достаточной точностью, если рассматривать их как точечные массы.
Объекты с двумя точечными массами с массами а также и векторы положения а также относительно некоторой инерциальной системы отсчета действуют гравитационные силы:
где - вектор относительного положения массы 1 по отношению к массе 2, выраженный как:
а также является единичным вектором в этом направлении и- длина этого вектора.
Разделив на их соответствующие массы и вычтя второе уравнение из первого, получим уравнение движения для ускорения первого объекта относительно второго:
( 1 )
где - гравитационный параметр и равен
Во многих приложениях можно сделать третье упрощающее предположение:
- 3. По сравнению с центральным телом масса вращающегося тела незначительна. Математически m 1 >> m 2 , поэтому α = G ( m 1 + m 2 ) ≈ Gm 1 .
Это предположение не является необходимым для решения упрощенной задачи двух тел, но оно упрощает вычисления, особенно со спутниками, вращающимися вокруг Земли, и планетами, вращающимися вокруг Солнца. Даже масса Юпитера меньше массы Солнца в 1047 раз [4], что составляет ошибку 0,096% в значении α. Известные исключения включают систему Земля-Луна (отношение масс 81,3), систему Плутон-Харон (отношение масс 8,9) и двойные звездные системы.
При этих предположениях дифференциальное уравнение для случая двух тел может быть полностью решено математически, и результирующая орбита, которая следует законам движения планет Кеплера, называется «орбитой Кеплера». Орбиты всех планет с высокой точностью соответствуют орбитам Кеплера вокруг Солнца. Небольшие отклонения вызваны гораздо более слабым гравитационным притяжением между планетами, а в случае Меркурия - общей теорией относительности . Орбиты искусственных спутников вокруг Земли в хорошем приближении являются орбитами Кеплера с небольшими возмущениями из-за гравитационного притяжения Солнца, Луны и сжатия Земли. В приложениях с высокой точностью, для которых уравнение движения должно быть интегрировано численно с учетом всех гравитационных и негравитационных сил (таких как давление солнечного излучения и атмосферное сопротивление ), концепции орбиты Кеплера имеют первостепенное значение и широко используются.
Кеплеровские элементы
Любую кеплеровскую траекторию можно определить шестью параметрами. Движение объекта, движущегося в трехмерном пространстве, характеризуется вектором положения и вектором скорости. Каждый вектор состоит из трех компонентов, поэтому общее количество значений, необходимых для определения траектории в пространстве, равно шести. Орбита обычно определяется шестью элементами (известными как кеплеровские элементы ), которые можно вычислить по положению и скорости, три из которых уже обсуждались. Эти элементы удобны тем, что их шесть, пять неизменны для невозмущенной орбиты (резкий контраст с двумя постоянно меняющимися векторами). Будущее местоположение объекта в пределах его орбиты может быть предсказано, а его новое положение и скорость могут быть легко получены с помощью элементов орбиты.
Два определяют размер и форму траектории:
- Большая полуось ()
- Эксцентриситет ()
Три определяют ориентацию плоскости орбиты :
- Наклон () определяет угол между плоскостью орбиты и плоскостью отсчета.
- Долгота восходящего узла () определяет угол между опорным направлением и восходящим пересечением орбиты на опорной плоскости (восходящий узел).
- Аргумент перицентра () определяет угол между восходящим узлом и перицентром.
И наконец:
- Истинная аномалия () определяет положение орбитального тела вдоль траектории, отсчитываемой от перицентра. Вместо истинной аномалии можно использовать несколько альтернативных значений, наиболее распространенным из которых являетсясредняя аномалия и, время, прошедшее с периапсиса.
Так как , а также представляют собой просто угловые измерения, определяющие ориентацию траектории в системе отсчета, они не являются строго необходимыми при обсуждении движения объекта в плоскости орбиты. Они были упомянуты здесь для полноты, но не требуются для приведенных ниже доказательств.
Математическое решение дифференциального уравнения ( 1 ) выше
Для движения под действием любой центральной силы, т. Е. Силы, параллельной r , удельный относительный угловой момент остается постоянным:
Поскольку векторное произведение вектора положения и его скорости остается постоянным, они должны лежать в одной плоскости, ортогональной к . Это означает, что вектор-функция представляет собой плоскую кривую .
Поскольку уравнение симметрично относительно его источника, его легче решить в полярных координатах. Однако важно отметить, что уравнение ( 1 ) относится к линейному ускорению. в отличие от угловых или радиальный ускорение. Поэтому при преобразовании уравнения нужно быть осторожным. Введение в декартову систему координати полярные единичные векторы в плоскости, ортогональной :
Теперь мы можем переписать вектор-функцию и его производные как:
(см. « Векторное исчисление »). Подставляя их в ( 1 ), находим:
Это дает необычное полярное дифференциальное уравнение:
( 2 )
Чтобы решить это уравнение, необходимо исключить все производные по времени. Это приносит:
( 3 )
Взяв производную по времени от ( 3 ), получаем
( 4 )
Уравнения ( 3 ) и ( 4 ) позволяют исключить производные по времени от. Чтобы исключить производные по времени от, цепное правило используется для поиска подходящих замен:
( 5 )
( 6 )
Используя эти четыре подстановки, можно исключить все производные по времени в ( 2 ), получив обыкновенное дифференциальное уравнение для как функция
( 7 )
Дифференциальное уравнение ( 7 ) можно решить аналитически путем подстановки переменной
( 8 )
Используя цепное правило для дифференциации, получаем:
( 9 )
( 10 )
Используя выражения ( 10 ) и ( 9 ) для а также получает
( 11 )
с общим решением
( 12 )
где е и- константы интегрирования, зависящие от начальных значений s и
Вместо использования постоянной интегрирования явно вводится соглашение о том, что единичные векторы определяющие систему координат в плоскости орбиты, выбираются так, чтобы принимает значение ноль, а е положительно. Это значит, что равен нулю в точке, где является максимальным и, следовательно, минимально. Определив параметр p как у одного есть это
Альтернативное происхождение
Другой способ решить это уравнение без использования полярных дифференциальных уравнений:
Определите единичный вектор такой, что а также . Следует, что
Теперь рассмотрим
(см. Векторное тройное произведение ). Заметь
Подстановка этих значений в предыдущее уравнение дает:
Объединение обеих сторон:
где c - постоянный вектор. Обозначение этого символа r дает интересный результат:
где угол между а также . Решение для r :
Заметь фактически являются полярными координатами вектор-функции. Делаем замены а также , мы снова приходим к уравнению
( 13 )
Это уравнение в полярных координатах для конического сечения с началом в фокальной точке. Аргумент называется «истинной аномалией».
Свойства уравнения траектории
Для это круг радиуса p .
Для это эллипс с
( 14 )
( 15 )
Для это парабола с фокусным расстоянием
Для это гипербола с
( 16 )
( 17 )
На следующем изображении показаны круг (серый), эллипс (красный), парабола (зеленый) и гипербола (синий).
Точка на горизонтальной линии, идущей вправо от точки фокусировки, - это точка с при котором расстояние до фокуса принимает минимальное значение перицентр. Для эллипса также существует апоцентр, для которого расстояние до фокуса принимает максимальное значение Для гиперболы диапазон для является
а для параболы диапазон равен
Используя цепное правило дифференцирования ( 5 ), уравнение ( 2 ) и определение p как получается, что радиальная составляющая скорости равна
( 18 )
и что тангенциальная составляющая (составляющая скорости, перпендикулярная к ) является
( 19 )
Связь полярного аргумента а время t немного отличается для эллиптических и гиперболических орбит.
Для эллиптической орбиты переходят к « эксцентрической аномалии » E, для которой
( 20 )
( 21 )
и следовательно
( 22 )
( 23 )
а угловой момент H равен
( 24 )
Интегрирование по времени t дает
( 25 )
в предположении, что время выбирается так, чтобы постоянная интегрирования была равна нулю.
Поскольку по определению p имеется
( 26 )
это можно написать
( 27 )
Для гиперболической орбиты используются гиперболические функции для параметризации
( 28 )
( 29 )
для чего есть
( 30 )
( 31 )
а угловой момент H равен
( 32 )
Интегрируя по времени t, получаем
( 33 )
т.е.
( 34 )
Чтобы узнать, какое время t соответствует определенной истинной аномалии вычисляется соответствующий параметр E, связанный со временем с помощью соотношения ( 27 ) для эллиптической орбиты и с соотношением ( 34 ) для гиперболической орбиты.
Отметим, что соотношения ( 27 ) и ( 34 ) определяют отображение между диапазонами
Некоторые дополнительные формулы
Для эллиптической орбиты из ( 20 ) и ( 21 ) получаем
( 35 )
и поэтому
( 36 )
Тогда из ( 36 ) следует, что
Из геометрической конструкции, определяющей эксцентрическую аномалию, ясно, что векторы а также находятся по одну сторону от оси x . Отсюда следует, что векторы а также находятся в одном квадранте. Следовательно, есть
( 37 )
и это
( 38 )
( 39 )
где ""- полярный аргумент вектора и n выбрано так, чтобы
Для численного расчета может использоваться стандартная функция ATAN2 (y, x) (или DATAN2 (y, x) с двойной точностью ), доступная, например, в языке программирования FORTRAN .
Обратите внимание, что это отображение между диапазонами
Для гиперболической орбиты из ( 28 ) и ( 29 ) получаем
( 40 )
и поэтому
( 41 )
В виде
и, как а также того же знака следует, что
( 42 )
Это соотношение удобно для перехода между «истинной аномалией» и параметром E , который связан со временем соотношением ( 34 ). Обратите внимание, что это отображение между диапазонами
и это можно вычислить, используя соотношение
Из соотношения ( 27 ) следует, что орбитальный период P для эллиптической орбиты равен
( 43 )
Поскольку потенциальная энергия, соответствующая силовому полю соотношения ( 1 ), равна
из ( 13 ), ( 14 ), ( 18 ) и ( 19 ) следует, что сумма кинетической и потенциальной энергии
для эллиптической орбиты
( 44 )
и из ( 13 ), ( 16 ), ( 18 ) и ( 19 ), что сумма кинетической и потенциальной энергии для гиперболической орбиты равна
( 45 )
Относительная инерциальная система координат
в орбитальной плоскости с по направлению к перицентру из ( 18 ) и ( 19 ) получаем, что компоненты скорости равны
( 46 )
( 47 )
См. Также Уравнение центра - Аналитические разложения.
Уравнение центра связывает среднюю аномалию с истинной аномалией для эллиптических орбит при небольшом числовом эксцентриситете.
Определение кеплеровской орбиты, соответствующей заданному начальному состоянию
Это « задача начального значения » для дифференциального уравнения ( 1 ), которое является уравнением первого порядка для 6-мерного «вектора состояния». когда написано как
( 48 )
( 49 )
При любых значениях исходного «вектора состояния» Кеплеровская орбита, соответствующая решению этой начальной задачи, может быть найдена с помощью следующего алгоритма:
Определите ортогональные единичные векторы через
( 50 )
( 51 )
с участием а также
Из ( 13 ), ( 18 ) и ( 19 ) следует, что, полагая
( 52 )
и определив а также такой, что
( 53 )
( 54 )
где
( 55 )
получается орбита Кеплера, которая для истинной аномалии имеет то же r , а также значения, определенные формулами ( 50 ) и ( 51 ).
Если эта орбита Кеплера, то также имеет такой же векторы для этой истинной аномалии как определенные ( 50 ) и ( 51 ) вектор состояния орбиты Кеплера принимает желаемые значения для истинной аномалии .
Стандартная инерционно неподвижная система координат в плоскости орбиты (с направленный от центра однородной сферы к перицентру), определяющий ориентацию конического сечения (эллипс, парабола или гипербола), затем можно определить с помощью соотношения
( 56 )
( 57 )
Отметим, что соотношения ( 53 ) и ( 54 ) имеют особенность при а также
т.е.
( 58 )
это тот случай, когда это круговая орбита, которая соответствует начальному состоянию
Оскулирующая орбита Кеплера
Для любого вектора состояния орбита Кеплера, соответствующая этому состоянию, может быть вычислена с помощью алгоритма, определенного выше. Сначала параметры определяются из а затем ортогональные единичные векторы в плоскости орбиты используя соотношения ( 56 ) и ( 57 ).
Если теперь уравнение движения имеет вид
( 59 )
где
это функция, отличная от
итоговые параметры
определяется все будет меняться со временем, в отличие от орбиты Кеплера, для которой только параметр будет отличаться
Вычисленная таким образом орбита Кеплера, имеющая тот же «вектор состояния», что и решение «уравнения движения» ( 59 ) в момент времени t , называется в это время «соприкасающейся».
Эта концепция, например, полезна в случае, если
где
представляет собой небольшую «возмущающую силу», вызванную, например, слабым гравитационным притяжением других небесных тел. В этом случае параметры оскулирующей орбиты Кеплера будут только медленно изменяться, и оскулирующая орбита Кеплера является хорошим приближением к реальной орбите в течение значительного периода времени до и после времени осколивания.
Эта концепция также может быть полезна для ракеты во время полета с двигателем, поскольку она затем сообщает, по какой орбите Кеплера ракета продолжит движение в случае отключения тяги.
Для орбиты, близкой к круговой, понятие « вектор эксцентриситета » определяется какявляется полезным. Из ( 53 ), ( 54 ) и ( 56 ) следует, что
( 60 )
т.е. - гладкая дифференцируемая функция вектора состояния также, если это состояние соответствует круговой орбите.
Смотрите также
- Проблема двух тел
- Гравитационная задача двух тел
- Проблема Кеплера
- Законы движения планет Кеплера
- Эллиптическая орбита
- Гиперболическая траектория
- Параболическая траектория
- Радиальная траектория
- Моделирование орбиты
Цитаты
- ^ Коперник. стр. 513–514
- ^ Gould, Алан (2016-09-24). «Иоганн Кеплер: его жизнь, его законы и времена» . НАСА . Проверено 3 декабря 2018 .
- ↑ Бейт, Мюллер, Уайт. стр 177–181
- ^ http://ssd.jpl.nasa.gov
Рекомендации
- Эль-Ясберг "Теория полета искусственных спутников Земли", Израильская программа научных переводов (1967)
- Бейт, Роджер; Мюллер, Дональд; Белый, Джерри (1971). Основы астродинамики . Dover Publications, Inc., Нью-Йорк. ISBN 0-486-60061-0.
- Коперник, Николай (1952), «Книга I, глава 4, Движение небесных тел является регулярным, круговым и вечным - или еще одним составным из круговых движений», О вращении небесных сфер , Великие книги западного мира , 16 , перевод Чарльза Гленна Уоллиса, Чикаго: Уильям Бентон, стр. 497–838.
Внешние ссылки
- Апплет JAVA, анимирующий орбиту спутника по эллиптической орбите Кеплера вокруг Земли с любым значением большой полуоси и эксцентриситета.