Search NASA⌕ Search

SEARCH · Search NASA

Results for “Material Response Modeling”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Thermomechanical characterization of Hastelloy-X under uniaxial cyclic loading

In most high-temperature engineering applications, components are subjected to complex combinations of thermal and mechanical loading during service. A number of viscoplastic constitutive models were proposed which potentially can provide mathematical descriptions of material response under such conditions. Implementation of these models into large finite element codes such as MARC has already resulted in much improved inelastic analysis capability for hot-section aircraft engine components. However, a number of questions remain regarding the validity of methods adopted in characterizing these constitutive models for particular high-temperature materials. One area of concern is that the majority of experimental data available for this purpose are determined under isothermal conditions. This is in contrast to service conditions which, as noted above, almost always involve some form of thermal cycling. The obvious question arises as to whether a constitutive model characterized using an isothermal data base can adequately predict material response under thermomechanical conditions. An experimental program was initiated within the HOST program to address this particular concern. The results of the most recent isothermal and thermomechanical experiments are described.

Ellis, J. R.↗
NASA NTRSDate not supplied

Identification and Study of Validation Level Test Cases for Computational Modeling of Non-Charring Ablators

Computational modeling of Thermal Protection System (TPS) materials, used for aerospace applications, provides numerous advantages in preliminary selection and design of a heatshield material and shape for atmospheric entry vehicles. However, to serve as a reliable tool for prediction of material thermal and ablative behavior, the modeling approach needs to be validated against real experimental and flight data, preferably at a range of applied conditions. The validation study is typically very complex as it requires reliable measured data not only for the material thermal response and surface recession, but also well characterized environmental conditions. The validation problem becomes even more complex when the material thermal response is dictated by multi-physics effects such as solid conduction, in-depth thermal decomposition, pyrolysis gas flow and chemical reactions. The multi-physics effects complicate not only the modeling effort, but also the experimental measurement for validation of various aspects of the highly coupled problem. In this study, an attempt is made to identify suitable experimental data that could serve as a source for validation of material thermal response modeling tools. To reduce the computational complexity, this study focuses only on non-charring ablators, where the material thermal response could be modeled with a single governing equation for solid conduction and the ablation is limited only to the surface of the material. With a well characterized and publicly available experimental data being sparse, the study is limited in presenting test cases for only three materials: camphor, graphite and FiberForm® in the sequence of increased modeling complexity. Graphite is a commonly used TPS material for aerospace applications, both for leading edges of high-speed vehicles and internal insulation of solid rocket motors. FiberForm® is a porous carbon pre-form used in preparation of the well known PICA material Tran et al. [1996]. Inclusion of camphor into the list is conditioned with the relative simplicity in modeling the material thermal and chemical response and the low-enthalpy flow environment. In addition, camphor has been used as a simple test material for study of flow transition behavior by Stock and Ginoux [1973] and assessment of a heatshield shape change at flight relevant conditions by Rotondi et al. [2022]. In this work, the identified experimental data was extracted from the public literature and test cases that yet have been published. As it was found from the review, not a single test case contains an exhaustive set of data that would validate every aspect of the material physics. However, in the data collected, various aspects of the material behavior can be still validated, such as surface and in-depth temperature, amount of recession and a shape change. The identified experimental data for each case is accompanied with a characterized flow environment and simulated boundary conditions predicted by a Data-Parallel Line Relaxation (DPLR) code Wright et al. [1998]. In addition, material thermal response numerical simulations in each test case are performed with Kentucky Aerothermodynamics and Thermal Response System (KATS-MR) Zibitsker et al. [2022] providing a comparative study and a sanity check for the proposed validation data. Sample results from the performed numerical study are shown below. Figure 1 shows distribution of surface heat flux and pressure values on a hemi-cylinder model made of FiberForm® and tested in HyMETS arc-jet facility. The results are shown for the high pressure condition among the two tests. Flow simulation was performed with DPLR code on a quarter of original geometry. In the figure, the quarter shape was mirrored across zx and xy planes to show the complete distribution. Figure 2 shows the material response results for the high pressure case (7500 Pa), simulated with KATS-MR and a comparison to the experimental data for the surface temperature and shape shape. The simulation was performed on a 2-D slice, extracted in the xy plane at the middle of the sample. Figure 3 shows the material response simulation for the low pressure case (3500 Pa) and a comparison to the experimental data for surface temperature and shape change.

