SEARCH · Engineering Papers

Results for “constitutive 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

NASA NTRS2017

Incorporation of Failure Into an Orthotropic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased use in the aerospace and automotive communities. The aerospace community has identified several key capabilities which are currently lacking in the available material models in commercial transient dynamic finite element codes. To attempt to improve the predictive capability of composite impact simulations, a next generation material model is being developed for incorporation within the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is used to allow for the uncoupling of the deformation and damage analyses. In the damage model, a semi-coupled approach is employed where the overall damage in a particular coordinate direction is assumed to be a multiplicative combination of the damage in that direction resulting from the applied loads in various coordinate directions. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure. Failure surfaces can be defined with any arbitrary shape, unlike traditional failure models where the mathematical functions used to define the failure surface impose a specific shape on the failure surface. In the current paper, the complete development of the failure model is described and the generation of a tabulated failure surface for a representative composite material is discussed.

Impact↗

NASA NTRS2017

Implementation of a Tabulated Failure Model Into a Generalized Composite Material Model Suitable for Use in Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased use in the aerospace and automotive communities. The aerospace community has identified several key capabilities which are currently lacking in the available material models in commercial transient dynamic finite element codes. To attempt to improve the predictive capability of composite impact simulations, a next generation material model is being developed for incorporation within the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is used to allow for the uncoupling of the deformation and damage analyses. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure. Failure surfaces can be defined with any arbitrary shape, unlike traditional failure models where the mathematical functions used to define the failure surface impose a specific shape on the failure surface. In the current paper, the complete development of the failure model is described and the generation of a tabulated failure surface for a representative composite material is discussed.

polymer matrix composites↗

NASA NTRS2022

Impact Simulation of Composite Panels for Aerospace Applications

The aerospace research community has been working towards the development of a comprehensive framework for the design of composite structures at the component, subassembly, assembly, and system levels. Such a framework requires accurate, efficient and easy-to-use mathematical and computational models that are in sync with physical reality. The research work in this paper is designed to help an aerospace structural engineer move towards this goal by simulating dynamic impact to composite panels. First, the existing orthotropic visco-elastic visco-plastic material model in LS-DYNA commercial finite element program, MAT_COMPOSITE_TABULATED_PLASTICITY_DAMAGE (MAT_213) is enhanced with the development and use of rate-dependent damage data. Using state-of-the-art experimental procedures, impact validation test data are created. Second, the versatility of the constitutive model is examined by its ability to simulate a low velocity (27 ft/s, 8 m/s) impact test, and two high velocity impact tests (385–417 ft/s, 117–127 m/s) that straddle the ballistic limit of the projectile and the composite panels. Third, these impact events are modeled using three different finite element types – under-integrated and fully-integrated solid elements, and fully-integrated thin shell elements. Insight into accuracy versus numerical efficiency is gained by comparing the performance of solid versus thin shell elements. The numerical results show that accurate predictions of these impact events can be obtained with reasonable effort and care.

Simulations↗

DOE OSTI2023

NEML2 - THE NEW ENGINEERING MATERIAL MODEL LIBRARY, VERSION 2

SF-22-157 NEML2 is an extension of an earlier open source code developed at ANL (SF-xxxx) called NEML. These libraries provide tools for constructing and using models of engineering materials, with a focus on high temperature materials for nuclear reactor applications. The software include tools for setting up material models, calibrating models, storing models for later use in other simulations, and running simple "material point" simulations of the material behavior. The main difference between NEML and NEML2 is the new software has been complete restructured to run on GPUs and other accelerators in order to improve the performance of methods used to train constitutive models against data and to improve the performance of the final models in production environments, like finite element software. While over time we expect NEML2 to become more-or-less feature complete with NEML we are filing a new software disclosure because (1) the library has been complete restructured and (2) we want to host NEML2 in a separate repository.

MESSNER, MARK↗

NASA NTRS2011

