Free energy of dislocations in a multi-slip geometry

The collective dynamics of dislocations is the underlying mechanism of plastic deformation in metallic crystals. Dislocation motion in metals generally occurs on multiple slip systems. The simultaneous activation of different slip systems plays a crucial role in crystal plasticity models. In this contribution, we study the energetic interactions between dislocations on different slip systems by deriving the free energy in a multi-slip geometry. In this, we restrict ourselves to straight and parallel edge dislocations. The obtained free energy has a long-range mean-field contribution, a statistical contribution and a many-body contribution. The many-body contribution is a local function of the total dislocation density on each slip system, and can therefore not be written in terms of the net dislocation density only. Moreover, this function is a strongly non-linear and non-convex function of the density on different slip systems, and hence the coupling between slip systems is of great importance.     

Inspection of free energy functions in gradient crystal plasticity

The dislocation density tensor computed as the cud of plastic distortion is regarded as a new constitutive variable in crystal plasticity. The dependence of the free energy function on the dislocation density tensor is explored starting from a quadratic ansatz. Rank one and logarithmic dependencies are then envisaged based on considerations from the statistical theory of dislocations. The rele- vance of the presented free energy potentials is evaluated from the corresponding analytical solutions of the periodic two-phase laminate problem under shear where one layer is a single crystal material undergoing single slip and the second one remains purely elastic.     

Continuum dynamics of the formation,migration and dissociation of self-locked dislocation structures on parallel slip planes

In continuum models of dislocations, proper formulations of short-range elastic interactions of dislocations are crucial for capturing various types of dislocation patterns formed in crystalline materials. In this article, the continuum dynamics of straight dislocations distributed on two parallel slip planes is modelled through upscaling the underlying discrete dislocation dynamics. Two continuum velocity field quantities are introduced to facilitate the discrete-to-continuum transition. The first one is the local migration velocity of dislocation ensembles which is found fully independent of the short-range dislocation correlations. The second one is the decoupling velocity of dislocation pairs controlled by a threshold stress value, which is proposed to be the effective flow stress for single slip systems. Compared to the almost ubiquitously adopted Taylor relationship, the derived flow stress formula exhibits two features that are more consistent with the underlying discrete dislocation dynamics: (i) the flow stress increases with the in-plane component of the dislocation density only up to a certain value, hence the derived formula admits a minimum inter-dislocation distance within slip planes; (ii) the flow stress smoothly transits to zero when all dislocations become geometrically necessary dislocations. A regime under which inhomogeneities in dislocation density grow is identified, and is further validated through comparison with discrete dislocation dynamical simulation results. Based on the findings in this article and in our previous works, a general strategy for incorporating short-range dislocation correlations into continuum models of dislocations is proposed.     

Study of size effects in thin films by means of a crystal plasticity theory based on DiFT

In a recent publication, we derived the mesoscale continuum theory of plasticity for multiple-slip systems of parallel edge dislocations, motivated by the statistical-based nonlocal continuum crystal plasticity theory for single-glide given by Yefimov et al. [2004b. A comparison of a statistical-mechanics based plasticity model with discrete dislocation plasticity simulations. J. Mech. Phys. Solids 52, 279-300]. In this dislocation field theory (DiFT) the transport equations for both the total dislocation density and geometrically necessary dislocation (GND) density on each slip system were obtained from the Peach-Koehler interactions through both single and pair dislocation correlations. The effect of pair correlation interactions manifested itself in the form of a back stress in addition to the external shear and the self-consistent internal stress. We here present the study of size effects in single crystalline thin films with symmetric double slip using the novel continuum theory. Two boundary value problems are analyzed: (1) stress relaxation in thin films on substrates subject to thermal loading, and (2) simple shear in constrained films. In these problems, earlier discrete dislocation simulations had shown that size effects are born out of layers of dislocations developing near constrained interfaces. These boundary layers depend on slip orientations and applied loading but are insensitive to the film thickness. We investigate the stress response to changes in controlled parameters in both problems. Comparisons with previous discrete dislocation simulations are discussed.     

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.     

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.     

Simulation of edge dislocation emission from a crack under in-plane shear

Considered is the tandem emission of dislocations and dislocation dipoles from a crack under in-plane shear in one slip system as well as multiple slip systems. Effective stress intensity factors are determined by considering zones of local distortion similar to that in macro-plasticity. The dislocation free zone (DFZ) is also obtained which is analogous to the core region in fracture mechanics. Studied are effects of dislocation emission or development of plastic zone in front of the crack tip on the potential crack propagation based on the strain energy density factor criterion.     

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.     

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.     

Molecular dynamics study of solute strengthening in Al/Mg alloys