ablation↗
NASA NTRS1977

Some fundamental aspects of the fatigue and fracture response of composite materials

A simple analytical model, based on the concept that the constraining plies surrounding cracks through the thickness of a given ply or plies control the spacing of those cracks by controlling the surrounding stress state, appears to agree well with experimental data. Such a model can be used to evaluate critically various constraint situations including the prediction of the ability of a given constraint layer to arrest flaw development. The analytical model, which uses an equilibrium element approach in the present case, is sensitive to lamina and laminate elastic properties, flaw size, stacking sequence, and an interlaminar stress transfer distance which is estimated from experimental observations. The analysis is easily extended to study flaw growth and plastic region formation.

Reifsnider, K. L.↗
DOE OSTI2025

High-Fidelity Modeling of Fuel-To-Coolant Thermomechanical Transport Behaviors Under Transient Conditions

This report summarizes the work completed under NEUP project number 21-24006. The objectives of this project are to advance the high-fidelity modeling capabilities and important phenomena that is important for high-burnup UO 2 and accident tolerant fuels (ATF) during transient conditions. Accurate modeling of the time-dependent phenomena that impact material performance must be used to determine the figures of merit and safety margin. Phenomena such as fuel fragmentation, cladding oxidation, pellet-clad interaction, clad ballooning, and clad rupture are examples that pose challenges to modeling during these transients. This project focused on the development of high-fidelity tightly coupled multiphysics models that can capture the time-dependent material response and associated thermal hydraulic conditions during these events. These models can then be validated against existing separate effects tests and in-pile integral experiments and will be used to model Transient Reactor Test facility (TREAT) loss-of-coolant accidents (LOCA) experiments. To achieve the project objective, we used a combination of NEAMS and NRC codes to model various LOCA test sets for the separate effects and in-pile integral experiments. BlueCRAB tool set, which can accurately predict material response at a sub-fuel pin level, as well as modeling the entire reactor system response to these events. Fuel performance was modeled using BISON (various versions) and FAST (version 1.2.1). BISON and FAST can model on a sub-fuel pin level the fuel performance under transient conditions. Both have simplified thermal hydraulic models that are capable of providing basic coolant boundary conditions. To better capture the thermomechanical interaction between the fuel, clad, and coolant, more sophisticated thermal hydraulic models are necessary.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Constitutive Model Development for Aging Polymer Encapsulants (ASC P&EM FY2021 L2 Milestone 7836)

