В связи с разработкой космогонической гипотезы О. Ю. Шмидта в 40-х гг. в СССР были выполнены многочисленные исследования финальных движений в задаче трёх тел; полученные в этих работах результаты имеют значение на неограниченном интервале времени. В США (1965) численным методом изучена эволюция орбит пяти внешних планет на интервале времени в 120 000 лет. Самым интересным результатом этой работы явилось открытие либрации Плутона относительно Нептуна, благодаря которой минимальное расстояние между этими планетами не может быть меньше 18 астрономических единиц, хотя в проекции на плоскость эклиптики орбиты Плутона и Нептуна пересекаются. В СССР выполнена обширная работа (1967) по применению теории вековых возмущений Лагранжа — Брауэра к изучению эволюции орбиты Земли на протяжении миллионов лет. Эта работа имеет важное значение для понимания изменения климата Земли в различные геологические эпохи.
Начало 20 в. было отмечено значительным прогрессом в разработке математических методов Н. м. Этот прогресс был связан прежде всего с работами французского математика А. Пуанкаре , русского математика А. М. Ляпунова и финского астронома К. Сундмана. Последнему удалось решить общую задачу трёх тел с помощью бесконечных степенных сходящихся рядов. Однако ряды Сундмана оказались совершенно непригодными для практического использования из-за их крайне медленной сходимости. Сходимость рядов в Н. м. тесно связана с так называемой проблемой малых делителей. Математические трудности этой проблемы в значительной степени преодолены в работах математиков школы А. Н. Колмогорова .
Развитие Н. м. в СССР тесно связано с деятельностью двух научных центров, возникших непосредственно после Великой Октябрьской социалистической революции: Теоретической астрономии института АН СССР в Ленинграде и кафедры небесной механики Московского университета (см. Астрономический институт : имени П. К. Штернберга ). В этих двух центрах сложились ленинградская и московская школы, которые определили развитие Н. м. в СССР. В Ленинграде вопросы Н. м. разрабатывались главным образом в связи с такими практическими задачами, как составление астрономических ежегодников, вычисление эфемерид малых планет и др. В Москве доминирующее влияние на протяжении многих лет имели космогонические проблемы, а также астродинамика.
Среди иностранных научных учреждений, ведущих исследования в области Н. м., видное место занимают: Вашингтонская морская обсерватория , Гринвичская астрономическая обсерватория , Бюро долгот в Париже, Астрономический институт в Гейдельберге и др.
Релятивистская Н.м. В середине 20 в. в связи с повышением точности оптических наблюдений небесных тел, развитием новых методов наблюдений (наблюдения доплеровского смещения, радиолокация и лазерная локация) и возможностью проведения экспериментов в Н. м. при помощи космических зондов и искусственных спутников всё большее значение приобретает учёт релятивистских эффектов в движении тел Солнечной системы. Эти проблемы решаются релятивистской Н. м., опирающейся на общую теорию относительности Эйнштейна (см. Тяготение ). Роль общей теории относительности для Н. м. не ограничивается учётом малых поправок к теориям движения небесных тел. С появлением общей теории относительности удалось дать объяснение явлению тяготения, и таким образом Н. м. как наука о гравитационном движении небесных тел по существу становится релятивистской.
Согласно основной идее общей теории относительности, свойства пространства событий реального мира определяются движением и распределением масс, а движение и распределение масс, в свою очередь, определяются метрикой пространства-времени . Эта взаимосвязь находит своё отражение в уравнениях поля — нелинейных уравнениях с частными производными, определяющих метрику поля. В теории тяготения Ньютона уравнения движения (законы механики Ньютона) постулируются отдельно от уравнений поля (линейные уравнения Лапласа и Пуассона для ньютонова потенциала). В общей же теории относительности уравнения движения тел содержатся в уравнениях поля. Однако строгое решение уравнений поля, представляющее интерес для Н. м., и вид строгих уравнений движения задачи n тел, даже для n = 2, в общей теории относительности не получены. Лишь для n = 1 удалось найти строгие решения уравнений поля: решение Шварцшильда для сферически симметричного неподвижного тела и решение Керра, описывающее поле вращающегося тела сферической структуры. Для решения задачи n тел (n > 2) приходится прибегать к приближённым методам и искать решение в виде рядов по степеням малых параметров. Таким параметром в случае движения тел Солнечной системы часто служит отношение квадрата характеристической скорости орбитального движения тел к квадрату скорости света. Вследствие малости этого отношения (около 10-8 ) в уравнениях движения и их решениях достаточно для всех практических приложений учитывать лишь члены первой степени относительно этого параметра.
Релятивистские эффекты в движении больших планет Солнечной системы могут быть получены с достаточной точностью на основе решения Шварцшильда. Основным эффектом при этом является вековое смещение перигелиев планет. В решении Шварцшильда имеется также релятивистский вековой член в движении узла орбиты, но выделить этот эффект в явном виде из наблюдений не удаётся. Частично этот вековой член учитывается в радиолокационном эффекте при радиолокации Меркурия и Венеры с Земли (радиолокационный эффект состоит в дополнительном по сравнению с ньютоновским запаздыванием сигнала при возвращении его на Землю). Этот эффект подтвержден экспериментально. Релятивистские эффекты в движении малых планет и комет выявить достаточно уверенно пока не удаётся из-за отсутствия хорошо разработанной ньютоновской теории движения этих объектов и недостаточного количества точных наблюдений.