An official website of the United States government
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
A Line-Free Discrete Dislocation Dynamics Method for Finite Domains
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
- Award ID(s):
- 1950027
- PAR ID:
- 10599384
- Publisher / Repository:
- Springer Nature Switzerland
- Date Published:
- ISBN:
- 978-3-031-50349-8
- Page Range / eLocation ID:
- 825 to 832
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Mascia, Corrado; Terracina, Andrea; Tesei, Alberto (Ed.)We study a model of dislocations in two-dimensional elastic media. In this model, the displacement satisfies the system of linear elasticity with mixed displacement-traction homogeneous boundary conditions in the complement of an open curve in a bounded planar domain, and has a specified jump, the slip, across the curve, while the traction is continuous there. The stiffness tensor is allowed to be anisotropic and inhomogeneous. We prove well-posedness of the direct problem in a variational setting, assuming the coefficients are Lipschitz continuous. Using unique continuation arguments, we then establish uniqueness in the inverse problem of determining the dislocation curve and the slip from a single measurement of the displacement on an open patch of the traction-free part of the boundary. Uniqueness holds when the elasticity operators admits a suitable decomposition and the curve satisfies additional geometric assumptions. This work complements the results in Arch. Ration. Mech. Anal., 236(1):71-111, (2020), and in Preprint arXiv:2004.00321, which concern three-dimensional isotropic elastic media.more » « less
-
We consider a model for elastic dislocations in geophysics. We model a portion of the Earth’s crust as a bounded, inhomogeneous elastic body with a buried fault surface, along which slip occurs. We prove well-posedness of the resulting mixed-boundary-value-transmission problem, assuming only bounded elastic moduli. We establish uniqueness in the inverse problem of determin- ing the fault surface and the slip from a unique measurement of the displacement on an open patch at the surface, assuming in addition that the Earth’s crust is an isotropic, layered medium with Lamé coefficients piecewise Lipschitz on a known partition and that the fault surface satisfies certain geo- metric conditions. These results substantially extend those of the authors in [Arch. Ration. Mech. Anal. 236, 71–111 (2020)].more » « less
-
ABSTRACT We investigate the appropriateness of calving or crevasse models from the literature using linear elastic fracture mechanics (LEFM). To this end, we compare LEFM model-predicted stress intensity factors (SIFs) against numerically computed SIFs using the displacement correlation method in conjunction with the finite element method. We present several benchmark simulations wherein we calculate the SIF at the tips of water-filled surface and basal crevasses penetrating through rectangular ice slabs under different boundary conditions, including grounded and floating conditions. Our simulation results indicate that the basal boundary condition significantly influences the SIF at the crevasse tips. We find that the existing calving models using LEFM are not generally accurate for evaluating SIFs in grounded glaciers or floating ice shelves. We also illustrate that using the ‘single edge crack’ weight function in the LEFM formulations may be appropriate for predicting calving from floating ice shelves, owing to the low fracture toughness of ice; whereas, using the ‘double edge crack’ or ‘central through crack’ weight functions is more appropriate for predicting calving from grounded glaciers. To conclude, we recommend using the displacement correlation method for SIF evaluation in real glaciers and ice shelves with complex geometries and boundary conditions.more » « less
-
Background. An assumption of Digital Image Correlation (DIC) is that the displacement field within each subset is relatively smooth, captured with reasonable accuracy by, for example, linear or quadratic shape functions. Although this assumption works well for many materials, it becomes problematic for heterogeneous materials, such as fiber networks, wherein the length scale of heterogeneity matches the size of a subset. Objective. Here we applied DIC to fibrous networks made of collagen, for which displacements at the scale of a subset are highly heterogeneous, but errors caused by the heterogeneity are difficult to quantify. We developed a method to quantify such errors. Methods. We began by generating a synthetic three-dimensional fiber network with structure matching that of gels made of fibrous collagen. We then formulated an algorithm to mimic the way in which a confocal microscope images the fibers at its focal plane, thereby generating synthetic images similar to those obtained in experiments. Displacement boundary conditions were applied to the synthetic fiber networks, and the resulting displacement fields were computed using a finite element solver. DIC was applied to the synthetic images, and displacements were compared to the data from the finite element method, enabling rigorous quantification of error. Results. Point-wise errors in the DIC-measured displacements were substantial, often exceeding 40%, but over regions far larger than the length scales of heterogeneity or the DIC subset size, errors were modest, e.g., ≤15%. Conclusions. Although DIC can accurately measure displacements of fiber networks at length scales larger than the subset window, quantification of mechanical behavior at the scale of material heterogeneity will require new methods to complement or replace the use of DIC.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Book Chapter:
- https://doi.org/10.1007/978-3-031-50349-8_71
