# Elastic Analysis for Nanocontact Problem with Surface Stress Effects under Shear Load

Authors: D. X. Lei; L. Y. Wang; Z. Y. Ou
Journal: Journal of Nanomaterials
 (2012)

---

## Abstract

Consideration of surface stress effects on the elastic field of nanocontact problem has extensive applications in several modern problems of solid mechanics. In this paper, the effects of surface stress on the contact problem at nanometers are studied in the frame of surface elasticity theory. Fourier integral transform method is adopted to derive the fundamental solution of the nanocontact problem under shear load. As two special cases, the deformations induced by a uniformly distributed shear load and a concentrated shear force are discussed in detail, respectively. The results indicate some interesting characteristics in nanocontact mechanics, which are distinctly different from those in macrocontact problem. At nanoscale, both the contact stresses and the displacements on the deformed surface transit continuously across the uniform distributed shear load boundary as a result of surface stress. In addition, the indent depth and the contact stress depend strongly on the surface stress for nanoindentation.


---

## Body

## 1. Introduction

Contact problems of an elastic half-space by applied loads on its plane surface are the subject of extensive literature [1]. The greater parts of these works of literature are concerned with the responses subjected to surface loads with only existing residual surface tension [2]. Corresponding problems of complete effects of surface stress (nonzero surface tension and surface elastic properties) have also attracted the attention of a number of authors and a considerable body of literature exists on this class of problems [3]. The surface of solids is a special region with very small thickness (a few times of atom-spacing). Since the equilibrium lattice spacing in the surface is different from that in the bulk, surface stress appears. For solids with large characteristic dimensions, the volume ratios of surface region to the bulk material is small; the effect of surface stress then can be neglected because of its relatively tiny contribution. However, for microsolids with large surface-to-bulk ratio the significance of surfaces is likely to be important, such as nanocontact problem [4]. This is extremely true for nanoscale materials or structures.When the characteristic size of an element or a solid approaches nanometers, its mechanical behavior displays remarkable size-dependent phenomena [5]. However, as there is no intrinsic length scale involved in the constitutive laws, the classical elastic theory cannot predict the size-dependent behavior of solids. Owing to the increasing ratio of surface-to-bulk volume, the effect of surface stress has been considered as one of the major factors contributing to the exceptional behaviors at nanoscale [2, 6–12].To account for the effect of surface stress in solid mechanics, Gurtin and Ian Murdoch [13, 14] developed a continuum theory of surface elasticity. For some elementary deformation modes, the prediction of surface elasticity showed a good agreement with directly atomic simulation. Therefore, the surface elasticity has been widely adopted to investigate the mechanical phenomena at nanoscale. Miller and Shenoy [15] studied the size-dependent elastic stiffness of structural elements such as nanobars, nanobeams, and nanoplates. Cammarata et al. [16] considered the size-dependent deformation in thin film with surface effect. Through atomic simulation, Shenoy [17] calculated the surface elastic constants of metallic fcc crystal surfaces. Gao et al. [18] developed a finite-element method to account for the effect of surface elasticity. Zhang and Wang [19, 20] investigated the effect of surface energy on the yield strength of nanoporous materials and nanowires. Wang and Feng [2] studied the response of a half-plane subjected to normal pressures with constant residual surface tension. Ou et al. [9–11] studied how the size dependence can be considered in the mechanical performance of nanoscale structures and devices. Zhao and Rajapakse [21] derived the fundamental solution of an elastic layer bonded to a rigid substrate with surface effects. Due to the complexity of the problem, however, until now there has been no systematic investigation into the complete effects of stress (nonzero surface tension and surface elastic properties) on an isotropic elastic bulk. In this paper, Fourier integral transform method is used to solve the nonclassical boundary-value problems with complete surface effects.In the present paper, the fundamental nanocontact problem of an isotropic elastic bulk with complete surface stress effects (nonzero surface tension and surface elastic properties) subjected to surface shear load is considered. The Fourier integral transform method is employed to obtain the fundamental solution of nanocontact problem with surface effects. This method allows us to easily extend our analysis to problems involving half-plane subjected to both concentrated shear load and uniformly distributed shear load on a finite region. The outline of the paper is organized as follows. The basic equations of surface elasticity are reviewed briefly in Section2. In Section 3, the Fourier integral transform method is adopted to solve the nanocontact problem with surface effects subjected to shear loads. The detailed results for the cases of uniformly distributed shear load and concentrated forces are discussed in Sections 4 and 5, respectively, and concluding remarks are presented in Section 6.


## 2. Basic Equations of Surface Elasticity

