# Contact Problem for an Elastic Layer on an Elastic Half Plane Loaded by Means of Three Rigid Flat Punches

Authors: T. S. Ozsahin; O. Taskıner
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

The frictionless contact problem for an elastic layer resting on an elastic half plane is considered. The problem is solved by using the theory of elasticity and integral transformation technique. The compressive loadsP and Q (per unit thickness in z direction) are applied to the layer through three rigid flat punches. The elastic layer is also subjected to uniform vertical body force due to effect of gravity. The contact along the interface between elastic layer and half plane is continuous, if the value of the load factor, λ, is less than a critical value, λcr. In this case, initial separation loads, λcr and initial separation points, xcr are determined. Also the required distance between the punches to avoid any separation between the punches and the elastic layer is studied and the limit distance between punches that ends interaction of punches is investigated for various dimensionless quantities. However, if tensile tractions are not allowed on the interface, for λ>λcr the layer separates from the interface along a certain finite region. Numerical results for distance determining the separation area, vertical displacement in the separation zone, contact stress distribution along the interface between elastic layer and half plane are given for this discontinuous contact case.


---

## Body

## 1. Introduction

Contact between deformable bodies abounds in industry and everyday life. Because of the industrial importance of the physical processes that take place during contact, a considerable effort has been made in their modeling, analysis, and numerical simulations.The range of application in contact mechanics starts with problems like foundations in civil engineering, where the lift off the foundation from soil due to eccentric forces acting on a building, railway ballasts, foundation grillages, continuous foundation beams, runaways, liquid tanks resting on the ground, and grain silos is considered. Furthermore, foundation including piles as supporting members or the driving of piles into the soil is of interest. Also classical bearing problem of steel constructions and the connection of structural members by bolts or screws are areas in which contact analysis enters the design process in civil engineering [1].A complete analysis of the interaction problem for elastic bodies generally requires the determination of stresses and strain within the individual bodies in contact, together with information regarding the distribution of displacements and stresses at the contact regions.The contact problem for an elastic layer has attracted considerable attention in the past because of its possible application to a variety of structures of practical interest. The layer usually rests on a foundation which may be either elastic or rigid and the body force due to gravity may be neglected [2–27] or effect of gravity may be taken into account [28–35]. Studies regarding the frictionless contact along the interface may be found in [2–15]. Also contact region may exhibit frictional characteristics, giving rise to normal and shear traction at the contact surface [16–27].While initial contact is determined by the geometric features of the bodies, the extent of the contact generally changes not only by the particular loads applied to the bodies but also with the elastic constants of the materials. Due to the bending of the layer under local compressive loads, in the absence of gravity effects, the contact area would decrease to a finite size which is independent of the magnitude of the applied load [2–7]. If the effect of gravity was taken into account, the normal stress along the layer subspace interface will be compressive and the contact is maintained through the frictionless interface unless the compressive applied load exceeds a certain critical value. When the magnitude of the compressive external load exceeds a certain value, a separation will take place between the layer and the foundation. The length of separation region along the interface is unknown, and the problem becomes a discontinuous contact problem [28–35].Interaction between an elastic medium and a rigid punch forms another group of contact problems. Rigid punches may be structural elements such as foundations, beams, and plates of finite or infinite extent resting on idealized linearly deformable elastic media. Here, the shape of the contact region may be known a priori and remains constant, or contact region may be changed due to the shape of punch profile [8–27]. The problem of flat-ended rigid punch has important applications in soil mechanics, particularly in estimating the safety of foundations. The application of the three punches for an elastic layer resting on an elastic half plane in soil mechanics is obvious; for example, the punches can be taken as foundations placed on layered soil. When the foundations are placed on soil, there is a possibility of pressure isobars of adjacent foundations overlapping each other. The soil is highly stressed in the zones of overlapping, or the difference of settlement between two adjacent foundations, commonly referred to as differential settlement may cause damage to the structure. It is possible to avoid overlapping of pressures or differential settlement by installing the foundations at considerable distance apart from each other.In this study, contact problem of the three punches for an elastic layer resting on an elastic half plane is considered according to the theory of elasticity with integral transformation technique. The compressive loadsP and Q (per unit thickness in z direction) are applied to the layer through three rigid flat punches. The width of midmost punch can be different from the other two punches and thickness of the layer is constant, h. The layer is subjected to homogeneous vertical body force due to gravity, ρ1g. All surfaces are frictionless. The layer remains in contact with the elastic half plane where the magnitude of the load factor λ is less than a critical value, λcr(λ=P/ρ1gh2). If λ>λcr, the contact is discontinuous and a separation takes place between the layer and the half plane. A numerical integration procedure is performed for the solution of the problems, and different parameters are researched for various dimensionless quantities for both continuous and discontinuous contact cases. Finally, numerical results are analyzed and conclusions are drawn.


## 2. General Expressions for Stresses and Displacements

