Search NASA⌕ Search

SEARCH · Search NASA

Results for “Brittle failure”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Brittle failure analysis and modeling of high-burnup PWR fuel cladding alloys

The aim of this research is the development of methods for predicting mechanical behavior and identification of limiting conditions to prevent brittle failure of high-burnup (HBU) pressure water reactor (PWR) fuel cladding alloys. A finite element (FE) model of the ring compression test (RCT) was created to analyze the failure behavior of zirconium-based alloys with radial hydrides during the RCT. An elastic-plastic material model describes the zirconium alloy. The stress-strain curve needed for the elastic-plastic material model was derived by inverse finite element analyses. Cohesive zone modeling is used to reproduce sudden load drops during RCT loading. Based on the failure mechanism in non-irradiated ZIRLO (R) claddings, a micro-mechanical model was developed that distinguishes between brittle failure along hydrides and ductile failure of the zirconium matrix. Two different cohesive laws representing these types of failure are present in the same cohesive interface. The key differences between these constitutive laws are the cohesive strength, the stress at which damage initiates, and the cohesive energy, which is the damage energy dissipated by the cohesive zone. Statistically generated matrix-hydride distributions were mapped onto the cohesive elements and simulations with focus on the first load drop were performed. Computational results are in good agreement with the RCT results conducted on high-burnup M5 (R) samples. It could be shown that crack initiation and propagation strongly depend on the specific configuration of hydrides and matrix material in the fracture area.

Simbruner, Kai↗
DOE OSTI2025

Modelling the brittle failure of graphite induced by the controlled impact of runaway electrons in DIII-D

The thermo-mechanical response of an ATJ graphite sample to controlled runaway electron (RE) dissipation, realized in DIII-D, is modelled with a novel work-flow that features the RE orbit code KORC, the Monte Carlo particle transport code Geant4 and the finite element multiphysics software COMSOL. KORC provides the RE striking positions and momenta, Geant4 calculates the volumetric energy deposition and COMSOL simulates the thermoelastic response. Brittle failure is predicted according to the maximum normal stress criterion, which is suitable for ATJ graphite owing to its linear elastic behavior up to fracture and its isotropic mechanical properties. Measurements of the conducted energy, damage topology, explosion timing and blown-off material volume, impose a number of empirical constraints that suffice to distinguish between different RE impact scenarios and to identify RE parameters which provide the best match to the observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Large deformation and brittle failure calculated using the dual-domain material point method

The dual domain material point (DDMP) method is explored as a candidate to be implemented in a general purpose code to perform simulations of materials with complex geometry that undergo large history-dependent deformation and failure. To test its candidacy, we study its mesh convergence, its sensitivity to mesh orientation, and its ability to handle softening and failure of a material. Simulations of large deformation and simulations of mechanical failure are performed using both DDMP and the material point method (MPM). When cell-crossing of material points is not an issue and when there are a sufficient number of material points in each computation cell, the numerical error decreases with the square of the cell size as expected for both MPM and DDMP. DDMP has reduced error compared with MPM when there are many instances of material points crossing cell boundaries due to the continuous nature of the modified gradient of the shape functions. Simulations of a specimen under tension are also performed where the background mesh is aligned and misaligned with the tension direction. MPM displays a significant mesh-dependent stress field, DDMP shows negligible mesh dependency. Despite a mesh orientation-dependent stress field from MPM, the critical tension and failure mode from both MPM and DDMP calculations have negligible mesh dependency when using a non-local failure model. If only the failure mode is important (i.e., local stresses are unimportant), MPM with a non-local failure model is a suitable method for modeling failure with small deformations. However, if local stresses are also important or if there are large deformations with many cell-crossings before failure, DDMP should be the method that is used. A needed improvement for DDMP is identified from our numerical simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A unified finite strain gradient-enhanced micropolar continuum approach for modeling quasi-brittle failure of cohesive-frictional materials

