Interaction between phase transformations and dislocations at the nanoscale. Part 1. General phase field approach

Thermodynamically consistent, three-dimensional (3D) phase field approach (PFA) for coupled multivariant martensitic transformations (PTs), including cyclic PTs, variant–variant transformations (i.e., twinning), and dislocation evolution is developed at large strains. One of our key points is in the justification of the multiplicative decomposition of the deformation gradient into elastic, transformational, and plastic parts. The plastic part includes four mechanisms: dislocation motion in martensite along slip systems of martensite and slip systems of austenite inherited during PT and dislocation motion in austenite along slip systems of austenite and slip systems of martensite inherited during reverse PT. The plastic part of the velocity gradient for all these mechanisms is defined in the crystal lattice of the austenite utilizing just slip systems of austenite and inherited slip systems of martensite, and just two corresponding types of order parameters. The explicit expressions for the Helmholtz free energy and the transformation and plastic deformation gradients are presented to satisfy the formulated conditions related to homogeneous thermodynamic equilibrium states of crystal lattice and their instabilities. In particular, they result in a constant (i.e., stress- and temperature-independent) transformation deformation gradient and Burgers vectors. Thermodynamic treatment resulted in the determination of the driving forces for change of the order parameters for PTs and dislocations. It also determined the boundary conditions for the order parameters that include a variation of the surface energy during PT and exit of dislocations. Ginzburg–Landau equations for dislocations include variation of properties during PTs, which in turn produces additional contributions from dislocations to the Ginzburg–Landau equations for PTs. A complete system of coupled PFA and mechanics equations is presented. A similar theory can be developed for PFA to dislocations and other PTs, like reconstructive PTs and diffusive PTs described by the Cahn–Hilliard equation, as well as twinning and grain boundaries evolution.     

Discrete Crystal Elasticity and Discrete Dislocations in Crystals

This article is concerned with the development of a discrete theory of crystal elasticity and dislocations in crystals. The theory is founded upon suitable adaptations to crystal lattices of elements of algebraic topology and differential calculus such as chain complexes and homology groups, differential forms and operators, and a theory of integration of forms. In particular, we define the lattice complex of a number of commonly encountered lattices, including body-centered cubic and face-centered cubic lattices. We show that material frame indifference naturally leads to discrete notions of stress and strain in lattices. Lattice defects such as dislocations are introduced by means of locally lattice-invariant (but globally incompatible) eigendeformations. The geometrical framework affords discrete analogs of fundamental objects and relations of the theory of linear elastic dislocations, such as the dislocation density tensor, the equation of conservation of Burgers vector, Kröner's relation and Mura's formula for the stored energy. We additionally supply conditions for the existence of equilibrium displacement fields; we show that linear elasticity is recovered as the Γ-limit of harmonic lattice statics as the lattice parameter becomes vanishingly small; we compute the Γ-limit of dilute dislocation distributions of dislocations; and we show that the theory of continuously distributed linear elastic dislocations is recovered as the Γ-limit of the stored energy as the lattice parameter and Burgers vectors become vanishingly small.     

New inroads in an old subject: Plasticity, from around the atomic to the macroscopic scale

Nonsingular, stressed, dislocation (wall) profiles are shown to be 1-d equilibria of a non-equilibrium theory of Field Dislocation Mechanics (FDM). It is also shown that such equilibrium profiles corresponding to a given level of load cannot generally serve as a travelling wave profile of the governing equation for other values of nearby constant load; however, one case of soft loading with a special form of the dislocation velocity law is demonstrated to have no ‘Peierls barrier’ in this sense. The analysis is facilitated by the formulation of a 1-d, scalar, time-dependent, Hamilton-Jacobi equation as an exact special case of the full 3-d FDM theory accounting for non-convex elastic energy, small, Nye-tensor-dependent core energy, and possibly an energy contribution based on incompatible slip. Relevant nonlinear stability questions, including that of nucleation, are formulated in a non-equilibrium setting. Elementary averaging ideas show a singular perturbation structure in the evolution of the (unsymmetric) macroscopic plastic distortion, thus pointing to the possibility of predicting generally rate-insensitive slow response constrained to a tensorial ‘yield’ surface, while allowing fast excursions off it, even though only simple kinetic assumptions are employed in the microscopic FDM theory. The emergent small viscosity on averaging that serves as the small parameter for the perturbation structure is a robust, almost-geometric consequence of large gradients of slip in the dislocation core and the persistent presence of a large number of dislocations in the averaging volume. In the simplest approximation, the macroscopic yield criterion displays anisotropy based on the microscopic dislocation line and Burgers vector distribution, a dependence on the Laplacian of the incompatible slip tensor and a nonlocal term related to a Stokes-Helmholtz-curl projection of an ‘internal stress’ derived from the incompatible slip energy.     

