Resumo Técnico: Sobre o Surgimento da Estrutura de Pêndulo na Locomoção com Múltiplos Contatos

Formulação do Problema

O Modelo de Pêndulo Invertido Linear (LIPM) é uma suposição onipresente no controle de locomoção de pernas, frequentemente imposta mais como uma restrição de modelagem do que surgindo naturalmente do objetivo do controlador. Embora existam extensões para altura variável ou momento angular, a forma de pêndulo é tipicamente adicionada "manualmente". Este trabalho investiga uma questão mais restrita: Qual função objetivo específica faz com que o padrão de força de pêndulo surja como solução ótima e onde essa imagem se desfaz?

Os autores analisam um pequeno Problema de Controle Ótimo Centróide (OCP) que penaliza a taxa de variação do momento angular ($\|\dot{H}_G\|$) ao lado de termos padrão para forças de contato e aceleração normal. O objetivo é determinar se a minimização da taxa de variação do momento angular leva naturalmente o sistema à variedade de pêndulo (onde $\dot{H}_G = 0$) e caracterizar os limites desse comportamento sob diferentes configurações de apoio (full-rank vs. $N=2$) e requisitos de tarefa.

Metodologia

O trabalho formula um OCP de Força Líquida (Net-Wrench OCP) (Definição 1), que minimiza uma função de custo composta pelos seguintes componentes:

1. Termo de Penalidade de Momento Angular: $\alpha \|\dot{H}_G\|^2$
2. Rastreamento de Tarefa: $\lambda \|\dot{H}_G - \dot{H}_{task}\|^2$ (para tarefas não relacionadas ao equilíbrio)
3. Termo de Penalidade de Aceleração Normal: $\beta (n^\top \ddot{c})^2$
4. Regularização de Força: $\gamma \sum \|f_i\|^2$

A análise ocorre em três etapas teóricas, validadas por experimentos numéricos em um quadrúpede de massa pontual e no robô Unitree Go1 no MuJoCo (tanto QP em Malha Aberta quanto MPC de Nível de Torque em Malha Fechada):

1. Análise do Apoio Full-Rank: Utilizando a Decomposição em Valores Singulares (SVD) da matriz Jacobiana de Momentos, os autores derivam o comportamento assintótico quando o peso do momento angular $\alpha \to \infty$.
2. Análise do Apoio $N=2$: Investigação do caso de rank deficiente (por exemplo, trote diagonal), onde o operador de momento tem rank 1, levando a resíduos inevitáveis e restrições do cone de atrito.
3. Análise Orientada por Tarefa: Investigação de como a adição de um termo de rastreamento de tarefa ($\lambda > 0$) desloca o ótimo da variedade de pêndulo.

Principais Contribuições e Resultados Teóricos

1. Surgimento da Estrutura de Pêndulo no Apoio Full-Rank

• Teorema 1: Para apoios full-rank, a solução ótima colapsa assintoticamente com uma taxa de $O(1/\alpha)$ em direção ao padrão de força de pêndulo ($\dot{H}_G \to 0$).
• Constante Geométrica: O fator de proporcionalidade é determinado pela SVD da matriz Jacobiana de Momentos, o que corresponde diretamente à geometria do polígono de sustentação.
• Conexão com Modelos Clássicos: Nesta variedade de pêndulo, o Ponto de Aplicação da Força Zero (ZMP) corresponde naturalmente ao ponto de apoio virtual, e a equação orbital da componente divergente do movimento (DCM) surge como consequência, não como suposição.

2. Limites Fundamentais no Apoio $N=2$ (Trote Diagonal)

• Teorema 5 (Limite Geométrico Inferior): Em apoios $N=2$, a matriz Jacobiana de Momentos tem rank 1. Uma taxa de variação de momento angular residual permanece ao longo da direção da linha que conecta os dois pés ($\hat{D}$), independentemente de qualquer redistribuição de força. Isso gera um limite inferior não nulo (piso) para $\|\dot{H}_G\|$.
• Teorema 6 (Quebra no Cone de Atrito): Uma segunda limitação surge do cone de atrito. Se a força necessária para compensar a componente de momento perpendicular estiver fora do cone de atrito, ocorre uma "quebra não suave" em uma aceleração horizontal crítica ($a_x^*$). Além desse ponto, o $\|\dot{H}_G\|$ alcançável excede estritamente o limite geométrico inferior.
• Implicação: No trote diagonal, a direção do pêndulo frequentemente cai fora do cone de atrito (por exemplo, $40,6^\circ$ contra $31^\circ$ para $\mu=0,6$), tornando a solução de pêndulo inadmissível e gerando um limite inferior independente de $\alpha$.

