Endowing variational phase-field fracture models with custom strength criteria

Este artigo propõe um novo arcabouço de fratura por campo de fase variacional que incorpora domínios elásticos arbitrários e critérios de resistência personalizados através da introdução de um potencial de dissipação dependente do estado, preservando, assim, a estrutura variacional enquanto permite o controle independente sobre a degradação elástica e a nucleação de trincas sob estados de tensão multiaxiais.

O essencial

Imagine que você esteja tentando prever como um pedaço de vidro ou concreto irá rachar quando você o aperta, estica ou torce. Por muito tempo, os cientistas usaram uma ferramenta matemática inteligente chamada método de "campo de fase" (phase-field) para simular isso. Pense neste método como um mapa meteorológico de alta tecnologia para rachaduras: em vez de desenhar uma linha nítida e irregular onde uma rachadura aparecerá, ele pinta uma zona suave e difusa que transita gradualmente de "sólido" para "quebrado".

No entanto, havia um problema com os mapas antigos. Eles eram como um terno de tamanho único. Eles assumiam que um material quebra da mesma forma, quer você o esteja puxando (tração) ou apertando (compressão). Na realidade, os materiais são exigentes. O concreto, por exemplo, odeia ser puxado, mas é bastante resistente quando apertado. Os modelos antigos não conseguiam distinguir facilmente essas diferentes "personalidades" de tensão sem quebrar as regras matemáticas que faziam o modelo funcionar.

A Nova Ideia: Um Terno Sob Medida

Os autores deste artigo propõem uma nova maneira de construir esses modelos. Eles chamam isso de "dotar modelos de fratura de campo de fase variacionais com critérios de resistência personalizados". Em termos simples, eles descobriram como dar a esses modelos de previsão de rachaduras um terno feito sob medida que se ajuste às regras específicas de qualquer material.

Aqui está como eles fizeram isso, usando uma analogia simples:

O Sistema de Duas Partes: O Paletó e o Escudo

Imagine que um material está vestindo duas camadas:

1. O Paletó (Energia Livre): Esta camada representa a rigidez do material. À medida que o material sofre danos (como um paletó ficando com buracos), ele fica mais fraco e menos rígido. Nos modelos antigos, o paletó e as regras para quando ele rasgava estavam colados. Se você mudasse o paletó, acidentalmente mudava as regras de quando ele rasgaria.
2. O Escudo (Potencial de Dissipação): Esta camada representa a resistência ou o "ponto de ruptura" do material. Ela decide exatamente quanta força é necessária para iniciar uma rachadura.

A Inovação:
Os autores perceberam que poderiam deixar o Escudo mudar sua forma dependendo de como você está empurrando ou puxando o material, sem bagunçar o Paletó.

• Modo Antigo: Se você quisesse que o material fosse mais forte em compressão do que em tração, teria que reescrever toda a matemática do paletó. Isso era confuso e muitas vezes quebrava a "estrutura variacional" (a lógica interna que mantém a simulação estável).
• Novo Modo: Eles tornaram o Escudo "dependente do estado". Isso significa que o Escudo pode parecer um círculo, uma oval ou um blob estranho, dependendo da direção da força.
  • Se você puxar o material, o Escudo pode ser pequeno (fácil de quebrar).
  • Se você apertar o material, o Escudo pode ser enorme (difícil de quebrar).
  • Crucialmente, o Paletó (rigidez) permanece exatamente o mesmo. Os dois agora são independentes.

O Mapa do "Domínio Elástico"

O artigo fala muito sobre o "domínio elástico". Imagine um mapa de uma zona segura. Enquanto as forças no material permanecerem dentro desta zona, o material estará seguro e não irá rachar.

• Nos modelos antigos, essa zona segura era sempre um círculo (ou semicírculo) perfeito e simétrico.
• Nos novos modelos, os autores podem desenhar essa zona segura de qualquer forma que desejarem.
  • Eles podem fazer uma Dupla-Elipse (como um formato de amendoim) para lidar com diferentes limites de estiramento versus compressão.
  • Eles podem fazer um cone Drucker-Prager (como um cone de sorvete) para modelar rochas e solos que se comportam de forma diferente sob pressão.
  • Eles podem fazer um formato Huber que permite que o material seja apertado sem rachar (não interpenetração), mas ainda assim quebre facilmente se for puxado.

