Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture modelling”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Fracture and crack growth in orthotropic laminates

A mathematical model based on the classical shear-lag assumptions is used to study the residual strength and fracture behavior of composite laminates with symmetrically placed buffer strips. The laminate is loaded by a uniform remote longitudinal tensile strain and has initial damage in the form of a transverse crack in the parent laminate between buffer strips. The crack growth behavior as a function of material properties, number of buffer-strip plies, spacing, width of buffer strips, longitudinal matrix splitting, and debonding at the interface is studied. Buffer-strip laminates are shown to arrest fracture and increase the residual strengths significantly over those of one material laminates, with S-glass being a more effective buffer strip material than Kevlar in increasing the damage tolerance of graphite/epoxy panels. For a typical graphite/epoxy laminate with S-glass buffer-strips, the residual strength is about 2.4 times the residual strength of an all graphite/epoxy panel with the same crack length. Approximately 50% of this increase is due to the S-glass/epoxy buffer-strips, 40% due to longitudinal splitting of the buffer strip interface and 10% due to bonding.

Goree, James G.↗
DOE OSTI2021

Mechanism of Microseismic Generation During Hydraulic Fracturing - With Evidence From HFTS 2 Observations

The objective of this study is to understand how microseismic events are generated during hydraulic fracturing, as well as the role of geomechanical conditions (i.e., stress and mechanical stratigraphy) in this process. In the industry, microseismic event clouds have been generally used as an “outer-boundary” of the “stimulated reservoir volume” (SRV). However, by comparing with other surveillance data (Low frequency Distributed Acoustic Sensing, or LF-DAS strain) in the Hydraulic Fracturing Test Site (HFTS) 2 experiment, we show that this assumption is fundamentally flawed. The HFTS 2 data has three unique observations that have not been commonly observed in other datasets: 1. Due to influence of offset pad depletion, microseismic data shows that hydraulic fractures from the child well can propagate over 3000 feet into the depleted low stress zone. 2. By comparing microseismic and horizontal fiber LF-DAS strain data, we observe that microseismic event cloud does not necessarily reflect the created hydraulic fracture volume. Particularly, the extent of microseismic event clouds near heel stages are much shorter than what is shown with LF-DAS strain data. 3. Microseismic event magnitudes are larger in the depleted regions. Through geomechanical analysis, we demonstrate that the “bedding-plane-slip” model is likely the mechanism for microseismic generation during hydraulic fracturing. Furthermore, this model successfully explains the above field observations from HFTS 2 experiment. We also provide a quantitative relationship connecting the microseismic event magnitude with fracture width increment and layer mechanical property contrast.

58 GEOSCIENCES↗
DOE OSTI2026

In-Situ Visualization of a Growing Brittle Crack in Aluminum Oxynitride Using Synchrotron X-Rays and the Double-Cleavage Drilled Compression Geometry

Brittle fracture is difficult to study in situ due to the speed of a growing crack and the often-catastrophic nature of failure in brittle materials. As a result, the influence of microstructural considerations, such as orientation, grain boundary locations, and strain field, on the crack path remains poorly understood. Presented in this study is a method addressing this knowledge gap, which utilizes the double-cleavage drilled compression geometry to achieve quasi-stable fracture in aluminum oxynitride (AlON). Synchrotron X-ray micro-computed tomography is used to characterize the crack shape and length, while high-energy diffraction microscopy provides information on the strains, orientations, and shapes of grains in the microstructure surrounding the crack tip. During testing, the crack grew in discrete and irregular jumps while the fracture toughness falls within reported ranges. The crack in AlON is found to have no greater tendency to crack intergranularly as compared to transgranularly, and grains which are cracked transgranularly do not display a trend in orientation or stress when compared to those around which the crack followed a grain boundary. The high resolution of the crack path and microstructural data provides a path forward for modeling and understanding 3D brittle fracture.

Gorske, Sara F.↗
DOE OSTI2024

Investigating the Crosslinking, Degradation, and Adhesion Behavior of Photovoltaic Encapsulants Under Thermal Accelerated Aging

