Search NASA⌕ Search

SEARCH · Search NASA

Results for “Brittle fracture”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

FFT-based solver for higher-order and multi-phase-field fracture models applied to strongly anisotropic brittle materials

This paper presents the application of a fast Fourier transform (FFT) based method to solve two phase field models designed to simulate crack growth of strongly anisotropic materials in the brittle regime. By leveraging the ability of the FFT-based solver to generate solutions with higher-order and global continuities, we design two simple algorithms to capture the complex fracture patterns (e.g. sawtooth, and curved crack growth) common in materials with strongly anisotropic surface energy via the multi-phase-field and high-order phase-field frameworks. A staggered operator-split solver is used where both the balance of linear momentum and the phase field governing equations are formulated in the periodic domain. The unit phase field of the initial failure region is prescribed by the penalty method to alleviate the sharp material contrast between the initial failure region and the base material. The discrete frequency vectors are generalized to estimate the second and fourth order gradients such that the Gibbs effect near shape interfaces or jump conditions can be suppressed. Furthermore, a preconditioner is adopted to improve the convergence rate of the iterative linear solver. Three numerical experiments are used to systematically compare the performance of the FFT-based method in the multi-phase-field and high-order phase-field frameworks.

42 ENGINEERING↗
DOE OSTI2025

Influence of build direction on the fracture mechanism of 3D printed octet lattices

Here, we investigate the effects of 3D printing build directions on the fracture properties of octet lattice metamaterials made of polylactic acid (PLA) and how these effects vary with the relative density of the lattices. Single-edge notch bend samples are 3D printed in two orthogonal build directions at various relative densities. Our results show that the work of fracture for octet lattices with build directions parallel to the crack plane is significantly higher than those with build directions perpendicular to the crack plane. We also observed that the ratio of specific work of fracture between the two build directions remains nearly constant across different relative densities. In contrast, the ratio of peak load between the two build directions decreases as relative density increases. Furthermore, the build direction dictates the fracture mechanism. While the perpendicular build direction predominantly results in a brittle (mode I) fracture, the parallel build direction leads to a complex, delamination-dominated (mode II type) fracture. This phenomenon is largely governed by the weak interfaces formed between the printed layers and their interaction with the lattice geometry. These results reveal that the build direction governs the fracture mechanism and work of fracture of these lattice metamaterials and is therefore an important design consideration.

3D printed PLA↗
DOE OSTI2021

Anisotropic thermal-conductivity degradation in the phase-field method accounting for crack directionality

Dynamic loads applied to metals may lead to brittle or ductile fracture depending on the imposed strain rates, material properties, and specimen geometry. With an increase in the applied velocity, brittle to ductile failure mode transition may also be observed. Furthermore, at high strain rates, ductile fracture may be preceded by shear bands, which are narrow bands of intense plastic deformation, typically accompanied by a significant rise in temperature. Reliable models are needed to predict the response of metals subject to dynamic loads. In particular, capturing the interplay between heat conduction and crack propagation using the phase-field fracture method is still an open research field. To accurately capture the heat transfer physics across crack surfaces, damage models degrading thermal-conductivity are necessary. While isotropic thermal-conductivity degradation models were proposed, they may lead to errors as they are indifferent to crack directionality. In our current work, a new anisotropic approach is proposed in which thermal-conductivity, which depends on the phase-field gradient, is degraded solely across the crack. Here, it is shown that this approach improves the near-field approximation of temperature and heat flux compared with isotropic degradation, when taking the discontinuous crack solutions as reference. Additionally, the anisotropic degradation approach is implemented in a unified dynamic fracture model and a couple of examples demonstrate the viability of this approach.

42 ENGINEERING↗
DOE OSTI2020

Scale bridging damage model for quasi-brittle metals informed with crack evolution statistics

