On the relation between plasticity, friction, and geometry

On the relation between plasticity, friction, and geometry


Plasticity refers to thermodynamically irreversible deformation associated with a change of configuration of materials. Friction is a phenomenological law that describes the forces resisting sliding between two solids or across an embedded dislocation. These two types of constitutive behaviors explain the deformation of a wide range of engineered and natural materials. Yet, they are typically described with distinct physical laws that cloud their inherent connexion. Here, I introduce a multiplicative form of kinematic friction that closely resembles the power-law flow of viscoplastic materials and that regularizes the constitutive behavior at vanishing velocity, with important implications for rupture dynamics. Using a tensor-valued state variable that describes the degree of localization, I describe a constitutive framework compatible with viscoplastic theories that captures the continuum between distributed and localized deformation and for which the frictional response emerges when the deformed region collapses from three to two dimensions.

I Introduction

Figure 1: Displacement fields in two-dimensional models of deformation due to increasingly localized deformation. We consider an embedded dislocation buried at 20 km depth, 20 km long in the down-dip direction, and dipping 45 from the vertical. The left panels are for a screw dislocation and the right ones for an edge dislocation. The top, middle and bottom panels are for shear zones of 10 km, 5 km, and 200 m, respectively. The contour lines are every 20 cm.

The deformation of solid media under stress is complex, often nonlinear, and involves a wide range of microphysical processes. Anelastic deformation represents the thermodynamically irreversible deformation of materials, implying that under this type of deformation materials do not recover their initial configuration if the applied stress is removed. Often, deformation includes both elastic (reversible) and anelastic components. Plastic flow refers to the distributed anelastic deformation of materials under stress and is observed in most engineered and natural materials, including metals, soils, and rocks. Friction describes the resistance to sliding of solid surfaces in contact. Being confined at a solid interface, friction represents a case of localized deformation.

The rheology of plastic flow has been classically described by the von Mises model[1]; [2] or its variants[3] that assumes no plastic flow when the stress is below a yield surface. Friction has also be described mathematically as a yield surface, implying that sliding may not occur unless the shear forces reach a critical value that depends on normal stress. These results were found experimentally by Da Vinci, re-discovered by Amontons[4], and confirmed by Coulomb[5]. These models afford a simple description of the phenomenology but they do not describe how fast the medium flows or breaks when placed under sufficient stress nor the complex dynamics or instabilities that can develop near the yield surface[6].

The limitations of the classic formulations have been remedied by the development of viscoplasticity[7]; [8] and the Dieterich-Ruina rate-and-state form of kinematic friction[9]; [10]. In the classic view of plastic flow and friction, anelastic deformation occurs only at the yield surface and the stress cannot exceed the yield surface (Fig. 1). In the modern view, the yield surface represents a reference point under which deformation is very slow and above which deformation is very fast. The concept of a yield surface is only useful if the anelastic deformation rate is superlinear with stress. The weak sensitivity of the frictional resistance to sliding velocity noted by Coulomb is captured by the logarithmic sensitivity to velocity in the modern formulation of rate-and-state friction (a wide range of slip velocity occurs for a narrow range of effective stress). For plastic flow, the strength is often rate dependent, taking a linear form for viscoelastic materials and often a power-law form for plastic materials. However, there is a continuum between linear, power-law, and exponential stress-strain-rate relationships embodied by processes such as diffusion creep, dislocation creep[11], and Peierls creep[12]; [13]. Like friction, plastic flow may also necessitate state variables for phenomena associated with transient creep[14] and work-hardening in general[15]. The viscoplasticity and rate-and-state friction theories have allowed modeling multiple cycles of ruptures with complex deformation histories[16]; [17].

As noted by others[18]; [19], the similitudes in the parallel development of the theories of plasticity and friction are striking. The mechanics of friction and plasticity are also deeply related. For example the micromechanics of plasticity can be explained as the motion of internal dislocations[20]; [21] and the rate- and state-dependence of friction can be explained by the plastic behavior of the contact asperities regulating the real area of contact[22]; [23]; [24]. In many cases, the frictional interface consists in a finite-width shear zone made of granular, damaged, or semi-brittle material and the assumption of a simple frictional behavior is for mathematical or modeling convenience only[25]; [26]. This suggests that these theories could be combined in a more general, unifying framework.

Plasticity and friction represent the end-members for distributed and localized deformation, respectively. However, they cannot capture the continuum between localized and distributed deformation. This can be particularly important when describing the inception of new frictional surfaces by the accumulation of damage or increased localization. During this process, the material first deforms plastically until a frictional interface has matured. When this point is reached, a new surface has emerged and friction can become the dominant deformation mechanism.

In geophysics, there is an impetus to combine long-term and short-term processes to allow faults to emerge from tectonic processes and resolve rupture dynamics on these new structures[27] as important feedbacks operate at these different time scales. That dynamic ruptures on frictional interfaces and plastic deformation operate at different time scales represents a significant numerical challenge. But there is no convenient theoretical framework to describe friction and plasticity in a unified way.

