The correlation of the indentation size effect measured with indenters of various shapes

Experimental results are presented which show that the indentation size effect for pyramidal and spherical indenters can be correlated. For a pyramidal indenter, the hardness measured in crystalline materials usually increases with decreasing depth of penetration, which is known as the indentation size effect. Spherical indentation also shows an indentation size effect. However, for a spherical indenter, hardness is not affected by depth, but increases with decreasing sphere radius. The correlation for pyramidal and spherical indenter shapes is based on geometrically necessary dislocations and work-hardening. The Nix and Gao indentation size effect model (J. Mech. Phys. Solids 46 (1998) 411) for conical indenters is extended to indenters of various shapes and compared to the experimental results.     

Non-local crystal plasticity model with intrinsic SSD and GND effects

A strain gradient-dependent crystal plasticity approach is presented to model the constitutive behaviour of polycrystal FCC metals under large plastic deformation. In order to be capable of predicting scale dependence, the heterogeneous deformation-induced evolution and distribution of geometrically necessary dislocations (GNDs) are incorporated into the phenomenological continuum theory of crystal plasticity. Consequently, the resulting boundary value problem accommodates, in addition to the ordinary stress equilibrium condition, a condition which sets the additional nodal degrees of freedom, the edge and screw GND densities, proportional (in a weak sense) to the gradients of crystalline slip. Next to this direct coupling between microstructural dislocation evolutions and macroscopic gradients of plastic slip, another characteristic of the presented crystal plasticity model is the incorporation of the GND-effect, which leads to an essentially different constitutive behaviour than the statistically stored dislocation (SSD) densities. The GNDs, by their geometrical nature of locally similar signs, are expected to influence the plastic flow through a non-local back-stress measure, counteracting the resolved shear stress on the slip systems in the undeformed situation and providing a kinematic hardening contribution. Furthermore, the interactions between both SSD and GND densities are subject to the formation of slip system obstacle densities and accompanying hardening, accountable for slip resistance. As an example problem and without loss of generality, the model is applied to predict the formation of boundary layers and the accompanying size effect of a constrained strip under simple shear deformation, for symmetric double-slip conditions.     

Crystal plasticity model with enhanced hardening by geometrically necessary dislocation accumulation

A strain gradient dependent crystal plasticity approach is used to model the constitutive behaviour of polycrystal FCC metals under large plastic deformation. Material points are considered as aggregates of grains, subdivided into several fictitious grain fractions: a single crystal volume element stands for the grain interior whereas grain boundaries are represented by bi-crystal volume elements, each having the crystallographic lattice orientations of its adjacent crystals. A relaxed Taylor-like interaction law is used for the transition from the local to the global scale. It is relaxed with respect to the bi-crystals, providing compatibility and stress equilibrium at their internal interface. During loading, the bi-crystal boundaries deform dissimilar to the associated grain interior. Arising from this heterogeneity, a geometrically necessary dislocation (GND) density can be computed, which is required to restore compatibility of the crystallographic lattice. This effect provides a physically based method to account for the additional hardening as introduced by the GNDs, the magnitude of which is related to the grain size. Hence, a scale-dependent response is obtained, for which the numerical simulations predict a mechanical behaviour corresponding to the Hall-Petch effect. Compared to a full-scale finite element model reported in the literature, the present polycrystalline crystal plasticity model is of equal quality yet much more efficient from a computational point of view for simulating uniaxial tension experiments with various grain sizes.     

Dislocation core spreading at interfaces between metal films and amorphous substrates

Experiments with transmission electron microscopy have shown that in a strong electron beam the contrast of dislocations may gradually disappear at an incoherent interface between a metal thin film and an amorphous substrate. There are reasons to believe that this phenomenon is caused by radiation-induced dislocation core spreading at the interface. A quantitative model accounting for this effect will be necessary for a better understanding of dislocation structures and plastic deformation in metal thin films. As a first step toward this objective, we develop a number of mathematical solutions for dislocation core spreading at an incoherent interface. For simplicity, we consider screw dislocations, and consider the interface to be characterized by a shear adhesive strength, τ₀, below which no core spreading occurs, and above which spreading takes place in a viscous manner. We determine the final equilibrium core width and the rate of core spreading for single or planar arrays of dislocations in a homogeneous bulk material or at the interface between a thin film and a semi-infinite substrate where the film and substrate may have the same, or different, elastic constants. Some of our solutions are analytic and others are based on an implicit finite difference method with a Gauss-Chebyshev quadrature scheme. The phenomenon of dislocation core spreading is expected to have a dramatic effect on the strength of crystalline films deposited on amorphous substrates.     

