Model Order Reduction of a Sliding Beam using a Global Basis: Formulation and Evaluation

Resumo Técnico: Redução de Ordem de Modelo de uma Viga Deslizante utilizando uma Base Global

Enunciado do Problema
Estruturas telescópicas, onipresentes em máquinas pesadas como guindastes e empilhadeiras, dependem de vigas deslizantes onde os comprimentos efetivos e as condições de contorno (BCs) mudam continuamente. A modelagem desses sistemas apresenta um desafio fundamental: as técnicas clássicas de Redução de Ordem de Modelo (MOR) assumem topologias fixas, computando formas modais uma única vez para uma configuração estática. Em vigas deslizantes, à medida que a estrutura se estende ou retrai, o comprimento efetivo e as BCs evoluem, tornando inválidas as coordenadas modais fixas. Atualizar a base modal durante a simulação para capturar essas mudanças normalmente requer a recomputação das coordenadas, o que anula os ganhos de eficiência computacional obtidos com a redução. Embora a PMOR (Redução de Ordem de Modelo Paramétrica) ofereça soluções, as abordagens locais existentes (que interpolam entre valores de parâmetros discretos) frequentemente resultam em modelos reduzidos que mudam com os parâmetros, complicando a interpretação física das coordenadas. Além disso, abordagens de PMOR globais, que projetam o sistema em um subespaço fixo, receberam pouca atenção no contexto da dinâmica de sistemas multibody com restrições para sistemas deslizantes.

Metodologia
Os autores propõem uma abordagem de PMOR global baseada em snapshots integrada dentro de um formalismo de multibody com restrições. A metodologia procede em duas etapas principais:

1. Construção de uma Base de Redução Global:

   • Snapshots: O sistema é discretizado usando elementos finitos de Euler–Bernoulli (EB) planares. Os snapshots são definidos como matrizes modais ($\Psi_i$) contendo um conjunto de modos próprios ($\phi$) computados para posições específicas do slider ($s$).
   • Parametrização: Dois parâmetros governam a geração dos snapshots: $n_{snap}$ (número de posições distintas do slider amostradas) e $r_{snap}$ (número de modos próprios retidos por snapshot). A BC deslizante é imposta algebricamente em nós específicos durante a computação desses modos próprios.
   • Compressão: Todos os snapshots são empilhados para formar uma matriz "pool" ($\Psi_{all}$). A Decomposição em Valores Singulares (SVD) via Decomposição de Valores Singulares (POD) é aplicada a este pool para extrair uma base de redução global compacta ($\Psi_{glob}$) consistindo nos vetores singulares esquerdos mais significativos.

2. Integração no Formalismo Multibody:

   • Imposição de Restrição: A viga deslizante é modelada usando uma estrutura de Lagrange Tipo 1 com restrições algébricas. O vetor de restrição $g$ impõe deslocamento zero na extremidade fixa e deslocamento transversal zero na posição do slider móvel. A posição do slider é dependente do tempo (reonômica) e holonômica.
   • Redução: As coordenadas completas do sistema $q$ são expressas como $q(t) = \Psi_{glob} p(t)$, onde $p$ são as coordenadas modais reduzidas. Esta projeção é aplicada às Equações Diferenciais Algébricas (DAEs).
   • Implementação: As matrizes de massa ($M_{red}$) e rigidez ($K_{red}$) reduzidas são constantes. No entanto, a matriz Jacobiana de restrição reduzida ($C_{q,red}$) e o vetor de força ($f_{red}$) são recomputados a cada passo de tempo para contabilizar o slider móvel e as cargas externas. A matriz Jacobiana de restrição é reduzida via $C_{q,red} = \Psi_{glob}^T C_q$.

Principais Contribuições

• Base Global em MBD com Restrições: O trabalho demonstra a integração de uma base de redução global em um formalismo de multibody com restrições para vigas deslizantes, uma combinação que não foi sistematicamente avaliada na literatura anterior.
• Formulação de Restrição Algébrica: Os autores formulam restrições deslizantes algébricas diretamente usando funções de forma de Hermite dentro da estrutura EB, garantindo continuidade $C^1$ conforme o slider se move entre os elementos.
• Avaliação Sistemática: Um estudo abrangente investiga o impacto da quantidade de snapshots ($n_{snap}$) e da contagem de modos próprios por snapshot ($r_{snap}$) na precisão, particularmente sob um orçamento fixo de modos globais ($r_{glob}$).
• Benchmarking: O método proposto é validado contra uma solução de referência de alta fidelidade usando a Formulação de Coordenadas Nodais Absolutas (ANCF) com uma implementação de junta deslizante, comparando tanto as respostas cinemáticas quanto a conservação de energia.

Resultos

• Precisão: Em um caso de teste desafiador envolvendo uma viga altamente flexível com carregamento dependente do tempo e movimento rápido do slider, a base de redução global (usando $n_{snap}=5, r_{snap}=6, r_{glob}=20$) manteve um erro de deslocamento RMS normalizado abaixo de 2% em comparação com a solução EB não reduzida. Especificamente, o erro foi de 1,82% para o deslocamento da viga e 1,57% para o deslocamento da ponta no caso de teste "evoluído".
• Conservação de Energia: A evolução da energia do modelo EB reduzido correspondeu ao benchmark ANCF, confirmando que a imposição da restrição algébrica não introduz nem dissipação nem injeção artificial de energia.
• Eficiência Computacional: A redução do espaço de coordenadas de 153 graus de liberdade para 20 resultou em uma redução no tempo de computação de aproximadamente 90% em comparação com a abordagem EB não reduzida. Comparado ao benchmark ANCF, a abordagem reduzida foi cerca de 50% mais rápida.
• Sensibilidade aos Parâmetros: O estudo revelou um compromisso na seleção de snapshots. Sob um orçamento fixo de modos globais, o equilíbrio entre o número de snapshots e os modos por snapshot é crítico. Poucos snapshots geram uma base mal informada, enquanto muitos snapshots (com poucos modos cada) podem levar a uma truncagem ineficiente e erros maiores. Um produto ótimo de $n_{snap} \times r_{snap} \approx 26$ foi identificado para o caso de teste específico.

Significância e Alegações
O artigo afirma que a base de redução global proposta aborda com sucesso o desafio da topologia evolutiva em vigas deslizantes sem sacrificar o significado físico das coordenadas modais. Ao projetar o sistema em um subespaço fixo derivado de um pool de snapshots, o método permite a simulação eficiente de sistemas multibody flexíveis com fronteiras móveis. Os autores afirmam que o método é adequado para simulação em tempo real de estruturas telescópicas, oferecendo um aumento significativo de velocidade mantendo uma precisão aceitável (erros < 2%). O trabalho estabelece uma base para a aplicação de PMOR dentro da dinâmica de multibody com restrições, preenchendo a lacuna entre as técnicas de redução clássicas e as complexidades de sistemas mecânicos dependentes de parâmetros. O trabalho futuro nota a extensão desta abordagem para formulações de referencial de corpo flutuante para lidar com grandes movimentos de corpo rígido.