The equivalent axisymmetric model for Berkovich indenters in power-law hardening materials

Nix and Gio [Nix, W.D., Gao, H.J., 1998. Indentation size effects in crystalline materials: a law for strain gradient plasticity. Journal of the Mechanics and Physics of Solids 46, 411–425] established an important relation between the micro-indentation hardness and indentation depth for axisymmetric indenters. For the Berkovich indenter, however, this relation requires an equivalent cone angle. Qin et al. [Qin, J., Huang, Y., Xiao, J., Hwang, K.C., 2009. The equivalence of axisymmetric indentation model for three-dimensional indentation hardness. Journal of Materials Research 24, 776–783] showed that the widely used equivalent cone angle from the criterion of equal base area leads to significant errors in micro-indentation, and proposed a new equivalence of equal cone angle for iridium. It is shown in this paper that this new equivalence holds for a wide range of plastic work hardening materials. In addition, the prior equal-base-area criterion does not hold because the Berkovich indenter gives much higher density of geometrically necessary dislocations than axisymmetric indenter. The equivalence of equal cone angle, however, does not hold for Vickers indenter.     

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.     

Analytical and experimental determination of the material intrinsic length scale of strain gradient plasticity theory from micro- and nano-indentation experiments

The enhanced gradient plasticity theories formulate a constitutive framework on the continuum level that is used to bridge the gap between the micromechanical plasticity and the classical continuum plasticity. They are successful in explaining the size effects encountered in many micro- and nano-advanced technologies due to the incorporation of an intrinsic material length parameter into the constitutive modeling. However, the full utility of the gradient-type theories hinges on one's ability to determine the intrinsic material length that scales with strain gradients, and this study aims at addressing and remedying this situation. Based on the Taylor's hardening law, a micromechanical model that assesses a nonlinear coupling between the statistically stored dislocations (SSDs) and geometrically necessary dislocations (GNDs) is used here in order to derive an analytical form for the deformation-gradient-related intrinsic length-scale parameter in terms of measurable microstructural physical parameters. This work also presents a method for identifying the length-scale parameter from micro- and nano-indentation experiments using both spherical and pyramidal indenters. The deviation of the Nix and Gao [Mech. Phys. Solids 46 (1998) 411] and Swadener et al. [J. Mech. Phys. Solids 50 (2002) 681; Scr. Mater. 47 (2002) 343] indentation size effect (ISE) models’ predictions from hardness results at small depths for the case of conical indenters and at small diameters for the case of spherical indenters, respectively, is largely corrected by incorporating an interaction coefficient that compensates for the proper coupling between the SSDs and GNDs during indentation. Experimental results are also presented which show that the ISE for pyramidal and spherical indenters can be correlated successfully by using the proposed model.     

Size effects in the conical indentation of an elasto-plastic solid

The size effect in conical indentation of an elasto-plastic solid is predicted via the Fleck and Willis formulation of strain gradient plasticity (Fleck, N.A. and Willis, J.R., 2009, A mathematical basis for strain gradient plasticity theory. Part II: tensorial plastic multiplier, J. Mech. Phys. Solids, 57, 1045–1057). The rate-dependent formulation is implemented numerically and the full-field indentation problem is analyzed via finite element calculations, for both ideally plastic behavior and dissipative hardening. The isotropic strain-gradient theory involves three material length scales, and the relative significance of these length scales upon the degree of size effect is assessed. Indentation maps are generated to summarize the sensitivity of indentation hardness to indent size, indenter geometry and material properties (such as yield strain and strain hardening index). The finite element model is also used to evaluate the pertinence of the Johnson cavity expansion model and of the Nix–Gao model, which have been extensively used to predict size effects in indentation hardness.     

Size-dependent sharp indentation—I: a closed-form expression of the indentation loading curve