In surface elasticity theory, a surface is regarded as a negligibly thin membrane that has material constants different from the bulk material and is adhered to the bulk without slipping. The equilibrium and constitutive equations in the bulk of material are the same as those in classical elastic theory, but the presence of surface stresses gives rise to a nonclassical boundary condition. For further details, the reader may refer to Gurtin and Ian Murdoch [13, 14]. Here only several basic equations of surface elasticity theory are reviewed.In the absence of body force, the equilibrium equations, constitutive law, and geometry relations in the bulk are as follows:(1)σij,j=0σij=2G(εij+v1-2vεkkδij),
where G and v are the shear modulus and Poisson’s ratio of the bulk material, σij and εij are the stress tensor and strain tensor in the bulk material, respectively. Throughout the paper, Einstein’s summation convention is adopted for all repeated Latin indices (1, 2, 3) and Greek indices (1, 2).The strain tensor is related to the displacement vectorui by
(2)εij=12(ui,j+uj,i).Assume that the surface of the material adheres perfectly to its bulk without slipping. Then the equilibrium conditions on the surface are expressed as(3)σβαnβ+σβα,βs=0,σijninj=σαβsκαβ,
where ni denotes the normal to the surface, καβ the curvature tensor of the surface, and σαβs the surface stress tensor.The surface stress tensor is related to the surface strain tensor by(4)σβαs=τsδβα+2(μs-τs)εβα+(λs+τs)εγγδβα,
where τs is the residual surface tension under unstrained conditions, μs and λs are surface Lamé constants which can be determined by atom simulations or experiments [17].


## 3. Nanocontact Model with Surface Effects

Considering a material occupying the upper half-planez>0, we refer to a Cartesian coordinate system (o-xyz), as shown in Figure 1, where the x-axis is along the surface, and the z-axis perpendicular to the surface. The plane-strain conditions are assumed to ε2i=0, and the contact is assumed to be frictionless. In this case, the boundary conditions (3) on the contact surfacez=0 are simplified to
(5)σzz=0σxz+s(x)=-(dτsdx+ksd2udx2),
where s(x) is the shear load applied on the surface, and ks=2μs+λs is a surface material constant.Figure 1
Schematic of contact problem under uniform shear load.For the considered plane problem, the equilibrium equations and Hooke’s law in the bulk reduce to(6)∂σxx∂x+∂σxz∂z=0,∂σzz∂z+∂σxz∂x=0,(7)εxx=12G[(1-v)σxx-vσzz],εzz=12G[(1-v)σzz-vσxx],εxx=σxz2G.The strains are related to the displacements by(8)εxx=∂u∂x,εzz=∂w∂z,εxz=12(∂u∂z+∂w∂x).
which satisfy the following compatibility condition:
(9)∂2εxx∂z2+∂2εzz∂x2=2∂2εxz∂x∂z.
As in classical theory of elasticity, the Airy stress function χ(x,z) is defined by
(10)σxx=∂2χ∂z2,σzz=∂2χ∂x2,σxz=-∂2χ∂x∂z.
Then the equilibrium equations in (6) are satisfied automatically, and the compatibility equation in (9) becomes
(11)(∂2∂x2+∂2∂z2)(∂2χ∂x2+∂2χ∂z2)=0.To solve the boundary-value problem, the Fourier integral transformation method is adopted to the coordinatex. Then, the Airy-stress function χ(x,z) and its Fourier transformation χ~(ξ,z) can be expressed as
(12)χ~(ξ,z)=12π∫-∞∞χ(x,z)eixξdξ,χ(ξ,z)=12π∫-∞∞χ~(x,z)e-ixξdξ.Substituting (12) into (11) and considering the condition that the stresses vanish at infinity, one obtains
(13)χ~(ξ,z)=(A+Bz)e-z|ξ|,
where A and B are generally functions of ξ as yet to be determined.Substituting (13) and (12) into (10), the stresses can be written as
(14)σxx=12π∫-∞∞ξ2[A(ξ)+(z-2|ξ|)B(ξ)]e-ixξ-z|ξ|dξ,σzz=-12π∫-∞∞ξ2[A(ξ)+zB(ξ)]e-ixξ-z|ξ|dξ,σxz=i2π∫-∞∞ξ[(1-|ξ|z)B(ξ)-|ξ|A(ξ)]e-ixξ-z|ξ|dξ.By substituting the stresses into (7) and using (8), the displacements are derived as
(15)w(x,z)=122πG∫-∞∞[|ξ|A(ξ)+(1-2v+z|ξ|)B(ξ)]×e-ixξ-z|ξ|dξ,u(x,z)=i22πG∫-∞∞[(2-v)|ξ|A(ξ)+(z|ξ|-2(1-v))×B(ξ)|ξ|]e-ixξ-z|ξ|dξ.Substituting (15) and (14) into the surface condition (5) leads to
(16)A(ξ)=0,B(ξ)=is~(ξ)ξ11+bξ.
Here
(17)b=ks(1-v)G,s~(ξ)=12π∫-∞∞s(x)eixξdx,
where b is a length parameter depending on the surface property and material elastic constants. It should be pointed out that this parameter indicates the thickness size of the zone where the surface effect is significant and plays a critical role in the surface elasticity. For metals, b is estimated on the order of nanometers [15, 19, 22].In what follows, the two special cases including a uniform distributed shear load over a finite region and a concentrated shear force will be dicussed in detail, respectively.


