Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic strain engineering”

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.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Spin-mediated shear oscillators in a van der Waals antiferromagnet

Understanding how microscopic spin configuration gives rise to exotic properties at the macroscopic length scale has long been pursued in magnetic materials. One seminal example is the Einstein-de Haas effect in ferromagnets, in which angular momentum of spins can be converted into mechanical rotation of an entire object. However, for antiferromagnets without net magnetic moment, how spin ordering couples to macroscopic movement remains elusive. Here we observed a seesaw-like rotation of reciprocal lattice peaks of an antiferromagnetic nanolayer film, whose gigahertz structural resonance exhibits more than an order-of-magnitude amplification after cooling below the Neel temperature. Using a suite of ultrafast diffraction and microscopy techniques, we directly visualize this spin-driven rotation in reciprocal space at the nanoscale. This motion corresponds to interlayer shear in real space, in which individual micro-patches of the film behave as coherent oscillators that are phase-locked and shear along the same in-plane axis. Using time-resolved optical polarimetry, we further show that the enhanced mechanical response strongly correlates with ultrafast demagnetization, which releases elastic energy stored in local strain gradients to drive the oscillators. Finally, our work not only offers the first microscopic view of spin-mediated mechanical motion of an antiferromagnet but it also identifies a new route towards realizing high-frequency resonators up to the millimetre band, so the capability of controlling magnetic states on the ultrafast timescale can be readily transferred to engineering the mechanical properties of nanodevices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Realizing mechanical frustration at the nanoscale using DNA origami

Structural designs inspired by physical and biological systems have been previously utilized to develop mechanical metamaterials with enhanced properties based on clever geometric arrangement of constituent building blocks. Here, we use the DNA origami method to realize a nanoscale metastructure exhibiting mechanical frustration, a counterpart of the well-known phenomenon of magnetic frustration. By selectively actuating reconfigurable struts, it adopts either frustrated or non-frustrated states, each characterized by distinct free energy profiles. While the non-frustrated state distributes the strain homogeneously, the frustrated mode concentrates it at a specific location. Molecular dynamics simulations reconcile the contrasting behaviors and provide insights into underlying mechanics. We explore the design space further by tailoring responses through structural modifications. Our work combines programmable DNA self-assembly with mechanical design principles to overcome engineering limitations encountered at the macroscale to design dynamic, deformable nanostructures with potential applications in elastic energy storage, nanomechanical computation, and allosteric mechanisms in DNA-based nanomachinery.

DNA nanostructures↗
DOE OSTI2024

Influence of chemical strains on the electrocaloric response, polarization morphology, tetragonality, and negative-capacitance effect of ferroelectric core-shell nanorods and nanowires

Using Landau-Ginzburg-Devonshire (LGD) approach, we proposed the analytical description of the influence of chemical strains on spontaneous polarization and the electrocaloric response in ferroelectric core-shell nanorods. We postulate that the nanorod core presents a defect-free single-crystalline ferroelectric material, and elastic defects are accumulated in the ultrathin shell, where they can induce tensile or compressive chemical strains. Finite-element modeling (FEM) based on the LGD approach reveals transitions of domain-structure morphology induced by chemical strains in the Ba Ti O 3 nanorods. Namely, tensile chemical strains induce and support the single-domain state in the central part of the nanorod, while the curled domain structures appear near the unscreened or partially screened ends of the rod. The vortexlike domains propagate toward the central part of the rod and fill it entirely, when the rod is covered by a shell with compressive chemical strains above some critical value. The critical value depends on the nanorod sizes, aspect ratio, and screening conditions at its ends. Both analytical theory and FEM predict that the tensile chemical strains in the shell increase the nanorod polarization, lattice tetragonality, and electrocaloric response well above the values corresponding to the bulk material. The physical reason for the increase is strong electrostriction coupling between the mismatch-type elastic strains induced in the core by chemical strains in the shell. Comparison with earlier XRD data confirmed an increase of the tetragonality ratio in tensile Ba Ti O 3 nanorods compared to the bulk material. Obtained analytical expressions, which are suitable for the description of strain-induced changes in a wide range of multiaxial ferroelectric core-shell nanorods and nanowires, can be useful for strain engineering of advanced ferroelectric nanomaterials for energy storage, harvesting, electrocaloric applications, and negative capacitance elements. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2023

