A multi-phase-field model for fiber-reinforced composite laminates based on puck failure theory

Este artigo propõe um modelo de campo de fase múltiplo baseado na teoria de falha de Puck, utilizando um método de sobreposição de malhas e duas variáveis independentes para prever com precisão os modos de falha e a interação entre danos em laminados de compósitos reforçados com fibras sob diversas condições de carregamento.

O essencial

Imagine que você está construindo um avião ou um carro de corrida usando materiais super leves e resistentes, como camadas de fibra de carbono coladas umas sobre as outras. Esses materiais são incríveis, mas têm um segredo: eles são como um "quebra-cabeça" de camadas. Se você puxar ou empurrar muito, elas podem se romper de maneiras muito diferentes, dependendo de como as fibras estão alinhadas em cada camada.

O problema é que prever exatamente onde e como esse material vai quebrar é um pesadelo para os engenheiros. Testar tudo fisicamente é caro, demorado e, às vezes, perigoso. É aí que entra este artigo.

Os autores criaram um super simulador de computador (um modelo matemático) que funciona como um "oráculo" para prever a falha desses materiais. Vamos descomplicar como eles fizeram isso usando analogias do dia a dia:

1. O Problema: O "Sanduíche" de Camadas

Pense no material composto como um sanduíche gigante feito de várias fatias de pão (as camadas).

• Em algumas fatias, as migalhas (fibras) estão alinhadas para a esquerda.
• Em outras, para a direita.
• Em outras, para cima e para baixo.

Se você apertar esse sanduíche, ele pode quebrar de duas formas principais:

1. Quebrando as fibras: Como se você estivesse tentando arrancar os fios de uma rede de pesca.
2. Quebrando a "cola" (matriz): Como se a massa entre as fibras estivesse rachando.

Antes, os computadores tinham dificuldade em simular isso porque tratavam o material como uma coisa só, sem distinguir entre a fibra e a cola.

2. A Solução: O "Duplo Detector de Fogo"

Os autores criaram um modelo que usa dois sensores independentes (chamados de "campos de fase") para monitorar o sanduíche:

• Sensor 1: Foca apenas nas fibras. Ele grita: "Ei, as fibras estão esticando demais!"
• Sensor 2: Foca apenas na cola (matriz). Ele grita: "A cola está rachando!"

Isso é como ter dois bombeiros diferentes em um prédio: um especializado em apagar incêndios na estrutura de aço (fibras) e outro especializado em apagar incêndios nas paredes de gesso (cola). Eles trabalham juntos, mas cada um sabe exatamente o que está acontecendo no seu domínio.

3. A Técnica Mágica: O "Papel Calco" (Mesh Overlay)

Aqui está a parte mais criativa. Como simular um sanduíche com 10 camadas sem deixar o computador lento?

• Eles usaram uma técnica chamada "Mesh Overlay" (sobreposição de malha).
• Imagine que você tem uma folha de papel quadriculado (a malha do computador). Em vez de desenhar 10 folhas de papel empilhadas (o que seria pesado), eles desenham uma única folha.
• Mas, em cima dessa mesma folha, eles desenham 10 "fantasmas" ou "papel calco".
  • No "papel calco 1", as fibras estão na horizontal.
  • No "papel calco 2", estão na vertical.
  • No "papel calco 3", estão em diagonal.
• Todos esses "papel calco" compartilham os mesmos pontos de conexão (nós), mas cada um calcula sua própria tensão e dano. É como se você tivesse várias camadas de vidro transparente empilhadas, e você pudesse ver a rachadura em cada uma delas sem precisar de 10 vidros físicos.

4. A Regra do Jogo: A Teoria de Puck

Para saber quando o sensor deve gritar, eles usaram uma regra famosa chamada Teoria de Puck.
Pense nisso como um "manual de instruções" que diz:

• "Se a fibra for puxada para a direita com força X, ela quebra."
• "Se a cola for espremida com força Y, ela racha."
• "Se houver uma mistura de torção e pressão, a falha será diferente."

O modelo segue esse manual rigorosamente para decidir quando o dano começa.

5. O Resultado: Previsão de Acidentes

Os autores testaram esse simulador em quatro situações diferentes, como se fossem testes de colisão em um videogame:

