Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture propagation”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2026

Subsurface mechanical damage of fused silica glass during grinding by various sub-aperture tools with and without ultrasonics

The subsurface mechanical damage (SSD) depth after grinding fused silica glass with a comprehensive set of sub-aperture fixed abrasive grinding tools [cup, wheel, belt, pad, and rotary face mill (with and without ultrasonics)] and process parameters has been statistically measured using the taper wedge technique and evaluated. Consistent with a previously reported grinding model [J. Non-Cryst. Solids 352, 5601–5617 (2006) Crossref , Materials Science and Technology of Optical Fabrication (Wiley & Sons, 2018)], based on the sliding indentation fracture by sliding particles or asperities where the normal load per particle determines the depth of the fracture (and ultimately the overall SSD depth), the dominant factor controlling SSD depth was found to be the abrasive size regardless of the tool type and process conditions. Compared to full aperture grinding methods, the overall SSD depth was higher using the sub-aperture tools, likely due to the higher effective pressure and higher load per particle distribution. Here, in addition to abrasive size, a significant reduction in SSD depth was achieved by: (1) reducing the load distribution on the abrasive particles via increase in contact area and/or decrease in mechanical loading; (2) using a more compliant host tool medium; and (3) in what we believe is a more novel way, using ultrasonics. Combining low abrasive size, larger contact area, and a compliant host, the 6 µm diamond in a resin matrix (Trizact) on a foam pad led to very low SSD depth (~ 4.6 µm), relatively fast grinding rate (186 mm 3 /h), and little or no figure degradation. This grinding tool/process is an attractive choice for final grind, resulting in significantly reduced polish out (i.e., “grey out”) times. With the rotary face mill tool, the use of ultrasonics consistently led to a SSD depth reduction (ranging from 17%–34%). A new fracture mechanics-based model, to the best of our knowledge, where the relevant normal load is parallel to the feed direction, has been developed to explain how ultrasonics leads to lower SSD depth. The key factors, supported by finite element stress analysis and load measurements, are (1) the initiation of fractures at higher z heights during the tool’s ultrasonic vertical oscillations, thus propagating less deep into workpiece; (2) reduction in load (and therefore reduction in fracture propagation distance) due to smaller tool-workpiece feed direction contact area (again caused by higher heights relative to depth of cut); (3) upward movement of the tool during oscillation leads to fracturing toward the surface instead into the depth; and finally (4) at tool’s lowest point of oscillation cycle, there may not be enough time for the fracture to propagate to its full length.

Optics and optical instruments↗
NASA NTRS1992

Static and dynamic strain energy release rates in toughened thermosetting composite laminates

In this work, the static and dynamic fracture properties of several thermosetting resin based composite laminates are presented. Two classes of materials are explored. These are homogeneous, thermosetting resins and toughened, multi-phase, thermosetting resin systems. Multi-phase resin materials have shown enhancement over homogenous materials with respect to damage resistance. The development of new dynamic tests are presented for composite laminates based on Width Tapered Double Cantilevered Beam (WTDCB) for Mode 1 fracture and the End Notched Flexure (ENF) specimen. The WTDCB sample was loaded via a low inertia, pneumatic cylinder to produce rapid cross-head displacements. A high rate, piezo-electric load cell and an accelerometer were mounted on the specimen. A digital oscilloscope was used for data acquisition. Typical static and dynamic load versus displacement plots are presented. The ENF specimen was impacted in three point bending with an instrumented impact tower. Fracture initiation and propagation energies under static and dynamic conditions were determined analytically and experimentally. The test results for Mode 1 fracture are relatively insensitive to strain rate effects for the laminates tested in this study. The test results from Mode 2 fracture indicate that the toughened systems provide superior fracture initiation and higher resistance to propagation under dynamic conditions. While the static fracture properties of the homogeneous systems may be relatively high, the apparent Mode 2 dynamic critical strain energy release rate drops significantly. The results indicate that static Mode 2 fracture testing is inadequate for determining the fracture performance of composite structures subjected to conditions such as low velocity impact. A good correlation between the basic Mode 2 dynamic fracture properties and the performance is a combined material/structural Compression After Impact (CAI) test is found. These results underscore the importance of examining rate-dependent behavior for determining the longevity of structures manufactured from composite materials.

Cairns, Douglas S.↗
NASA NTRS1991

Dynamic delamination crack propagation in a graphite/epoxy laminate

