Tensile Stresses in the Coating with Interlayer under Normal and Tangential Loading

In this paper, we obtained the analytic solution of the three-dimensional problem of elasticity concerning non-homogeneous half-space, in which the surface is under normal and tangential loading applied in a circular area. Half-space is composed of the homogeneous body, as well as coatings containing two homogeneous layers: interlayer and top layer. The analysis was carried out in cases where the Young’s modulus of the intermediate layer differed from the Young’s modulus of the substrate and the top layer. We have concentrated on first principal stress analysis. The obtained results can serve as an indication for the design of the composition of coatings in the context of tensile stress.


Introduction
Coatings improve a surface's tribological properties, particularly in terms of frictionally cooperative elements. Therefore, the appropriate selection of coatings may cause a decrease in the friction and reduction in the thermal or chemical adverse influence of an environment. Coatings are used when the mechanical contact of elements is unavoidable. Notable examples of coating usage are coatings of switches and plugs (in order to provide good electrical conductivity), tools, bearings, and optical laminas in surfaces subjected to mechanical contact.
In order to evaluate the mechanical properties of a coating, the stress field in the considered medium has to be determined. This stress field arises from contact pressure being applied to its surface. It is crucial to take into account the normal and tangential loading. The formation of cracks should be expected in places of maximum tensile stress. In homogeneous coatings, the maximum value of the σ 1 occurs on the surface in the neighbourhood of the edge of the contact area. This can be observed in articles presenting analytical [1,2] and numerical solutions [3,4] and both ways [5], as well as those describing numerical solutions confirmed by experiments [6] (similar to homogeneous half-space as can be seen in articles from several decades ago [7][8][9]) or the surroundings of the interface between the coating and the base [4,6,10].
Interlayers are used in order to increase the bonding force between the coating and base [3,[11][12][13][14]. The occurrence of these layers may also lead to a reduction in tensile stresses near the interface. The two-dimensional problem of the loading of inhomogeneous halfspace, composed of two homogeneous layers and a homogeneous body, was considered in [3]. Calculations were carried out using the finite element method. It was emphasised that the location of the maximum value of the characteristics of the stress tensor depends on the relations of Young's moduli of the components considered half-space, the ratios of the thickness of the components to the characteristic dimension of the contact area, and the coefficient of friction. In article [12], the loading of the coating with an interlayer was considered and the Huber von Mises stresses were calculated using the results of FEM calculations. On this basis, the authors concluded about the appearance of cracks. The axisymmetric problem of loading a body's surface with double-layer coating was considered in [13,14]. The shape of the indenter (its radius) and the thickness of the coating were analyzed in [15] to obtain the Huber von Mises stress distribution in the medium, and also thus for the moment of initiation of the destruction process.
The scientific research described in the literature was carried out using numerical methods, particularly the finite element method [3,4,6,12].
An analytical solution to the problem of punch movement on a multi-layer substrate was presented in [16]. The results concern the influence of punch speed, coefficient of friction, dimensions and force value on stress distribution. While designing the structure of coatings, a very important task is to ensure appropriate tribological properties. An equally important aspect is to determine the stress value that arises as a result of surface loading [6].
In the present work, we obtain an analytical solution of the three-dimensional problem of elasticity for half-space with a coating composed of two homogeneous layers. The external surface of the non-homogeneous half-space being considered is put under normal and tangential load, applied in a circular area. We focused on the first principal stress analysis.

Statement of the Problem
The subject matter of the considered problem is a homogeneous half-space with the coating of thickness, H. The half-space is composed of two homogeneous layers: the interlayer (thickness H 1 ) and the top layer (thickness H-H 1 ) ( Figure 1). The mechanical properties of the body are described by Young's modulus E and by Poisson's ratio, which is constant in each component of the considered inhomogeneous half-space. Between components of this medium, the ideal mechanical contact condition was satisfied. The surface of the considered half-space is under normal load and related with it by the Amontons-Coulomb's law tangential load: t = f p , where f -coefficient of friction. Loading is applied in a circular area Ω ( Figure 1) and has elliptical distribution.
where: x, y, z-Cartesian dimensionless coordinates related to radius a of loading circle.
area, and the coefficient of friction. In article [12], the loading of the coating with an inter layer was considered and the Huber von Mises stresses were calculated using the result of FEM calculations. On this basis, the authors concluded about the appearance of cracks The axisymmetric problem of loading a body's surface with double-layer coating was con sidered in [13,14]. The shape of the indenter (its radius) and the thickness of the coating were analyzed in [15] to obtain the Huber von Mises stress distribution in the medium and also thus for the moment of initiation of the destruction process. The scientific research described in the literature was carried out using numerica methods, particularly the finite element method [3,4,6,12].
An analytical solution to the problem of punch movement on a multi-layer substrat was presented in [16]. The results concern the influence of punch speed, coefficient of fric tion, dimensions and force value on stress distribution. While designing the structure o coatings, a very important task is to ensure appropriate tribological properties. An equally important aspect is to determine the stress value that arises as a result of surface loading [6].
In the present work, we obtain an analytical solution of the three-dimensional prob lem of elasticity for half-space with a coating composed of two homogeneous layers. Th external surface of the non-homogeneous half-space being considered is put under norma and tangential load, applied in a circular area. We focused on the first principal stres analysis.

Statement of the Problem
The subject matter of the considered problem is a homogeneous half-space with th coating of thickness, H. The half-space is composed of two homogeneous layers: the inter layer (thickness H1) and the top layer (thickness H-H1) ( Figure 1). The mechanical prop erties of the body are described by Young's modulus E and by Poisson's ratio, which i constant in each component of the considered inhomogeneous half-space. Between com ponents of this medium, the ideal mechanical contact condition was satisfied. The surfac of the considered half-space is under normal load and related with it by the Amontons Coulomb's law tangential load: t = fp, where f-coefficient of friction. Loading is applied in a circular area  ( Figure 1) and has elliptical distribution.

( )
where: x, y, z-Cartesian dimensionless coordinates related to radius a of loading circle.

Method of the Solution
Elasticity Equations (2)-(4) were solved using a two-dimensional Fourier integral transform.
Inserting Equations (12)-(19) into Equations (5)-(9), we obtain two independent systems of linear equations to determine functions a j (ξ, η) and b j (ξ, η). These systems contain 10 and 5 equations, respectively. After solving them and calculating the a j (ξ, η) and b j (ξ, η) functions, we use the inverse Fourier transform. As a result, we obtained stress tensor components in two-dimensional integral form. Integration was carried out using polar coordinates which were introduced in the integral transform plane. Integrals along the angular coordinate were calculated analytically. Integrals along the radial coordinate in internal points of considered half-space were calculated numerically using the Gaussian quadrature. When calculating the integrals which described stresses at the surface of considered inhomogeneous half-space (z = h), we consider the asymptotic behaviour of the integrand while s→∞.

Results Analysis
Analysis of the obtained relationship indicated that the analytical solution depends on eight dimensionless parameters: E 1 /E 0 , E 2 /E 0 , ν 0 , ν 1 , ν 2 , h, h 1 , f. Analysing interlayer influence, we consider cases in which: During the calculations we take, respectively: E1 = E 0 /2; E 1 = (E 0 + E 2 )/2; E 1 = 2E 2 . Participation of the interlayer's thickness in the thickness of the whole coating is in the range h 1 /h = 0.05 ÷ 0.25. The other parameters take the following values: E 2 /E 0 = 2 or E 2 /E 0 = 4; h = 0.4 or h = 0.8; f = 0; 0.1; 0.25; 0.5; ν 0 = ν 1 = ν 2 = 1/3. In the presented paper, we focus on the tensile stresses applied on the surface of inhomogeneous half-space and the interfaces between its components. Figures 2-5 shows the dimensionless parameter σ 1 /σ 1 S in the function of the relative thickness of interlayer h 1 /h, where: σ 1 -maximum tensile stress on the surface of considered half-space, and σ 1 S -maximum tensile stress in an adequate problem in which the coating does not contain an interlayer. The maximum increase in the tensile stress occurs when the Young modulus of the interlayer is less than the Young modulus of the base (E 1 /E 0 = 0.5). This stress increases simultaneously with the parameter h 1 /h. The interlayer described by the Young modulus E 1 = (E 0 + E 2 )/2 causes an insignificant increase in σ 1 stress on the surface. Only introducing the interlayer described by the Young modulus E 1 > E 2 causes a negligible decrease in the tensile stress on the surface. The difference between σ 1 and σ 1 S stresses decreases simultaneously when the coefficient of friction increases.

Results Analysis
Analysis of the obtained relationship indicated that the analytical solution depends on eight dimensionless parameters: E1/E0, E2/E0, 0, 1, 2, h, h1, f. Analysing interlayer influence, we consider cases in which: (1) E1 < E0; (2) E0 < E1 < E2; (3) E1 > E2. During the calculations we take, respectively: E1 = E0/2; E1 = (E0 + E2)/2; E1 = 2E2. Participation of the interlayer's thickness in the thickness of the whole coating is in the range h1/h = 0,05 ÷ 0,25. The other parameters take the following values: E2/E0 = 2 or E2/E0 = 4; h = 0,4 or h = 0,8; f = 0; 0,1; 0,25; 0,5; 0 = 1 = 2 = 1/3. In the presented paper, we focus on the tensile stresses applied on the surface of inhomogeneous half-space and the interfaces between its components. Figures 2-5 shows the dimensionless parameter 1/1 S in the function of the relative thickness of interlayer h1/h, where: 1-maximum tensile stress on the surface of considered half-space, and 1 S -maximum tensile stress in an adequate problem in which the coating does not contain an interlayer. The maximum increase in the tensile stress occurs when the Young modulus of the interlayer is less than the Young modulus of the base (E1/E0 = 0.5). This stress increases simultaneously with the parameter h1/h. The interlayer described by the Young modulus E1 = (E0 + E2)/2 causes an insignificant increase in 1 stress on the surface. Only introducing the interlayer described by the Young modulus E1 > E2 causes a negligible decrease in the tensile stress on the surface. The difference between 1 and 1 S stresses decreases simultaneously when the coefficient of friction increases.     Tensile stresses may occur in a medium with a homogeneous coating on the interlayer between the coating and base [4,5,10]. Figure 6 shows the correlation between maximum tensile stress on the surface and analogous stress on the interface between coating and base in a medium with a homogeneous coating in terms of thickness. Tensile stress on the interface takes significant values relative to the corresponding stresses on the surface at a relatively small coefficient of friction and parameter h > 0,5 ( Figure 6). That is why we analyse the influence of the interlayer to stress distribution near the interface for these parameters (Figure 7). This figure shows the ratio of the stresses 1/1 I in the function of the thickness of interlayer h1/h (1-maximum tensile stress on the interfaces between components of the considered medium with double-layered coatings, 1 I -maximum tensile stress on the interface between coating and base in the medium with the homogeneous coating ( Figure 6)). If the 1 stress take larger values on the interface between base and interlayer, then the graph is a continuous line; when larger values occur on the top layerinterlayer interface, the line is dashed. Introducing the interlayer, which is described by a Tensile stresses may occur in a medium with a homogeneous coating on the interlayer between the coating and base [4,5,10]. Figure 6 shows the correlation between maximum tensile stress on the surface and analogous stress on the interface between coating and base in a medium with a homogeneous coating in terms of thickness. Tensile stress on the interface takes significant values relative to the corresponding stresses on the surface at a relatively small coefficient of friction and parameter h > 0.5 ( Figure 6). That is why we analyse the influence of the interlayer to stress distribution near the interface for these parameters (Figure 7). This figure shows the ratio of the stresses σ 1 /σ 1 I in the function of the thickness of interlayer h 1 /h (σ 1 -maximum tensile stress on the interfaces between components of the considered medium with double-layered coatings, σ 1 I -maximum tensile stress on the interface between coating and base in the medium with the homogeneous coating ( Figure 6)). If the σ 1 stress take larger values on the interface between base and interlayer, then the graph is a continuous line; when larger values occur on the top layer-interlayer interface, the line is dashed. Introducing the interlayer, which is described by a Young modulus lower than base E 1 < E 0 (E 1 = E 0 /2) and larger than top layers modulus E 1 > E 2 (E 1 = 2E 2 ), results in the increase in stresses on the interface relative to corresponding stresses in the medium with homogeneous coating. A decrease in the stresses occurs when the Young modulus of the interlayer has a value between the top coat and base moduli:

Conclusions
In the presented work, the analytical solution to the three-dimensional problem of elasticity was considered. The research object is the half-space whose surface is under normal and tangential load. The half-space comprises a homogeneous body and coating, including the top layer and interlayer.

Conclusions
In the presented work, the analytical solution to the three-dimensional problem of elasticity was considered. The research object is the half-space whose surface is under normal and tangential load. The half-space comprises a homogeneous body and coating, including the top layer and interlayer.
Based on the stress analysis, introducing an interlayer with a Young modulus of E 1 < E 0 causes a significant increase in the σ 1 stress on the surface of the inhomogeneous half-space. When E 1 > E 2 , the tension stresses on the interface increase. Therefore, the choice of interlayer in which Young modulus is in range (E 0 , E 2 ) seems to be optimal. Introducing this layer may considerably reduce the tensile stresses near the interface between components of the considered medium. The increase in the stresses on the surface is insignificant. The test results can help in the selection of the composition of the cover in terms of minimizing the occurrence of tensile stress in it.