A 2D discrete model with a phase-field regularisation applied to dynamic crack propagation

Quando um material quebradiço como o vidro ou o aço temperado se quebra, ele não simplesmente parte ao meio. Em vez disso, uma rachadura percorre o sólido, frequentemente ramificando-se em uma teia complexa de fraturas que estilhaça o objeto em muitos pedaços. Prever exatamente para onde essas rachaduras viajarão e quão rápido elas se moverão é um dos desafios mais difíceis da engenharia. Se um modelo usado para simular essa destruição for falho, ele pode sugerir que uma ponte falhará em uma direção segura quando, na verdade, ela colapsará catastroficamente. Durante décadas, cientistas tentaram construir modelos computacionais que mimetizem a estrutura atômica dos materiais, tratando-os como uma grade de pequenos pontos conectados por molas. Embora esses modelos sejam poderosos, eles sofrem de uma falha peculiar: o caminho que a rachadura percorre frequentemente depende de como a grade do computador é rotacionada, e não da física do material em si. É como se uma rachadura em um pedaço real de vidro mudasse subitamente de direção simplesmente porque o engenheiro girou a tela do computador.

Para resolver este problema, Matías Braun, um pesquisador da Universidad Nebrija, na Espanha, desenvolveu uma nova maneira de simular a fratura dinâmica. Sua abordagem combina uma grade tradicional de pontos com um método de "campo de fase" (phase-field), uma técnica que trata a rachadura não como uma linha nítida que salta de ponto em ponto, mas como uma zona de dano suave e difusa que se espalha pelo material. Ao fundir essas duas ideias, Braun criou um modelo que permite que as rachaduras viajem em qualquer direção, independentemente de como a grade de computador subjacente está orientada. Ele testou este novo método em três cenários distintos: uma placa de aço atingida por um impacto de alta velocidade, uma placa com entalhe puxada até quebrar e um cilindro de metal explodindo de dentro para fora. Em todos os casos, o novo modelo produziu resultados que coincidiram com experimentos do mundo real e eliminaram o viés artificial causado pela grade do computador, oferecendo uma visão mais clara e confiável de como os materiais realmente falham.

O cerne do problema que Braun abordou reside em como os computadores lidam com o mundo invisível dos átomos. Em seu trabalho anterior, ele utilizou um modelo discreto onde o material era representado por uma rede triangular, um padrão de pontos conectados por linhas. Quando o material era tensionado, as conexões entre os pontos se quebravam, simulando uma rachadura. No entanto, como os pontos estavam arranjados em um padrão fixo, a rachadura era forçada a seguir as linhas da grade. Se a grade fosse rotacionada, a rachadura mudaria seu caminho, mesmo que o próprio material não tivesse mudado. Este é um problema significativo porque, na realidade, uma rachadura deve seguir o caminho de menor resistência determinado pelas propriedades do material, e não o ângulo arbitrário de uma malha de computador.

A solução de Braun foi introduzir uma regularização de campo de fase. Em vez de forçar a quebra de conexões específicas entre os pontos, ele atribuiu um valor a cada ponto que representa o quão danificado ele está. Esse valor varia de zero, significando que o material está perfeitamente intacto, a um, significando que está completamente quebrado. A rachadura torna-se uma região onde esse valor aumenta gradualmente, criando uma transição suave em vez de uma linha irregular. Isso permite que a rachadura se curve e se ramifique em qualquer direção, ignorando efetivamente a estrutura rígida da grade. Crucialmente, Braun calibrou este sistema usando duas propriedades fundamentais do material: sua resistência à tração, que é a força necessária para iniciar uma rachadura, e sua energia de fratura, que é a energia necessária para manter a rachadura crescendo. Ao ancorar a simulação nessas constantes físicas reais, o modelo garante que os resultados dependam do material, não da matemática.

