In-situ high-energy X-ray diffraction and crystal plasticity modeling to predict the evolution of te
Explore the source record for details and available documents.
Ferreri, Nicholas C.↗
SEARCH · Search NASA
Results for “anisotropic plasticity”
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.
Cross-slip and easy-glide CRSS in titanium: Theoretical predictions and in-situ TEM measurements
This study investigates the mechanics of prismatic and first-order pyramidal $\langle$a$\rangle$ slip in titanium (Ti), elucidating the physics of easy-glide and cross-slip through a combination of theory and experiments. Screw-character prismatic (Pr) dislocations in Ti are of particular interest because their complex cores can be stable or unstable, leading to activation by either cross-slip or planar glide. To investigate these mechanisms, site-specific micro tensile samples were prepared using focused ion beam (FIB) milling and mounted on a push-to-pull (PTP) device for in-situ transmission electron microscopy (TEM) tensile testing. The in-situ experiments provide direct observations of the onset of dislocation motion and the precise determination of the critical resolved shear stress (CRSS) for the activated mechanisms, and its evolution with load cycling. A comprehensive theory has been developed to predict the CRSS values for easy glide, cross-slip, and multiplication of dislocations. Predicted critical stresses for pyramidal (π)-to-Pr and reverse cross-slip agree closely with the experimental measurements. The latter cross-slip stress is a factor of two higher than that of unobstructed planar slip. The model accounts for overlapping dislocation cores and employs a Wigner-Seitz based cell to evaluate misfit energies. By combining ab initio density functional theory (DFT) with anisotropic elasticity, the framework identifies minimum energy pathways for dislocation glide, which can be intermittent and zig-zag. A simplified expression utilizing (π) and Pr Schmid factor ratios is proposed for critical stress corresponding to (π)-to-Pr cross-slip transition. The results are strongly dependent on crystal orientation, underscoring non-Schmid behavior. Overall, this study explores key critical stress parameters essential for informing higher-scale simulations of plasticity in Ti.
Critical stress↗
Effects of rhodium doping on dislocation nucleation in a [001] SrNi 2 P 2 single crystal under spherical nanoindentation
Nanoindentation was performed on SrNi 2 P 2 single crystals and their dilute solid solutions, Sr(Ni 1–x Rh x ) 2 P 2 (x = 0 ~ 0.055), under [001] loading to investigate the effects of elemental doping on Young’s modulus and dislocation nucleation stress. The results show that Young’s modulus and the dislocation nucleation stress decrease as the Rh content increases, and their Rh content dependence also varies with the Rh content. Electron diffraction analysis and the anisotropic lattice distortion calculations revealed that a local structural transition from the orthorhombic superstructure to the tetragonal structure and local compressive residual stresses could be the main reasons for the decrease in Young’s modulus and dislocation nucleation stress, respectively. The overlap of stress–strain fields between Rh atoms and the associated local structural transition contribute to their Rh content dependence. Here, the results also indicate that elemental doping does not necessarily strengthen materials if the incipient plasticity is controlled by dislocation nucleation, unlike conventional solid solution strengthening.
36 MATERIALS SCIENCE↗
Can changes in deformation regimes be inferred from crystallographic preferred orientations in polar ice?
Creep due to ice flow is generally thought to be the main cause for the formation of crystallographic preferred orientations (CPOs) in polycrystalline anisotropic ice. However, linking the development of CPOs to the ice flow history requires a proper understanding of the ice aggregate's microstructural response to flow transitions. In this contribution the influence of ice deformation history on the CPO development is investigated by means of full-field numerical simulations at the microscale. We simulate the CPO evolution of polycrystalline ice under combinations of two consecutive deformation events up to high strain, using the code VPFFT (visco-plastic fast Fourier transform algorithm) within ELLE. A volume of ice is first deformed under coaxial boundary conditions, which results in a CPO. The sample is then subjected to different boundary conditions (coaxial or non-coaxial) in order to observe how the deformation regime switch impacts the CPO. The model results indicate that the second flow event tends to destroy the first, inherited fabric with a range of transitional fabrics. However, the transition is slow when crystallographic axes are critically oriented with respect to the second imposed regime. Therefore, interpretations of past deformation events from observed CPOs must be carried out with caution, particularly in areas with complex deformation histories.
54 ENVIRONMENTAL SCIENCES↗
Visualization and validation of twin nucleation and early-stage growth in magnesium
The abrupt occurrence of twinning when Mg is deformed leads to a highly anisotropic response, making it too unreliable for structural use and too unpredictable for observation. Here, we describe an in-situ transmission electron microscopy experiment on Mg crystals with strategically designed geometries for visualization of a long-proposed but unverified twinning mechanism. Combining with atomistic simulations and topological analysis, we conclude that twin nucleation occurs through a pure-shuffle mechanism that requires prismatic-basal transformations. Also, we verified a crystal geometry dependent twin growth mechanism, that is the early-stage growth associated with instability of plasticity flow, which can be dominated either by slower movement of prismatic-basal boundary steps, or by faster glide-shuffle along the twinning plane. The fundamental understanding of twinning provides a pathway to understand deformation from a scientific standpoint and the microstructure design principles to engineer metals with enhanced behavior from a technological standpoint.
36 MATERIALS SCIENCE↗
A general solution for accelerating screw dislocations in arbitrary slip systems with reflection symmetry
Solutions to the differential equations of linear elasticity in the continuum limit in arbitrary crystal symmetry are known only for steady-state dislocations of arbitrary character, i.e. line defects moving at constant velocity. Troubled by singularities at certain ‘critical’ velocities (typically close to certain sound speeds), these dislocation fields are thought to be too idealized, and divergences are usually attributed to neglecting the finite size of the core and to the restriction to constant velocity. In the isotropic limit, accelerating pure screw and edge dislocations were studied some time ago. A generalization to anisotropic crystals has been attempted for pure screw and edge dislocations only for some special cases. This work aims to fill the gap of deriving a general anisotropic solution for pure screw dislocations applicable to slip systems featuring a reflection symmetry, a prerequisite to studying pure screw dislocations without mixing with edge dislocations. Finally, further generalizations to arbitrary mixed dislocations as well as regularizations of the dislocation core are beyond the scope of this paper and are left for future work.
42 ENGINEERING↗
A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism
We present a macroscale constitutive model that couples magnetism with thermal, elastic, plastic, and damage effects in an Internal State Variable (ISV) theory. Previous constitutive models did not include an interdependence between the internal magnetic (magnetostriction and magnetic flux) and mechanical fields. Although constitutive models explaining the mechanisms behind mechanical deformations caused by magnetization changes have been presented in the literature, they mainly focus on nanoscale structure–property relations. A fully coupled multiphysics macroscale ISV model presented herein admits lower length scale information from the nanoscale and microscale descriptions of the multiphysics behavior, thus capturing the effects of magnetic field forces with isotropic and anisotropic magnetization terms and moments under thermomechanical deformations. For the first time, this ISV modeling framework internally coheres to the kinematic, thermodynamic, and kinetic relationships of deformation using the evolving ISV histories. For the kinematics, a multiplicative decomposition of deformation gradient is employed including a magnetization term; hence, the Jacobian represents the conservation of mass and conservation of momentum including magnetism. The first and second laws of thermodynamics are used to constrain the appropriate constitutive relations through the Clausius–Duhem inequality. The kinetic framework employs a stress–strain relationship with a flow rule that couples the thermal, mechanical, and magnetic terms. Experimental data from the literature for three different materials (iron, nickel, and cobalt) are used to compare with the model’s results showing good correlations.
36 MATERIALS SCIENCE↗
A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism
We present a macroscale constitutive model that couples magnetism with thermal, elastic, plastic, and damage effects in an Internal State Variable (ISV) theory. Previous constitutive models did not include an interdependence between the internal magnetic (magnetostriction and magnetic flux) and mechanical fields. Although constitutive models explaining the mechanisms behind mechanical deformations caused by magnetization changes have been presented in the literature, they mainly focus on nanoscale structure–property relations. A fully coupled multiphysics macroscale ISV model presented herein admits lower length scale information from the nanoscale and microscale descriptions of the multiphysics behavior, thus capturing the effects of magnetic field forces with isotropic and anisotropic magnetization terms and moments under thermomechanical deformations. For the first time, this ISV modeling framework internally coheres to the kinematic, thermodynamic, and kinetic relationships of deformation using the evolving ISV histories. For the kinematics, a multiplicative decomposition of deformation gradient is employed including a magnetization term; hence, the Jacobian represents the conservation of mass and conservation of momentum including magnetism. The first and second laws of thermodynamics are used to constrain the appropriate constitutive relations through the Clausius–Duhem inequality. The kinetic framework employs a stress–strain relationship with a flow rule that couples the thermal, mechanical, and magnetic terms. Experimental data from the literature for three different materials (iron, nickel, and cobalt) are used to compare with the model’s results showing good correlations.
36 MATERIALS SCIENCE↗
Dislocation drag and its influence on elastic precursor decay
Plastic deformation is mediated by the creation and movement of dislocations, and at high stress the latter is dominated by dislocation drag from phonon wind. Here, by simulating a 1-D shock impact problem we analyze the importance of accurately modeling dislocation drag and dislocation density evolution in the high stress regime. Dislocation drag is modeled according to a first-principles derivation as a function of stress and dislocation character, and its temperature and density dependence are approximated to the extent currently known. Much less is known about dislocation density evolution, leading to far greater uncertainty in these model parameters. In studying anisotropic fcc metals with character dependent dislocations, the present work generalizes similar earlier studies by other authors.
36 MATERIALS SCIENCE↗
A new woven composite constitutive model validated by shock wave experiments
In this paper, we present results of plate impact simulations of shock compressed woven glass fiber-reinforced plastic (GRP) performed using the Arbitrary Lagrangian–Eulerian three-dimensional finite element code. A hyperelastic large-strain-based empirical Continuum Damage Mechanics (CDM) formulation is employed to describe damage initiation and growth in the shock-compressed GRP. The model parameters calibration scheme utilizes the Velocity Interferometer System for Any Reflector normal particle velocity measurements at the free surface of the GRP target plates. The impact velocity in the experiments ranged from 8.5 to 418 m/s. The finite element model considered planar 0°/90° bidirectional plies with an individual ply thickness of 0.68 mm, stacked to reach a total laminate thickness of 6.8 mm. The anisotropic elastic strains were estimated from the experimentally determined tetragonal symmetry stiffness matrix for the GRP. The strain-based damage model captures several salient features observed in the measured free surface particle wave profiles, including the shock rise time, onset of Elastic—Elastic Cracking, and the shape of the nonlinear portion of the experimental particle velocity profiles. Furthermore, the CDM model predicts the dominant damage mode to be matrix microcracking due to shear and the associated bulk expansion (bulking) under the global compressive loading in the plate impact configuration.
42 ENGINEERING↗
Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime
Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.
36 MATERIALS SCIENCE↗
Mechanics of the Ring Tension Test (RTT): A Finite Element-based Investigation
The ring tension test (RTT) is an experimental method for determining mechanical behavior in a material’s circumferential or hoop direction. It is a crucial test for testing anisotropic materials with tube geometry, such as nuclear fuel cladding or irradiated pipes. Several RTT configurations exist, each with their own advantages and disadvantages. However, this test is significantly more complex than traditional tensile testing and can be especially sensitive to small differences and inconsistencies in the test setup and geometry, ultimately affecting the derived mechanical properties. Previous research has focused on method development, and little work has been done on understanding the subtle differences between an ideal test and experiments, specifically when the tests are performed on highly irradiated materials in hot cells. In this work, a finite element-based investigation of the RTT is conducted. Two promising test configurations are investigated, comparing their ability to determine accurate material strengths through plastic deformation. Several non-ideal conditions and uncontrollable effects which are likely to occur during experimental testing such as machining tolerances, variations of specimen geometry from nominal dimensions, rotation of specimens and fixturing, and other test setup discrepancies are studied. The sensitivity of measured strengths to these conditions is presented. A mechanics-based approach to describing and correcting raw data to determine actual strengths is also included for one of the configurations, resulting in a robust correction method with highly accurate material strength measurements. Based on these analyses, the hemicylindrical mandrel configuration is recommended with a gauge region oriented at a 45° angle.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Room Temperature Quasi-static Characterization and Constitutive Model Parametrization of Flexible Polyurethane Foams of Different Densities Loaded in Different Orientations
This report describes the efforts to characterize and model General Plastics TF6070 and EF4000 flexible polyurethane foams under room temperature, large deformation quasi-static cyclic mechanical loading conditions. Densities from three to fifteen pound per cubic foot (PCF) are examined, which correspond to relative densities of approximately 4 to 20%. These foams are open cell structured and flexible at room temperature with a glass transition transition less than -30°C, and they fully recover their original shape when unloaded (at room temperature). Uniaxial compression tests were conducted with accompanying lateral image series for Digital Image Correlation (DIC) analysis with the goal of extracting transverse strain responses. Due to difficulties with DIC analysis at large strains, lateral strains were instead extracted for each test via edge tracking. The experimental results exhibit a nonlinear elastic response and anisotropic material behavior (particularly for the lower densities). Some hysteresis is observed that is different between the first and subsequent cycles of deformation indicating both a small degree of permanent damage (reduced stiffness during reloading) and viscoelasticity. These inelastic mechanisms are not considered in the modeling and calibration in this report. This work considers only the rate independent, room temperature foam behavior. Individual foam densities were calibrated for loading in two directions, parallel and perpendicular to the foam bubble rise direction, since the mechanical behavior is different in these two directions. The Flex Foam constitutive model was used for all parameterizations despite the fact that the model is isotropic. A review of the constitutive model is given as well as necessary data reduction procedures to parameterize it for each foam density and orientation are discussed. Finally, two different parameterizations are developed that take the undeformed foam density as an input that span all densities considered. These two parameterized models represent foams loaded in the rise and transverse directions respectively. We summarize the assumptions and limitations of the parameterizations provided in this report to guide analysis choices with them. All parameterizations presented herein have the following traits, room temperature, rate independent, damage-free, and non-dissipative . Isotropy (even if they are representing anisotropic data). Supplied Sierra Solid Mechanics Flex Foam Model Inputs are in units: pounds, inches, Celsius, and seconds
36 MATERIALS SCIENCE↗
Vortex phases and glassy dynamics in the highly anisotropic superconductor HgBa 2 CuO 4+δ
We present an extensive study of vortex dynamics in a high-quality single crystal of HgBa 2 CuO 4+ δ , a highly anisotropic superconductor that is a model system for studying the effects of anisotropy. From magnetization M measurements over a wide range of temperatures T and fields H , we construct a detailed vortex phase diagram. We find that the temperature-dependent vortex penetration field H p ( T ), second magnetization peak H smp ( T ), and irreversibility field H irr ( T ) all decay exponentially at low temperatures and exhibit an abrupt change in behavior at high temperatures T / T c >~0.5. By measuring the rates of thermally activated vortex motion (creep) S ( T , H ) = | d ln M ( T , H )/ d ln t |, we reveal glassy behavior involving collective creep of bundles of 2D pancake vortices as well as temperature- and time-tuned crossovers from elastic (collective) dynamics to plastic flow. Based on the creep results, we show that the second magnetization peak coincides with the elastic-to-plastic crossover at low T , yet the mechanism changes at higher temperatures.
36 MATERIALS SCIENCE↗
Evaluation of the ring tension test (RTT) for robust determination of material strengths
The ring tension test (RTT) is an experimental method for determining mechanical behavior in a material’s circumferential or hoop direction. It is a crucial test for testing anisotropic materials with tube geometry, such as nuclear fuel cladding or irradiated pipes. Several RTT configurations exist, each with their own advantages and disadvantages. However, this test is significantly more complex than traditional tensile testing and can be especially sensitive to small differences and inconsistencies in the test setup and geometry, ultimately affecting the derived mechanical properties. Previous research has focused on method development, and little work has been done on understanding the subtle differences between an ideal test and experiments, specifically when the tests are performed on highly irradiated materials in hot cells. In this work, a finite element-based investigation of the RTT is conducted to fill this gap. The two most used test configurations are investigated, comparing their ability to determine accurate material strengths through plastic deformation. Several non-ideal conditions and uncontrollable effects which are likely to occur during experimental testing such as machining tolerances, variations of specimen geometry from nominal dimensions, rotation of specimens and fixturing, and other test setup discrepancies are studied. The sensitivity of measured strengths to these conditions is presented. A mechanics-based approach to describing and correcting raw data to determine actual strengths is also included for one of the configurations, resulting in a robust correction method with highly accurate material strength measurements. Finally, based on these analyses, the hemicylindrical mandrel configuration is recommended with a gauge region oriented at a 45° angle.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)
The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.
36 MATERIALS SCIENCE↗
A publicly available PyTorch-$\mathrm{ABAQUS}$ $\mathrm{UMAT}$ deep-learning framework for level-set plasticity
Here this paper introduces a publicly available PyTorch-ABAQUS deep-learning framework of a family of plasticity models where the yield surface is implicitly represented by a scalar-valued function. In particular, our focus is to introduce a practical framework that can be deployed for engineering analysis that employs a user-defined material subroutine (UMAT/VUMAT) for ABAQUS, which is written in FORTRAN. To accomplish this task while leveraging the back-propagation learning algorithm to speed up the neural-network training, we introduce an interface code where the weights and biases of the trained neural networks obtained via the PyTorch library can be automatically converted into a generic FORTRAN code that can be a part of the UMAT/VUMAT algorithm. To enable third-party validation, we purposely make all the data sets, source code used to train the neural-network-based constitutive models, and the trained models available in a public repository. Furthermore, the practicality of the workflow is then further tested on a dataset for anisotropic yield function to showcase the extensibility of the proposed framework. A number of representative numerical experiments are used to examine the accuracy, robustness and reproducibility of the results generated by the neural network models.
36 MATERIALS SCIENCE↗
Computationally efficient method for determining limiting velocities of edge dislocations in anisotropic crystals
The continuum-limit theory of dislocations in crystals predicts divergences in the elastic energy at crystal-geometry dependent limiting velocities vL, which separate subsonic, transsonic, and supersonic dislocation glide regimes and are therefore import for material strength models at high strain rates. Although it is known how to calculate those limiting velocities, there is one special case - edge dislocations with reflection symmetry, but non-vanishing elastic constants c16 or c26 - where previous methods have been notoriously slow. In this letter, we address this deficiency by deriving a computationally efficient method for determining the limiting velocities of edge dislocations with reflection symmetry which is two orders of magnitude faster than the previous method.
36 MATERIALS SCIENCE↗