Dislocation microstructures and strain-gradient plasticity with one active slip plane

We study dislocation networks in the plane using the vectorial phase-field model introduced by Ortiz and coworkers, in the limit of small lattice spacing. We show that, in a scaling regime where the total length of the dislocations is large, the phase field model reduces to a simpler model of the strain-gradient type. The limiting model contains a term describing the three-dimensional elastic energy and a strain-gradient term describing the energy of the geometrically necessary dislocations, characterized by the tangential gradient of the slip. The energy density appearing in the strain-gradient term is determined by the solution of a cell problem, which depends on the line tension energy of dislocations. In the case of cubic crystals with isotropic elasticity our model shows that complex microstructures may form in which dislocations with different Burgers vector and orientation react with each other to reduce the total self-energy.     

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

Analysis of dislocation nucleation from a crystal surface based on the Peierls-Nabarro dislocation model

Dislocation nucleation from a stressed crystal surface is analyzed based on the Peierls-Nabarro dislocation model. The variational boundary integral approach is used to obtain the profiles of the embryonic dislocations in various three-dimensional nucleation configurations. The stress-dependent activation energies required to activate dislocations from their stable to unstable saddle point configurations are determined. Compared to previous analyses of this type of problem based on continuum elastic dislocation theory, the present analysis eliminates the uncertain core cutoff parameter by allowing for the existence of an extended dislocation core as the embryonic dislocation evolves. Moreover, atomic information can be incorporated to reveal the dependence of the nucleation process on the profile of the atomic interlayer potential as compared to continuum elastic dislocation theory in which only elastic constants and Burgers vector are relevant. Finally, the presented methodology can also be readily used to study dislocation nucleation from the surface heterogeneities such as cracks, steps, and quantum structures of electronic devices.     

The Force on a Dislocation Near a Weakly Bonded Interface

The physical system being studied here is an elastic dislocation in proximity to a material interface which has limited strength in shear. The configuration is relevant to elastic strain relaxation in bonded thin-film semiconductor material structures. A plane strain boundary value problem is formulated and an exact solution for the force on the dislocation versus distance from the interface is obtained. If the dislocation is close enough to the interface, its stress field causes irreversible slip across the interface. This slip, in turn, induces an attractive force on the dislocation. For the case of an edge dislocation with Burgers vector normal to the interface, it is found that the dislocation has a stable equilibrium position at a small distance from the interface. This revised version was published online in August 2006 with corrections to the Cover Date.     

A crystal plasticity model that incorporates stresses and strains due to slip gradients

This work is concerned with incorporating the kinematic and stress effects of excess dislocations in a constitutive model for the elastoplastic behavior of crystalline materials. The foundation of the model is a three term multiplicative decomposition of the deformation gradient in which the two classical terms of plastic and elastic deformation are included along with an additional term for long range strain due to the collective effects of excess dislocations. The long range strain is obtained from an assumed density of Volterra edge dislocations and is directly related to gradients in slip. A new material parameter emerges which is the size the region about a continuum point that contributes to long range strains.Using Hookean elasticity, the stress at a point is linearly related to the sum of the elastic plus the long range strain fields. However, the driving force for slip is postulated to be due only to the elastic stress so that the long range stress is a back stress in the constitutive relationship for plastic deformation. A consistent balance of the total deformation rate with the three proposed mechanisms of deformation leads to a set of differential equations that can be solved for the elastic stress, rotation and pressure which then implicitly defines the material state and equilibrium stress. Results from the simulation of a tapered tensile specimen demonstrate that the constitutive model exhibits isotropic and kinematic type hardening effects as well as changes in the pattern of plastic deformation and necking when compared to a material without slip gradient effects.     

