Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CRACKING”

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 OSTI2016

A method to generate conformal finite-element meshes from 3D measurements of microstructurally small fatigue-crack propagation [A method to generate conformal finite-element meshes from 3D measurements of microstructurally small fatigue-crack propagation: 3D Meshes of Microstructurally Small Crack Growth]

In an effort to reproduce computationally the observed evolution of microstructurally small fatigue cracks (MSFCs), a method is presented for generating conformal, finite-element (FE), volume meshes from 3D measurements of MSFC propagation. The resulting volume meshes contain traction-free surfaces that conform to incrementally measured 3D crack shapes. Grain morphologies measured using near-field high-energy X-ray diffraction microscopy are also represented within the FE volume meshes. Proof-of-concept simulations are performed to demonstrate the utility of the mesh-generation method. The proof-of-concept simulations employ a crystal-plasticity constitutive model and are performed using the conformal FE meshes corresponding to successive crack-growth increments. Although the simulations for each crack increment are currently independent of one another, they need not be, and transfer of material-state information among successive crack-increment meshes is discussed. The mesh-generation method was developed using post-mortem measurements, yet it is general enough that it can be applied to in-situ measurements of 3D MSFC propagation.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Assessing Gaseous Hydrogen Assisted Fatigue Crack Growth Susceptibility of Pipeline Steel Weld Fusion Zones and Heat Affected Zones

The objective of this work was twofold: (1) measure reliable fatigue crack growth relationships for X65 steel and its girth weld in high-pressure hydrogen gas to enable structural integrity assessments of hydrogen pipelines, and (2) evaluate the hydrogen accelerated fatigue crack growth susceptibility of the weld fusion zone and heat-affected zone relative to the base metal. Fatigue crack growth relationships (da/dN versus ΔK) were measured for girth welded X65 pipeline steel in high pressure hydrogen gas (21 MPa) and in air. Hydrogen assisted fatigue crack growth was observed for the base metal (BM), fusion zone (FZ), and heat-affected zone (HAZ), and was manifested through crack growth rates reaching nearly an order of magnitude acceleration over rates in air. Addditionally, at higher ΔK values, crack growth rates of BM, FZ, and HAZ were coincident; however, at lower ΔK, the fatigue crack growth relationships exhibited some divergence with the fusion zone having the highest crack growth rates. These relative fatigue crack growth rates in the BM, FZ, and HAZ were provisional, however, since both crack closure and residual stress contributed to the crack-tip driving force in specimens extracted from the HAZ. Despite the relatively high applied R-ratio (R = 0.5), crack closure was detected in the heat affected zone tests, in contrast to the absence of crack closure in the base metal tests. Crack closure corrections were performed using the adjusted compliance ratio method and the effect of residual stress on K max was determined by the crack-compliance method. Crack-tip driving forces that account for closure and residual stress effects were quantified as a weighted function of ΔK and K max (i.e., K norm ), and the resulting da/dN versus K norm relationships showed that the HAZ exhibited higher hydrogen accelerated fatigue crack growth rates than the BM at lower K norm values.

36 MATERIALS SCIENCE↗
DOE OSTI2025

In Situ Measurement of Stress Corrosion Crack Tip Properties During Environmentally Induced Cracking in Stainless Steel 304L

The influence of the external cathodic area, stress intensity, and electrochemical polarization level on the resultant crack growth rate and crack tip pH and potential was studied for stainless steel 304L in a sodium chloride environment. Before testing, the validity of in situ crack tip property measurements using microprobes and drilled holes was supported through modeling efforts and confirmed with experimental observations. Here, it was found that increasing the R-ratio during corrosion fatigue increased the crack growth rate and decreased the measured crack tip pH of the 304L sample. Decreasing the electrochemical potential (i.e., cathodic polarization) of the bulk sample increased the crack tip pH and stopped measured crack extension under constant K conditions. Conversely, increasing the potential of the bulk sample decreased the crack tip pH and induced cracking under constant K conditions. Decreasing the cathodic area external to the crack increased the measured pH at the crack tip and decreased the crack growth rate under fatigue conditions. For the given configuration, environment, and alloy studied, fatigue crack growth is not sustained by cathodic reactions solely occurring in the crack wake and notch. These results have implications for the ability to predict atmospheric crack growth rates from laboratory, full immersion experiments.