Para provar que este novo método funcionava, Braun realizou três testes rigorosos. O primeiro foi o teste Kalthoff–Winkler, um experimento clássico onde uma placa de aço com um entalhe pré-cortado é atingida por um projétil em alta velocidade. Neste cenário, uma rachadura dispara do entalhe em um ângulo específico. Modelos anteriores frequentemente previam o ângulo errado dependendo de como a grade era rotacionada. Nas simulações de Braun, a rachadura viajou consistentemente em um ângulo de aproximadamente 70 graus, correspondendo perfeitamente às observações experimentais. Este resultado mante-se verdadeiro quer a grade estivesse alinhada com a placa ou rotacionada em 15 ou 30 graus. A velocidade da rachadura também permaneceu estável, viajando a cerca de 70 por cento da velocidade do som no material, independentemente da orientação da grade.

O segundo teste envolveu uma placa com entalhe de borda única sob tensão, uma configuração usada para estudar como as rachaduras crescem lentamente sob um puxão constante. Aqui, o objetivo era ver se a rachadura viajaria em linha reta do entalhe até a borda oposta. Modelos de grade tradicionais frequentemente forçavam a rachadura a fazer zigue-zagues ou desviar com base no alinhamento da grade. O modelo de Braun, no entanto, produziu uma rachadura reta para cada orientação da grade. A força necessária para quebrar a placa também foi consistente, variando em menos de dois por cento, independentemente de como a grade era rotacionada. Isso demonstrou que o modelo poderia prever com precisão tanto o caminho quanto a resistência do material sem a influência artificial da malha do computador.

O teste final e mais complexo foi a fragmentação de um cilindro espesso sujeito a uma explosão interna. Neste cenário, nenhuma rachadura única é guiada por um entalhe; em vez disso, muitas rachaduras devem nuclear espontaneamente e correr para fora, estilhaçando o cilindro em dezenas de pedaços. Este é um evento caótico onde a simetria da grade poderia facilmente ditar o padrão dos fragmentos. Quando Braun executou a simulação com a grade centrada no eixo do cilindro, as rachaduras formaram um padrão simétrico. No entanto, quando ele deslocou a grade ligeiramente para que ela não estivesse perfeitamente centralizada, as rachaduras distribuíram-se irregularmente e o número de fragmentos mudou. Isso confirmou que o modelo era sensível à configuração física, e não apenas à simetria da grade. A simulação previu entre 12 e 19 fragmentos, um número que caiu dentro da faixa de outros métodos avançados, mas foi inferior a alguns, provavelmente porque a zona de dano suave no modelo impediu a formação de rachaduras que estavam muito próximas umas das outras.

Ao longo destes testes, os pesquisadores rastrearam cuidadosamente a energia do sistema. Eles compararam a energia perdida conforme o material se quebrava contra a energia necessária para criar as novas superfícies de rachadura. Em alguns casos, houve uma pequena discrepância entre a energia calculada pelo modelo de dano e a energia realmente liberada pelas conexões quebradas. Esta é uma limitação conhecida da abordagem híbrida, onde a força motriz para a rachadura é calculada de forma diferente da energia das ligações quebradas. No entanto, o balanço de energia mecânica permaneceu preciso dentro de meio por cento, confirmando que a simulação era fisicamente sólida. Os resultados sugerem que, embora o modelo não seja perfeito, ele representa um passo significativo à frente na compreensão da fratura dinâmica.

As implicações deste trabalho estendem-se para além de apenas melhores gráficos de computador. Ao remover a dependência da orientação da grade, os engenheiros podem confiar que suas simulações de testes de colisão, falhas estruturais ou testes de materiais estão refletindo a verdadeira física do material. O modelo consegue unir a natureza discreta dos átomos e o comportamento contínuo dos objetos sólidos. Mostra que, com as ferramentas matemáticas certas, é possível simular o processo caótico e violento de quebra sem deixar que a grade do computador dite o resultado. O estudo conclui que esta abordagem oferece uma maneira robusta de prever como os materiais falham sob condições extremas, proporcionando um caminho mais claro para o design de estruturas mais seguras e resilientes.