Mode I and mixed mode crack-tip fields in strain gradient plasticity

Strain gradients develop near the crack-tip of Mode I or mixed mode cracks. A finite strain version of the phenomenological strain gradient plasticity theory of Fleck-Hutchinson (2001) is used here to quantify the effect of the material length scales on the crack-tip stress field for a sharp stationary crack under Mode I and mixed mode loading. It is found that for material length scales much smaller than the scale of the deformation gradients, the predictions converge to conventional elastic-plastic solutions. For length scales sufficiently large, the predictions converge to elastic solutions. Thus, the range of length scales over which a strain gradient plasticity model is necessary is identified. The role of each of the three material length scales, incorporated in the multiple length scale theory, in altering the near-tip stress field is systematically studied in order to quantify their effect.     

A finite deformation theory of strain gradient plasticity

Plastic deformation exhibits strong size dependence at the micron scale, as observed in micro-torsion, bending, and indentation experiments. Classical plasticity theories, which possess no internal material lengths, cannot explain this size dependence. Based on dislocation mechanics, strain gradient plasticity theories have been developed for micron-scale applications. These theories, however, have been limited to infinitesimal deformation, even though the micro-scale experiments involve rather large strains and rotations. In this paper, we propose a finite deformation theory of strain gradient plasticity. The kinematics relations (including strain gradients), equilibrium equations, and constitutive laws are expressed in the reference configuration. The finite deformation strain gradient theory is used to model micro-indentation with results agreeing very well with the experimental data. We show that the finite deformation effect is not very significant for modeling micro-indentation experiments.     

A unified treatment of strain gradient plasticity

A theoretical framework is presented that has potential to cover a large range of strain gradient plasticity effects in isotropic materials. Both incremental plasticity and viscoplasticity models are presented. Many of the alternative models that have been presented in the literature are included as special cases. Based on the expression for plastic dissipation, it is in accordance with Gurtin (J. Mech. Phys. Solids 48 (2000) 989; Int. J. Plast. 19 (2003) 47) argued that the plastic flow direction is governed by a microstress q_ij and not the deviatoric Cauchy stress σ_ij′ that has been assumed by many others. The structure of the governing equations is of second order in the displacements and the plastic strains which makes it comparatively easy to implement in a finite element programme. In addition, a framework for the formulation of consistent boundary conditions is presented. It is shown that there is a close connection between surface energy of an interface and boundary conditions in terms of plastic strains and moment stresses. This should make it possible to study boundary layer effects at the interface between grains or phases. Consistent boundary conditions for an expanding elastic-plastic boundary are as well formulated. As examples, biaxial tension of a thin film on a thick substrate, torsion of a thin wire and a spherical void under remote hydrostatic tension are investigated.     

Grain boundary interface mechanics in strain gradient crystal plasticity

Interactions between dislocations and grain boundaries play an important role in the plastic deformation of polycrystalline metals. Capturing accurately the behaviour of these internal interfaces is particularly important for applications where the relative grain boundary fraction is significant, such as ultra fine-grained metals, thin films and micro-devices. Incorporating these micro-scale interactions (which are sensitive to a number of dislocation, interface and crystallographic parameters) within a macro-scale crystal plasticity model poses a challenge. The innovative features in the present paper include (i) the formulation of a thermodynamically consistent grain boundary interface model within a microstructurally motivated strain gradient crystal plasticity framework, (ii) the presence of intra-grain slip system coupling through a microstructurally derived internal stress, (iii) the incorporation of inter-grain slip system coupling via an interface energy accounting for both the magnitude and direction of contributions to the residual defect from all slip systems in the two neighbouring grains, and (iv) the numerical implementation of the grain boundary model to directly investigate the influence of the interface constitutive parameters on plastic deformation. The model problem of a bicrystal deforming in plane strain is analysed. The influence of dissipative and energetic interface hardening, grain misorientation, asymmetry in the grain orientations and the grain size are systematically investigated. In each case, the crystal response is compared with reference calculations with grain boundaries that are either ‘microhard’ (impenetrable to dislocations) or ‘microfree’ (an infinite dislocation sink).     

Generalizing J 2 flow theory: Fundamental issues in strain gradient plasticity