Degradation of photovoltaic (PV) module encapsulant characteristics that lead to mechanical embrittlement and delamination remains a cause of failure in solar installations. A multiscale reliability model connecting the encapsulant mechanical and fracture properties to the degraded molecular structure and interfacial bonding to adjacent solar cell and glass substrates was previously published. The model, developed primarily for poly(ethylene-co-vinyl acetate) acetate (EVA) encapsulants, remains to be experimentally validated. Determining the degradation and crosslinking kinetics of alternative encapsulants, such as polyolefin elastomer (POE) and EVA/POE/EVA composites (EPE), can generalize the model. In this work, we subject fully cured EVA, POE, and EPE encapsulants to accelerated thermal aging to determine how high temperatures impact reaction kinetics. An increase in gel content (crosslinking) and decrease in crystallinity of the encapsulants under hot-aerobic (90 degrees C, 22% RH) and hot-anaerobic (90 degrees C, sealed in N 2 air) aging were observed, even in the absence of UV and crosslinking initiators. Fourier transform infrared spectroscopy (FTIR)-attenuated total reflectance analysis showed insignificant encapsulant degradation, demonstrating the critical role of UV and moisture in accelerating degradation. Adhesion testing performed on coupon-level specimens (cell/encapsulant/glass laminates) showed decreases in adhesion energy, Gc , from 5000 h of hot-dry (90 degrees C, ~1% RH) and hot-humid (90 degrees C, 60% RH) aging. POE coupons demonstrated the best stability, followed by EPE then EVA. For EVA and POE, hot-humid aged coupons experienced a larger decrease in Gc due to enhanced hydrolytic degradation. Hot-dry aging condition demonstrated that thermal degradation of the interface could be significant even if the encapsulant experiences negligible degradation in the absence of UV and elevated humidity.

accelerated aging↗
NASA NTRS1989

Fatigue crack growth model RANDOM2 user manual, appendix 1

The FORTRAN program RANDOM2 is documented. RANDOM2 is based on fracture mechanics using a probabilistic fatigue crack growth model. It predicts the random lifetime of an engine component to reach a given crack size. Included in this user manual are details regarding the theoretical background of RANDOM2, input data, instructions and a sample problem illustrating the use of RANDOM2. Appendix A gives information on the physical quantities, their symbols, FORTRAN names, and both SI and U.S. Customary units. Appendix B includes photocopies of the actual computer printout corresponding to the sample problem. Appendices C and D detail the IMSL, Ver. 10(1), subroutines and functions called by RANDOM2 and a SAS/GRAPH(2) program that can be used to plot both the probability density function (p.d.f.) and the cumulative distribution function (c.d.f.).

Boyce, Lola↗
NASA NTRS1992

A creep cavity growth model for creep-fatigue life prediction of a unidirectional W/Cu composite

A microstructural model was developed to predict creep-fatigue life in a (0)(sub 4), 9 volume percent tungsten fiber-reinforced copper matrix composite at the temperature of 833 K. The mechanism of failure of the composite is assumed to be governed by the growth of quasi-equilibrium cavities in the copper matrix of the composite, based on the microscopically observed failure mechanisms. The methodology uses a cavity growth model developed for prediction of creep fracture. Instantaneous values of strain rate and stress in the copper matrix during fatigue cycles were calculated and incorporated in the model to predict cyclic life. The stress in the copper matrix was determined by use of a simple two-bar model for the fiber and matrix during cyclic loading. The model successfully predicted the composite creep-fatigue life under tension-tension cyclic loading through the use of this instantaneous matrix stress level. Inclusion of additional mechanisms such as cavity nucleation, grain boundary sliding, and the effect of fibers on matrix-stress level would result in more generalized predictions of creep-fatigue life.

Kim, Young-Suk↗
NASA NTRS2010

Fracture Mechanics Analyses of Reinforced Carbon-Carbon Wing-Leading-Edge Panels

Fracture mechanics analyses of subsurface defects within the joggle regions of the Space Shuttle wing-leading-edge RCC panels are performed. A 2D plane strain idealized joggle finite element model is developed to study the fracture behavior of the panels for three distinct loading conditions - lift-off and ascent, on-orbit, and entry. For lift-off and ascent, an estimated bounding aerodynamic pressure load is used for the analyses, while for on-orbit and entry, thermo-mechanical analyses are performed using the extreme cold and hot temperatures experienced by the panels. In addition, a best estimate for the material stress-free temperature is used in the thermo-mechanical analyses. In the finite element models, the substrate and coating are modeled separately as two distinct materials. Subsurface defects are introduced at the coating-substrate interface and within the substrate. The objective of the fracture mechanics analyses is to evaluate the defect driving forces, which are characterized by the strain energy release rates, and determine if defects can become unstable for each of the loading conditions.

Raju, Ivatury S.↗
NASA NTRS1985

The accuracy of approximate solutions in the analysis of fracture of composites

This paper concerns the accuracy of three related mathematical models (developed by Hedgepeth, Eringen and Sendeckyj and Jones) used in the stress analysis and in fracture studies of continuous-fiber composites. These models have particular application in the investigation of fiber and matrix stresses in unidirectional composites in the region near a crack tip. The interest in such models is motivated by the desire to be able to simplify the equations of elasticity to the point that they can be solved in a relatively easy manner.

