Instabilities in Dynamic Anti-plane Sliding of an Elastic Layer on a Dissimilar Elastic Half-Space

The stability of dynamic anti-plane sliding at an interface between an elastic layer and an elastic half-space with dissimilar elastic properties is studied. Friction at the interface is assumed to follow a rate- and state-dependent law, with a positive instantaneous dependence on slip velocity and a rate weakening behavior in the steady state. The perturbations at the interface are of the form exp(ikx ₁+pt), where k is the wavenumber, x ₁ is the coordinate along the interface, p is the time response to the perturbation and t is time. A key feature of the problem is that interfacial waves both in freely slipping contact as well as in bonded contact exist for the problem. Attention is focused on the role of the interfacial waves on slip stability. Instabilities are plotted in the $\operatorname{Re} (pL/V_{o})$ versus $\operatorname{Im} (p/|k|c_{s})$ plane, where L is a length scale in the friction law, V _o is the unperturbed slip velocity and c _s is the shear wave speed of the layer. Stability of both rapid and slow slip is studied. The results show one mechanism by which instabilities occur is the destabilization by friction of the interfacial waves in freely slipping contact/bonded contact. This occurs even in slow sliding, thus confirming that the quasi-static approximation is not valid for slow sliding. The effect of material properties and layer thickness on the stability results is studied.     

Dissipative interface waves and the transient response of a three-dimensional sliding interface with Coulomb friction

We investigate the linearized response of two elastic half-spaces sliding past one another with constant Coulomb friction to small three-dimensional perturbations. Starting with the assumption that friction always opposes slip velocity, we derive a set of linearized boundary conditions relating perturbations of shear traction to slip velocity. Friction introduces an effective viscosity transverse to the direction of the original sliding, but offers no additional resistance to slip aligned with the original sliding direction. The amplitude of transverse slip depends on a nondimensional parameter η=c_sτ₀/μv₀, where τ₀ is the initial shear stress, 2v₀ is the initial slip velocity, μ is the shear modulus, and c_s is the shear wave speed. As η→0, the transverse shear traction becomes negligible, and we find an azimuthally symmetric Rayleigh wave trapped along the interface. As η→∞, the inplane and antiplane wavesystems frictionally couple into an interface wave with a velocity that is directionally dependent, increasing from the Rayleigh speed in the direction of initial sliding up to the shear wave speed in the transverse direction. Except in these frictional limits and the specialization to two-dimensional inplane geometry, the interface waves are dissipative. In addition to forward and backward propagating interface waves, we find that for η>1, a third solution to the dispersion relation appears, corresponding to a damped standing wave mode. For large-amplitude perturbations, the interface becomes isotropically dissipative. The behavior resembles the frictionless response in the extremely strong perturbation limit, except that the waves are damped. We extend the linearized analysis by presenting analytical solutions for the transient response of the medium to both line and point sources on the interface. The resulting self-similar slip pulses consist of the interface waves and head waves, and help explain the transmission of forces across fracture surfaces. Furthermore, we suggest that the η→∞ limit describes the sliding interface behind the crack edge for shear fracture problems in which the absolute level of sliding friction is much larger than any interfacial stress changes.     

Analysis of dynamic instability of interfacial slip waves based on the surface impedance tensor

A new method relying on the Stroh formulism and the theory of the surface impedance tensor was developed to investigate the dynamic instability of interfacial slip waves. The concept of the surface impedance tensor was extended to the case where the wave speed is of a complex value, and the boundary conditions at the frictionally contacting interface were expressed by the surface impedance tensor. Then the boundary value problem was transformed to searching for zeroes of a complex polynomial in the unit circle. As an example, the steady frictional sliding of an elastic half-space in contact with a rigid flat surface was considered in details. A quartic complex characteristic equation was derived and its solution behavior in the unit circle was discussed. An explicit expression for the instability condition of the interfacial slip waves was presented.     

Computation of interface wave motions by reciprocity considerations

