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.
more »
« less
Theoretical development of continuum dislocation dynamics for finite-deformation crystal plasticity at the mesoscale
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.
more »
« less
- Award ID(s):
- 1663311
- NSF-PAR ID:
- 10188228
- Date Published:
- Journal Name:
- Journal of the mechanics and physics of solids
- Volume:
- 139
- ISSN:
- 0022-5096
- Page Range / eLocation ID:
- 103926
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)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.more » « less
-
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.more » « less
-
more » « less
Abstract The deformation of crystalline materials by dislocation motion takes place in discrete amounts determined by the Burgers vector. Dislocations may move individually or in bundles, potentially giving rise to intermittent slip. This confers plastic deformation with a certain degree of variability that can be interpreted as being caused by stochastic fluctuations in dislocation behavior. However, crystal plasticity (CP) models are almost always formulated in a continuum sense, assuming that fluctuations average out over large material volumes and/or cancel out due to multi-slip contributions. Nevertheless, plastic fluctuations are known to be important in confined volumes at or below the micron scale, at high temperatures, and under low strain rate/stress deformation conditions. Here, we develop a stochastic solver for CP models based on the residence-time algorithm that naturally captures plastic fluctuations by sampling among the set of active slip systems in the crystal. The method solves the evolution equations of explicit CP formulations, which are recast as stochastic ordinary differential equations and integrated discretely in time. The stochastic CP model is numerically stable by design and naturally breaks the symmetry of plastic slip by sampling among the active plastic shear rates with the correct probability. This can lead to phenomena such as intermittent slip or plastic localization without adding external symmetry-breaking operations to the model. The method is applied to body-centered cubic tungsten single crystals under a variety of temperatures, loading orientations, and imposed strain rates.
-
more » « less
Abstract The severe difficulty to resolve simultaneously both the macroscopic deformation process and the dislocation dynamics on the atomic scale limits our understanding of crystal plasticity. Here we use colloidal crystals, imaged on the single particle level by high-speed three-dimensional (3D) confocal microscopy, and resolve in real-time both the relaxation of the epitaxial misfit strain and the accompanying evolution of dislocations. We show how dislocation interactions give rise to the formation of complex dislocation networks in 3D and to unexpectedly sharp plastic relaxation. The sharp relaxation is facilitated by attractive interactions that promote the formation of new dislocations that are more efficient in mediating strain. Dislocation networks form fragmented structures, as dislocation growth is blocked by either attractive interactions, which result in the formation of sessile dislocation junctions, or by repulsion from perpendicular segments. The strength of these blocking mechanisms decreases with the thickness of the crystal film. These results reveal the critical role of dislocation interactions in plastic deformation of thin films and can be readily generalized from the colloidal to the atomic scale.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1016/j.jmps.2020.103926