Prediction of nonlinear specific heat during single crystal HMX phase transition

Here, we develop a thermodynamically consistent chemo-thermo-mechanical model for the β → δ phase transition of energetic HMX crystals. In contrast to previous models, which either considered specific heat to be a constant or utilized a calibrated function, this model provides novel expressions for the specific heats at constant volume and constant elastic strains derived directly from continuum mechanics. In addition, the model provides a novel prediction for the critical temperature at which the chemical heating rate achieves its extremum for Arrhenius kinetics. The numerical solution predicts highly nonlinear specific heat behavior including order of magnitude changes.

42 ENGINEERING↗
DOE OSTI2021

Estimation of the orientation of stress in the Earth’s crust without earthquake or borehole data

Abstract Mechanical stress acting in the Earth’s crust is a fundamental property that is important for a wide range of scientific and engineering applications. The orientation of maximum horizontal compressive stress can be estimated by inverting earthquake source mechanisms and measured directly from borehole-based measurements, but large regions of the continents have few or no observations. Here we present an approach to determine the orientation of maximum horizontal compressive stress by measuring stress-induced anisotropy of nonlinear susceptibility, which is the derivative of elastic modulus with respect to strain. Laboratory and Earth experiments show that nonlinear susceptibility is azimuthally dependent in an anisotropic stress field and is maximum in the orientation of maximum horizontal compressive stress. We observe this behavior in the Earth—in Oklahoma and New Mexico, U.S.A, where maximum nonlinear susceptibility coincides with the orientation of maximum horizontal compressive stress measured using traditional methods. Our measurements use empirical Green’s functions and solid-earth tides and can be applied at different temporal and spatial scales.

58 GEOSCIENCES↗
DOE OSTI2022

High Velocity, Low‐Voltage Collective In‐Plane Switching in (100) BaTiO 3 Thin Films

Abstract Ferroelectrics are being increasingly called upon for electronic devices in extreme environments. Device performance and energy efficiency is highly correlated to clock frequency, operational voltage, and resistive loss. To increase performance it is common to engineer ferroelectric domain structure with highly‐correlated electrical and elastic coupling that elicit fast and efficient collective switching. Designing domain structures with advantageous properties is difficult because the mechanisms involved in collective switching are poorly understood and difficult to investigate. Collective switching is a hierarchical process where the nano‐ and mesoscale responses control the macroscopic properties. Using chemical solution synthesis, epitaxially nearly‐relaxed (100) BaTiO 3 films are synthesized. Thermal strain induces a strongly‐correlated domain structure with alternating domains of polarization along the [010] and [001] in‐plane axes and 90° domain walls along the [011] or [01] directions. Simultaneous capacitance–voltage measurements and band‐excitation piezoresponse force microscopy revealed strong collective switching behavior. Using a deep convolutional autoencoder, hierarchical switching is automatically tracked and the switching pathway is identified. The collective switching velocities are calculated to be ≈500 cm s −1 at 5 V (7 kV cm −1 ), orders‐of‐magnitude faster than expected. These combinations of properties are promising for high‐speed tunable dielectrics and low‐voltage ferroelectric memories and logic.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Control of the Folding Dynamics of Self-Reconfiguring Magnetic Microbots Using Liquid Crystallinity

Reconfigurable microdevices are being explored in a range of contexts where their life-like abilities to move and change shape are important. While much work has been done to control the motion of these microdevices by engineering their geometry and composition, little is known about their dynamics in complex fluid environments with non-Newtonian rheology. Here in this paper, it is shown how the actuation dynamics of reconfigurable microdevices made by assembly of patchy magnetic microcubes, which are referred to as “microbots,” can be modulated by their interactions with the anisotropic viscoelastic environment of a liquid crystal (LC). The free energy arising from the elastic strain of LC and formation of topological defects around the microbots influences their folding dynamics, which can be tuned by tailoring both the far-field orientation of the LC and the local ordering of the LC at the microbot surfaces. These findings represent a first step toward establishing a general set of design rules to control the dynamics of microbots using complex anisotropic fluids.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Preliminary Experimental Results in Support of the Development of EPP and SMT Design Methods and Viscoplastic Model for A709

