Modeling Fatigue Crack Growth Using Symbolic Regression
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Fatigue analysis tools can vary across industries. For example, automotive engineers often use the Wöhler (S-N) approach to design for safe-life, while aerospace engineers prioritize damage tolerance and inspection intervals, relying instead on crack growth models such as Paris’ law. Although both approaches may deal with the control of cracks in similar materials, their analysis tools and material characterizations are fundamentally distinct. This divide mirrors the classic split between stress-based strength analysis and linear elastic fracture mechanics. However, modern nonlinear models that incorporate material softening, such as cohesive laws, blur this boundary and capture fracture behaviors across scales. This presentation describes the CF23 fatigue model, which uses cohesive softening to link S-N crack initiation with crack propagation rates. CF23 spans the full fatigue spectrum, from initial propagation transients to steady-state growth and threshold conditions, offering a unified framework that bridges Wöhler and Paris-based methodologies. Example applications include fatigue crack propagation transients in adhesive interfaces and skin/stiffener separation.
This paper presents recent developments in advanced analysis methods for the computation of stress site damage. The method of solution is based on the p-version of the finite element method. Its implementation was designed to permit extraction of linear stress intensity factors using a superconvergent extraction method (known as the contour integral method) and evaluation of the J-integral following an elastic-plastic analysis. Coarse meshes are adequate for obtaining accurate results supported by p-convergence data. The elastic-plastic analysis is based on the deformation theory of plasticity and the von Mises yield criterion. The model problem consists of an aluminum plate with six equally spaced holes and a crack emanating from each hole. The cracks are of different sizes. The panel is subjected to a remote tensile load. Experimental results are available for the panel. The plasticity analysis provided the same limit load as the experimentally determined load. The results of elastic-plastic analysis were compared with the results of linear elastic analysis in an effort to evaluate how plastic zone sizes influence the crack growth rates. The onset of net-section yielding was determined also. The results show that crack growth rate is accelerated by the presence of adjacent damage, and the critical crack size is shorter when the effects of plasticity are taken into consideration. This work also addresses the effects of alternative stress-strain laws: The elastic-ideally-plastic material model is compared against the Ramberg-Osgood model.
The present report follows on the cohesive fatigue damage model methodology proposed in NASA-TP-2018-219838. In that report, an empirical function describing the incremental damage due to cyclic loading was used to calculate fatigue damage within a cohesive formulation. The form of the function was developed such that, when integrated at a constant stress amplitude from no damage to failure, it produces a life versus load response that is consistent with an S-N diagram. Therefore, the parameters of the damage model could be obtained by fitting the model predictions to an S-N diagram. The finite element analyses performed demonstrate that the cohesive fatigue accumulation function provides a link between the S-N diagram that describes crack initiation, and the Paris law that characterizes the rate of crack propagation. However, when the model was proposed, it was not known whether the form of the damage accumulation function associated with a desired S-N diagram is unique and, if not, if the link between S-N and the Paris law is unique and independent of the fatigue function selected. In the effort described herein, several alternative forms of the damage function that reproduce the desired features of S-N diagrams were found and evaluated. The effects of each of these functions on the predicted parameters of the Paris law and the propagation threshold are discussed. The results indicate that the predicted exponent m of the Paris law is indeed independent of the damage accumulation function. However, different functions predict different values for the pre-factor C of the Paris law. Therefore, the proper damage accumulation function must be selected by comparison with experiments. One of the new damage accumulation functions proposed herein was found to be particularly useful for analysis because of the ease with which the model parameters can be determined with a minimal amount of experimental information. The effectiveness of the proposed methodology and damage function was demonstrated by conducting analyses of a double cantilever beam test, a mixed-mode bending test, and a three-point bending test of a skin/doubler specimen. The results indicate that the same set of model parameters can provide accurate predictions of the rate of fatigue crack propagation for a variety of material interfaces, mode mixities, load levels, and stress ratios.
After the rollout of Space Shuttle Discovery in April 2005 in preparation for return-to-flight, there was a failure of the Orbiter (OV-103) helium signature leak test in the gaseous hydrogen (GH2) system. Leakage was attributed to the Flow Control Valve (FCV) in Main Engine 3. The FCV determined to be the source of the leak for OV-103 is designated as LV-58. The nitrile/Buna N rubber O-ring seal was removed from LV-58, and failure analysis indicated radial cracks providing leak paths in one quadrant. Cracks were eventually found in 6 of 9 FCV O-rings among the three Shuttle Orbiters, though none were as severe as those for LV-58, OV-103. Testing by EM10 at MSFC on all 9 FCV O- rings included: laser dimensional, Shore A hardness and properties from a dynamic mechanical analyzer (DMA) and an Instron tensile machine. The following test data was obtained on the cracked quadrant of the LV-58, OV-103 O-ring: (1) the estimated compression set was only 9.5%, compared to none for the rest of the O-ring; (2) Shore A hardness for the O.D. was higher by almost 4 durometer points than for the rest of the O-ring; and (3) DMA data showed that the storage/elastic modulus E was almost 25% lower than for the rest of the O-ring. Of the 8 FCV O-rings tested on an Instron, 4 yielded tensile strengths that were below the MIL spec requirement of 1350 psi-a likely influence of rubber cracking. Comparisons were made between values of modulus determined by DNA (elastic) and Instron (Young s). Each nitrile/Buna N O-ring used in the FCV conforms to the MIL-P-25732C specification. A number of such O-rings taken from shelf storage at MSFC and Kennedy Space Center (KSC) were used to generate a reference curve of DMA glass transition temperature (Tg) vs. shelf storage time ranging from 8 to 26 years. A similar reference curve of TGA onset temperature (of rubber weight loss) vs. shelf storage time was also generated. The DMA and TGA data for the used FCV O-rings were compared to the reference curves. Correlations were also made between the DMA modulus (at 22 C) and Shore A hardness for all 9 of the FCV O-rings used among the three Shuttle Orbiters. The radial cracking in the FCV O-rings was determined to be due to ozone attack, as nitrile/Buna N rubber is susceptible to such attack. Nitrile/Buna N material under MIL-P25732C should be used in a hydraulic fluid environment to help protect it from cracking. However, the FCV O-rings were used in an air only environment. The FCV design has as much as a 9-mil gap that allows the O.D. of the O-ring to be directly exposed to ozone, pressurized air and some elevated temperatures, accelerating the weathering process that leads to O-ring cracking. Space Shuttle flights will likely not continue past 2010. Therefore, Shuttle management decided to continue using the nitrile/Buna N material for the FCVs, but have each O-ring replaced after 3 years to minimize any chances for crack initiation.
Volume IV of the Report contains appendices that were not cited in Volume I. These consist of documents produced by NASA and other organizations, which were provided to the Columbia Accident Investigation Board in support of its inquiry into the February 1, 2003 destruction of the Space Shuttle Columbia. The documents are compiled in this volume in the interest of establishing a complete record, but they do not necessarily represent the views of the Board. Volume I contains the Board's findings, analysis, and recommendations. The documents in Volume III through V are also contained in their original color format on the DVD disc in the back of Volume II. Volumes II – VI contain a number of conclusions and recommendations, several of which were adopted by the Board in Volume I. The other conclusions and recommendations drawn in Volumes II – VI do not necessarily reflect the opinion of the Board, but are included for the record. When there is conflict, Volume I takes precedence.
Following the successful Artemis I skip reentry on December 11, 2022, unexpected char liberation from the Orion heat- shield prompted formation of an Anomaly Response Team to determine root cause and establish flight rationale and corrective actions for subsequent missions. Through a comprehensive investigation including detailed hardware analysis, modeling & simulation, and ground testing, the team determined that the low permeability Avcoat experienced extreme internal gas pressure buildup during the skip entry that could not adequately outgas, leading to crack formation and char liberation. Key discoveries included significant performance differences between permeable and impermeable heatshield regions and successful replication of the anomaly through ground testing. Based on extensive ground testing and analysis, Artemis II flew a modified trajectory without the skip entry to minimize crack-inducing conditions, and successfully splashed down on April 10, 2026 with significantly reduced char loss. For Artemis III and beyond, a corrective action was implemented to produce a more permeable version of Avcoat that will substantially reduce internal pressure accumulation during entry.
The present study investigates the mixed mode interlaminar fracture toughness of six polymer matrix composite materials at both room temperature and cryogenic liquid nitrogen temperature. The single leg bending test employed gives a mix of Mode I (opening) and Mode II (in-plane shearing) cracking. The tested thermoplastic-matrix composites were shown to have total mixed mode toughnesses 3.1–4.3 times higher than those of representative thermoset resin composites.
This project, jointly sponsored by Rocketdyne and CSTAR, involves the development of laser joining of materials which have heretofore been impractical to bond. Of particular interest are joints between stainless steel and copper and also aluminum 6061 to aluminum 2219. CSTAR has a unique opportunity in this area since both the process and development and diagnostics are of interest to industry. Initial results using the pulse tailored laser welding technique developed in CLA for joining crack sensitive materials have proven promising for the aluminum joints based upon metallurgical and electronic microprobe analysis. A declaration of success requires additional mechanical testing. A CW technique has been applied to the stainless-copper joining with some preliminary success. These joints are of significant interest for aeronautics and rocket propulsion applications and the project is expected to continue.
A 350 micron thick 7 wt.% Y 2 O 3 -ZrO 2 (7YSZ) ceramic Thermal Barrier Coating (TBC) was manufactured in an argon atmosphere and with a strictly controlled substrate temperature on flat IN 600 samples and on IN 100 aircraft turbine blades. This new atmosphere and temperature controlled spraying (ATCS) allows the reduction of the residual stresses and an improvement in the microstructure and in the mechanical properties of the coating. The good performance of these new TBC's was assessed in thermal fatigue and in thermal shock tests carried out in air at 1100°C and at 1280°C respectively. The oxidation rate of the metal-ceramic interface was measured at 1100°C. Failure of the TBC's was induced only by the deterioration (like oxidation, cracking, ...) of the uncoated parts of the samples, as testing temperature was too severe for the base metals. Spalling was thereafter driven by compressive stresses near the aluminum oxide layer at the metal-ceramic interface. The results of these tests underline the good quality of the ceramic coating manufactured by the new ATCS technology.
Atmospheric Plasma Spray (APS) and Vacuum Plasma Spray (VPS) techniques were used to develop Ti/2 vol.% hBN coatings, for extreme space environments and tested aboard the International Space Station as part of the MISSE-17 (Materials International Space Station Experiments) program. The coatings were exposed to atomic oxygen, space radiation, and low-orbit thermal cycling. VPS coatings showed a 56% increase in microhardness, a 26% rise in elastic modulus, minimal porosity and crack density changes compared to APS coatings. The change in mechanical properties is attributed to the formation of TiO, TiO₂ and TiN from nitrogen retention, alongside radiation-induced dislocations, which enhanced surface hardening. The oxidation of titanium led to the formation of TiO and TiO₂, while boron nitride was retained and underwent transmutation in VPS coatings. XPS and EDS analyses confirmed the enhanced space-environment resistance of VPS coatings, making them ideal for long-term spacecraft protection in lunar and Martian conditions.
A novel benzoxazine made from furfurylamine, paraformaldehyde, and piceol, Bz-FA-HA, is assessed for applications in fiber reinforced (FR) composites. Piceol is a biobased phenolic compound derived from the roots of Norwegian spruce trees and contains a methyl ketone group at the para position. Bz-FA-HA is a liquid at room temperature, has a viscosity of < 1 Pa.s at temperatures above 90 °C, and a Tonset of cure at 132 °C. The carbonyl is found to react with the furan ring, yielding a crosslinking reaction, when cured above 180 °C as indicated by differential scanning calorimetry and thermogravimetric analysis coupled with Fourier transform infrared spectroscopy. Poly(Bz-FA-HA) has a Tg > 350 °C, attributed to the crosslinking reaction. Furthermore, the storage modulus is > 3 GPa, regardless of cure temperature. Poly(Bz-FA-HA) has a char yield at 800 °C of 65.0 % (62.3 % at 1000 °C), and a Tonset of decomposition of 357 °C in nitrogen. The resulting carbon formed during pyrolysis shrinks during the carbonization reaction and scanning electron microscopy imaging shows a cross-section with micro cracks. The high processability, advanced mechanical properties, and exceptional char yield make it a promising candidate as the matrix for FR composites.
The National Aeronautics and Space Administration (NASA) is visiting the moon again. This time, the objective is to explore establishing a permanent lunar base. To achieve both longterm human habitation on the moon and future deep space travel, it is crucial to make the most of the in-situ resources and build autonomous systems on the moon to support the construction of a lunar habitat. NASA’s In-situ Resource Utilization (ISRU) program aims to minimize the need to ship heavy prefabricated structures, reducing cost and enhancing sustainability. Here, we developed an economical extrusion method for printing lunar regolith-based composites using shape memory vitrimer as a binder. A rheological study is conducted to determine the extrudability of the composite with different regolith weight percentages. Several characterizations were conducted on the composites. The as-printed composites exhibited compressive and flexural strengths of 73.32 MPa and 156.59 MPa, respectively, and good impact tolerance. The composite maintained 57.92% of its mechanical properties even after the second crack healing cycle. The composites also exhibit shape fixity ratio of 90.02% and shape recovery ratio of 83.46%. The simple synthesis method, sustainability, and good thermomechanical properties make the 3D printed composite an ideal material for lunar construction applications.
Friction stir welded (FSW) aluminum blanks offer larger preform sizes than available from the mill in single width plates. Large blanks made from aerospace-grade aluminum alloys are essential for making large domes and cones for launch vehicles and crew vehicles. Previous research efforts have demonstrated >5m FSW blanks can be formed into large domes or cones using spin forming, but the subsequent solution heat treatment consistently led to abnormal grain growth (AGG). The AGG phenomenon occurs because the FSW blanks must be spin formed in the -O annealed temper to withstand the high strain process without cracking. Subsequently, the solution heat treatment process leads to the onset of AGG, whereby a small number of grains grow much larger than the surrounding matrix. Ultimately, the large grains and bimodal grain-size distribution leads to unacceptable ductility in the short-transverse direction of the part. Alternatively, explosive hydroforming is a nascent process that can be applied to large and thick aerospace-grade aluminum alloy blanks in the -T3 temper and only require a post-forming artificial aging process to a -T8. Avoiding the solution heat treatment (SHT) process mitigates AGG, and also eliminates the need for very rare drop-furnace quench baths that can accommodate a part size larger than 5m in diameter. This presentation summarizes a Lockheed Martin-led effort that demonstrated explosive hydroforming of a 1.9m diameter, 39mm thick 2050 aluminum lithium FSW blank. Details will be provided on the process flow, mechanical properties, and metallographic analyses. This work was performed in partnership with NASA Langley Research Center (LaRC) under a NASA Announcement of Collaboration Opportunities program.
In January of 1999 the NM DS-2 Mars microprobe will be launched to impact on Mars in December. The technical objectives of the missions are to demonstrate: key technologies, a passive atmospheric entry, highly integrated microelectronics which can withstand both low temperatures and high decelerations, and the capability to conduct in-situ, surface and subsurface science data acquisition. The scientific objectives are to determine if ice is present below the Martian surface, measure the local atmospheric pressure, characterize the thermal properties of the martian subsurface soil, and to estimate the vertical temperature gradient of the Martian soil. The battery requirements are 2-4 cell batteries, with voltage of 6-14 volts, capacity of 550 mAh at 80C, and 2Ah at 25C, shelf life of 2.5 years, an operating temperature of 60C and below, and the ability to withstand shock impact of 80,000 g's. The technical challenges and the approach is reviewed. The Li-SOCL2 system is reviewed, and graphs showing the current and voltage is displayed, along with the voltage over discharge time. The problems encountered during the testing were: (1) impact sensitivity, (2) cracking of the seals, and (3) delay in voltage. A new design resulted in no problems in the impact testing phase. The corrective actions for the seal problems involved: (1) pre weld fill tube, (2) an improved heat sink during case to cover weld and (3) change the seal dimensions to reduce stress. To correct the voltage delay problem the solutions involved: (1) drying the electrodes to reduce contamination by water, (2) assemblage of the cells within a week of electrode manufacture, (3) ensure electrolyte purity, and (4) provide second depassivation pulse after landing. The conclusions on further testing were that the battery can: (1) withstand anticipated shock of up to 80,000 g, (2) meet the discharge profile post shock at Mars temperatures, (3) meet the required self discharge rate and (4) meet environmental requirements.
The objective of this research work has been to provide analytical background and support to the ongoing experimental program at NASA, White Sands Test Facility, involving testing composite overwrapped pressure vessels (COPV) for impact damage and cyclic pressurization. Preliminary theoretical basis, including the governing equations for a shallow shell subjected to internal pressure, has been established. Effects of the Griffith type cracks on the structural integrity of the cylindrical vessel were evaluated by methods of Fracture Mechanics. The results indicate that the effective mass of the pressure vessel is an important factor influencing the response to impact events. We also have found that the material properties of the target, contained in the constitutive equations of the composite attached to the Aluminum liner, dominate the impact event in the low velocity range, the material properties become less important, while the target mass distribution and the impactor mass become more significant as the velocity of the impactor increases. Therefore, at high-velocity impact it is not only the kinetic energy of the impactor but also its mass which has a significant effect on the dynamics of the event, and consequently on the induced damage. This work also suggests a methodology for an assessment of the rate of loading effects on the degradation of the material toughness associated with a high-velocity impact where the rate effects become significant. To model the rate dependence of the material response a viscoelastic-plastic constitutive equations were assumed, and on this basis predictions are made regarding the rate dependent material resistance curve. Other dynamic phenomena associated with the impact event have been treated in the framework of the Computational Mechanics using the courtesy of Prof. P. Guebelle and his graduate student at University of Illinois at Urbana-Champaign who have an access to a super-fast computer located on their campus. Finally, the guidelines for a follow-up research program are provided in the body of this report. They address three major areas: theoretical research, numerical studies, and further experimental work.
Additive manufacturing has proven useful for introducing oxide particles into superalloys, to provide oxide dispersion strengthening (ODS) at high temperatures. This study screened how such an approach can influence microstructure and failure modes for superalloy 718. The objective of this study was to compare tensile and creep failure modes for additively manufactured superalloy 718 having introduced oxide particles. Fine yttrium-oxide (yttria) particles were introduced into 718 powder using two different powder mixing methods. Additive manufacturing by laser powder bed fusion was then used to prepare specimen blanks for each case, oriented parallel and transverse to the building direction. These were subsequently given consistent stress relief, hot isostatic pressurization, and final heat treatments, then subjected to tensile and creep tests at varied temperatures. Additively manufactured 718 without ODS (“no-ODS”) had nearly equiaxed grains at 50 μm in width. Roll-mixed ODS material had a bimodal grain size distribution, with fine grains at 9 μm and coarse grains at 66 μm in average width. The fine grains were elongated about 5x in the building direction. Acoustic-mixed ODS materials had grains at 10 μm in width that were elongated about 10x in the build direction. Compared to no-ODS, roll-mixed and acoustic mixed ODS materials did not show improved tensile strengths at room temperature, 760 °C, and 1093 °C. The tensile failure strains parallel to the building direction of roll-mixed ODS material were lowest in these tests. Flattened clumps of yttria had formed transverse to the building direction in roll-mixed material, which were easily cracked. Creep rupture response at 760 °C in the building direction was highest for acoustic-mixed material test. Creep tests at 1093 °C showed varied results, with low rupture lives associated with enhanced cavitation at grain boundaries and surface oxidation. Transverse to the building direction, all three materials had low failure strains in tests at 760 °C. The elongated grain boundaries failed readily in all transverse specimens tested at 760 °C and 1093 °C.