It has not been a simple matter to obtain a sound extension of the classical J₂ flow theory of plasticity that in- corporates a dependence on plastic strain gradients and that is capable of capturing size-dependent behaviour of metals at the micron scale. Two classes of basic extensions of clas- sical J₂ theory have been proposed: one with increments in higher order stresses related to increments of strain gradi- ents and the other characterized by the higher order stresses themselves expressed in terms of increments of strain gra- dients. The theories proposed by Muhlhaus and Aifantis in 1991 and Fleck and Hutchinson in 2001 are in the first class, and, as formulated, these do not always satisfy ther- modynamic requirements on plastic dissipation. On the other hand, theories of the second class proposed by Gudmundson in 2004 and Gurtin and Anand in 2009 have the physical deficiency that the higher order stress quantities can change discontinuously for bodies subject to arbitrarily small load changes. The present paper lays out this background to the quest for a sound phenomenological extension of the rate- independent J₂ flow theory of plasticity to include a de- pendence on gradients of plastic strain. A modification of the Fleck-Hutchinson formulation that ensures its thermo- dynamic integrity is presented and contrasted with a compa- rable formulation of the second class where in the higher or- der stresses are expressed in terms of the plastic strain rate. Both versions are constructed to reduce to the classical J₂ flow theory of plasticity when the gradients can be neglected and to coincide with the simpler and more readily formulated J₂ deformation theory of gradient plasticity for deformation histories characterized by proportional straining.     

Energetic dislocation interactions and thermodynamical aspects of strain gradient crystal plasticity theories

This paper focuses on the unification of two frequently used and apparently different strain gradient crystal plasticity frameworks: (i) the physically motivated strain gradient crystal plasticity models proposed by Evers et al. [2004a. Non-local crystal plasticity model with intrinsic SSD and GND effects. Journal of the Mechanics and Physics of Solids 52, 2379-2401; 2004b. Scale dependent crystal plasticity framework with dislocation density and grain boundary effects. International Journal of Solids and Structures 41, 5209-5230] and Bayley et al. [2006. A comparison of dislocation induced back stress formulations in strain gradient crystal plasticity. International Journal of Solids and Structure 43, 7268-7286; 2007. A three dimensional dislocation field crystal plasticity approach applied to miniaturized structures. Philosophical Magazine 87, 1361-1378] (here referred to as Evers-Bayley type models), where a physical back stress plays the most important role and which are further extended here to deal with truly large deformations, and (ii) the thermodynamically consistent strain gradient crystal plasticity model of Gurtin (2002-2008) (here referred to as the Gurtin type model), where the energetic part of a higher order micro-stress is derived from a non-standard free energy function. The energetic micro-stress vectors for the Gurtin type models are extracted from the definition of the back stresses of the improved Evers-Bayley type models. The possible defect energy forms that yield the derived physically based micro-stresses are discussed. The duality of both type of formulations is shown further by a comparison of the micro-boundary conditions. As a result, this paper provides a direct physical interpretation of the different terms present in Gurtin's model.     

A gradient theory of small-deformation, single-crystal plasticity that accounts for GND-induced interactions between slip systems

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 densities of (geometrically necessary) edge and screw dislocations, densities that describe the accumulation of dislocations, and densities that characterize forest hardening. The form of the forest densities is based on an explicit kinematical expression for the normal Burgers vector on a slip plane.     

A homogenization theory of strain gradient single crystal plasticity and its finite element discretization

In this study, a homogenization theory based on the Gurtin strain gradient formulation and its finite element discretization are developed for investigating the size effects on macroscopic responses of periodic materials. To derive the homogenization equations consisting of the relation of macroscopic stress, the weak form of stress balance, and the weak form of microforce balance, the Y-periodicity is used as additional, as well as standard, boundary conditions at the boundary of a unit cell. Then, by applying a tangent modulus method, a set of finite element equations is obtained from the homogenization equations. The computational stability and efficiency of this finite element discretization are verified by analyzing a model composite. Furthermore, a model polycrystal is analyzed for investigating the grain size dependence of polycrystal plasticity. In this analysis, the micro-clamped, micro-free, and defect-free conditions are considered as the additional boundary conditions at grain boundaries, and their effects are discussed.     

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.     

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.     

Nanoindentation and the indentation size effect: Kinetics of deformation and strain gradient plasticity

