2017 есть дефекты Небесная механика Банникова Е

2017 ( есть дефекты) Небесная механика Банникова Е. Ю.

Законы Кеплера • Первый закон Кеплера… и длина эллипса Параметрическое уравнение эллипса x acos y bsin Длина эллипса 2 L rd E(e) 1 e /2 sin 2 d 2 0 0 2 y d 4 a. E(e) x 2 - полный эллиптический интеграл 2 -го рода 2 Длина окружности E(0) dl / L 2 a

Гравитационный потенциал Материальной точки • • GM r Сферической оболочки GM , r a r Sph GM a , r a Внутри сферической оболочки пробная частица находится в невесомости. Теорема Ньютона: обобщение на эллипсоидальный слой

Притяжение пробной частицы внутри сферы: элементарные соображения Телесный угол d d. S 2 r 1 d. S 2 r 2 2 1 2 21 F dm 1 r 2 d. S r 2 F 2 r 1 dm 2 r 1 d. S 2 1 1 2 2 M S F 2 r r 2 Const

Гравитационный потенциал шара Потенциал сферы a GM 4 G 2 r r Sph r a GM 4 G a, r a a шар сфера внутр шар G 4 r 0 dr GM r (3 23 a внутр шар GM r внеш внутр сфера r 2 шар , внешн сфера a r a 2 a r dr r 2) r

Гравитационный потенциал шара Теорема Дирихле 4 G , внутри объема вне объема 1 out r r a 2 2 r r r 0 r , out 0, c 2 GM r a out c 1 r out c 2 GM r 1 2 out 2 inner r 2 4 r r r 3 r 2 G G c 1 r c 2

Гравитационный потенциал шара Учитываем сшивку на границе потенциала и силы Потенциал шара out r a inner r a GM , r a r шар GM (3 a 2 r 2), r a 2 a 3

Гравитационный потенциал • Интегрирование по объему • Суммирование по элементарным составляющим • Теорема Дирихле Сфера, цилиндр, шар….

Гравитационный потенциал эллипсоида Теорема Лапласа Однородные софокусные эллипсоиды притягивают внешнюю точку с силами, одинаково направленными, а по величине пропорциональными их массам Fx M Fx/ M / Теорема Ляпунова Шар обладает минимальной потенциальной энергий

Задача Эйлера о двух неподвижных массах • Гравитационный потенциал сжатого сфероида эквивалентен потенциалу стержня мнимой длины. • Метод эквигравитирующих стержней

Разложение потенциала в ряд Лапласа

Задача многих тел • Произвольная инерциальная с. к. Rj R i mj R 3 G Ri N j 1 ij Порядок системы 6(n+1) (0) Ri R (t 0 ) i R Ri (t 0 ) i (0) Первые интегралы: i i mi R a, mi R a t b i R m R i i Etot T U i закон движения центра масс i I закон сохранения момента количества движения закон сохранения энергии В скалярном виде 10 первых интегралов в произвольной инерциальной с. к.