This SAND report fulfills the completion requirements for the ASC Physics and Engineering Modeling Level 2 Milestone 7836 during Fiscal Year 2021. The Sandia Simplified potential energy clock (SPEC) non-linear viscoelastic constitutive model was developed to predict a whole host of polymer glass physical behaviors in order to provide a tool to assess the effects of stress on these materials over their lifecycle. Polymer glasses are used extensively in applications such as electronics packaging, where encapsulants and adhesives can be critical to device performance. In this work, the focus is on assessing the performance of the model in predicting material evolution associated with long-term physical aging, an area that the model has not been fully vetted in. These predictions are key to utilizing models to help demonstrate electronics packaging component reliability over decades long service lives, a task that is very costly and time consuming to execute experimentally. The initiating hypothesis for the work was that a model calibration process can be defined that enables confidence in physical aging predictions under ND relevant environments and timescales without sacrificing other predictive capabilities. To test the hypothesis, an extensive suite of calibration and aging data was assembled from a combination of prior work and collaborating projects (Aging and Lifetimes as well as the DoD Joint Munitions Program) for two mission relevant epoxy encapsulants, 828DGEBA/DEA and 828DGEBA/T403. Multiple model calibration processes were developed and evaluated against the entire set of data for each material. A qualitative assessment of each calibration's ability to predict the wide range of aging responses was key to ranking the calibrations against each other. During this evaluation, predictions that were identified as non-physical, i.e., demonstrated something that was qualitatively different than known material behavior, were heavily weighted against the calibration performance. Thus, unphysical predictions for one aspect of aging response could generate a lower overall rating for a calibration process even if that process generated better quantitative predictions for another aspect of aging response. This insurance that all predictions are qualitatively correct is important to the overall aim of utilizing the model to predict residual stress evolution, which will depend on the interplay amongst the different material aging responses. The DSC-focused calibration procedure generated the best all-around aging predictions for both materials, demonstrating material models that can qualitatively predict the whole host of different physical aging responses that have been measured. This step forward in predictive capability comes from an unanticipated source, utilization of calorimetry measurements to specify model parameters. The DSC-focused calibration technique performed better than compression-focused techniques that more heavily weigh measurements more closely related to the structural responses to be predicted. Indeed, the DSC-focused calibration procedure was only possible due to recent incorporation of the enthalpy and heat capacity features into SPEC that was newly verified during this L2 milestone. Fundamentally similar aspects of the two material model calibrations as well as parametric studies to assess sensitives of the aging predictions are discussed within the report. A perspective on the next steps to the overall goal of residual stress evolution predictions under stockpile conditions closes the report.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modular Linear Thermoviscoelastic Model

Time-dependent, viscoelastic responses of materials like polymers and glasses have long been studied. As such, a variety of models have been put forth to describe the behavior including simple rheological models (e.g. Maxwell, Kelvin), linear "fading memory" theories, and hereditary integral based linear thermal viscoelastic approaches as well as more recent nonlinear theories that are either integral, fictive temperature, or differential internal state variable based. The current work details a new LINEAR_THERMOVISCOELASTIC model that has been added to LAME. This formulation represents a viscoelastic theory that neglects some of the phenomenological details of the PEC/SPEC models in favor of efficiency and simplicity. Furthermore, this new model is a first step towards developing "modular" viscoelastic capabilities akin to those available with hardening descriptions for plasticity models in LAME. Specifically, multiple different (including user-defined) shift-factor forms are implemented with each being easily selected via parameter specification rather than requiring distinct material models.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Multi-Phase Modeling Framework Suitable for Dynamic Applications

Under dynamic loading conditions and the associated extreme conditions many metals will undergo phase transformations. The change in crystal structure associated with solid–solid phase transformations can significantly alter the subsequent mechanical response of the material. For the interpretation of experiments involving dynamic loading it is beneficial to have a modeling framework that captures key features of the material response while remaining relatively simple. We introduce a candidate framework and apply it to the metal tin to highlight a range of behaviors that are captured by the model. We also discuss potential extensions to capture additional behaviors that could be important for certain materials and loading scenarios. The model is useful for analysis of results from dynamic experiments and offers a point of departure for more complex model formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Physically realistic, parametric model for excitonic critical point parabolic band oscillators

Critical point parabolic band (CPPB) oscillators are often useful to model the optical response of semiconductor materials, such as hybrid organic–inorganic lead halide-based perovskites, to incident photons in the form of the complex dielectric function (ε=ε1+iε2) spectra. Some models of ε using CPPB oscillators are not guaranteed Kramers–Kronig (KK) consistent (and therefore not physically realistic), may have excess or arbitrary parameter values, or may require prohibitively long computational time when used to fit ellipsometric spectra. For excitonic CPPBs, commonly used to describe the optical response of hybrid organic–inorganic lead halide-based perovskite materials, a physically realistic, parametric model of ε is developed from the KK relationship between ε1 and ε2 for a number of CPPB oscillators with an Urbach tail below the lowest direct transition. This parametric model is shown to produce the same line shape reported from previous works accurately and more quickly than other available KK-consistent CPPB models.

