Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Material Response 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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Evolution of global and local deformation in additively manufactured octet truss lattice structures

Additively manufactured lattice truss structures, often referred to as architected cellular materials, present significant advantages over conventional structures due to their unique characteristics such as high strength-to-weight ratios and surface area-to-volume ratios. These geometrically complex structures, however, come with concomitant challenges for qualification and inspection. In this study, compression testing interrupted with micro-computed tomography inspection was conducted to monitor the evolution of global and local deformation throughout the loading process of 304L stainless steel octet truss lattice structures. Both two- and three-dimensional image analysis techniques were leveraged to characterize geometric heterogeneities resulting from the laser powder bed fusion manufacturing process as well as track the structure throughout deformation. Variations from model-predicted behavior resulting from these heterogeneities are considered relative to the predicted and actual responses of the structures during compression to better understand, model, and predict the octet truss lattice structure compression response.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Polysulfone shock compressed above the decomposition threshold: Velocimetry and modeling of two-wave structures

Polysulfone was shock loaded to pressures of 14.6–26.2 GPa in a series of gas gun-driven plate-impact experiments measuring material response with embedded electro-magnetic particle velocity gauges and optical velocimetry. The embedded electro-magnetic particle velocity gauges did not show a distinct two-wave structure but did show rounding that suggested a reaction but not a distinct separation of the reactants and product waves. In contrast, the transmission experiments fielded with optical velocimetry, with product pressures ranging from 21.0 to 26.2 GPa, showed well-defined two-wave structures due to shock-driven chemical decomposition of the polymer to products at a higher density. Distinct two-wave structures have not previously been published in literature during polymer compression; here, we observed these two-wave structures at both the polymer/lithium fluoride and polymer/polymethylmethacrylate interfaces. Hydrodynamic simulations of the experiments were performed using a newly constructed SESAME equation of state (EOS) for the polymer, a thermochemical EOS for the decomposition products, and an Arrhenius reaction rate model for the kinetics of the reaction. Finally, simulation results also demonstrated two-wave structures but were unable to quantitatively reproduce either the embedded gauge or optical velocimetry data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Thixotropy, antithixotropy, and viscoelasticity in hysteresis

Thixotropy, antithixotropy, and viscoelasticity are three types of time-dependent dynamics that involve fundamentally different underlying physical processes. Here, we show that the three dynamics exhibit different signatures in hysteresis by examining the fingerprints of the simplest thixotropic kinetic model, a new antithixotropic model that we introduce here, and the Giesekus model. We start by showing that a consistent protocol to generate hysteresis loops is a discrete shear-rate controlled ramp that begins and ends at high shear rates, rather than at low shear rates. Using this protocol, we identify two distinguishing features in the resulting stress versus shear rate loops. The first is the direction of the hysteresis loops: clockwise for thixotropy, but counterclockwise for viscoelasticity and antithixotropy. A second feature is achieved at high ramping rates where all responses lose hysteresis: the viscoelastic response shows a stress plateau at low shear rates due to lack of stress relaxation, whereas the thixotropic and antithixotropic responses are purely viscous with minimal shear thinning or thickening. We establish further evidence for these signatures by experimentally measuring the hysteresis of Laponite suspensions, carbon black suspensions, and poly(ethylene oxide) solutions, each representing a historically accepted example of each class of material behavior. The signatures measured in experiments are consistent with those predicted by the three models. This study reveals different fingerprints in hysteresis loops associated with thixotropy, antithixotropy, and viscoelasticity, which may be helpful in distinguishing the three time-dependent responses.

Mechanics↗
DOE OSTI2023

Cost‐Responsive Optimization of Nickel Nanoparticle Synthesis

