Уравнения движения планет. Решение задачи Кеплера в пакете Mathcad

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

где

- радиус-векторы первого и второго тела, соответственно,

- потенциальная взаимодействия тел, g - гравитационная постоянная. Введем вектор, направленный от первого тела ко второму телу

Тогда в системе отсчета с началом координат в центре масс рассматриваемой системы тел

Величину, определяемую в соответствии с (7), принято называть приведенной массой. Функция (6) формально совпадает с функцией Лагранжа одной материальной точки с массой m, движущейся в потенциале

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

, создаваемом неподвижным центром с массой m₁+m₂. Отметим, что если масса одного из взаимодействующих тел значительно меньше массы другого тела, последнее можно рассматривать как неподвижный притягивающий центр, и найденная зависимость

будет описывать траекторию движения более легкого тела.

В противном случае, решив задачу о движении тела с массой m в потенциале

по зависимости

, в соответствии с (4), (5) находят траектории каждой частицы

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

)

в полном соответствии с законом всемирного тяготения Ньютона.

Отметим два важных свойства силы тяготения, вытекающих из (10):

1) сила зависит только от расстояния между телами;

2) сила направлена по прямой, проходящей через центры взаимодействующих тел.

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

Сохранение момента импульса, в свою очередь, означает, что траектория движения тела в центральном поле лежит в плоскости, которой перпендикулярен вектор

. Кроме того, движение тела ограничивается условиями сохранения полной энергии

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

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

Введя обозначение

и сократив общие множители, запишем выражения (13), (14), составляющие систему ДУ второго порядка, в виде

Предваряя численное решение системы уравнений (15), (16), проведем обезразмеривание этих уравнений. Если в качестве единиц измерения расстояния и времени выбрать радиус орбиты R и период обращения Т, соответствующие движению тела по окружности, то можно ввести безразмерные переменные

. Выполнив в (15), (16) замену переменных x® X, y® Y, t® t , получаем:

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

и скоростью тела

соотношением

При движении в гравитационном поле по окружности центростремительное ускорение обусловлено гравитационной силой. Следовательно,

Выражение (21), являясь общим условием любой круговой орбиты, позволяет найти зависимость периода движения от радиуса орбиты. Период движения

поэтому, подставив в (22) выражение (21), получим

Подставляя выражение (23) в (17), (18), получаем окончательно обезразмеренную систему уравнений

Из уравнений (24), (25) видна их универсальность - они не зависят ни от периода обращения тела вокруг центра поля, ни от радиуса орбиты. Следовательно, величина

, входящая в (17) и (18), одинакова для всех тел, совершающих движение в гравитационном поле по замкнутым траекториям, что является доказательством справедливости третьего закона Кеплера.

При решении системы дифференциальных уравнений будем считать, что в начальный момент времени тело находилось в точке с радиус-вектором

, скорость тела была направлена вертикально вверх,

(рис. 2). Так как система уравнений (24), (25) является безразмерной, необходимо также привести к безразмерному виду начальные условия. Выполнив, как и выше, замену переменных

, приводим начальные условия к следующему виду:

Однако использовать конкретные числовые значения R, T, M для проверки законов Кеплера не требуется, так как безразмерные начальные условия также обладают известным универсализмом. Для того чтобы это показать, найдем безразмерную скорость тела, движущегося в гравитационном поле по окружности. Подставив (21), (23) в (27), получаем

Следовательно, для получения орбит, отличных от круговых, достаточно задавать значения начальной скорости, отличные от 2p.

Знаете ли Вы, что такое мысленный эксперимент, gedanken experiment?
Это несуществующая практика, потусторонний опыт, воображение того, чего нет на самом деле. Мысленные эксперименты подобны снам наяву. Они рождают чудовищ. В отличие от физического эксперимента, который является опытной проверкой гипотез, "мысленный эксперимент" фокуснически подменяет экспериментальную проверку желаемыми, не проверенными на практике выводами, манипулируя логикообразными построениями, реально нарушающими саму логику путем использования недоказанных посылок в качестве доказанных, то есть путем подмены. Таким образом, основной задачей заявителей "мысленных экспериментов" является обман слушателя или читателя путем замены настоящего физического эксперимента его "куклой" - фиктивными рассуждениями под честное слово без самой физической проверки.
Заполнение физики воображаемыми, "мысленными экспериментами" привело к возникновению абсурдной сюрреалистической, спутанно-запутанной картины мира. Настоящий исследователь должен отличать такие "фантики" от настоящих ценностей.

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

Это мы видим на примере СТО и ОТО, превратившихся в своеобразный вид религии, управляющей наукой и общественным мнением. Никакое количество фактов, противоречащих им, не может преодолеть формулу Эйнштейна: "Если факт не соответствует теории - измените факт" (В другом варианте " - Факт не соответствует теории? - Тем хуже для факта").

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

Эксперимент на то и эксперимент, что он есть не изощрение мысли, а проверка мысли. Непротиворечивая внутри себя мысль не может сама себя проверить. Это доказано Куртом Гёделем.

Понятие "мысленный эксперимент" придумано специально спекулянтами - релятивистами для шулерской подмены реальной проверки мысли на практике (эксперимента) своим "честным словом". Подробнее читайте в FAQ по эфирной физике.