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 343 records · Page 19

A predictive model for failure properties of thermoset resins

A predictive model for the three-dimensional failure behavior of engineering polymers has been developed in a recent NASA-sponsored research program. This model acknowledges the underlying molecular deformation mechanisms and thus accounts for the effects of different chemical compositions, crosslink density, functionality of the curing agent, etc., on the complete nonlinear stress-strain response including yield. The material parameters required by the model can be determined from test-tube quantities of a new resin in only a few days. Thus, we can obtain a first-order prediction of the applicability of a new resin for an advanced aerospace application without synthesizing the large quantities of material needed for failure testing. This technology will effect order-of-magnitude reductions in the time and expense required to develop new engineering polymers.

Caruthers, James M.↗

Comprehensive Micromechanics-Analysis Code - Version 4.0

Version 4.0 of the Micromechanics Analysis Code With Generalized Method of Cells (MAC/GMC) has been developed as an improved means of computational simulation of advanced composite materials. The previous version of MAC/GMC was described in "Comprehensive Micromechanics-Analysis Code" (LEW-16870), NASA Tech Briefs, Vol. 24, No. 6 (June 2000), page 38. To recapitulate: MAC/GMC is a computer program that predicts the elastic and inelastic thermomechanical responses of continuous and discontinuous composite materials with arbitrary internal microstructures and reinforcement shapes. The predictive capability of MAC/GMC rests on a model known as the generalized method of cells (GMC) - a continuum-based model of micromechanics that provides closed-form expressions for the macroscopic response of a composite material in terms of the properties, sizes, shapes, and responses of the individual constituents or phases that make up the material. Enhancements in version 4.0 include a capability for modeling thermomechanically and electromagnetically coupled ("smart") materials; a more-accurate (high-fidelity) version of the GMC; a capability to simulate discontinuous plies within a laminate; additional constitutive models of materials; expanded yield-surface-analysis capabilities; and expanded failure-analysis and life-prediction capabilities on both the microscopic and macroscopic scales.

Arnold, S. M.↗

nn-PINNs: Non-Newtonian physics-informed neural networks for complex fluid modeling

Time- and rate-dependent material functions in non-Newtonian fluids in response to different deformation fields pose a challenge in integrating different constitutive models into conventional computational fluid dynamic platforms. Considering their relevance in many industrial and natural settings alike, robust data-driven frameworks that enable accurate modeling of these complex fluids are of great interest. The main goal is to solve the coupled Partial Differential Equations (PDEs) consisting of the constitutive equations that relate the shear stress to the deformation and fully capture the behavior of the fluid under various flow protocols with different boundary conditions. Here, in this work, we present non-Newtonian physics-informed neural networks (nn-PINNs) for solving systems of coupled PDEs adopted for complex fluid flow modeling. The proposed nn-PINN method is employed to solve the constitutive models in conjunction with conservation of mass and momentum by benefiting from Automatic Differentiation (AD) in neural networks, hence avoiding the mesh generation step. nn-PINNs are tested for a number of different complex fluids with different constitutive models and for several flow protocols. These include a range of Generalized Newtonian Fluid (GNF) empirical constitutive models, as well as some phenomenological models with memory effects and thixotropic timescales. nn-PINNs are found to obtain the correct solution of complex fluids in spatiotemporal domains with good accuracy compared to the ground truth solution. We also present applications of nn-PINNs for complex fluid modeling problems with unknown boundary conditions on the surface, and show that our approach can successfully recover the velocity and stress fields across the domain, including the boundaries, given some sparse velocity measurements.

42 ENGINEERING↗

Probabilistic Mesomechanical Fatigue Model

A probabilistic mesomechanical fatigue life model is proposed to link the microstructural material heterogeneities to the statistical scatter in the macrostructural response. The macrostructure is modeled as an ensemble of microelements. Cracks nucleation within the microelements and grow from the microelements to final fracture. Variations of the microelement properties are defined using statistical parameters. A micromechanical slip band decohesion model is used to determine the crack nucleation life and size. A crack tip opening displacement model is used to determine the small crack growth life and size. Paris law is used to determine the long crack growth life. The models are combined in a Monte Carlo simulation to determine the statistical distribution of total fatigue life for the macrostructure. The modeled response is compared to trends in experimental observations from the literature.

Tryon, Robert G.↗

Computation of Effective Mechanical Properties and Mechanical Erosion Modeling of TPS Materials

