Non-convex rate dependent strain gradient crystal plasticity and deformation patterning

A rate dependent strain gradient crystal plasticity framework is presented where the displacement and the plastic slip fields are considered as primary variables. These coupled fields are determined on a global level by solving simultaneously the linear momentum balance and the slip evolution equation, which is derived in a thermodynamically consistent manner. The formulation is based on the 1D theory presented in Yalcinkaya et al. (2011), where the patterning of plastic slip is obtained in a system with non-convex energetic hardening through a phenomenological double-well plastic potential. In the current multi-dimensional multi-slip analysis the non-convexity enters the framework through a latent hardening potential presented in Ortiz and Repettto (1999) where the microstructure evolution is obtained explicitly via a lamination procedure. The current study aims the implicit evolution of deformation patterns due to the incorporated physically based non-convex potential.     

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.     

Micromechanical modelling of the effect of plastic deformation on the mechanical behaviour in pseudoelastic shape memory alloys

Except for the recoverable strain induced by phase transformation, NiTi alloys are very ductile even in the martensite phase. The purpose of the present paper is to study the influence of permanent deformation, which results from plastic deformation of martensite, on the mechanical behaviour of pseudoelastic NiTi alloys. Based on phenomenological theory of martensitic transformation and crystal plasticity, a new three dimensional micromechanical model is proposed by coupling both the slip and twinning deformation mechanisms. The present model is implemented as User MATerial subroutine (UMAT) into ABAQUS/Standard to study the influences of plastic deformation on the stress and strain fields, and on the evolution of martensite transformation. Results show that with the increasing of plastic deformation the residual strain increases and the phase transformation stress–strain curves from the martensite to austenite become steeper and less obvious. Both characteristics, stabilisation of martensite and impedance of the reverse transformation, due to plastic deformation are captured.     

A Dedicated DIC Methodology for Characterizing Plastic Deformation in Single Crystals

This paper focuses on the development of an appropriate Digital Image Correlation (DIC) methodology based on Image Registration and dedicated for characterizing the plastic deformation in single crystals. A pure nickel single crystal specimen is plastically deformed in tension and investigated by DIC technique. Based on the measured kinematic fields, the proposed method enables to identify the slip activity on the crystal surface and to locate precisely the slip band interfaces at microscale which behave as kinematic discontinuities. The computed displacement data are projected on a well-defined physical basis containing slip details, then the strain fields can be derived directly from a set of analytical functions. The possible errors in displacement induced by this projection approach are evaluated. Finally, some results of the evaluated strain fields are presented. It demonstrates that the developed DIC methodology allows quantitative characterization of a heterogeneous deformation process and promotes further relationships to be established between slip activity and strain field evolution in single crystals.     

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.     

On dislocation patterning: Multiple slip effects in the rate equation approach

Strain localization and dislocation pattern formation are typical features of plastic deformation in metals and alloys. Glide and climb dislocation motion, along with accompanying production/annihilation processes, lead to the occurrence of instabilities of initially uniform dislocation distributions. These instabilities result to the development of various types of dislocation microstructures (dislocation cells, slip and kink bands, persistent slip bands, labyrinth structures, etc.), depending on the externally applied loading and the intrinsic lattice constraints. The term “dislocation patterning” was introduced over 20 years ago by the third author and a corresponding “gradient dislocation dynamics” framework was suggested to describe such phenomena. In the W–A model proposed at that time by the last two authors, it was shown how coupled nonlinear evolution equations of the reaction-diffusion type for the forest (immobile) and gliding (mobile) dislocation densities can generate dislocation microstructures which correspond to walls perpendicular to the slip direction for Cu-crystals oriented for single slip under cyclic loading conditions. This model is adapted to the multiple slip case here. Weakly nonlinear analysis predicts that dislocation patterns should correspond to domains of walls perpendicular to each slip direction and separated by domain walls in the same orientations. This result is confirmed by numerical analysis and experimental observations. The present model generalizes the original W–A model to the case of multiple slip and considers also explicitly gradient effects by allowing for non-uniform dislocation velocities and internal stress effects.     

A dislocation density based crystal plasticity finite element model: Application to a two-phase polycrystalline HCP/BCC composites

We present a multiscale model for anisotropic, elasto-plastic, rate- and temperature-sensitive deformation of polycrystalline aggregates to large plastic strains. The model accounts for a dislocation-based hardening law for multiple slip modes and links a single-crystal to a polycrystalline response using a crystal plasticity finite element based homogenization. It is capable of predicting local stress and strain fields based on evolving microstructure including the explicit evolution of dislocation density and crystallographic grain reorientation. We apply the model to simulate monotonic mechanical response of a hexagonal close-packed metal, zirconium (Zr), and a body-centered cubic metal, niobium (Nb), and study the texture evolution and deformation mechanisms in a two-phase Zr/Nb layered composite under severe plastic deformation. The model predicts well the texture in both co-deforming phases to very large plastic strains. In addition, it offers insights into the active slip systems underlying texture evolution, indicating that the observed textures develop by a combination of prismatic, pyramidal, and anomalous basal slip in Zr and primarily {110}〈111〉 slip and secondly {112}〈111〉 slip in Nb.     

