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.     

Crystal plasticity model with enhanced hardening by geometrically necessary dislocation accumulation

A strain gradient dependent crystal plasticity approach is used to model the constitutive behaviour of polycrystal FCC metals under large plastic deformation. Material points are considered as aggregates of grains, subdivided into several fictitious grain fractions: a single crystal volume element stands for the grain interior whereas grain boundaries are represented by bi-crystal volume elements, each having the crystallographic lattice orientations of its adjacent crystals. A relaxed Taylor-like interaction law is used for the transition from the local to the global scale. It is relaxed with respect to the bi-crystals, providing compatibility and stress equilibrium at their internal interface. During loading, the bi-crystal boundaries deform dissimilar to the associated grain interior. Arising from this heterogeneity, a geometrically necessary dislocation (GND) density can be computed, which is required to restore compatibility of the crystallographic lattice. This effect provides a physically based method to account for the additional hardening as introduced by the GNDs, the magnitude of which is related to the grain size. Hence, a scale-dependent response is obtained, for which the numerical simulations predict a mechanical behaviour corresponding to the Hall-Petch effect. Compared to a full-scale finite element model reported in the literature, the present polycrystalline crystal plasticity model is of equal quality yet much more efficient from a computational point of view for simulating uniaxial tension experiments with various grain sizes.     

On large-strain finite element solutions of higher-order gradient crystal plasticity

A finite-strain higher-order gradient crystal plasticity model accounting for the backstress effect originating from the existence of geometrically necessary dislocations (GNDs) is applied to plane strain finite element analysis. Different element types are tested to seek out an element formulation that is reliable and useful for solving problems involving severe plastic deformation. In the present finite element formulation, the GND density rates are chosen to be additional nodal degrees of freedom. Different orders of shape functions are employed for the interpolation of displacement rates and GND density rates. Their effects on solutions are examined in detail by considering three boundary value problems: a simple shear of a constrained layer (a film), a compression problem with loading surfaces impenetrable to dislocations, and a tension problem involving shear band formation. In all the cases, the formulation in which eight-node elements with reduced integration and four-node elements with full integration are used respectively for displacement rates and the GND density rates gives reasonable solutions. In addition to the discussion on the choice of finite elements, detailed behavior in gradient-dependent solids, such as the accumulation of GND density and the distribution of backstress on each slip system, is investigated by utilizing the reliable computational results obtained.     

Continuum thermodynamic models for crystal plasticity including the effects of geometrically-necessary dislocations

The purpose of this work is the formulation of constitutive models for the inelastic material behaviour of single crystals and polycrystals in which geometrically necessary dislocations (GNDs) may develop and influence this behaviour. To this end, we focus on the dependence of the development of such dislocations on the inhomogeneity of the inelastic deformation in the material. More precisely, in the crystal plasticity context, this is a relation between the density of GNDs and the inhomogeneity of inelastic deformation in glide systems. In this work, two models for GND density and its evolution, i.e., a glide-system-based model, and a continuum model, are formulated and investigated. As it turns out, the former of these is consistent with the original two-dimensional GND model of Ashby (Philos. Mag. 21 (1970) 399), and the latter with the more recent model of Dai and Parks (Proceedings of Plasticity ’97, Neat Press, 1997, p. 17). Since both models involve a dependence of the inelastic state of a material point on the (history of the) inhomogeneity of the glide-system inelastic deformation, their incorporation into crystal plasticity modelling necessarily implies a corresponding non-local generalization of this modelling. As it turns out, a natural quantity on which to base such a non-local continuum thermodynamic generalization, i.e., in the context of crystal plasticity, is the glide-system (scalar) slip deformation. In particular, this is accomplished here by treating each such slip deformation as either (1), a generalized “gradient” internal variable, or (2), as a scalar internal degree-of-freedom. Both of these approaches yield a corresponding generalized Ginzburg-Landau- or Cahn-Allen-type field relation for this scalar deformation determined in part by the dependence of the free energy on the dislocation state in the material. In the last part of the work, attention is focused on specific models for the free energy and its dependence on this state. After summarizing and briefly discussing the initial-boundary-value problem resulting from the current approach as well as its algorithmic form suitable for numerical implementation, the work ends with a discussion of additional aspects of the formulation, and in particular the connection of the approach to GND modelling taken here with other approaches.     

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.     

