Heterogeneous crystal and poly-crystal plasticity modeling from a transformation field analysis within a regularized Schmid law

The regularized Schmid law (RSL) has recently been proposed as a plastic flow criterion for poly-crystals under the crude assumptions of either uniform stress or uniform strain. We first reconsider this law for application to heterogeneous intra-crystalline plasticity, with reference to a Homogeneous Equivalent Super-Crystal. We then extend the modeling to poly-crystals with the goal to account for both stress and strain heterogeneities within as well as between grains. The transformation field analysis (TFA) is used as the homogenization procedure. This TFA is known to be accurate for materials that can be described as assemblies of plastically homogeneous domains. Otherwise, the estimates of the material effective behavior that result from its application are too stiff. Because stress and strain fields are almost everywhere uniform in laminates, we consider crystal slip organizations into multi-laminate structures. It is demonstrated that laminate layers either parallel to slip planes or normal to slip directions do not contribute to the over-stiffness due to the TFA. Thus, hierarchical multi-laminate (HML) structures are introduced where the successive laminate orientations are taken parallel to the crystal slip planes. It is shown that a conveniently weighted superposition of all the possible plane hierarchies cancels out most of the undesirable TFA contributions to the overall stiffness estimates. A relevant extension to poly-crystal plasticity of this (RSL-TFA-HML) modeling is presented.     

A finite deformation theory of higher-order gradient crystal plasticity

For higher-order gradient crystal plasticity, a finite deformation formulation is presented. The theory does not deviate much from the conventional crystal plasticity theory. Only a back stress effect and additional differential equations for evolution of the geometrically necessary dislocation (GND) densities supplement the conventional theory within a non-work-conjugate framework in which there is no need to introduce higher-order microscopic stresses that would be work-conjugate to slip rate gradients. We discuss its connection to a work-conjugate type of finite deformation gradient crystal plasticity that is based on an assumption of the existence of higher-order stresses. Furthermore, a boundary-value problem for simple shear of a constrained thin strip is studied numerically, and some characteristic features of finite deformation are demonstrated through a comparison to a solution for the small deformation theory. As in a previous formulation for small deformation, the present formulation applies to the context of multiple and three-dimensional slip deformations.     

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 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.     

X-ray microdiffraction and strain gradient crystal plasticity studies of geometrically necessary dislocations near a Ni bicrystal grain boundary

We compare experimental measurements of inhomogeneous plastic deformation in a Ni bicrystal with crystal plasticity simulations. Polychromatic X-ray microdiffraction, orientation imaging microscopy and scanning electron microscopy, were used to characterize the geometrically necessary dislocation distribution of the bicrystal after uniaxial tensile deformation. Changes in the local crystallographic orientations within the sample reflect its plastic response during the tensile test. Elastic strain in both grains increases near the grain boundary. Finite element simulations were used to understand the influence of initial grain orientation and structural inhomogeneities on the geometrically necessary dislocations arrangement and distribution and to understand the underlying materials physics.     

The structure and energetics of, and the plasticity caused by, Eshelby dislocations

The structure of coaxial, or Eshelby, dislocations are computed using isotropic elasticity for arrays of up to 500 dislocations. The energies of these arrays are determined in order to predict the lowest energy configuration and multiple meta-stable configurations are often found. The energy from these elasticity predictions shows good agreement with molecular statics simulations of aluminum. From these simulations, the torque-twist curves are predicted and compared with molecular dynamics simulations.     

A theory of grain boundaries that accounts automatically for grain misorientation and grain-boundary orientation

This work is an attempt to answer the question:

Is there a physically natural method of characterizing the possible interactions between the slip systems of two grains that meet at a grain boundary—a method that could form the basis for the formulation of grain-boundary conditions?

Here we give a positive answer to this question based on the notion of a Burgers vector as described by a tensor field G on the grain boundary [Gurtin, M.E., Needleman, A., 2005. Boundary conditions in small-deformation single-crystal plasticity that account for the Burgers vector. J. Mech. Phys. Solids 53, 1-31]. We show that the magnitude of G can be expressed in terms of two types of moduli: inter-grain moduli that characterize slip-system interactions between the two grains; intra-grain moduli that for each grain characterize interactions between any two slip systems of that grain.We base the theory on microscopic force balances derived using the principle of virtual power, a version of the second law in the form of a free-energy imbalance, and thermodynamically compatible constitutive relations dependent on G and its rate. The resulting microscopic force balances represent flow rules for the grain boundary; and what is most important, these flow rules account automatically—via the intra- and inter-grain moduli—for the relative misorientation of the grains and the orientation of the grain boundary relative to those grains.     

Direct calculations of material parameters for gradient plasticity

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

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.     

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 theory for grain boundaries with strain-gradient plasticity

In this work, the effect of the material microstructural interface between two materials (i.e., grain boundary in polycrystalls) is adopted into a thermodynamic-based higher order strain gradient plasticity framework. The developed grain boundary flow rule accounts for the energy storage at the grain boundary due to the dislocation pile up as well as energy dissipation caused by the dislocation transfer through the grain boundary. The theory is developed based on the decomposition of the thermodynamic conjugate forces into energetic and dissipative counterparts which provides the constitutive equations to have both energetic and dissipative gradient length scales for the grain and grain boundary. The numerical solution for the proposed framework is also presented here within the finite element context. The material parameters of the gradient framework are also calibrated using an extensive set of micro-scale experimental measurements of thin metal films over a wide range of size and temperature of the samples.     

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.     

Micromechanical modeling of the elasto-viscoplastic behavior of semi-crystalline polymers