The current theoretical study deals with computation of Stoneley waves along a solid–solid interface and Scholte waves (also called Scholte-Gogoladze) along a solid–liquid interface by reciprocity considerations. Closed-form solutions of the wave motions generated by a time-harmonic line load applied in two bonded elastic half-spaces of different material properties are derived in a simple manner. In order to perform direct applications of reciprocity theorems, we introduce in this article new expressions for the displacements of free interface waves. Reciprocity relations between an actual state, interface wave motion generated by a time-harmonic line load, and a virtual state, an appropriately chosen free wave traveling along the interface, are derived. Scattered amplitudes of Stoneley waves and Scholte waves due to the load are thus computed. To show application of the obtained results, scattering of Stoneley wave by a delamination at the interface is then studied.     

Effects of interfacial friction on flaw tolerant adhesion between two dissimilar elastic solids

Flaw tolerance refers to a state in which a pre-existing crack-like flaw does not propagate even as the material is stretched to failure near its theoretical strength. Such an optimal scenario can be achieved when the characteristic length scale is reduced to below a critical value. So far, the critical conditions to achieve flaw tolerance have been discussed mostly for homogeneous materials or for two dissimilar materials in frictionless or perfectly bonded adhesion. In this paper, we consider the role of friction in flaw tolerant adhesion between two dissimilar elastic solids. We adopt a frictional contact model in which slip is allowed wherever the shear stress along the interface reaches a threshold value defined as the friction strength. The critical length scale for flaw tolerance is derived analytically for a penny-shaped crack and for an external circular crack. Compared to the cases of frictionless contact, we find that interfacial friction can reduce the critical length scales for flaw tolerance by up to 12.5%.     

General results for complete contacts subject to oscillatory shear

In this paper, we consider the general interfacial characteristics of a square elastic block, pressed onto an elastically similar half-plane by a constant normal force, and subjected to oscillatory shear. It is found that there is a critical coefficient of friction, 0.543, above which the contact is permanently stuck along its entire length for a shearing force below about 55% of that needed to cause sliding. For shearing forces above this, the contact interface will either shakedown to a fully adhered state (depending on the degree of reversal of the shear loading) or will exhibit cyclic slip at an interior point. If the coefficient of friction is below 0.543, the application of normal load alone will produce equal and opposite slip zones attached to the contact edges. The subsequent imposition of a shear force causes the leading edge slip zone to increase in length while the presence of residual slipping tractions at the trailing edge causes the trailing edge to lift off. Under oscillatory loading, the contact edges cycle between slip and separation over a minute region while an interior point may exhibit cyclic slip if the loading history is sufficiently demanding. The results found are of practical relevance to the study of fretting fatigue of complete contacts, such as some types of spline joint.     

Uniqueness of Stoneley waves in pre-stressed incompressible elastic media

The main aim of this paper is to prove, for the general case, the uniqueness of Stoneley waves propagating along the bonded interface of two pre-stressed incompressible elastic half-spaces. In order to do that the authors have used the complex function method. By this approach, it is shown that the secular equation of Stoneley waves in pre-stressed incompressible elastic half-spaces has at most one solution in the complex plane. This says that if a Stoneley wave exists, then it is unique.     

Transition from stick to slip in Hertzian contact with “Griffith” friction: The Cattaneo–Mindlin problem revisited

Classically, the transition from stick to slip is modelled with Amonton–Coulomb law, leading to the Cattaneo–Mindlin problem, which is amenable to quite general solutions using the idea of superposing normal contact pressure distributions – in particular superposing the full sliding component of shear with a corrective distribution in the stick region. However, faults model in geophysics and recent high-speed measurements of the real contact area and the strain fields in dry (nominally flat) rough interfaces at macroscopic but laboratory scale, all suggest that the transition from ‘static’ to ‘dynamic’ friction can be described, rather than by Coulomb law, by classical fracture mechanics singular solutions of shear cracks. Here, we introduce an ‘adhesive’ model for friction in a Hertzian spherical contact, maintaining the Hertzian solution for the normal pressures, but where the inception of slip is given by a Griffith condition. In the slip region, the standard Coulomb law continues to hold. This leads to a very simple solution for the Cattaneo–Mindlin problem, in which the “corrective” solution in the stick area is in fact similar to the mode II equivalent of a JKR singular solution for adhesive contact. The model departs from the standard Cattaneo–Mindlin solution, showing an increased size of the stick zone relative to the contact area, and a sudden transition to slip when the stick region reaches a critical size (the equivalent of the pull-off contact size of the JKR solution). The apparent static friction coefficient before sliding can be much higher than the sliding friction coefficient and, for a given friction fracture “energy”, the process results in size and normal load dependence of the apparent static friction coefficient. Some qualitative agreement with Fineberg's group experiments for friction exists, namely the stick–slip boundary quasi-static prediction may correspond to the arrest of their slip “precursors”, and the rapid collapse to global sliding when the precursors arrest front has reached about half the interface may correspond to the reach of the “critical” size for the stick zone.     