Por Que Isso Importa (Segundo o Artigo)

Os autores testaram seu novo método com várias "receitas" diferentes (modelos M1 a M5). Eles simularam um disco de material sendo puxado e empurrado de diferentes ângulos.

1. Flexibilidade: Eles mostraram que poderiam criar um modelo onde o material quebra facilmente ao ser puxado, mas é muito forte ao ser apertado, e vice-versa, tudo isso mantendo a matemática limpa e estável.
2. Independência: Eles provaram que você pode ajustar a "rigidez" (o quanto ele dobra) e a "resistência" (quando ele quebra) separadamente. Antes, mudar um frequentemente forçava a mudança do outro.
3. Precisão: As simulações mostraram que as rachaduras começavam exatamente onde o mapa da "zona segura" personalizada dizia que deveriam, correspondendo a diferentes condições de carregamento complexas (como torcer e apertar ao mesmo tempo).

A Conclusão

Este artigo não afirma que vai curar doenças ou construir novas pontes imediatamente. Em vez disso, ele fornece um novo kit de ferramentas matemáticas mais flexível. Ele permite que cientistas construam simulações de computador que respeitem as regras específicas e peculiares de diferentes materiais (como concreto, rocha ou tecido biológico) sem quebrar as leis fundamentais da física que tornam as simulações confiáveis. É como atualizar de um mapa genérico e pré-fabricado para um GPS que pode desenhar rotas personalizadas para qualquer terreno que você apresentar a ele.

Resumo técnico

Resumo Técnico: Dotação de Modelos de Fratura por Campo de Fase Variacional com Critérios de Resistência Customizados

Definição do Problema
Embora a modelagem de fratura por campo de fase tenha se tornado uma teoria consolidada para capturar padrões complexos de fissuração, uma limitação significativa permanece na capacidade de prescrever arbitrariamente critérios de resistência e descrever consistentemente a evolução do domínio elástico correspondente sob estados de tensão multiaxiais. As abordagens existentes geralmente se enquadram em três categorias: (a) divisão da energia livre, que mantém a estrutura variacional, mas carece de flexibilidade para definir domínios elásticos arbitrários; (b) modificação direta das forças de condução ou resistência do campo de fase, o que oferece versatilidade, mas sacrifica a estrutura variacional; e (c) acoplamento com deformações plásticas ou não lineares, que preserva a consistência variacional, mas introduz variáveis internas adicionais e complexidade de formulação. Crucialmente, na maioria dos métodos existentes, o comportamento de degradação elástica e o critério de resistência estão fortemente acoplados, impedindo o ajuste independente da degradação de rigidez e da resistência do material.

Metodologia
Os autores propõem uma estratégia fundamentalmente diferente, formulada dentro do arcabouço energético de Materiais de Padrão Generalizado (GSM). A inovação central é a introdução de um potencial de dissipação dependente do estado. Diferente dos modelos padrão, onde o potencial de dissipação é constante ou depende apenas da variável de campo de fase, esta abordagem permite que o potencial de dissipação dependa do estado atual do material (especificamente, dos invariantes de deformação ou tensão).

Componentes metodológicos principais incluem:

