Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Damage evolution”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Efforts to stabilize composite localization elements

Localization finite elements seek to provide a robust framework for modeling ductile failure. They utilize the same constitutive model as the bulk material through the introduction of a length scale in a specialized deformation gradient that regularizes displacement discontinuity. Similar to many other elements, localization elements exhibit locking and associated pressure oscillations under incompressible plastic flow, which is a critical issue when attempting to model pressure-driven damage evolution. These issues can be drastically improved through what are essentially reduced integration techniques for the Jacobian and pressure, but there seem to be pressure-related instabilities that persist and are specific to localization elements. This memo summarizes recent efforts to mitigate and understand this problem, mostly for the 12-node composite wedge localization element in particular. At this point, it remains unclear whether the pressure fields within any localization element can be sufficiently stabilized in order to properly model failures that include softening or fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Transmission Line Based Sensor Design for FMCW Radar Processing of Collision Damage

This report presents a framework for transmission line ( T-L ) sensor design in Frequency Modulated Continuous Wave (FMCW) radar applications, with particular emphasis on modeling collision damage effects. Classical electromagnetic (EM) propagation theory is combined with impact mechanics to develop a model for analyzing transmission line behavior as a sensor under high-velocity impact conditions. The framework includes detailed mathematical derivations, spatial and temporal damage evolution models, and practical implementation considerations.

42 ENGINEERING↗
DOE OSTI2020

A thermodynamics-based damage model for the non-linear mechanical behavior of SiC/SiC ceramic matrix composites in irradiation and thermal environments

A damage model is developed and validated with experimental data for the non-linear mechanical behavior of SiC/SiC composite materials in nuclear applications. Cyclic thermal and mechanical loading associated with neutron irradiation effects of these composites leads to wide-spread and progressive micro-cracking that leads to loss of thermal conductivity and further enhancement of thermo-mechanical damage. A physics-based model of wide-spread micro-cracking is developed within the thermodynamic framework of continuum damage mechanics. Evolution equations for damage parameters that describe the growth of continuum damage are developed, where the material variables are obtained from experiments. The model novelty is in coupling mechanical, thermal, and irradiation damage through a consistent thermodynamic framework, including loss of thermal conductivity due to the evolution of mechanically induced micro-cracks. A number of thermo-mechanical experiments were conducted to confirm model assumptions. The model is shown to be validated with out-of-pile experiments, and then implemented using commercial finite element code COMSOL to the fuel cladding problem with normal and high radiation dose cases.

Materials Science↗
DOE OSTI2020

Damage initiation and evolution in Al-Si layered microstructures under shock loading conditions at atomic scales

Designing materials with microstructural features like multiphase interfaces have shown significant promise for usage in the next generation defense and nuclear applications. The dynamic response of these interfaces in extreme environments is observed to vary with the deformability of individual phases, which could subsequently alter the favored damage nucleation sites related to spall failure. Here, we investigate the role of spacing of interfaces in a nanocrystalline layered Aluminum-Silicon system on the shock wave propagation behavior, and microstructural and defect evolution using classical molecular dynamics simulations. The simulations are carried out with different hypothetical Al/Si microstructures including variations in the distribution of Si grains as layers in a nanocrystalline Al matrix. The molecular dynamics simulations suggest that spall failure is preferentially initiated at Al/Si interfaces. In addition, for the same system volume and same concentration of Si, a higher number of layers is marked by an increase in shock wave velocity and reduced resistance to spall failure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Role of Co-injected Helium on Swelling and Cavity Evolution at High Damage Levels in Ferritic-Martensitic Steels