“Frictionless” and “frictional” ThermoElastic Dynamic Instability (TEDI) of sliding contacts

Recently, we found that a new form of coupled instability, named ThermoElastic Dynamic Instability (TEDI), can occur by interaction between frictional heating and the natural dynamic modes of sliding bodies. This is distinct from the classical dynamic instabilities (DI) which is produced by an interaction between the frictional forces at the sliding interface and the natural modes of vibration of the bodies if the friction coefficient is sufficiently high, and also from ThermoElastic Instability (TEI), which is due to the interaction of frictional heating and thermal expansion, leading for example to low pitched brake noise above some critical speed. This result was relative to an highly idealized system, comprising an elastic layer sliding over a rigid plane including both dynamic and thermoelastic effects, but neglecting shear waves at the interface due to frictional tractions (from which the denomination “frictionless TEDI”). We demonstrate here that including these shear waves destabilizes both the shear and dilatational vibration modes of the system at arbitrarily small friction coefficients and speeds, where DI and TEI are predicted to be stable. A detailed study of the new modes and transient simulations show that for low pressures and high speed, the system tends towards the results of the previous model (“frictionless TEDI”), i.e. the tendency to a state in which the layer bounces over the plane, with alternating periods of sliding contact and separation. In the case of low speeds and high pressures, viceversa, the system is dominated by the modes near the resonance of the shear and dilatational modes, with a resulting complex behaviour, but generally leading to stick-slip regimes, reducing the jumping mode of “frictionless TEDI”, because stick reduces or stops frictional heating production.     

Mechanics of advancing pin-loaded contacts with friction

This paper considers finite friction contact problems involving an elastic pin and an infinite elastic plate with a circular hole. Using a suitable class of Green's functions, the singular integral equations governing a very general class of conforming contact problems are formulated. In particular, remote plate stresses, pin loads, moments and distributed loading of the pin by conservative body forces are considered. Numerical solutions are presented for different partial slip load cases. In monotonic loading, the dependence of the tractions on the coefficient of friction is strongest when the contact is highly conforming. For less conforming contacts, the tractions are insensitive to an increase in the value of the friction coefficient above a certain threshold. The contact size and peak pressure in monotonic loading are only weakly dependent on the pin load distribution, with center loads leading to slightly higher peak pressure and lower peak shear than distributed loads. In contrast to half-plane cylinder fretting contacts, fretting behavior is quite different depending on whether or not the pin is allowed to rotate freely. If pin rotation is disallowed, the fretting tractions resemble half-plane fretting tractions in the weakly conforming regime but the contact resists sliding in the strongly conforming regime. If pin rotation is allowed, the shear traction behavior resembles planar rolling contacts in that one slip zone is dominant and the peak shear occurs at its edge. In this case, the effects of material dissimilarity in the strongly conforming regime are only secondary and the contact never goes into sliding. Fretting tractions in the forward and reversed load states show shape asymmetry, which persists with continued load cycling. Finally, the governing integro-differential equation for full sliding is derived; in the limiting case of no friction, the same equation governs contacts with center loading and uniform body force loading, resulting in identical pressures when their resultants are equal.     

A new secular equation for slip waves along the interface of two dissimilar anisotropic elastic half-spaces in sliding contact

We consider wave propagation along the interface of two dissimilar anisotropic elastic half-spaces that are in sliding contact. A new secular equation is obtained that covers all special cases in one equation. One special case is when a Rayleigh wave (called the R$R$-wave) can propagate in both half-spaces with the same wave speed. Another special case is when a slip wave (called the S$S$-wave) can propagate in each of the half-spaces with the same wave speed. If a Rayleigh wave and a slip wave can propagate in one of the half-spaces it is called the RS-wave. In this case an interfacial slip wave exists in which the other half-space is at rest unless an RS-wave can also propagate in the other half-space. The results for general anisotropic elastic materials are applied to orthotropic materials.     

