Студопедия Главная Случайная страница Обратная связь

Разделы: Автомобили Астрономия Биология География Дом и сад Другие языки Другое Информатика История Культура Литература Логика Математика Медицина Металлургия Механика Образование Охрана труда Педагогика Политика Право Психология Религия Риторика Социология Спорт Строительство Технология Туризм Физика Философия Финансы Химия Черчение Экология Экономика Электроника

Физика. Квантовая физика




ЗАДАЧА ОБ ОТНОСИТЕЛЬНОМ ДВИЖЕНИИ N ГРАВИТИРУЮЩИХ ТЕЛ

 

Кафедра астрономии

Базей Александр Анатольевич

Одесса 2001


ЗАДАЧА N ТЕЛ.

Задача N тел заключается в изучении движения конечного числа материальных точек под действием их взаимного притяжения по закону Ньютона.

Рассмотрим движение материальных точек вблизи центрального тела, масса которого много больше массы каждой точки. Пример такой модели – Солнечная система. Специфика Солнечной системы, в частности, в том, что {1}

v в самом Солнце сосредоточено 99% массы всей системы;

v планеты достаточно удалены друг от друга.

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

Переход к системе отсчета, связанной с центром Солнца

В силу сказанного, целесообразно начало системы отсчета связать с центром Солнца. Пусть имеется N гравитирующих тел:

m0 x0, h0, z0;

m1 x1, h1, z1;

. . . . .

mN-1 xN-1, hN-1, zN-1.

Расстояние между телами - rij, .

В инерциальной системе отсчета уравнения движения k-го тела:

[1]

Штрих у верхнего индекса суммирования, как и ранее, означает, что в сумме опускаются слагаемые с одинаковыми индексами.

Каждое слагаемое равно притяжению к каждому отдельному телу на основании закона Ньютона. Выделим притяжение Солнца:

[2]

Теперь выпишем уравнения движения Солнца:

[3]

Перейдем в уравнениях движения от сил к ускорениям. Для этого разделим систему [2] на mk, а систему [3] на m0. Полученные уравнения вычтем: [2]-[3].

Таким образом, здесь записаны дифференциальные уравнения движения относительно центра Солнца. Тем самым совершен переход к системе отсчета, связанной с центром Солнца. Сделано это путем вычитания уравнений движения Солнца из уравнений движения k–го тела. {2} Формально уравнения перехода можно записать так

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

[4]

полное ускорение k-го ускорение k–й планеты ускорение Солнца

ускорение тела Солнцем от i–х планет от i–х планет

[4]

Если в правой части в системе [4] стоят нули, (массы всех тел, за исключением m0 и mk, равны нулю) то уравнения обращаются в уравнения задачи 2-х тел. В задаче 2-х тел имеет место кеплерово движение, обычно называемое невозмущенным, а всякое отклонение от такого движения называют возмущением. {3}

Системе уравнений [4] можно придать изящную форму введением некоторой функции. Надо сконструировать такую функцию, чтобы ее производные по соответствующим координатам давали правые части уравнений движения [4]. {4}

Легко проверить, что такая функция

, [5]

где ri – гелиоцентрические расстояния возмущающих планет с координатами (xi, yi, zi),

xk, yk, zk – координаты возмущаемой планеты,

rki – взаимные расстояния.

Окончательно получим

[6]

Rk, наличие которого справа вызывает возмущения, называется возмущающей или пертурбационной функцией. Производные Rk по координатам точки mk дают компоненты ускорения, испытываемого телом mk в его относительном движении со стороны всех других тел, за исключением тела mk, принятого за начало координат. {5}

Поскольку уравнения системы [6] являются основой теории движения планет, то их иногда называют планетной формой уравнений относительного движения. При изучении гелиоцентрического движения планет за центральное тело принимают Солнце. В этом случае возмущающая функция мала, так как в нее входят множителями планетные массы; следовательно, правые части уравнений [6] будут малы, и их влияние действительно целесообразно рассматривать как возмущения. Отсюда, в частности, трудности теории движения Луны.

Система [6] порядка 6N-6 может рассматриваться как результат исключения при помощи интегралов сохранения импульса (формулы [3] и [4] из главы “ЗАДАЧА О ДВИЖЕНИИ N ГРАВИТИРУЮЩИХ ТЕЛ” (введение)) из системы [1] порядка 6N шести величин . {6}

 