A micromechanically based constitutive model for the elasto-viscoplastic deformation and texture evolution of semi-crystalline polymers is developed. The model idealizes the microstructure to consist of an aggregate of two-phase layered composite inclusions. A new framework for the composite inclusion model is formulated to facilitate the use of finite deformation elasto-viscoplastic constitutive models for each constituent phase. The crystalline lamellae are modeled as anisotropic elastic with plastic flow occurring via crystallographic slip. The amorphous phase is modeled as isotropic elastic with plastic flow being a rate-dependent process with strain hardening resulting from molecular orientation. The volume-averaged deformation and stress within the inclusions are related to the macroscopic fields by a hybrid interaction model. The uniaxial compression of initially isotropic high density polyethylene (HDPE) is taken as a case study. The ability of the model to capture the elasto-plastic stress-strain behavior of HDPE during monotonic and cyclic loading, the evolution of anisotropy, and the effect of crystallinity on initial modulus, yield stress, post-yield behavior and unloading-reloading cycles are presented.     

A theory of strain-gradient plasticity for isotropic, plastically irrotational materials. Part I: Small deformations

This study develops a small-deformation theory of strain-gradient plasticity for isotropic materials in the absence of plastic rotation. The theory is based on 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 , the gradient of the plastic strain-rate, and

•

the free energy ψ to depend on the Burgers tensor .

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 the plastic strain. The microstresses are strictly dissipative when ψ is independent of the Burgers tensor, but when ψ depends on G the gradient microstress is partially energetic, and this, in turn, leads to a back stress and (hence) to Bauschinger-effects in the flow rule. It is further shown that dependencies of the microstresses on lead to strengthening and weakening effects in the flow rule.Typical macroscopic boundary conditions are supplemented by nonstandard microscopic boundary conditions associated with flow, and, as an aid to numerical solutions, a weak (virtual power) formulation of the nonlocal flow rule is derived.     

A one-dimensional theory of strain-gradient plasticity: Formulation, analysis, numerical results

This study develops a one-dimensional theory of strain-gradient plasticity based on: (i) a system of microstresses consistent with a microforce balance; (ii) a mechanical version of the second law that includes, via microstresses, work performed during viscoplastic flow; (iii) a constitutive theory that allows

•

the free-energy to depend on the gradient of the plastic strain, and

•

the microstresses to depend on the gradient of the plastic strain-rate.

The constitutive equations, whose rate-dependence is of power-law form, are endowed with energetic and dissipative gradient length-scales L and l, respectively, and allow for a gradient-dependent generalization of standard internal-variable hardening. The microforce balance when augmented by the constitutive relations for the microstresses results in a nonlocal flow rule in the form of a partial differential equation for the plastic strain. Typical macroscopic boundary conditions are supplemented by nonstandard microscopic boundary conditions associated with flow, and properties of the resulting boundary-value problem are studied both analytically and numerically. The resulting solutions are shown to exhibit three distinct physical phenomena:

(i)

standard (isotropic) internal-variable hardening;

(ii)

energetic hardening, with concomitant back stress, associated with plastic-strain gradients and resulting in boundary layer effects;

(iii)

dissipative strengthening associated with plastic strain-rate gradients and resulting in a size-dependent increase in yield strength.

     

Geometrically necessary dislocations in viscoplastic single crystals and bicrystals undergoing small deformations

In this study we develop a gradient theory of small-deformation single-crystal plasticity that accounts for geometrically necessary dislocations (GNDs). The resulting framework is used to discuss grain boundaries. The grains are allowed to slip along the interface, but growth phenomenona and phase transitions are neglected. The bulk theory is based on the introduction of a microforce balance for each slip system and includes a defect energy depending on a suitable measure of GNDs. The microforce balances are shown to be equivalent to nonlocal yield conditions for the individual slip systems, yield conditions that feature backstresses resulting from energy stored in dislocations. When applied to a grain boundary the theory leads to concomitant yield conditions: relative slip of the grains is activated when the shear stress reaches a suitable threshold; plastic slip in bulk at the grain boundary is activated only when the local density of GNDs reaches an assigned threshold. Consequently, in the initial stages of plastic deformation the grain boundary acts as a barrier to plastic slip, while in later stages the interface acts as a source or sink for dislocations. We obtain an exact solution for a simple problem in plane strain involving a semi-infinite compressed specimen that abuts a rigid material. We view this problem as an approximation to a situation involving a grain boundary between a grain with slip systems aligned for easy flow and a grain whose slip system alignment severely inhibits flow. The solution exhibits large slip gradients within a thin layer at the grain boundary.     

Size effects and dislocation patterning in two-dimensional bending

We perform atomistic Monte Carlo simulations of bending a Lennard-Jones single crystal in two dimensions. Dislocations nucleate only at the free surface as there are no sources in the interior of the sample. When dislocations reach sufficient density, they spontaneously coalesce to nucleate grain boundaries, and the resulting microstructure depends strongly on the initial crystal orientation of the sample. In initial yield, we find a reverse size effect, in which larger samples show a higher scaled bending moment than smaller samples for a given strain and strain rate. This effect is associated with source-limited plasticity and high strain rate relative to dislocation mobility, and the size effect in initial yield disappears when we scale the data to account for strain rate effects. Once dislocations coalesce to form grain boundaries, the size effect reverses and we find that smaller crystals support a higher scaled bending moment than larger crystals. This finding is in qualitative agreement with experimental results. Finally, we observe an instability at the compressed crystal surface that suggests a novel mechanism for the formation of a hillock structure. The hillock is formed when a high angle grain boundary, after absorbing additional dislocations, becomes unstable and folds to form a new crystal grain that protrudes from the free surface.