Plastic deformation modelling of tempered martensite steel block structure by a nonlocal crystal plasticity model

The plastic deformations of tempered martensite steel representative volume elements with different martensite block structures have been investigated by using a nonlocal crystal plasticity model which considers isotropic and kinematic hardening produced by plastic strain gradients. It was found that pronounced strain gradients occur in the grain boundary region even under homogeneous loading. The isotropic hardening of strain gradients strongly influences the global stress–strain diagram while the kinematic hardening of strain gradients influences the local deformation behaviour. It is found that the additional strain gradient hardening is not only dependent on the block width but also on the misorientations or the deformation incompatibilities in adjacent blocks.     

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.     

A semi-implicit integration scheme for rate independent finite crystal plasticity

An efficient new numerical integration scheme is presented for rate independent crystal plasticity theory. A key feature of this approach is the ability to identify active slip systems prior to determining their shearing rate. Options are described for various cases of slip system hardening, including self hardening and latent hardening. Alternatives for the constitutive update are explored, including hyperelasticity based on the multiplicative decomposition of the deformation gradient as well as application of the consistency condition in a much more efficient hypoelastic formulation. Several conclusions are drawn concerning the influences of elastic and plastic properties on the activation of slip systems and their subsequent shearing rates. Key among these is the fact that, once activated, shearing rates are independent of the levels of shear flow resistance on the slip systems, provided that the plastic hardening moduli are much less in magnitude than the elastic moduli, as is usually the case. Determination of active slip systems and their shearing rates depend on the degree of elastic anisotropy of the crystal, but not on the magnitude of elastic stiffness.     

Mechanism-based strain gradient crystal plasticity—I. Theory

We have been developing the theory of mechanism-based strain gradient plasticity (MSG) to model size-dependent plastic deformation at micron and submicron length scales. The core idea has been to incorporate the concept of geometrically necessary dislocations into the continuum plastic constitutive laws via the Taylor hardening relation. Here we extend this effort to develop a mechanism-based strain gradient theory of crystal plasticity. In this theory, an effective density of geometrically necessary dislocations for a specific slip plane is introduced via a continuum analog of the Peach-Koehler force in dislocation theory and is incorporated into the plastic constitutive laws via the Taylor relation.     

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.     

考虑运动硬化的DD3镍基合金力学行为研究

详细介绍了镍基合金的晶体塑性本构模型,在Asaro大变形晶体塑性框架下,详细介绍了镍基合金的晶体塑性本构模型,在Asaro大变形晶体塑性框架下,引入了运动硬化规律,考虑了温度和应变率对晶体塑性变形的影响,通过针对每个滑移系考虑屈服准则和流动规律建立了晶体塑性模型. 对积分过程进行了推导,通过编写ABAQUS材料用户子程序(UMAT), 实现本构模型的有限元积分算法. 在此基础上模拟了DD3镍基单晶合金在单轴拉伸和循环载荷下的响应,并与实验数据进行了对比. 利用该模型可以很好地模拟镍基单晶所具有的各向异性特性,体现了镍基单晶在循环载荷作用下的拉-压不对称性.      

Grain boundary interface mechanics in strain gradient crystal plasticity

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

Dual-phase steel sheets under cyclic tension–compression to large strains: Experiments and crystal plasticity modeling

