Au fondement théorique de la mécanique céleste se trouve l'application des lois de Newton couplée à la conservation de l'énergie mécanique, du moment cinétique et, dans le cas newtonien pur, du vecteur de Laplace-Runge-Lenz. Lorsque l'on considère un système isolé à deux corps, les équations différentielles du mouvement sont intégrables et admettent des solutions analytiques exactes correspondant aux coniques : ellipses pour les orbites liées, paraboles ou hyperboles pour les trajectoires de passage. Cette intégrabilité exceptionnelle disparaît dès que l'on introduit un troisième corps interactif, marquant la transition vers le problème à trois corps, qui s'est révélé non intégrable en général. Cette limitation analytique a ouvert la voie à l'étude des systèmes dynamiques non linéaires, à la théorie des bifurcations et à la naissance de la mécanique chaotique, dont Poincaré a posé les bases en démontrant la sensibilité extrême aux conditions initiales et l'existence de structures complexes dans l'espace des phases.

Pour traiter les configurations réelles, où de multiples perturbations gravitationnelles et non gravitationnelles s'accumulent, la mécanique céleste a développé un arsenal de méthodes d'approximation. La théorie des perturbations, initiée par Laplace, Lagrange et poursuivie par Delaunay, Brouwer ou encore Lie, permet de calculer les déviations orbitales engendrées par des corps secondaires, l'aplatissement planétaire, les effets de marée, la pression de radiation ou la traînée atmosphérique. Ces développements en séries asymptotiques ou en variables canoniques ont conduit à la découverte des résonances orbitales, phénomènes clés pour expliquer la structure des anneaux planétaires, la migration des planètes dans les disques protoplanétaires, ou encore la stabilité relative du Système solaire. La formulation hamiltonienne et lagrangienne (La Mécanique) y occupe une place centrale, offrant un cadre géométrique symplectique qui préserve les invariants du mouvement et facilite l'analyse qualitative des trajectoires, des points de Lagrange aux variétés invariantes.

L'avènement de l'informatique a profondément renouvelé la discipline à partir du milieu du XX^e siècle. Les intégrations numériques à haute précision, couplées à des algorithmes symplectiques, des méthodes de régularisation et des techniques de filtrage d'erreurs, permettent désormais de simuler l'évolution de systèmes à N corps sur des échelles de temps cosmologiques. Ces simulations ont confirmé la nature chaotique sous-jacente à de nombreux systèmes astronomiques, montrant que la stabilité à long terme n'est pas une propriété déterministe absolue mais probabiliste, avec des horizons de prédictibilité limités par l'amplification exponentielle des incertitudes. 

Parallèlement, la mécanique céleste classique a dû être enrichie par la relativité générale, qui introduit des corrections mesurables dans les champs gravitationnels intenses ou à grande vitesse, comme la précession du périhélie de Mercure, la lenteur des horloges en orbite, ou l'évolution des pulsars binaires par émission d'ondes gravitationnelles, phénomènes aujourd'hui exploités pour tester les théories de la gravitation.

Aujourd'hui, la mécanique céleste demeure un champ interdisciplinaire en pleine expansion, à la croisée de l'astrophysique, des mathématiques appliquées et de l'ingénierie spatiale. Elle est indispensable à la conception des trajectoires de missions interplanétaires, au maintien de constellations de satellites, au suivi des débris orbitaux et à l'évaluation des risques d'impact d'astéroïdes géocroiseurs. Les recherches actuelles étudient la dynamique des disques circumstellaires, la formation hiérarchique des architectures planétaires, les mécanismes de capture et d'éjection dans les amas denses, ainsi que l'intégration de l'apprentissage automatique pour accélérer les calculs, identifier des résonances cachées ou optimiser des manoeuvres à faible poussée.

Malgré les avancées considérables, des questions théoriques majeures persistent, notamment la caractérisation rigoureuse de la stabilité à long terme des systèmes à N corps, la modélisation unifiée des effets non gravitationnels à différentes échelles, et l'articulation cohérente entre dynamique classique, relativiste et, dans des contextes extrêmes, quantique. La mécanique céleste continue ainsi d'incarner une entreprise scientifique ancienne mais toujours vivante : traduire, par le calcul et la validation observationnelle, l'architecture dynamique de l'univers.