Characterization of yielding behavior of polycrystalline metals with single crystal plasticity based on representative characteristic length

Single crystal plasticity based on a representative characteristic length is proposed and introduced into a homogenization approach based on finite element analyses, which are applied to characterization of distinctive yielding behaviors of polycrystalline metals, yield-point elongation, and grain size strengthening. The computational manner for an implicit stress update is derived with the framework of a standard multi-surface plasticity at finite strain, where the evolution of the characteristic lengths are numerically converted from the accumulated slips of all of slip systems by exploiting the mathematical feature of the characteristic length as the intermediate function of the plastic internal variables. Furthermore, a constitutive model for a single crystal reproduces the stress–strain curve divided into three parts. Using two-scale finite element analysis, the macroscopic stress–strain response with yield-point elongation under a situation of low dislocation density is reproduced. Finally, the grain size effect on the yield strength is analyzed with modeling of the grain boundary in the context of the proposed constitutive model and is discussed from both macroscopic and microscopic views.     

Reverse plasticity in single crystal silicon nanospheres

Nanoparticles in the range of 20 to 100 nm in size can be deposited, isolated, and individually probed for their mechanical properties. With a hypersonic plasma particle deposition technique, this has been successfully accomplished for silicon and titanium. We have already shown that silicon nanoparticles are superhard in the 30 to 50 GPa range after work hardening (Gerberich, W.W., Mook, W.M., Perrey, C.R., Carter, C.B., Baskes, M.I., Mukherjee, R., Gidwani, A., Heberlein, J., McMurry, P.H., Girshick, S.L., 2003a. Superhard silicon nanospheres. J. Mech. Phys. Solids 51, 979). At the same time when small nanospheres are compressed, a fraction of the plastic strain is reversed after unloading. Initially, the amount of reverse dislocation motion was small but appeared to accelerate once a threshold strain was reached. The cumulative reverse plastic strain from repeated loading of the same nanosphere appeared to increase from less than 0.04 to approximately 0.4 as cumulative strain increased from 0.2 to 0.6. For large strains then, it appears that a greater amount of plastic strain is recovered after unloading. This can at least partially be understood in terms of the enormous back stress developed at the small scale when dislocations are only a few nm apart. As the ramifications to nanoscopic features on MEMS, micromachines and magnetic recording devices is considerable, it is desirable to understand if a length scale can be developed for such phenomena. In terms of classic dislocation theory an attempt is made. Problems and prospects are discussed with regards to predictive models for hardness and reverse plasticity.     

Atomistic simulations of elastic deformation and dislocation nucleation during nanoindentation

Nanoindentation experiments have shown that microstructural inhomogeneities across the surface of gold thin films lead to position-dependent nanoindentation behavior [Phys. Rev. B (2002), to be submitted]. The rationale for such behavior was based on the availability of dislocation sources at the grain boundary for initiating plasticity. In order to verify or refute this theory, a computational approach has been pursued. Here, a simulation study of the initial stages of indentation using the embedded atom method (EAM) is presented. First, the principles of the EAM are given, and a comparison is made between atomistic simulations and continuum models for elastic deformation. Then, the mechanism of dislocation nucleation in single crystalline gold is analyzed, and the effects of elastic anisotropy are considered. Finally, a systematic study of the indentation response in the proximity of a high angle, high sigma (low symmetry) grain boundary is presented; indentation behavior is simulated for varying indenter positions relative to the boundary. The results indicate that high angle grain boundaries are a ready source of dislocations in indentation-induced deformation.     

Mechanism-based strain gradient crystal plasticity—I. Theory

