Search NASA⌕ Search

SEARCH · Search NASA

Results for “Cohesive zone”

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 OSTI2021

Simulations of crack extensions in arc-shaped tension specimens of uncharged and hydrogen-charged 21-6-9 austenitic stainless steels using cohesive zone modeling with varying cohesive parameters

Crack extensions in uncharged and hydrogen-charged side-grooved A(T) specimens of conventionally forged 21-6-9 austenitic stainless steels are simulated using the cohesive zone modeling approach. Two-dimensional plane strain finite element analyses with fixed cohesive parameters are first conducted to fit the experimental load-displacement curves. Similar analyses using varying cohesive parameters as functions of the crack extension are then conducted to fit the experimental load-crack extension and crack extension-displacement curves. The computational results with varying cohesive parameters can fit very well the experimental data. The computational results also indicate that the average cohesive energy for the hydrogen-charged A(T) specimen is lower than that for the uncharged A(T) specimen.

42 ENGINEERING↗
DOE OSTI2024

Unified cohesive zone model (UCZM) for fracturing and fragmenting solids

Here, a Unified Cohesive Zone Model (UCZM), which inherits most of the advantages while overcoming the shortcomings of existing Cohesive Zone Models (CZMs), is proposed. Similar to the traditional extrinsic CZM approach, UCZM dynamically inserts the cohesive elements into the system based on local material states (e.g., stress, strain). However, the transition from continua to discontinua is smoothly achieved, thereby eliminating the “time-discontinuous” issue seen in the extrinsic CZM. Moreover, within the novel UCZM framework, the point of transition from continua to discontinua is controllable through the introduction of crack initialization criteria. As a result, the UCZM allows any material models (e.g., elastic, plastic, damage models) for continuum solids and for discrete fracture behavior to work together. In essence, both an enhanced extrinsic cohesive zone model and an intrinsic cohesive zone model can be represented by the proposed unified model. The proposed UCZM has been verified through different numerical examples. The work demonstrates that the UCZM is a highly effective approach for modeling fracture and fragmentation processes in solids.

42 ENGINEERING↗
DOE OSTI2024

Bayesian calibration and uncertainty quantification of a rate-dependent cohesive zone model for polymer interfaces

In this work we present a rate-dependent cohesive zone model for the fracture of polymeric interfaces and performs a Bayesian calibration, an uncertainty quantification, and a sensitivity analysis for the model. The proposed cohesive zone model accounts for both reversible elastic and irreversible rate-dependent separation sliding deformation at the interface. The viscous dissipation due to the irreversible opening at the interface is modeled using elastic-viscoplastic kinematics that incorporates the effects of strain rate. Inverse calibration of parameters for such complex models through trial and error is challenging due to the large number of parameters of the model. Moreover, the calibrated parameter values are often non-unique and uncertain when the available experimental data is limited. To tackle this challenge, we employ a Bayesian calibration approach to identify parameters from experimental data, the resulting parameters significantly enhance the accuracy of the model. To quantify the uncertainty associated with the inverse parameter estimation, a modular Bayesian approach is employed to calibrate the unknown model parameters, accounting for the parameter uncertainty of the cohesive zone model. The advantages of the Bayesian calibration over a deterministic parameter fit are demonstrated. Further, to quantify the model uncertainties, such as incorrect assumptions or missing physics, a discrepancy function is introduced, which significantly improves the model’s prediction. Finally, the total uncertainty of the model is quantified in a predictive setting. A sensitivity analysis is performed to assess how changes in the input variables of the model affect the peak load, facilitating the identification of a concise set of highly influential parameters. The present approach can be used for calibration and uncertainty quantification for other complex computational mechanics models. It should also facilitate the designing of interface materials under uncertainty.

42 ENGINEERING↗
DOE OSTI2026

Finite deformation implementation of a mixed-mode single-integral type cohesive zone with reorienting surfaces of separation