A numerical study of crack tip constraint in ductile single crystals

In this work, the effect of crack tip constraint on near-tip stress and deformation fields in a ductile FCC single crystal is studied under mode I, plane strain conditions. To this end, modified boundary layer simulations within crystal plasticity framework are performed, neglecting elastic anisotropy. The first and second terms of the isotropic elastic crack tip field, which are governed by the stress intensity factor K and T-stress, are prescribed as remote boundary conditions and solutions pertaining to different levels of T-stress are generated. It is found that the near-tip deformation field, especially, the development of kink or slip shear bands, is sensitive to the constraint level. The stress distribution and the size and shape of the plastic zone near the crack tip are also strongly influenced by the level of T-stress, with progressive loss of crack tip constraint occurring as T-stress becomes more negative. A family of near-tip fields is obtained which are characterized by two terms (such as K and T or J and a constraint parameter Q) as in isotropic plastic solids.     

Wedge indentation into elastic-plastic single crystals, 1: Asymptotic fields for nearly-flat wedge

Asymptotic stress and deformation fields under the contact point singularities of a nearly-flat wedge indenter and of a flat punch are derived for elastic ideally-plastic single crystals with three effective in-plane slip systems that admit a plane strain deformation state. Face-centered cubic (FCC), body-centered cubic (BCC), and hexagonal-close packed (HCP) crystals are considered. The asymptotic fields for the flat punch are analogous to those at the tip of a stationary crack, so a potential solution is that the deformation field consists entirely of angular constant stress plastic sectors separated by rays of plastic deformation across which stresses change discontinuously. The asymptotic fields for a nearly-flat wedge indenter are analogous to those of a quasistatically growing crack tip fields in that stress discontinuities can not exist across sector boundaries. Hence, the asymptotic fields under the contact point singularities of a nearly-flat wedge indenter are significantly different than those under a flat punch. A family of solutions is derived that consists entirely of elastically deforming angular sectors separated by rays of plastic deformation across which the stress state is continuous. Such a solution can be found for FCC and BCC crystals, but it is shown that the asymptotic fields for HCP crystals must include at least one angular constant stress plastic sector. The structure of such fields is important because they play a significant role in the establishment of the overall fields under a wedge indenter in a single crystal. Numerical simulations—discussed in detail in a companion paper—of the stress and deformation fields under the contact point singularity of a wedge indenter for a FCC crystal possess the salient features of the analytical solution.     

Interaction between phase transformations and dislocations at the nanoscale. Part 1. General phase field approach

Thermodynamically consistent, three-dimensional (3D) phase field approach (PFA) for coupled multivariant martensitic transformations (PTs), including cyclic PTs, variant–variant transformations (i.e., twinning), and dislocation evolution is developed at large strains. One of our key points is in the justification of the multiplicative decomposition of the deformation gradient into elastic, transformational, and plastic parts. The plastic part includes four mechanisms: dislocation motion in martensite along slip systems of martensite and slip systems of austenite inherited during PT and dislocation motion in austenite along slip systems of austenite and slip systems of martensite inherited during reverse PT. The plastic part of the velocity gradient for all these mechanisms is defined in the crystal lattice of the austenite utilizing just slip systems of austenite and inherited slip systems of martensite, and just two corresponding types of order parameters. The explicit expressions for the Helmholtz free energy and the transformation and plastic deformation gradients are presented to satisfy the formulated conditions related to homogeneous thermodynamic equilibrium states of crystal lattice and their instabilities. In particular, they result in a constant (i.e., stress- and temperature-independent) transformation deformation gradient and Burgers vectors. Thermodynamic treatment resulted in the determination of the driving forces for change of the order parameters for PTs and dislocations. It also determined the boundary conditions for the order parameters that include a variation of the surface energy during PT and exit of dislocations. Ginzburg–Landau equations for dislocations include variation of properties during PTs, which in turn produces additional contributions from dislocations to the Ginzburg–Landau equations for PTs. A complete system of coupled PFA and mechanics equations is presented. A similar theory can be developed for PFA to dislocations and other PTs, like reconstructive PTs and diffusive PTs described by the Cahn–Hilliard equation, as well as twinning and grain boundaries evolution.     

High strain-rate modeling of the interfacial effects of dispersed particles in high strength aluminum alloys

