Search NASA⌕ Search

SEARCH · Search NASA

Results for “cohesive models”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Mesh Convergence Requirements for Composite Damage Models

The ability of the finite element method to accurately represent the response of objects with intricate geometry and loading renders the finite element method as an extremely versatile analysis technique for structural analysis. Finite element analysis is routinely used in industry to calculate deflections, stress concentrations, natural frequencies, buckling loads, and much more. The method works by discretizing complex problems into smaller, simpler approximations that are valid over small uniform domains. For common analyses, the maximum size of the elements that can be used is often be determined by experience. However, to verify the quality of a solution, analyses with several levels of mesh refinement should be performed to ensure that the solution has converged. In recent years, the finite element method has been used to calculate the resistance of structures, and in particular that of composite structures. A number of techniques such as cohesive zone modeling, the virtual crack closure technique, and continuum damage modeling have emerged that can be used to predict cracking, delaminations, fiber failure, and other composite damage modes that lead to structural collapse. However, damage models present mesh refinement requirements that are not well understood. In this presentation, we examine different mesh refinement issues related to the representation of damage in composite materials. Damage process zone sizes and their corresponding mesh requirements will be discussed. The difficulties of modeling discontinuities and the associated need for regularization techniques will be illustrated, and some unexpected element size constraints will be presented. Finally, some of the difficulties in constructing models of composite structures capable of predicting transverse matrix cracking will be discussed. It will be shown that to predict the initiation and propagation of transverse matrix cracks, their density, and their saturation may require models that are significantly more refined than those that have been contemplated in the past.

Davila, Carlos G.↗
DOE OSTI2025

Interfacial Strength in Hierarchical Carbon Fiber Composites: Interplay of Interphase Modulus and Roughness

A facile, direct deposition approach that exploits van der Waals interactions between carbonaceous materials is utilized to create unidirectional hybrid carbon fiber composites. Two small molecule crosslinkers, a trifunctional aromatic (TL) and a difunctional aliphatic (DL) acyl chloride, are first utilized to create a crosslinked interphase with a softer and stiffer modulus respectively. TL crosslinked interphase with a higher modulus improved the tensile strength by 50%, despite non-covalent linking between fiber and matrix, elucidating the critical role of the interphase in alleviating modulus mismatch between the high modulus carbon fiber and the rubbery matrix. Fractional quantities of carbon nanotubes are additionally dispersed in the small molecule crosslinkers which behaved as a dispersant, helping introduce nanoasperities on the carbon fiber surface. Strong “pi-pi” interactions between CNTs and CF contributed to tensile properties, which are increased by 66% compared to the control. A cohesive zone model suggests that a stiffer interphase is better able to exploit surface heterogeneities and roughness on the fiber, synergistically enhancing interfacial strength.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗
NASA NTRS1993

A thermodynamic analysis of propagating subcritical cracks with cohesive zones

The results of the so-called energetic approach to fracture with particular attention to the issue of energy dissipation due to crack propagation are applied to the case of a crack with cohesive zone. The thermodynamic admissibility of subcritical crack growth (SCG) is discussed together with some hypotheses that lead to the derivation of SCG laws. A two-phase cohesive zone model for discontinuous crack growth is presented and its thermodynamics analyzed, followed by an example of its possible application.

Allen, David H.↗
NASA NTRS2009

The Development of Directional Decohesion Finite Elements for Multiscale Failure Analysis of Metallic Polycrystals

Atomistic simulations of intergranular fracture have indicated that grain-scale crack growth in polycrystalline metals can be direction dependent. At these material length scales, the atomic environment greatly influences the nature of intergranular crack propagation, through either brittle or ductile mechanisms, that are a function of adjacent grain orientation and direction of crack propagation. Methods have been developed to obtain cohesive zone models (CZM) directly from molecular dynamics simulations. These CZMs may be incorporated into decohesion finite element formulations to simulate fracture at larger length scales. A new directional decohesion element is presented that calculates the direction of Mode I opening and incorporates a material criterion for dislocation emission based on the local crystallographic environment to automatically select the CZM that best represents crack growth. The simulation of fracture in 2-D and 3-D aluminum polycrystals is used to illustrate the effect of parameterized CZMs and the effectiveness of directional decohesion finite elements.

Saether, Erik↗
NASA NTRS2014

