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.     

An elasto-plastic theory of dislocation and disclination fields

A linear theory of the elasto-plasticity of crystalline solids based on a continuous representation of crystal defects – dislocations and disclinations – is presented. The model accounts for the translational and rotational aspects of lattice incompatibility, respectively associated with the presence of dislocations and disclinations. The defects content relates to the incompatible plastic strain and curvature tensors. The stress state is described by using the conjugate variables to strain and curvature, i.e., the stress and couple-stress tensors. Defect motion is described by two transport equations. A dynamic interplay between dislocations and disclinations results from a disclination-induced source term in the transport of dislocations. Thermodynamic guidance provides the driving forces conjugate to dislocation and disclination velocity in a continuous context, as well as admissible constitutive relations for the latter. When dislocation and disclination velocity vanish, the model reduces to deWit’s elasto-static theory of crystal defects. It also reduces to Acharya’s linear elasto-plastic theory for dislocation fields when the disclination density is ignored. The theory is intended for use in instances where rotational defects matter, such as grain boundaries. To illustrate its applicability, a finite high-angle tilt boundary is modeled using a disclination dipole and its behavior under tensile loading normal to the boundary is shown.     

Continuum representation of systems of dislocation lines: A general method for deriving closed-form evolution equations

Plasticity is governed by the evolution of, in general anisotropic, systems of dislocations. We seek to faithfully represent this evolution in terms of density-like variables which average over the discrete dislocation microstructure. Starting from T. Hochrainer's continuum theory of dislocations (CDD) (Hochrainer, 2015), we introduce a methodology based on the ‘Maximum Information Entropy Principle’ (MIEP) for deriving closed-form evolution equations for dislocation density measures of different order. These equations provide an optimum representation of the kinematic properties of systems of curved and connected dislocation lines with the information contained in a given set of density measures. The performance of the derived equations is benchmarked against other models proposed in the literature, using discrete dislocation dynamics simulations as a reference. As a benchmark problem we study dislocations moving in a highly heterogeneous, persistent-slip-band like geometry. We demonstrate that excellent agreement with discrete simulations can be obtained in terms of a very small number of averaged dislocation fields containing information about the edge and screw components of the total and excess (geometrically necessary) dislocation densities. From these the full dislocation orientation distribution which emerges as dislocations move through a channel-wall structure can be faithfully reconstructed.     

Dislocation Patterns and Work-Hardening in Crystalline Plasticity

We propose a continuum model for the evolution of the total dislocation densities in fcc crystals, in the framework of rate-independent plasticity. The basic physical features which are taken into account are: (i) the role of dislocations in hardening; (ii) the relations between the slip velocity and dislocation mobility; (iii) the energetics of self and mutual interactions between dislocations; (iv) nonlocal effects in the interaction between dislocations. A set of reaction–diffusion equations is obtained, with mobilities depending on the slip velocities, which is able to describe the formation of dislocation walls and cells. To this effect, the results of numerical simulations in two special cases are presented. Mathematics Subject Classifications (2000) 74C15, 74C20.     

Viscoplasticity of metals: Comments on statistical approaches to dislocation reactions

Dislocations, i.e., line defects in the crystal structure, are at the origin of viscoplastic deformation in metals. In the literature, the statistical approaches to the dislocation dynamics in terms of dislocation densities are often used. There, the key ingredient for modeling the ubiquitous strain hardening is the interplay between the mobile and forest (i.e., immobile) dislocations, which is captured in terms of reaction-type contributions in the evolution of the dislocation densities. In this paper, we demonstrate that a certain class of such models is in conflict with thermodynamic principles. The origin of this conflict is identified. Specifically, the absence of the reversal processes for any of the reactions is problematic and sharply contrasts to usual chemical reactions. Possible solutions for restoring the thermodynamic admissibility are discussed.     

An analysis of dislocation nucleation near a free surface