The ASME code qualification effort for Alloy 709 (A709) is currently underway to qualify it for Class A construction in the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code, Section III, Division 5. The United States (U.S.) Department of Energy (DOE) national laboratories are collaborating in the advanced materials development initiative to investigate the mechanical performance of A709 in support of its code qualification. As part of the A709 code qualification effort, this report summarizes ORNL’s initial experimental findings that support the integration of A709 into the elastic-perfectly plastic (EPP) strain limits code case. It also covers thermomechanical fatigue testing conducted to develop viscoelastic material models, along with the preliminary results of creep-fatigue experiments at 816°C using the Simplified Model Test (SMT) method.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Report on FY 2021 testing in support of integrating Alloy 800H and 2.25Cr-1Mo for EPP analysis

The use of simplified Elastic-Perfectly Plastic (EPP) design methods avoids the complexities and limitations of the elastic analysis approach in evaluating strain limits and creep-fatigue damage for Class A metallic coolant boundary components in ASME Boiler and Pressure Vessel Code, Section III, Division 5, Subsection HB, Subpart B. No stress classification and linearization are required in the EPP methods. Further, the strain limits and creep-fatigue damage evaluation procedures based on elastic analysis results were developed for the reactor operating conditions where uncoupled, rate-independent plasticity and stationary creep can be used to approximate the deformation in the component adequately. But at very high temperatures, the deformation is rate-dependent, and it can only be approximated adequately by unified viscoplastic model. Thus, for example, the elastic analysis approach cannot be used to evaluate strain limits and creep-fatigue damage for Alloy 617 above 650 °C. The EPP methods remove the temperature restriction. They can be used to evaluate strain limits and creep-fatigue damage for the full range of temperatures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Prediction of carbon nanostructure mechanical properties and the role of defects using machine learning

Graphene-based nanostructures hold immense potential as strong and lightweight materials, however, their mechanical properties such as modulus and strength are difficult to fully exploit due to challenges in atomic-scale engineering. This study presents a database of over 2,000 pristine and defective nanoscale CNT bundles and other graphitic assemblies, inspired by microscopy, with associated stress–strain curves from reactive molecular dynamics (MD) simulations using the reactive INTERFACE force field (IFF-R). These 3D structures, containing up to 80,000 atoms, enable detailed analyses of structure-stiffness-failure relationships. By leveraging the database and physics- and chemistry-informed machine learning (ML), accurate predictions of elastic moduli and tensile strength are demonstrated at speeds 1,000 to 10,000 times faster than efficient MD simulations. Hierarchical Graph Neural Networks with Spatial Information (HS-GNNs) are introduced, which integrate chemistry knowledge. HS-GNNs as well as extreme gradient boosted trees (XGBoost) achieve forecasts of mechanical properties of arbitrary carbon nanostructures with only 3 to 6% mean relative error. The reliability equals experimental accuracy and is up to 20 times higher than other ML methods. Predictions maintain 8 to 18% accuracy for large CNT bundles, CNT junctions, and carbon fiber cross-sections outside the training distribution. The physics- and chemistry-informed HS-GNN works remarkably well for data outside the training range while XGBoost works well with limited training data inside the training range. The carbon nanostructure database is designed for integration with multimodal experimental and simulation data, scalable beyond 100 nm size, and extendable to chemically similar compounds and broader property ranges. The ML approaches have potential for applications in structural materials, nanoelectronics, and carbon-based catalysts.

Winetrout, Jordan J.↗
DOE OSTI2024

Advances in Molecular Beam Epitaxy Growth of Ultra-Wide Bandgap Ga2O3 Based Alloys