Investigating Some Technical Issues on Cohesive Zone Modeling of Fracture

This study investigates some technical issues related to the use of cohesive zone models (CZMs) in modeling fracture processes. These issues include: why cohesive laws of different shapes can produce similar fracture predictions; under what conditions CZM predictions have a high degree of agreement with linear elastic fracture mechanics (LEFM) analysis results; when the shape of cohesive laws becomes important in the fracture predictions; and why the opening profile along the cohesive zone length needs to be accurately predicted. Two cohesive models were used in this study to address these technical issues. They are the linear softening cohesive model and the Dugdale perfectly plastic cohesive model. Each cohesive model constitutes five cohesive laws of different maximum tractions. All cohesive laws have the same cohesive work rate (CWR) which is defined by the area under the traction-separation curve. The effects of the maximum traction on the cohesive zone length and the critical remote applied stress are investigated for both models. For a CZM to predict a fracture load similar to that obtained by an LEFM analysis, the cohesive zone length needs to be much smaller than the crack length, which reflects the small scale yielding condition requirement for LEFM analysis to be valid. For large-scale cohesive zone cases, the predicted critical remote applied stresses depend on the shape of cohesive models used and can significantly deviate from LEFM results. Furthermore, this study also reveals the importance of accurately predicting the cohesive zone profile in determining the critical remote applied load.

Wang, John T.↗

NASA NTRS2012

Dynamic Load Measurement of Ballistic Gelatin Impact Using an Instrumented Tube

Bird strikes are a common problem for the aerospace industry and can cause serious damage to an aircraft. Ballistic gelatin is frequently used as a surrogate for actual bird carcasses in bird strike tests. Numerical simulations of these tests are used to supplement experimental data, therefore it is necessary to use numerical modeling techniques that can accurately capture the dynamic response of ballistic gelatin. An experimental technique is introduced to validate these modeling techniques. A ballistic gelatin projectile is fired into a strike plate attached to a 36 in. long sensor tube. Dynamic load is measured at two locations relative to the strike plate using strain gages configured in a full Wheatstone bridge. Data from these experiments are used to validate a gelatin constitutive model. Simulations of the apparatus are analyzed to investigate its performance.

Seidt, J. D.↗

DOE OSTI2024

A non-isothermal breakage-damage model for plastic-bonded granular materials incorporating temperature, pressure, and rate dependencies

Plastic-bonded granular materials (PBM) are widely used in industrial sectors, including building construction, abrasive applications, and defense applications such as plastic-bonded explosives. The mechanical behavior of PBM is highly nonlinear, irreversible, rate dependent, and temperature sensitive governed by various micromechanical attributions such as grain crushing and binder damage. This paper presents a thermodynamically consistent, microstructure-informed constitutive model to capture these characteristic behaviors of PBM. Key features of the model include a breakage internal variable to upscale the grain-scale information to the continuum level and to predict grain size evolution under mechanical loading. In addition, a damage internal state variable is introduced to account for the damage, deterioration, and debonding of the binder matrix upon loading. Temperature is taken as a fundamental external state variable to handle non-isothermal loading paths. The proposed model is able to capture with good accuracy several important aspects of the mechanical properties of PBM, such as pressure-dependent elasticity, pressure-dependent yield strength, brittle-to-ductile transition, temperature dependency, and rate dependency in the post-yielding regime. Furthermore, the model is validated against multiple published datasets obtained from confined and unconfined compression tests, covering various PBM compositions, confining pressures, temperatures, and strain rates.

Breakage↗

NASA NTRS1990

Inelastic deformation under nonisothermal loading