Computationally efficient methods for bridging length scales, from highly resolved micro/meso-scale models that can explicitly model crack growth, to macro-scale continuum models that are more suitable for modeling large parts, have been of interest to researchers for decades. In this work, an improved brittle damage model is presented for the simulation of dynamic fracture in continuum scale quasi-brittle metal components. Crack evolution statistics, including the number, length, and orientation of individual cracks, are extracted from high-fidelity, finite discrete element method (FDEM) simulations and used to generate effective material moduli that reflect the material’s damaged state over time. This strategy allows for the retention of small-scale physical behaviors such as crack growth and coalescence in continuum scale hydrodynamic simulations. However, the high-fidelity simulations required to generate the crack statistics are computationally expensive. Thus, steps were taken to produce a flexible constitutive model to reduce the number of costly high-fidelity simulations needed to produce accurate results. A new stress based degradation criterion is introduced for the degradation of individual material zones. This allows for the development of a heterogeneous damage distribution within the bulk material. Then a flow stress model is added to the hydrodynamic simulation to account for plasticity in quasi-brittle materials. As a result, the effective moduli model can be applied to a larger range of materials. The effective moduli constitutive model is used to simulate beryllium flyer plate experiments. The results from the continuum scale simulations using statistics from a single high-fidelity simulation are found to be in excellent agreement with numerical and experimental velocity interferometer data. The same set of crack statistics are used to extrapolate the results of a higher rate flyer plate case using the effective moduli model. In conclusion, the extension of this model to higher rate cases shows promise for further reducing the number of costly high-fidelity simulations needed to generate crack statistics.

36 MATERIALS SCIENCE↗
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↗
DOE OSTI2025

Effect of Loading Frequency on Fatigue Behavior of Additively Manufactured Short Fiber Thermoplastics via Infrared Thermography

Fracture, Fatigue, Failure, Damage Evolution and Thermomechanics & Infrared Imaging, Volume 4 of the Proceedings of the 2025 SEM Annual Conference & Exposition on Experimental and Applied Mechanics, the fourth volume of five from the Conference, brings together contributions to this important area of research and engineering. The collection presents early findings and case studies on a wide range of areas, including: • 3D Deformation and Failure • In-situ Techniques and Microscale Effects on Mechanical Behavior • Vibration and Cyclic Fatigue • Integration of Models and Experiments • Interface-Mediated Damage and Failure • Fracture and Fatigue in Brittle Materials • Material Characterizations Using Thermography • Fatigue, Damage & Fracture Evaluation Using Infrared Thermography

Pathak, Pharindra [Auburn University]↗
DOE OSTI2024

Stress Dependency of Brittle Creep in Granite: Insights into Source Mechanisms and Parameters

Creep in rocks refers to the gradual deformation of rock material over time under the influence of constant stress. Characterizing these deformations is of great importance for engineering design, geotechnical assessment, mining operations, geological studies, and understanding natural hazards. While laboratory experiments and a variety of numerical approaches have offered explanations for microcrack interaction and damage accumulation under the three stages of creep (primary, secondary and tertiary) in conventional creep experiments, the micromechanisms of the fractures produced in brittle creep and its dependency on the applied stress have not been explored in detail. The present study focused on investigating the fracturing mechanisms that occur during creep-induced fracturing at different stress levels and estimation of the source parameters and energy budget components. A series of uniaxial compression creep experiments have been conducted at different stress level ratios (70%, 75%, 80% and 85%), to the unconfined compressive strength (UCS) of double-flawed Barre granite specimen. Creep measurements were complemented with the Acoustic Emission (AE) measurements. The creep-induced fractures were classified into double-couple (DC), compensated linear vector dipole (CLVD) and isotropic (ISO) components using the AE moment tensor decomposition technique. The results show that non-double-couple sources dominated during creep at all the specified stress levels; however, their proportions decreased as the stress level was increased. The source parameters estimation indicated a significant increase in the magnitude of the events and the radiated seismic energy with increasing levels of stress and a slight increase in the evaluated source radius and stress drop with increasing stress levels. Furthermore, this study contributes to the existing knowledge of creep-induced fracturing by providing insights into the fracturing mechanisms and the radiated seismic energy produced, which can be helpful for the development of improved models and strategies for rock engineering and geoscience applications.