Abstract Early‐stage cost evaluation during catalyst development holds the potential to accelerate the commercialization and deployment of advanced catalytic materials for sustainable chemical processes. The modeling and assessment of manufacturing costs as early as the laboratory synthesis scale, for example, focusing on materials costs and synthesis performance metrics, can support the development of an experimental–economic feedback loop that enables rapid insight into cost drivers associated with catalyst synthesis and highlights areas that require focused research and development effort. Ultimately, this feedback loop supports the realization of an economic understanding of the overall synthetic process and highlights opportunities to reduce costs, serving as the foundation for the scale‐up of catalyst manufacturing. Herein, a case study is presented utilizing CatCost, a free and publicly available estimation tool for the evaluation of catalyst manufacturing costs, to perform a cost‐responsive optimization of the synthesis of nickel nanoparticles (Ni NPs). It is demonstrated that reagent substitutions with more cost‐effective analogs, coupled with stoichiometric optimization, afford a 58% reduction in raw materials cost without changing the product yield or properties.

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

Modeling fluvial erosion on regional to continental scales

The fluvial system is a major concern in modeling landform evolution in response to tectonic deformation. Three stream bed types (bedrock, coarse-bed alluvial, and fine-bed alluvial) differ in factors controlling their occurrence and evolution and in appropriate modeling approaches. Spatial and temporal transitions among bed types occur in response to changes in sediment characteristics and tectonic deformation. Erosion in bedrock channels depends upon the ability to scour or pluck bed material; this detachment capacity is often a power function of drainage area and gradient. Exposure of bedrock in channel beds, due to rapid downcutting or resistant rock, slows the response of headwater catchments to downstream baselevel changes. Sediment routing through alluvial channels must account for supply from slope erosion, transport rates, abrasion, and sorting. In regional landform modeling, implicit rate laws must be developed for sediment production from erosion of sub-grid-scale slopes and small channels.

Howard, Alan D.↗
DOE OSTI2025

MPACT Safeguards Modeling: FY25 Update

Sandia National Laboratories develops and maintains several open-source software packages to support material accountancy analyses. This includes the Material Accountancy Performance Indicator Toolkit (MAPIT), the Fissile Facility Flow Modeler (F3M) and the Separation and Safeguards Performance Model Library (SSPM-L). MAPIT is responsible for performing statistical safeguards analyses on bulk and itemized data from nuclear fuel cycle facilities and can operate on real or synthetic data. MAPIT is the only open-source software for such analyses. F3M is a library of modules, built in MATLAB Simulink, that contain pre made blocks to represent different generic fuel cycle processes. These blocks can be used together in a modular fashion to represent and simulate nuclear fuel cycle processes with the goal of improving facility-level accountancy during the design phase. F3M is also an open-source library. Finally, the SSPM-L library is a series of completed models built from F3M. The library includes facility models such as a generic PUREX facility and a fuel fabrication facility. The SSPM-L library is not open source, but is available to collaborators with a relevant use case. These tools include modeling and simulation pipelines to simulate nuclear fuel cycle facilities and the underlying software needed to simulate measurement uncertainty and perform statistical analyses. Together, these tools can perform end-to-end nuclear material accountancy analyses. This report documents the various improvements made to these tools in FY25. Specifically, we added new statistical test, new statistical modeling capabilities, new fuel cycle facility models, and launched a new open-source model component library.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Deep glassy state dynamic data challenge glass models: Elastic models

The idea of an “ideal” glass transition temperature has persisted at least since the work from Kauzmann when it was observed that the entropy of glass-forming liquids extrapolated to below that of the crystal, thus suggesting the need for a phase transition at a finite or non-zero absolute temperature. This thermodynamic paradox was also found to be related to the observations of a diverging of the extrapolated viscosity or relaxation times at a temperature near to this ideal glass transition. Recently, however, we have carried out experiments using both an ancient amber material and an ultra-stable amorphous fluoropolymer that challenge the ideas of the divergence of the viscosity or relaxation times at this ideal glass transition. In the present manuscript we have evaluated two theories of the glass transition that are based on ideas related to elasticity of the amorphous glass-forming material. We find that the models from both J. Dyre and his group and of K.S. Schweizer and his group not only show non-diverging behavior but are also in some agreement with the new data, though still showing somewhat slower dynamics than those observed in the experiments. The work shows that the data are good enough to distinguish between theories and it is suggested that other mechanisms may be needed to fully describe the non-diverging responses of the ultra-stable glasses.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The Role of Unit-Cell Topology in Modulating the Compaction Response of Additively Manufactured Cellular Materials using Simulations and Validation Experiments

