SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

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

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

NASA NTRS1993

A bilinear failure criterion for mixed-mode delamination

Many different failure criteria have been suggested for mixed-mode delamination toughness, but few sets of mixed-mode data exist that are consistent over the full range of Mode 1 opening load to Mode 2 shear load range. The mixed-mode bending (MMB) test was used to measure the delamination toughness of a brittle epoxy composite, a state-of-the-art toughened epoxy composite, and a tough thermoplastic composite over the full mixed-mode range. To gain insight into the different failure responses of the different materials, the delamination fracture surfaces were also examined. An evaluation of several failure criteria that have been reported in the literature was performed, and the range of responses modeled by each criterion was analyzed. A bilinear failure criterion was introduced based on a change in the failure mechanism observed from the delamination surfaces. The different criteria were compared to the failure response of the three materials tested. The responses of the two epoxies were best modeled with the new bilinear failure criterion. The failure response of the tough thermoplastic composite could be modeled well with the bilinear criterion but could also be modeled with the more simple linear failure criterion. Since the materials differed in their mixed-mode failure response, mixed-mode delamination testing will be needed to characterize a composite material. This paper presents consistent sets of mixed-mode data, provides a critical evaluation of the mixed-mode failure criteria, and should provide general guidance for selecting an appropriate criterion for other materials.

Reeder, James R.↗

NASA NTRS1998

Global Failure Modes in High Temperature Composite Structures

Composite materials have been considered for many years as the major advance in the construction of energy efficient aerospace structures. Notable advances have been made in understanding the special design considerations that set composites apart from the usual "isotropic" engineering materials such as the metals. As a result, a number of significant engineering designs have been accomplished. However, one shortcoming of the currently favored composites is their relatively unforgiving behavior with respect to failure (brittleness) under seemingly mild impact conditions and large efforts are underway to rectify that situation, much along the lines of introducing thermoplastic matrix materials. Because of their relatively more pronounced (thermo) viscoelastic behavior these materials respond with "toughness" in fracture situations. From the point of view of applications requiring material strength, this property is highly desirable. This feature impacts several important and distinct engineering problems which have been' considered under this grant and cover the 1) effect of impact damage on structural (buckling) stability of composite panels, the 2) effect of time dependence on the progression of buckling instabilities, and the 3) evolution of damage and fracture at generic thickness discontinuities in structures. The latter topic has serious implications for structural stability problems (buckling failure in reinforced shell structures) as well as failure progression in stringer-reinforced shell structures. This grant has dealt with these issues. Polymer "toughness" is usually associated with uncrosslinked or thermo-plastic polymers. But, by comparison with their thermoset counterparts they tend to exhibit more pronounced time dependent material behavior; also, that time dependence can occur at lower temperatures which places restriction in the high temperature use of these "newer and tougher" materials that are not quite so serious with the thermoset matrix materials. From a structural point of view the implications of this material behavior are potentially severe in that structural failure characteristics are no longer readily observed in short term qualification tests so characteristic for aerospace structures built from typical engineering metals.

Knauss, W. G.↗

NASA NTRS2011

Statistical Models of Fracture Relevant to Nuclear-Grade Graphite: Review and Recommendations

The nuclear-grade (low-impurity) graphite needed for the fuel element and moderator material for next-generation (Gen IV) reactors displays large scatter in strength and a nonlinear stress-strain response from damage accumulation. This response can be characterized as quasi-brittle. In this expanded review, relevant statistical failure models for various brittle and quasi-brittle material systems are discussed with regard to strength distribution, size effect, multiaxial strength, and damage accumulation. This includes descriptions of the Weibull, Batdorf, and Burchell models as well as models that describe the strength response of composite materials, which involves distributed damage. Results from lattice simulations are included for a physics-based description of material breakdown. Consideration is given to the predicted transition between brittle and quasi-brittle damage behavior versus the density of damage (level of disorder) within the material system. The literature indicates that weakest-link-based failure modeling approaches appear to be reasonably robust in that they can be applied to materials that display distributed damage, provided that the level of disorder in the material is not too large. The Weibull distribution is argued to be the most appropriate statistical distribution to model the stochastic-strength response of graphite.

