Search NASA⌕ Search

SEARCH · Search NASA

Results for “material response simulation”

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 145 records · Page 8

Numerical Simulation of Rocket Exhaust Interaction with Lunar Soil

This technology development originated from the need to assess the debris threat resulting from soil material erosion induced by landing spacecraft rocket plume impingement on extraterrestrial planetary surfaces. The impact of soil debris was observed to be highly detrimental during NASA s Apollo lunar missions and will pose a threat for any future landings on the Moon, Mars, and other exploration targets. The innovation developed under this program provides a simulation tool that combines modeling of the diverse disciplines of rocket plume impingement gas dynamics, granular soil material liberation, and soil debris particle kinetics into one unified simulation system. The Unified Flow Solver (UFS) developed by CFDRC enabled the efficient, seamless simulation of mixed continuum and rarefied rocket plume flow utilizing a novel direct numerical simulation technique of the Boltzmann gas dynamics equation. The characteristics of the soil granular material response and modeling of the erosion and liberation processes were enabled through novel first principle-based granular mechanics models developed by the University of Florida specifically for the highly irregularly shaped and cohesive lunar regolith material. These tools were integrated into a unique simulation system that accounts for all relevant physics aspects: (1) Modeling of spacecraft rocket plume impingement flow under lunar vacuum environment resulting in a mixed continuum and rarefied flow; (2) Modeling of lunar soil characteristics to capture soil-specific effects of particle size and shape composition, soil layer cohesion and granular flow physics; and (3) Accurate tracking of soil-borne debris particles beginning with aerodynamically driven motion inside the plume to purely ballistic motion in lunar far field conditions. In the earlier project phase of this innovation, the capabilities of the UFS for mixed continuum and rarefied flow situations were validated and demonstrated for lunar lander rocket plume flow impingement under lunar vacuum conditions. Applications and improvements to the granular flow simulation tools contributed by the University of Florida were tested against Earth environment experimental results. Requirements for developing, validating, and demonstrating this solution environment were clearly identified, and an effective second phase execution plan was devised. In this phase, the physics models were refined and fully integrated into a production-oriented simulation tool set. Three-dimensional simulations of Apollo Lunar Excursion Module (LEM) and Altair landers (including full-scale lander geometry) established the practical applicability of the UFS simulation approach and its advanced performance level for large-scale realistic problems.

Liever, Peter↗

Motion Imagery and Robotics Application (MIRA): Standards-Based Robotics

This technology development originated from the need to assess the debris threat resulting from soil material erosion induced by landing spacecraft rocket plume impingement on extraterrestrial planetary surfaces. The impact of soil debris was observed to be highly detrimental during NASA s Apollo lunar missions and will pose a threat for any future landings on the Moon, Mars, and other exploration targets. The innovation developed under this program provides a simulation tool that combines modeling of the diverse disciplines of rocket plume impingement gas dynamics, granular soil material liberation, and soil debris particle kinetics into one unified simulation system. The Unified Flow Solver (UFS) developed by CFDRC enabled the efficient, seamless simulation of mixed continuum and rarefied rocket plume flow utilizing a novel direct numerical simulation technique of the Boltzmann gas dynamics equation. The characteristics of the soil granular material response and modeling of the erosion and liberation processes were enabled through novel first principle-based granular mechanics models developed by the University of Florida specifically for the highly irregularly shaped and cohesive lunar regolith material. These tools were integrated into a unique simulation system that accounts for all relevant physics aspects: (1) Modeling of spacecraft rocket plume impingement flow under lunar vacuum environment resulting in a mixed continuum and rarefied flow; (2) Modeling of lunar soil characteristics to capture soil-specific effects of particle size and shape composition, soil layer cohesion and granular flow physics; and (3) Accurate tracking of soil-borne debris particles beginning with aerodynamically driven motion inside the plume to purely ballistic motion in lunar far field conditions.

Martinez, Lindolfo↗

Modeling the Proton Radiation Belt With Van Allen Probes Relativistic Electron-Proton Telescope Data