Towards that objective, I highlight two important results. In Section II, I show that the kinematics of distributed deformation and localized deformation can be captured by a single mathematical expression[28]. This highlights a continuum of kinematic descriptions between distributed and localized deformation. In Section III, I show that the modern rheological laws that describe plastic flow have a striking similarity with the ones describing frictional dynamics, despite some subtle but important differences. To do so, I introduce the multiplicative form of rate-and-state friction. In Section IV, I conjecture a constitutive framework that reduces to viscoplasticity in some situations and to kinematic friction framework in others. The model predicts that frictional behavior is an end-member rheological response due to the collapse of deformation region from three to two dimensions.

In this paper, I do not discuss a new microphysical model of plasticity or friction. Such important models shed light on the constitutive equations of anelasticity and can be found in the work of various authors[29]; [30]; [31]; [32]; [33]; [34]. Instead, I discuss the relation between the macroscopic representations of plasticity and friction in a continuum that are directly incorporated in analytic or numerical models of time-dependent deformation.

Ii Kinematic continuum between distributed and localized deformation

In this Section, I consider the kinematics of deformation of plastic and brittle materials and how anelastic deformation induces strain in the surrounding medium by elastic coupling. The static or quasi-static equilibrium attained in the presence of anelastic deformation can be modeled using transformation strain[15]; [35], which can be due to thermal stress, multi-phase flow, or plastic deformation[26]. Steketee[25]; [36] showed that the Schmidt tensor, a particular form of transformation strain, can capture the motion of dislocations. However, because of the fundamental difference in nature between distributed and localized deformation, these cases were historically treated distinctly.

Figure 2: Displacement fields due to increasingly localized deformation. We consider the plastic deformation of a cuboid representative volume element (rectangle) buried at 20 km depth, 20 km long in the down-dip direction, and dipping 45 from the vertical. The left panels are for out-of-plane plastic strain and the right ones for in-plane plastic strain. The top, middle and bottom panels are for shear zones of 10 km, 5 km, and 200 m, respectively. The contour lines are every 20 cm. With increasing localization of plastic deformation, the deformation converges to the cases of a screw dislocation and an edge dislocation for the left and right panels, respectively.

Here, I show how various degrees of localization can be described within a single framework. To do so, I consider the plastic deformation of a shear zone confined in a larger continuum and subjected to the transformation strain components . The region that deforms anelastically takes the form of a cuboid with dimension , , and . For convenience, let us imagine a reference system aligned with the cuboid so that the subscript indices match those of the basis vectors. The resulting deformation can be evaluated by solving the governing equation for elasticity using equivalent body forces to incorporate the plastic strain[26]; [28]; [37]. For illustration purposes, I consider a half space with a free surface and the cases of in-plane and anti-plane deformation.

I subject the shear zone to increasing plastic strain in a cuboid of decreasing width . For the sake of simplicity, I adopt a two-dimensional model where and I choose a grid size of 250 m for all calculations, which also represents the dimension of a representative volume element. I reduce the plastic strain to the nontrivial component for the case of in-plane strain and to the nontrivial component for the case of anti-plane strain. The cumulative displacement across the shear zone due to plasticity alone, i.e., irreversible, is


for in-plane strain and


for anti-plane strain. In these simulations, I keep and constant and I reduce incrementally. The resulting deformations are shown in Fig. 2. When the width of the shear zone is lower than the grid size, the shear zone virtually collapses from three to two dimensions, i.e., from a volume to a surface. The solutions converge towards the one for an edge dislocation for in-plane strain and towards the one for a screw dislocation for anti-plane strain. Using in-plane strain as an example, the limit case of localized deformation is obtained for


That is, the motion of dislocation is associated with infinite strain. For brittle deformation, plastic strain is infinite, the support is zero, but the product is finite and the cumulative plastic displacement is called slip[18].

These results indicate that brittle deformation can be thought of as an end-member of plasticity for infinite strain over a zero volume, at least in terms of the kinematics of deformation. These results have been exploited to simulate fault slip in the framework of transformation strain using numerical methods[25]; [26] and to derive analytic expressions for the deformation of elasto-plastic media for various degrees of localization, from distributed to localized[28]; [37]. In reality, even the dislocations associated with the motion of lattice defects in crystal structures have a finite width commensurate to the atomic scale (Å), so plastic strain is always bounded.

Unified kinematic models of anelastic deformation from distributed to localized have been explored[26]; [28] and form a motivation to explore the possibility of a unification of the constitutive laws.

Iii The multiplicative form of kinematic friction and its relation to viscoplasticity

