Unveiling the Multiphysics Complexity: An Isogeometric Framework for Inducing Bifurcation and Tracing Post-Buckling Paths in Electroelastic Thin Shells

Resumo Técnico: Revelando a Complexidade Multifísica: Uma Estrutura Isogeométrica para Induzir Bifurcação e Rastrear Caminhos de Pós-Flambagem em Cascas Eletroelásticas Finas

Definição do Problema
Elastômeros dielétricos (DEs) são polímeros eletroativos capazes de grandes deformações sob campos elétricos, tornando-os promissores para atuadores suaves, sensores e colhedores de energia. No entanto, simular com precisão seu complexo acoplamento multifísico não linear permanece um desafio. Especificamente, os métodos de elementos finitos (FEM) convencionais baseados em polinômios de Lagrange lutam com os requisitos de continuidade $C^1$ da teoria de cascas de Kirchhoff–Love, que é essencial para modelar estruturas finas sem graus de liberdade rotacionais. Além disso, as estruturas numéricas existentes frequentemente carecem de robustez para rastrear bifurcações de quebra de simetria e caminhos de equilíbrio de pós-flambagem induzidos por carregamento eletromecânico acoplado, particularmente ao lidar com grandes deformações, incompressibilidade de material e não linearidades geométricas.

Metodologia
Os autores propõem uma estrutura abrangente de análise isogeométrica (IGA) baseada em superfícies de subdivisão Catmull–Clark para enfrentar esses desafios. A metodologia integra os seguintes componentes:

1. Formulação de Casca de Kirchhoff–Love Isogeométrica: A estrutura utiliza uma teoria de casca Kirchhoff–Love não linear adaptada para elastômeros dielétricos. Ela incorpora cinemática geometricamente não linear, acoplamento eletromecânico induzido pela tensão de Maxwell, incompressibilidade do material e pré-tensão inicial.
2. Superfícies de Subdivisão: Superfícies de subdivisão Catmull–Clark são empregadas para discretizar a geometria da superfície média e os campos de deslocamento. Esta escolha garante a continuidade global $C^1$ exigida pela teoria de Kirchhoff–Love, mesmo em malhas não estruturadas com vértices extraordinários, proporcionando simultaneamente representações geométricas suaves.
3. Modelagem Constitutiva e Decomposição de Tensão: A tensão total é decomposta em tensão mecânica hiperelástica, tensão de Maxwell (induzida eletricamente) e um termo de pressão hidrostática que impõe a incompressibilidade. A formulação impõe explicitamente a condição de tensão plana ($S_{33}=0$) através de um procedimento de condensação estática, derivando operadores tangentes consistentes analiticamente para eliminar o componente de tensão de espessura como uma variável independente.
4. Algoritmo de Rastreamento de Bifurcação: Para capturar instabilidades de quebra de simetria e rastrear caminhos de pós-flambagem, os autores empregam um algoritmo Newton–Raphson em estágios combinado com a continuação de arco esférico de Crisfield. Uma característica fundamental é o uso de perturbação de automodo: quando um autovalor zero é detectado na matriz de rigidez tangente, a solução é perturbada ao longo do respectivo autovetor para quebrar a simetria e guiar o solver para o ramo instável de pós-bifurcação.
5. Protocolo de Carregamento em Estágios: A implementação numérica segue um caminho de carregamento de três estágios para garantir convergência robusta: (1) carregamento de pressão da configuração não deformada até um nível intermediário; (2) carregamento de voltagem com pressão fixa; e (3) carregamento adicional de pressão com voltagem fixa.

Principais Contribuições
O artigo delineia quatro contribuições primárias:

1. Formulação de Casca Não Linear: Desenvolvimento de uma formulação de casca Kirchhoff–Love específica para elastômeros dielétricos que trata consistentemente a cinemática de deformação finita, o acoplamento eletromecânico, a tensão de Maxwell e a incompressibilidade do material.
2. Estrutura de Rastreamento de Bifurcação: Um framework numérico sistemático combinando continuação de arco esférico em estágios com perturbação de automodo para detectar robustamente instabilidades de quebra de simetria e atravessar ramos de equilíbrio instáveis em cascas eletroelásticas finas.
3. Discretização $C^1$-Contínua: A aplicação de superfícies de subdivisão Catmull–Clark para satisfazer os requisitos de continuidade da teoria de Kirchhoff–Love, garantindo campos de deformação suaves mesmo durante severa localização e autocontato em cenários de pós-flambagem.
4. Validação e Demonstração: Validação numérica abrangente contra soluções analíticas e benchmarks, demonstrando a capacidade do método em capturar grandes deformações, início de bifurcação, troca de modo e respostas de pós-flambagem.

Resultados
A estrutura proposta foi validada através de três exemplos numéricos:

• Cascas Esféricas Electroelásticas: O comportamento de inflação de uma membrana esférica foi simulado sob vários potenciais elétricos. Os resultados numéricos mostraram excelente concordância com soluções analíticas, capturando com precisão o efeito de amolecimento eletrostático (um deslocamento para a esquerda e para baixo na curva pressão–deslocamento com o aumento da voltagem).
• Placas Circulares Pré-tensionadas: O estudo investigou a interação entre pré-tensão mecânica e eletro-amolecimento. Os resultados indicaram que, embora os campos elétricos reduzam a rigidez estrutural e a capacidade de carga, o aumento da pré-tensão mecânica efetivamente suprime o efeito de eletro-amolecimento, aumentando a estabilidade da membrana.
• Membranas Toroidais: Este exemplo focou em quebra de simetria e comportamento de pós-flambagem. A estrutura rastreou com sucesso caminhos de equilíbrio complexos, revelando que os campos elétricos atuam como agentes desestabilizadores, deslocando pontos limites para pressões e deslocamentos menores. As simulações demonstraram que voltagens mais altas induzem abaulamento localizado severo e autocontato, suprimindo o ramo de equilíbrio simétrico e levando a modos de deformação assimétricos. O método lidou com sucesso com a transição de estados axissimétricos para não axissimétricos e capturou a evolução dos caminhos de pós-flambagem.

Significância e Alegações
O artigo afirma que sua estrutura isogeométrica fornece uma ferramenta robusta e precisa para analisar o comportamento não linear complexo de cascas finas de elastômero dielétrico. Ao superar as limitações de continuidade dos FEM convencionais e a dificuldade de rastrear caminhos instáveis, o método permite a simulação de instabilidades de quebra de simetria e respostas de pós-flambagem que eram anteriormente difíceis de resolver em geometrias de cascas finas sob carregamento eletromecânico combinado. Os autores posicionam este trabalho como um passo em direção ao design e otimização de atuadores suaves, sensores e colhedores de energia, observando que a estrutura é particularmente adequada para simular fenômenos não lineares envolvendo modos de deformação suaves e padrões de instabilidade complexos. Os autores notam modestamente que trabalhos futuros focarão na extensão da estrutura para incluir efeitos dinâmicos, viscoelasticidade e laminados multicamadas, mas não reivindicam aplicação imediata a essas áreas no estudo atual.