Here, in this work, a novel framework for modeling quasi-brittle crack propagation and shear band dominated failure of cohesive-frictional materials like concrete, mortar, rock, tough ceramics, energetic materials, but also granular materials like sands or powders in terms of a unified continuum approach is proposed. It is based on a combination of the gradient-enhanced continuum with gradients of internal variables for representing quasi-brittle cracking, and the micropolar continuum, accounting for the deformation of the microstructure. For developing the gradient-enhanced micropolar framework, the set of balance equations and the kinematic relations are derived, and the constitutive relations are established in a general manner. The framework is formulated in a geometrically exact setting, based on the thermodynamically sound theory of hyperelasto-plasticity, and the numerical implementation by means of the finite element method is discussed. For assessing the approach, realizations of this new approach in terms of constitutive models for particular materials are developed. They are applied to numerical benchmark examples, investigating various loading conditions, and the obtained results are validated by means of a comparison with experiments from the literature.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Probability of brittle failure

A methodology was developed for collecting statistically representative data for crack initiation and arrest from small number of test specimens. An epoxy (based on bisphenol A diglycidyl ether and polyglycol extended diglycyl ether and cured with diethylene triamine) is selected as a model material. A compact tension specimen with displacement controlled loading is used to observe multiple crack initiation and arrests. The energy release rate at crack initiation is significantly higher than that at a crack arrest, as has been observed elsewhere. The difference between these energy release rates is found to depend on specimen size (scale effect), and is quantitatively related to the fracture surface morphology. The scale effect, similar to that in statistical strength theory, is usually attributed to the statistics of defects which control the fracture process. Triangular shaped ripples (deltoids) are formed on the fracture surface during the slow subcritical crack growth, prior to the smooth mirror-like surface characteristic of fast cracks. The deltoids are complementary on the two crack faces which excludes any inelastic deformation from consideration. Presence of defects is also suggested by the observed scale effect. However, there are no defects at the deltoid apexes detectable down to the 0.1 micron level.

Kim, A.↗
DOE OSTI2023

A robust interface finite element formulation for modeling brittle material failure problems

Abstract Failure of many brittle materials and structures can be modeled using interface‐oriented finite elements combined with intrinsic cohesive zone models. The discontinuous Galerkin (DG) finite element method provides an innovative framework for modeling brittle crack propagation with zero‐thickness interface elements, which can accommodate extrinsic cohesive laws to avoid the artificial compliance required in intrinsic cohesive models. However, robust formulations and implementations of DG methods are critical in alleviating the well‐known convergence issues for both crack nucleation and propagation with reduced instability. This paper presents a robust interface element formulation by modifying the incomplete interior penalty Galerkin (IIPG) method, which successfully avoids the initial element interface penetration across elements that occurs prior to crack nucleation, and thereby greatly reduces the instability issue as cracks open. We further verified and validated our implementation by using a bar tension test and a beam fracturing benchmark. The robustness of our proposed interface element method was demonstrated by a micromechanics fiber/matrix debonding problem with 64 fibers embedded in a bulk matrix.

Liu, Ruijie↗
NASA NTRS1985

Tectonic determinations of lithospheric thickness on Ganymede and Callisto

The concept of the Maxwell time of a viscoelastic material (4.5) is used in conjunction with calculated thermal profiles to evaluate the significance of tectonic estimates of lithospheric thickness. Thermal lithospheric thicknesses provide fundamental constraints on planetary thermal histories that complement the constraints provided by dateable surface deposits of endogenic origin. Lithospheric constraints are of particular value on the icy satellites where our understanding of both rheology and surface ages is considerably poorer than it is for the terrestrial planets. Certain extensional tectonic features can and have been used to estimate lithospheric thicknesses on Ganymede and Callisto. These estimates, however, refer to the depth of the elastic lithosphere defined by the zone of brittle failure. The relation between the elastic lithosphere and the thermal lithosphere (generally defined by the zone of conductive heat transport) is not straightforward, because the depth of brittle failure depends not only on the thermal profile, but also on rheology and strain rate (or the characteristic time over which stresses build towards failure). Characteristic time considerations are not trivial in this context because stresses generating brittle failure on the icy satellites may be produced by impacts, with characteristic times of seconds to days, or by geologic processes with time scales of hundreds of millions of years.