The goal of this presentation is to provide a general overview of the multi-scale modeling formulation to determine if there is additional surface recession in Thermal Protection Systems (TPS) materials as a result of mechanical erosion due to high shear conditions during atmospheric entry. This modeling process is performed at different scales by leveraging two computational frameworks developed at NASA: the Porous Microstructure Analysis (PuMA) software, and the Porous material Analysis Toolbox based on OpenFOAM (PATO). PuMA specializes in computing effective macro-scale material properties by performing material response simulations on 3D digital micro-scale representations of porous micro-structures. The modeling of TPS materials at the micro-scale is essential to understand how they behave as part of a heat shield assembly. The first part of the presentation will detail the implementation of PuMA’s cell-centered finite volume elasticity solver, which allows the computation of macro-scale effective mechanical properties of heterogeneous and anisotropic materials such as fibrous and woven TPS composites. These homogenized mechanical properties are used by PATO’s mechanical erosion model to predict the TPS material’s recession at a larger scale. This work will also provide some examples of multi-scale analysis from the fiber level up to the unit cell. The second part of the presentation will focus on the macro-scale approach to determine if erosion at the heat shield’s surface occurs due to mechanical and thermal loads experienced during atmospheric entry. To accomplish this, a solid mechanics module was integrated within PATO enabling it to model the potential mechanical erosion in three steps: first, after obtaining the effective mechanical properties with PuMA, the implemented stress analysis solver computes the stress and the displacement fields for the TPS material using the wall shear stress tensor, computed using a CFD solver, as boundary conditions; then, regions on the surface where the stress meets the failure criteria are identified; finally, the failed material is removed and the mesh is redistributed accordingly. The outcome is a model capable of predicting the total recession in the material due to surface chemistry and mechanical erosion.

Mechanical Properties↗

The effect of gravity induced free convection upon the melting phenomena of a finite paraffin slab for thermal control

Two theoretical models were developed to predict the thermal response of the phase change material to a given hot plate temperature. A two-dimensional pure conduction model was developed to predict the melting of the phase change material when heat transfer was a function of conduction. A combined conduction-convection model, also two-dimensional, was developed to predict the phase change phenomena when heat transfer was a function of conduction and gravity-induced free convection. Both models were solved using explicit finite difference approximations on a digital computer. The experimental equipment consisted of a rectangular cell utilizing a heat chamber, an expansion chamber, and a test chamber; a sixteen channel multipoint recorder, and a fluid flow system. The recorder monitered hot and cold plate temperatures and interior node temperatures at two second intervals. A comparison of theoretical temperature profiles and experimental temperature profiles is presented for six runs at various angles of inclination of the test cell with respect to the horizontal direction. A detailed discussion of results is presented.

Bain, R. L.↗

Micromechanics Analysis Code Post-Processing (MACPOST) User Guide

As advanced composite materials have gained wider usage. the need for analytical models and computer codes to predict the thermomechanical deformation response of these materials has increased significantly. Recently, a micromechanics technique called the generalized method of cells (GMC) has been developed, which has the capability to fulfill this -oal. Tc provide a framework for GMC, the Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC) has been developed. As MAC/GMC has been updated, significant improvements have been made to the post-processing capabilities of the code. Through the MACPOST program, which operates directly within the MSC/PATRAN graphical pre- and post-processing package, a direct link between the analysis capabilities of MAC/GMC and the post-processing capabilities of MSC/PATRAN has been established. MACPOST has simplified the production, printing. and exportation of results for unit cells analyzed by MAC/GMC. MACPOST allows different micro-level quantities to be plotted quickly and easily in contour plots. In addition, meaningful data for X-Y plots can be examined. MACPOST thus serves as an important analysis and visualization tool for the macro- and micro-level data generated by MAC/GMC. This report serves as the user's manual for the MACPOST program.

Goldberg, Robert K.↗

Impact Response of Thermoplastic Composites – Experiments and Modeling

This study investigates the impact response of thermoplastic composite materials. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK reinforced with T700G (T700/LMPAEK ). Several High Energy Dynamic Impact (HEDI) tests were performed wherein an aluminum projectile impacted T700/LMPAEK panels at different velocities. In the impact experiment studied herein, the projectile impacted the T700/LMAPEK panel at 73.2 m/s and rebounded. The ensuing damage to the T700/LMPAEK panel was characterized via ultrasonic C-scans, visualizing damage modes including matrix cracking and delaminations. Since the projectile rebounded, the fiber breakage damage mode was not evident. The impact event was modeled using a combined continuum damage mechanics and cohesive zone analysis approach with the commercially available hydrocode LS-DYNA®. The dynamic deformation and damage of the continuum, i.e., the plies in the T700/LMPAEK laminate, was modeled using a Deformation Gradient Decomposition (DGD) material model, i.e., MAT299 available within LS-DYNA. The interlaminar regions were treated as zero-thickness cohesive zones, wherein the dynamic delamination mechanics was modeled using a traction-separation formulation. Mixed-mode delamination growth was captured using the Benzegaggh-Kenane (B-K) traction-separation interpolation law. Strain-rate sensitivity of interlaminar fracture toughness was included using approximations based on published literature. The outcomes of the modeling effort are discussed herein which demonstrate good correlations with test data regarding panel deflection, panel-projectile interactions, damage modes and damage extent.

