Rapid Sliding Indentation with Friction of a Pre-Stressed Thermoelastic Material

A rigid indentor travels with a constant speed over the surface of an isotropic thermoelastic half-space. Friction exists between the indentor and half-space, and the latter is initially in equilibrium at a uniform temperature under a uniform normal pre-stress. This pre-stress, below but near yield, is assumed to produce deformations that dominate the additional deformations due to indentation. Thus, the process is treated as one of small deformations superposed upon (relatively) large. The governing equations for the superposed deformation are those of nonisotropic coupled thermoelasticity. A steady-state two-dimensional study uses robust asymptotic analytical solutions to reduce the associated mixed boundary value problem to a classical singular integral equation which can be solved analytically. The solutions show that the pre-stress-induced de facto nonisotropy alters the values of the rotational and dilatational wave and Rayleigh speeds in the half-space and, in the case of a compressive pre-stress, generates a second, lower, critical speed. In addition, pre-stress generates noncritical sliding speeds at which the friction-dependent integral equation eigenvalue changes sign. For purposes of illustration, expressions for the half-space surface temperature change and its average over the contact zone, the equations necessary to determine contact zone size and location, the resultant moment on the indentor, and the maximum compressive stress on the contact zone are presented for a parabolic indentor. This revised version was published online in August 2006 with corrections to the Cover Date.     

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.     

Rapid sliding indentation with friction on a transversely isotropic thermoelastic half-space

A rigid insulated die slides at a constant sub-critical speed on a transversely isotropic half-space in the presence of friction. In a two-dimensional analysis of the dynamic steady-state, the coupled equations of thermoelasticity are invoked. All elements of the Coulomb friction model are strictly enforced, thus giving rise to auxiliary conditions, including two unilateral constraints.Robust asymptotic forms of an exact solution to a related problem with unmixed boundary conditions lead to analytical solutions for the sliding indentation problem. The solution expressions, abetted by calculations for zinc, show the role of frictional heating on the half-space surface. The effects of friction and sliding speed on contact zone size and location and average contact zone temperature are also studied.The analysis is aided by factoring procedures that simplify the complicated forms that arise in anisotropic elasticity. A scheme that renders expressions for roots of certain irrational functions analytic to within a single quadrature also plays a role.     

Steady sliding frictional contact problem for a 2d elastic half-space with a discontinuous friction coefficient and related stress singularities

The steady sliding frictional contact problem between a moving rigid indentor of arbitrary shape and an isotropic homogeneous elastic half-space in plane strain is extensively analysed. The case where the friction coefficient is a step function (with respect to the space variable), that is, where there are jumps in the friction coefficient, is considered. The problem is put under the form of a variational inequality which is proved to always have a solution which, in addition, is unique in some cases. The solutions exhibit different kinds of universal singularities that are explicitly given. In particular, it is shown that the nature of the universal stress singularity at a jump of the friction coefficient is different depending on the sign of the jump.     

Application of a contact model due to L. A. galin to problems of contact interaction between stringers and massive deformable bodies

L. A. Galin’s contact model for a narrow beam bending on an elastic half-space and Melan’s contact model for a stringer are used to consider two problems of contact interaction between one or two identical symmetrically loaded stringers with small rectangular cross-sections and an elastic half-space. The basic characteristics of these problems are expressed by explicit formulas, and the results of their numerical analysis are given as well.     

Bonded contact of a flexible elliptical disk with a transversely isotropic half-space

The article concerns the problem of bonded contact of a thin, flexible elliptical disk with a transversely isotropic half-space. Three different cases of loading have been considered: (a) the disk is loaded by a transverse force, whose line of action passes through the center of the disk and lies in the plane of the disk; (b) the disk is subjected to a rotation by a torque, whose axis is perpendicular to the surface of the half-space; (c) the half-space with the bonded disk is under uniform stress field at infinity. The problem corresponding to all three cases is reduced, in a unified manner, to a set of coupled two-dimensional integral equations. Closed-form solutions for these equations have been obtained by using Galins theorem.     

Rolling contact of a rigid sphere/sliding of a spherical indenter upon a viscoelastic half-space containing an ellipsoidal inhomogeneity