corrosion fatigue↗
DOE OSTI2019

The influence of plasticity-induced crack closure on creep-fatigue crack growth in two heat-resistant steels

Finite element simulations are performed to predict plasticity-induced crack closure and creep-fatigue crack growth rates in two heat-resistant steels, the austenitic 20Cr-25Ni (Alloy 709) and martensitic 9Cr-1Mo (modified 1Cr-9Mo) steels. In a creep-fatigue cycle, total crack growth rate is usually computed by the addition of the fatigue crack growth rate during cyclic loading and creep crack growth rate during hold time. Two-dimensional finite element analyses of compact tension specimens are performed to simulate crack growth under cyclic and time-dependent loading conditions, with the consideration of elastic-plastic and creep deformations of the material at the crack tip. The simulations quantify the effect of hold time on crack opening load induced by the combined action of plasticity-induced crack closure and creep-induced stress relaxation at the crack tip. It is shown that increasing hold time during a creep-fatigue cycle results in a decrease of crack-tip opening load, thus increasing crack growth rate during the next cyclic loading. Experimental results of creep-fatigue crack growth rates in Alloy 709 are presented. The finite element simulations predicted crack growth rates under creep-fatigue loading in agreement with experimental values for both Alloy 709 and modified 9Cr-1Mo steels.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Direct observations and characterization of crack closure during microstructurally small fatigue crack growth via in-situ high-energy X-ray characterization

This study provides direct observation of the crack closure mechanism of a naturally occurring, tortuous, 3D microstructurally small fatigue crack (SFC) in additively manufactured Inconel 718. In-situ non-destructive characterization is performed using high-energy X-ray diffraction techniques to capture the evolution of the 3D microstructure and micromechanical response in the vicinity of the crack front. Based on the stress state of twelve grains analyzed at the crack tip, the crack closure events of the SFC front was found to be spatially heterogeneous with respect to loading progression governed by the local stress state of the grains (specifically the stress reversal from compression to tension). From this analysis, three grains that displayed different degrees of crack closure were further investigated, based on the orientation of the crack relative to the grains and the associated modality of crack growth. The stress normal to the crack plane and the associated degree of Mode I crack behavior were correlated with events of the crack opening earlier during the loading cycle. This was further corroborated by diffraction spot spreading analysis that quantified the crystallographic lattice distortion caused by the opening crack. Additionally, the detailed characterization of the opening behavior of the grains located at the crack tip and their associated states of stress elucidates the mechanism governing crack closure and will inform future modeling efforts of this phenomenon.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Effect of micromorphology at the fatigue crack tip on the crack growth in electron beam welded Ti-6Al-4V joint

In this paper, we describe experiments on welded joints of Ti-6Al-4V alloy specimens exhibiting fatigue characteristics in the base metal (BM), hot affected zone (HAZ) and fuse zone (FZ). The effect of micromorphology on crack propagation at the tip of the fatigue crack in joints formed by electron beam welding was investigated using an optical microscope, transmission electron microscope and other methodologies. The results demonstrated that the fatigue crack originated in and propagated along α/β boundaries in the BM. In the HAZ, the fatigue crack occurred at the boundary between martensite laths, and propagated through most irregular-equiaxed α phases and a few martensite laths. In the FZ, the fatigue crack originated at the boundaries between the fine crushing phases among martensite laths, and propagated along a majority of α/β boundaries and several narrow martensite laths. The electron beam welded joint of Ti-6Al-4V alloy showed instances of zigzag fatigue cracks that increased in degree from lowest in the HAZ, moderate in the FZ to greatest in the BM. Conversely, fatigue crack growth rate (FCGR) was greatest in the HAZ, less in the FZ and slowest in the BM. - Highlights: •Ti-6Al-4V welded joint exhibits different fatigue characteristics. •The fatigue crack propagates along α/β boundaries in the BM. •The fatigue crack propagates through α phases and martensite laths in the HAZ. •The fatigue crack propagates along α/β boundaries and martensite laths in the FZ. •Fatigue crack growth rate is fastest in the HAZ, less in the FZ, slowest in the BM.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A mesh-independent method for planar three-dimensional crack growth in finite domains: A mesh-independent method for crack growth in 3D finite domains

