Abstract

This paper presents the analysis of a layered elastic half space under the action of axisymmetric surface loading and the influence of the surface energy effects. The boundary value problems for the bulk and the surface are formulated based on classical linear elasticity and a complete Gurtin-Murdoch constitutive relation. An analytical technique using Love’s representation and the Hankel integral transform is employed to derive an integral-form solution for both displacement and stress fields. An efficient numerical quadrature is then applied to accurately evaluate all involved integrals. Selected numerical results are presented to portray the influence of various parameters on elastic fields. Numerical results indicate that the surface stress displays a significant influence on both displacement and stress fields. It is also found that the layered half space becomes stiffer with the presence of surface stresses. In addition, unlike the classical elasticity solution, size-dependent behavior of elastic fields is noted. The present analytical solutions provide fundamental understanding of the influence of surface energy on layered elastic materials. It can also be used as a benchmark solution for the development of numerical techniques such as FEM and BEM, for analysis of more complex problems involving a layered medium under the influence of surface energy effects.

1. Introduction

Nanotechnology has received wide attention in recent years due to its vast applications in various disciplines such as biology, chemistry, physics, medicines, material sciences, and engineering. In the fields of material sciences and engineering, studies related to mechanical behavior of nanostructured materials have become a subject of numerous investigations due to the fact that understanding fundamental aspects of their behaviors at nanoscale level is important for optimum design of nanosized devices and structures. There are two approaches that have commonly been employed to theoretically investigate mechanical behaviors of materials at nanoscale, namely, atomistic simulation and modified continuum-based model. Atomistic modeling techniques require a very large computational effort, although they are considered very accurate. A modified continuum-based model then becomes an attractive alternative in obtaining first-approximation to predict mechanical behaviors of nanostructured materials. Due to their high surface-to-volume ratio, nanoscale elements usually exhibit high influence of surface/interface free energy, which is the energy associated with atoms at or near a free surface (e.g., see [1]); consequently, their mechanical behavior becomes size-dependent [2]. Thus, surface energy effects, which are generally ignored in conventional continuum mechanics problems, need to be taken into account in modified continuum-based simulation for nanoscale systems. A theoretical framework based on continuum mechanics concepts was proposed by Gurtin and Murdoch [3, 4] to take into consideration the influence of surface energy effects. In their model, an elastic surface was formed as a mathematical layer of zero thickness perfectly bonded to the underlying bulk material without slipping. Several studies were carried out to verify that modified continuum-based simulations with surface energy effects and size-dependency can be employed to model nanostructured elements with acceptable accuracy. For instance, Miller and Shenoy [5] examined the size-dependent behavior of nanostructured elements (i.e., bar, beam, and plate) by adopting the Gurtin-Murdoch model and found that their results were in a good agreement with those obtained from direct atomistic simulations. Dingreville et al. [6] developed a continuum framework to incorporate the surface free energy in the framework of continuum mechanics and demonstrated that overall mechanical behaviors of nanostructured elements such as particles, wires, and films were found to be size-dependent. There also exist other continuum-based theories that have been developed to take into account the size-dependent material behaviors at the nanoscale level such as the strain gradient elasticity theory by Mindlin [7]. The theory proposed by Mindlin has not been widely adopted in the modeling of nanoscale systems since it involves several additional material parameters and higher-order governing equations. Simplified versions of Mindlin’s theory have then been proposed, and analytical solutions to various continuum mechanics problems were presented based on its simplified versions (e.g., see [810]).

Over the last two decades, several researchers have investigated a variety of continuum mechanics problems by adopting the Gurtin-Murdoch theory of surface elasticity. For example, Huang and Yu [11] studied an elastic half plane under surface loading with consideration of surface energy effects. An elastic layer with finite thickness, subjected to surface loading under plane-strain and axisymmetric conditions, was also considered by Zhao and Rajapakse [12]. Intarit et al. [13] derived fundamental solutions of an elastic half plane under internal loading and dislocations. An elastic half plane under surface shear loading was also investigated by Lei et al. [14]. Recently, a nanocontact problem of layered viscoelastic solids with surface energy effects was presented by Abdel Rahman and Mahmoud [15]. All these studies, however, considered the surface stress tensor as a 2D quantity with its out-of-plane components being neglected. Wang et al. [16] showed that the out-of-plane terms of the surface displacement gradient could be significant even in the case of small deformations particularly for curved and rotated surfaces. The complete version of Gurtin-Murdoch model, with consideration of the out-of-plane term, has later been employed to examine various continuum mechanics problems, for example, problems related to an internally loaded elastic layer under plane-strain condition [17] and axisymmetric loading [18], respectively; contact problem [19]; nanoindentation [20, 21]; nanobeams [22]; nanoplate [23]; and nanosized cracks [24, 25]. In addition, the influence of surface energy is also significant in problems related to soft elastic solids [26].