58 GEOSCIENCES↗
DOE OSTI2021

Nanomechanical investigation of the interplay between pore morphology and crack orientation of amorphous silica

Porous amorphous silica (a-SiO 2 ) is of both fundamental and practical interests, as they exhibit a large specific surface area and tunable porous network. However, the brittle nature of a-SiO 2 and the presence of pre-existing cracks at both micro- and nano-scales lead to complex mechanical behavior. Here, we systematically investigate the effects of pre-existing crack and its orientation on the mechanical properties of a-SiO 2 with varying pore shapes using reactive molecular dynamics simulations. We demonstrate that pore shape will primarily influence the Young’s modulus (E) and critical energy release rate (G IC ). We further investigate the impact of pore shape and crack orientation by local characterization of the structural parameters. By defining the high stress and inter-mediate regions, the overall mechanical properties are found to be greatly influenced by the pore shape which can be reflected through the spatial distribution of von Mises stress. Overall, G IC is found to increase with the increase of ligament length (also known as pore wall thickness). Meanwhile, the effect of the pre-existing crack on the crack propagation process is confirmed by analyzing the density distribution evolution. These results highlight the interplay between pore morphology and crack orientation in controlling the fracture behaviors in brittle porous materials.

42 ENGINEERING↗
DOE OSTI2024

Strain Field Around Individual Dislocations Controls Failure

Understanding material failure on a fundamental level is a key aspect in the design of robust structural materials, especially for metals and alloys capable to undergo plastic deformation. In the last decade, significant progress is made in quantifying the stresses associated with failure in both experiments and simulations. Nonetheless, the processes occurring on the most essential level of individual dislocations that govern semi-brittle and ductile fracture are still experimentally not accessible, limiting the failure prediction capabilities. Therefore, in the present work, a one-of-a-kind nanoscale fracture experiment is conducted on a single crystalline Cr bending beam in situ in the transmission electron microscope and for the first time quantify the transient strains around individual dislocations, as well as of the whole dislocation network during crack opening. The results reveal the importance of both pre-existing and newly emitted dislocations for crack-tip shielding via their intrinsic strain field and provide guidelines to design more damage tolerant materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ultra-grain refinement creates FCC pure cobalt with high strength and high ductility

Although pure cobalt is generally known to have a hexagonal close-packed (HCP) structure at room temperature, we show that its high-temperature face-centered cubic (FCC) phase can be strongly stabilized through grain refinement, resulting in FCC pure cobalt at room temperature. Ultrafine-grained (UFG) FCC cobalt exhibits a hierarchical microstructure consisting of dense stacking fault (SF) networks in the dominant FCC grains and numerous SFs and thin FCC layers within a few HCP plates. This unique microstructure leads to a high tensile strength exceeding 1 GPa, together with a tensile elongation of over 35%, thereby surpassing the well-known strength–ductility trade-off of pure metals. In-situ synchrotron X-ray diffraction revealed that the UFG FCC cobalt exhibited a markedly enhanced deformation-induced FCC→HCP martensitic transformation, which provided sustained strain hardening through the transformation-induced plasticity (TRIP) effect. Furthermore, ultra-grain refinement dramatically suppressed premature void and crack formation, causing a transition in the fracture mode from brittle to ductile. These findings advance the fundamental understanding of the phase stability and TRIP-assisted deformation in elemental cobalt and offer new guidelines for the microstructure-driven design of high-performance cobalt-based structural alloys.

Metastability↗
DOE OSTI2025

Impact of post-irradiation annealing on mechanical performance of irradiated 718 alloy

