A comparison of nonlocal continuum and discrete dislocation plasticity predictions

Discrete dislocation simulations of two boundary value problems are used as numerical experiments to explore the extent to which the nonlocal crystal plasticity theory of Gurtin (J. Mech. Phys. Solids 50 (2002) 5) can reproduce their predictions. In one problem simple shear of a constrained strip is analyzed, while the other problem concerns a two-dimensional model composite with elastic reinforcements in a crystalline matrix subject to macroscopic shear. In the constrained layer problem, boundary layers develop that give rise to size effects. In the composite problem, the discrete dislocation solutions exhibit composite hardening that depends on the reinforcement morphology, a size dependence of the overall stress-strain response for some morphologies, and a strong Bauschinger effect on unloading. In neither problem are the qualitative features of the discrete dislocation results represented by conventional continuum crystal plasticity. The nonlocal plasticity calculations here reproduce the behavior seen in the discrete dislocation simulations in remarkable detail.     

Dislocation-based micropolar single crystal plasticity: Comparison of multi- and single criterion theories

Two new formulations of micropolar single crystal plasticity are presented within a geometrically linear setting. The construction of yield criteria and flow rules for generalized continuum theories with higher-order stresses can be done in one of two ways: (i) a single criterion can be introduced in terms of a combined equivalent stress and inelastic rate or (ii) or individual criteria can be specified for each conjugate stress/inelastic kinematic rate pair, a so-called multi-criterion theory. Both single and multi-criterion theories are developed and discussed within the context of dislocation-based constitutive models. Parallels and distinctions are made between the proposed theories and some of the alternative generalized crystal plasticity models that can be found in the literature. Parametric numerical simulations of a constrained thin film subjected to simple shear are conducted via finite element analysis using a simplified 2-D version of the fully 3-D theory to highlight the influence of specific model components on the resulting deformation under both loading and unloading conditions. The deformation behavior is quantified in terms of the average stress-strain response and the local shear strain and geometrically necessary dislocation density distributions. It is demonstrated that micropolar single crystal plasticity can qualitatively capture the same range of behaviors as slip gradient-based models, while offering a simpler numerical implementation and without introducing plastic slip rates as generalized traction-conjugate velocities subject to an additional microforce balance.     

Non-local crystal plasticity model with intrinsic SSD and GND effects

A strain gradient-dependent crystal plasticity approach is presented to model the constitutive behaviour of polycrystal FCC metals under large plastic deformation. In order to be capable of predicting scale dependence, the heterogeneous deformation-induced evolution and distribution of geometrically necessary dislocations (GNDs) are incorporated into the phenomenological continuum theory of crystal plasticity. Consequently, the resulting boundary value problem accommodates, in addition to the ordinary stress equilibrium condition, a condition which sets the additional nodal degrees of freedom, the edge and screw GND densities, proportional (in a weak sense) to the gradients of crystalline slip. Next to this direct coupling between microstructural dislocation evolutions and macroscopic gradients of plastic slip, another characteristic of the presented crystal plasticity model is the incorporation of the GND-effect, which leads to an essentially different constitutive behaviour than the statistically stored dislocation (SSD) densities. The GNDs, by their geometrical nature of locally similar signs, are expected to influence the plastic flow through a non-local back-stress measure, counteracting the resolved shear stress on the slip systems in the undeformed situation and providing a kinematic hardening contribution. Furthermore, the interactions between both SSD and GND densities are subject to the formation of slip system obstacle densities and accompanying hardening, accountable for slip resistance. As an example problem and without loss of generality, the model is applied to predict the formation of boundary layers and the accompanying size effect of a constrained strip under simple shear deformation, for symmetric double-slip conditions.     

Plasticity size effects in voided crystals

The shear and equi-biaxial straining responses of periodic voided single crystals are analysed using discrete dislocation plasticity and a continuum strain gradient crystal plasticity theory. In the discrete dislocation formulation, the dislocations are all of edge character and are modelled as line singularities in an elastic material. The lattice resistance to dislocation motion, dislocation nucleation, dislocation interaction with obstacles and annihilation are incorporated through a set of constitutive rules. Over the range of length scales investigated, both the discrete dislocation and strain gradient plasticity formulations predict a negligible size effect under shear loading. By contrast, under equi-biaxial loading both plasticity formulations predict a strong size dependence with the flow strength approximately scaling inversely with the void spacing. Excellent agreement is obtained between predictions of the two formulations for all crystal types and void volume fractions considered when the material length scale in the non-local plasticity model is chosen to be (about 10 times the slip plane spacing in the discrete dislocation models).     

Plane constrained shear of single crystal strip with two active slip systems

Within continuum dislocation theory the plane constrained shear of a single crystal strip with two active slip systems is considered. An analytical solution is found for symmetric double slip which exhibits the energetic and dissipative thresholds for dislocation nucleation, the Bauschinger translational work hardening, and the size effects. Comparison with discrete dislocation simulations shows good agreement between the discrete and continuum approaches. Numerical procedures in the general case of non-symmetric double slip are proposed.     

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.     

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.     

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.     

Multiscale modeling of the plasticity in an aluminum single crystal