The influence of helium co-injection at rates from 0 to 4 appm He/dpa on swelling in ferritic-martensitic alloys T91 and HT9 was explored. Irradiations with 5.0 MeV Fe++ ions and degraded He++ ions were performed at 445°C up to damage levels of 150 dpa and helium co-injection rates of 0, 0.02, 0.2 and 4 appm He/dpa in T91, and at 460°C to a damage level of 188 dpa and helium co-injection rates of 0, 0.06 and 4 appm He/dpa in HT9. Helium was observed to enhance cavity nucleation at low damage levels, resulting in the maximum swelling at the highest helium co-injection rate. As the damage level was increased, the helium content at which swelling was maximized shifted to lower helium concentrations, ultimately resulting in the highest swelling occurring with zero helium by 150 dpa. Finally, this behavior was due to the helium-stabilized bubble microstructure that increased the cavity sink strength and reduced both cavity growth rate and swelling relative to the helium-free condition.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Genetic programming for interpretable, data-driven continuum damage models.

The damage mechanisms that lead to failure in engineering alloys have been studied extensively, but converting this knowledge into constitutive models that are suitable for engineering-scale analysis remains a challenge. Evolution laws for continuum damage have been developed in the past and have proven effective but suffer from many non-physical assumptions that inhibit the overall accuracy of the model. Further, the assumptions inherent in these existing models prevent them from being applicable to a broad class of materials. At the same time, computational models of fine-scale damage mechanisms continue to advance making it tractable to generate large training data sets through computer simulation. Data-driven machine learning approaches can leverage these data sets to avoid making limiting assumptions, and instead produce models directly from the results of microstructural simulations and/or experiments. Many of these machine learning approaches are rapid and accurate, but they offer little to no insight into the underlying relationships among state variables being discovered. Conversely, genetic programming symbolic regression (GPSR) is a machine learning method that produces analytic expressions relating the state variables, allowing maximal insight and interpretability. To that end, we propose using GPSR as a data-driven method of obtaining microstructurally informed continuum damage models. Data is generated using microstructural simulations of damage evolution, parameterized over microstructural statistics (i.e., pore shape) and nominally applied deformations. Analytic expressions for damage evolution are obtained from the data using GPSR, and these expressions are then utilized within a continuum constitutive model. Overall, this approach is a promising method of automatically obtaining analytic relations describing constitutive phenomena in a material.

Buche, Michael Robert↗
DOE OSTI2021

Attaining regularization length insensitivity in phase-field models of ductile failure

In this work, a cohesive phase-field model of ductile fracture in a finite-deformation setting is presented. The model is based on a free-energy function in which both elastic and plastic work contributions are coupled to damage. Using a strictly variational framework, the field evolution equations, damage kinetics, and flow rule are jointly derived from a scalar least-action principle. Particular emphasis is placed on the use of a rational function for the stress degradation that maintains a fixed effective strength with decreasing regularization length. The model is employed to examine crack growth in pure mode-I problems through the generation of crack growth resistance (J-R) curves. In contrast to alternative models, the current formulation gives rise to J-R curves that are insensitive to the regularization length. Numerical evidence suggests convergence of local fields with respect to diminishing regularization length as well.

42 ENGINEERING↗
DOE OSTI2024

Sound speed measurements in shock compressed cemented tungsten carbide: Evolution of elastic moduli with damage at pressures to 100 GPa

The motivation of the present study is to gain insights into the evolution of elastic properties of cemented tungsten carbides (WC) shock compressed to 100 GPa. Seven plate impact experiments—two front surface impact and five release wave overtake—are conducted to make simultaneous measurements of Hugoniot states and longitudinal sound speeds in shocked WC with 3.7wt.% cobalt binder. The sound speeds along with estimates for bulk sound speeds, obtained using the Birch–Murnaghan EoS, are analyzed to determine the elastic moduli—longitudinal, bulk, and shear—as a function of Hugoniot stress. The longitudinal and bulk sound speeds at Hugoniot states of interest are found to increase linearly with longitudinal stress. Consistent with the increase in sound speeds, the longitudinal and bulk moduli also increase with Hugoniot stress; however, the increase in longitudinal modulus is modest when compared to predictions of theoretical models that account for pressure and temperature dependence of elastic moduli, but with no damage. The shear moduli remain nearly constant at ~318 GPa over the range of Hugoniot states investigated. These values are, however, much lower than those predicted by the Steinberg–Guinan model with no damage. Poisson’s ratio decreases initially from its ambient value of 0.208 to ~0.199 for Hugoniot stress ≤10 GPa indicating consolidation of the WC microstructure with low initial stress; however, with an increase in Hugoniot stress to ~100 GPa, Poisson’s ratio increases to ~0.317, indicating degradation of shear moduli with increasing stress. The product of density and Grüneisen parameter (ρΓ), after an initial spike, remains nearly constant for volumetric strains ≥0.07. The maximum average temperature rise is estimated to be ~286°C at the highest Hugoniot stress employed in the study.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mechanism-Based Energy Regularization in Computational Modeling of Quasibrittle Fracture

