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 55 records · Page 3

Multiscale Modeling of Damage Processes in Aluminum Alloys: Grain-Scale Mechanisms

This paper has two goals related to the development of a physically-grounded methodology for modeling the initial stages of fatigue crack growth in an aluminum alloy. The aluminum alloy, AA 7075-T651, is susceptible to fatigue cracking that nucleates from cracked second phase iron-bearing particles. Thus, the first goal of the paper is to validate an existing framework for the prediction of the conditions under which the particles crack. The observed statistics of particle cracking (defined as incubation for this alloy) must be accurately predicted to simulate the stochastic nature of microstructurally small fatigue crack (MSFC) formation. Also, only by simulating incubation of damage in a statistically accurate manner can subsequent stages of crack growth be accurately predicted. To maintain fidelity and computational efficiency, a filtering procedure was developed to eliminate particles that were unlikely to crack. The particle filter considers the distributions of particle sizes and shapes, grain texture, and the configuration of the surrounding grains. This filter helps substantially reduce the number of particles that need to be included in the microstructural models and forms the basis of the future work on the subsequent stages of MSFC, crack nucleation and microstructurally small crack propagation. A physics-based approach to simulating fracture should ultimately begin at nanometer length scale, in which atomistic simulation is used to predict the fundamental damage mechanisms of MSFC. These mechanisms include dislocation formation and interaction, interstitial void formation, and atomic diffusion. However, atomistic simulations quickly become computationally intractable as the system size increases, especially when directly linking to the already large microstructural models. Therefore, the second goal of this paper is to propose a method that will incorporate atomistic simulation and small-scale experimental characterization into the existing multiscale framework. At the microscale, the nanoscale mechanics are represented within cohesive zones where appropriate, i.e. where the mechanics observed at the nanoscale can be represented as occurring on a plane such as at grain boundaries or slip planes at a crack front. Important advancements that are yet to be made include: 1. an increased fidelity in cohesive zone modeling; 2. a means to understand how atomistic simulation scales with time; 3. a new experimental methodology for generating empirical models for CZMs and emerging materials; and 4. a validation of simulations of the damage processes at the nano-micro scale. With ever-increasing computer power, the long-term ability to employ atomistic simulation for the prognosis of structural components will not be limited by computation power, but by our lack of knowledge in incorporating atomistic models into simulations of MSFC into a multiscale framework.

Hochhalter, J. D.↗

New Developments in the Embedded Statistical Coupling Method: Atomistic/Continuum Crack Propagation

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain has been enhanced. The concurrent MD-FEM coupling methodology uses statistical averaging of the deformation of the atomistic MD domain to provide interface displacement boundary conditions to the surrounding continuum FEM region, which, in turn, generates interface reaction forces that are applied as piecewise constant traction boundary conditions to the MD domain. The enhancement is based on the addition of molecular dynamics-based cohesive zone model (CZM) elements near the MD-FEM interface. The CZM elements are a continuum interpretation of the traction-displacement relationships taken from MD simulations using Cohesive Zone Volume Elements (CZVE). The addition of CZM elements to the concurrent MD-FEM analysis provides a consistent set of atomistically-based cohesive properties within the finite element region near the growing crack. Another set of CZVEs are then used to extract revised CZM relationships from the enhanced embedded statistical coupling method (ESCM) simulation of an edge crack under uniaxial loading.

Saether, E.↗

Fracture-Based Mesh Size Requirements for Matrix Cracks in Continuum Damage Mechanics Models

This paper evaluates the ability of progressive damage analysis (PDA) finite element (FE) models to predict transverse matrix cracks in unidirectional composites. The results of the analyses are compared to closed-form linear elastic fracture mechanics (LEFM) solutions. Matrix cracks in fiber-reinforced composite materials subjected to mode I and mode II loading are studied using continuum damage mechanics and zero-thickness cohesive zone modeling approaches. The FE models used in this study are built parametrically so as to investigate several model input variables and the limits associated with matching the upper-bound LEFM solutions. Specifically, the sensitivity of the PDA FE model results to changes in strength and element size are investigated.

