Search for: All records

For the past century, dislocations have been understood to be the carriers of plastic deformation in crystalline solids. However, their collective behavior is still poorly understood. Progress in understanding the collective behavior of dislocations has primarily come in one of two modes: the simulation of systems of interacting discrete dislocations and the treatment of density measures of varying complexity that are considered as continuum fields. A summary of contemporary models of continuum dislocation dynamics is presented. These include, in order of complexity, the two-dimensional statistical theory of dislocations, the field dislocation mechanics treating the total Kröner–Nye tensor, vector density approaches that treat geometrically necessary dislocations on each slip system of a crystal, and high-order theories that examine the effect of dislocation curvature and distribution over orientation. Each of theories contain common themes, including statistical closure of the kinetic dislocation transport equations and treatment of dislocation reactions such as junction formation. An emphasis is placed on how these common themes rely on closure relations obtained by analysis of discrete dislocation dynamics experiments. The outlook of these various continuum theories of dislocation motion is then discussed. 

An accurate description of the evolution of dislocation networks is an essential part of discrete and continuum dislocation dynamics models. These networks evolve by motion of the dislocation lines and by forming junctions between these lines via cross slip, annihilation and junction reactions. In this work, we introduce these dislocation reactions into continuum dislocation models using the theory of de Rham currents. We introduce dislocations on each slip system as potentially open lines whose boundaries are associated with junction points and, therefore, still create a network of collectively closed lines that satisfy the classical relations and for the dislocation density tensor and the plastic distortion . To ensure this, we leverage Frank’s second rule at the junction nodes and the concept of virtual dislocation segments. We introduce the junction point density as a new state variable that represents the distribution of junction points within the crystal containing the dislocation network. Adding this information requires knowledge of the global structure of the dislocation network, which we obtain from its representation as a graph. We derive transport relations for the dislocation line density on each slip system in the crystal, which now includes a term that corresponds to the motion of junction points. We also derive the transport relations for junction points, which include source terms that reflect the topology changes of the dislocation network due to junction formation. 

The sensitivity of recrystallization kinetics in metals to the heterogeneity of microstructure and deformation history is a widely accepted experimental fact. However, most of the available recrystallization models employ either a mean field approach or use grain-averaged parameters, and thus neglecting the mesoscopic heterogeneity induced by prior deformation. In the present study, we investigate the impact of deformation-induced dislocation (subgrain) structure on the kinetics of recrystallization in metals using the phase-field approach. The primary focus here is upon the role of dislocation cell boundaries. The free energy formulation of the phase-field model accounts for the heterogeneity of the microstructure by assigning localized energy to the resulting dislocation microstructure realizations generated from experimental data. These microstructure realizations are created using the universal scaling laws for the spacing and the misorientation angles of both the geometrically necessary and incidental dislocation boundaries. The resulting free energy is used into an Allen-Cahn based model of recrystallization kinetics, which are solved using the finite element method. The solutions thus obtained shed light on the critical role of the spatial heterogeneity of deformation in the non-smooth growth of recrystallization nuclei and on the final grain structure. The results showed that, in agreement with experiment, the morphology of recrystallization front exhibits protrusions and retrusions. By resolving the subgrain structure, the presented algorithm paves the way for developing predictive kinetic models that fully account for the deformed state of recrystallizing metals. 

Continuum dislocation dynamics models of mesoscale plasticity consist of dislocation transport-reaction equations coupled with crystal mechanics equations. The coupling between these two sets of equations is such that dislocation transport gives rise to the evolution of plastic distortion (strain), while the evolution of the latter fixes the stress from which the dislocation velocity field is found via a mobility law. Earlier solutions of these equations employed a staggered solution scheme for the two sets of equations in which the plastic distortion was updated via time integration of its rate, as found from Orowan’s law. In this work, we show that such a direct time integration scheme can suffer from accumulation of numerical errors. We introduce an alternative scheme based on field dislocation mechanics that ensures consistency between the plastic distortion and the dislocation content in the crystal. The new scheme is based on calculating the compatible and incompatible parts of the plastic distortion separately, and the incompatible part is calculated from the current dislocation density field. Stress field and dislocation transport calculations were implemented within a finite element based discretization of the governing equations, with the crystal mechanics part solved by a conventional Galerkin method and the dislocation transport equations by the least squares method. A simple test is first performed to show the accuracy of the two schemes for updating the plastic distortion, which shows that the solution method based on field dislocation mechanics is more accurate. This method then was used to simulate an austenitic steel crystal under uniaxial loading and multiple slip conditions. By considering dislocation interactions caused by junctions, a hardening rate similar to discrete dislocation dynamics simulation results was obtained. The simulations show that dislocations exhibit some self-organized structures as the strain is increased. 

The equations of dislocation transport at finite crystal deformation were developed, with a special emphasis on a vector density representation of dislocations. A companion thermodynamic analysis yielded a generalized expression for the driving force of dislocations that depend on Mandel (Cauchy) stress in the reference (spatial) configurations and the contribution of the dislocation core energy to the free energy of the crystal. Our formulation relied on several dislocation density tensor measures linked to the incompatibility of the plastic distortion in the crystal. While previous works develop such tensors starting from the multiplicative decomposition of the deformation gradient, we developed the tensor measures of the dislocation density and the dislocation flux from the additive decomposition of the displacement gradient and the crystal velocity fields. The two-point dislocation density measures defined by the referential curl of the plastic distortion and the spatial curl of the inverse elastic distortion and the associate dislocation currents were found to be more useful in deriving the referential and spatial forms of the transport equations for the vector density of dislocations. A few test problems showing the effect of finite deformation on the static dislocation fields are presented, with a particular attention to lattice rotation. The framework developed provides the theoretical basis for investigating crystal plasticity and dislocation patterning at the mesoscale, and it bears the potential for realistic comparison with experiments upon numerical solution.