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.     

Characterization of yielding behavior of polycrystalline metals with single crystal plasticity based on representative characteristic length

Single crystal plasticity based on a representative characteristic length is proposed and introduced into a homogenization approach based on finite element analyses, which are applied to characterization of distinctive yielding behaviors of polycrystalline metals, yield-point elongation, and grain size strengthening. The computational manner for an implicit stress update is derived with the framework of a standard multi-surface plasticity at finite strain, where the evolution of the characteristic lengths are numerically converted from the accumulated slips of all of slip systems by exploiting the mathematical feature of the characteristic length as the intermediate function of the plastic internal variables. Furthermore, a constitutive model for a single crystal reproduces the stress–strain curve divided into three parts. Using two-scale finite element analysis, the macroscopic stress–strain response with yield-point elongation under a situation of low dislocation density is reproduced. Finally, the grain size effect on the yield strength is analyzed with modeling of the grain boundary in the context of the proposed constitutive model and is discussed from both macroscopic and microscopic views.     

Crystal plasticity modeling of slip activity in Ti–6Al–4V under high cycle fatigue loading

Deformation micromechanisms of a Ti–6Al–4V alloy under fatigue loading at room temperature are studied using a three-dimensional crystal plasticity constitutive model. The model employs a minimum set of fitting parameters based on experimental data for Ti–6Al–4V. Single slip is strongly favored through a softening law that affects mainly the driving force for slip on the first activated slip system. Cyclic deformation behavior at the macroscopic scale and at the local scale of grains is analyzed through the simulation of 20 cycles of fatigue on a polycrystalline structure of 900 randomly oriented grains. The progressive activation of slip (basal, prismatic, and pyramidal) is analyzed and compared to experimental observations.     

A three-dimensional crystal plasticity model for duplex Ti–6Al–4V

A rate dependent crystal plasticity model for the α/β Ti–Al alloy Ti–6Al–4V with duplex microstructure is developed and presented herein. Duplex Ti–6Al–4V is a dual-phase alloy consisting of an hcp structured matrix primary α-phase and secondary lamellar α + β domains that are composed of alternating layers of secondary α laths and bcc structured residual β laths. The model accounts for distinct three-dimensional slip geometry for each phase, anisotropic and length scale dependent slip system strengths, the non-planar dislocation core structure of prismatic screw dislocations in the primary α-phase, and crystallographic texture. The model is implemented in the general purpose finite element code (ABAQUS, 2005. Ver 6.5, Hibbitt, Karlsson, and Sorensen, Inc., Pawtucket, RI) via a UMAT subroutine.     

Modelling the behaviour of polycrystalline austenitic steel with twinning-induced plasticity effect

A micromechanical model using the scale transition method in elastoviscoplasticity has been developed to describe the behaviour of those austenitic steels that display a TWIP effect. A physically based constitutive equation at the grain scale is proposed considering two inelastic strain modes: crystallographic slip and twinning. The typical organizations of microtwins observed in electron microscopy are considered, and the twin–slip as well as the twin–twin interactions are accounted for. The parameters for slip are first fitted on the uniaxial tensile response obtained at intermediate temperatures (when twinning is inhibited). Then, the parameters associated with twinning are identified using the stress–strain curve at room temperature. The simulated results in both macro and micro scales are in good agreement with experimentally obtained results.     

A multiscale approach for modeling scale-dependent yield stress in polycrystalline metals

Modeling of scale-dependent characteristics of mechanical properties of metal polycrystals is studied using both discrete dislocation dynamics and continuum crystal plasticity. The initial movements of dislocation arc emitted from a Frank-Read type dislocation source and bounded by surrounding grain boundaries are examined by dislocation dynamics analyses system and we find the minimum resolved shear stress for the FR source to emit at least one closed loop. When the grain size is large enough compared to the size of FR source, the minimum resolved shear stress levels off to a certain value, but when the grain size is close to the size of the FR source, the minimum resolved shear stress shows a sharp increase. These results are modeled into the expression of the critical resolved shear stress of slip systems and continuum mechanics based crystal plasticity analyses of six-grained polycrystal models are made. Results of the crystal plasticity analyses show a distinct increase of macro- and microscopic yield stress for specimens with smaller mean grain diameter. Scale-dependent characteristics of the yield stress and its relation to some control parameters are discussed.     

A finite element formulation based on non-associated plasticity for sheet metal forming