Consider a frictionless elastic layer of thicknessh lying on an elastic half plane. The geometry and coordinate system are shown in Figure 1. The governing equations are (1a)μk∇2uk+2μk(κk-1)∂∂x(∂uk∂x+∂vk∂y)=0,(1b)μk∇2vk+2μk(κk-1)∂∂y(∂uk∂x+∂vk∂y)=ρkg,vkvkvkvkvkvkvkvkvkvkvkvkvkvkvk(k=1,2)where ρkg is the intensity of the body force acting vertically in which ρk and g are mass density and gravity acceleration. uk and vk are the x and y components of the displacement vector, μk and κk represent shear modulus and elastic constant of the layer and the half plane, respectively.κk=(3-γk)/(1+γk) for plane stress and κk=(3-4γk) for plane strain. γk is the Poisson ratio (k=1,2). Subscript 1 indicates the elastic layer and subscript 2 indicates the elastic half plane.Figure 1
Elastic layer resting on an elastic half plane loaded by means of three rigid flat punches.u
p and vp represent the displacements for the case in which gravity forces are considered. uh and vh are the displacements when the gravity forces are ignored, and total field of displacements may be expressed as(2a)u=up+uh,(2b)v=vp+vh.Observing thatx=0 is a plane of symmetry, it is sufficient to consider the problem in the region 0≤x≤∞ only. Using the symmetry consideration, the following expressions may be written: (3a)u1(x,y)=-u1(-x,y),(3b)v1(x,y)=v1(-x,y),(3c)u1(x,y)=2π∫0∞ϕ1(α,y)sin(αx)dα,(3d)v1(x,y)=2π∫0∞Ψ1(α,y)cos(αx)dα,where ϕ1 and Ψ1 functions are inverse Fourier transforms of  u1and v1 respectively. Taking necessary derivatives of (3c) and (3d), substituting them into (1a) and (1b), and solving the second-order differential equations, the following equations may be obtained for displacements: (4a)u1h(x,y)=2π∫0∞[(A+By)e-αy+(C+Dy)eαy]sin(αx)dα,(4b)v1h(x,y)=2π∫0∞[[A+(κ1ακ+y)B]e-αy+[-C+(κ1α-y)D]eαy]cos(αx)dα.Using Hooke’s law and (4a) and (4b), stress components which do not include the gravity force may be expressed as follows: (5a)σx1h(x,y)=4μ1π∫0∞[[α(A+By)-(3-κ12)B]e-αy+++++[α(C+Dy)+(3-κi2)D]eαy]×cos(αx)dα,(5b)σy1h(x,y)=4μ1π∫0∞[-[α(A+By)+(κ1+12)B]e-αy+[-α(C+Dy)+(κ1+12)D]eαy]×cos(αx)dα,(5c)τxy1h(x,y)=4μ1π∫0∞[-[α(A+By)+(κi-12)B]e-αy+[α(C+Dy)-(κi-12)D]eαy]×sin(αx)dα.For the case in which gravity force exist, particular part of the displacement components corresponding toρ1g, the following expressions are obtained, that is, special solution of the Navier equations for a layer with a height h: (6a)u1p=3-κ18μ1ρ1gh2x,(6b)v1p=ρ1gy2μ1[(κ1-1)(κ1+1)(y-h)-(κ1+1)8h],(6c)σy1p=ρ1g(y-h),(6d)σx1p=ρ1g(y-h2)(1+κ2)(1+κ2)+2μ2(1+κ1),(6e)τxy1p=0.Considering the orthogonal axes shown in Figure1, displacements will be zero for y=-∞, and if μ2, ν2 are the elastic constants of the half plane, then the homogenous field of displacements and stresses of the elastic half plane may be obtained as (7a)u2h(x,y)=2π∫0∞[(E+Fy)eαy]sin(αx)dα,(7b)v2h(x,y)=2π∫0∞[[-E+(κ2α-y)F]eαy]cos(αx)dα,(7c)σx2h(x,y)=4μ2π∫0∞[[α(E+Fy)+(3-κ22)F]eαy]×cos(αx)dα,(7d)σy2h(x,y)=4μ2π∫0∞[[+(κ2+12)F]-α(E+Fy)+(κ2+12)F]eαy]cos(αx)dα,(7e)τxy2h(x,y)=4μ2π∫0∞[[α(E+Fy)-(κ2-12)F]eαy]×sin(αx)dα,subscript 2 indicates the elastic half plane. Note that the body force acting in the foundation is neglected since it does not disturb the contact pressure distribution. A, B, C, D, E, and F are the unknown constants, which will be determined from the boundary and continuity conditions at y=0 and y=h.


## 3. Case of Continuous Contact