Stress analysis of a layered elastic medium under applied surface loading has a rich history (e.g., see [2730]) due to its close relevance to various engineering applications, such as characterization of mechanical properties of layered materials: for example, protective coatings, multilayer capacitors, and layered composite materials; analysis and design of pavement and foundations; and in situ testing of soils and rocks and so forth. A review of literature indicates that studies related to a layered elastic medium with consideration of surface energy effects based on the Gurtin-Murdoch theory are very limited. This class of problems has extensive applications in the study of nanocoatings and nanoscale surface layers that are used in electronic devices, tribological and biomaterial applications, advanced industrial materials, communication devices, and so forth. The main objective of this paper is to present analytical solutions to a layered elastic half space under axisymmetric surface loading by adopting the complete Gurtin-Murdoch theory of surface elasticity. The boundary value problems of a layered elastic half space under axisymmetric surface loading involving nonclassical boundary conditions due to surface stress influence are formulated by employing Love’s strain potential and the Hankel integral transform. Selected numerical results for displacements and stresses due to applied vertical and radial loading are presented to portray the influence of layer thickness, surface material parameters, and size-dependency on elastic fields. The present fundamental solution is useful in the development of boundary integral equation methods for the investigation of more complicated problems such as nanoindentation and contact problems involving a layered elastic medium. In addition, the present numerical results can also be employed as a benchmark solution in the development of numerical techniques such as finite element and boundary element methods for analysis of a variety of problems with the influence of surface energy such as nanoscale problems and soft elastic solids.

2. Governing Equations and General Solutions

Consider a layered elastic half space consisting of two elastic materials with different properties perfectly bonded together, in which the upper material is an elastic layer of finite thickness and subjected to axisymmetric vertical and radial surface loads denoted by and , respectively, as shown in Figure 1. According to the Gurtin-Murdoch surface elasticity theory, both materials consist of two parts, the bulk material and the surface, which is a zero-thickness layer perfectly bonded to the bulk material without slipping. The equilibrium equations, the constitutive equations, and the strain-displacement relationship of the bulk material under axisymmetric deformations are the same as those in the classical elasticity theory, which are given, respectively, bywhere denote the components of stress tensors; denote the components of strain tensors; and denote the components of displacement tensors, respectively. In addition, and are Lamé constants of a bulk material. On the surface, the equilibrium conditions in terms of the generalized Young-Laplace equation [31], the surface constitutive relations, and the strain-displacement relationship can be expressed, respectively, as [3, 4]where the superscript “” is used to denote the quantities corresponding to the surface; and are surface Lamé constants; is the residual surface stress (or surface tension) under unstrained conditions; and denotes the prescribed traction on the surface. Equation (4) can be viewed as the out-of-plane contribution of the preexisting surface tension in the deformed configuration whereas the surface gradient of the displacement acts as the out-of-plane component of the unit vector tangent to the surface in the deformed state. This term has been ignored in several previous studies even though the contribution of could be significant even in the case of small deformations (e.g., see [17, 18, 20]).

The general solutions for the stresses and displacements in the bulk material under axisymmetric deformations can be expressed by using Love’s strain potential and Hankel integral transform as [32]whereIn addition, , , , and are the arbitrary functions to be determined from the appropriate boundary conditions.

3. Solutions for Axisymmetric Surface Loading