In this paper, a closed-form expression of the size-dependent sharp indentation loading curve has been proposed based on dimensional analysis and the finite deformation Taylor-based nonlocal theory (TNT) of plasticity (Int. J. Plasticity 20 (2004) 831). The key issue is to link the results of FEM based on TNT plasticity with those obtained using conventional FEM by taking as the effective strain gradient, η, that presented in the work of Nix and Gao (J. Mech. Phys. Solids 46 (1998) 411), thus avoiding large-scale finite element computations using strain gradient plasticity theories. Two experiments carried out on 316 stainless-steel and pure titanium have been used to verify the effectiveness of the present analytical model; the results demonstrate that the present analytical expression of the size-dependent indentation loading curve corresponds very well to the experimental indentation loading curve. The empirical constant, α, in the Taylor model estimated from the experimental data has the correct order of magnitude. Also, the results presented in this part can be further applied to establish an analytical framework to extract the plastic properties of metallic materials with sharp indentation on a small scale where the size effect caused by geometrically necessary dislocations is significant. This will be discussed in detail in the second part of the paper.     

The finite deformation theory of Taylor-based nonlocal plasticity

Recent experiments have shown that metallic materials display significant size effect at the micron and sub-micron scales. This has motivated the development of strain gradient plasticity theories, which usually involve extra boundary conditions and possibly higher-order governing equations. We propose a finite deformation theory of nonlocal plasticity based on the Taylor dislocation model. The theory falls into Rice's theoretical framework of internal variables [J Mech Phys Solids 19 (1971) 433], and it does not require any extra boundary conditions. We apply the theory to study the micro-indentation hardness experiments, and it agrees very well with the experimental data over a wide range of indentation depth.     

The role of dissipation and defect energy in variational formulations of problems in strain-gradient plasticity. Part 1: polycrystalline plasticity

A general set of flow laws and associated variational formulations are constructed for small-deformation rate-independent problems in strain-gradient plasticity. The framework is based on the thermodynamically consistent theory due to Gurtin and Anand (J Mech Phys Solids 53:1624–1649, 2005), and includes as variables a set of microstresses which have both energetic and dissipative components. The flow law is of associative type. It is expressed as a normality law with respect to a convex but otherwise arbitrary yield function, or equivalently in terms of the corresponding dissipation function. Two cases studied are, first, an extension of the classical Hill-Mises or J ₂ flow law and second, a form written as a linear sum of the magnitudes of the plastic strain and strain gradient. This latter form is motivated by work of Evans and Hutchinson (Acta Mater 57:1675–1688, 2009) and Nix and Gao (J Mech Phys Solids 46:411–425, 1998), who show that it leads to superior correspondence with experimental results, at least for particular classes of problems. The corresponding yield function is obtained by a duality argument. The variational problem is based on the flow rule expressed in terms of the dissipation function, and the problem is formulated as a variational inequality in the displacement, plastic strain, and hardening parameter. Dissipative components of the microstresses, which are indeterminate, are absent from the formulation. Existence and uniqueness of solutions are investigated for the generalized Hill-Mises and linear-sum dissipation functions, and for various combinations of defect energy. The conditions for well-posedness of the problem depend critically on the choice of dissipation function, and on the presence or otherwise of a defect energy in the plastic strain or plastic strain gradient, and of internal-variable hardening.     

Discrete dislocation plasticity analysis of the wedge indentation of films

The plane strain indentation of single crystal films on a rigid substrate by a rigid wedge indenter is analyzed using discrete dislocation plasticity. The crystals have three slip systems at ±35.3^° and 90^° with respect to the indentation direction. The analyses are carried out for three values of the film thickness, 2, 10 and , and with the dislocations all of edge character modeled as line singularities in a linear elastic material. The lattice resistance to dislocation motion, dislocation nucleation, dislocation interaction with obstacles and dislocation annihilation are incorporated through a set of constitutive rules. Over the range of indentation depths considered, the indentation pressure for the 10 and thick films decreases with increasing contact size and attains a contact size-independent value for contact lengths . On the other hand, for the films, the indentation pressure first decreases with increasing contact size and subsequently increases as the plastic zone reaches the rigid substrate. For the 10 and thick films sink-in occurs around the indenter, while pile-up occurs in the film when the plastic zone reaches the substrate. Comparisons are made with predictions obtained from other formulations: (i) the contact size-independent indentation pressure is compared with that given by continuum crystal plasticity; (ii) the scaling of the indentation pressure with indentation depth is compared with the relation proposed by Nix and Gao [1998. Indentation size effects in crystalline materials: a law for strain gradient plasticity. J. Mech. Phys. Solids 43, 411-423]; and (iii) the computed contact area is compared with that obtained from the estimation procedure of Oliver and Pharr [1992. An improved technique for determining hardness and elastic-modulus using load and displacement sensing indentation experiments, J. Mater. Res. 7, 1564-1583].     