In the present paper, a finite element formulation based on non-associated plasticity is developed. In the constitutive formulation, isotropic hardening is assumed and an evolution equation for the hardening parameter consistent with the principle of plastic work equivalence is introduced. The yield function and plastic potential function are considered as two different functions with functional form as the yield function of Hill [Hill, R., 1948. Theory of yielding and plastic flow of anisotropic metals. Proc. Roy. Soc. A 193, 281–297] or Karafillis–Boyce associated model [Karafillis, A.P. Boyce, M., 1993. A general anisotropic yield criterion using bounds and a transformation weighting tensor. J. Mech. Phys. Solids 41, 1859–1886]. Algorithmic formulations of constitutive models that utilize associated or non-associated flow rule coupled with Hill or Karafillis–Boyce stress functions are derived by application of implicit return mapping procedure. Capabilities in predicting planar anisotropy of the Hill and Karafillis–Boyce stress functions are investigated considering material data of Al2008-T4 and Al2090-T3 sheet samples. The accuracy of the derived stress integration procedures is investigated by calculating iso-error maps.     

Geometrically necessary dislocation structure organization in FCC bicrystal subjected to cyclic plasticity

Crystal plasticity finite element analysis of cyclic deformation of compatible type FCC bicrystals are performed. The model specimen used in the analysis is a virtual FCC bicrystal with an isotropic elastic property; therefore, the effect of constraint due to elastic incompatibility does not appear. The results of the analysis show the strain-amplitude-dependence of both the organization of the GND structure and the stress–strain behavior. The calculated stress–strain curve with the largest strain amplitude shows additional cyclic hardening. The microscopic mechanisms of the strain-amplitude-dependent organization of the GND structure and additional cyclic hardening behavior are discussed in terms of the activation of secondary slip system(s). Finally, the effects of the elastic anisotropy, the lattice friction stress and the interaction between dislocations are also argued.     

Discrete dislocation plasticity analysis of the grain size dependence of the flow strength of polycrystals

The grain size dependence of the flow strength of polycrystals is analyzed using plane strain, discrete dislocation plasticity. Dislocations are modeled as line singularities in a linear elastic solid and plasticity occurs through the collective motion of large numbers of dislocations. Constitutive rules are used to model lattice resistance to dislocation motion, as well as dislocation nucleation, dislocation annihilation and the interaction with obstacles. The materials analyzed consist of micron scale grains having either one or three slip systems and two types of grain arrangements: either a checker-board pattern or randomly dispersed with a specified volume fraction. Calculations are carried out for materials with either a high density of dislocation sources or a low density of dislocation sources. In all cases, the grain boundaries are taken to be impenetrable to dislocations. A Hall–Petch type relation is predicted with Hall–Petch exponents ranging from ≈0.3 to ≈1.6 depending on the number of slip systems, the grain arrangement, the dislocation source density and the range of grain sizes to which a Hall–Petch expression is fit. The grain size dependence of the flow strength is obtained even when no slip incompatibility exists between grains suggesting that slip blocking/transmission governs the Hall–Petch effect in the simulations.     

Debonding failure and size effects in micro-reinforced composites

Failure in micro-reinforced composites is investigated numerically using the strain-gradient plasticity theory of Gudmundson [Gudmundson, P., 2004. A unified treatment of strain gradient plasticity. Journal of the Mechanics and Physics of Solids 52 (6) 1379–1406] in a plane strain visco-plastic formulation. Bi-axially loaded unit cells are used and failure is modeled using a cohesive zone at the reinforcement interface. During debonding a sudden stress drop in the overall average stress–strain response is observed. Adaptive higher-order boundary conditions are imposed at the reinforcement interface for realistically modeling the restrictions on moving dislocations as debonding occurs. It is found that the influence of the imposed higher-order boundary conditions at the interface is minor. If strain-gradient effects are accounted for a void with a smooth shape develops at the reinforcement interface while a smaller void having a sharp tip nucleates if strain-gradient effects are excluded. Using orthogonalization of the plastic strain gradient with three corresponding material length scales it is found that, the first length scale dominates the evaluated overall average stress–strain response, the second one only has a small effect and the third one has an intermediate effect. Finally, studies of reinforcement having elliptical cross-sections show rather significant gradients of stress which is not seen for the corresponding circular cross-sections. Also, an increased drop in the overall load carrying capacity is observed for cross-sections elongated perpendicular to the principal tensile direction compared to the corresponding circular cross-sections.     

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

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

Effects of spatial grain orientation distribution and initial surface topography on sheet metal necking