A comparison between crystal and isotropic strain gradient plasticity theories with accent on the role of the plastic spin

In the small deformation range, we consider crystal and isotropic “higher-order” theories of strain gradient plasticity, in which two different types of size effects are accounted for: (i) that dissipative, entering the model through the definition of an effective measure of plastic deformation peculiar of the isotropic hardening function and (ii) that energetic, included by defining the defect energy (i.e., a function of Nye's dislocation density tensor added to the free energy; see, e.g., [Gurtin, M.E., 2002. A gradient theory of single-crystal viscoplasticity that accounts for geometrically necessary dislocations. J. Mech. Phys. Solids 50, 5–32]). In order to compare the two modellings, we recast both of them into a unified deformation theory framework and apply them to a simple boundary value problem for which we can exploit the Γ-convergence results of [Bardella, L., Giacomini, A., 2008. Influence of material parameters and crystallography on the size effects describable by means of strain gradient plasticity. J. Mech. Phys. Solids 56 (9), 2906–2934], in which the crystal model is made isotropic by imposing that any direction be a possible slip system. We show that the isotropic modelling can satisfactorily approximate the behaviour described by the isotropic limit obtained from the crystal modelling if the former constitutively involves the plastic spin, as in the theory put forward in Section 12 of [Gurtin, M.E., 2004. A gradient theory of small-deformation isotropic plasticity that accounts for the Burgers vector and for dissipation due to plastic spin. J. Mech. Phys. Solids 52, 2545–2568]. The analysis suggests a criterium for choosing the material parameter governing the plastic spin dependence into the relevant Gurtin model.     

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.     

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.     

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.     

Nucleation of partial dislocations at a crack and its implication on deformation mechanisms of nanostructured metals

Nucleation of partial dislocations at a crack is analyzed based a multiscale model that incorporates atomic information into continuum-mechanics approach. The crack and the slip plane as the extension of the crack are modeled as a surface of displacement discontinuities embedded in an elastic medium. The atomic potential between the adjacent atomic layers along the slip plane is assumed to be the generalized stacking fault energy, which is obtained based on atomic calculations. The relative displacements along the slip plane, corresponding to the configurations of partial dislocations and stacking faults, are solved through the variational boundary integral method. The energetics of partial dislocation nucleation at the crack in FCC metals Al and Cu are comparatively studied for their distinctive difference in the intrinsic stacking fault energy. Compared with nucleation of perfect dislocations in previous studies, several new features have emerged. Among them, the critical stress and activation energy for nucleation of partial dislocations are markedly lowered. Depending on the value of stacking fault energy and crack configuration, the saddle-point configurations of partial dislocations can be vastly different in terms of the nucleation sequence and the size of the stacking fault. These findings have significant implication for identifying the fundamental dislocation and grain-boundary-mediated deformation mechanisms, which may account for the limiting strength and the high strain rate sensitivity of nanostructured metals.     

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.     

Energy dissipation mechanisms in ductile fracture

The objective is to investigate energy dissipation mechanisms that operate at different length scales during fracture in ductile materials. A dimensional analysis is performed to identify the sets of dimensionless parameters which contribute to energy dissipation via dislocation-mediated plastic deformation at a crack tip. However, rather than using phenomenological variables such as yield stress and hardening modulus in the analysis, physical variables such as dislocation density, Burgers vector and Peierls stress are used. It is then shown via elementary arguments that the resulting dimensionless parameters can be interpreted in terms of competitions between various energy dissipation mechanisms at different length scales from the crack tip; the energy dissipations mechanisms are cleavage, crack tip dislocation nucleation and also dislocation nucleation from a Frank-Read source. Therefore, the material behavior is classified into three groups. The first two groups are the well-known intrinsic brittle and intrinsic ductile behavior. The third group is designated to be extrinsic ductile behavior for which Frank-Read dislocation nucleation is the initial energy dissipation mechanism. It is shown that a material is predicted to exhibit extrinsic ductility if the dimensionless parameter bρ_disl^1/2 (b is Burgers vector, ρ_disl is dislocation density) is within a certain range defined by other dimensionless parameters, irrespective of the competition between cleavage and crack tip dislocation nucleation. The predictions compare favorably to the documented behavior of a number of different classes of materials.