Impact modeling↗

Bayesian calibration of strength model parameters from Taylor impact data

Materials strength plays a key role in determining the mechanical response of engineered structures. As such, accurate strength models are crucial in simulations involving complex loading conditions, particularly when deformation in the plastic regime is deemed important. In this work, a Gaussian process based surrogate for a finite element simulation of a Taylor impact test is developed and used for Bayesian calibration of the Preston–Tonks–Wallace strength model. Here, the surrogate model is shown to closely approximate the salient features of the Taylor cylinder deformation and is validated against simulation output before being used in the strength model calibration routine. The results show that Taylor impact test data can be used in the calibration of constitutive equations through the use of a combination of data science techniques, namely Gaussian processes and Bayesian inference.

36 MATERIALS SCIENCE↗

Inelastic peridynamic model for molecular crystal particles

The peridynamic theory of solid mechanics is applied to modeling the deformation and fracture of micrometer-sized particles made of organic crystalline material. A new peridynamic material model is proposed to reproduce the elastic–plastic response, creep, and fracture that are observed in experiments. The model is implemented in a three-dimensional, meshless Lagrangian simulation code. In the small deformation, elastic regime, the model agrees well with classical Hertzian contact analysis for a sphere compressed between rigid plates. Under higher load, material and geometrical nonlinearity is predicted, leading to fracture. Finally, the material parameters for the energetic material CL-20 are evaluated from nanoindentation test data on the cyclic compression and failure of micrometer-sized grains.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Directed Energy Deposition Process Modeling, Validation, and Process-Informed Optimization

The directed energy deposition (DED) process, one of the most popular additive manufacturing techniques in use today, involves various complex physical mechanisms that are not yet well understood. In this regard, computational tools show promise for elucidating the manufacturing process and enabling nondestructive performance evaluations of manufactured parts. To better control and optimize the DED process?thereby improving the manufactured product? the present work develops and demonstrates a novel artificial intelligence (AI)-based process control and optimization technique. Specifically, a geometry-free thermo-mechanical model with adaptive subdomain construction is developed to accurately capture the material?s thermo-mechanical response under cyclical reheating and high cooling rates [1]. The model is demonstrated and validated with experimental measurements, in light of various geometries and processing parameters. Moreover, based on this thermo-mechanical model, a physics-informed reduced-order model is developed to enable quick predictions of the temperature field at every time step. Furthermore, an AI-based controller is developed that can adapt to the ever-changing system states by automatically adjusting the manufacturing process parameters. This entire work is based on the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) [2] and its recent integration with Libtorch (the C++ frontend of PyTorch [3]). The combined development of the adaptive material deposition scheme, thermo-mechanical model, associated reduced-order model, and AI-based process controller carries great potential for improving advanced manufacturing processes.