Sliding Contact with Friction at Arbitrary Constant Speeds on a Pre-Stressed Highly Elastic Half-Space

Resisted by Coulomb friction, a rigid indentor slides at a constant arbitrary speed on a generalized neo-Hookean half-space under pre-stress. A dynamic steady-state situation in plane strain is assumed, and is treated as the superposition of contact-triggered infinitesimal deformations upon finite deformations due to pre-stress. Exact solutions are presented for both deformations, and the infinitesimal component exhibits the anisotropy typically induced by pre-stress, and wave speeds that are sensitive to pre-stress. In view of the unilateral constraints of contact, these and other critical speeds define the sliding speed ranges for physically-acceptable solutions. In particular, a Rayleigh speed is the upper bound for subsonic sliding. Solutions are further constrained by the unilateral requirement that contact zone shear must oppose indentor/half-space slip. The generic parabolic indentor is used for illustration, and it is found that traction continuity at the contact zone leading edge is lost for supersonic sliding and at the single sliding speed allowed in the frictionless limit in the trans-sonic range. A range of acceptable pre-stresses is also identified; for pre-stresses that lie out of range, either a negative Poisson effect occurs, or the Rayleigh wave disappears, thereby precluding sliding in the subsonic range. This revised version was published online in August 2006 with corrections to the Cover Date.     

一般各向异性单侧接触界面上波的反射和折射

研究简谐弹性波在一般各向异性介质单侧接触界面上的反射和折射问题．利用Fouier分析方法将非线性Coulomb摩擦接触边界波动问题化为一组代数方程．给出了确定局部分离、滑移和粘着区的思路和方法及各区域的解；讨论了出现界面局部分离和滑移的条件．对特定材料组合情况进行了详细数值计算，给出了界面力、相对滑移速度、张开位移、高频谐波的反射折射系数等特征参量；考察了平面和反平面波动的耦合及整体滑移等．其中关于高频谐波的结果可对已有实验结果给出很好的定性解释．在大多数情况下，即使对摩擦系数无穷大的粘滞接触界面，分离区端部也总是存在一个很小的滑移区。     

Re-polarization of elastic waves at a frictional contact interface—II. Incidence of a P or SV wave

This is Part II of a two-part paper which analyses the re-polarization of elastic waves at a frictional contact interface between two solids. The re-polarization of SH waves was solved in Part I by the use of the Fourier analysis. Here, in Part II, we consider the re-polarization of P or SV waves. It is assumed that the two solids are pressed together and, at the same time, loaded by anti-plane and in-plane shearing traction. If the incident wave is sufficiently strong, localized separation and slip may take place at the interface. As a result, the incident in-plane wave is re-polarized at the interface so that the anti-plane waves (SH waves) are induced. Using the method similar to that of Part I and considering the boundary conditions involving separation and slip, we manage to reduce the problem to a set of algebraic equations coupled with simple integral equations. An iterative method is developed based on the solution to the perfectly bonded interface. The locations and sizes of the separation and slip zones, the interface traction, the slip velocities, the global sliding velocities and the energy dissipation and partition are displayed for the case of two identical materials. It is found that the separation zones and the gaps are independent of the induced waves.     

A modified torsional kolsky bar for investigating dynamic friction

This paper introduces an experiment to investigate dry sliding resistance of frictional interfaces at normal pressures up to 100 MPa, slip speeds up to 10 m/s and slip distances of approximately 10 mm. This new apparatus involves a novel modification of the conventional torsional Kolsky bar apparatus, employed extensively in the past for investigating high strain rate behavior of engineering materials. The new experimental configuration represents a significant improvement over conventional tribology experiments because it uses elastic torsional waves with a superimposed static compressive force to control the interfacial traction. Moreover, the apparatus allows critical frictional parameters such as the interfacial sliding resistance, slip speeds and slip without the use of transducers at the frictional interface. The usefulness of the device is demonstrated by presenting results of high-speed friction on 6061-T6 Al/1018 steel and Carpenter Hampden tool steel/7075-T6 Al tribo pairs.     

A three-layer structure model for the effect of a soft middle layer on Love waves propagating in layered piezoelectric systems