• Separação de Mecanismos: A densidade de energia livre ($\psi$) é responsável exclusivamente por codificar a degradação elástica (perda de rigidez), enquanto o potencial de dissipação dependente do estado ($\phi$) codifica o critério de resistência e a forma do domínio elástico. Esse desacoplamento permite que essas duas propriedades distintas do material sejam controladas independentemente.
• Dependência da Energia de Fratura: A energia de fratura $G_c$ é tratada como uma função da direção de carregamento (representada por um ângulo $\vartheta$ no plano de deformação volumétrica-desviatória ou de tensão). Essa dependência é incorporada via uma função de fratura $G_f(\varepsilon, \alpha)$ dentro do potencial de dissipação.
• Consistência Variacional: Ao derivar as equações de evolução a partir de um balanço de energia global e de uma condição de estabilidade (Formulação Energética), o modelo proposto retém a estrutura variacional do problema, preservando vantagens teóricas e computacionais associadas.
• Implementação Numérica: Um esquema de solução escalonado (staggered) é empregado, resolvendo os problemas de equilíbrio e dano alternadamente. A função de fratura dependente do estado é congelada no passo de carga convergido anterior durante as iterações para aproximar a energia dissipada dependente do estado.

Contribuições Principais e Modelos de Materiais
O artigo demonstra a flexibilidade deste arcabouço construindo e analisando cinco modelos de materiais específicos (M1–M5) que combinam diferentes leis de degradação e critérios de resistência:

1. M1 (AT1 Padrão): Degradação elástica total (FED) com um critério de resistência isotrópico padrão.
2. M2 (Dupla Elipse): FED com um critério de resistência de dupla elipse ajustável, permitindo limites volumétricos e desviatórios arbitrários.
3. M3 (Drucker–Prager): FED com um critério de Drucker–Prager, adequado para materiais sensíveis à pressão.
4. M4 (Drucker–Prager com Não-Interpenetração): Degradação elástica parcial (PED) onde o módulo de compressibilidade volumétrica não é degradado, mas o critério de resistência de Drucker–Prager é recuperado.
5. M5 (Huber com Não-Interpenetração): PED combinado com um critério do tipo Huber, permitindo o ajuste independente dos limites volumétricos de tração e de compressão, além de garantir a não-interpenetração sob compressão.

Resultados
O artigo apresenta resultados analíticos e numéricos para um problema de disco biaxial sob várias direções de carregamento ($\vartheta$):

• Respostas Homogêneas: Soluções analíticas confirmam que os modelos reproduzem com sucesso evoluções distintas dos domínios elásticos. Por exemplo, M3 e M4 compartilham a mesma forma de domínio elástico, mas exibem respostas mecânicas diferentes devido à presença ou ausência de degradação do módulo de compressibilidade. O M5 demonstra um comportamento de degradação não simétrico, distinguindo entre deformações volumétricas de tração e compressão.
• Nucleação de Fissuras: Simulações numéricas de nucleação de fissuras sob carregamento multiaxial mostem que os modelos preveem com precisão os pontos de nucleação correspondentes às superfícies de resistência prescritas. Os resultados ilustram que o arcabouço proposto pode capturar comportamentos de nucleação de modo misto (ex: expansão isotrópica, tração uniaxial, cisalhamento puro, compressão uniaxial) consistentes com os critérios de resistência específicos definidos para cada modelo.
• Independência de Características: Os resultados validam que a degradação elástica (governada pela energia livre) e o critério de resistência (governada pelo potencial de dissipação) podem ser ajustados de forma independente, uma característica não facilmente alcançável em abordagens anteriores.

Significância e Alegações
Os autores alegam que este trabalho abre caminho para desenvolvimentos futuros ao fornecer um arcabouço flexível e variacionalmente consistente para incorporar critérios de resistência arbitrários em modelos de fratura por campo de fase. A principal significância reside na capacidade de:

1. Preservar a estrutura variacional do problema ao introduzir critérios de resistência complexos e dependentes do estado.
2. Desacoplar a degradação de rigidez da resistência, permitindo a modelagem de materiais de forma mais realista e controlável de maneira independente.
3. Construir sistematicamente modelos para materiais complexos (ex: materiais quase-frágeis, friccionais ou sensíveis à pressão) sem recorrer a termos de fonte não-variacionais ou variáveis internas excessivas.

O artigo conclui modestamente, observando que, embora o arcabouço lide bem com a nucleação e a evolução estável, capturar eventos de propagação de fissura instáveis exigirá a integração futura de estratégias de regularização temporal, tais como termos inerciais ou viscosos. O presente trabalho estabelece a base teórica para essas extensões.