Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ductile damage”

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 OSTI2026

Effect of macroscopic surface defects on dynamic damage: An experimental and numerical study

This study examines the impact of macroscopic surface defects on the dynamic ductile damage behavior of polycrystalline metals using plate-impact experiments. Defects of various shapes (flat, round, and point) were manufactured on the free or impact surfaces of annealed copper specimens. The experiments were diagnosed with photon Doppler velocimetry measurements and soft recovery techniques. The experimental results revealed that defect shape and location significantly affect velocity–time profiles and void distribution. In order to understand the dynamics of shock propagation and corresponding ductile damage evolution within the target specimens, the impact experiments were modeled with a calibrated Tepla model, a dynamic ductile damage model for polycrystalline metals [Nguyen et al., Int. J. Solids Struct. 329, 113833 (2026)]. Overall, our resulting simulated velocity showed good agreement with measured velocity, and our simulated porosity distributions qualitatively matched experimental data. Based on our simulation results, defects on the free surface were found to distort rarefaction waves and therefore the corresponding spall planes. On the other hand, defects on the impact surface generated a delayed shock when the flyer plate and the defected target area were in contact, leading to distortion of the spall plane. The distortion of the spall plane resulted in a non-uniform distribution of voids within the defected specimens. Lastly, we examined the role of local work hardening due to the defect manufacturing process on the velocity and porosity distribution, using Tepla simulations with a simplified representation of local hardening. This investigation highlights the importance of defect geometry, location, and local hardening associated with defect generation in dynamic ductile damage processes.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Reformulating a Gurson-based dynamic damage model and demonstrating improved predictive power and numerical robustness

The original Tepla (TEnsile PLAsticity) ductile damage model, based on the Gurson yield surface, has long been used to model damage evolution and material failure under dynamic loading. Unfortunately, Tepla suffered from mesh sensitivity, numerical instability, and limited predictive capability. Here, we theoretically reformulate Tepla to address these issues. We especially focus on the prediction of porosity, which is the key state variable used for modeling ductile damage, as compared to more easily measured surface velocities, which are at best an indirect measure of damage. Key model changes include separating the viscosity during volumetric void growth from underlying shear strength behavior and switching to an iterative bisection solver. The new Tepla is then calibrated on incipient spall experiments on half-hard copper and tantalum, which demonstrate its ability to simultaneously fit the model to recovered porosity distributions and measured surface velocities, a stringent test. Improved numerical behavior, such as greatly reduced mesh sensitivity, is also shown in those simulations. Finally, the new Tepla model is applied to several high-explosive loaded, sweeping wave experiments, showing the ability of the model to predict behavior on tests with significantly different loading conditions and histories than the calibration data.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling detonation with CartaBlanca simulations

The accurate modeling of high explosive (HE) detonation and the con- sequent large solid deformation, failure, plastic flow, porosity growth, and shock wave propagation is important because simulations can capture spatial and temporal features that experimental diagnostics cannot capture. However, the simulation of the explosive event poses challenges to a computational scientist. These include the accurate modeling of ductile damage, crack formation, plastic deformation, as well as physical and nu- merical instabilities. The material response can be history-dependent and subject to large material deformation. Our research simulates the impact of a high explosive (Detasheet) onto a tantulum metal plate. We performed the simulations using CartaBlanca at different mesh resolutions. We decided the study would be impactful if we perform the simulations with the Material Point Method. Differences were observed at the different mesh resolutions in velocity and nodal stress magnitude, so increased mesh resolutions may be required. In addition, the Discontinuous Galerkin method would be needed to account for the large gas expansion.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1973

The effect of yield strength and ductility to fatigue damage

The cumulative damage of aluminium alloys with different yield strength and various ductility due to seismic loads was studied. The responses of an idealized beam with a centered mass at one end and fixed at the other end to El Centro's and Taft's earthquakes are computed by assuming that the alloys are perfectly elastoplastic materials and by using numerical technique. Consequently, the corresponding residual plastic strain can be obtained from the stress-strain relationship. The revised Palmgren-Miner cumulative damage theorem is utilized to calculate the fatigue damage. The numerical results show that in certain cases, the high ductility materials are more resistant to seismic loads than the high yield strength materials. The results also show that if a structure collapse during the earthquake, the collapse always occurs in the very early stage.

Yeh, H. Y.↗
DOE OSTI2024

Damage evolution and ductile fracture of commercially-pure titanium sheets subjected to simple tension and cyclic bending under tension

