Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “crack nucleation”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Crack nucleation from non-metallic inclusions in aluminum alloys described by peridynamics simulations

Conventional engineering methods oftentimes have challenges in the quantification of crack nucleation processes from manufacturing defects that are relevant for engineering component lifing. In this work, we present peridynamic simulation framework for the description of the crack nucleation process from cluster of non-metallic inclusions in aluminum alloys. Our non-local simulation framework characterizes crack nucleation process as multiple micro-crack nucleation events from individual inclusions, and eventually one micro-crack dominates. We define individual stages of the crack nucleation process, i.e., nucleation, micro-crack, technical, and crack initiation, that allow a quantification and meta model development of both the individual stages and the entire crack nucleation process.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Crack nucleation at forging flaws studied by non-local peridynamics simulations

We present a computational study and framework that allows us to study and understand the crack nucleation process from forging flaws. Forging flaws may be present in large steel rotor components commonly used for rotating power generation equipment including gas turbines, electrical generators, and steam turbines. The service life of these components is often limited by crack nucleation and subsequent growth from such forging flaws, which frequently exhibit themselves as non-metallic oxide inclusions. The fatigue crack growth process can be described by established engineering fracture mechanics methods. However, the initial crack nucleation process from a forging flaw is challenging for traditional engineering methods to quantify as it depends on the details of the flaw, including flaw morphology. We adopt the peridynamics method to describe and study this crack nucleation process. For a specific industrial gas turbine rotor steel, we present how we integrate and fit commonly known base material property data such as elastic properties, yield strength, and S-N curves, as well as fatigue crack growth data into a peridynamic model. The obtained model is then utilized in a series of high-performance two-dimensional peridynamic simulations to study the crack nucleation process from forging flaws for ambient and elevated temperatures in a rectangular simulation cell specimen. The simulations reveal an initial local nucleation at multiple small oxide inclusions followed by micro-crack propagation, arrest, coalescence, and eventual emergence of a dominant micro-crack that governs the crack nucleation process. The dependence on temperature and density of oxide inclusions of both the details of the microscopic processes and cycles to crack nucleation is also observed. The results are compared with fatigue experiments performed with specimens containing forging flaws of the same rotor steel.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Crack nucleation and growth during dynamic indentation of chemically-strengthened glass

Dynamic point impact loading is a primary cause of fracture of glass screens on mobile devices. An improved understanding of crack initiation and evolution under high-speed indentation could contribute to the development of materials with better performance, but experimental observations are challenging due to the short timescales and limited depth-of-field of optical microscopy. To address this need, we have observed fracture of a chemically-strengthened glass during dynamic indentation using in situ x-ray phase-contrast imaging (XPCI). Median cracks initiate below the surface of the glass, at a depth approximately corresponding to the depth at which the surface residual compressive stress diminishes to zero. These cracks initially propagate at ~10ms –1 , rapidly accelerating to > 100ms –1 . We also observe some evidence for rate-dependent behavior, in that indentation at the lowest rates studied here (below about ~0.15ms –1 ) fails to initiate cracks regardless of the depth of the indentation (up to 18μm), while indentation at higher rates produces median cracks that either arrest or cause complete fracture, depending on the depth of indentation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development of Graphite Thermal and Mechanical Modeling Capabilities in Grizzly

Nuclear-grade graphites are used extensively in the core designs of multiple types of advanced nuclear reactors. In the reactor environment, graphite is exposed over long durations to extreme conditions, including high temperatures, radiation and potentially molten salt and oxygen. Exposure to these conditions can cause several degradation mechanisms in graphite, including nonuniform volumetric strains induced by irradiation and thermal expansion, which lead to stresses that can compromise the performance of graphite components. Evaluating component integrity, predicting component performance over the reactor lifetime, and developing design standards all require robust tools for predicting fracture initiation and propagation in graphite structural components in nuclear reactors. This report documents progress in an ongoing effort to develop modeling and simulation tools in the Grizzly code for predicting the performance of graphite exposed to reactor conditions. Recent developments include a set of thermal and mechanical models that now include the IG-110, NBG-18, and H-451 graphite grades. Improvements have also been made to a nonlinear damaged plasticity model applicable to predicting damage under tension and compression to quasibrittle materials, including graphite. In addition, enhancements have been made to the extended finite element method implementation targeted at simulating graphite fracture. These include new capabilities for crack nucleation in the interior of a solid body, improved treatment of crack nucleation on free surfaces, and more robust modeling approaches for crack growth approaching free surfaces or other cracks.