Here, the effect of post-irradiation annealing on solution-annealed 718 alloy was investigated using advanced mechanical testing, fractography, scanning electron microscopy, and transmission electron microscopy. Specimens were extracted from a proton beam window operated at the Spallation Neutron Source, irradiated with 940 MeV protons to a maximum dose of approximately 9.7 displacements per atom (dpa) at a calculated temperature not exceeding 110 °C while in service. Helium and hydrogen concentrations reached about 1700 and 6900 atomic parts per million (appm), respectively. Despite irradiation and high tensile strength (yield stress over 1 GPa), the material exhibited significant ductility. Annealing at 500 °C, 700 °C, and 900 °C for 30 min resulted in an appreciable decrease in yield strength and an increase in ductility for annealing treatments at 500 °C and 900 °C relative to the strength and ductility of the as-irradiated material. The presence of helium and hydrogen led to cavity formation and cleavage-like brittle features on fractured surfaces; however, high-magnification imaging revealed the presence of small-scale ductile dimples, indicating that the fracture mechanism remained mixed. The annealed specimens retained total elongation levels of 14–33 %, and there was no sudden drop in ductility after heat treatments. The ductility level in the irradiated and annealed material is notable despite the presence of helium and hydrogen.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Toughening materials: enhancing resistance to fracture

It has been said that ‘God invented plasticity, but the Devil invented fracture!’ Both mechanisms represent the two prime modes of structural failure, respectively, plastic collapse and the rupture/breaking of a component, but the concept of developing materials with enhanced resistance to fracture can be difficult. This is because fracture resistance invariably involves a compromise—between strength and ductility, between strength and toughness—fundamentally leading to a ‘conflict’ between nano-/micro-structural damage and the mechanisms of toughening. Here, we examine the two major classes of such toughening: (i) intrinsic toughening, which occurs ahead of a crack tip and is motivated by plasticity—this is the principal mode of fracture resistance in ductile materials, and (ii) extrinsic toughening, which occurs at, or in the wake of, a crack tip and is associated with crack-tip shielding—this is generally the sole mode of fracture resistance in brittle materials. We briefly examine how these distinct mechanistic processes have been used to toughen synthetic materials—intrinsically in gradient materials and in multiple principal-element metallic alloys with the example of metallic glasses and high-entropy alloys, and extrinsically in ceramics with the example of ceramic-matrix composites—in comparison to Nature which has been especially adept in creating biological/natural materials which are toughened by one or both mechanistic classes, despite often consisting of constituents with meagre mechanical properties. The success of Nature has been driven by its ability to cultivate the development of materials with multiple length-scale hierarchical structures that display ingenious gradients and structural adaptability, a philosophy which we need to emulate and more importantly learn to synthesize to make structural materials of the future with unprecedented combinations of mechanical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A chain stretch-based gradient-enhanced model for damage and fracture in elastomers

Similar to quasi-brittle materials, it has been recently shown that elastomers can exhibit a macroscopically diffuse damage zone that accompanies the fracture process. In this study, we introduce a stretch-based gradient-enhanced damage (GED) model that allows the fracture to localize and also captures the development of a physically diffuse damage zone. This capability contrasts with the paradigm of the phase field method for fracture, where a sharp crack is numerically approximated in a diffuse manner. Capturing fracture localization and diffuse damage in our approach is achieved by considering nonlocal effects that encompass network topology, heterogeneity, and imperfections. These considerations motivate the use of a statistical damage function dependent upon the nonlocal deformation state. From this model, fracture toughness is realized as an output. While GED models have been classically utilized for damage modeling of structural engineering materials (e.g., concrete), they face challenges when trying to capture the cascade from damage to fracture, often leading to damage zone broadening (de Borst and Verhoosel, 2016). This deficiency contributed to the popularity of the phase-field method over the GED model for elastomers and other quasi-brittle materials. Other groups have proceeded with damage-based GED formulations that prove identical to the phase-field method (Lorentz et al., 2012), but these inherit the aforementioned limitations. To address this issue in a thermodynamically consistent framework, we implement two modeling features (a nonlocal driving force bound and a simple relaxation function) specifically designed to capture the evolution of a physically meaningful damage field and the simultaneous localization of fracture, thereby overcoming a longstanding obstacle in the development of these nonlocal strain- or stretch-based approaches. Here, we discuss several numerical examples to understand the features of the approach at the limit of incompressibility, and compare them to the phase-field method as a benchmark for the macroscopic response and fracture energy predictions.

