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.     

Multiscale modeling of crack initiation and propagation at the nanoscale

Fracture occurs on multiple interacting length scales; atoms separate on the atomic scale while plasticity develops on the microscale. A dynamic multiscale approach (CADD: coupled atomistics and discrete dislocations) is employed to investigate an edge-cracked specimen of single-crystal nickel, Ni, (brittle failure) and aluminum, Al, (ductile failure) subjected to mode-I loading. The dynamic model couples continuum finite elements to a fully atomistic region, with key advantages such as the ability to accommodate 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. An ad hoc computational technique is also applied to dissipate localized waves formed during crack advance in the atomistic zone, whereby an embedded damping zone at the atomistic/continuum interface effectively eliminates the spurious reflection of high-frequency phonons, while allowing low-frequency phonons to pass into the continuum region.The simulations accurately capture the essential physics of the crack propagation in a Ni specimen at different temperatures, including the formation of nano-voids and the sudden acceleration of the crack tip to a velocity close to the material Rayleigh wave speed. The nanoscale brittle fracture happens through the crack growth in the form of nano-void nucleation, growth and coalescence ahead of the crack tip, and as such resembles fracture at the microscale. When the crack tip behaves in a ductile manner, the crack does not advance rapidly after the pre-opening process but is blunted by dislocation generation from its tip. The effect of temperature on crack speed is found to be perceptible in both ductile and brittle specimens.     

A single theory for some quasi-static,supersonic, atomic,and tectonic scale applications of dislocations

We describe a model based on continuum mechanics that reduces the study of a significant class of problems of discrete dislocation dynamics to questions of the modern theory of continuum plasticity. As applications, we explore the questions of the existence of a Peierls stress in a continuum theory, dislocation annihilation, dislocation dissociation, finite-speed-of-propagation effects of elastic waves vis-a-vis dynamic dislocation fields, supersonic dislocation motion, and short-slip duration in rupture dynamics.     

Dislocation-density function dynamics – An all-dislocation,full-dynamics approach for modeling intensive dislocation structures

It has long been recognized that a successful strategy for computational plasticity will have to bridge across the meso scale in which the interactions of high quantities of dislocations dominate. In this work, a new meso-scale scheme based on the full dynamics of dislocation-density functions is proposed. In this scheme, the evolution of the dislocation-density functions is derived from a coarse-graining procedure which clearly defines the relationship between the discrete-line and density representations of the dislocation microstructure. Full dynamics of the dislocation-density functions are considered based on an “all-dislocation” concept in which statistically stored dislocations are preserved and treated in the same way as geometrically necessary dislocations. Elastic interactions between dislocations in a 3D space are treated in accordance with Mura's formula for eigen stress. Dislocation generation is considered as a consequence of dislocations to maintain their connectivity, and a special scheme is devised for this purpose. The model is applied to simulate a number of intensive microstructures involving discrete dislocation events, including loop expansion and shrinkage under applied and self stress, dipole annihilation, and Orowan looping. The scheme can also handle high densities of dislocations present in extensive microstructures.     

Finite element approximation of field dislocation mechanics

A tool for studying links between continuum plasticity and dislocation theory within a field framework is presented. A finite element implementation of the geometrically linear version of a recently proposed theory of field dislocation mechanics (J. Mech. Phys. Solids 49 (2001) 761; Proc. Roy. Soc. 459 (2003) 1343; J. Mech. Phys. Solids 52 (2004) 301) represents the main idea behind the tool. The constitutive ingredients of the theory under consideration are simply elasticity and a specification of dislocation velocity and nucleation. The set of equations to be approximated are non-standard in the context of solid mechanics applications. It comprises the standard second-order equilibrium equations, a first-order div-curl system for the elastic incompatibility, and a first-order, wave-propagative system for the evolution of dislocation density. The latter two sets of equations require special treatment as the standard Galerkin method is not adequate, and are solved utilizing a least-squares finite element strategy. The implementation is validated against analytical results of the classical elastic theory of dislocations and analytical results of the theory itself. Elastic stress fields of dislocation distributions in generally anisotropic media of finite extent, deviation from elastic response, yield-drop, and back-stress are shown to be natural consequences of the model. The development of inhomogeneity, from homogeneous initial conditions and boundary conditions corresponding to homogeneous deformation in conventional plasticity, is also demonstrated. To our knowledge, this work represents the first computational implementation of a theory of dislocation mechanics where no analytical results, singular solutions in particular, are required to formulate the implementation. In particular, a part of the work is the first finite element implementation of Kröner's linear elastic theory of continuously distributed dislocations in its full generality.     

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.     

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.     

Atomistic/continuum simulation of interfacial fracture Part II: Atomistic/dislocation/continuum simulation

