Modelado Matemático del Problema de los Tres Cuerpos: De las Ecuaciones a las Órbitas Periódicas
El problema de los tres cuerpos consiste en calcular las trayectorias de tres objetos masivos bajo la ley de gravitación universal de Newton. Las condiciones iniciales especifican las posiciones y velocidades de los cuerpos en el espacio 3D, lo que lleva a un sistema de 18 ecuaciones diferenciales de segundo orden no lineales.
En el problema completo, las coordenadas y velocidades de cada cuerpo están ligadas a aceleraciones que dependen de sus distancias mutuas. Las masas normalizadas (relación con la máxima) simplifican los cálculos. Aquí está el sistema en forma de derivadas primeras:
X1' = VX1
Y1' = VY1
Z1' = VZ1
X2' = VX2
Y2' = VY2
Z2' = VZ2
X3' = VX3
Y3' = VY3
Z3' = VZ3
VX1' = -m2*(X1-X2)/((X1-X2)^2+(Y1-Y2)^2+(Z1-Z2)^2)^1.5-m3*(X1-X3)/((X1-X3)^2+(Y1-Y3)^2+(Z1-Z3)^2)^1.5
VY1' = -m2*(Y1-Y2)/((X1-X2)^2+(Y1-Y2)^2+(Z1-Z2)^2)^1.5-m3*(Y1-Y3)/((X1-X3)^2+(Y1-Y3)^2+(Z1-Z3)^2)^1.5
VZ1' = -m2*(Z1-Z2)/((X1-X2)^2+(Y1-Y2)^2+(Z1-Z2)^2)^1.5-m3*(Z1-Z3)/((X1-X3)^2+(Y1-Y3)^2+(Z1-Z3)^2)^1.5
VX2' = -m3*(X2-X3)/((X2-X3)^2+(Y2-Y3)^2+(Z2-Z3)^2)^1.5-m1*(X2-X1)/((X1-X2)^2+(Y1-Y2)^2+(Z1-Z2)^2)^1.5
VY2' = -m3*(Y2-Y3)/((X2-X3)^2+(Y2-Y3)^2+(Z2-Z3)^2)^1.5-m1*(Y2-Y1)/((X1-X2)^2+(Y1-Y2)^2+(Z1-Z2)^2)^1.5
VZ2' = -m3*(Z2-Z3)/((X2-X3)^2+(Y2-Y3)^2+(Z2-Z3)^2)^1.5-m1*(Z2-Z1)/((X1-X2)^2+(Y1-Y2)^2+(Z1-Z2)^2)^1.5
VX3' = -m1*(X3-X1)/((X1-X3)^2+(Y1-Y3)^2+(Z1-Z3)^2)^1.5-m2*(X3-X2)/((X3-X2)^2+(Y3-Y2)^2+(Z3-Z2)^2)^1.5
VY3' = -m1*(Y3-Y1)/((X1-X3)^2+(Y1-Y3)^2+(Z1-Z3)^2)^1.5-m2*(Y3-Y2)/((X3-X2)^2+(Y3-Y2)^2+(Z3-Z2)^2)^1.5
VZ3' = -m1*(Z3-Z1)/((X1-X3)^2+(Y1-Y3)^2+(Z1-Z3)^2)^1.5-m2*(Z3-Z2)/((X3-X2)^2+(Y3-Y2)^2+(Z3-Z2)^2)^1.5Variantes Simplificadas: Problemas Planar y Restringido
El problema planar (P3BP) se restringe al espacio 2D, reduciendo las ecuaciones a 12. Aquí, los cuerpos permanecen en el plano definido por las condiciones iniciales.
Para el problema planar:
X1' = VX1
X2' = VX2
X3' = VX3
Y1' = VY1
Y2' = VY2
Y3' = VY3
VX1' = m2*(X2-X1)/((X1-X2)^2+(Y1-Y2)^2)^1.5+m3*(X3-X1)/((X1-X3)^2+(Y1-Y3)^2)^1.5
VX2' = m3*(X3-X2)/((X2-X3)^2+(Y2-Y3)^2)^1.5+m1*(X1-X2)/((X2-X1)^2+(Y2-Y1)^2)^1.5
VX3' = m1*(X1-X3)/((X3-X1)^2+(Y3-Y1)^2)^1.5+m2*(X2-X3)/((X3-X2)^2+(Y3-Y2)^2)^1.5
VY1' = m2*(Y2-Y1)/((X1-X2)^2+(Y1-Y2)^2)^1.5+m3*(Y3-Y1)/((X1-X3)^2+(Y1-Y3)^2)^1.5
VY2' = m3*(Y3-Y2)/((X2-X3)^2+(Y2-Y3)^2)^1.5+m1*(Y1-Y2)/((X2-X1)^2+(Y2-Y1)^2)^1.5
VY3' = m1*(Y1-Y3)/((X3-X1)^2+(Y3-Y1)^2)^1.5+m3*(Y2-Y3)/((X3-X2)^2+(Y3-Y2)^2)^1.5El problema restringido (PR3BP) asume una masa pequeña para el tercer cuerpo. Ecuaciones para su movimiento:
X' = VX
Y' = VY
VX' = X+2*Y4-m1*(X+m2)/((X+m2)^2+Y^2)^1.5-m2*(X-m1)/((X-m1)^2+Y^2)^1.5
VY' = Y-2*Y3-m1*Y/((X+m2)^2+Y^2)^1.5-m2*Y/((X-m1)^2+Y^2)^1.5El problema restringido circular (CR3BP) modela una nave espacial en un sistema de dos cuerpos con el más pequeño en una órbita circular. Se usa el parámetro μ = m2 / (m1 + m2):
X' = VX
Y' = VY
z' = VZ
VX' = 2*VY+X-(1-mu) * (X+mu)/((X+mu)^2+Y^2+Z^2)^1.5-mu * (X-1+mu)/((X-1+mu)^2+Y^2+Z^2)^1.5
VY' = -2*VX+Y-(1-mu) * Y/((X+mu)^2+Y^2+Z^2)^1.5-mu * Y/((X-1+mu)^2+Y^2+Z^2)^1.5
VZ' = -(1-mu) * Z/((X+mu)^2+Y^2+Z^2)^1.5-mu * Z/((X-1+mu)^2+Y^2+Z^2)^1.5Estos sistemas se resuelven mediante métodos numéricos en Mathematica, Python (SciPy) o MATLAB.
Soluciones Analíticas: Euler y Lagrange
Las únicas soluciones analíticas conocidas son las soluciones colineales de Euler (1767) y las soluciones triangulares de Lagrange (1772). En el caso de Lagrange, los cuerpos comienzan en los vértices de un triángulo equilátero con velocidades coincidentes. Las trayectorias son elipses rotadas 120°; los cuerpos mantienen la configuración triangular.
Estas soluciones son inestables: pequeñas desviaciones en las condiciones iniciales llevan al caos. Las implementaciones prácticas requieren precisión hasta 10^{-10}.
Órbitas Periódicas en el Problema Planar
Henri Poincaré demostró la existencia de infinitas soluciones periódicas para masas iguales. Búsquedas computacionales han revelado familias con simetría.
- Familia Broucke-Enò-Hadjidimitriou (1975): Órbitas de 8 y 16 símbolos. Ejemplo R4 (período 5.4):
```
T0 = 0, T1 = 5.4
X1(T0) = 0.8733047091, X2(T0) = -0.6254030288, X3(T0) = -0.2479016803
Y1(T0) = 0, Y2(T0) = 0, Y3(T0) = 0
VX1(T0) = 0, VX2(T0) = 0, VX3(T0) = 0
VY1(T0) = 1.010776444, VY2(T0) = -1.683353346, VY3(T0) = 0.6725769022
```
- Familia Shuvakov-Dmitrashinovich (2013): Dragonfly (período 22):
```
T0 = 0, T1 = 22
X1(T0) = -1, X2(T0) = 1, X3(T0) = 0
Y1(T0) = 0, Y2(T0) = 0, Y3(T0) = 0
VX1(T0) = 0.08058, VX2(T0) = 0.08058, VX3(T0) = -0.16116
VY1(T0) = 0.58884, VY2(T0) = 0.58884, VY3(T0) = -1.17768
```
Otras: Butterfly, Moth, Goggles, Yin-Yang.
- Familia Shin (2016): Oval, cara de gato (período 5.1).
Para reproducirlas, usa Runge-Kutta de orden 8 con paso adaptativo.
Puntos Clave
- El problema completo requiere 18 ecuaciones; las simplificaciones lo reducen a 4–6 para simulaciones prácticas.
- Las soluciones analíticas de Euler/Lagrange son inestables y solo aplican en condiciones ideales.
- Las órbitas periódicas se encuentran mediante métodos numéricos; las familias Broucke y Shuvakov contienen cientos de variantes.
- Las condiciones iniciales deben especificarse con alta precisión para la periodicidad.
- CR3BP es crítico para misiones de NASA/ESA: maniobras en sistemas Tierra-Luna y Sol-Júpiter.
— Editorial Team
Aún no hay comentarios.