36 MATERIALS SCIENCE↗
DOE OSTI2019

An anisotropic damage model based on dislocation-mediated nucleation of cracks under high-rate compression

In this report, we developed a thermodynamically-consistent, rate-dependent micromechanics model for brittle damage nucleated by dislocation plasticity applicable for large deformations. Dislocation substructure evolution was used to inform a nucleation criterion for a microcrack. Under global compression, the sliding of a microcrack induces formation of wing cracks. Effective stress drives dynamic growth of these cracks under a 3D stress state, resulting in an anisotropic material stiffness. The model was also advanced to predict grain size dependence of a polycrystalline solid. Internal variables were constrained based on the laws of thermodynamics. Material constants were calibrated for polycrystalline beryllium to demonstrate the applicability of the model to simulate dynamic failure under compression. We demonstrate the versatility of the model to capture brittle to ductile transition governed by temperature and strain rate. The predictive capability of the model to simulate failure stress and failure strain is compared with dynamic and quasistatic data on beryllium.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A phase-field model of quasi-brittle fracture for pressurized cracks: Application to UO 2 high-burnup microstructure fragmentation

Here, in this paper, we present a phase-field model of quasi-brittle fracture with pressurized cracks, with dedicated applications for polycrystalline materials. The model is formulated as a minimization problem within the variational framework. The external work done by pressure on the crack surfaces is included in the objective function. Several careful modeling choices lead to a regularization-length-independent critical strength. The model is constructed to give a softening response with an underlying linear traction-separation law. The pressure-dependent softening response and the regularization of the prescribed pressure are demonstrated with a (quasi) one-dimensional numerical analysis. In a two-dimensional numerical analysis under plane strain assumptions, the critical stress corresponding to crack propagation (as predicted by our quasi-brittle fracture model) is compared with linear elastic fracture mechanics (LEFM) analytical solutions. Our model is further utilized to simulate fission-gas-induced fragmentation of the UO 2 high-burnup structure (HBS). Simulation results show that pressurized bubbles can cause crack nucleation and propagation, and that the bubble size and the surrounding external pressure affect the critical pressure corresponding to crack nucleation. Simulations of a partial HBS at different recrystallization stages show that different grain structures (due to recrystallization) also influence crack paths and fragmentation morphology.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A robust interface finite element formulation for modeling brittle material failure problems

Abstract Failure of many brittle materials and structures can be modeled using interface‐oriented finite elements combined with intrinsic cohesive zone models. The discontinuous Galerkin (DG) finite element method provides an innovative framework for modeling brittle crack propagation with zero‐thickness interface elements, which can accommodate extrinsic cohesive laws to avoid the artificial compliance required in intrinsic cohesive models. However, robust formulations and implementations of DG methods are critical in alleviating the well‐known convergence issues for both crack nucleation and propagation with reduced instability. This paper presents a robust interface element formulation by modifying the incomplete interior penalty Galerkin (IIPG) method, which successfully avoids the initial element interface penetration across elements that occurs prior to crack nucleation, and thereby greatly reduces the instability issue as cracks open. We further verified and validated our implementation by using a bar tension test and a beam fracturing benchmark. The robustness of our proposed interface element method was demonstrated by a micromechanics fiber/matrix debonding problem with 64 fibers embedded in a bulk matrix.

Liu, Ruijie↗
DOE OSTI2019

Direct measurements of slip irreversibility in a nickel-based superalloy using high resolution digital image correlation

