A comprehensive anisotropic phase-field fracture formulation with a novel positive/negative projection scheme

Resumo Técnico: Uma Formulação Abrangente de Fratura por Campo de Fase Anisotrópica com um Novo Esquema de Projeção Positivo/Negativo

Definição do Problema
A teoria de campo de fase (PF - phase-field) surgiu como uma abordagem eficiente para simular padrões de fratura complexos sem caminhos predefinidos. No entanto, estender esses modelos para materiais frágeis anisotrópicos apresenta desafios significativos. Embora os modelos de PF isotrópicos frequentemente dependam de decomposições de energia de deformação (como divisões espectrais ou divisões volumétrico-desviatóricas) para evitar o crescimento de trincas não físico sob compressão, esses métodos frequentemente falham em manter a consistência termodinâmica ou a ortogonalidade de energia quando aplicados à elasticidade anisotrópica. Consequentemente, os esquemas de projeção positivo/negativo (PN) existentes frequentemente assumem que a tensão e a deformação são coaxiais — uma suposição válida para materiais isotrópicos, mas que não se mantém em geral para sólidos anisotrópicos, onde as direções dos principais tensões e deformações divergem. Consequentemente, há uma falta de um arcabouço rigorosamente consistente que simultaneamente lide com a elasticidade anisotrópica, garanta o equilíbrio termodinâmico e suprima o dano sob estados compressivos sem o uso de aproximações heurísticas.

Metodologia
Os autores desenvolveram um arcabouço de campo de fase abrangente para materiais frágeis anisotrópicos, implementado no software de elementos finitos comercial Abaqus via sub-rotinas de elemento de usuário (UEL) e de material de usuário (UMAT). O núcleo da metodologia envolve três componentes primários:

1. Arcabouço Variacionalmente Consistente: A formulação adere às leis de balanço fundamentais e às primeira e segunda leis da termodinâmica. Utiliza uma estratégia de solução monolítica para resolver simultaneamente as equações acopladas de deslocamento e campo de fase, garantindo uma convergência robusta.
2. Elasticidade Anisotrópica: O modelo incorpora a simetria de material trigonal (representativa de cristais do tipo quartzo) através da rotação do tensor de rigidez do sistema de coordenadas do cristal para o quadro global de elementos finitos. Isso permite capturar a rigidez e a resposta mecânica dependentes da orientação.
3. Novo Esquema de Projeção Positivo/Negativo Anisotrópico (PNP): Para abordar a não-coaxialidade de tensão e deformação em materiais anisotrópicos, os autores propõem um novo esquema PNP baseado na norma de energia complementar.
   • Estratégia de Mapeamento: O método mapeia os tensores de tensão e deformação anisotrópicos reais para um espaço isotrópico auxiliar e fictício usando tensores de transformação derivados de equivalência de energia.
   • Projeção: Neste espaço isotrópico auxiliar, onde tensão e deformação são coaxiais, resolve-se um Problema de Complementaridade Linear (LCP) para decompor a deformação em partes positiva (tração) e negativa (compressão). Isso garante que a decomposição seja ortogonal na norma de energia e termodinamicamente consistente.
   • Retro-transformação: As quantidades projetadas são mapeadas de volta para o espaço anisotrópico real. A função de degradação do campo de fase é aplicada apenas à contribuição de tração (positiva), enquanto a contribuição de compressão (negativa) permanece intacta, efetivamente suprimindo o crescimento de trincas sob compressão.

Principais Contribuições

• Divisão Anisotrópica Termodinamicamente Consistente: O artigo introduz uma extensão rigorosa do esquema PNP de Wu e Cervera para a elasticidade anisotrópica. Ao utilizar um mapeamento de equivalência de energia para um espaço isotrópico auxiliar, a formulação garante que a energia motriz da trinca esteja consistentemente associada à deformação de tração, mesmo quando tensão e deformação não são coaxiais.
• Implementação Monolítica no Abaqus: Os autores fornecem uma implementação prática usando um UEL para o sistema acoplado deslocamento-dano e um UMAT para visualização de pós-processamento. Esta abordagem aproveita a robustez do solver do Abaqus enquanto permite um comportamento constitutivo anisotrópico complexo.
• Verificação Abrangente: O arcabouço é verificado através de uma hierarquia de testes, incluindo testes de elemento único sob tração, compressão e cisalhamento; testes de patch com malhas distorcidas para garantir a robustez da formulação do elemento; e múltiplos benchmarks de fratura quase-estática (tração de entalhe único, compressão, cisalhamento e flexão em três pontos).

Resultados
Os resultados numéricos demonstram a eficácia da formulação proposta:

• Benchmarks Isotrópicos: O modelo reproduz soluções analíticas e resultados da literatura para benchmarks isotrópicos padrão (ex: tração de entalhe único e flexão em três pontos) com alta precisão. Ele suprime com sucesso o dano sob compressão pura e exibe respostas simétricas para carregamento de cisalhamento positivo e negativo.
• Comportamento Anisotrópico: Para materiais trigonais (cristais de areia de sílica), o modelo captura respostas distintas dependentes da orientação.
  • Rigidez e Resistência: As simulações mostram cargas de pico mais altas e respostas mais rígidas quando o carregamento se alinha com a direção "forte" cristalográfica (eixo c) em comparação com a direção "fraca" (eixo a).
  • Trajetórias de Trinca: O modelo prevê trajetórias de trinca dependentes da orientação. Na orientação forte, as trincas desviam da direção do entalhe inicial devido aos termos de acoplamento anisotrópico (ex: $C_{14}$), curvando-se para uma direção energeticamente preferida antes de reorientar próximo a bordas livres. Na orientação fraca, as trincas propagam-se de forma mais linear.
  • Compressão e Cisalhamento: A divisão baseada em LCP evita com sucesso a evolução de dano sob carregamento compressivo para ambas as orientações. Sob cisalhamento, o modelo captura a assimetria na resposta entre carregamento de cisalhamento positivo e negativo para a orientação forte, uma consequência direta da simetria trigonal e do acoplamento constitutivo, mantendo a simetria para a orientação fraca.
• Estudos Paramétricos: O estudo confirma que o refinamento da malha, a escala de comprimento de regularização e a taxa de liberação de energia crítica ($G_c$) influenciam a carga de pico e o comportamento de amolecimento pós-pico, consistentes com as expectativas da teoria de campo de fase.

Significância
O artigo afirma que este trabalho fornece um "arcabouço abrangente e inovador" para quantificar o comportamento de fratura de materiais frágeis anisotrópicos. Sua principal significância reside em preencher a lacuna entre a consistência termodinâmica rigorosa e a modelagem prática de dependência de direção. Ao estender o esquema PNP para a elasticidade anisotrópica através de um mapeamento de equivalência de energia, a formulação supera as limitações de divisões heurísticas anteriores que assumiam a coaxialidade tensão-deformação. Os autores afirmam que o modelo é capaz de representar fenômenos de fratura em sólidos frágeis dependentes de direção, capturando não apenas rigidez e cargas de pico dependentes da orientação, mas também trajetórias de trinca complexas impulsionadas pela anisotropia. O trabalho é apresentado como um passo em direção a simulações mais precisas de materiais como o quartzo e outras estruturas cristalinas, embora os autores notem limitações atuais à simetria trigonal, geometrias 2D e energia de fratura independente de orientação, sugerindo estas como alvos para extensões futuras.