Time-Optimal Switching Surfaces for Triple Integrator under Full Box Constraints

Resumo Técnico: Superfícies de Comutação Otimizadas no Tempo para Triplo Integrador sob Restrições de Caixa Completas

Enunciado do Problema
Este artigo aborda o problema de controle ótimo no tempo para um sistema de triplo integrador sujeito a restrições de caixa completas. As dinâmicas do sistema envolvem entradas de controle limitadas por jerk ($u$) e restrições de estado na posição ($x_1$), velocidade ($x_2$) e aceleração ($x_3$). O objetivo é transferir o sistema de um estado inicial arbitrário $\mathbf{x}_0$ para um estado terminal fixo $\mathbf{x}_f$ no tempo mínimo $t_f$.

Embora o controle ótimo no tempo para duplos integradores seja bem compreendido, o caso do triplo integrador apresenta desafios significativos, particularmente ao lidar com:

1. Restrições assimétricas: Onde os limites superior e inferior para o jerk e os estados não são iguais em magnitude.
2. Condições de contorno não estacionárias: Estados iniciais e terminais arbitrários.
3. Restrições de posição ativas: Cenários onde a trajetória deve tocar ou deslizar ao longo das fronteiras de restrição de posição ($x_3 = \bar{x}_3$ ou $x_3 = \underline{x}_3$).

Soluções existentes, como o algoritmo comercial Ruckig (versão comunitária), frequentemente falham sob restrições assimétricas ou restrições de posição ativas. Abordagens analíticas anteriores ou ignoraram restrições de estado, assumiram restrições simétricas, ou falharam em fornecer uma caracterização completa das superfícies de comutação e das condições ativas para restrições de posição.

Metodologia
Os autores utilizam o Princípio do Máximo de Pontryagin (PMP) para caracterizar a lei de controle ótima, que segue uma estrutura Bang-Bang-Singular (BBS). O núcleo da metodologia envolve uma análise geométrica do espaço de estados para derivar superfícies de comutação otimizadas no tempo.

1. Lei de Comutação Aumentada (ASL): A trajetória ótima é representada como uma sequência de arcos (arcos bang-bang não restritos e arcos restritos onde os estados são fixos) e "marcadores tangentes". Um marcador tangente ocorre quando a restrição de posição torna-se ativa tangencialmente (ou seja, $x_3 = \bar{x}_3$, $x_2 = 0$ e $x_1 \geq 0$).
2. Construção da Superfície de Comutação:
   • Sem Restrições de Posição: Os autores constroem variedades de comutação por integração reversa a partir do estado terminal $\mathbf{x}_f$. Eles definem superfícies correspondentes a ASLs específicas (por exemplo, $\bar{0}\bar{0}\bar{0}$, $\bar{0}\bar{1}\bar{0}\bar{0}$) que particionam o espaço de estados 3D.
   • Com Restrições de Posição: O artigo deriva as condições específicas sob as quais as restrições de posição tornam-se ativas. Introduz novas superfícies de comutação (por exemplo, $\bar{0}(\bar{3}, 2)\bar{0}\bar{0}\bar{0}$) que contabilizam trajetórias tangentes à fronteira de posição. O Teorema 2 fornece uma classificação completa de todas as formas possíveis de ASL para o triplo integrador com restrições de caixa completas, estabelecendo que as trajetórias ótimas consistem em sequências de arcos não restritos e marcadores tangentes.
3. Implementação Algorítmica:
   • Algoritmo 1: Determina a direção relativa entre $\mathbf{x}_0$ e $\mathbf{x}_f$ analisando projeções de 2ª ordem, lidando com casos onde as superfícies de comutação são multivaloradas.
   • Algoritmo 2: Resolve o problema sem restrições de posição conectando estados via controles ótimos de 2ª ordem e resolvendo equações polinomiais para as durações dos arcos.
   • Algoritmo 3: Estende a solução para restrições de caixa completas. Se a solução não restrita violar as restrições de posição, o algoritmo conecta recursivamente o estado a um marcador tangente (derivado do Teorema 2) e resolve os subproblemas.
   • O sistema de equações para as durações dos arcos é resolvido usando bases de Gröbner, reduzindo o problema a uma única equação polinomial de grau $\leq 6$.

Principais Contribuições

• Caracterização Completa das Superfícies de Comutação: Pela primeira vez, o artigo fornece uma caracterização analítica completa das superfícies de comutação otimizadas no tempo para o triplo integrador sob restrições de caixa completas, incluindo a estrutura geométrica de soluções com restrições de posição ativas.
• Derivação das Condições Ativas: O artigo deriva explicitamente as condições sob as quais as restrições de posição tornam-se ativas (a condição de "marcador tangente"), uma lacuna anteriormente ausente na literatura.
• Teorema 2 (Classificação ASL): Estabelece o conjunto completo de Leis de Comutação Aumentada admissíveis para este problema, provando que as trajetórias ótimas seguem sequências específicas envolvendo marcadores tangentes.
• Algoritmo Eficiente: Um algoritmo é proposto que planeja trajetórias otimizadas no tempo para estados de contorno arbitrários e restrições assimétricas.

Resultados
O método proposto foi avaliado contra três baselines: Ruckig (versão comunitária), CasADi (baseado em IPOPT) e SCP (Programação Convexa Sequencial).

• Taxa de Sucesso: O método proposto alcançou uma taxa de sucesso de 100% na busca por soluções viáveis tanto para casos simétricos (Caso 1) quanto assimétricos com restrições de posição ativas (Caso 2). Em contraste, o Ruckig falhou no Caso 2 devido à assimetria, e o SCP/CasADi mostraram taxas de sucesso significativamente menores (34,1% e 97,7% respectivamente no Caso 2).
• Eficiência Computacional: O tempo computacional médio por trajetória foi de aproximadamente 10 $\mu$s (especificamente $2,38 \times 10^{-6}$ s para o Caso 1 e $1,25 \times 10^{-5}$ s para o Caso 2). Isso representa uma redução de 5 ordens de grandeza em comparação com baselines baseadas em otimização (CasADi: ~0,28s; SCP: ~7,28s).
• Optimalidade: O método encontrou consistentemente soluções globalmente ótimas. Em casos de teste onde métodos baseados em otimização convergiram para ótimos locais (ASLs subótimos), o método proposto reduziu o tempo terminal em até 90,9% (Caso 1) e 73,4% (Caso 2).

Significância
Os autores afirmam que este trabalho fornece a primeira caracterização completa de superfícies de comutação e um algoritmo eficiente para o problema de controle ótimo no tempo de 3ª ordem sob condições terminais não estacionárias e restrições de posição ativas. Ao afastar-se da otimização iterativa em direção a uma solução geométrica e analítica, o artigo demonstra que o controle de alta precisão e otimizado no tempo pode ser alcançado com tempos de computação na escala de microssegundos, tornando-o altamente adequado para aplicações industriais em tempo real, como movimento robótico e usinagem por controle numérico, onde as restrições são frequentemente assimétricas e ativas.