A conventional theory of mechanism-based strain gradient plasticity

There exist two frameworks of strain gradient plasticity theories to model size effects observed at the micron and sub-micron scales in experiments. The first framework involves the higher-order stress and therefore requires extra boundary conditions, such as the theory of mechanism-based strain gradient (MSG) plasticity [J Mech Phys Solids 47 (1999) 1239; J Mech Phys Solids 48 (2000) 99; J Mater Res 15 (2000) 1786] established from the Taylor dislocation model. The other framework does not involve the higher-order stress, and the strain gradient effect come into play via the incremental plastic moduli. A conventional theory of mechanism-based strain gradient plasticity is established in this paper. It is also based on the Taylor dislocation model, but it does not involve the higher-order stress and therefore falls into the second strain gradient plasticity framework that preserves the structure of conventional plasticity theories. The plastic strain gradient appears only in the constitutive model, and the equilibrium equations and boundary conditions are the same as the conventional continuum theories. It is shown that the difference between this theory and the higher-order MSG plasticity theory based on the same dislocation model is only significant within a thin boundary layer of the solid.     

The indenter tip radius effect in micro- and nanoindentation hardness experiments.

Nix and Gao established an important relation between the microindentation hardness and indentation depth. Such a relation has been verified by many microindentation experiments (indentation depths in the micrometer range), but it does not always hold in nanoindentation experiments (indentation depths approaching the nanometer range). Indenter tip radius effect has been proposed by Qu et al. and others as possibly the main factor that causes the deviation from Nix and Gao's relationship. We have developed an indentation model for micro- and nanoindentation, which accounts for two indenter shapes, a sharp, conical indenter and a conical indenter with a spherical tip. The analysis is based on the conventional theory of mechanism-based strain gradient plasticity established from the Taylor dislocation model to account for the effect of geometrically necessary dislocations. The comparison between numerical result and Feng and Nix's experimental data shows that the indenter tip radius effect indeed causes the deviation from Nix-Gao relation, but it seems not be the main factor. The project supported by the National Natural Science Foundation of China (10121202) and the Ministry of Education of China (20020003023)     

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.     

Numerical implementation of a thermo-elastic–plastic constitutive equation in consideration of transformation plasticity in welding

Finite element analysis of welding processes, which entail phase evolution, heat transfer and deformations, is considered in this paper. Attention focuses on numerical implementation of the thermo-elastic–plastic constitutive equation proposed by Leblond et al. [J. Mech. Phys. Solids 34(4) (1986a) 395; J. Mech. Phys. Solids 34(4) (1986b) 411] in consideration of the transformation plasticity. Based upon the multiplicative decomposition of deformation gradient, hyperelastoplastic formulation is borrowed for efficient numerical implementation, and the algorithmic consistent moduli for elastic–plastic deformations including transformation plasticity are obtained in the closed form. The convergence behavior of the present implementation is demonstrated via a couple of numerical examples.     

The coupled effects of plastic strain gradient and thermal softening on the dynamic growth of voids

This paper examines the combined effects of temperature, strain gradient and inertia on the growth of voids in ductile fracture. A dislocation-based gradient plasticity theory [J. Mech. Phys. Solids 47 (1999) 1239, J. Mech. Phys. Solids 48 (2000) 99] is applied, and temperature effects are incorporated. Since a strong size-dependence is introduced into the dynamic growth of voids through gradient plasticity, a cut-off size is then set by the stress level of the applied loading. Only those voids that are initially larger than the cut-off size can grow rapidly. At the early stages of void growth, the effects of strain gradients greatly increase the stress level. Therefore, thermal softening has a strong effect in lowering the threshold stress for the unstable growth of voids. Once the voids start rapid growth, however, the influence of strain gradients will decrease, and the rate of dynamic void growth predicted by strain gradient plasticity approaches that predicted by classical plasticity theories.     

