Search NASA⌕ Search

SEARCH · Search NASA

Results for “strain simulation”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Measurement and analysis of critical CTOA for an aluminum alloy sheet

The stable tearing behavior of thin sheets of 2024-T3 aluminum alloy was investigated for middle crack tension, M(T), and compact tension, C(T), specimens. The surface crack-tip opening angle (CTOA), applied loads, crack extension, and local displacements were measured. A critical CTOA fracture criterion was incorporated into a two-dimensional, elastic plastic finite element code and used to simulate the experimental fracture behavior. The CTOA measurements and observations of the fracture surfaces have shown that large values for surface CTOA were observed for small crack extensions (less than the sheet thickness); substantial tunneling of the crack was associated with small crack extensions; crack tunneling in the M(T) specimen was less than that observed in the C(T) configuration; for larger crack extensions, the measured CTOA values were determined to be approximately 6 degrees for both the M(T) and C(T) configuration; and for larger crack extensions, crack tunneling remained constant. The two-dimensional finite element predictions of fracture behavior assumed a constant critical CTOA value of 6 degrees and accounted for local crack tip constraint with a plane strain core of elements ahead of the crack tip. The plane strain core extended 5 mm above the crack plane. The simulations were within +/- 4 percent of the maximum applied load for the C(T) tests within 2 percent for the M(T) tests.

Dawicke, D. S.↗
DOE OSTI2023

Ion-beam radiation-induced Eshelby transformations: The mean and variance in hydrostatic and shear residual stresses

Ion beam plays a pivotal role in ion implantations and the fabrication of nanostructures. However, there lacks a quantitative model to describe the residual stresses associated with the ion-beam radiation. Radiation-induced residual stress/transformation strain have been mostly recognized in the hydrostatic sub strain space. Here, in this study, we use molecular dynamics (MD) simulations to show that the response of a material to irradiation is generally anisotropic that depends on the ion-beam direction, and should be described using tensorial quantities. We demonstrate that accelerator-based ion beam irradiation, combined with the intrinsic lattice anisotropy and externally induced anisotropy (such as anisotropic mechanical loadings), causes radiation-actuated shear transformation strains in addition to hydrostatic expansion. We map out these complex correlations for several materials. Radiation-induced defects are shown to be responsible for residual shear stresses in the manner of Eshelby inclusion transformation. We propose such tensorial response model should be considered for accurate nanoscale fabrication using ion-beam irradiation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1978

Analysis and Test of Deep Flaws in Thin Sheets of Aluminum and Titanium. Volume 2: Crack Opening Displacement and Stress-Strain Data

Static fracture tests were performed on surface flawed specimens of aluminum and titanium alloys. A simulated proof overload cycle was applied prior to all of the cyclic tests. Variables included in each test series were flaw shapes and thickness. Additionally, test temperature was a variable for the aluminum test series. The crack opening displacement and stress-strain data obtained are presented.

Finger, R. W.↗
DOE OSTI2021

Nonlinear convergence in contact mechanics: Immersed boundary finite volume

In this report we present an immersed boundary finite volume (IBM) method for simulating quasistatic contact mechanics of linearly elastic domains at small strains. In IBM, all external boundaries and internal contacts of an object are represented by embedded surfaces inside a Cartesian mesh, which need not conform to the grid lines. The contact constraints consist of the non-penetrability condition and Coulomb’s friction law, which are discretized using special interpolation stencils and enforced via penalty parameters. The resulting nonlinear system depends on displacement unknowns only. To solve it, we use the Newton method but find that it diverges frequently. To understand the divergence pattern, we analyze a simplified 2-cell problem and show that the global convergence of Newton cannot be ensured for any choice of penalty parameters. We thus propose a modified Newton solver, which guarantees convergence for the 2-cell problem and is numerically verified to converge for all the challenging simulations considered herein. While both 1 st - and 2 nd -order variants of IBM, in displacement unknowns, are proposed, the modified Newton solver applies only to the 1 st -order variant.

42 ENGINEERING↗
NASA NTRS1992

Laboratory evaluation of footings for lunar telescopes

Presented here are the results of laboratory experiments with diffferent footing shapes for lunar telescopes. These experiments used a variety of soils including some to simulate regolith response. Based on what is known of regolith and regolith-structure interaction, a shallow-multiple-contact points footing foundation can be adequately designed to support lunar telescopes. Plane-strain load-displacement tests were conducted with different footings and different lunar simulants in a deep transparent plexiglass container. The model footings considered include the rectangular, hemispherical, and spudcan designs. Simulants used to reproduce the mechanical properties of the lunar regolith were fly ash, crushed basalt with and without glass, and a processed lunar simulant. Load-displacement curves were obtained for the different footings in Ottawa sand and in the crushed basalt with glass. The spudcan footing was found to be self-digging and yet stiff, thus providing excellent lateral stability in a large variety of soils.