An empirical model of the proton radiation belt is constructed from data taken during 2013-2017 by the Relativistic Electron-Proton Telescopes on the Van Allen Probes satellites. The model intensity is a function of time, kinetic energy in the range 18-600 megaelectronvolts, equatorial pitch angle, and L shell of proton guiding centers. Data are selected, on the basis of energy deposits in each of the nine silicon detectors, to reduce background caused by hard proton energy spectra at low L. Instrument response functions are computed by Monte Carlo integration, using simulated proton paths through a simplified structural model, to account for energy loss in shielding material for protons outside the nominal field of view. Overlap of energy channels, their wide angular response, and changing satellite orientation require the model dependencies on all three independent variables be determined simultaneously. This is done by least squares minimization with a customized steepest descent algorithm. Model uncertainty accounts for statistical data error and systematic error in the simulated instrument response. A proton energy spectrum is also computed from data taken during the 8 January 2014 solar event, to illustrate methods for the simpler case of an isotropic and homogeneous model distribution. Radiation belt and solar proton results are compared to intensities computed with a simplified, on-axis response that can provide a good approximation under limited circumstances.

Selesnick, R. S.↗

Optimal fabrication processes for unidirectional metal-matrix composites: A computational simulation

A method is proposed for optimizing the fabrication process of unidirectional metal matrix composites. The temperature and pressure histories are optimized such that the residual microstresses of the composite at the end of the fabrication process are minimized and the material integrity throughout the process is ensured. The response of the composite during the fabrication is simulated based on a nonlinear micromechanics theory. The optimal fabrication problem is formulated and solved with non-linear programming. Application cases regarding the optimization of the fabrication cool-down phases of unidirectional ultra-high modulus graphite/copper and silicon carbide/titanium composites are presented.

Saravanos, D. A.↗

Optimal fabrication processes for unidirectional metal-matrix composites - A computational simulation

A method is proposed for optimizing the fabrication process of unidirectional metal matrix composites. The temperature and pressure histories are optimized such that the residual microstresses of the composite at the end of the fabrication process are minimized and the material integrity throughout the process is ensured. The response of the composite during the fabrication is simulated based on a nonlinear micromechanics theory. The optimal fabrication problem is formulated and solved with nonlinear programming. Application cases regarding the optimization of the fabrication cool-down phases of unidirectional ultra-high modulus graphite/copper and silicon carbide/titanium composites are presented.

Saravanos, D. A.↗

Quasi-Static and Dynamic Analysis of Composite Panels: Characterization and Validation

The predictive capability of a progressive damage model to simulate the quasi-static and dynamic crushing of composite specimens is evaluated in this article. The material model is commercially available within the ABAQUS finite element software package, i.e., “Damage for Fiber reinforced Composites”, and is often utilized to simulate damage progression and failure in static and dynamic failure scenarios such as crushing of a composite plate. The composite specimens used in this study consist of quasi-isotropic C-shaped and corrugated composite panels comprised of carbon fiber/epoxy braided fabric. The C-shaped panel specimen is used to calibrate the model and optimize material/model properties through comparison with experimental results. Then, using the same model parameters, the quasi-static crush response of the corrugated panel is predicted. The predictive capability of the model is further demonstrated by simulating the dynamic crushing response of both C-shaped and corrugated panels using the same material properties as in the quasi-static crush case with only minor adjustments for the dynamic elastic properties (i.e., rate dependent modulus). Results show that the in-built material model available in ABAQUS can successfully reproduce experimental results for both panels (C-shaped and corrugated) subjected to both quasi-static and dynamic loading scenarios. A small-scale parametric study on physical (experimentally measurable) and nonphysical (purely mathematical) model parameters was also conducted. It can be concluded that, achieving successful simulation results requires a more in-depth understanding of the influence of model parameter variation on the mechanical response of the composite together with the strategies and challenges of the utilized modeling methodology (e.g., FEA mesh density, etc.).

Naghipour Ghezeljeh, Paria↗

Simulating the Impact Response of Composite Airframe Components

In 2010, NASA Langley Research Center obtained residual hardware from the US Army's Survivable Affordable Repairable Airframe Program (SARAP). The hardware consisted of a composite fuselage section that was representative of the center section of a Black Hawk helicopter. The section was fabricated by Sikorsky Aircraft Corporation and designated the Test Validation Article (TVA). The TVA was subjected to a vertical drop test in 2008 to evaluate a tilting roof concept to limit the intrusion of overhead mass items, such as the rotor transmission, into the fuselage cabin. As a result of the 2008 test, damage to the hardware was limited primarily to the roof. Consequently, when the post-test article was obtained in 2010, the roof area was removed and the remaining structure was cut into six different types of test specimens including: (1) tension and compression coupons for material property characterization, (2) I-beam sections, (3) T-sections, (4) cruciform sections, (5) a large subfloor section, and (6) a forward framed fuselage section. In 2011, NASA and Sikorsky entered into a cooperative research agreement to study the impact responses of composite airframe structures and to evaluate the capabilities of the explicit transient dynamic finite element code, LS-DYNA®, to simulate these responses including damage initiation and progressive failure. Finite element models of the composite specimens were developed and impact simulations were performed. The properties of the composite material were represented using both a progressive in-plane damage model (Mat 54) and a continuum damage mechanics model (Mat 58) in LS-DYNA. This paper provides test-analysis comparisons of time history responses and the location and type of damage for representative I-beam, T-section, and cruciform section components.