Nemeth, Noel N.↗

NASA NTRS1989

The nature of cometary materials

Because cometary surfaces are likely to be far colder and of a different composition than planetary surfaces, there are some new considerations that must be examined in regards to placing instrumented packages or sample return devices on their surfaces. The qualitative analysis of the problem of attaching hardware to a comet and not being ejected back into space can be divided into two parts. The first problem is to pierce the mantle and obtain access to the icy core. Drilling through the mantle requires that the drilling forces be reacted. Reacting such forces probably requires attachment to the icy core below. Therefore, some kinetic impact piercing device is likely to be required as the first act of attachment. The second problem for a piercing device to overcome is the force produced by the impact kinetic energy that tries to eject the piercing device back into space. The mantle and icy core can absorb some of the impact kinetic energy in the form of fracture formation and friction energy. The energy that is not absorbed in these two ways is stored by the core as elastic deformation of the mantle and icy core. It is concluded that because the cometary materials are almost certainly brittle and the icy core is likely to be self lubricating, the elastic rebound and gas pressure expulsion forces must be counteracted by forces greater than those that may be provided by a piercing device or its capture devices (barbs).

Stephens, James↗

DOE OSTI2026

Influence of preheating on tungsten damage under pulsed hydrogen plasma exposure

Preheated tungsten plate behavior under multiple pulsed hydrogen plasma irradiations in the PW-7 plasma accelerator was investigated. The preheating temperature was 973 K, which is above the ductile to brittle transition temperature (DBTT ≈ 673 K). The plate was exposed to three consecutive series of 13 plasma pulses each. Preheating sufficiently suppressed surface cracking at the early irradiation stage: after 13 pulses, cracks were rare, localized, and did not form extensive fractures. During subsequent irradiation series, the area fraction of fine cracks increased by ∼ 2% due to accumulation of dislocations around the pre-existing defects. However, their width remained nearly constant, averaging 0.2–0.3 μm, with rare cracks up to 1 μm in size. Repeated plasma exposure resulted in pronounced surface roughening associated with multiple melting cycles, crater formation and blister rupture, with roughness reaching ∼ 204 nm and a recorded sample mass loss of ∼ 2.8 mg. These results indicate that controlling tungsten preheating well above DBTT effectively reduces damage. However, the surface condition during initial irradiation significantly influences its subsequent morphological evolution. Preheating to 773 K, selected regarding the minimum DBTT value, has previously been demonstrated to be an effective approach for limiting crack growth. A further increase in temperature of approximately ∼ 1000 K did not qualitatively improve results: crack still formed and widths did not decrease. At the same time, the absence of further crack coarsening indicates that the limiting effect of preheating has been reached. Considering the increased energy costs, further increasing the preheating temperature seems unjustified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Integrated design of aluminum-enriched high-entropy refractory B2 alloys with synergy of high strength and ductility

Refractory high-entropy alloys (RHEAs) are promising high-temperature structural materials. Their large compositional space poses great design challenges for phase control and high strength-ductility synergy. The present research pioneers using integrated high-throughput machine learning with Monte Carlo simulations supplemented by ab initio calculations to effectively navigate phase selection and mechanical property predictions, developing single-phase ordered B2 aluminum-enriched RHEAs (Al-RHEAs) demonstrating high strength and ductility. These Al-RHEAs achieve remarkable mechanical properties, including compressive yield strengths up to 1.7 gigapascals, fracture strains exceeding 50%, and notable high-temperature strength retention. They also demonstrate a tensile yield strength of 1.0 gigapascals with a ductility of 9%, albeit with B2 ordering. Furthermore, we identify valence electron count domains for alloy ductility and brittleness with the explanation from density functional theory and provide crucial insights into elemental influence on atomic ordering and mechanical performance. The work sets forth a strategic blueprint for high-throughput alloy design and reveals fundamental principles governing the mechanical properties of advanced structural alloys.