The discrete crack mechanics (DCM) method is a dislocation-based crack modeling technique where cracks are constructed using Volterra dislocation loops. The method allows for the natural introduction of displacement discontinuities, avoiding numerically expensive techniques. Mesh dependence in existing computational modeling of crack growth is eliminated by utilizing a superposition procedure. The elastic field of cracks in finite bodies is separated into two parts: the infinite-medium solution of discrete dislocations and an finite element method solution of a correction problem that satisfies external boundary conditions. In the DCM, a crack is represented by a dislocation array with a fixed outer loop determining the crack tip position encompassing additional concentric loops free to expand or contract. Solving for the equilibrium positions of the inner loops gives the crack shape and stress field. The equation of motion governing the crack tip is developed for quasi-static growth problems. Additionally, convergence and accuracy of the DCM method are verified with two- and three-dimensional problems with well-known solutions. Crack growth is simulated under load and displacement (rotation) control. In the latter case, a semicircular surface crack in a bent prismatic beam is shown to change shape as it propagates inward, stopping as the imposed rotation is accommodated.

42 ENGINEERING↗
DOE OSTI2023

Inversion-based correction of Double-Torsion (DT) subcritical crack growth tests for crack profile geometry

Because of its simplicity and the ability to produce a stable, slow-propagating crack, the Double-Torsion (DT) method has been used widely for investigating the critical and subcritical propagation of a slow-propagating tensile (mode-I) crack. However, to determine the complex relationship between the crack velocity $\mathcal{v_c}$ vs. the strain energy release rate $\mathscr{G}$ (or the stress intensity factor K) from laboratory measurements, several corrections must be made to account for the impact of sample and crack geometry. Particularly, DT test typically produces a crack with a curved edge profile instead of a straight line, causing the local $\mathcal{v_c}$ and $\mathscr{G}$ vary along the crack front. The experimentally measured $\mathcal{v_c}$ and $\mathscr{G}$ data merely reflect collective, averaged behavior of the crack. This makes inversion for the intrinsic, “true” crack growth kinetics necessary, based upon the knowledge of the crack geometry. Simple and effective correction methods have been proposed and validated for the slow, chemical-reaction-controlled part (Region I) of the $\mathcal{v_c}$-$\mathscr{G}$ curve. However, reliable methods for the highly nonlinear, transport-dominated part (Region II) and its sudden transition to the dynamic propagation part (Region III) are still lacking. Here we propose a method for determining the intrinsic $\mathcal{v_c}$-$\mathscr{G}$ relationship across all three Regions based upon DT test data, using a simple model function and its numerical inversion. The performance of this approach is examined and demonstrated using both synthetic and laboratory data for subcritical crack growth in soda lime glass.

42 ENGINEERING↗
DOE OSTI2021

Hydro-chemo-mechanical phase field formulation for corrosion induced cracking in reinforced concrete