Size-effects in plane strain sheet-necking

A finite strain generalization of the strain gradient plasticity theory by Fleck and Hutchinson (J. Mech. Phys. Solids 49 (2001a) 2245) is proposed and used to study size effects in plane strain necking of thin sheets using the finite element method. Both sheets with rigid grips at the ends and specimens with shear free ends are analyzed. The strain gradient plasticity theory predicts delayed onset of localization when compared to conventional theory, and it depresses deformation localization in the neck. The sensitivity to imperfections is analyzed as well as differently hardening materials.     

A variational model for fracture mechanics: Numerical experiments

In the variational model for brittle fracture proposed in Francfort and Marigo [1998. Revisiting brittle fracture as an energy minimization problem. J. Mech. Phys. Solids 46, 1319-1342], the minimum problem is formulated as a free discontinuity problem for the energy functional of a linear elastic body. A family of approximating regularized problems is then defined, each of which can be solved numerically by a finite element procedure. Here we re-formulate the minimum problem within the context of finite elasticity. The main change is the introduction of the dependence of the strain energy density on the determinant of the deformation gradient. This change requires new, more general existence and Γ-convergence results. The results of some two-dimensional numerical simulations are presented, and compared with corresponding simulations made in Bourdin et al. [2000. Numerical experiments in revisited brittle fracture. J. Mech. Phys. Solids 48, 797-826] for the linear elastic model.     

A comparison between crystal and isotropic strain gradient plasticity theories with accent on the role of the plastic spin

In the small deformation range, we consider crystal and isotropic “higher-order” theories of strain gradient plasticity, in which two different types of size effects are accounted for: (i) that dissipative, entering the model through the definition of an effective measure of plastic deformation peculiar of the isotropic hardening function and (ii) that energetic, included by defining the defect energy (i.e., a function of Nye's dislocation density tensor added to the free energy; see, e.g., [Gurtin, M.E., 2002. A gradient theory of single-crystal viscoplasticity that accounts for geometrically necessary dislocations. J. Mech. Phys. Solids 50, 5–32]). In order to compare the two modellings, we recast both of them into a unified deformation theory framework and apply them to a simple boundary value problem for which we can exploit the Γ-convergence results of [Bardella, L., Giacomini, A., 2008. Influence of material parameters and crystallography on the size effects describable by means of strain gradient plasticity. J. Mech. Phys. Solids 56 (9), 2906–2934], in which the crystal model is made isotropic by imposing that any direction be a possible slip system. We show that the isotropic modelling can satisfactorily approximate the behaviour described by the isotropic limit obtained from the crystal modelling if the former constitutively involves the plastic spin, as in the theory put forward in Section 12 of [Gurtin, M.E., 2004. A gradient theory of small-deformation isotropic plasticity that accounts for the Burgers vector and for dissipation due to plastic spin. J. Mech. Phys. Solids 52, 2545–2568]. The analysis suggests a criterium for choosing the material parameter governing the plastic spin dependence into the relevant Gurtin model.     

Well-posedness of a model of strain gradient plasticity for plastically irrotational materials

The initial boundary value problem corresponding to a model of strain gradient plasticity due to [Gurtin, M., Anand, L., 2005. A theory of strain gradient plasticity for isotropic, plastically irrotational materials. Part I: Small deformations. J. Mech. Phys. Solids 53, 1624–1649] is formulated as a variational inequality, and analysed. The formulation is a primal one, in that the unknown variables are the displacement, plastic strain, and the hardening parameter. The focus of the analysis is on those properties of the problem that would ensure existence of a unique solution. It is shown that this is the case when hardening takes place. A similar property does not hold for the case of softening. The model is therefore extended by adding to it terms involving the divergence of plastic strain. For this extended model the desired property of coercivity holds, albeit only on the boundary of the set of admissible functions.     

Small scale,grain size and substrate effects in nano-indentation experiment of film–substrate systems

The mechanical properties of film–substrate systems have been investigated through nano-indentation experiments in our former paper (Chen, S.H., Liu, L., Wang, T.C., 2005. Investigation of the mechanical properties of thin films by nano-indentation, considering the effects of thickness and different coating–substrate combinations. Surf. Coat. Technol., 191, 25–32), in which Al–Glass with three different film thicknesses are adopted and it is found that the relation between the hardness H and normalized indentation depth h/t, where t denotes the film thickness, exhibits three different regimes: (i) the hardness decreases obviously with increasing indentation depth; (ii) then, the hardness keeps an almost constant value in the range of 0.1–0.7 of the normalized indentation depth h/t; (iii) after that, the hardness increases with increasing indentation depth. In this paper, the indentation image is further investigated and finite element method is used to analyze the nano-indentation phenomena with both classical plasticity and strain gradient plasticity theories. Not only the case with an ideal sharp indenter tip but also that with a round one is considered in both theories. Finally, we find that the classical plasticity theory can not predict the experimental results, even considering the indenter tip curvature. However, the strain gradient plasticity theory can describe the experimental data very well not only at a shallow indentation depth but also at a deep depth. Strain gradient and substrate effects are proved to coexist in film–substrate nano-indentation experiments.     

Modeling of polycrystals using a gradient crystal plasticity theory that includes dissipative micro-stresses

This study investigates thermodynamically consistent dissipative hardening in gradient crystal plasticity in a large-deformation context. A viscoplastic model which accounts for constitutive dependence on the slip, the slip gradient as well as the slip rate gradient is presented. The model is an extension of that due to Gurtin (Gurtin, M. E., J. Mech. Phys. Solids, 52 (2004) 2545–2568 and Gurtin, M. E., J. Mech. Phys. Solids, 56 (2008) 640–662)), and is guided by the viscoplastic model and algorithm of Ekh et al. (Ekh, M., Grymer, M., Runesson, K. and Svedberg, T., Int. J. Numer. Meths Engng, 72 (2007) 197–220) whose governing equations are equivalent to those of Gurtin for the purely energetic case. In contrast to the Gurtin formulation and in line with that due to Ekh et al., viscoplasticity in the present model is accounted for through a Perzyna-type regularization. The resulting theory includes three different types of hardening: standard isotropic hardening is incorporated as well as energetic hardening driven by the slip gradient. In addition, as a third type, dissipative hardening associated with plastic strain rate gradients is included. Numerical computations are carried out and discussed for the large strain, viscoplastic model with non-zero dissipative backstress.     

