MathCAD   КМ   ТПОИ   3GL   к экономической информатике   4GL - визуальным средам

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

Поршнев С.В.

Уравнения движения планет

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

, (1)

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

. (2)

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

. (3)

Из (2), (3) находим:

, (4)

. (5)

Подставляя (4), (5) в (1), получаем

, (6)

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

. (7)

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

В противном случае, решив задачу о движении тела с массой m в потенциале по зависимости , в соответствии с (4), (5) находят траектории каждой частицы .

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

, (8)

получим уравнение движения тела

, (9)

которое при m1>> m2 принимает вид

(10)

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

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

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

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

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

, (11)

где .

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

(12а)

и величины

. (12б)

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

pic0002.jpg (4662 bytes)

Рис.1

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

, (13)

. (14)

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

, (15)

. (16)

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

, (17)

. (18)

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

. (19)

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

, (20)

откуда находим

. (21)

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

, (22)

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

. (23)

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

, (24)

. (25)

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

pic0003.jpg (3778 bytes)

Рис.2

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

, (26)

, (27)

где T определяется выражением (23).

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

. (28)

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

на верх   MathCAD   КМ   ТПОИ   3GL   к экономической информатике   4GL - визуальным средам

Знаете ли Вы, в чем ложность понятия "физический вакуум"?

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

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

Понять абсурдность этого понятия легче всего обратившись к истокам его возникновения. Рождено оно было Полем Дираком в 1930-х, когда стало ясно, что отрицание эфира в чистом виде, как это делал великий математик, но посредственный физик Анри Пуанкаре, уже нельзя. Слишком много фактов противоречит этому.

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

Однако такая постановка является внутренне противоречивой (виртуальные частицы ненаблюдаемы и их по произволу можно считать в одном случае отсутствующими, а в другом - присутствующими) и противоречащей релятивизму (то есть отрицанию эфира, так как при наличии таких частиц в вакууме релятивизм уже просто невозможен). Подробнее читайте в FAQ по эфирной физике.

НОВОСТИ ФОРУМА

Форум Рыцари теории эфира


Рыцари теории эфира
 10.11.2021 - 12:37: ПЕРСОНАЛИИ - Personalias -> WHO IS WHO - КТО ЕСТЬ КТО - Карим_Хайдаров.
10.11.2021 - 12:36: СОВЕСТЬ - Conscience -> РАСЧЕЛОВЕЧИВАНИЕ ЧЕЛОВЕКА. КОМУ ЭТО НАДО? - Карим_Хайдаров.
10.11.2021 - 12:36: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от д.м.н. Александра Алексеевича Редько - Карим_Хайдаров.
10.11.2021 - 12:35: ЭКОЛОГИЯ - Ecology -> Биологическая безопасность населения - Карим_Хайдаров.
10.11.2021 - 12:34: ВОЙНА, ПОЛИТИКА И НАУКА - War, Politics and Science -> Проблема государственного терроризма - Карим_Хайдаров.
10.11.2021 - 12:34: ВОЙНА, ПОЛИТИКА И НАУКА - War, Politics and Science -> ПРАВОСУДИЯ.НЕТ - Карим_Хайдаров.
10.11.2021 - 12:34: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от Вадима Глогера, США - Карим_Хайдаров.
10.11.2021 - 09:18: НОВЫЕ ТЕХНОЛОГИИ - New Technologies -> Волновая генетика Петра Гаряева, 5G-контроль и управление - Карим_Хайдаров.
10.11.2021 - 09:18: ЭКОЛОГИЯ - Ecology -> ЭКОЛОГИЯ ДЛЯ ВСЕХ - Карим_Хайдаров.
10.11.2021 - 09:16: ЭКОЛОГИЯ - Ecology -> ПРОБЛЕМЫ МЕДИЦИНЫ - Карим_Хайдаров.
10.11.2021 - 09:15: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от Екатерины Коваленко - Карим_Хайдаров.
10.11.2021 - 09:13: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от Вильгельма Варкентина - Карим_Хайдаров.
Bourabai Research - Технологии XXI века Bourabai Research Institution