Finite compressibility and strain hardening in elasto-plastic models of amorphous matter

Título: Compressibilidade Finita e Endurecimento por Deformação em Modelos Elasto-Plásticos de Matéria Amorfica

1. Problema e Contexto

O comportamento plástico e de escoamento (yielding) de sólidos amorfos (como ligas metálicas amorfas, vidros moles e materiais granulares) é governado por transformações de cisalhamento locais. Após uma transformação, o material circundante sofre incrementos de tensão descritos pelo problema da inclusão de Eshelby. Um resultado fundamental da solução de Eshelby é que a tensão resultante dentro e ao redor da inclusão depende criticamente da compressibilidade do material, geralmente expressa pela razão entre o módulo de compressão ($K$) e o módulo de cisalhamento ($\mu$), ou pelo coeficiente de Poisson ($\nu$).

Apesar de trabalhos seminais (como o de Argon) terem destacado a importância da compressibilidade na barreira de energia para transformações de cisalhamento, a maioria dos tratamentos recentes de matéria amorfa baseados em modelos de automatos elasto-plásticos ignora a compressibilidade finita, tratando frequentemente o material como incompressível. Além disso, existe uma correlação empírica conhecida: ligas metálicas com alta razão $K/\mu$ tendem a ser mais dúcteis. O objetivo deste trabalho é preencher essa lacuna, investigando quantitativamente como a variação de $K/\mu$ afeta o comportamento macroscópico de escoamento e endurecimento em sistemas amorfos.

2. Metodologia

Os autores utilizam um modelo de automato elasto-plástico mesoscópico baseado em um mosaico (tiling) de sítios idênticos que podem sofrer transformações plásticas discretas.

• Discretização da Elasticidade: O modelo emprega uma discretização em grade escalonada (staggered grid) para resolver as equações de elasticidade linear. Esta abordagem, comum em sismologia, resolve problemas patológicos encontrados em esquemas de diferenças centradas de segunda ordem para termos de acoplamento de cisalhamento ($K_{12}$ e $K_{21}$).
• Parâmetro de Controle: A única variável de ajuste material é a compressibilidade adimensional $K/\mu$. O modelo introduz um parâmetro $\alpha = 2K/(K+\mu)$, que varia de 0 (comprimível máximo) a 2 (incompressível).
• Protocolos de Carregamento: O sistema é submetido a dois tipos de deformação:
  1. Cisalhamento Cíclico: Amplitude $\Gamma$, variando de baixas (regime elástico) a altas amplitudes.
  2. Cisalhamento Direto (Forward Steady Shear): Deformação monotônica até atingir um estado estacionário.
• Mecanismo de Escoamento: Cada sítio possui uma função de energia potencial por partes quadráticas. Quando a tensão elástica local atinge um limiar (definido como unidade), o sítio sofre uma transformação plástica (mudança de vale de potencial), gerando uma redistribuição de tensão em todo o sistema via o kernel de Eshelby.

3. Contribuições Principais e Resultados

A. O Papel da Tensão de Retorno de Eshelby ($\sigma_0$)

O estudo identifica que a tensão de retorno (back stress) gerada por uma transformação plástica, denotada por $\sigma_0$, é o fator determinante para o comportamento macroscópico.

• $\sigma_0$ é negativo e sua magnitude aumenta com o aumento de $K/\mu$ (ou seja, sistemas menos compressíveis geram uma tensão de retorno mais forte).
• A relação entre a tensão de retorno e os parâmetros do sistema é dada por $\sigma_0 \approx 1 - 2\alpha A_0$, onde $A_0$ é uma constante geométrica da rede.

B. Cisalhamento Cíclico e Transição de Histerese

Os autores analisam a transição entre três regimes observados em cisalhamento cíclico:

1. Regime Elástico (ER): Sem rearranjos plásticos.
2. Regime Plástico Reversível (RPR): Ciclos de histerese onde o sistema retorna à configuração microscópica original após um número inteiro de ciclos.
3. Regime Difusivo (DR): Trajetórias difusivas e irreversibilidade.

Resultados Chave:

• Limiar de Histerese ($\Gamma_0$): A amplitude crítica $\Gamma_0$ que marca a transição do regime elástico para o plástico reversível é determinada por $\Gamma_0 = (1 + \sigma_0)/2$.
• Efeito da Compressibilidade: Como $\sigma_0$ torna-se mais negativo com o aumento de $K/\mu$, o valor de $\Gamma_0$ diminui. Isso significa que sistemas menos compressíveis (maior $K/\mu$) exibem histerese e plasticidade reversível em amplitudes de deformação menores.
• Endurecimento Emergente: O modelo demonstra uma transição clara no comportamento de endurecimento:
  • Para $\Gamma < \Gamma_0$: Comportamento de endurecimento isotrópico (o intervalo elástico cresce simetricamente).
  • Para $\Gamma > \Gamma_0$: Comportamento de endurecimento cinemático (o intervalo elástico desloca-se no espaço de tensões sem mudar de tamanho).
  • Essa transição ocorre exatamente no ponto crítico $\Gamma_0$.

C. Cisalhamento Direto (Steady Shear)

No regime de cisalhamento direto, o sistema atinge um estado estacionário com uma tensão de escoamento $\sigma_y$.

• Dependência de $K/\mu$: A tensão de escoamento $\sigma_y$ diminui à medida que $K/\mu$ aumenta. Sistemas menos compressíveis fluem a tensões menores.
• Colapso de Dados: As distribuições de tensão $P(\sigma)$ em estado estacionário para diferentes $K/\mu$ colapsam em uma única curva universal quando as tensões são reescaladas usando $\sigma_0$. Isso permite derivar uma expressão analítica para $\sigma_y$ em termos de $\sigma_0$.
• Mecanismo: A maior amplitude de redistribuição de tensão em sistemas menos compressíveis (devido ao maior $\alpha$) faz com que os sítios atinjam o limiar de escoamento mais cedo, reduzindo a tensão macroscópica necessária para manter o fluxo.

4. Significado e Implicações

• Explicação da Ductilidade: O trabalho fornece uma explicação mecânica fundamental para a observação experimental de que ligas metálicas com alta razão $K/\mu$ são mais dúcteis. No modelo, a maior ductilidade (escoamento a tensões mais baixas e maior capacidade de deformação plástica reversível) surge diretamente da dependência da tensão de retorno de Eshelby com o coeficiente de Poisson, sem a necessidade de parâmetros ad-hoc de "volume livre" ou mecanismos microscópicos complexos.
• Validação de Modelos de Campo Médio: Os resultados mostram que características macroscópicas complexas, como a transição de endurecimento isotrópico para cinemático e a dependência da tensão de escoamento com a compressibilidade, podem emergir de modelos simples baseados apenas na redistribuição elástica de Eshelby. Isso valida o uso de abordagens de campo médio que focam na distribuição de tensões locais.
• Aplicabilidade: Os resultados são relevantes para uma vasta gama de materiais amorfos, incluindo ligas metálicas amorfas, vidros coloidais, emulsões e espumas, onde a razão $K/\mu$ pode variar significativamente dependendo da composição, história de resfriamento (quench) ou processamento mecânico.

Em resumo, o artigo demonstra que a compressibilidade finita é um parâmetro crucial que controla a amplitude da redistribuição de tensões após eventos de escoamento local, ditando diretamente a tensão de escoamento macroscópica, a amplitude crítica para histerese e o tipo de endurecimento observado em materiais amorfos.