This paper describes a study into the evolution of damage in commercial-purity titanium (CP–Ti) sheets subjected to cyclic bending under tension (CBT) and uniaxial tension (simple tension, ST). Sections were taken from sheets that were strained to various levels under both methods and were imaged using X-ray computed tomography (XCT) to reveal insights into the size, density, and distribution of microvoids in the sheets. Digital image correlation (DIC) was also used to observe the surface strain of CBT and ST sheets during testing. These results showed that elongation to failure (ETF) for CBT deformation is about 2.5× higher than for ST, local longitudinal strains in the bulk of the CBT sample are around 1.3× higher than the peak strain in the ST necked region, and 1.7× higher in the localized region of CBT failure. It was found that the volume-density of voids in both sheets followed a similar exponential increase with strain, reaching nearly 45× higher in the CBT than the ST sheets before the onset of failure, mainly due to the higher strain levels achieved. A higher volume density of voids developed in the center of sheets at high levels of strain, for both processes, with the density in the sheet center becoming approximately double that found near the edges. Scanning electron microscopy (SEM) was used to examine the fracture surface of CBT and ST sheets after failure. The observations are presented and discussed highlighting the CBT process as effective to delay ductile fracture of the CP-Ti sheets.

Commercially pure titanium↗
DOE OSTI2026

Assessing the accuracy of time-fraction and ductility exhaustion approaches for creep-fatigue damage prediction through feature-test validation

Determining creep-fatigue damage formation is critical for elevated temperature components integrity. This study evaluates creep-fatigue assessments with emphasis on differences between creep damage models. Evaluated are the time-fraction model and the classical and stress-modified ductility exhaustion creep damage models. This work extends the domain of stress-modified ductility exhaustion models by introducing and validating such formalism to Ni-based alloys. The fidelity of the assessments was benchmarked against uniaxial creep-fatigue and multiaxial feature tests of Alloy 617. For uniaxial specimens, best estimate predictions rank ductility exhaustion as the most accurate and precise and time-fraction as markedly conservative. For feature tests, ductility exhaustion predictions are within < 4.0 times difference, whereas time-fraction underpredicts life by factors of 7–14. The observations suggest ductility exhaustion models as alternative to time-fraction models in design codes for situations requiring characterization of the design margin. Further work in employing such models to assess other relevant phenomena (e.g., stress relaxation cracking) is discussed.

creep-fatigue↗
DOE OSTI2024

Modeling plasticity-mediated void growth at the single crystal scale: A physics-informed machine learning approach

Modeling the evolution of voids during plastic flow as well as their effects on plastic dissipation is critical for both component manufacturing and lifetime estimation purposes. To this end, we propose a rate-dependent constitutive model to homogenize the effects of semi-randomly distributed voids on single crystal plasticity whilst capturing void interaction and plastic anisotropy. Here, this present work focuses on the case of face centered cubic crystals to introduce an anisotropic gauge function applicable within the crystal plasticity formalism. The approach combines analytical methods to describe the micromechanics of the system in combination with symbolic regression to capture analytically intractable mechanisms from data. The hybrid framework uses a physics-informed genetic programming-based symbolic regression algorithm to solve a multiform optimization problem simultaneously producing a new gauge function and a new strain rate equation. This is also a multi-objective optimization problem with many competing objectives. A new search and selection step is introduced to the genetic algorithm that promotes convergence toward a global solution that better satisfies all the objectives. Overall, the symbolic equations produced leverage data-driven methods to achieve greater accuracy than comparable alternatives on an analytically intractable problem while maintaining model transparency.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A phase-field fracture formulation for generalized standard materials: The interplay between thermomechanics and damage

Accurately modeling fracture of ductile materials poses open challenges in the field of computational mechanics due to the multiphysics nature of their failure processes. Integrating the interplay between thermodynamics and damage into ductile fracture models is vital for predicting critical failure modes. Here, in this paper, we develop a versatile phase-field (PF) framework for modeling ductile fracture, taking into account finite-strain elasto-plasticity. The framework stems from a variational formulation of constitutive relations for generalized standard materials (GSMs), whose response is described by a Helmholtz free energy and a dissipation pseudo-potential. Its variational structure is based on a minimum principle for a functional that expresses the sum of power densities for reversible and irreversible processes. By minimizing this functional with a constraint on a von Mises yield function, we derive the evolution equation for the equivalent plastic strain and an associative flow rule. This constrained optimization problem is analytically solved for a wide class of thermo-viscoplasticity models. The key innovations of the current work include (i) a cubic plastic degradation function that accounts for a non-vanishing damage-dependent yield stress, (ii) closed-form expressions of the Helmholtz free energy and dissipation pseudo-potential for three thermo-viscoplasticity models, (iii) an extended Johnson–Cook plasticity model with a nonlinear hardening law, and (iv) a plastic work heat source that depends on the plastic degradation function and a variable Taylor–Quinney (TQ) coefficient. The capabilities of the proposed framework are tested with the aid of four ductile fracture problems, including the Sandia Fracture Challenge. In each of these problems, we examine the evolution of relevant field variables such as the PF order parameter, the equivalent plastic strain, the temperature, and the internal power dissipation density, in addition to the overall structural response quantified by the force–displacement curve. These numerical studies demonstrate that the proposed framework effectively represents ductile fracture, yielding computational results that exhibit good agreement with experimental data.