The interfacial effects of dispersed particles on the dynamic deformation of high strength aluminum alloys have been investigated using an eigenstrain-based formulation coupled with dislocation-density based crystalline plasticity and a microstructurally based finite element framework. This accounts for the unrelaxed plastic strains associated with the interfacial behavior of dispersed particles, such as Orowan looping. Particle spacing had a significant effect on the distribution of plastic shear slip, with localization occurring between the particles for smaller particle spacing. The eigenstress field associated with larger particles led to longer-range interaction of pressure fields, which can promote void coalescence for nucleated voids at the particle-matrix interface. Grain orientation also had a significant effect on the behavior associated with the particles, with plastic shear slip localizing at the particle-matrix interfaces for low angle grain-boundary (GB) misorientations, and at GBs and GB junctions for high angle GB misorientations.     

Self-consistent modeling of large plastic deformation, texture and morphology evolution in semi-crystalline polymers

A self-consistent model for semi-crystalline polymers is proposed to study their constitutive behavior, texture and morphology evolution during large plastic deformation. The material is considered as an aggregate of composite inclusions, each representing a stack of crystalline lamellae with their adjacent amorphous layers. The deformation within the inclusions is volume-averaged over the phases. The interlamellar shear is modeled as an additional slip system with a slip direction depending on the inclusion's stress. Hardening of the amorphous phase due to molecular orientation and, eventually, coarse slip, is introduced via Arruda-Boyce hardening law for the corresponding plastic resistance. The morphology evolution is accounted for through the change of shape of the inclusions under the applied deformation gradient. The overall behavior is obtained via a viscoplastic tangent self-consistent scheme. The model is applied to high density polyethylene (HDPE). The stress-strain response, texture and morphology changes are simulated under different modes of straining and compared to experimental data as well as to the predictions of other models.     

Dislocation pile-ups in bicrystals within continuum dislocation theory

Within continuum dislocation theory the plastic deformation of bicrystals under a mixed deformation of plane constrained uniaxial extension and shear is investigated with regard to the nucleation of dislocations and the dislocation pile-up near the phase boundaries of a model bicrystal with one active slip system within each single crystal. For plane uniaxial extension, we present a closed-form analytical solution for the evolution of the plastic distortion and of the dislocation network in the case of symmetric slip planes (i.e. for twins), which exhibits an energetic as well as a dissipative threshold for the dislocation nucleation. The general solution for non-symmetric slip systems is obtained numerically. For a combined deformation of extension and shear, we analyze the possibility of linearly superposing results obtained for both loading cases independently. All solutions presented in this paper also display the Bauschinger effect of translational work hardening and a size effect typical to problems of crystal plasticity.     

Cylindrical void in a rigid-ideally plastic single crystal III: Hexagonal close-packed crystal

The analytical solution is derived for the plane strain stress field around a cylindrical void in a hexagonal close-packed single crystal with three in-plane slip systems oriented at the angle π/3 with respect to one another. The critical resolved shear stress on each slip system is assumed to be equal. The crystal is loaded by both internal pressure and a far-field equibiaxial compressive stress. The deformation field takes the form of angular sectors, called slip sectors, within which only one slip system is active; the boundaries between different sectors are radial lines. The stress fields are derived by enforcing equilibrium and a rigid, ideally plastic constitutive relationship, in the spirit of anisotropic slip line theory. The results show that each slip sector is divided into smaller regions denoted as stress sectors and the stress state valid within each stress sector is derived. It is shown that stresses are unique and are continuous within stress sectors and across stress sector boundaries, but the gradient of stresses is not continuous across the boundaries between stress sectors. The solution shows self-similarity in that the stresses over the entire domain can be determined from the stresses within a small region adjacent to the void by invoking certain scaling and symmetry properties. In addition, the stress state exhibits periodicity along logarithmic spirals which emanate from the void. The results predict that the mean value of in-plane pressure required to activate plastic deformation around a void in a single crystal can be higher than that necessary for a void in an isotropic material and is sensitive to the orientation of the slip systems relative to the void.     

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.     

On the strain hardening and texture evolution in high manganese steels: Experiments and numerical investigation

We present a systematic investigation on the strain hardening and texture evolution in high manganese steels where twinning induced plasticity (TWIP) plays a significant role for the materials' plastic deformation. Motivated by the stress–strain behavior of typical TWIP steels with compositions of Fe, Mn, and C, we develop a mechanistic model to explain the strain-hardening in crystals where deformation twinning dominates the plastic deformation. The classical single crystal plasticity model accounting for both dislocation slip and deformation twinning are then employed to simulate the plastic deformation in polycrystalline TWIP steels. While only deformation twinning is activated for plasticity, the simulations with samples composed of voronoi grains cannot fully capture the texture evolution of the TWIP steel. By including both twinning deformation and dislocation slip, the model is able to capture both the stress–strain behaviors and the texture evolution in Fe–Mn–C TWIP steel in different boundary-value problems. Further analysis on the strain contributions by both mechanisms suggests that deformation twinning plays the dominant role at the initial stage of plasticity in TWIP steels, and dislocation slip becomes increasingly important at large strains.     

