Efficient strain-space hyperreduction in large-deformation solid mechanics

Resumo Técnico: Hiperredução Eficiente em Espaço de Deformação para Mecânica de Sólidos de Grandes Deformações

Definição do Problema
Técnicas de Redução de Ordem de Modelo (MOR) que operam no espaço de deformação demonstraram recentemente um desempenho superior na homogeneização computacional, particularmente para Elementos de Volume Representativo (RVEs). No entanto, sua aplicação tem sido amplamente restrita a este contexto específico. Estender essas técnicas para problemas gerais de mecânica de sólidos de grandes deformações apresenta desafios significativos. Diferente das formulações baseadas em deslocamento, os métodos baseados em deformação não podem impor diretamente condições de contorno (BCs) de Dirichlet, pois as deformações nos pontos de integração na fronteira não são explicitamente conhecidas, e os valores de contorno são prescritos como deslocamentos, não como deformações. Além disso, as abordagens padrão de espaço de deformação frequentemente dependem de BCs periódicas, que não são aplicáveis a problemas de macroescala com valores de contorno arbitrários.

Metodologia
Os autores propõem um framework generalizado para aplicar MOR de espaço de deformação a problemas com BCs de Dirichlet arbitrárias e parametrizadas. A metodologia central envolve três componentes principais:

1. Tratamento de Contorno via Levantamento (Lifting): Para satisfazer BCs de Dirichlet não nulas, os autores empregam uma estratégia de levantamento. Um campo de referência consistente com o contorno é computado offline para valores de deslocamento unitários em cada fronteira independente. O campo solução real é então decomposto neste campo de referência (escalonado pelos valores de contorno reais) e um "campo de flutuação". O campo de flutuação desaparece em todas as fronteiras de Dirichlet por construção, permitindo que a base reduzida (derivada de Decomposição de Modo Próprio, POD) satisfaça as BCs a priori. Esta abordagem transforma BCs inhomogêneas em homogêneas para o solver reduzido.

2. Formulações de Espaço de Deformação: O artigo generaliza três técnicas específicas de hiperredução para este novo cenário:

   • ECM (Método de Cubatura Empírica) de Espaço de Deformação: Reformulado para realizar a integração reduzida diretamente nos pontos de Gauss, contornando operações ao nível do elemento.
   • E3C (Cubatura de Cluster Empiricamente Corrigida) de Espaço de Deformaçao: Generalizada para problemas mecânicos arbitrários. O E3C otimiza os pontos de integração e os pesos no espaço de deformação, permitindo que sejam distintos dos pontos de Gauss da malha original.
   • EMSL (Amostragem de Material e Linearização Empírica) de Espaço de Deformação: Estendido além dos RVEs. O EMSL lineariza o comportamento do material em torno de deformações médias médias por cluster. Ao contrário de outros métodos que integram aproximadamente uma lei de material exata, o EMSL integra exatamente uma aproximação linear por cluster, resultando em um problema afim por passo de carga que não requer um solver Newton-Raphson iterativo durante a fase online.

3. Linha de Base de Comparação: Estes métodos de espaço de deformação são testados contra o método padrão de Amostragem e Pesagem Conservador de Energia (ECSW) no espaço de deslocamento, que serve como uma linha de base robusta e preservadora de estrutura.

Principais Contribuições

• Generalização de MOR de Espaço de Deformação: O trabalho estende com sucesso as técnicas de MOR de espaço de deformação (ECM, E3C, EMSL) da homogeneização computacional para a mecânica de sólidos de grandes deformações com BCs de Dirichlet arbitrárias.
• Novo Tratamento de Contorno: A introdução de uma formulação de campo de flutuação baseada em levantamento permite a imposição de condições de contorno de Dirichlet arbitrárias no espaço de deformação sem exigir dados de deformação de contorno.
• Benchmarking Abrangente: Os autores fornecem uma comparação rigorosa entre a hiperredução de espaço de deformação versus espaço de deslocamento em dois exemplos hiperelásticos envolvendo estimativa de parâmetros inversos, variando parâmetros de material e trajetórias de carga complexas (tração, cisalhamento e carga mista).

Resultos
Os métodos foram testados em dois problemas de placas hiperelásticas com furos elípticos:

1. Exemplo 1 (Deformação Moderada): Uma placa sob carregamento de tração com até 20% de deformação.
2. Exemplo 2 (Grande Deformação e Carga Complexa): A mesma geometria sujeita a tração (até 50%), cisalhamento (até 25%) e carga mista, com uma malha mais fina e maior custo computacional.

Achados de Desempenho:

• Precisão vs. Tempo de Execução: Em ambos os exemplos, os métodos de espaço de deformação superaram consistentemente o ECSW de espaço de deslocamento no compromisso entre precisão e tempo de execução.
• Acelerações (Speedups):
  • E3C: Alcançou acelerações de aproximadamente 10.000 vezes (Exemplo 1) e 30.000 vezes (Exemplo 2) em relação ao modelo de ordem total, mantendo alta precisão (ex: <1% de erro). O E3C geralmente forneceu a maior precisão para um determinado orçamento computacional entre os métodos de espaço de deformação quando suficientes pontos de integração foram utilizados.
  • EMSL: Demonstrou eficiência excepcional de tempo de execução, alcançando acelerações de até 100.000 vezes (Exemplo 1) e 300.000 vezes (Exemplo 2). O EMSL destacou-se quando os orçamentos online e offline eram muito limitados, muitas vezes superando o E3C em velocidade, mas às vezes ao custo de uma precisão ligeiramente inferior em regimes altamente não lineares em comparação ao E3C com muitos pontos de integração.
  • ECM: Desempenhou melhor que o ECSW, mas geralmente ficou atrás do E3C e do EMSL nas métricas específicas testadas.
• Custos Offline: Embora o E3C possua uma fase de treinamento offline mais cara devido à otimização iterativa, os autores observam que isso é gerenciável e significativamente menor que o tempo necessário para gerar dados de treinamento. O EMSL possui o menor custo offline.

Significância e Alegações
O artigo afirma que a metodologia proposta justifica novas pesquisas em técnicas de MOR de espaço de deformação além da homogeneização computacional. Os autores destacam que estas técnicas oferecem:

• Eficiência Superior: Reduções dramáticas no tempo de execução online (ordens de magnitude) mantendo alta precisão.
• Não Intrusividade: Os métodos requerem apenas dados de deslocamento, deformação e tensão (facilmente extraídos de solvers black-box) e uma rotina de material. Eles contornam operações ao nível do elemento e não requerem conhecimento dos detalhes de discretização (ex: tipos de elementos) ou vetores de resíduo.
• Aplicabilidade: O tratamento bem-sucedido de BCs de Dirichlet arbitrárias remove uma barreira importante para a aplicação dessas técnicas eficientes em problemas estruturais de macroescala.

Os autores concluem que, embora restem desafios para problemas altamente não lineares (como elastoplasticidade) e para a estimativa de deformações por cluster no EMSL, o desempenho demonstrado sugere que estes métodos podem permitir novas aplicações de engenharia em controle, otimização e ajuste de modelos de materiais, onde o custo computacional é atualmente proibitivo.