The objective of this paper is to evaluate, both experimentally and analytically, the appropriate forms of the hardening evolution equations in unified constitutive models for conditions involving nonisothermal loading. Critical experiments were performed for a cast nickel-base superalloy by using variable temperature tensile, creep, and cyclic tests in the 538-982 C temperature range. These experimental results were compared with both isothermal data and predictions of the Bodner-Partom-Bodner, 1987 elastic-viscoplastic theory to assess the effects of thermal history on constitutive behavior. The results indicate that the hardening evolution equations based on isothermal data are applicable for nonisothermal loading of these precipitation strengthened alloys. Additional thermal history effect terms in the hardening evolution equations were not required beyond those accounting for the variation of material constants with temperature. Using material constants determined solely from isothermal data, the inelastic deformation behavior of B1900 + Hf subject to thermomechanical loading were adequately predicted by the Bodner-Partom model.

Chan, K. S.↗

DOE OSTI2025

Probabilistic Predictions for Fastener Failure in the Sandia Mechanics Challenge Using the Discrete-Direct Uncertainty Quantification Approach

This paper documents the blind and post-blind analysis predictions for the 2023 Sandia Mechanics Challenge (SMC), which involved predicting the behavior of a threaded fastener joint structure subjected to shock loading. Utilizing repeat sets of fastener calibration data from various experimental configurations including tension, double shear, and joint tension, we developed a library of calibrated models which were propagated through the application model using the Discrete-Direct (DD) uncertainty quantification (UQ) approach. Although the initial blind predictions did not incorporate spare-sample processing to quantify fastener failure probabilities, the analyses yielded reasonable conclusions aligned with experimental results. In the post-blind analysis phase, we focused on enhancing the fidelity of the aluminum constitutive model and innovating the DD approach to obtain probabilistic predictions for fastener failure, particularly when quantities of interest (QoIs) approach their bounds. The improved aluminum model captures the behavior of the cantilever under shock loading more accurately, predicting both partial and complete cracks, although it tends to underpredict failure propagation. The enhanced DD approach facilitates probabilistic predictions that reflect the interdependent failure mechanisms of the fasteners and the cantilever, revealing that while certain fasteners are more likely to fail, the failure does not necessarily follow a progressive pattern. Overall, the post-blind analyses significantly improved the predictive capabilities of the model, providing valuable insights into the SMC application and establishing a robust foundation for informed engineering decisions. The methodology demonstrates a cost-effective and extensible approach suitable for a wide range of applications, highlighting the importance of uncertainty quantification to provide context for engineering decision making.

42 ENGINEERING↗

NASA NTRS1991

Thermal barrier coating life prediction model development, phase 2

The objective of this program was to generate a life prediction model for electron-beam-physical vapor deposited (EB-PVD) zirconia thermal barrier coating (TBC) on gas turbine engine components. Specific activities involved in development of the EB-PVD life prediction model included measurement of EB-PVD ceramic physical and mechanical properties and adherence strength, measurement of the thermally grown oxide (TGO) growth kinetics, generation of quantitative cyclic thermal spallation life data, and development of a spallation life prediction model. Life data useful for model development was obtained by exposing instrumented, EB-PVD ceramic coated cylindrical specimens in a jet fueled burner rig. Monotonic compression and tensile mechanical tests and physical property tests were conducted to obtain the EB-PVD ceramic behavior required for burner rig specimen analysis. As part of that effort, a nonlinear constitutive model was developed for the EB-PVD ceramic. Spallation failure of the EB-PVD TBC system consistently occurred at the TGO-metal interface. Calculated out-of-plane stresses were a small fraction of that required to statically fail the TGO. Thus, EB-PVD spallation was attributed to the interfacial cracking caused by in-plane TGO strains. Since TGO mechanical properties were not measured in this program, calculation of the burner rig specimen TGO in-plane strains was performed by using alumina properties. A life model based on maximum in-plane TGO tensile mechanical strain and TGO thickness correlated the burner rig specimen EB-PVD ceramic spallation lives within a factor of about plus or minus 2X.

Meier, Susan Manning↗

NASA NTRS2016

Incorporation of Damage and Failure into an Orthotropic Elasto-Plastic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

