Skip to main content

ORIGINAL RESEARCH article

Front. Mech. Eng., 22 May 2020
Sec. Tribology
This article is part of the Research Topic Friction and Multi-Field Problems in Sliding Contacts View all 7 articles

Stress Distributions in the Cattaneo–Mindlin Problem on a Contact With Slip and Adhesion of Two Cylindrical Bodies

\nTaras Klimchuk
Taras Klimchuk1*Volodymyr OstrykVolodymyr Ostryk2
  • 1Department of Theoretical and Applied Mechanics, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
  • 2Institute of Applied Physics, National Academy of Sciences of Ukraine, Sumy, Ukraine

Two-dimensional Cattaneo–Mindlin problem on a contact with slip and adhesion between two elastic bodies of the same material is considered in classical statement for circular cylinders as well as in more general statements: for non-circular cylinders with higher order of touch and cylindrical bodies with a wavy periodic surface. Stress distributions are found both in the contact area and inside the bodies. The question of how the friction affects a formation of adhesion zone in the contact area and distribution of maximum tangential stress within the bodies is investigated.

Introduction

If elastic bodies are made of the same material and their sizes as well as the radii of curvature of their surfaces are much larger than the size of contact area, then a mutual slip of contacting surfaces will not occur in the case of normal contact. It means that friction forces do not appear in the contact area (Johnson, 1985). If a shear loading is applied to such contacted bodies, then depending on its magnitude, compared with a normal loading, two cases are possible. The first one is full mutual slip of surfaces of the bodies in whole contact area, which occurs under ultimate shear loading. In the second case, surfaces of the bodies clutch in the middle of contact area under the friction influence and slide against each other at its edges. In the last case, we have a contact with slip and adhesion. The size of adhesion zone is unknown in advance and depends, in particular, on a friction coefficient.

Contact with slip and adhesion between elastic bodies of the same material, first considered independently by Cattaneo (1938) and Mindlin (1949), is commonly called the Cattaneo–Mindlin problem (Johnson, 1985). In a two-dimensional case, this problem involves the contact of elastic cylindrical bodies under conditions of plane strain. In a three-dimensional case, it is about the contact of elastic spheres or bodies, which have smooth convex surfaces and touch each other previously at a point. The Cattaneo (1938) and Mindlin (1949) have found distribution of tangential stress in the contact area. The stress within the bodies, in particular, a parameter corresponding to the von Mises yield criterion was investigated in Dini and Hills (2004).

The generalization of the Cattaneo–Mindlin problem on the cases of higher-order touch, multiple contact areas and periodic contact are presented in Ciavarella (1998a,b), Jager (1998), Block and Keer (2008), Papangelo and Ciavarella (2015), and Papangelo et al. (2015). The stress distributions inside the contacting bodies were not studied in these papers. Periodic contact with slip and adhesion problems of plane strain between two elastic bodies with the same material (one is a half-space and another has a complex shape surface) were considered in Goryacheva et al. (2012), Chumak et al. (2014), Goryacheva and Martynyak (2014), and Malanchuk et al. (2017).

This article deals with two-dimensional Cattaneo–Mindlin problem as in classical statement for circular cylinders and also in two more general statements: for non-circular cylinders with higher order of touch and cylindrical bodies with a wavy periodic surface. In addition to the distribution of contact stresses, the distribution of stresses is found inside the bodies, in particular, the distribution of maximum tangential stress. Previously, a study of distribution of maximum tangential stress in elastic bodies, which are under conditions of smooth and sliding contact, was performed by Belyaev (1924), Sneddon (1946, 1948), Kuznetsov (1978), Kuznetsov and Gorokhovsky (1978), Johnson (1985), Ostryk (2015a,b), Klimchuk and Ostrik (2017), and Klimchuk and Ostryk (2018). Results of such studies enable to reveal zones with high values of tangential stresses, where a plastic flow of material appears according to the Tresca yield criterion.

Contact of Circular Cylinders

