Active learning of continuum-scale mechanics by guiding discrete dislocation dynamics simulations
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Dislocations are the primary carriers of plasticity in metallic materials. Understanding the basic mechanisms for dislocation movement is paramount to predicting the material mechanical response. Relying on atomistic simulations, we observe a transition from non-Arrhenius to Arrhenius behavior in the rate of an edge dislocation overcoming the long-range elastic interaction with a prismatic loop in tungsten. Close to the critical resolved shear stress, the process shows a non-Arrhenius behavior at low temperatures. However, as the temperature increases, the activation entropy starts to dominate, leading to a traditional Arrhenius-like behavior. Here, we have computed the activation entropy analytically along the minimum energy path following Schoeck’s method, which captures the cross-over between anti-Arrhenius and Arrhenius domains. Also, the Projected Average Force Integrator (PAFI), another simulation method to compute free energies along an initial transition path, exhibits considerable concurrence with Schoeck’s formalism. We conclude that entropic effects need to be considered to understand processes involving dislocations bypassing elastic barriers close to the critical resolved shear stress. More work needs to be performed to fully understand the discrepancies between Schoeck’s and PAFI compared to molecular dynamics.
In this article, we present the creep characteristics of two reduced activation ferritic-martensitic steels of identical starting compositions formed by different fabrication routes: a nanostructured ferritic alloy commonly referred to as a castable nanostructured alloy (CNA) and a sintered nanostructured alloy (SNA) variant. Through a series of nanoindentation experiments spanning a temperature range of 25 °C to 650 °C, with a maximum load of 100 mN, we find creep behaviors in the cast and sintered materials to be remarkably similar. The creep stress exponent (n) for CNA and SNA were found to be in the range of 8–35 and the activation volume was ∼14–42b 3 , underscoring a dominance of dislocation-mediated mechanisms in both alloys. Notably, we observed a decline in the creep stress exponent with increasing temperature, attributable to the heightened influence of thermally activated dislocations. This phenomenon suggests a potential transition in the deformation mechanism towards a thermally activated dislocation climb process, significantly impacting the observed creep behavior.
Order strengthening in nickel-based superalloys is associated with the extra stress required for dislocations to bypass the $\gamma'$ precipitates distributed in the $\gamma$ matrix. Depending on the operating conditions and microstructure, a rich variety of bypass mechanism has been identified, with various shearing and Orowan looping processes gradually giving way to climb bypass as the operating conditions change from the low/intermediate temperatures and high stress regime, to the high temperature and low stress regime. When anti phase boundary (APB) shearing and Orowan looping mechanisms operate, the classical picture is that, at for a given volume fraction, the bypass mechanism changes from shearing to looping with increased particle size and within a broad coexistence size window. Another possibility, which is supported by indirect experimental evidence, is that a third “hybrid” transition mechanism may operate. Here, in this paper we use discrete dislocation dynamics (DDD) simulations to study dislocation bypass mechanisms in Ni-based superalloys. We develop a new method to compute generalized stacking fault forces in DDD simulations, based on a concept borrowed from complex analysis and known as the winding number of a closed curve about a point. We use this method to study the mechanisms of bypass of a square lattice of spherical $\gamma'$ precipitates by $a/2\langle110\rangle$ {111} edge dislocations, as a function of the precipitates volume fraction and size. We show that not only the hybrid mechanism is possible, but also that it operates as the transition mechanism between the shearing and looping regimes over a wide range of precipitates volume fraction and radii. Based on our simulation results, we propose a simple model for the strength of this mechanism. We also consider the effects of a $\gamma$/$\gamma'$ lattice misfit on the bypass mechanisms, which we approximate by an additional precipitate stress computed according to Eshelby’s inclusion theory. We show that in the shearing and hybrid looping-shearing regimes, a lattice misfit generally results in an increased bypass stress. For sufficiently high lattice misfit, the critical bypass configuration in attractive dislocation-precipitates interactions changes dramatically, and the bypass stress is controlled by the pinning of the trailing dislocation on the exit side of the precipitates, similar to what has been reported in the high-temperature creep literature.
Monolithic integration of lattice-mismatched semiconductor materials opens up access to a wide range of bandgaps and new device functionalities. However, it is inevitably accompanied by defect formation. A thorough analysis of how these defects propagate and interact with interfaces is critical to understanding their effects on device parameters. In this study, we present a comprehensive study of dislocation networks in the GaSb/GaAs heteroepitaxial system using transmission electron microscopy (TEM). Specifically, the sample analyzed is a GaSb film grown on GaAs using dislocation–reduction strategies such as interfacial misfit array formation and introduction of a dislocation filtering layer. Using various TEM techniques, it is shown that such an analysis can reveal important information on the dislocation behavior including filtering mechanism, types of dislocation reactions, and other interactions with interfaces. A novel method that enables plan-view imaging of deeply embedded interfaces using TEM and a demonstration of independent imaging of different dislocation types are also presented. While clearly effective in characterizing dislocation behavior in GaSb/GaAs, we believe that the methods outlined in this article can be extended to study other heteroepitaxial material systems.
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 was 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. Finally, the simulations show that dislocations exhibit some self-organized structures as the strain is increased.
We report the mechanical behavior of materials operating under high temperatures is strongly influenced by creep mechanisms such as dislocation climb, which is controlled by the diffusion of vacancies. However, atomistic simulations of these mechanisms have traditionally been impractical due to the long time scales required. To overcome these time scale challenges, we use Parallel Trajectory Splicing (ParSplice), an accelerated molecular dynamics method, to simulate dislocation climb in nickel. We focus on modeling the activity of a vacancy near a jog on an edge dislocation in order to observe vacancy pipe diffusion and vacancy absorption at the jog. From rigorously constructed trajectories encompassing more than 2000 vacancy absorption events over a simulation time of more than 4μs at 900 K, a comprehensive sampling of available atomistic mechanisms is collated and analyzed further with molecular statics calculations. We estimate average rates for pipe diffusion and vacancy absorption into the jog using data from the dynamic and static calculations, finding very good agreement. Our results strongly suggest that the dominant mechanism for vacancy absorption by jogs is via biased diffusion to the dislocation core followed by fast pipe diffusion to the jog.
Additively manufactured (AM) 316 L stainless steel (SS) after hot isostatic pressing (HIP) was found to exhibit superior resistance to irradiation-assisted stress corrosion cracking (IASCC) in high-temperature water, as compared to wrought 316 L SS. The well-accepted IASCC factors of radiation-induced segregation (RIS) and radiation hardening are not accurate descriptions of IASCC susceptibility in this case. A decreased strain localization along grain boundaries, caused by dislocation channel broadening, was confirmed to suppress crack initiation. A unique distribution of irradiation-induced defects in HIP AM 316 L SS eased dislocation cross-slip compared to those in the wrought counterpart, thus increasing the channel width near the grain boundaries. Here for the first time, this study highlights the importance of dislocation channel broadening as a potential mechanism to further improve the IASCC resistance of 316 L SS and provides direct experimental evidence based on commercial-grade materials.
Refractory superalloys (RSAs) are a new class of high-temperature structural materials consisting of a nanometer-sized mixture of two coherent A2 (disordered BCC) and B2 (ordered) phases. While these RSAs typically exhibit superior mechanical properties at elevated temperatures, they suffer from limited plasticity at lower temperatures, often attributed to the continuous ordered B2 matrix in these two-phase alloys. To knowledgably control the microstructure and mechanical properties of these alloys, the properties of the constitutive phases in RSAs should be known. Currently this information is absent. Here, in the present work, the microstructure and mechanical properties (at 20–1200 °C) of a candidate Al 0.5 NbTa 0.8 Ti 1.5 V 0.2 Zr RSA and two single-phase alloys, which compositions correspond to the compositions of the A2 and B2 phases in the selected RSA, are reported. Interestingly, and quite unexpectedly, the intrinsic mechanical behavior of the constituent A2 and B2 phases, deciphered for the first time in the present study, revealed that they are softer and substantially more ductile than the parent two-phase RSA. Furthermore, the ordered B2 phase exhibits a deformation twinning mechanism leading to twinning induced plasticity (TWIP) coupled with a low APB energy and activity of 1/2<111> dislocations. These mechanisms are thought to be responsible for the high plasticity of the ordered B2 phase. The results presented here question the current assumption that the ordered B2 phase in these RSAs has limited plasticity and strongly indicate that the nature of the B2/A2 interface boundaries and coherency strengthening play an important role in controlling the mechanical properties of the two-phase RSA.
A long standing challenge in computational materials science is to establish a quantitative connection between the macroscopic properties of plastic deformation with the microscopic mechanisms of dislocations in crystalline materials. Although the discrete dislocation dynamics (DDD) simulation method has been developed for several decades with the goal of addressing this challenge, a one-to-one comparison between the DDD predictions on single crystal stress–strain curves and experimental measurements under identical conditions has not been possible to date. Such a comparison is an essential step towards establishing a dislocation-physics based theory of plasticity and a multiscale framework of the plastic behaviors of crystalline materials. Here we provide direct comparisons between the stress–strain curves of Cu single crystals under high strain rate loading in the [0 0 1] and [0 1 1] directions obtained from miniaturized desktop Kolsky bar experiments and those from DDD simulations under identical loading conditions. With an appropriate set of parameters, DDD simulations can produce stress–strain curves that are in reasonable agreement with the experimental results. However, the dislocation mobility values needed to achieve this agreement are an order of magnitude lower than expected based on previous measurements and atomistic simulations. We hypothesize that this discrepancy could be caused by drag forces from jogs and point defects produced during the plastic deformation. Cross-slip of screw dislocations is also found to be necessary to capture the experimental stress–strain behavior, especially for the [0 1 1] loading direction. Finally, this work provides an example of how direct comparisons between DDD simulations and experimental measurements can provide new insight into the fundamental mechanisms of plastic deformation.
Using a nondestructive, in situ near-infrared defect imaging technique, we obtained time-resolved images of dark line defect (DLD) growth in the cavities of several high-power 790-nm diode lasers and calculated individual DLD growth rates. A clear longitudinal spatial dependence in the DLD growth rate was observed, in which the DLD growth rate was faster near the anti-reflective (AR) facet. Longitudinal spatial hole burning simulation results indicate that the total photon density increases toward the AR facet, while the carrier density decreases toward the AR facet. Here, these results suggest that photon absorption plays a key role in enhancing DLD growth beyond the contribution of the traditional carrier capture mechanism. Asymmetric DLD growth that was preferentially directed toward the high-reflectivity facet was also observed, providing further support for the role of photon absorption in the growth process. To account for these findings, we propose a photon-accelerated mechanism for the enhancement of DLD growth in which the slower carrier capture process in the traditional recombination-enhanced dislocation glide mechanism is replaced by a faster photon absorption process.
Deformation at high strain rates often results in high stresses on many engineering materials, potentially leading to catastrophic failure without proper design. High-strain-rate mechanical testing is thus needed to improve the design of future structural materials for a wide range of applications. Although several high-strain-rate mechanical testing techniques have been developed to provide a fundamental understanding of material responses and microstructural evolution under high-strain-rate deformation conditions, these tests are often very time consuming and costly. In this work, we utilize a high-strain-rate nanoindentation testing technique and system in combination with transmission electron microscopy to reveal the deformation mechanisms and dislocation substructures that evolve in pure metals from low (10 –2 s –1 ) to very high indentation strain rates (10 4 s –1 ), using face-centered cubic aluminum and body-centered cubic molybdenum as model materials. The results help to establish the conditions under which micro- and macro-scale tests can be compared with validity and also provide a promising pathway that could lead to accelerated high-strain-rate testing at substantially reduced costs.
The precipitation-strengthened Nickel alloy Inconel® 740H® (IN740H) exhibits increased ductility at higher applied strain rates during quasi-static tensile tests at an elevated temperature of 760°C. The examination of fracture surfaces in this context reveals a noteworthy transition of underlying fracture mechanisms from transgranular to intergranular fracture as the applied strain rate decreases from 1×10 -3 / s to 0.83×10 -4 / s . To thoroughly understand the mechanical response of IN740H under these conditions, this study develops a crystal plasticity finite element (CPFE) model. Further, this model incorporates various deformation mechanisms including dislocation slips, climb, and grain boundary sliding, which are relevant to the test conditions. The model is calibrated using data from both tensile tests at different strain rates and creep tests across a broad stress range at 760°C, enabling the accurate determination of model parameters for each mechanism. Simulation results well captured the experimental observations of different failure modes. At higher strain rates, the model shows a dominance of dislocation slip leading to heterogeneous plastic deformation and formation of transgranular shear bands causing the failure, while at lower strain rates, an increased activity of grain boundary sliding causes grain boundaries crack leading to intergranular failure.
Medium- and high-entropy alloys based on the CrCoNi-system have been shown to display outstanding strength, tensile ductility and fracture toughness (damage-tolerance properties), especially at cryogenic temperatures. Here we examine the J Ic and (back-calculated) K JIc fracture toughness values of the face-centered cubic, equiatomic CrCoNi and CrMnFeCoNi alloys at 20 K. At flow stress values of ~1.5 GPa, crack-initiation K JIc toughnesses were found to be exceptionally high, respectively 235 and 415 MPa√m for CrMnFeCoNi and CrCoNi, with the latter displaying a crack-growth toughness K ss exceeding 540 MPa√m after 2.25 mm of stable cracking, which to our knowledge is the highest such value ever reported. Characterization of the crack-tip regions in CrCoNi by scanning electron and transmission electron microscopy reveal deformation structures at 20 K that are quite distinct from those at higher temperatures and involve heterogeneous nucleation, but restricted growth, of stacking faults and fine nanotwins, together with transformation to the hexagonal closed-packed phase. The coherent interfaces of these features can promote both the arrest and transmission of dislocations to generate respectively strength and ductility which strongly contributes to sustained strain hardening. Indeed, we believe that these nominally single-phase, concentrated solid-solution alloys develop their fracture resistance through a progressive synergy of deformation mechanisms, including dislocation glide, stacking-fault formation, nano-twinning and eventually in situ phase transformation, all of which serve to extend continuous strain hardening which simultaneously elevates strength and ductility (by delaying plastic instability), leading to truly exceptional resistance to fracture.
The ability of a body-centered-cubic metal to deform plastically is limited by the thermally activated glide motion of screw dislocations, which are line defects with a mobility exhibiting complex dependence on temperature, stress, and dislocation segment length. We derive an analytical expression for the velocity of dislocation glide, based on a statistical mechanics argument, and identify an apparent phase transition marked by a crossover temperature above which the activation energy for glide effectively halves, changing from the formation energy of a double kink to that of a single kink. The analysis is in quantitative agreement with direct kinetic Monte Carlo simulations.
Abstract Silicon-based layered nanocomposites, comprised of covalent-metal interfaces, have demonstrated elevated resistance to radiation. The amorphization of the crystalline silicon sublayer during irradiation and/or heating can provide an additional mechanism for accommodating irradiation-induced defects. In this study, we investigated the mechanical strength of irradiated Si-based nanocomposites using atomistic modeling. We first examined dose effects on the defect evolution mechanisms near silicon-gold crystalline and amorphous interfaces. Our simulations reveal the growth of an emergent amorphous interfacial layer with increasing dose, a dominant factor mitigating radiation damage. We then examined the effect of radiation on the mechanical strength of silicon-gold multilayers by constructing yield surfaces. These results demonstrate a rapid onset strength loss with dose. Nearly identical behavior is observed in bulk gold, a phenomenon that can be rooted to the formation of radiation-induced stacking fault tetrahedra which dominate the dislocation emission mechanism during mechanical loading. Taken together, these results advance our understanding of the interaction between radiation-induced point defects and metal-covalent interfaces.
Severe plastic deformation of metals is known to lead to superior properties that cannot be achieved by any traditional metallurgic process. Origin of the superior properties is perceived to be closely associated with grain refinement, a fundamental process during the severe plastic deformation, which is essentially the formation of new grain boundaries. However, the atomistic mechanism of grain boundary formation remains largely obscure. Here, by using in-situ transmission electron microscopy and molecular dynamic simulation, we reveal, for the first time at atomic level, the shear-induced low-angle grain boundary (LAGB) formation processes in Au nanocrystal. In this work, we discover the LAGB formation is accomplished through inward propagation of nanotwins accompanied by dislocations gliding on the twin boundaries, a nanotwin-mediated dislocation slip mechanism, which shows reversible characteristic under reciprocating shear load and is affected by the nanocrystal microstructure and orientation. Our result unveils unprecedented atomistic insights on shear driven grain refinement towards nanostructure of superior properties.
Sputter deposited nickel-molybdenum-tungsten (Ni-Mo-W) thin films possess a beneficial suite of properties that stem from the extremely fine growth twins that form during the deposition process. Previously these materi- als were only characterized in tension, however, in this study in situ micropillar compression and post-mortem microstructural analysis of nanotwinned Ni 84 Mo 11 W 5 micropillars were employed to measure the compressive response and elucidate the attendant deformation mechanisms. The pillars exhibit Hookean behavior to com- pressive strengths of 3-3.5 GPa and the onset of non-linear plastic flow was manifest by discrete strain bursts and highly localized shear bands. Plastic deformation was concentrated at the top of the pillar, while the bulk of the micropillar was nominally unaffected. Post-mortem investigations indicate that at sufficiently high stresses shear banding is triggered, resulting in intense and highly localized plastic deformation that led to the formation of twin-free nanocrystalline grains within highly deformed shear bands. By contrast, the regions adjacent to the shear bands were unaffected. The absence of detwinning and dislocation glide mechanisms was unexpected and in direct contrast to what has been observed in nanotwinned Cu-Al. Post-mortem observations of the Ni-Mo-W micropillars suggest that the ultrafine twins create a unique form of dislocation starvation and source-limited plasticity. The ultrahigh compressive strength is governed by the triggering of shear bands rather than the activa- tion of dislocation glide. Finally, the specialized nature of plasticity in nanotwinned Ni-Mo-W is clear, even though the precise trigger for shear band formation in nanotwinned Ni-Mo-W remains to be identified.