A mesoscale investigation of strain rate effect on dynamic deformation of single-crystal copper

A combined finite element (FE) simulation and discrete dislocation dynamics (DD) approach has been developed in this paper to investigate the dynamic deformation of single-crystal copper at mesoscale. The DD code yields the plastic strain based on the slip of dislocations and serves as a substitute for the 3D constitutive form used in the usual FE computation, which is implemented into ABAQUS/Standard with a user-defined material subroutine. On the other hand, the FE code computes the displacement and stress field during the dynamic deformation. The loading rate effects on the yield stress and the deformation patterning of single-crystal copper are investigated. With the increasing of strain rate, the yield stress of single-crystal copper increases rapidly. A critical strain rate exists in each single-crystal copper block for the given size and dislocation sources, below which the yield stress is relatively insensitive to the strain rate. The dislocation patterning changes from non-uniform to uniform under high-strain-rate. The shear stresses in the bands are higher than that in the neighboring regions, which are formed shear bands in the crystal. The band width increases with the strain rate, which often take places where the damage occurs.     

Incorporation of twinning into a crystal plasticity finite element model: Evolution of lattice strains and texture in Zircaloy-2

A crystal plasticity finite element code is developed to model lattice strains and texture evolution of HCP crystals. The code is implemented to model elastic and plastic deformation considering slip and twinning based plastic deformation. The model accounts for twinning reorientation and growth. Twinning, as well as slip, is considered to follow a rate dependent formulation. The results of the simulations are compared to previously published in situ neutron diffraction data. Experimental results of the evolution of the texture and lattice strains under uniaxial tension/compression loading along the rolling, transverse, and normal direction of a piece of rolled Zircaloy-2 are compared with model predictions. The rate dependent formulation introduced is capable of correctly capturing the influence of slip and twinning deformation on lattice strains as well as texture evolution.     

Thermodynamically consistent phase field approach to dislocation evolution at small and large strains

A thermodynamically consistent, large strain phase field approach to dislocation nucleation and evolution at the nanoscale is developed. Each dislocation is defined by an order parameter, which determines the magnitude of the Burgers vector for the given slip planes and directions. The kinematics is based on the multiplicative decomposition of the deformation gradient into elastic and plastic contributions. The relationship between the rates of the plastic deformation gradient and the order parameters is consistent with phenomenological crystal plasticity. Thermodynamic and stability conditions for homogeneous states are formulated and satisfied by the proper choice of the Helmholtz free energy and the order parameter dependence on the Burgers vector. They allow us to reproduce desired lattice instability conditions and a stress-order parameter curve, as well as to obtain a stress-independent equilibrium Burgers vector and to avoid artificial dissipation during elastic deformation. The Ginzburg–Landau equations are obtained as the linear kinetic relations between the rate of change of the order parameters and the conjugate thermodynamic driving forces. A crystalline energy coefficient for dislocations is defined as a periodic step-wise function of the coordinate along the normal to the slip plane, which provides an energy barrier normal to the slip plane and determines the desired, mesh-independent height of the dislocation bands for any slip system orientation. Gradient energy contains an additional term, which excludes the localization of a dislocation within a height smaller than the prescribed height, but it does not produce artificial interface energy. An additional energy term is introduced that penalizes the interaction of different dislocations at the same point. Non-periodic boundary conditions for dislocations are introduced which include the change of the surface energy due to the exit of dislocations from the crystal. Obtained kinematics, thermodynamics, and kinetics of dislocations at large strains are simplified for small strains and rotations, as well.     

Numerical investigation of the interaction between the martensitic transformation front and the plastic strain in austenite

Phase-field simulations of the martensitic transformation (MT) in an austenitic matrix which has already undergone the plastic deformation are carried out. For this purpose the elasto-plastic phase-field approach of incoherent MT developed in a previous work [Kundin et al., 2011. A phase-field model for incoherent martensitic transformations including plastic accommodation processes in the austenite. J. Mech. Phys. Solids 59, 2082–2012] is used. The evolution equation for the dislocation density field is extended by taking into account the thermal and athermal annihilation of the dislocations in the austenitic matrix and the athermal annihilation at the transformation front. It is shown that the plastic deformation in the austenite caused by the MT interacts with the dislocation field and the MT front that leads to an inhomogeneous increasing of the total dislocation density. During the phase transformation one part of the dislocations in the austenite is inherited by the martensitic phase and this inheritance depends on the kinetics and the crystallography of MT. Another part of dislocations annihilates at the transformation front and decreases the dislocation density in the growing martensite. Based on the simulation results the specific type of phenomenological dependency between the inherited dislocations, the martensite phase fraction and the plastic deformation is proposed.