The cross-sectional warping problem for hyperelastic beams: An efficient formulation in Voigt notation

1. Problema e Motivação

A teoria clássica de vigas é frequentemente limitada a pequenas deformações elásticas e à suposição de que as seções transversais permanecem planas (não deformáveis). No entanto, muitas aplicações modernas, como robótica macia, materiais metamateriais e reforços de concreto, envolvem comportamentos não lineares, grandes deformações e materiais hiperelásticos.

Modelos tridimensionais (3D) de elementos finitos são precisos, mas computacionalmente caros para simulações estruturais de vigas. Alternativas unidimensionais (1D) baseadas em teoria de vigas são eficientes, mas muitas vezes falham ao capturar a deformação da seção transversal (warping) e comportamentos constitutivos não lineares complexos.

O Problema de Warping da Seção Transversal (CSWP - Cross-Sectional Warping Problem), introduzido anteriormente por Arora et al., oferece uma abordagem computacional para derivar relações constitutivas não lineares gerais para vigas, resolvendo um problema de valor de contorno 2D na seção transversal. No entanto, a formulação original e suas extensões utilizam uma formulação energética conjugada baseada no Tensor de Primeira Piola-Kirchhoff (PK1) e no gradiente de deformação. Isso exige a manipulação de tensores de segunda ordem não simétricos e envolve uma linearização complexa com tensores constitutivos de quarta ordem, o que aumenta o custo computacional e a complexidade da implementação numérica.

2. Metodologia Proposta

O artigo propõe uma reformulação totalmente material do CSWP para vigas 3D hiperelásticas e deformáveis por cisalhamento. Os pontos centrais da metodologia são:

• Formulação Material Simétrica: Em vez de usar o tensor PK1, os autores utilizam o Tensor de Segunda Piola-Kirchhoff (PK2) e o Tensor de Deformação de Green-Lagrange (GL). Ambos são simétricos, o que é uma vantagem fundamental.
• Notação de Voigt: Devido à simetria dos tensores PK2 e GL, a formulação pode ser expressa eficientemente na notação de Voigt (vetores e matrizes). Isso reduz tensores de segunda e quarta ordem para representações matriciais-vetoriais, simplificando drasticamente a implementação numérica e reduzindo o custo computacional.
• Formulação Variacional: O problema é formulado como um problema de minimização de energia com restrições. A função de warping ($u$) é determinada minimizando a energia de deformação total sobre a seção transversal, sujeita a restrições de movimento rígido (fixação do centro de massa e rotação) impostas via multiplicadores de Lagrange ($\lambda, \mu$).
• Discretização Isogeométrica: O problema é discretizado utilizando funções de base B-spline (Análise Isogeométrica - IGA), aproveitando a biblioteca de código aberto NLIGA.
• Cálculo de Sensibilidades e Rigidez: Para obter a matriz de rigidez da viga ($C_b$), são necessárias as sensibilidades da função de warping em relação às medidas de deformação da viga ($\varepsilon, \kappa$). Os autores utilizam um método adjunto derivado diretamente do sistema não linear discretizado para calcular essas sensibilidades de forma eficiente, evitando a necessidade de formulações de elementos finitos adicionais complexas.

3. Principais Contribuições

1. Nova Formulação PK2: Introdução de uma formulação energeticamente conjugada do CSWP baseada em tensões PK2 e deformações GL, adequada para vigas 3D hiperelásticas.
2. Operadores de Deformação-Deslocamento (B e $\bar{B}$): Derivação e uso de operadores diferenciais específicos para resolver o CSWP na notação de Voigt e determinar a rigidez tangente da viga.
3. Implementação Numérica Detalhada: Descrição completa da implementação numérica, incluindo a linearização do sistema e a construção da matriz de rigidez tangente.
4. Verificação Numérica: Validação rigorosa comparando a nova formulação "PK2" com a formulação original "PK1" (de Arora et al.) e com resultados estabelecidos na literatura.
5. Reprodutibilidade: Publicação de um repositório de código aberto (GitHub) contendo a implementação em MATLAB (baseada na biblioteca NLIGA) e todos os exemplos numéricos, permitindo que outros pesquisadores reproduzam e estendam os resultados.

4. Resultados

Os resultados numéricos demonstraram que a nova formulação é equivalente à formulação PK1 existente, mas com vantagens computacionais:

• Equivalência: As formulações PK1 e PK2 produziram resultados idênticos para deslocamentos de warping, forças resultantes, momentos e matrizes de rigidez.
• Convergência: O erro de resíduo ($L_2$-norm) e o número de iterações de Newton foram idênticos para ambas as formulações, confirmando a consistência matemática.
• Validação: O modelo foi validado em casos de torção pura, carregamento multiaxial e cisalhamento puro, utilizando diferentes modelos de materiais hiperelásticos (Saint Venant-Kirchhoff, Neo-Hooke e Mooney-Rivlin) e geometrias de seção transversal (quadrada e circular).
• Visualização: As distribuições de deslocamento fora do plano e tensões de von Mises mostraram-se consistentes com a literatura e com as expectativas físicas.

5. Significado e Impacto

Este trabalho preenche uma lacuna importante na literatura ao fornecer uma abordagem unificada e computacionalmente eficiente para obter modelos constitutivos de vigas que incluem materiais hiperelásticos e efeitos inelásticos (como inchamento e viscoelasto-plasticidade), mantendo a conformidade com a cinemática do contínuo 3D.

A principal vantagem da formulação proposta é a eficiência computacional e a facilidade de implementação devido ao uso da notação de Voigt e de tensores simétricos. Isso torna o método mais acessível para integração em frameworks de simulação existentes e abre caminho para extensões futuras para problemas mais complexos, como plasticidade e multiphysics (ex.: quimioelasticidade). A disponibilidade do código fonte fomenta a adoção mais ampla e o desenvolvimento de novos modelos constitutivos não lineares para vigas.