A Shifted Cohesive-Zone Method for Non-Interface-Fitted Meshes with Applications to Crystal Plasticity

1. Formulação do Problema

A simulação precisa de problemas de mecânica dos sólidos dominados por interfaces (por exemplo, contornos de grão, trincas, interfaces de fase) em microestruturas complexas é desafiadora. Os Modelos de Zona Coesiva (CZM) tradicionais exigem malhas ajustadas à interface (onde as arestas da malha se alinham perfeitamente com as interfaces do material).

• Limitações das Malhas Ajustadas à Interface: A geração de malhas quadrilaterais ou hexaédricas de alta qualidade que se conformem a microestruturas 3D complexas é computacionalmente cara e frequentemente resulta em elementos distorcidos (Jacobianos negativos) ou requer ferramentas de malhagem complexas.
• Limitações dos Métodos Existentes Não Ajustados à Interface: Métodos como XFEM ou Métodos de Elementos Finitos Imersos (IFEM) frequentemente exigem o enriquecimento do espaço de aproximação localmente, levando a espaços de funções não uniformes e complexidade de implementação.
• Lacuna: Há necessidade de um método que permita o uso de malhas padrão, não conformes (especialmente quadrilaterais/hexaédricas) para problemas de interface coesiva, mantendo a precisão e evitando o ônus da geração de malhas.

2. Metodologia: O Método de Zona Coesiva Deslocado (SCZM)

Os autores estendem o Método de Fronteira Deslocada (SBM) para formulações de zona coesiva. A filosofia central é impor condições de interface em uma interface substituta ($\Gamma_h$) definida pelas facetas da malha, em vez da verdadeira interface física ($\Gamma$).

Componentes Teóricos Chave:

• Construção da Interface Substituta: A verdadeira interface é aproximada pela coleção de facetas da malha ($\Gamma_h$) que intersectam a verdadeira interface.
• Formulação Fraca Deslocada:
  • A lei de tração-separação é imposta em $\Gamma_h$ em vez de $\Gamma$.
  • Correção Geométrica: Uma expansão de transporte de superfície de primeira ordem é utilizada para mapear a integral da verdadeira interface para a interface substituta. Isso inclui um fator de correção de área geométrica ($|n \cdot n_h|$) para contabilizar o ângulo entre a normal verdadeira ($n$) e a normal da faceta ($n_h$).
  • Correção Direcional: Uma adição crítica ao SBM padrão. Como a normal da faceta difere da normal verdadeira, surge uma componente de tração tangencial. Os autores derivam um termo de correção direcional ($\sigma \tau_h$) para contabilizar esse desalinhamento, garantindo consistência no vetor de tração.
• Aproximação do Salto de Deslocamento: O salto de deslocamento ($\llbracket u \rrbracket$) necessário para a lei coesiva é aproximado na localização substituta usando uma expansão de Taylor de primeira ordem: $\llbracket u \rrbracket|_{\Pi(x)} \approx \llbracket u \rrbracket + \llbracket \nabla u \rrbracket \cdot d$, onde $d$ é o vetor de distância até a verdadeira interface.
• Modelos Constitutivos: O arcabouço suporta modelos dependentes da história, especificamente Plasticidade Cristalina (viscoplasticidade de cristal único com grande deformação com endurecimento Voce) acoplada a várias leis de tração-separação (exponencial, bilinear de modo misto e modelos acoplados 3D).

Implementação e Algoritmos:

• Software: Implementado no framework de código aberto MOOSE, acoplado ao NEML2 para modelagem constitutiva.
• Classificação de Pontos: Um algoritmo de classificação de pontos aprimorado por PCA é desenvolvido para determinar a qual grão um elemento pertence. Em vez de lançamento de raios fixo, raios são lançados ao longo da componente principal de variância mínima da geometria da fronteira, reduzindo significativamente as verificações de interseção.
• Construção da Substituta: Elementos são atribuídos a um único ID de grão com base na fração volumétrica dominante.
• Verificação: Um procedimento de transferência de solução projeta os resultados do SCZM em uma malha ajustada à interface gerada (IFM) para visualização e comparação.

