The investigation of the motion of planets, the Moon, and the Sun based on a new principle of interaction

Full Text

Введение. Закон всемирного тяготения был сформулирован Ньютоном с целью разработки математической теории движения небесных тел. На основе закона тяготения Ньютона была разработана небесная механика, которая до середины XIX века вполне удовлетворительно объясняла движение небесных тел в Солнечной системе. В начале 1850 г. французским ученым Леверье (U.J.J. Le Verrier) была доказана невозможность представления наблюдений прохождения Меркурия по диску Солнца на основе ньютоновской динамики любой системой оскулирующих элементов и масс известных ему планет [1, 2]. В конце XIX столетия, наряду с установленными эмпирическими трудностями, стали все более выявляться и обсуждаться проблемы, связанные с феноменологическим характером ньютоновской теории. Закон всемирного тяготения Ньютона, лежащий в основе небесной механики, не связан с какимлибо конкретным представлением о природе или механизме взаимного притяжения между телами. Различные затруднения вызывали у ученых сомнения в точности закона Ньютона и делались многочисленные попытки внести поправки в формулу закона тяготения. Эти попытки улучшения закона всемирного тяготения оставались безуспешными, и лишь создание общей теории относительности, казалось бы, устранило все имеющие трудности, связанные с прогнозированием движения небесных тел, объяснив при этом невязки теоретических расчетов с наблюдениями в смещении перигелия Меркурия. Однако общая теория относительности не свободна от недостатков. Важнейшим из них является вопрос о природе гравитации. В рамках общей теории относительности он так же, как и в теории гравитации Ньютона, рассматривается чисто феноменологически. Другим недостатком общей теории относительности является существенное усложнение дифференциальных уравнений в задаче n тел. Решение этой задачи приходится искать в виде рядов по степеням малых параметров, при этом учитываются лишь начальные члены разложения и совершенно не исследуется вопрос о сходимости рядов. Основным недостатком как ньютоновской динамики, так и теории тяготения Эйнштейна является, по мнению автора, наделение массы свойством, порождающим поле тяготения. В первом случае предполагается, что масса обладает свойством притяжения других материальных тел, во втором случае масса наделяется свойством искривлять окружающее пространство. Все вышеперечисленные недостатки как ньютоновской динамики, так и общей теории относительности отражаются на степени точности и достоверности исследований, проводимых на их основе. 1. Вывод дифференциальных уравнений движения материальных тел, основанных на новом принципе взаимодействия. Возникает вопрос: возможно ли создание математической модели, в которой присутствует другое объяснение природы гравитации. На основании многочисленных фундаментальных исследований, проведенных ранее, можно ответить, что в рамках существующих моделей этот вопрос не может быть решен положительно. По-видимому, первопричину гравитации следует искать не в наличии массы в материи, а в свойствах взаимодействия окружающего пространства с движущейся ма119 З а у с а е в А. Ф. терией. Под окружающим пространством можно понимать физический вакуум, свойства которого в настоящее время пока находятся в стадии формирования. Физический вакуум стал предметом изучения физики благодаря трудам известных ученых: П. Дирака, Р. Фейнмана, Де Ситтера, Д. Казимира, Я. Б. Зельдовича и др. [3-6]. Ученые считают, что физический вакуум является особым состоянием материи, претендующим на первооснову мира. Физический вакуум генетически предшествует полю и веществу, он порождает их, поэтому вся Вселенная живет по его законам, которые науке предстоит еще открыть. На основании построения простейшей модели, описывающей взаимодействие небесных тел с окружающим пространством, автором данной статьи получены дифференциальные уравнения движения n материальных тел в барицентрической системе координат. Для вывода этих уравнений допускается ряд упрощений: пространство, в котором движутся материальные тела, имеющие сферическую форму с равномерно распределенной плотностью, обладает свойством сжатия при изменении пространственной плотности. В любой заданный момент времени каждая фиксированная область пространства имеет свою определенную плотность. Очевидно, что в пространственных областях, где находятся материальные тела, плотность пространства выше по сравнению с окрестностью, окружающей это тело. При движении материальных тел происходит изменение состояния всей системы, в местах, где находились материальные тела, образуются пустоты, т. е. области, плотности которых значительно меньше плотности окружающей среды. Реакция пространства на образование пустот должна приводить к выравниванию плотностей или заполнению пустот окружающей средой. Таким образом, пространство сжимается, заполняя пустоты, образованные движущимися материальными телами. Вывод дифференциальных уравнений движения основан на следующей идее. В каждый фиксированный момент времени материальное тело занимает в пространстве определенный объем. При перемещении тела пространство, занимаемое им в предыдущий момент времени, освобождается. При этом освободившееся пространство заполняется окружающей его средой, тем самым происходит сжатие окружающего пространства на величину объема, освобожденного движущемся объектом. Различные материальные тела занимают в пространстве определенные объемы. Однако вытесняемый ими объем, как правило, не равен фактическому объему и находится в прямой зависимости от плотности материального тела, т. е. для всякого материального тела существует предельная плотность, для которой можно рассчитать вытесняемый этим телом объем. Для сферически симметричных тел этот объем можно вычислить по формуле 4 3 V0 = πr0 . 3 где r0 назовем эффективным радиусом материального тела, который равен радиусу этого тела, сжатого до предельной плотности. На произвольном расстоянии r от центра материального тела пространство будет сжиматься на величину объема 4 3 4 3 4 3 πr - πr1 = πr0 3 3 3 120 (1) Исследование движения планет, Луны и Солнца, основанное на новом принципе взаимодействия или, сокращая соотношение (1) на постоянный множитель, получим 3 3 r3 - r1 = r0 , (2) где r1 - неизвестное расстояние, подлежащее определению. Из выражения (2) найдем r1 по формуле r1 = 3 3 r3 - r0 . (3) Из соотношения (2) находим величину сжатия на произвольном расстоянии r: 3 r0 (4) r - r1 = 2 2. r + rr1 + r1 С учетом (3) выражение (4) запишем в виде 3 r0 r - r1 = . 3 3 r2 + r 3 r3 - r0 + 3 (r3 - r0 )2 Определим теперь T0 (время сжатия пространства от r = r0 до r1 = 0), полагая, что r функционально зависит от ускорения w(r) и времени t по закону r3 w(r)t2 r= 0 = . (5) r2 2 Из (5) следует, что ускорение за единицу времени обратно пропорционально квадрату расстояния. Выразим эту зависимость соотношением w(r) = µ , r2 где µ - постоянная величина, значение которой находим из соотношения 2 µ = r0 w(r0 ). Для нахождения T0 , используем соотношение dr = w(r)t. dt Представим его в виде µt dr = 2. dt r Интегрируя последнее выражение, получим 3 2 r0 µT0 = . 3 2 (6) 3 2r0 . 3µ (7) Выразив T0 из (6), находим 2 T0 = 121 З а у с а е в А. Ф. Из (1), (2) следует, что относительно материального тела сжатое пространство можно представить в виде последовательности концентрических сфер, объемы которых равны объему данного тела. Допуская, что сжатие пространства от r до r1 между каждой последовательной парой сфер происходит за одинаковое время, получим 3 r0 r2 + r 3 r3 - r0 + 3 3 3 (r3 - r0 )2 = 2 w(r)T0 . 2 (8) Из выражений (7) и (8) находим w(r) = 3µ r2 3 +r r3 3 - r0 + 3 3 (r3 - r0 )2 . (9) Таким образом, получено соотношение для нахождения ускорения на произвольном расстоянии от сферически симметричного (с равномерно распределенной плотностью) материального тела. Дифференциальные уравнения движения в классической ньютоновской задаче n тел в прямоугольных координатах с началом в центре масс всей системы n материальных точек имеют следующий вид [7, 8]: d2 X = dt2 Xi - X Yi - Y d2 Y = k 2 mi , , 3 2 dt ∆i ∆3 i i d2 Z Zi - Z = , k 2 mi dt2 ∆3 i k 2 mi i (10) i где ∆2 = (Xi - X)2 + (Yi - Y )2 + (Zi - Z)2 . Здесь через X, Y , Z обозначены i барицентрические координаты возмущаемого тела, а через mi , Xi , Yi , Zi - массы и барицентрические координаты возмущающих тел. Для получения дифференциальных уравнений движения, основанных на новом принципе взаимодействия, в барицентрической системе координат достаточно умножить ускорение возмущаемого тела на направляющие косинусы данного ускорения. С учетом (9), дифференциальные уравнения движения возмущаемого тела в барицентрической системе координат примут следующий вид [9, 10]: d2 X = dt2 d2 Y = dt2 d2 Z = dt2 i 2 Xi - X 3w0i r0i , ∆i 2 + ∆ 3 ∆3 - r 3 + 3 (∆3 - r 3 )2 ∆i i i 0i i 0i 2 Yi - Y 3w0i r0i , ∆i ∆2 + ∆ 3 ∆3 - r3 + 3 (∆3 - r3 )2 i Zi - Z ∆i ∆2 + ∆ i i i i i i 0i i 2 3w0i r0i 3 3 ∆3 - r0i + i (11) 0i , 3 3 (∆3 - r0i )2 i где r0i - эффективный радиус i-того тела, w0i - соответствующее ускорение для i-того тела на расстоянии r0i от центра массы. 122 Исследование движения планет, Луны и Солнца, основанное на новом принципе взаимодействия Дифференциальные уравнения (11) имеют более сложный вид по сравнению с уравнениями (10). Однако они значительно проще дифференциальных уравнений движения, учитывающих релятивистские эффекты [11-13]. Под эффективным радиусом материального тела будем понимать такой радиус сферы, которому соответствует сжатие пространства на величину фактического объема. Из физических соображений ясно, что тела различной плотности и одинакового размера не могут иметь равные эффективные радиусы. Их значения можно определить путем согласования решения уравнений (11) с наблюдениями. Следует отметить, что система дифференциальных уравнений (11) не содержит масс тел и силовых взаимодействий. Дифференциальные уравнения (11) сохраняют принцип дальнодействия, так как в этих уравнениях отсутствуют силовые взаимодействия между телами и процесс сжатия пространства относительно всех материальных тел происходит одновременно. Таким образом, дифференциальные уравнения (10), (11) описывают движение материальных тел в барицентрической системе координат. Уравнения (10) - обычные дифференциальные уравнения движения в ньютоновской форме. Вывод дифференциальных уравнений (11) основан на новом принципе взаимодействия материальных тел, в которых отсутствует явно понятие массы и силы. Понятие ускорения определяется через закон изменения величины радиус-вектора при сжатии пространства в процессе движения материального тела. При r0i → 0, т. е. при замене материальных тел конечных размеров материальными телами с бесконечно малыми линейными размерами, уравнения (11) обращаются в уравнения (10). 2. Численное интегрирование уравнений движения больших планет, Луны и Солнца на интервале времени с 1600 по 2200 гг. В настоящее время одной из высокоточных численных теорий движения больших планет является планетная теория DE405/LE405, созданная американскими учеными Ньюхалом, Стендишем, Вильямсом (X. X. Newhall, E. M. Standish, Jr., and J. G. Williams) [11, 14], на интервале времени с 2305424.5 J.D. (1599 12 05) до 2525008.5 J.D. (2201 02 20). Полученные ими эфемериды для внутренних планет полностью согласованы с оптическими и радиолокационными наблюдениями. С целью проверки эффективности различных математических моделей, описывающих движение больших планет, автором данной работы проведено сопоставление результатов численного интегрирования дифференциальных уравнений (10), (11) на интервале времени с 1600 по 2200 гг. с DE405/LE405. ˙ ˙ ˙ Начальные данные координат X, Y , Z и скоростей X, Y , Z, фактических радиусов планет, Луны и Солнца r и ускорений w приведены в табл. 1. В правую часть уравнений (11) в качестве постоянных входят эффективные радиусы r0i и ускорения w0i планет и Солнца, значения которых зависят от плотности тел. Очевидно, чем меньше плотность небесного тела, тем больше различаются его эффективный и фактический радиус. На основании тестовых расчетов путем согласования решения уравнений (11) с данными, полученными на основании DE405/LE405, были найдены эффективные радиусы и ускорения для всех планет Луны и Солнца. Численное интегрирование уравнений движения (10), (11) было проведено модифицированным методом Эверхарта 27 порядка с шагом интегрирования 3 дня [16-18]. Результаты численного интегрирования представлены в табл. 2-9, в ко123 З а у с а е в А. Ф. Таблица 1 Начальные данные (координаты (а.е.), скорости (a.е./сут), радиус (км) и ускорение (км/c2 )) для планет, Луны и Солнца на момент времени T = 2440400.5 J.D. (1969 06 28.0) [The initial data (coordinates (au), speeds (au/day), radius (km), and acceleration (km/sec2 )) for the planets, the Moon and the Sun on the time T = 2440400.5 J.D. (1969 06 28.0)] Меркурий [Mercury] X = 0.361762714604 Y = -0.090781967730 Z = -0.085714983182 ˙ ˙ ˙ X = 0.00336749391398 Y = 0.02489452044680 Z = 0.01294630068860 r = 2.424 · 103 w = 0.369066656 · 10-2 Венера [Venus] X = 0.612751941342 Y = -0.348365368495 Z = -0.195278288980 ˙ ˙ ˙ X = 0.0109520683617 Y = 0.0156176843653 Z = 0.00633110555360 3 -2 r = 6.100 · 10 w = 0.899027143 · 10 Земля [Earth] X = 0.120527237123 Y = -0.925814243017 Z = -0.401527009924 ˙ ˙ ˙ X = 0.01680396477149 Y = 0.00175034387379 Z = 0.000759242499157 3 -2 r = 6.20315 · 10 w = 0.979831333 · 10 Луна [Moon] X = 0.119719059795 Y = -0.927808873018 Z = -0.402614272585 ˙ = 0.01740504958809 ˙ = 0.00158289841317 ˙ X Y Z = 0.00067368035418 r = 1.738 · 103 w = 0.163 · 10-2 Марс [Mars] X = -0.110186074283 Y = -1.327599456133 Z = -0.605889132614 ˙ ˙ ˙ X = 0.0144816530597 Y = 0.00024246311776 Z = -0.000281520734247 3 -2 r = 4.812 · 10 w = 0.469169773 · 10 Юпитер [Jupiter] X = -5.379706898836 Y = -0.830480581460 Z = -0.224828700228 ˙ ˙ ˙ X = 0.00109201154301 Y = -0.00651811656579 Z = -0.00282078316536 r = 7.144 · 104 w = 2.484123914 · 10-2 Сатурн [Saturn] X = 7.894392441979 Y = 4.596477801627 Z = 1.558697573530 ˙ ˙ ˙ X = -0.00321755523930 Y = 0.00433580985896 Z = 0.00192864656566 r = 6.044 · 104 w = 1.038017148 · 10-2 Уран [Uranus] X = -18.265398306822 Y = -1.161944505518 Z = -0.250103483937 ˙ ˙ ˙ X = 0.000221188417418 Y = -0.00376247593285 Z = -0.00165101470307 4 -2 r = 2.486 · 10 w = 0.939014231 · 10 Нептун [Neptune] X = -16.055042583768 Y = -23.942181216179 Z = -9.400156723549 ˙ ˙ ˙ X = 0.00264277104336 Y = -0.00149831445536 Z = -0.000679041903018 r = 2.650 · 104 w = 1.000508651 · 10-2 Плутон [Pluto] X = -30.483319603999 Y = -0.872478355496 Z = 8.911563040990 ˙ ˙ ˙ X = 0.000322210447723 Y = -0.00314357030215 Z = -0.00107794882974 3 -2 r = 2.900 · 10 w = 0.1167183115 · 10 Солнце [Sun] X = 0.0045025081562339 Y = 0.0007670747009324 Z = 0.0002660568051770 ˙ ˙ ˙ X = -0.3517482096 · 10-6 Y = 0.51776253996 · 10-5 Z = 0.22291018544 · 10-5 r = 6.9596405726 · 105 w = 273.9920358 · 10-3 124 Исследование движения планет, Луны и Солнца, основанное на новом принципе взаимодействия Таблица 2 Элементы орбит внутренних планет, полученные по данным DE405/LE405 на дату 1602 09 04 [The elements of the orbits of the inner planets, obtained according to the DE405/LE405 on the date 1602 09 04] Object M a e ω Ω i Mercury Venus Earth + Moon Mars 293.1782 108.2540 246.8564 293.1929 0.387099 0.723331 1.000004 1.523694 0.205553 0.006953 0.016900 0.092958 27.9534 53.8602 106.2233 283.4517 48.8676 77.8905 355.2989 50.8750 7.0244 3.3961 0.0456 1.8777 M -mean anomaly (degrees); a - orbit major semiaxis (au); e - orbit eccentricity; ω - orbit ascending node-perihelion angle (degrees); Ω - longitude of orbit ascending node (degrees); i - orbit inclination to the ecliptic plane (degrees) Таблица 3 Расхождения в элементах орбит численного интегрирования уравнений движения (10) и (11) с данными табл. 2 на дату 1602 09 04 [Differences in the orbital elements of the numerical integration of the motion equations (10) and (11) with the data of table 2 on the date 1602 09 04] Object ∆M ∆a ∆e ∆ω ∆Ω ∆i Mercury / Equ. (11) Mercury / Equ. (10) -0.0003 0.0137 -0.000001 -0.000001 0 0 0.0006 0.0443 -0.0002 -0.0001 0 0 Venus / Equ. (11) Venus / Equ. (10) -0.0018 -0.0010 0 0 0 0 0.0070 0.0092 0 0 0 0 Earth + Moon / Equ. (11) Earth + Moon / Equ. (10) 0.0001 0.0002 0 0 0 0 0.0042 0.0048 -0.0009 -0.0010 0 0 Mars / Equ. (11) Mars / Equ. (10) 0.0001 0.0001 0 0 0 0 0.0012 0.0013 0 0 0 0 Таблица 4 Элементы орбит внутренних планет, полученные по данным DE405/LE405 на дату 1702 09 24 [The elements of the orbits of the inner planets, obtained according to the DE405/LE405 on the date 1702 09 24] Object M a e ω Ω i Mercury Venus Earth + Moon Mars 83.4455 336.6162 264.5779 2.9941 0.387097 0.723327 1.000020 1.523726 0.205573 0.006946 0.016821 0.093150 28.2393 54.0330 107.0936 284.2193 48.7423 77.6097 354.8071 50.5817 7.0184 3.3955 0.0323 1.8697 Таблица 5 Расхождения в элементах орбит численного интегрирования уравнений движения (10) и (11) с данными табл. 4 на дату 1702 09 24 [Differences in the orbital elements of the numerical integration of the motion equations (10) and (11) with the data of table 4 on the date 1702 09 24] Object ∆M ∆a ∆e ∆ω ∆Ω ∆i Mercury / Equ. (11) Mercury / Equ. (10) 0 0.0002 0 0.000001 0 0 0.0001 0.0319 0 0.0001 0 0 Venus / Equ. (11) Venus / Equ. (10) -0.0009 -0.0004 0 0 0 0 0.0046 0.0063 0 0 0 0 Earth + Moon / Equ. (11) Earth + Moon / Equ. (10) 0.0001 0.0002 0 0 0 0 0.0030 0.0028 -0.0010 -0.0010 0 0 Mars / Equ. (11) Mars / Equ. (10) 0 0.0001 0 0 0 0 -0.0011 -0.0010 0 0 0 0 125 З а у с а е в А. Ф. Таблица 6 Элементы орбит внутренних планет, полученные по данным DE405/LE405 на дату 2069 06 28 [The elements of the orbits of the inner planets, obtained according to the DE405/LE405 on the date 2069 06 28] Object M a e ω Ω i Mercury Venus Earth + Moon Mars 0.2858 33.0207 172.6113 359.2689 0.387099 0.723332 0.999991 1.523765 0.205631 0.006779 0.016691 0.093561 29.2845 55.2393 286.0928 286.8867 48.2828 76.5939 177.0134 49.5031 6.9965 3.3926 0.0155 1.8398 Таблица 7 Расхождения в элементах орбит численного интегрирования уравнений движения (10) и (11) с данными табл. 6 на дату 2069 06 28 [Differences in the orbital elements of the numerical integration of the motion equations (10) and (11) with the data of table 6 on the date 2069 06 28] Object ∆M ∆a ∆e ∆ω ∆Ω ∆i -6 Mercury / Equ. (11) Mercury / Equ. (10) -0.0001 -0.0039 0 0 10 0 0 -0.0019 0 0.0001 0 0 Venus / Equ. (11) Venus / Equ. (10) 0.0001 -0.0001 0 0 0 0 -0.0014 -0.0021 0 0 0 0 Earth + Moon / Equ. (11) Earth + Moon / Equ. (10) 0 0 0 0 0 0 -0.0001 -0.0003 -0.0008 -0.0007 0 0 Mars / Equ. (11) Mars / Equ. (10) 0 0 10-6 10-6 0 0 -0.0003 -0.0004 0 0 0 0 Таблица 8 Элементы орбит внутренних планет, полученные по данным DE405/LE405 на дату 2200 04 10 [The elements of the orbits of the inner planets, obtained according to the DE405/LE405 on the date 2200 04 10] Object M a e ω Ω i Mercury Venus Earth + Moon Mars 358.8282 241.9619 91.7078 190.2019 0.387099 0.723335 1.000000 1.523747 0.205658 0.0066686 0.016646 0.093485 29.6552 55.6271 288.0876 287.8232 48.1184 76.2302 175.5105 49.1155 6.9887 3.3912 0.0325 1.8291 Таблица 9 Расхождения в элементах орбит численного интегрирования уравнений движения (10) и (11) с данными табл. 8 на дату 2200 04 10 [Differences in the orbital elements of the numerical integration of the motion equations (10) and (11) with the data of table 8 on the date 2200 04 10] Object ∆a ∆e ∆ω ∆Ω ∆i Mercury / Equ. (11) Mercury / Equ. (10) 0 0 0 -10-6 0 0 0.0001 -0.0275 -0.0001 -0.0001 0 0 Venus / Equ. (11) Venus / Equ. (10) -0.0010 0.0005 0 0 0 0 -0.0043 -0.0057 0 0 0 0 Earth + Moon / Equ. (11) Earth + Moon / Equ. (10) -0.0001 -0.0002 0 0 0 0 -0.0011 -0.0015 -0.0008 -0.0008 0 0 Mars / Equ. (11) Mars / Equ. (10) 126 ∆M -0.0001 -0.0001 0 0 10-6 10-6 -0.0008 -0.0009 -0.0001 -0.0001 0 0 Исследование движения планет, Луны и Солнца, основанное на новом принципе взаимодействия торых приняты следующие обозначения: M - средняя аномалия; ω - аргумент перигелия; Ω - долгота восходящего узла; i - наклонение; ∆M , ∆a, ∆e, ∆ω, ∆Ω, ∆i - отклонения элементов орбит, вычисленных на основе решения уравнений (10) и (11), от элементов, найденных по данным DE405/LE405. Все эти величины выражены в градусах и долях градусов; большая полуось a - в астрономических единицах. В табл. 2, 4, 6, 8 приведены элементы орбит больших планет на четыре момента оскуляции, полученные с помощью использования банка данных DE405/LE405 [14]. В табл. 3, 5, 7, 9 показаны расхождения в элементах орбит численного интегрирования уравнений движения (10) и (11) с данными табл. 2, 4, 6, 8 соответственно. При этом данные первой строки табл. 3, 5, 7, 9 найдены путем решения дифференциальных уравнений (11), полученных на основе нового принципа взаимодействия, а данные второй строки получены путем решения классических ньютоновских дифференциальных уравнений (10). Из сопоставления орбитальных элементов, приведенных в табл. 3, 5, 7, следует, что элементы орбит, найденные путем решения дифференциальных уравнений (11) и полученные с помощью использования банка данных DE405/LE405, отличаются незначительно. Максимальное различие в элементах орбит наблюдается в аргументах перигелиев в 1602 г. Для Меркурия расхождения в значениях аргумента перигелия составляет 0.0006◦ , что соответствует отклонению векового смещения от значений, полученных на основе банка данных DE405/LE405 - 0.59 . Для Венеры вековое отклонение составляет 6.87 , для барицентра Земля + Луна - 4.12 и для Марса - 1.18 . Эти расхождения связаны с различной скоростью векового движения аргументов перигелиев у различных планет. В табл. 10 представлены вековые смещения долгот перигелиев внутренних планет, вычисленных на основании решения дифференциальных уравнений (10) и полученных с использованием банка данных DE405/LE405. Как видно из табл. 10, вековое смещение долготы перигелия Меркурия в результате учета релятивистских эффектов и смещения полученного на основании решения уравнений (11) отличаются друг от друга на 0.3 , что находится в пределах вычислительной погрешности. Скорости движения перигелиев Венеры, Земли и Марса отличаются друг от друга более значительно. Однако вследствие того, что погрешности смещений перигелиев внутренних планет составляют для Меркурия ±0.44 , для Венеры, Земли и Марса - ±5.28 , ±1.79 и ±0.27 соответственно [13], нет возможности однозначно ответить на вопрос, какие смещения перигелиев планет, приведенные в табл. 10 для внутренних планет, являются истинными. Отличия остальных элементов орбит: больших полуосей a, эксцентриситетов e, долгот восходящих узлов Ω и наклонений i (см. табл. 2-9) у всех Таблица 10 Вековое движение перигелиев внутренних планет [The secular movement of the perihelion of the inner planets] Solution Mercury Venus Earth + Moon Mars DE405/LE405 Equ. (11) 43.25 ± 0.44 42.95 8.06 ± 5.28 2.18 5.01 ± 1.79 0.60 1.07 ± 0.27 0.12 127 З а у с а е в А. Ф. планет на всем интервале интегрирования с 1602 по 2200 гг. незначительные. Сравнение элементов орбит внутренних планет, вычисленных с использованием дифференциальных уравнений (11) и DE405/LE405 (см. табл. 2- 9), указывает на вполне удовлетворительное согласие, так как максимальная разность между элементами орбит находится в пределах вычислительных погрешностей. Возникает естественный вопрос: выдвигались ли подобные гипотезы ранее, делались ли предположения, что причины проявления гравитации являются неотъемлемым свойством окружающего пространства, а не массы сосредоточенной в теле? Подобные объяснения причины гравитации были высказаны Б. Риманом [19], а впоследствии А. Пуанкаре [20]. Рассуждая о природе тяготения, Б. Риман пишет [19, c. 472]: «Существующую (. . . ) в каждой точке пространства определенную по величине и направлению силу ускорения я пытаюсь объяснить движением некоей субстанции, наполняющей все бесконечное пространство. Эту субстанцию можно представить себе как физическое пространство, точки которого движутся в геометрическом пространстве. На основании этого допущения все воздействия весомых тел на весомые тела передаются в пустом пространстве посредством названной субстанции. Таким образом, формы движения, лежащие в существе света и теплоты, посылаемых небесными телами, суть не что иное, как формы движения этой субстанции. Но названные явления, именно тяготение и распространение света сквозь пустое пространство, - единственные, которые должны были бы быть объяснены только движением этой субстанции». А. Пуанкаре высказывает схожие с Б. Риманом суждения о причине тяготения. Говоря о тяготении, А. Пуанкаре пишет [20, c. 503]: «В то время как в ньютоновской механике количество энергии движущегося тела зависит от инерции тела, находящегося в движении, здесь энергия зависит от того, что называют инерцией эфира по отношению к электромагнитным силам. Инерция эфира возрастает вместе со скоростью и становится бесконечно большой, когда скорость электрона приближается к скорости света. Таким образом, кажущаяся масса электрона возрастает со скоростью (. . . ). При этом новом представлении постоянной массы материи не существует. Инерцией обладает не материя, а эфир; он один оказывает сопротивление движению, так что можно было бы сказать: нет материи, есть только дыры в эфире». Таким образом, размышляя о природе тяготения, как Б. Риман, так и А. Пуанкаре приходят к тому, что причину тяготения следует искать не в наличии массы в материальном теле, а в свойстве бесконечного пространства, которое у Б. Римана называется движущейся субстанцией, а у А. Пуанкаре свойством эфира. В предложенной нами модели считается, что изменением движения материальных тел в каждый момент времени управляет состояние пространства. Относительно любой точки пространства можно сказать, что в данный момент времени она имеет определенное состояние вследствие того, что пространство рассматривается как одно целое. В этой связи понятие одновременности для всех точек пространства будет связано с понятием состояния пространства в данный момент времени, которое является вполне определенным. Силы взаимодействия между материальными телами заменяются изменением состояния пространства. В этом случае отсутствует необходимость в релятивистском подходе при исследовании движения материальных тел, 128 Исследование движения планет, Луны и Солнца, основанное на новом принципе взаимодействия если считать, что изменения состояния всех точек пространства происходит одновременно. На основании проведенных исследований можно сделать следующие выводы: дифференциальные уравнения движения (11) вполне удовлетворительно описывают движение больших планет, Луны и Солнца на интервале времени 600 лет; они значительно проще дифференциальных уравнений, учитывающих релятивистские эффекты; кроме того, по затратам машинного времени более чем в 2 раза эффективнее последних.

About the authors

Anatoliy F Zausaev

Samara State Technical University

Email: zausaevaf@mail.ru
(Dr. Phys. & Math. Sci.; zausaevaf@mail.ru), Professor, Dept. of Applied Mathematics & Computer Science 244, Molodogvardeyskaya st., Samara, 443100, Russian Federation

References

Supplementary files