Journal of Nanomaterials

Volume 2012, Article ID 505034, 7 pages

http://dx.doi.org/10.1155/2012/505034

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

School of Science, Lanzhou University of Technology, Lanzhou 730050, China

Received 24 February 2012; Accepted 15 May 2012

Academic Editor: Sergio J. Mejía-Rosales

Copyright © 2012 D. X. Lei et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### 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.

#### 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 Section 2. 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: where and are the shear modulus and Poisson’s ratio of the bulk material, and 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 vector by

Assume that the surface of the material adheres perfectly to its bulk without slipping. Then the equilibrium conditions on the surface are expressed as where denotes the normal to the surface, the curvature tensor of the surface, and the surface stress tensor.

The surface stress tensor is related to the surface strain tensor by where is the residual surface tension under unstrained conditions, and 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-plane , we refer to a Cartesian coordinate system (), as shown in Figure 1, where the -axis is along the surface, and the -axis perpendicular to the surface. The plane-strain conditions are assumed to , and the contact is assumed to be frictionless. In this case, the boundary conditions (3) on the contact surface are simplified to where is the shear load applied on the surface, and is a surface material constant.

For the considered plane problem, the equilibrium equations and Hooke’s law in the bulk reduce to

The strains are related to the displacements by which satisfy the following compatibility condition: As in classical theory of elasticity, the Airy stress function is defined by Then the equilibrium equations in (6) are satisfied automatically, and the compatibility equation in (9) becomes

To solve the boundary-value problem, the Fourier integral transformation method is adopted to the coordinate . Then, the Airy-stress function and its Fourier transformation can be expressed as

Substituting (12) into (11) and considering the condition that the stresses vanish at infinity, one obtains where and are generally functions of as yet to be determined.

Substituting (13) and (12) into (10), the stresses can be written as

By substituting the stresses into (7) and using (8), the displacements are derived as

Substituting (15) and (14) into the surface condition (5) leads to Here where 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, 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 load acts over the region , one has Therefore is given by Then the stresses are obtained as The corresponding displacements are expressed as

On the contact surface , the stresses is given by

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, when , 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 where and .

Assuming that the origin has no displacement in the direction, that is, , one obtains

If the shear displacement is specified to be zero at a distance on the contact surface, that is, , the displacement on the surface is derived as

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, , becomes significant at nanoscale. For the illustration of the effects of surface energy, is taken as 0, 0.2, 0.5, and 2 in our calculations. Figures 2 and 3 show the distribution of the stresses and on the contact surface, where the solution of is consistent with the classical elastic result. When the loading size is comparable to the parameter , 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 , 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 (), the normal stress is decreasing monotonically with respect to , and the inverse is observed outside the loading region (). It is also found in Figure 3 that shear stress changes smoothly across the loading boundary , which is different from a rapid jump from zero outside the load region to 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.

Due to the different surface stress value, the surface displacement on the contact surface is plotted in Figure 4 with 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 . With the increase of surface stress, the indent depth decreases continuously. The horizontal displacement is displayed in Figure 5, where we set . and . It is seen that the slope of the deformed surface for is continuous everywhere. However, the classical elasticity theory predicted unreasonably that the gradient of the deformed surface is infinite at the load boundary, as seen from the curve of .

#### 5. Concentrated Shear Load

Now the extreme situation that a point force acts on the surface at the origin is considered by setting as . Thus, one has

Substituting (26) into (14) and (15), the corresponding stress and displacement distributions are given as follows, respectively,

However, when , which implies that the surface influence is ignored, the stresses of the half-plane can be obtained as which are consistent with those in classical contact mechanics [1].

On the contact surface, the stresses and displacements are where

The variation of the distribution of stresses on the contact surface subjected to a concentrated shear force with various values of is shown in Figures 6 and 7. In classical elasticity without surface effect (i.e., ), the stresses tend to infinity at the loading point , while it approaches finite values when surface stresses are present. And the stresses and at decrease with the increase of surface stress characterized by . 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., ), in which the indent depth at 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.

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.

#### Acknowledgments

The supports from the National Natural Science Foundation (Grant no. 11062004), the Natural Science Foundation of Gansu province (1010RJZA021) and the Doctoral Foundation of Lanzhou University of Technology are acknowledged.

#### References

