A variational critical-state theory of friction

1. Problema e Contexto

O artigo aborda a complexa mecânica das gouges de falha (camadas de material granular esmagado e compactado entre superfícies de falhas geológicas). Este material desempenha um papel fundamental em processos naturais como terremotos, deslizamentos de terra e erupções vulcânicas.

• Desafio Principal: Caracterizar o comportamento mecânico da gouge sob cisalhamento é um problema aberto na física das fontes sísmicas. O material sofre grandes deformações, exibe inelasticidade, dependência da taxa de deformação (rate-dependence) e mudanças volumétricas (dilatação ou compactação).
• Limitações Atuais: Modelos existentes frequentemente dependem de formulações empíricas (como as leis de atrito rate-and-state) ou modelos de plasticidade com superfícies de escoamento prescritas a priori, o que limita a adaptabilidade a grandes deformações e a conexão física direta com a microestrutura.
• Objetivo: Desenvolver uma estrutura teórica rigorosa, baseada em princípios variacionais e cinemática finita, para descrever a resposta de camadas granulares sob cisalhamento e dilatação, conectando a mecânica dos solos (modelo Cam-Clay) à física de falhas.

2. Metodologia

Os autores propõem uma formulação variacional de cinemática finita para materiais rígido-viscoplásticos sensíveis à pressão.

• Fundamentação Teórica:

  • Utilizam a decomposição multiplicativa do gradiente de deformação ($F = F^e F^p$).
  • Assumem elasticidade rígida (deformações elásticas desprezíveis comparadas às deformações plásticas finitas típicas de falhas), o que simplifica o tensor de deformação elástica para uma pura rotação.
  • A evolução da plasticidade é derivada do princípio da dissipação máxima, garantindo objetividade material e consistência termodinâmica.
  • O modelo é especializado para o modelo Cam-Clay (um modelo de estado crítico amplamente utilizado em geotecnia), adaptado para grandes deformações.

• Formulação Variacional:

  • Define-se uma função de dissipação que depende do tensor de taxa de deformação, da pressão de consolidação e de potenciais cinéticos.
  • As relações tensão-deformação efetivas são obtidas minimizando essa função de dissipação, resultando em uma estrutura potencial que define tanto a lei de atrito quanto a evolução das variáveis de estado internas (como o volume/consolidação).

• Aplicação a Camadas Granulares:

  • O modelo é aplicado a uma camada de gouge com espessura inicial $h_0$, submetida a cisalhamento simples e variação de volume.
  • Derivam-se soluções analíticas fechadas para o caso de cisalhamento simples sob tensão normal constante, assumindo que a taxa de cisalhamento é muito maior que a taxa de dilatação.

• Validação:

  • Realizam testes numéricos de ponto material (usando um esquema de Euler implícito) para simular ensaios de oedometer (consolidação) e ensaios de cisalhamento (taxa constante e "velocity-step").
  • Comparam os resultados numéricos e as soluções analíticas com dados experimentais de laboratório e previsões teóricas existentes.

3. Principais Contribuições

1. Estrutura Variacional Unificada: Apresentam uma formulação que deriva a lei de atrito efetiva e a evolução do estado de consolidação a partir de um único funcional de dissipação, eliminando a necessidade de postular leis constitutivas ad hoc.
2. Conexão com o Estado Crítico: Adaptam o modelo Cam-Clay (originalmente para argilas) para falhas geológicas, unificando compressibilidade, escoamento e endurecimento/amacimento em um framework de grandes deformações.
3. Soluções Analíticas Fechadas: Derivam expressões explícitas para a tensão de cisalhamento e a evolução da espessura da camada em função do deslizamento, capturando o comportamento transitório e o estado estacionário.
4. Interpretação Física das Variáveis de Estado: Identificam a variável de estado (comumente usada em leis empíricas) como o estado efetivo de consolidação da camada, fornecendo uma interpretação física clara baseada na mudança de volume.
5. Mecanismo de Transição de Estabilidade: Demonstram como a razão entre a pressão de pré-consolidação e a tensão de confinamento controla a transição entre endurecimento por taxa (rate-strengthening, deslizamento estável) e amolecimento por taxa (rate-weakening, instabilidade sísmica).

4. Resultados Chave

• Resposta de Consolidação: O modelo reproduz corretamente o endurecimento do material com o aumento da consolidação (maior pré-consolidação), alinhando-se com ensaios de oedometer.
• Comportamento em Cisalhamento (Velocity-Step):
  • O modelo captura a resposta não monotônica típica: um salto instantâneo na tensão de cisalhamento seguido por uma evolução gradual para um novo estado estacionário após uma mudança brusca na taxa de deslizamento.
  • Preveem corretamente a dilatação em taxas de carregamento mais rápidas e a compactação em taxas mais lentas.
• Transição de Comportamento Reológico:
  • Identificaram que o parâmetro chave $\hat{\sigma}_n$ (razão entre tensão de confinamento e pressão de pré-consolidação) determina o comportamento.
  • Quando a camada está altamente consolidada (alta pré-consolidação relativa), o material tende a dilatar e amolecer (rate-weakening), favorecendo a nucleação de terremotos.
  • Quando menos consolidada, o material compacta e endurece (rate-strengthening), favorecendo o deslizamento estável (creep).
• Validação Experimental: As previsões numéricas e analíticas mostram excelente concordância qualitativa e quantitativa com dados experimentais de laboratório sobre camadas granulares e com as leis de atrito rate-and-state empíricas.

5. Significado e Implicações

Este trabalho oferece uma fundação racional para as leis empíricas de atrito rate-and-state, que são amplamente utilizadas para modelar a nucleação e propagação de terremotos, mas cujos parâmetros são frequentemente inferidos empiricamente sem uma base física microestrutural clara.

• Avanço Teórico: Ao conectar a mecânica de solos (estado crítico) com a física de falhas, o modelo permite prever o comportamento de falhas sob diferentes condições de carregamento e história de consolidação, algo difícil com modelos puramente empíricos.
• Previsibilidade de Terremotos: A capacidade do modelo de capturar a transição entre deslizamento estável e instável (sismogênico) baseada em variáveis físicas mensuráveis (consolidação e confinamento) oferece novos insights sobre onde e por que os terremotos nucleiam.
• Futuro: A estrutura variacional proposta facilita a extensão do modelo para incluir efeitos de fluidos porosos (pressão de poros) e aquecimento por atrito, tornando-o uma ferramenta promissora para simulações mais realistas de dinâmica de falhas.

Em resumo, o artigo substitui a descrição fenomenológica do atrito de falhas por uma teoria constitutiva rigorosa baseada em princípios termodinâmicos e de estado crítico, melhorando a compreensão dos mecanismos físicos que governam a sismicidade.