Science & Technology - Other Topics↗

NASA NTRS1994

A new method for failure prediction of SR-200 beryllium sheet structures

Contemporary applications of failure criteria frequently incorporate two-dimensional or simplified three-dimensional methodologies for prediction of stresses. Motivation behind the development of a new multi-dimensional failure criterion is due mainly to the lack of a sufficiently accurate mathematical tool that accounts for the behavior of brittle material with anisotropic properties. Such a criterion should be able to provide a reliable maximum load estimate so that design of the structure is not penalized in terms of excessive weight requirements. The failure criterion developed is represented by a fracture surface in a six-dimensional stress space. The criterion is applied for failure prediction of SR-200 beryllium sheet structures, a non-homogeneous orthotropic material used widely in space applications. Two experiments are used to verify the criterion.

Papados, P. P.↗

NASA NTRS1994

Processing and Mechanical Properties of NiAl-Based In-Situ Composites

In-situ composites based on the NiAl-Cr eutectic system were successfully produced by containerless processing and evaluated. The NiAl-Cr alloys had a fibrous microstructure while the NiAl-(Cr,Mo) alloys containing 1 at. percent or more molybdenum exhibited a lamellar structure. The NiAl-28Cr-6Mo eutectic displays promising high temperature strength while still maintaining a reasonable room temperature fracture toughness when compared to other NiAl-based materials. The Laves phase NiAlTa was used to strengthen NiAl and very promising creep strengths were found for the directionally solidified NiAl-NiAlTa eutectic. The eutectic composition was found to be near NiAl-15.5Ta (at. percent) and well aligned microstructures were produced at this composition. An off-eutectic composition of NiAl-14.5Ta was also processed, consisting of NiAl dendrites surrounded by aligned eutectic regions. The room temperature toughness of these two phase alloys was similar to that of polycrystalline NiAl even with the presence of the brittle Laves phase NiAlTa. Polyphase in-situ composites were generated by directional solidification of ternary eutectics. The systems investigated were the Ni-Al-Ta-X (X=Cr, Mo, or V) alloys. Ternary eutectics were found in each of these systems and both the eutectic composition and temperature were determined. Of these ternary eutectics, the one in the NiAl-Ta-Cr system was found to be the most promising. The fracture toughness of the NiAl-(Cr,Al)NiTa-Cr eutectic was intermediate between those of the NiAl-NiAlTa eutectic and the NiAl-Cr eutectic. The creep strength of this ternary eutectic was similar to or greater than that of the NiAl-Cr eutectic.

Johnson, David Ray↗

DOE OSTI2020

Multiscale Thermomechanical Damage Model with Internal State Variables for Ceramic Matrix Composites

A length scale-dependent matrix damage model has been developed to simulate the response of ceramic matrix composites (CMCs). This model captures damage initiation and propagation in the brittle matrix material by employing internal state variable (ISV) theory within a multiscale modeling framework to obtain damaged matrix stress-strain constitutive relationships at each length scale. Material degradation due to matrix cracking is captured by a damage variable determined using fracture mechanics and the self-consistent scheme. An additional state variable is introduced to capture the nucleation and growth of matrix porosity, a key mechanism contributing to matrix nonlinearity. Porosity occurs as a result of material diffusion around grain boundaries and the growth of the porosity state variable is related to the material entropy dissipation and the volumetric strain. The nonlinear predictive capabilities of the material model are demonstrated for monolithic silicon carbide, unidirectional (UD) carbon fiber/silicon carbide matrix (C/SiC) CMC, and five-harness satin (5HS) woven C/SiC CMC.

Skinner, Travis D.↗

DOE OSTI2023

From Well Log to Formation Model: A Novel Laboratory Calibrated Methodology with Demonstration

