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.     

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.     

On the formulations of higher-order strain gradient crystal plasticity models

Recently, several higher-order extensions to the crystal plasticity theory have been proposed to incorporate effects of material length scales that were missing links in the conventional continuum mechanics. The extended theories are classified into work-conjugate and non-work-conjugate types. A common feature of the former is that existence of higher-order stresses work-conjugate to gradients of plastic strain is presumed and an extended principle of virtual work involving such an additional virtual work contribution is formulated. Meanwhile, in the latter type, the higher-order stress quantities do not appear explicitly. Instead, rates of crystallographic slip are influenced by back stresses that arise in response to spatial gradients of the geometrically necessary dislocation densities. The work-conjugate type and the non-work-conjugate type of theories have different theoretical backgrounds and very unlike mathematical representations. Nevertheless, both types of theories predict the same kind of material length scale effects. We have recently shown that there exists some equivalency between the two approaches in the special situation of two-dimensional single slip under small deformation. In this paper, the discussion is extended to a more general situation, i.e. the context of multiple and three-dimensional slip deformations.     

Effect of yield surface vertex on crack-tip fields in mode III

An asymptotic crack-tip analysis of stress and strain fields is carried out for an antiplane shear crack (Mode III) based on a corner theory of plasticity. Because of the nonproportional loading history experienced by a material element near the crack tip in stable crack growth, classical flow theory may predict an overly stiff response of the elastic plastic solid, as is the case in plastic buckling problems. The corner theory used here accounts for this anomalous behavior. The results are compared with those of a similar analysis based on the J₂ flow theory of plasticity.     

A study on the bending stiffness of single-walled carbon nanotubes and related issues