A study of the indentation size effect (ISE) in aluminum and alpha brass is presented. The study employs rate effects to examine the fundamental mechanisms responsible for the ISE. These rate effects are characterized in terms of the rate sensitivity of the hardness, , where H is the hardness and is an effective strain rate in the plastic volume beneath the indenter. can be measured using indentation creep, load relaxation, or rate change experiments. The activation volume V∗, calculated based on which can traditionally be used to compare rate sensitivity data from a hardness test to conventional uniaxial testing, is calculated. Using materials with different stacking fault energy and specimens with different levels of work hardening, we demonstrate how increasing the dislocation density affects V∗; these effects may be taken as a kinetic signature of dislocation strengthening mechanisms. We noticed both H and exhibit an ISE. The course of V∗ vs. H as a result of the ISE is consistent with the course of testing specimens with different level of work hardening. This result was observed in both materials. This suggests that a dislocation mechanism is responsible for the ISE. When the results are fitted to a strain gradient plasticity model, the data at deep indents (microhardness and large nanoindentation) exhibit a straight-line behavior closely identical to literature data. However, for shallow indents (nanoindentation data), the slope of the line severely changes, decreasing by a factor of 10, resulting in a “bilinear behavior”.     

A nonlocal strain gradient plasticity theory for finite deformations

Strain gradient plasticity for finite deformations is addressed within the framework of nonlocal continuum thermodynamics, featured by the concepts of (nonlocality) energy residual and globally simple material. The plastic strain gradient is assumed to be physically meaningful in the domain of particle isoclinic configurations (with the director vector triad constant both in space and time), whereas the objective notion of corotational gradient makes it possible to compute the plastic strain gradient in any domain of particle intermediate configurations. A phenomenological elastic–plastic constitutive model is presented, with mixed kinematic/isotropic hardening laws in the form of PDEs and related higher order boundary conditions (including those associated with the moving elastic/plastic boundary). Two fourth-order projection tensor operators, functions of the elastic and plastic strain states, are shown to relate the skew-symmetric parts of the Mandel stress and back stress to the related symmetric parts. Consistent with the thermodynamic restrictions therein derived, the flow laws for rate-independent associative plasticity are formulated in a six-dimensional tensor space in terms of symmetric parts of Mandel stresses and related work-conjugate generalized plastic strain rates. A simple shear problem application is presented for illustrative purposes.     

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.     

The bimodal theory of plasticity: A form-invariant generalisation

The bimodal plasticity model of fibre-reinforced materials is currently available and applicable only in association with thin-walled fibrous composites containing a family of straight fibres which are conveniently assumed parallel with the x₁-axis of an appropriately chosen Cartesian co-ordinate system. Based on reliable experimental evidence, the model suggests that plastic slip in the composite operates in two distinct modes; the so-called matrix dominated mode (MDM) which depends on a matrix yield stress, and the fibre dominated mode (FDM) which depends also on the fibre yield stress. Each mode is activated by different states of applied stress, has its own yield surface (or surfaces) in the stress space and has its own segment on the overall yield surface of the composite. This paper employs theory of tensor representations and produces a form-invariant generalisation of both modes of the model. This generalisation furnishes the model with direct applicability to relevant plasticity problems, regardless of the shape of the fibres or the orientation of the co-ordinate system. It thus provides a proper mathematical foundation that underpins important physical concepts associated with the model while it also elucidates several technical relevant issues. A most interesting of those issues is the revelation that activation of the MDM plastic regime is possible only if the applied stress state allows the fibres to act like they are practically inextensible. Moreover, activation of the more dominant, between the two MDM plastic slip branches is possible only if conditions of material incompressibility hold, in addition to the implied condition of fibre inextensibility. A direct mathematical connection is thus achieved between basic, experimentally verified concepts of the bimodal plasticity model and a relevant mathematical model originated earlier from the theory of ideal fibre-reinforced materials. An additional issue of discussion involves the number of independent yield stress parameters that the bimodal theory needs to take into consideration. Moreover, an analytical expression is provided of a relatively simple mathematical surface that possesses all known features of the FDM yield surface; currently captured with the aid of both experimental and computational means. The present study is guided by the existing relevant experimental evidence which, however, is principally associated with the plastic behaviour of solids reinforced by strong fibres. Nevertheless, several of the outlined developments are expected to be applicable to composite materials containing a single family of more compliant or even weak fibres.     

Continuum theory of dislocations revisited

