Modified augmented Lagrangian preconditioning for mixed-dimensional beam-solid coupling

Resumo Técnico: Pré-condicionamento de Lagrange Aumentado Modificado para Acoplamento Misturado de Dimensões Viga-Sólido

Enunciado do Problema
O artigo aborda os desafios numéricos associados ao acoplamento de dimensões mistas de corpos sólidos tridimensionais (3D) e vigas de Kirchhoff–Love de um único eixo (1D) livres de torção incorporadas. Esta abordagem de modelagem é essencial para simular materiais reforçados com fibras, como compósitos de fibras curtas ou tecidos biológicos, onde as fibras são significativamente esbeltas em comparação com a matriz sólida. O acoplamento é imposto via multiplicadores de Lagrange para garantir restrições posicionais exatas entre a linha central da viga e o sólido circundante.

Esta formulação leva a um sistema de equações lineares com uma estrutura de ponto de sela indefinida. O sistema apresenta dificuldades específicas:

1. Indefinição: A presença de um bloco zero na diagonal (associado aos multiplicadores de Lagrange) impede o uso de solvers padrão de definição positiva.
2. Subproblemas Singulares: O subproblema da viga, quando considerado isoladamente (como exigido por pré-condicionadores de bloco), é um problema de Neumann puro devido à ausência de condições de contorno de Dirichlet explícitas nas fibras incorporadas. Isso torna a matriz de rigidez da viga singular.
3. Sensibilidade aos Parâmetros: O condicionamento do sistema é fortemente influenciado pelos parâmetros de modelagem, incluindo a razão de rigidez entre viga e sólido ($E = E_B/E_S$), a razão de volume viga-sólido ($V$) e a razão de malha.
4. Limitações dos Métodos Existentes: Trabalhos anteriores dos autores utilizaram regularização por penalidade, o que reduz o tamanho do sistema, mas introduz má condição e inconsistência variacional. Por outro lado, abordagens diretas de multiplicadores de Lagrange exigem estratégias de pré-condicionamento especializadas que ainda não foram totalmente exploradas para este contexto específico de elasticidade de valores vetoriais e dimensões mistas.

Metodologia
Os autores propõem uma estratégia de pré-condicionamento de bloco de Lagrange Aumentado Modificado (MAL) para resolver o sistema de ponto de sela resultante. A metodologia procede da seguinte forma:

1. Formulação Aumentada: O sistema original de multiplicadores de Lagrange é reformulado em um sistema de Lagrange aumentado. Isso envolve a adição de termos de penalidade ($\epsilon M^T W^{-1} M$ e $\epsilon D^T W^{-1} D$) aos blocos de rigidez do sólido e da viga. Crucialmente, esta aumentação:

   • Preserva a imposição exata das restrições de acoplamento (ao contrário dos métodos de penalidade pura).
   • Regulariza o subproblema singular da viga, tornando a matriz de rigidez da viga invertível mesmo para pequenos parâmetros de penalidade.
   • Mantém a estrutura de ponto de sela do sistema global.

2. Construção do Pré-condicionador de Bloco: Partindo de um pré-condicionador de Lagrange aumentado ideal, os autores derivam variantes triangulares de bloco práticas. Estes pré-condicionadores permitem o tratamento independente dos blocos de sólido, viga e complemento de Schur:

   • Bloco do Sólido: Pré-condicionado usando métodos de Multigrid Algébrico (AMG) (especificamente Agregação Suavizada e Agregação Simples).
   • Bloco da Viga: Como a matriz aumentada da viga é bloco-diagonal (para fibras independentes) e invertível, ela é resolvida usando métodos diretos (fatoração LU).
   • Complemento de Schur: O núcleo do método envolve a aproximação do complemento de Schur aumentado. Os autores introduzem três variantes distintas:
     • Variante I: Aproxima o complemento de Schur usando inversos aproximados esparsos (SPAI) para o bloco da viga e uma aproximação diagonal para o bloco do sólido, assumindo um pequeno parâmetro de penalidade ($\epsilon \ll E_B$).
     • Variante II: Utiliza uma aproximação diagonal baseada na matriz de escala $W$ (especificamente $\hat{S} \approx -\frac{1}{\epsilon}W$), requerendo um parâmetro de penalidade cuidadosamente ajustado para agrupar os autovalores.
     • Variante III: Introduz parâmetros de penalidade separados para as contribuições do sólido ($\epsilon_S$) e da viga ($\epsilon_B$). Ao definir $\epsilon_S \to 0$ e ajustar $\epsilon_B$, esta variante visa desacoplar os problemas de escala, utilizando uma aproximação de complemento de Schur diagonal.

