skip to main content


Title: On the computational solution of vector-density based continuum dislocation dynamics models: a comparison of two plastic distortion and stress update algorithms
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
Award ID(s):
1663311
NSF-PAR ID:
10283948
Author(s) / Creator(s):
Date Published:
Journal Name:
International journal of plasticity
Volume:
138
ISSN:
0749-6419
Page Range / eLocation ID:
102943
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The inherent inconsistency in identifying the phase field in the phase field crystal theory with the material mass and, simultaneously, with material distortion is discussed. In its current implementation, elastic relaxation in the phase field crystal occurs on a diffusive time scale through a dissipative permeation mode. The very same phase field distortion that is included in solid elasticity drives diffusive motion, resulting in a non physical relaxation of the phase field crystal. We present two alternative theories to remedy this shortcoming. In the first case, it is assumed that the phase field only determines the incompatible part of the elastic distortion, and therefore one is free to specify an additional compatible distortion so as to satisfy mechanical equilibrium at all times (in the quasi static limit). A numerical solution of the new model for the case of a dislocation dipole shows that, unlike the classical phase field crystal model, it can account for the known law of relative motion of the two dislocations in the dipole. The physical origin of the compatible strain in this new theory remains to be specified. Therefore, a second theory is presented in which an explicit coupling between independent distortion and phase field accounts for the time dependence of the relaxation of fluctuations in both. Preliminary details of its implementation are also given. 
    more » « less
  2. 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
  3. 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
  4. 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
  5. A phase-field model for thermomechanically-induced fracture in NiTi at the single crystal level, i.e., fracture under loading paths that may take advantage of either of the functional properties of NiTi–superelasticity or shape memory effect–, is presented, formulated within the kinematically linear regime. The model accounts for reversible phase transformation from austenite to martensite habit plane variants and plastic deformation in the austenite phase. Transformation-induced plastic deformation is viewed as a mechanism for accommodation of the local deformation incompatibility at the austenite–martensite interfaces and is accounted for by introducing an interaction term in the free energy derived based on the Mori–Tanaka and Kröner micromechanical assumptions and the hypothesis of martensite instantaneous growth within austenite. Based on experimental observations suggesting that NiTi fractures in a stress-controlled manner, damage is assumed to be driven by the elastic energy, i.e., phase transformation and plastic deformation are assumed to contribute in crack formation and growth indirectly through stress redistribution. The model is restricted to quasistatic mechanical loading (no latent heat effects), thermal loading sufficiently slow with respect to the time rate of heat transfer by conduction (no thermal gradients), and a temperature range below 𝑀𝑑, which is the temperature above which the austenite phase is stable, i.e., stress-induced martensitic transformation is suppressed. The numerical implementation of the model is based on an efficient scheme of viscous regularization in both phase transformation and plastic deformation, an explicit numerical integration via a tangent modulus method, and a staggered scheme for the coupling of the unknown fields. The model is shown able to capture transformation-induced toughening, i.e., stable crack advance attributed to the shielding effect of inelastic deformation left in the wake of the growing crack under nominal isothermal loading, actuation-induced fracture under a constant bias load, and crystallographic dependence on crack pattern. 
    more » « less