Single-walled carbon nanotubes (SWCNTs) are usually modeled as elastic tubes and their bending stiffness D is often related to their axial stretching modulus E (Young's modulus) as in mechanics of materials (i.e. D=EI where I is the moment of inertia of the tube). However, recent studies show that large discrepancies may exist when this relationship is used to predict Young's modulus of carbon nanotubes (CNTs) through bending dominated deformations. In the present paper, the bending stiffness of SWCNTs and some related issues are investigated by the combined use of the molecular-mechanics (M-M) model and the deformation mapping technique. Based on the analysis results, the contradictions mentioned above can be explained well. Furthermore, an analytical expression for the bending stiffness of SWCNTs is also presented. It shows that the bending stiffness of a SWCNT is approximately proportional to the cube of its radius which agrees well with the existing molecular dynamics simulation and continuum theory based results.     

A finite-deformation,gradient theory of single-crystal plasticity with free energy dependent on the accumulation of geometrically necessary dislocations

This paper develops a gradient theory of single-crystal plasticity based on a system of microscopic force balances, one balance for each slip system, derived from the principle of virtual power, and a mechanical version of the second law that includes, via the microscopic forces, work performed during plastic flow. When combined with thermodynamically consistent constitutive relations the microscopic force balances become nonlocal flow rules for the individual slip systems in the form of partial differential equations requiring boundary conditions. Central ingredients in the theory are geometrically necessary edge and screw dislocations together with a free energy that accounts for work hardening through a dependence on the accumulation of geometrically necessary dislocations.     

Size effect in quasibrittle failure: Analytical model and numerical simulations using cohesive zone model

The recent rewriting of the Ba?ant’s size effect law (Morel, 2008) which has suggested the existence of an additional asymptotic regime for intermediate structure sizes is now compared to numerical simulations of fracture of geometrically similar notched structures of different sizes extending over 2.4 decades. The quasibrittle fracture behavior is simulated through cohesive zone model (bilinear softening) using a constant set of cohesive parameters whatever the specimen size D is. The R-curves resulting from the load–displacement responses are estimated and appear as size-independent. On this basis, the different asymptotic regimes expected for the size effect on fracture properties at peak load such as the relative crack length, the resistance to crack growth and the nominal strength are shown in fair agreement with the size effect observed on the results obtained from numerical simulations.     

A gradient theory of single-crystal viscoplasticity that accounts for geometrically necessary dislocations

This study develops a gradient theory of single-crystal plasticity that accounts for geometrically necessary dislocations. The theory is based on classical crystalline kinematics; classical macroforces; microforces for each slip system consistent with a microforce balance; a mechanical version of the second law that includes, via the microforces, work performed during slip; a rate-independent constitutive theory that includes dependences on a tensorial measure of geometrically necessary dislocations. The microforce balances are shown to be equivalent to nonlocal yield conditions for the individual slip systems. The field equations consist of the yield conditions coupled to the standard macroscopic force balance; these are supplemented by classical macroscopic boundary conditions in conjunction with nonstandard boundary conditions associated with slip. As an aid to solution, a weak (virtual power) formulation of the nonlocal yield conditions is derived. To make contact with classical dislocation theory, the microstresses are shown to represent counterparts of the Peach-Koehler force on a single dislocation.     

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.     

A simple theory of strain gradient plasticity based on stress-induced anisotropy of defect diffusion

A new strain gradient plasticity theory is formulated to accommodate more than one material length parameter. The theory is an extension of the classical J₂ flow theory of metal plasticity to the micron scale. Distinctive features of the proposed theory as compared to other existing theories are the simplicities of mathematical formulation, numerical implementation and physical interpretation.     

Size-dependent rigidities of nanosized torsional elements

A theory for the prediction of the size dependence of torsional rigidities of nanosized structural elements is developed. It is shown that, to a very good approximation, the torsional rigidity (D) of a nanosized bar differs from the prediction of standard continuum mechanics (D_c) as (D−D_c)/D_c=Ah₀/a where A is a non-dimensional constant, a is the size scale of the cross-section of the bar and h₀ is a material length equal to the ratio of the surface elastic constant to the bulk elastic constant. The theory developed is compared with direct atomistic calculations (“numerical experiment”) of the torsional rigidity bars made of several FCC metals modeled using the embedded atom method. Very good agreement is obtained between theory and simulation. The framework presented here can aid the development of design methodologies for nanoscale structural elements without the need for full scale atomistic simulations.     

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.     

Size effects in generalised continuum crystal plasticity for two-phase laminates

The solutions of a boundary value problem are explored for various classes of generalised crystal plasticity models including Cosserat, strain gradient and micromorphic crystal plasticity. The considered microstructure consists of a two-phase laminate containing a purely elastic and an elasto-plastic phase undergoing single or double slip. The local distributions of plastic slip, lattice rotation and stresses are derived when the microstructure is subjected to simple shear. The arising size effects are characterised by the overall extra back stress component resulting from the action of higher order stresses, a characteristic length l_c describing the size-dependent domain of material response, and by the corresponding scaling law lⁿ as a function of microstructural length scale, l. Explicit relations for these quantities are derived and compared for the different models. The conditions at the interface between the elastic and elasto-plastic phases are shown to play a major role in the solution. A range of material parameters is shown to exist for which the Cosserat and micromorphic approaches exhibit the same behaviour. The models display in general significantly different asymptotic regimes for small microstructural length scales. Scaling power laws with the exponent continuously ranging from 0 to −2 are obtained depending on the values of the material parameters. The unusual exponent value −2 is obtained for the strain gradient plasticity model, denoted “curl H^p” in this work. These results provide guidelines for the identification of higher order material parameters of crystal plasticity models from experimental data, such as precipitate size effects in precipitate strengthened alloys.     

Scale transition of a higher order plasticity model – A consistent homogenization theory from meso to macro

Standard plasticity models cannot capture the microstructural size effect associated with grain sizes, as well as structural size effects induced by external boundaries and overall gradients. Many higher-order plasticity models introduce a length scale parameter to resolve the latter limitation – microstructural influences are not explicitly account for. This paper adopts two distinct length scales in the formulation, i.e. an intrinsic length scale (l) governing micro-processes such as dislocation pile-up at internal boundaries, as well as the characteristic grain size (L), and aims to unravel the interaction between these two length scales and the characteristic specimen size (H) at the macro level. At the meso-scale, we adopt the strain gradient plasticity model developed in Gurtin (2004) [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] which accounts for the direct influence of grain boundaries. Through a novel homogenization theory, the plasticity model is translated consistently from meso to macro. The two length scale parameters (l and L) manifest themselves naturally at the macro scale, hence capturing both types of size effects in an average sense. The resulting (macro) higher-order model is thermodynamically consistent to the meso model, and has the same structure as a micromorphic continuum. Finally, we consider a bending example for the two limiting cases – microhard and microfree conditions at grain boundaries – and illustrate the excellent match between the meso and homogenized solutions.     

A generalized approach to formulate the consistent tangent stiffness in plasticity with application to the GLD porous material model

It has been shown that the use of the consistent tangent moduli is crucial for preserving the quadratic convergence rate of the global Newton iterations in the solution of the incremental problem. In this paper, we present a general method to formulate the consistent tangent stiffness for plasticity. The robustness and efficiency of the proposed approach are examined by applying it to the isotropic material with J₂ flow plasticity and comparing the performance and the analysis results with the original implementation in the commercial finite element program ABAQUS. The proposed approach is then applied to an anisotropic porous plasticity model, the Gologanu–Leblond–Devaux model. Performance comparison between the consistent tangent stiffness and the conventional continuum tangent stiffness demonstrates significant improvement in convergence characteristics of the overall Newton iterations caused by using the consistent tangent matrix.     

Finite viscoplasticity of amorphous glassy polymers in the logarithmic strain space

The paper outlines a new constitutive model and experimental results of rate-dependent finite elastic–plastic behavior of amorphous glassy polymers. In contrast to existing kinematical approaches to finite viscoplasticity of glassy polymers, the formulation proposed is constructed in the logarithmic strain space and related to a six-dimensional plastic metric. Therefore, it a priori avoids difficulties concerning with the uniqueness of a plastic rotation. The constitutive framework consists of three major steps: (i) A geometric pre-processing defines a total and a plastic logarithmic strain measures determined from the current and plastic metrics, respectively. (ii) The constitutive model describes the stresses and the consistent moduli work-conjugate to the logarithmic strain measures in an analogous structure to the geometrically linear theory. (iii) A geometric post-processing maps the stresses and the algorithmic tangent moduli computed in the logarithmic strain space to their nominal, material or spatial counterparts in the finite deformation space. The analogy between the formulation of finite plasticity in the logarithmic strain space and the geometrically linear theory of plasticity makes this framework very attractive, in particular regarding the algorithmic implementation. The flow rule for viscoplastic strains in the logarithmic strain space is adopted from the celebrated double-kink theory. The post-yield kinematic hardening is modeled by different network models. Here, we compare the response of the eight chain model with the newly proposed non-affine micro-sphere model. Apart from the constitutive model, experimental results obtained from both the homogeneous compression and inhomogeneous tension tests on polycarbonate are presented. Besides the load–displacement data acquired from inhomogeneous experiments, quantitative three-dimensional optical measurements of the surface strain fields are carried out. With regard to these experimental data, the excellent predictive quality of the theory proposed is demonstrated by means of representative numerical simulations.     

A critical state sand model with elastic–plastic coupling

Soil elastic moduli are highly pressure-dependent. Experimental findings have indicated that the elastic shear modulus of sands depends on p^χ, where p is mean principal effective stress and χ is a non-dimensional parameter. χ practically remains unchanged for shear strains less than 10⁻⁵ where the mechanical behavior is purely elastic. However, experiments have revealed that the emergence of plasticity for shear strains larger than 10⁻⁵ provokes a gradual increase in χ. Technically, this observation is an elastic–plastic coupling effect in which plasticity causes to change the elastic characteristics. Here, this issue is considered in hyper-elasticity framework in conjunction with a critical state compatible bounding surface plasticity platform for granular soils. To this aim, constitutive equations linking χ to a proper kinematic hardening parameter are presented. Then, using the proposed approach, a hyper-elastic theory is modified to consider the mentioned elastic–plastic coupling effect in the whole domain of the elastoplastic behavior. Adopting the improved hyper-elasticity necessitates the modification of a number of basic plasticity platform elements. In this regard, dilatancy and plastic hardening modulus of the bounding surface platform are modified. Successful performance of the modified constitutive model is presented against experimental data of loading/unloading triaxial tests.     

Wedge indentation of thin films modelled by strain gradient plasticity

A plane strain study of wedge indentation of a thin film on a substrate is performed. The film is modelled with the strain gradient plasticity theory by Gudmundson [Gudmundson, P., 2004. A unified treatment of strain gradient plasticity. Journal of the Mechanics and Physics of Solids 52, 1379–1406] and analysed using finite element simulations. Several trends that have been experimentally observed elsewhere are captured in the predictions of the mechanical behaviour of the thin film. Such trends include increased hardness at shallow depths due to gradient effects as well as increased hardness at larger depths due to the influence of the substrate. In between, a plateau is found which is observed to scale linearly with the material length scale parameter. It is shown that the degree of hardening of the material has a strong influence on the substrate effect, where a high hardening modulus gives a larger impact on this effect. Furthermore, pile-up deformation dominated by plasticity at small values of the internal length scale parameter is turned into sink-in deformation where plasticity is suppressed for larger values of the length scale parameter. Finally, it is demonstrated that the effect of substrate compliance has a significant effect on the hardness predictions if the effective stiffness of the substrate is of the same order as the stiffness of the film.     

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.