Abstract Quasibrittle materials are featured by a strain-softening constitutive behavior under many loading scenarios, which could eventually lead to localization instability. It has long been known that strain localization would result in spurious mesh sensitivity in finite element (FE) simulations. Previous studies have shown that, for the case of fully localized damage, the mesh sensitivity can be mitigated through energy regularization of the material constitutive law. However, depending on the loading configuration and structural geometry, quasibrittle structures could exhibit a complex damage process, which involves both localized and diffused damage patterns at different stages of loading. This study presents a generalized energy regularization method that considers the spatial and temporal evolution of damage pattern. The method introduces a localization parameter, which describes the local damage pattern. The localization parameter governs the energy regularization of the constitutive model, which captures the transition from diffused to localized damage during the failure process. The method is cast into an isotropic damage model, and is further extended to rate-dependent behavior. The energy regularization scheme is directly incorporated into the kinetics of damage growth. The model is applied to simulate static and dynamic failures of ceramic specimens. It is shown that the present model is able to effectively mitigate the spurious mesh sensitivity in FE simulations of both types of failure. The present analysis demonstrates the essential role of mechanism-based energy regularization of constitutive relation in FE simulations of quasibrittle fracture.

Mechanics↗
DOE OSTI2022

Effect of Stress on Irradiation Responses of Highly Oriented Pyrolytic Graphite

The effect of stress on irradiation responses of highly oriented pyrolytic graphite (HOPG) was studied by combing molecular dynamics (MD) simulation, proton irradiation, and Raman characterization. MD simulations of carbon knock-on at energies < 60 eV were used to obtain average threshold displacement energies (E¯d) as a function of strain ranging from 0 to 10%. Simulations at a higher irradiation energy of 2–5 keV were used to study the effect of strain on damage cascade evolution. With increasing tensile strain, E¯d was reduced from 35 eV at 0% strain to 31 eV at 10% strain. The strain-reduced E¯d led to a higher damage peak and more surviving defects (up to 1 ps). Furthermore, high strains induced local cleavage around the cavities, as one additional mechanism of damage enhancement. Experimentally, HOPG film was folded, and the folded region with the maximum tensile stress was irradiated by a 2 MeV proton beam. Raman characterization showed significantly enhanced D to G modes in comparison to the stress-free irradiation. Based on the strain dependence of E¯d and the Kinchin–Pease model, a formula for displacement estimation under different tensile strains is proposed. The stress effects need to be considered in graphite applications in a reactor’s harsh environment where both neutron damage and stress are present.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Photoluminescence probing of light absorption centers at silica laser damage

We use photoluminescence (PL) imaging to study damage growth precursors within laser damage sites on the surface of silica. Damage site evolution is induced by multiple shots of UV nanosecond pulsed laser at various energy densities and monitored throughout the early stages of growth. Wide-field PL imaging rapidly locates microscopic light absorption centers within the silica damage site. Our quantitative analysis shows that damage sites with strong local PL intensity show a higher probability of growth upon subsequent laser pulses. Scanning electron microscopy (SEM) paired with a study of PL spectrum shows that the strong PL intensity appears from the subsurface fractures with high defect density, which provides a local light absorption center leading to significant damage growth. We believe that this result offers an efficient optical damage mitigation strategy by providing a rapid and non-destructive optical inspection approach.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A novel peridynamics-based approach to predict pharmaceutical tablet robustness