Let us consider two elastic cylinders, which have radii R1 and R2 as well as the same elastic constants (G is a shear modulus and ν a Poisson's ratio). Cylinders touch previously along the common generatrix (Figure 1A). They are being compressed by normal forces of intensity P uniformly distributed through the axial coordinate and form a contact area −axa (Figure 1B). After that, cylinders are applied by uniformly distributed tangential forces of intensity Q, perpendicular to generatrix (Figure 1C). Under the friction influence, the contact area divides into the adhesion zone cxc in the middle as well as slip zones −ax < c, c < xa at its edges. Boundary points of cylinders have the same tangential displacements in the adhesion zone. Tangential and normal stresses are subjected to the Amonton friction law in slip zones. Thus, they are proportional with a friction coefficient μ0. The plane strain independent of the axial coordinate takes place within the cylinders. Owing to the assumption that a contact area is small (a < < R1, R2), the problem is formulated for two elastic half-planes. Surfaces of the bodies near the contact area are approximated by parabolic cylinders y=x2/(2R1), y=-x2//(2R2).

FIGURE 1
www.frontiersin.org

Figure 1. Scheme of contact: (A) bodies prior to loading and (B) the first stage and (C) the second stage of loading.

On the first stage of normal loading, boundary conditions of the problem have the form

(uy(2)-uy(1))|y=0=x22R*,  σy(1)|y=0=σy(2)|y=0,-axa, σy(1)|y=0=σy(2)|y=0=0, |x|>a,τxy(1)|y=0=τxy(2)|y=0=0, -<x<, 1R*=1R1+1R2    (1)

where R* is a combined radius of curvature. According to the last condition, shear stresses are absent in the contact area due to the problem symmetry with respect to axis Oy. With the use of a solution of this problem (Johnson, 1985), the contact pressure function

p(x)=-12Gσy(1)|y=0=-12Gσy(2)|y=0, -axa    (2)

has a following form:

p(x)=Pπa2Ga2-x2.    (3)

A half-width of contact area is determined by an equality

a=2(1-ν)PR*πG,    (4)

which follows from an equilibrium condition

-aap(x)dx=P2G.    (5)

On the second stage of loading by tangential forces Q, considered in Cattaneo (1938) and Mindlin (1949), boundary conditions are the following:

(ux(2)-ux(1))|y=0=0, -cxc,  τxy(1)|y=0=μ0σy(1)|y=0,  c<|x|a,(uy(2)-uy(1))|y=0=x22R*,  σy(1)|y=0=σy(2)|y=0,  τxy(1)|y=0=τxy(2)|y=0,  -axa,σy(1)|y=0=σy(2)|y=0=τxy(1)|y=0=τxy(2)|y=0=0,  |x|>a.    (6)

Normal stress σy(1)|y=0 in the contact area remains the same as on the first stage. It is determined by equalities (2) and (3). According to the Cattaneo–Mindlin solution, a dimensionless function of tangential contact tractions

q(x)=-12Gτxy(1)|y=0=-12Gτxy(2)|y=0,  -axa    (7)

in slip and adhesion zones has the forms

q(x)=μ0Pπa2G(a2-x2-c2-x2), -cxc,q(x)=μ0p(x), c<|x|a.    (8)

The relative size of adhesion zone

ca=1-Qμ0P    (9)

is determined by an equilibrium condition

-aaq(x)dx=Q2G.    (10)

Figure 2 shows the distribution of normalized normal p̄(x)=(2Ga/P) and tangential q̄(x)=(2Ga/P)q(x) contact tractions in the case μ0 = 0.25 for different values of ratio Q/P. The distribution of contact pressure p̄(x) is permanent, and the distribution of tangential tractions q̄(x) has characteristic breaks at points x = ±c of transition from the adhesion zone to slip zones. For values Q/P = 0.05, 0.1, 0.15, and 0.2, the relative size c/a of the adhesion zone is 0.894, 0.775, 0.633, and 0.447, respectively. In the case of Q/P = 0, tangential tractions do not appear (q(x) ≡ 0) in the contact area (corresponding to a smooth contact). In the case of Q/P = μ0, we have a sliding contact. Herewith, the adhesion zone disappears (c/a = 0).

FIGURE 2
www.frontiersin.org

Figure 2. Distribution of normal p̄(x) and tangential q̄(x) contact tractions (μ0 = 0.25, ν = 0.3, Q/P = 0.05, 0.1, 0.15, 0.2, 0.25).

Stresses at each point (x, y) of elastic half-plane (cylinder 1) are found with use of the Kolosov formulas (Muskhelishvili, 1966):

12G(σx+σy)=-4Φ(z),  12G(σy-σx+2iτxy)                            =2[Φ¯(z)+Φ(z)+(z-z¯)Φ(z)],    (11)

where Φ(z) is the Kolosov–Muskhelishvili potential, an analytic function of complex variable z = x + iy, presented by the Cauchy type integral

Φ(z)=12πi-aap(s)-iq(s)s-zds.    (12)

According to formulas (3) and (8), density of integral (12) has the form

p(s)iq(s)=Pπa2G{ (1iμo)a2s2+iμ0c2s2,csc,(1iμo)a2s2,as<c,     c<sa.     (13)

Taking into account (Equation 13), we write the function Φ(z) from Equation (12) as follows:

Φ(z)=Pπa2G12πi [ (1iμ0)aaa2s2szds               +iμ0ccc2s2szds ].    (14)

With the use of a value of integral

-aaa2-s2s-zds=π(z2-a2-z),    (15)

obtained by the contour integration method (Muskhelishvili, 1966), finally, we find that

Φ(z)=-iP2πa2G[(1-μ0)z2-a2+iμ0z2-c2-z].    (16)

With the use of formula (11), the combinations of stresses are determined:

σx+σy=4Pπa2[ (1iμ0)z2a2+iμ0z2c2z ],σxσy+2iτxy=4Pπa2 [ μ0(z2a2z2c2) +y((1iμ0)zz2a2+iμ0zz2c21) ].    (17)

To find separate components of stress tensor from Equation (17), we set distances r1±, r2±, ρ and angles ϑ1±, ϑ2±, ϑ by the following dependences:

za=r1±eiϑ1±,   z±c=r2±eiϑ2±,   z=ρeiϑ.    (18)

Then we get

σx+σy=4Pπa2[r1+r1(sinϑ1++ϑ12μ0cosϑ1++ϑ12)+r2+r2μ0cosϑ2++ϑ22y],σxσy+2iτxy=4Pπa2[μ0r1+r1eiϑ1++ϑ12+μ0r2+r2eiϑ2++ϑ22+(1iμ0)ρyr1+r1ei(ϑϑ1++ϑ12)+iμ0ρyr2+r2ei(ϑϑ2++ϑ22)y].    (19)

Equality (19) yields the following stresses:

σx=2Pπa2{ r1+r1(sinϑ1++ϑ122μ0cosϑ1++ϑ12) +2μ0r2+r2cosϑ2++ϑ22+ρ2ar1+r1 [ cos(ϑϑ1++ϑ12) +μ0sin(ϑϑ1++ϑ12) ]sin(ϑ1+ϑ1)μ0ρ2cr2+r2sin(ϑϑ1++ϑ12)sin(ϑ2+ϑ2)1ar1+r1sin(ϑ1+ϑ1) },σy=2Pπa2{ r1+r1sinϑ1++ϑ12 +μ0ρ2cr2+r2sin(ϑϑ2++ϑ22)sin(ϑ2+ϑ2)ρ2ar1+r1 [ cos(ϑϑ1++ϑ12) +μ0sin(ϑϑ1++ϑ12) ]sin(ϑ1+ϑ1) },τxy=2Pπa2{ μ0r1+r1sinϑ1++ϑ12 +μ0r2+r2sinϑ2++ϑ22+ρ2ar1+r1 [ sin(ϑϑ1++ϑ12) μ0cos(ϑϑ1++ϑ12)]sin(ϑ1+ϑ1)+μ0ρ2cr2+r2cos(ϑϑ2++ϑ22)sin(ϑ2+ϑ2) }.    (20)

In addition, we find the principal normal σ1, σ2, σ3 and tangential τ1, τ2, τ3 stresses (Timoshenko and Goodyear, 1979)

σ1,2=12(σx+σy)±τ1,σ3=ν(σx+σy),τ1        =12|σy-σx+2iτxy|,τ2,3=12|σ2,1-σ3|    (21)

as well as maximum tangential stress

τmax=max(τ1,τ2,τ3).    (23)

Figures 3, 4 present the level lines of maximum tangential stress τmax, calculated for Poisson's ratio ν = 0.3. The solid lines correspond to τmax, and the dashed ones point to τ1, if τ1 < τmax = τ2. Figure 3A refers to a smooth contact (Q/P = 0), Figure 3B0 = 0.25, Q/P = 0.15), and Figure 4A0 = 0.5, Q/P = 0.3) refer to a contact with slip and adhesion, and Figure 3C0 = Q/P = 0.25) and Figure 4B0 = Q/P = 0.5) refer to a sliding contact. Distribution of stress τmax for a contact with slip and adhesion is transitional between the same distribution for the smooth and sliding contacts. In this case, the largest values of stress τmax increase, while the force ratio Q/P moves from 0 to μ0. If a friction coefficient μ = 0.25, then the value max τmax = max τ1 is reached at point x = ± x0, y = ± y0 inside each of contacting bodies. This point is shifting in a direction of shear loading and closer to the surface of the body in comparison with a smooth contact (Table 1). If a friction coefficient μ0 = 0.5, then with increase in Q/P, this point comes first to the surface of the body; and second, while moving to a sliding contact, the value τmax becomes the largest at all points of contact area (maxτmax=maxτ1=0.52Pπa, −1 < x/a < 1, Figure 4B). In addition, under sliding contact conditions for μ0 = 0.25 (as shown in Figure 3C), the largest value τmax=0.3312Pπa, slightly bigger than the value maxτ1=0.3302Pπa inside a half-plane, is reached within the slip zone at points x/a = ± 0.243 by stresses τ2 and τ3, respectively. For μ0 = 0.5 and Q/P = 0.3 (Figure 4A), the largest value τmax=τ1=0.3872Pπa is the same at all points of slip zones (0.633a < |x| < a).

FIGURE 3
www.frontiersin.org

Figure 3. Distribution of maximum tangential stress τmax in the cases of (A) smooth contact (μ0 = 0, Q/P = 0), (B) contact with slip and adhesion (μ0 = 0.25, Q/P = 0.15), and (C) sliding contact (μ0 = Q/P = 0.25).

FIGURE 4
www.frontiersin.org

Figure 4. Distribution of maximum tangential stress τmax in the cases of (A) contact with slip and adhesion (μ0 = 0.5, Q/P = 0.3) and (B) sliding contact (μ0 = Q/P = 0.5).

TABLE 1
www.frontiersin.org

Table 1. Points with the largest value of maximum tangential stress.

A more detailed numerical analysis shows that for friction coefficients from μ0 < 0.25 and for any ratio of forces in range 0 ≤ Q/P ≤ μ0, the value max τmax is reached at point inside a body. For each μ0 ≥ 0.25, there is such a transient value (Q/P)* that the value max τmax is realized inside a half-plane while 0 ≤ Q/P < (Q/P)* and at the boundary of a half-plane in slip zones while (Q/P)*Q/P ≤ μ0. Values (Q/P)* are given in Table 2 for various friction coefficients. The largest value τmax = τ2, 3 is reached at certain symmetrical points x = ±x1 inside the right and left slip zones, while 0.25 ≤ μ0 < 0.5. The largest value τmax = τ1 is constant along both slip zones while μ0 ≥ 0.5.

TABLE 2
www.frontiersin.org

Table 2. Transient values (Q/P)* for various friction coefficients μ0.

To summarize these results, we have obtained in analytical form the largest values τmax in slip zones. With the use of relations (17), (21), and (22), the following expressions are found for the principal tangential stresses:

τ1=μ02Pπa1-(ca)2,τ2,3     =Pπa2|(1-2ν)(a2-x2+sgnx·μ0x2-c2)     ±μ0a2-c2|,      y=0,c<|x|<a.    (23)

Here, we find that

maxc<|x|<aτmax|y=o={ τ2,3|y=0,x=±x1=Pπa2a2c2[ (12ν)1+μ02+μ0 ],μ0μ0*,τ1|y=0,c<| x |<a=μ02Pπa1(c/a)2,μ0μ0*,x1=μ02a2+c21+μ02μ0*=12ν2ν(1ν).    (24)

In particular, μ0*=0.436 for ν = 0.3. Values maxc<|x|<aτmax|y=0 in Equation (24) become the largest values τmax, if they are bigger than the values max τ1 inside a half-plane. This is illustrated in Figures 3C, 4A,B.

We note that in the case μ0 = 0.5, Q/P = 0.15, the distribution of maximum tangential stress is almost the same as for μ0 = 0.25, Q/P = 0.15 (shown in Figure 3B). This fact indicates that the distribution of maximum tangential stresses is not almost affected by the distribution of tangential tractions in the contact area, but only the main vector of tractions is significant.

Touch of Higher Order

Consider the Cattaneo–Mindlin problem for non-circular cylinders, if their prior contact along the common generatrix has higher order. In this case, we set surfaces of bodies near the contact area by equations y=B1|x|n, y=-B2|x|n (B1 > 0, B2 > 0; n = 2, 3, …). If n = 2, then we have the classic Cattaneo–Mindlin problem for circular cylinders.

Further, we use singular integral relations (Johnson, 1985):

ux(1,2)x|y=0=-(1-2ν)p(x)±2(1-ν)1πa-aq(s)s-xds,uy(1,2)x|y=0=(1-2ν)q(x)±2(1-ν)1πa-ap(s)s-xds,    (25)

where functions p(x) and q(x) are determined by Equations (2) and (7), respectively, while −axa and equal to zero at |x| > a.

The boundary conditions for the first stage of normal loading are set in Equation (1), except the first condition, which should be written as follows:

(uy(2)-uy(1))|y=0=B|x|n,-axa,B=B1+B2.    (26)

Satisfying differentiated condition (26) by the second relation (25) with q(x) ≡ 0, we obtain a singular integral equation

1πa-ap(s)s-xds=-nBsgnx4(1-ν)|x|n-1,-a<x<a.    (27)

With the solution of Equation (27), bounded at the edges of interval −a < x < a, we find in the form (Gakhov, 1963)

p(x)=nB4(1-ν)πa2-x2In(x),In(x)=a-a|s|n-1sgnsa2-s2dss-x.    (28)

The recurrent relation holds for an integral in Equation (28),

In(x)=x2In-2(x)+2Jn-2,n=3,4,,Jn=0asna2-s2ds=n-1na2Jn-2=anJn,J2m=(2m-1)!!(2m)!!π2,J2m-1=(2m-2)!!(2m-1)!!,m=1,2,,    (29)

with initial values

I1(x)=2a2-x21na+a2-x2|x|,I2(x)=π.

Applying a recurrent formula (29) several times, we find

I2m(x)=2(J2m-2+J2m-4x2++J2xn-4)+I2(x)x2(m-1),I2m-1(x)=2(J2m-3+J2m-5x2++J1x2m-4)+I1(x)x2(m-1).

Therefore, on the basis of Equation (28), we obtain the distribution of contact pressure

p(x)=P2πa2Gp~(x,a),p~(x,a)=1-(-1)n2n!!(n-1)!!xn-1an-21na+a2-x2|x|+a2-x2k=1[n/2][n/2]k[(n-1)/2]k(xa)2k-2,-axa,    (30)

where [n/2] is an entire part of n/2, [n]k = n(n − 1)…(nk + 1). Dependence of half-width a of the contact area on the force P

a=(1-ν)PnJnBGn    (31)

is found from an equilibrium condition (5).

In the case of even n, formulas (30) and (31) have been obtained by Shtayerman (1949).

Figure 5 shows the dimensionless contact pressure p̄(x)=(2Ga/P)p(x) distribution on half-width of contact area (0 ≤ xa) for different values of indicator n. It can be noticed that with increase in n, the contact pressure decreases inside the contact area but increases at its edge.

FIGURE 5
www.frontiersin.org

Figure 5. Distribution of contact pressure p̄(x) for different n.

On the second stage, elastic bodies are loaded by tangential forces Q. Contact area −axa divides into the adhesion zone −cxc and the slip zones −ax < −c, c < xa. In boundary condition (6), the third condition should be replaced by Equation (26) and the normal stress σy(1)|y=0 in the contact area should take the form of Equations (2) and (30). Satisfying differentiated first condition (6) by the first relation (25), we obtain an integral equation with respect to the function of tangential contact tractions q(x) in the adhesion zone

-aaq(s)s-xds=0,-c<x<c.    (32)

Because the normal and tangential tractions in slip zones are connected by equality

q(x)=μ0p(x),c<|x|a,    (33)

Equation (32) is written as follows:

-ccq(s)s-xds=-μ0(ca+-a-c)p(s)s-xds,-c<x<c.    (34)

The substitution of unknown function

q0(x)=q(x)-μ0p(x),-c<x<c    (35)

yields the following equation:

-ccq0(s)s-xds=-μ0-aap(s)s-xds,-c<x<c.    (36)

Taking into account (Equation 27), integral Equation (36) takes the form

1π-ccq0(s)s-xds=μ0nBsgnx4(1-ν)|x|n-1,-c<x<c    (37)

and differs from the integral Equation (27) in both an integration interval (a is replaced by c) and the multiplier −μ0 at the right part. Therefore, with respect to Equation (30), we obtain its bounded solution in the form

q0(x)=-μ0P2πa2G(ca)n-2p~(x,c).    (38)

With the use of Equations (35) and (38), the function of tangential contact tractions is found in the adhesion zone

q(x)=q0(x)+μ0p(x)=μ0P2πa2G[p~(x,a)-(ca)n-2p~(x,c)].    (39)

With the use of an equilibrium condition (10), the relative size of adhesion zone is determined

ca=1-Qμ0Pn.    (40)

Stresses in half-plane y ≥ 0 are found with use of formula (11), in which the potential density Φ(z) from Equation (12) is given as follows:

p(s)​​​iq(s) ​= ​P2πa2G{ (1iμ0)p˜(s,a)+iμ0(ca)n2p˜(s,c),csc,(1iμ0)p˜(s,a),as<c,c<sa.    (41)

Taking into account (Equation 41), the function Φ(z) from Equation (12) is presented as

Φ(z)=P2πa2G[(1-iμ0)Φ~(z,a)+iμ0(ca)n-2Φ~(z,c)],                                  Φ~(z,a)=12πi-aap~(s,a)s-zds.    (42)

Further, to find the function Φ~(z,a) in Equation (42), we consider two cases.

1. The case of even n. Using the second formula of Equation (30), we have that

Φ˜(z,a)=12πi( nn1U0(z)+n(n2)(n1)(n3)U1(z)a2+ ...+n!!(n1)!!Un/21(z)an2 ),Uk(z)=aas2ka2s2szds,    k=0,1, ... .     (43)

The following recurrent formula is true for integral Uk(z) in Equation (43):

Uk(z)=z2Uk-1(z)+2zVk-1,   k=1,2,,Vk=0as2ka2-s2ds=J2k+22k+1,   k=0,1,.    (44)

According to Equation (15), the integral U0(z) in Equation (43), which is an initial value for formula (44), can be presented as

U0(z)=π(z2-a2-z).    (45)

Consistently applying recurrent formula (44) and taking into account Equation (45), we find

Uk(z)=πz [ (2k3)!!(2k)!! a2k+(2k5)!!(2k2)!!a2k2z2+...+(1)!!2!!a2z2k2+z2k1(z2a2z) ],  k=0,1,... .    (46)

Substituting Uk(z)(k = 0, 1, .. . , n/2-1) from Equation (46) into the first formula of Equation (43), after transformations, we get

   Φ~(z,a)=-i2k=1n\2(ankz2-a2+cnkz)(za)2k-2,ank=[n/2]k[(n-1)/2]k, cnk=m=kn/2(2m-2k-3)!!(2m-2k)!!anm,    (47)

assuming that (−)!! = 1, (−3)!! = −1, and

  Φ˜z(z,a)=i2k=1n/2(bnkzz2a2(2k1)cnk)(za)2k2,bnk=[n/2]k[(n3)/2k, n=2,4,...; k=1,2,...,n/2.    (48)

In particular, for the case of n = 4 we have

a41=43, a42=83, b42=4, b41=-8, c41=0, c42=-83,  Φ~(z,a)=-2i3[(1+2z2a2)z2-a2-2z3a2].    (49)

On the basis of Equations (42), (47), and (48), we find the stresses with use of formula (11).

1. The case of odd n. Similar to the case of even n, we obtain

   Φ˜(z,a)=ia2 [ n!!(n1)!!( 1πk=1n1Jnk1nk(za)k1   + i(za)n1lnia+z2a2z )   +1az2a2k=1(n1)/2(ank+cnkza)(za)2k2 ],  cnk=m=k(n1)/2(2m2k3)!!(2m2k)!!anm,Φ˜z(z,a)=i2 [ n!!(n1)!!1πk=2n1k1nkJnk1(za)k2 az2a2k=1(n3)/2bnk(za)2k1  +n!!(n3)!!(ilnia+z2a2z+az2a2)(za)n2+k=1(n1)/2(2k1)cnk(za)2k2 ],J2k=2J2k, J2k1=0.    (50)

values τ*=πa2Pmaxτmax=πa2Pmaxτ1 calculated for ν = 0.3 in the case of smooth contact (Q/P = 0) are given in Table 3 as well as relative coordinates ±x¯0=±x0/a,y¯0=y0/a of points, in which these values are reached. With increase in n, these points approach the edges of contact area.

TABLE 3
www.frontiersin.org

Table 3. Points with values of τ* for different n.

The distribution of maximum tangential stress τmax in the case n = 4, ν = 0.3 is given in Figure 6. Figure 6A refers to a smooth contact (Q/P = 0), Figure 6B to a contact with slip and adhesion (μ0 = 0.25, Q/P = 0.15), and Figure 6C to a sliding contact (μ0 = Q/P = 0.25). As in Contact of Circular Cylinders, solid lines correspond to τmax and dashed ones to τ1, if τ1 < τmax = τ2. By comparing with a second order touch (n = 2, Figure 3), the stress concentration is observed closer to the edges of contact area.

FIGURE 6
www.frontiersin.org

Figure 6. Distribution of maximum tangential stress τmax in the cases of (A) smooth contact (μ0 = 0, Q/P = 0), (B) contact with slip and adhesion (μ0 = 0.25, Q/P = 0.15), and (C) sliding contact (μ0 = Q/P = 0.25) for n = 4.

Periodic Case

Assume that surfaces of elastic cylindrical bodies are l-periodic along the axis Ox and touch previously through the lines x = kl(k = 0, ±1, ±2, …). Surfaces are set mathematically by equations y = -B1(cos2πxl-1), y = B2(cos2πxl-1) (B1 > 0, B2 > 0). In this variant of the Cattaneo–Mindlin problem, the bodies are compressed at infinity by the normal forces of intensity p, creating the contact areas −a + klxa + kl. After this, the tangential forces of intensity q are applied at infinity and contact areas divide into adhesion −c + klxc + kl and slip c + kl < |x| ≤ a + kl(k = 0, ±1, ±2, …) zones.

Singular integral relation (25) in the periodic case takes the form

ux(1,2)x|y=0=-(1-2v)p(x)±2(1-v)1l-aaq(s)cotπl(s-x)ds,uy(1,2)x|y=0=-(1-2v)q(x)±2(1-v)1l-aap(s)cotπl(s-x)ds,    (51)

where functions p(x), q(x) are l-periodic and defined by equalities (2) and (7).

On the first stage of normal loading, boundary conditions are as follows:

(uy(2)-uy(1))|y=0=B(cos2πxl+1), σy(1)|y=0=σy(2)|y=0,-a+klxa+kl,σy(1)|y=0=σy(2)|y=0=0, a+kl<x<l-a+kl,τxy(1)|y=0=τxy(2)|y=0=0, -<x<,B=B1+B2,  k=0,±1,  ±2,.    (52)

Satisfying differentiated first condition (52) by the second relation (51) with q(x) ≡ 0, we obtain a singular integral equation.

-aap(s)cotπl(s-x)ds=-πB1-vsin2πxl,  -a<x<a.    (53)

The solution of Equation (53), bounded at the ends of interval −a < x < a, has the form (Shtayerman, 1949)

p(x)=2πB(1-v)lcosπxl.sinπl(a-x)sinπl(a+x),         -axa.    (54)

With the use of an equilibrium condition,

-aap(x)dx=pl2G,    (55)

a half-width of each contact area is determined

a=lπarcsin(1-v)plπBG.    (56)

The boundary conditions on the second stage of tangential loading are as follows:

(ux(2)-ux(1))|y=0=0,  -c+klxc+kl,   τxy(1)|y=0                                         =μ0σy(1)|y=0, c+kl<|x|a+kl,(uy(2)-uy(1))|y=0=B(cos2πxl+1),  σy(1)|y=0                                         =σy(2)|y=0,  τxy(1)|y=0=τxy(2)|y=0,                                         -a+klxa+kl,                  σy(1)|y=0=σy(2)|y=0=τxy(1)|y=0=τxy(2)|y=0=0,a                                        +kl<x<l-a+kl,  k=0,  ±1,±2,.    (57)

Here, the normal stress σy(1)|y=0 is known and presented by formulas (2) and (54).

Using the first boundary condition (57) and the first relation (51), in similar way to Touch of Higher Order, we obtain the integral equation with respect to l-periodic function q0(x) from Equation (35)

-ccq0(s)cotπl(s-x)ds=μ0πB1-vsin2πxl,   -c<x<c.    (58)

The last one is similar to Equation (53). Its bounded solution has the form

q0(x)=-μ02πB(1-v)lcosπxl.sinπl(c-x)sinπl(c+x),           -cxc.    (59)

From Equations (35) and (59), we find l-periodic function of tangential contact tractions

           q(x)=μ0p(x),  c<| x |a,q(x)=μ02πB(1v)lcosπxl.( sinπl(cx)sinπl(c+x) sinπl(ax)sinπl(a+x) ),cxc.     (60)

From an equilibrium condition,

-aaq(x)=ql2G    (61)

the relative width of adhesion zones is determined as

ca=1πarcsin(1-Qμ0Psinπal).    (62)

The stresses in half-plane y ≥ 0 are found with use of formula (11), where the potential Φ(z) has the form

    Φ(z)=12ilaa[ p(s)iq(s) ]cotπl(sz)ds,p(s)iq(s)={     (1iμ0)p(s)iq0(s),csc,(1iμ0)p(s),as<c,c<sa,    (63)

or

Φ(z)=1-iμ02il-aap(s)cotπl(s-z)ds-12l-ccq0(s)cotπl(s-z)ds.    (64)

Calculating the integral

-aacosπsl.sinπl(a-s)sinπl(a+s)cotπl(s-z)ds=lcosπzl.(sinπl(z-a)sinπl(z+a)-sinπzl)    (65)

and taking into account (Equation 56), we obtain

Φ(z)=ipcos(πz/l)Gsin2(πa/l) [ (1iμ0)sinπl(za)sinπl(z+a) +iμ0sinπl(zc)sinπl(z+c)sinπzl ],Φ¯(z)+Φ(z)=2μ0pcos(πz/l)Gsin2(πa/l)( sinπl(za)sinπl(z+a) sinπl(zc)sinπl(z+c) ),Φ(z)=iπpcos(πz/l)Glsin2(πa/l) [ (1iμ0)[cos(2πz/l)+sin2(πa/l)]sin(πz/l)sin[(π/l)(za)]sin[(π/l)(z+a)] +iμ0[cos(2πz/l)+sin2(πc/l)]sin(πz/l)sin[(π/l)(zc)]sin[(π/l)(z+c)] ].    (66)

Substituting expressions from Equation (66) into Equation (11), the stresses are determined.

In Figure 7, the distribution of maximum tangential stress τmax (ν = 0.3, l/a = 8) is shown in cases (A) μ0 = 0.25, Q/P = 0.15; (B) μ0 = 0.5, Q/P = 0.3. The solid lines correspond to τmax, and the dashed and dash-dotted lines correspond to τ1, if τ1 < τmax = τ2 and τ1 < τmax = τ3, respectively. The maximum tangential stress τmax = τ1 reaches the largest value within a half-plane [case (A)] or at half-plane boundary in points x ± c[c/a = 0.623 case (B)]. In the last case, the maximum tangential stress slightly decreases (within 0.6%) along the slip zones. The presented distribution differs a little from the stress distribution in the corresponding cases for a single contact area (l/a = ∞, Contact of Circular Cylinders, Figures 3B, 4A) and coincides with them at l/a ≥ 20.

FIGURE 7
www.frontiersin.org

Figure 7. Distribution of maximum tangential stress τmax in the periodic case (l/a = 8) for (A) μ0 = 0.25, Q/P = 0.15 and (B) μ0 = 0.5, Q/P = 0.3.

Conclusions

With the use of numerical analysis of obtained stress distributions in the Cattaneo–Mindlin problem, it is established that the point, in which maximum tangential stress reaches its largest value, is located inside each of contacting elastic bodies for a friction coefficient μ < 0.25. With increase in a friction coefficient and a shear loading, this point moves gradually to the boundary of the body (in the contact area). In the case of a touch of higher order, this point moves closer to the far edge (x = a) of contact area. At the same time, the local increase of stress occurs near the edge x = −a. Comparing with the case of a single contact area, in the periodic problem, the distribution of maximum tangential stress is almost unchanged, if a contact area does not exceed a quarter of period (l/a ≥ 8).

Data Availability Statement

The datasets generated for this study are available on request to the corresponding author.

Author Contributions

VO set a problem and performed the analytical calculations. TK created the code and performed the numerical calculations. VO and TK analyzed the results.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Belyaev, M. M. (1924). “Local stresses during compression of elastic bodies,” in Engineering Structures and Structural Mechanics (Saint Petersburg: The Way), 30–43.

Google Scholar

Block, J. M., and Keer, L. M. (2008). Periodic contact problems in plane elasticity. J. Mech. Mater. Struct. 3, 1207–1237. doi: 10.2140/jomms.2008.3.1207

CrossRef Full Text | Google Scholar

Cattaneo, C. (1938). Sul contatto di due corpi elastici: distribuzione locale degli sforzi. Rendiconti dell' Accademia nazionale dei Lincei. 27, 474–478.

Google Scholar

Chumak, K., Malanchuk, N., and Martynyak, R. (2014). Partial slip contact problem for solids with regular surface texture assuming thermal insulation or thermal permeability of interface gaps. Int. J. Mech. Sci. 84, 138–146. doi: 10.1016/j.ijmecsci.2014.04.015

CrossRef Full Text | Google Scholar

Ciavarella, M. (1998a). The generalized Cattaneo partial slip plane contact problem. I – theory. Int. J. Solids Struct. 35, 2349–2362. doi: 10.1016/S0020-7683(97)00154-6

CrossRef Full Text | Google Scholar

Ciavarella, M. (1998b). The generalized Cattaneo partial slip plane contact problem. II – examples. Int. J. Solids Struct. 35, 2363–2378. doi: 10.1016/S0020-7683(97)00155-8

CrossRef Full Text | Google Scholar

Dini, D., and Hills, D. A. (2004). Bounded asymptotic solutions for incomplete contacts in partial slip. Int. J. Solids Struct. 41, 7049–7062. doi: 10.1016/j.ijsolstr.2004.05.058

CrossRef Full Text | Google Scholar

Gakhov, F. D. (1963). Boundary Value Problems. Moscow: Fizmatgiz.

Google Scholar

Goryacheva, I. G., Malanchuk, N. I., and Martynyak, R. M. (2012). Contact interaction of bodies with a periodic relief during partial slip. J. Appl. Math. Mech. 76, 695–709. doi: 10.1016/j.jappmathmech.2013.03.002

CrossRef Full Text | Google Scholar

Goryacheva, I. G., and Martynyak, R. M. (2014). Contact problems for textured surfaces involving frictional effects. Proc. Inst. Mech. Eng. Part J. 228, 707–716. doi: 10.1177/1350650114528318

CrossRef Full Text | Google Scholar

Jager, J. (1998). A new principle in contact mechanics. ASME J. Tribol. 120, 677–684. doi: 10.1115/1.2833765

CrossRef Full Text | Google Scholar

Johnson, K. L. (1985). Contact Mechanics. Cambridge: Cambridge University Press. doi: 10.1017/CBO9781139171731

CrossRef Full Text | Google Scholar

Klimchuk, T. V., and Ostrik, V. I. (2017). Frictional contact between an elastic strip and a semi-infinite punch with rounded edge. Acta Mechanica 228, 3619–3631. doi: 10.1007/s00707-017-1866-8

CrossRef Full Text | Google Scholar

Klimchuk, T. V., and Ostryk, V. I. (2018). Smooth contact of a semi-infinite punch with rounded edge and an elastic strip. J. Math. Sci. 231, 650–664. doi: 10.1007/s10958-018-3842-9

CrossRef Full Text | Google Scholar

Kuznetsov, Y.A., and Gorokhovsky, G. A. (1978). The effect of roughness on the stress state of frictionally contacting bodies. Appl. Mech. 14, 62–68. doi: 10.1007/BF00885748

CrossRef Full Text | Google Scholar

Kuznetsov, Y. A. (1978). The superposition principle in the solution of contact – friction stress problems. Wear 50, 183–189. doi: 10.1016/0043-1648(78)90255-7

CrossRef Full Text | Google Scholar

Malanchuk, N., Slobodyan, B., and Martynyak, R. (2017). Frictional contact of two solids with a periodically grooved surface in the presence of an ideal gas in interface gaps. J. Theor. Appl. Mech. 55, 1181–1192. doi: 10.15632/jtam-pl.55.4.1181

CrossRef Full Text | Google Scholar

Mindlin, R. D. (1949). Compliance of elastic bodies in contact. J. Appl. Mech. 16, 259–268.

Google Scholar

Muskhelishvili, N. I. (1966). Some Basic Problems of the Mathematical Theory of Elasticity. Moscow: Nauka.

Google Scholar

Ostryk, V. I. (2015a). Contact Mechanics. Kyiv: Kyiv University Press.

Google Scholar

Ostryk, V. I. (2015b). Smooth contact of a polynomial profile stamp with an elastic half-plane. Bull. Kyiv Univ. 1, 41–44.

Google Scholar

Papangelo, A., and Ciavarella, M. (2015). Cattaneo – Mindlin plane problem with Griffith friction. Wear 342, 398–407. doi: 10.1016/j.wear.2015.10.005

CrossRef Full Text | Google Scholar

Papangelo, A., Ciavarella, M., and Barber, J. R. (2015). Fracture mechanics implications for apparent static friction coefficient in contact problems involving slip-weakening laws. Proc. R Soc. A 471:271. doi: 10.1098/rspa.2015.0271

CrossRef Full Text | Google Scholar

Shtayerman, I. Y. (1949). Contact Problem of the Theory of Elasticity. Moscow: Gostekhizdat.

Google Scholar

Sneddon, I. N. (1946). Boussinesq's problem for a flat-ended cylinder. Proc. Camb. Philos. Soc. 42, 29–39. doi: 10.1017/S0305004100022702

CrossRef Full Text | Google Scholar

Sneddon, I. N. (1948). Boussinesq's problem for a rigid cone. Proc. Camb. Philos. Soc. 44, 492–507. doi: 10.1017/S0305004100024518

CrossRef Full Text | Google Scholar

Timoshenko, S. P., and Goodyear, J. N. (1979). Theory of Elasticity. Moscow: Nauka.

Google Scholar

Keywords: smooth contact, sliding contact, slip and adhesion, plain strain, maximum tangential stress, Cattaneo-Mindlin problem, periodic problem, cylindrical bodies

Citation: Klimchuk T and Ostryk V (2020) Stress Distributions in the Cattaneo–Mindlin Problem on a Contact With Slip and Adhesion of Two Cylindrical Bodies. Front. Mech. Eng. 6:22. doi: 10.3389/fmech.2020.00022

Received: 12 January 2020; Accepted: 09 April 2020;
Published: 22 May 2020.

Edited by:

Elena Torskaya, Institute for Problems in Mechanics (RAS), Russia

Reviewed by:

Sergey Lychev, Institute for Problems in Mechanics (RAS), Russia
Antonio Papangelo, Politecnico di Bari, Italy

Copyright © 2020 Klimchuk and Ostryk. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Taras Klimchuk, tarasyniv@ukr.net

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.