Chua, Koon M.↗
NASA NTRSDate not supplied

Particle Impact Simulation and Ignition Prediction

An experimentally calibrated tool is needed to predict if a system is susceptible to failure by particle impact ignition (PI) based on use conditions, materials, and flow geometry. This tool will accelerate new components, evaluating existing hardware, and help disposition anomalies. Conduct particle impact testing with in-situ diagnostics and complementary simulations on subset of key engineering materials (IN718, M400, 316L, 6061, Ti64, Zr) to develop a proof-of-concept predictive tool for assessing the risk of PI for idealized geometries (spherical particles) in realistic environments. Assess particle/target interactions (coefficient of restitution, ignition, kindling) using instrumented particle impact rigs while systematically varying key parameters (materials, particle size, environment, target configuration). Determine key field variables (temperature, strain, stress) in particle impacts using Multiphysics finite element and hydrocode simulations validated through comparison with experimental measurements and observations. Synthesize experiments and simulations into constitutive models for PI that can be integrated with existing computational fluid dynamics (CFD) and Debris Transport Analysis (DTA) tools in future efforts

particle impact↗
NASA NTRS2024

Particle Impact Simulation and Ignition Prediction

An experimentally calibrated tool is needed to predict if a system is susceptible to failure by particle impact ignition (PI) based on use conditions, materials, and flow geometry. This tool will accelerate new components, evaluating existing hardware, and help disposition anomalies. - Conduct particle impact testing with in-situ diagnostics and complementary simulations on subset of key engineering materials (IN718, M400, 316L, 6061, Ti64, Zr) to develop a proof-of-concept predictive tool for assessing the risk of PI for idealized geometries (spherical particles) in realistic environments. - Assess particle/target interactions (coefficient of restitution, ignition, kindling) using instrumented particle impact rigs while systematically varying key parameters (materials, particle size, environment, target configuration). - Determine key field variables (temperature, strain, stress) in particle impacts using Multiphysics finite element and hydrocode simulations validated through comparison with experimental measurements and observations. - Synthesize experiments and simulations into constitutive models for PI that can be integrated with existing computational fluid dynamics (CFD) and Debris Transport Analysis (DTA) tools in future efforts.

Jonathan Tylka↗
NASA NTRS1999

Exercise Countermeasures for Bone Loss During Space Flight: A Method for the Study of Ground Reaction Forces and their Implications for Bone Strain

Effective countermeasures to prevent loss of bone mineral during long duration space flight remain elusive. Despite an exercise program on MIR flights, the data from LeBlanc et al. (1996) indicated that there was still a mean rate of loss of bone mineral density in the proximal femur of 1.58% per month (n=18, flight duration 4 - 14.4 months). The specific mechanisms regulating bone mass are not known, but most investigators agree that bone maintenance is largely dependent upon mechanical demand and the resultant local bone strains. A plausible hypothesis is that bone loss during space flight, such as that reported by LeBlanc et al. (1996), may result from failure to effectively load the skeleton in order to generate localized bone strains of sufficient magnitude to prevent disuse osteoporosis. A variety of methods have been proposed to simulate locomotor exercise in reduced gravity. In such simulations, and in an actual microgravity environment, a gravity replacement load (GRL) must always be added to return the exercising subject to the support surface and the resulting skeletal load is critically dependent upon the magnitude of the GRL. To our knowledge, GRLs during orbital flight have only been measured once (on STS 81) and it is likely that most or all prior treadmill exercise in space has used GRLs that were less than one body weight. McCrory (1997) has shown that subjects walking and running in simulated zero-G can tolerate GRLs of 1 if an appropriate harness is used. Several investigators have attempted to measure in vivo strains and forces in the bones of humans, but have faced ethical and technical limitations. The anteromedial aspect of the tibial midshaft has been a common site for the placement of strain gauges; one reason to measure strains in the anterior tibia is that this region is surgically accessible. Aamodt et al. (1997) were able to measure strains on the lateral surface of the proximal femur only because their experimental subjects were already scheduled for hip surgery. Lu et al. (1997) used an instrumented massive proximal femoral prosthesis along with electromyographic measurements to demonstrate that femoral forces depend on muscular activity. These analyses of in vivo bone mechanics are valuable. The invasive nature of the procedures involved, however, limits both the number of subjects and the number of strain gauge locations. Further, the results of these studies may be confounded by the inclusion of subjects with pathological conditions. Gross et al. (1992) measured strain at three locations on the equine third metacarpal and used those data to construct a computer model of the internal strain environment of the bone. An analogous placement of multiple gauges in living humans would be difficult and potentially hazardous because of the depth of soft tissue overlying the tibia and femur.