The finite element method is used to numerically simulate localized necking in AA6111-T4 under stretching. The measured EBSD data (grain orientations and their spatial distributions) are directly incorporated into the finite element model and the constitutive response at an integration point is described by the single crystal plasticity theory. We assume that localized necking is associated with surface instability, the onset of unstable growth in surface roughening. It is demonstrated that such a surface instability/necking is the natural outcome of the present approach, and the artificial initial imperfection necessitated by the macroscopic M–K approach [Marciniak and Kuczynski (1967). Int. J. Mech. Sci. 9, 609–620] is not relevant in the present analysis. The effects of spatial orientation distribution, material strain rate sensitivity, texture evolution, and initial surface topography on necking are discussed. It is found that localized necking depends strongly on both the initial texture and its spatial orientation distribution. It is also demonstrated that the initial surface topography has only a small influence on necking.     

Path-dependence of the forming limit stresses in a sheet metal

The effect of changing strain paths on the forming limit stresses of sheet metals is investigated using the Marciniak–Kuczyński model and a phenomenological plasticity model with non-normality effects [Kuroda, M., Tvergaard, V., 2001. A phenomenological plasticity model with non-normality effects representing observations in crystal plasticity. J. Mech. Phys. Solids 49, 1239–1263]. Forming limits are simulated for linear stress paths and two types of combined loading: a combined loading consisting of two linear stress paths in which unloading is included between the first and second loadings (combined loading A), and combined loading in which the strain path is abruptly changed without unloading (combined loading B). The forming limit stresses calculated for combined loading A agree well with those calculated for the linear stress paths, while the forming limit curves in strain space depend strongly on the strain paths. The forming limit stresses calculated for the combined loading B do not, however, coincide with those calculated for the linear stress paths. The strain-path dependence of the forming limit stress is discussed in detail by observing the strain localization process.     

Distribution based model for the grain boundaries in polycrystalline plasticity

This contribution focuses on the development of constitutive models for the grain boundary region between two crystals, relying on the dislocation based polycrystalline model documented in (Evers, L.P., Parks, D.M., Brekelmans, W.A.M., Geers, M.G.D., 2002. Crystal plasticity model with enhanced hardening by geometrically necessary dislocation accumulation. J. Mech. Phys. Solids 50, 2403–2424; Evers, L.P., Brekelmans, W.A.M., Geers, M.G.D., 2004a. Non-local crystal plasticity model with intrinsic SSD and GND effects. J. Mech. Phys. Solids 52, 2379–2401; Evers, L.P., Brekelmans, W.A.M., Geers, M.G.D., 2004b. Scale dependent crystal plasticity framework with dislocation density and grain boundary effects. Int. J. Solids Struct. 41, 5209–5230). The grain boundary is first viewed as a geometrical surface endowed with its own fields, which are treated here as distributions from a mathematical point of view. Regular and singular dislocation tensors are introduced, defining the grain equilibrium, both in the grain core and at the boundary of both grains. Balance equations for the grain core and grain boundary are derived, that involve the dislocation density distribution tensor, in both its regular and singular contributions. The driving force for the motion of the geometrically necessary dislocations is identified from the pull-back to the lattice configuration of the quasi-static balance of momentum, that reveals the duality between the stress and the curl of the elastic gradient. Criteria that govern the flow of mobile geometrically necessary dislocations (GNDs) through the grain boundary are next elaborated on these bases. Specifically, the sign of the projection of a lattice microtraction on the glide velocity defines a necessary condition for the transmission of incoming GNDs, thereby rendering the set of active slip systems for the glide of outgoing dislocations. Viewing the grain boundary as adjacent bands in each grain with a constant GND density in each, the driving force for the grain boundary slip is further expressed in terms of the GND densities and the differently oriented slip systems in each grain. A semi-analytical solution is developed in the case of symmetrical slip in a bicrystal under plane strain conditions. It is shown that the transmission of plastic slip occurs when the angle made by the slip direction relative to the grain boundary normal is less than a critical value, depending on the ratio of the GND densities and the orientation of the transmitted dislocations.     

Yielding and intrinsic plasticity of Ti–Zr–Ni–Cu–Be bulk metallic glass