Figure 3: Representative elements for viscoplasticity and kinematic friction. A) The dynamics of viscoplastic flow is represented by a stress-strain-rate constitutive relationship, the simplest one being a spring-dashpot model. B) The dynamics of frictional sliding relies on a velocity-stress constitutive relationship that incorporates the orientation of the representative surface element. The simplest model is a spring-slider.

Viscoplasticity describes the rate-dependent deformation of plastic materials under stress[38]. The model describes the response of a representative volume element (Fig. 3a). The simplest model for viscoplasticity is a combination of a spring and a (potentially nonlinear) dashpot. Broadly speaking, the rate of plastic strain can be described as


where is the plastic strain rate tensor, is the stress tensor, and the are state variables following the evolution laws


The inclusion of state variables and the cumulative strain in the constitutive behavior allows for the transient and work-hardening effects and makes the rate of deformation non-conservative, i.e., dissipative and path dependent.

Modern friction describes the rate of sliding across a representative surface element upon application of shear and normal tractions[9]; [10]. The simplest constitutive model is to combine a spring and a slider (Fig. 3b). Rate-and-state friction is a phenomenological law that describes the frictional resistance to sliding including healing and weakening. This framework has some known shortcomings[39]; [40]; [41], but it has nonetheless been successful at explaining many dynamic mechanical systems over a wide range of materials and boundary conditions[42]; [43]; [44]; [45]; [46]; [47]; [48]; [49]; [50]; [51]. A friction law can be written as follows


where v is the sliding velocity, the are state variables with the evolution laws


and n is the local normal direction of the surface where frictional sliding occurs. The notable difference between plasticity and friction is the role of geometry. The constitutive equation for viscoplasticity applies everywhere in the medium. In constrast, the traction and the slip vector cannot be established without a proper representation of the surface in terms of its orientation and position. Indeed, frictional sliding may be associated with the transformation strain


that explicitly incorporates the normal vector of the frictional surface. Any constitutive law for anelastic deformation that does not explicitly include geometry cannot fully represent friction.

iii.1 The multiplicative and additive forms of kinematic friction

Figure 4: A finite shear zone model for frictional sliding. The relative velocity across frictional interface is accommodated by plastic strain within the shear zone.

With this important difference in mind, I now describe some similarities in form between the constitutive laws for viscoplasticity and kinematic friction. I consider a thought experiment where the frictional interface is modeled by a narrow shear zone of width undergoing plastic deformation (Fig. 4) with the following rheology


where is a reference strain rate, is the norm of the shear stress in the direction of sliding, is the average grain size in a gouge layer, is a reference value, and is a reference stress. I assume that grain size is a state variable that varies with


where the first term in the right-hand side gives rise to healing and the second term accounts for dynamic recrystallization[52]. Classically, viscoplastic rheologies are formulated with a different nomenclature, but I make this choice of parameterization for the sake of the argument. The power exponent is assumed in the range , corresponding to low-temperature creep. The grain size exponent is assumed in the range , the bounds corresponding to negligible and significant dependence on grain size, respectively.

The above formulation can be shown to be mathematically equivalent to rate-and-state friction with the aging law with the following few assumptions. First, to ease the comparison, I turn the strain rate into a velocity using and . Using the state variable , with , I obtain


Then, I recast the stress-strain-rate relationship in terms of frictional resistance


As the sensitivity of stress to velocity is only weak, I consider a Taylor series expansion using around , expanding the relationship (12) accordingly around . Defining friction as and keeping only the linear terms, I obtain


The correspondence with the classic nomenclature of rate-and-state friction is immediate if we note , , and . To go further let us consider the evolution law (10). Turning the strain rate into a velocity, I obtain


or, simply,


which is the aging law of rate-and-state friction[10]. In this formulation the state variable is interpreted as the average grain size scaled by a reference sliding velocity and the reference grain size coincides with the width of the gouge layer, which also matches with the characteristic weakening distance. The direct effect coefficient is related to the reciprocal of the power exponent, consistent with values of of the order of . The parameter is controlled by the dependence to grain size.

I note that there is no laboratory experiment that describes a combination of power-law flow with a high stress exponent with a grain-size dependence such as in (9), although it complies to the case of dislocation glide[53] with . Nonetheless, the plastic model for friction encapsulated in equations (9) and (10) is validated a posteriori by a mathematical equivalence with rate-and-state friction, an empirical law tuned to a number of laboratory experiments. Admittedly, the proof is only mathematical, and therefore the evolution law (10) may not represent grain size but any other physical property that varies with the same rate.

I call (11) and (12) the multiplicative or power-law form of rate-and-state friction, and the Dieterich-Ruina formulation (13) the additive form of rate-and-state friction. The multiplicative form is more appealing than the additive one because the latter is ill-posed for vanishing velocities. Indeed, (13) incorrectly predicts a negative norm of the traction vector for sufficiently small velocities and an infinite value for a vanishing velocity. The viscoplastic model of frictional resistance also brings new insight. Indeed, some laboratory measurements support that the critical slip distance correlates with the width of the frictional interface[54]; [55].