Croft, S. K.↗
NASA NTRS1995

Analysis of the progressive failure of brittle matrix composites

This report investigates two of the most common modes of localized failures, namely, periodic fiber-bridged matrix cracks and transverse matrix cracks. A modification of Daniels' bundle theory is combined with Weibull's weakest link theory to model the statistical distribution of the periodic matrix cracking strength for an individual layer. Results of the model predictions are compared with experimental data from the open literature. Extensions to the model are made to account for possible imperfections within the layer (i.e., nonuniform fiber lengths, irregular crack spacing, and degraded in-situ fiber properties), and the results of these studies are presented. A generalized shear-lag analysis is derived which is capable of modeling the development of transverse matrix cracks in material systems having a general multilayer configuration and under states of full in-plane load. A method for computing the effective elastic properties for the damaged layer at the global level is detailed based upon the solution for the effects of the damage at the local level. This methodology is general in nature and is therefore also applicable to (0(sub m)/90(sub n))(sub s) systems. The characteristic stress-strain response for more general cases is shown to be qualitatively correct (experimental data is not available for a quantitative evaluation), and the damage evolution is recorded in terms of the matrix crack density as a function of the applied strain. Probabilistic effects are introduced to account for the statistical nature of the material strengths, thus allowing cumulative distribution curves for the probability of failure to be generated for each of the example laminates. Additionally, Oh and Finney's classic work on fracture location in brittle materials is extended and combined with the shear-lag analysis. The result is an analytical form for predicting the probability density function for the location of the next transverse crack occurrence within a crack bounded region. The results of this study verified qualitatively the validity of assuming a uniform crack spacing (as was done in the shear-lag model).

Thomas, David J.↗
DOE OSTI2022

Dynamic properties of FeCrMnNi, a high entropy alloy

The goal of this paper was to assess the dynamic properties of a new class of materials, High entropy alloys (HEAs). Specifically, gas gun experiments coupled with recovery were performed on a specific HEA, FeCrMnNi, to measure its equation of state (EOS) and its spall strength, which is a measure of the stress required to nucleate voids under shock loading conditions. While there has been a plethora of work performed to assess the mechanical properties of these HEAs under uniaxial stress conditions as a function of strain rate, investigation of its properties in dynamic extremes remains rare. The current work fills this gap in knowledge for this novel class of materials. Our results show that the measured EOS for this material at one velocity was in reasonable agreement with an estimated Hugoniot. Furthermore, while the spall strength for this material was found to be ~1.9 GPa, with some variation based on sample location, the main failure mode was brittle. This brittle failure leading to formation of cracks in the material is different than the ductile failure observed in pure iron and its alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A novel post-processing method for progressive failure analysis of brittle composite compression

Finite element analysis of brittle materials in axial compression typically uses element deletion to allow continued global deformation post-element-failure. However, element deletion produces cyclic load-displacement curves that underestimate energy absorption and are not representative of a continuum system. Two key observations support the conclusion that results from an appropriately discretized model can be an adequate representation of a continuum system. Specifically, the frequency of the oscillations in the load-displacement curve is directly dependent upon element length in the loading direction, and the peak amplitudes of oscillations are mesh size independent. A method of post-processing the analysis results, by connecting the peak amplitudes of oscillations, is proposed and applied to a series of continuous carbon fiber composite crush tubes. The load-displacement curve, stable crushing load, and specific energy absorption of the post-processed results compare well to an experimental study of crush tubes with similar layups.

Materials Science↗
NASA NTRS1983

Experimental deformation of polycrystalline H2O ice at high pressure and low temperature - Preliminary results