Fatigue crack nucleation in crystalline materials typically develops due to highly localized cyclic slip. During a fatigue cycle, reverse slip differs locally from slip in the forward direction particularly in precipitate-containing materials such as superalloys. In this paper we report the first direct measurements of irreversibility at the scale of individual slip bands by high-resolution digital image correlation (DIC) in a polycrystalline nickel-based superalloy. Quantitative measurements of the slip irreversibility are challenging for regions of material that have a size that captures the microstructure and its variability. High spatial resolution at the nanometer scale during experimental measurements is needed to observe slip localization during deformation. Moreover, large fields are also needed to obtain the material response over statistically representative populations of microstructural configurations. Recently, high resolution scanning electron microscope (SEM) digital image correlation (DIC) has been extended for quantitative analysis of discontinuities induced by slip events using the Heaviside-DIC method. Therefore, this novel method provides quantitative measurements of slip localization at the specimen surface. In this paper, the Heaviside-DIC method is used to measure slip irreversibility and plastic strain accumulation in a nickel-based superalloy. The method detects bands with high levels of irreversibility early in cycling that ultimately form fatigue cracks upon further cycling. The local microstructural configurations that induce large amounts of plasticity and slip irreversibility are correlated to crack nucleation locations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Four-dimensional microstructurally small fatigue crack growth in cyclically loaded nickel superalloy specimen

Here, we nondestructively image a low-solvus high refractory (LSHR) nickel base specimen microstructure during cyclic loading to characterize crack propagation during the microstructurally small crack (MSC) growth regime. Micro-computed tomography (uCT) is used to characterize the material density and identify crack interfaces, while high-energy diffraction microscopy (HEDM) provides snapshots of the microstructural orientation field changes. From these measurements, data are computationally reconstructed and coregistered to describe the specimen’s configuration before load, before crack nucleation, and throughout the crack propagation/loading process. An MSC forms at a FIB-cut notch, which was previously placed to concentrate stress at the center of the specimen gauge, and propagates along a facet parallel to the local {1 1 1} crystallographic plane. Analysis of available crystallographic slip directions and measured point-to-point disorientation axes indicate probable single and coplanar slip activity along FCC octahedral planes within this region. Similar analyses of other crack propagation facets reveal generally slower growth along directions favoring zig-zag cross and codirectional cross slip, which are often disfavored at room temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Influence of Nominal Composition on the Microstructure, Tensile Properties, and Weldability of Cast Monel Alloys

Cast Monel alloys are used in many industrial applications that require a combination of good mechanical properties and excellent resistance to corrosion. Despite relative widespread use, there has been limited prior research investigating the fundamental composition–structure–property relationships. Here in this work, microstructural characterization, thermal analysis, electron probe microanalysis, tensile testing, and Varestraint testing were used to assess the effects of variations in nominal composition on the solidification path, microstructure, mechanical properties, and solidification cracking susceptibility of cast Monel alloys. It was found that Si segregation caused the formation of silicides at the end of solidification in grades containing at least 3 wt pct Si. While increases to Si content led to significant improvements in strengthening due to the precipitation of β 1 -Ni 3 Si, the silicide eutectics acted as crack nucleation sites during tensile loading which severely reduced ductility. The solidification cracking susceptibility of low-Si Monel alloys was found to be relatively low. However, increases to Si concentration and the onset of associated eutectic reactions increased the solidification temperature range and drastically reduced cracking resistance. Increases in the Cu and Mn concentrations were found to reduce the solubility limit of Si in austenite which promoted additional eutectic formation and exacerbated the reductions in ductility and/or weldability.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Local structural ordering determines the mechanical damage tolerance of amorphous grain boundary complexions

Amorphous grain boundary complexions act as toughening features within a microstructure because they can absorb dislocations more efficiently than traditional grain boundaries. This toughening effect should be a strong function of the local internal structure of the complexion, which has recently been shown to be determined by grain boundary crystallography. To test this hypothesis, molecular dynamics are used here to simulate dislocation absorption and damage nucleation for complexions with different distributions of structural short-range order. Here, the complexion with a more disordered structure away from the dislocation absorption site is actually found to better resist crack nucleation, as damage tolerance requires delocalized deformation and the operation of shear-transformation zones through the complexion thickness. The more damage tolerant complexion accommodates plastic strain efficiently within the entire complexion, providing the key mechanistic insight that local patterning and asymmetry of structural short-range order controls the toughening effect of amorphous complexions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