I emphasize that (11) is not formally a viscoplastic law because the dynamic variables include the shear and normal traction instead of the deviatoric stress and the pressure. That is, the geometric layout of the shear zone has been incorporated in (11) to conform to the general case (6). However, the similarity between (11) and the power-law form of viscoplasticity is striking and hints at the existence of a constitutive relationship that incorporates plasticity and friction as end-members.

iii.2 Comparison with laboratory data

Figure 5: Healing and the additive and multiplicative forms of rate-and-state friction. The static coefficient of friction data for quartz sandstone from[56] for hold tests. The hold time corresponds to the duration of stick. The data are reduced with the additive model (13) (solid lines) and the multiplicative model (12) (dashed profiles). Both models explain the data well within their scatter.

The choice of a constitutive framework may have important implications on the predictions of failure modes, from laboratory samples to large-scale geodynamics. Therefore, one should proceed carefully when introducing a different friction law. The additive form (13) of rate-and-state friction was derived empirically by Dieterich[9] and Ruina[10] based partially on hold experiments that allowed re-strengthening at zero slip rate[56]. Here, I revisit the data collected by Dieterich[56] to assess the merit of using (11) to explain these observations. In Fig. 5, I show the static coefficient of friction of five laboratory experiments conducted at different normal stresses for a wide range of hold times. While different, the predictions of the additive (solid profiles) and multiplicative (dashed profiles) forms of rate-and-state friction are both well within the data scatter and have virtually equal merit in reducing the laboratory data.

The apparent linear scaling of the static coefficient of friction with the logarithm of hold time () motivated Dieterich for his now classic formulation. The multiplicative form of rate-and-state friction (11) implies a different scaling, namely a linear dependence between the logarithm of the coefficient of friction with the logarithm of hold time (). Some broader considerations justify the latter as a better choice.

Tarantola[57] showed that the natural distance between two physical quantities and that cannot be negative, which he refers to as Jeffreys quantities, is . This metric is insensitive to the choice of physical unit, such that is the same whether and are expressed in Pa or MPa. As Jeffreys quantities live in the interval , their associated random variable cannot be normally distributed and their non-informative prior cannot be uniformly distributed, both precluding using simple least squares for parameter estimation. The yield stress and the amplitude of the traction vector, which are equated through the consistency relation , are Jeffrey parameters. Therefore, it is preferable to consider the metric rather than to describe the difference between any two friction coefficients. The same argument applies for the velocity dependence , as velocity and its conjugate, slowness, are also Jeffrey’s parameters. Consequently, a power-law form for the yield stress dependence on the velocity and the state variable is better suited than a logarithmic form.

iii.3 The grain-size evolution of Hall and Parmentier

Several evolution laws have been proposed in the field of viscoplasticity that may inspire new developments in kinematic friction. Hall and Parmentier[52] introduce an evolution law for grain size where the thermally activated grain growth reduces the interfacial energy associated with grain boundaries and grain size reduction occurs by sub-grain rotation. They write


where is temperature dependent and controls the healing rate, is a non-dimensional parameter controlling the rate of grain size reduction, and is a non-dimensional material parameter in the range 1 to 4. To reduce the evolution law (16) to a reference form, I conduct the following change of variable. First, using the geometric setting of Fig. 4, I approximate the strain rate by the ratio of the cumulative plastic velocity across the shear zone scaled by the width


where I used a similar scaling for the reference strain rate. Then, I introduce a new state variable that scales linearly with grain size


Finally, I introduce a critical distance


With these changes, the coupled constitutive equations (9) and (16) become


This formulation closely resembles the power-law form of rate-and-state friction. Indeed, the grain-size evolution law (16) takes the form of the aging law (15) of rate-and-state friction for . However, if the steady-state grain size obeys the piezometric relationship[58], then grain-size evolution only leads to velocity-strengthening behavior (e.g., ).

Because the mechanisms of grain-size evolution during plastic flow and the evolution of effective contact area during frictional sliding may be different, the similarities between these evolution laws indicate that multiple mechanisms of material reconfiguration at the microscopic level may take a similar form mathematically and sometimes produce the same macroscopic constitutive behavior.

iii.4 The grain-size evolution of Rozel, Ricard and Bercovici

Rozel, Ricard and Bercovici[59] employ basic non-equilibrium thermodynamics to propose a general equation for the mean grain size evolution under the assumption that the whole grain size distribution remains self-similar. They obtain the evolution law of the form


where the healing term is similar to the formulation of Hall and Parmentier[52], but the rate of grain size reduction is proportional to the shear stress. Assuming the geometrical setting of Fig. 4 and the constitutive equation (9), we have (recall the weak dependence of frictional strength on sliding velocity). Therefore, the evolution law (21) can be simplified to