ВОПРОСЫ.

1. В чем заключается специфика Солнечной системы с точки зрения задачи N тел?

2. Как совершить переход от инерциальной системы отсчета к системе отсчета, связанной с центром Солнца?

3. Что называют возмущением в движении тела?

4. Как строится возмущающая функция?

5. Что дают производные возмущающей функции по координатам данного тела?

6. Каков порядок системы дифференциальных уравнений, описывающих движение в задаче N тел относительно центрального тела и почему?

Возмущенное движение

Перепишем систему уравнений движения k-й планеты массой mk без индекса k:

[7]

где ,

x, y, z – координаты планеты,

r – расстояние этой планеты от Солнца,

R - соответствующая возмущающая функция.

Для Солнечной системы, как правило, полагают, что основная плоскость это плоскость эклиптики, а ось ОХ направлена в точку весны.

Заметим, что система [7], описывающая возмущенное движение – система неоднородных дифференциальных уравнений, а система [8] описывающая невозмущенное движение – система однородных дифференциальных уравнений:

[8]

Эта система уже решена, проинтегрирована, и в результате получена траектория движения планеты – орбита, которая описывается кеплеровыми элементами орбиты - постоянными величинами:

[9]

Решение уравнений [7] уже не приведет к кеплерову движению по коническому сечению. Но отказываться от простого и хорошо изученного кеплерова движения нежелательно. В таком случае предпочитают считать, что планета движется по коническому сечению, но само это сечение непрерывно изменяется. Плоскость орбиты изменяется: она поворачивается и покачивается. Сама орбита как бы «дышит», вытягивается или сокращается, поворачивается в своей плоскости, оставаясь, однако, в любой момент коническим сечением.

Поэтому решение уравнений [7] можно искать в той же форме, считая элементы орбиты не постоянными, а надлежащим образом выбранными функциями времени. Полученные таким образом элементы называются мгновенными элементами орбиты. {1} Совокупность этих элементов дает мгновенную орбиту планеты для момента t. {2} Знание мгновенной орбиты позволяет вычислять координаты планеты для любого момента t по формулам кеплерова движения [9].

Для невозмущенного движения элементы орбиты постоянны. Поэтому

, , .

Для возмущенного движения элементы орбиты зависят от времени. Поэтому

,

,

.

Поскольку мы имеем только 3 условия [9], которым должны удовлетворять шесть функций (система [7] – 6-го порядка, [9] определяет 3 величины), то мгновенные элементы определяются не однозначно – через точку в пространстве можно провести множество орбит с самыми разными элементами. Их нужно подчинить еще трем дополнительным условиям.

Потребуем, чтобы не только координаты , но и их производные выражались через мгновенные элементы орбиты. {3} То есть требуем

, , , или

,

, [10]

,

тогда элементы полностью определяют движение по возмущенной орбите. Такие мгновенные элементы, однозначно определяемые формулами [9] и [10], называются оскулирующими элементами орбиты.

Оскулирующими элементами для момента t называются такие элементы орбиты, которые дают положение и скорость для этого момента по формулам невозмущенного движения. {4}

Пусть Р – положение планеты в момент t0. Если в этот момент возмущающие планеты прекращают свое действие на тело Р, то его орбита будет эллипсом PQ с постоянными элементами . Если известны координаты и вектор скорости, то имеем необходимое количество уравнений для определения элементов орбиты

[11]

Действительная же орбита рассматриваемой планеты, когда учтено действие возмущающих планет, отличается от эллипса, это траектория РР’. Но в точке Р в момент t0 совпадают координаты и вектор скорости планеты на возмущенной и невозмущенной орбитах.

Невозмущенная орбита, имеющая общую точку и вектор скорости в момент t0 с возмущенной орбитой называется оскулирующей в момент t0, а постоянные - оскулирующими элементами орбиты для момента t0. {5} («оскулировать» - значит соприкасаться)

Иными словами, оскулирующая орбита представляет собой такую орбиту, по которой бы начала двигаться в некоторый момент планета, если бы все возмущения в этот момент внезапно исчезли.