Dynamic delamination crack propagation in a (90/0) 5s Graphite/Epoxy laminate with an embedded interfacial crack was investigated experimentally using high speed photography. The dynamic motion was produced by impacting the beamlike laminate specimen with a silicon rubber ball. The threshold impact velocities required to initiate dynamic crack propagation in laminates with varying initial crack positions were determined. The crack propagation speeds were estimated from the photographs. Results show that the through the thickness position of the embedded crack can significantly affect the dominant mechanism and the threshold impact velocity for the onset of crack movement. If the initial delamination is placed near the top of bottom surface of the laminate, local buckling of the delaminated plies may cause instability of the crack. If the initial delamination lies on the midplane, local buckling does not occur and the initiation of crack propagation appears to be dominated by Mode II fracture. The crack propagation and arrest observed was seen to be affected by wave motion within the delamination region.

Grady, J. E.↗
DOE OSTI2020

A numerical-homogenization based phase-field fracture modeling of linear elastic heterogeneous porous media

Most porous media, such as geomaterials and biomaterials are highly heterogeneous in nature, and they contain large variations of microscopic pore structures, such as pore sizes, pore distribution, and pore shapes. The oscillation of microscopic structures is a substantial challenge in theoretical characterization and is usually ignored in continuous modeling. However, mechanical behavior of porous media such as deformation and failure, are essentially impacted by the microscopic heterogeneity which needs to be considered in modeling a porous media. Here, this research proposes a numerical modeling framework with a capability to investigate the effect of microscopic heterogeneity on the macroscopic fracture behavior in porous media by using a numerical homogenization technique, combined with the phase-field fracture modeling method. This numerical modeling strategy computes a homogenized elasticity tensor based on microscopic heterogeneous pore structures heterogenous porous domain by solving boundary value problems at microscopic domain. The strain energy and subsequent propagation of macroscopic fractures will be updated using homogenized stiffness information. Using this numerical scheme, the microscopic pore structure’s impact on the fracture behavior through the homogenized elastic tensor will be taken into account. This multiscale technique is benchmarked against classical problems. Finally, the results highlight the importance of the underlying pore structure and reveal that both fracture strength and propagation path can be influenced by the microscopic heterogeneity.

36 MATERIALS SCIENCE↗
NASA NTRS1985

Dynamic fracture mechanics

Dynamic fracture and crack propagation concepts for ductile materials are reviewed. The equations for calculating dynamic stress integrity and the dynamic energy release rate in order to study dynamic crack propagation are provided. The stress intensity factor versus crack velocity relation is investigated. The uses of optical experimental techniques and finite element methods for fracture analyses are described. The fracture criteria for a rapidly propagating crack under mixed mode conditions are discussed; crack extension and fracture criteria under combined tension and shear loading are based on maximum circumferential stress or energy criteria such as strain energy density. The development and use of a Dugdale model and finite element models to represent crack and fracture dynamics are examined.

Kobayashi, A. S.↗
NASA NTRS1972

Investigation of the effect of vacuum environment on the fatigue and fracture behavior of 7075-T6.

Axial-load fatigue-life, fatigue-crack propagation, and fracture-toughness experiments were conducted on sheet specimens made of 7075-T6 aluminum alloy. These experiments were conducted at air pressures ranging from 101 kN/sq m to 7 micronewtons/sq m to determine the effect of air pressure on fatigue behavior. Analysis of the results from the fatigue-life experiments indicated that for a given stress level, the lower the air pressure was the longer the fatigue life. At a pressure of 7 micronewtons/sq m, fatigue lives were 15 to 30 times longer than at 101 kN/sq m. Analysis of the results from the fatigue-crack-growth experiments indicates that at low values of stress-intensity range, the fatigue-crack-growth rates were approximately twice as high at atmospheric pressure as in vacuum. However, at higher values of stress-intensity range, the fatigue-crack-growth rates were nominally the same in vacuum and at atmospheric pressure.

Hudson, C. M.↗
DOE OSTI2022

Scaling Behavior of Thermally Driven Fractures in Deep Low‐Permeability Formations: A Plane Strain Model With 1‐D Heat Conduction

