Beñat Gurrutxaga-Lerma

A dynamic discrete dislocation plasticity method for the simulation of plastic relaxation under shock loading

In this article, it is demonstrated that current methods of modelling plasticity as the collective motion of discrete dislocations, such as two-dimensional discrete dislocation plasticity (DDP), are unsuitable for the simulation of very high strain rate processes (106 s−1 or more) such as plastic relaxation during shock loading. Current DDP models treat dislocations quasi-statically, ignoring the time-dependent nature of the elastic fields of dislocations. It is shown that this assumption introduces unphysical artefacts into the system when simulating plasticity resulting from shock loading. This deficiency can be overcome only by formulating a fully time-dependent elastodynamic description of the elastic fields of discrete dislocations. Building on the work of Markenscoff & Clifton, the fundamental time-dependent solutions for the injection and non-uniform motion of straight edge dislocations are presented. The numerical implementation of these solutions for a single moving dislocation and for two annihilating dislocations in an infinite plane are presented. The application of these solutions in a two-dimensional model of time-dependent plasticity during shock loading is outlined here and will be presented in detail elsewhere.

The role of the mobility law of dislocations in the plastic response of shock loaded pure metals

This article examines the role that the choice of a dislocation mobility law has in the study of plastic relaxation at shock fronts. Five different mobility laws, two of them phenomenological fits to data, and three more based on physical models of dislocation inertia, are tested by employing dynamic discrete dislocation plasticity (D3P) simulations of a shock loaded aluminium thin foil. It is found that inertial laws invariably entail very short acceleration times for dislocations changing their kinematic state. As long as the mobility laws describe the same regime of terminal speeds, all mobility laws predict the same degree of plastic relaxation at the shock front. This is used to show that the main factor affecting plastic relaxation at the shock front is in fact the speed of dislocations.

Elastodynamic image forces on dislocations.

The elastodynamic image forces on edge and screw dislocations in the presence of a planar-free surface are derived. The explicit form of the elastodynamic fields of an injected, quiescent screw dislocation are also derived. The resulting image forces are affected by retardation effects: the dislocations experience no image force for a period of time defined by the arrival and reflection at the free surface of the dislocation fields. For the case of injected, stationary dislocations, it is shown that the elastodynamic image force tends asymptotically to the elastotatic prediction. For the case of injected, moving dislocations, it is shown that the elastodynamic image force on both the edge and the screw dislocations is magnified by inertial effects, and becomes increasingly divergent with time; this additional effect, missing in the elastostatic description, is shown to be substantial even for slow moving dislocations. Finally, it is shown that the elastodynamic image force of an edge dislocation moving towards the surface at the Rayleigh wave speed becomes repulsive, rather than attractive; this is suggestive of instabilities at the core of the dislocation, and likely resonances with the free surface.

The Role of Homogeneous Nucleation in Planar Dynamic Discrete Dislocation Plasticity

© 2015 by ASME. Homogeneous nucleation of dislocations is the dominant dislocation generation mechanism at strain rates above 10 8 s -1 ; at those rates, homogeneous nucleation dominates the plastic relaxation of shock waves in the same way that Frank-Read sources control the onset of plastic flow at low strain rates. This article describes the implementation of homogeneous nucleation in dynamic discrete dislocation plasticity (D3P), a planar method of discrete dislocation dynamics (DDD) that offers a complete elastodynamic treatment of plasticity. The implemented methodology is put to the test by studying four materials - Al, Fe, Ni, and Mo - that are shock loaded with the same intensity and a strain rate of 10 10 S -1 . It is found that, even for comparable dislocation densities, the lattice shear strength is fundamental in determining the amount of plastic relaxation a material displays when shock loaded.

Dynamic Discrete Dislocation Plasticity

Abstract This chapter concerns with dynamic discrete dislocation plasticity (D3P), a two-dimensional method of discrete dislocation dynamics aimed at the study of plastic relaxation processes in crystalline materials subjected to weak shock loading. Traditionally, the study of plasticity under weak shock loading and high strain rate has been based on direct experimental measurement of the macroscopic response of the material. Using these data, well-known macroscopic constitutive laws and equations of state have been formulated. However, direct simulation of dislocations as the dynamic agents of plastic relaxation in those circumstances remains a challenge. In discrete dislocation dynamics (DDD) methods, in particular the two-dimensional discrete dislocation plasticity (DDP), the dislocations are modeled as discrete discontinuities in an elastic continuum. However, current DDP and DDD methods are unable to adequately simulate plastic relaxation because they treat dislocation motion quasi-statically, thus neglecting the time-dependent nature of the elastic fields and assuming that they instantaneously acquire the shape and magnitude predicted by elastostatics. This chapter reproduces the findings by Gurrutxaga-Lerma et al. (2013) , who proved that under shock loading, this assumption leads to models that invariably break causality, introducing numerous artifacts that invalidate quasi-static simulation techniques. This chapter posits that these limitations can only be overcome with a fully time-dependent formulation of the elastic fields of dislocations. In this chapter, following the works of Markenscoff & Clifton (1981) and Gurrutxaga-Lerma et al. (2013) , a truly dynamic formulation for the creation, annihilation, and nonuniform motion of straight edge dislocations is derived. These solutions extend the DDP framework to a fully elastodynamic formulation that has been called dynamic discrete dislocation plasticity (D3P). This chapter describes the several changes in paradigm with respect to DDP and DDD methods that D3P introduces, including the retardation effects in dislocation interactions and the effect of the dislocation’s past history. The chapter then builds an account of all the methodological aspects of D3P that have to be modified from DDP, including mobility laws, generation rules, etc. Finally, the chapter explores the applications D3P has to the study of plasticity under shock loading.

Elastodynamic image forces on screw dislocations in the presence of phase boundaries.

The elastodynamic image forces acting on straight screw dislocations in the presence of planar phase boundaries are derived. Two separate dislocations are studied: (i) the injected, non-moving screw dislocation and (ii) the injected (or pre-existing), generally non-uniformly moving screw dislocation. The image forces are derived for both the case of a rigid surface and of a planar interface between two homogeneous, isotropic phases. The case of a rigid interface is shown to be solvable employing Heads image dislocation construction. The case of the elastodynamic image force due to an interface is solved by deriving the reflected waves contribution to the global solution across the interface. This entails obtaining the fundamental solution (Greens function) for a point unit force via Cagniards method, and then applying the convolution theorem for a screw dislocation modelled as a force distribution. Complete, explicit formulae are provided when available. It is shown that the elastodynamic image forces are generally affected by retardation effects, and that those acting on the moving dislocations display a dynamic magnification that exceed the attraction (or repulsion) predicted in classical elastostatic calculations.