Abstract
Sequences of dynamic instabilities are analyzed for a single degree of freedom elastic system which slides along a surface having frictional resistance depending on slip rate and slip rate history, in the manner of Dieterich, Ruina and others. The system is represented as a rigid block in contact with a fixed surface and having a spring attached to it whose opposite end is forced to move at a uniform slow speed. The resulting “stick‐slip” motions are well understood in the classical case for which there is an abrupt drop from “static” to “sliding” frictional resistance. We analyze them here on the basis of more accurate frictional constitutive models. The problem has two time scales, an inertial scale set by the natural oscillation periodTof the analogous frictionless system asT/2π and a state relaxation scaleL/Voccurring in evolution, over a characteristic slip distanceL, of frictional stress τ towards a “steady state” value τss(V) associated with slip speedV. We show that τ ≃ τss(V) during motions for which acceleration a satisfiesaL/V2≪ 1, and that this condition is met during an inertia controlled instability in typical circumstances for which the unstable slip is much greater thanL. SinceV/ais of orderT/2π during inertia controlled motion, one hasL/VT/2π, whereasL/VT/2π during much of the essentially quasi‐static “stick” part of the cycle when there is a sufficiently small imposed velocity at the load point. Thus the physically irrelevant time scale (L/Vduring inertial controlled motion,T/2π during quasi‐static motion) is much shorter than the relevant scale, which is troublesome from a numerical point of view as it is the shorter time scale which constrains allowable step size. We propose efficient numerical procedures to deal with such response, in which the full equations with inertia and state relaxation are solved only in a transition regime when both time scales are significant. We show results for several friction laws, all having history dependence based on a single evolving state variable and all having properties that ∂τ/∂V> 0 for instantaneous changes inV, that τ evolves towards τss(V) as exp (−δ/L) with ongoing slip δ whenV= const, and that dτss(V)/dV< 0 except possibly at highV. During the dynamic instabilities we find that motion continues at a nearly steady state condition, τ≃τss(V), until dynamic overshoot becomes so significant that “arrest” begins. In the arrest stage,Vdrops rapidly to very much lower values (never zero in our models) under nearly fixed state conditions, and then the long quasi‐static “stick” phase of the motion begins again.