in the context of a shear zone, where for the law of Rozel, Ricard and Bercovici[59] and for the law of Hall and Parmentier[52]. The weakening term is now similar to the one in (16), except for a factor . To reduce the evolution law (22) to a reference form, I use the state variable (18) and introduce the scaling parameter


With these changes, the coupled constitutive equations (12) and (22) become


Because of the choice of parameterization, eq. (24) reduces to the evolution law of Hall and Parmentier[52] and the aging law of rate-and-state friction for specific values of and . For example, the formulation simplifies to the aging law of rate-and-state friction for and .

These two examples demonstrate that largely different assumptions about the microphysics of the state evolution may lead to similar mathematical expressions for the evolution law. This may explain why the rate-and-state framework has been so successful at capturing deformation in the laboratory for so many different types of materials and stress regimes.

iii.5 Implications for rupture dynamics

Figure 6: Earthquake cycle simulations with the additive (13) and multiplicative (12) forms of rate-and-state friction. A) Simulations with the additive form, which is insensitive to the value of the static friction coefficient. B, C, D) Simulations with the multiplicative form of rate-and-state friction with a static friction coefficient of , 0.2, and 0.1, respectively. Rupture dynamics with the multiplicative form is sensitive to the strength of the fault, which among other characteristics affects the nucleation, the source-time function, and the static and dynamic stress drops of each event.

The dynamics of frictional instabilities depends on the details of the constitutive framework. To illustrate the importance of understanding the frictional resistance near zero slip speed, I explore the predictions of the additive and multiplicative forms of rate-and-state friction with numerical models of mode III ruptures in anti-plane strain. Particularly, I show that the rupture modes (the peak velocity, stress drop, stability) strongly depend on the static coefficient of friction when the multiplicative form of rate-and-state friction is assumed.

I consider a vertical frictional interface embedded in an elastic layer with a free surface. Instead of solving for the inertial contribution directly, I use radiation damping as an approximation within the boundary integral method, as in previous work[60]; [16]; [61]. I consider a 35 km vertical fault that is velocity weakening between 5 and 15 km depth with and velocity strengthening elsewhere with , all under a uniform normal stress of MPa. The system is loaded at a uniform rate of m/s for s and the initial condition is chosen for steady state at that slip velocity.

The results shown in Fig. 6 compare dynamic rupture simulations assuming the additive (top panel) and the multiplicative (three bottom panels) forms of rate-and-state friction. For a strong fault with , the results from both forms are similar in terms of recurrence times, overall number of instabilities, peak velocities and stress drops. As the static coefficient of friction decreases, the peak velocity and stress drop decrease and the details of the rupture evolution change significantly, except for the first event in the simulation that is strongly influenced by the initial condition. The importance of the static coefficient of friction comes about from the relative amplitude of the yield stress and the stress drop during failure and this effect should be exacerbated by seismic radiations.

The additive form of rate-and-state friction has been interpreted in terms of statistical physics[62]; [43] as a probability for forward jumps of the microasperities forming the real area of contact proportional to . When backward jumps are taken into account, the probability of forward jumps becomes proportional to , leading to a term in instead of , as in[62]; [43]


The multiplicative form (12) assumes the probability of forward jumps of the form , i.e., a type of power-law distribution. Further laboratory experiments may be needed to constrain the constitutive behavior at small slip velocity.

Iv A unifying constitutive framework for plastic flow and friction

Figure 7: A constitutive framework for plastic flow that incorporates the role of geometry. A) The spatial velocity gradient L describes how the plastic velocity v changes spatially. The velocity gradient is controlled by the stress tensor and the geometry of internal deformation, here modeled as a cuboid of dimension , , and with face areas , , and . B) The constitutive laws of viscoplasticity, corresponding to a stress-strain-rate relationship, are obtained for . C) The constitutive laws for kinematic friction is obtained when the deformed region collapses from a volume to a surface, here with , leading to a velocity-traction relationship.

The unified kinematic description of distributed to localized deformation described in Section II and the similarities between the power-law viscoplastic laws and rate-and-state friction discussed in Section III hint at the existence of a constitutive framework that would simplify to plasticity in some situations and to friction in others.

In this section, I describe a unifying constitutive framework where the continuum between plastic flow and friction is only controlled by geometrical factors. When friction is interpreted as localized plastic flow, the constitutive relationship is anisotropic. Therefore, I introduce a heuristic that controls plastic anisotropy based on the geometry of internal deformation within a representative volume element. As internal tractions are modulated by the area on which they are applied, I consider the ratio of internal surface areas

A (26)

as a dimensionless second-order tensor state variable of the representative volume element. In order to keep aspect ratios free variables, I consider a cuboid element of dimension , , and with face areas , , and (Fig. 7a). I define the weighted moment density tensor

M (27)

Each element of M is proportional to a component of the force acting in opposite directions on faces and . I then define the weighted mean moment density


and the weighted mean moment density tensor as


Most naturally, the definition of the deviatoric weighted moment density tensor follows as


