Log-linear Dynamic Inversion for Thrusting Spacecraft on SE2(3)

1. Problema e Contexto

O crescimento acelerado de constelações de satélites (como Starlink e OneWeb) exige um aumento significativo em manobras orbitais. Para garantir a segurança e a eficiência de combustível, os sistemas de controle modernos dependem cada vez mais de manobras de empuxo contínuo e precisas.

No entanto, os modelos clássicos de dinâmica relativa, como as equações de Hill–Clohessy–Wiltshire (HCW) e Tschauner–Hempel/Yamanaka–Ankersen (TH/YA), são baseados em linearizações locais em torno de uma órbita de referência. Esses modelos falham em capturar o acoplamento não linear entre translação e rotação que ocorre durante manobras com empuxo, onde a aceleração é comandada no corpo e rotacionada para o referencial inercial. Isso resulta em perda de precisão e validade durante manobras de longa duração ou acopladas à atitude.

O objetivo deste trabalho é fornecer uma descrição geométrica robusta das manobras de espaçonaves com empuxo, superando as limitações das linearizações locais tradicionais.

2. Metodologia

Os autores propõem uma formulação baseada na teoria de sistemas invariantes, utilizando o Grupo de Lie $SE_2(3)$ (o grupo das isometrias espaciais diretas duplas), que estende o grupo $SE(3)$ ao incluir a velocidade translacional como parte do elemento do grupo.

• Representação do Estado: O estado da espaçonave (posição, velocidade e atitude) é mapeado em um elemento $X \in SE_2(3)$.
• Erro Invariante: Em vez de subtrair vetores no espaço euclidiano (que não tem significado físico em variedades não lineares), o erro de rastreamento é definido como um erro invariante à esquerda: $\eta = X^{-1}\bar{X}$, onde $\bar{X}$ é a trajetória de referência. O erro logarítmico $\xi = \text{Log}(\eta)^\vee$ vive na álgebra de Lie associada ($\mathfrak{se}_2(3)$).
• Dinâmica Misto-Invariante: A dinâmica do sistema é expressa na forma $\dot{X} = (M - C)X + X(N + C)$. Os autores identificam que, exceto pelo termo de gravidade $g(p)$, o sistema seria estritamente "grupo-afim". A gravidade é a única obstrução à afinidade do grupo.
• Abordagem de Controle:
  1. Limitação da Não-Linearidade: Demonstra-se que o desajuste de gravidade ($g(p) - g(\bar{p})$) atua como uma perturbação que pode ser limitada, permitindo o uso de ferramentas de estabilidade linear sem compensação de gravidade.
  2. Inversão Dinâmica: Propõe-se uma lei de controle de inversão dinâmica que elimina completamente o termo de gravidade, resultando em uma dinâmica de erro estritamente log-linear.

3. Principais Contribuições

O artigo apresenta três contribuições fundamentais:

1. Formulação Geométrica: Demonstra que a dinâmica translacional e rotacional de uma espaçonave com empuxo pode ser embutida no grupo de Lie $SE_2(3)$ como um sistema misto-invariante, identificando a gravidade como o único obstáculo à afinidade do grupo.
2. Limite de Não-Linearidade da Gravidade: Deriva um limite matemático rigoroso para o erro induzido pela gravidade. Mostra-se que, para erros de rastreamento limitados e altitudes orbitais acima de um mínimo, o desajuste de gravidade atua como uma perturbação linear proporcional ao erro de posição. Isso permite o uso de ferramentas padrão de estabilidade e alcançabilidade (como LMIs) sem necessidade de compensação de gravidade.
3. Lei de Controle de Inversão Dinâmica: Apresenta uma lei de controle que cancela exatamente o termo de gravidade, transformando a dinâmica de erro fechada em um sistema estritamente log-linear e independente da trajetória de referência.

4. Resultados

Os resultados foram validados através de simulações numéricas em uma órbita altamente elíptica (tipo Molniya), um caso de teste rigoroso devido à grande variação no raio orbital e no gradiente de gravidade.

• Validação da Dinâmica Log-Linear: A evolução do erro calculada diretamente pela dinâmica log-linear (Eq. 18) foi comparada com a integração clássica das equações de Newtonianas. Mesmo com separações de posição crescendo de 219 m para 2.407 km e erros de velocidade de 1,97 km/s, a discrepância entre os dois métodos foi de apenas $1,6 \times 10^{-8}\%$.
• Verificação do Limite de Gravidade: O termo real de desajuste de gravidade foi comparado com o limite teórico derivado. O limite global calculado (11,6 m/s²) cobriu com segurança o desajuste máximo observado (2,85 m/s²) ao longo de toda a órbita, confirmando a validade da aproximação mesmo em órbitas elípticas extremas.
• Precisão: Os resíduos de posição e velocidade permaneceram abaixo de $4 \times 10^{-4}$ m e $4 \times 10^{-7}$ m/s, respectivamente, indicando que as diferenças são devidas apenas a tolerâncias de integração numérica e não a erros de modelagem.

5. Significado e Impacto

Este trabalho é significativo por várias razões:

• Unificação Geométrica: Clarifica como a dinâmica de espaçonaves com empuxo se encaixa no quadro de sistemas invariantes, oferecendo uma representação estruturalmente consistente para erros de rastreamento.
• Superioridade sobre Modelos Clássicos: Diferente das equações HCW/TH/YA, que exigem re-linearização frequente e falham em grandes desvios ou acoplamentos atitude-translação, a abordagem $SE_2(3)$ é globalmente válida (dentro do raio de injetividade do log) e lida naturalmente com o acoplamento.
• Aplicabilidade Prática: A capacidade de tratar a gravidade como uma perturbação limitada ou eliminável via controle permite o uso de técnicas de controle linear robustas (como LQR ou análise de conjuntos alcançáveis via LMIs) em cenários complexos de manobras de aproximação e rendezvous.
• Futuro: O trabalho abre caminho para o desenvolvimento de planos de manobra de $\Delta v$ impulsivos baseados em trajetórias de referência com empuxo dentro de um único quadro geométrico unificado.

Em resumo, o artigo fornece uma base teórica sólida e validada numericamente para o controle de espaçonaves com empuxo, substituindo aproximações locais por uma formulação geométrica global e aproximadamente linear.