Steel corrosion induced cracking in reinforced concrete structures is often caused by volumetric straining at the steel-concrete interface. This volumetric strain is mainly controlled by the rate of generation of corrosion product, which depends on the diffusion of Cl{sup −} and dissolved O{sub 2} as well as the chemical reaction at the steel-concrete interface. While many numerical modeling techniques have been proposed to tackle this complex phenomenon, most have been limited to using staggered solution schemes, where the diffusion, chemical reaction, mechanical strain, and crack fields are solved for separately. Unfortunately, such approaches introduce numerical errors that can lead to unrealistic predictions. Therefore, a coupled solution, is needed to remove the formulation induced errors and understand this complex phenomenon. To provide a unifying solution framework for corrosion induced cracking in reinforced concrete, a phase field formulation is presented and implemented with FEM algorithms. The governing equations are derived from the system's free energy and Fick's law of diffusion. The proposed formulation was first validated in 1D by comparing the modeling results with the semi-analytical solution approximated via Laurent series to avoid singularities in the complex domain. 2D plane strain solutions to the embedded rebar undergoing corrosion induced cracking with pre-cracks are also presented and compared with the approximated solutions from linear elastic fracture mechanics theory (LEFM). The modeling of crack propagation enabled by the proposed formulation is also presented along with parametric studies to reveal the roles of the fracture toughness and permeability of the concrete on the corrosion induced cracking. A case study was also conducted on H-pile steel with concrete jackets to demonstrate the feasibility of modeling corrosion induced fracture in reinforced concrete structures with complex steel-concrete interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2018

An arc-length method for controlled cohesive crack propagation using high-order XFEM and Irwin’s crack closure integral

Numerical modeling of cohesive crack growth in quasi-brittle materials is challenging, primarily due to the combination of (i) nonlinearity associated with the fracture process zone (FPZ), (ii) arbitrary directions to which a crack may propagate, and (iii) snap-back or snap-through instabilities encountered in the response of the structure. To address these challenges, here we propose a novel arc-length method that can follow the equilibrium path of cohesive crack propagation. The proposed approach is based on the extended finite element method (XFEM) with scalar high-order enrichment functions and Irwin’s crack closure integral, which allows for direct control of the applied loads necessary to propagate cohesive cracks. This is achieved by augmenting a constraint equation written in terms of stress intensity factors (SIFs), and expressed explicitly in terms of the enriched degrees of freedom, which is an attractive feature achieved with Irwin’s integral, since SIFs can be written in closed-form. Note that singular enrichments are active in an unstable crack propagation state and automatically vanish in stable crack configurations. Furthermore, to propagate cracks in arbitrary directions, we employ a maximum circumferential stress criterion implemented by (i) direct usage of the SIFs, and by (ii) a new stress-based nonlocal implementation of this principle. Various benchmark problems including pure mode I and mixed-mode fracture are solved to demonstrate the predictive capability of the present framework for cohesive crack modeling.

42 ENGINEERING↗
DOE OSTI2018

Predicting the 3D fatigue crack growth rate of small cracks using multimodal data via Bayesian networks: In-situ experiments and crystal plasticity simulations

Small crack propagation accounts for most of the fatigue life of engineering structures subject to high cycle fatigue loading conditions. Determining the fatigue crack growth rate of small cracks propagating into polycrystalline engineering alloys is critical to improving fatigue life predictions, thus lowering cost and increasing safety. In this work, cycle-by-cycle data of a small crack propagating in a beta metastable titanium alloy is available via phase and diffraction contrast tomography. Crystal plasticity simulations are used to supplement experimental data regarding the micromechanical fields ahead of the crack tip. Experimental and numerical results are combined into a multimodal dataset and sampled utilizing a non-local data mining procedure. Furthermore, to capture the propensity of body-centered cubic metals to deform according to the pencil-glide model, a non-local driving force is postulated. The proposed driving force serves as the basis to construct a data-driven probabilistic crack propagation framework using Bayesian networks as building blocks. The spatial correlation between the postulated driving force and experimental observations is obtained by analyzing the results of the proposed framework. Results show that the above correlation increases proportionally to the distance from the crack front until the edge of the plastic zone. Moreover, the predictions of the propagation framework show good agreement with experimental observations. Finally, we studied the interaction of a small crack with grain boundaries (GBs) utilizing various slip transmission criteria, revealing the tendency of a crack to cross a GB by propagating along the slip directions minimizing the residual Burgers vector within the GB.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Determination of Constraint-Independent Crack Tip Opening Angle for Stable Crack Growth in High-Strength Ductile Steels