Таким образом, оскулирующая орбита в момент t0 это та кеплерова орбита, которая соответствует положению и скорости рассматриваемой планеты в этот момент. Моменту t1, следующему за t0, будет соответствовать новая оскулирующая орбита с элементами , связанными с новыми координатами и скоростью планеты в ее действительном движении. И так далее, - оскулирующая орбита может быть вычислена для любого момента времени, каково бы ни было движение планеты, если известны ее положение и скорость в этот момент времени.

Разности называются возмущениями в элементах орбиты за время t1-t0.

Возмущения в координатах определяются аналогично.

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

 

ВОПРОСЫ.

1. Дайте понятие мгновенных элементов орбиты.

2. Как определяется мгновенная орбита?

3. В чем отличие мгновенных и оскулирующих элементов орбиты?

4. Что называется оскулирующими элементами орбиты?

5. Что называется оскулирующей орбитой?

Отличается ли скорость планеты на оскулирующей и на мгновенной орбите в момент оскуляции?

Отличается ли скорость планеты на истинной и на мгновенной орбите в момент оскуляции?

 

Уравнения Лагранжа

Далее будем искать вариацию оскулирующей орбиты. Перепишем систему уравнений [11] несколько иначе:

[12]

[13]

где - сопряженные элементы орбиты:

, - среднее суточное движение,

.

Такое разделение элементов орбиты можно заметить, например, в задаче 2-х тел при вычислении координат:

.

В возмущенном движении - функции времени. Продифференцируем [13] еще раз по времени:

,

ускорение в ускорение в

возмущенном невозмущенном

движении движении

,

.

Сравнивая полученное уравнение с системой [7], приходим к выводу

[14]

Теперь продифференцируем систему [12] по времени:

Но , поскольку по определению, движение по оскулирующей орбите происходит с той же скоростью, что и по истинной траектории в момент оскуляции. {1} Поэтому

[15]

В системы [14] и [15] входят вариации элементов орбиты . В совокупности имеем 6 уравнений и 6 неизвестных вариаций i=1,2,3. Формально поставленная задача решена – найдена вариация оскулирующей орбиты. Но реальное ее решение значительно сложнее.

 

Для решения уравнений систем [14] и [15] относительно упростим задачу, введя новые коэффициенты.

Производные возмущающей функции по координатам заменим на производные по элементам орбиты . Будем проводить преобразования первых уравнений систем [14] и [15]. Действия с другими уравнениями аналогичны.

Домножим [14] ; [15] и вычтем из первого произведения второе.

,

,

.

Под знаками сумм записаны скобки Пуассона:

, , тогда

. [16]

Аналогично -

,

.

Сложим эти три уравнения

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

[17]

Выражения вида , определяемые равенством [17], называются скобками Лагранжа относительно . {2}

Вспомним, что возмущающая функция R зависит от координат x, y, z рассматриваемой планеты и xk, yk, zk возмущающих планет или, что то же самое, R зависит от их элементов орбиты . Следовательно

.

Итак,

. [18]

Здесь записано 3 уравнения (i=1,2,3). Аналогичные уравнения для элементов орбиты группы :

. [19]

Системы [18] и [19] – уравнения Лагранжа для вариации элементов орбиты . Далее будет показано, что скобки Лагранжа являются функциями только шести элементов орбиты и явно не зависят от времени. Таким образом, имеется 6 уравнений [18] и [19], причем каждое является линейным относительно неизвестных . Если известны значения скобок Лагранжа и производные возмущающей функции по соответствующим элементам орбиты, то легко можно найти вариации элементов орбиты.

Если - оскулирующие элементы орбиты известны на какую-либо эпоху t, то в эпоху t+Dt элементы орбиты:

,

- в этом смысл уравнений Лагранжа. {3}

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

 

ВОПРОСЫ.

1. Отличается ли скорость планеты на оскулирующей и на истинной орбите в момент оскуляции?

2. Что называется скобками Лагранжа?

3. В чем смысл уравнений Лагранжа?

Свойства скобок Лагранжа

В этом пункте под элементом орбиты будем подразумевать любой элемент орбиты. Тогда, по определению, скобки Лагранжа

.

Свойство I. - это очевидно из определения.