Boundary value problem of a layered elastic half space subjected to axisymmetric normal and tangential traction, denoted by and , respectively, applied at its surface as shown in Figure 1 is considered in this section. To solve this problem, the layered half space is divided into two domains. Domain “1” represents the upper layer and domain “2” represents the underlying half space. The general solutions of the bulk material in domain “1” are given by (6) whereas those of domain “2” can also be obtained from (6) by replacing the arbitrary functions to with the arbitrary functions to , respectively. Note that and are imposed to ensure the regularity of the solutions at infinity for domain “2”; in addition, the subscript is used to denote the quantities corresponding to domains “1” and “2,” respectively. The solutions of to can be determined by solving the following boundary and continuity conditions:where is a surface material constant. It should be noted that (8) are nonclassical boundary conditions obtained from (2) to (5). In view of (6) together with the assumption that the surface residual stress is constant, the following six linear algebraic equations are established to solve for the arbitrary functions to :where the following nondimensional quantities in the above equations are defined thus: ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ; and . In addition, the functions and are obtained from the surface loads and , respectively, asin which ;  ; and . The arbitrary functions to for given functions of the applied surface loads and can then be obtained by solving the linear equation system (9). Substitution of the arbitrary functions to into (6) yields the displacement and stress fields at an arbitrary point of the layered elastic half space shown in Figure 1.

4. Numerical Results and Discussion

The complete solutions of displacements and stresses corresponding to the boundary value problems of a layered elastic medium under axisymmetric surface loading shown in Figure 1 are given by (6) with the arbitrary functions to obtained by solving (9). The solutions of to are obtained explicitly for applied normal and tangential traction, p(r) and , respectively, but their expressions for each loading case are not presented here for brevity. Given the complexity of the arbitrary functions to , closed-form solutions to the displacement and stress fields cannot be obtained. Therefore it is essential to determine all elastic fields by numerically evaluating the semi-infinite integrals appearing in (6). It is found that those semi-infinite integrals with respect to ξ  can be accurately evaluated by employing the numerical quadrature scheme based on 21-point Gauss-Kronrod rule [33]. The accuracy of the present solution is first verified by comparing with the existing solution given by Gerrard [27], who presented the classical solutions (without the influence of surface energy effects) of a layered elastic half space subjected to axisymmetric surface loading. Table 1 presents a comparison of normalized displacements at the surface and normalized stresses at the interface along the radial direction of a layered elastic half space under uniformly distributed normal traction , acting over a circular area of radius at the surface. The comparison of surface displacements and stresses at the interface of the layered half space under linearly distributed shear traction applied over a circular area of radius at the surface is also shown in Table 2. In addition, with Poisson’s ratio and are considered for the numerical results presented in both tables. The solutions for normalized displacements and stresses from the present study are obtained by setting the parameters associated with the surface energy effects to be zero; that is, and . It is evident that excellent agreement between the two solutions is observed in both displacements and stresses shown in Tables 1 and 2.

Numerical results for vertical and radial displacements and vertical and shear stresses corresponding to a layered elastic half space with the influence of surface energy effects subjected to axisymmetric surface loading as shown in Figure 1 are presented next. Two cases of axisymmetric surface loading, namely, the vertical loading and the radial loading, are considered in the numerical study. The vertical loading denotes the case where uniformly distributed normal traction is applied over a circular area of normalized radius . The radial loading represents the case where the layered half space is subjected to linearly distributed tangential traction over a circular area of normalized radius , where is the maximum traction at the edge of the loading region. The functions defined as shown in (10) are given, respectively, for the vertical loading and the radial loading as follows:In addition, the numerical results presented hereafter correspond to the case where the material for the upper layer (domain “1”) is whereas is chosen for the underlying half space (domain “2”), respectively. The material properties for both domains are given in Table 3 [5].

Figure 2 presents radial variations of nondimensional displacements at the top surface and nondimensional stresses at the interface of a layered elastic half space under the vertical loading for different values of normalized thickness of the top layer . Note that the stress profiles in all figures presented in this section are computed at the interface at the bulk material of the underlying half space. Figure 2(a) shows radial profiles of vertical and radial surface displacements for various values of whereas the profiles of normal and shear stresses at the interface are illustrated in Figure 2(b). The classical solutions, also presented in these figures for comparison, are obtained by setting the parameters associated with the surface energy effects to be zero; that is, and . It is evident from Figure 2 that although the results from the present study and the classical solution display similar trends for both displacements and stresses at all values of , the surface energy effect renders the layered medium stiffer. The present solution yields lower surface displacements and stresses at the interface. The influence of surface energy is however less significant in the interface stresses, especially in the case of the shear stress. It is also found that the magnitude of all displacements and stresses decrease with increasing the normalized thickness of the layer since the upper layer is stiffer than the underlying half space (Lamé constants of Si [] are higher than those of Al []). In addition, as the layer thickness increases both vertical and radial surface displacements move towards the homogeneous half-space solutions presented by Intarit [34], and both solutions are virtually identical when .