A material model which incorporates several key capabilities which have been identified by the aerospace community as lacking in the composite impact models currently available in LS-DYNA(Registered Trademark) is under development. In particular, the material model, which is being implemented as MAT 213 into a tailored version of LS-DYNA being jointly developed by the FAA and NASA, incorporates both plasticity and damage within the material model, utilizes experimentally based tabulated input to define the evolution of plasticity and damage as opposed to specifying discrete input parameters (such as modulus and strength), and is able to analyze the response of composites composed with a variety of fiber architectures. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. The capability to account for the rate and temperature dependent deformation response of composites has also been incorporated into the material model. For the damage model, a strain equivalent formulation is utilized to allow for the uncoupling of the deformation and damage analyses. In the damage model, a diagonal damage tensor is defined to account for the directionally dependent variation of damage. However, in composites it has been found that loading in one direction can lead to damage in multiple coordinate directions. To account for this phenomena, the terms in the damage matrix are semi-coupled such that the damage in a particular coordinate direction is a function of the stresses and plastic strains in all of the coordinate directions. The onset of material failure, and thus element deletion, is being developed to be a function of the stresses and plastic strains in the various coordinate directions. Systematic procedures are being developed to generate the required input parameters based on the results of experimental tests.

Finite Element Method↗

NASA NTRS1976

Development of a solid propellant viscoelastic dynamic model

The results of a one year study to develop a dynamic response model for the Space Shuttle Solid Rocket Motor (SRM) propellant are presented. An extensive literature survey was conducted, from which it was concluded that the only significant variables affecting the dynamic response of the SRM propellant are temperature and frequency. Based on this study, and experimental data on propellants related to the SRM propellant, a dynamic constitutive model was developed in the form of a simple power law with temperature incorporated in the form of a modified power law. A computer program was generated which performs a least-squares curve-fit of laboratory data to determine the model parameters and it calculates dynamic moduli at any desired temperature and frequency. Additional studies investigated dynamic scaling laws and the extent of coupling between the SRM propellant and motor cases. It was found, in agreement with other investigations, that the propellant provides all of the mass and damping characteristics whereas the case provides all of the stiffness.

Hufferd, W. L.↗

NASA NTRS1990

Application of thermal life prediction model to high-temperature aerospace alloys B1900+Hf and Haynes 188

The results of the application of a newly proposed thermomechanical fatigue (TMF) life prediction method to a series of laboratory TMF results on two high-temperature aerospace engine alloys are presented. The method, referred to as TMF/TS-SRP, is based on three relatively recent developments: the total strain version of the method of Strainrange Partitioning (TS-SRP), the bithermal testing technique for characterizing TMF behavior, and advanced viscoplastic constitutive models. The high-temperature data reported in a companion publication are used to evaluate the constants in the model and to provide the TMF verification data to check its accuracy. Predicted lives are in agreement with the experimental lives to within a factor of approximately 2.

Halford, Gary R.↗

NASA NTRS1992

Application of a thermal fatigue life prediction model to high-temperature aerospace alloys B1900 + Hf and Haynes 188

The results of the application of a newly proposed thermomechanical fatigue (TMF) life prediction method to a series of laboratory TMF results on two high-temperature aerospace engine alloys are presented. The method, referred to as TMF/TS-SRP, is based on three relatively recent developments: the total strain version of the method of Strain Range Partitioning (TS-SRP), the bithermal testing technique for characterizing TMF behavior, and advanced viscoplastic constitutive models. The high-temperature data reported in a companion publication are used to evaluate the constants in the model and to provide the TMF verification data to check its accuracy. Predicted lives are in agreement with the experimental lives to within a factor of approximately 2.

Halford, G. R.↗

DOE OSTI2026

Calibration of the Nonlinear Viscoelastic SPECtacular Model for 828-CTBN/DEA/GMB and 828/Z/AlOx