1. Teste de Esticão (Coupon): Puxar uma tira de material até quebrar.
2. Buraco no Meio (Open Hole): Um material com um furo no centro (como uma janela de avião) sendo puxado. O modelo previu exatamente onde as rachaduras se formariam ao redor do buraco.
3. Tensão Compacta: Um teste para ver como o material resiste a trincas que já existem.
4. Dupla Entalhe: Um material com dois cortes, simulando pontos fracos.

O que eles descobriram?
O simulador conseguiu prever com muita precisão:

• Qualitativamente: "Ah, a rachadura vai começar aqui e se espalhar para lá, formando um padrão X." (Eles compararam as imagens do computador com fotos reais de testes e eram muito parecidas).
• Quantitativamente: "O material vai aguentar até 25.000 Newtons de força antes de quebrar." (Os números batiam muito bem com os testes reais).

Resumo Final

Em termos simples, os autores criaram um simulador de "vidro mágico" para materiais compostos.

• Eles não precisam mais construir e destruir centenas de peças físicas para ver o que acontece.
• O computador consegue "ver" através das camadas, entendendo que a fibra e a cola têm personalidades diferentes.
• Isso ajuda engenheiros a projetar aviões e carros mais seguros, sabendo exatamente onde o material vai falhar antes mesmo de construí-lo.

É como ter uma bola de cristal que não apenas prevê o futuro, mas mostra exatamente como o material vai se despedaçar, permitindo que os designers mudem o projeto para evitar o desastre.

Resumo técnico

Aqui está um resumo técnico detalhado do artigo "A Multi-Phase-Field Model for Fiber-Reinforced Composite Laminates Based on Puck Failure Theory", apresentado em português:

Título: Um Modelo de Multi-Campo de Fase para Lâminas de Compósitos Reforçados com Fibras Baseado na Teoria de Falha de Puck

1. Problema e Motivação

Os compósitos reforçados com fibras (FRCs) são amplamente utilizados em setores como aeroespacial e automotivo devido à sua alta relação resistência/peso. No entanto, prever a fratura nesses materiais é um desafio significativo devido à complexidade dos mecanismos de falha, que incluem trincas na matriz, falha das fibras, delaminação e interações complexas entre essas modos.
A literatura existente sobre métodos de campo de fase (Phase-Field - PF) para compósitos apresenta limitações:

• Muitos modelos são restritos a lâminas unidirecionais ou sequências simples (ex: cross-ply).
• Dificuldade em capturar a falha sob carregamento compressivo devido a divisões de energia inadequadas.
• Falta de critérios robustos para a nucleação de trincas (iniciação da falha) que considerem a anisotropia e os diferentes modos de falha (fibra vs. matriz).
• A maioria dos modelos não consegue simular a propagação de trincas no plano (in-plane) para laminados com múltiplas camadas e sequências complexas sem custo computacional proibitivo (como em modelos 3D completos).

2. Metodologia

O artigo propõe um modelo numérico inovador que combina a Teoria de Falha de Puck com uma abordagem de Multi-Campo de Fase em duas dimensões, utilizando uma técnica de sobreposição de malha (mesh overlay).

• Formulação Variacional Multi-Campo de Fase:

  • O modelo utiliza duas variáveis de campo de fase independentes: $d_f$ (falha dominada por fibras) e $d_m$ (falha dominada por matriz/inter-fibra).
  • A energia total do sistema é decomposta em energia de deformação associada às fibras e à matriz, além da energia de superfície (fratura).
  • São introduzidos dois comprimentos característicos ($\ell_f$ e $\ell_m$) e dois tensores estruturais para descrever com precisão os modos de dano anisotrópicos.
  • A degradação da rigidez é governada por funções de degradação baseadas no modelo AT2.

• Critério de Iniciação (Teoria de Puck):

  • A iniciação do dano não é baseada apenas na energia, mas acionada pelos critérios de falha de Puck, que distinguem entre falha de fibra (tração/compressão) e falha de matriz (Modos A, B e C).
  • Fatores de exposição ($F_f, F_m$) são calculados com base no estado de tensão local. Quando atingem 1, a energia de acionamento para o campo de fase correspondente é ativada.

