An implicit octree-based adaptive Material Point Method

Resumo Técnico: Um Método de Ponto de Material Adaptativo Baseado em Octree Implícito

Definição do Problema
O Método de Ponto de Material (MPM) é um arcabouço robusto para simular grandes deformações, particularmente em interações entre estruturas contínuas e rígidas. No entanto, as implementações padrão do MPM dependem de malhas de fundo estruturadas, o que exige alta resolução em todo o domínio para capturar com precisão as interfaces em evolução entre corpos rígidos e contínuos deformáveis. Esse requisito leva a sistemas numéricos extensos e tempos de execução computacional excessivos, especialmente em três dimensões. Além disso, esquemas de integração temporal implícitos, que são vantajosos para problemas de engenharia de longa duração (ex: interação solo-estrutura), sofrem com a baixa escalabilidade de solvers lineares esparsos à medida que o número de graus de liberdade aumenta. O método também enfrenta instabilidades numéricas bem conhecidas: a instabilidade de travessia de célula (cell-crossing instability) e o problema do "pequeno corte" (small-cut problem), onde elementos são parcialmente cortados por fronteiras, levando a sistemas mal condicionados.

Metodologia
Este artigo apresenta o primeiro GIMPM (Método de Ponto de Material de Interpolação Generalizada) baseado em octree implícito. O núcleo da abordagem é uma malha de fundo adaptativa que refina e coarsifica (torna mais grossa) automaticamente com base na posição do corpo rígido que interage com o contínuo.

• Discretização por Octree: A malha de fundo utiliza uma estrutura de octree para concentrar o esforço computacional espacialmente. Elementos são refinados apenas em regiões que exigem alta fidelidade (próximo à interface), enquanto elementos grosseiros são mantidos em outros locais. Isso introduz nós pendentes (hanging nodes) (nós na face ou aresta de um vizinho mais grosseiro), que quebram a continuidade $C^0$ padrão.
• Formulação de Base GIMP Discreta: Para lidar com a malha de octree não conforme, os autores formulam uma base GIMP discreta. Em vez de construir funções de peça por peça complexas para todas as possíveis variações de interseção, as funções de base são computadas via uma regra de quadratura de ponto de Gauss única sobre a interseção do domínio do ponto de material e os elementos de fundo. Essa abordagem incorpora naturalmente as equações de restrição necessárias para ligar os nós pendentes aos seus nós de vértice de suporte, mantendo a continuidade $C^1$ das funções de base.
• Estabilização por Fantasma (Ghost Stabilisation): O problema do pequeno corte, que causa má condicionação em domínios não ajustados (unfitted), é mitigado pela extensão das técnicas de estabilização fantasma para malhas de octree. Isso envolve a adição de uma forma bilinear fantasma à matriz de rigidez através de "faces fantasmas" (interfaces entre elementos ativos e inativos ou de fronteira), garantindo estabilidade sem a necessidade de técnicas de agregação que poderiam complicar o arcabouço.
• Estratégia Adaptativa: A adaptividade é impulsionada pela posição do corpo rígido. A estratégia envolve três etapas: identificar elementos que interceptam o corpo rígido, refinar elementos dentro de uma distância especificada ($dx_{region}$) e refinar os pontos de material (GIMPs) com base no tamanho local do elemento. Os dados do material são transferidos de GIMPs grosseiros para os refinados, embora os autores notem que, para propriedades de material dependentes da posição, essa transferência depende de coordenadas de referência em vez de um arcabouço energético.

Principais Contribuições

1. Primeiro GIMPM de Octree Implícito: O artigo introduz o primeiro GIMPM utilizando uma malha de fundo de octree, permitindo a solução eficiente de problemas de interação contínuo-estrutura de grande deformação.
2. Formulação de Base Discreta: Um novo método de integração discreta para funções de base GIMP em malhas não conformes, evitando a necessidade de acoplamento de malha hierárquica ou definições complexas de peças.
3. Estabilização Fantasma Estendida: A adaptação da estabilização fantasma para malhas de octree com nós pendentes para resolver a instabilidade de pequeno corte.
4. Adaptatividade Automática: Uma estratégia sistemática que adapta tanto a malha computacional quanto a discretização do ponto de material com base na interface de contato, eliminando a necessidade de algoritmos separados de desrefinamento.

Resultos
O método foi validado através de cinco exemplos numéricos:

• Estudos de Convergência: Uma coluna sob peso próprio e uma coluna sob tensão de contato constante demonstraram convergência das soluções de tensão com o refinamento da malha, validando a formulação de nó pendente e a abordagem de penalidade de contato.
• Contato Dinâmico: Uma esfera rolando ladeira abaixo validou a implementação de fricção e a formulação de contato na presença de nós pendentes, mostrando excelente concordância com soluções analíticas para vários coeficientes de fricção.
• Aplicações Geotécnicas:
  • Teste de Penetração de Cone (CPT): A malha de octree adaptativa reproduziu com sucesso os resultados experimentais de CPT. O estudo destacou que, embora o refinamento aumente a rigidez neste setup específico (devido à definição do módulo de Young nos centros dos GIMPs), os resultados convergiram com o refinamento.
  • Penetração de Âncora: Uma âncora arrastada através de um leito de areia foi simulada e comparada com dados experimentais e uma simulação anterior de malha estruturada.

Desempenho e Impacto Ambiental
No caso da penetração da âncora, a abordagem de octree adaptativa alcançou resultados comparáveis ou superiores a uma malha conforme de alta resolução.

• Velocidade: O método adaptativo alcançou um aceleração de até 5,5 vezes em comparação com a simulação de malha estruturada.
• Eficiência: A redução nos graus de liberdade e nos pontos de material levou a uma diminuição significativa nos recursos computacionais.
• Pegada de Carbono: O artigo relata uma redução de 21 vezes nas emissões de CO2 para a simulação adaptativa em comparação com o equivalente de malha estruturada, calculada com base no uso de recursos de Computação de Alto Desempenho (HPC).

Significância
O artigo afirma que este trabalho representa um passo significativo para a eficiência do MPM implícito para problemas de grande deformação. Ao integrar a adaptatividade de octree com o GIMPM, o método concentra a resolução computacional estritamente na interface em evolução, reduzindo o número de graus de liberdade e pontos de material. Isso não apenas acelera as simulações, mas também oferece benefícios ambientais substanciais ao reduzir o consumo de energia necessário para a modelagem geotécnica de alta fidelidade. Os autores enfatizam que o método requer o mínimo de entrada do usuário (definindo o tamanho mínimo do elemento e a região de refinamento) e fornece um arcabouço robusto para interações complexas entre solo e estruturas rígidas.