Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic model”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A parametric study of slow dynamic nonlinear elasticity with comparisons to models

Several phenomenological models that aspire to quantitative description of anomalous nonlinear mesoscopic elasticity are reviewed and compared with laboratory measurements. This class of nonlinearity, best known perhaps for slow dynamics and aging, is seen widely in imperfectly consolidated granular solids but is not well understood. Typical slow dynamic tests show that a modest conditioning oscillatory "pump" strain depresses material stiffness, which then recovers like the logarithm of time after conditioning ceases. Several phenomenological models based on physical arguments have been proposed that predict the material stiffness response to arbitrary pump strain histories during conditioning and recovery. Approximate closed form and numerical solutions to the models are presented that predict the quantitative influence of three key pump parameters: the pump's strain amplitude, the pump's strain rate, and the pump’s duration. Laboratory measurements on Berea sandstone, concrete and a confined single aluminum bead find that slow dynamic responses are linear in pump strain and independent of pump frequency. Measurements also show that, after pump-off, stiffness recovers over times far longer than the pump duration. These observations and others are compared to model predictions. One of the considered models, based on a picture of fast brittle damage and slow healing, successfully matches all these behaviors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Integrating Ultra-Coarse-Grained Protein Models into Accessible Workflows for Multiscale Molecular Dynamics

To capture protein conformational transitions using molecular dynamics (MD), several simulation resolutions covering different spatial and temporal scales are typically needed. All-atom (AA) simulations provide fine resolution, but are computationally infeasible for large systems over longer durations. Coarse-grained (CG) and ultra-coarse-grained (UCG) models have a lower resolution and computational cost while still being able to conserve essential protein features. Prior work on a Multiscale Machinelearned Modeling Infrastructure (MuMMI) combined both AA and CG simulations to study RAS-RAF protein interactions, leveraging CG models for longer time scales and using AA to investigate unusual conformations in greater detail. However, MuMMI is still resource-intensive, and this study aims to maximize exploration of the protein conformational space while reducing computational cost. In this paper, we build on prior work that integrates UCG models based on heterogeneous elastic network modeling (hENM) into the MuMMI workflow. We demonstrate that UCG models enable accurate sampling of protein conformations, focusing on simulating RAS-RAF protein interactions. Using higher-resolution CG Martini simulation data, we can automatically refine intramolecular interactions in UCG models. We present a scalable Python package that uses fluctuations observed in higher-resolution CG Martini simulations to estimate bond coefficients of the UCG model. We built novel machine learning-based backmapping methods to recover more detailed CG Martini structures from UCG structures, using diffusion models to learn the mapping between scales. Finally, we present UCG-mini-MuMMI, an accessible and less compute-intensive version of MuMMI as a resource for the scientific community. Incorporating UCG models into MD studies is applicable to a broad range of systems and proteins, and our study offers insights into the advantages and limitations of these methods.

Chemical structure↗
DOE OSTI2015

Probing the folded state and mechanical unfolding pathways of T4 lysozyme using all-atom and coarse-grained molecular simulation

The Bacteriophage T4 Lysozyme (T4L) is a prototype modular protein comprised of an N-terminal and a C-domain domain, which was extensively studied to understand the folding/unfolding mechanism of modular proteins. To offer detailed structural and dynamic insights to the folded-state stability and the mechanical unfolding behaviors of T4L, we have performed extensive equilibrium and steered molecular dynamics simulations of both the wild-type (WT) and a circular permutation (CP) variant of T4L using all-atom and coarse-grained force fields. Our all-atom and coarse-grained simulations of the folded state have consistently found greater stability of the C-domain than the N-domain in isolation, which is in agreement with past thermostatic studies of T4L. While the all-atom simulation cannot fully explain the mechanical unfolding behaviors of the WT and the CP variant observed in an optical tweezers study, the coarse-grained simulations based on the Go model or a modified elastic network model (mENM) are in qualitative agreement with the experimental finding of greater unfolding cooperativity in the WT than the CP variant. Interestingly, the two coarse-grained models predict different structural mechanisms for the observed change in cooperativity between the WT and the CP variant—while the Go model predicts minor modification of the unfolding pathways by circular permutation (i.e., preserving the general order that the N-domain unfolds before the C-domain), the mENM predicts a dramatic change in unfolding pathways (e.g., different order of N/C-domain unfolding in the WT and the CP variant). Based on our simulations, we have analyzed the limitations of and the key differences between these models and offered testable predictions for future experiments to resolve the structural mechanism for cooperative folding/unfolding of T4L.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1988

