A variational critical-state theory of friction

TL;DR

This paper introduces a variational, finite-kinematics model for fault gouge behavior in earthquakes, capturing complex inelastic and rate-dependent responses, and connects it to empirical friction laws.

Contribution

It develops a novel variational framework for fault gouge, incorporating rate- and state-dependent behavior with explicit solutions and experimental validation.

Findings

Model accurately predicts shear response of granular layers.

Connections established between the model and empirical friction laws.

Validated against standard laboratory shear experiments.

Abstract

Friction plays a fundamental role in many natural processes, including earthquakes, landslides, and volcanic eruptions. Earthquakes occur when highly compressed fault surfaces accumulate large enough shear stresses, causing the faults to move relative to one another, or slip. The slip is accommodated within a thin layer of comminuted granular material -- called fault gouge -- between the fault surfaces. As a result, characterizing the mechanical behavior of fault gouge in response to shear is a major open problem in earthquake source physics. Modeling gouge is complicated by large deformations, inelasticity, rate dependence, and volumetric changes. As such, researchers typically rely on empirical formulations to capture the effective response. Here, we systematically develop a variational, finite-kinematics framework for fault gouge. We first describe a general theory for a rigid-viscoplastic, pressure-sensitive material, where the plasticity evolution follows from the principle of maximum dissipation. Then, we specialize the governing equations for a Cam-Clay material within a shearing and dilating layer. We rely on convexity considerations and experimental observations from consolidation tests of granular layers to calibrate the model and develop explicit solutions for the rate- and state- dependent response of the model to shear tests under constant compressive normal stress and prescribed shearing rate. To validate the model, we select common rate functions and compare numerical material point tests and theoretical solutions to standard laboratory experiments of shearing granular layers. Lastly, we discuss connections of the model to empirical rate-and-state friction laws.