This paper describes a numerical, hierarchical multiscale modeling methodology involving two distinct bridges over three different length scales that predicts the work hardening of face centered cubic crystals in the absence of physical experiments. This methodology builds a clear bridging approach connecting nano-, micro- and meso-scales. In this methodology, molecular dynamics simulations (nanoscale) are performed to generate mobilities for dislocations. A discrete dislocations numerical tool (microscale) then uses the mobility data obtained from the molecular dynamics simulations to determine the work hardening. The second bridge occurs as the material parameters in a slip system hardening law employed in crystal plasticity models (mesoscale) are determined by the dislocation dynamics simulation results. The material parameters are computed using a correlation procedure based on both the functional form of the hardening law and the internal elastic stress/plastic shear strain fields computed from discrete dislocations. This multiscale bridging methodology was validated by using a crystal plasticity model to predict the mechanical response of an aluminum single crystal deformed under uniaxial compressive loading along the [4 2 1] direction. The computed strain-stress response agrees well with the experimental data.     

Plastic size effect analysis of lamellar composites using a discrete dislocation plasticity approach

Plastic size effect analysis of lamellar composites consisting of elastic and elastic-plastic layers is performed using a discrete dislocation plasticity approach, which is based on applying periodic homogenization to the superposition method for discrete dislocation plasticity. In this approach, the decomposition of displacements into macro and perturbed components circumvents the calculation of superposing displacement fields induced by dislocations in an infinitely homogeneous medium, resulting in two periodic boundary value problems specialized for the analysis of representative volume elements. The present approach is verified by analyzing a model lamellar composite that includes edge dislocations fixed at interfaces. The plastic size effects due to dislocation pile-ups at interfaces are also analyzed. The analysis shows that, strain hardening in elastic-plastic layers arises depending on two factors, namely the thickness and stiffness of elastic layers; and the gap between slip planes in adjacent elastic-plastic layers. In the case where the thickness of elastic layers is several dozen nm, strain hardening in elastic-plastic layers is restrained as the gap of the slip planes decreases. This particular effect is attributed to the long range stress due to pile-ups in adjacent elastic-plastic layers.     

A discrete dislocation analysis of bending

Bending of a strip in plane strain is analyzed using discrete dislocation plasticity where the dislocations are modeled as line defects in a linear elastic medium. At each stage of loading, superposition is used to represent the solution in terms of the infinite medium solution for the discrete dislocations and a complementary solution that enforces the boundary conditions, which is non-singular and obtained from a linear elastic, finite element solution. The lattice resistance to dislocation motion, dislocation nucleation and dislocation annihilation are incorporated into the formulation through a set of constitutive rules. Solutions for cases with multiple slip systems and with a single slip system are presented. The bending moment versus rotation relation and the evolution of the dislocation structure are outcomes of the boundary value problem solution. The effects of slip geometry, obstacles to dislocation motion and specimen size on the moment versus rotation response are considered. Also, the evolution of the dislocation structure is studied with emphasis on the role of geometrically necessary dislocations. The dislocation structure that develops leads to well-defined slip bands, with the slip band spacing scaling with the specimen height.     

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

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.     

Direct calculations of material parameters for gradient plasticity

Length scale parameters introduced in gradient theories of plasticity are calculated in closed form with a continuum dislocation based theory. The similarity of the governing equations in both models for the evolution of plastic deformation of a constrained thin film makes it possible to identify parameters of the gradient plasticity theory with the dislocation based model. A one-to-one identification is not possible given that gradient plasticity does not account for individual dislocations. However, by comparing the mean plastic deformation across the film thickness we find that the length scale parameter, l, introduced in the gradient plasticity theory depends on the geometry as well as material constants.     

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.     

Micro-plasticity of surface steps under adhesive contact: Part I—Surface yielding controlled by single-dislocation nucleation

Mechanics of nano- and meso-scale contacts of rough surfaces is of fundamental importance in understanding deformation and failure mechanisms of a solid surface, and in engineering fabrication and reliability of small surface structures. We present a micro-mechanical dislocation model of contact-induced deformation of a surface step or ledge, as a unit process model to construct a meso-scale model of plastic deformations near and at a rough surface. This paper (Part I) considers onset of contact-induced surface yielding controlled by single-dislocation nucleation from a surface step. The Stroh formalism of anisotropic elasticity and conservation integrals are used to evaluate the driving force on the dislocation. The driving force together with a dislocation nucleation criterion is used to construct a contact-strength map of a surface step in terms of contact pressure, step height, surface adhesion and lattice resistance. Atomistic simulations of atomic surface-step indentation on a gold (1 0 0) surface have been also carried out with the embedded atom method. As predicted by the continuum dislocation model, the atomistic simulations also indicate that surface adhesion plays a significant role in dislocation nucleation processes. Instabilities due to adhesion and dislocation nucleation are evident. The atomistic simulation is used to calibrate the continuum dislocation nucleation criterion, while the continuum dislocation modeling captures the dislocation energetics in the inhomogeneous stress field of the surface-step under contact loading. Results show that dislocations in certain slip planes can be easily nucleated but will stay in equilibrium positions very close to the surface step, while dislocations in some other slip planes easily move away from the surface into the bulk. This phenomenon is called contact-induced near-surface dislocation segregation. As a consequence, we predict the existence of a thin tensile-stress sub-layer adjacent to the surface within the boundary layer of near-surface plastic deformation. In the companion paper (Part II), we analyze the surface hardening behavior caused by multiple dislocations.