Cyclic uniaxial and biaxial hardening of type 304 stainless steel modeled by the viscoplasticity theory based on overstress

The isotropic theory of viscoplasticity based on overstress does not use a yield surface or a loading and unloading criterion. The inelastic strain rate depends on overstress, the difference between the stress and the equilibrium stress, and is assumed to be rate dependent. Special attention is paid to the modeling of elastic regions. For the modeling of cyclic hardening, such as observed in annealed Type 304 stainless steel, and additional growth law for a scalar quantity which represents the rate independent asymptotic value of the equilibrium stress is added. It is made to increase with inelastic deformation using a new scalar measure which differentiates between nonproportional and proportional loading. The theory is applied to correlate uniaxial data under two step amplitude loading including the effect of further hardening at the high amplitude and proportional and nonproportional cyclic loadings. Results are compared with corresponding experiments.

Yao, David↗
DOE OSTI2020

Numerical Modeling of Anisotropic Elastic-Wave Sensitivity Propagation for Optimal Design of Time-Lapse Seismic Surveys

Reliable subsurface time-lapse seismic monitoring is crucial for many geophysical applications, such as enhanced geothermal system characterization, geologic carbon utilization and storage, and conventional and unconventional oil/gas reservoir characterization, etc. In this work, we develop an elastic-wave sensitivity propagation method for optimal design of cost-effective time-lapse seismic surveys considering the fact that most of subsurface geologic layers and fractured reservoirs are anisotropic instead of isotropic. For anisotropic media, we define monitoring criteria using qP- and qS-wave sensitivity energies after decomposing qP- and qS-wave components from the total elastic-wave sensitivity wavefield using a hybrid time- and frequency-domain approach. Geophones should therefore be placed at locations with significant qP- and qS-wave sensitivity energies for cost-effective time-lapse seismic monitoring in an anisotropic geology setting. Our numerical modeling results for a modified anisotropic Hess model demonstrate that, compared with the isotropic case, subsurface anisotropy changes the spatial distributions of elastic-wave sensitivity energies. Consequently, it is necessary to consider subsurface anisotropies when designing the spatial distribution of geophones for cost-effective time-lapse seismic monitoring. This finding suggests that it is essential to use our new anisotropic elastic-wave sensitivity modeling method for optimal design of time-lapse seismic surveys to reliably monitor the changes in subsurface reservoirs, fracture zones or target monitoring regions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Examination of the collision force method for analyzing the responses of simple containment/deflection structures to impact by one engine rotor blade fragment

An approximate collision analysis, termed the collision-force method, was developed for studying impact-interaction of an engine rotor blade fragment with an initially circular containment ring. This collision analysis utilizes basic mass, material property, geometry, and pre-impact velocity information for the fragment, together with any one of three postulated patterns of blade deformation behavior: (1) the elastic straight blade model, (2) the elastic-plastic straight shortening blade model, and (3) the elastic-plastic curling blade model. The collision-induced forces are used to predict the resulting motions of both the blade fragment and the containment ring. Containment ring transient responses are predicted by a finite element computer code which accommodates the large deformation, elastic-plastic planar deformation behavior of simple structures such as beams and/or rings. The effects of varying the values of certain parameters in each blade-behavior model were studied. Comparisons of predictions with experimental data indicate that of the three postulated blade-behavior models, the elastic-plastic curling blade model appears to be the most plausible and satisfactory for predicting the impact-induced motions of a ductile engine rotor blade and a containment ring against which the blade impacts.

Zirin, R. M.↗
DOE OSTI2022

Efficient machine-learning model for fast assessment of elastic properties of high-entropy alloys

We combined descriptor-based analytical models for stiffness-matrix and elastic-moduli with mean-field methods to accelerate assessment of technologically useful properties of high-entropy alloys, such as strength and ductility. Model training for elastic properties uses Sure-Independence Screening (SIS) and Sparsifying Operator (SO) method yielding an optimal analytical model, constructed with meaningful atomic features to predict target properties. Computationally inexpensive analytical descriptors were trained using a database of elastic properties determined from density functional theory for binary and ternary subsets of Nb-Mo-Ta-W-V refractory alloys. The optimal Elastic-SISSO models, extracted from an exponentially large feature space, give an extremely accurate prediction of target properties, similar to or better than other models, with some verified from existing experiments. Here we also show that electronegativity variance and elastic-moduli can directly predict trends in ductility and yield strength of refractory HEAs, and reveals promising alloy concentration regions.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A discrete element method representation of an anisotropic elastic continuum

