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.     

An alternative treatment of phenomenological higher-order strain-gradient plasticity theory

Phenomenological higher-order strain-gradient plasticity is here presented through a formulation inspired by previous work for strain-gradient crystal plasticity. A physical interpretation of the phenomenological yield condition that involves an effect of second gradient of the equivalent plastic strain is discussed, applying a dislocation theory-based consideration. Then, a differential equation for the equivalent plastic strain-gradient is introduced as an additional governing equation. Its weak form makes it possible to deduce and impose extra boundary conditions for the equivalent plastic strain. A connection between the present treatment and strain-gradient theories based on an extended virtual work principle is discussed. Furthermore, a numerical implementation and analysis of constrained simple shear of a thin strip are presented.     

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 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 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 theory of strain-gradient plasticity for isotropic, plastically irrotational materials. Part II: Finite deformations

This paper generalizes to finite deformations our companion paper [Gurtin, M.E., Anand, L., 2004. A theory of strain-gradient plasticity for isotropic, plastically irrotational materials. Part I: Small deformations. Journal of the Mechanics and Physics of Solids, submitted]. Specifically, we develop a gradient theory for finite-deformation isotropic viscoplasticity in the absence of plastic spin. The theory is based on the Kröner–Lee decomposition F = F^eF^p of the deformation gradient into elastic and plastic parts; a system of microstresses consistent with a microforce balance; a mechanical version of the second law that includes, via microstresses, work performed during viscoplastic flow; a constitutive theory that allows:

• the microstresses to depend on D^p, the gradient of the plastic stretching,

• the free energy ψ to depend on the Burgers tensor G = F^pCurlF^p.

The microforce balance when augmented by constitutive relations for the microstresses results in a nonlocal flow rule in the form of a tensorial second-order partial differential equation for F^p. The microstresses are strictly dissipative when ψ is independent of the Burgers tensor, but when ψ depends on G the microstresses are partially energetic, and this, in turn, leads to backstresses and (hence) Bauschinger-effects in the flow rule. The typical macroscopic boundary conditions are supplemented by nonstandard microscopic boundary conditions associated with viscoplastic flow, and, as an aid to numerical solution, a weak (virtual power) formulation of the nonlocal flow rule is derived. Finally, the dependences of the microstresses on D^p are shown, analytically, to result in strengthening and possibly weakening of the body induced by viscoplastic flow.     

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.     

Modeling the evolution of crystallographic dislocation density in crystal plasticity

Dislocations are the most important material defects in crystal plasticity, and although dislocation mechanics has long been understood as the underlying physical basis for continuum crystal plasticity formulations, explicit consideration of crystallographic dislocation mechanics has been largely absent in working constitutive models. Here, dislocation density state variables evolve from initial conditions according to equations based on fundamental concepts in dislocation mechanics such as the conservation of Burgers vector in multiplication and annihilation processes. The model is implemented to investigate the polyslip behavior of single-crystal aluminum. The results not only capture the mechanical stress/strain response, but also detail the development of underlying dislocation structure responsible for the plastic behavior.     

A comparison of nonlocal continuum and discrete dislocation plasticity predictions

Discrete dislocation simulations of two boundary value problems are used as numerical experiments to explore the extent to which the nonlocal crystal plasticity theory of Gurtin (J. Mech. Phys. Solids 50 (2002) 5) can reproduce their predictions. In one problem simple shear of a constrained strip is analyzed, while the other problem concerns a two-dimensional model composite with elastic reinforcements in a crystalline matrix subject to macroscopic shear. In the constrained layer problem, boundary layers develop that give rise to size effects. In the composite problem, the discrete dislocation solutions exhibit composite hardening that depends on the reinforcement morphology, a size dependence of the overall stress-strain response for some morphologies, and a strong Bauschinger effect on unloading. In neither problem are the qualitative features of the discrete dislocation results represented by conventional continuum crystal plasticity. The nonlocal plasticity calculations here reproduce the behavior seen in the discrete dislocation simulations in remarkable detail.     

Multiscale modeling of the plasticity in an aluminum single crystal