The strengthening of Al by Mg solute atoms is investigated using molecular dynamics (MD) studies of single dislocations moving through a field of randomly placed solutes. The MD method permits explicit treatment of “core” effects, dislocation pinning and deceleration, and dislocation unpinning by thermal activation, all under an applied load. Choice of an appropriate MD simulation cell size is assessed using analytic concepts developed by Labusch. The interaction energy of a single Mg atom with straight edge and screw dislocations is computed and compared with continuum models. Using the single Mg energies, a one-dimensional energy landscape for the motion of a straight edge dislocation through a random field of Mg solutes is computed. The minima in this landscape match well with those found in the MD simulations at zero temperature. The stress to unpin a straight edge dislocation trapped in a local energy minimum generated by the solutes is then predicted semi-analytically using the energy landscape, and good agreement is obtained with the MD results. At temperatures of 300 and 500 K, the thermally activated rate of unpinning vs. stress and temperature is calculated semi-analytically, and agreement with the full MD results is again obtained with the fitting of a single attempt frequency in a transition state model. The agreement of the semi-analytical models provides a basis for calculating yield stress vs. strain rate and temperature, resulting from statistical pinning, for the case of non-interacting dislocations on a single slip system, and for extending the analysis to study dynamic strain aging effects resulting from diffusion of Mg atoms around a pinned dislocation.     

A multi-scale computational model of crystal plasticity at submicron-to-nanometer scales

Plastic flow in crystal at submicron-to-nanometer scales involves many new interesting problems. In this paper, a unified computational model which directly combines 3D discrete dislocation dynamics (DDD) and continuum mechanics is developed to investigate the plastic behaviors at these scales. In this model, the discrete dislocation plasticity in a finite crystal is solved under a completed continuum mechanics framework: (1) an initial internal stress field is introduced to represent the preexisting stationary dislocations in the crystal; (2) the external boundary condition is handled by finite element method spontaneously; and (3) the constitutive relationship is based on the finite deformation theory of crystal plasticity, but the discrete plastic strains induced by the slip of the newly nucleated or propagating dislocations are calculated by dislocation dynamics methodology instead of phenomenological evolution equations used in conventional crystal plasticity. These discrete plastic strains are then localized to the continuum material points by a Burgers vector density function proposed by us. Various processes, such as loop dislocation evolution, dislocation junction formation etc., are simulated to verify the reliability of this computational model. Specifically, a uniaxial compression test for micro-pillars of Cu is simulated by this model to investigate the ‘dislocation starvation hardening’ observed in the recent experiment.     

Martensitic phase transition involving dislocations

A model of solid–solid phase transition involving dislocations in crystals is proposed within the nonlinear continuum dislocation theory (CDT). The co-existence of phases having piecewise constant plastic slip in laminates is possible for the two-well free energy density. The jumps of the plastic slip across the phase interfaces determine the surface dislocation densities at those incoherent boundaries. The number of phase interfaces should be determined by comparing the energy of dislocation arrays and the relaxed energy minimized among uniform plastic slips.     

Studying the effect of grain boundaries in dislocation density based crystal-plasticity finite element simulations

A dislocation density based constitutive model for the face centered cubic crystal structure has been implemented into a crystal-plasticity finite element framework and extended to consider the mechanical interaction between mobile dislocations and grain boundaries by the authors [Ma, A., Roters, F., Raabe, D., 2006a. A dislocation density based constitutive model for crystal-plasticity FEM including geometrically necessary dislocations. Acta Materialia 54, 2169–2179; Ma, A., Roters, F., Raabe, D., 2006b. On the consideration of interactions between dislocations and grain boundaries in crystal-plasticity finite element modeling – theory, experiments, and simulations. Acta Materialia 54, 2181–2194]. The approach to model the grain boundary resistance against slip is based on the introduction of an additional activation energy into the rate equation for mobile dislocations in the vicinity of internal interfaces. This energy barrier is derived from the assumption of thermally activated dislocation penetration events through grain boundaries. The model takes full account of the geometry of the grain boundaries and of the Schmid factors of the critically stressed incoming and outgoing slip systems. In this study we focus on the influence of the one remaining model parameter which can be used to scale the obstacle strength of the grain boundary.     

Continuum approximation of the Peach-Koehler force on dislocations in a slip plane

We derive a continuum model for the Peach-Koehler force on dislocations in a slip plane. To represent the dislocations, we use the disregistry across the slip plane, whose gradient gives the density and direction of the dislocations. The continuum model is derived rigorously from the Peach-Koehler force on dislocations in a region that contains many dislocations. The resulting continuum model can be written as the variation of an elastic energy that consists of the contribution from the long-range elastic interaction of dislocations and a correction due to the line tension effect.     

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

Strain hardening in 2D discrete dislocation dynamics simulations: A new ‘2.5D’ algorithm

The two-dimensional discrete dislocation dynamics (2D DD) method, consisting of parallel straight edge dislocations gliding on independent slip systems in a plane strain model of a crystal, is often used to study complicated boundary value problems in crystal plasticity. However, the absence of truly three dimensional mechanisms such as junction formation means that forest hardening cannot be modeled, unless additional so-called ‘2.5D’ constitutive rules are prescribed for short-range dislocation interactions. Here, results from three dimensional dislocation dynamics (3D DD) simulations in an FCC material are used to define new constitutive rules for short-range interactions and junction formation between dislocations on intersecting slip systems in 2D. The mutual strengthening effect of junctions on preexisting obstacles, such as precipitates or grain boundaries, is also accounted for in the model. The new ‘2.5D’ DD model, with no arbitrary adjustable parameters beyond those obtained from lower scale simulation methods, is shown to predict athermal hardening rates, differences in flow behavior for single and multiple slip, and latent hardening ratios. All these phenomena are well-established in the plasticity of crystals and quantitative results predicted by the model are in good agreement with experimental observations.