Свойство II. - это следует из определения скобок Пуассона.

Свойство III. - скобки Лагранжа явно от времени не зависят. {1} Докажем это свойство.

Элементы орбиты характеризуют движение по коническому сечению. Движение по произвольной траектории в общем случае нельзя описывать с помощью кеплеровых элементов орбиты. В данном случае описание произвольного движения производится с помощью медленно изменяющихся кеплеровых элементов орбиты (поскольку произвольная траектория незначительно отличается от конического сечения из-за малости возмущающих сил). Для доказательства этого свойства необходимо воспользоваться уравнениями невозмущенного движения [8]. Тогда координаты x, y, z из [8] однозначно выражаются через кеплеровы элементы орбиты.

Введем обозначение . Тогда эти уравнения примут вид

.

Продифференцируем скобку Пуассона по времени:

, где . Подставим эту производную

, и прибавим ноль.

Получим

.

А так как скобка Лагранжа это сумма трех скобок Пуассона, то складывая

,

что и требовалось доказать.

Явная независимость скобок Лагранжа от времени позволяет после их вычисления придать времени произвольное значение. Это свойство широко используется.

 

ВОПРОСЫ.

1. Назовите свойства скобок Лагранжа?

Вычисление скобок Лагранжа

Для вычисления скобок Лагранжа потребуется вычислить частные производные координат по элементам орбиты.

Пусть x, y, z – эклиптические координаты возмущаемой планеты. Положение плоскости оскулирующей орбиты и направление на перигелий относительно эклиптической системы отсчета определяется элементами . Введем орбитальные прямоугольные координаты планеты . Ось направлена на перигелий орбиты, ось . Направляющие косинусы этих осей относительно эклиптической системы отсчета x, y, z были получены в главе «ЗАДАЧА 2-х ТЕЛ» и определены там равенствами [39]. Там же они обозначены как px,py,pz,qx,qy,qz. Найдем теперь направляющие косинусы прямой, направленной на полюс орбиты С.

Из сферической геометрии известно

 
 


 

Величины дуг ,

,

,

тогда из сферического треугольника CNХ, где следует

или

.

Из сферического треугольника CNY, где ,

или

.

Дуга CZ соединяет полюсы больших кругов – эклиптики и орбиты, угол между которыми i

или

Итак, направляющие косинусы всех трех осей орбитальной системы отсчета относительно эклиптической

[20]

Производные направляющих косинусов по элементам орбиты легко вычисляются

 

  px py pz qx qy qz
W - py px - qy qx
w qx qy qz - px - py - pz
i rxsinw rysinw rzsinw rxcosw rycosw rzcosw

 

ВЫЧИСЛЕНИЕ СКОБОК ЛАГРАНЖА ТИПА

По определению

.

Здесь координаты зависят от элементов орбиты , но если перейти к орбитальной системе отсчета, то не зависят от (в орбитальной системе отсчета очевидно ).

Теперь

раскрываем скобки и приводим подобные слагаемые

 

 

 

Сомножитель при сумме – постоянная величина (закон сохранения момента импульса, «ЗАДАЧА 2-х ТЕЛ» - [37]):

, где .

( , поскольку слева стоят орбитальные координаты)

Итак

[21]

 

Для дальнейших вычислений найдем некоторые вспомогательные величины. Пусть g - любой из 6-ти элементов орбиты. Поскольку , то , и, подставляя производные из таблицы

аналогично

[22]

Пользуясь этими равенствами, вычислим

.

Легко проверить, что , и по свойству направляющих косинусов .

[23]

Продифференцируем равенство [23] по d, где d – тоже любой из 6-ти элементов орбиты:

[24]

Перепишем [23]

и снова продифференцируем по g:

[25]

Вычтем [25]-[24]:

[26]

 

Теперь вернемся к формуле [21], с учетом полученного:

[27]

В частности, если , то

,

.

Если , то

.

В правой части перейдем к другим элементам орбиты:

, где n - среднее движение («ЗАДАЧА 2-х ТЕЛ» - [59]: ),

Похожим образом вычисляются скобки Лагранжа типа и . Значения всех скобок Лагранжа приведены в следующей таблице (первым в скобках Лагранжа выбирается элемент, стоящий в заголовке строки, вторым – стоящий в заголовке столбца):

  a e c W w i