36 MATERIALS SCIENCE↗

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-degree 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 International Planetary Probe Workshop (IPPW), NASA Headquarters postulated that a common aeroshell 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-degrees-of-freedom 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 Thermal Protection System (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. Presentation: We will present an overview of the study scope, highlights of the trade studies and design driver analyses, and the final recommendations of a common probe design and assembly. We will also indicate limitations that the common probe design may have for the different destinations. Finally, recommended qualification approaches for missions will be presented.

Hwang, H. H.↗

A phoenics model for the triaxial loading of an initially cylindrical mass of rate-type material with provisions for bulging and yield

This work traces the response of a granular material via the Ten Coefficient Truesdell rate-type constituitive model into the simplest meaningful loading: the triaxial test configuration. A functional relation has been posed for computing the rather peculiar relation between average applied stress and average porosity. Using that relation an attack has been mounted on the dilemma that exists between dynamic and constitutive use of the pressure variable; that is relating dynamic pressure, thermodynamic pressure, stress deviator and higher stress invariants. The resolution was as a linear superposition with a one-way feedback, in that while the dynamic component could not effect the constituitive component, the converse was not true since density appears in the momentum transport relation.

French, K. W., Jr.↗

A constitutive model for glass–ceramic materials

Glass–ceramics have received recent attention for use in glass–ceramic to metal hermetic seals. Due to their heterogeneous microstructure, these materials exhibit a number of advantageous responses over conventional glass based seals. Key amongst them is the possibility of a controllable thermal strain response and apparent coefficient of thermal expansion which may be used to minimize thermally induced residual stresses for aforementioned seals. Furthermore, these behaviors result from an inorganic glass matrix and variety of crystalline ceramic phases including silica polymorph(s) that may undergo reversible solid-to-solid transformations with associated inelastic strain. Correspondingly, these materials exhibit complex thermomechanical responses associated with multiple inelastic mechanisms (viscoelasticity and phase transformation). While modeling these behaviors is essential for developing and analyzing the corresponding applications, no such model exists. Therefore, in this work a three-dimensional continuum constitutive model for glass–ceramic materials combining these various inelastic mechanisms is developed via an internal state variable approach. A corresponding fully implicit three dimensional numerical formulation is also proposed and implemented. The model is used to simulate existing experiments and validate the proposed formalism. As an example, the simple seal problem of a glass–ceramic seal inside a concentric metal shell is explored. Finally, the impact of the cooling rates, viscoelastic shift factors, and inelastic strain on final residual stress state are all investigated and the differing contributions highlighted.

36 MATERIALS SCIENCE↗

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↗

Strain measurements in a rotary engine housing

The development of structural design tools for Rotary Combustion Engines (RCE) using Finite Element Modeling (FEM) requires knowledge about the response of engine materials to various service conditions. This paper describes experimental work that studied housing deformation as a result of thermal, pressure and mechanical loads. The measurement of thermal loads, clamping pressure, and deformation was accomplished by use of high-temperature strain gauges, thermocouples, and a high speed data acquisition system. FEM models for heat transfer stress analysis of the rotor housing will be verified and refined based on these experimental results.

Lee, C. M.↗

Conjugate (solid/fluid) computational fluid dynamics analysis of the space shuttle solid rocket motor nozzle/case and case field joints

Three-dimensional, conjugate (solid/fluid) heat transfer analyses of new designs of the Solid Rocket Motor (SRM) nozzle/case and case field joints are described. The main focus was to predict the consequences of multiple rips (or debonds) in the ambient cure adhesive packed between the nozzle/case joint surfaces and the bond line between the mating field joint surfaces. The models calculate the transient temperature responses of the various materials neighboring postulated flow/leakpaths into, past, and out from the nozzle/case primary O-ring cavity and case field capture O-ring cavity. These results were used to assess if the design was failsafe (i.e., no potential O-ring erosion) and reusable (i.e., no excessive steel temperatures). The models are adaptions and extensions of the general purpose PHOENICS fluid dynamics code. A non-orthogonal coordinate system was employed and 11,592 control cells for the nozzle/case and 20,088 for the case field joints are used with non-uniform distribution. Physical properties of both fluid and solids are temperature dependent. A number of parametric studies were run for both joints with results showing temperature limits for reuse for the steel case on the nozzle joint being exceeded while the steel case temperatures for the field joint were not. O-ring temperatures for the nozzle joint predicted erosion while for the field joint they did not.

Doran, D.↗

A sensitivity analysis of twinning crystal plasticity finite element model using single crystal and poly crystal Zircaloy

The popularity of crystal plasticity finite element method (CPFEM) models is increasing due to their ability to predict the mechanical response of crystalline materials such as metals and metal alloys more accurately than traditional continuum mechanics models. This is since the crystal plasticity models consider the effect of atomic structure, microstructural morphology, and properties of individual grains. These CPFEM models use a large number of material parameters in order to capture the mesoscale physics which comes with the downside of the tedious calibration process. In this paper, a CPFEM code was developed to include the twinning induced grain reorientation and subsequent crystallographic slip for HPC material. The developed code is incorporated in a large-scale, parallelized nonlinear solver WARP3D. Further, a sensitivity analysis with respect to 22 material parameters was then conducted using single crystal and polycrystal representative volume element (RVE) of Zircaloy material. Loading was applied along five different crystallographic orientations for single crystal RVE and along three directions namely, rolling (RD), transverse (TD), and normal (ND) direction for polycrystal RVE. Results obtained from the sensitivity analysis were used for the calibration of material parameters for Zircaloy. Finally, developed code along with calibrated material parameters was used to investigate the effect of the hydride phase formation in Zircaloy which is a typical case observed for nuclear applications. It was found that the volume fraction of the hydride phase has a significant impact on the mechanical properties of Zircaloy.

36 MATERIALS SCIENCE↗