36 MATERIALS SCIENCE↗
NASA NTRS1972

Effect of cryogenic irradiation on NERVA structural alloys

Several alloys (Hastelloy X, AISI 347, A-286 bolts, Inconel 718, Al 7039-T63 and Ti-5Al-2.5Sn ELI) were irradiated in liquid nitrogen (140 R) to neutron fluences between 10 to the 17th power and 10 to the 19th power nvt (E greater than 1.0 Mev). After irradiation, tensile properties were obtained in liquid nitrogen without permitting any warmup except for some specimens which were annealed at 540 R. The usual trend of radiation damage typical for materials irradiated at and above room temperature was observed, such as an increase in strength and decrease in ductility. However, the damage at 140 R was greater because this temperature prevented the annealing of radiation-induced defects which occurs above 140 R.

Dixon, C. E.↗
DOE OSTI2025

The Ductility of 49Fe-49Co-2V Soft Magnetic Alloy Bar: Surface Effects and Test Methods

The tensile ductility of 49Fe-49Co-2V (Hiperco® 50A) bar was investigated in both as-received and heat-treated conditions. The as-received/machined specimens exhibit very low ductility compared to samples where heat treatment was the final step prior to testing. Microstructural characterization showed that internal residual strain from bar processing and, most importantly, surface machining damage, cause lower elongation in the as-received material. Because fracture of this intermetallic alloy initiates at the surface, it is particularly susceptible to surface machining damage, i.e., the near-surface region has already exhausted most of its ability to accumulate tensile strain. During heat treatment, the internal residual strain and near-surface machining damage are eliminated and ductility is improved, despite a higher degree of crystallographic ordering in the heat-treated condition (which typically lowers ductility). Furthermore, if machining is again performed after heat treatment, the material again exhibits brittle behavior, even with only light touch-up machining passes. Here, in this work, methods of tensile strain measurement were investigated, namely conventional knife-edge extensometry and noncontact digital image correlation (DIC) on heat-treated material. For clip-on knife-edge extensometry, the range of failure strain was 2.5-5.5% for heat-treated Hiperco. For noncontact methods, ductility up to 7% was observed. The results highlight the tendency for the alloy to fail at surface imperfections, even those produced by application of the extensometer itself. Noncontact laser extensometry is recommended for determining the intrinsic ductility of the alloy. A method of laser surface modification was developed which increased ductility by ~ 100% compared to unmodified samples. The high cooling rates achieved during laser surface processing can bypass the ordering reaction and produce a ductile disordered structure at the surface that exhibits ductile fracture characteristics.

EBSD↗
NASA NTRS1980

Cyclic plasticity and failure of structural components

An analytical approach for low-cycle fatigue prediction is presented. The approach combines a cyclic plasticity model with the finite element method and a damage accumulation criterion for ductile metals. The cyclic plasticity model is based on the concept of the combination of several yield surfaces. The surfaces are related to the material uniaxial stress-strain curve idealized by piecewise linear segments. The damage criterion is based on the Coffin-Manson formulae modified for the mean stress variation effect. It is extended to the multiaxial varying stress-strain field and applied for both the crack initiation and the crack growth processes. The stable slow crack growth rate is approximated by the damage accumulation gradient computed from the cracked finite element models. This procedure requires fatigue testing data of only smooth specimens under constant strain amplitudes. The present approach is illustrated by numerical examples of an aircraft wing stiffened panel subjected to compression, which causes material yielding and residual tension.

Kalev, I.↗
DOE OSTI2025

Experimental data for damage mechanics simulation challenge