A preliminary study is carried out of involving 70 constant strain deformation tests on pure polycrystalline H2O ice under conditions covering most of the stability field of ice I sub h. Brittle failure of I sub h is found to be promoted by lower P, lower T, and higher strain rates. Ductile flow is found to be promoted by higher P, higher T, and lower strain rates. The brittle failure of ice I sub h is found to be most unusual. The fracture strength is a positive function of P only below 50 MPa. At pressures greater than this, the fracture strength is independent of P, and the fracture plane lies approximately 45 deg from the load axis. It is believed that existing extrapolation based on existing experimental data to Ganymede and Callisto may be badly in error.

Durham, W. B.↗
NASA NTRS2009

Multiscale Computer Simulation of Tensile and Compressive Strain in Polymer- Coated Silica Aerogels

While the low thermal conductivities of silica aerogels have made them of interest to the aerospace community as lightweight thermal insulation, the application of conformal polymer coatings to these gels increases their strength significantly, making them potentially useful as structural materials as well. In this work we perform multiscale computer simulations to investigate the tensile and compressive strain behavior of silica and polymer-coated silica aerogels. Aerogels are made up of clusters of interconnected particles of amorphous silica of less than bulk density. We simulate gel nanostructure using a Diffusion Limited Cluster Aggregation (DLCA) procedure, which produces aggregates that exhibit fractal dimensions similar to those observed in real aerogels. We have previously found that model gels obtained via DLCA exhibited stress-strain curves characteristic of the experimentally observed brittle failure. However, the strain energetics near the expected point of failure were not consistent with such failure. This shortcoming may be due to the fact that the DLCA process produces model gels that are lacking in closed-loop substructures, compared with real gels. Our model gels therefore contain an excess of dangling strands, which tend to unravel under tensile strain, producing non-brittle failure. To address this problem, we have incorporated a modification to the DLCA algorithm that specifically produces closed loops in the model gels. We obtain the strain energetics of interparticle connections via atomistic molecular statics, and abstract the collective energy of the atomic bonds into a Morse potential scaled to describe gel particle interactions. Polymer coatings are similarly described. We apply repeated small uniaxial strains to DLCA clusters, and allow relaxation of the center eighty percent of the cluster between strains. The simulations produce energetics and stress-strain curves for looped and nonlooped clusters, for a variety of densities and interaction parameters.

Good, Brian↗
DOE OSTI2022

Dynamic properties of FeCrMnNi, a High Entropy Alloy

The goal of this paper was to assess the dynamic properties of a new class of materials, High entropy alloys (HEAs). Specifically, gas gun experiments coupled with recovery were performed on a specific HEA, FeCrMnNi, to measure its equation of state (EOS) and its spall strength, which is a measure of the stress required to nucleate voids under shock loading conditions. While there has been a plethora of work performed to assess the mechanical properties of these HEAs under uniaxial stress conditions as a function of strain rate, investigation of its properties in dynamic extremes remains rare. The current work fills this gap in knowledge for this novel class of materials. Our results show that the measured EOS for this material at one velocity was in reasonable agreement with an estimated Hugoniot. Furthermore, while the spall strength for this material was found to be ~1.9 GPa, with some variation based on sample location, the main failure mode wasbrittle. This brittle failure leading to formation of cracks in the material is different than the ductile failure observed in pure iron and its alloys. DOE/NV/03624--1298

36 MATERIALS SCIENCE↗
DOE OSTI2020

High-burnup Experiments in Reactivity Initiated Accidents (HERA)

HERA is a joint experimental program (JEEP) operating within the Nuclear Energy Agency’s (NEA’s) framework for irradiation experiments (FIDES). HERA is dedicated to the understanding of light water reactor (LWR) fuel performance at high burnup under reactivity-initiated accidents (RIA). During RIAs the fuel cladding can be breached by one of two mechanisms; pellet-cladding mechanical interaction or high temperature balloon and rupture. Pellet-cladding mechanical interaction (PCMI) failure occurs in the early phase of an RIA transient where the fuel pellet expands rapidly into the cladding prior to any meaningful heat transfer from the pellet to the cladding. This results in high stresses in the cladding resulting in crack propagation through the base Zircaloy metal and axially along the cladding tube. Failure by PCMI is a brittle failure mode and occurs at low cladding strains, approximately 2% hoop strain. If the cladding survives the PCMI phase of the transient, it may still fail due to ballooning and rupture induced by extended exposure to high temperature. Thermal transport from the fuel pin to the coolant is significantly reduced if a transition to film boiling on the cladding surface occurs. This result in both a rapid increase in fuel temperature, which drives fission gas release, and a spike in cladding temperature, which results in loss of mechanical strength. If the internal pressure of the fuel rod is greater than the system pressure the cladding can inelastically deform and eventually rupture. Failure by swelling and rupture occurs at larger strains, generally greater than 5% hoop strain. While PCMI failures are generally seen as more limiting, the dominance of one failure mode over the other can be affected by both the cladding material conditions, and the transient evolution.