Peterman, M.↗
DOE OSTI2022

Development of integrated screening, cultivar optimization, and verification research (DISCOVR): A coordinated research-driven approach to improve microalgal productivity, composition, and culture stability for commercially viable biofuels production

To address major knowledge gaps and barriers to the commercial development of algal biomass for biofuels and co-products, a collaborative consortium, Development of Integrated Screening, Cultivar Optimization, and Verification Research (DISCOVR), was established in 2016. Funded by the U.S. Department of Energy (DOE) Bioenergy Technologies Office (BETO), this consortium constitutes a partnership between four DOE national laboratories - Pacific Northwest National Laboratory (PNNL), Los Alamos National Laboratory (LANL), the National Renewable Energy Laboratory (NREL), and Sandia National Laboratories (SNL) - and the Arizona Center for Algae Technology and Innovation (AzCATI) at Arizona State University. To address the barriers of strain selection for achieving high seasonal productivities with a suitable composition and culture resilience, a tiered strain down-selection pipeline is implemented. At Tier I, the temperature and salinity tolerance of strains is determined in flask cultures; at Tier II, the areal biomass productivity and composition is determined in climate-simulation photobioreactors; at Tier III, the productivity and culture stability are determined in outdoor raceways. The top performing strains move forward to long-term testing at the algae testbed site at AzCATI to generate annual biomass productivity data. Concurrent to the strain down-selection in the DISCOVR pipeline, hypotheses for increasing biomass productivity, shifting biomass composition to enhance intrinsic value, and improving culture stability and resistance to pests are also tested. Techno-economic analyses are carried out to determine whether promising findings from laboratory studies or proposed modifications in outdoor pond cultivation conditions translate into reductions in the minimum biomass selling price (MBSP). Notably, in the three years following the launch of DISCOVR, annual biomass productivity has increased from 11.7 to 17.6 g m -2 day -1 , resulting in an MBSP decrease from 824 to 611 $ ton -1 .

09 BIOMASS FUELS↗
DOE OSTI2023

Microstructure characterization and elastic-plastic self-consistent simulation studies of anisotropic deformation of β -tin

Uniaxial compression tests of 99.9% pure polycrystalline $β$-tin (Sn) were conducted at various strain rates (10 -3 /s, 10 -1 /s and 1/s) and temperatures (294 K and 193 K) to understand the effects of these variables on the stress–strain response, microstructure evolution, and stress relaxation behavior. Multiple Sn specimens were subjected to complex compressive loading/unloading/reloading paths at different strain rates. Specimens initially compressed at higher strain rates showed strain rate-dependent texture evolution and more pronounced relaxation upon reloading, as compared to those pre-strained at lower rates at room temperature (RT). Compression tests conducted at low temperature (193 K) revealed increased strength, similar to the enhanced strength observed when the strain rate was increased at RT. Neutron diffraction was employed to characterize the initial and final bulk textures of the RT specimens. Electron backscatter diffraction was utilized to examine the crystallographic grain orientation and morphology, thereby identifying the signatures of dislocation-mediated deformation, grain refinement, recrystallization, and twinning behaviors. Elastic-plastic self-consistent simulations were performed to investigate the deformation modes responsible for the strain rate-dependent macroscopic stress–strain response and texture evolution. Effects of crystallographic orientation on stress relaxation behavior was also examined. In conclusion, the model predictions are in reasonable agreement with experimental observations.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Study of strain boundary conditions and GaAs buffer sizes in InGaAs quantum dots

NEMO 3-D has been developed for the simulation of electronic structure in self-assembled InGaAs quantum dots on GaAs substrates. Typical self-assembled quantum dots in that material system contain about 0.5 to 1 million atoms. Effects of strain by the surrounding GaAs buffer modify the electronic structure inside the quantum dot significantly and a large GaAs buffer must be included in the strain and electronic structure.

transport tunneling quantum↗
DOE OSTI2026

A finite viscoelastic constitutive model for low to high strain rate response of elastomers with application of strain rate-induced glass transition