Influence of material parameters and crystallography on the size effects describable by means of strain gradient plasticity

In the context of single-crystal strain gradient plasticity, we focus on the simple shear of a constrained strip in order to study the effects of the material parameters possibly involved in the modelling. The model consists of a deformation theory suggested and left undeveloped by Bardella [(2007). Some remarks on the strain gradient crystal plasticity modelling, with particular reference to the material length scales involved. Int. J. Plasticity 23, 296–322] in which, for each glide, three dissipative length scales are considered; they enter the model through the definition of an effective slip which brings into the isotropic hardening function the relevant plastic strain gradients, averaged by means of a p-norm. By means of the defect energy (i.e., a function of Nye's dislocation density tensor added to the free energy; see, e.g., Gurtin [2002. A gradient theory of single-crystal viscoplasticity that accounts for geometrically necessary dislocations. J. Mech. Phys. Solids 50, 5–32]), the model further involves an energetic material length scale. The application suggests that two dissipative length scales may be enough to qualitatively describe the size effect of metals at the microscale, and they are chosen in such a way that the higher-order state variables of the model be the dislocation densities. Moreover, we show that, depending on the crystallography, the size effect governed by the defect energy may be different from what expected (based on the findings of [Bardella, L., 2006. A deformation theory of strain gradient crystal plasticity that accounts for geometrically necessary dislocations. J. Mech. Phys. Solids 54, 128–160] and [Gurtin et al. 2007. Gradient single-crystal plasticity with free energy dependent on dislocation densities. J. Mech. Phys. Solids 55, 1853–1878]), leading mostly to some strengthening. In order to investigate the model capability, we also exploit a Γ-convergence technique to find closed-form solutions in the “isotropic limit”. Finally, we analytically show that in the “perfect plasticity” case, should the dissipative length scales be set to zero, the presence of the sole energetic length scale may lead, as in standard plasticity, to non-uniqueness of solutions.