## 4. Uniform Distributed Shear Load

In the case a uniform shear loads acts over the region |x|≤a, one has
(18)s~(ξ)=2πssin⁡(aξ)ξ.
Therefore B(ξ) is given by
(19)B(ξ)=i2πssin⁡(aξ)ξ211+bξ.
Then the stresses are obtained as
(20)σxx=2sπ∫0∞[zξ-2ξ(1+bξ)]sin⁡(aξ)sin⁡(xξ)e-zξdξ,σzz=-2sπ∫0∞[z1+bξ]sin⁡(aξ)sin⁡(xξ)e-zξdξ,σxz=-2sπ∫0∞[1-zξξ(1+bξ)]sin⁡(aξ)sin⁡(xξ)e-zξdξ.
The corresponding displacements are expressed as
(21)u(x,z)=-sπG∫0∞[(2-v)zξ-2(1-v)ξ2(1+bξ)]×sin⁡(aξ)cos⁡(xξ)e-zξdξ,w(x,z)=sπG∫0∞1-2v+zξξ2(1+bξ)sin⁡(aξ)sin⁡(xξ)e-zξdξ.On the contact surfacez=0, the stresses is given by
(22)σxx=-4sπ∫0∞(bat+1)-1sin⁡(t)tsin⁡(xat)dt,σxz=-2sπ∫0∞(bat+1)-1sin⁡(t)tcos⁡(xat)dt.A closed-form solution for the elastic field of the half-plane cannot be obtained due to the complexity of the integrals involved in the solution. However, whenb=0, that is, the surface influence is ignored in (22), the stresses of the half-plane are in agreement with those in the classical elastic results [1], and are expressed as
(23)σxx=s2π[4ln⁡(r1r2)-(cos⁡(2θ1)-cos⁡(2θ2))],σxz=-s2π[2(θ2-θ1)+(sin⁡(2θ1)-sin⁡(2θ2))],
where tanθ1,2=z/(x∓a) and r1,2=[(x∓a)2+z2]1/2.Assuming that the origin has no displacement in thez direction, that is, w(0,0)=0, one obtains
(24)w(x,0)=(1-2v)asπG∫0∞(bat+1)-1sin⁡(t)t2sin⁡(xat)dt.If the shear displacementu is specified to be zero at a distance l on the contact surface, that is, u(l,0)=0, the displacement on the surface is derived as
(25)u(x,0)=2(1-v)asπG∫0∞(bat+1)-1sin⁡(t)t2×[cos⁡(xat)-cos⁡(lat)]dt.It is instructive to examine the influence of the surface stress on the stresses and displacements of the contact surface and compare them with those in classical contact problem. Surface stress can be obtained from experiments [2] or atomic simulations [15]. According to these results, the influence of surface stresses, characterized by the intrinsic length, b=ks(1-v)/G, becomes significant at nanoscale. For the illustration of the effects of surface energy, b/a is taken as 0, 0.2, 0.5, and 2 in our calculations. Figures 2 and 3 show the distribution of the stresses σxx and σxz on the contact surface, where the solution of b/a=0 is consistent with the classical elastic result. When the loading size a is comparable to the parameter b, that is, of the order of nanometers, the influence of surface stress is evidently significant. It can be seen from Figure 2 that the normal stress transits continuously across the loading boundary x=±a, which is opposed to a singularity predicted by classical elasticity. In addition, the maximum normal stress in the bulk increases with an increase in surface stress. In the loading region (|x/a|<1.0), the normal stress is decreasing monotonically with respect to x, and the inverse is observed outside the loading region (|x/a|>1.0). It is also found in Figure 3 that shear stress changes smoothly across the loading boundary x=±a, which is different from a rapid jump from zero outside the load region to s inside the load region predicted by the classical contact problem. In addition, the maximum shear stress in the bulk decreases continuously with the increase of surface stress.Figure 2
Distribution of contact normal stress under a uniform shear load.Figure 3
Distribution of contact shear stress under a uniform shear load.Due to the different surface stress value, the surface displacement on the contact surface is plotted in Figure4 with m1=(1-2v)s/πG which also shows that the gradient of the surface displacement is continuous everywhere on the deformed surface, though the classical elasticity gives only discontinuous displacement gradient at x=±a. With the increase of surface stress, the indent depth decreases continuously. The horizontal displacement is displayed in Figure 5, where we set l=5a. and m2=2(1-v)s/πG. It is seen that the slope of the deformed surface for a>0 is continuous everywhere. However, the classical elasticity theory predicted unreasonably that the gradient of the deformed surface is infinite at the load boundaryx=±a, as seen from the curve of b/a=0.Figure 4
Surface displacement under a uniform shear load.Figure 5
Surface indentation under a uniform shear load.