In this paper, the frictionless rolling contact problem between a rigid sphere and a viscoelastic half-space containing one elastic inhomogeneity is solved. The problem is equivalent to the frictionless sliding of a spherical tip over a viscoelastic body. The inhomogeneity may be of spherical or ellipsoidal shape, the later being of any orientation relatively to the contact surface. The model presented here is three dimensional and based on semi-analytical methods. In order to take into account the viscoelastic aspect of the problem, contact equations are discretized in the spatial and temporal dimensions. The frictionless rolling of the sphere, assumed rigid here for the sake of simplicity, is taken into account by translating the subsurface viscoelastic fields related to the contact problem. Eshelby's formalism is applied at each step of the temporal discretization to account for the effect of the inhomogeneity on the contact pressure distribution, subsurface stresses, rolling friction and the resulting torque. A Conjugate Gradient Method and the Fast Fourier Transforms are used to reduce the computation cost. The model is validated by a finite element model of a rigid sphere rolling upon a homogeneous vciscoelastic half-space, as well as through comparison with reference solutions from the literature. A parametric analysis of the effect of elastic properties and geometrical features of the inhomogeneity is performed. Transient and steady-state solutions are obtained. Numerical results about the contact pressure distribution, the deformed surface geometry, the apparent friction coefficient as well as subsurface stresses are presented, with or without heterogeneous inclusion.     

Two-dimensional generalized thermoelasticity problem for a half-space subjected to ramp-type heating

In this paper, we will consider a half-space filled with an elastic material, which has constant elastic parameters. The governing equations are taken in a unified system from which the field equations for coupled thermoelasticity as well as for generalized thermoelasticity can be easily obtained as particular cases. A linear temperature ramping function is used to more realistically model thermal loading of the half-space surface. The medium is assumed initially quiescent. Laplace and Fourier transform techniques are used to obtain the general solution for any set of boundary conditions. The general solution obtained is applied to a specific problem of a half-space subjected to ramp-type heating. The inverse Fourier transforms are obtained analytically while the inverse Laplace transforms are computed numerically using a method based on Fourier expansion techniques. Some comparisons have been shown in figures to estimate the effect of the ramping parameter of heating with different theories of thermoelasticity.     

On the stability of a plate under longitudinal compression on a two-layer half-space with lower layer prestressed by gravity forces

In the plane (plane strain) and axially symmetric statements, we study the problem of stability, under the action of longitudinal compressing forces, of an infinite elastic plate in two-sided contact with an elastic half-space. The upper layer of finite depth is described by the usual equations of linear theory of elasticity; the lower layer, which is geometrically nonlinear, incompressible, and infinite in depth, is prestressed by gravity forces. The total adhesion between the layer of finite depth and the lower half-space is realized. It is also assumed that the same adhesion takes place between the upper layer of the half-space and the plate with the contact tangential stresses taken into account.The results can be used to calculate the working capacity of coated bodies and layered composites and in problems of geophysics.The problem of stability of an infinite elastic plate under longitudinal compression under conditions of two-sided contact with an elastic base was studied earlier in the monograph [1] (Fuss-Winkler base) and in [2–4].     

Action of an elliptic punch on a transversally isotropic half-space

The 3D contact problem on the action of a punch elliptic in horizontal projection on a transversally isotropic elastic half-space is considered for the case in which the isotropy planes are perpendicular to the boundary of the half-space. The elliptic contact region is assumed to be given (the punch has sharp edges). The integral equation of the contact problem is obtained. The elastic rigidity of the half-space boundary characterized by the normal displacement under the action of a given lumped force significantly depends on the chosen direction on this boundary. In this connection, the following two cases of location of the ellipse of contact are considered: it can be elongated along the first or the second axis of Cartesian coordinate system on the body boundary. Exact solutions are obtained for a punch with base shaped as an elliptic paraboloid, and these solutions are used to carry out the computations for various versions of the five elastic constants. The structure of the exact solution is found for a punch with polynomial base, and a method for determining the solution is proposed.     

Contact analysis in presence of spherical inhomogeneities within a half-space