To model material ductile failure and crack propagation, cohesive zone elements can be embedded along potential fracture paths in a finite element simulation. When damage criteria are met, elements in the mesh decohere, simulating the formation and propagation of a crack. In this paper, we present a novel computational algorithm based on finite deformation theory, essential to modeling crack initiation and growth in solids undergoing large deformations. This new algorithm was formulated within a Lagrangian frame of reference to extend previous cohesive zone algorithms to include modeling crack growth in finite deformation contexts. The local coordinate system, necessary for defining an embedded cohesive zone, is constructed based upon the current configuration and is updated within the nonlinear iteration process, thereby resulting in the convergence of the solution for a growing crack in a large deformation quasi-static setting. The model’s accuracy was demonstrated by comparing finite element model simulation results with the analytic case of a constant surface separation, as shown in the verification examples. The power and efficacy of the algorithm to capture large deformations during crack growth were then demonstrated with a double cantilever beam example case. It indicates that the model can be applied to a variety of physical circumstances for predicting crack initiation and growth with delamination and fracture.

42 ENGINEERING↗
DOE OSTI2025

A cohesive zone treatment for the material point method involving problems of large deformation and damage

A new algorithm is described that permits the use of cohesive zones in the material point method for problems involving large deformation and fracture. In contrast to previous cohesive zone implementations, this method does not utilize massless surface-element particles. Instead, cohesive tractions are computed using the shape function mappings from a reference grid configuration in combination with explicitly defined particle surface normals and surface positions. These normals and relative surface positions are updated each time step according to particle deformation. The tractions are converted to cohesive forces using the nodal areas and mapped back to particles using the same reference shape function mappings. These forces are then remapped by conventional particle-to-grid interpolation as external forces using the current-configuration shape-function mappings. This allows highly compliant cohesive zones to function over jump displacements larger than a grid cell. Upon damage, these interfaces can revert to conventional multi-field contact surfaces. This approach is general and readily applies to two and three dimensions as well as being compatible with damage-field gradient partitioning offering exceptional computational flexibility. The framework for this method enables other capabilities, such as improved contact precision using explicitly defined surface normals and positions, and a method to mitigate spurious material damage at weak discontinuities between stiff brittle materials and soft or compliant materials.

Cohesive zone↗
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↗
DOE OSTI2024

Poromechanical cohesive interface element with combined Mode I-II cohesive zone elastoplasticity for simulating fracture in fluid-saturated porous media

A combined Mode I-II cohesive zone (CZ) elasto-plastic constitutive model, and a two-dimensional (2D) cohesive interface element (CIE) are formulated and implemented at small strain within an ABAQUS User Element (UEL) for simulating 2D crack nucleation and propagation in fluid-saturated porous media. Here, the CZ model mitigates problems of convergence for the global Newton-Raphson solver within ABAQUS, which when combined with a viscous stabilization procedure allows for simulation of post-peak response under load control for coupled poromechanical finite element analysis, such as concrete gravity dam stability analysis. Verification examples are presented, along with a more complex ambient limestone-concrete wedge fracture experiment, water-pressurized concrete wedge experiment, and concrete gravity dam stability analyses. A calibration procedure for estimating the CZ parameters is demonstrated with the limestone-concrete wedge fracture process. For the water-pressurized concrete wedge fracture experiment it is shown that the inherent time-dependence of the poromechanical CIE analysis provides a good match with experimental force versus displacement results at various crack mouth opening rates, yet misses the pore water pressure evolution ahead of the crack tip propagation. This is likely a result of the concrete being partially-saturated in the experiment, whereas the finite element analysis assumes fully water saturated concrete. For the concrete gravity dam analysis, it is shown that base crack opening and associated water uplift pressure leads to a reduced Factor of Safety, which is confirmed by separate analytical calculations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Mortar-based Cohesive Zone Model with Application to TRISO Particle Debonding

The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering-scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. This report outlines work towards the development of an integrated, mortar-based cohesive zone framework with mechanical contact. We model debonding of tri-structural isotropic (TRISO) layers using a mortar approach, extending the framework used for mechanical contact, using weighted quantities. A bilinear mixed mode traction approach, already used in Jiang et al. (2021), is applied analogously herein to model debonding between TRISO particle layers. The examples include the debonding of the inner pyrolytic carbon (IPyC) from the silicon carbide (SiC) layers, as well as partial and full debonding of the buffer from the IPyC layers. Both 2D and 3D examples are utilized to demonstrate the superior numerical performance of the new mortar approach, which allows for modeling arbitrarily meshed interfaces (i.e., non-matching discretizations).

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Fracture Analysis of Cohesive Zone Models for Modeling Residual Stress Induced Delamination in Composite Structures

