Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fracture Models”

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 1,099 records · Page 61

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↗

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.↗

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.↗

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.↗

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↗

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↗

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↗

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↗

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↗

Statistical analyses of reservoir and fracturing parameters for a multifractured shale oil reservoir in Mississippi

The use of multifractured horizontal wells has improved the efficiency of hydrocarbon extraction from shale gas and oil plays. It is highly desirable to estimate the characteristics of the reservoir and well fracturing through production data analysis. Production rate transient data from 16 wells in the Mississippi sections were analyzed to estimate local reservoir permeability and hydraulic fracture parameters. Key factors affecting well productivity have been identified for optimizing future well completion. Based on the theory of distance of investigation during transient linear flow, the reservoir matrix permeability of TMS was estimated to be between 53 nd and 210 nd, averaging at 116 nd. The matrix permeability follows the lognormal distribution and is considered homogeneous according to the coefficient of variation. The fracture half-length was estimated using the matrix permeability data from the rate transient analysis. The fracture half-length was found to have a mean value of 234 ft with a standard deviation of 66 ft in a normal distribution. The fracture conductivity was back-calculated by matching pseudosteady production rate data to Guo et al's (SPE Reserve Eval Eng.12, 2009, 879) productivity model for boundary-dominated flow. The fracture conductivity was estimated to range from 0.4 md-ft to 3.2 md-ft, averaging at 1.4 md-ft. Based on Guo et al's (SPE Reserve Eval Eng.12, 2009, 879) well productivity model applied to TMS condition, well productivity can be improved significantly by increasing fracture length and conductivity. Future TMS wells should be completed with conductivity values of greater than 10 md-ft.

04 OIL SHALES AND TAR SANDS↗

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↗

Simulating the Fracture of Notched Mortar Beams through Extended Finite-Element Method and Peridynamics

This paper simulates fracture in notched mortar beams under three-point bending using extended finite element method (XFEM) and peridynamics. A three-phase microstructure (i.e., cement paste, aggregates, and paste-aggregate interface) is used for constitutive modeling of the mortar to obtain the elastic properties for simulation. In the XFEM approach, the simulated homogenized elastic modulus is used along with the total fracture energy of the cement mortar in a damage model to predict the fracture response of the mortar including crack propagation and its fracture parameters (Mode I stress intensity factor, K IC and critical crack tip opening displacement, CTOD C ). The damage model incorporates a maximum principal stress-based damage initiation criteria and a traction-separation law for damage evolution. In the peridynamics approach, a bond-based model involving a prototype microelastic brittle (PMB) material model is used. The elastic properties and fracture energy release rates are used as inputs in the PMB model, along with the choice of peridynamic horizon size. Comparison with experimental fracture properties ( K IC , CTOD C ) as well as crack propagation paths from digital image correlation show that both the approaches yield satisfactory results, particularly for K IC and crack extension. Furthermore, both these methods can be adopted for fracture simulation of cement-based materials.

42 ENGINEERING↗

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↗

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↗

The effect of thickness on fatigue crack propagation in 7475-T731 aluminum alloy sheet

Tests were conducted on three thicknesses of 7475-T731 aluminum alloy sheet to investigate the effect of thickness on fatigue crack propagation under constant amplitude loading conditions and on retardation following a single peak overload. Constant amplitude loading tests were performed at stress ratios of 0.05 and 0.75 to obtain data for conditions with crack closure and without crack closure, respectively. At both stress ratios a thickness effect was clearly evident, with thicker specimens exhibiting higher growth rates in the transition from plane strain to plane stress region. The effect of thickness for a stress ratio of 0.05 corresponded well with the fracturing mode transitions observed on the specimens. A model based on the strain energy release rate which accounted for the fracture mode transition was found to correlate the thickness effects well. The specimens tested at the stress ratio of 0.75 did not make the transition from tensile mode to shear mode, indicating that another mechanism besides crack closure or fracture mode transition was active.

Daiuto, R. A.↗

Capacitor Test, Evaluation. and Modeling Within NASA Electronic Parts and Packaging (NEPP) Program. "Why Ceramic Capacitors Fracture During Manual Soldering and How to Avoid Failures"

Presentation discusses: (1) Why Multi-Layer Ceramic Capacitors(MLCCs) crack during manual soldering? Workmanship and parts issues. (2) Do existing qualification requirements assure crack-free soldering? MIL-spec Thermal Shock (TS) testing. MIL-spec Resistance to Soldering Heat (RSH) test. (3) What test can assure reliable soldering? Mechanical characteristics of ceramics. Comparison of three TS techniques: LND, TSD, and IWT. (4) Simulation of TS conditions.

Teverovsky, Alexander↗