On the convexity of phase-field fracture formulations: Analytical study and comparison of various degradation functions

Efficient and accurate fracture modeling is of great importance in applications where catastrophic outcomes under extreme scenarios are possible. The phase-field (PF) approach to fracture received significant attention over the past decade, due to its capability to capture complicated fracture patterns (e.g., crack merging and branching). Specifically, crack initiation and propagation are modeled via minimization of the total energy functional, which is regularized with the aid of a phase field. Despite the promising results and modeling capabilities of the PF method in many applications, the solution of fracture problems remains computationally challenging mainly due to the non-convexity of the total energy functional with respect to the combined unknown (phase field and displacement) fields. Understanding the effects of their coupling on convexity is crucial in order to address frequently encountered hurdles in fracture modeling (e.g., inefficient solvers and non-physical crack nucleation). In this paper, we develop convexity criteria for a wide class of PF fracture formulations. For this class of formulations, the second variation of the total energy functional is expressed in terms of Hessian matrices (evaluated at individual material points). Depending on the choice of geometric crack functions and degradation functions, we classify the formulations into three categories and analytically study each one separately. To study the sign of the second variation, we derive inequalities which are satisfied at material points when the Hessian matrix is locally positive semi-definite. These inequalities provide objective criteria for comparing degradation functions. Finally, the applicability of the proposed convexity criteria is demonstrated in the context of a one-dimensional problem, solved using a conventional monolithic solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Resolving electrochemically triggered topological defect dynamics and structural degradation in layered oxides

Understanding topological defects-controlled structural degradation of layered oxides—a key cathode material for high-performance lithium-ion batteries—plays a critical role in developing next-generation cathode materials. Here, by constructing a nanobattery in an electron microscope enabling atomic-scale monitoring of electrochemcial reactions, we captured the electrochemically driven atomistic dynamics and evolution of dislocations—a most important topological defect in material. We deciphered how dislocations nucleate, move, and annihilate within layered cathodes at the atomic scale. Specifically, we found two types of dislocation configurations, i.e., single dislocations and dislocation dipoles. Both pure dislocation glide/climb and mixed motions were captured, and the dislocation glide and climb velocities were first experimentally measured. Moreover, dislocation activity-mediated structural degradation such as crack nucleation, phase transformation, and lattice reorientation was unraveled. Our work provides deep insights into the atomistic dynamics of electrochemically driven dislocation activities in layered oxides.

Wang, Chunyang (ORCID:0000000184613952)↗
DOE OSTI2020

Role of Material Condition in Precursor Corrosion Damage and Stress Corrosion Crack Initiation Behavior of Alloy 600 in PWR Primary Water

Stress corrosion crack (SCC) initiation testing has been performed on a 15% cold-worked UNS N06600 (Alloy 600) heat in the mill-annealed (MA), solution annealed (SA), and thermally treated (TT) conditions to assess the role of grain boundary (GB) carbides on stress-assisted intergranular attack (IGA) and short crack nucleation and growth. The SCC initiation tests were conducted at constant load equivalent to materials yield stress in 360oC simulated PWR primary water with dissolved hydrogen at the Ni/NiO stability line. Results revealed the highest SCC initiation susceptibility in the Alloy 600 MA material followed by the TT and SA materials, suggesting the GB carbide distribution did not have a controlling effect on SCC initiation resistance. Quantitative assessments were conducted on IGA and short cracks to help understand this phenomenon and the role of GB carbides on precursor damage development that lead to differences in macroscopic SCC initiation behavior are discussed.

Zhai, Ziqing↗
DOE OSTI2025

Towards robust surrogate models: Benchmarking machine learning approaches to expediting phase field simulations of brittle fracture