- K. L. Johnson,
*Contact Mechanics*, Cambridge University Press, Cambridge, Mass, USA, 1985. - G. F. Wang and X. Q. Feng, “Effects of surface stresses on contact problems at nanoscale,”
*Journal of Applied Physics*, vol. 101, Article ID 013510, 2007. View at Google Scholar - M. E. Kassner, S. Nemat-Nasser, Z. Suo et al., “New directions in mechanics,”
*Mechanics of Materials*, vol. 37, no. 2-3, pp. 231–259, 2005. View at Publisher · View at Google Scholar · View at Scopus - Y. Huang, F. Zhang, K. C. Hwang, W. D. Nix, G. M. Pharr, and G. Feng, “A model of size effects in nano-indentation,”
*Journal of the Mechanics and Physics of Solids*, vol. 54, no. 8, pp. 1668–1686, 2006. View at Publisher · View at Google Scholar · View at Scopus - R. Dingreville, J. Qu, and M. Cherkaoui, “Surface free energy and its effect on the elastic behavior of nano-sized particles, wires and films,”
*Journal of the Mechanics and Physics of Solids*, vol. 53, no. 8, pp. 1827–1854, 2005. View at Publisher · View at Google Scholar · View at Scopus - A. J. Kulkarni and M. Zhou, “Surface-effects-dominated thermal and mechanical responses of zinc oxide nanobelts,”
*Acta Mechanica Sinica*, vol. 22, pp. 217–224, 2006. View at Publisher · View at Google Scholar - B. Movchan and I. S. Jones, “Asymptotic and numerical study of a surface breaking crack subject to a transient thermal loading,”
*Acta Mechanica Sinica*, vol. 22, pp. 22–27, 2006. View at Publisher · View at Google Scholar - G. F. Wang and X. Q. Feng, “Interface effects on effective elastic moduli of nanocrystalline materials,”
*Materials Science and Engineering A*, vol. 363, pp. 1–8, 2003. View at Publisher · View at Google Scholar - Z. Y. Ou, G. F. Wang, and T. J. Wang, “Effect of residual surface tension on the stress concentration around a nanosized spheroidal cavity,”
*International Journal of Engineering Science*, vol. 46, pp. 475–485, 2008. View at Publisher · View at Google Scholar - Z. Y. Ou, G. F. Wang, and T. J. Wang, “Elastic fields around a nanosized spheroidal cavity under arbitrary uniform remote loadings,”
*European Journal of Mechanics A/Solids*, vol. 28, pp. 110–120, 2009. View at Publisher · View at Google Scholar - Z. Y. Ou, G. F. Wang, and T. J. Wang, “An analytical solution for the elastic fields near spheroidal nano-inclusions,”
*Acta Mechanica Sinica*, vol. 25, pp. 821–830, 2009. View at Publisher · View at Google Scholar - Z. Y. Ou and S. D. Pang, “A screw dislocation interacting with a coated nano-inhomogeneity incorporating interface stress,”
*Materials Science and Engineering A*, vol. 528, no. 6, pp. 2762–2775, 2011. View at Publisher · View at Google Scholar · View at Scopus - M. E. Gurtin and A. Ian Murdoch, “A continuum theory of elastic material surfaces,”
*Archive for Rational Mechanics and Analysis*, vol. 57, no. 4, pp. 291–323, 1975. View at Publisher · View at Google Scholar · View at Scopus - M. E. Gurtin and A. Ian Murdoch, “Surface stress in solids,”
*International Journal of Solids and Structures*, vol. 14, no. 6, pp. 431–440, 1978. View at Google Scholar · View at Scopus - R. E. Miller and V. B. Shenoy, “Size-dependent elastic properties of nanosized structural elements,”
*Nanotechnology*, vol. 11, no. 3, pp. 139–147, 2000. View at Publisher · View at Google Scholar · View at Scopus - R. C. Cammarata, K. Sieradzki, and F. Spaepen, “Simple model for interface stresses with application to misfit dislocation generation in epitaxial thin films,”
*Journal of Applied Physics*, vol. 87, no. 3, pp. 1227–1234, 2000. View at Google Scholar · View at Scopus - V. B. Shenoy, “Atomistic calculations of elastic properties of metallic fcc crystal surfaces,”
*Physical Review B*, vol. 71, Article ID 094104, 2005. View at Google Scholar - W. Gao, S. Yu, and G. Huang, “Finite element characterization of the size-dependent mechanical behaviour in nanosystems,”
*Nanotechnology*, vol. 17, no. 4, pp. 1118–1122, 2006. View at Publisher · View at Google Scholar · View at Scopus - W. X. Zhang and T. J. Wang, “Effect of surface energy on the yield strength of nanoporous materials,”
*Applied Physics Letters*, vol. 90, Article ID 063104, 2007. View at Google Scholar - W. X. Zhang, T. J. Wang, and X. Chen, “Effect of surface stress on the asymmetric yield strength of nanowires,”
*Journal of Applied Physics*, vol. 103, Article ID 123527, 2008. View at Google Scholar - X. J. Zhao and R. K. N. D. Rajapakse, “Analytical solutions for a surface-loaded isotropic elastic layer with surface energy effects,”
*International Journal of Engineering Science*, vol. 47, pp. 1433–1444, 2009. View at Publisher · View at Google Scholar - F. Q. Yang and, “Size-dependent effective modulus of elastic composite materials: spherical nanocavities at dilute concentrations,”
*Journal of Applied Physics*, vol. 95, pp. 3516–3520, 2004. View at Publisher · View at Google Scholar