with the apparent shear stress


Both plastic flow and kinematic friction can be formulated with a single constitutive law at the spatial scale of the representative volume element following


where is the spatial anelastic velocity gradient, v is the anelastic velocity field, and the anelastic strain rate is given by


The change of scale to the dimension of internal deformation requires the condition number of the state variable tensor


that quantifies the degree of localization within the representative volume element.

As an example, using the framework of equations (26-34), I consider two end-member cases. First, I assume that the cuboid is a cube (Fig. 7b), such as and . The state variable becomes and the weighted moment density tensor . The weighted mean moment density reduces to and the weighted mean moment density tensor simplifies to


i.e., the isotropic pressure tensor. Finally, the deviatoric moment tensor reduces to


i.e., the deviatoric stress tensor, which plays a central role in the theory of viscoplasticity. With , the proposed constitutive law (32) reduces to the typical viscoplastic rheology


as the symmetry of the deviatoric stress implies in this case.

Now I consider a second case where the cuboid flattens to a surface (Fig. 7c) with , leading to and . In practice, may be finite, but I consider the limit case. The state variable simplifies to


and the weighted moment density tensor follows as


i.e., only the components corresponding to the traction vector are non zero and the tensor is not symmetric. The weighted mean moment density reduces to , where is the normal stress to the surface . The weighted mean moment density tensor becomes


where is the normal stress in this case. Finally, the deviatoric moment density tensor reads


The norm of the deviatoric moment density tensor, , corresponds to the norm of the shear tractions that are parallel to the surface , as is required to describe the local frictional resistance. The deviatoric moment density tensor is not symmetric either, allowing for simple shear. For the case where the cuboid region collapses to a surface with , only two nontrivial components of the constitutive relationship remain


that correspond to the strike-slip and dip-slip component of the velocity gradient. Here, is the rake angle, with and . The spatial velocity gradient is not symmetric, creating simple shear, compatible with fault slip. Both and are finite because the internal deformation is averaged over the dimension of the representative volume element. The internal deformation is given by , which approaches infinity with increasing localization. As described in Section III, the formulation (42) corresponds to the multiplicative form of rate-and-state friction in three dimensions.

In the proposed constitutive framework, whether the deformation simplifies to plasticity (distributed) or friction (localized) depends on a pair of geometrical factors in three dimensions, either the pair and , the pair and , or the pair and . In two dimensions, if is infinite, the behavior will be controlled by the aspect ratio .

The comparison of these two end-members, one corresponding to the distributed macroscopic stress and the other to the tractions on a surface, shows that friction can emerge from anisotropic plasticity in the special case of large aspect ratios in the internal geometry of the representative element.

V Conclusions

The kinematics of plastic deformation in a continuum can incorporate all degrees of localization, from distributed to localized with a single mathematical expression. The constitutive law for kinematic friction closely resembles the power laws used in viscoplastic theories, in particular when the multiplicative form of rate-and-state friction is employed.

A unifying constitutive framework can be identified, in which viscoplasticity and kinematic friction represent end-members of the rheological behavior for distributed and localized deformation, respectively. The behavior is controlled by the aspect ratio of the deformation region. It is also likely that different microphysical processes become activated depending of the aspect ratio.

While I hope that this discussion will provide novel insight into the place of friction in the general framework of anelastic deformation, I have left out some important issues such as the physical origin of strain localization, in particular, what controls the thickness of a shear zone or fault gouge, and the true nature of the state variable used in (11) and (12), which is left to future work.