Data-driven approaches have the potential to make modeling complex, nonlinear physical phenomena significantly more computationally tractable. For example, computational modeling of fracture is a core challenge where machine learning techniques have the potential to provide a much needed speedup that would enable progress in areas such as multi-scale modeling and uncertainty quantification. Currently, phase field modeling (PFM) of fracture is one such approach that offers a convenient variational formulation to model crack nucleation, branching and propagation. To date, machine learning techniques have shown promise in approximating PFM simulations. While standard fracture benchmarks represent realistic scenarios frequently observed in practice, they typically do not provide sufficiently challenging tests for data-driven methods. Here, to address this gap, we introduce a challenging dataset based on PFM simulations designed to benchmark and advance ML methods for fracture modeling. This dataset includes three energy decomposition methods, two boundary conditions, and 1000 random initial crack configurations for a total of 6000 simulations. Each sample contains 100 time steps capturing the temporal evolution of the crack field. Alongside this dataset, we also implement and evaluate Physics Informed Neural Networks (PINN), Fourier Neural Operators (FNO), and UNet models as baselines, and explore the impact of ensembling strategies on prediction accuracy. With this combination of our dataset and baseline models drawn from the literature we aim to provide a standardized and challenging benchmark for evaluating machine learning approaches to solid mechanics. Our results highlight both the promise and limitations of popular current models, and demonstrate the utility of this dataset as a testbed for advancing machine learning in fracture mechanics research.

Benchmark dataset↗
DOE OSTI2021

Thermomechanical Microstructural Predictions of Fracture Nucleation of Zircaloy-4 Alloys With δ and ε Hydride Distributions

Abstract A crystalline dislocation-density formulation that was incorporated with a nonlinear finite-element (FE) method was utilized to understand and to predict the thermomechanical behavior of an hexagonal closest packed (h.c.p.) zircaloy system with hydrides with either face-centered cubic (f.c.c.) or body-centered cubic (b.c.c.) hydrides. This formulation was then used with a recently developed fracture methodology that is adapted for finite inelastic strains and multiphase crystalline systems to understand how different microstructurally based fracture modes nucleate and propagate. The interrelated microstructural characteristics of the different crystalline hydride and matrix phases with the necessary orientation relationships (ORs) have been represented, such that a detailed physical understanding of fracture nucleation and propagation can be predicted for the simultaneous thermomechanical failure modes of hydride populations and the matrix. The effects of volume fraction, morphology, crystalline structure, and orientation and distribution of the hydrides on simultaneous and multiple fracture modes were investigated for radial, circumferential, and mixed distributions. Another key aspect was accounting for temperatures changes due to the effects of thermal conduction and dissipated plastic work and their collective effects on fracture. For hydrided aggregates subjected to high temperatures, thermal softening resulted in higher ductility due to increased dislocation-density activity, which led to higher shear strain accumulation and inhibited crack nucleation and growth. The predictions provide validated insights into why circumferential hydrides are more fracture-resistant than radial hydrides for different volume fractions and thermomechanical loading conditions.

Engineering↗
DOE OSTI2024

Poromechanical cohesive interface element with combined Mode I-II cohesive zone elastoplasticity for simulating fracture in fluid-saturated porous media

A combined Mode I-II cohesive zone (CZ) elasto-plastic constitutive model, and a two-dimensional (2D) cohesive interface element (CIE) are formulated and implemented at small strain within an ABAQUS User Element (UEL) for simulating 2D crack nucleation and propagation in fluid-saturated porous media. Here, the CZ model mitigates problems of convergence for the global Newton-Raphson solver within ABAQUS, which when combined with a viscous stabilization procedure allows for simulation of post-peak response under load control for coupled poromechanical finite element analysis, such as concrete gravity dam stability analysis. Verification examples are presented, along with a more complex ambient limestone-concrete wedge fracture experiment, water-pressurized concrete wedge experiment, and concrete gravity dam stability analyses. A calibration procedure for estimating the CZ parameters is demonstrated with the limestone-concrete wedge fracture process. For the water-pressurized concrete wedge fracture experiment it is shown that the inherent time-dependence of the poromechanical CIE analysis provides a good match with experimental force versus displacement results at various crack mouth opening rates, yet misses the pore water pressure evolution ahead of the crack tip propagation. This is likely a result of the concrete being partially-saturated in the experiment, whereas the finite element analysis assumes fully water saturated concrete. For the concrete gravity dam analysis, it is shown that base crack opening and associated water uplift pressure leads to a reduced Factor of Safety, which is confirmed by separate analytical calculations.

97 MATHEMATICS AND COMPUTING↗