Amorphous elastomers exhibit significant rate-stiffening and unique viscous flow characteristics across a wide range of strain rates, often undergoing glass transition above a strain rate threshold. We have developed a thermodynamically-consistent and micromechanically-inspired constitutive model for soft elastomeric materials to capture the rate-dependent stress-strain behavior and hysteresis when subjected to low to high strain rates. Here, our proposed constitutive model encapsulates the viscous flow of materials through molecular motion at low strain rates and local rearrangement and alignment of the molecules trying to overcome the intermolecular resistance at high strain rates, essentially covering the glass transition. We applied our constitutive model to uniaxial compression experiments performed at low and high strain rates for polyborosiloxane (PBS) to identify the material parameters, and subsequently, performed numerical simulations of single and multi-cycle compression, stress relaxation, and small amplitude oscillatory tension-compression. Our analyses indicate that the model predicts higher total energy dissipation with increasing strain rate; however, dissipation associated with molecular relaxation decreases (forming a cusp) because, beyond a crossover strain rate, intermolecular rearrangement and alignment become dominant, which is consistent with the onset of the glass transition. For cyclic loading-unloading, we observed that dissipation over a cycle remains constant at low strain rates but decreases non-monotonically at high strain rates before becoming constant, with the peak stress over the cycle becoming higher, which can be interpreted as more loading being carried elastically by the polymer network as the intermolecular rearrangement process occurs. Additionally, our model was able to predict the qualitative nature of the storage modulus and loss modulus in the limit of small strain over a wide range of frequency sweeps.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Crack Growth Prediction Methodology for Multi-Site Damage: Layered Analysis and Growth During Plasticity

A finite element program has been developed to perform quasi-static, elastic-plastic crack growth simulations. The model provides a general framework for mixed-mode I/II elastic-plastic fracture analysis using small strain assumptions and plane stress, plane strain, and axisymmetric finite elements. Cracks are modeled explicitly in the mesh. As the cracks propagate, automatic remeshing algorithms delete the mesh local to the crack tip, extend the crack, and build a new mesh around the new tip. State variable mapping algorithms transfer stresses and displacements from the old mesh to the new mesh. The von Mises material model is implemented in the context of a non-linear Newton solution scheme. The fracture criterion is the critical crack tip opening displacement, and crack direction is predicted by the maximum tensile stress criterion at the crack tip. The implementation can accommodate multiple curving and interacting cracks. An additional fracture algorithm based on nodal release can be used to simulate fracture along a horizontal plane of symmetry. A core of plane strain elements can be used with the nodal release algorithm to simulate the triaxial state of stress near the crack tip. Verification and validation studies compare analysis results with experimental data and published three-dimensional analysis results. Fracture predictions using nodal release for compact tension, middle-crack tension, and multi-site damage test specimens produced accurate results for residual strength and link-up loads. Curving crack predictions using remeshing/mapping were compared with experimental data for an Arcan mixed-mode specimen. Loading angles from 0 degrees to 90 degrees were analyzed. The maximum tensile stress criterion was able to predict the crack direction and path for all loading angles in which the material failed in tension. Residual strength was also accurately predicted for these cases.

James, Mark Anthony↗
DOE OSTI2024

Anomalous temperature dependence of elastic limit in metallic glasses

Abstract Understanding the atomistic mechanisms of inelastic deformation in metallic glasses (MGs) remains challenging due to their amorphous structure, where local carriers of plasticity cannot be easily defined. Using molecular dynamics (MD) simulations, we analyzed the onset of inelastic deformation in CuZr MGs, specifically the temperature dependence of the elastic limit, in terms of localized shear transformation (ST) events. We find that although the ST events initiate at lower strain with increasing temperature, the elastic limit increases with temperature in certain temperature ranges. We explain this anomalous behavior through the framework of an energy-strain landscape (ESL) constructed from high-throughput strain-dependent energy barrier calculations for the ST events identified in the MD simulations. The ESL reveals that the anomalous behavior is caused by the transition of ST events from irreversible to reversible with increasing temperature. An analytical formulation is developed to predict this transition and the temperature dependence of the elastic limit.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Multiscale analysis of large twist ferroelectricity and swirling dislocations in bilayer hexagonal boron nitride