In this work, we develop a physically-based crystal plasticity model for the prediction of cyclic tension–compression deformation of multi-phase materials, specifically dual-phase (DP) steels. The model is elasto–plastic in nature and integrates a hardening law based on statistically stored dislocation density, localized hardening due to geometrically necessary dislocations (GNDs), slip-system-level kinematic backstresses, and annihilation of dislocations. The model further features a two level homogenization scheme where the first level is the overall response of a two-phase polycrystalline aggregate and the second level is the homogenized response of the martensite polycrystalline regions. The model is applied to simulate a cyclic tension–compression–tension deformation behavior of DP590 steel sheets. From experiments, we observe that the material exhibits a typical decreasing hardening rate during forward loading, followed by a linear and then a non-linear unloading upon the load reversal, the Bauschinger effect, and changes in hardening rate during strain reversals. To predict these effects, we identify the model parameters using a portion of the measured data and validate and verify them using the remaining data. The developed model is capable of predicting all the particular features of the cyclic deformation of DP590 steel, with great accuracy. From the predictions, we infer and discuss the effects of GNDs, the backstresses, dislocation annihilation, and the two-level homogenization scheme on capturing the cyclic deformation behavior of the material.     

Size effects and idealized dislocation microstructure at small scales: Predictions of a Phenomenological model of Mesoscopic Field Dislocation Mechanics: Part I

A Phenomenological Mesoscopic Field Dislocation Mechanics (PMFDM) model is developed, extending continuum plasticity theory for studying initial-boundary value problems of small-scale plasticity. PMFDM results from an elementary space-time averaging of the equations of Field Dislocation Mechanics (FDM), followed by a closure assumption from any strain-gradient plasticity model that attempts to account for effects of geometrically necessary dislocations (GNDs) only in work hardening. The specific lower-order gradient plasticity model chosen to substantiate this work requires one additional material parameter compared to its conventional continuum plasticity counterpart. The further addition of dislocation mechanics requires no additional material parameters. The model (a) retains the constitutive dependence of the free-energy only on elastic strain as in conventional continuum plasticity with no explicit dependence on dislocation density, (b) does not require higher-order stresses, and (c) does not require a constitutive specification of a ‘back-stress’ in the expression for average dislocation velocity/plastic strain rate. However, long-range stress effects of average dislocation distributions are predicted by the model in a mechanistically rigorous sense. Plausible boundary conditions (with obvious implication for corresponding interface conditions) are discussed in some detail from a physical point of view. Energetic and dissipative aspects of the model are also discussed. The developed framework is a continuous-time model of averaged dislocation plasticity, without having to rely on the notion of incremental work functions, their convexity properties, or their minimization. The tangent modulus relating stress rate and total strain rate in the model is the positive-definite tensor of linear elasticity, and this is not an impediment to the development of idealized microstructure in the theory and computations, even when such a convexity property is preserved in a computational scheme. A model of finite deformation, mesoscopic single crystal plasticity is also presented, motivated by the above considerations.Lower-order gradient plasticity appears as a constitutive limit of PMFDM, and the development suggests a plausible boundary condition on the plastic strain rate for this limit that is appropriate for the modeling of constrained plastic flow in three-dimensional situations.     

Strain gradient crystal plasticity effects on flow localization

In metal grains one of the most important failure mechanisms involves shear band localization. As the band width is small, the deformations are affected by material length scales. To study localization in single grains a rate-dependent crystal plasticity formulation for finite strains is presented for metals described by the reformulated Fleck–Hutchinson strain gradient plasticity theory. The theory is implemented numerically within a finite element framework using slip rate increments and displacement increments as state variables. The formulation reduces to the classical crystal plasticity theory in the absence of strain gradients. The model is used to study the effect of an internal material length scale on the localization of plastic flow in shear bands in a single crystal under plane strain tension. It is shown that the mesh sensitivity is removed when using the nonlocal material model considered. Furthermore, it is illustrated how different hardening functions affect the formation of shear bands.     

Strain gradient crystal plasticity analysis of a single crystal containing a cylindrical void

