Peridynamic modeling of the crack velocity dependence via an incubation time fracture criterion

Este estudo utiliza uma abordagem peridinâmica com um critério de fratura baseado em tempo de incubação para modelar os experimentos de Ravi-Chandar e Knauss sobre Homalite-100, revelando que as variações no fator de intensidade de tensão do Modo-I em velocidades de trinca constantes e o início da ramificação microscópica em velocidades mais altas fornecem novas perspectivas sobre a natureza da dependência da velocidade da trinca na fratura dinâmica.

O essencial

Imagine que você está assistindo a uma trinca correndo através de um pedaço de plástico frágil, como uma folha de Homalite-100. Nos velhos tempos da física, os cientistas pensavam que, se você soubesse a velocidade com que a trinca estava se movendo, poderia calcular exatamente quanto "esforço" (ou pressão) a estava empurrando para frente. Era como pensar que, se um carro está indo a 60 milhas por hora, o motor deve estar produzindo exatamente 100 cavalos de potência. Simples, não é?

Mas experimentos na década de 1980 mostraram que isso não era verdade. Às vezes, a trinca estava se movendo na mesma velocidade exata, mas a pressão que a empurrava era drasticamente diferente. Era como se dois carros estivessem ambos fazendo 60 milhas por hora, mas um tivesse um motor minúsculo e o outro um propulsor de foguete. Os cientistas ficaram perplexos: Por que a mesma velocidade tem "empurrões" diferentes?

Este artigo é uma história de detetive onde os autores usam um novo tipo de simulação computacional para resolver esse mistério.

A Ferramenta do Detetive: Peridinâmica

A maioria dos modelos computacionais de trincas é como uma corrente de dominós. Se um dominó cai, ele empurra o próximo. Mas, se um dominó está faltando (uma trinca), a corrente se quebra e a matemática fica travada.

Os autores usaram um método chamado Peridinâmica. Pense nisso não como uma corrente, mas como um enxame de abelhas. Cada abelha pode falar com todas as outras abelhas dentro de uma certa distância, mesmo que haja um vão no meio. Se uma abelha voar para longe (uma trinca se formar), as outras abelhas apenas param de falar com ela, mas o resto do enxame continua se movendo perfeitamente bem. Isso permite que o computador lide com quebras e trincas sem ficar confuso.

O Ingrediente Secreto: O "Tempo de Incubação"

A verdadeira descoberta neste artigo é como eles decidiram quando uma trinca deve realmente se romper.

Da maneira antiga, se a pressão ficasse alta o suficiente, o material se quebrava instantaneamente. Mas os autores usaram uma regra chamada Critério do Tempo de Incubação.

Imagine que você está tentando quebrar um galho seco. Você não apenas puxa e ele se quebra instantaneamente. Você puxa, segura ali por uma fração de segundo enquanto as fibras se esticam e enfraquecem, e então ele se quebra. Essa fração de segundo é o "tempo de incubação".

Os autores programaram seu enxame de computadores para lembrar dos últimos microssegundos de pressão. O material só se quebra se a média da pressão durante esse curto período de "incubação" for alta o suficiente. Isso leva em conta o fato de que os materiais precisam de um pouquinho de tempo para "decidir" se quebrar.

O Que Eles Encontraram

Eles executaram simulações das placas de plástico sendo puxadas para fora, assim como os experimentos reais. Aqui está o que eles descobriram:

1. O Quebra-Cabeça Velocidade vs. Pressão: Assim como nos experimentos reais, seu computador mostrou que, para a mesma velocidade de trinca, a pressão (Fator de Intensidade de Tensão) não era um único número. Era uma faixa. Às vezes era baixa, às vezes alta.
2. O Efeito "Ramificação Microscópica": Quando a trinca se movia lentamente, ela ia em linha reta. Mas quando acelerava (acima de 400 metros por segundo), começava a ficar instável. Ela começou a brotar pequenas trincas laterais microscópicas, como um galho de árvore se dividindo em galhinhos.
   • A Analogia: Imagine um corredor fazendo um sprint. Em uma trote constante, ele corre em linha reta. Mas quando ele corre a toda velocidade, ele começa a oscilar e fazer zig-zags ligeiros para manter o equilíbrio.
   • O Resultado: Essas pequenas "oscilações" (ramificações microscópicas) fizeram a leitura da pressão subir e descer selvagemente. Isso explicou por que a pressão não era única para uma dada velocidade; a trinca estava mudando fisicamente sua forma ligeiramente enquanto corria.

A Conclusão

O artigo conclui que a razão pela qual vemos diferentes valores de pressão para a mesma velocidade de trinca é porque a trinca não é uma linha suave e perfeita. É uma coisa caótica e viva que flutua.

• Em velocidades mais baixas: A trinca é estável, e a pressão é relativamente estável.
• Em velocidades mais altas: A trinca começa a "ramificar-se microscopicamente" (brotar pequenas trincas laterais). Esse caos faz com que a pressão salte para cima e para baixo, criando a dispersão vista nos experimentos.

Ao usar esse "enxame de abelhas" (Peridinâmica) combinado com o "período de espera" (Tempo de Incubação), os autores recriaram com sucesso a relação bagunçada e não única entre velocidade de trinca e pressão que experimentos do mundo real haviam mostrado por décadas. Eles provaram que o "ruído" nos dados não é um erro; é uma característica física real de como trincas de movimento rápido se comportam.

Resumo técnico

Resumo Técnico: Modelagem Peridinâmica da Dependência da Velocidade de Trinca via Critério de Fratura com Tempo de Incubação

