ORIGINAL RESEARCH article
Volume 6 - 2020 | https://doi.org/10.3389/fmech.2020.00022
Stress Distributions in the Cattaneo–Mindlin Problem on a Contact With Slip and Adhesion of Two Cylindrical Bodies
- 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.
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 −a ≤ x ≤ a (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 c ≤ x ≤ c in the middle as well as slip zones −a ≤ x < c, c < x ≤ a 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 , .
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
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
has a following form:
A half-width of contact area is determined by an equality
which follows from an equilibrium condition
Normal stress 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
in slip and adhesion zones has the forms
The relative size of adhesion zone
is determined by an equilibrium condition
Figure 2 shows the distribution of normalized normal and tangential contact tractions in the case μ0 = 0.25 for different values of ratio Q/P. The distribution of contact pressure is permanent, and the distribution of tangential tractions 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. Distribution of normal and tangential 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):
where Φ(z) is the Kolosov–Muskhelishvili potential, an analytic function of complex variable z = x + iy, presented by the Cauchy type integral
According to formulas (3) and (8), density of integral (12) has the form
Taking into account (Equation 13), we write the function Φ(z) from Equation (12) as follows:
With the use of a value of integral
obtained by the contour integration method (Muskhelishvili, 1966), finally, we find that
With the use of formula (11), the combinations of stresses are determined:
To find separate components of stress tensor from Equation (17), we set distances , , ρ and angles , , ϑ by the following dependences:
Then we get
Equality (19) yields the following stresses:
In addition, we find the principal normal σ1, σ2, σ3 and tangential τ1, τ2, τ3 stresses (Timoshenko and Goodyear, 1979)
as well as maximum tangential stress
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 3B (μ0 = 0.25, Q/P = 0.15), and Figure 4A (μ0 = 0.5, Q/P = 0.3) refer to a contact with slip and adhesion, and Figure 3C (μ0 = Q/P = 0.25) and Figure 4B (μ0 = 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 (, −1 < x/a < 1, Figure 4B). In addition, under sliding contact conditions for μ0 = 0.25 (as shown in Figure 3C), the largest value , slightly bigger than the value 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 is the same at all points of slip zones (0.633a < |x| < a).
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. 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).
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.
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:
Here, we find that
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 , (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):
where functions p(x) and q(x) are determined by Equations (2) and (7), respectively, while −a ≤ x ≤ a 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:
Satisfying differentiated condition (26) by the second relation (25) with q(x) ≡ 0, we obtain a singular integral equation
With the solution of Equation (27), bounded at the edges of interval −a < x < a, we find in the form (Gakhov, 1963)
The recurrent relation holds for an integral in Equation (28),
with initial values
Applying a recurrent formula (29) several times, we find
Therefore, on the basis of Equation (28), we obtain the distribution of contact pressure
where [n/2] is an entire part of n/2, [n]k = n(n − 1)…(n − k + 1). Dependence of half-width a of the contact area on the force P
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 distribution on half-width of contact area (0 ≤ x ≤ a) 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.
On the second stage, elastic bodies are loaded by tangential forces Q. Contact area −a ≤ x ≤ a divides into the adhesion zone −c ≤ x ≤ c and the slip zones −a ≤ x < −c, c < x ≤ a. In boundary condition (6), the third condition should be replaced by Equation (26) and the normal stress 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
Because the normal and tangential tractions in slip zones are connected by equality
Equation (32) is written as follows:
The substitution of unknown function
yields the following equation:
Taking into account (Equation 27), integral Equation (36) takes the form
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
With the use of Equations (35) and (38), the function of tangential contact tractions is found in the adhesion zone
With the use of an equilibrium condition (10), the relative size of adhesion zone is determined
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:
Taking into account (Equation 41), the function Φ(z) from Equation (12) is presented as
Further, to find the function in Equation (42), we consider two cases.
1. The case of even n. Using the second formula of Equation (30), we have that
The following recurrent formula is true for integral Uk(z) in Equation (43):
According to Equation (15), the integral U0(z) in Equation (43), which is an initial value for formula (44), can be presented as
Consistently applying recurrent formula (44) and taking into account Equation (45), we find
Substituting Uk(z)(k = 0, 1, .. . , n/2-1) from Equation (46) into the first formula of Equation (43), after transformations, we get
assuming that (−)!! = 1, (−3)!! = −1, and
In particular, for the case of n = 4 we have
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
values calculated for ν = 0.3 in the case of smooth contact (Q/P = 0) are given in Table 3 as well as relative coordinates of points, in which these values are reached. With increase in n, these points approach the edges of contact area.
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. 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.
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 = y = (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 + kl ≤ x ≤ a + kl. After this, the tangential forces of intensity q∞ are applied at infinity and contact areas divide into adhesion −c + kl ≤ x ≤ c + kl and slip c + kl < |x| ≤ a + kl(k = 0, ±1, ±2, …) zones.
Singular integral relation (25) in the periodic case takes the form
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:
Satisfying differentiated first condition (52) by the second relation (51) with q(x) ≡ 0, we obtain a singular integral equation.
The solution of Equation (53), bounded at the ends of interval −a < x < a, has the form (Shtayerman, 1949)
With the use of an equilibrium condition,
a half-width of each contact area is determined
The boundary conditions on the second stage of tangential loading are as follows:
Here, the normal stress 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)
The last one is similar to Equation (53). Its bounded solution has the form
From Equations (35) and (59), we find l-periodic function of tangential contact tractions
From an equilibrium condition,
the relative width of adhesion zones is determined as
The stresses in half-plane y ≥ 0 are found with use of formula (11), where the potential Φ(z) has the form
Calculating the integral
and taking into account (Equation 56), we obtain
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. 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.
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.
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.
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
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
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
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
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, firstname.lastname@example.org