Jackson, Karen E.↗

Lunar powder simulator under lunarlike conditions - Dielectric properties

A powered basalt simulating lunar regolith material was subjected to temperature variations similar to those found in the lunar surface during a lunation. To obtain the dielectric response under simulated lunar conditions, all measurements were made in appropriate vacuums which, in addition, allowed for an evaluation of moisture effects in the powdered sample. It is found that during approximately seven tenths of a lunation the dielectric response of a 5- to 10-cm surface layer should be fairly constant. During the remaining three tenths, it should undergo rapid changes. Changes in the dielectric properties of the surface layer should be controlled by temperature. In the subjacent layers such changes are thought to be controlled by the increasing density of the regolith.

Alvarez, R.↗

Cyclic arc plasma tests of RSI materials using a preheater

The results of a test program are reported in which a preheater was used with an arc plasma stream to study the thermal response of samples of candidate reusable surface insulation materials for the space shuttle. The preheater simulated the shuttle temperature history during the first and last portions of the test cycle, which could not be simulated by the air arc plasma flow. Pre- and post-test data taken for each of the materials included magnified views, optical properties, and chemical analyses. The test results indicate that the mullite base samples experience higher surface temperatures than the other materials at heating rates greater than 225 kw/sq m. The ceramic fibrous mullite and silica coatings show noncatalytic wall behavior. Internal temperature response data for the materials are compared and correlated with analytical predictions.

Stewart, D. A.↗

Development and Characterization of a Rate-Dependent Three-Dimensional Macroscopic Plasticity Model Suitable for Use in Composite Impact Problems

Several key capabilities have been identified by the aerospace community as lacking in the material/models for composite materials currently available within commercial transient dynamic finite element codes such as LS-DYNA. Some of the specific desired features that have been identified include the incorporation of both plasticity and damage within the material model, the capability of using the material model to analyze the response of both three-dimensional solid elements and two dimensional shell elements, and the ability to simulate the response of composites composed with a variety of composite architectures, including laminates, weaves and braids. In addition, a need has been expressed to have a material model that utilizes tabulated experimentally based input to define the evolution of plasticity and damage as opposed to utilizing discrete input parameters (such as modulus and strength) and analytical functions based on curve fitting. To begin to address these needs, an orthotropic macroscopic plasticity based model suitable for implementation within LS-DYNA has been developed. Specifically, the Tsai-Wu composite failure model has been generalized and extended to a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The coefficients in the yield function are determined based on tabulated stress-strain curves in the various normal and shear directions, along with selected off-axis curves. Incorporating rate dependence into the yield function is achieved by using a series of tabluated input curves, each at a different constant strain rate. The non-associative flow-rule is used to compute the evolution of the effective plastic strain. Systematic procedures have been developed to determine the values of the various coefficients in the yield function and the flow rule based on the tabulated input data. An algorithm based on the radial return method has been developed to facilitate the numerical implementation of the material model. The presented paper will present in detail the development of the orthotropic plasticity model and the procedures used to obtain the required material parameters. Methods in which a combination of actual testing and selective numerical testing can be combined to yield the appropriate input data for the model will be described. A specific laminated polymer matrix composite will be examined to demonstrate the application of the model.

Impact↗

The Development of a CO2 Test Capability in the NASA JSC ARCJet for Mars Reentry Simulation