• Técnica de Sobreposição de Malha (Mesh Overlay):

  • Para simular laminados com diferentes orientações de fibras sem criar uma malha 3D complexa, cada camada (ply) é representada por uma malha separada (malha "dummy") que compartilha os mesmos nós da malha base.
  • O deslocamento é acoplado entre todas as camadas (garantindo compatibilidade), enquanto as tensões e o dano são desacoplados e calculados individualmente para cada orientação de fibra.
  • Limitação: O modelo é 2D e, portanto, ignora a delaminação (falha interlaminar), assumindo interfaces perfeitas entre as camadas.

• Implementação:

  • O modelo foi implementado como um elemento de usuário (UEL) no software comercial Abaqus, utilizando o solucionador Newton-Raphson.

3. Principais Contribuições

1. Integração Puck-Phase-Field: É uma das primeiras aplicações da teoria de falha de Puck dentro de um framework de campo de fase para capturar a nucleação e propagação de trincas em laminados complexos.
2. Modelo 2D Multi-Camadas: Desenvolvimento de um método eficiente computacionalmente (2D) capaz de simular a propagação de trincas no plano para laminados com múltiplas sequências de empilhamento, contornando a necessidade de modelos 3D pesados.
3. Separação de Modos de Falha: Uso de dois campos de fase independentes permite analisar a interação e competição entre falha de fibra e falha de matriz, algo difícil em modelos de dano contínuo tradicionais (CDM).
4. Validação Exaustiva: O modelo foi validado contra quatro conjuntos de dados experimentais distintos, cobrindo diferentes tipos de carregamento e geometrias.

4. Resultados e Validação

O modelo foi testado em quatro cenários de referência, comparando resultados numéricos com dados experimentais da literatura:

• Testes de Coupon (Tração e Compressão):

  • Simulou laminados com ângulos de 45°, 60° e 75°.
  • O modelo capturou corretamente o padrão de dano (trincas de matriz seguidas por falha de fibra) e as curvas força-deslocamento.
  • Houve uma boa concordância quantitativa nas tensões máximas, embora o modelo tenha ligeiramente superestimado a carga máxima em alguns casos devido à ausência de plasticidade.

• Tração com Orifício (Open Hole Tension - OHT):

  • Sequência [0°/90°]: O modelo previu com precisão o dano de matriz ao redor do orifício e a subsequente falha de fibra na direção perpendicular à carga, alinhando-se bem com os resultados experimentais de pico de carga.
  • Sequência [45°/-45°]: Capturou o fenômeno de "blindagem de trinca" (crack shielding) e a propagação assimétrica, demonstrando capacidade de lidar com comportamentos complexos de dano.

• Tração Compacta (Compact Tension):

  • Comparação com laminados de blocos cruzados ([0°/90°]).
  • O modelo reproduziu bem o comportamento linear inicial e a propagação de trincas. As discrepâncias no final da curva foram atribuídas à delaminação (não modelada) e ao puxamento de blocos de fibra nos experimentos reais.

• Tração com Dupla Entalhe (Double-Notched Tension):

  • Testado em sequências cruzadas, [45/0/-45] e quasi-isotrópicas.
  • O modelo demonstrou capacidade de simular o desenvolvimento de danos subcríticos e a interação entre modos de falha de fibra e matriz, replicando qualitativamente os padrões de fratura observados experimentalmente.

5. Significância e Conclusões

O trabalho apresenta uma ferramenta computacional robusta para a análise de falhas em compósitos laminados.

• Vantagens: Oferece uma alternativa computacionalmente eficiente aos modelos 3D, permitindo a análise de propagação de trincas no plano para laminados complexos. A separação dos campos de fase permite uma compreensão física mais profunda dos mecanismos de falha.
• Limitações: O modelo é estritamente 2D e não captura delaminação ou efeitos de tensões fora do plano. A não inclusão de plasticidade limita a precisão em regimes de pós-pico para alguns materiais.
• Impacto: O artigo fornece códigos e dados abertos, promovendo a transparência e permitindo que pesquisadores independentes validem e expandam o modelo. A abordagem proposta é um avanço significativo na modelagem de dano em compósitos, unindo a precisão da teoria de Puck com a capacidade de simulação de fratura do método de campo de fase.