This paper describes a numerical, hierarchical multiscale modeling methodology involving two distinct bridges over three different length scales that predicts the work hardening of face centered cubic crystals in the absence of physical experiments. This methodology builds a clear bridging approach connecting nano-, micro- and meso-scales. In this methodology, molecular dynamics simulations (nanoscale) are performed to generate mobilities for dislocations. A discrete dislocations numerical tool (microscale) then uses the mobility data obtained from the molecular dynamics simulations to determine the work hardening. The second bridge occurs as the material parameters in a slip system hardening law employed in crystal plasticity models (mesoscale) are determined by the dislocation dynamics simulation results. The material parameters are computed using a correlation procedure based on both the functional form of the hardening law and the internal elastic stress/plastic shear strain fields computed from discrete dislocations. This multiscale bridging methodology was validated by using a crystal plasticity model to predict the mechanical response of an aluminum single crystal deformed under uniaxial compressive loading along the [4 2 1] direction. The computed strain-stress response agrees well with the experimental data.     

A model framework for anisotropic damage coupled to crystal (visco)plasticity

We develop a model framework for anisotropic damage coupled to crystal (visco)plasticity, which is based on the concept of a fictitious (undamaged) configuration. The theoretical setting is that of finite strains, which is natural when studying crystal inelasticity even in the case of actual small strains. It turns out that the evolution law for damage, which reflects degradation in the slip planes and which is the key new relation, bears strong resemblance with the inelastic flow rule. Some numerical results showing qualitatively the anisotropic development of damage concludes the paper.     

A crystal plasticity model that incorporates stresses and strains due to slip gradients

This work is concerned with incorporating the kinematic and stress effects of excess dislocations in a constitutive model for the elastoplastic behavior of crystalline materials. The foundation of the model is a three term multiplicative decomposition of the deformation gradient in which the two classical terms of plastic and elastic deformation are included along with an additional term for long range strain due to the collective effects of excess dislocations. The long range strain is obtained from an assumed density of Volterra edge dislocations and is directly related to gradients in slip. A new material parameter emerges which is the size the region about a continuum point that contributes to long range strains.Using Hookean elasticity, the stress at a point is linearly related to the sum of the elastic plus the long range strain fields. However, the driving force for slip is postulated to be due only to the elastic stress so that the long range stress is a back stress in the constitutive relationship for plastic deformation. A consistent balance of the total deformation rate with the three proposed mechanisms of deformation leads to a set of differential equations that can be solved for the elastic stress, rotation and pressure which then implicitly defines the material state and equilibrium stress. Results from the simulation of a tapered tensile specimen demonstrate that the constitutive model exhibits isotropic and kinematic type hardening effects as well as changes in the pattern of plastic deformation and necking when compared to a material without slip gradient effects.     

Mechanism-based strain gradient crystal plasticity—II. Analysis

In part I of this series (Mechanism-based strain gradient crystal plasticity—I. Theory. J. Mech. Phys. Sol. (2005), accepted for publication), we have proposed a theory of mechanism-based strain gradient crystal plasticity (MSG-CP) to model the effect of inherent anisotropy of a crystal lattice on size-dependent non-uniform plastic deformation at micron and submicron length scales. In the present paper, several example problems are investigated to show how crystal anisotropy is reflected by the MSG-CP theory.     

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 finite element modelling of low cycle fatigue in fcc metals

A new dislocation-based model for low cycle fatigue in fcc metals at a length scale smaller than the feature size of the dislocation structures is presented. It uses the crystal plasticity finite element method and dislocation densities as internal variables. Equations for the dipole distance distribution, for the double cross slip mechanism and a new dislocation multiplication law are introduced, which can predict the emergence of vein and channel structures starting from a randomly perturbed dislocation distribution. The characteristics of these structures in copper and aluminium, as well as the mechanical properties, are compared with experiments. Compared with existing density-based theories, the capability to reproduce dislocation patterning is a significant step forward.     

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.     

Viscoplastic flow rule for dislocation-mediated plasticity from systematic coarse-graining

The plastic response of metals is determined by the collective, coarse-grained dynamics of dislocations, rather than by the dynamics of individual dislocations. The evolution equations at both levels are quite different, for example considering their dependence on the applied mechanical load. On the one hand, the relation between the configurational force and dislocation velocity for individual dislocations is linear. On the other hand, in phenomenological crystal plasticity models, the relation between load and plastic slip is highly non-linear and often taken of power-law form. In this work, it is shown that this difference is justified and a consequence of emergent effects. Previously, an expression for the macroscopic dislocation flux was derived by systematic coarse graining (Kooiman et al., 2015). This expression has been evaluated numerically in this paper. The resulting relation between dislocation flux and applied mechanical load is found to be of power-law form with an exponent 3.7, while the underlying Discrete Dislocation Dynamics has a linear flux–load relation.