Physics↗
NASA NTRS1980

Erosion of composite materials

A model for describing the response of uncoated and coated fiber reinforced composites subjected to repeated impingements of liquid (rain) droplets is presented. The model is based on the concept that fatigue is the dominant factor in the erosion process. Algebraic expressions are provided which give the incubation period, the rate of mass loss past the incubation period, and the total mass loss of the material during rain impact. The influence of material properties on erosion damage and the protection offered by different coatings are discussed and the use of the model in the design in the design of structures and components is illustrated.

Springer, G. S.↗
NASA NTRS2007

Materials Database Development for Ballistic Impact Modeling

A set of experimental data is being generated under the Fundamental Aeronautics Program Supersonics project to help create and validate accurate computational impact models of jet engine impact events. The data generated will include material property data generated at a range of different strain rates, from 1x10(exp -4)/sec to 5x10(exp 4)/sec, over a range of temperatures. In addition, carefully instrumented ballistic impact tests will be conducted on flat plates and curved structures to provide material and structural response information to help validate the computational models. The material property data and the ballistic impact data will be generated using materials from the same lot, as far as possible. It was found in preliminary testing that the surface finish of test specimens has an effect on measured high strain rate tension response of AL2024. Both the maximum stress and maximum elongation are greater on specimens with a smoother finish. This report gives an overview of the testing that is being conducted and presents results of preliminary testing of the surface finish study.

Pereira, J. Michael↗
NASA NTRS2002

3D Finite Element Analysis of Particle-Reinforced Aluminum

Deformation in particle-reinforced aluminum has been simulated using three distinct types of finite element model: a three-dimensional repeating unit cell, a three-dimensional multi-particle model, and two-dimensional multi-particle models. The repeating unit cell model represents a fictitious periodic cubic array of particles. The 3D multi-particle (3D-MP) model represents randomly placed and oriented particles. The 2D generalized plane strain multi-particle models were obtained from planar sections through the 3D-MP model. These models were used to study the tensile macroscopic stress-strain response and the associated stress and strain distributions in an elastoplastic matrix. The results indicate that the 2D model having a particle area fraction equal to the particle representative volume fraction of the 3D models predicted the same macroscopic stress-strain response as the 3D models. However, there are fluctuations in the particle area fraction in a representative volume element. As expected, predictions from 2D models having different particle area fractions do not agree with predictions from 3D models. More importantly, it was found that the microscopic stress and strain distributions from the 2D models do not agree with those from the 3D-MP model. Specifically, the plastic strain distribution predicted by the 2D model is banded along lines inclined at 45 deg from the loading axis while the 3D model prediction is not. Additionally, the triaxial stress and maximum principal stress distributions predicted by 2D and 3D models do not agree. Thus, it appears necessary to use a multi-particle 3D model to accurately predict material responses that depend on local effects, such as strain-to-failure, fracture toughness, and fatigue life.

Shen, H.↗
NASA NTRS2019

Common Probe Design Study and Follow-On Activities

