- Award ID(s):
- 1663311
- NSF-PAR ID:
- 10283948
- 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
-
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
-
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
-
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.
-
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
-
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