3. Desvio Orientado por Tarefa

• Teorema 7: Quando um termo de tarefa ($\lambda > 0$) é adicionado, a taxa ótima de variação do momento angular torna-se uma interpolação suave: $\dot{H}_G^* = \frac{\lambda}{\alpha + \lambda} \dot{H}_{task}^*$.
• Diagrama de Fase: O sistema transita suavemente de um regime "dominado pelo equilíbrio" (tipo pêndulo, $\lambda/\alpha \ll 1$) para um regime "dominado pela tarefa" (não tipo pêndulo, $\lambda/\alpha \gg 1$). Isso explica por que o mesmo robô pode parecer semelhante ao LIPM em equilíbrio estático, mas ser não tipo pêndulo durante manobras ágeis, sem trocar o controlador.

Resultados Experimentais

Os autores validaram essas afirmações usando o Unitree Go1 e um modelo de massa pontual:

• Verificação de Taxa (Massa Pontual): Com o aumento de $\alpha$, o desvio do padrão de força de pêndulo diminuiu por um fator de 227, com uma inclinação Log-Log de $-1$, confirmando a taxa $O(1/\alpha)$ (Teorema 1).
• Constante Geométrica (Go1 N=4): Em um apoio de quatro pés, a constante de escala empírica ($K_e = 8,4$) concordou dentro de $\sim 16\%$ com a previsão analítica derivada dos polígonos de sustentação dos pés ($K_a = 9,8$). Os valores singulares da matriz Jacobiana de Momentos coincidiram com os polígonos de sustentação físicos até quatro casas decimais.
• Verificação de Piso (Go1 N=2): No trote diagonal, $\|\dot{H}_G\|$ saturou em um limite geométrico inferior previsto de $0,2835 \, \text{m}^2/\text{s}^2$ e concordou com o limite analítico até quatro algarismos.
• Realidade em Malha Fechada: Em uma MPC de Nível de Torque em Malha Fechada, o limite inferior persistiu, mas foi mais alto ($\approx 1,71 \, \text{m}^2/\text{s}^2$) do que a previsão estática. Os autores atribuem essa lacuna a efeitos dinâmicos (impactos de inércia, transições de contato, erros de servomecanismo) não capturados na análise quase-estática.
• Identidade ZMP: Os experimentos numéricos confirmaram que, na variedade de pêndulo, o ZMP baseado em força converge para o ponto de apoio virtual, validando a relação teórica entre o LIPM e o ótimo do OCP.

Significado e Afirmações

O trabalho afirma tornar a relação entre LIPM e controle ótimo explícita e verificável. Ele argumenta que:

1. O LIPM é uma consequência e não apenas uma escolha: O padrão de pêndulo surge naturalmente quando a taxa de variação do momento angular é minimizada em apoios full-rank, desde que o peso $\alpha$ seja suficientemente alto.
2. Os limites são previsíveis: O colapso da suposição de pêndulo em apoios $N=2$ não é uma falha do controlador, mas uma restrição fundamental geométrica e de atrito. O trabalho fornece expressões fechadas para os limites inferiores e a aceleração crítica na qual as restrições de atrito se tornam ativas.
3. Visão unificada de Equilíbrio e Tarefa: A razão $\lambda/\alpha$ fornece um parâmetro contínuo que explica a transição entre comportamento focado no equilíbrio (semelhante ao LIPM) e comportamento focado na tarefa (ágil).

Os autores permanecem modestos quanto aos resultados em malha fechada, observando que, embora o comportamento qualitativo (a existência de um limite inferior) se mantenha em simulações dinâmicas, a lacuna quantitativa entre a teoria estática e a realidade em malha fechada é significativa e provavelmente devida a dinâmicas não modeladas (compliance, deslizamento, taxas de controle finitas). Eles não afirmam resolver essas lacunas dinâmicas, mas sim definir os limites estáticos dentro dos quais essas dinâmicas operam.