The Common Probe Study was funded by the NASA's Planetary Science Division in the Science Mission Directorate in 2018 to investigate the feasibility of a common aeroshell design for atmospheric probe missions at Venus, Jupiter, Saturn, Uranus, and Neptune. The study involved 4 NASA Centers: Ames Research Center, Goddard Space Flight Center, Langley Research Center, and the Jet Propulsion Laboratory. The common aeroshell design that was studied was a 400 kg, 1.5 m diameter, 45-degree sphere cone shape with a high density heatshield material (Heatshield for Extreme Entry Environments Technology, or HEEET) and a parachute system to extract the descent vehicle. This size of aeroshell could accommodate a descent vehicle of 0.75 m diameter, which could encompass both Tier 1 and Tier 2 science instruments at each of the 5 destinations. Study methodology: First, a notional payload of instruments for each destination was defined based on the top priority measurements indicated by the Planetary Science Decadal Survey. Steep and shallow entry flight path angles (EFPA) were defined for each planet based on qualification and operational g-load limits for current, state-of-the-art instruments. Interplanetary trajectories were then identified that bounded the EFPA range.Next, 3-DoF simulations for entry trajectories were run using the entry state vectors from the interplanetary trajectories. Conical ribbon parachutes were sized based on heatshield separation dynamics. Aero-heating correlations were used to generate stagnation point convective and radiative heat flux profiles. High fidelity thermal response models for various TPS materials were used to size stagnation point thicknesses, with margins based on previous studies. Backshell TPS masses were assumed based on scaled heat fluxes from the heatshield and also from previous mission concepts.Based on these analyses, we have found that the common design is applicable for atmospheric probe missions for 4 out of the 5 destinations. Because of the unique gravity well for Jupiter, the entry environments are more severe resulting in heat loads an order of magnitude higher than for the other destinations.The next step is to determine what follow-on activities NASA should engage in. A questionnaire for the atmospheric probe community has been developed, with a focus on what size of aeroshell should be further analyzed (smaller or same diameter), and what incentives would make using such an aeroshell, if assembled and available, desirable to mission proposers.Preliminary results from this questionnaire will be presented.

Hwang, Helen H.↗
NASA NTRS2018

A Common Probe Design for Multiple Planetary Destinations

Atmospheric probes have been successfully flown to planets and moons in the solar system to conduct in-situ measurements. They include the Pioneer Venus multi-probes, the Galileo Jupiter probe, and Huygens probe. Probe mission concepts to five destinations, including Venus, Jupiter, Saturn, Uranus, and Neptune, have all utilized similar-shaped aeroshells and concept of operations, namely a 45 deg sphere cone shape with high density heatshield material and parachute system for extracting the descent vehicle from the aeroshell. Each concept designed its probe to meet specific mission requirements and to optimize mass, volume, and cost. At the 2017 IPPW, NASA Headquarters postulated that a common aero-shell design could be used successfully for multiple destinations and missions. This "common probe" design could even be assembled with multiple copies, properly stored, and made available for future NASA missions, potentially realizing savings in cost and schedule and reducing the risk of losing technologies and skills difficult to sustain over decades. Thus the NASA Planetary Science Division funded a study to investigate whether a common probe design could meet most, if not all, mission needs to the five planetary destinations with extreme entry environments. The Common Probe study involved four NASA Centers and addressed these issues, including constraints and inefficiencies that occur in specifying a common design.Study methodology: First, a notional payload of instruments for each destination was defined based on priority measurements from the Planetary Science Decadal Survey. Steep and shallow entry flight path angles (EFPA) were defined for each planet based on qualification and operational g-load limits for current, state-of-the-art instruments. Interplanetary trajectories were then identified for a bounding range of EFPA. Next, 3-DoF simulations for entry trajectories were run using the entry state vectors from the interplanetary trajectories. Aeroheating correlations were used to generate stagnation point convective and radiative heat flux profiles for several aeroshell shapes and entry masses. High fidelity thermal response models for various TPS materials were used to size stagnation point thicknesses, with margins based on previous studies. Backshell TPS masses were assumed based on scaled heat fluxes from the heatshield and also from previous mission concepts.

Hwang, Helen H.↗
DOE OSTI2024

MatCal Users Guide: Release 1.3.0