a
e
c
W
w
i

Здесь .

 

Таким образом, все шесть уравнений Лагранжа [18] и [19] в развернутом виде:

 

,

,

,

,

,

.

Подставляя сюда значения скобок Лагранжа:

,

,

, [28]

,

,

.

Систему уравнений Лагранжа [28] надо решить относительно :

,

,

, [29]

,

,

.

Отметим, что время входит в уравнения Лагранжа только посредством производных пертурбационных функций R. {1}

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

При изучении Солнечной системы часто встречаются движения по слабонаклоненным орбитам. Чем ближе к нулю наклонение i, тем менее определенным становится положение восходящего узла орбиты. При i=0 положение узла становится полностью неопределенным. Вследствие этого долгота восходящего узла W находится из наблюдений с ошибкой, неограниченно возрастающей по мере приближения i к нулю.

Чтобы избежать возникающих отсюда неудобств при определении положения перигелия, вместо аргумента перигелия w пользуются элементом орбиты , называемым долготой перигелия. {2} Кроме того, Вместо элемента вводят среднюю долготу в эпоху . Тогда уравнения Лагранжа выглядят немного иначе.

 

ВОПРОСЫ.

1. Входит ли в уравнения Лагранжа время?

2. Что называется долготой перигелия и с какой целью она вводится?

Возмущенное движение в Солнечной системе

Уравнения Лагранжа были широко использованы для аналитического изучения движений в Солнечной системе.

Решение системы [29] можно искать методом последовательных приближений из-за малости планетных масс по сравнению с солнечной. Тогда решение будет иметь вид

,

, [30]

Ряды расположены по целым степеням масс, через a0, e0,... обозначены оскулирующие элементы орбиты в некоторый начальный момент времени t=0. Через обозначены члены рядов n-й степени относительно масс. Эти величины называются возмущениями n-го порядка. {1}

Для оценки важности отдельных возмущений употребляется несколько характеристик.

Прежде всего – это порядок возмущения. Чем ниже порядок, тем больше (при прочих равных условиях) влияние рассматриваемого возмущения. {2}

С другой стороны все слагаемые в решениях типа [30] делятся на периодические, вековые и смешанные возмущения.

Кроме того, возмущения характеризуются своей степенью, рангом и классом. Эти характеристики мы здесь не рассматриваем.

При изучении движения планет приходится иметь дело главным образом с периодическими возмущениями 1-го порядка. В зависимости от величины периода, возмущения делятся на:

v короткопериодические – их периоды меньше или равны периодам возмущающей или возмущаемой планет, амплитуды таких членов и соответствующие возмущения обычно невелики; {3}

v долгопериодические – их наличие часто обусловлено приблизительной соизмеримостью средних угловых движений n возмущающей и возмущаемой планет, величины их обычно на 2 порядка больше величин короткопериодических членов, и поэтому им соответствуют очень большие возмущения. {4}

С долгопериодическими возмущениями астрономы впервые встретились при попытке объяснить движения Юпитера и Сатурна их взаимным притяжением. Под влиянием долгопериодических неравенств долгота Сатурна может изменяться на 50’.

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

Наличие вековых членов первого порядка обусловлено той частью возмущающей функции, которая при разложении в ряд не зависит от времени. Если при решении уравнений Лагранжа возникают вековые члены в возмущении данного элемента орбиты, то он будет монотонно изменяться с малой постоянной скоростью (вековым образом) без ограничений. {5}

Теорема Лапласа-Лагранжа. Если средние движения n планет несоизмеримы, то большие полуоси а не имеют вековых возмущений 1-го порядка. {6}

Эта теорема связана с вопросом об устойчивости Солнечной системы. При рассмотрении возмущений 2-го и более высоких порядков в изменении большой полуоси а появляются вековые члены, и следовательно, область изменения а уже не будет ограниченной для всех значений времени – устойчивость, в прежнем смысле, будет отсутствовать. На первый взгляд кажется, что об устойчивости Солнечной системы говорить нельзя. Но здесь надо вспомнить, что при решении уравнений Лагранжа возмущающая функция R раскладывается в ряд, а значит промежуток времени, на котором это разложение справедливо – ограничен. Поэтому ясно, что возмущения высоких порядков получены с определенной точностью для ограниченных интервалов времени. Таким образом, об устойчивости Солнечной системы этот анализ ничего не говорит. {7}