This paper presents a fast method of solving contact problems when one of the mating bodies contains multiple heterogeneous inclusions, and numerical results are presented for soft or stiff inhomogeneities. The emphasis is put on the effects of spherical inclusions on the contact pressure distribution and subsurface stress field in an elastic half-space. The computing time and allocated memory are kept small, compared to the finite element method, by the use of analytical solution to account for the presence of inhomogeneities. Eshelby’s equivalent inclusion method is considered in the contact solver. An iterative process is implemented to determine the displacements and stress fields caused by the eigenstrains of all spherical inclusions. The proposed method can be seen as an enrichment technique for which the effect of heterogeneous inclusions is superimposed on the homogeneous solution in the contact algorithm. 3D and 2D Fast Fourier Transforms are utilized to improve the computational efficiency. Configurations such as stringer and cluster of spherical inclusions are analyzed. The effects of Young’s modulus, Poisson’s ratio, size and location of the inhomogeneities are also investigated. Numerical results show that the presence of inclusions in the vicinity of the contact surface could significantly changes the contact pressure distribution. From a numerical point of view the role of Poisson’s ratio is found very important. One of the findings is that a relatively ‘soft’ and nearly incompressible inclusion – for example a cavity filled with a liquid – can be more detrimental for the stress state within the matrix than a very hard inclusion with a classical Poisson’s ratio of 0.3.     

Steady and transient Green’s functions for anisotropic conduction in an exponentially graded solid

The problem of the determination of Green’s function in conduction for a rectilinearly anisotropic solid with an exponential grading along a certain direction is studied. Domains of an unbounded space and a half-space, either three-dimensional or two-dimensional, are considered. Along the boundary of the domain, homogeneous boundary conditions of the first and second kinds are imposed. We find interestingly that, under this specific type of grading, the Green’s functions permit an algebraic decomposition, which will in turn greatly simplify the formulation. The method of Fourier transform is employed for the Green’s function for a half-space or a half-plane. Although the derivation process is quite tedious, we show analytically that the inverse transform can be found exactly and their resulting expressions are surprisingly neat and compact. In addition, both steady-state and transient-state field solutions are considered. By taking Laplace transform with respect to the time variable, we show that the mathematical frameworks for the steady-state and transient-state Green’s functions are entirely analogous. Thereby, the transient-state Green’s function is readily obtained by taking Laplace inverse transform back to the time domain. These derived fundamental solutions will serve as benchmark results for modeling some inhomogeneous materials. In the absence of grading term, we have verified analytically that our solutions agree exactly with previously known Green’s functions for homogeneous media.     

Three-dimensional analytical elasticity solution for loaded functionally graded coated half-space

In the present paper, a three-dimensional problem of elasticity of normal and tangential loading of surface of the functionally graded coated half-space is considered. In case when Poisson's ratio is constant and the Young's modulus is a power or exponential function of the distance from the surface of the half-space, analytical solution using Fourier transform is obtained. Stress field due to Hertz contact pressure in an elliptical region are studied as a function of the parameter b/a (where a and b are axes of the contact ellipse) and coating thickness.     

Indentation of an elastic layer by a rigid cylinder

The Green’s functions for the indentation of an elastic layer resting on or bonded to a rigid base by a line load are found efficiently and accurately by a combination of contour integration with a series expansion for small arguments. From the form of the equations it is clear that the function is oscillatory when the layer is free to slip over the base, but for the bonded layer, the function simply decays to zero after a single overshoot.The deformation due to pressure distributions of the form of the product of a polynomial with an elliptical (“Hertzian”) term is calculated and the coefficients chosen to match the indentation shape to that of a cylindrical indenter. The resulting pressure distributions behave much as in Johnson’s approximate theory, becoming parabolic instead of elliptical as the ratio b/d of contact width to layer thickness increases, or, for the bonded incompressible (ν = 1/2) layer, becoming bell-shaped for very large b/d.The relation between the approach δ and the contact width b curves has been investigated, and some anomalies in published asymptotic equations noted and, perhaps, resolved.A noticeable feature of our method is that, unlike previous solutions in which the full mixed boundary value problem (given indenter shape / stress-free boundary) has been solved, the bonded incompressible solid causes no problems and is handled just as for lower values of Poisson’s ratio.     

Plane contact problem for periodic laminated composite involving frictional heating

Summary The stationary problem of a rigid thermally insulated punch sliding over the boundary surface of a periodic two-layered thermoelastic half-space is considered. The heat generated in the contact area is assumed to be caused by frictional forces. The problem is formulated within the framework of thermoelasticity with microlocal parameters, and it is reduced to a system of two integral equations, which is solved numerically. The effects connected with the composite structure are analyzed.     

The effect of fractional parameter on a perfect conducting elastic half-space in generalized magneto-thermoelasticity