Coupled atomistic/dislocation/continuum simulation of interfacial fracture is performed in this paper. The model consists of a nanoscopic core made by atomistic assembly and a surrounding elastic continuum with discrete dislocations. Atomistic dislocations nucleate from the crack tip and move to the continuum layer where they glide according to the dislocation dynamics curve. An atoms/continuum averlapping belt is devised to facilitate the transition between the two scales. The continuum constraint on the atomic assembly is imposed through the mechanics atmosphere along the overlapping belt. Transmissions of mechanics parameters such as displacements, stresses, masses and momenta across the belt are realized. The present model allows us to explore interfacial fracture processes under different mode mixity. The effect of atomistic zigzag interface on the fracture process is revealed: it hinders dislocation emission from the crack tip, especially under high mode mixity. The project supported by the National Natural Science Foundation of China     

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.     

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.     

Multi-scale plasticity modeling: Coupled discrete dislocation and continuum crystal plasticity

A hierarchical multi-scale model that couples a region of material described by discrete dislocation plasticity to a surrounding region described by conventional crystal plasticity is presented. The coupled model is aimed at capturing non-classical plasticity effects such as the long-range stresses associated with a density of geometrically necessary dislocations and source limited plasticity, while also accounting for plastic flow and the associated energy dissipation at much larger scales where such non-classical effects are absent. The key to the model is the treatment of the interface between the discrete and continuum regions, where continuity of tractions and displacements is maintained in an average sense and the flow of net Burgers vector is managed via “passing” of discrete dislocations. The formulation is used to analyze two plane strain problems: (i) tension of a block and (ii) crack growth under mode I loading with various sizes of the discrete dislocation plasticity region surrounding the crack tip. The computed crack growth resistance curves are nearly independent of the size of the discrete dislocation plasticity region for region sizes ranging from to . The multi-scale model can reduce the computational time for the mode I crack analysis by a factor of 14 with little or no loss of fidelity in the crack growth predictions.     

A statistical analysis of the elastic distortion and dislocation density fields in deformed crystals

The statistical properties of the elastic distortion fields of dislocations in deforming crystals are investigated using the method of discrete dislocation dynamics to simulate dislocation structures and dislocation density evolution under tensile loading. Probability distribution functions (PDF) and pair correlation functions (PCF) of the simulated internal elastic strains and lattice rotations are generated for tensile strain levels up to 0.85%. The PDFs of simulated lattice rotation are compared with sub-micrometer resolution three-dimensional X-ray microscopy measurements of rotation magnitudes and deformation length scales in 1.0% and 2.3% compression strained Cu single crystals to explore the linkage between experiment and the theoretical analysis. The statistical properties of the deformation simulations are analyzed through determinations of the Nye and Kröner dislocation density tensors. The significance of the magnitudes and the length scales of the elastic strain and the rotation parts of dislocation density tensors are demonstrated, and their relevance to understanding the fundamental aspects of deformation is discussed.     

A discrete mechanics approach to dislocation dynamics in BCC crystals

A discrete mechanics approach to modeling the dynamics of dislocations in BCC single crystals is presented. Ideas are borrowed from discrete differential calculus and algebraic topology and suitably adapted to crystal lattices. In particular, the extension of a crystal lattice to a CW complex allows for convenient manipulation of forms and fields defined over the crystal. Dislocations are treated within the theory as energy-minimizing structures that lead to locally lattice-invariant but globally incompatible eigendeformations. The discrete nature of the theory eliminates the need for regularization of the core singularity and inherently allows for dislocation reactions and complicated topological transitions. The quantization of slip to integer multiples of the Burgers’ vector leads to a large integer optimization problem. A novel approach to solving this NP-hard problem based on considerations of metastability is proposed. A numerical example that applies the method to study the emanation of dislocation loops from a point source of dilatation in a large BCC crystal is presented. The structure and energetics of BCC screw dislocation cores, as obtained via the present formulation, are also considered and shown to be in good agreement with available atomistic studies. The method thus provides a realistic avenue for mesoscale simulations of dislocation based crystal plasticity with fully atomistic resolution.     

Mechanical annealing under low-amplitude cyclic loading in micropillars

Mechanical annealing has been demonstrated to be an effective method for decreasing the overall dislocation density in submicron single crystal. However, simultaneously significant shape change always unexpectedly happens under extremely high monotonic loading to drive the pre-existing dislocations out of the free surfaces. In the present work, through in situ TEM experiments it is found that cyclic loading with low stress amplitude can drive most dislocations out of the submicron sample with virtually little change of the shape. The underlying dislocation mechanism is revealed by carrying out discrete dislocation dynamic (DDD) simulations. The simulation results indicate that the dislocation density decreases within cycles, while the accumulated plastic strain is small. By comparing the evolution of dislocation junction under monotonic, cyclic and relaxation deformation, the cumulative irreversible slip is found to be the key factor of promoting junction destruction and dislocation annihilation at free surface under low-amplitude cyclic loading condition. By introducing this mechanics into dislocation density evolution equations, the critical conditions for mechanical annealing under cyclic and monotonic loadings are discussed. Low-amplitude cyclic loading which strengthens the single crystal without seriously disturbing the structure has the potential applications in the manufacture of defect-free nano-devices.     

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.