С другой стороны, есть основания полагать, что вековые возмущения это периодические возмущения с очень большими периодами -продолжительностью несколько тысяч лет. Вековой член lt может оказаться приближением к долгопериодическому члену sin lt, когда lt мало. Это справедливо для любого элемента орбиты.

 

ВОПРОСЫ.

1. Что называется возмущениями n–го порядка?

2. Что больше при прочих равных условиях: влияние возмущения 1-го или 2-го порядка?

3. Дайте понятие короткопериодических возмущений.

4. Дайте понятие долгопериодических возмущений, когда они обычно возникают?

5. Дайте понятие вековых возмущений.

6. Сформулируйте теорему Лапласа-Лагранжа.

7. Известно, что при рассмотрении возмущений 2-го и высших порядков появляются вековые возмущения больших полуосей орбит планет. Что можно сказать об устойчивости Солнечной системы?

 

ВОЗМУЩЕНИЯ. ГЕОМЕТРИЧЕСКИЕ СООБРАЖЕНИЯ.

Возмущающая сила

Напомним, что кеплерово движение в задаче 2-х тел, обычно называют невозмущенным, а всякое отклонение от такого движения называют возмущением.

Предположим, что тело m движется из точки P0 в направлении Q0 со скоростью V0 и что кроме притяжения тела S не имеется сил, действующих на него. Тогда согласно результатам задачи 2-х тел оно будет двигаться по коническому сечению С0, элементы которого однозначно определены начальными условиями.

Предположим, что когда m приходит в точку P1, то оно подвергается мгновенному импульсу в направлении P1Q1; это положение и новый вектор скорости определяют новое коническое сечение С1, по которому будет двигаться тело. Придя в точку Р2, тело m подвергается действию мгновенного импульса в направлении P2Q2, и это определит новую коническую траекторию С2. Итак тело будет двигаться по коническим сечениям, время от времени изменяющимся под влиянием возмущающих импульсов.

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

Для нахождения действия возмущающей силы на элементы орбиты ее удобно разложить на три прямоугольных составляющих. Это можно сделать несколькими способами. Мы разложим возмущающую силу (рис.5) на три составляющих так:

S – ортогональную составляющую, перпендикулярную к плоскости орбиты и положительную, если она направлена в сторону северного полюса эклиптики;

T – тангенциальную составляющую, направленную по касательной к траектории и положительную, если она действует по направлению движения; (лежит в плоскости орбиты)

N – нормальную составляющую, перпендикулярную к касательной к траектории и положительную, если она направлена внутрь орбиты. (лежит в плоскости орбиты)

 

 

ВОПРОСЫ.

1. Как направлена ортогональная составляющая возмущающей силы?

2. Как направлена тангенциальная составляющая возмущающей силы?

3. Как направлена нормальная составляющая возмущающей силы?

 

Действия составляющих возмущающей силы

Для краткости и определенности изложения рассмотрим возмущения на примере системы Земля-Луна-Солнце. Далее:

возмущаемое тело – Луна;

возмущающее тело – Солнце.

Это случай, который рассматривал Ньютон по существу такими же методами, как и применяемые здесь ниже. Рассмотрим характер возмущений, возникающих от одних положительных составляющих; отрицательные составляющие изменяют элементы в обратную сторону.

В рассмотрении находятся независимые элементы орбиты a,e,T0,w,W,i.

 

ВОЗДЕЙСТВИЕ ОРТОГОНАЛЬНОЙ СОСТАВЛЯЮЩЕЙ S

Поскольку S перпендикулярна плоскости орбиты, то она может изменять только положение плоскости орбиты в пространстве, но не величину и форму самой орбиты. {1} S не изменяет a,e,T0. Если w отсчитывается от неподвижной линии в плоскости орбиты, то и аргумент перигея остается неизменным, но обычно w считается от восходящего узла орбиты. Поэтому, если отрицательную скорость возрастания долготы узла W умножить на cosi, то результат даст скорость возрастания w вследствие изменения начала, от которого он отсчитывается (рис.6).

