A N A L Y T I C A L A N D F I N I T E E L E M E N T S O L U T I O N S F O R A N E L A S T I C M A T R I X W I T H T W O - D I M E N S I O N A L N A N O S C A L E I N H O M O G E N E I T I E S By Lian Tian B. S., University of Science and Technology of China, P. R. China, 2003 A THESIS SUBMITTED IN PARTIAL FULLFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF APPLIED SCIENCE in THE FACULTY OF GRADUATE STUDIES (MECHANICAL ENGINEERING) THE UNIVERSITY OF BRITISH COLUMBIA June 2006 © Lian Tian, 2006 Abstract A new frontier of research at the interface of material science and mechanics of solids has emerged with the current focus on the development of nanomaterials such as nanotubes, nanowires and nanoparticles. As the dimensions of a structure approach the nanoscale, its properties can be size-dependent. The classical continuum theory, however, does not admit an intrinsic size, and is not applicable to the analysis of nanoscale materials and structures. Mechanics of nanomaterials-based composites can be understood by incorporating the effects of surface and interfacial energy. A fundamental problem in the study of behaviour of such materials is the examination of size-dependent elastic field of an elastic matrix with nanoscale inhomogeneities. Classical inhomogeneity problems have a rich history since the celebrated work of Eshelby. However, the classical solutions cannot be applied to study nanoscale inhomogeneity problems and new solutions accounting for surface/interface energy have to be derived. This thesis therefore presents an analytical scheme and a finite element formulation to study the size-dependent elastic field of an elastic matrix containing two-dimensional nanoscale inhomogeneities. The Gurtin-Murdoch surface/interface elasticity model is applied to incorporate the surface/interface energy effects. By using the complex potential technique of Muskhelishvili, a closed-form analytical solution is obtained for the elastic field of a nanoscale circular inhomogeneity in an infinite matrix under arbitrary remote loading and a uniform eigenstrain. In the case of an elliptical inhomogeneity, the analytic potential functions are obtained approximately. A new finite element formulation that takes into account the surface stress effects is presented. Elastic field is found to depend on the characteristic dimensions of the inhomogeneity, surface elastic constants and surface residual stress. A striking feature of the new solutions is the existence of singular elastic fields below some dimensions of the inhomogeneity. This phenomenon requires careful further investigation. Eshelby tensor of a nanoscale circular inhomogeneity in an infinite matrix due to a uniform eigenstrain is uniform but becomes size-dependent; however, the tensor is size-dependent and non-uniform in the case of an elliptical inhomogeneity. i i Table of Contents Abstract ii Table of Contents iii List of Tables v List of Figures . vi List of Symbols x Acknowledgements xi 1. Introduction 1 1.1 Nanotechnology 1 1.2 Nanoscale Mechanics 3 1.3 Review of Surface Elasticity Model ...5 1.4 Review of Classical and Modern Inhomogeneity Problems 10 1.5 Scope of the Current Work 11 2. Elastic Field of an Infinite Matrix with a Nanoscale Circular Inhomogeneity...14 2.1 Problem Description and Basic Equations 14 2.2 Formulation of Inhomogeneity Problem 15 2.3 Numerical Results for Elastic Field Around a Circular Hole 22 2.4 Numerical Results for Elastic Field of a Circular Inhomogeneity 29 3. Elastic Field of an Infinite Matrix with a Nanoscale Elliptical Inhomogeneity...33 3.1 Problem Description 33 3.2 Solution of Elastic Field 34 3.3 Numerical Results for Elastic Field Around an Ell iptical Hole 42 3.4 Numerical Results for Elastic Field of an Ell iptical Inhomogeneity 50 4. 2-D Finite Element Analysis of Nanoscale Inhomogeneity Problems 56 4.1 Finite Element Formulation 56 4.2 Finite Element Simulation of Inhomogeneity Problems 62 4.2.1 Elastic Field Around a Hole in a Finite Plate 62 4.2.2 Elastic Field of a Finite Plate with Two Circular Inhomogeneities 72 4.2.3 Elastic Field Around a Circular Hole in a Finite Anisotropic Plate 76 5. Summary and Conclusions '. 79 i i i 5.1 Summary and Major Findings 79 5.2 Suggestions for Future Work 81 Bibliography 83 Appendices 88 Appendix A 88 Appendix B 92 Appendix C 92 Appendix D 94 iv List of Tables Table 2.1 Surface elastic constants. Units are N/m. 22 v List of Figures Figure 1.1 Potential benefits of nanoscience & nanotechnology 3 Figure 2.1 A nanoscale circular inhomogeneity in an infinite matrix 14 Figure 2.2 Variation of stress concentration factor at 0 = 0 24 Figure 2.3 Nondimensional hoop stress along the hole surface (R0 =5nm) 25 Figure 2.4 Nondimensional radial displacement along the hole surface 26 Figure 2.5 Variation of nondimensional hoop and radial stresses along the x i -direction for different Ks (R0= 5nm; d=0) 27 Figure 2.6 Variation of nondimensional hoop and radial stress components along the xj-direction for a hole with residual surface stress ( r ° ^ 0 ; d=0)........29 Figure 2.7 Variation of nondimensional interfacial hoop stress at 0 = 0 with loading ratio or inhomogeneity radius (solid lines for matrix and dash lines for inhomogeneity) 30 Figure 2.8 Nondimensional hoop stress along the inhomogeneity/matrix interface (solid lines for matrix and dash lines for inhomogeneity) 31 Figure 2.9 Variation of nondimensional interfacial hoop stress due to an eigenstrain E* with inhomogeneity radius (solid lines and left Y-axis for matrix, and dash lines and right Y-axis for inhomogeneity). 32 Figure 3.1 A nanoscale elliptical inhomogeneity in an infinite matrix 33 Figure 3.2 Conformal mapping from z-plane to £ -plane 36 Figure 3.3 Variation of stress concentration factor at 6 = 0 with hole size R and Ks (a/b=\.5) under uniaxial loading <jyy=cr0 44 Figure 3.4 Variation of stress concentration factor at 9 = 0 with hole size R and Ks (a/b=3) under uniaxial loading CT^=CT 0 46 Figure 3.5 Variation of nondimensional tangential stress along the hole surface with i?=6nm (alb=\.5) 47 Figure 3.6 Variation of nondimensional normal displacement along the hole surface for different Ks (a/b=\.5; i?=6nm) 48 v i Figure 3.7 Variation of stress components, CT22 and cr,,, along the xi-direction for different Ksunder uniaxial loading cr^ =tr0 (a/b=l.5; i? = 6nm)...49 Figure 3.8 Variation of stress components, a22 anda u , along the xi-direction for a hole with residual surface stress under uniaxial loading ay°y=aQ ( f l/ft=1.5;rV0) 50 Figure 3.9 Variation of the nondimensional tangential stress at 6 = 0 with the inhomogeneity size R and Ks (alb-\.5) under uniaxial loading (cr^, = <T0): solid lines for matrix and dash lines for inhomogeneity 51 Figure 3.10 Variation of the nondimensional tangential stress at 0 = 0 with the inhomogeneity size R and K (a/b=3) under uniaxial loading (<j™y =cr0): solid lines for matrix and dash lines for inhomogeneity 52 Figure 3.11 Variation of nondimensional tangential stress along the interface with i?=5nm and a/b=\.5: solid lines for matrix and dash lines for inhomogeneity 53 Figure 3.12 Variation of the nondimensional tangential stress at 6 = 0 on the interface with the inhomogeneity size R and Ks (a/b=l .5) under a uniform dilatational eigenstrain e* (solid lines for matrix and dash lines for inhomogeneity) 54 Figure 3.13 Varia t ion of nondimensional strains along the interface in an inhomogeneity for different values of Ks and R under a uniform dilatational eigenstrain 55 Figure 4.1 Nanoscale arbitrarily shaped inhomogeneities in a matrix material 57 Figure 4.2 Surface coordinates (n, t) and Cartesian coordinates (x, y) 59 Figure 4.3 A nanoscale elliptical hole in a finite plate 63 Figure 4.4 Comparison of analytical and finite element hoop stress solutions (dash lines for analytical solution of infinite plate and solid lines for finite plate with 2i?//=0.05 and 2R/h=0.05) 64 vu Figure 4.5 Variation of nondimensional hoop stress at 9 = 0 of a finite plate with R 65 Figure 4.6 Variation of nondimensional hoop stress at hole surface with 9 66 Figure 4.7 Variation of nondimensional hoop stress at 9 = 0 with 2R/1 (2R/h=0.\; R=5nm for Ks>0 and i?=30nm forKs < 0) 66 Figure 4.8 Variation of nondimensional tangential stress of an elliptical hole with R and 0 (a/b=\.5) 67 Figure 4.9 Variation of nondimensional hoop stress at 9 - 0 with 2a// (2b/h=0.\; i?=5nm for Ks > 0 and #=30nm forKs < 0) 68 Figure 4.10 A nanoscale circular hole in a semi-infinite plate 69 Figure 4.11 Variation of nondimensional hoop stress at 9 = ± nj2 with hole radius for two different D/R values (B/R=2 and 1.1) 70 Figure 4.12 Variation of nondimensional hoop stress on hole surface with angle 0 (D/R=2 and 1.1) 71 Figure 4.13 Variation of nondimensional hoop stress with D/R (R=5nm for Ks > 0 and i?=22nm for Ks < 0) 72 Figure 4.14 Two circular inhomogeneities in a finite plate 73 Figure 4.15 Variation of nondimensional interfacial hoop stress with radius and angle (D/R=2; Solid lines for matrix and dash lines for inhomogeneity) 74 Figure 4.16 Variation of nondimensional interfacial hoop stress with radius and angle (D/R=\.\; Solid lines for matrix and dash lines for inhomogeneity) 75 Figure 4.17 Variation of nondimensional hoop stress with D/R (7?=5nm for A^s > 0 and R=7nm for Ks < 0; Solid lines for 0 - n and dash lines forf9 = 0)...76 Figure 4.18 Variation of nondimensional hoop stress on hole surface at 9 = 0 with R (2#//=0.05 and 2R/h=0.05) 77 Vll l Figure 4.19 Variation of nondimensional hoop stress along hole surface with 6 (2R/l=0.05 and 2R/h=0.05) 77 Figure 4.20 Variation of nondimensional normal displacement with 0 (2i?//=0.05 and 2R/h=0.05; Solid lines for anisotropic case and dash lines for isotropic case) 78 ix List of Symbols a Length of semi-major axis of ellipse b Length of semi-minor axis of ellipse d Loading ratio E Young's modulus Ks Surface elastic modulus defined as 2/us + A s - r° Ro Radius or Curvature Radius R Radius or Average of lengths of semi-major and semi-minor axes of ellipse SaPyS Surface elastic constant tensor Displacement component in the /-direction U B Bulk elastic strain energy Surface elastic strain energy w Potential energy of external loads 7 Surface free energy Surface stress tensor under unstrained condition r° Residual surface stress under unstrained condition aa Stress component Lame constants Surface Lame constants ' MM Lame constants of matrix Lame constants of inhomogeneity V Poisson's ratio Poisson's ratio of matrix Poisson's ratio of inhomogeneity n Total potential energy X A c k n o w l e d g e m e n t s I wish to express my greatest appreciation and gratitude to my advisor, Professor Nimal Rajapakse, for suggesting the topic of this thesis, and for his patient guidance and continuous support which enabled me to continue this work and complete the thesis. I also gratefully acknowledge the support from Professor Nimal Rajapakse's research grant from the Natural Sciences and Engineering Research Council of Canada. I would also like to express my sincere thanks to the members of my thesis committee, Professor Mohamed Gadala and Professor Mu Chiao, for their kind help and for spending their time reading my thesis and make valuable comments. I also want to thank my parents, my sisters and my friends for their constant support, encouragement and help. xi Chapter 1 INTRODUCTION 1.1 Nanotechnology In 1959, Nobel Laureate physicist Richard Feynman gave a now famous lecture, "There is plenty of room at the bottom" [1]. He stimulated his audience with the possibility that in future one could manipulate and control, at the atomic level and create new materials, structures and devices. Nanotechnology is concerned with the development of such nanomaterials, nanostructures and nanodevices. Feynman's lecture laid the foundation stone to conceptualize and develop the research base for nanotechnology. As science and technology advance, there has been a great demand for new materials with better properties for applications in various advanced technology developments. As a result of several decades of research and development, nanoscale materials and structures can now be fabricated and instruments such as scanning tunneling microscopes and atomic force microscopes have been invented for nanostructure measurement, characterization and manipulation. In 1999, another Nobel Laureate chemist Richard Smalley presented his views on nanotechnology: "We are about to be able to build things that work on the smallest possible length scales, atom by atom with the ultimate level of finesse. These little nanothings, and the technology that assembles and manipulates them - nanotechnology - will revolutionize our industries, and our lives [2]". Nanoscience and nanotechnology are now attracting considerable interest and investment in order to develop revolutionary new applications in a wide range of disciplines. So what are exactly 'nano', nanoscale materials and structures, etc? 'Nano' means a billionth, and a nanometer (nm) is 10"9 m which is roughly four times the diameter of an individual atom. Ten nanometers is about 1,000 times smaller than the diameter of a human hair [3]. Nanostructured or nanoscale materials mean materials that have at least one of the overall dimensions in the nanometer range (about 1 nm to 100 nm), such as 1 nanoparticles, nanowires and nanofilms. They are intermediate in size between atomic/molecular and microscopic structures. By assembling these nanoscale materials and structures, it is possible to manufacture nanocomposites, nanodevices and nanosystems with desired properties and functions. What are the potential advantages of nanoscale materials and structures? How can nanoscience and nanotechnology affect our lives? Nanomaterials and nanodevices are based on fundamentally new molecular organization and they exhibit novel physical, chemical, and biological properties and phenomena. The assembly of the nanometer-scale components into large structures can produce nanocomposites, nanodevices and nanosystems with unique properties and functions. As shown in Fig. 1.1, the impact of nanoscience and nanotechnology is pervasive [3,4]. In information technology, nanoscale materials and structures have already found their way into information processing devices: faster sensors for signal acquisition and smaller size devices for denser information storage, etc. For medicine and patient care applications, it may be possible to produce new nanoparticle-based systems for drug delivery and nanosensors for early detection of diseases. In aeronautics and space exploration, potential applications include nanoinstrumentation for microspacecrafts and thermal barriers and wear-resistant nanostructured coatings. For environment and energy, nanotechnology can be used to develop semiconductor nanoparticles for efficient solar cells, reactive metal fine powders for groundwater decontamination, and nanoscale catalysis for increasing the efficiency of chemical reactions and combustion. In materials and manufacturing applications, nanostructured ceramics can be both harder and less brittle than currently available ceramics and nanoscale reinforcements in polymer composites can have greatly improved mechanical properties. Lighter, stronger and programmable materials, wear-resistant tires, tougher coatings, flame-retardant plastics, self-repairing materials and biologically inspired materials can be imagined. It is quite apparent that nanoscale materials and structures can provide mankind with significant benefits covering a wide range of applications. 2 Information Technology Medicine & Health Aeronautics & Space exploration Nanoscience & Nanotechnology Environment & Energy Materials & Manufacturing Other Applications Fig. 1.1 Potential benefits of nanoscience & nanotechnology. 1.2 Nanoscale Mechanics Before larger scale industrial application of nanoscale materials and nanostructures, it is necessary to have a comprehensive fundamental understanding of all aspects of their behaviour so that functionally and economically better designs can be achieved. Nanoscale mechanics, which deals with the study of mechanical properties and mechanical response of nanoscale materials and structures, is an important part of nanoscience and nanotechnology. Many applications of nanoscale materials and structures require knowledge of their deformation and stress field, strength, and fracture. For example, a carbon nanotube used as tip of a scanning probe microscope can be resistant to damage due to its flexibility [5]; Nanocomposites are expected to have higher stiffness and strength than classical composites due to nanoscale effects such as surface energy [3]. There are challenging issues in inelastic deformation and fracture of materials in the development of nanodevices, biopolymers, and hybrid bio-abio systems [6]. Different from continuum mechanics, which deals with bulk materials (usually applicable at length scales greater than 20 microns) and molecular dynamics which deal with behaviour at the molecular levels, nanoscale mechanics is concerned with objects with such characteristic lengths that neither the atomistic nor continuum models are applicable. In fact, nanoscale mechanics is now progressing in the direction of extending the concepts and methods of traditional continuum mechanics and bridging those with the effects at the molecular levels [7]. Such approaches are now known as multi-scale 3 modeling and involve consideration and combining of effects that occur at different length scales. As concepts in continuum mechanics are well known in the mechanics community, it is useful to summarize the key features of the molecular dynamics/atomistic models. Different from the continuum mechanics in which the material is considered as a continuum, the atomistic models deal with the motion of atoms and simulate the behaviour of a nanoscale object by a cluster of atoms. The motion of atoms is governed by Newton's second law in which the atomic forces are derived from the interatomic potentials. The key issue is to determine the effective interatomic potentials efficiently and accurately. Molecular dynamic (MD) models are currently limited in applications due to heavy computing requirements. An interesting class of problems in nanoscale mechanics deals with nanoscale inhomogeneities in materials. In practical terms, these include defects and reinforcements. Nanoparticles, nanowires, nanotubes, etc used in nanocomposites or quantum dots and wires used in semiconductor applications are good examples of nanoscale inhomogeneities [8]. It is well known that elastic field of quantum dots in a semiconductor device or nanoscale inhomogeneities in a composite can significantly influence the properties of a device or composite. Natural biological materials (e.g., nacre, bone), which are nanocomposites of soft proteins and brittle minerals, are very hard and tough, and by studying their behaviour at nanoscale it would be possible to develop bio-inspired materials in the laboratory [9] As in the case of classical inhomogeneity problems, the mechanical field and effective material properties need to be determined for nanoscale inhomogeneity problems. The classical concepts of continuum mechanics need to be modified to account for effects that exist at the nanoscale. A significant such effect is the surface effect. For nanoscale problems, the ratio of surface to volume is high, and the effect of surface energy cannot be neglected. Several questions arise here. How can the surface energy 4 effects be included in the modeling of the inhomogeneity problems? How can atomistic structure and properties be included in such models? Is the behaviour of a nanoscale inhomogeneity vastly different from a macroscale inhomogeneity? These are a few questions that motivated the present study and to conduct research to examine the mechanical field of nanoscale inhomogeneities. 1.3 Review of Surface Elasticity Model The surface/interface effects can be easily seen at the atomic scale, and this has been clearly pointed out and explained by Streitz et al. [10] and Dingreville et al. [11]. Due to different local environment, atoms at or near a free surface or interface have different equilibrium positions than do atoms in the bulk of a material. As a result, the energy of these atoms is, in general, different from that of the atoms in the bulk. The excess energy associated with the surface/interface atoms is called surface/interfacial free energy. The properties of a material, which are sensitive to the atomic positions or energies, are affected by the surface/interfacial free energy associated with these atoms at or near the surface/interface. For a material in which the number of atoms near the surface/interface is small compared to the total number of atoms, the surface/interface effects are insignificant and often can be neglected. For nanoscale structures, however, the surface/interface effects can be important due to the high ratio of surface/interface area to volume. Nanoscale materials may exhibit novel phenomena such as higher elastic modulus and mechanical strength compared to the conventional materials due to the surface/interface energy effects. The ratio of surface free energy y (J/m ) and Young's modulus E (J/m3), y/E, is dimensional (m) and points to some other inherent parameter of a material [12]. This intrinsic length scale is usually small, in the nanometer range or even smaller. When a material or a constituent of a composite has one characteristic length comparable to the intrinsic scale, the surface free energy can play an important role on the properties of the material, and thus the properties become size-dependent. Many experiments have reported such size-dependent behavior. For example, atomic 5 force microscopy is used to determine the Young's modulus, strength and toughness of nanorods and nanotubes [13] and elastic modulus of nanowires [14]. Unlike bulk material elements, the measured effective elastic properties are highly dependent on the size of the nanostructure. Similar behavior is observed for nanocomposites. Singh et al. [15] investigated the toughness of nanoparticle-reinforced composites with varied particle sizes. More recently, Cadek et al. [16] measured the tensile modulus of nanotube-reinforced polymer composites. The above experimental studies reported that the mechanical properties depend on the size of the reinforcements and the surface area per unit volume. To explain the size-dependent behavior at the nanoscale, several researchers have performed atomistic computer simulations. Sun and Zhang used a semi-continuum model [17] and Liang et al. used the embedded-atom-method inter-atomic potentials [18] to study elastic constants of plate-like nanomaterials and nanowires respectively, and reported the size-dependency of the elastic constants. Ji and Gao [9] explained the high toughness and strength of biological nano-composites through the virtual internal bond model and demonstrated the important role of embedded nanoscale minerals in such materials. Many simulation methods such as ab initio molecular dynamics, Monte Carlo simulations, classical empirical inter-atomic potentials have been developed to compute the properties of nanoscale materials and structures. These computer simulation methods, though very powerful, are often too complex and can model materials only in small scale as they need substantial computing power and are therefore not attractive in engineering applications. Extension of continuum concepts to nanoscale is considered an attractive option by engineers because of the complexities associated with the molecular dynamics/atomistic models and the associated high computing cost. The surface (or interface)/inter-phase energy based concepts have been introduced to extend the classical continuum theories. It is assumed that the two homogeneous phases are separated by a layer which is defined as an inter-phase or a dividing surface/interface. In the concept of 6 inter-phase, a system is considered to be made of three phases—the two bulk phases and an inter-phase. The inter-phase has a finite volume and the boundaries of the inter-phase are usually chosen to be at locations at which the properties of atoms are no longer varying significantly with position [11]. It is pointed out that, in the analytical or numerical analysis of boundary-value problems of composites, it is convenient to replace the inter-phase model by a proper surface/interface model [19]. In the concept of a surface/interface model, the two bulk phases are considered to be separated by a single dividing surface, and the surface is assumed to be a mathematical layer of zero-thickness. A surface stress tensor exists on the surface layer owing to the surface/interfacial free energy. The displacement and/or stress discontinuities are assumed at the interface. In the case of nanoscale materials and structures, the surface/interface stress model has become an attractive option because the surface/interface contribution to the property of the material can be accounted and the surface/interface stress can properly represent the effect of the surface/interfacial energy. Gibbs [20] first introduced the concepts of surface/interface energy and surface stress. In the Gibbs's formalism of thermodynamics [20], a quantity, y, called surface free energy, is defined to represent the reversible work per unit area due to creating a new surface. Gibbs [20] also pointed out that in the case of solids, there is a different type of quantity, called the surface stress that is associated with the reversible work per unit area needed to elastically stretch a pre-existing surface. From the thermodynamics of solid surfaces, the relationship between the surface stress and the surface free energy has been derived as [21,22] ' °ap =Y5ap+drld£aP . 0-1) where oap and eap denote the surface stress and strain, respectively, and 8tj is the Kronecker delta. It should be noted that the surface stress tensor is a 2D quantity in the tangent plane of the surface and the strain normal to the surface is excluded in Eq. (1.1). Thus, the Greek indices take the value of 1 or 2. In the right-hand side of Eq. (1.1), the first term has no explicit variation of the strain, but the second term indicates a variation 7 of the surface free energy with respect to surface strain. Nix and Gao [23] presented an atomistic interpretation of the interface stress and showed that Eq. (1.1) is an expression in the Eulerian frame of reference and the first term does not appear in the embedded Lagrangian coordinates. Eq. (1.1) can be rewritten by introducing a set of surface elastic constants. Miller and Shenoy [24] and Shenoy [25] suggested a linear constitutive equation of the following form: where r0ap is the surface stress when the bulk is unstrained, and SaPy5 is the fourth order surface elastic constant tensor. In general, there are at most six independent elastic constants for a crystal surface. The number of independent elastic constants can be reduced when a surface possesses geometric symmetry [25]. For an isotropic surface, Gurtin and Murdoch [26] and Gurtin et al. [27] proposed the following generic and simple expression for surface stress-strain relation. °* = A * + (*S + *°KSafi + 2(MS ~r°)sap, (1.3) where As and jus are the surface Lame constants and r° is the residual surface stress under unstrained conditions. The surface elastic constants are quite different from the bulk elastic constants and have to be known before applying the surface stress model. With the assumption of isotropy, Gurtin and Murdoch [28] computed some sample values of the surface moduli. By using molecular dynamics, Dingreville et al. [11] presented the calculation of the surface elastic constants with surface energy. A systematic study of the surface elastic constants has been performed by Shenoy [25] and their values can be calculated from atomistic simulations. It is found that the surface elastic constants tensor SaPyS need not be positive definite. Therefore, the quadratic from SapyS£ap£yS may be negative, but it does not violate the basic thermodynamic postulates. Shenoy [25] pointed out that a 8 surface region with a special atomic structure cannot exist without the bulk and the total energy (bulk + surface) still satisfies the positive definite condition. At present, the determination of surface elastic constants by experiments seems to be a challenging task and requires further efforts to address some fundamental challenges [25]. In general, the surface/interface properties are anisotropic and depend on the crystallographic direction of the surface/interface, and thus surfaces/interfaces in solids usually have anisotropic stresses [25,27,29]. However, as pointed out by Weissmuller and Cahn [29], "obtaining empirical database to determine fully the interface stress would be an enormous undertaking out of proportion to its usefulness". Moreover, one can assume that the interface has isotropic stress and it is still meaningful to use suitable averages of the interface stress [29,30]. Therefore, the isotropic assumption of surface/interface stress is considered acceptable in the study of surface/interface effects and to illustrate the key features of physical behaviour [31-33]. It is noted that the surface/interface stress model is only an idealization of a complex problem. In fact, the surface/interface region contains several layers of atoms which have different properties from that of atoms in the bulk. Therefore, the surface/interface stress model is valid if and only if the bulk material is much larger than the size of several layers of atoms [11]. If the bulk material contains only a small number of atoms which is comparable to that of the atoms at or near the surface/interface, the validity of macroscopic thermodynamic quantities such as surface free energy is questionable and the applicability of surface/interface stress model can be challenged [11]. The presence of surface/interface stresses results in a jump of traction from one bulk phase to the other and results in a non-classical stress boundary condition. This boundary condition and the surface stress-strain relations together with the equilibrium equations of bulk body form a coupled system of field equations. The displacement at the interface is assumed to be continuous or discontinuous in the literature. For the first kind 9 of deformation, the two homogeneous bulk phases are perfectly bonded, but for the second kind, slip or imperfect bonding exists at the interface. 1.4 Review of Classical and Modern Inhomogeneity Problems Elastic materials containing single or multiple holes, inhomogeneities and eigenstrain regions are commonly known as inhomogeneity problems in mechanics literature. Such problems have a wide range of practical applications including the determination of effective properties of composites, mechanical fields due to phase transformations in materials, stress concentration factors, etc. Inhomogeneity problems have been studied extensively in micromechanics for several decades since the celebrated work of Eshelby [34-36]. The books by Mura [37], Nemat-Nasser and Hori [38], and Marko and Preziosi [39] have provided compressive reviews of the micromechanics of heterogeneous materials. The above references deal exclusively with the classical inhomogeneity problems based on continuum mechanics. Nowadays, an inhomogeneity can be fabricated with dimensions of the order of 1-lOOnm by taking advantage of nanotechnology. Research on modern inhomogeneity problems is emerging and experimental, theoretical and computational tools are used to examine this class of problems. An interesting research topic in nanomechanics is the study of elastic field of nanoscale inhomogeneities in an elastic matrix. Due to high surface/interface to volume ratio of nanoscale inhomogeneities, it is necessary to incorporate the surface/interface effects which are neglected in the classical solution. The surface stress model has been employed to study modern inhomogeneity probelms. Yang [40] obtained closed-form solutions of the effective shear moduli and bulk modulus of composite with dilute nanocavities when considering only constant surface stress (or surface free energy), and found that the solutions are a function of the surface free energy, size of the nanocavity and the surface stress effect only influences the elastic field around the hole surface. By considering the surface stress without a residual stress, Duan et al. [33] examined the effective bulk and shear moduli of a solid containing nano-inhomogeneities and showed that the effective properties depend on the surface elastic 10 constants. Sharma et al. [31] presented a closed-form expression of the size-dependent elastic state of a spherical nano-inhomogeneity under far-field triaxial loading and a dilatational eigenstrain by applying the Gurtin-Murdoch model with both residual surface stress and strain-dependent surface stress. Sharma and Ganti [32] studied Eshelby's tensor of a nano-inhomogeneity undergoing a dilatational eigenstrain and pointed out that only inhomogeneities with a constant curvature admit a uniform elastic state. Duan et al. [41] extended the Eshelby formalism for a nanoscale spherical inhomogeneity subjected to an arbitrary uniform eigenstrain (not dilatational) by considering the surface/interface effects and demonstrated that the Eshelby and stress concentration tensors are, in general, not uniform inside the inhomogeneity but are position-dependent and size-dependent. While the present thesis was in the final completion, Gao et al. [42] reported the extension of the classical finite element method (FEM) by introducing surface elements to take into account the surface effect and then applied the method to investigate the interaction between two nanovoids and effective moduli of nanoporous materials. It is noted that all the above studies assumed that both bulk and surfaces/interfaces are elastically isotropic. 1.5 Scope of the Current Work Based on the above introduction and literature survey, it is found that nanoscale materials and structures have attracted significant interest due to their potential use in the development of materials with novel and improved properties. Inhomogeneity problems at the nanoscale have become an important research topic with significant impact on advanced materials development based on nanotechnology. Recent research has incorporated the surface/interface effects by applying the surface/interface stress model to study the mechanical properties and behavior of an inhomogeneity-matrix system. However, studies on practically useful inhomogeneity shapes (elliptic and spheroidal inhomogeneities), multiple inhomogeneities, non-symmetric loading, etc have not yet appeared in the literature. Therefore, this thesis tries to fill a part of this gap by focusing on the study of two-dimensional inhomogeneity problems at the nanoscale. The main objective of the thesis is to provide a comprehensive theoretical analysis of a two-11 dimensional elastic matrix containing a dilute distribution of nanoscale circular/elliptical inhomogeneities under arbitrary loading and investigate the elastic field of the inhomogeneity and matrix. A finite element scheme is also presented to analyze finite plate with an inhomogeneity and more complex geometries such as two interacting inhomogeneities problems. Chapter 2 presents the details of the derivation of an analytical solution for elastic field of a nanoscale circular inhomogeneity embedded in an infinite isotropic elastic matrix. Perfect bonding is assumed at the inhomogeneity-matrix interface. The Gurtin-Murdoch surface stress model is applied to take into account the surface stress effects. Muskhelishvili's complex potential approach is extended for the first time to solve this class of problems with a non-classical boundary condition. The closed-form solution derived in this thesis corresponds to a far-field uniform biaxial traction and a uniform eigenstrain in the inhomogeneity. Selected numerical results are presented for the elastic field around a circular hole and the inhomogeneity-matrix interface. Chapter 3 considers an inhomogeneity having the shape of an ellipse. Similar to the case of circular inhomogeneity, displacement continuity at the interface is assumed and the Gurtin-Murdoch model is used. By using the conformal mapping and the complex variable technique of Muskhelishvili, infinite series expressions of analytic potential functions are obtained approximately. Selected numerical results are then presented for the elastic field around an elliptical hole and an inhomogeneity-matrix interface. In Chapter 4, a finite element scheme is developed to study more complex inhomogeneity problems. The present approach, although developed independently of Gao. et al. [42], is conceptually identical. However, this thesis presents a more careful validation of the F E M scheme by using the analytical solution for elliptical inhomogeneity and examining other special features such as the instability of elastic field for negative value of surface elastic constants. So in Chapter 4, the finite element formulation is first presented to incorporate the surface/interface stress effect. The finite element solution is then verified by considering the stress field around a circular hole in a 12 large plate under tension and comparing the results with the theoretical solution derived in Chapter 2. The elastic field around a circular/elliptical hole in a finite plate is studied by using the finite element scheme by taking into account the surface stress effects. To study the interaction between inhomogeneities, the elastic field of a finite plate with two nanoscale circular inhomogeneities is also investigated. Finally, an example of anisotropic case is presented. Chapter 5 presents a summary of the thesis, major findings of the current study and suggestions for future work. 13 Chapter 2 ELASTIC FIELD OF AN INFINITE MATRIX WITH A NANOSCALE CIRCULAR INHOMOGENEITY 2.1 Problem Description and Basic Equations In this chapter, a two-dimensional problem of the elastic field solution for an infinite matrix containing nanoscale circular inhomogeneities at dilute distribution is considered. The surface effects need to be incorporated due to the nanoscale dimensions of the inhomogeneities. The inhomogeneities are assumed to be so far apart that their interaction can be neglected. Thus the problem is simplified to an infinite matrix containing a single nanoscale circular inhomogeneity. In the current study, the matrix is subjected to uniform far-field tractions CT ,^<7"2 and cr,"' and a uniform eigenstrain s* is prescribed in the inhomogeneity (Fig. 2.1). The matrix and inhomogeneity materials are assumed to be linearly elastic, homogeneous and isotropic with Lame constants X^, //M and \ , /ax, respectively. Note that the subscripts M and I are used to identify quantities associated with the matrix and inhomogeneity respectively. The inhomogeneity, with its center at the origin of the coordinate system, has radius R0 as shown in Fig. 2.1. 1 Fig. 2.1 A nanoscale circular inhomogeneity in an infinite matrix. 14 For plane problems of classical elasticity, Muskhelishvili [43] provided an analytical solution by using the complex potential functions method. Muskhelishvili's approach is used to study the current problem. The displacement and stress components in Cartesian and polar coordinates, (x\, xi, X?) and (r, 6 , X3) respectively, can be expressed in terms of two analytic functions <j>(z) and y/(z) as [43]: 2ju(ur+iue) = Q-w[K</>(z)-z^(z)-y/(z)\, (2.1) CTrr+oeg=2[$(z) + <t>(z)\, (2.2) ic r , =<f(z) + Jlz~)-e2e[zf(z) + y,Xz)], (2.3) <T - 1 C 2JX{U\ +'m2) = K(j)(z)-zcf) (z)-y/(z), (2.4) cTn+c722 = 2^(z) + r/'(z)j, (2.5) CT22 - c r n +2icr12 =2[zf(z) + y/\z)\, (2.6) where w, and ai} are displacement and stress components respectively, z = x, + ix 2 = re'9, K = 3 - 4v for plane stain and (3 - v) /(l + v) for plane stress, and ju and v are the shear modulus and Poisson's ratio, respectively. The overbar in Eqs. (2.1)-(2.6) represents the complex conjugate, and the prime denotes differentiation with respect to the variable z. 2.2 Formulation of Inhomogeneity Problem Assume that there is perfect bonding at the inhomogeneity-matrix interface, then the displacements are continuous at the interface: (Ur + ™8)M = ("r + iUd), + ("r + 'XU9)| » » t T = /?„ , (2.7) where the last term is the displacement induced by the prescribed uniform dilatational eigenstrain £*, i.e., e*u = e*22 = £*, and (ur + iue)* = RQ£\ r=Kn The surface traction on the inhomogeneity-matrix interface is discontinuous due to the surface stress effect. The inhomogeneity and matrix phases are assumed to be separated by a single dividing interface layer with zero thickness and surface stress exists 15 on this layer owing to the surface/interfacial free energy. This interface layer has different elastic properties from that of the bulk. For isotropic bulk and interface, the following field equations and constitutive relations can be established based on the theory proposed by Gurtin and Murdoch [26] and Gurtin et al. [27]: In the bulk (matrix and inhomogeneity): < . = 0 , a* =13^+2^, (2.8) On the surface/interface: fc«J+<,,=0,. (2.9) . t > J = < ^ > (2-10) < = ^ a + 2 ^ S - r 0 ) ^ a + ( A s + r ° ) ^ a , (2.11) where superscripts B and S are used to denote the quantities corresponding to bulk (matrix and inhomogeneity) and surface/interface; cr,y and £ y denote stress and strain, respectively; A and ju are the bulk Lame constants; Sy is the Kronecker delta; nj is the normal vector on the surface; Xs and /us are the surface Lame constants; r° is the residual surface stress under unstrained conditions; kap is the curvature tensor of the surface/interface and [J*]] = (*)M - (*), denotes the jump across the matrix-inhomogeneity interface. It should be noted that the surface stress tensor is a 2D quantity and the strain normal to the surface is excluded in Eq. (2.11). Thus, the Greek indices take the value of 1 or 2, while Latin subscripts adopt values from 1 to 3. In the bulk, the equilibrium equations and constitutive relation are the same as those in the classical elasticity. On the surface, the surface stress [Eq. (2.11)] is a sum of a residual stress and a linear function of surface strain, and the surface Lame constants are different from that of the bulk. Eq. (2.9) and (2.10) are the equilibrium equations in the tangential plane and normal direction of the interface respectively. The presence of surface stress results in a non-classical stress boundary condition. When no surface stress is considered, these equations reduce to the classical traction continuity equations. In the (r, 6, X3) coordinate [X3 is the direction perpendicular to the (r, #)-plane], Eqs. (2.9) and (2.10) can be written as: 16 On the surface/interface: In 9 -direction: In x3-direction: dcr00 do R0d9 dx3 8<J33 der, 03 dx3 RQd9 = 0, = 0. In r -direction: R, (2.12) (2.13) (2.14) For plane problems, cr®, = <J39 = <r|3 = 0 and the derivatives with respect to X3 are zero. Thus, Eq. (2.13) is automatically satisfied. Eqs. (2.12) and (2.14) can be expressed in the following complex variable form: .dtrl R0d9 (2.15) The left-hand side of Eq. (2.15) can be written in terms of potential functions by using Eq. (2.3). For the right-hand side, the surface stress is: a% = T0+2(jus-T0)eee+(As+T0)(£33 + egg). (2.16) The elastic strain sgg at the surface can be obtained from the following equations: 1 £ee-£rr+i£re= — [z<t>'(z) + V/'(z)F ae (2.17) (2.18) r=Rn Here, Q = X + p. for plane strain and Q = ju(3X + 2ju)/(X + 2/J.) for plane stress. Therefore, . (2.19) 1 2Q <f>(z) + *'(z)]+ -±-[zf{z)+ ¥{z)Yd + -L \zf{z)+ <p (z)\-2'e r=Hn Note that the strain egg is continuous at the interface because of the continuous displacement at the interface. In the following derivation, the strain ege is calculated from the matrix. Note that e33 =0 for plane strain and s33 =v(err + £gg)/(v for plane stress. Because of the discontinuity of the strain s33 at the interface for plane stress, the mean strain is used, i.e., 17 *33 = i [ ( % ) M • (2- 2 0) (^33)M c a n ^ e obtained by using Eq. (2.17). Due to the eigenstrain effect, (£33), is, ( % ) . = ^ l k + * J ; + 2 * ' ] , ' (2.21) where (£,.,. +£ee)i is the elastic strain in the inhomogeneity, which can be obtained from Eq. (2.17). From the above discussion, we have to obtain the complex potential functions which should satisfy the boundary conditions. To simplify the expression in the following derivation, introduce a non-dimensional complex variable £ such that: z = m(t) = R0£, S = r,ee, r0=r/R0. (2.22) Note that at the interfacer0 = 1, and <j)(^) - </)(m(^)) and y/(%) -y/(m(^)). The complex potentials ^ M ( £ ) , ^ ( £ ) , V M (£ ) a n c * corresponding to the matrix and inhomogeneity are now expanded into the following Laurent series form: 0M($) = A$ + f i A n r , vM(£) = B€ + f t B n r , (2.23) n=l «=1 A(£) = E ^ " . ^(^ = 1 ^ . (2.24) n=l n=l Note that the constant terms have been omitted in Eqs. (2.23) and (2.24) since they represent the rigid body displacements and have no effect on the stress distribution. In Eq. (2.23), A and B are real and complex numbers respectively characterizing the remote stress field. In view of Eqs. (2.5) and (2.6), A A ?R O O , OO ' ^ A 00 0 0 , / ^ ' CO r \ r \ C r i l + ° " 2 2 = — ' °"22 _ C r i l + 2 l 0 " l 2 (2.25) where cr,™ , cr™ and cr,™ are the far-field stresses which in the case of our study are assumed to be: (TJJ = o"0 , cr,™ = d<r0 , cr,™ = 0. (2.26) Here, d is a real number characterizing the loading ratio cr,™/cr™,. In the present case, both A and B are real numbers. Note that d=\ and d=0 refer to biaxial and uniaxial loadings, respectively, while d=-l represents pure shear loading. 18 The following set of equations is obtained by substituting Eqs. ( 2 . 1 ) , ( 2 . 2 3 ) and ( 2 . 2 4 ) into Eq. ( 2 . 7 ) and equating the coefficients of em < ?. 0 = - 2 F 2 , -^{KUAX-B)=-3F3-GX, MM • ^ ^ . = - ( " + 3 ) ^ 3 - ^ , , ' ( " = 1 , 2 , 3 . . . ) , ( 2 . 2 7 ) MM ^(KMA-A-B])=K1F]-F]+2MIR0S\ MM MM ^ ( ^ - 5 „ + 2 ) = / c , F „ + 2 , ( « = 1 , 2 , 3 . . . ) . ( 2 . 2 8 ) MM Substitution of Eqs. ( 2 . 3 ) , ( 2 . 1 6 ) , ( 2 . 2 3 ) and ( 2 . 2 4 ) into Eq. ( 2 . 1 5 ) yields: 2A + B} - ( F , + ^ ) = r ° +4A5/J}R0£' - A , 5 , - A , 5 , + 2 ( A , + A 3 ) A 2 ^ + A 4 A 5 ( F , +F}), 2B2 - 2F2 = A , 5 2 + 4 A 4 A 5 F 2 , - AX - B + 3 F 3 + G , = -AXB- 2AXA\ + 3A]B3 + ( A , + A 3 ) A 2 I , - 3 A 4 A 5 F 3 , -AL+B3-F3=ALB+2AXA\ - 3 A , 5 3 - ( A , + A 3 ) A 2 v 4 , + 3 A 4 A 5 F 3 , - I„ + l + (n + 3)FN+3 + GN+X = -Ax(n + \)(n + 2)I„+1 + A 2 ( A , + A 3 ){n + \)AN+X + A ,(n + 3)BN+3 -(n + 3)A4A5FN+3, (n=\, 2 , 3 . . . ) , - {n +1)AN+X + BN+3 - FN+3 = A , (n + \){n + 2)AN+X - A2 ( A , + A 3 )(n + \)AN+] -Ax(n + 3)Bn+3+(n + 3)A4A5FN+3, (n=\, 2 , 3 . . . ) , ( 2 . 2 9 ) where 4MMRO QM QI A 3 = 0 for plane strain and A j t £ — H M _ f o r plane stress, 4MMRo ^ M " 1 A 5 = 0 for plane strain and ^ + T — ^ — for plane stress. ( 2 . 3 0 ) 1 9 Note that all the parameters A , , A 3 and A 5 are defined as some quantity divided by the inhomogeneity radius R0. These quantities only depend on the material properties of the bulk and interface, and hence represent the intrinsic length scales. The solution of Eqs. (2.27)-(2.29) yields, A A - ^ 2 - c i + 3 K A 5 c , ) B ^ A + ] = 0 ; ( n = 1 ) 2 j 3 } j ( 2 3 1 ) d] + KMC2 +d2+ 3A 4 A 5 rf 3 A 7 r ° +2(2A 5 A 7 - c 7 ) / / , V +2(A 7 c 4 + A 6 A 8 c 7 ) ^ " = A 1 ' # 2 = U > A 6 c 7 + A 7 c , BiJc^c5-c3+3A4A5c6)B^ ^ 3 = 0 , ( ^ 1 , 2 , 3 . : . ) , (2.32) dx +KMC2 +d2 + 3A 4 A 5 c / 3 F = - A 6 r 0 + 2 A 6 ( c 1 A 8 - c 4 ) ^ - 2 ( 2 A 5 A 6 + c , ) / / , V f = Q 2 ( A 6 c 7 + A 7 c , ) ' 2 ' ^ 3 = ^ . ( " C ' + C 2 7 4 ^ 5 ^ ' ^ - = 0 , ("=1,2,3...), (2.33) /e, a, + KMC2 +d2+ 3 A 4 A 5 a 3 G , = - A 6 (- c}+c2-c3+ 3 A 4 A 5 c 6 )fl _ B] 3A 6 (- e,+c2 - c 4 - c5)ff y dx +KMC2 +d2 + 3A 4 A 5 c / 3 /r, rf, + /vMc 2 + rf2 + 3A 4 A 5 rf 3 G n + , =0, (w=l,2,3...), (2.34) where A 6 = ^ i / ^ M » A 7 = ( ^ i - 1 ) / 2 , A 8 = ( * r M - l ) / 2 , c , = l + 2A, , c 2 = A 6 ( l + 3A,), c3 = A 6 ( l - A , - A 6 ) / * r , , c 4 = A 2 ( A , + A 3 ) - l , c5 = A 6 [ l - A 2 ( A , + A 3 ) + 2A, + i r M A , ] , c 6 = A 6 ( A 6 - l ) / / c , , c 7 = l + A 4 A 5 , rf, =1 + 4A, +A 2 (A , + A 3 ) , d2 = A 6 [ l + A, + A 2 ( A , + A3) +KMA6]/^ , rf3=A6(l + / c M A 6 ) / / c I . (2.35) The complete elastic field is explicitly given by Eqs. (2.1)-(2.6) together with Eqs. (2.23), (2.24) and (2.30)-(2.35). After some manipulation, the stresses along the interface are: In the matrix: 2/4-.B, (5-35 3 )cos20 ( T e 0 = — + " » 20 2 A + B, (-4A, - 5 + 35 3) cos 20 CTRR= L + - ! 3- , Ro Ro (-2 A, + B.+ 3B3) sin 20 O"r0 = „ » U - i 0 ) ^0 In the inhomogeneity: _ 2F, | (12F3+G|)cos2.9 _ 2FX G,cos2r9 ^ i?0 /?0 rr RQ R0 _^ ________ ( 2 3 7 ) Note that due to the surface/interface stress effects, A , ( A 3 , A 5 ) * 0 and r° * 0 , and the stress state depends on the size of the inhomogeneity R0 as the coefficients A \, B\, 5 3 , F\, F), and G\ in Eqs. (2.31)-(2.34) are nonlinear functions of R0. When R0 is quite larger than the intrinsic scales, these parameters A , , A 3 and A 5 are very closed to zero, and the surface effects will only depend on constant surface stress r ° . When R0 is comparable to the intrinsic length, A , , A 3 and A 5 are large and the surface effects would be important. When no surface/interface stresses exist, i.e., A, = A 3 = A 5 = 0 and r° = 0 , the above results reduce to the classical solution in which the elastic state is size-independent. The shear and radial stresses are equal on either side of the interface in the classical case and, A 6 ( l + O ( l - r f ) c r o s i n 2 0 2(1 + K M A 6 ) _ = A 6 ( l + A g ) ( l + flr)o-0-4/v?' A 6 ( l + y M )(l-fl> 0 cos26> 2 ( A 6 + A 7 ) 2(l + * - M A 7 ) Hoop stress on the interface are different and given by, In the matrix: a = ( A 6 + 2 A 7 - A 6 A g ) ( l + <0o-0+4;/,g* | (4 -3A 6 +/c M A 6 ) ( l -^)cr 0 cos26? 2 ( A 6 + A 7 ) 2(1 + * M A 6 ) t K ' } In the inhomogeneity: _ = ( A 6 + A 6 A g ) ( l + rf)o-n-4/vr* ( A6(\ + KM){\-d)oocos20 2(A 6 + A 7 ) 2{\ + KMA,) 21 As can be seen from Eqs. (2.38)-(2.40), the stress state is independent of RQ in the classical solution. Note that the above solution reduces to the special case that the inhomogeneity is a hole by setting the shear modulus of the inhomogeneity //, to be zero. 2.3 Numerical Results for Elastic Field Around a Circular Hole Selected numerical results are presented in this and next section based on the surface elastic constants obtained from past studies. Experiments have been performed to determine the surface stress which has an order of 1 N/m [14,44,45,etc], but the surface elastic constants are difficult to be measured and no results are available at present. The embedded atom method was used by Miller and Shenoy [24] and Shenoy [25] to determine the surface elastic constants. Their results indicated that the surface elastic constants depend on the material type and the surface crystal orientation as shown in Table 2.1 for isotropic surface. Although the surface properties are generally anisotropic, it is assumed that isotropic case is sufficient to illustrate the main features of the size-dependent response. The plane strain case in which the surface effect is represented by the parameters A^ s and r° is investigated in the numerical study without loss of any generality. While the following numerical results are not to provide very accurate values, the study is to show the main behavior with surface effects. In the calculations, unless specified otherwise, KS = ± 1 0 N / m , (As + r ° ) = ±10N/m and r° is between - lN/m and lN/m. Table 2.1 Surface elastic constants. Units are N/m. Surface r° KS A s + r ° A l [100] -5.4251 3.4939 0.5689 -7.9253 4.0628 A l [111] -0.3760 6.8511 0.9108 5.1882 7.7619 N i [111] -0.6729 -1.8585 -0.1153 -3.0730 -1.9738 In this section, an infinite plane of aluminum containing a circular hole under far-field loading is considered. The bulk elastic constants for aluminum are: X_ = 58.17 GPa , MM =26.13 GPa [46]. 22 Based on the analysis presented in the previous section, hoop stress (plane strain) at 9 = 0 on hole surface is given by: Note that the first term is the classical elasticity result and the last two terms represent the surface stress effects and contain non-linear terms of R0. The first two terms are linear with respect to the loading magnitude while the third term is independent of external loading and linear with respect to the residual surface stress r ° . When d=\, i.e., under radially symmetric loading, the result is the same as that obtained by Sharma and Ganti [32]. Following the classical definition, a stress concentration factor can be defined for a hole by normalizing hoop stress at 0 = 0 by the remote loading magnitude when T° = 0. The effect of the surface elastic constant, Ks = 2/us +AS - r ° , is first studied by setting r° = 0. In this case based on Eq. (2.41), hoop stress concentration factor (<rgg /cr 0) is independent of the loading value <70 . Fig. 2.2 (a) shows the stress concentration factor for various values of far-field loading ratio d for a hole with radius = 5nm. Stress concentration factor varies linearly with d and is slightly increased or decreased compared to the classical result ( r ° = 0, Ks = 0) depending on whether A ' s is negative or positive. The influence of Ks appears to be more prominent when d is negative. Fig. 2.2 (b) shows the stress concentration factor for various values of the radius of the hole. The classical solution in which r° = 0 and Ks =0 is, as expected, independent of the radius, while the surface stress effects cause the stress concentration factor to be highly size-dependent especially when the radius is less than lOnm. The stress concentration factor increases or decreases rapidly when the cavity radius is less than 1 Onm depending on whether A^s is negative or positive. Surface stress effect is negligible when the hole radius is over 15nm and the stress concentration factor is equal to the classical elasticity solution. Similar behavior is observed for other values of d as well. crdg=(3-d)cr0- (2.41) 23 I i 1 ' 1 ' 1 • 1 -1.0 -0.5 0.0 0.5 1.0 d (a) Radius of the hole RQ = 5nm i • 1 i 1 • 1 • 1 • 1 0 10 20 30 40 50 R0(nm) (b) Loading ratio d=0 Fig. 2.2. Variation of stress concentration factor at 9 - 0 . Variation of nondimensional hoop stress [ ogg (R0 ,0)/<jo] along the hole surface is shown in Fig. 2.3 with hole radius of 5nm when r° = 0. The results are shown only from 9 = 0 to O-K/2 since the elastic field is symmetric with respect to both x and y coordinates. Fig. 2.3 (a) shows the result for the case of uniaxial loading in ^-direction (rf=0) when r° = 0. At 9 = 0, the nondimensional stress is reduced for positive KS and increased for negative KS compared to the classical result; Opposite behavior is 24 observed at 0 = 7r/2. Similar result can be seen for d--\ which represents pure shear loading [Fig. 2.3 (b)]. Under biaxial loading (d=\), the elastic filed is radially symmetric and the hoop stress does not vary along the hole. Positive or negative K5 decreases or increases the nondimensional hoop stress and the surface effects can be seen from Fig.2.2 (a) aXd=\. For other values of loading ratio d, the loading can be seen as a combination of shear and uniaxial loadings when d is negative or biaxial and uniaxial loadings when d is positive. 0.0 0.3 0.6 0.9 1.2 1.5 0.0 0.3 0.6 0.9 1.2 1.5 6>(rad) e(rad) (a) Loading ratio d=0 (b) Loading ratio d=-\ Fig. 2.3. Nondimensional hoop stress along the hole surface (RQ = 5nm). Fig. 2.4 shows the nondimensional radial displacement (normalized by cr0R0/juM ) along the hole with R0 =5nm when r° = 0. Under uniaxial loading in y-direction [Fig. 2.4 (a)], the nondimensional displacement is increased for positive Ks and decreased for negative Ks. When CT0 > 0, a positive Ks causes the hole to shrink while a negative Ks causes the hole to expand; opposite phenomena are observed when cr0 < 0. Under biaxial loading (d=l), the displacement does not vary along the hole and the surface effects are the similar to that under uniaxial loading. Under pure shear loading (d=-\), different behavior is observed as shown in Fig. 2.4 (b). When the classical displacement is positive, the displacement is decreased for positive Ks and increased for negative Ks. When the classical displacement is negative, an opposite surface effect is 25 seen for the same sign of Ks. As a result, no matter whether cr0 > 0 or <r0 < 0, a positive Ks maintains the circular shape of the hole while a negative Ks distort its shape. These behaviors imply that there exists local hardening or softening in the vicinity of the hole and this is consistent with other studies on the elastic moduli of nanoplate and nanobeam [24,47] and the surface stress effect on a solid containing a nano-scale spherical cavity [33,48]. 0.0 0.3 0.6 0.9 1.2 1.5 e(rad) (b) Loading ratio oN-1 Fig. 2.4 Nondimensional radial displacement along the hole surface. 26 Figure 2.5 shows the variation of nondimensional hoop stress and radial stress along positive xi-direction. The surface stress effect is evident near the hole surface but diminishes quite rapidly as xi increases especially in the case of hoop stress. Figure 2.5 (b) shows that when compared to the classical solution, the nondimensional radial stress is larger near the hole surface and slightly smaller far from the hole for positive values of Ks. Opposite behavior is noted for negative values of A^ s . (a) Hoop stress (b) Radial stress Fig. 2.5 Variation of nondimensional hoop and radial stresses along the xi-direction for different Ks (R0 = 5nm; d=0). Behaviour similar to that shown in Fig 2.2- 2.5 is observed for the case of plane stress. For plane stress, A 3 in Eq. (2.30) is nonzero and depends on the surface parameter (/Is +T°). In the numerical calculation, ± 1 0 N / m is used for this surface parameter. It is noted that the stresses are very close (the difference is less than 4%) to plane strain solutions with A^s equal to lON/m or -lON/m, and the difference in displacements is little larger (more than 4%). Consider next only the influence of residual stress r° (the surface stress is a constant value r°) on the elastic field of a plane containing a hole. With A, = A 3 = 0, the hoop stress at the hole surface is crgd = [(1 + d) + 2(1 - d) cos 26}j0 - r°/R0, (2.42) and the radial displacement at the hole surface is 27 r MM 2M] (2.43) M 8 2 Note that Eqs. (2.42) and (2.43) apply to both plane stress and plane strain cases. Two second terms on the right hand sides of Eqs. (2.42) and (2.43) are independent of the magnitude of remote loading and correspond to stress due to the residual surface stress. As cr0 would appear in the denominator of the second term on the right hand side of Eq. (2.42), a stress concentration factor of the classical form is not defined whenr 0 ^ 0 . It can be seen that the effect of r° is independent of the loading ratio d and the angle 0. It is clear from Eq. (2.42) that circumferential dependence of hoop stress along the hole surface is given by the classical elasticity solution minus a constant term that is linearly proportional to r° and inversely proportional to hole radius. The radial displacement is also linearly proportional to r° [Eq. (2.43)] and it is obvious that the residual stress r° can increase or decrease the radial displacement depending on the sign of r ° . Therefore, the residual stress also causes local hardening or softening around the hole and this is also observed for a solid with a spherical hole [40], To show the effect of r° on the stresses along jq-axis, let agg and <rfr denote hoop and radial stresses corresponding to the classical elasticity solution respectively and let cTg0 and o~fr denote hoop and radial stresses due to the residual surface stress. Fig. 2.6 shows the variation of crgg and ocn normalized by cr0 and GW and asrr normalized by r°/R0 along the positive xi-direction under uniaxial loading. Again, the residual surface stress shows significant influence on stress field in the vicinity of the hole surface. Its effect is negligible at a distant greater than four times the hole radius. Note that crgg and asrr along x\ are proportional to - — .o ( .. Y 2 _o f .. Y 2 and — respectively when Ks = 0. Hence, the nondimensional stress components due to r° shown in Fig.2.6 are independent of the radius Ro. 28 (a) Hoop stress (b) Radial stress Fig. 2 . 6 Variation of nondimensional hoop and radial stress components along the . xi-direction for a hole with residual surface stress ( r ° * 0 ; d=0). 2.4 Numerical Results for Elastic field of a Circular Inhomogeneity In recent years, quantum dot and wire structures have attracted considerable attention due to their potential application in nanotechnology [ 8 ] . It is also well known that mechanical and opto-electronic properties of nanocomposites are significantly influenced by the elastic field of the inhomogeneities. It is therefore important to understand the elastic field of a nanoscale inhomogeneity in a matrix material. To show the surface/interface effect on the elastic field, a matrix-inhomogeneity system made out of InAs/GaAs is considered in this section. The bulk Lame constants used are: A , = 5 0 . 6 6 G P a , / J , = 1 9 . 0 G P a for InAs, and ^ = 6 4 . 4 3 G P a , juM = 3 2 . 9 G P a for GaAs [ 4 9 ] . Consider first the case of a GaAs plane subjected to far-field loading with no eigenstrain in the InAs inhomogeneity. Hoop stress at the point 0 = 0 on the interface is investigated. Fig. 2 . 7 (a) shows the nondimensional interfacial hoop stress of the inhomogeneity and the matrix for various values of d. The results are similar to that of a circular hole [Fig. 2 . 2 (a)] and the effect of Ks is slightly more prominent in the inhomogeneity than in the matrix. When d=-l and Ks = l O N / m , the differences between the present and the classical results are 7 % for the inhomogeneity and 4 . 5 % for the matrix. Fig. 2 . 7 (b) shows the nondimensional hoop stress for different values of the 2 9 inhomogeneity radius with loading ratio d=0. Similar to the case of a circular hole, the surface stress effect is significant when the inhomogeneity radius is less than lOnm. Similar behavior can be observed for other values of d as well. 1 1 1 1 1 ' 1 ' 1 -1.0 -0.5 0.0 0.5 1.0 d (a) Inhomogeneity radius R0 = 5 nm x°=0, A:S=-10N/m — : T°=O, KS=O • -^-•:T°=0,A?S=10N/m A i • 1 • 1 • 1 • 1 ' 1 0 5 10 15 20 25 #0(nm) (b) Loading ratio a=0 Fig. 2.7 Variation of nondimensional interfacial hoop stress at 9 = 0 with loading ratio or inhomogeneity radius (solid lines for matrix and dash lines for inhomogeneity). Fig. 2.8 shows the nondimensional hoop stress along the inhomogeneity-matrix interface with R0 = 5nm. The surface stress effects are similar to the case of a circular hole. For ^ =0 and d=-l, a positive (negative) A^s decreases (increases) the stress at 6 - 0 30 and increases (decreases) the stress at 0 = x/2 in both matrix and inhomogeneity. Under biaxial loading (d=l), the stress is radially symmetric and increases or decreases depending on the sign of Ks. The solutions for displacements are also very similar to the case of a circular hole. 0.0 -0.3 -T°=0, i ^ - l O N / m -T°=0, KS=0 -T°=0, ^ l O N / m - * - T ° = 0 , A: s=-10N/m —T°=O , A : S = O - ^ - T ° = 0 , K S =10N/m 0.0 0.3 0.6 0.9 6 (rad) (a) Loading ratio d-0 1.5 1.5-1 1.0 0.5H © D 0.0 o -0.5--1.0--1.5 — T°=0, A^=-10N/m — x°=0, A^=0 <*-T°=0, A^=10N/m — - x ° = 0 , A ^ - l O N / m - - - T°=0, A^=0 • ^ - x 0 = 0 , A ^ l O N / m 0.0 0.3 1.2 1.5 0.6 0.9 9 (rad) (b) Loading ratio d=-\ Fig. 2.8 Nondimensional hoop stress along the inhomogeneity/matrix interface (solid lines for matrix and dash lines for inhomogeneity). The influence of eigenstrain e* in the inhomogeneity is now considered in the absence of far-field loading. In this case, the stress field is radially symmetric and there is 31 no shear stress. It is noted that Eshelby's tensor is uniform in this case but also size-dependent due to surface stress effect. Fig. 2.9 shows the nondimensional hoop stress <jgd jfJLx£* at the inhomogeneity-matrix interface. The nondimensional hoop stress is tensile in the matrix, while it is compressive in the inhomogeneity. In the classical case, hoop stress in the inhomogeneity and matrix have the same absolute value but opposite signs. Size-dependent behavior of hoop stress is clearly evident for an inhomogeneity with a radius smaller than 15nm. 2.6 2.54 U) 2.4 2.3 2.2 - T 1 • 1 a \ - a -a - T °=0 ,A : S =-10N /m T°=O,A:S=O -A -T°=0, A^=10N/m > v. / - a - - x °=0, /T^-lON/m - - - T°=0, KS=0 • i i - A - x°=0,A:S=10N/m i 1 1 1 1 • i i i i i . -2.1 -2.2 -2.3 1 -2.4 t> -2.5 -2.6 25 30 0 5 10 15 20 fl0(nm) Fig. 2.9 Variation of nondimensional interfacial hoop stress due to an eigenstrain s* with inhomogeneity radius (solid lines and left Y-axis for matrix, and dash lines and right Y -axis for inhomogeneity). Not presented here are the results for an inhomogeneity/matrix system subjected to both far-field loading and an eigenstrain in the inhomogeneity. However, it can be seen from Eqs. (2.31)-(2.34) that the elastic field components are linear functions of eigenstrain s* and loading ratio d (or A and B) and the superposition principle can be employed to obtain the solution. The effect of r° is not discussed here for an eigenstrain. It is clear from Eqs. (2.36) and (2.37) that the effect of r° is very similar to the case of a circular hole. In addition, the solution for plane stress is very close to that of plane strain for a prescribed eigenstrain. 32 Chapter 3 ELASTIC FIELD OF AN INFINITE MATRIX WITH A NANOSACLE ELLIPTICAL INHOMOGENEITY 3.1 Problem Description In this chapter, a two-dimensional problem of the elastic field solution for a matrix material containing nanoscale elliptical inhomogeneities at dilute distribution is considered. As in Chapter 2, the surface stress effects need to be incorporated and the interaction between inhomogeneities is neglected. This problem can be idealized as an infinite plane containing a single nanoscale elliptical inhomogeneity. The case of an elliptical inhomogeneity is more practically useful than a circular inclusion. On the other hand, the solution for an elliptical inhomogeneity cannot be expressed in closed-form and the salient feature of the elastic field can only be studied by numerical means. ° 2 2 t 1 „ ^21 I Fig. 3.1. A nanoscale elliptical inhomogeneity in an infinite matrix. In the current study, the plane is subjected to uniform far-field tractions cr,", cr"2 and cr," and a prescribed uniform eigenstrain E* in the inhomogeneity (Fig. 3.1), and the matrix and inhomogeneity materials are assumed to be linearly elastic, homogeneous and 33 isotropic with Lame constants , /JU and \ , , respectively. The matrix occupies a region denoted by SM, and the inhomogeneity, with the center at the origin of the coordinate system, occupies a region denoted by S\. The ellipse r represents the matrix-inhomogeneity interface. 3.2 Solution of Elastic Field The complex potential function method of Muskhelishvili [43] is employed to study the two-dimensional inhomogeneity problem under consideration. For plane problems, the displacement and stress components in Cartesian coordinates (x\, xj, X3) can be expressed in terms of two analytic functions <p(z) and y/(z) as shown in Eqs. (2.4)-(2.6). At the interface r , the boundary displacements and tractions can be written in the normal-tangential coordinates (n, t) as: 2M« B +i" , ) = l ^ ) - ^ ? f > ) - ^ ) ^ i a J (3-1) am-xam =jXz) + Jw-[zf(z) + y,Xz)}i2ia, (3.2) where t is the unit tangent, and n is the outward unit normal at the interface which in complex form is e i a (where a is the angle between the normal direction n and the positive x\-axis). Assume that there is perfect bonding at the inhomogeneity-matrix interface, then the displacements are continuous at the interface: k + l ^ L = (Un+™, ),+("„ 0 n r > (3-3) where the last term is the displacement induced by the prescribed uniform dilatational eigenstrain e*, i.e., £•*, =e*22 = s*, and: {un+iu,X =ze'e-'a, on T . (3.4) The surface traction on the inhomogeneity-matrix interface is discontinuous due to the surface stress effect. For isotropic bulk and interface, the field equations and constitutive relations based on the theory proposed by Gurtin and Murdoch [26] and Gurtin et al. [27] are given in Eqs. (2.8)-(2.11). In the («, t, x 3) coordinates [X3 is the 34 direction perpendicular to the (n, r)-plane], the equilibrium equations (2.9) and (2.10) can be written as: On the surface/interface r : In/-direction: |x„B, J + ^ - + - ^ - = 0 , (3.5) dt dx. In ^-direction: ] + ^ 2 L + ^ L = 0 , (3.6) dx3 dt Indirection: [ o r „ B J = ^ , (3.7) * o where R0 is the curvature radius. For plane problems, cr„B3 = cr3, = af3 = 0 and the derivatives with respect to x3 are zero. Thus Eq. (3.6) is automatically satisfied. Eqs. (3.5) and (3.7) can be expressed in the following complex variable form: k . - ^ , J = 4 + i ^ - 0-8) The left-hand side of Eq. (3.8) can be obtained from Eq. (3.2). For the right-hand side, the surface normal stress in the tangential direction is: al =T° + 2<7/ s -r°K,+(A s +r°)(*3 3+*„), on T . (3.9) Special attention is required when calculating the strain £33 at the interface, because the strain at either side of the interface can be different even though the displacement is continuous as assumed. Thus the interface has associated with it two interface stresses. Here, the average of the two interface stresses is taken as this is consistent with the case of a spherical inhomogeneity in which the interface stress is continuous and only one interface stress or the average of the two is used in Eq. (3.8) [32]. In the case of a hole, there is only one interface stress, or more exactly, the surface stress, and only this stress appears in the right-hand side of Eq. (3.8). Tangential elastic stain at the surface en can be obtained from the following equations: etl+em=£u+e22=^\z) + <p\z)), on T, (3.10) 35 EU-Em+2\sm={s22-en+2iEuya =-[z(f\Z) + y/\z)Ya, on T. (3.11) Here, Q = X + p for plane strain and p(3A + 2fx)j{X + 2p) for plane stress. Therefore, e„=j^{z) + f { 7 ^ ^ (3.12) In the matrix, the elastic strain sn is also the actual strain. In the inhomogeneity, due to the eigenstrain effect, {EII)1=(EII)\+ E* , where {£„)] is the elastic strain in the inhomogeneity which can be obtained form Eq. (3.12). For the other strain in Eq. (3.8), E 3 2 =0 for plane strain in both matrix and inhomogeneity, and for plane stress (E33 ) m = v M {snn + EU )EM / (v M -1) in the matrix and O33)i = v i [i£nn + £„ X+2£' ] / ( v i -1) i n t h e inhomogeneity. Here (em + e„ )eM and ( £ m +£„)* are the elastic strain in the matrix and inhomogeneity respectively which can be obtained by using Eq. (3.10). Following Muskhelishvili [43] and England [50], introduce the following mapping function to simplify the geometry of the problem (Fig. 3.2): z - plane L\ - plane Fig. 3.2. Conformal mapping from z-plane to t\ -plane. If] I v z = m(t) = R({ + -) , V ^ = -£ 1 1 + (L)2 1/2 = C + in = re" (3.13) 36 Here, a + b \-alb ri r f . . . R = , m = and / = va -b , (3.14) 2 \ + a/b where a and 6 are the length of semi-axes of the ellipse and 0 < m < 1. When m=0 the ellipse becomes a circle and in the limit m=\ it becomes a crack. The mapping function transforms region SM into the exterior region of the unit circle (|df| = 1) and region S i into an annular region between the unit circle and a circle of radius |£| = 4m . Here the region S i is imagined to be cut along the line L = {(x,,0): - / < x, < /} which is transformed into the circle with radius = 4m . Assume that there are no singularity points in the region S i , then ^,(z) and i//, (z) must be holomorphic in region S i and <r\(z) = <r\{z), ^,(z) = i / , (z) , zeL. (3.15) Consequently, the conditions (3.15) ensure that 0,(z) and y,(z) are analytic functions throughout the region S i . Now take ^(£) = 0(m(d;)) , and = y/(m(£,)) in the mapped plane . Therefore, the conditions (3.15) become: M4) = 6<Z), V,(£) = V.(I)> Vdj:|d;| = v ^ . (3.16) The complex potentials 0M(£), ^ , (£ ) , y M ( £ ) a n a " V1O5) corresponding to the matrix and inhomogeneity are now expanded into the following Laurent series form. <t>u(S) = AS + fjA„r , ¥u{$) = B$ + fjBnr, (3-17) <PX (£) = £ for + F„r ) , ^ , (£) = £ fee;" + )• (3-18) n=\ n=l Note that the constant terms have been omitted in Eqs. (3.17) and (3.18) since they represent the rigid body displacements and have no effect on the stress distribution. A and B are given constants characterizing the remote stress field. In view of Eqs. (2.5) and (2.6), 37 CO 00 c r u +CT 2 2 = 2(A + A) CO 00 '22 w 11 crn +2icr12 = — , (3.19) R - „ .* R where cr," , cr"2 and c,™ are the far-field stresses. . According to England [50], the imaginary part of .4 is related to the rotation cox at infinity and, coa=(\ + K)lm(A)/2ju. (3.20) In the current problem, the rotation is zero, so that: A = A. (3.21) Following Stagni [51] and Shen et al. [52], introduce an auxiliary functions Q(£) and a new auxiliary function ©(£) such that, m (£) (3.22) ©(£) = Q ( £ ) m ( | ) - 0" (£)[m(l/1)J. (3.23) The auxiliary functions QM(%), © M ( £ ) , Q,(£)and ©,(£) corresponding to the matrix and inhomogeneity can also be expanded into Laurent series as: R. (3.24) (3.25) n.K) = A.+i;(^"+^r). ©.(£) = o 0 +X(o„r+^-") Eq. (3.24) combined with Eqs. (3.17), (3.22) and (3.23) leads to the following relations: C = mA + B, C 0 = 0 , (3.26) D = 0, D0=B, £>,=0 , D2=-Cm-C_+A + mA\, D3 = -2C2 + 2mA2, Dn+3 = -(n + 2)C„ + 2 + mnCn + m(n + 2)An+2 -nA„, (n=\, 2, 3...). (3.27) Eq. (3.25) combined with Eqs. (3.18) and (3.23) gives the following relations between the coefficients On, Pn and En, F„, Ln and Mn. 38 O0 = -3mL3 + Lx + 3E3 - mEx, 0= -{n + 3)mLn+3 +(n + 1)Z„+1 +(n + 3)En+3 - (n + \)mEn+x, {n=\, 2, 3...), (3.28) PX = -2mL2 + 2E2, P2 = -mlx -Mx+Ex + mFx, P3 = -2M2 + 2mF2, PN+3 = nmMn -(n + 2)Mn+2 -nFn + (n + 2)mFn+2, (n=\,2,3...). (3.29) Furthermore, the conditions (3.16) combined with Eqs. (3.18), (3.22) and (3.25) yield the following relations. Fn=m"En, Hn=m"Gn, Mn = m"LN +nm"-\\-m2)En. (3.30) Using Eqs. (3.1), (3.4) and (3.22), the continuity of displacement across the interface, Eq. (3.3), can be expressed as: ~ k^M (£) - " (^1)] = ^ - VA - ^ Ml)]+ * • (3-31) Noting that ^ = ei0 on the interface, substitute Eqs. (3.17), (3.18), (3.24) and (3.25) into Eq. (3.31) and equate the coefficients of ew9 to the following relations. LQ=0, A3(KMA-Cx)- KXEX -Mx -\-2fj.xRs*, A3(-CN+X)= >cxEn+x-Mn+X, (n=\, 2,3,...), (3.32) A3(icMAX -C)= KXFx -Lx + 2juxRme*, A 3 K M 4 + I = K , ^ , + I - A . + I » («=1,2,3, . . . ) . (3.33) where A 3 =//,///M . To solve the present problem, it is required to obtain the coefficients including AN, C „ , D N , E„, Ln, ON and P N . Once these are known, the other coefficients ( B N , Fn, HN and Mn) can be determined. Based on the above analysis, DN can be expressed in terms of AN and C„ from Eq. (3.27); En and Ln can be expressed in terms of AN and CN from Eqs. (3.30), (3.32) and (3.33); ON and PN can be expressed in terms of En and Ln from Eqs. (3.28)-(3.30), and thus in terms of AN and CN. The only unknown coefficients are therefore AN and CN. Additional relations between the unknowns can be obtained from the boundary condition (3.8) as follows. 39 In the £ -plane, following England [50], e2i« = p i « » l (£) . I = e i e on T , (3.34) and the derivatives with respect to the tangential direction t can be expressed as: d _d_dz_ _ J L J _ & _ ^ \ _ ^ _ Yt~~dz~dt a i a 7 _ d j dz_dt d£ dj^dt' — = ie 5r •a/ = -ie (3.35) (3.36) In the £ -plane, Eq. (3.2) becomes: r ... LV m (£) w (£) •e2l(?, o n f . Multiplying the above expression by the (non-vanishing) factor m (£) , and eliminating y/\4) by using Eq. (3.22), yields, _ (3.37) m ^ ) K „ - i c T j = ^ ( ^ ) - e 2 i e Q ^ ) . The right-hand side of Eq. (3.8) is written as: <rs„ , ^ | [ ( cT , f ) M + (^) 1 ] | . a ^ ^ ^ + K ^ ) , ] R0 dt dt (3,38) Eq. (3.38) are derived in Appendix A by substituting the interface stresses (3.9) into the above equation and multiplying the resulting expression by a factor Consequently, by combining Eqs. (3.17), (3.18), (3.24), (3.25), (3.37) and (A.5), Eq. (3.8) becomes: \(n+])d X n(Cn -En-Mn + A - Ce 2 i* - £ n(An -Fn-Ln Y 1 1 m<£)lR m(4)/R\ j\go +K + +h])e~i0 +(c?2 +h2)eW +(c?3 + h3)e2W + (ft + K )e™ + (Ss + "s)*-3'9 + (ft + K)*™ + ^ + ***** 40 CO 00 00 00 + f X e i ( " + 3 ) ( ? + i X e - i ( " + 3 ) * +__U„Q'{"+4)0 +J_Vnc-'{"+4)e (3.39) n=\ «=1 n=\ n=\ J Here, the coefficients gQ , gn («=1,2 6), //„ («=1,2,...,8), 5 „ , 7/„ , <7„ and Vn are defined in Appendix A. It can be seen that these coefficients are functions of the nondimensional parameters A, and A 4 or A 5 and A 6 which are defined as the ratios of the intrinsic lengths to the characteristic length of the material (the inhomogeneity size), and therefore, these parameters and the residual stress r° represent the surface effects. Next, employing a method similar to that used by Shen et al. [53], (£)/R can be expanded into an infinite series of the following form (see Appendix C). 1 1 -(l + Z>*sin2#)"3/2, b* = 4m (3.40) k=\ J - l k=\ Eq. (3.41) can be rewritten as: (1 + Z>*sin2f?r3/2 ~ where ~ J . X-1 T S„i2ke , -\2k6\ , j e i 2 J g + e - i 2 J g _ ^ e i 2 ( J - l ) g + e . 2 ( W ) g ' (\-me2,e)(\-me-2'9) 1 ( l - m e 2 , 0 ) ( l - w e - 2 , e ) J r , V" 1 r rA2k6 . -\2k6\ fo+2-f2k(Q + e ) k=\ (3.41) (3.42) f0=(\ + m2)l0-2ml 2 ' (3.43) f2k = (\ + m2)I2k-ml2k_2-mI2k+2, (fc=l, 2,..., J-l), flJ ~ ^2J ~ mh.J-2 • Take the denominator ( l -me 2 ' e ) ( l -we" 2 ' e ) into the first term of the right-hand 1 side of Eq. (3.39), i.e. (see Appendix D): , and expand it into power series of the following form 41 CO 00 m'tf)/R (\-me2'd)(\-me-2ie) (\-m2)2 T ^m n e i 2 " 9 +l + ^ ( m n + ( l - m 2 ) » m " ) e (3.44) Thereafter, using Eqs. (3.40), (3.42) and (3.44), Eq. (3.39) becomes ±n(Cn -En -MBy™ +A-Ce2>°-±n{An -Fn-Ln>* x /o+Z/2t(e i 2 M+e" i 2 W) to+A+(ft+^)e- i e+(g 2+^)e i e+(g 3+^)e k=\ + (ft + K )e2i" + (ft + h5 )e"3ie + (g 6 + h6 )e3i* + V ~ 4 i * + /\e oo oo + 1 v ( " + 3 ) *+X -i(n+4)0 (3.45) By equating the coefficients of c'"6 in Eq. (3.45), the relationships between the unknown coefficients can be obtained. This yields a sufficient number of equations to solve for the unknown coefficients. Depending on the level of accuracy required, different values of J and the number of the coefficients in the power series, i.e., n, are chosen. 3.3 Numerical Results for Elastic Field Around an Elliptical Hole Selected numerical results for plane strain case are presented in this and ensuing sections. The surface/interface effects are represented by the residual surface stress, r ° , and parameter, J^ s (or r ° , A, and A 6 ) , and the values of these parameters for numerical calculation are same as in Chapter 2, i.e., Ks - ± 1 0 N / m and r° is between -1 N/m and 1 N/m. The number of the terms in the infinite series representation of the complex potential functions is chosen so that the error in the numerical calculation is maintained below 1%. This is achieved by increasing the number of terms in the series representation until the difference between two consecutive sums is less than 1%. Accuracy of the numerical calculations is checked by setting the surface elastic constants and residual 42 surface stress to negligibly small values and comparing the resulting numerical solutions with the classical elasticity solution for an elliptic hole in an infinite plane subjected to remote uni-axial tension. It is found that the two solutions agree very closely. In this section, an infinite plane of aluminum containing an elliptical hole under far-field loading is considered. The effect of the surface elastic constant, Ks = 2jus + ?L - r°, is first studied by setting r° = 0. In this case, the stress concentration factor, <j„ /cr 0 , at 6 = 0 is independent of the magnitude of the applied loading. Fig. 3.3 shows the stress concentration factor for various hole sizes R [= (a + b) 12] under uniaxial loading with alb =1.5. As expected the classical solution in which r° =0 and Ks = 0 is independent of the hole size, while the surface stress effects cause the stress concentration factor to be size-dependent especially when R is less than 20nm. The stress concentration factor increases or decreases when Ks < 0 or Ks > 0 . The difference between the classical and current results can reach 27% when R = lnmfor A^s >0. The surface stress effects are negligible when R is over 40nm and the stress concentration factor is equal to the classical elasticity solution. This behavior is similar to that of the circular case and of the spherical case [31]. However, for Ks <0, the stress is found to become highly oscillatory and singular at some values for R less than 6nm. It is noted that this unstable phenomenon does not happen for A" s > 0 in which case the stress becomes smaller as R decreases. 43 4.4 4.2 4.0 3.8 -€ 3.6 3.4 3.2 3.0 2.8 : -T°=0,A:S=-10N/m - x°=0, A^=0 -T°=0, /CS=10N/m - i 1 r-0 10 20 30 40 50 60 70 80 90 100 fl(nm) Fig. 3 . 3 . Variation of stress concentration factor at 6 - 0 with hole size R and Ks (alb-\.5) under uniaxial loading cr ,^ = cr0. It is proposed to further investigate this phenomenon by considering the case of a circular hole under uniaxial or biaxial loading. Closed form of analytical solution for the plane strain case is: a o a l + d | l ~ d c o z 2 0 | ( A , A 2 - l X l + ^ ) + ^ 0 / ^ o f ( A , A 2 - l X l - ^ ) c o „ 2 g ; 'o 1 + 2A 1 + 4 A , + A , A 2 <jrr_\ + d ^{AiA2-\\\ + d)+T°/cr0R0 i {\A2+2Aj\-d)cQz2e^ 1 + 2 A , 1 + 4 A , + A , A 2 ( 3 . 4 6 ) ( 3 . 4 7 ) ur (KM+\){l + d) ^{AlA2-\\l + d)+T0/cr0R0 (l-d) cr0RjjuM 8 2 ( 1 + 2 A , ) i ( y M + l X l + 2 A , X l - r f ) cos 20 1 + 4 A , + A , A 2 ± ( A , A , - I X I - < 0 „ „ cos26>- 4 V 1 2 — ^ ^cos20, 1+ 4A, + A , A 2 ( 3 . 4 8 ) °ORO/MM 4 1 - ' , . s i n 2 g - ^ 1 - y " » - ' ' X U 2 A l ) s . n 2 g _ i ( A l A ; - l X l - r f ) s i n 2 , i 1 + 4 A , + A , A 2 1 + 4 A , + A , A 2 ( 3 . 4 9 ) 4 4 Here, cr0 is the loading stress and R0 is the radius of the circular hole and d = 0 for uniaxial loading and d- \ for biaxial loading. When AT s <0, A , (A, = Ks/4juMR0 for circular case, see Chapter 2) is negative. Therefore, the denominators containing A, in Eqs. (3.46)-(3.49) can be zero for certain radii resulting in singular stress and displacement fields. The corresponding radii are very small, less than lnm, for the values of Ks corresponding to aluminum. A similar phenomenon exist for the elliptical case, however, there appear to be many hole sizes R which induce singular stresses. Some of these R values can be relatively large (>5nm) depending on the geometry of the ellipse, surface elastic constant Ks and elastic properties of the matrix material. Such unstable behavior should not occur in experiments involving real materials because of the plastic properties of nanoscale materials. This implies that Gurtin-Murdoch model has some restrictions when applied to solve situations involving negative Ks values. Fig. 3.4 shows the stress concentration factor for various hole sizes R with alb =3 under uniaxial loading. The result is similar to the case of alb =1.5. However, the effect of the surface elastic constants Ks is more noticeable. The stress concentration factor shows high size-dependency when R is less than 40nm. The difference between the classical and the current results can reach 40% when R = 1 nm for Ks > 0. The stress for Ks < 0 is unstable below 15nm. It is noted that the results are similar for other values of alb, and as the value of alb increases, the effect of the surface elastic constant Ks becomes more pronounced and the value of R below which the stress is unstable when Ks < 0 increases. 45 8 0 10 20 30 40 50 60 70 80 90 100 «(nm) Fig. 3.4. Variation of stress concentration factor at 9 = 0 with hole size R and KS (a/b=3) under uniaxial loading o™y = a0. Fig. 3.5 shows the nondimensional tangential stress [<J„(R,9)/<rQ] along the hole with A*=6nm and a/b=\.5. The result is shown from 9 = 0 to 9 = n/2 due to symmetry. Under biaxial loading as shown in Fig. 3.5 (a), the nondimensional tangential stress is decreased for positive Ks and increased for negative Ks for the whole surface. Different behavior is observed under uniaxial loading as shown in Fig. 3.5 (b). At 9 = 0 , the nondimensional tangential stress is reduced for positive A^s and increased for negative Ks, but opposite behavior can be seen at 9 = 71/2. Similar result as uniaxial loading can be observed for shear loading. It is noted that these results are similar to the case of a circular hole. 46 0(rad) 0 (rad) (a) Under biaxial loading cr™ = cr™ = <r0 (b) Under uniaxial laoding cryy=cr0 Fig. 3.5. Variation of nondimensional tangential stress along the hole surface withi?=6nm (a/b=\.5). Fig. 3.6 shows the nondimensional normal displacement (normalized by cvQR/juM ) along the hole surface for a/b=l.5 and Z?=10nm under biaxial and uniaxial loading. It can be seen that positive values of surface elastic constant Ks cause the hole to shrink when compared to the classical case (no surface stress effect) while a negative Ks causes the hole to expand more. Similar behavior can be observed in the case of uniaxial loading. Under pure shear loading (figure not shown), the absolute value of normal displacement is decreased or increased when Ks > 0 or Ks < 0. This is similar to the case of a circular hole. As a result, when Ks is positive the shape of the hole is relatively unchanged while for a negative Ks more shear distortion of the hole takes place. Therefore, the surface elastic constant (Ks) causes local hardening or softening around the hole. A larger positive Ks introduces a larger hardening zone, while smaller negative value of Ks (larger absolute value) introduces a larger softening zone. 47 0.0 0.3 0.6 0.9 1.2 1.5 0(rad) (a) Under biaxial loading (<j_ = cr" = <r0) 0(rad) (b) Under uniaxial loading (cr^, = cr0) Fig. 3.6. Variation of nondimensional normal displacement along the hole surface for different Ks (a/b=\.5; R=6nm). Fig.3.7 shows the variation of stresses cr22 and an along the positive xi-direction. The surface stress effect is significant only near the hole surface but diminishes quite rapidly as x\ increases especially in the case of a22 . Fig. 3.7 (b) shows that, when compared to the classical solution, the nondimensional stress cru/cr0 is higher near the 48 hole surface but slightly smaller far from the hole for positive values of Ks. Opposite behavior is noted for negative values of Ks. The behavior is similar to the case of a circular hole (Chapter 2). (a) Variation of stress cr22 (b) Variation of stress cr,, Fig. 3.7. Variation of stress components, cr22 and cr,,, along the xi-direction for different Ks under uniaxial loading ay=a0 (alb=\.S; R = 6nm). It is noted that similar phenomena are observed for the case of plane stress and other values of alb. For plane stress, A 4 in Eq. (A.2) is nonzero and depends on the surface parameter (A s + r ° ) . When ± 1 0 N / m are used in the numerical calculation for this surface parameter, it can be seen that the results of stresses and displacements for plane stress are very close to that of plane strain for the same parameter K5 (10 N/m or -10 N/m) and the same value of alb. Consider next the influence of residual stress r° on the elastic field of a plane containing an elliptical hole by setting Ks - 0 (or A, =0). As expected, the effect of r° is similar to that in the case of circular hole as r° bahaviors like a loading and is independent of the applied remote loading. The numerical solutions for stress do not show any instability problem with respect to r° . This behaviour can be easily confirmed 49 from Eq. (2.42) or Eq. (3.46) for the circular case as the denominators cannot be equal to zero. To show the effect of residual stress, r° , let cr2c2 and denote the stress components corresponding to the classical elasticity solution respectively, and let cr22 and of, denote the stress components due to the residual surface stress r ° . Fig. 3.8 shows the variation of a_\ and cr,^ normalized by cr0 and cr22 and erf, normalized by r°/R along the positive xpdirection under uniaxial loading when alb=\.5. The residual surface stress shows a significant influence on stress field only in the vicinity of the hole surface. Its effect is negligible at a distant greater than four times the major semi-axis. Note that nondimensional stress components due to r° shown in Fig. 3.8 are independent of R. (a) Variation of stress <r22 (b) Variation of stress cr,, Fig. 3.8. Variation of stress components, cr22 and cr,,, along the xi-direction for a hole with residual surface stress under uniaxial loading cr^ = cr0 (a/6=1.5;r° * 0). 3.4 Numerical Results for Elastic Field of an Elliptical Inhomogeneity This section presents results for elastic field of an infinite plane containing a nanoscale elliptical inhomogeneity in the case of plane strain. To show the 50 surface/interface effect on the elastic field, a matrix-inhomogeneity system made out of InAs/GaAs is considered. Consider first the case of a GaAs plane subjected to far-field loading with no eigenstrain in the InAs inhomogeneity. The tangential stress o~„ at the point 6 = 0 on the interface is investigated. Figs. 3.9 and 3.10 show the nondimensional tangential stress for various values of R with a/b=\.5 and a/b=3, respectively. The nondimensional tangential stresses in both the inhomogeneity and matrix are increased and decreased when Ks < 0 and Ks > 0 respectively. The effect of Ks is slightly more prominent in the inhomogeneity than in the matrix. For example, when /?=lnm and Ks =10N/m in the case of a/b=\.5, the differences between the present and the classical results are 25.5% for the inhomogeneity and 17.5% for the matrix. Similar to the case of an elliptical hole (Figs. 3.3 and 3.4), the effect of Ks is more pronounced as the value of alb increases. Stresses in both the matrix and inhomogeneity for Ks < 0 become unstable below a certain value of R depending on the value of alb. 1.6 r o D ° 1.0 - •• - : T°=0, A ^ - l O N / m , - - - :T°=0,A: s=0 ^ - : T°=(), /Ts=10N/m • 0.6 0 10 20 30 40 50 fi(nm) Fig. 3.9. Variation of the nondimensional tangential stress at 0 = 0 with the inhomogeneity size R and Ks (a/b=\.5) under uniaxial loading (ayy = cr0): solid lines for matrix and dash lines for inhomogeneity. 51 1.6-1.4-- » -: x°=0, A^=- lON/m 1.2- — , — : T°=O,A:S=O 1.0-- A - , - A - : x°=0, K M o N / m 0.8-0.6-i. • i • i • i • i • i • i • i 1 i 1 i • 0 10 20 30 40 50 60 70 80 90 100 R0(nm) Fig. 3.10. Variation of the nondimensional tangential stress at 6 = 0 with the inhomogeneity size R and Ks (a/b=3) under uniaxial loading (cr^ = cr0): solid lines for matrix and dash lines for inhomogeneity. Fig. 3.11 shows the nondimensional tangential stress along the inhomogeneity-matrix interface when A=5nm and a/b=\.5. Under biaxial loading, the nondimensional stresses in both matrix and inhomogeneity are increased or decreased when Ks < 0 or A^s > 0. Under uniaxial loading, however, the nondimensional stresses in both matrix and inhomogeneity are decreased at 6 = 0 and increased at 0 = n/2 whenA^s >0; however, opposite behavior is observed for Ks < 0. The behavior of the tangential stress under pure shear loading is similar to the case of uniaxial loading. Therefore, the surface effect is the same as the cases of an elliptical hole and of a circular inhomogeneity. The results of displacements, which are not shown here, is also very similar to the case of an elliptical hole. 52 1.5-1.2-x°=0, A^ =-10N/m . : x °=o, ^=0 0.9-0.6-\ -:T°=0,A^=10N/m 0.3-0.0--0.3- 1 i 1 i i > i • i e(Rad) (b) Under uniaxial loading <j™y = <r0 Fig. 3.11. Variation of nondimensional tangential stress along the interface with i?=5nm and a/b=\.5: solid lines for matrix and dash lines for inhomogeneity. The influence of eigenstrain, e*, in the inhomogeneity is now considered in the absence of far-field loading. Fig. 3.12 shows the nondimensional tangential stress cr,,/ Hi8* o n m e interface between the inhomogeneity and the matrix for various values of R when a/b=l.5. The results are similar to Fig. 3.9. The size-dependent 53 behaviour of the stress field is evident for an inhomogeneity with R smaller than 15nm. Again, the effect of Ks is slightly more pronounced in the inhomogeneity than in the matrix and stress is unstable for Ks < 0 when R is below 3nm. 10 20 30 40 3.3 ' c o 3.2 2.9 -1.8 " co ^T-2.0 D -2.2 -2.4 — i — 10 i — 20 30 40 tf(nm) 50 V —•—T°=0,.R:s=-10N/m T°=0, ^=0 —A-x°=0,A : s =10N/m - • -T°=0,A^=-10N /m T°=0, A^=0 - A - t0=0, Ks=WWm 50 Fig. 3.12. Variation of the nondimensional tangential stress at 0 = 0 on the interface with the inhomogeneity size R and Ks (a/b=\.5) under a uniform dilatational eigenstrain e* (solid lines for matrix and dash lines for inhomogeneity). N o w consider Eshelby's problem in the presence of surface/interface effects. Eshelby [34] obtained a uniform elastic field for an ellipsoidal inhomogeneity under a uniform eigenstrain in the classical case. As shown in Fig. 3.13, the nondimensional normal strains in the xi-direction and X2-direction are no longer uniform and vary along the interface (in the inhomogeneity, £ = £e+e', where e is the actual strain and se is the elastic strain). As the inhomogeneity becomes smaller (e.g., /?=10nm), the effect of Ks and the non-uniformity of the elastic field are more prominent. The non-uniformity of strain field is more obvious as alb increases and more surface stress effects can be observed for el2 than £,e,. Similar behavior is observed for shear strain. Therefore, the Eshelby's tensor in the current case is size-dependent and non-uniform for an elliptical 54 i n h o m o g e n e i t y , w h i c h i s c o n s i s t e n t w i t h the r e su l t o b t a i n e d b y S h a r m a a n d G a n t i [32 ] f o r a n i n h o m o g e n e i t y w i t h n o n - c o n s t a n t c u r v a t u r e . S i m i l a r r e su l t s are o b s e r v e d f o r o the r v a l u e s o f alb. T h e e f f e c t o f r ° i s n o t d i s c u s s e d h e r e , b u t i t c a n be e x p e c t e d that t h i s e f f e c t i s v e r y s i m i l a r to the case o f a n e l l i p t i c a l h o l e . I n a d d i t i o n , the n u m e r i c a l r e su l t s f o r p l a n e stress case are v e r y s i m i l a r to that o f p l a n e s t r a i n . - T°=0, >TS=-1 ON/m, R= 1 Onm T°=0, /T S =-1 -T°=0 , K%=0^ - 0, A?S - ON/m, fl=50nm -0.500 -0.501 - x°=0, fC^= 1 ON/m, fl=50nm -T°=0 , A^=l ON/m, R=\ Onm 0.0 0.5 1.0 1.5 2.0 2.5 3.0 6>(rad) ( a ) a / 6 = 1 . 5 -0.76 T°=0, K*= 1 ON/m, R= 1 Onm T°=0, A^=10N/m, /?=50nni/M -T"=0, £ S =-10N/m, fl=10nm 0.0 0.5 1.0 1.5 2.0 2.5 3.0 0(rad) (b ) alb=3 -0.10 -0.1 H -0.12 ^-0.13 00-0.14 -0.15-1 -0.16 -t°=0, KS=-10N/m, /?=10nm -T°=0, KS=-10N/m, «=50nm -T°=0 ,K S =0 -x°=0, ATS=10N/m,/?=50nm ^ -T°=0, KS=\ON/m, R=\Onm 0.24 0.08 4 fl=50nm\ T ° = 0 , / : S = 1 0 N / m x°=0, KS=1 ON/m, R=\Onm 0.0 0.5 1.0 1.5 2.0 2.5 3.0 0(rad) (c )a f/6=1 .5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 0(rad) (d) alb=3 F i g . 3 .13 . V a r i a t i o n o f n o n d i m e n s i o n a l s t ra ins a l o n g the i n t e r f a c e i n a n i n h o m o g e n e i t y f o r d i f f e r e n t v a l u e s o f Ks a n d R u n d e r a u n i f o r m d i l a t a t i o n a l e i g e n s t r a i n . 55 Chapter 4 2-D FINITE ELEMENT ANALYSIS OF NANOSCALE INHOMOGENEITY PROBLEMS 4.1 Finite Element Formulation In this chapter, the classical finite element method (FEM) is extended to develop a new formulation that takes into account the surface stress effects. The finite element scheme is thereafter applied to study the elastic field of a finite plate with a single nanoscale inhomogeneity or two interacting nanoscale inhomogeneities. The advantage of the finite element method when compared to the analytical methods used in Chapters 2 and 3 is that it can be easily applied to solve more complex nanoscale mechanics problems such as arbitrarily shaped inhomogeneities, anisotropic matrix and inhomogeneity materials, and unit cells containing multiple inhomogeneities. This section presents the detailed derivation of a displacement-based, two-dimensional finite element formulation for a Gurtin-Murdoch type continuum with surface stress effects by using variational methods. The problem under consideration is schematically shown in Fig. 4.1 and a Cartesian coordinate (x, y) [z is the third axis perpendicular to the (x, _y)-plane] is used in the analysis. The matrix and inhomogeneities materials have no body forces and are linearly elastic, homogeneous and anisotropic. As the present study is concerned with two-dimensional problems, the surface stress constitutive equations are isotropic. The inhomogeneities are in the nanometer range and can be arbitrarily shaped. Perfect bonding condition is assumed at the inhomogeneity-matrix interface. The matrix-inhomogeneity system is subjected to external loads on the outer surface S of the matrix materials and/or a prescribed eigenstrain in the inhomogeneities. 56 Fig. 4.1 Nanoscale arbitrarily shaped inhomogeneities in a matrix material. The total potential energy Tl of the system under consideration consists of the bulk elastic strain energy UB of the material system, the potential energy of external loads W and the surface elastic strain energy Us. Therefore, U = UB +US +W , (4.1) where f/B and W can be calculated from the classical equations as: U»= ^ ' a ^ d V - r ^ ^ d V , (4.2) and w = -j{u}T {r}ds, (4.3) Here, the superscript B denotes quantities corresponding to the bulk (both the matrix and inhomogeneities); the superscript T denotes the transpose of a matrix or vector; VM and VX denote the volume of the matrix and inhomogeneities respectively; and {T} and {W} are the surface traction vector due to applied loads and the surface displacement vector respectively. The surface/interface elastic strain energy Us associated with the surfaces stresses can be expressed as, L / s = j f a % d e a p d r , (4.4) r+s where the superscript S denotes quantities corresponding to the surface. 57 The integral with respect to dr is taken over the matrix-inhomogeneity interfaces r and the boundary surface of the matrix 5 when the boundary surface elastic strain energy has to be accounted. In the Gurtin-Murdoch surface stress model, it can be assumed that there exists a mathematical surface layer of zero-thickness on which a surface stress tensor exists. Therefore, the integral with respect to the surface/interface in Eq. (4.4) is similar to Eq. (4.3). Introduce the element shape function matrix [N(x, y)] such that M = [iV]{z7}, (4.5) where {u} is the displacement vector at a general point within an element, and {u} is the nodal displacement vector. In two-dimensional problems, these two displacement vectors can be written as, {u}= {ux,uy}T and {u}={(ux)] (uy\ (ux)2 {uy)2---{ux)n (uy)n}T , (4.6) where the subscript ' n ' denotes number of the nodes per element. Differentiating Eq. (4.5) with respect to the coordinates, the corresponding element strain vector can be expressed as, {u} = [B]{u}, (4.7) where [B] is a matrix containing derivatives of the element shape functions . For simplicity, the plane shown in Fig. 4.1 is assumed to have unit thickness in the z-direction in the evaluation of all integrals. Integrating Eq. (4.2) with respect to the strains by using linear stress-strain relations for the matrix and inhomogeneity materials results in, UB = J i M T M M W V + J i M T [ D \ {s}dV - \{ef [D\ {a }dV , (4.8) where [Z)]M and [D\ are the 2-D elasticity matrices of the matrix and inhomogeneity materials respectively, and {e*} is the prescribed uniform eigenstrain in the inhomogeneity. For linear elastic orthotropic material, the elasticity matrix in plane strain case is given by: w = d{u}_ 3X: etc,. 58 [D] = c c 12 22 0 c , 6 6 (4.9) while in plane stress case, [D} = C, c 2 C 3 3 C n -c c C 3 3 0 12 c c ^ 13'-' 23 C 3 3 r ^ 2 2 c 2 ^23 c ^33 0 0 0 c66 (4.10) where C 0 is the 3D elasticity matrix component in the Voigt notation. Substitution of Eq. (4.7) into Eq. (4.8) yields: c/ B = \i{uY[BY[D}u[B}{u}dV+ fcmBY{D}M«}dV-jfc}WM,{*>-Substitution of Eq. (4.5) into Eq. (4.3) yields: W = - J M « } ) T {T}dS = -l{uYMT {T}dS . (4.12) To obtain the surface elastic strain energy due to surface stress components in the conventional finite element form, the surface stress tensor has to be expressed in the (x, y) coordinate. Fig. 4.2 shows the surface normal-tangential coordinate system (n, t) and the Cartesian coordinates (x, y). The angle between the normal direction n and the positive x-axis is denoted by 6. Fig .4.2 Surface coordinates (n, t) and Cartesian coordinates (x, y): 59 Note that in two-dimensional problems only the surface normal stress component in the /-direction, of , enters into the calculation of the energy and the surface stresses in the (x, y) coordinates can be expressed in terms of of as, { ^ < ^ J = {P2 q2 -pqfaf,, (4.13) where p = sin 0, q = cos 0. For an isotropic surface, the surface stress constitutive relation can be expressed as [26]: * i = M a , + 2l" S - r° )sap + (As + r° \n8ap . (4.14) and the various quantities in Eq. (4.14) are defined in Chapter 2. Using Eq. (4.14), *?,=T° + K \ + ( A S + T % , (4.15) where A^s = 2ps + As - r ° , and £zz=0 for plane strain and £ _ = - { c u e x x + C 2 i £ y y ) / C 3 3 for plane stress. Standard tensor transformation rules yield, £u = P2£xx + <?£yy ~ 2Pq£*y = {?* <f ~ P^xx £x* Yty I > (4-16) Substituting Eq. (4.15) into Eq. (4.13) and expressing all the strains in the (x, y) coordinates yields R <*% <}T=[^lsK £yy / J (4.17) where Ksp4+Vxsp2 Ksp2q2+V2sp2 -Ksp3q Ksp2q2+V,sq2 Ksq4+V2sq2 -Kspq3 -KSp3q-V?pq -K'pq'-Vfpq Ksp2q2 (4.18) Here, V* = V2 = 0 for plane strain which results in a symmetric matrix [ D \ , and F, s =-(;ts +T°)CU/C33 and V2 = -(/I s +T0)C23/C33 for plane stress which results in an unsymmetric matrix [ D \ . By using Eq. (4.17) and Eq. (4.7), Eq. (4.4) can be expressed as: 60 r+s r+s = ^{uY[Bj[Dl[B]{u}dr+ \{uY[BY{p 2q 2-pq} Tr0dr. (4.19) r+s r+s Now invoke the stationary condition of n , i.e., <5Q = 0 , with respect to the unknown nodal displacement vector. In view of Eq. (4.1), the variations of the three terms in FI can be expressed as, 5U»={fiiY l[Bj[D}M[B\iV {u}+{SuY \[Bf[D\[B}iV {u} l[Bf[D\{e}dV dW = -{SuY j[Nj{T}dV SU* = {SUY{ j[BY[D\[B}ir {u}+{6uY\ \[BJ{p 2 q 2 -pqfv'dT \r+s \r+s (4.20) (4.21) (4.22) where [ D \ =J^D\ + [D j ) and is a symmetric matrix. Finally, the equilibrium equation can be expressed as, MM-I/-} . where [K]=j[BY[DUBW+l[BY[D\lB}iV+ J [ ^ M d T , (4.23) (4.24) r+s is the total stiffness matrix obtained by assembling the element stiffness matrices , and {/}= \M{T}dV+ \[BY[D\{e}dV- j[Bf{p2 f-pq?T°dr, (4.25) s V] r+s is the total equivalent nodal load vector obtained by assembling element load vectors. From Eqs. (4.23)-(4.25), it can be seen that presence of surface stresses changes the stiffness matrices and nodal load vectors of the elements that have a side on the surface/interface on which surfaces stresses exist. Moreover, the presence of a non-zero residual surface stress r° produces an additional term for the nodal force vector. Residual 61 surface stress also has an effect on the stiffness of an element as the matrix [ D \ is a function of r ° . In the absence of surface stress effects, the above equations reduce to the classical elasticity case. In the following numerical study, the effect of r° is not presented and KS = ± 1 0 N / m , and the plane strain case is considered without loss of any generality. 4.2 Finite Element Simulation of Inhomogeneity Problems A computer code based on the finite element formuation presented in the preceeding section was developed. The commercial software package Hypermesh is used to generate a two-dimensional finite element mesh and 8-node isoparametric elements are used to discretize the domain. The mesh information is then read by a F O R T R A N based finite element program developed by the author. A l l area and line integrals are computed by using 3x3 and 3 Gauss integration schemes respectively. The linear equation system is solved by an IMSL subroutine. Stresses in an element are first calculated at the 2 x 2 Gauss integration points by using basic relations in elasticity and then extrapolated to the node points by using a bilinear extrapolation. The final nodal stress is the average of the stresses of the node calculated from different elements associated with the node. 4.2.1 Elastic Field Around a Hole in a Finite Plate A finite plate of isotropic aluminum containing a nanoscale elliptical hole is subjected to uniform remote traction in the ^-direction as shown in Fig. 4.3. The properties of aluminum are given in Chpater 2. The hole with a major axis length equal to 2a and minor axis length equal to 2b is located at the middle of the plate. The origin of the Cartesian coordinate system (x, y) is at the center of the hole and the normal-tangential coordinate system («, /) at a point on the surface of the hole is as shown in Fig. 4.3. First consider the special case of a circular hole (a=b=R). Surface stresses are present around the surface of the hole and do not exist along the outer boundary of the plate. 62 "O"o y n d > x h a ± Fig. 4.3. A nanoscale elliptical hole in a finite plate. Figure 4.4 shows the variation of nondimensional hoop stress with the hole radius for 2A7/=0.05 and 2A7/z=0.05. This case approximates the behaviour of a circular hole in an infinite plane and the corresponding analytical solution is also shown in Fig. 4.4. The agreement between the two solutions is very good and confirms the high accuracy of the present finite element solutions. An interesting observation from Fig. 4.4 is the unstable beahviour of the solution below some hole radius when Ks <0. Similar behaviour was noted earlier in the case of the analytical solution of an elliptical hole. It is noted that stiffness matrix of the finite element mesh becomes non-positive definite below a certain R value when Ks < 0. The reason is that, forKs < 0, the surface elastic energy is negative and it could dominate over the bulk elastic energy for small values of R. As found in Chapter 3 using the analytical solution, hoop stress can also be singular in the case of an infinite plane with a circular hole but for R less than lnm. However, in the case of finite element analysis, the plane has finite dimensions and this causes the unstable behaviour of the solution to occur at a R value greater than lnm. Such unstable behaviour does not happen for A^s > 0 in which case hoop stress becomes smaller as R decreases. For A^s > 0, the total stiffness matrix remain positive definite for all values of R. 63 *—,--•>-: T^O.A^-lON/m — , :x°=0,^=0 -:T°=0, J^=10N/m i ' 1 ' 1 1 1 ' 1 1 0 10 20 30 40 50 R(nm) Fig. 4.4 Comparison of analytical and finite element hoop stress solutions (dash lines for analytical solution of infinite plate and solid lines for finite plate with 2J?//=0.05 and 2R/h=0.05). Consider next the case of a circular hole in a finite plate. The distance between the center of the hole and the top and bottom edges is kept at ten times the radius of the hole, i.e., 2R/h=0.\, and hoop stress at 0 = 0 is shown in Fig. 4.5 for two values of the length ratio (2R/l=0.\ and 0.4). The results are similar to the case of an infinite plane with a circular hole [Fig. 2.2(b) in Chapter 2] and Fig. 4.4, but the stresses are due to the plate boundary effects. As in the case of Fig. 4.4, hoop stress becomes highly oscillatory and singular at some radii when Ks < 0. The radius below which this unstable behaviour occurs increases as the plate becomes smaller. 64 4.0 3.2 J T°=0, A:S=-10N/m t°=0, A^=0 t°=0, A^lON/m x°=0,^S=-10N/m x°=0,KS=0 —*— T°=0, /TS=10N/m 2.4 2.5 0 10 20 30 40 50 0 10 20 30 40 R(nm) (b) 2A7/=0.4and 2R/h=0.\ 50 60 fl(nm) (a)2A7/=0.1 and 2R/h=0A Fig. 4.5. Variation of nondimensional hoop stress at 0 = 0 of a finite plate with R. Fig. 4.6 shows the variation of nondimensional hoop stress on hole surface with angle 0. The solution behaviour is very similar to the infinite plane solution in Chapter 2. The maximum value of hoop stress increases as the plate dimensions get smaller. Note that although the R values corresponding to K =10 and K = -10 are different, hoop stress solutions for the two cases are quite close to each other. Fig. 4.7 shows the nondimensional hoop stress at 0 = 0 for various values of 2R/l. In this figure, A=30nm for A^s <0 as the stress is unstable for R less than 30nm when 2R/l = 0.5 . As 2R/t increases, hoop stress corresponding to the classical and current cases increase nonlinearly. 65 2 ^ 1 T ^ O . ^ - l O N / m l i°=0, A^=0 T°=0, A ^ l O N / m -1 0.0 0.3 0.6 0.9 1.2 1.5 0(rad) (a) 2J?//=0.1 and 2R/h=0A (R=5nm for Ks >0 and#=15nmfor Ks <0) 0.0 0.3 0.6 0.9 1.2 1.5 0(rad) (b) 2/J//=0.4 and 2R/h=0.\ (#=5nm for A: s > 0 and i?=25nm for Ks < 0) Fig. 4.6. Variation of nondimensional hoop stress at hole surface with G. Fig. 4.7 Variation of nondimensional hoop stress at 6 = 0 with 2R/1 (2R/h=0A; R=5nm for / : s > 0 and i?=30nm forKs < 0). Consider now the case of an elliptical hole with a/b=\.5 (Fig. 4.3). Fig. 4.8 shows the variation of the nondimensional tangential stress on the hole surface with R [=(a+b)/2] and 0 for 26//? = 0.1 and 2a//= 0.1 and 0.4. The results are similar to the case of an infinite plane containing an elliptic hole with alb=\.5 [Figs. 3.3 and 3.5(b) in Chapter 3]. 66 As in the case of analytical solution, hoop stress becomes unstable for R less than 14nm and 24nm for 2a/1 = 0.1 and 0.4 respectively. Note that stress instability starts when R is less than 6nm for an infinite plane with an identical elliptical hole (Chapter 3). Fig. 4.9 shows the nondimensional hoop stress at 9 = 0 for various values of 2a/l and the behavior is similar to the case of a circular hole. 4.4 4.0 = 3.6 32 2.8 • f f ' k —•—,- - ° - : T°=0, A^=-10N/ra , : x°=0, A^=0 A -: x°=0, A^=10N/m 10 20 30 R(nm) 40 50 o D (a) Variation with R: Dash line for infinite plate and solid line for finite plate with 2b/h=0A and 2a//=0.1 - ~ — T°=0, /T^ - lON/m T° = O , A : S = O T°=0, K S=10N/m 20 30 40 50 60 fl(nm) (b) Variation with R: 2b/h=0.\ and 2a/l=0A 0.6 0.9 0(rad) (c) Variation with 9 (2b/h=0A and 2a//=0.1; R=5nmforKs >0 and A=15nmforA^s <0) 0.6 0.9 0(rad) (d) Variation with 9 (2b/h=0.1 and 2a//=0.4; tf=5nm for Ks > 0 andA=25nmforA: s <0) Fig. 4.8 Variation of nondimensional tangential stress of an elliptical hole with R and 9 (a/b=\.5). 67 0.1 0.2 0.3 0.4 0.5 lall Fig. 4.9 Variation of nondimensional hoop stress at 0 = 0 with lall (2b/h=Q.\;R=5mn for Ks > 0 and £=30nm foriC s < 0). A remotely loaded semi-infinite plate of isotropic aluminum with a circular hole as shown in Fig. 4.10 is considered. The objective is to consider the interaction between the hole and free surface. In the finite element analysis, the distances between the center of the hole and the left, right and bottom edges are set to ten times the radius of the hole, i.e., h/R=l0 and l/2R=\0 to simulate the case of a semi-infinite medium. The influence of free surface is controlled by the distance between the center of the hole and top surface which is denoted by D. The origins of the Cartesian coordinate system (x, y) and polar coordinate system (r, 9) are located at the center of the hole. 68 y D r x R h 0 0 Fig. 4.10 A nanoscale circular hole in a semi-infinite plate. Fig. 4.11 shows the variation of nondimensional hoop stress at 0 = 7r/2 6 = -nj'2 with hole radius R for D/R=2 and 1.1. The stress field is unstable for R less than 17nm and 22nm for D/R=2 and 1.1 respectively when Ks is negative. As expected hoop stress at 0 = -7r/2 is smaller than that at 0 - n/2 . A substantial increase in hoop stress is noted at 0 = nj2 when D//?=l.l. Fig. 4.12 shows the variation of nondimensional hoop stress along the hole surface. Hoop stress is not symmetric with respect to 0 = 0 in the present case and the dependence on 6 becomes quite complex when D/i?=l.l. The surface stress effect is much more pronounced on the portion of the hole surface that is closer to the free surface when D//?=l.l. 69 (a)D/J?=2(0 = - ; r /2) (b) D/R=2 (9 = TT/2) 1="- 1 • 1 • 1 • 1 • 1 ' 1 ' 1 • ' P i , • 1 • 1 • 1 • 1 r—1 • 1 • 1 0 10 20 30 40 50 60 70 0 10 20 30 40 50 60 70 80 (^nm) tf(nm) (c)D/R= 1.1 (0 = -x/2) (d)D/7?=l.l (0 = n/2) Fig. 4.11 Variation of nondimensional hoop stress at 9 = ± nj2 with hole radius for two different D/R values (D/R=2 and 1.1). 70 ^(rad) 6>(rad) (a) D/R=2; R=5nm for Ks > 0 and (b) D/R=\. 1; tf=5nm for Ks > 0 and i?=17nmforA:s <0 #=22nm for ATS < 0 Fig. 4.12 Variation of nondimensional hoop stress on hole surface with angle 6 (D/R=2 and 1.1). Fig. 4.13 shows the variation of nondimensional hoop stress at 6 = 7v/2 and 0 = -7r/2 with D/R. As D/R decreases, hoop stress increases rapidly at both points but the influence at 6 = K/2 is higher. The surface stress effect is more pronounced for positive Ks values at 0 = /r/2 . When ^ s = 1 0 N / m and i?=5nm, the increases of nondimensional hoop stress due to surface stress effect are 6.7% for D/R=2 and 8.4% for D/i?=l.l at 0 = -7c/2, and the changes are 8.7% for D/R=2 and 31.2% for D/R=\.l at 0 = TT/2. 71 (a)0 = - ;r/2 (b)d = x/2 Fig. 4.13 Variation of nondimensional hoop stress with D/R (i?=5nm for Ks > 0 and R=22nm for A' 5 < 0). 4.2.2 Elastic Field of a Finite Plate with Two Circular Inhomogeneities This section considers a finite plate of GaAs, containing two circular InAs inhomogeneities with identical radii as shown in Fig. 4.14. Both the matrix and inhomogeneities materials are assumed to be isotropic and their properties are given in Chapter 2. The plate is subjected to uniform remote tension in the vertical direction. The two inhomogeneities are at a distance 2D between their centers, and are located symmetrically in the plate. Let l/R=\0 and h/R=20, so that the effect of the plate edges on the elastic field in and around the inhomogeneities is very small when compared to the interaction between the two inhomogeneities if the distance between them is small. Only one quarter of the problem has to be considered in the finite element analysis and the Cartesian coordinate system (x, y) and polar coordinate system (r, 6) are defined as shown in Fig. 4.14. 72 1^0 y l c r 0 Fig. 4.14 Two circular inhomogeneities in a finite plate. Fig. 4.15 and 4.16 shows the variation of nondimensional hoop stress in the matrix and inhomogeneity at two locations (9 = 0 and 9 = n ) on their interface for different values of radius R and angle 6 when D/R=\.\ and 2.0 respectively. As D/R decreases, the interaction is more significant. This can be seen by comparing Fig. 4.15 and Fig. 4.16. The general variation of hoop stress is similar to the case of a single inclusion presented in Chapter 2. Hoop stresses in both the matrix and inhomogeneity are smaller at 9 = 0 when compared to 9 = n and so is the surface stress effect. The stress is unstable for R less than 7nm for Ks = - l O N / m . 73 o ^ 1 . 0 D 0.8 . - : T ° = 0 , A:S=-10N/m - : T°=0, Ks=0 ^-:T 0 =0, /T^lON/m 20 R(nm) (a) 0 = 0 30 1.44 1.2 o ^ 1 . 0 to 0.84 40 T0=0,.Ks=-10N/m : T°=0, ^ = 0 : x°=0, ^=10N/m 10 20 #(nm) (b) 0 = 7t 30 40 1.0 1.5 2.0 £>(rad) (c) Variation along the interface (7?=5nm for Ks > 0 and R=6nm for < 0) Fig. 4.15 Variation of nondimensional interfacial hoop stress with radius and angle (D/R=2; Solid lines for matrix and dash lines for inhomogeneity). 74 l , 1 , , . p = — , , , 1 . 0.0 0.5 1.0 1.5 2.0 2.5 3.0 6V(rad) (c) Variation along the interface (A=5nm for Ks > 0 and R=7nm fovKs < 0) Fig. 4.16 Variation of nondimensional interfacial hoop stress with radius and angle (D/R=l.\; Solid lines for matrix and dash lines for inhomogeneity). Fig. 4.17 shows the variation of nondimensional hoop stress in the matrix and inhomogeneity at two locations on their interface (8 = 0 and 6 = n) with D/R. Note that R=5nm for Ks = 10N / m and R-7 nm for KS = -10N / m . Hoop stresses at 6 = 0 do not change much with the D/R and the difference between the solutions with or without surface stress effects is similar to the case of a single inhomogeneity. However, more influence of the interaction between the inhomogeneities is noted at 0 - n but this effect 75 diminishes rapidly with increasing D/R and the interaction between the inhomogeneities vanishes when D/R > 2 but substantial influence can be seen at 6 = n when DIR less than 1.5. The surface stress effect is a little more pronounced at 0 = n than at 0 = 0. D/R D / f t (a) Matrix (b) Inhomogeneity Fig. 4.17 Variation of nondimensional hoop stress with D/R (A=5nm for Ks > 0 and R=7nm for Ks < 0; Solid lines for 6 = n and dash lines for 6 = 0). 4.2.3 Elastic Field Around a Circular Hole in a Finite Anisotropic Plate Consider a circular hole (a=b=R) embedded in the middle of a finite plate under uniform remote loading in the ^-direction (Fig. 4.3). The plate under consideration is FCC A l which has cubic symmetry. The coordinate directions are assumed to be the crystallographic directions and the three independent elastic constants are: C„ =118.10GPa, C 1 2 = 62.293GPa and C 6 6 = 36.706GPa [24]. The surface properties are assumed to be isotropic. Therefore, the new F E M formulation developed in section 4.1 is also applicable in the case of plane strain state. Figs. 4.18 and 4.19 shows the variation of nondimensional hoop stress with the hole radius and angle respectively for 2A7/=0.05 and 2A7//=0.05. The classical nondimensional hoop stress at 9 = 0 is a little far from 3 due to the anisotropic properties of Aluminum. The surface stress effects are similar to that in the isotropic case (Figs. 4.4 and 4.6). The hoop stress Is also found to become unstable for the radius less than lOnm when the surface elastic modulus,./^, is 76 negative. Fig. 4.20 shows the variation of nondimensional normal displacement along the hole surface for both isotropic and anisotropic cases when 2R/l=0.Q5 and 2R/h=0.05. In both cases, the nondimensional displacements are increased or decreased when Ks < 0 or Ks>0. T°=0, A^=-10N/m T°=0, A^=0 T°=0, / T ^ l O N / m 0 40 50 10 20 30 R(nm) Fig. 4.18 Variation of nondimensional hoop stress on hole surface at 6 = 0 with R (2/?//=0.05 and 2R/h=0.05). 2 0 •1 -x°=0, K 5 ^ - 1 0 N / m •T°=0, K?=0 - x ° = 0 , A : S = 1 0 N / m 0.0 0.3 0.6 0.9 1.2 1.5 0(rad) Fig. 4.19 Variation of nondimensional hoop stress along hole surface with 0 (2M=0.05 and 2MJ=0.05). 77 0.0 0.3 0.6 0.9 1.2 1.5 0(rad) Fig. 4.20 Variation of nondimensional normal displacement with 6 (2R/l=0.05 and 2R/h=0.05; Solid lines for anisotropic case and dash lines for isotropic case). 78 Chapter 5 S U M M A R Y A N D C O N C L U S I O N S 5.1 Summary and Major Findings The major findings and conclusions of the current study are given below. ( 1 ) Based on the Gurtin-Murdoch surface/interface elasticity model, a closed-form analytical solution is obtained for the elastic field of an infinite isotropic elastic matrix containing a nanoscale circular inhomogeneity under arbitrary remote loading or a uniform eigenstrain in the inhomogeneity. The extension of the complex potential function method of Muskhelishvili is demonstrated for the first time for nanoscale inhomogeneity problems. The new solution reduces to the classical elasticity solution in the absence of surface stress effects. The stress state shows strong dependency on inhomogeneity radius, surface elastic constants and residual surface stress when the inhomogeneity radius is less than lOnm. Hoop stress around matrix-inhomogeneity interface can be increased or decreased due to surface stress effects and can be singular for an inhomogeneity with radius less than lnm when the surface stress modulus, Ks, is negative. The Eshelby tensor for a uniform eigenstrain is uniform but becomes size-dependent for a nanoscale circular inhomogeneity. The effect of surface stresses becomes negligible at a distance greater than three to four times the radius of the inhomogeneity. The circular nanoscale inhomogeneity solution presented in this thesis is a benchmark solution in nanomechanics and can be used in the validation of numerical methods such as the finite element method and in the study of effective properties of nanocomposites with dilute concentrations of reinforcing particles. (2) Following the successful derivation of the closed-form analytical solution for a circular inhomogeneity, the complex potential function method is extended to consider the case of a nanoscale elliptical inhomogeneity. The elliptical inhomogeneity problem is fundamentally important to the manufacture of nanocomposites and quantum dots as it allows some basic understanding of the optimum shape of inhomogeneities for enhanced properties and load transfer. The analytic potential functions, which are expressed by 79 infinite power series, can only be obtained approximately for an elliptical inhomogeneity. The elastic field is therefore investigated through a comprehensive numerical study. As in the case of a circular inhomogeneity, the elastic state is size-dependent, and hoop stress around a hole or an inhomogeneity interface is decreased or increased depending on whether Ks is positive or negative. An interesting feature of the elliptical inhomogeneity solution is the instability of the elastic field below a certain value of the sum of the half-lengths of the major and minor axes of the inhomogeneity for negative A^ s values. Such unstable phenomenon does not occur when A^s is positive. The reason for this behaviour is the non-positive definite nature of the elastic energy associated with the surface stresses when A^ sis negative and its dominance over the bulk strain energy for very small dimensions of the inhomogeneity. In such situations, the Gurtin-Murdoch model essentially breaks down and atomistic or other models need to be chosen. The Eshelby tensor for a uniform eigenstrain is size-dependent and no longer uniform in the case of an elliptical inhomogeneity. (3) The classical finite element method (FEM) is extended in this thesis to develop a new formulation for the analysis of two-dimensional problems involving an elastic matrix with arbitrary shaped nanoscale inhomogeneities. The formulation is based on the Gurtin-Murdoch continuum model and the principle of minimum total potential energy. The influence of the surface stresses appears in the formulation through a term corresponding to surface elastic strain energy due to surface stresses. It is found that the stiffness matrix associated with the surface stresses can be positive definite or non-positive definite depending on the values of the surface elastic constants. This can lead to a total stiffness matrix that is non-positive definite and the resulting elastic field shows unstable and singular behaviour. Unstable behaviour of the elastic field is noted in the finite element analysis for all negative values of Ks when the characteristic dimension of the inhomogeneity is below a certain value. As the dimensions of the matrix domain are increased, the finite element solutions become very close to the analytical solution for an infinite matrix. This confirms the high accuracy of the present finite element scheme. 80 (4) The stress field of a finite plate with a nanoscale circular or elliptical hole under remote tension is studied by using the F E M . As the ratio of the width of the plate to the hole radius decreases, the surface stress effects become more pronounced, and the hole size below which the elastic field becomes unstable for negative Ks also increases. Surface stress effects are more significant in the case of an elliptical hole in a finite plate and further increase as the ellipse becomes flatter. The stress field of a semi-infinite plate with a circular hole near the free edge is obtained by using the F E M . When the distance between the centre of the hole and the free edge is greater than three times the radius of the hole, the edge effect is quite small. As the hole approaches the edge, the edge effect and the surface stress effects result in a significant increase of hoop stress. The stress field of a finite plate with two identical circular inhomogeneities is also studied by using the F E M . The general trend of the solution is similar to that of a single circular inhomogeneity. As the distance between the two inhomogeneities decreases, their interaction becomes more significant and the surface stress effect becomes more pronounced. Unstable behaviour of the elastic field is also noted in this case when the inhomogeneity radius is below a certain value and K5 is negative. The elastic field of a circular hole in an anisotropic matrix is shown and the result is found to be similar to that of the isotropic case. 5.2 Suggestions for Future Work Current understnading of mechanics of nanomaterials is very limited. The study of nanoscale inhomogeneity problems has important implications in the development of nanomaterials with superior properties. The presnt study has examined some basic problems related to nanoscale inhomogeneities and further studies are necessary to advance the development of nanomaterials and the scientific base for nanotechnolgy. It is therefore recommeded that following studies be undertaken to advance the current state-of-the-art of mechanics of nanomaterials. (1) Isotropic behavior of matrix, inhomogeneity and surface/interface is assumed in the analytical solutions derived in this thesis. However, anisotropic behaviour is more practically useful in nanotechnology applications. Benchmark analytical solutions for a 81 circular (2-D) and spherical (3-D) inhomogeneity in an anisotropic matrix will be useful to gain some fundamental understanding of the behaviour of nanocomposites and other nanomaterials. The extension and numerical implementation of the new F E M formulation developed in the current study to include a range of anisotropic material models for both bulk and surface responses and multiple inhomogeneities of different shapes will be useful to the study of nanomaterials and calculation of effective properties of nanocomposites. (2) Inhomogeneities are assumed to be perfectly bonded to the matrix in the present study. It is useful to examine the cases of slip, twist and wrinkling of the inhomogeneity surface/interface to gain some fundamental understanding of the failure modes of nanocomposites and other nanomaterials. In addition, the matrix and inhomogeneity materials are assumed to be perfectly elastic in the present study but plastic deformations can play a significant role at the nanoscale. Consideration of plasticity effects may be useful in examining the instability of elastic field that is observed for negative values of Ks when the inhomogeneity characteristic dimensions are below a certain value. Such unusual behaviour may not be present in the case of an elasto-plastic model. (3) This thesis establishes the necessary framework to conduct a finite element study of a unit cell of a nanocomposite material. A comprehensive study of a unit cell taking into account different inhomogeneity shapes and sizes, anisotropic behaiviour of the matrix/inhomogeneities, interface effects, inhomogeneity densities and distributions, etc. will be very useful to the study of effective properties of nanocomposites. (4) Experimental studies and atomistic simulations are recommended to validate the applicability of the Gurtin-Murdoch model and the analytical and F E M solutions obtained in the thesis. 82 B I B L I O G R A P H Y [I] Feynman, R.P., "There is plenty of room at the bottom", In: Gilbert, H.D. (Eds), Miniaturization, Reinhold, New York, 1961. [2] Smalley, R. E., "Congressional hearings", Summer 1999. [3] "Nanotechnology - a revolution in the making -- vision for R & D in the next decade", A report of the interagency working group on the nanoscience, engineering, and technology, presented to the OSTP committee on technology, March 10, 1999. [4] Klabunde, K.J . , "Introduction to nanotechnology", In: Klabunde, K.J . , (Eds), Nanoscale Materials in Chemistry, John Wiley & Sons, New York, 2001. [5] Dai, H. , Hafner, J.H., Rinzler, A .G. , Colbert, D.T. and Smalley, R.E., "Nanotubes as nanoprobes in scanning probe microscopy", Nature, Vol . 384 (14), pp. 147-150, 1996. [6] Kassner, M.E . et al, "New directions in mechanics", Mech. Mater., Vol . 37(2-3), pp. 231-259,2005. [7] Kouris, D. and Gao, H. , "Nanomechanics of surfaces and interfaces", ASME J. Appl. Mech., Vol 69(4), pp. 405-406, 2002. [8] Maranganti, R. and Sharma, P., " A review of strain field calculations in embedded quantum dots and wires", In: Reith, M . and Schommers, W., (Eds.), Handbook of Theoretical and Computational Nanotechnology, American Scientific Publishers, Stevenson Ranch, California (Chapter 118), 2006. [9] Ji, B. and Gao, H. , " A study of fracture mechanisms in biological nano-composites via the virtual internal bond model", Mater. Sci. Eng. A, Vol . 366(1), 96-103, 2004. [10] Streitz, F.H., Cammarata, R.C. and Sieradzki, K. , "Surface-stress effects on elastic properties. I. Thin metal films", Phys. Rev. B, Vol . 49(15), pp. 10699-10706, 1994. [II] Dingreville, R., Qu, J. and Cherkaoui, M . , "Surface free energy and its effect on the elastic behavior of nano-sized particles, wires and films", J. Mech. Phys. Solids, Vol . 53(8), pp. 1827-1954, 2005. [12] Yakobson, B.I., "Nanomechanics", In: William, A . et al. (Eds), Handbook of Nanoscience, Engineering and Technology (Chapter 17), CRC Press, Boca Raton, 2003. 83 [13] Wong, E.W., Sheehan, P.E. and Lieber, C M . , "Nanobeam mechanics: elasticity, strength, and toughness of nanorods and nanotubes", Science, Vol . 277(5334), pp. 1971-1975, 1997. [14] Cuenot, S., Fretigny, C , Demoustier-Champagne S. and Nysten, B., "Surface tension effect on the mechanical properties of nanomaterials measured by atomic force microscopy", Phys. Rev. B, Vol . 69(16), pp. 165410-1-165410-5, 2004. [15] Singh, R.P., Zhang, M . and Chan, D., "Toughening of a brittle thermosetting polymer: effects of reinforcement particle size and volume fraction", J. Mater. Sci., Vol . 37(4), pp. 781-788, 2002. [16] Cadek, M . et al., "Reinforcement of polymers with carbon nanotubes: the role of nanotube surface area", Nano Lett., Vol . 4(2), pp. 353-356, 2004. [17] Sun, C T . and Zhang, H. , "Size-dependent elastic moduli of platelike nanomaterials", J. Appl. Phys., Vol . 93(2), pp. 1212-1218, 2002. [18] Liang, H. , Upmanyu, M . and Huang, H. , "Size-dependent elasticity of nanowires: nonlinear effects", Phys. Rev. B, Vol . 71(24), pp. 241403-1-241403-4, 2005. [19] Rubin, M.B. and Benveniste, Y . , " A cosserat shell model for interphases in elastic media", J. Mech. Phys. Solids, Vol . 52(5), pp. 1023-1052, 2004. [20] Gibbs, J.W., The Scientific Papers of J. Willard Gibbs. Vol . 1. Longmans Green, London, 1906. [21] Cammarata, R . C , "Surface and interface stress effects in thin films", Prog. Surf. 5c/., Vol . 46(1), 1-38, 1994. [22] Shuttleworth, R., "The surface tension of solids", Proc. Phys. Soc. A, Vol . 63, pp. 444.457, 1950. [23] Nix, W.D. and Gao, H. , "An atomistic interpretation of interface stress", Scr. Mater., Vol . 39(12), pp. 1653-1661, 1998. [24] Miller, R.E. and Shenoy, V .B . , "Size-dependent elastic properties of nanosized structural elements". Nanotechnology, Vol . 11(3), pp. 139-147, 2000. [25] Shenoy, V .B . , "Atomistic calculations of elastic properties of metallic fee crystal surfaces", Phys. Rev. B, Vol . 71(9), pp. 094104-1-094104-11, 2005. [26] Gurtin, M.E . and Murdoch, A.I., " A continuum theory of elastic material surfaces", Arch. Rat. Mech. Anal., Vol . 57(4), pp. 291-323, 1975. 84 [27] Gurtin, M.E. , Weissmuller, J. and Larche, F., " A general theory of curved deformable interfaces in solids at equilibrium", Philos. Mag. A, Vol . 78(5), pp. 1093-1109,1998. [28] Gurtin, M.E . and Murdoch, A.I., "Surface stress in solids", Int. J. Solids Struct., Vol. 14(6), pp. 431-440, 1978. [29] Weissmuller, J. and Cahn, J.W., "Mean Stresses in microstructures due to interface stresses: A generalization of a capillary equation for solids", Acta Mater., Vol. 45(5), pp. 1899-1906, 1997. [30] Ibach, H. , "The role of surface stress in reconstruction, epitaxial growth and stabilization of mesoscopic structures", Surf. Sci. Rep., Vol . 29(5-6), pp. 195-263, 1997. [31] Sharma, P., Ganti, S. and Bhate, N . , "Effect of surfaces on the size-dependent elastic state of nano-inhomogeneities", Appl. Phys. Lett., Vol . 82(4), pp. 535-537, 2003. [32] Sharma, P. and Ganti, S., "Size-dependent Eshelby's tensor for embedded nano-inclusions incorporating surface/interface energies", ASME J. Appl. Mech., Vol . 71(5), pp. 663-671,2004. [33] Duan, H.L., Wang, J., Huang, Z.P. and Karihaloo, B.L. , "Size-dependent effective elastic constants of solids containing nano-inhomogeneities with interface stress", J. Mech. Phys. Solids, Vol . 53(7), pp. 1574-1596, 2005. [34] Eshelby, J.D., "The determination of the elastic field of an ellipsoidal inclusion and related problems", Proc. R. Soc. London, Ser. A, Vol . 241, 376-396, 1957. [35] Eshelby, J.D., "The elastic field outside an ellipsoidal inclusion", Proc. R. Soc. London, Ser. A, Vol . 252, pp. 561-569,1959. [36] Eshelby, J.D., "Elastic inclusions and inhomogeneities", In: Sneddon, I.N. and Hi l l , R. (Eds), Progress in Solid Mechanics 2, North Holland, Amsterdam, pp. 89-140,1961. [37] Mura, T., Micromechanics of Defects in Solids, Martinus Nijhoff, Dordrecht, Netherlands, 1987. [38] Nemat-Nasser, S. and Hori, M . , Micromechanics: Overall Properties of Heterogeneous Solids, Elsevier, New York, 1999. [39] Markov, K. and Preziosi, L., Heterogeneous Media: Micromechanics Modeling Methods and Simulations, Birkhauser Verlag, Switzerland, 2000. 85 [40] Yang, F., "Size-dependent effective modulus of elastic composite materials: spherical nanocavities at dilute concentrations", J. Appl. Phys., Vol . 95 (7), pp. 3516-3520,2004. [41] Duan, H.L., Wang, J., Huang, Z.P. and Karihaloo, B.L. , "Eshelby formalism for nano-inhomogeneities", Proc.R. Soc. A, Vol . 461(2062), pp. 3335-3353, 2005. [42] Gao, W., Yu, S.W. and Huang, G.Y., "Finite element characterization of the size-dependent mechanical behavior in nanosystems", Nanotechnology, Vol . 17(4), pp. 1118-1122,2006. [43] Muskhelishvili, N.I., Some Basic Problems of the Mathematical Theory of Elasticity. P. Noodhoff Ltd., Groningen, The Netherlands, 1963. [44] Ruud, J.A., Witvrouw, A. and Spaepen, F., "Bulk and interface stresses in silver-nickel multilayered thin films", J. Appl. Phys., Vol . 74(4), pp. 2517-2523, 1993. [45] Josell, D., Bonevitch, J.E., Shao, I. and Cammarata, R.C., "Measuring the interface stress: Silver/nickel interfaces", J. Mater. Res., Vol . 14(11), pp. 4358-4365, 1999. [46] Meyers, M . A . and Chawla, K . K . , Mechanical Behavior of Materials. Prentice Hall, Englewood Cliffs, NJ , pp. 92, 1999. [47] Zhou, L .G. and Huang, H. , "Are surfaces elastically softer or stiffer?", Appl. Phys. Lett., Vol . 84(11), pp. 1940-1942, 2004. [48] L i . Z.R., Lim, C.W. and He. L .H. , "Stress concentration around a nano-scale spherical cavity in elastic media: effect of surface stress", Eur. J. Mech. A/Solids, Vol . 25(2), pp. 260-270, 2006. [49] Sharma, P. and Ganti, S., "Interfacial elasticity corrections to size-dependent strain-state of embedded quantum dots", Phys. Stat. Sol. B, Vol . 234(3), pp. R10-R12, 2002. [50] England, A .H . , Complex Variable Methods in Elasticity. Wiley-Interscience, London, 1971. [51] Stagni, L, "Elastic field perturbation by an elliptic inhomogeneity with a sliding interface", J. Appl. Math. Phys. (ZAMP), Vol . 42(6), pp. 811-820, 1991. [52] Shen, H. , Schiavone, P., Ru, C Q . and Mioduchowski, A . , "Interfacial thermal stress analysis of an elliptic inclusion with a compliant interphase layer in plane elasticity", Int. J. Solids Struct., Vol . 38(42-43), pp. 7587-7606, 2001. 86 Shen, H. , Schiavone, P., Ru, C Q . and Mioduchowski, A. , "An elliptic inclusion with imperfect interface in anti-plane shear", Int. J. Solids Struct., Vol . 37(33), pp. 4557-4575, 2000. 87 APPENDICES Appendix A T T • • 1 (\-m2)IR , . ,. , Using the expression for the curvature — = — (see Appendix B) and Eq. (3.9), the interface stress in the matrix can be expressed as: .2 2 u It R2 + A, 0 M ( £ ) < f + 0 M ( £ ) c f R + ( A , + A 4 ) A ; R R R , o n T , (A.1) where A , = , KS =2p.s+As-r°, A2 A (A.2) By using Eqs. (3.9), (3.35) and (3.36), the derivative term in Eq. (3.38) is written as: ii^km\^m\4)/R dt 1 _ M ^ [ _ ^ + (1 + w 2 ) 5 3 _ w | 5 ] _ M M ( | ) [ 2 ( 1 + W 2 ) ^ _ 4 W ] \m{g)lR\ [ R R A , e M ( ^ ) r _ R [_ m ^ + (i + m i ) ^ _ m ^ ] + AIQM(£) [2(1 + m 2 ) £ - * _ Am] R m ( | ) I i i?2 + A 2 ( A , + A 4 ) 7? tf2 on r. (A.3) In the inclusion, considering the effect of eigenstrain, the actual strain is regarded as the sum of elastic strain and eigenstrain. However, the calculation is the same as that in the matrix and the result is similar. In Eqs. (A.1)-(A.3), (j)u and 0 M are replaced by <f)x 88 and 0 , , respectively; A, and A 2 by A 6 and A 7 , respectively; A 4 by A 5 ; r° by r° + 8(A6 + A^^Re*, and others keep the same . Here, /L s +r° .v/, 4//,A v, -1 A 6 = A 7 (A.4) m'(£)m (£) R Multiplying Eq. (3.38) by a factor ^ ^ and using Eqs. (3.17), (3.18), m (3.24), (3.25), (A.l) and (A.3) yields: R 0 dt -4i6 , 4 i 0 + __ <?„ei("+3)(? + __ r„e- (" + 3 ) e + __ /7„e i (" + 4 ) < ? +]T K„e" •i(n+4)9 «=1 n=1 , on r, (A.5) where g 0 = j r ° ( l - m 4 ) + A 2 ( A , + A 4 ) ( 2 - 2 m V + 4A 1m(£> 0-D 0) + A 2 ( A , + AA)[(5m-m3)A] + ( -3w-m 3 )J 1 ]+12A 2 (A, + AA)m2(A3 + A\) + A, (1 - m 2)(D 2 + D2) + 4A,m(-D 4 + 7J4), g, =14A 2(A, +A4)wv42 +6A 2 (A, + AA)m2A2 - 2 0 A 2 ( A , + A4)m24, + 2A,Z>3 + 3A]mD3 - 5A,m£>5, g 2 = - 6 A 2 ( A , + A 4 ) w 2 ^ 2 +A 2 (A , +A 4 ) ( -10w-4m 3 ) I 2 +20A 2 (A, +A4)m24, - 3A,m£>3 - 2A,w 2 A + 5A,/wZ)s, g 3 = - i r 0 w ( l - ^ 2 ) + A 2 ( A , + A 4 ) ( w - m 3 ) ^ + A,(3 + w2)7J0 + A 2 ( A , + A 4)[(-3 + m2)A] + 2m2A] + (27m + 3m3) A, -30m 2 ^ 5 ] + 2A,w(D 2 + A ) + A,(3 + m 2 )D 4 - 6 A , w D 6 , gA = -±T°m(\-m2) + A2(A]+A,)(3m3-3m)A-Al(l + 3m2)D0 + A 2 (A, + A 4 ) [ - 2 m 2 ^ + (1 + m2 )A\ + (-21m - 9m3 ) I 3 + 30m2% \ -2A]mD2-2A]mD2+Al(-\-3m2)D4+6AlmD6, 89 g5 - A 2 ( A , + A 4)[(-8 - 2m2)A2 + (44m + 8ro 3 )4 - 42m2A6] + A,w£>3 +A,(4 + 2w 2 )D 5 -7A]mD1, g6 = A 2 ( A , + A4)[(4 + 6m 2)Z 2 +(-36w-16w 3 )4, + 42m2A6\ - A,mZ>3 + A,( -2 - 4w2)7J5 + 7A,mZ)7, S„ = A 2 (A, + A 4)[/«(« - n2 )A„ + (* + 2)[n + 2 + (2» + 3)m2 ]An+2 ] + A 2 ( A , + A4){(« + 4)[w + m3 -(n + 5)(2m + m3)]An+4 + (n + 6)(n + 7)m2An+6} + A,(n -l)m7Jn+3 + A,[3 + m2 - (n + 5)(1 + m2)]Dn+5 + A, (n + 7)mDn+1, T= A 2 (A, + A 4 )[(«2 - n)mAn - (n + 2)(2nm2 +m2+n + 4)An+2 ] + A 2 ( A , + A 4)[(« + 4)(«m 3 + 2m3 + 2nm + \\m)An+4 -(n + 6)(n + 7)m2An+6\ + A, (1 - n)mDn+3 + A, [(n + 5)(1 + m2) - (1 + 3w2)]Z)„+ 5 - A, (/, + 7)mD„ + 7 , (A.6) /*0 = i [ r ° + 8(A 5 + A 6 ) M l R £ * | l - m 4 ) + A6[4m(O0 -O0) + (1 -m2)(P2 + P2) + 4m(P< -/>)] +A 7 (A 5 + A 6)[(l -5m2)E, + (1 + 3m2)Ex + \2m(Ei - £ 3 ) j + A 7 ( A 5 + A 6 ) [ ( 5 m - m 3 ) ^ + ( -3m-m 3 ) ^ + 12m2(F3 - F 3 ) ] , /T, = -5A 6 mO| + 3A 6 m^ + 2 A 6 ^ + 2A 6 P 3 + 3A6m?] - 5A6mP5 +A1(A5+A6)[6mE2+(4 + \0m2)E2-20mE4] +A 7 (A 5 +A 6)(l4mF 2 + 6m2F2-20m2F4), h2 = 5A 6mO, - 2A6m2P{ - 3A6m?] - 3A6mP3 - 2A6m2P] + 5A6mP5 - A 7 ( A 5 + A j ( l 4 m 2 £ 2 +6m£ 2 - 2 0 m £ 4 ) +A1(A5+A6)[-6m2F2+(-\0m-4m3)F2+20m2F4], h3 = - j m(l - m 2 | r ° + 8(A 6 + A5)julR£,]+ A 6 (3 + m2)O0 - 6A6m02 + 2A6m(P2 +P2) + A6(3 + m2)P4 - 6A6mP6 +A 7 (A 5 + A 6)[2m£, - (m + m 3)f, + (3 + 9 m 2 ) £ 3 -30mE5] +A 7 (A 5 + A 6)[(-3 + m2)F, + 2 m 2 ^ + (27m + 3m 3)F 3 -30m 2F 5], h4=-±m(\- m2 \T° + 8(A 6 + A,)^Rs'\- A 6 ( l + 3m 2)O 0 + 6A6m02 90 - 2A6m(F> + P2) - A 6 (1 + 3m 2)F 4 + 6A6mP6 +A 7 (A 5 + A6)[(3m3 -m)F, -2mE, - (3 +27 m2)E3 +30mE5] +A 7 (A 5 + A 6 ) [ -2m 2 F, + (\ +m2)F]-(2\m + 9m3)F3+30m2F5], h5 = A6(4 + 2m2)0]-7A6m03 + A6mP, + A6mP3 + A 6 (4 + 2m2)P5 -7A6mP1 - A 7 (A 5 + A6)[(6m + 4m3 )E2 - (16 + 36m2 ) £ 4 + 42mF6 ] +A 7 (A 5 + A 6)[(-8 + 2m2)F2 + (44m + 8m 3)F 4 -42m2F6], h6 = -A6(2 + 4m2)0] + 7A6m<93 - A 6m/^ - A6mP3 - A 6 (2 + 4m2)P5 + 7A 6mP 7 +A 7(A 6+A 5)[(2m + 8m 3 )F 2 - (8 + 4 4 m 2 ) F 4 + 4 2 m £ 6 ] +A 7 (A 6 +A5)[(4 + 6m 2)F 2 - (36m + 16m3)F4 + 42m 2 F 6 ], h, = A 6 (5 + 3m 2)0, - 8 A 6 m 0 4 + A6(5 + 3m2)P6 - 8A 6 mP 8 +A 7 (A 6 +A5)[-(15m + 9m 3)F 3 + (25 + 55m2)ES -SdmE,j +A 7 (A 6 +A5)[-(15 + 9m 2)F 3 + (65m +15m3)F5 -56m 2 F 7 ] , \ = -A6(3 + 5m2)02 + 8A6m<94 - A 6 (3 + 5m 2)P 6 +8A 6mP 8 +A 7 (A 6 +A 5 )[(9m +15m3 )E3 -(15 + 65m2 )E5 + 56m£ 7 ] +A 7 (A 6 +A 5 )[(9 + 15m 2)F 3-(55m + 25m 3 )F 5 +56m 2 F 7 j , U„ = nA6mOn - A 6 [n + 3 + (n + 5)m2 ]o„+2 + (n + 8)A 6mO„ + 4 + «A 6mF„ + 4 - A6[n + 3 + (n + 5)m 2]P„+ 6 - (« + 8)A 6mP n + 8 - n(n + 1)A7 (A 6 +A 5 )m 2 F n + 1 +A 7 (A 6 +A 5)(« + 3)[(2n + 3)m + (n + 5)m2 ] F „ + 3 - A 7 ( A 6 +A5)(n + 5)[w + 3 + (2n + 13)m2 ]En+5 +(n + 7)(n + 8)A 7 (A 6 +A5)mE, - n(n + 1)A7 (A 5 +A 5)mF„ + 1 +A 7 (A 6 +A5)(n + 3)[n + 3 + (2n + 5)m2 ]F n + 3 + A 7 (A 6 + A 5 ) { - (n + 5)[(2* +1 l)m + (n + 5)m3 ]F„ + 5 + (H + 7)(« + 8)m 2 F„ + 7 }, F„ = -nA6mOn + A 6 [* + 3 + (n + 3)m2 ]o„+2 - ( * + 8)A 6mO n + 4 ~ nA6mP„+4 + A6[n + 5 + ( / i + 3)m2]/>B+6 - ( w + 8)A 6mP„ + 8 + n(n + 1)A7 (A 6 +A 5 )m 2 F„ + 1 - A 7 (A 6 +A5)(n + 3)[(2n + 5)m + (n + 3)m2 ] F „ + 3 - A 7 ( A 6 +A5)(n + 5)[w + 5 + (2n +1 l)m 2 ]F„ + 5 - (« + 7)(« + 8)A 7 (A 6 +A 5)m£„ 91 n+3 + n{n + 1)A 7 (A 6 +A 5)mF n + 1 - A 7 (A 6 +A5)(w + 3)[n + 5 + (2n + 3)m2 ]F„ . +A 7 (A 6 +A5){(« + 5)[(2n +13)m + (w + 3)m 3 ]F„ + 5 - ( « + 7)(w + 8)m 2F n + 7}. (A.7) Appendix B In the (x\, ^-coordinate system, the equation of ellipse is: 2 2 Xj^ __ _ , J. 1.1 ' a~ b where a and & are the semi-axes of the ellipse. Write the variable xj and its derivatives with respect to x\ as a function of x\ and X2 in the first quadrant: (B.l) x2=b,\—Y V o 6 2 x, a x, a2 x 2 Thereafter, the curvature can be written as: J _ a2 x\ (l + x 2 2 J b* x ,2 ^ 2 ' 1+ 4 2 v a x2 j Note that z + z z-z a = Ri\ + m), b = R(\-m), z = R 2 * 2 / The curvature is expressed in terms of the complex variable £ as follows: 1 a ( Z - z ) ' 4 - ( z + z ) 2 3/2 (1-m 2 ) 1 R m\4)lR 3 -(B.2) (B.3) (B.4) (B.5) Appendix C For any integer k, I" (l + b" sin 2 e)V1 sm(k6)dG = 0, | " (l + 6* sin 2 r?)"3/2 cos(2& +1)0 dd = 0, (C.l) 92 Therefore, (l + b' sin 2 6>)~3/2 cos[2(ifc + 1)]0 d6 = f (l + sin2 0)"3'2 [cos lk9 cos 20 - sin Ike sin 20] J0 = { ' ( l + 6* sin 2 0)~3/2 cos 2*0 cos 20^0—^- | + sin 2 0) (^l + fc* sin 2 0 ) c o s 2 £ 0 J0 I* + sin2 f?) [cos2(£ + l)0 + cos2(£-l)0]o'0 4k + 2k_l^1 + y s . n 2 ^"3/2 c o s ( 2 ^ ^ } d Q (C2) 6 Now define r ( l + 6*s in 2 0) 3 / 2 which leads to the relation I2k = I_2k , and Eq. (C.2) becomes 1 Ant \_3/2 I2k =— £ (l + 6*sin 20) cos(2ke)de, (C.3) _ 1 + 2*:/ \ 4fc + 2ta' Using the Fourier series in the complex form, (l + b' sin 2 0)"3/2 =I0+±Ilk(e'2*e + e-2**). (C.5) For a large integer the right-hand side of Eq. (C.5) approaches a geometric series. To find the ratio of this geometric series, assume that for large k, ^ P = 7- (C6) *2k For large k, Eq. (C.4) can be approximately reduced to: n-r 1 — _______ j U ~ l2(k+\) + h{k-\) ^ ^ l2k • Combined with Eq. (C.6), the above expression leads to 0 ^ 2 + l - i ^ ^ . (C.7) b Noting the expression for b" in Eq. (3.40), the solution of Eq. (C.7) is 93 2 + b' b* '2 + b'V V b* j •1 = m < 1, (C.8) which implies that the geometric series is convergent. Now rewrite Eq. (C.5) as: (l + b* sin 2 d)-yi = I0 + £ I2k(e'2k° + e ^ ) + £ 7 2 ( J + t ) [ e " ™ ' + e " * ™ ' ] and consider the third term on the right-hand side as a geometric series for large integer J (approximately). Therefore, J-I (l + b' sin 2 G)V2 = 70 + 7 2 t (e i 2 W + e"'2*")+ /. = / o+E / 2 , ( e ' 2 " + e - ' 2 - ) + / : • + \-meae I-me" 28 k=\ e i 2 J g + e - i 2 J g _ w | e . 2 ( J - l ) 0 + e i 2 ( W ) g -' l -me i 2 ")( l-me- i 2 ( ? ) " • (C9) Appendix D Consider the following equations (\S\ = 1): 1 TT^ ( + i 7 f l V 1 • 1 1-me Therefore ±i20 n=0 m (£)//? 1-me" 20 = X ( m e - 2 ^ . (D . l ) n=0 ( l -we i 2 e ) ( l -we" i 2 e ) = E ( - e - 2 ^ | ; ( ^ 2 ^ | ; ( m e - - ) " n=0 n=0 n=0 n=0 n=l 1-m2 1 f^-i / _ i ? « V 1 n=0 l - m 2 X ^ e - - ^ ^ X ( - e " - ) " + l + X (me- ) " (l-m 2) 2 I t J [w" + (l _ m2 \im" Yln6 +1 + J (me129)" (D.2) 94
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Analytical and finite element solutions for an elastic...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Analytical and finite element solutions for an elastic matrix with two-dimensional nanoscale inhomogeneities Tian, Lian 2006
pdf
Notice for Google Chrome users:
If you are having trouble viewing or searching the PDF with Google Chrome, please download it here instead.
If you are having trouble viewing or searching the PDF with Google Chrome, please download it here instead.
Page Metadata
Item Metadata
Title | Analytical and finite element solutions for an elastic matrix with two-dimensional nanoscale inhomogeneities |
Creator |
Tian, Lian |
Publisher | University of British Columbia |
Date Issued | 2006 |
Description | A new frontier of research at the interface of material science and mechanics of solids has emerged with the current focus on the development of nanomaterials such as nanotubes, nanowires and nanoparticles. As the dimensions of a structure approach the nanoscale, its properties can be size-dependent. The classical continuum theory, however, does not admit an intrinsic size, and is not applicable to the analysis of nanoscale materials and structures. Mechanics of nanomaterials-based composites can be understood by incorporating the effects of surface and interfacial energy. A fundamental problem in the study of behaviour of such materials is the examination of size-dependent elastic field of an elastic matrix with nanoscale inhomogeneities. Classical inhomogeneity problems have a rich history since the celebrated work of Eshelby. However, the classical solutions cannot be applied to study nanoscale inhomogeneity problems and new solutions accounting for surface/interface energy have to be derived. This thesis therefore presents an analytical scheme and a finite element formulation to study the size-dependent elastic field of an elastic matrix containing two-dimensional nanoscale inhomogeneities. The Gurtin-Murdoch surface/interface elasticity model is applied to incorporate the surface/interface energy effects. By using the complex potential technique of Muskhelishvili, a closed-form analytical solution is obtained for the elastic field of a nanoscale circular inhomogeneity in an infinite matrix under arbitrary remote loading and a uniform eigenstrain. In the case of an elliptical inhomogeneity, the analytic potential functions are obtained approximately. A new finite element formulation that takes into account the surface stress effects is presented. Elastic field is found to depend on the characteristic dimensions of the inhomogeneity, surface elastic constants and surface residual stress. A striking feature of the new solutions is the existence of singular elastic fields below some dimensions of the inhomogeneity. This phenomenon requires careful further investigation. Eshelby tensor of a nanoscale circular inhomogeneity in an infinite matrix due to a uniform eigenstrain is uniform but becomes size-dependent; however, the tensor is size-dependent and non-uniform in the case of an elliptical inhomogeneity. |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Date Available | 2010-01-16 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0080742 |
URI | http://hdl.handle.net/2429/18383 |
Degree |
Master of Applied Science - MASc |
Program |
Mechanical Engineering |
Affiliation |
Applied Science, Faculty of Mechanical Engineering, Department of |
Degree Grantor | University of British Columbia |
GraduationDate | 2006-11 |
Campus |
UBCV |
Scholarly Level | Graduate |
AggregatedSourceRepository | DSpace |
Download
- Media
- 831-ubc_2006-0677.pdf [ 4.1MB ]
- Metadata
- JSON: 831-1.0080742.json
- JSON-LD: 831-1.0080742-ld.json
- RDF/XML (Pretty): 831-1.0080742-rdf.xml
- RDF/JSON: 831-1.0080742-rdf.json
- Turtle: 831-1.0080742-turtle.txt
- N-Triples: 831-1.0080742-rdf-ntriples.txt
- Original Record: 831-1.0080742-source.json
- Full Text
- 831-1.0080742-fulltext.txt
- Citation
- 831-1.0080742.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
data-media="{[{embed.selectedMedia}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
https://iiif.library.ubc.ca/presentation/dsp.831.1-0080742/manifest