Gallium oxide (Ga2O3) is an emerging ultra-wide bandgap semiconductor material that has attracted attention for its potential to outperform existing SiC and GaN based devices operating at high breakdown voltages and high temperature. Isovalent alloying of In and Al in Ga2O3 provides the ability to engineer bandgap energy and strain of the material. Alloying with Al increases the bandgap energy and the theoretically achievable Baliga's figure of merit, a key measure of a material's ultimate performance limits for high power switching devices. Alloying with In introduces compressive strain and can be used to counteract the tensile strain of Al incorporation. The resulting (AlxGa1-x-yIny)2O3 alloy can be lattice-matched to commercially available Ga2O3 wafers and has a tunable bandgap energy greater than that of Ga2O3, 4.76 eV. Such lattice-matched material can be grown arbitrarily thick without the detrimental effects of elastic strain and relaxation, making it suitable for high voltage diodes and transistors. However, efforts to synthesize isovalent alloys are complicated by their tendency to phase separate into corundum Al2O3 or bixbyite In2O3. Literature reports of the quaternary (AlxGa1-x-yIny)2O3 are limited to <1% unintentional indium incorporation in In-catalyzed (AlxGa1-x)2O3. The primary limitation to quaternary growth is the limited incorporation of indium at elevated growth temperatures. This limited incorporation is due to both the volatility of indium oxide and Al and Ga cation exchange reactions which replace indium in In2O3. We report on the development of a novel high-throughput molecular beam epitaxy (MBE) technique to screen the growth conditions for the ternary alloy (InyGa1-y)2O3, and the application of these findings to the first successful synthesis of phase pure monoclinic (AlxGa1-x-yIny)2O3 by MBE. By leveraging the unique sub-oxide chemistry of Ga2O3 and in-situ monitoring of crystal properties by reflection high-energy electron diffraction (RHEED), a cyclical growth and etch-back method is developed and applied to rapidly characterize the (InyGa1-y)2O3 growth space. This cyclical method provides approximately 10x increase in experimental throughput and up to 46x improvement in Ga2O3 substrate utilization. Appropriate growth conditions for monoclinic (InyGa1- y)2O3 are identified by machine learning analysis of RHEED patterns and targeted growths are characterized ex-situ to confirm improved In incorporation. These growth conditions are then combined with established (AlxGa1-x)2O3 growth conditions to grow quaternary (AlxGa1-x-yIny)2O3 with Al mole fractions ranging from 1.4% - 24.4% and In mole fractions ranging from 3.1% to 15.5%. The chemical and optical properties of the alloys are investigated by XRD, XPS, and spectroscopic ellipsometry. A lattice-matched (AlxGa1-x-yIny)2O3 alloy is examined by 4D-STEM and the chemical and physical uniformity of Al and In incorporation are discussed.

alloy↗
DOE OSTI2023

Toughening Brittle Bio‐P3HB with Synthetic P3HB of Engineered Stereomicrostructures

Abstract Poly(3‐hydroxybutyrate) (P3HB), a biologically produced, biodegradable natural polyester, exhibits excellent thermal and barrier properties but suffers from mechanical brittleness, largely limiting its applications. Here we report a mono‐material product design strategy to toughen stereoperfect, brittle bio or synthetic P3HB by blending it with stereomicrostructurally engineered P3HB. Through tacticity ([ mm ] from 0 to 100 %) and molecular weight ( M n to 788 kDa) tuning, high‐performance synthetic P3HB materials with tensile strength to ≈30 MPa, fracture strain to ≈800 %, and toughness to 126 MJ m −3 (>110× tougher than bio‐P3HB) have been produced. Physical blending of the brittle P3HB with such P3HB in 10 to 90 wt % dramatically enhances its ductility from ≈5 % to 95–450 % and optical clarity from 19 % to 85 % visible light transmittance while maintaining desirably high elastic modulus (>1 GPa), tensile strength (>35 MPa), and melting temperature (160–170 °C). This P3HB‐toughening‐P3HB methodology departs from the traditional approach of incorporating chemically distinct components to toughen P3HB, which hinders chemical or mechanical recycling, highlighting the potential of the mono‐material product design solely based on biodegradable P3HB to deliver P3HB materials with diverse performance properties.

Zhang, Zhen↗
DOE OSTI2023

Toughening Brittle Bio‐P3HB with Synthetic P3HB of Engineered Stereomicrostructures