Следовательно, при рассмотрении возмущений, вызываемых ортогональной составляющей S, достаточно рассмотреть изменения в W и i.

Пусть P0Q0 (рис.7) лежит в плоскости невозмущенной орбиты. Предположим, что когда Луна находится в P0, происходит мгновенный импульс P0S0 (S>0); тогда она будет двигаться по направлению P0P1 и новые узел и наклонение будут W1 и i1. Очевидно, что i1 > i0 и W1 < W0. Предположим, что когда Луна приходит в P1, имеет место новый мгновенный импульс P1S1 (S<0). Обозначим новые узел и наклонение W2 и i2, и очевидно, что i2 < i1 и W2 < W1. Если P0W1 = W1P1, P0S0 = P1S1 и скорость Луны в P0 равна скорости в P1, то i0 = i2. Общий результат – обратное движение узла и неизменное наклонение – восходящий узел все время отступает.

Из рис.8 для нисходящего узла орбиты видно, что положительная ортогональная составляющая S перед пересечением эклиптики увеличивает узел и уменьшает наклонение орбиты; а после пересечения эклиптики S<0 снова приводит к увеличению узла и увеличению наклонения орбиты.

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

Итак, положительная S заставляет узлы двигаться вперед, если Луна находится в I и II четверти, и отступать, если Луна находится в III и IV четверти. Аналогично можно сделать выводы и для наклонения орбиты:

  a e w W i
S>0 0, если отсчитывать от неподвижной линии в плоскости орбиты I, II четв. - увеличение; III, IV четв. - уменьшение I, IV четв. - увеличение; II, III четв. - уменьшение

 

ВОЗДЕЙСТВИЕ ТАНГЕНЦИАЛЬНОЙ СОСТАВЛЯЮЩЕЙ T

Тангенциальная составляющая возмущающей силы не изменяет долготы восходящего узла W и наклонения i, поскольку лежит в плоскости орбиты. {2}

Большая полуось определяется из интеграла энергии . Для эллиптической орбиты a>0, а так как положительная Т увеличивает скорость V и не изменяет мгновенно расстояние r, то положительная Т увеличивает большую полуось а в любой части орбиты Луны. {3}

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

Из интеграла энергии , , ;

заметим, что ,

.

Подставляя, - чем больше скорость, тем больше изменение а под действием тангенциальной составляющей возмущающей силы, то есть - в перигее орбиты Луны. {4}

 

Тангенциальная составляющая изменяет величину скорости, но не изменяет направления движения. Фокус F2 остается неизменным, поскольку там находится Земля, r2 остается постоянным (рис.9). Но, согласно результатам предыдущего рассмотрения, большая полуось возрастает. Так как r2=2a-r1 в то время как r2 неизменно, то увеличивается r1, фокус F1 перемещается в положение F1’. Линия апсид вращается вперед от AП к A’П’. Поэтому – положительная тангенциальная составляющая возмущающей силы заставляет линию апсид вращаться вперед в течение первой половины периода, и назад в течение второй половины периода; в перигее и апогее линия апсид неподвижна.

Отметим, что поскольку большая полуось Луны изменяется сильнее всего в перигее, то наиболее быстро линия апсид вращается вблизи KL.

 

Эксцентриситет орбиты определяется, как (рис.10). Когда Луна находится в перигее при T>0 большая полуось а - растет, r2 – неизменно:

2a=r1+r2. Следовательно, r1 - растет, причем и новый эксцентриситет:

.

Вблизи перигея положительная тангенциальная составляющая возмущающей силы увеличивает эксцентриситет орбиты.

Пусть теперь Луна находится в апогее, (рис.11) при T>0 большая полуось а по-прежнему растет, r2 – неизменно, поэтому r1 – растет: . Поэтому расстояние F1F2 уменьшается на , большая ось возрастает на эту величину

.

Вблизи апогея положительная тангенциальная составляющая возмущающей силы уменьшает эксцентриситет орбиты.

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

  a e w W i
T>0 всегда увеличение MПN - увеличение; NAM - уменьшение I, II четв. - увеличение; III, IV четв. - уменьшение

 

 

ВОЗДЕЙСТВИЕ НОРМАЛЬНОЙ СОСТАВЛЯЮЩЕЙ N