The Atmospheric Reentry Materials and Structures Evaluation Facility (ARMSEF) located at NASA Johnson Space Center is used for simulating the extreme environment experienced upon reentry for the development and certification of thermal protection systems (TPS). The facility supports a large variety of programs and was heavily leveraged for the certification and operational support of the TPS for the Orbiter and, more recently, the development of the heat shield for CEV. This paper will provide more detail into the heritage of the facility. Unique attributes of the facility include a modular aerodynamically stabilized arc heater and independently controlled O2 and N2 for the test gases. When combining the O2 and N2 in a 23:77 mass ratio respectively the Earth s atmosphere is accurately simulated and via modification of this ratio the investigation of the effects of atomic oxygen on a material s response is possible. In the summer of 2010 a development effort was started to add CO2 as a third independently controlled test gas such that, when combined with N2, opens up the possibility of accurately simulating a Martian reentry environment. This paper will discuss the test facility, especially the arc heater, in more detail. Initial testing involved relatively low concentrations of CO2 combined with N2 for the primary purpose of gathering data to answer two pressing safety concerns. The first being the rate of production of carbon monoxide (CO) within the ejector vacuum system. The main concern was that CO can be flammable and possibly explosive at high enough concentrations and pressures. The hazard control during the development phase involved the use of injecting N2 inside the test chamber diffuser to dilute and reduce the concentration of any and all CO present. A residual gas analyzer (RGA) was used to determine the relative amount of CO in the exhaust gas and provide a conversion rate of CO2 to CO. This paper will discuss in more detail the results of the RGA data and the calculated conversion rate. The second safety concern addressed is the possible formation of hydrogen cyanide (HCN) and cyanide (CN). HCN would primarily be present in the cooling water while the CN would most probably condense onto the interior surfaces of the test chamber. Water samples and wipes of the test chamber surfaces were analyzed by an industrial hygienist for the presence of HCN and CN. His paper will discuss these results in more detail. Throughout this development effort measurements of the CO2:N2 flowfield were made with heat flux and pressure probes and with laser induced fluorescence (LIF) of the atomic oxygen. This paper will discuss these results

DelPapa, Steven V.↗

Significance of DSMC Computed Aerothermal Environments in the Rarefied Regime for Atmospheric Entry Material Response

During Mars atmospheric entry, the Mars Science Laboratory (MSL) was protected by a 4.5 meters diameter ablative heatshield assembled in 113 tiles. The heatshield was made of NASA's flagship ablative material, the Phenolic Impregnated Carbon Ablator (PICA). Prior work compared the traditional one-dimensional and three-dimensional material response models at different locations in the heatshield. It was observed that the flow was basically one-dimensional in the nose and flank regions, but three-dimensional flow effects were observed in the outer flank. The objective of this work is to study the effects of the aerothermal environment on the material response. We extend prior work by computing aerothermal environments using the direct simulation Monte Carlo (DSMC) code SPARTA and the CFD code Data Parallel Line Relaxation (DPLR). SPARTA is used to compute environment in the rarefied regime prior to 48.4s of entry where the Knudsen number is such that the Navier-Stokes equations can be inaccurate. Similarly to previous work, the DPLR software is used to compute the hypersonic environment for laminar then turbulent boundary layer assumptions from 48.4 s up to 100 s after Entry Interface (EI) along the MSL 08-TPS-02/01a trajectory. We observe that extending the aerothermal environments to times prior to 48.4 s modifies the thermal response of the heat shield at the surface and in-depth; however the effects on the recession are minimal. Additionally, using the assumption of a turbulent boundary layer versus a laminar one leads to higher surface and in-depth temperatures, larger recession, and a displacement of the peak heating and peak recession location.

Aerothermal↗

Simulating Damage Due to a Lightning Strike Event: Effects of Temperature Dependent Properties on Interlaminar Damage

A multidirectional, carbon fiber-epoxy, composite panel is subjected to a simulated lightning strike, within a finite element method framework, and the effect of material properties on the failure (delamination) response is investigated through a detailed numerical study. The numerical model of the composite panel consists of individual homogenized plies with user-defined, cohesive interface elements between them. Lightning strikes are simulated as an assumed combination of excessive heat and high pressure loadings. It is observed that the initiation and propagation of lightning-induced delamination is a significant function of the temperature dependency of interfacial fracture toughness. This dependency must be defined properly in order to achieve reliable predictions of the present lightning-induced delamination in the composite panel.

computational mechanics↗

An Investigation of High-Cycle Fatigue Models for Metallic Structures Exhibiting Snap-Through Response