3. Contribuições Principais

1. Primeira Aplicação à Plasticidade Cristalina: Demonstra o SCZM para modelos constitutivos dependentes da história, um passo significativo além das aplicações anteriores de elasticidade linear.
2. Formulação Fraca Simplificada: Deriva uma forma fraca reduzida na interface substituta que simplifica a implementação de resíduos não lineares e matrizes tangentes consistentes.
3. Correção Direcional: Introduz um termo de correção necessário para malhas não ajustadas à interface para lidar com o desalinhamento normal, prevenindo oscilações de tensão e erros.
4. Classificação Consciente da Geometria: Desenvolve um algoritmo baseado em PCA que acelera a construção de domínios substitutos, tornando o método escalável para geometrias 3D complexas.
5. Disponibilidade de Código Aberto: Fornece a primeira implementação dentro do framework MOOSE, permitindo integração com módulos de física existentes.

4. Resultados Numéricos e Verificação

O método foi validado através de uma série progressiva de benchmarks:

• Solução Fabricada (MMS):
  • Demonstrou convergência de primeira ordem para casos gerais (limitada pela aproximação geométrica da interface).
  • Demonstrou convergência de segunda ordem para casos simplificados onde os termos de correção de campo desaparecem.
  • Abordagens padrão não ajustadas à interface (sem deslocamento) falharam em convergir.
• Benchmarks de Interface Única:
  • Elasticidade: Os resultados do SCZM corresponderam às soluções de referência ajustadas à interface (IFM) para forças de reação, liberação de energia e campos de tensão.
  • Plasticidade Cristalina: O SCZM com correção direcional correspondeu estreitamente às soluções de referência IFM para forças de reação e distribuições de tensão. A abordagem "ingênua" (sem deslocamento ou correção) mostrou erros significativos.
• RVEs Policristalinos (2D e 3D):
  • 2D: O SCZM em malhas quadrilaterais não conformes convergiu para a solução IFM conforme a resolução da malha aumentou. Capturou com precisão curvas globais de força-deslocamento, evolução de dano e concentrações locais de tensão (por exemplo, em junções triplos).
  • 3D: Aplicado a um RVE 3D de 15 grãos com elementos hexaédricos. O SCZM reproduziu os padrões de deformação, forças de reação e evolução de dano da solução IFM com alta fidelidade, apesar da geometria complexa dos grãos.
• Desempenho: O algoritmo de classificação aprimorado por PCA mostrou tempo de execução quase constante independentemente da complexidade da fronteira, oferecendo acelerações significativas em relação ao lançamento de raios por força bruta.

5. Significado e Conclusão

• Eficiência: O SCZM elimina a necessidade de geração de malhas conformes à interface, demorada, particularmente para microestruturas 3D complexas onde malhas hexaédricas são preferidas para precisão, mas difíceis de gerar.
• Precisão: Ao incorporar correções geométricas e direcionais, o SCZM alcança precisão comparável aos métodos ajustados à interface sem o ônus da malhagem.
• Versatilidade: O método é robusto em diferentes tipos de elementos (triangulares, quadrilaterais, hexaédricos) e modelos de material (elasticidade linear a plasticidade cristalina complexa).
• Trabalho Futuro: Os autores sugerem estender o método para incluir correções de lacuna geométrica de ordem superior (inspiradas no Gap-SBM) e acoplamento com refinamento adaptativo de malha (baseado em octree) para maior eficiência.

Em resumo, este artigo apresenta um arcabouço robusto, preciso e eficiente para simular falha de interface coesiva em microestruturas complexas usando malhas não conformes padrão, preenchendo uma lacuna crítica entre modelagem constitutiva avançada (plasticidade cristalina) e capacidades práticas de malhagem.