Recently, Sherief et al. (Int. J. Solids Struct. 47:269–275, 2010) proposed a model in generalized thermoelasticity based on the fractional order time derivatives. The propagation of electro-magneto-thermoelastic disturbances in a perfectly conducting elastic half-space is investigated in the context of the above fractional order theory of generalized thermoelasticity. There acts an initial magnetic field parallel to the plane boundary of the half-space. Normal mode analysis together with the eigenvalue approach technique is used to solve the resulting non-dimensional coupled governing equations of the problem. The obtained solution is then applied to two specific problems for the half-space, whose boundary is subjected to (i) thermally isolated surfaces subjected to time-dependent compression and (ii) a time-dependent thermal shock and zero stress. The effects of fractional parameter and magnetic field on the variations of different field quantities inside the half-space are analyzed graphically.     

Dynamic frictional indentation of an elastic half-plane by a rigid punch

Dynamic rigid indentation of a linearly elastic half-plane in the presence of Coulomb friction is studied in this paper. A rigid punch, which is either wedge- or parabolic-shaped, is rapidly driven into the deformable body so that stress waves are generated. The contact region is assumed to extend at a constant sub-Rayleigh speed (this situation can be achieved by conveniently specifying the kinetic and geometric characteristics of indentor), whereas, due to symmetry, friction acts in opposing directions on opposite sides of the indentor. As the present exact analysis shows, this sign reversal of the tangential traction along the half-plane surface creates an extra stress-singularity at the changeover point of the boundary conditions (due to symmetry, this point here coincides with the point where the indentor apex makes contact with the half-plane surface). The study exploits the problem's self-similarity by utilizing homogeneous-function techniques previously used by L.M. Brock, along with the Riemann-Hilbert problem analysis. Representative numerical results are given for the wedge indentation case.     

Contact interaction of a rigid heat-conducting punch with an elastic layer involving nonstationary frictional heating

A mathematical formulation is given and a solution is found to the quasistatic contact problem of thermoelasticity for a rigid heat-conducting punch moving over an elastic layer with fixed base. The interaction is accompanied by heating due to frictional forces obeying Amonton’s law. The problem is reduced to a system of integral equations with time-varying limits of integration. The structure of these equations depends on the type of thermal and physical conditions on the contact surface. An algorithm is proposed for the numerical solution of this kind of equations. The variation in the contact pressure and contact area with time is studied __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 12, pp. 35–46, December 2005.     

Mechanics of non-slipping adhesive contact on a power-law graded elastic half-space

In previous work about axisymmetric adhesive contact on power-law graded elastic materials, the contact interface was often assumed to be frictionless, which is, however, not always the case in practical applications. In order to elucidate the effect of friction and the coupling between normal and tangential deformations, in the present paper, the problem of a rigid punch with a parabolic shape in non-slipping adhesive contact with a power-law graded half-space is studied analytically via singular integral equation method. A series of closed-form analytical solutions, which include the frictionless and homogeneous solutions as special cases, are obtained. Our results show that, compared with the frictionless case, the interfacial friction tends to reduce the contact area and the indentation depth during adhesion. The magnitude of the coupling effect depends on both the Poisson ratio and the gradient exponent of the half-space. This effect vanishes for homogeneous incompressible as well as for linearly graded materials but becomes significant for auxetic materials with negative Poisson’s ratio. Furthermore, influence of mode mixity on the adhesive behavior of power-law graded materials, which was seldom touched in literature, is discussed in details.     

Mechanics of indentation for piezoelectric thin films on elastic substrate

Frictionless normal indentation problem of rigid flat-ended cylindrical, conical and spherical indenters on piezoelectric film, which is either in frictionless contact with or perfectly bonded to an elastic half-space (substrate), is investigated. Both conducting and insulating indenters are considered. With Hankel transform, the general solutions of the homogeneous governing equations for the piezoelectric layer and the elastic half-space are presented. Using the boundary conditions for a vertical point force or a point electric charge, and the boundary conditions on the film/substrate interface, the Green’s functions can be obtained by solving sets of simultaneous linear algebraic equations. The solution of the indentation problem is obtained by integrating these Green’s functions over the contact area with unknown surface tractions or electric charge distribution, which will be determined from the boundary conditions on the contact surface between the indenter and the film. The solution is expressed in terms of dual integral equations that are converted to a Fredholm integral equation of the second kind and solved numerically. Numerical examples are also presented. The comparison between two film/substrate bonding conditions is made. It shows that the indentation rigidity of the film/substrate system is lower when the film is in frictionless contact with the substrate. The effects of the Young’s modulus and Poisson’s ratio of the elastic substrate, indenter electrical condition and indenter prescribed electric potential on the indentation responses are presented.