This work demonstrates how the characterization and modeling of both elastic and creep properties are essential to describe zones in layered rock formations as either low-stress targets for stimulation or high-stress barriers to fracture growth. Prediction of fracture height is critical for designing stimulation operations in oil and gas wells. Ideally, fractures are placed in target zones which will produce hydrocarbons and should not propagate into zones expected to be unproductive or to produce unwanted fluids such as water which in turn must be treated and/or disposed. The essential task in designing stimulation plans is predicting which zones have low horizontal stresses and which will be high-stress barriers to fracture growth. Despite this importance, there are gaps in current knowledge and a complete workflow from laboratory characterization to a finite element model which includes time dependent rock deformation is required. While the research and methodology presented here also have application to CO 2or hydrogen storage, wastewater injection, and geothermal applications, the focus will be on hydrocarbon extraction. This thesis presents the results of a characterization-to-prediction workflow for the Caney shale, which is an emerging hydrocarbon resource in Oklahoma, USA. It begins with an investigation to enable critical evaluation of the Caney zonation into nominally “brittle” and “ductile” zones based on properties observed from well logs. It shows none of the zones are consistently “brittle” or “ductile” mechanical behavior based on the variety of definitions of these terms. However, the nominally ductile zones are weaker and more prone to creep. A laboratory investigation of samples including strength, elastic, and creep properties, is then used in a finite element model of stress evolution. The model includes both elastic deformation and viscoplastic creep. Results predict the least creep-prone layers to have the lowest horizontal stresses, therefore comprising hydraulic fracturing targets. The most creep-prone layers attain a horizontal stress similar to the vertical stress and therefore are predicted to be high stress barriers to hydraulic fracture stimulation. In addition to defining stimulation target intervals, the model shows how as tectonic strain rate increases, there is a transition from creep-dominated stresses to stresses dominated by elasticity.

Benge, Margaret↗

DOE OSTI2024

PDMATLAB2D: A Peridynamics MATLAB Two-dimensional Code

PDMATLAB2D is a meshfree peridynamics implementation in MATLAB suitable for simulation of two-dimensional fracture problems. The purpose of this code is twofold. First, it provides an entry-level peridynamics computational tool for educational and training purposes. Second, it serves as an accessible and easily modifiable computational tool for peridynamics researchers who would like to adapt the code for a multitude of peridynamics simulation scenarios. The current version of the code implements a bond-based brittle elastic peridynamic model and a critical stretch criterion for bond breaking. However, the code is designed to be extendable for other peridynamic models and computational features. In this paper, we provide an overview of the code structure and functions with illustrative examples. In conclusion, due to the integrated computation and postprocessing MATLAB capabilities, PDMATLAB2D can serve as an effective testbed for testing new constitutive models and advanced numerical features for peridynamics computations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Impact of supercritical carbon dioxide on the frictional strength of and the transport through thin cracks in shale

Understanding the mechanical and transport behavior of thin (i.e. small aperture) cracks slipping under supercritical carbon dioxide (sc-CO 2 ) conditions is essential to evaluate the integrity of sealing formations with buoyant sc-CO 2 below and the success of waterless fracturing. The two major items of interest in this work are frictional strength and permeability change of the crack. We used a triaxial cell that permits in situ visualization to conduct and monitor slippage along the faces of narrow cracks subjected to triaxial stresses. Such cracks are analogs to small geological faults. We tested carbonate-rich, 1-inch diameter Wolfcamp shale samples that are saw cut 30° to vertical to create a thin crack. Friction coefficients ranged from about 0.6 to 0.8 consistent with expectations for brittle rocks. The sc-CO 2 generally did not alter friction coefficient over the time scale of experiments. From a transport perspective, saturating cracks with sc-CO 2 substantially decreased permeability of the crack by 26%–52%, while slip resulted in a variety of permeability responses. Overall, the combined impact of sc-CO 2 saturation and slip reduced fault permeability for all tests. Our observations support the notion that the sealing capacity of some caprocks improves when saturated with sc-CO 2 and that some slip of small fractures is not necessarily detrimental to caprock integrity.

58 GEOSCIENCES↗

NASA NTRS2002

