High-order isogeometric phase-field fracture modeling: a numerical-experimental study

Quando engenheiros projetam pontes, aeronaves ou até mesmo ferramentas cotidianas, eles devem antecipar como os materiais falharão. A forma mais dramática de falha é a fratura, onde uma peça sólida se divide subitamente. Por décadas, cientistas confiaram em duas formas principais para estudar isso: testes físicos em laboratório e simulações computacionais. Testes físicos são provas incontestáveis de como um material se comporta, mas são caros, demorados e destrutivos. As simulações computacionais oferecem uma alternativa mais rápida e barata, permitindo que projetistas testem milhares de variações antes que um único pedaço de metal seja cortado. No entanto, um desafio persistente tem sido o abismo entre os dois. Muitos modelos computacionais são construídos para descrever materiais que quebram como vidro, enquanto o mundo real é cheio de metais que dobram e esticam antes de quebrar. Esse descompasso frequentemente força os engenheiros a escolher entre um modelo que é matematicamente simples, mas fisicamente impreciso, ou um que é complexo, mas difícil de executar.

Uma equipe de pesquisadores de universidades da Itália e da Alemanha explorou recentemente um meio-termo surpreendente. Eles investigaram se um modelo computacional projetado especificamente para materiais frágeis, semelhantes ao vidro, poderia prever com sucesso o comportamento de um metal dúctil, maleável. O objeto de estudo foi uma liga de alumínio específica, comumente usada na construção e manufatura, conhecida por sua capacidade de se deformar significativamente antes de quebrar. Os pesquisadores buscaram verificar se poderiam usar uma estrutura matemática "frágil" para simular um evento real "dúctil", focando especificamente em como as rachaduras começam e se propagam através de uma amostra de metal com entalhes pré-cortados. O objetivo deles não era recriar cada detalhe minúsculo do dobramento do metal, mas ver se o modelo poderia prever com precisão duas coisas críticas: a força máxima que o metal poderia suportar antes de falhar e o caminho exato que a rachadura percorreria ao rasgar o material.

O experimento começou com testes físicos em corpos de prova de metal em formato de "osso de cachorro", um formato padrão para testar a resistência do material. Os pesquisadores utilizaram uma liga de alumínio conhecida como EN AW-6060, que contém silício e magnésio. Eles prepararam duas versões do corpo de prova: uma com 4 milímetros de espessura e outra com 6 milímetros. Cada peça apresentava dois pequenos cortes triangulares, ou entalhes, posicionados ligeiramente fora do centro. Esses entalhes atuam como pontos fracos onde uma rachadura provavelmente começará. A equipe colocou esses corpos de prova em uma poderosa máquina de teste que os puxava até que quebrassem. Eles observaram que a espessura do metal alterava dramaticamente o resultado. O espécime mais espesso, de 6 milímetros, desenvolveu uma rachadura que viajou diagonalmente de um entalhe ao outro, conectando-os em uma única linha nítida. Em contraste, o espécime mais fino, de 4 milímetros, desenvolveu duas rachaduras separadas que se curvaram para longe uma da outra e nunca se encontraram. Essa diferença de comportamento é um exemplo clássico de como as restrições físicas da espessura de um material podem alterar seu modo de falha, um fenômeno bem conhecido na engenharia, mas difícil de capturar em modelos computacionais simples.

Para simular esses eventos, os pesquisadores recorreram a um método computacional sofisticado chamado análise isogeométrica, que utiliza curvas matemáticas avançadas para representar a forma do objeto de maneira mais precisa do que os métodos tradicionais baseados em grades. Eles empregaram um tipo específico de modelo conhecido como modelo de fratura de campo de fase (phase-field). Nessa abordagem, uma rachadura não é representada como uma linha aguda e irregular que o computador precisa rastrear manualmente. Em vez disso, a rachadura é tratada como uma zona suave e difusa onde o material está perdendo gradualmente sua resistência. Essa variável de "campo de fase" atua como um interruptor de intensidade (dimmer), reduzindo a capacidade do material de se manter unido de totalmente intacto para completamente quebrado. Os pesquisadores escolheram uma versão deste modelo que foi originalmente projetada para materiais frágeis, o que significa que assume que o material quebra sem dobrar ou esticar significativamente. Esta foi uma escolha ousada, pois o alumínio que estavam testando é dúctil, o que significa que ele estica e absorve energia através da deformação plástica antes de finalmente romper.