A method for modeling cubically anisotropic elasticity within the discrete element method is presented. The discrete element method (DEM) is an approach originally intended for modeling granular materials (sand, soil, and powders); however, recent developments have usefully extended it to model stochastic mechanical processes in monolithic solids which, to date, have been assumed to be elastically isotropic. The method presented here for efficiently capturing cubic elasticity in DEM is an important prerequisite for further extending DEM to capture the influence of elastic anisotropy on the mechanical response of polycrystals, composites, etc. The system demonstrated here uses a directionally assigned stiffness in the bonds between adjacent elements and includes separate schemes for achieving anisotropy with Zener ratios greater and smaller than one. Finally, the model framework is presented along with an analysis of the accessible space of elastic properties that can be modeled and an artificial neural network interpolation scheme for mapping input parameters to model elastic behavior.

42 ENGINEERING↗
DOE OSTI2013

Folding model analysis of pion elastic and inelastic scattering from {sup 6}Li and {sup 12}C

{pi}{sup {+-}}-Nucleus scattering cross sections are calculated applying the Watanabe superposition model with a phenomenological Woods-Saxon potential. The phenomenological potential parameters are searched for {pi}{sup {+-}} scattering from {sup 6}Li and {sup 12}C to reproduce not only differential elastic cross sections but also inelastic and total and reaction cross sections at pion kinetic energies from 50 to 672 MeV. The optical potentials of {sup 6}Li and {sup 12}C are calculated in terms of the alpha particle and deuteron optical potentials. Inelastic scattering has been analyzed using the distorted waves from elastic-scattering data. The values of deformation lengths thus obtained compare very well with the ones reported earlier.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2004

Stability Analysis of the Slowed-Rotor Compound Helicopter Configuration

The stability and control of rotors at high advance ratio are considered. Teetering, articulated, gimbaled, and rigid hub types are considered for a compound helicopter (rotor and fixed wing). Stability predictions obtained using an analytical rigid flapping blade analysis, a rigid blade CAMRAD II model, and an elastic blade CAMRAD II model are compared. For the flapping blade analysis, the teetering rotor is the most stable, 5howing no instabilities up to an advance ratio of 3 and a Lock number of 18. With an elastic blade model, the teetering rotor is unstable at an advance ratio of 1.5. Analysis of the trim controls and blade flapping shows that for small positive collective pitch, trim can be maintained without excessive control input or flapping angles.

Johnson, Wayne↗
NASA NTRS2007

Stability Analysis of the Slowed-Rotor Compound Helicopter Configuration

The stability and control of rotors at high advance ratio are considered. Teetering, articulated, gimbaled, and rigid hub types are considered for a compound helicopter (rotor and fixed wing). Stability predictions obtained using an analytical rigid flapping blade analysis, a rigid blade CAMRAD II model, and an elastic blade CAMRAD II model are compared. For the flapping blade analysis, the teetering rotor is the most stable, showing no instabilities up to an advance ratio of 3 and a Lock number of 18. A notional elastic blade model of a teetering rotor is unstable at an advance ratio of 1.5, independent of pitch frequency. Analysis of the trim controls and blade flapping shows that for small positive collective pitch, trim can be maintained without excessive control input or flapping angles.

Floros, Matthew W.↗
DOE OSTI2023

Beyond the Stokes approximation: shallow visco-elastic ice-sheet models

Abstract The hypothesis that ice-sheet evolution is only controlled by the long-term non-Newtonian viscous behavior of ice has been challenged by observations indicating that effects like brittle failure, stick-slip sliding, tides and wave action may affect ice-sheet evolution on sub-daily timescales. Over these timescales, the quasi-static-creep approximation is no longer appropriate and elastic effects become important. Simulating elastic effects in ice-sheet models over relevant timescales, however, remains challenging. Here, we show that by including a visco-elastic rheology and reintroducing the oft neglected acceleration term back into the ice-sheet stress balance, we can create a visco-elastic system where the velocity is locally determined and information propagates at the elastic wave speed. Crucially, the elastic wave speed can be treated like an adjustable parameter and set to any value to reproduce a range of phenomena, provided the wave speed is large compared to the viscous velocity. We illustrate the system using three examples. The first two examples demonstrate that the system converges to the steady-state viscous and elastic limits. The third example examines ice-shelf rifting and iceberg calving. This final example hints at the utility of the visco-elastic formulation in treating both long-term evolution and short-term environmental effects.

58 GEOSCIENCES↗
DOE OSTI2019

Dynamics, Radiation, and Overall Energy Budget of Earthquake Rupture With Coseismic Off-Fault Damage