The crack tip opening angle (CTOA) is one of fracture toughness parameters that has been used for decades in describing large stable crack growth in thin-walled aerospace structures under the low-constraint conditions. Recently, the pipeline industry has developed a growing interest in using the CTOA parameter to serve as the minimum required fracture toughness to arrest dynamic crack propagation in modern gas transmission pipelines made of high-strength ductile steel. To meet this industrial need, the CTOA test standard ASTM E3039 was therefore developed for measuring the constant critical CTOA. ASTM E3039 recommends a drop weight tearing test (DWTT) specimen with a shallow crack for standard CTOA testing, but its CTOA may depend on the low constraint condition of the DWTT specimen at the crack tip. Verifying the constraint independence of the DWTT-measured CTOA thus becomes indispensable for applying CTOA toughness to the running fracture control in the pipeline design. For this purpose, the present paper evaluates critical CTOA values in a set of fracture toughness tests on single-edge notched bend (SENB) specimens with shallow and deep cracks, based on four CTOA estimation models. Among these, the Ln(P)-LLD linear fit model is similar to that recommended by ASTM E3039 for CTOA calculation. Fracture test data for X80 pipeline steel and HY80 structural steel were considered in the CTOA evaluation. The results showed that the four CTOA models were able to determine a constraint independent CTOA value for stable crack growth in the SENB specimens. As a result, a single, reliable, constant CTOA value could be determined regardless of the specimen geometry or the crack-tip constraint conditions. Therefore, the CTOA measured using ASTM E3039 is constraint-independent and transferable to use in cases of actual cracks propagating in gas transmission pipelines.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Instability in dynamic fracture and the failure of the classical theory of cracks

Cracks, the major vehicle for material failure, undergo a micro-branching instability at ~40% of their sonic limiting velocity in three dimensions. Recent experiments showed that in thin systems cracks accelerate to nearly their limiting velocity without micro-branching, until undergoing an oscillatory instability. Despite their fundamental importance, these dynamic instabilities are not explained by the classical theory of cracks, which is based on linear elasticity and an extraneous local symmetry criterion to predict crack paths. We develop a two-dimensional theory for predicting arbitrary paths of ultrahigh-speed cracks, which incorporates elastic nonlinearity without extraneous criteria. We show that cracks undergo an oscillatory instability controlled by small-scale, near crack-tip, elastic nonlinearity. This instability occurs above an ultrahigh critical velocity and features an intrinsic wavelength proportional to the ratio of the fracture energy to the elastic modulus, in quantitative agreement with experiments. As a result, this ratio emerges as a fundamental scaling length assumed to play no role in the classical theory of cracks, but shown here to strongly influence crack dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

An Enhanced Virtual Crack Closure Technique for Stress Intensity Factor Calculation along Arbitrary Crack Fronts and the Application in Hydraulic Fracturing Simulation

The virtual crack closure technique (VCCT) is widely used for calculating energy release rates along crack fronts and modeling the propagation of cracks in solid materials. Although the VCCT formulation for smooth crack fronts has been sufficiently addressed in the literature, the application of VCCT to a nonsmoothed crack front with sharp corners warrants further investigation. In this study, we describe an enhanced VCCT to calculate energy release rates and stress intensity factors for cracks with arbitrary shapes in 3D domains discretized on structured grids. The formulations of the enhanced VCCT were developed and implemented into a multiphysics simulation environment capable of simulating crack propagation in the framework of linear elastic fracture mechanics. Comparisons with existing analytical/numerical solutions and other VCCT approaches were performed to verify the enhanced VCCT in terms of SIF calculation along nonsmoothed crack fronts. We then applied the enhanced VCCT to a hydraulically driven penny-shaped fracture problem to further demonstrate its capability to simulate nonsmoothed fracture propagation.

02 PETROLEUM↗