(HERA)↗
DOE OSTI2023

High-burnup Experiments in Reactivity Initiated Accidents (HERA)

HERA is a joint experimental program (JEEP) operating within the Nuclear Energy Agency’s (NEA’s) framework for irradiation experiments (FIDES). HERA is dedicated to the understanding of light water reactor (LWR) fuel performance at high burnup under reactivity-initiated accidents (RIA). During RIAs the fuel cladding can be breached by one of two mechanisms; pellet-cladding mechanical interaction or high temperature balloon and rupture. Pellet-cladding mechanical interaction (PCMI) failure occurs in the early phase of an RIA transient where the fuel pellet expands rapidly into the cladding prior to any meaningful heat transfer from the pellet to the cladding. This results in high stresses in the cladding resulting in crack propagation through the base Zircaloy metal and axially along the cladding tube. Failure by PCMI is a brittle failure mode and occurs at low cladding strains, approximately 2% hoop strain. If the cladding survives the PCMI phase of the transient, it may still fail due to ballooning and rupture induced by extended exposure to high temperature. Thermal transport from the fuel pin to the coolant is significantly reduced if a transition to film boiling on the cladding surface occurs. This result in both a rapid increase in fuel temperature, which drives fission gas release, and a spike in cladding temperature, which results in loss of mechanical strength. If the internal pressure of the fuel rod is greater than the system pressure the cladding can inelastically deform and eventually rupture. Failure by swelling and rupture occurs at larger strains, generally greater than 5% hoop strain. While PCMI failures are generally seen as more limiting, the dominance of one failure mode over the other can be affected by both the cladding material conditions, and the transient evolution.

(HERA)↗
DOE OSTI2021

Deformation and fracture behavior of a laser powder bed fusion processed stainless steel: In situ synchrotron x-ray computed microtomography study

Tensile plastic deformation and fracture behavior of a laser powder bed fusion processed 316L stainless steel alloy was investigated in situ using high-resolution synchrotron x-ray computed microtomography (sXCT). Three different cases were studied: (1) specimen A1, near full-density with initial porosity of 0.04%, (2) specimen B1, initial lack-of-fusion (LOF) porosity of about 3% with the major axis of defects aligned nearly perpendicular to the tensile loading direction (LD), and (3) specimen B, similar to B1 but with the major axis of LOF defects aligned nearly parallel to LD. The effects of type, density, and orientation of defects in as-printed samples on plastic deformation and failure modes were studied by analyzing sXCT image sequences and quantitative defect properties. Specifically, the evolutions of porosity, size, aspect ratio, and angle of gas pores and LOF defects were characterized as a function of applied strain. Moreover, interactions between as-printed defects and deformation-induced crack propagation were investigated. The tensile behavior of specimen A1 was comparable to that of a commercial wrought plate. Furthermore, the quantitative in-situ sXCT results showed a void-growth phenomenon in A1 but without noticeable coalescence, resulting in a ductile cup-cone failure. Specimen B1 showed damage-driven brittle failure with limited plasticity due to early and rapid propagation of open voids. On the other hand, specimen B, with a similar amount of LOF defects as B1, showed deceivingly strong and ductile tensile behavior similar to A1 without any major crack propagations due to the favorable orientation relationship between LOF defects and the LD.

36 MATERIALS SCIENCE↗