skip to main content

Title: Meso-scale dislocations and friction of shape-complementary soft interfaces
The interface between two surfaces patterned with complementary shapes such as arrays of ridge–channel structures or pillars accommodates relative misorientation and lattice mismatch by spontaneous production of dislocation arrays. Here, we show that the relative sliding of such an interface is accomplished by dislocation glide on the interfacial plane. An exception is the singular case where the lattices are perfectly matched across the sample dimension, in which case sliding is accompanied by motion of edge-nucleated defects. These are meso-scale analogues of molecular sliding friction mechanisms between crystalline interfaces. The dislocations, in addition to the long-range elastic energy associated with their Burgers vectors, also cause significant out-of-plane dilation, which props open the interface locally. For this reason, the sliding friction is strongly pressure dependent; it also depends on the relative orientation of the patterns. Sliding friction can be strongly enhanced compared with a control, showing that shape-complementary interfaces can be engineered for strongly enhanced pressure- and orientation-dependent frictional properties in soft solids.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Journal of The Royal Society Interface
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Localized frictional sliding on faults in the continental crust transitions at depth to distributed deformation in viscous shear zones. This brittle‐ductile transition (BDT), and/or the transition from velocity‐weakening (VW) to velocity‐strengthening (VS) friction, are controlled by the lithospheric thermal structure and composition. Here, we investigate these transitions, and their effect on the depth extent of earthquakes, using 2D antiplane shear simulations of a strike‐slip fault with rate‐and‐state friction. The off‐fault material is viscoelastic, with temperature‐dependent dislocation creep. We solve the heat equation for temperature, accounting for frictional and viscous shear heating that creates a thermal anomaly relative to the ambient geotherm which reduces viscosity and facilitates viscous flow. We explore several geotherms and effective normal stress distributions (by changing pore pressure), quantifying the thermal anomaly, seismic and aseismic slip, and the transition from frictional sliding to viscous flow. The thermal anomaly can reach several hundred degrees below the seismogenic zone in models with hydrostatic pressure but is smaller for higher pressure (and these high‐pressure models are most consistent with San Andreas Fault heat flow constraints). Shear heating raises the BDT, sometimes to where it limits rupture depth rather than the frictional VW‐to‐VS transition. Our thermomechanical modeling framework can be used to evaluate lithospheric rheology and thermal models through predictions of earthquake ruptures, postseismic and interseismic crustal deformation, heat flow, and the geological structures that reflect the complex deformation beneath faults.

    more » « less
  2. The γ-surface represents the energetic cost associated with relative, rigid body sliding of crystal planes and contains useful information related to plastic deformation of the respective crystal. Here, we present γ-surfaces for the most active glide planes of the energetic molecular crystal cyclotetramethylene-tetranitramine in the monoclinic β phase, i.e., (101) and (011), at pressures up to 15 GPa. We observe the existence of stable staking faults in both planes and at all pressures and report the increase in the stacking fault energy with pressure. We also report the energetic barriers for sliding along minimum energy paths in various directions contained in these planes as well as the critical resolved shear stress at which the crystal becomes unstable in the absence of crystal defects. [100] traces of the γ-surface for multiple planes such as (001), (010), and (021) are further evaluated in view of the previously reported importance of this slip direction for dislocation cross-slip. It is observed that increasing the pressure does not modify the topology of the γ-surface in an essential way, which implies that although barriers for slip increase, the general phenomenology of dislocation motion is not modified qualitatively by the pressure. The energy barriers increase faster with pressure in the (011) plane, and hence, it is implied that the (101) plane is the most active glide plane at high pressures. The results are generally relevant for studies of plastic deformation in this molecular crystal.

    more » « less
  3. Nearly all frictional interfaces strengthen as the logarithm of time when sliding at ultra-low speeds. Observations of also logarithmic-in-time growth of interfacial contact area under such conditions have led to constitutive models that assume that this frictional strengthening results from purely time-dependent, and slip-insensitive, contact-area growth. The main laboratory support for such strengthening has traditionally been derived from increases in friction during “load-point hold” experiments, wherein a sliding interface is allowed to gradually self-relax down to subnanometric slip rates. In contrast, following step decreases in the shear loading rate, friction is widely reported to increase over a characteristic slip scale, independent of the magnitude of the slip-rate decrease—a signature of slip-dependent strengthening. To investigate this apparent contradiction, we subjected granite samples to a series of step decreases in shear rate of up to 3.5 orders of magnitude and load-point holds of up to 10,000 s, such that both protocols accessed the phenomenological regime traditionally inferred to demonstrate time-dependent frictional strengthening. When modeling the resultant data, which probe interfacial slip rates ranging from 3 . μ m · s − 1 . to less than 10 − 5 μ m · s − 1 , we found that constitutive models where low slip-rate friction evolution mimics log-time contact-area growth require parameters that differ by orders of magnitude across the different experiments. In contrast, an alternative constitutive model, in which friction evolves only with interfacial slip, fits most of the data well with nearly identical parameters. This leads to the surprising conclusion that frictional strengthening is dominantly slip-dependent, even at subnanometric slip rates. 
    more » « less
  4. Understanding the orientation dynamics of anisotropic colloidal particles is important for suspension rheology and particle self-assembly. However, even for the simplest case of dilute suspensions in shear flow, the orientation dynamics of non-spherical Brownian particles are poorly understood. Here we analytically calculate the time-dependent orientation distributions for non-spherical axisymmetric particles confined to rotate in the flow–gradient plane, in the limit of small but non-zero Brownian diffusivity. For continuous shear, despite the complicated dynamics arising from the particle rotations, we find a coordinate change that maps the orientation dynamics to a diffusion equation with a remarkably simple ratio of the enhanced rotary diffusivity to the zero shear diffusion: $D_{eff}^{r}/D_{0}^{r}=(3/8)(p-1/p)^{2}+1$ , where $p$ is the particle aspect ratio. For oscillatory shear, the enhanced diffusion becomes orientation dependent and drastically alters the long-time orientation distributions. We describe a general method for solving the time-dependent oscillatory shear distributions and finding the effective diffusion constant. As an illustration, we use this method to solve for the diffusion and distributions in the case of triangle-wave oscillatory shear and find that they depend strongly on the strain amplitude and particle aspect ratio. These results provide new insight into the time-dependent rheology of suspensions of anisotropic particles. For continuous shear, we find two distinct diffusive time scales in the rheology that scale separately with aspect ratio $p$ , as $1/D_{0}^{r}p^{4}$ and as $1/D_{0}^{r}p^{2}$ for $p\gg 1$ . For oscillatory shear flows, the intrinsic viscosity oscillates with the strain amplitude. Finally, we show the relevance of our results to real suspensions in which particles can rotate freely. Collectively, the interplay between shear-induced rotations and diffusion has rich structure and strong effects: for a particle with aspect ratio 10, the oscillatory shear intrinsic viscosity varies by a factor of ${\approx}2$ and the rotational diffusion by a factor of ${\approx}40$ . 
    more » « less
  5. Basal slip along glaciers and ice streams can be significantly modified by external time-dependent forcing, although it is not clear why some systems are more sensitive to tidal stresses. We have conducted a series of laboratory experiments to explore the effect of time varying load point velocity on ice-on-rock friction. Varying the load point velocity induces shear stress forcing, making this an analogous simulation of aspects of ice stream tidal modulation. Ambient pressure, double-direct shear experiments were conducted in a cryogenic servo-controlled biaxial deformation apparatus at temperatures between −2°C and −16°C. In addition to a background, median velocity (1 and 10 μm/s), a sinusoidal velocity was applied to the central sliding sample over a range of periods and amplitudes. Normal stress was held constant over each run (0.1, 0.5 or 1 MPa) and the shear stress was measured. Over the range of parameters studied, the full spectrum of slip behavior from creeping to slow-slip to stick-slip was observed, similar to the diversity of sliding styles observed in Antarctic and Greenland ice streams. Under conditions in which the amplitude of oscillation is equal to the median velocity, significant healing occurs as velocity approaches zero, causing a high-amplitude change in friction. The amplitude of the event increases with increasing period (i.e. hold time). At high normal stress, velocity oscillations force an otherwise stable system to behave unstably, with consistently-timed events during every cycle. Rate-state friction parameters determined from velocity steps show that the ice-rock interface is velocity strengthening. A companion paper describes a method of analyzing the oscillatory data directly. Forward modeling of a sinusoidally-driven slider block, using rate-and-state dependent friction formulation and experimentally derived parameters, successfully predicts the experimental output in all but a few cases. 
    more » « less