This work comprises Earth Observatory of Singapore contribution no. 174. This research is supported by the National Research Foundation of Singapore under the NRF Fellowship scheme (National Research Fellow Awards No. NRF-NRFF2013-04) and by the Earth Observatory of Singapore, the National Research Foundation, and the Singapore Ministry of Education under the Research Centres of Excellence initiative.



  1. MT Huber. The specific shear strain work as criterion of material strength. Czasopismo Techniczne, Lwów, 1904.
  2. R v Mises. Mechanik der festen körper im plastisch-deformablen zustand. Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, Mathematisch-Physikalische Klasse, 1913(4):582–592, 1913.
  3. D. C. Drucker and W. Prager. Soil mechanics and plastic analysis or limit design. Quarterly of applied mathematics, 10(2):157–165, 1952.
  4. G. Amontons. De la resistance causée dans les machines, tant par les frottemens des parties qui les composent, que par roideur des cordes qu’on y employe, & la maniere de calculer l’un & l’autre. Mémoires de l’Académie royale des science, 1699.
  5. CA Coulomb. The theory of simple machines. Mem. Math. Phys. Acad. Sci, 10(161-331):4, 1785.
  6. C. H. Scholz. The mechanics of earthquakes and faulting. 439 pp., Cambidge Univ. Press, New York, NY, 1990.
  7. EN da C Andrade. On the viscous flow in metals, and allied phenomena. Proc. Royal Soc. London (A), pages 1–12, 1910.
  8. F. H. Norton. The creep of steel at high temperatures. Number 35. McGraw-Hill Book Company, Incorporated, 1929.
  9. J. H. Dieterich. Time-dependent friction and the mechanics of stick-slip. Pure Appl. Geophys., 116(4-5):790–806, 1978.
  10. A. Ruina. Slip instability and state variable friction laws. J. Geophys. Res., 88:10,359–10,370, 1983.
  11. S.-I. Karato and H. Jung. Effects of pressure on high-temperature dislocation creep in olivine. Philosophical Magazine, 83(3):401–414, 2003.
  12. R Peierls. The size of a dislocation. Proc. Phys. Soc., 52(1):34, 1940.
  13. FRN Nabarro. Dislocations in a simple cubic lattice. Proc. Phys. Soc., 59(2):256, 1947.
  14. S. Masuti, S. Barbot, S. Karato, L. Feng, and P. Banerjee. Upper mantle water stratification inferred from the 2012 Mw 8.6 Indian Ocean earthquake. Nature, 538:373–377, 2016.
  15. S. Nemat-Nasser. Plasticity. A Treatise on Finite Deformation of Heterogeneous Inelastic Materials. Cambridge Monographs on Mechanics. Cambridge University Press, 2004.
  16. V. Lambert and S. Barbot. Contribution of viscoelastic flow in earthquake cycles within the lithosphere-asthenosphere system. Geophys. Res. Lett., 43(19):142–154, 2016.
  17. Kali L Allison and Eric M Dunham. Earthquake cycle simulations with rate-and-state friction and power-law viscoelasticity. Tectonophysics, 2017.
  18. R. Michalowski and Mroz Z. Associated and non-associated sliding rules in contact friction problems. Archiv. Mech., 30:259–276, 1978.
  19. Curnier A. A theory of friction. Int. J. Solids Struct., 20(7):637–647, 1984.
  20. V. Volterra. Sur l’équilibre des corps élastiques multiplement connexes. Annales scientifiques de l’École Normale Supérieure, 24:401–517, 1907.
  21. E. Orowan. Zur kristallplastizität. iii. Zeitschrift für Physik A Hadrons and Nuclei, 89(9):634–659, 1934.
  22. F. P. Bowden and D. Tabor. The friction and lubrication of Solids, Part I. Clarendon Press, Oxford, 1950.
  23. E. Rabinowicz. The intrinsic variables affecting the stick-slip process. Proc. Phys. Soc. London, 71(4):668–675, 1958.
  24. F. P. Bowden and D. Tabor. The friction and lubrication of Solids, Part II. Clarendon Press, Oxford, 1964.
  25. J. A. Steketee. On Volterra’s Dislocations in a Semi-Infinite Elastic Medium. Can. J. Phys., 36:192–205, 1958.
  26. S. Barbot and Y. Fialko. Fourier-domain Green’s function for an elastic semi-infinite solid under gravity, with applications to earthquake and volcano deformation. Geophys. J. Int., 182(2):568–582, 2010.
  27. R. Herrendörfer, Y. Van Dinther, T. Gerya, and L. A. Dalguer. Earthquake supercycle in subduction zones controlled by the width of the seismogenic zone. Nature Geoscience, 8(6):471–474, 2015.
  28. S. Barbot, J. DP Moore, and V. Lambert. Displacement and stress associated with distributed anelastic deformation in a half-space. Bull. Seism. Soc. Am., 107(2):821–855, 2017.
  29. J. H. Dieterich and B. D. Kilgore. Direct observation of frictional contacts: New insights for sliding memory effects. Pure Appl. Geophys., 143:283–302, 1994.
  30. Norman H. Sleep. Physical basis of evolution laws for rate and state friction. Geochemistry, Geophysics, Geosystems, 6(11), 2005.
  31. AR Niemeijer and CJ Spiers. A microphysical model for strong velocity weakening in phyllosilicate-bearing fault gouges. J. Geophys. Res., 112(B10), 2007.
  32. T. Putelat, J. HP Dawes, and J. R Willis. On the microphysical foundations of rate-and-state friction. J. Mech. and Phys. Solids, 59(5):1062–1075, 2011.
  33. H. Perfettini and A. Molinari. A micromechanical model of rate and state friction: 1. static and dynamic sliding. J. Geophys. Res., 122(4):2590–2637, 2017.
  34. A. Molinari and H. Perfettini. A micromechanical model of rate and state friction: 2. effect of shear and normal stress changes. J. Geophys. Res., 122(4):2638–2652, 2017.
  35. S. Nemat-Nasser and M. Hori. Micromechanics: overall properties of heterogeneous materials. Elsevier, 2 edition, 1999.
  36. J. A. Steketee. Some geophysical applications of the elasticity theory of dislocations. Can. J. Phys., 36:1168–1198, 1958.
  37. S. Barbot. Displacement and stress induced by anelastic deformation confined in a tetrahedral volume. Bull. Seism. Soc. Am., submitted, 2018.
  38. J. L. Chaboche. A review of some plasticity and viscoplasticity constitutive theories. Int. J. Plasticity, 24:1642–1693, 2008.
  39. M. Nakatani. Conceptual and physical clarification of rate and state friction: Frictional sliding as a thermally activated rheology. J. Geophys. Res., 106(B7):13347–13380, 2001.
  40. E. Bayart, A. M. Rubin, and C. Marone. Evolution of fault friction following large velocity jumps. Eos Trans. AGU, 87(52), 2006.
  41. P. Segall. Earthquake and volcano deformation. Princeton University Press, Princeton, NJ, 2010.
  42. S. T. Tse and J. R. Rice. Crustal earthquake instability in relation to the depth variation of frictional slip properties. J. Geophys. Res., 91(B9):9452–9472, 1986.
  43. N. Lapusta, J. R. Rice, Y. BenZion, and G. Zheng. Elastodynamics analysis for slow tectonic loading with spontaneous rupture episodes on faults with rate- and state-dependent friction. J. Geophys. Res., 105(B10):23765–23789, October 2000.
  44. JM Carlson and AA Batista. Constitutive relation for the friction between lubricated surfaces. Phys. Rev. E, 53(4):4153, 1996.
  45. A. Lemaître. Rearrangements and dilatancy for sheared dense materials. Phys. Rev. Lett., 89(19):195503, 2002.
  46. IS Aranson, LS Tsimring, and VM Vinokur. Stick-slip friction and nucleation dynamics of ultrathin liquid films. Phys. Rev. B, 65(12):125402, 2002.
  47. M. Urbakh, J. Klafter, D. Gourdon, and J. Israelachvili. The nonlinear nature of friction. Nature, 430(6999):525–528, 2004.
  48. N. Lapusta and S. Barbot. Models of earthquakes and aseismic slip based on laboratory-derived rate and state friction laws. In A. Bizzarri and H. S. Bhat, editors, The mechanics of Faulting: From Laboratory to Real Earthquakes, pages 153–207. Research Signpost, Trivandrum, Kerala, India, 2012.
  49. S. Barbot, N. Lapusta, and J. P. Avouac. Under the hood of the earthquake machine: Towards predictive modeling of the seismic cycle. Science, 336(6082):707–710, 2012.
  50. D. M. Veedu and S. Barbot. The Parkfield tremors reveal slow and fast ruptures on the same asperity. Nature, 532(7599):361–365, 2016.
  51. Y. Kaneko, S. B Nielsen, and B. M Carpenter. The onset of laboratory earthquakes explained by nucleating rupture on a rate-and-state fault. J. Geophys. Res., 121(8):6071–6091, 2016.
  52. C. E Hall and EM Parmentier. Influence of grain size evolution on convective instability. Geochemistry, Geophysics, Geosystems, 4(3), 2003.
  53. J. M Warren and G. Hirth. Grain size sensitive deformation mechanisms in naturally deformed peridotites. Earth Planet. Sci. Lett., 248(1):438–450, 2006.
  54. C. Marone and B. Kilgore. Scaling of the critical slip distance for seismic faulting with shear strain in fault zones. Nature, 362:618–620, 1993.
  55. C. J. Marone. Laboratory-derived friction laws and their application to seismic faulting. An. Rev. Earth Planet. Sc., 26:643–696, 1998.
  56. J. H. Dieterich. Time-dependent friction in rocks. J. Geophys. Res., 77:3690–3697, 1972.
  57. Albert Tarantola. Elements for physics: quantities, qualities, and intrinsic theories. Springer Science & Business Media, 2006.
  58. D. Van der Wal, P. Chopra, M. Drury, and J. F. Gerald. Relationships between dynamically recrystallized grain size and deformation conditions in experimentally deformed olivine rocks. Geophys. Res. Lett., 20(14):1479–1482, 1993.
  59. A. Rozel, Y. Ricard, and D. Bercovici. A thermodynamically self-consistent damage equation for grain size evolution during dynamic recrystallization. Geophys. J. Int., 184(2):719–728, 2011.
  60. Q. Qiu, E. M Hill, S. Barbot, J. Hubbard, W. Feng, E. O Lindsey, L. Feng, K. Dai, S. V Samsonov, and P. Tapponnier. The mechanism of partial rupture of a locked megathrust: The role of fault morphology. Geology, 44(10):875–878, 2016.
  61. A. Goswami and S. Barbot. Slow-slip events in semi-brittle serpentinite fault zones. Scientific Reports, 2018, in press.
  62. J. R. Rice and Y. Ben-Zion. Slip complexity in earthquake fault models. Proc. Nat. Natl. Acad. Sci., 93(9):3811–3818, 1996.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description