Leone, Frank A.↗

Microwave Signatures of Melting/Refreezing Snow: Observations and Modeling Using Dense Medium Radiative Transfer Theory

Microwave brightness temperatures of snow covered terrains can be modeled by means of the Dense Radiative Transfer Medium Theory (DMRT). In a dense medium, such as snow, the assumption of independent scattering is no longer valid and the scattering of correlated scatterers must be considered. In the DMRT, this is done considering a pair distribution function of the particles position. In the electromagnetic model, the snowpack is simulated as a homogeneous layer having effective permittivity and albedo calculated through the DMRT. In order to account for clustering of snow crystals, a model of cohesive particles can be applied, where the cohesion between the particles is described by means of a dimensionless parameters called stickiness (z), representing a measure of the inversion of the attraction of the particles. The lower the z the higher the stickiness. In this study, microwave signatures of melting and refreezing cycles of seasonal snowpacks at high altitudes are studied by means of both experimental and modeling tools. Radiometric data were collected 24 hours per day by the University of Michigan Tower Mounted Radiometer System (TMRS). The brightness temperatures collected by means of the TMRS are simulated by means of a multi-layer electromagnetic model based on the dense medium theory with the inputs to the model derived from the data collected at the snow pits and from the meteorological station. The paper is structured as follows: in the first Section the temperature profiles recorded by the meteorological station and the snow pit data are presented and analyzed; in the second Section, the characteristics of the radiometric system used to collect the brightness temperatures are reported together with the temporal behavior of the recorded brightness temperatures; in the successive Section the multi-layer DMRT-based electromagnetic model is described; in the fourth Section the comparison between modeled and measured brightness temperatures is discussed. We dedicate the last Section to the conclusions and future works.

Tedesco, Marco↗

Modeling of Response of Thermoplastic Composites Under High Energy Impact

This work investigates the response of a thermoplastic panel under high energy ballistic impact. The composite material system is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK3). Several high energy dynamic impact tests were performed where an aluminum projectile impacted various panels at different impact velocities. For this work, the test performed at an impact velocity of 73.2 m/s was investigated. The projectile rebounded upon impact, and the damage to the panel was characterized through ultrasonic and CT-scanning. Response of the panel under impact and the associated damage is discussed. Numerical simulation of the impact process was performed using the LS-DYNA® computer code [1] with a new material model MAT_COMPOSITE_DAMAGE_DGD (MAT 299), a modification over MAT 261 material model in terms of an improved representation of damage coupling and additional capability for strain rate effects. Delaminations were modeled using cohesive zone contact formulation in LS-DYNA®, wherein the mixed-mode fracture is captured via the Benzegaggh-Kenane (B-K) mode-dependent fracture energy interpolation law. Simulation of this impact test demonstrates a good match of results in terms of panel deflection and panel-projectile interactions. The overall extent of delaminations was under-predicted by the simulation. This work forms the basis to model tests performed at higher velocities, involving projectile penetration of panel.

Composite Materials↗

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.↗

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.↗

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↗

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↗

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↗

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)↗

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↗

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↗

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.↗

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.↗

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↗

Interaction of Delaminations and Matrix Cracks in a CFRP Plate, Part II: Simulation Using an Enriched Shell Finite Element Model

Numerical simulations are presented of a recently developed test which creates multiple delaminations in a CFRP laminate specimen that grow and interact via transverse matrix cracks. A novel shell element enriched with the Floating Node Method, and a damage algorithm based on the Virtual Crack Closure Technique, were used to successfully simulate the tests. Additionally, a 3D high mesh fidelity model based on cohesive zones and continuum damage mechanics was used to simulate the tests and act as a representative of other similar state-of-the-art high mesh fidelity modeling techniques to compare to the enriched shell element. The enriched shell and high mesh fidelity models had similar levels of accuracy and generally matched the experimental data. With runtimes of 36 minutes for the shell model and 55 hours for the high mesh fidelity model, the shell model is 92 times faster than the high- fidelity simulation.

Mark W McElroy↗