A variationally consistent beam-to-beam point coupling formulation for geometrically exact beam theories

Título

Uma formulação variacionalmente consistente de acoplamento pontual entre vigas para teorias de vigas geometricamente exatas

1. Problema e Motivação

Estruturas esbeltas do tipo viga são comuns em aplicações de engenharia (estruturas reforçadas com fibras, sistemas de treliça, mecanismos flexíveis, filamentos biológicos). Frequentemente, essas estruturas interagem em locais discretos através de juntas ou conectores.
O principal desafio abordado neste trabalho é a modelagem precisa dessas interações quando vigas diferentes são descritas por formulações numéricas distintas. Isso inclui diferenças em:

• Teorias de vigas subjacentes (ex: Simo-Reissner vs. Kirchhoff-Love).
• Esquemas de interpolação (ex: Hermite vs. Lagrange).
• Parametrizações de rotação (ex: vetores de rotação vs. tangente + torção).

Abordagens tradicionais de acoplamento direto (conectar graus de liberdade nodais) são limitadas a nós de elementos finitos e tornam-se complexas ou inviáveis quando as parametrizações das orientações diferem. Além disso, a imposição de restrições via métodos de penalidade pode levar a matrizes de sistema indefinidas e problemas de convergência.

2. Metodologia

O artigo propõe uma formulação de acoplamento pontual entre vigas independente da formulação, baseada na Teoria de Vigas Geometricamente Exatas (discretizada pelo Método dos Elementos Finitos).

Conceitos Chave:

• Variáveis Cinemáticas: O acoplamento é definido exclusivamente em termos da cinemática da seção transversal: posição do centróide ($r$) e orientação ($\Lambda \in SO(3)$).
• Medidas de Deformação Generalizadas:
  • Posicional: Introduz uma medida de deformação espacial que calcula a diferença entre a distância relativa atual e a distância relativa de referência, considerando um "push-forward" médio das duas seções para garantir simetria. Isso permite acoplar seções com centróides não coincidentes na configuração de referência.
  • Rotacional: Utiliza um vetor de rotação relativa baseado em rotações efetivas ($\hat{\Lambda}_i = \Lambda_i (\Lambda_i^0)^T$), garantindo que a restrição seja nula na configuração de referência livre de tensões.
• Imposição de Restrições:
  • Método dos Multiplicadores de Lagrange: Fornece uma estrutura variacional consistente, gerando um sistema de ponto de sela.
  • Regularização de Penalidade: Permite condensar os multiplicadores de Lagrange, eliminando variáveis adicionais do sistema, mas requer cuidados na escolha do parâmetro de penalidade para evitar má condicionamento.
• Propriedades Fundamentais: A formulação foi projetada para satisfazer:
  1. Objetividade: Invariância sob movimentos rígidos superpostos.
  2. Simetria: A formulação não depende da ordem de indexação das vigas (crucial para a regularização de penalidade).
  3. Consistência: A restrição é nula na configuração de referência livre de tensões.

Implementação Discreta:

A formulação é apresentada de forma geral, permitindo a integração direta em qualquer modelo de elementos finitos de vigas geometricamente exatas, independentemente das variáveis primárias ou esquemas de interpolação. As equações de acoplamento são mapeadas para os graus de liberdade discretos dos elementos finitos através de matrizes de forma.

3. Principais Contribuições

1. Generalidade: Uma abordagem unificada que permite acoplar vigas com formulações, interpolações e parametrizações de rotação completamente diferentes.
2. Flexibilidade Geométrica: O ponto de acoplamento não precisa estar em um nó de elemento finito; pode estar em uma posição arbitrária ao longo do elemento. O método lida corretamente com offsets (deslocamentos) iniciais entre os centróides das seções.
3. Consistência Variacional: Derivação rigorosa das contribuições virtuais de trabalho e das linearizações consistentes (matrizes tangentes) necessárias para esquemas de solução não linear (Newton-Raphson).
4. Validação de Propriedades: Demonstração matemática e numérica da objetividade, simetria e consistência com a configuração de referência.

4. Resultados Numéricos

O método foi validado através de quatro exemplos numéricos utilizando o código de pesquisa 4C e o pré-processador BeamMe:

• Estrutura em L (Offset Zero e Não-Zero):
  • Validou-se que a solução com multiplicadores de Lagrange coincide exatamente com uma conexão nodal direta (para offset zero).
  • Para offset não-zero, a formulação modela corretamente o acoplamento rígido, enquanto a regularização de penalidade converge para a solução de referência à medida que o parâmetro de penalidade aumenta.
• Duas Vigas Cruzadas (Convergência Espacial):
  • Analisou-se a convergência para malhas com número par e ímpar de elementos.
  • Observou-se que, quando o ponto de acoplamento coincide com um nó (malha par), a ordem de convergência é ótima ($O(h^2)$). Quando o ponto está no meio de um elemento (malha ímpar), a ordem cai para $O(h)$, o que é uma propriedade inerente às formulações de vigas exatas, não um defeito do acoplamento.
  • Confirmou-se a objetividade: a energia elástica total permanece constante durante rotações rígidas do sistema.
• Dupla Hélice (Dinâmica e Heterogeneidade):
  • Simulação dinâmica de uma estrutura de hélice dupla com 10 vigas conectadoras.
  • Demonstrou-se a capacidade de acoplar vigas Simo-Reissner e Kirchhoff-Love, com diferentes ordens de interpolação (Lagrange e Hermite) e parametrizações de rotação, em uma única simulação dinâmica.
• Cilindro Enrolado com Fios (Flambagem):
  • Modelagem de um cilindro formado por fibras circulares e axiais interligadas em 160 pontos.
  • O método capturou com sucesso o comportamento de flambagem não linear e a resposta força-deslocamento, demonstrando robustez em estruturas complexas com muitas conexões.

5. Significado e Conclusão

Este trabalho fornece um framework robusto e independente da formulação para o acoplamento de modelos de vigas geometricamente exatas. A principal vantagem é a capacidade de simular sistemas complexos onde diferentes tipos de elementos finitos (necessários para diferentes escalas ou comportamentos físicos) precisam interagir rigidamente sem a necessidade de compatibilidade de malhas ou parametrizações.

A formulação elimina as limitações de acoplar apenas em nós e resolve problemas de convergência associados a métodos de penalidade mal condicionados, oferecendo uma solução matematicamente consistente para a análise de estruturas esbeltas interconectadas. O trabalho abre caminho para futuras extensões, como acoplamentos linha-a-linha.