Radial profiles of normalized surface displacements and normalized stresses at the interface of the layered elastic medium under the vertical loading are shown in Figure 3 to demonstrate the influence of the residual surface stress on elastic fields. The values of the residual surface stress in the underlying half space are varied (i.e., ) whereas other material parameters associated with both upper layer and underlying half space given in Table 3 remain unchanged. In addition, the normalized thickness of is considered in the numerical results shown in this figure. Once again, the influence of the surface stress is clearly observed from the displacement and stress solutions presented in Figure 3. The values of all displacements and stresses from the present study are substantially reduced from their classical elasticity counterparts as the value of the residual surface stress increases.

The next numerical results are presented to demonstrate the size-dependent behavior of the present solution when the influence of surface energy effects is considered. Figure 4 shows radial variations of vertical and radial surface displacements and the vertical and shear stresses at the interface of the layered half space under the vertical loading for different values of the normalized radius of loading area (i.e., ). In addition, the thicknesses of the top layer and the circular loading area are varied while their ratio is maintained at . Note that the solution when corresponds to the case where the thickness of the layer is equal to the characteristic length (). The corresponding nondimensional solution for the classical elasticity case is also shown, and it is size-independent. The size-dependency of the present solution is clearly observed in all displacement and stress profiles. It is evident from the numerical results presented in Figure 4 that the present solution accounting for surface energy effects approaches the classical solution as the loading radius increases. This is consistent with the fact that a larger loading area would produce higher displacements and stresses.

The final set of the numerical results correspond to the case where the layered elastic half space is subjected to the radial loading, in which the tangential traction is applied linearly distributed over a circular area of normalized radius . Figure 5 presents radial profiles of nondimensional displacements at the top surface and nondimensional stresses at the interface for different values of . It is evident from Figure 5 that both displacements and stresses of the layered half space under the radial loading depend more significantly on surface energy effects for all values of when compared to the results presented in Figure 2 under the vertical loading case. The presence of surface stresses significantly lowers the magnitude of all displacements and stresses shown in Figure 5. In addition, all displacements and stresses are reduced as the normalized thickness of the layer increases. Once again, both vertical and radial surface displacements are practically the same as the half-space solutions given by Intarit [34] when similar to what was observed in the vertical loading case.

5. Conclusions

An analytical treatment of a layered elastic half space under axisymmetric surface loading, taking into account the influence of surface energy effects, is presented in this paper. The boundary value problem corresponding to a layered elastic half space subjected to axisymmetric normal and tangential traction is formulated based on the complete Gurtin-Murdoch theory of surface elasticity and the application of Hankel integral transforms. An efficient numerical quadrature scheme is developed to evaluate all involved integrals. Selected numerical results for radial profiles of displacements and stresses are presented to portray the influence of various parameters on elastic fields. The numerical results indicate that the surface energy effects play an important role in both stress and displacement fields of a layered elastic medium. Surface residual stress reduces the surface displacements and stresses at the interface of the layered medium. Unlike the classical elasticity solution, the present study shows substantial size-dependency of elastic fields. Both displacement and stress fields show strong dependency on the loading area, and the magnitudes of nondimensional displacements and stresses increase with increasing the radius of the loading area. The application of the present fundamental solution for nanoindentation and contact problems is currently underway and will be presented in a future publication. In addition, the present solution can also be used as a benchmark solution for assessing the accuracy of numerical models such as the finite element and boundary element methods, which can be used to investigate more complicated problems in the presence of surface energy effects.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The work presented in this paper was supported by the Thailand Research Fund (Grants no. BRG 5880017 and PHD/0096/2554). The support is gratefully acknowledged.