Задача двух тел: Произвольная инерциальная с. к. R R Gm R 1 2 2 R 1 3 12 R R R 2 Gm 2 1 1 R 12 3 R(0) 1, 2 0 ) R (t 1, 2 (0) R 1, 2 R (t ) 1, 2 m Барицентрическая с. к. i 3 m 1 2 1 G ( 3 m 1 m 2 2) 1 3 G (1 m 2 2 m 1 m 2 2) 3 2 i i 0 0

Задача двух тел Относительная система координат r r G(M m) r 3 r (0) r(t ) 0 (0) r r(t ) 0 0 Первые интегралы: r r I момент на единицу массы Etot T U Порядок системы =6, но 4 первых (в скалярах) интегралов. Не хватает…. .

Задача двух тел. Интеграл Лапласа r I 0 r 3 G(M m) Интеграл Лапласа r I r r вектор Лапласа Уравнения связи между первыми интегралами I 0 Вектор момента и вектор Лапласа перпендикулярны 2 2 E I 2 tot 5 независимых первых интегралов => задача двух тел в относит с. к. (система 6 -го порядка) сводится к одному уравнению

Задача двух тел. Орбитальная с. к. I h r M r r r (r I ) r r r cos m n x p r 1 ecos rsin r I 2 r p I 2 e

Задача двух тел

Уравнение Кеплера r 2 I n n tg n (t dn ) 3/2 n 2 (1 ecos ) p 0 1 E e tg 2 1 2 e Уравнение Кеплера E esin. E n(t ) T 2 a 3/2 G(M m)

Смещение перигелия Меркурия Ньютоновское приближение rg GM 2 c 2 2 GMm Em du 2 u 2 u I dn I 2 Релятивистская задача 2 rgmc 2 mc du 3 2 u r 2 u 22 (1 d I I g n E 2 ) mc 2 4 Максимальное смещение перигелия наблюдается для Меркурия и составляе 43’’ за 100 лет.

Задача трех тел r r 1 2 r 2 1 3 r 1 r 23 U x 2 ny x y 2 nx U y U z z 1, 2 Gm 1, 2 n 2 U(x, y, z) 2 (x 2 y 2 ) 1 2 r r 2 1 V 2 2 U CJ C J - интеграл Якоби (интеграл относительной энергии)

Задача трех тел n 2 U(x, y, z) 2 (x V 2 U C 2 J Поверхности нулевой скорости 2 1 n (x y ) 2 r 1 2 r 2 CJ - интеграл Якоби 2 2 2 CJ 2 2 ) y 1 2 r 1 r 2

Кривые нулевой скорости m=0. 04; x 0=1. 179; y 0=0; z 0=0. 0; Vx 0=0. 0; Vy 0=-0. 238; Vz 0=0. 0; m=0. 04; x 0=1. 12; y 0=0; z 0=0. 0; Vx 0=0. 0; Vy 0=-0. 238; Vz 0=0. 0;

Точки Лагранжа L 4, L 5 r 1 2 r 1 x 1/ 2 2 y 2 2 1 L 1, L 2 L 3 1/3 2 r 1 2 3 2 7 r 1 12 1 1 3/

Семейство Хильды Резонанс 3: 2 L 3, L 4, L 5 – афелии астероидов

Янус и Эпиметей

Метод Лагранжа оскулирующих элементов r r 3 r F возм F 0 На малом интервале – невозмущенное кеплеровское движения, соответствующее разным начальным условиям. Возмущенная орбита является огибающей семейства невозмущенных Планетные уравнения Лагранжа d. Ei j dt (E , U i возм ) Uвозм - возмущающий потенциал i E (i, , , e, p, )

Метод Лагранжа оскулирующих элементов Разложение, усреднение U…. Вековые (секулярные) возмущения a 0 a e 0 e A(cos 0 0 Bt cos ) 0 Ct Астероид, возмущенный Юпитером

Метод Лагранжа или метод оскулирующих элементов d. Ei i (t, E j) dt Для двух планет E E 0 A 0(t t 0 ) Ak 1, k 2 k 1 n k 2 n 2 1 1 k 1, k 2 1 sin(k M 2 Bk 1, k 2 2 k. M ) k 1 n k 2 n 2 1 cos(k M 1 1 2 2 k. M ) k 1, 2 – собственная частота движения планеты и частота возмущающей силы (m Периодические возмущения : Amlp T A(k 1 , k 2 ) k 1 n k 2 n 2 1 2 k 1 n k 2 n 2 1

Метод Лагранжа или метод оскулирующих элементов • Резонансные возмущения Пример – астероид на резонансной орбите с Юпитером a. J 5. 20 a. A 3. 27 TA TJ 0. 5 2 n J n A

Метод Лагранжа или метод оскулирующих элементов • Возмущения кеплеровских элементов Юпитер-Сатурн Пример Юпитер-Сатурн TJ 2 TS 5 Эксцентриситет орбиты каждой из планет периодически изменяется с периодом 70 100 лет Наклонение – 51 000 лет

Метод Лагранжа или метод оскулирующих элементов • Короткопериодические возмущения k 1 n k 2 n 2 k n 1 1 1 T Amlp A(k 1 , k 2 ) k 1 n k 2 n 2 1 2 k 1 n 1 Периоды возмущений порядка орбитального периода, амплитуда мала. • Долгопериодические возмущения (при малых k 1, k 2) k 1 n k 2 n 2 0 1 n 1 k 2 n k 1 2 Отношение средних движений =простой дроби - резонансное состояние

Спутник-пастух колец Сатурна Открытие в 1990 г. при анализе изображений Вояджер-2 (1981 г. ) 1991 г. – официально назван в честь бога пастухов Расположен внутри люка Энке и движется почти в плоскости экватора Сатурна

Спутник-пастухи Сатурна Прометей (внутренний) Пандора (внешний) Cassini image

Спутники-пастухи Cassini's confirmation that a small moon orbits within the Keeler gap in Saturn's rings is made all the more exciting by this image, in which the disk of the 7 kilometer-wide body is resolved for the first time. The new body, provisionally named S/2005 S 1, was first seen in a time-lapse sequence of images taken on May 1, 2005, as Cassini began its climb to higher elevations in orbit around. In the vicinity of the little moon, the Keeler gap edges bear striking similarities to the scalloped edges of the 322 kilometer-wide Encke gap, where the small moon Pan (25 kilometers across) resides. From the size of the waves seen in the scalloped edges of the Encke gap, imaging scientists were able to estimate the mass of Pan. They expect to do the same eventually with S/2005 S 1. This image was obtained with the Cassini spacecraft narrow-angle camera on May 2, 2005, at a distance of about 594, 000 kilometers from Saturn. https: //www. nasa. gov/mission_pages/cassini/multimedia/pia 06237. html

Гиперион (спутник перевертыш) Двуликий… Япет Спин-орбитальный резонанс Cassini images

Астероид Круитни Орбитальный резонанс с Землей 1: 1 Проекция на эклиптику Сопутствующая ситема к.

Сечение Пуанкаре Сечение (отображение) Пуанкаре – отображение проекции траектории на выделенную плоскость фазового пространства. Для траектории на плоскости: 4 -х мерное фазовое пространство ( x, y, x, y) Одну из переменных можно исключить, воспользовавшись интегралом Якоби или полной энергией => 3 D фазовое пр-во ( x, y, x ) Выделяем одну из плоскостей, например, y=0. Т. о. получаем проекцию на фазовую плоскость ( x, x )

Сечение Пуанкаре 2 Потенциал Хенона-Хейлиса U(x, y) 1 x y 2 2 2 2 x 2 E =1/40 tot N=50 – число частиц 150 периодов 3 y y 3 Траектория одной из частиц 150 периодов

Сечение Пуанкаре Потенциал Хенона-Хейлиса E =1/40 tot 1 2 2 y y 3 2 x U(x, y) x 2 3 N=50 – число частиц a few 100 periods

Сечение Пуанкаре Потенциал Хенона-Хейлиса E =1/20 tot 1 2 2 U(x, y) x y 2 2 y y 3 2 x 2 3

Сечение Пуанкаре Потенциал Хенона-Хейлиса E =1/12 tot 1 2 2 U(x, y) x y 2 2 y y 3 2 x 2 3

Сечение Пуанкаре Потенциал Хенона-Хейлиса E =1/8 tot 1 2 2 U(x, y) x y 2 2 y y 3 2 x 2 3

Сечение Пуанкаре Потенциал Хенона-Хейлиса E =1/6 tot 1 2 2 U(x, y) x y 2 2 y y 3 2 x 2 3