Molecular dynamics analyses of defect-free aluminum single crystals subject to bending are carried out to investigate dislocation nucleation from free surfaces. A principal aim of the analyses is to provide background for the development of dislocation nucleation criteria for use in discrete dislocation plasticity calculations. The molecular dynamics simulations use an embedded atom potential for aluminum. Bending is imposed on a strip by specifying a linear variation of displacement rate on opposite edges. The overall bending response is determined and the character of the dislocations nucleated is identified. It is found that the stress magnitudes at the instant of dislocation nucleation are nearly an order of magnitude smaller than for homogeneous bulk dislocation nucleation. The characterization of dislocation nucleation in terms of various phenomenological nucleation criteria is explored, in particular: (i) a critical resolved shear stress; (ii) the onset of an elastic instability; and (iii) a critical stress-gradient criterion. It is found that dislocation nucleation is not well-represented by a critical value of the resolved shear stress but is reasonably well-represented by the critical stress-gradient criterion.     

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.     

A dynamic flow rule for viscoplasticity in polycrystalline solids under high strain rates

For visco-plasticity in polycrystalline solids under high strain rates, we introduce a dynamic flow rule (also called the micro-force balance) that has a second order time derivative term in the form of micro-inertia. It is revealed that this term, whose physical origin is traced to dynamically evolving dislocations, has a profound effect on the macro-continuum plastic response. Based on energy equivalence between the micro-part of the kinetic energy and that associated with the fictive dislocation mass in the continuous dislocation distribution (CDD) theory, an explicit expression for the micro-inertial length scale is derived. The micro-force balance together with the classical momentum balance equations thus describes the viscoplastic response of the isotropic polycrystalline material. Using rational thermodynamics, we arrive at constitutive equations relating the thermodynamic forces (stresses) and fluxes. A consistent derivation of temperature evolution is also provided, thus replacing the empirical route. The micro-force balance, supplemented with the constitutive relations for the stresses, yields a locally hyperbolic flow rule owing to the micro-inertia term. The implication of micro-inertia on the continuum response is explicitly demonstrated by reproducing experimentally observed stress–strain responses under high strain-rate loadings and varying temperatures. An interesting finding is the identification of micro-inertia as the source of oscillations in the stress–strain response under high strain rates.     

Couple stresses in crystalline solids: origins from plastic slip gradients, dislocation core distortions, and three-body interatomic potentials

In the presence of plastic slip gradients, compatibility requires gradients in elastic rotation and stretch tensors. In a crystal lattice the gradient in elastic rotation can be related to bond angle changes at cores of so-called geometrically necessary dislocations. The corresponding continuum strain energy density can be obtained from an interatomic potential that includes two- and three-body terms. The three-body terms induce restoring moments that lead to a couple stress tensor in the continuum limit. The resulting stress and couple stress jointly satisfy a balance law. Boundary conditions are obtained upon stress, couple stress, strain and strain gradient tensors. This higher-order continuum theory was formulated by Toupin (Arch. Ration. Mech. Anal. 11 (1962) 385). Toupin's theory has been extended in this work to incorporate constitutive relations for the stress and couple stress under multiplicative elastoplasticity. The higher-order continuum theory is exploited to solve a boundary value problem of relevance to single crystal and polycrystalline nano-devices. It is demonstrated that certain slip-dominated deformation mechanisms increase the compliance of nanostructures in bending-dominated situations. The significance of these ideas in the context of continuum plasticity models is also dwelt upon.     

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.     

Line-integral representations for the elastic displacements,stresses and interaction energy of arbitrary dislocation loops in transversely isotropic bimaterials

The elastic displacements, stresses and interaction energy of arbitrarily shaped dislocation loops with general Burgers vectors in transversely isotropic bimaterials (i.e. joined half-spaces) are expressed in terms of simple line integrals for the first time. These expressions are very similar to their isotropic full-space counterparts in the literature and can be easily incorporated into three-dimensional (3D) dislocation dynamics (DD) simulations for hexagonal crystals with interfaces/surfaces. All possible degenerate cases, e.g. isotropic bimaterials and isotropic half-space, are considered in detail. The singularities intrinsic to the classical continuum theory of dislocations are removed by spreading the Burgers vector anisotropically around every point on the dislocation line according to three particular spreading functions. This non-singular treatment guarantees the equivalence among different versions of the energy formulae and their consistency with the stress formula presented in this paper. Several numerical examples are provided as verification of the derived dislocation solutions, which further show significant influence of material anisotropy and bimaterial interface on the elastic fields and interaction energy of dislocation loops.     

Modeling of thermally activated dislocation glide and plastic flow through local obstacles