Elastomers↗
DOE OSTI2023

Fracture of PC and PMMA cylinders in the Taylor tests

The present paper analyzes failure observed in polymer rods collided with a steel anvil under the Taylor cylinder impact experiments. The experimental and numerical analysis demonstrates that association of the failure patterns with the loading mode depends on the rod material. Two materials are considered in the paper, specifically, polymethylmethacrylate (PMMA) and polycarbonate (PC). Damage in PMMA rods at the Taylor tests is shown to be dominated by the shear load resulting in plastic deformation. In contrast, the major mechanism of damage in PC rods is governed by omni-directional (pressure) mode of load. Analysis of the damage in sections of the rods associated with the glass transition is also addressed in the present theoretical study. Failure in these polymers follows the traditional fracture wave mechanism for brittle glass-type materials due to load in the corresponding mode. Numerical analysis employs a two-phase model with improved inter-phase thermal balance accounting for the inter-phase mass, momentum, and heat exchanges during the glass phase transition, while damage being accumulated in time. In conclusion, the simulation results demonstrate a good correlation with the experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Irradiation damage concurrent challenges with RAFM and ODS steels for fusion reactor first-wall/blanket: a review

Reduced activation ferritic martensitic (RAFM) and oxide dispersion strengthened (ODS) steels are the most promising candidates for fusion first-wall/blanket (FW/B) structures. The performance of these steels will deteriorate during service due to neutron damage and transmutation-induced gases, such as helium/hydrogen, at elevated operating temperatures. Here, after highlighting the operating conditions of fusion reactor concepts and a brief overview, the main irradiation-induced degradation challenges associated with RAFM/ODS steels are discussed. Their long-term degradation scenarios such as (a) low-temperature hardening embrittlement (LTHE)—including dose-temperature dependent yield stress, tensile elongations, necking ductility, test temperature effect on hardening, Charpy impact ductile-to-brittle transition temperature and fracture toughness, (b) intermediate temperature cavity swelling, (c) the effect of helium on LTHE and cavity swelling, (d) irradiation creep and (e) tritium management issues are reviewed. The potential causes of LTHE are discussed, which highlights the need for advanced characterisation techniques. The mechanical properties, including the tensile/Charpy impact of RAFM and ODS steels, are compared to show that the current generation of ODS steels also suffers from LTHE, and shows irradiation hardening up to high temperatures of ~400 °C–500 °C. To minimise this, future ODS steel development for FW/B-specific application should target materials with a lower Cr concentration (to minimise α'), and minimise other elements that could form embrittling phases under irradiation. RAFM steel-designing activities targeting improvements in creep and LTHE are reviewed. The need to better understand the synergistic effects of helium on the thermo-mechanical properties in the entire temperature range of FW/B is highlighted. Because fusion operating conditions will be complex, including stresses due to the magnetic field, primary loads like coolant pressure, secondary loads from thermal gradients, and due to spatial variation in damage levels and gas production rates, an experimentally validated multiscale modelling approach is suggested as a pathway to future reactor component designing such as for the fusion neutron science facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Jeankempite, Ca 5 (AsO 4 ) 2 (AsO 3 OH) 2 (H 2 O) 7 , a new arsenate mineral from the Mohawk Mine, Keweenaw County, Michigan, USA

Abstract Jeankempite, Ca 5 (AsO 4 ) 2 (AsO 3 OH) 2 (H 2 O) 7 , is a new mineral species (IMA2018-090) discovered amongst coatings of arsenate minerals on oxidised copper arsenides from the Mohawk No. 2 mine, Mohawk, Keweenaw County, Michigan, USA. The new mineral occurs as lamellar bundles of colourless to white plates up to 1 mm wide and is visually indistinguishable from guérinite, with which it forms intergrowths. Jeankempite is transparent to translucent with a waxy lustre and white streak, is non-fluorescent under longwave and shortwave ultraviolet illumination, has a Mohs hardness of ~1.5 and brittle tenacity with uneven fracture. Crystals are flattened on {01 $\bar{1}$ } and exhibit perfect cleavage on {01 $\bar{1}$ }. Optically, jeankempite is biaxial (+), α = 1.601(2), β = 1.607(2), γ = 1.619(2) (white light); 2V meas. = 72(2)° and 2V calc. = 71.0°. The empirical formula is (Ca 4.97 Na 0.013 Mg 0.017 )(As 3.99 S 0.01 ) 4 O 23 H 16 , based on 23 O and 16 H atoms per formula unit. Thermogravimetric analysis indicates that jeankempite undergoes four weight losses totalling 16.82%, close to the expected loss of 16.30%, corresponding to eight H 2 O. Jeankempite is triclinic, P $\bar{1}$ , a = 6.710(6), b = 14.901(14), c = 15.940(15) Å, α = 73.583(12)°, β = 81.984(12)°, γ = 82.754(12)°, V = 1507(2) Å 3 and Z = 3. The final structure was refined to R 1 = 0.0591 for 2781 reflections with I obs > 3σ I . The crystal structure of jeankempite is built from a network of edge- and vertex-sharing CaO 6 , CaO 7 and AsO 4 polyhedra, and we hypothesise that the new mineral has formed due to a topotactic reaction brought on by dehydration of preexisting guérinite.

Mineralogy↗
DOE OSTI2024

Optimal Geometry for Focused Ion Beam-Milled Samples for Direct-Pull Micro-Tensile Testing Performed In Situ in a Scanning Electron Microscope

A thorough procedure was developed to efficiently manufacture dogbone samples using focused ion beam (FIB) milling for micro-tensile testing. A Bruker PI 89 PicoIndenter, Billerica, MA, USA, was used as a case study, although the analysis and results are applicable to other micro-mechanical testing systems capable of mounting a standard, Ø12.7 mm × Ø3.2 mm pin, scanning electron microscopy (SEM) pin stub (Ted Pella, Redding, CA, USA). Nine dogbones were made from an Fe-45Cu alloy additively manufactured using powder-fed laser-directed energy deposition (DED-LB). Testing showed that fracture was confined to the gauge section for all dogbones and that the fracture mode, ductile vs. brittle, was entirely dependent on the grain orientation relative to the loading direction. The analysis showed that the measured plastic strain to failure can vary from >11% (optimal geometry) to <1% (non-optimal geometry) in micro-tensile testing of high-tensile-strength (>1 GPa) metallic materials. Subsequently, a finite element analysis (FEA) was conducted to identify the improved dogbone geometries. A total of ten thousand dogbone geometries were tested, and their dimensions were defined by a set of four adjustable parameters (corner radius, load surface angle, load surface length, and dogbone head length). The gauge width and gauge length were fixed to 4 µm and 10 µm, respectively. Three-dimensional surface plots of the stress concentration as a function of two parameters were used to identify the optimal ranges of parameter values. The addition of maximum width and length constraints, measuring 25 µm and 30 µm, respectively, allowed us to identify an optimal geometry at load surface angles of 30° and 45°. Their respective dimensions (corner radius, load surface length, and dogbone head length) are, in µm, 12, 6, and 7 and 10, 7, and 7. Testing these two optimal geometries with a range of gauge lengths from 4 to 20 µm showed that smaller gauge lengths only slightly reduced the detrimental stress concentration outside the gauge section. However, smaller gauge lengths will notably improve the FIB surface polishing step as tapering is reduced with smaller dogbone lengths.

Chemistry↗