Нормальная составляющая возмущающей силы не изменяет долготы восходящего узла W и наклонения i, поскольку лежит в плоскости орбиты. {5}

Из интеграла энергии следует, что большая полуось зависит только от величины скорости, а не от ее направления. Так как N действует перпендикулярно к касательной, то она изменяет только направление скорости и, значит, не влияет на величину большой полуоси а. {6}

 

PT – касательная к орбите (рис.12). Действием нормальной составляющей она изменяется в PT’. По свойству эллипса радиусы к двум фокусам образуют равные углы с касательной, поэтому r1 перейдет r1’; и так как N не влияет на длину большой полуоси, то длина r1=r1’.

Следовательно, если Луна находится в области LAK, то положительная нормальная составляющая возмущающей силы будет вращать линию апсид вперед (w возрастает), а если в области KПL, то назад (w уменьшается). В точках K и L нормальная составляющая не изменяет направления линии апсид.

 

По определению . После действия нормальной составляющей эксцентриситет становится равным (рис.12), так как большая полуось не изменяется. Сравнивая треугольник PF1F2 и треугольник PF1’F2, видим, что сторона PF2 общая, PF1=PF1’. Отсюда следует, меньше , и значит : e’<e – эксцентриситет уменьшается.

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

  a e w W i
N>0 I, II четв. - уменьшение; III, IV четв. - увеличение в интервале LAK - увеличение; в интервале KПL - уменьшение

 

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

 

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

 

ВОПРОСЫ.

1. Какие элементы орбиты изменяет ортогональная составляющая возмущающей силы и почему?

2. Какие элементы орбиты изменяет тангенциальная составляющая возмущающей силы и почему?

3. Как тангенциальная составляющая изменяет большую полуось орбиты?

4. В какой части орбиты влияние тангенциальной составляющей возмущающей силы наиболее эффективно?

5. Какие элементы орбиты изменяет нормальная составляющая возмущающей силы и почему?

6. Изменяет ли нормальная составляющая величину большой полуоси и почему?

Негравитационные возмущения

ВОЗМУЩАЮЩЕЕ ДЕЙСТВИЕ СОПРОТИВЛЯЮЩЕЙСЯ СРЕДЫ

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

наклонение i и долгота восходящего узла W не изменяются;

большая полуось а непрерывно уменьшается;

линия апсид (w) подвергается периодическим изменениям, вращаясь назад в течение первой половины обращения и настолько же вперед в течение второй половины;

эксцентриситет периодически уменьшается и увеличивается, с другой стороны сопротивление среды при больших скоростях зависит от V2, поэтому изменение значительно больше в перигее, чем в апогее, и все действие в течение периода заключается в уменьшении эксцентриситета. {1}

 

ВОЗМУЩЕНИЯ, ВОЗНИКАЮЩИЕ ИЗ-ЗА СЖАТИЯ ЦЕНТРАЛЬНОГО ТЕЛА

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

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

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

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

Описанные возмущения тем больше, чем больше сжатие планеты и чем ближе спутник. В случае Луны возмущения от сжатия Земли очень малы. Наиболее яркий пример возмущений этого рода в Солнечной системе мы находим в орбите 5-го спутника Юпитера. Сжатие Юпитера настолько велико и орбита спутника так мала, что его линия апсид продвигается приблизительно на 900О в год.

 

Кроме того, возможно воздействие реактивной составляющей в случае, например, испарения газа с поверхности тела (кометы). Направление и величина возмущающей силы в этом случае зависит от многих факторов.

Световое и корпускулярное излучение Солнца направлено по радиус-вектору в сторону, противоположную направлению силы тяжести. Поэтому можно говорить об уменьшении притяжения из-за светового и корпускулярного давления Солнца. Для учета этого факта эквивалентно уменьшают массу Солнца. В случае рассмотрения движения пылевых частиц и очень мелких тел эту силу надо учитывать обязательно. {2}

 

 

ЭФФЕКТ ПОЙНТИНГА-РОБЕРТСОНА







Дата добавления: 2014-10-29; просмотров: 194. Нарушение авторских прав

Studopedia.info - Студопедия - 2014-2017 год . (0.05 сек.) русская версия | украинская версия