A fracture study of coupon-scale composite cylinders with embedded defects was conducted with an objective to assess and validate a modeling approach using two available cohesive material models. The study included experimental and simulation evaluations of initiation of crack growth and progression. Interrupted thermal experiments used acoustic emissions monitoring to identify the onset of crack progression during each cooling interval and ultrasonic scanning provided images of defect growth. Verification, validation, and uncertainty quantification (VVUQ) processes were performed in the assessment of the simulation predicted temperature at which crack propagation begins (quantity of interest). The Sobol sensitivity analysis identified the hoop direction elastic modulus in the carbon fiber reinforced polymer (CFRP) plies as the most influential parameter for simulations using both cohesive models, accounting for at least 70% of the variation in the temperature at crack propagation. The UQ temperature range for the Tvergaard-Hutchinson model was higher (more conservative) than the experimental acoustic measurement indicators of crack progression, while the temperature range for the Thouless-Parmigiani model enveloped the experimental data points for the primary defect size of 0.75 x 1 in. The simulations could not capture the stable crack growth indicated in the experiments. This is likely due to the models’ inability to represent anisotropic fracture toughness attributed to the structure of the orthotropic fiber weave in a woven composite laminate.

42 ENGINEERING↗
DOE OSTI2026

Non-conformal interface-cohesive modeling with the shifted boundary method

The accurate simulation of boundary- and interface-dominated problems on complex geometries remains challenging when boundary- or interface-fitted meshes are difficult to generate, particularly for curved boundaries, polycrystalline microstructures, and dense interface networks. The Shifted Boundary Method (SBM) alleviates this meshing burden by shifting the enforcement of boundary conditions from the true boundary to a nearby surrogate boundary and recovering the effect of the true boundary through geometric correction terms, thereby enabling standard finite element spaces on non-boundary-fitted meshes. In this report, we develop a general shiftedboundary and shifted-interface framework within the open-source MOOSE framework. We first present a general SBM implementation for complex geometries on non-boundary-fitted meshes. We then adopt the Shifted Interface Method (SIM) for internal interfaces and develop a unified shifted-interface treatment in which the interface law is enforced on a surrogate interface and the effect of the true interface is recovered through shifted jumps, fluxes, and tractions. This perspective brings scalar thermal-contact and vector-valued cohesive-zone mechanics into a single framework, the latter realized as the Shifted Cohesive Zone Method (SCZM) and coupled with history-dependent constitutive models from NEML2. We further extend the MOOSE mesh infrastructure to support cohesive-zone calculations on distributed meshes. The framework is verified and demonstrated through three progressive studies: Poisson’s equation on a smoothed starshaped domain, a manufactured thermal-contact problem on a non-interface-fitted mesh, and a two-dimensional polycrystalline representative volume element combining crystal plasticity with cohesive grain-boundary interfaces. Across these studies, the shifted formulations reproduce boundary- and interface-fitted reference solutions with high fidelity, indicating that the proposed framework provides an accurate and efficient route to boundary- and interface-dominated simulations on arbitrary geometries without requiring fitted meshes.

Yang, Cheng-Hau↗
DOE OSTI2024

On plasticity-enhanced interfacial toughness in bonded joints