With its atomically thin structure and intrinsic ferroelectric properties, heterodeformed bilayer hexagonal boron nitride (hBN) has gained prominence in next-generation non-volatile memory applications. However, studies to date have focused almost exclusively on small-twist bilayer hBN, leaving the question of whether ferroelectricity can persist under small heterostrain and large heterodeformation entirely unexplored. In this work, we establish the crystallographic origin of ferroelectricity in bilayer hBN configurations heterodeformed relative to high-symmetry configurations such as AA-stacking and 21.786789° twisted configurations (Σ7), using Smith normal form bicrystallography. We then demonstrate out-of-plane ferroelectricity in bilayer hBN across configurations vicinal to both the AA and Σ7 stackings. Atomistic simulations reveal that AA-vicinal systems support ferroelectricity under both small twist and small strain, with polarization switching in the latter governed by the deformation of swirling dislocations rather than the straight interface dislocations seen in the former. For Σ7-vicinal systems, where existing interatomic potentials underperform particularly under extreme out-of-plane compression, we develop a density-functional-theory-informed continuum framework—the bicrystallography-informed frame-invariant multiscale (BFIM) model, which captures out-of-plane ferroelectricity in heterodeformed configurations vicinal to Σ7 stacking. Interface dislocations in these large heterodeformed bilayer configurations exhibit markedly smaller Burgers vectors compared to interface dislocations in small-twist and small-strain bilayer hBN. The BFIM model reproduces experimental results and provides a powerful, computationally efficient framework for predicting ferroelectricity in large-unit-cell heterostructures where atomistic simulations are prohibitively expensive.

Ahmed, Md Tusher [Univ. of Illinois at Urbana-Cham↗
NASA NTRS2004

Using the GeoFEST Faulted Region Simulation System

GeoFEST (the Geophysical Finite Element Simulation Tool) simulates stress evolution, fault slip and plastic/elastic processes in realistic materials, and so is suitable for earthquake cycle studies in regions such as Southern California. Many new capabilities and means of access for GeoFEST are now supported. New abilities include MPI-based cluster parallel computing using automatic PYRAMID/Parmetis-based mesh partitioning, automatic mesh generation for layered media with rectangular faults, and results visualization that is integrated with remote sensing data. The parallel GeoFEST application has been successfully run on over a half-dozen computers, including Intel Xeon clusters, Itanium II and Altix machines, and the Apple G5 cluster. It is not separately optimized for different machines, but relies on good domain partitioning for load-balance and low communication, and careful writing of the parallel diagonally preconditioned conjugate gradient solver to keep communication overhead low. Demonstrated thousand-step solutions for over a million finite elements on 64 processors require under three hours, and scaling tests show high efficiency when using more than (order of) 4000 elements per processor. The source code and documentation for GeoFEST is available at no cost from Open Channel Foundation. In addition GeoFEST may be used through a browser-based portal environment available to approved users. That environment includes semi-automated geometry creation and mesh generation tools, GeoFEST, and RIVA-based visualization tools that include the ability to generate a flyover animation showing deformations and topography. Work is in progress to support simulation of a region with several faults using 16 million elements, using a strain energy metric to adapt the mesh to faithfully represent the solution in a region of widely varying strain.

Geophyical Finite Element Simulation Tool (GeoFEST↗
DOE OSTI2023

Ab-initio simulations of coherent phonon-induced pumping of carriers in zirconium pentatelluride

Abstract Laser-driven coherent phonons can act as modulated strain fields and modify the adiabatic ground state topology of quantum materials. Here we use time-dependent first-principles and effective model calculations to simulate the effect of the coherent phonon induced by strong terahertz electric field on electronic carriers in the topological insulator ZrTe 5 . We show that a coherent A 1g Raman mode modulation can effectively pump carriers across the band gap, even though the phonon energy is about an order of magnitude smaller than the equilibrium band gap. We reveal the microscopic mechanism of this effect which occurs via Landau–Zener-Stückelberg tunneling of Bloch electrons in a narrow region in the Brillouin zone center where the transient energy gap closes when the system switches from strong to weak topological insulator. The quantum dynamics simulation results are in excellent agreement with recent pump-probe experiments in ZrTe 5 at low temperature.

36 MATERIALS SCIENCE↗
NASA NTRS1978

Computer-aided analysis and design of the shape rolling process for producing turbine engine airfoils

Mild steel (AISI 1018) was selected as model cold-rolling material and Ti-6Al-4V and INCONEL 718 were selected as typical hot-rolling and cold-rolling alloys, respectively. The flow stress and workability of these alloys were characterized and friction factor at the roll/workpiece interface was determined at their respective working conditions by conducting ring tests. Computer-aided mathematical models for predicting metal flow and stresses, and for simulating the shape-rolling process were developed. These models utilize the upper-bound and the slab methods of analysis, and are capable of predicting the lateral spread, roll-separating force, roll torque and local stresses, strains and strain rates. This computer-aided design (CAD) system is also capable of simulating the actual rolling process and thereby designing roll-pass schedule in rolling of an airfoil or similar shape. The predictions from the CAD system were verified with respect to cold rolling of mild steel plates. The system is being applied to cold and hot isothermal rolling of an airfoil shape, and will be verified with respect to laboratory experiments under controlled conditions.

Lahoti, G. D.↗