An official website of the United States government
A method for solving general boundary-value problems involving discrete dislocations is introduced. Plastic flow emerges from the motion of dislocations in an incremental fashion. At each increment, the displacement, strain and stress fields in the body are obtained by superposition of the infinite medium fields associated with individual dislocations and an image field that enforces boundary conditions. Dislocations are represented as monopoles and dislocation events are treated as a transportation map problem. Long-range interactions are accounted for through linear elasticity with a core regularization procedure. At the current state of development of the method, no ad hoc short-range interactions are included. An approximate loop nucleation model is used for large-scale computations. The image problem is solved using a finite element formulation with the following features: (i) a single Cholesky decomposition of the global stiffness matrix, (ii) a consistent enforcement of traction and displacement boundary conditions, and (iii) image force interpolation using an efficient BB-tree algorithm. To ensure accuracy, we explore stable time steps and employ monopole splitting techniques. Special attention is given to the interaction of curved dislocations with arbitrary domain boundaries and free surfaces. The capabilities of the framework are illustrated through a wire torsion problem.
more »
« less
A discrete dislocation analysis of size-dependent plasticity in torsion
A method for solving three dimensional discrete dislocation plasticity boundary-value problems using a monopole representation of the dislocations is presented. At each time step, the displacement, strain and stress fields in a finite body are obtained by superposition of infinite body dislocation fields and an image field that enforces the boundary conditions. The three dimensional infinite body fields are obtained by representing dislocations as being comprised of points, termed monopoles, that carry dislocation line and Burgers vector information. The image fields are obtained from a three dimensional linear elastic finite element calculation. The implementation of the coupling of the monopole representation with the finite element method, including the interaction of curved dislocations with free surfaces, is presented in some detail because it differs significantly from an implementation with a line based dislocation representation. Numerical convergence and the modeling of dislocation loop nucleation for large scale computations are investigated. The monopole discrete dislocation plasticity framework is used to investigate the effect of size and initial dislocation density on the torsion of wires with diameters varying over three orders of magnitude. Depending on the initial dislocation source density and the wire diameter, three regimes of torsion–twist response are obtained: (i) for wires with a sufficiently small diameter, plastic deformation is nucleation controlled and is strongly size dependent; (ii) for wires with larger diameters dislocation plasticity is dislocation interaction controlled, with the emergence of geometrically necessary dislocations and dislocation pile-ups playing a key role, and is strongly size dependent; and (iii) for wires with sufficiently large diameters plastic deformation becomes less heterogeneous and the dependence on size is greatly diminished.
more »
« less
- Award ID(s):
- 1950027
- PAR ID:
- 10599383
- Publisher / Repository:
- Elsevier
- Date Published:
- Journal Name:
- Journal of the Mechanics and Physics of Solids
- Volume:
- 190
- Issue:
- C
- ISSN:
- 0022-5096
- Page Range / eLocation ID:
- 105709
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Plasticity in body centered cubic (BCC) crystals is shown to be controlled by slow screw dislocation motion, owing to the thermally-activated process of kink pair nucleation and migration. Through three dimensional discrete dislocation dynamics simulations, this work unravels the mystery of how such slow screw dislocation behavior contributes to extremely rapid strain bursts in submicron BCC tungsten (W) pillars, which is typical of BCC metals. It is found that strain bursts are dominated by the motion of non-screw dislocations at low strain rate, but are more influenced by screw dislocations at high strain rate. The total, and partial strain burst magnitude due to screw dislocations alone, are found to exhibit rate dependence following a power law statistics with exponent of 0.65. Similar power law statistics are also obeyed for the standard deviation of the corresponding plastic strain rate. The role of screw dislocations is attributed to the changing nature of dislocation source operation at different strain rates. The corresponding spatial distribution of plastic deformation is also discussed based on the uniqueness of the simulation method in reproducing the distribution of slipped area and plastic strain with very high spatial resolution.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
-
A novel dislocation-density-based crystal plasticity model for nanocrystalline face-centered cubic metals is developed based on the thermally-activated mechanism of dislocations depinning from grain boundaries. Dislocations nucleated from grain boundary dislocation sources are assumed to be the primary carriers of plasticity in the nanocrystals. The evolution of the dislocation density thereby involves a competition between the nucleation of dislocations from grain boundary defect structures, such as ledges, and the absorption of dislocations into the grain boundary via diffusion processes. This model facilitates the simulation of plastic deformation in nanocrystalline metals, with consideration of the initial microstructure resulting from a particular processing method, to be computed as a direct result of dislocation-mediated plasticity only. The exclusion of grain boundary-mediated plasticity mechanisms in the formulation of the crystal plasticity model allows for the exploration of the fundamental role dislocations play in nanocrystalline plasticity. The combined effect of average grain size, grain size distribution shape, and initial dislocation density on the mechanical performance and strain-rate sensitivity are explored with the model. Further, the influence of the grain boundary diffusivity on post-yielding strain-hardening behavior is investigated to discern the impact that the choice of processing route has on the resulting deformation response of the material.more » « less
-
This work presents a multiscale study of the uniaxial compression of Si pillars, with diameters ranging from 50 nm to 360 nm, using the Concurrent Atomistic-Continuum (CAC) method. The simulations reproduce the brittle and ductile deformation behaviors of Si pillars observed in experiments. For defect-free Si pillars compressed by a perfectly smooth flat punch with a repulsive force field to reflect an assumed rigid indenter, dislocations are nucleated from the corner of the bottom surface for pillars with diameters of 100 nm and below, while for pillars with diameters of 220 nm and above, dislocations nucleate from the top surface; multiple slip systems are activated in all pillars except for the pillar with a diameter of 50 nm. A strong size effect is thus demonstrated with regard to the nucleation of dislocations. Another key finding is the critical role of defects on the indenter surface. For a perfectly flat indenter, all the defect-free Si pillars with diameters ranging from 50 nm to 360 nm exhibit ductile deformation. By contrast, for an indenter with surface steps, all pillars with diameters of 100 nm and above deform in a brittle manner. These surface steps cause sequential nucleation of dislocations and activation of two slip systems, leading to dislocation intersection and formation of a sessile Lomer lock. Continued pileups of dislocations against the Lomer lock lead to the initiation of a crack at the intersection. The deformation mechanism underlying the crack formation is thus demonstrated.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1016/j.jmps.2024.105709