Nonlinear viscoelastic model calibrations are developed for two filled epoxies: an adduct of Epon 828 DGEBA and CTBN cross-linked with diethanolamine (DEA) filled with glass microballoons, called 828-CTBN/DEA/GMB, and Epon 828 DGEBA cross-linked with Ancamine Z and filled with alumina, called 828/Z/AlOx. The thermo-mechanical responses of these materials are represented using the SPECtacular constitutive model, an extension of the Simplified Potential Energy Clock (SPEC) model. The model couples

Cundiff, Kenneth Noel [Sandia National Laboratorie↗

NASA NTRS1996

Computational Materials Research

Computational Materials aims to model and predict thermodynamic, mechanical, and transport properties of polymer matrix composites. This workshop, the second coordinated by NASA Langley, reports progress in measurements and modeling at a number of length scales: atomic, molecular, nano, and continuum. Assembled here are presentations on quantum calculations for force field development, molecular mechanics of interfaces, molecular weight effects on mechanical properties, molecular dynamics applied to poling of polymers for electrets, Monte Carlo simulation of aromatic thermoplastics, thermal pressure coefficients of liquids, ultrasonic elastic constants, group additivity predictions, bulk constitutive models, and viscoplasticity characterization.

Hinkley, Jeffrey A.↗

DOE OSTI2020

Modeling The Effects of Loading Scenario and Thermal Expansion Coefficient on Potential Failure of Cryo-compressed Hydrogen Vessels

A multiscale thermomechanical model for a simplified Type-3 cryogenic compressed hydrogen (H2) storage vessel is developed in this paper. The model accounts for the temperature-dependent elastic-plastic behavior of the vessel carbon/epoxy composite overwrap and aluminum alloy liner. The homogenized thermo-elastic-plastic behavior for the individual laminas of the vessel layup is obtained by an incremental Eshelby-Mori-Tanka approach associated with a micromechanical failure criterion to predict lamina failure while a standard elastic-plastic constitutive model is used to describe the behavior of a typical aluminum alloy assumed for the liner. The vessel response to external loadings is achieved by a finite element method. Four loading scenarios representing four thermomechanical cycles applied to the vessel are analyzed to evaluate constituent and lamina stresses as well as the associate failure criterion during the cycle according to these scenarios. The model can provide helpful guidance to mitigate thermal stresses by an adequate selection of loading scenario, optimizing the layup and by tailoring thermomechanical properties of the resin matrix.

36 MATERIALS SCIENCE↗

NASA NTRS2005

Measurement and Prediction of the Thermomechanical Response of Shape Memory Alloy Hybrid Composite Beams

Previous work at NASA Langley Research Center (LaRC) involved fabrication and testing of composite beams with embedded, pre-strained shape memory alloy (SMA) ribbons within the beam structures. That study also provided comparison of experimental results with numerical predictions from a research code making use of a new thermoelastic model for shape memory alloy hybrid composite (SMAHC) structures. The previous work showed qualitative validation of the numerical model. However, deficiencies in the experimental-numerical correlation were noted and hypotheses for the discrepancies were given for further investigation. The goal of this work is to refine the experimental measurement and numerical modeling approaches in order to better understand the discrepancies, improve the correlation between prediction and measurement, and provide rigorous quantitative validation of the numerical analysis/design tool. The experimental investigation is refined by a more thorough test procedure and incorporation of higher fidelity measurements such as infrared thermography and projection moire interferometry. The numerical results are produced by a recently commercialized version of the constitutive model as implemented in ABAQUS and are refined by incorporation of additional measured parameters such as geometric imperfection. Thermal buckling, post-buckling, and random responses to thermal and inertial (base acceleration) loads are studied. The results demonstrate the effectiveness of SMAHC structures in controlling static and dynamic responses by adaptive stiffening. Excellent agreement is achieved between the predicted and measured results of the static and dynamic thermomechanical response, thereby providing quantitative validation of the numerical tool.

Davis, Brian↗