Reliability Analysis of Brittle Material Structures - Including MEMS(?) - With the CARES/Life Program

Brittle materials are being used, or considered, for a wide variety of high tech applications that operate in harsh environments, including static and rotating turbine parts. thermal protection systems, dental prosthetics, fuel cells, oxygen transport membranes, radomes, and MEMS. Designing components to sustain repeated load without fracturing while using the minimum amount of material requires the use of a probabilistic design methodology. The CARES/Life code provides a general-purpose analysis tool that predicts the probability of failure of a ceramic component as a function of its time in service. For this presentation an interview of the CARES/Life program will be provided. Emphasis will be placed on describing the latest enhancements to the code for reliability analysis with time varying loads and temperatures (fully transient reliability analysis). Also, early efforts in investigating the validity of using Weibull statistics, the basis of the CARES/Life program, to characterize the strength of MEMS structures will be described as as well as the version of CARES/Life for MEMS (CARES/MEMS) being prepared which incorporates single crystal and edge flaw reliability analysis capability. It is hoped this talk will open a dialog for potential collaboration in the area of MEMS testing and life prediction.

Nemeth, Noel N.↗

DOE OSTI2021

Recovering the effects of subgrid heterogeneity in simulations of radionuclide transport through fractured media

Groundwater flow and contaminant transport through fractured media can be simulated using Discrete Fracture Network (DFN) models which provide a natural description of structural heterogeneity. However, this approach is computationally expensive, with the large number of intersecting fractures necessitated by many real-world applications requiring modeling simplifications to be made for calculations to be tractable. Upscaling methods commonly used for this purpose can result in some loss of local-scale variability in the groundwater flow velocity field, resulting in underestimation of particle travel times, transport resistance and retention in transport calculations. In this paper, a transport downscaling algorithm to recover the transport effects of heterogeneity is tested on a synthetic Brittle Fault Zone model, motivated by the problem of large safety assessment calculations for geological repositories of spent nuclear fuel. We show that the variability in the local-scale velocity field which is lost by upscaling can be recovered by sampling from a library of DFN transport paths, accurately reproducing DFN transport statistic distributions and radionuclide breakthrough curves in an upscaled model.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1993

Tension fracture of laminates for transport fuselage. Part 2: Large notches

Tests were conducted on over 200 center-crack specimens to evaluate: (a) the tension-fracture performance of candidate materials and laminates for commercial fuselage applications; and (b) the accuracy of several failure criteria in predicting response. Crack lengths of up to 12 inches were considered. Other variables included fiber/matrix combination, layup, lamination manufacturing process, and intraply hybridization. Laminates fabricated using the automated tow-placement process provided significantly higher tension-fracture strengths than nominally identical tape laminates. This confirmed earlier findings for other layups, and possibly relates to a reduced stress concentration resulting from a larger scale of repeatable material inhomogeneity in the tow-placed laminates. Changes in material and layup result in a trade-off between small-notch and large-notch strengths. Toughened resins and 0 deg-dominate layups result in higher small-notch strengths but lower large-notch strengths than brittle resins, 90 deg and 45 deg dominated layups, and intraply S2-glass hybrid material forms. Test results indicate that strength-prediction methods that allow for a reduced order singularity of the crack-tip stress field are more successful at predicting failure over a range of notch sizes than those relying on the classical square-root singularity. The order of singularity required to accurately predict large-notch strength from small-notch data was affected by both material and layup. Measured crack-tip strain distributions were generally higher than those predicted using classical methods. Traditional methods of correcting for finite specimen width were found to be lacking, confirming earlier findings with other specimen geometries. Fracture tests of two stiffened panels, identical except for differing materials, with severed central stiffeners resulted in nearly identical damage progression and failure sequences. Strain-softening laws implemented within finite element models appear attractive to account for load redistribution in configured structure due to damage-induced crack tip softening

Walker, Tom H.↗

NASA NTRS2008

Impact Resistance of Environmental Barrier Coated SiC/SiC Composites