Bulk metallic glass with composition Ti₄₀Zr₂₅Ni₈Cu₉Be₁₈ exhibits considerably high compressive yield stress, significant plasticity (with a concomitant vein-like fracture morphology) and relatively low density. Yielding and intrinsic plasticity of this alloy are discussed in terms of its thermal and elastic properties. An influence of normal stresses acting on the shear plane is evidenced by: (i) the fracture angle (<45°) and (ii) finite-element simulations of nanoindentation curves, which require the use of a specific yield criterion, sensitive to local normal stresses acting on the shear plane, to properly match the experimental data. The ratio between hardness and compressive yield strength (constraint factor) is analyzed in terms of several models and is best adjusted using a modified expanding cavity model incorporating a pressure-sensitivity index defined by the Drucker–Prager yield criterion. Furthermore, comparative results from compression tests and nanoindentation reveal that deformation also causes strain softening, a phenomenon which is accompanied with the occurrence of serrated plastic flow and results in a so-called indentation size effect (ISE). A new approach to model the ISE of this metallic glass using the free volume concept is presented.     

Effects of texture on shear band formation in plane strain tension/compression and bending

In this study, effects of typical texture components observed in rolled aluminum alloy sheets on shear band formation in plane strain tension/compression and bending are systematically studied. The material response is described by a generalized Taylor-type polycrystal model, in which each grain is characterized in terms of an elastic–viscoplastic continuum slip constitutive relation. First, a simple model analysis in which the shear band is assumed to occur in a weaker thin slice of material is performed. From this simple model analysis, two important quantities regarding shear band formation are obtained: i.e. the critical strain at the onset of shear banding and the corresponding orientation of shear band. Second, the shear band development in plane strain tension/compression is analyzed by the finite element method. Predictability of the finite element analysis is compared to that of the simple model analysis. Third, shear band developments in plane strain pure bending of a sheet specimen with the typical textures are studied. Regions near the surfaces in a bent sheet specimen are approximately subjected to plane strain tension or compression. From this viewpoint, the bendability of a sheet specimen may be evaluated, using the knowledge regarding shear band formation in plane strain tension/compression. To confirm this and to encompass overall deformation of a bent sheet specimen, including shear bands, finite element analyses of plane strain pure bending are carried out, and the predicted shear band formation in bent specimens is compared to that in the tension/compression problem. Finally, the present results are compared to previous related studies, and the efficiency of the present method for materials design in future is discussed.     

Towards a homogenized plasticity theory which predicts structural and microstructural size effects

We study an idealized bending problem where two types of size effects are present – one induced by the non-uniform (macro) deformation, the other due to the (internal) resistance at grain boundaries. Classical models are not able to capture either of the two types of size dependent behavior. A remedy is to adopt a gradient crystal plasticity formulation which allows one to study the direct influence of different microstructural properties on the material response. However, it is computationally expensive to do so for a typical engineering problem since the discretization has to be done at a sub-granular level. In this paper, a homogenization theory is proposed such that the small deformation gradient crystal plasticity framework by Cermelli and Gurtin [Cermelli, P., Gurtin, M.E., 2002. Geometrically necessary dislocations in viscoplastic single crystals and bicrystals undergoing small deformations. Int. J. Solids Struct. 39, 6281–6309] translates from the micro to macro level consistently. Microstructural properties thus propagate naturally to the macro scale and the homogenized solutions compare well with the fine scale analyses for the two limit cases – microhard and microfree conditions. Three length scale parameters, i.e. the intrinsic length scale, grain size and the foil thickness, manifest themselves in the homogenized solution, thus capturing both types of size effects. We further discuss on the interplay and competition between the two size effects.     

Modeling aspects of low plastic strain amplitude multiaxial cyclic plasticity in conventional and ultrafine grain nickel

Aspects of the cyclically saturated responses of initially annealed, conventional grain size (average grain diameter of approximately 50 μm) and electrodeposited, ultrafine grain (grains from 20 to 500 nm) nickel to reversed proportional and 90° out-of-phase axial-torsional, strain-controlled cycling at a nominally constant equivalent inelastic strain amplitude of approximately 100 μ strain are reported. An anisotropic, axial-torsional subspace version of Abdel-Karim and Ohno’s kinematic hardening model is presented. Within the framework of conventional small strain, rate-independent plasticity, this approach is used to model the responses. An anhysteretic, phenomenolically based, magnetomechanical model is coupled to the rate-independent plasticity model to include the cyclic magnetostriction response. The kinematic hardening parameter determination scheme, using the proportional path responses, is described. The model correlations achieved are presented and the ability of the resulting models to capture the 90° out-of-phase responses is examined. The model parameter sets, as determined from the proportional responses, require small changes to result in more accurate correlation of the 90° out-of-phase responses and the implications of this are discussed. The relative values of the model parameters between the two materials reflect the initial microstructures. Persistent mean stresses associated with mean total strains imposed are successfully modeled for the proportional strain path responses but not for the 90° out-of-phase responses.