Additive manufacturing has enabled a transformational ability to create cellular structures (or foams) with tailored topology. Compared to their monolithic polymer counterparts, cellular structures are potentially suitable for systems requiring materials with high specific energy-absorbing capability to provide enhanced damping. In this work, we demonstrate the utility of controlling unit-cell topology with the intent of obtaining a desired stress–strain response and energy density. Using mesoscale simulations that resolve the unit-cell sub-structures, we validate the role of unit-cell topology in selectively activating a buckling mode and thereby modulating the characteristic stress–strain response. Simulations incorporate a linear viscoelastic constitutive model and a hyperelastic model for simulating large deformation of the polymer under both tension and compression. Simulated results for nine different cellular structures are compared with experimental data to gain insights into three different modes of buckling and the corresponding stress–strain response.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multiphysics Modeling Framework for Soft PVC Gel Sensors with Experimental Comparisons

Polyvinyl chloride (PVC) gels have recently been found to exhibit mechanoelectrical transduction or sensing capabilities under compressive loading applications. This phenomenon is not wholly understood but has been characterized as an adsorption-like phenomena under varying amounts and types of plasticizers. A different polymer lattice structure has also been tested, thermoplastic polyurethane, which showed similar sensing characteristics. This study examines mechanical and electrical properties of these gel sensors and proposes a mathematical framework of the underlying mechanisms of mechanoelectrical transduction. COMSOL Multiphysics is used to show solid mechanics characteristics, electrostatic properties, and transport of interstitial plasticizer under compressive loading applications. The solid mechanics takes a continuum mechanics approach and includes a highly compressive Storakers material model for compressive loading applications. The electrostatics and transport properties include charge conservation and a Langmuir adsorption migration model with variable diffusion properties based on plasticizer properties. Results show both plasticizer concentration gradient as well as expected voltage response under varying amounts and types of plasticizers. Experimental work is also completed to show agreeance with the modeling results.

42 ENGINEERING↗
NASA NTRS2004

The Influence of Mesh Density on the Impact Response of a Shuttle Leading-Edge Panel Finite Element Simulation

A study was performed to examine the influence of varying mesh density on an LS-DYNA simulation of a rectangular-shaped foam projectile impacting the space shuttle leading edge Panel 6. The shuttle leading-edge panels are fabricated of reinforced carbon-carbon (RCC) material. During the study, nine cases were executed with all possible combinations of coarse, baseline, and fine meshes of the foam and panel. For each simulation, the same material properties and impact conditions were specified and only the mesh density was varied. In the baseline model, the shell elements representing the RCC panel are approximately 0.2-in. on edge, whereas the foam elements are about 0.5-in. on edge. The element nominal edge-length for the baseline panel was halved to create a fine panel (0.1-in. edge length) mesh and doubled to create a coarse panel (0.4-in. edge length) mesh. In addition, the element nominal edge-length of the baseline foam projectile was halved (0.25-in. edge length) to create a fine foam mesh and doubled (1.0- in. edge length) to create a coarse foam mesh. The initial impact velocity of the foam was 775 ft/s. The simulations were executed in LS-DYNA version 960 for 6 ms of simulation time. Contour plots of resultant panel displacement and effective stress in the foam were compared at five discrete time intervals. Also, time-history responses of internal and kinetic energy of the panel, kinetic and hourglass energy of the foam, and resultant contact force were plotted to determine the influence of mesh density. As a final comparison, the model with a fine panel and fine foam mesh was executed with slightly different material properties for the RCC. For this model, the average degraded properties of the RCC were replaced with the maximum degraded properties. Similar comparisons of panel and foam responses were made for the average and maximum degraded models.

Jackson, Karen E.↗
DOE OSTI2021

A predictive discrete-continuum multiscale model of plasticity with quantified uncertainty