An elastic layer with a height ofh resting on an elastic half plane, shown in Figure 1, is analyzed for unit thickness in z direction. Widths of punches at both sides are similar, (c-b)/h, and each of these punches transmits a concentrated load of Q to the elastic layer. The width of the midmost punch is different, 2a/h and it is subjected to a concentrated load, P. All surfaces are frictionless. Particularly, the initial separation load (λcr) and point (xcr) where the layer separated from the elastic half plane and the variation of the stress distribution between elastic layer and elastic half plane is examined depending on material properties, width of punches, and magnitude of the external loads, P and Q. Due to the different settlement of punches, a separation takes place between punch I and elastic layer, if punches are close enough. Therefore, the critical distance between the punches indicating the initiation of separation between the punch I and the elastic layer is researched and also limit distance between the punches where the interaction of punches ends is investigated.If load factor(λ) is sufficiently small, then the contact along the layer-subspace, y=0,0<x<∞, will be continuous, and A, B, C, D, E, and F must be determined from the following boundary and continuity conditions:(8a)σy1(x,h)={-P(x),0<x<a-Q(x),b<x<c0,a<x<b,c<x<∞,(8b)τxy1(x,h)=0,0<x<∞,(8c)τxy1(x,0)=0,0<x<∞,(8d)τxy2(x,0)=0,0<x<∞,(8e)σy1(x,0)=σy2(x,0),0<x<∞,(8f)∂∂x[v2(x,0)-v1(x,0)]=0,0<x<∞,(8g)∂∂x[v1(x,h)]=0,0<x<a,(8h)∂∂x[v1(x,h)]=0,b<x<c,in which subscripts 1 and 2 indicate relation to the elastic layer and the elastic half plane, respectively.P(x) is the unknown contact pressure under punch I and Q(x) is the unknown contact pressure under punch II, which have not been determined yet. If a separation occurs between the elastic layer and elastic half plane, this will give rise to a discontinuous contact position and the following results for former solution will no longer be valid and new solution will be attained for the latter case.Equilibrium conditions of the problem may be expressed as(9a)∫0aP(x)dx=P2,(9b)∫bcQ(x)dx=Q.Displacement and stress expressions (4a), (4b), (5a)–(5c), (6a)–(6e), and (7a)–(7e) are substituted into boundary conditions (8a)–(8f), and unknown constants A, B, C, D, E,andF are determined in terms of unknown functions P(x) and Q(x). By making use of (8g) and (8h), after some simple manipulations, one may obtain the following singular integral equations for P(x) and Q(x) [36, 37]: (10a)-1πμ1∫0a[k1(x,t)+(1+κ1)4(1t+x+1t-x)]P(t)dt-1πμ1∫bc[k1(x,t)+(1+κ1)4(1t+x+1t-x)]×Q(t)dt=0,0<x<a,(10b)-1πμ1∫0a[k1(x,t)+(1+κ1)4(1t+x+1t-x)]P(t)dt-1πμ1∫bc[k1(x,t)+(1+κ1)4(1t+x+1t-x)]×Q(t)dt=0,b<x<c,where(11)k1(x,t)=∫0∞{{α38(1+κ1)×[[(1+κ2)+μ2μ1(1+κ1)]+e-2αh[4αh(1+κ2)-2μ2μ1(1+κ1)]+e-4αh[-1-κ2+μ2μ1(1+κ1)]]}(Δ)-1-(1+κ1)4{α3/8(1+κ1)}×{sinα(t+x)-sinα(t-x)}dα,
in which
(12)Δ=α3(2{(-1-κ2-μ2μ1(1+κ1))+e-4αh[-1-κ2+μ2μ1(1+κ1)]+e-2αh[4αhμ2μ1(-1-κ1)(1+κ2)(2+4α2h2)+4αhμ2μ1(-1-κ1)]})-1.If evaluated values ofA, B, C, and D in terms of P(x) and Q(x) are substituted into (5b), the expression of the contact stress between elastic layer and half plane may be obtained as
(13)σy1(x,0)=-ρ1gh-1π∫0ak2(x,t)P(t)dt-1π∫bck2(x,t)Q(t)dt,0<x<∞,
in which ρ1 and g are mass density and gravity acceleration, respectively, where
(14)k2(x,t)=∫0∞{α3μ2μ1(1+κ1)×[e-3αh(-1+αh)+e-αh(1+αh)]α3μ2μ1(1+κ1)}(Δ)-1×{cosα(t+x)+cosα(t-x)}dα.To simplify the numerical analysis, the following dimensionless quantities are introduced:(15a)x1=ar1,(15b)t1=as1,(15c)x2=c-b2r2+c+b2,(15d)t2=c-b2s2+c+b2,(15e)g1(s1)=P(as1)P/h,(15f)g2(s2)=Q(((c-b)/2)s2+(c+b)/2)P/h,(15g)α=wh,(15h)λ=Pρ1gh2.Substituting from (15a)–(15h), (9a), (9b) and (10a), (10b) may be expressed as(16a)∫01g1(s1)ahds1=12,(16b)∫-11g2(s2)c-b2hds2=QP,(16c)-1π∫01g1(s1)ahds1[m1(r1,s1)+(1+κ1)4++++++++++×(1a(s1+r1)-1a(s1-r1))]-1π∫-11g2(s2)c-b2hds2×[m2(r1,s2)+(1+κ1)4×(1(((c-b)/2)s2+(c+b)/2)+ar1-1(((c-b)/2)s2+(c+b)/2)-ar1)]=0,0<r1<1,(16d)-1π∫01g1(s1)ahds1×[m3(r2,s1)+(1+κ1)4+++×(1as1+(((c-b)/2)r2+(c+b)/2)+++++-1as1-(((c-b)/2)r2+(c+b)/2))]-1π∫-11g2(s2)c-b2ds2×[1(c-b/2)(s2-r2)m4(r2,s2)+(1+κ1)4+++×(1((c-b)/2)(s2+r2)+(c+b)+++++-1((c-b)/2)(s2-r2))]=0,-1<r2<1,(16e)σy1(x,0)P/h=-1λ-1π∫01m5(r1,s1)g1(s1)ahds1-1π∫-11m6(r2,s2)g2(s2)c-b2hds2,where(17a)m1(r1,s1)=k1(x1,t1),(17b)m2(r1,s2)=k1(x1,t2),(17c)m3(r2,s1)=k1(x2,t1),(17d)m4(r2,s2)=k1(x2,t2),(17e)m5(r1,s1)=k2(x1,t1),(17f)m6(r2,s2)=k2(x2,t2).The index of the integral equations in (16a)–(16e) is +1; so the functions g1(s1) and g2(s2) may be expressed in the following forms:(18a)g1(s1)=G1(s1)(1-s12)-1/2,(0<s1<1),(18b)g2(s2)=G2(s2)(1-s22)-1/2,(-1<s2<1),where G1(s1) is bounded in [0,1] and G2(s2) is bounded in [-1,1]. Then using appropriate Gauss-Chebyshev integration formula given in [36], (16a)–(16d) are replaced by the following algebraic equations:(19a)∑i=1nπWiahG1(s1i)=12,(19b)∑i=1nπWic-b2hG2(s2i)=QP,(19c)-∑i=1nWiG1(s1i)ah[(1a(s1i+r1j)-1a(s1i-r1j))m1(r1j,s1i)+(1+κ1)4+++++++++×(1a(s1i+r1j)-1a(s1i-r1j))]-∑i=1nWiG2(s2i)c-b2h×[(1a(s1i+r1j)-1a(s1i-r1j))m2(r1j,s2i)+(1+κ1)4×(1(((c-b)/2)s2i+(c+b)/2)+ar1j++++-1(((c-b)/2)s2i+(c+b)/2)-ar1j)]=0,(j=1,…,n-1),(19d)-∑i=1nWiG1(s1i)ah×[1as1i-((c-b/2)r2j+(c+b/2))m3(r2j,s1i)+(1+κ1)4×(1as1i+(((c-b)/2)r2j+(c+b)/2)xxxxx×-1as1i-(((c-b)/2)r2j+(c+b)/2))]-∑i=1nWiG2(s2i)c-b2h×[1((c-b)/2)(s2i-r2j)m4(r2j,s2i)+(1+κ1)4+++×(1((c-b)/2)(s2i+r2j)+(c+b)+++++-1((c-b)/2)(s2i-r2j))]=0,(j=1,…,n-1),(19e)W1=Wn=12n-2,Wi=1n-1,(i=2,…,n-1),(19f)s1i=cos(i-12n-1π),s2i=cos(i-1n-1π),(i=1,…,n),(19g)r1j=cos(2j-14n-2π),r2j=cos(2j-12n-2π),(j=1,…,n-1).The unknownsG1(s1i) and G2(s2i)(i=1,…,n) are determined from the system (19a)–(19d). By using (18a), (18b), substituting the results into (16e), and using Gauss-Chebyshev integration formula, the contact stress σy1(x,0)h/P is evaluated. It should be observed that the integral equations (16c) and (16d) are valid provided the contact stress σy1(x,0)h/P is compressive everywhere; that is, 0<x<∞. The critical load factor, λcr and the corresponding location of interface separation, xcr can be determined through the use of the following condition for various dimensionless quantities:
(20)σy1(x,0)P/h=0.


## 4. Case of Discontinuous Contact

Since the interface cannot carry tensile tractions forλ >λcr, there will be separation between the elastic layer and the elastic half plane in the neighborhood of x=xcr on the contact plane y=0, as shown in Figure 1. Assuming that the separation region is described by d<x<e, y=0, where d and e are unknowns and functions of λ, boundary and continuity conditions for the discontinuous contact case are defined as follows:(21a)σy1(x,h)={-P(x),0<x<a-Q(x),b<x<c0,a<x<b,c<x<∞,(21b)τxy1(x,h)=0,0<x<∞,(21c)τxy1(x,0)=0,0<x<∞,(21d)τxy2(x,0)=0,0<x<∞,(21e)σy1(x,0)=σy2(x,0),0<x<∞,(21f)∂∂x[v2(x,0)-v1(x,0)]={φ(x),d<x<e0,0<x<d,e<x<∞,(21g)σy1(x,0)=σy2(x,0)=0,d<x<e,(21h)∂∂x[v1(x,h)]=0,0<x<a,(21i)∂∂x[v1(x,h)]=0,c<x<b.After utilizing the boundary and continuity conditions defined in (21a)–(21f), new values for the constants A, B, C, D, E, and F which appear in (4a), (4b), (5a)–(5c), and (7a)–(7e) may be obtained in terms of new unknown functions P(x), Q(x), and φ(x). Unknown functions are then determined from the conditions (21h)-(21i) which have not yet been satisfied. These conditions give the following system of singular integral equations:(22a)-1πμ1∫0a[k1(x,t)+(1+κ1)4(1t+x-1t-x)]P(t)dt-1πμ1∫bc[k1(x,t)+(1+κ1)4(1t+x-1t-x)]Q(t)dt+1π∫dek2(x,t)φ(t)dt=0,0<x<a,(22b)-1πμ1∫0a[k1(x,t)+(1+κ1)4(1t+x-1t-x)]P(t)dt-1πμ1∫bc[k1(x,t)+(1+κ1)4(1t+x-1t-x)]Q(t)dt+1π∫dek2(x,t)φ(t)dt=0,b<x<c,(22c)1π∫0ak2(x,t)P(t)dt+1π∫bck2(x,t)Q(t)dt-μ1π∫de[k3(x,t)-4μ2/μ1(1+κ2)+μ2/μ1(1+κ1)×(1t+x+1t-x)4μ2/μ1(1+κ2)+μ2/μ1(1+κ1)]φ(t)dt-ρ1gh=0,d<x<e,where kernels k1(x,t) and k2(x,t) are given by (11) and (14) and
(23)k3(x,t)=∫0∞{{-2α3μ2μ1{e-2αh(2+4α2h2)-e-4αh-1}}×(Δ)-1+4μ2/μ1(1+κ2)+μ2/μ1(1+κ1){-2α3μ2μ1{e-2αh(2+4α2h2)-e-4αh-1}}}×{sinα(t+x)+sinα(t-x)}dα,
in which Δ is given by (12).The index of integral equations (22a) and (22b) is +1. On the other hand, the index of the singular integral equation (22c) is −1 due to the physical requirement of smooth contact at the end points d and e [37]. Thus, in solving the problem the two conditions which would account for the unknowns d and e are the consistency condition of integral equation (22c) and the single-valuedness condition
(24)∫deφ(x)dx=0.Defining the following dimensionless quantities(25a)x3=e-d2r3+e+d2,t3=e-d2s3+e+d2,(25b)g3(s3)=μ1φ(((e-d)/2)s3+(e+d)/2)P/h;by making use of (15a)–(15h), the integral equations (22a)–(22c) may be expressed as follows:(26a)-1π∫01g1(s1)ahds1[(1a(s1+r1)-1a(s1-r1))m1*(r1,s1)+(1+κ1)4++++++++++×(1a(s1+r1)-1a(s1-r1))]-1π∫-11g2(s2)c-b2hds2×[1((c-b/2)s2+(c+b/2))-ar1m2*(r1,s2)+(1+κ1)4×(1(((c-b)/2)s2+(c+b)/2)+ar1-1(((c-b)/2)s2+(c+b)/2)-ar1)]+1π∫-11g3(s3)m3*(r1,s3)e-d2hds3=0,0<r1<1,(26b)-1π∫-11g1(s1)ahds1×[-1as1-((c-b/2)r2+(c+b/2))m4*(r2,s1)+(1+κ1)4×(1as1+(((c-b)/2)r2+(c+b)/2)-1as1-(((c-b)/2)r2+(c+b)/2))]-1π∫-11g2(s2)c-b2hds2[m5*(r2,s2)+(1+κ1)4++++++++++++++×(1((c-b)/2)(s2+r2)+(c+b)+++++++++++++++-1((c-b)/2)(s2-r2)(1+κ1)4)]+1π∫-11g3(s3)m6*(r2,s3)e-d2hds3=0,-1<r2<1,(26c)-1π∫-11g1(s1)m7*(r3,s1)ahds1-1π∫-11g2(s2)m8*(r3,s2)c-b2hds2-1π∫-11g3(s3)e-d2hds3×[m9*(r3,s3)-4μ2/μ1(1+κ2)+μ2/μ1(1+κ1)×(1((e-d)/2)(s3+r3)+(e+d)+1((e-d)/2)(s3-r3))]-1λ=0,-1<r3<1,where(27a)m1*(r1,s1)=k1(x1,t1),(27b)m2*(r1,s2)=k1(x1,t2),(27c)m3*(r1,s3)=k2(x1,t3),(27d)m4*(r2,s1)=k1(x2,t1),(27e)m5*(r2,s2)=k1(x2,t2),(27f)m6*(r2,s3)=k2(x2,t3),(27g)m7*(r3,s1)=k2(x3,t1),(27h)m8*(r3,s2)=k2(x3,t2),(27i)m9*(r3,s3)=k3(x3,t3).Similar to (16a), (16b), additional condition (24) may be expressed as
(28)∫-11g3(s3)ds3=0.To solve the system of integral equations, it is found to be more convenient to assume that (26c) as well as (26a) and (26b) has an index +1 [29]; consequently, the function gi(si)  (i=1,…,3) may be expressed in the form (29a)g1(s1)=G1(s1)(1-s12)-1/2,0<s1<1,(29b)gi(si)=Gi(si)(1-si2)-1/2,-1<si<1,(i=2,3),where Gi(si) is a bounded function. In order to insure smooth contact at the end points of the separation region, we then impose the following conditions on G3(s3):
(30)G3(-1)=0,G3(1)=0.Equations (26a)–(26c), (16a), (16b), and (28) can easily be reduced to the following system of linear algebraic equations by employing the appropriate Gauss-Chebyshev integration formula [36]: (31a)-∑i=1nWiG1(s1i)ah[1a(s1i+r1j)m1*(r1j,s1i)+(1+κ1)4+++++++++×1a(s1i+r1j)-1a(s1i-r1j)]-∑i=1nWiG2(s2i)c-b2h×[-1((c-b/2)s2i+(c+b/2))-ar1j)m2*(r1j,s2i)+(1+κ1)4×(1(((c-b)/2)s2i+(c+b)/2)+ar1j-1(((c-b)/2)s2i+(c+b)/2)-ar1j)]+∑i=2n-1WiG3(s3i)e-d2hm3*(r1j,s3i)=0,(j=1,…,n-1),(31b)-∑i=1nWiG1(s1i)ah×[1as1i-(((c-b)/2)r2j+(c+b)/2)m4*(r2j,s1i)+(1+κ1)4+++×(1as1i+(((c-b)/2)r2j+(c+b)/2)+++++-1as1i-(((c-b)/2)r2j+(c+b)/2))]-∑i=1nWiG2(s2i)c-b2h[1(c-b/2)(s2i-r2j)m5*(r2j,s2i)+(1+κ1)4+++++++++++++×(1((c-b)/2)(s2i-r2j)+(c+b)+++++++++++++++±1((c-b)/2)(s2i-r2j))]+∑i=2n-1WiG3(s3i)e-d2hm6*(r2j,s3i)=0(j=1,…,n-1),(31c)-∑i=1nWiG1(s1i)ahm7*(r3j,s1i)-∑i=1nWiG2(s2i)c-b2hm8*(r3j,s2i)-∑i=2n-1WiG3(s3i)e-d2h×[1(e-d/2)(s3i-r3j)m9*(r3j,s3i)-4μ2/μ1(1+κ2)+μ2/μ1(1+κ1)+++×(1((e-d)/2)(s3i-r3j)+(e+d)++++++1((e-d)/2)(s3i-r3j))]-1λ=0,(j=1,…,n-1),(31d)∑i=1nπWiahG1(s1i)=12,(31e)∑i=1nπWic-b2hG2(s2i)=QP,(31f)∑i=2n-1πWiG3(s3i)=0,where Wi, si, and rj are given by (19e)–(19g) (r2=r3,s2=s3). It was shown in [36] that the consistency condition is automatically satisfied if the Gauss-Chebyshev integration formula is used for solving integral equations. Thus, (31a)–(31d) and (31e)-(31f) give 3n equation for 3n unknowns G1(si), G2(si), G3(sj)  (i=1,…,n),  (j=2,…,n-1), d and e. The equation system is nonlinear in d and e; so an interpolation scheme is required for the solution. Selected values of d and e are substituted into (31a)–(31d), and G1(s), G2(s), and G3(s) are obtained which must satisfy (31e), (31f) at the same time for known λ>λcr. If (31e), (31f) are not satisfied, then solution must be repeated with new values of d and e until the (31e), (31f) are satisfied at the same time.It should be noted that (22c) gives the σy1(x,0)h/P outside as well as inside the separation region (e,f). Thus, once the functions G1(s), G2(s), and G3(s) and the constants d and e are determined, contact stress σy1(x,0)h/P may be easily evaluated. The displacement component v*(x,0) in the separation region (d,e), referring to (21f) and (25b), may be obtained from(32a)v*(x,0)=v2(x,0)-v1(x,0)=∫dxg3(t)dt,d<x<e
or
(32b)μ1P/hv*(x,0)=e-d2h∫-1r3G3(s3)ds3,-1<r3<1,
where
(32c)x=e-d2r3+e+d2.Also using appropriate Gauss-Chebyshev integration formula and taking+1 as the index of (32b), the following expression may be written for the vertical displacement in the separation region: (33)μ1P/hv*(x,0)=e-d2h∑i=2k-1WiG3(s3i),(k=2,…,n-1),
where Wi and si are given by (19e)–(19g).


## 5. Results and Discussion

Some of the calculated results obtained from the solution of the continuous and discontinuous contact problems for various dimensionless quantities such asμ2/μ1, a/h, (c-b)/h, Q/P, λ, and (b-a)/h are presented in Figures 2, 3, 4, 5, and 6 and Tables 1, 2, 3, 4, and 5. First separation point between elastic layer and elastic half plane is determined and first separation region is investigated. Contact pressure σy1(x,0)h/P is presented. Depending on load factor λ, possibilities of other separation regions between elastic layer and elastic half plane are determined. Also possibility of separation between first punch and elastic layer is researched. Besides, the distance (b-a)/h, that ends the interaction of punches is examined. It is assumed that Q/h≥P/h.Table 1
Variation of minimum value of distance between two punches(b-a)/h to avoid separation under first punch with a/h and (c-b)/h for various values of Q/P (μ2/μ1=6.48).
(
c
-
b
)
/
h
a
/
h
(
b
-
a
)
/
h
2
P
4
P
6
P
8
P
10
P
12
P
0.25
0.5
0.2773
0.5465
0.7263
0.8564
0.9569
1.0383
0.75
0.4645
0.7676
0.9485
1.075
1.1716
1.250
1.0
0.603
0.9107
1.089
1.2138
1.3105
1.3893
0.5
0.5
0.185
0.4403
0.617
0.7466
0.8474
0.9291
0.75
0.3617
0.6597
0.8405
0.9675
1.065
1.144
1.0
0.4983
0.8035
0.983
1.1086
1.206
1.2856
0.75
0.5
0.108
0.3432
0.516
0.645
0.746
0.8285
0.75
0.2715
0.5613
0.742
0.8699
0.9685
1.0485
1.0
0.4045
0.7073
0.8875
1.015
1.1134
1.4445
1.0
0.5
0.0357
0.2504
0.4195
0.548
0.6497
0.7332
0.75
0.191
0.4694
0.6502
0.7797
0.8798
0.9613
1.0
0.320
0.6188
0.8005
0.930
1.030
1.1123Table 2
Variation of distance between two punches(b-a)/h that ends interaction of punches with Q/P (μ2/μ1=2.75, a/h=0.25, and (c-b)/h=0.5).
Q
(
b
-
a
)
/
h
Punch I
Punch II
λ
c
r
right
(
x
c
r
right
-
a
)
/
h
λ
c
r
left
=
λ
c
r
right
P
10.5585
96.6748
2.091
96.6175
2.092
2
P
9.5874
96.0544
2.091
48.3612
2.092
4
P
8.6032
94.2026
2.094
24.1974
2.092
6
P
8.1313
91.6959
2.100
16.1360
2.092
8
P
7.8582
88.6769
2.111
12.1037
2.091
10
P
7.6742
85.2495
2.128
9.6838
2.091
12
P
7.5378
81.4924
2.153
8.0704
2.091
14
P
7.4308
81.4924
2.192
6.9178
2.090Table 3
Variation of distance between two punches(b-a)/h that ends interaction of punches with μ2/μ1 (Q=6P, a/h=0.25, and (c-b)/h=0.5).
μ
2
/
μ
1
(
b
-
a
)
/
h
Punch I
Punch II
λ
c
r
right
(
x
c
r
right
-
a
)
/
h
λ
c
r
left
=
λ
c
r
right
0.15
14.8404
164.5971
5.309
36.6977
4.795
0.36
12.1858
156.6651
3.678
30.6651
3.596
0.61
10.8897
143.5372
3.066
26.8566
3.027
1.65
9.2688
114.5372
2.418
20.5185
2.402
2.75
8.1313
91.6959
2.100
16.1360
2.092
6.48
6.4981
65.5235
1.812
11.4136
1.793Table 4
Variation of distance between two punches(b-a)/h that ends interaction of punches with a/h and (c-b)/h (Q=4P, μ2/μ1=0.36).
a
/
h
(
c
-
b
)
/
h
(
b
-
a
)
/
h
Punch I
Punch II
λ
c
r
right
(
x
c
r
right
-
a
)
/
h
λ
c
r
left
=
λ
c
r
right
(
b
-
x
c
r
left
)
/
h
=
(
x
c
r
right
-
c
)
/
h
0.25
1.5
10.8345
150.2590
3.775
54.1240
3.3816
0.25
1.0
11.8775
164.4214
3.655
49.0259
3.4495
0.25
0.5
12.7008
169.3000
3.628
45.9626
3.5968
0.50
1.0
11.6754
173.4645
3.523
49.0235
3.4494
0.75
1.5
10.4960
176.6011
3.421
54.1289
3.3810Table 5
Variation of load factor values with distance between two punches(b-a)/h (Q=2P, μ2/μ1=1.65, a/h=0.5, and (c-b)/h=1).
(
b
-
a
)
/
h
Punch I
Punch II
λ
c
r
right
x
c
r
right
λ
c
r
left
x
c
r
left
λ
c
r
right
x
c
r
right
0.5
30.6814
4.294
1
32.4279
4.806
3
34.0783
6.808
5
28.4331
3.167
28.4331
3.167
34.2756
8.803
6
32.0727
4.184
32.0727
4.184
34.3055
9.801
7
33.8531
5.211
33.8531
5.211
34.3175
10.800
7.9446
123.1714
2.835
34.3182
6.148
34.3182
11.743Figure 2
Variation of minimum value of distance between two punches(b-a)/h to avoid separation under first punch with Q/P for various values of μ2/μ1 (a/h=0.5, (c-b)/h=1).Figure 3
Contact stress distribution between elastic layer and elastic half plane for the cases of continuous (λ<λcr) and discontinuous contact (λ>λcr) (Q=2P, μ2/μ1=2.75, a/h=0.25, (b-a)/h=1.5, and (c-b)/h=0.5).Figure 4
Contact stress distribution between elastic layer and elastic half plane for the case of discontinuous contact (μ2/μ1=0.36, a/h=0.125, (b-a)/h=0.5, (c-b)/h=0.5, and λ=75>λcr).Figure 5
Contact stress distribution between elastic layer and elastic half plane for the case of discontinuous contact (Q=2P, a/h=0.5, (b-a)/h=2, (c-b)/h=1, and λ=100>λcr).Figure 6
Separation displacementv*(x,0) between elastic layer and elastic half space as a function of x for various values of second punch load Q(μ2/μ1=6.48, a/h=0.75, (b-a)/h=1, (c-b)/h=1.5, and λ=50>λcr), (*: [35]).Figure2 shows limit distance between punches that initiates separation under first punch for various values of material constant, μ2/μ1 with Q/P. The distance (b-a)/h increases, maintaining continuous contact between first punch and elastic layer with Q/P. Besides, for bigger values of μ2/μ1, limit distance that initiates separation under first punch decreases. In such a case, elastic half plane gets stiffer and it becomes easy to separate first punch from the elastic layer.Variation of critical distance between punches witha/h and (c-b)/h for various values of Q/P is presented in Table 1. For fixed values of second punch width, (c-b)/h, an increase in first punch width requires longer distance between punches to avoid separation under first punch. On the contrary, for fixed values of first punch width, a/h, an increase in second punch width decreases (b-a)/h and punches can be placed closer to each other without separation.Interaction between punches ends for a definite value of(b-a)/h. Tables 2–4 show the critical value of the distance that ends interaction between punches with dimensionless quantities μ2/μ1, a/h, (c-b)/h, and Q/P. In such a case, there is no need to consider punches together. Also these tables show the values of the load factor that cause separation between elastic layer and elastic half plane, λcr. For λ=λcr, σy1(x,0)h/P is zero. Contact between punches and elastic layer is continuous.Table3 shows the critical distance between punches that ends interaction of punches with elastic constant, μ2/μ1. For small values of μ2/μ1; that is, it is easy to bend elastic layer, interaction between punches ends in a longer distance. Initial separation point xcr between elastic layer and elastic half plane from the origin x=0 decreases with an increase in μ2/μ1. Critical load factor also decreases in this situation.Distance(b-a)/h that ends interaction between punches increases with a decrease in second punch width while first punch width is fixed. If both first and second widths are increased, interaction of punches ends in a shorter distance. This situation is presented in Table 4. In this case, critical load factor, λcr increases but initial separation point, xcr decreases with increment in punch widths. Separation also occurs at the right-hand side of the second punch.For fixed values ofa/h=0.5, (c-b)/h=1, Q=2P, and μ2/μ1=1.65, variations of critical load factor, λcr and initial separation point, xcr are given in Table 5. For small values of (b-a)/h, initial separation point between elastic layer and elastic half plane appears at the right-hand side of second punch. If (b-a)/h increases, in this case separation takes place between two punches. Keeping on increasing, the distance (b-a)/h ends the interaction of punches. For (b-a)/h=7.9446, there is no need to consider punches together.In Figures3–5, the normalized contact stress distribution σy1(x,0)h/P at the interface of elastic layer and elastic half plane is given for the problems described in Section 3, and Section 4. Different scales have been used for continuous and discontinuous contact cases in order to include the entire pressure distribution and to give sufficient details in compact forms.Variation of the contact stress distributionσy1(x,0)h/P with load factor λ=P/ρ1gh2 for fixed values of Q=2P, μ2/μ1=2.75, a/h=0.25, (b-a)/h=1.5, and (c-b)/h=0.5 is presented in Figure 3, λ<λcr and λ>λcr show contact stress distribution for the continuous and discontinuous contact cases, respectively. It can be seen from the figure that another separation region is possible if distance between punches (b-a)/h or load factor λ is increased. Contact pressure has peaks around the edges of the rigid punches.In Figure4, the variation of the normalized contact stress σy1(x,0)h/P with Q/P is given for the discontinuous contact case (λ=75>λcr). Figure 4 shows that both contact stress and separation zone (e-d)/h increase with an increase in Q/P. Variation of the normalized contact stress σy1(x,0)h/P for the discontinuous contact case shows three different regions between elastic layer and elastic half plane. These are continuous contact region, separation zone, and also continuous contact region, where the effects of external load Q/h and P/h decrease and disappear infinitely.Variation of contact stressσy1(x,0)h/P with μ2/μ1=1.65 is shown in Figure 5. Separation zone increases as elastic half plane gets stiffer than the elastic layer. In this case, peak value of the contact stress also increases.Results calculated from (33) giving the displacement v*(x,0) in the separation region d<x<e as a function of x are shown in Figure 6. Also results are compared with those of [35]. It is seen that separation region d<x<e and displacement values increase with an increase in Q/P ratio as this is the case in [35]. Separation zone is nearly the same, but displacement values are smaller than those of [35] for fixed values of μ2/μ1=6.48, a/h=0.75, (b-a)/h=1, and (c-b)/h=1.5.


## 6. Conclusion

In this paper, continuous and discontinuous contact problems for an elastic layer resting on an elastic half plane loaded by means of three rigid flat punches are considered. Numerical procedures developed in this study can be used to find approximate solutions to problems of engineering interest.Numerical results show that the punch widths, elastic constants, external loads, distance between punches play a very important role in the formation of the continuous and discontinuous contact areas, initial separation point and load, separation displacement, limit distance that ends the interaction of punches, critical distance that causes separation under first punch, and the contact pressure distribution.From the study, the following conclusions may be drawn:(i)
it is easy to separate first punch from elastic layer ifQ/P or a/h increases. This results also decrease in μ2/μ1 or (c-b)/h.(ii)
Increment inQ/P, a/h or (c-b)/h ends interaction of punches in a shorter distance. Also an increase in μ2/μ1 causes the same result.(iii)
First separation between elastic layer and elastic half plane occurs between punches or in the regionc<x<∞ depending on λ. If (b-a)/h is big enough, interaction of punches disappears.(iv)
Size of separation region is not affected much withQ/P, but separation displacement is increased with an increase in Q/P.



---

*Source: 137427-2013-04-22.xml*