Goree, J. G.↗
NASA NTRS1985

The accuracy of approximate solutions in the analysis of fracture of composites

This paper concerns the accuracy of three related mathematical models (developed by Hedgepeth, Eringen and Sendeckyj and Jones) used in the stress analysis and in fracture studies of continuous-fiber composites. These models have particular application in the investigation of fiber and matrix stresses in unidirectional composites in the region near a crack tip. The interest in such models is motivated by the desire to be able to simplify the equations of elasticity to the point that they can be solved in a relatively easy manner.

Goree, J. G.↗
DOE OSTI2026

Machine Learning an Ab-Initio Based Bond-Order Potential for Bismuthene

Bismuthene is a heavy 2D material whose strong spin–orbit coupling and recently observed single-element ferroelectricity have intensified interest in its structural, vibrational, and transport properties. Accurate modeling of these behaviors requires a short-range interatomic potential that can reproduce the underlying bonding physics at a fraction of the computational cost of first-principles methods. However, such a potential is currently unavailable. Here, in this work, we construct a Tersoff bond-order potential for β-bismuthene using a reinforcement-learning framework that integrates a continuous Monte Carlo Tree Search with a simplex-based local optimizer. The optimized parameter sets reproduce first-principles lattice constants, cohesive energy, the equation of state, elastic constants, and phonon dispersion. We validate the models by performing thermal-conductivity calculations and uniaxial fracture simulations our findings confirm the reliability of the resulting models across multiple thermomechanical regimes. Comparison of the three best solutions reveals how differences in pairwise interactions, angular terms, and bond-order behavior govern phonon features and mechanical responses. We demonstrate an interpretable and computationally efficient potential for bismuthene and demonstrate a general reinforcement-learning strategy for developing bond-order models in emerging 2D materials.

deformation↗
NASA NTRSDate not supplied

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials (3MDCP)

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Woven↗
NASA NTRSDate not supplied

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Discrete Element Method↗
DOE OSTI2024

Utah FORGE 5-2615: Determination and Analysis of Thermo-poromechanical Response of Fractured Rock - 2024 Annual Workshop Presentation

This is a presentation on the Determination and Modeling-Informed Analysis of Thermo-poromechanical Response of Fractured Rock for Application to FORGE by the University of Oklahoma, presented by Ahmad Ghassemi. This video presentation discusses how to improve understanding and control of coupled thermoporomechanical (or thermo-hydro-mechanical-THM) processes in reservoir development, and to demonstrate its role in interpretations of the fracture closure pressure. This presentation was featured in the Utah FORGE R&D Annual Workshop on August 13-15, 2024.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Simulating Fracture Sealing by Granular LCM Particles in Geothermal Drilling

Lost circulation occurs when the returned fluid is less than what is pumped into the well due to loss of fluid to pores or fractures. A lost-circulation event is a common occurrence in a geothermal well. Typical geothermal reservoirs are often under-pressured and have larger fracture apertures. A severe lost-circulation event is costly and may lead to stuck pipe, well instability, and well abandonment. One typical treatment is adding lost-circulation materials (LCMs) to seal fractures. Conventional LCMs fail to properly seal fractures because their mechanical limit is exceeded at elevated temperatures. In this paper, parametric studies in numerical simulations are conducted to better understand different thermal effects on the sealing mechanisms of LCMs. The computational fluid dynamics (CFDs) and the discrete element method (DEM) are coupled to accurately capture the true physics of sealing by granular materials. Due to computational limits, the traditional Eulerian–Eulerian approach treats solid particles as a group of continuum matter. With the advance of modern computational power, particle bridging is achievable with DEM to track individual particles by modeling their interactive forces between each other. Particle–fluid interactions can be modeled by coupling CFD algorithms. Fracture sealing capability is investigated by studying the effect of four individual properties including fluid viscosity, particle size, friction coefficient, and Young’s modulus. It is found that thermally degraded properties lead to inefficient fracture sealing.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Hybrid fiber metal composite laminate interlaminar reinforcement through metal interlocks