Multiscale models of materials, consisting of upscaling discrete simulations to continuum models, are unique in their capability to simulate complex materials behavior. The fundamental limitation in multiscale models is the presence of uncertainty in the computational predictions delivered by them. In this work, a sequential multiscale model has been developed, incorporating discrete dislocation dynamics (DDD) simulations and a strain gradient plasticity (SGP) model to predict the size effect in plastic deformations of metallic micro-pillars. The DDD simulations include uniaxial compression of micro-pillars with different sizes and over a wide range of initial dislocation densities and spatial distributions of dislocations. An SGP model is employed at the continuum level that accounts for the size-dependency of flow stress and hardening rate. Sequences of uncertainty analyses have been performed to assess the predictive capability of the multiscale model. The variance-based global sensitivity analysis determines the effect of parameter uncertainty on the SGP model prediction. The multiscale model is then constructed by calibrating the continuum model using the data furnished by the DDD simulations. A Bayesian calibration method is implemented to quantify the uncertainty due to microstructural randomness in discrete dislocation simulations (density and spatial distribution of dislocations) on the macroscopic continuum model prediction (size effect in plastic deformation). Here, the outcomes of this study indicate that the discrete-continuum multiscale model can accurately simulate the plastic deformation of micro-pillars, despite the significant uncertainty in the DDD results. Additionally, depending on the macroscopic features represented by the DDD simulations, the SGP model can reliably predict the size effect in plasticity responses of the micropillars with below 10% of error.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Microstructure dependent shear localization in 316L SS

Shear localization is characterized by the concentration of plastic deformation within narrow bands of material and constitutes the precursor to failure in materials subjected to high strain rates. These bands have been observed experimentally, but classical continuum material models have not been able to fully capture their development. The aim of this research is to investigate the role of microstructural differences in the dynamic shear response of wrought and additively manufactured (AM) 316L SS. Compact Forced Simple Shear specimens for Digital Image Correlation (CFSS-DIC) of wrought and additively manufactured 316L SS were tested on a Split Hopkinson Pressure Bar (SHPB) in conjunction with infrared thermometry measurements. Results indicate that the microstructures of wrought and AM 316L SS exhibit differences that affect the shear response of this materials. In conclusion, CFSS-DIC samples tested on the SHPB did not exhibit recrystallization inside the shear region, fact confirmed by thermometry measurements.

36 MATERIALS SCIENCE↗
NASA NTRS2013

Impact Testing of Aluminum 2024 and Titanium 6Al-4V for Material Model Development

One of the difficulties with developing and verifying accurate impact models is that parameters such as high strain rate material properties, failure modes, static properties, and impact test measurements are often obtained from a variety of different sources using different materials, with little control over consistency among the different sources. In addition there is often a lack of quantitative measurements in impact tests to which the models can be compared. To alleviate some of these problems, a project is underway to develop a consistent set of material property, impact test data and failure analysis for a variety of aircraft materials that can be used to develop improved impact failure and deformation models. This project is jointly funded by the NASA Glenn Research Center and the FAA William J. Hughes Technical Center. Unique features of this set of data are that all material property data and impact test data are obtained using identical material, the test methods and procedures are extensively documented and all of the raw data is available. Four parallel efforts are currently underway: Measurement of material deformation and failure response over a wide range of strain rates and temperatures and failure analysis of material property specimens and impact test articles conducted by The Ohio State University; development of improved numerical modeling techniques for deformation and failure conducted by The George Washington University; impact testing of flat panels and substructures conducted by NASA Glenn Research Center. This report describes impact testing which has been done on aluminum (Al) 2024 and titanium (Ti) 6Al-4vanadium (V) sheet and plate samples of different thicknesses and with different types of projectiles, one a regular cylinder and one with a more complex geometry incorporating features representative of a jet engine fan blade. Data from this testing will be used in validating material models developed under this program. The material tests and the material models developed in this program will be published in separate reports.

Pereira, J. Michael↗
NASA NTRS1983

Solar maximum mission: Ground support programs at the Harvard Radio Astronomy Station

Observations of the spectral characteristics of solar radio bursts were made with new dynamic spectrum analyzers of high sensitivity and high reliability, over the frequency range 25-580 MHz. The observations also covered the maximum period of the current solar cycle and the period of international cooperative programs designated as the Solar Maximum Year. Radio data on shock waves generated by solar flares were combined with optical data on coronal transients, taken with equipment on the SMM and other satellites, and then incorporated into computer models for the outward passage of fast-mode MHD shocks through the solar corona. The MHD models are non-linear, time-dependent and for the most recent models, quasi-three-dimensional. They examine the global response of the corona for different types of input pulses (thermal, magnetic, etc.) and for different magnetic topologies (for example, open and closed fields). Data on coronal shocks and high-velocity material ejected from solar flares have been interpreted in terms of a model consisting of three main velocity regimes.