We have been developing the theory of mechanism-based strain gradient plasticity (MSG) to model size-dependent plastic deformation at micron and submicron length scales. The core idea has been to incorporate the concept of geometrically necessary dislocations into the continuum plastic constitutive laws via the Taylor hardening relation. Here we extend this effort to develop a mechanism-based strain gradient theory of crystal plasticity. In this theory, an effective density of geometrically necessary dislocations for a specific slip plane is introduced via a continuum analog of the Peach-Koehler force in dislocation theory and is incorporated into the plastic constitutive laws via the Taylor relation.     

High strain gradient plasticity associated with wedge indentation into face-centered cubic single crystals: Geometrically necessary dislocation densities

Experimental studies on indentation into face-centered cubic (FCC) single crystals such as copper and aluminum were performed to reveal the spatially resolved variation in crystal lattice rotation induced due to wedge indentation. The crystal lattice curvature tensors of the indented crystals were calculated from the in-plane lattice rotation results as measured by electron backscatter diffraction (EBSD). Nye's dislocation density tensors for plane strain deformation of both crystals were determined from the lattice curvature tensors. The least L²-norm solutions to the geometrically necessary dislocation densities for the case in which three effective in-plane slip systems were activated in the single crystals associated with the indentation were determined. Results show the formation of lattice rotation discontinuities along with a very high density of geometrically necessary dislocations.     

Direct calculations of material parameters for gradient plasticity

Length scale parameters introduced in gradient theories of plasticity are calculated in closed form with a continuum dislocation based theory. The similarity of the governing equations in both models for the evolution of plastic deformation of a constrained thin film makes it possible to identify parameters of the gradient plasticity theory with the dislocation based model. A one-to-one identification is not possible given that gradient plasticity does not account for individual dislocations. However, by comparing the mean plastic deformation across the film thickness we find that the length scale parameter, l, introduced in the gradient plasticity theory depends on the geometry as well as material constants.     

Higher-order stress and grain size effects due to self-energy of geometrically necessary dislocations

The higher-order stress work-conjugate to slip gradient in single crystals at small strains is derived based on the self-energy of geometrically necessary dislocations (GNDs). It is shown that this higher-order stress changes stepwise as a function of in-plane slip gradient and therefore significantly influences the onset of initial yielding in polycrystals. The higher-order stress based on the self-energy of GNDs is then incorporated into the strain gradient plasticity theory of Gurtin [2002. A gradient theory of single-crystal viscoplasticity that accounts for geometrically necessary dislocations. J. Mech. Phys. Solids 50, 5-32] and applied to single-slip-oriented 2D and 3D model crystal grains of size D. It is thus found that the self-energy of GNDs gives a D^-1-dependent term for the averaged resolved shear stress in such a model grain under yielding. Using published experimental data for several polycrystalline metals, it is demonstrated that the D^-1-dependent term successfully explains the grain size dependence of initial yield stress and the dislocation cell size dependence of flow stress in the submicron to several-micron range of grain and cell sizes.     

Modelling the torsion of thin metal wires by distortion gradient plasticity

Under small strains and rotations, we apply a phenomenological higher-order theory of distortion gradient plasticity to the torsion problem, here assumed as a paradigmatic benchmark of small-scale plasticity. Peculiar of the studied theory, proposed about ten years ago by Morton E. Gurtin, is the constitutive inclusion of the plastic spin, affecting both the free energy and the dissipation. In particular, the part of the free energy, called the defect energy, which accounts for Geometrically Necessary Dislocations, is a function of Nye's dislocation density tensor, dependent on the plastic distortion, including the plastic spin. For the specific torsion problem, we implement this distortion gradient plasticity theory into a Finite Element (FE) code characterised by implicit (Backward Euler) time integration, numerically robust and accurate for both viscoplastic and rate-independent material responses. We show that, contrariwise to other higher-order theories of strain gradient plasticity (neglecting the plastic spin), the distortion gradient plasticity can predict some strengthening even if a quadratic defect energy is chosen. On the basis of the results of many FE analyses, concerned with (i) cyclic loading, (ii) switch in the higher-order boundary conditions during monotonic plastic loading, (iii) the use of non-quadratic defect energies, and (iv) the prediction of experimental data, we mainly show that (a) including the plastic spin contribution in a gradient plasticity theory is highly recommendable to model small-scale plasticity, (b) less-than-quadratic defect energies may help in describing the experimental results, but they may lead to anomalous cyclic behaviour, and (c) dissipative (unrecoverable) higher-order finite stresses are responsible for an unexpected mechanical response under non-proportional loading.     