The pharmaceutical drug product development process can be greatly accelerated through the use of modeling and simulation techniques to predict the manufacturability and performance of a given formulation. The anticipation and possible mitigation of tablet damage due to manufacturing stresses represents a specific area of interest in the pharmaceutical industry for predicting formulation and tableting performance. While the finite element method (FEM) has been extensively used for predicting the mechanical behavior of powder material in the compaction processes, a shortcoming of the approach is the inherent difficulty to predict discontinuities (e.g., damage or cracking) within a tablet as FEM is a continuum-based approach. In this work, we propose a novel method utilizing peridynamics (PD), a numerical method that can capture discontinuities such as tablet fracture, to predict the evolution of damage and breakage in pharmaceutical tablets. The approach links (1) the finite element method – to elucidate the behavior of powders during die compaction – with (2) the peridynamics modeling technique – to model the discontinuous nature of damage and predict tablet breakage during the critical stages of unloading and ejection from the compression die. This short communication presents a proof of concept including a workflow to calibrate the linked FEM-PD simulation models. Further, it demonstrates promising results from a preliminary experimental validation of the approach. Following further development, this approach could be used to guide the optimization of compression processes through targeted changes to formulation material properties, compression process conditions, and/or tooling geometries to deliver improved process efficiency and tablet robustness.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling creep behavior in ceramic matrix composites

Here, in this work, a three-dimensional viscoplasticity formulation with progressive damage is developed and used to investigate the complex time-dependent constituent load transfer and progressive damage behavior in ceramic matrix composites (CMCs) subjected to creep. The viscoplasticity formulation is based on Hill's orthotropic plastic potential, an associative flow rule, and the Norton-Bailey creep power law with Arrhenius temperature dependence. A fracture mechanics-informed isotropic matrix damage model is used to account for CMC brittle matrix damage initiation and propagation, in which two scalar damage variables capture the effects of matrix porosity as well as matrix property degradation due to matrix crack initiation and propagation. The Curtin progressive fiber damage model is utilized to simulate progressive fiber failure. The creep-damage formulation is subsequently implemented as a constitutive model in the generalized method of cells (GMC) micromechanics formulation to simulate time-dependent deformation and material damage under creep loading conditions. The developed framework is used to simulate creep of single fiber SiC/SiC microcomposites. Simulation results are in excellent agreement with experimental and numerical data available in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A thermodynamically consistent finite strain phase field approach to ductile fracture considering multi-axial stress states

Phase field models for ductile fracture have gained significant attention in the last two decades due to their ability in implicitly tracking the nucleation and propagation of cracks. However, most crack phase field formulations for elastoplastic solids focus only on the effects of plastic deformation, and do not consider the different multi-axial stress states that may arise in practical designs. In this work, a thermodynamically consistent phase field approach coupled with finite strain plasticity, considering multi-axial stress states is presented. In order to account for the coupling between plasticity and stress states, the Stress-Weighted Ductile Fracture Model (SWDFM) is utilized. The SWDFM represents a criterion for predicting ductile crack initiation under both monotonic and cyclic loadings based on histories of an internal plastic variable, stress triaxiality, and the Lode angle parameter. The excellent performance of the SWDFM for predicting ductile crack initiation motivates for its incorporation into a phase field approach for predicting both crack initiation and propagation through degradation of the fracture toughness. Moreover, based on the second law of thermodynamics, exact requirements are imposed on the rate at which the fracture toughness can evolve. Here, a novel function for degrading the plastic yield surface during the evolution of damage is introduced. This function, in line with experimental observations, leads to an accumulation of plastic deformation in damaged regions of a solid, and avoids numerical instabilities arising from concentrations of large plastic deformations in severely damaged regions. For validating the proposed model, results of computational simulations are compared to data from selected tests considering different multi-axial stress states. Comparisons of the numerical results with data from laboratory experiments demonstrate the capabilities of the proposed framework.

