Stability analysis of admittance control using asymmetric stiffness matrix

Resumo Técnico: Análise de Estabilidade de Controle de Admitância com Matriz de Rigidez Assimétrica

1. Problema e Motivação

Em tarefas de automação que envolvem contato rico (como montagem de precisão "peg-in-hole"), o ajuste da rigidez do sistema de controle é um fator crítico. Métodos de controle de impedância e admitância são comumente usados para garantir complacência ativa.

• Limitação Atual: A maioria das análises de estabilidade baseia-se em matrizes de rigidez simétricas e na abordagem de passividade. Embora a introdução de uma matriz de rigidez assimétrica possa expandir o intervalo de seleção de parâmetros e melhorar a capacidade de indução de movimento do robô, ela viola a condição de passividade.
• O Desafio: A parte assimétrica da matriz gera um campo de força rotacional (campo de "curl"), que pode excitar oscilações em espiral (spiral oscillations) e aumentar a energia do sistema, tornando impossível garantir a estabilidade apenas através de métodos baseados em passividade. Até este estudo, não havia uma prova formal de estabilidade para esses sistemas com matrizes assimétricas.

2. Metodologia

Os autores propõem uma análise de estabilidade baseada na localização de raízes (root locus) e na análise de autovalores do sistema, em vez de depender exclusivamente da passividade.

• Modelo Matemático:
  • O sistema utiliza um modelo de admitância planar ($F_{ext} = M\ddot{x} + D\dot{x} + Kx$).
  • A matriz de rigidez $K$ é decomposta em uma parte simétrica ($K_s$) e uma parte assimétrica ($K_a$).
  • A matriz de amortecimento $D$ é assumida como diagonal.
• Análise de Autovalores:
  • Os autores derivaram as equações de estado do sistema e analisaram os autovalores da matriz do sistema ($A$).
  • Demonstraram que a presença de componentes imaginários nos autovalores da matriz de rigidez é a causa direta das oscilações em espiral.
• Derivação de Condições de Estabilidade:
  • Utilizando o método de root locus, os autores mapearam o comportamento dos polos do sistema à medida que o coeficiente de amortecimento ($d$) varia.
  • Derivaram desigualdades matemáticas que definem os limites de parâmetros para garantir que todos os autovalores sejam reais e negativos, assegurando a convergência do sistema.
  • A condição crítica para evitar oscilações em espiral é que a parte assimétrica ($k_{am}$) satisfaça a desigualdade:
    $$|k_{am}| < \sqrt{k_{sm}^2 + \frac{1}{4}(k_{xm} - k_{ym})^2}$$
    Onde $k_{sm}$ é o componente simétrico não diagonal e $k_{am}$ é o componente assimétrico.
  • Foi estabelecido um limite inferior para o coeficiente de amortecimento ($d$) necessário para estabilizar o sistema, mesmo na presença de assimetria.

3. Principais Contribuições

1. Identificação da Causa das Oscilações: Clarificaram que as oscilações em espiral são excitadas quando a matriz de rigidez é assimétrica e possui autovalores complexos (com parte imaginária).
2. Condição de Ausência de Oscilação: Demonstraram que a condição para a ausência de oscilações em espiral é equivalente à condição de que os autovalores da matriz de rigidez sejam números reais.
3. Análise de Estabilidade sem Passividade: Provaram que, embora a passividade não possa ser garantida com matrizes assimétricas, a estabilidade (convergência) pode ser garantida ajustando os parâmetros de amortecimento com base na análise de autovalores e root locus.
4. Fórmulas de Projeto: Forneceram equações explícitas para calcular o amortecimento mínimo necessário ($d$) em função dos parâmetros de rigidez assimétrica e simétrica.

4. Resultados (Simulação e Experimento)

• Simulações:
  • Compararam a resposta temporal de sistemas com e sem oscilações em espiral.
  • Quando a condição de autovalores reais não foi satisfeita, o sistema apresentou oscilações persistentes mesmo com alto amortecimento ($\zeta > 1$), falhando em atingir o amortecimento crítico.
  • Quando a condição foi satisfeita (autovalores reais), o sistema convergiu sem oscilações em espiral, atingindo o amortecimento crítico.
  • Estudos de energia mostraram que, embora a passividade local possa ser violada (a energia pode aumentar temporariamente durante fases de baixa velocidade devido ao campo de curl), a energia total converge se os parâmetros de amortecimento estiverem dentro do limite derivado.
• Experimentos:
  • Um braço robótico foi programado para mover-se verticalmente contra uma placa metálica, seguindo a superfície no eixo Y.
  • Matriz Simétrica: O robô apresentou um erro de posicionamento de 8 mm e uma força de contato menor devido a deslocamentos indesejados no eixo Z.
  • Matriz Assimétrica (Triangular): Com parâmetros ajustados conforme a teoria proposta, o erro de posicionamento foi reduzido para 6 mm e a força de contato foi mantida mais alta, pois o movimento foi confinado principalmente ao eixo Y, demonstrando a eficácia da rigidez assimétrica para tarefas específicas.

5. Significado e Conclusão

Este trabalho é significativo porque expande o horizonte de projeto para controle de robôs em contato.

• Flexibilidade de Projeto: Permite que engenheiros utilizem matrizes de rigidez assimétricas para otimizar o desempenho em tarefas complexas (como alinhamento e montagem) sem sacrificar a estabilidade, desde que os parâmetros de amortecimento sejam calculados rigorosamente.
• Superação de Limitações: Oferece uma alternativa viável à restrição de passividade, que limita o espaço de parâmetros em sistemas multivariáveis.
• Futuro: Os autores sugerem que a teoria pode ser generalizada para sistemas com 3 ou mais graus de liberdade e estendida para controle de rigidez variável, possivelmente utilizando conceitos de "tanques de energia" para análise de estabilidade.

Em resumo, o artigo fornece a base teórica e prática para o uso seguro e estável de matrizes de rigidez assimétricas em controle de admitância, transformando uma fonte potencial de instabilidade em uma ferramenta de controle mais versátil.