The strength limit in a bio-inspired metallic nanocomposite

Large-scale molecular dynamics simulations are performed to investigate the plastic deformation behavior of a bio-inspired metallic nanocomposite which consists of hard nanosized Ni platelets embedded in a soft Al matrix. The investigation is restricted to an idealized nanocomposite structure with regular platelet distributions in a quasi-two-dimensional geometry under quasi-static loading conditions. This restriction enables us to study size dependent material properties over a wide range of length scales with a fully atomistic resolution of the material and thus without any a priori assumptions of the deformation processes. The simulation results are analyzed with respect to the prevailing deformation mechanisms and their influence on the mechanical properties of the nanocomposite with various geometrical variations. It is found that interfacial sliding contributes significantly to the plastic deformation despite a strong bonding across the interface. Critical for the strength of the nanocomposite is the geometric confinement of dislocation processes in the plastic phase, which strongly depends on the length scale and the morphology of the nanostructure. However, for the smallest structural scales, the softening caused by interfacial sliding prevails, giving rise to a maximum strength.     

Mechanical properties of nanostructure of biological materials

Natural biological materials such as bone, teeth and nacre are nanocomposites of protein and mineral with superior strength. It is quite a marvel that nature produces hard and tough materials out of protein as soft as human skin and mineral as brittle as classroom chalk. What are the secrets of nature? Can we learn from this to produce bio-inspired materials in the laboratory? These questions have motivated us to investigate the mechanics of protein-mineral nanocomposite structure. Large aspect ratios and a staggered alignment of mineral platelets are found to be the key factors contributing to the large stiffness of biomaterials. A tension-shear chain (TSC) model of biological nanostructure reveals that the strength of biomaterials hinges upon optimizing the tensile strength of the mineral crystals. As the size of the mineral crystals is reduced to nanoscale, they become insensitive to flaws with strength approaching the theoretical strength of atomic bonds. The optimized tensile strength of mineral crystals thus allows a large amount of fracture energy to be dissipated in protein via shear deformation and consequently enhances the fracture toughness of biocomposites. We derive viscoelastic properties of the protein-mineral nanostructure and show that the toughness of biocomposite can be further enhanced by the viscoelastic properties of protein.     

Dislocation shielding of a cohesive crack

Dislocation interaction with a cohesive crack is of increasing importance to computational modelling of crack nucleation/growth and related toughening mechanisms in confined structures and under cyclic fatigue conditions. Here, dislocation shielding of a Dugdale cohesive crack described by a rectangular traction-separation law is studied. The shielding is completely characterized by three non-dimensional parameters representing the effective fracture toughness, the cohesive strength, and the distance between the dislocations and the crack tip. A closed form analytical solution shows that, while the classical singular crack model predicts that a dislocation can shield or anti-shield a crack depending on the sign of its Burgers vector, at low cohesive strengths a dislocation always shields the cohesive crack irrespective of the Burgers vector. A numerical study shows the transition in shielding from the classical solution of Lin and Thomson (1986) in the high strength limit to the solution in the low strength limit. An asymptotic analysis yields an approximate analytical model for the shielding over the full range of cohesive strengths. A discrete dislocation (DD) simulation of a large (>10³) number of edge dislocations interacting with a cohesive crack described by a trapezoidal traction-separation law confirms the transition in shielding, showing that the cohesive crack does behave like a singular crack at very high cohesive strengths (∼7 GPa), but that significant deviations in shielding between singular and cohesive crack predictions arise at cohesive strengths around 1GPa, consistent with the analytic models. Both analytical and numerical studies indicate that an appropriate crack tip model is essential for accurately quantifying dislocation shielding for cohesive strengths in the GPa range.     

A deformation theory of strain gradient crystal plasticity that accounts for geometrically necessary dislocations