A unified phenomenological model is developed to study the dislocation glide through weak obstacles during the first stage of plastic deformation in metals. This model takes into account both the dynamical responses of dislocations during the flight process and thermal activations while dislocations are bound by obstacle arrays. The average thermal activation rate is estimated using an analytical model based on the generalized Friedel relations. Then, the average flight velocity after an activation event is obtained numerically by discrete dislocation dynamics (DD). To simulate the dynamical dislocation behavior, the inertia term is implemented into the equation of dislocation motion within the DD code. The results from the DD simulations, coupled with the analytical model, determine the total dislocation velocity as a function of the stress and temperatures. By choosing parameters typical of the face centered cubic metals, the model reproduces both obstacle control and drag control motion in low and high velocity regimes, respectively. As expected by other string models, dislocation overshoots of obstacles caused by the dislocation inertia at the collisions are enhanced as temperature goes down.     

Continuum framework for dislocation structure,energy and dynamics of dislocation arrays and low angle grain boundaries

We present a continuum framework for dislocation structure, energy and dynamics of dislocation arrays and low angle grain boundaries that are allowed to be nonplanar or nonequilibrium. In our continuum framework, we define a dislocation density potential function on the dislocation array surface or grain boundary to describe the orientation dependent continuous distribution of dislocations in a very simple and accurate way. The continuum formulations incorporate both the long-range dislocation interaction and the local dislocation line energy, and are derived from the discrete dislocation model. The continuum framework recovers the classical Read–Shockley energy formula when the long-range elastic fields of the low angle grain boundaries are canceled out. Applications of our continuum framework in this paper are focused on dislocation structures on static planar and nonplanar low angle grain boundaries and misfitting interfaces. We present two methods under our continuum framework for this purpose, including the method based on the Frank׳s formula and the energy minimization method. We show that for any (planar or nonplanar) low angle grain boundary, the Frank׳s formula holds if and only if the long-range stress field in the continuum model is canceled out, and it does not necessarily hold for a total energy minimum dislocation structure.     

Shocks and slip systems: Predictions from a mesoscale theory of continuum dislocation dynamics

Exploring a recently developed mesoscale continuum theory of dislocation dynamics, we derive three predictions about plasticity and grain boundary formation in crystals. (1) There is a residual stress jump across grain boundaries and plasticity-induced cell walls as they form, which self-consistently acts to attract neighboring dislocations; residual stress in this theory appears as a remnant of the driving force behind wall formation under both polygonization and plastic deformation. We derive the predicted asymptotic late-time dynamics of the grain-boundary formation process. (2) During grain boundary formation at high temperatures, there is a predicted cusp in the elastic energy density. (3) In early stages of plasticity, when only one type of dislocation is active (single-slip), cell walls do not form in the theory; instead we predict the formation of a hitherto unrecognized jump singularity in the dislocation density.     

Multiscale plasticity modeling: coupled atomistics and discrete dislocation mechanics

A computational method (CADD) is presented whereby a continuum region containing dislocation defects is coupled to a fully atomistic region. The model is related to previous hybrid models in which continuum finite elements are coupled to a fully atomistic region, with two key advantages: the ability to accomodate discrete dislocations in the continuum region and an algorithm for automatically detecting dislocations as they move from the atomistic region to the continuum region and then correctly “converting” the atomistic dislocations into discrete dislocations, or vice-versa. The resulting CADD model allows for the study of 2d problems involving large numbers of defects where the system size is too big for fully atomistic simulation, and improves upon existing discrete dislocation techniques by preserving accurate atomistic details of dislocation nucleation and other atomic scale phenomena. Applications to nanoindentation, atomic scale void growth under tensile stress, and fracture are used to validate and demonstrate the capabilities of the model.     

A note on material forces in finite inelasticity

In this contribution we aim to elaborate material forces in the context of multiplicative elasto-plasticity, which is considered as a representative and general framework for finite inelasticity. The comparison of different representations of the balance of linear momentum enables us to identify relevant Eshelbian stress tensors and corresponding volume forces. These material, or rather configurational, forces incorporate dislocation density tensors due to the general incompatibility of the underlying intermediate configuration. As an interesting application, the celebrated Peach–Koehler force, driving single dislocations in the context of finite-deformation inelasticity, allows representation in terms of the derived configurational volume forces.