Enunciado do Problema
Este estudo aborda uma questão central não resolvida na mecânica da fratura dinâmica: a dependência do fator de intensidade de tensão (SIF) instantâneo do Modo-I ($K_I$) em relação à velocidade de propagação da trinca ($\dot{a}$). A Mecânica da Fratura Elástica Linear Clássica (LEFM) e critérios dinâmicos baseados em Freund e Rosakis frequentemente sugerem uma relação única entre $K_I$ e $\dot{a}$. No entanto, dados experimentais de Ravi-Chandar e Knauss em placas de polímero Homalite-100 frágil demonstram uma dependência não única, onde uma única velocidade de trinca corresponde a uma faixa de valores de SIF. Além disso, em velocidades mais altas, ocorre ramificação microscópica (micro-branching), levando a uma dispersão significativa nas medições de SIF. O desafio reside em reproduzir numericamente essa não unicidade e as flutuações associadas sem depender de propriedades materiais dependentes da taxa ajustadas empiricamente.

Metodologia
Os autores empregam uma estrutura peridinâmica não local para simular a propagação dinâmica de trincas em Homalite-100, replicando as condições experimentais de Ravi-Chandar e Knauss. A inovação metodológica central é a integração do Critério de Fratura com Tempo de Incubação (ITFC) na formulação peridinâmica.

• Estrutura Peridinâmica: As simulações utilizam um modelo peridinâmico baseado em estado ordinário implementado no software de código aberto Peridigm. O material é discretizado em pontos interagindo dentro de um horizonte esférico ($\delta$), permitindo o tratamento natural de descontinuidades como trincas sem operadores diferenciais espaciais.
• Critério de Fratura com Tempo de Incubação (ITFC): Diferentemente dos critérios clássicos que assumem falha instantânea, o ITFC postula que a fratura requer um tempo finito ($\tau$) para que processos preparatórios (por exemplo, microfissuração) ocorram. O critério avalia a tensão média no tempo sobre o período de incubação $\tau$ e um comprimento característico $d$.
  • O comprimento característico $d$ é derivado do critério quasi-estático de Irwin para garantir consistência com o SIF estático crítico ($K_{Ic}$).
  • A implementação peridinâmica adapta o critério de fratura de Tensão Remota (RS). A falha da ligação é determinada pela integração da tensão normal da ligação sobre o tempo de incubação $\tau$ (Eq. 11).
  • Para gerenciar a complexidade computacional e focar na ramificação microscópica em vez da ramificação macroscópica completa, o critério é avaliado apenas em planos paralelos ($0^\circ$) e perpendiculares ($90^\circ$) ao eixo principal da trinca.
• Configuração da Simulação: O modelo simula um corpo de prova de 300 mm $\times$ 500 mm com um pré-encaixe de 50 mm. O carregamento dinâmico é aplicado por meio de um pulso de força volumétrica trapezoidal com um tempo de subida de 25 $\mu$s. O SIF dinâmico é calculado usando um contorno retangular móvel para o integral J dinâmico, garantindo que o caminho de integração envolva a ponta principal da trinca enquanto exclui as ramificações microscópicas em desenvolvimento.

Principais Resultados
As simulações numéricas reproduzem com sucesso os fenômenos observados experimentalmente em relação à relação $K_I$-$\dot{a}$:

1. Dependência Não Única $K_I$-$\dot{a}$: Os resultados confirmam que, para uma velocidade de propagação de trinca quase constante, o SIF instantâneo exibe flutuações significativas. Consequentemente, o modelo prevê uma faixa de valores de SIF para uma dada velocidade, espelhando a não unicidade observada nos experimentos.
2. Estabilização de Velocidade e Transientes: Enquanto os experimentos mostram que a velocidade da trinca se estabiliza quase instantaneamente, as simulações exibem uma fase inicial de aceleração que dura aproximadamente 30–40 $\mu$s antes de atingir um estado estacionário. Durante essa fase transitória, o histórico do SIF mostra variações pronunciadas.
3. Ramificação Microscópica e Dispersão: Em taxas de carregamento mais altas (resultando em velocidades > 400 m/s), as simulações exibem o início da ramificação microscópica. Esse fenômeno correlaciona-se com um aumento marcado na amplitude das oscilações do SIF. A dispersão numérica nos valores de SIF em altas velocidades (0,53–1,25 MPa$\sqrt{m}$) alinha-se razoavelmente bem com as faixas experimentais (0,45–1,5 MPa$\sqrt{m}$).
4. Validação do Mecanismo: O estudo demonstra que a dispersão observada no SIF é fisicamente justificada pela dinâmica da ponta da trinca, especificamente a fase de aceleração e o início da ramificação microscópica, e não sendo um artefato de erro de medição.

Significado e Alegações
O artigo alega que a implementação peridinâmica do ITFC fornece uma abordagem valiosa para avaliar a tenacidade à fratura dinâmica e modelar a propagação de trincas em uma ampla faixa de taxas de carregamento. O significado principal reside na capacidade do modelo de capturar a sensibilidade à taxa da fratura frágil e a não unicidade da relação $K_I$-$\dot{a}$ sem introduzir parâmetros materiais dependentes da taxa arbitrários.

Os autores observam modestamente que, embora o modelo capture qualitativamente as tendências experimentais, discrepâncias quantitativas permanecem, particularmente em relação ao momento da estabilização da velocidade e à magnitude precisa das flutuações do SIF. Eles atribuem a previsão bem-sucedida da dependência não única à inclusão do tempo de incubação, que accounts pela história temporal da tensão, e à natureza não local da peridinâmica, que lida naturalmente com os campos singulares e descontinuidades associados à fratura dinâmica. O estudo conclui que a interação entre a fase de aceleração da trinca, as flutuações do SIF e o início da ramificação microscópica são os mecanismos governantes por trás do comportamento $K_I$-$\dot{a}$ observado experimentalmente.