Here, the performance and reliability of many structures and components depend on the integrity of interfaces between dissimilar materials. Interfacial toughness Γ is the key material parameter that characterizes resistance to interfacial crack growth, and Γ is known to depend on many factors including temperature. For example, previous work showed that the toughness of an epoxy/aluminum interface decreased 40 % as the test temperature was increased from -60 °C to room temperature (RT). Interfacial integrity at elevated temperatures is of considerable practical importance. Recent measurements show that instead of continuing to decrease with increasing temperature, Γ increases when test temperature is above RT. Cohesive zone finite element calculations of an adhesively bonded, asymmetric double cantilever beam specimen of the type used to measure Γ suggest that this increase in toughness may be a result of R-curve behavior generated by plasticity-enhanced toughening during stable subcritical crack growth with interfacial toughness defined as the critical steady-state limit value. In these calculations, which used an elastic-perfectly plastic epoxy model with a temperature-dependent yield strength, the plasticity-enhanced increase in Γ above its intrinsic value Γ o depended on the ratio of interfacial strength σ* to the yield strength σyb of the bond material. There is a nonlinear relationship between Γ/Γ o and σ*/σ yb with the value Γ/Γ o increasing rapidly above a threshold value of σ*/σ yb . The predicted increase in toughness can be significant. For example, there is nearly a factor of two predicted increase in Γ/Γ o during micrometer-scale crack-growth when σ*/σ yb = 2 (a reasonable choice for σ*/σ yb ). Furthermore, contrary to other reported results, plasticity-enhanced toughening can occur prior to crack advance as the cohesive zone forms and the peak stress at the tip of the original crack tip translates to the tip of the fully formed cohesive zone. These results suggest that plasticity-enhanced toughening should be considered when modeling interfaces at elevated temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Transient creep‐fatigue crack growth in creep‐brittle materials: Application to Alloy 718

Abstract This study presents incremental finite element computations of creep‐fatigue crack growth in Alloy 718 at in air. Alloy 718 is representative of creep‐brittle materials, in which viscoplastic deformation is restricted near the tip of a growing crack. The computations predict crack growth using a unique combination of an irreversible cohesive zone formulation and a strain gradient viscoplastic material model based on the Kocks–Mecking formulation. Cohesive zone damage parameters are estimated using sustained loading and constant‐amplitude cyclic loading experiments. Computations of crack extension under three different waveforms containing overloads all predict post‐overload retardation. The amount of retardation depends strongly on the overload ratio, consistent with experiments in the literature using similar waveforms. Analysis of the crack‐tip fields demonstrates retardation is associated with unloading in the highly deformed material near the advancing crack tip. Dynamic recovery and geometrically necessary dislocations are shown to significantly influence post‐overload crack extension.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Numerical modeling and experimental validation of low velocity impact of woven GFRP/CFRP composites

Low-velocity impact of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) composite laminates was studied experimentally and numerically. Hybrid laminates containing blocked layers of GFRP/CFRP/GFRP with all plies oriented at 0° were investigated. Relatively high impact energies were used to obtain full perforation of the laminate in a low-velocity impact setup. Numerical simulations were carried out using the in-house transient dynamics finite element code, Sierra/SM, developed at Sandia National Laboratories. A three-dimensional continuum damage model was used to describe the response of a woven composite ply. Two methods for handling delamination were considered and compared: (1) cohesive zone modeling and (2) continuum damage mechanics. The reduced model size achieved by omission of the cohesive zone elements produced acceptable results at reduced computational cost. Further, the comparison between different modeling techniques can be used to inform modeling decisions relevant to low velocity impact scenarios. The modeling was validated by comparing with the experimental results and showed good agreement in terms of predicted damage mechanisms and impactor velocity and force histories.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fragmentation analysis of a bar with the Lip-field approach

The Lip-field approach was introduced in Moës and Chevaugeon (2021) as a new way to regularize softening material models. It was tested in 1D quasistatic in Moës and Chevaugeon (2021) and 2D quasistatic in Chevaugeon and Moës (2021): this paper extends it to 1D dynamics, on the challenging problem of dynamic fragmentation. The Lip-field approach formulates the mechanical problem to be solved as an optimization problem, where the incremental potential to be minimized is the non-regularized one. Spurious localization is prevented by imposing a Lipschitz constraint on the damage field. Here, the displacement and damage field at each time step are obtained by a staggered algorithm, that is the displacement field is computed for a fixed damage field, then the damage field is computed for a fixed displacement field. Indeed, these two problems are convex, which is not the case of the global problem where the displacement and damage fields are sought at the same time. The incremental potential is obtained by equivalence with a cohesive zone model, which makes material parameters calibration simple. A non-regularized local damage equivalent to a cohesive zone model is also proposed. It is used as a reference for the Lip-field approach, without the need to implement displacement jumps. These approaches are applied to the brittle fragmentation of a 1D bar with randomly perturbed material properties to accelerate spatial convergence. Both explicit and implicit dynamic implementations are compared. Favorable comparison to several analytical, numerical and experimental references serves to validate the modeling approach.

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