Modeling Fatigue Crack Growth Using Symbolic Regression
Explore the source record for details and available documents.
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.
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.
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.
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.
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.
A novel, in-situ, scanning electron microscope (SEM) mechanical testing capability for materials at the microscale which provides experimental validation to a machine learning (ML) toolset for full-field validation of physics-based micromechanics models is being developed by researchers at NASA Glenn Research Center. These are enabling technologies for the integration of multiscale digital twins for materials into system level models which will result in the improved performance, material discovery, reduced production cost and time, rapid characterization, and prognostic structural health monitoring (SHM) for materials and structures for extreme environments in support of NASA space exploration missions. In order to bridge the material structure-to-system gap for digital twins, physics-based models must be experimentally validated at multiple length scales. Seminal microscale experiments, conducted at the Air Force Research Laboratory (AFRL), were limited to transverse compression of single-layer, unidirectional thermoset polymer matrix composite (PMC) micropillar specimens [1]. The early phases of the current project followed those initial results and setup to reproduce the compression testing of PMC material on the custom-built piezoelectric actuated micromechanical testing rig built by MicroTesting Solutions LLC. In this work, samples of thermoplastic PMC material were first machined into 3 mm cubes, and then further machining and final milling was done using a Focused Ion Beam (FIB). The initial experiment was done on a pillar roughly 20 µm x 20 µm x 40 µm tall. Additional pillars were milled with final sizes ranging from 20 µm x 20 µm x 40 µm tall to 40 µm x 40 µm x 65 µm tall. A speckle pattern for in-situ full-field measurements using Digital Image Correlation (DIC) was applied with platinum, which was coated on the surface, and then the FIB was used to mill away some of the coating to produce an irregular pattern of Pt on the pillar surface. The samples were loaded into the custom testing rig and placed into the SEM and loaded under compression until failure. Images were collected in the SEM during testing. Post-processing of the images was conducted using DIC to obtain full-field displacement and strain measurements elucidating the role of the matrix as well as fiber-fiber interaction at the microscale within the composite subjected to compression loading well into the non-linear regime of the material. Moreover, the evolution of fiber-matrix debonding and matrix cracking is observed in-situ at the microscale. This data, along with images segmented with a newly developed ML toolset [2], was used to create and validate physics-based micromechanics models. An image of the failed micropillar is shown in Figure 1. The techniques developed in the initial compression experiment was tailored to the validation needs of the models and expanded to include different sized samples as well as possibly tension and fatigue.
Zeolites are a class of crystalline aluminosilicate materials that form the backbone of the chemical process industry worldwide. They are used primarily as adsorbents and catalysts and support to a significant extent the positive balance of trade realized by the chemical industry in the United States (around $19 billion in 1991). The magnitude of their efforts can be appreciated when one realizes that since their introduction as 'cracking catalysts' in the early 1960's, they have saved the equivalent of 60 percent of the total oil production from Alaska's North Slope. Thus the performance of zeolite catalysts can have a profound effect on the U.S. economy. It is estimated that a 1 percent increase in yield of the gasoline fraction per barrel of oil would represent a savings of 22 million barrels of crude oil per year, representing a reduction of $400 million in the United States' balance of payments. Thus any activity that results in improvement in zeolite catalyst performance is of significant scientific and industrial interest. In addition, due to their 'stability,' uniformity, and, within limits, their 'engineerable' structures, zeolites are being tested as potential adsorbents to purify gases and liquids at the parts-per-billion levels needed in today's electronic, biomedical, and biotechnology industries and for the environment. Other exotic applications, such as host materials for quantum-confined semiconductor atomic arrays, are also being investigated. Because of the importance of this class of material, extensive efforts have been made to characterize their structures and to understand their nucleation and growth mechanisms, so as to be able to custom-make zeolites for a desired application. To date, both the nucleation mechanics and chemistry (such as what are the 'key' nutrients) are, as yet, still unknown for many, if not all, systems. The problem is compounded because there is usually a 'gel' phase present that is assumed to control the degree of supersaturation, and this gel undergoes a continuous 'polymerization' type reaction during nucleation and growth. Generally, for structure characterization and diffusion studies, which are useful in evaluating zeolites for improving yield in petroleum refining as well as for many of the proposed new applications (e.g., catalytic membranes, molecular electronics, chemical sensors) large zeolites (greater than 100 to 1000 times normal size) with minimum lattice defects are desired. Presently, the lack of understanding of zeolite nucleation and growth precludes the custom design of zeolites for these or other uses. It was hypothesized that the microgravity levels achieved in an orbiting spacecraft could help to isolate the possible effects of natural convection (which affects defect formation) and minimize sedimentation, which occurs since zeolites are twice as dense as the solution from which they are formed. This was expected to promote larger crystals by allowing growing crystals a longer residence time in a high-concentration nutrient field. Thus it was hypothesized that the microgravity environment of Earth orbit would allow the growth of large, more defect-free zeolite crystals in high yield.