The paper continues the discussion of continuum theory of dislocations suggested by Berdichevsky and Sedov (PMM 31(6): 981–1000, 1967). The major new points are: the choice of energy, the variational form of the governing dynamical equations, the variational principle for the final plastic state.     

Formulation of strain gradient plasticity with interface energy in a consistent thermodynamic framework

In this work, the strain gradient formulation is used within the context of the thermodynamic principle, internal state variables, and thermodynamic and dissipation potentials. These in turn provide balance of momentum, boundary conditions, yield condition and flow rule, and free energy and dissipative energies. This new formulation contributes to the following important related issues: (i) the effects of interface energy that are incorporated into the formulation to address various boundary conditions, strengthening and formation of the boundary layers, (ii) nonlocal temperature effects that are crucial, for instance, for simulating the behavior of high speed machining for metallic materials using the strain gradient plasticity models, (iii) a new form of the nonlocal flow rule, (iv) physical bases of the length scale parameter and its identification using nano-indentation experiments and (v) a wide range of applications of the theory. Applications to thin films on thick substrates for various loading conditions and torsion of thin wires are investigated here along with the appropriate length scale effect on the behavior of these structures. Numerical issues of the theory are discussed and results are obtained using Matlab and Mathematica for the nonlinear ordinary differential equations (NODE) which constitute the boundary value problem.This study reveals that the micro-stress term has an important effect on the development of the boundary layers and hardening of the material at both hard and soft interface boundary conditions in thin films. The interface boundary conditions are described by the interfacial length scale and interfacial strength parameters. These parameters are important to control the size effect and hardening of the material. For more complex geometries the generalized form of the boundary value problem using the nonlocal finite element formulation is required to address the problems involved.     

Parameter identification of advanced plastic strain rate potentials and impact on plastic anisotropy prediction

In the work presented in this paper, several strain rate potentials are examined in order to analyze their ability to model the initial stress and strain anisotropy of several orthotropic sheet materials. Classical quadratic and more advanced non-quadratic strain rate potentials are investigated in the case of FCC and BCC polycrystals. Different identifications procedures are proposed, which are taking into account the crystallographic texture and/or a set of mechanical test data in the determination of the material parameters.     

On the homogenization of metal matrix composites using strain gradient plasticity

The homogenized response of metal matrix composites（MMC） is studied using strain gradient plasticity.The material model employed is a rate independent formulation of energetic strain gradient plasticity at the micro scale and conventional rate independent plasticity at the macro scale. Free energy inside the micro structure is included due to the elastic strains and plastic strain gradients. A unit cell containing a circular elastic fiber is analyzed under macroscopic simple shear in addition to transverse and longitudinal loading. The analyses are carried out under generalized plane strain condition. Micro-macro homogenization is performed observing the Hill-Mandel energy condition,and overall loading is considered such that the homogenized higher order terms vanish. The results highlight the intrinsic size-effects as well as the effect of fiber volume fraction on the overall response curves, plastic strain distributions and homogenized yield surfaces under different loading conditions. It is concluded that composites with smaller reinforcement size have larger initial yield surfaces and furthermore,they exhibit more kinematic hardening.     

Phase-field slip-line theory of plasticity

A variational approach to determine the deformation of an ideally plastic substance is proposed by solving a sequence of energy minimization problems under proper conditions to account for the irreversible character of plasticity. The flow is driven by the local transformation of elastic strain energy into plastic work on slip surfaces, once that a certain energetic barrier for slip activation has been overcome. The distinction of the elastic strain energy into spherical and deviatoric parts is used to incorporate in the model the idea of von Mises plasticity and isochoric plastic strain. This is a “phase field model” because the matching condition at the slip interfaces is substituted by the evolution of an auxiliary phase field that, similar to a damage field, is unitary on the elastic phase and null on the yielded phase. The slip lines diffuse in bands, whose width depends upon a material length-scale parameter.Numerical experiments on representative problems in plane strain give solutions with noteworthy similarities with the results from classical slip-line field theory, but the proposed model is much richer because, accounting for elastic deformations, it can describe the formation of slip bands at the local level, which can nucleate, propagate, widen and diffuse by varying the boundary conditions. In particular, the solution for a long pipe under internal pressure is very different from the one obtainable from the classical macroscopic theory of plasticity. For this case, the location of the plastic bands may be an insight to explain the premature failures that are sometimes encountered during the manufacturing process. This practical example enhances the importance of this new theory based on the mathematical sciences.