Maxwell, A.↗
NASA NTRS2002

Finite-Strain Fractional-Order Viscoelastic (FOV) Material Models and Numerical Methods for Solving Them

Fraction-order viscoelastic (FOV) material models have been proposed and studied in 1D since the 1930's, and were extended into three dimensions in the 1970's under the assumption of infinitesimal straining. It was not until 1997 that Drozdov introduced the first finite-strain FOV constitutive equations. In our presentation, we shall continue in this tradition by extending the standard, FOV, fluid and solid, material models introduced in 1971 by Caputo and Mainardi into 3D constitutive formula applicable for finite-strain analyses. To achieve this, we generalize both the convected and co-rotational derivatives of tensor fields to fractional order. This is accomplished by defining them first as body tensor fields and then mapping them into space as objective Cartesian tensor fields. Constitutive equations are constructed using both variants for fractional rate, and their responses are contrasted in simple shear. After five years of research and development, we now possess a basic suite of numerical tools necessary to study finite-strain FOV constitutive equations and their iterative refinement into a mature collection of material models. Numerical methods still need to be developed for efficiently solving fraction al-order integrals, derivatives, and differential equations in a finite element setting where such constitutive formulae would need to be solved at each Gauss point in each element of a finite model, which can number into the millions in today's analysis.

Freed, Alan D.↗
NASA NTRS2010

Micromechanics Modeling of Composites Subjected to Multiaxial Progressive Damage in the Constituents

The high-fidelity generalized method of cells composite micromechanics model is extended to include constituent-scale progressive damage via a proposed damage model. The damage model assumes that all material nonlinearity is due to damage in the form of reduced stiffness, and it uses six scalar damage variables (three for tension and three for compression) to track the damage. Damage strains are introduced that account for interaction among the strain components and that also allow the development of the damage evolution equations based on the constituent material uniaxial stress strain response. Local final-failure criteria are also proposed based on mode-specific strain energy release rates and total dissipated strain energy. The coupled micromechanics-damage model described herein is applied to a unidirectional E-glass/epoxy composite and a proprietary polymer matrix composite. Results illustrate the capability of the coupled model to capture the vastly different character of the monolithic (neat) resin matrix and the composite in response to far-field tension, compression, and shear loading.

Bednarcyk, Brett A.↗
NASA NTRS1987

Nonisothermal elasto-visco-plastic response of shell-type structures

A mathematical model and solution methodologies for analyzing structural response of thin, metallic shell-type structures under large transient, cyclic or static thermomechanical loads is disussed. Among the system responses, which are associated with these load conditions, are thermal buckling and creep buckling. Thus, geometric as well as material-type nonlinearities (of high order) can be anticipated and have been considered in the development of the mathematical model. Furthermore, this was accommodated in the solution procedures. A complete true ab-inito rate theory of kinematics and kinetics for continuum and curved thin structures, without any restriction on the magnitude of the strains or the deformation, was formulated. The time dependence and large strain behavior are incorporated through the introduction of the time rates of the metric and curvature in two coordinate systems, a fixed (spatial) one and a convected (material) coordinate system. The relations between the time derivative and the covariant derivatives (gradients) have been developed for curved space and motion, so that the velocity components supply the connection between the equations of motion and the time rate of change of the metric and curvature tensors.

Simitses, G. J.↗
NASA NTRS1987

Analysis of shell-type structures subjected to time-dependent mechanical and thermal loading

A general mathematical model and solution methodologies for analyzing structural response of thin, metallic shell-type structures under large transient, cyclic, or static thermomechanical loads are developed. Among the system responses, which are associated with these load conditions, are thermal buckling, creep buckling and ratcheting. Thus, geometric as well as material type nonlinearities (of high order) can be anticipated and must be considered in the development of the mathematical model. Furthermore, this must also be accommodated in the solution procedures.

Simitses, G. J.↗