Abstract Poly(3‐hydroxybutyrate) (P3HB), a biologically produced, biodegradable natural polyester, exhibits excellent thermal and barrier properties but suffers from mechanical brittleness, largely limiting its applications. Here we report a mono‐material product design strategy to toughen stereoperfect, brittle bio or synthetic P3HB by blending it with stereomicrostructurally engineered P3HB. Through tacticity ([ mm ] from 0 to 100 %) and molecular weight ( M n to 788 kDa) tuning, high‐performance synthetic P3HB materials with tensile strength to ≈30 MPa, fracture strain to ≈800 %, and toughness to 126 MJ m −3 (>110× tougher than bio‐P3HB) have been produced. Physical blending of the brittle P3HB with such P3HB in 10 to 90 wt % dramatically enhances its ductility from ≈5 % to 95–450 % and optical clarity from 19 % to 85 % visible light transmittance while maintaining desirably high elastic modulus (>1 GPa), tensile strength (>35 MPa), and melting temperature (160–170 °C). This P3HB‐toughening‐P3HB methodology departs from the traditional approach of incorporating chemically distinct components to toughen P3HB, which hinders chemical or mechanical recycling, highlighting the potential of the mono‐material product design solely based on biodegradable P3HB to deliver P3HB materials with diverse performance properties.

Chemistry↗
DOE OSTI2023

NanoMechanics: Elasticity and Friction in Nano‐Objects

A large scientific and technological effort is underway to investigate the properties of two‐dimensional (2D) materials to become building blocks in integrated nano‐electronic and photonic circuits, composites, coatings, energy harvesting nano‐ systems, nano‐sensors, and nano‐electro‐mechanical systems (NEMS).While several experiments and calculations have revealed exciting novel phenomena in these nanostructures, many scientific and technological questions remain open. In particular, understanding and controlling the structure and mechanical properties of 2D materials at the interface with a solid surface is of key importance to enable the aforementioned nanotechnologies, as well as to fulfill the potential of strain engineering to tailor their properties. The vision of this DoE research program was to investigate novel mechanical properties and phenomena in 1D and 2D materials with the overarching goal of defining a new basic understanding of mechanical behavior in nano and quantum systems. The group of the PI has developed in the last years several studies on the mechanical properties of Carbon nanotubes and oxide nanobelts, more recently the PI has focused her attention on the properties of two‐dimensional materials, such as graphene and MoS 2 , which are a few‐atomic‐layer thick films and hold a great potential for technological and energy applications. The most studied 2D material is graphene, existing as a single layer of graphite or a few‐layer thick epitaxial graphene film. Graphene possesses a large in‐plane Young’s modulus as well as high intrinsic carrier mobility, and high in‐plane thermal conductivity. Besides graphene, also 2D films of graphene oxide (GO), hexagonal Boron Nitride (h‐BN), and transition metal dichalcogenides such as MoS 2 exhibit unique and excellent properties and hold great promise for nanotechnology applications. This research program was aimed at developing a new basic understanding of the mechanical properties of 1D and 2D materials at the interface with solid substrates, with a focus on the role of defects, materials structure, and substrate interaction. Ultimately, we aimed at developing the basic knowledge and experimental tools for manipulating 1D and 2D materials’ structure, stacking, substrate interaction, defects, and number of layers to realize an entirely new class of ultra‐hard, ultra‐thin, and ultra‐light active materials with ad‐ hoc mechanical properties for a variety of applications.

2D Materials↗
DOE OSTI2025

Calibrating a finite-strain phase-field model of fracture for bonded granular materials with uncertainty quantification

To study the mechanical behavior of mock high explosives, an experimental and simulation program was developed to calibrate, with quantified uncertainty, a material model of the bonded granular material Idoxuridine and nitroplasticized Estane-5703. This paper reports on the efficacy of such a framework as a generalizable methodology for calibrating material models against experimental data with uncertainty quantification. Additionally, this paper studies the effect of two manufacturing temperatures and three initial granular configurations on the unconfined compressive behavior of the resulting bonded granular materials. In each of these cases, the same calibration framework was used; in that, hundreds of high-fidelity direct numerical simulations using a new, graphics processing unit-enabled, high-performance finite element method software, Ratel, were run to calibrate a finite-strain phase-field fracture model against experimental data. It was found that manufacturing temperature influenced the elastic response of the mock high explosives, with higher temperatures yielding a stiffer response. By contrast, it was found that the initial configuration of the grains had a negligible impact on the overall behavior of the mock high explosives though it remains possible that local damage accumulation within the specimens could be altered by the initial configurations. Overall, the calibration framework was successful at creating well-calibrated models, showing its usefulness as an engineering and scientific tool.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Achieving ultrahigh modulus of resilience and enhanced thermal stability in ZnO x /SU-8 interpenetrating network polymer nanocomposite nanopillars