A equipe executou suas simulações usando as mesmas dimensões e propriedades de material dos testes físicos. Eles ajustaram um parâmetro fundamental em seu modelo, conhecido como comprimento interno, que controla a largura da zona de dano difusa ao redor da rachadura. Ao ajustar cuidadosamente esse valor, descobriram que podiam replicar os diferentes comportamentos observados no laboratório. Quando definiram o comprimento interno para um valor maior, a simulação produziu um padrão de rachadura que conectava os dois entalhes diagonalmente, correspondendo ao espécime mais espesso de 6 milímetros. Quando reduziram o comprimento interno, a simulação mostrou duas rachaduras curvas e separadas que não se conectavam, espelhando perfeitamente o comportamento do espécime mais fino de 4 milímetros. Esta foi uma descoberta significativa porque sugeriu que uma estrutura matemática única e relativamente simples poderia capturar fenômenos físicos complexos simplesmente ajustando um parâmetro, sem a necessidade de modelar explicitamente a complicada física do estiramento do metal.

Os resultados não foram apenas visualmente semelhantes; eles também foram numericamente precisos. O modelo computacional previu a força máxima que o metal suportaria antes de quebrar com uma precisão notável. Para o espécime mais espesso, a força prevista estava dentro de 1 por cento da força real medida. Para o espécime mais fino, a previsão estava dentro de 3,5 por cento. Esses números são excepcionalmente próximos para uma simulação que, intencionalmente, ignorou a complexa deformação plástica que ocorre nos metais reais. Os pesquisadores observaram que, embora o modelo não pudesse reproduzir o amolecimento lento e gradual do metal após atingir sua força máxima, ele se destacou ao prever o momento da falha e o caminho que a rachadura tomaria. Isso torna o modelo uma ferramenta poderosa para engenheiros que precisam conhecer os limites de segurança de um projeto antes de ele ser construído.

O estudo conclui que o uso de um modelo de fratura frágil para simular materiais dúctis pode ser uma estratégia válida e conservadora para o design de engenharia. Ao dizer "conservadora", os autores querem dizer que o modelo tende a prever a falha em cargas ligeiramente menores do que as que o material realmente suportaria, o que é uma abordagem segura para aplicações críticas de segurança. O sucesso desta abordagem depende da formulação matemática específica utilizada, que permite que a energia seja dissipada precocemente no processo, mimetizando a absorção de energia de um material dúctil mesmo sem calcular explicitamente o estiramento. Os pesquisadores reconhecem que este é um compromisso fenomenológico; o modelo não conta toda a história do que acontece dentro do metal após ele ceder, mas conta a parte mais importante: quando e como ele irá quebrar.

Este trabalho abre uma nova perspectiva sobre como podemos usar simulações computacionais para compreender o mundo físico. Sugere que nem sempre precisamos do modelo mais complexo para obter a resposta correta. Às vezes, um modelo mais simples, aplicado com as ferramentas matemáticas certas e uma compreensão cuidadosa de suas limitações, pode fornecer previsões precisas e confiáveis. Para a liga de alumínio testada, os pesquisadores demonstraram que um modelo de fratura frágil poderia prever com sucesso tanto a força necessária para quebrar a peça quanto o intrincado caminho que a rachadura seguiria, simplesmente ajustando a escala da zona de dano. Esta descoberta oferece um caminho prático para engenheiros que projetam com metais, permitindo-lhes usar simulações de ordem superior eficientes para garantir segurança e desempenho sem se perderem no imenso custo computacional de modelar cada detalhe microscópico da deformação plástica. O estudo é um testemunho do poder de combinar testes laboratoriais rigorosos com métodos numéricos avançados para resolver desafios de engenharia do mundo real.