3. Implementação Numérica: Os sistemas lineares são resolvidos usando GMRES. A implementação utiliza os frameworks de software Trilinos e 4C, aproveitando o paralelismo de memória distribuída (MPI).

Contribuições Principais

• Novo Framework de Pré-condicionamento: O artigo apresenta a primeira aplicação de pré-condicionadores de bloco de Lagrange aumentado modificado especificamente adaptados para o acoplamento viga-sólido de dimensões mistas com multiplicadores de Lagrange.
• Tratamento de Subproblemas Singulares: A estratégia de aumentação regulariza efetivamente o subproblema singular da viga sem recorrer a técnicas complexas de projeção de nulidade ou métodos de deflação, permitindo o uso de solvers diretos padrão no bloco da viga.
• Variantes de Complemento de Schur: A derivação sistemática e a comparação de três estratégias de aproximação de complemento de Schur (incluindo a introdução de parâmetros de penalidade independentes para os domínios sólido e viga na Variante III).
• Análise de Escalabilidade: Avaliação abrangente de escalabilidade forte e fraca em clusters de memória distribuída para problemas com milhões de graus de liberdade.

Resultados
Experimentos numéricos em elementos de volume representativos (RVEs) e uma aplicação de placa sanduíche composta híbrida produziram os seguintes achados:

• Robustez: Os pré-condicionadores propostos demonstram robustez em relação à razão de rigidez ($E$), razão de volume ($V$) e refinamento de malha. A Variante III, em particular, mostrou quase-robustez de parâmetro para contrastes de rigidez moderados a altos.
• Independência de Malha: Os métodos alcançaram contagens de iteração quase independentes da malha (tipicamente entre 10 e 50 iterações para RVEs, e cerca de 20 para a placa composta) através de diferentes densidades de discretização.
• Escalabilidade: Os métodos exibiram escalabilidade forte e fraca favorável em clusters paralelos (até 512 processadores). As Variantes II e III mostraram eficiência paralela superior (até 97% de eficiência de escalabilidade fraca para o setup) em comparação com a Variante I, principalmente devido ao custo computacional da aproximação do complemento de Schur na Variante I.
• Aplicação de Engenharia: Na aplicação à placa sanduíche composta híbrida, o pré-condicionador ILU ingênuo falhou em convergir. Em contraste, os pré-condicionadores MAL propostos (especialmente a Variante IIIc) alcançaram convergência com baixas contagens de iteração e tempos de solução significativamente reduzidos (até 4x mais rápidos que a variante de pior desempenho) em comparação com abordagens padrão.

Significância
O artigo afirma que os pré-condicionadores de bloco de Lagrange aumentado modificado propostos fornecem uma alternativa robusta e eficiente tanto para formulações baseadas puramente em penalidade quanto para abordagens padrão de multiplicadores de Lagrange para simulações de grande escala de modelos de dimensões mistas em mecânica dos sólidos e estrutural. Ao permitir a imposição exata de restrições enquanto regularizam subproblemas singulares, o método facilita a simulação de estruturas reforçadas com fibras complexas, apresentando altos contrastes de rigidez e frações de volume variáveis. Os autores concluem que esta abordagem é particularmente adequada para aplicações de engenharia envolvendo fibras incorporadas, oferecendo um caminho viável para simulações de alta fidelidade e escaláveis de materiais compósitos. O trabalho futuro nota a extensão deste framework para formulações de vigas mais complexas (ex: Simo–Reissner) que incluem efeitos de torção e cisalhamento.