We propose a deformation theory of strain gradient crystal plasticity that accounts for the density of geometrically necessary dislocations by including, as an independent kinematic variable, Nye's dislocation density tensor [1953. Acta Metallurgica 1, 153-162]. This is accomplished in the same fashion as proposed by Gurtin and co-workers (see, for instance, Gurtin and Needleman [2005. J. Mech. Phys. Solids 53, 1-31]) in the context of a flow theory of crystal plasticity, by introducing the so-called defect energy. Moreover, in order to better describe the strengthening accompanied by diminishing size, we propose that the classical part of the plastic potential may be dependent on both the plastic slip vector and its gradient; for single crystals, this also makes it easier to deal with the “higher-order” boundary conditions. We develop both the kinematic formulation and its static dual and apply the theory to the simple shear of a constrained strip (example already exploited in Shu et al. [2001. J. Mech. Phys. Solids 49, 1361-1395], Bittencourt et al. [2003. J. Mech. Phys. Solids 51, 281-310], Niordson and Hutchinson [2003. Euro J. Mech. Phys. Solids 22, 771-778], Evers et al. [2004. J. Mech. Phys. Solids 52, 2379-2401], and Anand et al. [2005. J. Mech. Phys. Solids 53, 1789-1826]) to investigate what sort of behaviour the new model predicts. The availability of the total potential energy functional and its static dual allows us to easily solve this simple boundary value problem by resorting to the Ritz method.     

Size and temperature effects on the fracture mechanisms of silicon nanowires: Molecular dynamics simulations

We present molecular dynamics simulations of [1 1 0]-oriented Si nanowires (NWs) under a constant strain rate in tension until failure, using the modified embedded-atom-method (MEAM) potential. The fracture behavior of the NWs depends on both temperature and NW diameter. For NWs of diameter larger than 4 nm, cleavage fracture on the transverse (1 1 0) plane are predominantly observed at temperatures below 1000 K. At higher temperatures, the same NWs shear extensively on inclined {1 1 1} planes prior to fracture, analogous to the brittle-to-ductile transition (BDT) in bulk Si. Surprisingly, NWs with diameter less than 4 nm fail by shear regardless of temperature. Detailed analysis reveals that cleavage fracture is initiated by the nucleation of a crack, while shear failure is initiated by the nucleation of a dislocation, both from the surface. While dislocation mobility is believed to be the controlling factor of BDT in bulk Si, our analysis showed that the change of failure mechanism in Si NWs with decreasing diameters is nucleation controlled. Our results are compared with a recent in situ tensile experiment of Si NWs showing ductile failure at room temperature.     

Strain hardening in 2D discrete dislocation dynamics simulations: A new ‘2.5D’ algorithm

The two-dimensional discrete dislocation dynamics (2D DD) method, consisting of parallel straight edge dislocations gliding on independent slip systems in a plane strain model of a crystal, is often used to study complicated boundary value problems in crystal plasticity. However, the absence of truly three dimensional mechanisms such as junction formation means that forest hardening cannot be modeled, unless additional so-called ‘2.5D’ constitutive rules are prescribed for short-range dislocation interactions. Here, results from three dimensional dislocation dynamics (3D DD) simulations in an FCC material are used to define new constitutive rules for short-range interactions and junction formation between dislocations on intersecting slip systems in 2D. The mutual strengthening effect of junctions on preexisting obstacles, such as precipitates or grain boundaries, is also accounted for in the model. The new ‘2.5D’ DD model, with no arbitrary adjustable parameters beyond those obtained from lower scale simulation methods, is shown to predict athermal hardening rates, differences in flow behavior for single and multiple slip, and latent hardening ratios. All these phenomena are well-established in the plasticity of crystals and quantitative results predicted by the model are in good agreement with experimental observations.     

Global load-sharing model for unidirectional hybrid fibre-reinforced composites

A promising strategy to increase the tensile failure strain of carbon fibre-reinforced composites is to hybridise carbon fibres with other, higher-elongation fibres. The resulting increase in failure strain is known as the hybrid effect. In the present article, a global load-sharing model for hybrid composites is developed and used to carry out a parametric study for carbon/glass hybrids. Hybrid effects of up to 15% increase in failure strain are predicted, corresponding reasonably well to literature data. Scatter in the carbon fibre strength is shown to be crucial for the hybrid effect, while the scatter in glass fibre strength is much less important. In contrast to reports in earlier literature, the ratio of failure strains of the two fibres has only a small influence on the hybrid effect. The results provide guidelines for designing optimal hybrid composites.