Como calcular trajetória elíptica na prática
Muita gente tenta resolver problema de trajetória elíptica sem entender o que realmente está acontecendo. O erro mais comum é partir direto para a fórmula do período orbital sem antes verificar se o corpo está realmente em órbita elíptica ou se é só um ajuste numérico que precisa de correção. Eu já vi engenheiros perderem meia manhã com isso. O que define uma trajetória elíptica não é apenas a excentricidade estar entre 0 e 1. É a energia orbital total do sistema. Se a energia específica for negativa, o corpo está preso em órbita fechada. Se for zero, é parabólica. Se for positiva, é hiperbólica. Esse é o primeiro filtro que você deve aplicar antes de qualquer coisa.
O que é trajetória elíptica mesmo
Uma trajetória elíptica é o caminho que um corpo segue quando gravita em torno de outro corpo, mantendo distância variável ao longo do tempo. O corpo mais massivo fica em um dos focos da elipse, não no centro. Isso parece óbvio, mas muita gente coloca o foco na origem e depois se perde quando a simulação não fecha. Os parâmetros orbitais básicos são semieixo maior (a), excentricidade (e), e semi-latus rectum (p = a(1 - e²)). A partir desses três, você consegue derivar o resto. O período orbital vem da terceira lei de Kepler: T = 2(a³/), onde é o produto gravitacional padrão do corpo central. Em sistemas terrestres, vale aproximadamente 3,986 × 10¹ m³/s².
O que as pessoas não explicam bem é que a trajetória elíptica não é algo que se resolve com uma única equação. Você trabalha com duas coordenadas: a verdadeira anomalia () e a anomalia excêntrica (E). A relação entre elas exige solução numérica porque a equação de Kepler, M = E - e·sen(E), é transcendental. Não tem forma fechada para isolar E. Na prática, eu uso o método de Newton-Raphson para resolver a equação de Kepler. Começo com E = M para excentricidades baixas e E = para excentricidades altas. Geralmente convergi em três iterações, às vezes quatro se a excentricidade passar de 0,8. Para excentricidades acima de 0,95, o Newton-Raphson perde estabilidade e eu jogo direto no método de biseção com refinamento por regula falsi.
Um detalhe que quase ninguém menciona: se você está calculando trajetória elíptica para propulsão de foguetes, o problema muda completamente. Aí entra a equação de Tsiolkovsky, que não tem relação com Kepler. Misturar os dois quadros teóricos é um erro clássico que eu vi acontecer em relatórios de missões simples.
👉 Clique no botão abaixo para saber mais sobre o assunto!
Problema real que eu enfrentei e como resolvi
Eu estava trabalhando numa simulação de transferência orbital para um satélite de órbita baixa quando me deparei com um problema estranho. A trajetória elíptica calculada dava resultados coerentes nas primeiras passagens, mas a cada órbita a energia do sistema aumentava cerca de 0,03%. Depois de cinquenta revoluções, o satélite estava em uma órbita claramente divergente, o que era fisicamente impossível no modelo ideal. O culpado era o integrador numérico. Eu estava usando Runge-Kutta de quarta ordem com passo fixo de 60 segundos. Para órbitas com excentricidade próxima de 0,7, o passo era grande demais perto do perigeu, onde a velocidade muda drasticamente. O integrador "perdía" energia na precisão numérica e o sistema ganhava erro acumulado.
A solução foi trocar para um integrador adaptativo (Dormand-Prince 4(5)) com tolerância relativa de 1×10¹ e tolerância absoluta de 1×10¹². O passo passou a variar entre 2 segundos no perigeu e 300 segundos no apogeu. O erro de energia caiu para da ordem de 1×10¹ por órbita. Nada gratuito, claro. O tempo de simulação triplicou, mas pelo menos os resultados não pareciam ficção científica. Outro problema prático que eu enfrento frequentemente é a conversão entre elementos orbitais keplerianos e estados de posição e velocidade. A transformação direta funciona bem para excentricidades moderadas. Quando e > 0,9, a singularidade perto da órbita retílinea faz os cálculos explodirem. A solução que eu adoto é usar elementos equânimes de Vallado nesses casos, que removem a singularidade sem alterar a física do problema.
Limitações e quando a trajetória elíptica não funciona
Vale ser honesto aqui: trajetória elíptica como modelo puro só vale para problemas de dois corpos. Assim que você introduz terceiro corpo significativo — Lua, Sol, pressão de radiação solar, arrasto atmosférico — o modelo kepleriano deixa de ser suficiente. Para órbitas baixas abaixo de 600 km, o arrasto atmosférico modifica o semieixo maior em metros por dia. Sua "elipse" se fecha cada vez menor até o corpo reentrar. Para órbitas geoestacionárias, as perturbações solares e lunares empurram a excentricidade e a inclinação gradualmente. Sem manobras de correção, um satélite em GEO tende a desenvolver excentricidade de 0,15 a 0,2 em alguns anos. A trajetória elíptica teórica continua válida em cada instante, mas os parâmetros mudam continuamente.
Se você precisa de precisão orbital de longo prazo, o modelo de trajetória elíptica isolado não serve. Você precisa de um modelo de forças perturbadoras. Herrick-Ablaut e métodos de integração numérica com forças não-conservativas são o próximo passo. Não tente forçar Kepler a resolver o que ele não resolve. Para quem quer implementar isso do zero, o Python com a biblioteca poliastro cobre a maioria dos casos práticos. Ela já lida com a equação de Kepler, conversões de coordenadas e integrações adaptativas. O código leva menos de duzentas linhas para um simulador básico funcional. Eu recomendo começar por aí em vez de reimplementar tudo manualmente.