## 5. Concentrated Shear Load

Now the extreme situation that a point forceS acts on the surface at the origin is considered by setting 2sa→S as a→0. Thus, one has
(26)B(ξ)=i12πSξ11+bξ.Substituting (26) into (14) and (15), the corresponding stress and displacement distributions are given as follows, respectively,
(27)σxx=Sπ∫0∞zξ-21+bξsin⁡(xξ)e-zξdξ,σzz=-Sπ∫0∞zξ1+bξsin⁡(xξ)e-zξdξ,σxz=-Sπ∫0∞1-zξ1+bξcos⁡(xξ)e-zξdξ,(28)u(x,z)=-SπG∫0∞(2-v)zξ-2(1-v)ξ(1+bξ)cos⁡(xξ)e-zξdξ,w(x,z)=SπG∫0∞1-2v+zξξ(1+bξ)sin⁡(xξ)e-zξdξ.However, whenb=0, which implies that the surface influence is ignored, the stresses of the half-plane can be obtained as
(29)σxx=-2Sπx3(x2+z2)2,σxz=-2Sπx2z(x2+z2)2
which are consistent with those in classical contact mechanics [1].On the contact surface, the stresses and displacements are(30)σxx=-G1-v∫0∞(bΛ2t+1)-1sin⁡(xΛ2t)dt,σxz=-G1-2v∫0∞(bΛ2t+1)-1cos⁡(xΛ2t)dt,(31)u(x,0)=Λ1∫0∞1t(bΛ2t+1)-1[cos⁡⁡(xΛ2t)-cos⁡⁡(lΛ2t)]dt,w(x,0)=Λ2∫0∞1t(bΛ2t+1)-1sin⁡(xΛ2t)dt,
where
(32)Λ1=2(1-v)SπG,Λ2=(1-2v)SπG.The variation of the distribution of stresses on the contact surface subjected to a concentrated shear force with various values ofb/Λ2is shown in Figures 6 and 7. In classical elasticity without surface effect (i.e., b/Λ2=0), the stresses tend to infinity at the loading point x=0, while it approaches finite values when surface stresses are present. And the stresses σxz and σxx at x=0 decrease with the increase of surface stress characterized by b/Λ2. Figure 8 shows the maximum surface displacement in the bulk decreases with the increase of surface stress. Figure 9 shows the surface indentation for several representative values of surface stresses. Different from the classical elastic solution without surface effects (i.e., b/Λ2=0), in which the indent depth at x=0 is maximum, the indent depth will decreases when the effect of surface stress is taken into account. The indent depth decreases as the surface stress increases.Figure 6
Distribution of the contact normal stress under a concentrated shear load.Figure 7
Distribution of the contact shear stress under a concentrated shear load.Figure 8
Surface displacement under a concentrated shear load.Figure 9
Surface indentation under a concentrated shear load.Though only two special cases have been considered in contact problems with surface stresses, based on the solution of a point force applied on a semi-infinite body, one can deal with various types of contact problems by solving an integral equation, as those have been done in the classical contact mechanics.


## 6. Conclusions

The two-dimensional nanocontact problems for elastic bulk materials subjected to surface shear loads are investigated with completed surface stress effect (nonzero surface tension and surface elastic properties). A set of analytical solutions are presented by using Fourier integral transform method. The closed-form solution can be derived for the case of a half-plane. For the two particular cases of a uniform distributed shear load and a concentrated shear force, the results are analyzed in detail and compared with the classical linear elastic solutions. It is found that the surface elasticity theory illuminates some interesting characteristics of nanocontact problems, which are distinctly different from the classical solutions of elasticity without surface effects. The contact stresses and the displacements of the deformed surface change smoothly across the loading boundary; moreover the stress and the indent depth show a significant dependence on the surface stress. Therefore the effects of surface stresses should be considered for nanocontact problems.



---

*Source: 505034-2012-06-28.xml*