Earthquake ruptures dynamically activate coseismic off–fault damage around fault cores. Systematic field observation efforts have shown the distribution of off–fault damage around main faults, while numerical modeling using elastic–plastic off–fault material models have demonstrated the evolution of coseismic off–fault damage during earthquake ruptures. Laboratory–scale microearthquake experiments have pointed out the enhanced high–frequency radiation due to the coseismic off–fault damage. However, the detailed off–fault fracturing mechanisms, subsequent radiation, and its contribution to the overall energy budget remain to be fully understood because of limitations of current observational techniques and model formulations. Here, we constructed a new physics–based dynamic earthquake rupture modeling framework, based on the combined finite–discrete element method, to investigate the fundamental mechanisms of coseismic off–fault damage, and its effect on the rupture dynamics, the radiation and the overall energy budget. We conducted a 2–D systematic case study with depth and showed the mechanisms of dynamic activation of the coseismic off–fault damage. We found the decrease in rupture velocity and the enhanced high–frequency radiation in near field due to the coseismic off–fault damage. We then evaluated the overall energy budget, which shows a significant contribution of the coseismic off–fault damage to the overall energy budget even at depth, where the damage zone width becomes narrower. Here, the present numerical framework for the dynamic earthquake rupture modeling thus provides new insights into earthquake rupture dynamics with the coseismic off–fault damage.

58 GEOSCIENCES↗
NASA NTRS1987

Effect of helicopter blade dynamics on blade aerodynamic and structural loads

The effect of rotor blade dynamics on aerodynamic and structural loads is examined for a conventional, main-rotor helicopter using a comprehensive rotorcraft analysis (CAMRAD) and flight-test data. The impact of blade dynamics on blade section lift-coefficient time histories is studied by comparing predictions from a rigid-blade analysis and an elastic-blade analysis with helicopter flight test data. The elastic blade analysis better predicts high-frequency behavior of section lift. In addition, components of the blade angle of attack such as elastic blade twist, blade flap rate, blade slope velocity, and inflow are examined as a function of blade mode. Elastic blade motion changed blade angle of attack by a few tenths of a degree, and up to the sixth rotor harmonic. A similar study of the influence of blade dynamics on bending and torsion moments was also conducted. A correlation study comparing predictions from several elastic-blade analyses with flight-test data revealed that an elastic-blade model consisting of only three elastic bending modes (first and second flap and first lag), and two elastic torsion modes was sufficient for good correlation.

Heffernan, Ruth M.↗
NASA NTRS1989

Results of including geometric nonlinearities in an aeroelastic model of an F/A-18

An integrated, nonlinear simulation model suitable for aeroelastic modeling of fixed-wing aircraft has been developed. While the author realizes that the subject of modeling rotating, elastic structures is not closed, it is believed that the equations of motion developed and applied herein are correct to second order and are suitable for use with typical aircraft structures. The equations are not suitable for large elastic deformation. In addition, the modeling framework generalizes both the methods and terminology of non-linear rigid-body airplane simulation and traditional linear aeroelastic modeling. Concerning the importance of angular/elastic inertial coupling in the dynamic analysis of fixed-wing aircraft, the following may be said. The rigorous inclusion of said coupling is not without peril and must be approached with care. In keeping with the same engineering judgment that guided the development of the traditional aeroelastic equations, the effect of non-linear inertial effects for most airplane applications is expected to be small. A parameter does not tell the whole story, however, and modes flagged by the parameter as significant also need to be checked to see if the coupling is not a one-way path, i.e., the inertially affected modes can influence other modes.

Buttrill, Carey S.↗
NASA NTRS1993

Prediction of stable tearing of 2024-T3 aluminum alloy using the crack-tip opening angle approach

In this study, the crack-tip opening angle (CTOA) approach was incorporated into a damage growth finite element program, MADGIC (Micromechanics Analysis and Damage Growth in Composites), and was used to predict stable tearing in a middle-crack tension 2024-T3 aluminum alloy specimen. The MADGIC code is a displacement based finite element program implemented with an incremental elastic-plastic algorithm used to model elastic-plastic behavior and a nodal splitting and nodal force relaxation algorithm used to generate crack surfaces. Predictions of the applied stress as a function of crack extension and applied stress as a function of load-line displacement were in good agreement with experiments and with similar predictions made using an existing finite element program, ZIP2D. In addition, path integrals, namely, the J-integral and T*-integral, were also evaluated and compared with the CTOA approach. There appears to be a weak relationship between the CTOA and the T*-integral evaluated on a specific integration path during crack extension beyond maximum applied stress. This study further verifies that the CTOA can be used as an effective elastic-plastic fracture mechanics parameter to predict crack growth.

Bakuckas, J. G., Jr.↗