Impact performance of 2D woven SiC/SiC composites coated with 225 and 525 microns thick environmental barrier coating (EBC) was investigated. The composites were fabricated by melt infiltration and the EBC was deposited by plasma spray. Impact tests were conducted at room temperature and at 1316 C in air using 1.59 mm diameter steel-balls at projectile velocities ranging from 110 to 375 m/s . Both microscopy and non-destructive evaluation (NDE) methods were used to determine the extent of damage in the substrate and coating with increasing projectile velocity. The impacted specimens were tensile tested at room temperature to determine their residual mechanical properties. At projectile velocities less than 125 m/s , no detectable internal damage was noticed in the MI SiC/SiC composites coated with 525 microns EBC. With increase in projectile velocity beyond this value, spallation of EBC layers, delamination of fiber plies, and fiber fracture were detected. At a fixed projectile velocity, the composites coated with 525 microns EBC showed less damage than those coated with 225 microns EBC. Both types of coated composites retained a large fraction of the baseline properties of the as-fabricated composites and exhibited non-brittle failure after impact testing. Furnace exposure of impacted specimens in a moisture environment at 1316 C for 500 h indicated that the through-the-thickness cracks in the coating and delamination cracks in the substrate generated after impact testing acted as conduits for internal oxidation.

Bhatt, Ramakrishna T.↗

NASA NTRS2013

Unit-Sphere Multiaxial Stochastic-Strength Model Applied to Anisotropic and Composite Materials

Models that predict the failure probability of brittle materials under multiaxial loading have been developed by authors such as Batdorf, Evans, and Matsuo. These "unit-sphere" models assume that the strength-controlling flaws are randomly oriented, noninteracting planar microcracks of specified geometry but of variable size. This methodology has been extended to predict the multiaxial strength response of transversely isotropic brittle materials, including polymer matrix composites (PMCs), by considering (1) flaw-orientation anisotropy, whereby a preexisting microcrack has a higher likelihood of being oriented in one direction over another direction, and (2) critical strength, or K (sub Ic) orientation anisotropy, whereby the level of critical strength or fracture toughness for mode I crack propagation, K (sub Ic), changes with regard to the orientation of the microstructure. In this report, results from finite element analysis of a fiber-reinforced-matrix unit cell were used with the unit-sphere model to predict the biaxial strength response of a unidirectional PMC previously reported from the World-Wide Failure Exercise. Results for nuclear-grade graphite materials under biaxial loading are also shown for comparison. This effort was successful in predicting the multiaxial strength response for the chosen problems. Findings regarding stress-state interactions and failure modes also are provided.

Nemeth, Noel, N.↗

DOE OSTI2023

Implementation of extrinsic cohesive zone model (ECZM) in 2D finite-discrete element method (FDEM) using node binding scheme

The combined finite-discrete element method (FDEM) has been widely used for rock fracturing simulations. Conventionally, FDEM is realized using the intrinsic cohesive zone model (ICZM); however, it has the drawback of artificial compliance and high computational expense. As a complement, the extrinsic cohesive zone model (ECZM) is seen to be realized in FDEM recently, whereas the node splitting scheme utilized is cumbersome. Here, within the framework of ICZM-based FDEM, we propose a node binding scheme to efficiently bind the pre-discretized finite elements and thus guarantee the continuum behavior of materials in the elastic stage. The yield surfaces, controlled by ECZM, are dynamically embedded by invoking the pre-inserted cohesive elements. The effectiveness and efficiency of the proposed approach are validated and tested by performing a suite of numerical experiments. Compared with ICZM-based FDEM, the proposed approach can correctly capture material deformation and reduce the computation cost. In contrast to the existing ECZM-based FDEM, the proposed approach can overcome the frequent and complex element topology updating. Finally, this work provides a novel perspective that fully inherits the advantages of both ICZM and ECZM, but circumvents their shortcomings, which guarantees a more efficient and effective simulation of brittle material evolution from continuum to discontinuum.

58 GEOSCIENCES↗