The modulus of resilience, a mechanical property that quantifies the maximum strain energy density a material can store during elastic deformation, is a crucial parameter for materials used in flexible displays, micro/nano-electro-mechanical system (M/NEMS) actuators, and ultra-sensitive pressure sensors. In this study, ZnO x /SU-8 nanocomposite nanopillars with a diameter of 300 nm, fully infiltrated with a uniformly distributed, interpenetrating amorphous ZnO x filler network, were synthesized via vapor-phase infiltration (VPI). In-situ uniaxial nano-compression tests revealed that the modulus of resilience of ZnO x /SU-8 reaches ∼ 12 MJ/m 3 , which is an ultrahigh value among all engineering materials with comparable strength. In addition, the synthesis fidelity, inorganic infiltration depth, and mechanical performance were all significantly improved compared to VPI-synthesized AlO x nanocomposites. Thermal stability, another key requirement for M/NEMS device materials operating under extreme environments, was also notably enhanced. Furthermore, partial crystallization of the amorphous ZnO x fillers during annealing contributed to an additional increase in modulus of resilience, reaching up to ∼ 13.9 MJ/m 3 . This work presents an effective fabrication strategy for producing nanostructured organic–inorganic hybrid nanocomposites with ultrahigh modulus of resilience and superior thermal stability, paving the way for their integration into next-generation flexible displays and high-performance M/NEMS devices working under harsh environments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A materials data framework and dataset for elastomeric foam impact mitigating materials

The availability of materials data for impact-mitigating materials has lagged behind applications-based data. For example, data describing on-field helmeted impacts are available, whereas material behaviors for the constituent impact-mitigating materials used in helmet designs lack open datasets. Here, we describe a new FAIR (findable, accessible, interoperable, reusable) data framework with structural and mechanical response data for one example elastic impact protection foam. The continuum-scale behavior of foams emerges from the interplay of polymer properties, internal gas, and geometric structure. This behavior is rate and temperature sensitive, therefore, describing structure-property characteristics requires data collected across several types of instruments. Data included are from structure imaging via micro-computed tomography, finite deformation mechanical measurements from universal test systems with full-field displacement and strain, and visco-thermo-elastic properties from dynamic mechanical analysis. These data facilitate modeling and design efforts in foam mechanics, e.g., homogenization, direct numerical simulation, or phenomenological fitting. The data framework is implemented using data services and software from the Materials Data Facility of the Center for Hierarchical Materials Design.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mesoscale simulations of pressure-shear loading of granular tungsten carbide

Numerical simulations of pressure-shear loading of a granular material are performed using the shock physics code CTH. A simple mesoscale model for the granular material is used that consists of a randomly packed arrangement of solid circular or spherical grains of uniform size separated by vacuum. The grain material is described by a simple shock equation of state, elastic perfectly plastic strength model, and fracture model with baseline parameters for WC taken from previous mesoscale modeling work. Simulations using the baseline material parameters are performed at the same initial conditions of pressure-shear experiments on dry WC powders. Except for some localized flow regions appearing in simulations with an approximate treatment of sliding interfaces among grains, the samples respond elastically during shear, which is in contrast to experimental observations. By extending the simulations to higher shear wave amplitudes, macroscopic shear failure of the simulated samples is observed with the shear strength increasing with increasing stress confinement. The shear strength is also found to be strongly dependent on the grain interface treatment and on the fracture stress of the grains, though the variation in shear strength due to fracture stress decreases with increasing stress confinement. At partial compactions, the transverse velocity histories show strain-hardening behavior followed by formation of a shear interface that extends through the transverse dimensions of the sample. Near full compaction, no strain hardening is observed and, instead, the sample transitions sharply from an elastic response to formation of an internal shear interface. Agreement with experiment is shown to worsen with increasing confinement stress with simulations overpredicting the shear strengths measured in experiment. The source of the disagreement can be ultimately attributed to the Eulerian nature of the simulations, which do not treat contact and fracture realistically.

36 MATERIALS SCIENCE↗