Adaptive Shape Functions and Internal Mesh Adaptation for Modelling Progressive Failure in Adhesively Bonded Joints

Enhanced finite elements are elements with an embedded analytical solution that can capture detailed local fields, enabling more efficient, mesh independent finite element analysis. The shape functions are determined based on the analytical model rather than prescribed. This method was applied to adhesively bonded joints to model joint behavior with one element through the thickness. This study demonstrates two methods of maintaining the fidelity of such elements during adhesive non-linearity and cracking without increasing the mesh needed for an accurate solution. The first method uses adaptive shape functions, where the shape functions are recalculated at each load step based on the softening of the adhesive. The second method is internal mesh adaption, where cracking of the adhesive within an element is captured by further discretizing the element internally to represent the partially cracked geometry. By keeping mesh adaptations within an element, a finer mesh can be used during the analysis without affecting the global finite element model mesh. Examples are shown which highlight when each method is most effective in reducing the number of elements needed to capture adhesive nonlinearity and cracking. These methods are validated against analogous finite element models utilizing cohesive zone elements.

computational mechanics↗
NASA NTRSDate not supplied

Heat Pipe Heat Exchanger for Nuclear Electric Propulsion Power Conversion System

Heat pipe reactors have been considered by the Space Nuclear Propulsion program for Nuclear Electric Propulsion (NEP) power conversion systems and will require the use of heat exchangers to transfer heat via heat pipes to the Brayton working fluid from the reactor. Sodium (Na) and lithium (Li) were considered as viable working fluids inside the heat pipes which were assumed to have the same geometry based on studies and information from the Los Alamos National Laboratory. The heat exchanger was assumed to be a rectangular duct with heat pipes serving as tubes from previous NEP work and recommendations. Based on this geometry, Zukauskas correlations were used to model the convective heat transfer and pressure losses. Parametric sizing of the reactor component involved operational limits-based heat pipe thermal hydraulic modeling in cohesion with required user input geometry for the in-core lattice and various subcomponents. This work considered various power conversion inlet temperatures (PCIT) of 1100 K, 1150 K, and 1200 K for Na heat pipes and 1100 K, 1150 K, 1200 K, and 1400 K for Li heat pipes based on recommendations from prior work. Using these different PCITs, the subsystem masses and pressure losses were determined and analyzed. Na showed a lower overall operating temperature and about a fifth of the maximum heat throughput capability than that of Li for the same geometry. Due to this, the entire Na-based subsystem ended up being three times more massive than the Li-based subsystem given five times the required number of heat pipes. At the low PCIT of 1100 K, the Na-based subsystem exhibited the lowest pressure losses given the large overall cross sectional flow area and relatively low frictional pressure losses. However, as the PCIT increased, the frictional pressure losses increased resulting in higher pressure losses at the 1200 K PCIT than Li-based subsystem. However, the Li-based subsystem exhibited the largest pressure losses of all analyzed cases at the 1400 K PCIT due to the low density of the Brayton working fluid at this temperature.

electric↗
DOE OSTI2026

Influence of Shear Strength Assumptions on BISON Debonding Simulations

Accurately predicting the thermomechanical response of buffer–IPyC debonding in TRISO fuel particles requires reliable mechanical property inputs for each coating layer, particularly the normal and shear strengths that influence interlayer delamination and stress concentrations. Micro tensile testing of AGR-2 fuel particles provided experimentally measured normal strengths for the buffer, IPyC, and buffer–IPyC interface; however, shear strength was not measured. As a result, BISON simulations of interface debonding must rely on assumed shear strength values, typically estimated as 20–40% of the measured ultimate tensile strength. This study evaluates how these assumed shear strength values influence cohesive zone model (CZM) predictions of buffer–IPyC separation in AGR 2 TRISO particles. Using micro tensile data from three AGR 2 compacts (2 1 3, 5 1 3, and 6 3 3), BISON simulations were performed with multiple shear strength assumptions to quantify their effect on radial and tangential stress evolution, debonding, and gap propagation. The results show that shear strength is a high sensitivity parameter: increasing the assumed shear strength significantly alters the stress distribution at the buffer–IPyC junction, shifts the predicted debonding location, and changes the extent of partial gap formation. While normal strength controls the initiation of interface separation, shear strength strongly influences the mode mixity of the failure process and the resulting stress concentrations transmitted to the IPyC and SiC layers. These findings highlight a critical gap in current TRISO mechanical characterization. Without experimentally measured shear strength, BISON simulations must rely on approximations that introduce uncertainty into predictions of coating layer integrity and fission product barrier performance. Future fuel qualification campaigns should therefore consider measurement of shear strength at the interlayer interfaces to reduce model uncertainty and improve the fidelity of TRISO fuel performance simulations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