Abstract Injection of cold fluids through/into deep formations may cause significant cooling, thermal stress, and possible thermal fracturing. In this study, the thermal fracturing of low‐permeability formations under one‐dimensional heat conduction was investigated using a plane strain model. Dimensionless governing equations, with dimensionless fracture length , aperture , spacing , time , and effective confining stress , were derived. Solution of single thermal fracture was derived analytically, while solution of multiple fractures with constant (or dynamic) spacing were obtained using the displacement discontinuity method (and stability analysis). For single fracture, increases nonlinearly with and then transitions to scaling law , indicating that late‐time fracture length increases linearly with the square root of cooling time. For constantly spaced fractures, deviates from the single‐fracture solution at a later for a larger , showing slower propagation under inter‐fracture stress interaction . For dynamically spaced fractures, fracture arrest induced by stress interaction was determined by the stability analysis; the fully transient solution provides evolution of dimensionless fracture length, spacing, aperture, and pattern; a similar scaling law, with , obtained shows the effect of both stress interaction and fracture arrest. The solution and scaling law provide fast predictions for all reservoir and cooling conditions using (single) model parameter . Application to a geothermal site with demonstrates that thermal fractures reach 0.67, 6.25, and 78.00 m in length, 0.49, 2.30, and 13.00 m in spacing, and 0.43, 2.09, and 12.19 mm in aperture at 1, 100, and 10,000 days.

Chen, Bin↗
NASA NTRS2001

Gear Crack Propagation Path Studies: Guidelines for Ultra-Safe Design

Design guidelines have been established to prevent catastrophic rim fracture failure modes when considering gear tooth bending fatigue. Analysis was performed using the finite element method with principles of linear elastic fracture mechanics. Crack propagation paths were predicted for a variety of gear tooth and rim configurations. The effects of rim and web thicknesses, initial crack locations, and gear tooth geometry factors such as diametral pitch, number of teeth, pitch radius, and tooth pressure angle were considered. Design maps of tooth/rim fracture modes including effects of gear geometry, applied load, crack size, and material properties were developed. The occurrence of rim fractures significantly increased as the backup ratio (rim thickness divided by tooth height) decreased. The occurrence of rim fractures also increased as the initial crack location was moved down the root of the tooth. Increased rim and web compliance increased the occurrence of rim fractures. For gears with constant pitch radii, coarser-pitch teeth increased the occurrence of tooth fractures over rim fractures. Also, 25 deg pressure angle teeth had an increased occurrence of tooth fractures over rim fractures when compared to 20 deg pressure angle teeth. For gears with constant number of teeth or gears with constant diametral pitch, varying size had little or no effect on crack propagation paths.

Lewicki, David G.↗
NASA NTRS1978

Correlations between ultrasonic and fracture toughness factors in metallic materials

A heuristic mathematical basis was proposed for the experimental correlations found between ultrasonic propagation factors and fracture toughness factors in metallic materials. A crack extension model was developed wherein spontaneous stress (elastic) waves produced during microcracking are instrumental in promoting the onset of unstable crack extension. Material microstructural factors involved in the process are measurable by ultrasonic probing. Experimental results indicate that ultrasonic attenuation and velocity measurements will produce significant correlations with fracture toughness properties and also yield strength.

Vary, A.↗
NASA NTRS1978

Correlations between ultrasonic and fracture toughness factors in metallic materials

A heuristic mathematical basis is proposed for the experimental correlations found between ultrasonic propagation factors and fracture toughness factors in metallic materials. A crack extension model is proposed wherein spontaneous stress (elastic) waves produced during microcracking are instrumental in promoting the onset of unstable crack extension. Material microstructural factors involved in this process are measurable by ultrasonic probing. Experimental results indicate that ultrasonic attenuation and velocity measurements will produce significant correlations with fracture toughness properties and also yield strength.

Vary, A.↗
NASA NTRS1974

Roche limit of a solid body

Tidal fission of both impacting and orbiting linear elastic solid bodies based on Kelvin's theory of earth tides is considered. It is shown that there can be more than one mutually exclusive modes of fracture - the particular mode in which a body fractures depending on its size and strength. The analysis gives a vivid picture of the propagation of the fracture with a decreasing distance from the planet. Expressions for the initiation and completion of fracture are obtained which are displayed graphically for a rigid body. The effect of elasticity on the breakup altitude is discussed. For orbiting solid bodies, the study gives the upper limit of the breakup altitude as 0.38R (where R is the radius of planet), which is much less than the value 1.44R used for such bodies in the past. The results presented include a previously given theory by Sekiguchi as a part. For the special case of a liquid body, comparison is made with Roche's calculation and the difference explained.

Aggarwal, H. R.↗
NASA NTRS1996

Gear crack propagation investigations