42 ENGINEERING↗
DOE OSTI2022

In-situ imaging of flexure-induced fracture in tape-laminate composites using high-resolution X-ray computed tomography

Here, this study presents development of an in-situ flexure test for imaging progressive inter- and intralaminar fracture in tape-laminate composites using X-ray computed tomography (CT). The intent of this test is to provide detailed experimental observations of ply-level damage that can be used to validate existing, and develop new, progressive damage analysis (PDA) tools. The test consists of a vertically mounted specimen which is flexed using two eccentric compressive loads using an in-situ uniaxial load stage. The flexure specimen contains a starter notch mid-span which promotes initiation of composite failure within the X-ray field of view. Specimens with two different laminate stacking sequences were tested, imaged, analyzed. For a quasi-isotropic laminate with large angle changes between adjacent plies, there was near simultaneous growth of transverse cracks and delaminations below the midplane of the laminate. For a laminate with small angles between adjacent plies, there was extensive formation of transverse crack networks which penetrated the laminate thickness without delamination growth. In addition to imaging fracture, the X-ray CT data from both specimen types were used to quantify the variability in thickness and analyze the local orientation of individual plies. Overall, the proposed test and the image-data analysis methodology provided an important insight into the fracture processes in tape laminates and highlights the inherent ply-level geometrical variabilities that should be accounted for in PDA simulations.

42 ENGINEERING↗
DOE OSTI2021

A mechanistic damage model for solid oxide fuel cell ceramic materials - Part I: Constitutive modeling

A multiscale mechanistic damage model for solid oxide fuel cell (SOFC) ceramics that combines micromechanics of stiffness reduction due to material porosity change and microcracking with a continuum damage mechanics (CDM) description for the evolution of microcracks up to fracture is developed in this work. The model also accounts for volumetric swelling of the anode due to redox cycling. Porosity change and swelling during redox combined with normal operating thermomechanical loads could increase the stress, strain and damage distributions in SOFC stacks leading them to failure. The proposed model involves three governing parameters: (i) the porosity in terms of the pore or void volume fraction, (ii) material swelling magnitude prescribed at a given loading step, and (iii) the damage variable that describes microcracking caused by thermomechanical loads and swelling. Pores and microcracks in the ceramics are modeled as randomly distributed ellipsoidal inclusions with negligible stiffness by an Eshelby-Mori-Tanaka formulation combined with an inclusion orientation distribution method. Microcracking damage evolution is described by a CDM formulation. Volumetric swelling is treated in a similar way to thermal expansion in the constitutive relations. After validation and numerical checks, the damage model is used to analyze a simplified SOFC positive-electrode/electrolyte/negative-electrode (PEN) structure subjected a redox cycle in addition to thermomechanical loads experienced during normal operation of a SOFC stack.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Dependence of simulated radiation damage on crystal structure and atomic misfit in metals

This study investigates the evolution of radiation damage in three metals in the low temperature and high radiant flux regime using molecular statics and a Frenkel pair accumulation method to simulate up to 2.0 displacements per atom. The metals considered include Fe, equiatomic CrCoNi, and a fictitious metal with similar bulk properties to the CrCoNi composed of a single atom type referred to as an A-atom. CrCoNi is found to sustain higher concentrations of dislocations than either the Fe or A-atom systems and more stacking faults than the A-atom system. The results suggest that the difference between the concentrations of vacancies and interstitials is substantially smaller for CrCoNi than the A-atom system, perhaps reflecting that the sink capture radius is smaller in CrCoNi due to the roughened potential energy landscape. A model that partitions the major contributions from defects to the stored energy is described, and serves to highlight a general need for higher fidelity approaches to point defect identification.

36 MATERIALS SCIENCE↗