A High-Performance Discrete-Element Framework for Simulating Flow and Jamming of Moisture Bearing Biomass Feedstocks

We developed and verified a high-performance open-source discrete element method (DEM) solver with simultaneously-supported feedstock-specific interaction models, including bonded-sphere, liquid bridge, cohesion, and non-linear contact models. Our solver uses parallel data structures on hybrid central and graphics processing unit (CPU/GPU) architectures, with favorable strong scaling performance observed for large problem sizes comprised of (100 M particles), and 4X single-node GPU speedup. The particles for corn stover feedstock were conceptualized and calibrated based on experimental measurements and results. Sensitivity analyses demonstrate that the mass flow rate from a wedge hopper is governed primarily by moisture content, friction coefficient, and cohesion energy density. The model is used to reproduce experimentally observed hopper jamming results, highlighting that the experimental no-flow trends can only be achieved by using non-spherical particles, liquid bridge and cohesion models, highlighting the importance of using concurrent feedstock specialized models for the effective representation of biomass material handling problems.

bioenergy↗
NASA NTRS2013

Model-based Systems Engineering: Creation and Implementation of Model Validation Rules for MOS 2.0

Model-based Systems Engineering (MBSE) is an emerging modeling application that is used to enhance the system development process. MBSE allows for the centralization of project and system information that would otherwise be stored in extraneous locations, yielding better communication, expedited document generation and increased knowledge capture. Based on MBSE concepts and the employment of the Systems Modeling Language (SysML), extremely large and complex systems can be modeled from conceptual design through all system lifecycles. The Operations Revitalization Initiative (OpsRev) seeks to leverage MBSE to modernize the aging Advanced Multi-Mission Operations Systems (AMMOS) into the Mission Operations System 2.0 (MOS 2.0). The MOS 2.0 will be delivered in a series of conceptual and design models and documents built using the modeling tool MagicDraw. To ensure model completeness and cohesiveness, it is imperative that the MOS 2.0 models adhere to the specifications, patterns and profiles of the Mission Service Architecture Framework, thus leading to the use of validation rules. This paper outlines the process by which validation rules are identified, designed, implemented and tested. Ultimately, these rules provide the ability to maintain model correctness and synchronization in a simple, quick and effective manner, thus allowing the continuation of project and system progress.

Model-based Systems Engineering (MBSE)↗
NASA NTRSDate not supplied

Challenges in Scaling Analysis from Coupons to Elements in 3D Woven Composites for a Launch Vehicle End-Ring Structure

The limitations of structural design with 3D woven composites based on homogenized material models characterized using coupon test data, representative unit cell analysis, or a combination of the two are discussed in the context of a space launch system (SLS)-scale end ring. Selected building block test configurations at the coupon, element, and subcomponent scales are described in detail. Highly instrumented tests and X-ray computed tomography inspections are used to characterize the structural behavior and associated damage progression with increasing load. It is shown that a single homogenized material constitutive relationship, applied at coupon and element scale, is not sufficient to model the structural response observed in both tests. In fact, substantially different stress-strain laws are required to predict the response of the coupon and element scale specimens, despite similar loading conditions. A model with discrete representation of the binder tows and cohesive interfaces to model matrix cracking shows much better agreement with the test data at the two length scales. Thus, evidence is provided showing that predictive structural analysis of 3D woven composites requires discrete representation of the tows, matrix cracking morphology, and their interactions in the structural configuration. With accurate knowledge of the element scale behavior, engineering approximations readily enable analysis scale-up to the subcomponent level.

Andrew C Bergan↗
NASA NTRS2018

A Dual Spring Modeling Approach for Static and Fatigue Failure Assessments of Carbon/Epoxy Composite Sub-Elements