Analytical and experimental studies were performed to investigate the effect of gear rim thickness on crack propagation life. The FRANC (FRacture ANalysis Code) computer program was used to simulate crack propagation. The FRANC program used principles of linear elastic fracture mechanics, finite element modeling, and a unique re-meshing scheme to determine crack tip stress distributions, estimate stress intensity factors, and model crack propagation. Various fatigue crack growth models were used to estimate crack propagation life based on the calculated stress intensity factors. Experimental tests were performed in a gear fatigue rig to validate predicted crack propagation results. Test gears were installed with special crack propagation gages in the tooth fillet region to measure bending fatigue crack growth. Good correlation between predicted and measured crack growth was achieved when the fatigue crack closure concept was introduced into the analysis. As the gear rim thickness decreased, the compressive cyclic stress in the gear tooth fillet region increased. This retarded crack growth and increased the number of crack propagation cycles to failure.

Lewicki, David G.↗
DOE OSTI2022

Unsupervised learning from three-component accelerometer data to monitor the spatiotemporal evolution of meso-scale hydraulic fractures

Enhanced geothermal systems can provide a substantial share of the global energy demand. There exist several hurdles in the engineering implementations of such geothermal systems. One such hurdle is the accurate monitoring of the fracture networks created in subsurface through hydraulic stimulation of these systems. Micro seismicity associated with the stimulation is the primary means to locate the event hypocenters for estimating the stimulated rock volume. Existing methods for location the hypocenters are restricted to only the highest amplitude impulsive signals that are simultaneously detected on several sensors. Consequently, a large portion (usually ~99%) of the measurements are left unused. In this paper, an unsupervised manifold-approximation followed by clustering of 3-component accelerometer data is used to analyze the seismicity recorded on a monitoring well. With this method, a larger portion of the measured signal is used for the monitoring of the hydraulic fracture network. We analyze the EGS Collab experiment 1 microseismic data, recorded at the Sanford Underground Research Facility, South Dakota. Using the data from a single three-component accelerometer, the polarization features viz. Azimuth, incidence, rectilinearity, and planarity are used as inputs for the unsupervised manifold approximation followed by clustering. Our study shows that density-based clusters in the projected 3D space correspond to distinct types of hydraulically fractured zones around the injection point. Finally, we show that the temporal evolution of these clusters can be used to track fracture creation and propagation.

58 GEOSCIENCES↗
DOE OSTI2024

Arbitrary Order Virtual Element Methods for High‐Order Phase‐Field Modeling of Dynamic Fracture

ABSTRACT Accurate modeling of fracture nucleation and propagation in brittle and ductile materials subjected to dynamic loading is important in predicting material damage and failure under extreme conditions. Phase‐field fracture models have garnered a lot of attention in recent years due to their success in representing damage and fracture processes in a wide class of materials and under a variety of loading conditions. Second‐order phase‐field fracture models are by far the most popular among researchers (and increasingly, among practitioners), but fourth‐order models have started to gain broader acceptance since their more recent introduction. The exact solution corresponding to these high‐order phase‐field fracture models has higher regularity. Thus, numerical solutions of the model equations can achieve improved accuracy and higher spatial convergence rates. In this work, we develop a virtual element framework for the high‐order phase‐field model of dynamic fracture. The virtual element method (VEM) can be regarded as a generalization of the classical finite element method. In addition to many other desirable characteristics, the VEM allows computing on polytopal meshes. Here, we use ‐conforming virtual elements and the generalized‐ time integration method for the momentum balance equation, and adopt ‐conforming virtual elements for the high‐order phase‐field equation. We verify our virtual element framework using classical quasi‐static benchmark problems and demonstrate its capabilities with the aid of numerical simulations of dynamic fracture in brittle materials.

42 ENGINEERING↗
NASA NTRS1989

Kinetics of fracture in Fe-3Si steel under mode I loading

This paper deals with experimental studies of fatigue crack propagation (FCP) in Fe-3Si steel. The FCP experiments were performed on single-edge crack specimens. Results show that an extensive damage zone consisting of slip bands surrounds and precedes the propagating crack. The system of the crack and the damage zone constitutes the crack layer (CL). The results demonstrate that fracture propagates by the translation, expansion, and distortion of the part of CL called the active zone. The contours of the damage distribution show that different loading stresses may cause significantly different shapes of damage. Application of the stability criteria is demonstrated. It was shown that the critical energy release rate depends on the history of loading.

Bessendorf, Michael H.↗