Aircraft and automobile industries are continually seeking high-performance and lightweight solutions. Fiber-reinforced composites in conjunction with metals are being considered in the form of hybrid materials. There is a continuous need for improvement of interfacial bonding between composite layers and metal constituents. In this study, an innovative hybrid fiber interlocking metal hooks laminate (FIMHL) system was developed in which fiber-reinforced thermoset composite, and metal sheets are mechanically bonded together using out of plane hooks stamped in the metal. Process optimization was performed to gain the full benefit of the through-thickness reinforcement. Microstructural analysis showed improved interaction between the metal hooks and the fiber layers which reduced porosity and resin richness. This was reflected in mechanical properties, as tensile and flexural strength of FIMHL was enhanced by 38.5 and 18.8%, respectively, after process optimization. Furthermore, normalized weight fraction (76.4%) properties of FIMHL also confirmed that the increase in mechanical properties of optimized panel were not only due to increased reinforcement (glass fiber + aluminum) volume fraction but also because of improved metal hooks and fibers interaction. There was no significant effect on lap shear and fracture toughness, as they depend on bend back behavior of the hooks. A bilinear traction-separation model was used to characterize mode-I interlaminar fracture toughness properties, and 4.5% variation was recorded when modeling and experimental values were compared.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Skin factor and potential formation damage from chemical and mechanical processes in a naturally fractured carbonate aquifer with implications to CO 2 sequestration

Here, in this study, we investigate formation damage due to acidization and water injection tests into the naturally fractured carbonate Middle Duperow Formation at Kevin Dome, Montana, potentially diminishing the chance of a future successful Geological Carbon Sequestration (GCS) project. Multiple well-test analytical models, correlated with core description and lithology data, are used to determine flow behavior and communication between the water injection interval and surrounding formations. An improved three-dimensional (3D) geologic model with dual-continuum matrix and fracture properties is constructed based on most recent seismic, core, and water sample measurements. Brine injection is simulated to verify the interpretation from the analytical models, followed by CO 2 injection simulation. Geochemical calculations are performed to understand the in-situ processes that led to formation damage. Our findings suggest: (1) there are two possible scenarios that could lead to a positive total effective skin factor and permeability decline: partial penetration and formation damage; (2) analytical models indicate a positive total skin factor, contradicting results of a previous study suggesting that the well was mildly stimulated; (3) numerical simulation supports the formation damage hypothesis by matching the pressure buildup observed during the latter two brine injection tests; (4) several mechanical and chemical processes may have occurred during injection to clog the matrix/fracture system: anhydrite fines migration and/or calcite precipitation. We then make preventative suggestions for future GCS projects into carbonate reservoirs and remediation recommendations for GCS operation at the Kevin Dome.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Develop a new integrated macro→micro←nano (MMN) multiscale modeling framework to optimize high strength aluminum alloys and processes for vehicle light-weighting​

Bending tests provide a means to study plane strain fracture performance of 6000 series aluminum alloys. Metrics from bending tests have been correlated with self-pierce riveting (SPR) performance of a high strength AA6111 automotive aluminum alloy in previous works. Using the ORNL HPC resources, this project developed an innovative macro→micro←nano (MMN) multiscale microstructure-based finite element (FE) code to further understanding of the relationship between microstructure and fracture properties of high-strength 6000-series alloys. This work started with microstructural characterization in both mesoscale and nanoscale and bending performance characterization of AA6111 HS2-T6 alloy at Ford, the MMN framework was applied to this alloy to simulate 3-point VDA bending. From the results of the macro-modeling of 3-point VDA bending, the critical region of fracture was identified, and the region geometry was used to construct the micro-model. The fracture criterion of micron-scale precipitates and aluminum matrix which contains submicron and nano particles (AL-SMP-NP) in the micro-model was calibrated and validated by comparing simulated and measured bending results. With the AL-SMP-NP fracture strain obtained, the fracture strain of Al-matrix containing nano particles (ALNP) will similarly be determined by a submicron scale-model using an edge-constrained FE modeling approach developed by Hu et al. With the ALNP fracture strain obtained, the fracture strain of Al-matrix containing no particles will similarly be determined by a nano scale-model using an edge-constrained FE modeling approach. After the MMN framework is built and fracture criterion calibrated, nano-model FE simulations with virtual microstructures was performed to obtain a reduced order model (ROM) of the fracture criterion of the ALNP as a function of volume fraction, size, and distribution of the nanoparticle. This nano→submicron→micro modeling part allows exploration of the influence of different material nanostructures from different process conditions on the bending properties within the multiscale bending simulation framework and the ROM of material bendability as a function of nanoparticle size and shape was established. This obtained reduced order model (ROM) could help guide the design and selection of materials to improve existing SPR process models that could replace trial-and-error rivet/die selection and help to design new rivet/die combinations capable of robustly joining new higher strength 6000 5 series alloys in automotive body structures. This would enable lightweighting of Ford vehicles leading to greater fuel efficiency and reduce manufacturing time and energy.

36 MATERIALS SCIENCE↗