A study is undertaken to develop a methodology for determining the suitability of various high-cycle fatigue models for metallic structures subjected to combined thermal-acoustic loadings. Two features of this problem differentiate it from the fatigue of structures subject to acoustic loading alone. Potentially large mean stresses associated with the thermally pre- and post-buckled states require models capable of handling those conditions. Snap-through motion between multiple post-buckled equilibrium positions introduces very high alternating stress. The thermal-acoustic time history response of a clamped aluminum beam structure with geometric and material nonlinearities is determined via numerical simulation. A cumulative damage model is employed using a rainflow cycle counting scheme and fatigue estimates are made for 2024-T3 aluminum using various non-zero mean fatigue models, including Walker, Morrow, Morrow with true fracture strength, and MMPDS. A baseline zero-mean model is additionally considered. It is shown that for this material, the Walker model produces the most conservative fatigue estimates when the stress response has a tensile mean introduced by geometric nonlinearity, but remains in the linear elastic range. However, when the loading level is sufficiently high to produce plasticity, the response becomes more fully reversed and the baseline, Morrow, and Morrow with true fracture strength models produce the most conservative fatigue estimates.

Przekop, Adam↗

Quasi-Static and Dynamic Analysis of Composite Panels: Challenges in Engineering-Level Commercial Software

In engineering applications, where analysis time and efficiency are essential, analytical models and commercially available finite element (FE) material models are the two mainstream methods used to analyze composite crush scenarios. A reliable numerical model addressing the crushing response of composite materials requires a nonlinear progressive damage analysis which includes appropriate failure mechanisms such as matrix cracking, fragmentation, delamination, and fiber breakage. One might be able to address all these interacting failure responses through implementation/validation of a user defined material model or by utilization of a commercially available model, the latter being the most efficient engineering approach. In this paper, the quasi-static and dynamic crushing response of composite specimens, consisting of a C-channel and a corrugated geometry manufactured with carbon fiber/epoxy matrix fabric, are modeled using the progressive composite damage model within ABAQUS. The novel part of this work, which has not been previously addressed in the literature, is the utilization of a detailed parametric study to optimize composite damage model parameters for a given quasi-static crush simulation in ABAQUS. C-channel experiments were used to calibrate the model and then refine/optimize a set of material parameters that best fit the experimental results. These identical parameters were then used to simulate the crush response of the corrugated channel specimens, thereby making it a pure prediction and not a double recalibration. Next, the dynamic response of both panels were predicted, further confirming the predictive capability of the developed numerical model. The effect of parameter variation on the crush response is fully documented for future use in engineering applications utilizing this specific ABAQUS damage model.

Paria Naghipour↗

Initial comparison of SSPM ground test results and flight data to NASCAP simulations

The satellite surface potential monitor (SSPM) has been developed for the P78-2 SCATHA (Spacecraft Charging At The High Altitudes) satellite to determine the response of selected spacecraft materials to charged-particle environmental fluxes. Since the monitor infers total surface voltages from a single point on the interior side of insulators, a ground simulation program was undertaken to develop analytical techniques to model the monitors and to obtain an experimental calibration of the relationship between the flight measurement techniques and actual measurements. The experimental testing was conducted using monoenergetic electron beams irradiating samples in the dark. An analytical computer model was developed in the NASCAP (NASA Charging Analyzer Program) code. The analytical model material properties for Kapton that controlled backscatter and secondary yield were adjusted to obtain a single set of values that produced reasonable fits for both voltages and currents. The analytical techniques developed in the ground technology investigation have been applied to space flight conditions. Predictions were compared with limited flight data. The agreement is very good indicating that the technique, and NASCAP, can be used to predict spacecraft material charging behavior. Details of the testing, the analytical modelling technique and flight data comparisons are presented.

Stevens, N. J.↗

Numerical simulation of experiments in the Giant Planet Facility

Utilizing a series of existing computer codes, ablation experiments in the Giant Planet Facility are numerically simulated. Of primary importance is the simulation of the low Mach number shock layer that envelops the test model. The RASLE shock-layer code, used in the Jupiter entry probe heat-shield design, is adapted to the experimental conditions. RASLE predictions for radiative and convective heat fluxes are in good agreement with calorimeter measurements. In simulating carbonaceous ablation experiments, the RASLE code is coupled directly with the CMA material response code. For the graphite models, predicted and measured recessions agree very well. Predicted recession for the carbon phenolic models is 50% higher than that measured. This is the first time codes used for the Jupiter probe design have been compared with experiments.

Green, M. J.↗

Probabilistic structural analysis methods and applications

An advanced algorithm for simulating the probabilistic distribution of structural responses due to statistical uncertainties in loads, geometry, material properties, and boundary conditions is reported. The method effectively combines an advanced algorithm for calculating probability levels for multivariate problems (fast probability integration) together with a general-purpose finite-element code for stress, vibration, and buckling analysis. Application is made to a space propulsion system turbine blade for which the geometry and material properties are treated as random variables.

Cruse, T. A.↗