While there are many computational approaches for simulating damage in rock and other materials, few have been ground truth tested with either known experimental data or with blind data sets. Here, in this work, we present a bench-mark laboratory data set for a damage mechanics challenge to compare computational approaches on damage evolution in brittle-ductile materials. The samples were fabricated through additive manufacturing to produce repeatable specimens designed to fail in controlled ways. The failure was induced in the samples using a 3-point bending test to produce different Modes such as Mode I and mixed Modes including I-II, I-III and I-II-III Modes to generate a calibration data set and a blind challenge data set. Data collected included spatial and temporal measurements from traditional digital load–displacement sensors, 2D digital image correlation measurement to map surface deformations, 3D X-ray microscopy to ground-truth the crack-failure geometry, and laser profilometry to capture surface roughness. The data sets are available, on a data repository, to the community to advance computational models to improve our ability to predict damage in brittle-ductile materials.

3-point bending↗
NASA NTRS1983

Particulate erosion mechanisms

Particulate damage and erosion of ductile metals are today plaguing design and field engineers in diverse fields of engineering and technology. It was found that too many models and theories were proposed leading to much speculation from debris analysis and failure mechanism postulations. Most theories of solid particle erosion are based on material removal models which do not fully represent the actual physical processes of material removal. The various mechanisms proposed thus far are: melting, low-cycle fatigue, extrusion, delamination, shear localization, adhesive material transfer, etc. The experimental data on different materials highlighting the observed failure modes of the deformation and cutting wear processes using optical and scanning electron microscopy are presented. The most important mechanisms proved from the experimental observations of the specimens exposed to both spherical and angular particles are addressed, and the validity of the earlier theories discussed. Both the initial stages of damage and advanced stages of erosion were studied to gain a fundamental understanding of the process.

Veerabhadrarao, P.↗
NASA NTRS2020

Analysis of Open Hole Compression Specimens Using the CompDam Continuum Damage Mechanics Model

This report documents a validation study on the standard Open Hole Compression (OHC) laminate test specimen conducted as part of the NASA Advanced Composites Project (ACP). Tests were conducted on OHC specimens with hard, quasi-isotropic, and soft layups using digital image correlation and X-Ray computed tomography to capture the structural response and damage evolution. Progressive damage models were constructed for use with the CompDam continuum damage mechanics code following the best practices established during the ACP. Detailed interrogation of the analysis results and comparison with experimental measurements provide a basis for assessing the capability of the modeling approach for OHC. The structural response is found to be captured well, with strengths predicted within 3\% of the experimental values for hard and quasi-isotropic laminates. In the soft laminate, the model predicts failure to be more brittle than the nonlinear, ductile response that was measured. Damage states extracted from the models at the same load level as test measurements are overlaid to show directly the similarities and differences between test and analysis results. Studying the damage evolution predicted by the analysis reveals that the failure process is a competition between fiber damage and delamination/sub-laminate buckling, with fiber damage dominating the collapse in the hard laminate and sub-laminate buckling governing in the soft laminate. Finally, a series of parametric studies varying numerical solution parameters (mesh size, mass scaling) and physical properties (fiber direction compressive strength and toughness) reveal sensitivities and deficiencies of the model. To the authors' knowledge, this study is the first for OHC specimens to include detailed evaluation of damage mode interactions, direct overlay of predicted and measured damage states, and sensitivity of the predicted results to difficult-to-measure fiber direction material properties.

Continuum damage mechanics↗
NASA NTRS1972

Correlation of irradiation data using activation fluences and irradiation temperature.

A new method of correlating radiation damage data is tested using actual measured data taken from the open literature. This method, the activation fluence method, was found to be as accurate as other contemporary models with which it was compared. The new method also has several advantages over the other methods. The method employs a new entity, the activation fluence (time-integrated specific activation rate), as the independent variables in a regression model. Temperature at which the irradiation takes place is also a variable. Although the method was tested for a specific type of damage (change in nil-ductility transition temperature for A302-B steel) it has no inherent restrictions and is limited only by the imagination of the user.

Lynch, J. H.↗
DOE OSTI2025

An FFT-based micromechanical model for gradient enhanced brittle fracture

Damage models incorporated within FFT-based micromechanical methods have received much attention recently because of the need to better understand and predict brittle and ductile fracture. An important aspect of a damage model is non-local regularization, which removes the mesh dependence of the predictions that otherwise become physically unacceptable upon grid refinement. In this work, the Helmholtz-type equation for non-local gradient regularization of a damage model on a distorted grid is solved using an FFT-based approach. Further, the resulting system of equations is solved using the Jacobi iterative method. The model is applied to simulate brittle fracture of an intermetallic. The influence of the time and space discretization, the length-scale parameter, and intermetallic crystallographic orientation on crack evolution is studied.

36 MATERIALS SCIENCE↗