A three-layer structure model is proposed for investigating the effect of a soft elastic middle layer on the propagation behavior of Love waves in piezoelectric layered systems, with "soft" implying that the bulk-shear-wave velocity of the middle layer is smaller than that of the upper sensitive layer. Dispersion equations are obtained for unelectroded and traction-free upper surfaces which, in the limit, can be reduced to those for classical Love waves. Systematic parametric studies are subsequently carried out to quantify the effects of the soft middle layer upon Love wave propagation, including its thickness, mass density, dielectric constant and elastic coefficient. It is demonstrated that whilst the thickness and elastic coefficient of the middle layer affect significantly Love wave propagation, its mass density and dielectric constant have negligible influence. On condition that both the thickness and elastic coefficient of the middle layer are vanishingly small so that it degenerates into an imperfectly bonded interface, the three-layer model is also employed to investigate the influence of imperfect interfaces on Love waves propagating in piezoelectric layer/elastic substrate systems. Upon comparing with the predictions obtained by employing the traditional shear-lag model, the present three-layer structure model is found to be more accurate as it avoids the unrealistic displacement discontinuity across imperfectly bonded interfaces assumed by the shearlag model, especially for long waves when the piezoelectric layer is relatively thin.     

水工圆形隧洞围岩衬砌摩擦滑动接触的新解法

在实际工程中,围岩和衬砌接触时,它们之间并非完全光滑,也并非可以承受任意大的摩擦力.如果围岩与衬砌之间的剪应力大于所能承受的最大静摩擦力,接触面间将发生切向滑动,定义接触面上产生最小滑动量的状态为衬砌的真实工作状态,这种接触即为摩擦滑动接触.以库仑摩擦模型模拟围岩和衬砌之间的摩擦滑动接触,在考虑支护滞后效应的前提下,利用平面弹性复变函数方法列出了应力边界条件、应力连续条件以及位移连续条件的方程, 再结合最优化理论,建立了具有一般性的摩擦滑动接触解法.在利用混合罚函数法求解最优化问题的过程中,减少了设计变量的个数,极大地简化了优化模型,提升了优化过程的迭代速度以及优化结果的精度.以此为基础,获得了围岩和衬砌相互作用下圆形水工隧洞的应力解析解.该方法可以求解光滑接触和完全接触两种极限情况,具有一般性.同时,利用一种精确的计算方法得到了不同情况下满足完全接触条件摩擦系数的阈值,还分析了衬砌和围岩边界上切向应力的变化规律.      

Fracture initiation in the contact region under shear

In the contact region between sliding elastic bodies, there are subregions where the interacting shores are bonded and subregions where they can slide along each other. It is convenient to interpret the latter as transverse shear cracks with slip resistance forces acting on their closed shores. In the end regions of such a crack, stress concentration may lead to fracture initiation in the contacting bodies. Experimental results and an analytic model of the phenomenon are given for a situation where the fracture intersects the contact plane tilted with respect to the direction of the loads.     

An example of stick–slip and stick–slip–separation waves

The dynamical problem of a brake-like mechanical system composed of an elastic cylindrical tube with Coulomb's friction in contact with a rigid and rotating cylinder is considered. This model problem enables us to give an example of non-trivial periodic solutions in the form of stick–slip or stick–slip–separation waves propagating on the contact surface. A semi-analytical analysis of stick–slip waves is obtained when the system of governing equations is reduced by condensation to a simpler system involving only the contact displacements. This reduced system, of only one space variable in addition to time, can be solved almost analytically and gives some interesting informations on the existence and the characteristics of stick–slip waves such as the wave numbers on the circumference, stick and slip proportions, wave celerities, tangential and normal forces. It is shown in particular that the stick–slip–separation solutions would occur for small normal pressures or high rotational speeds. Since the analytical discussion becomes cumbersome in this case, a second approach based on numerical analysis by the finite element method is performed. The existence and the characteristics of stick–slip and stick–slip–separation waves are discussed numerically.     

On formulas for the velocity of Stoneley waves propagating along the loosely bonded interface of two elastic half-spaces

Formulas for the velocity of Stoneley waves propagating along the loosely bonded interface of two isotropic elastic half-spaces are derived using the complex function method. The derivation also shows that if a Stoneley wave exists, then it is unique. By using the obtained formulas, we can easily reproduce the numerical results previously obtained by Murty [G. S. Murty, Phys. Earth Planet. Interiors 11 (1975), 65–79.] by directly solving the secular equation.