Statistical failure prediction models for brittle materials
Previously cited in issue 13, p. 2107, Accession no. A82-30104
SEARCH · Search NASA
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.
Previously cited in issue 13, p. 2107, Accession no. A82-30104
MCFET computer program, named for combined Monte Carlo Finite Element Technique, implements specialized microstructural-lattice simulation technique. Developed to simulate microstructural evolution in material systems in which modulated phases occur and directionality of modulation influenced by internal and external stresses. Written in FORTRAN.
Micromechanical Combined Stress Analysis (MICSTRAN) computer code provides materials engineers with easy-to-use personal-computer-based software tool to calculate overall properties of composite, given properties of fibers and matrix. Computes overall thermoelastic parameters and stresses by micromechanical analysis. Written in FORTRAN 77.
Radiation from galactic cosmic rays (GCR) and solar particle events (SPE) is a serious hazard to humans and electronic instruments during space travel, particularly on prolonged missions outside the Earth s magnetic fields. Galactic cosmic radiation (GCR) is composed of approx. 98% nucleons and approx. 2% electrons and positrons. Although cosmic ray heavy ions are 1-2% of the fluence, these energetic heavy nuclei (HZE) contribute 50% of the long-term dose. These unusually high specific ionizations pose a significant health hazard acting as carcinogens and also causing microelectronics damage inside spacecraft and high-flying aircraft. These HZE ions are of concern for radiation protection and radiation shielding technology, because gross rearrangements and mutations and deletions in DNA are expected. Calculations have shown that HZE particles have a strong preference for interaction with light nuclei. The best shield for this radiation would be liquid hydrogen, which is totally impractical. For this reason, hydrogen-containing polymers make the most effective practical shields. Shielding is required during missions in Earth orbit and possibly for frequent flying at high altitude because of the broad GCR spectrum and during a passage into deep space and LunarMars habitation because of the protracted exposure encountered on a long space mission. An additional hazard comes from solar particle events (SPEs) which are mostly energetic protons that can produce heavy ion secondaries as well as neutrons in materials. These events occur at unpredictable times and can deliver a potentially lethal dose within several hours to an unshielded human. Radiation protection for humans requires safety in short-term missions and maintaining career exposure limits within acceptable levels on future long-term exploration missions. The selection of shield materials can alter the protection of humans by an order of magnitude. If improperly selected, shielding materials can actually increase radiation damage due to penetration properties and nuclear fragmentation. Protecting space-borne microelectronics from single event upsets (SEUs) by transmitted radiation will benefit system reliability and system design cost by using optimal shield materials. Long-term missions on the surface of the Moon or Mars will require the construction of habitats to protect humans during their stay. One approach to the construction is to make structural materials from lunar or Martian regolith using a polymeric material as a binder. The hydrogen-containing polymers are considerably more effective for radiation protection than the regolith, but the combination minimizes the amount of polymer to be transported. We have made composites of simulated lunar regolith with two different polymers, LaRC-SI, a high-performance polyimide thermoset, and polyethylene, a thermoplastic.
The constitutive properties of a cross-linked single-walled carbon nanotube material are predicted with a multi-scale model. The material is modeled as a transversely isotropic solid using concepts from equivalent-continuum modeling. The elastic constants are determined using molecular dynamics simulation. Some parameters of the molecular force field are determined specifically for the cross-linker from ab initio calculations. A demonstration of how the cross-linked nanotubes may affect the properties of a nanotube/polyimide composite is included using a micromechanical analysis.
In the present study, the assessment and evaluation of various acoustic tile designs were conducted using three-dimensional finite element analysis, which included static analysis, thermal analysis and modal analysis of integral and non-integral tile design options. Various benchmark specimens for acoustic tile designs, including CMC integral T-joint and notched CMC plate, were tested in both room and elevated temperature environment. Various candidate ceramic matrix composite materials were used in the numerical modeling and experimental study. The research effort in this program evolved from numerical modeling and concept design to a combined numerical analysis and experimental study. Many subjects associated with the design and performance of the acoustic tile in jet engine exhaust nozzle have been investigated.
The objectives of this program are the investigation of fundamental approaches to high temperature crack initiation life prediction, identification of specific modeling strategies and the development of specific models for component relevant loading conditions. A survey of the hot section material/coating systems used throughout the gas turbine industry is included. Two material/coating systems will be identified for the program. The material/coating system designated as the base system shall be used throughout Tasks 1-12. The alternate material/coating system will be used only in Task 12 for further evaluation of the models developed on the base material. In Task II, candidate life prediction approaches will be screened based on a set of criteria that includes experience of the approaches within the literature, correlation with isothermal data generated on the base material, and judgements relative to the applicability of the approach for the complex cycles to be considered in the option program. The two most promising approaches will be identified. Task 3 further evaluates the best approach using additional base material fatigue testing including verification tests. Task 4 consists of technical, schedular, financial and all other reporting requirements in accordance with the Reports of Work clause.
This study demonstrates that coupling of a material thermal response code and a flow solver with nonequilibrium gas/surface interaction for simulation of charring carbon ablators can be performed using an implicit approach. The material thermal response code used in this study is the three-dimensional version of Fully Implicit Ablation and Thermal response program, which predicts charring material thermal response and shape change on hypersonic space vehicles. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation method. Coupling between the material response and flow codes is performed by solving the surface mass balance in flow solver and the surface energy balance in material response code. Thus, the material surface recession is predicted in flow code, and the surface temperature and pyrolysis gas injection rate are computed in material response code. It is demonstrated that the time-lagged explicit approach is sufficient for simulations at low surface heating conditions, in which the surface ablation rate is not a strong function of the surface temperature. At elevated surface heating conditions, the implicit approach has to be taken, because the carbon ablation rate becomes a stiff function of the surface temperature, and thus the explicit approach appears to be inappropriate resulting in severe numerical oscillations of predicted surface temperature. Implicit coupling for simulation of arc-jet models is performed, and the predictions are compared with measured data. Implicit coupling for trajectory based simulation of Stardust fore-body heat shield is also conducted. The predicted stagnation point total recession is compared with that predicted using the chemical equilibrium surface assumption
The ability of a material model to capture in-plane matrix mode I and mode II crack growth is an essential component for modeling ply level damage evolution in composite structures. Previous studies using a continuum damage mechanics (CDM) approach have shown success in satisfying benchmark solutions for mode I and II crack growth. However, success was shown using a fiber-aligned meshing strategy, which encourages matrix cracks to propagate in a single band of elements, along the fiber direction. Generating a fiber-aligned mesh becomes a highly involved process for laminates including off-axis (non 0° or 90°) plies. The objective of this study is to quantify the effect of non-fiber aligned mesh discretization on predictions of in-plane matrix crack propagation. The approach taken incrementally varies the mesh orientation angle relative to the fiber orientation; more specifically, misaligned meshes are used to quantify the effect of element angle orientation relative to the initial crack orientation on the energy released during matrix crack propagation simulations using a CDM method. CDM solutions obtained with the misaligned meshes are evaluated against known benchmarks for mode I and II matrix crack growth. The CDM solutions reveal a near-polynomial trend of increased predicted failure stress with increased mesh misalignment angle; hence implying a potential relationship between element orientation angle and apparent fracture toughness.
During the entry of Mars Science Laboratory (MSL) the heatshield was equipped with the instrumentation suite Mars Entry, Descent, and Landing Instruments (MEDLI). In-depth thermocouple (TC) data was used to reconstruct the surface heating[1,2]. Discrepancies between MEDLI’s flight heating data and predicted thermal response and recession[3] led to a campaign to characterize and test Mars2020 and MEDLI2 (Mars Entry, Descent, and Landing Instrument 2)[4] flight-lot thermal protection materials (TPS) at relevant ground test conditions. This paper covers the development and validation of flight-lot material response models for the heatshield material Phenolic Impregnated Carbon Ablator (PICA) and the backshell material SLA-561V (Super Lightweight Ablator). Virgin and char properties such as thermal conductivity and emissivity were updated as a function of temperature and pressure in the MEDLI2 specific material response models. Other material properties were characterized and determined to be similar to prior TPS material measurements used for the heritage models. MEDLI2 also investigated the impact of the protective coating that was applied to MSL and Mars2020’s PICA heatshield to reduce the particulate matter shedding. Evidence from arc jet experiments on coated PICA showed that although there were no adverse effects on the performance of the TPS, there were differences in the in-depth temperatures and reduced recession compared to material response predictions of uncoated PICA.[5,6]. Fully Implicit Ablation and Thermal response (FIAT)[7] material response simulations were conducted with the heritage and MEDLI2 flight-lot material models. FIAT-predicted in-depth temperature responses and recession predictions were compared to flight-lot certification ground-test arc jet test data and MEDLI2 thermocouple data.
The postbuckling failure of square composite plates with central holes is analyzed numerically and experimentally. The particular plates studies have stacking sequences of: (+ and - 45/0/90)(sub 2S); (+ and - 45/0(sub 2))(sub 2S); (+ and - 45/0(sub 6))(sub S); and (+ and - 45)(sub 4S). A simple plate geometry, one with a hole diameter to plate width ratio of 0.3 is compared. Failure load, failure mode, and failure location are predicted numerically by using the finite element method. Predictions are compared with experimental results. In numerical failure analysis the interlaminar shear stresses, as well as the inplane stresses are taken into account. An issue addressed in this study is the possible mode shape change of the plate during loading. It is predicted that the first three laminates fail due to excessive stresses in the fiber direction, and more importantly, that the load level is independent of whether the laminate is deformed in a one-half or two-half wave configuration. It is predicted that the fourth laminate fails due to excessive inplane shear stress. Interlaminar shear failure is not predicted for any laminates. For the first two laminates the experimental observations correlated well with the predictions. Experimentally, the third laminate failed along the side support due to interlaminar shear strength S(sub 23). The fourth experimental laminate failed due to inplane shear in the location predicted, however material softening resulted in a different failure load from predictions.
Recently applications have exposed polymer matrix composite materials to very high strain rate loading conditions, requiring an ability to understand and predict the material behavior under these extreme conditions. In this second paper of a two part report, a three-dimensional composite micromechanical model is described which allows for the analysis of the rate dependent, nonlinear deformation response of a polymer matrix composite. Strain rate dependent inelastic constitutive equations utilized to model the deformation response of a polymer are implemented within the micromechanics method. The deformation response of two representative laminated carbon fiber reinforced composite materials with varying fiber orientation has been predicted using the described technique. The predicted results compare favorably to both experimental values and the response predicted by the Generalized Method of Cells, a well-established micromechanics analysis method.
Materials and devices used on the exterior of spacecraft in low Earth orbit (LEO) are subjected to environmental threats that can cause degradation in material properties, possibly threatening spacecraft mission success. These threats include: atomic oxygen (AO), ultraviolet and x-ray radiation, charged particle radiation, temperature extremes and thermal cycling, micrometeoroid and debris impacts, and contamination. Space environmental threats vary greatly based on spacecraft materials, thicknesses and stress levels, and the mission environment and duration. For more than a decade the Materials International Space Station Experiment (MISSE) has enabled the study of the long duration environmental durability of spacecraft materials in the LEO environment. The overall objective of MISSE is to test the stability and durability of materials and devices in the space environment in order to gain valuable knowledge on the performance of materials in space, as well as to enable lifetime predictions of new materials that may be used in future space flight. MISSE is a series of materials flight experiments, which are attached to the exterior of the International Space Station (ISS). Individual experiments were loaded onto suitcase-like trays, called Passive Experiment Containers (PECs). The PECs were transported to the ISS in the Space Shuttle cargo bay and attached to, and removed from, the ISS during extravehicular activities (EVAs). The PECs were retrieved after one or more years of space exposure and returned to Earth enabling post-flight experiment evaluation. MISSE is a multi-organization project with participants from the National Aeronautics and Space Administration (NASA), the Department of Defense (DoD), industry and academia. MISSE has provided a platform for environmental durability studies for thousands of samples and numerous devices, and it has produced many tangible impacts. Ten PECs (and one smaller tray) have been flown, representing MISSE 1 through MISSE 8, yielding long-duration space environmental performance and durability data that enable material validation, processing recertification and space qualification; improved predictions of materials and component lifetimes in space; model verification and development; and correlation factors between space-exposure and ground-facilities enabling more accurate in-space performance predictions based on ground-laboratory testing. A few of the many experiment results and observations, and their impacts, are provided. Those highlighted include examples on improved understanding of atomic oxygen scattering mechanisms, LEO coating durability results, and polymer erosion yields and their impacts on spacecraft design. The MISSE 2 Atomic Oxygen Scattering Chamber Experiment discovered that the peak flux of scattered AO was determined to be 45 deg from normal incidence, not the model predicted cosine dependence. In addition, the erosion yield (E(sub y)) of Kapton H for AO scattered off oxidized-Al is 22% of the E(sub y) of direct AO impingement. These results were used to help determine the degradation mechanism of a cesium iodide detector within the Hubble Space Telescope Cosmic Origins Spectrograph Experiment. The MISSE 6 Indium Tin Oxide (ITO) Degradation Experiment measured surface electrical resistance of ram and wake ITO coated samples. The data confirmed that ITO is a stable AO protective coating, and the results validated the durability of ITO conductive coatings for solar arrays for the Atmosphere-Space Transition 2 Explorer program. The MISSE 2, 6 and 7 Polymer Experiments have provided LEO AO Ey data on over 120 polymer and composites samples. The flight E(sub y) values were found to range from 3.05 x 10(exp -26) cu cm/atom for the AO resistant polymer CORIN to 9.14 x 10(exp -26) cu cm/atom for polyoxymethylene (POM). In addition, flying the same polymers on different missions has advanced the understanding of the AO E(sub y) dependency on solar exposure for polymers containing fluorine. The MISSE polymer results are highly requested and have impacted spacecraft design for WorldView-2 & -3, the Global Precipitation Measurement-Microwave Imager, and other spacecraft. The flight data has enabled the development of an Atomic Oxygen Erosion Predictive Tool that allows the erosion prediction of new and non-flown polymers. The data has also been used to develop a new NASA Technical Standards Handbook "Spacecraft Polymers Atomic Oxygen Durability Handbook." Many intangible benefits have also been derived from MISSE. For example, over 40 students have collaborated on Glenn's MISSE experiments, which have resulted in greater than $80K in student scholarships and awards in national and international science fairs. Students have also given presentations and won poster competition awards at international space conferences.
This paper describes a broad effort that is aimed at understanding the fundamental mechanisms of crack growth and using that understanding as a basis for designing materials and enabling predictions of fracture in materials and structures that have small characteristic dimensions. This area of research, herein referred to as Damage Science, emphasizes the length scale regimes of the nanoscale and the microscale for which analysis and characterization tools are being developed to predict the formation, propagation, and interaction of fundamental damage mechanisms. Examination of nanoscale processes requires atomistic and discrete dislocation plasticity simulations, while microscale processes can be examined using strain gradient plasticity, crystal plasticity and microstructure modeling methods. Concurrent and sequential multiscale modeling methods are being developed to analytically bridge between these length scales. Experimental methods for characterization and quantification of near-crack tip damage are also being developed. This paper focuses on several new methodologies in these areas and their application to understanding damage processes in polycrystalline metals. On-going and potential applications are also discussed.
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.
The first two years of a two-phase program aimed at improving the high temperature crack initiation life prediction technology for gas turbine hot section components are discussed. In Phase 1 (baseline) effort, low cycle fatigue (LCF) models, using a data base generated for a cast nickel base gas turbine hot section alloy (B1900+Hf), were evaluated for their ability to predict the crack initiation life for relevant creep-fatigue loading conditions and to define data required for determination of model constants. The variables included strain range and rate, mean strain, strain hold times and temperature. None of the models predicted all of the life trends within reasonable data requirements. A Cycle Damage Accumulation (CDA) was therefore developed which follows an exhaustion of material ductility approach. Material ductility is estimated based on observed similarities of deformation structure between fatigue, tensile and creep tests. The cycle damage function is based on total strain range, maximum stress and stress amplitude and includes both time independent and time dependent components. The CDA model accurately predicts all of the trends in creep-fatigue life with loading conditions. In addition, all of the CDA model constants are determinable from rapid cycle, fully reversed fatigue tests and monotonic tensile and/or creep data.
Heat transport programs adaptable for absorption analysis. Lightweight sandwich panel specimen used for comparison of water absorption measurements with program predictions. In program model, moisture -- like heat in heat-transport problem moves through variety of materials and structures along complex paths.
A model to calculate fatigue life is developed based on the assumption that fatigue life is entirely composed of crack growth from an initial microstructural inhomogeneity. Specifically, growth is considered to start from either an ellipsoidal void, a cracked particle, or a debonded particle. The capability of predicting fatigue life from material microstructure is based on linear elastic fracture mechanics principles, the sizes of the crack-initiating microstructural inhomogeneities, and an initiation parameter that is proportional to the cyclic plastic zone size. A key aspect of this modeling approach is that it is linked with a general purpose probability program to analyze the effect of the distribution of controlling microstructural features within the material. This enables prediction of fatigue stress versus life curves for various specimen geometries using distributional statistics obtained from characterizations of the microstructure. Results are compared to experimental fatigue data from an aluminum alloy.