The effects of void size and hardening in a hexagonal close-packed single crystal containing a cylindrical void loaded by a far-field equibiaxial tensile stress under plane strain conditions are studied. The crystal has three in-plane slip systems oriented at the angle 60° with respect to one another. Finite element simulations are performed using a strain gradient crystal plasticity formulation with an intrinsic length scale parameter in a non-local strain gradient constitutive framework. For a vanishing length scale parameter the non-local formulation reduces to a local crystal plasticity formulation. The stress and deformation fields obtained with a local non-hardening constitutive formulation are compared to those obtained from a local hardening formulation and to those from a non-local formulation. Compared to the case of the non-hardening local constitutive formulation, it is shown that a local theory with hardening has only minor effects on the deformation field around the void, whereas a significant difference is obtained with the non-local constitutive relation. Finally, it is shown that the applied stress state required to activate plastic deformation at the void is up to three times higher for smaller void sizes than for larger void sizes in the non-local material.     

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.     

Bauschinger and size effects in thin-film plasticity due to defect-energy of geometrical necessary dislocations

The Bauschinger and size effects in the thinfilm plasticity theory arising from the defect-energy of geometrically necessary dislocations (GNDs) are analytically investigated in this paper. Firstly, this defect-energy is deduced based on the elastic interactions of coupling dislocations (or pile-ups) moving on the closed neighboring slip plane. This energy is a quadratic function of the GNDs density, and includes an elastic interaction coefficient and an energetic length scale L. By incorporating it into the work- conjugate strain gradient plasticity theory of Gurtin, an energetic stress associated with this defect energy is obtained, which just plays the role of back stress in the kinematic hardening model. Then this back-stress hardening model is used to investigate the Bauschinger and size effects in the tension problem of single crystal Al films with passivation layers. The tension stress in the film shows a reverse dependence on the film thickness h. By comparing it with discrete-dislocation simulation results, the length scale L is determined, which is just several slip plane spacing, and accords well with our physical interpretation for the defect- energy. The Bauschinger effect after unloading is analyzed by combining this back-stress hardening model with a friction model. The effects of film thickness and pre-strain on the reversed plastic strain after unloading are quantified and qualitatively compared with experiment results.     

Creep,plasticity, and fatigue of single crystal superalloy

Single crystal components in gas turbine engines are subject to such extreme temperatures and stresses that life prediction becomes highly inaccurate resulting in components that can only be shown to meet their requirements through experience. Reliable life prediction methodologies are required both for design and life management. In order to address this issue we have developed a thermo-viscoplastic constitutive model for single crystal materials. Our incremental large strain formulation additively decomposes the inelastic strain rate into components along the octahedral and cubic slip planes. We have developed a crystallographic-based creep constitutive model able to predict sigmoidal creep behavior of Ni base superalloys. Inelastic shear rate along each slip system is expressed as a sum of a time dependent creep component and a rate independent plastic component. We develop a new robust, computationally efficient rate-independent crystal plasticity approach and combined it with creep flow rule calibrated for Ni-based superalloys. The transient variation of each of the inelastic components includes a back stress for kinematic hardening and latent hardening parameters to account for the stress evolution with inelastic strain as well as the evolution for dislocation densities. The complete formulation accurately predicts both monotonic and cyclic tests at different crystallographic orientations for constant and variable temperature conditions (low cycle fatigue (LCF) and thermo-mechanical fatigue (TMF) tests). Based on the test and modeling results we formulate a new life prediction criterion suitable for both LCF and TMF conditions.     

A modified model for simulating latent hardening during the plastic deformation of rate-dependent FCC polycrystals

A strain hardening model for the plastic deformation of rate-dependent FCC crystals is proposed based on experimental observations previously reported for single crystals. This model, which is an extension of that employed by et al. [1983], includes both the self-hardening and latent hardening of the slip systems. The differential hardening of the latent systems is assumed to arise from the interaction between glide dislocations and forests. With this hardening model and a rate-sensitive crystal plasticity theory, the deformation behavior of FCC polycrystals can be predicted from the deformation response of the constituent single crystals. As examples, the uniaxial tensile behaviour of pure aluminum and copper polycrystals is simulated using the extended model, and the results are compared with published experimental data. The effects of latent hardening on polycrystal deformation, especially on flow stress and the formation of tensile textures, are discussed.