A dual spring model is developed for the static and fatigue damage predictions of laminates interface in composite structures. Stress concentrations can be induced by the defects formed in the fabrication or service process. A conventional S-N based fatigue damage model may not be accurate to predict the fatigue life of a structure with high stress concentration. With the dual spring model, static delamination failure can be simulated using springs of cohesive type material model while fatigue delamination development can be predicted using linear springs, where the crack driving force is computed based on virtual crack closure technique (VCCT). A Paris law type fatigue growth law with its mode mixity is applied for fatigue crack growth prediction. After verified using benchmark examples, including Double Cantilever Beam (DCB), End-Notched Flexure (ENF) and Mix-Mode Bending (MMB), the proposed dual spring model is applied in the static and fatigue damage prediction of NASA/Boeing sub-elements and UTC sub-elements.

Cui, Xiaodong↗
NASA NTRS1990

A phenomenological model of coating/substrate adhesion and interfacial bimetallic peeling stress in composite mirrors

Adhesion and interfacial stress between metal films and structural composite material substrates is discussed. A theoretical and conceptual basis for selecting coating materials for composite mirror substrates is described. A phenomenological model that interrelates cohesive tensile strength of thin film coatings and interfacial peeling stresses is presented. The model serves as a basis in determining gradiated materials response and compatibility of composite substrate and coating combinations. Parametric evaluation of material properties and geometrical factors such as coating thickness are used to determine the threshold stress levels for maintaining adhesion at the different interfaces.

Mcelroy, Paul M.↗
NASA NTRS1985

Thermoplastic matrix composite processing model

The effects the processing parameters pressure, temperature, and time have on the quality of continuous graphite fiber reinforced thermoplastic matrix composites were quantitatively accessed by defining the extent to which intimate contact and bond formation has occurred at successive ply interfaces. Two models are presented predicting the extents to which the ply interfaces have achieved intimate contact and cohesive strength. The models are based on experimental observation of compression molded laminates and neat resin conditions, respectively. Identified as the mechanism explaining the phenomenon by which the plies bond to themselves is the theory of autohesion (or self diffusion). Theoretical predictions from the Reptation Theory between autohesive strength and contact time are used to explain the effects of the processing parameters on the observed experimental strengths. The application of a time-temperature relationship for autohesive strength predictions is evaluated. A viscoelastic compression molding model of a tow was developed to explain the phenomenon by which the prepreg ply interfaces develop intimate contact.

Dara, P. H.↗
NASA NTRS2020

Evaluation of Fatigue Damage Accumulation Functions for Delamination Initiation and Propagation

The present report follows on the cohesive fatigue damage model methodology proposed in NASA-TP-2018-219838. In that report, an empirical function describing the incremental damage due to cyclic loading was used to calculate fatigue damage within a cohesive formulation. The form of the function was developed such that, when integrated at a constant stress amplitude from no damage to failure, it produces a life versus load response that is consistent with an S-N diagram. Therefore, the parameters of the damage model could be obtained by fitting the model predictions to an S-N diagram. The finite element analyses performed demonstrate that the cohesive fatigue accumulation function provides a link between the S-N diagram that describes crack initiation, and the Paris law that characterizes the rate of crack propagation. However, when the model was proposed, it was not known whether the form of the damage accumulation function associated with a desired S-N diagram is unique and, if not, if the link between S-N and the Paris law is unique and independent of the fatigue function selected. In the effort described herein, several alternative forms of the damage function that reproduce the desired features of S-N diagrams were found and evaluated. The effects of each of these functions on the predicted parameters of the Paris law and the propagation threshold are discussed. The results indicate that the predicted exponent m of the Paris law is indeed independent of the damage accumulation function. However, different functions predict different values for the pre-factor C of the Paris law. Therefore, the proper damage accumulation function must be selected by comparison with experiments. One of the new damage accumulation functions proposed herein was found to be particularly useful for analysis because of the ease with which the model parameters can be determined with a minimal amount of experimental information. The effectiveness of the proposed methodology and damage function was demonstrated by conducting analyses of a double cantilever beam test, a mixed-mode bending test, and a three-point bending test of a skin/doubler specimen. The results indicate that the same set of model parameters can provide accurate predictions of the rate of fatigue crack propagation for a variety of material interfaces, mode mixities, load levels, and stress ratios.

Cohesive elements↗
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↗