Any continuum mechanics model will require three components: (1) a discretized geometry of the boundary value problem being studied, (2) the partial differential equations to be solved, and (3) the initial conditions and boundary conditions for the problem. To describe material behavior in these computational models, material models contribute to (2) the underlying equations and, occasionally, to (3) the initial conditions for the simulation. These material models can exhibit a mathematical form that is empirically based, based on first principles, or developed from both empirical observations and known physics. In general, these models are meant to represent a class of materials with well understood behavior. As a result, material models have parameters that must be tuned or calibrated so that the model response matches characterization data available for the specific material it is intended to represent when used to simulate a specific system. For simple models, such as isotropic, linear elastic materials in solid mechanics, this calibration process can be a simple analytical calculation directly extracting the parameters from experimental measurements. For complex models that have many inputs and require many characterization datasets to adequately identify the material behavior, the model calibration process can require an inverse problem approach where an optimization is performed to tune the model parameters to the available data.

36 MATERIALS SCIENCE↗
NASA NTRS2010

Effect of Microscopic Damage Events on Static and Ballistic Impact Strength of Triaxial Braid Composites

The reliability of impact simulations for aircraft components made with triaxial-braided carbon-fiber composites is currently limited by inadequate material property data and lack of validated material models for analysis. Methods to characterize the material properties used in the analytical models from a systematically obtained set of test data are also lacking. A macroscopic finite element based analytical model to analyze the impact response of these materials has been developed. The stiffness and strength properties utilized in the material model are obtained from a set of quasi-static in-plane tension, compression and shear coupon level tests. Full-field optical strain measurement techniques are applied in the testing, and the results are used to help in characterizing the model. The unit cell of the braided composite is modeled as a series of shell elements, where each element is modeled as a laminated composite. The braided architecture can thus be approximated within the analytical model. The transient dynamic finite element code LS-DYNA is utilized to conduct the finite element simulations, and an internal LS-DYNA constitutive model is utilized in the analysis. Methods to obtain the stiffness and strength properties required by the constitutive model from the available test data are developed. Simulations of quasi-static coupon tests and impact tests of a represented braided composite are conducted. Overall, the developed method shows promise, but improvements that are needed in test and analysis methods for better predictive capability are examined.

Littell, Justin D.↗
NASA NTRS1999

Strain Rate Dependent Deformation and Strength Modeling of a Polymer Matrix Composite Utilizing a Micromechanics Approach

Potential gas turbine applications will expose polymer matrix composites to very high strain rate loading conditions, requiring an ability to understand and predict the material behavior under extreme conditions. Specifically, analytical methods designed for these applications must have the capability of properly capturing the strain rate sensitivities and nonlinearities that are present in the material response. The Ramaswamy-Stouffer constitutive equations, originally developed to analyze the viscoplastic deformation of metals, have been modified to simulate the nonlinear deformation response of ductile, crystalline polymers. The constitutive model is characterized and correlated for two representative ductile polymers. Fiberite 977-2 and PEEK, and the computed results correlate well with experimental values. The polymer constitutive equations are implemented in a mechanics of materials based composite micromechanics model to predict the nonlinear, rate dependent deformation response of a composite ply. Uniform stress and uniform strain assumptions are applied to compute the effective stresses of a composite unit cell from the applied strains. The micromechanics equations are successfully verified for two polymer matrix composites. IM7/977-2 and AS4/PEEK. The ultimate strength of a composite ply is predicted with the Hashin failure criteria that were implemented in the composite micromechanics model. The failure stresses of the two composite material systems are accurately predicted for a variety of fiber orientations and strain rates. The composite deformation model is implemented in LS-DYNA, a commercially available transient dynamic explicit finite element code. The matrix constitutive equations are converted into an incremental form, and the model is implemented into LS-DYNA through the use of a user defined material subroutine. The deformation response of a bulk polymer and a polymer matrix composite are predicted by finite element analyses. The results compare reasonably well to experimental values, with some discrepancies. The discrepancies are at least partially caused by the method used to integrate the rate equations in the polymer constitutive model.

Goldberg, Robert K.↗