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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS2017

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↗

NASA NTRSDate not supplied

Impact of Thermoplastic Composites: Testing and Modeling

Thermoplastic Composites (TPCs) are being increasingly considered for aerospace applications given their faster manufacturing cycles and lower cost. The NASA High-Rate Composite Aircraft Manufacturing (HiCAM) project aims to evaluate and mature thermoplastic manufacturing technologies to achieve a 4X – 6X increase in the production rate of Next Generation Single Aisle (NGSA) commercial aircraft. Within the purview of the NASA HiCAM project, this work investigates the response of thermoplastic composites under High Energy Dynamic Impact (HEDI) conditions. HEDI tests were conducted on panels fabricated with a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK) material system. A blunt metal projectile was chosen as the impactor which impacted the T700/LMPAEK panels over a range of impact velocities imparted by a gas gun test setup. As expected, lower velocity impacts caused the projectile to rebound, whereas higher velocity impacts resulted in the projectile penetrating the test panels. For the test cases in which the projectile rebounded, the predominant damage modes recorded with Ultrasonic (UT) scans included interlaminar delaminations which exhibited a “rotating fan” like structure when viewed in the through-thickness direction. Additionally, the dynamic deflection of the center point of the back face of the test panels and the velocity of the projectile was also measured during the test. Finite element models were developed to model a test case in which the projectile rebounded. These models were developed using LS-DYNA® wherein the primary objective was to evaluate the ability of the material model MAT299 to capture the dynamic response and damage modes in the T700/LMPAEK test panels. MAT299 is a Deformation Gradient Decomposition (DGD) based Continuum Damage Mechanics (CDM) material model. The interlaminar delaminations were modeled using the cohesive contact formulation available within LS-DYNA, wherein the mixed-mode fracture is captured via the Benzegaggh-Kenane (B-K) mode-dependent fracture energy interpolation law. The above-mentioned modeling methods have mostly been applied to thermoset composites and adapting them to TPCs involves addressing challenges associated with appropriately representing material behavior. Therefore, the proposed paper would discuss the systematic approach taken to adapt thermoset composites modeling practices for applications to TPCs, while appropriately addressing material behavior. The proposed paper would show that the predicted back face center point deflection correlated reasonably well with experimental data and the peak deflection was predicted to be within 5% of the experimental measurements. The projectile rebound velocity, while predicted to be higher than the experimental measurement, was within reasonable bounds. Additionally, it will be shown that the predicted delamination shapes were in reasonable agreement with experimental data.

Composite Materials↗

NASA NTRS2014

Challenges in Modelling of Lightning-Induced Delamination; Effect of Temperature-Dependent Interfacial Properties

Lightning is a major cause of damage in laminated composite aerospace structures during flight. Due to the dielectric nature of Carbon fiber reinforced polymers (CFRPs), the high energy induced by lightning strike transforms into extreme, localized surface temperature accompanied with a high-pressure shockwave resulting in extensive damage. It is crucial to develop a numerical tool capable of predicting the damage induced from a lightning strike to supplement extremely expensive lightning experiments. Delamination is one of the most significant failure modes resulting from a lightning strike. It can be extended well beyond the visible damage zone, and requires sophisticated techniques and equipment to detect. A popular technique used to model delamination is the cohesive zone approach. Since the loading induced from a lightning strike event is assumed to consist of extreme localized heating, the cohesive zone formulation should additionally account for temperature effects. However, the sensitivity to this dependency remains unknown. Therefore, the major focus point of this work is to investigate the importance of this dependency via defining various temperature dependency profiles for the cohesive zone properties, and analyzing the corresponding delamination area. Thus, a detailed numerical model consisting of multidirectional composite plies with temperature-dependent cohesive elements in between is subjected to lightning (excessive amount of heat and pressure) and delamination/damage expansion is studied under specified conditions.

interlaminar damage↗

NASA NTRSDate not supplied

Modeling of Response of Thermoplastic Composites Under High Energy Impact

This work investigates the response of thermoplastic composites under impact loading. Panels fabricated with the material system under consideration – carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK) – were subject to impact with an Aluminum projectile, at low and high velocities. At lower velocities, the projectile rebounded from the panel after impact and caused damage e.g. delaminations. Finite element models were developed to model this test using LS-DYNA®, with material model MAT299, a continuum damage model based on the LaRC04 criteria, to capture the dynamic response and damage modes in the T700/LMPAEK panels. Interlaminar delaminations were modeled using the cohesive contact formulation available within LS-DYNA. Previous work by the authors showed that delaminations were consistently underpredicted in size, relative to the test. It was surmised that this was due to effects of strain rate on intralaminar and interlaminar response of the T700/LMPAEK material system. Therefore, the current work is focused on incorporating strain rate effects in modeling the interlaminar response. Rate dependence was incorporated using ‘one-way tie-break contact’ in LS-DYNA, wherein, the Mode-II fracture toughness of the LMPAEK resin, G IIc , was assumed to be an exponential function of the strain rate. In particular, G IIc , was modeled to exponentially decrease at higher strain rates. Models with this form of rate dependence favorably correlated with experimental data regarding delamination extent and rebound velocity, relative to prior models without strain rate effects.

Composite Materials↗

NASA NTRS2002

Nanostructured Composites: Effective Mechanical Property Determination of Nanotube Bundles

Carbon nanotubes naturally tend to form crystals in the form of hexagonally packed bundles or ropes that should exhibit a transversely isotropic constitutive behavior. Although the intratube axial stiffness is on the order of 1 TPa due to a strong network of delocalized bonds, the intertube cohesive strength is orders of magnitude less controlled by weak, nonbonding van der Waals interactions. An accurate determination of the effective mechanical properties of nanotube bundles is important to assess potential structural applications such as reinforcement in future composite material systems. A direct method for calculating effective material constants is developed in the present study. The Lennard-Jones potential is used to model the nonbonding cohesive forces. A complete set of transverse moduli are obtained and compared with existing data.

Saether, E.↗

NASA NTRSDate not supplied

Scalability of Cohesive Fatigue Analyses Using Explicit Solvers

A cohesive fatigue law has been integrated into a constitutive material model compatible with an explicit finite element solver. The cohesive fatigue model response is based on engineering approximations of the endurance limit and the Goodman diagram. This approach can predict stress-life diagrams for crack initiation, the Paris law regime, and transient effects of crack initiation and stable tearing. Simplified cyclic loading is utilized so that the applied load(or displacement) corresponds to the peak load of a fatigue cycle. Loads are held constant during fatigue while damage develops with increasing solution increments. An automatically-calculated ratio of fatigue cycles per solution increment controls the rate of damage growth, ensuring that damage growth is modeled with a sufficient minimum number of increments and damage growth advances to a minimum desired extent within the explicit analysis step time. The compatibility with an explicit finite element solver enables the analysis of structures that are computationally intractable for implicit finite element solvers. Scalability studies are conducted for geometrically nonlinear problems involving fiber-reinforced composite structures that exhibit fatigue damage growth of interacting matrix cracks and delaminations

Frank A Leone↗

DOE OSTI2024

Deformation and Delamination of Dynamically Bulging Bilayered Films

We employ finite element simulations to study the dynamic bulging of films under shock tube conditions and compare these responses with quasi-static bulging. Both scenarios subject the film to time-varying pressure. Two film types are considered: uniform material films and bilayers with a relatively soft interior layer and a hard exterior layer. The interface between the materials in the bilayer cases is modeled as a cohesive zone that follows a bilinear normal traction-separation cohesive law. Additionally, bilayer films with preexisting cracks in the interior layer are analyzed. Parametric studies are also conducted to explore the effects of loading rate and specimen size. The mechanical response of all materials is assumed to be strain-rate independent, focusing solely on the inertial effects in the response of the films under dynamic loading, and the absence of such effects under quasi-static loading. The results indicate that quasi-static bulging exhibits a predictable plate bending like deformation. In contrast, depending on the loading rate and specimen size, dynamic bulging can give rise to elastic waves, and a sequence of deformation processes including initial uniform acceleration, followed by bending-like deformation, and finally through-thickness biaxial stretching. Dynamic loading also leads to larger bulge formation and greater deformation than quasi-static loading. A key finding is the role of preexisting cracks in driving delamination under dynamic loading, which is absent in quasi-static loading. This work enhances our understanding of dynamic bulging and motivates further research on characterizing delamination in layered films using dynamic bulge tests.

36 MATERIALS SCIENCE↗

NASA NTRS2010

Computer Models Simulate Fine Particle Dispersion

Through a NASA Seed Fund partnership with DEM Solutions Inc., of Lebanon, New Hampshire, scientists at Kennedy Space Center refined existing software to study the electrostatic phenomena of granular and bulk materials as they apply to planetary surfaces. The software, EDEM, allows users to import particles and obtain accurate representations of their shapes for modeling purposes, such as simulating bulk solids behavior, and was enhanced to be able to more accurately model fine, abrasive, cohesive particles. These new EDEM capabilities can be applied in many industries unrelated to space exploration and have been adopted by several prominent U.S. companies, including John Deere, Pfizer, and Procter & Gamble.

Source record↗

NASA NTRS2012

Structural Integrity Testing Method for PRSEUS Rod-Wrap Stringer Design

NASA Langley Research Center and The Boeing Company are developing an innovative composite structural concept, called PRSEUS, for the flat center section of a future environmentally friendly hybrid wing body (HWB) aircraft. The PRSEUS (Pultruded Rod Stitched Efficient Unitized Structure) concept uses dry textile preforms for the skins, frames, and stiffener webs. The highly loaded stiffeners are made from precured unidirectional carbon/epoxy rods and dry fiber preforms. The rods are wrapped with the dry fiber preforms and a resin infusion process is used to form the rod-wrap stiffeners. The structural integrity of the rod-wrap interface is critical for maintaining the panel s high strength and bending rigidity. No standard testing method exists for testing the strength of the rod-wrap bondline. Recently, Boeing proposed a rod push-out testing method and conducted some preliminary tests using this method. This paper details an analytical study of the rod-wrap bondline. The rod-wrap interface is modeled as a cohesive zone for studying the initiation and growth of interfacial debonding during push-out testing. Based on the correlations of analysis results and Boeing s test data, the adequacy of the rod-wrap testing method is evaluated, and potential approaches for improvement of the test method are proposed.

Wang, John T.↗

NASA NTRSDate not supplied

Analysis of Crack Propagation Transients in Accelerated Fatigue Characterization of Adhesives

The characterization of the rates of crack propagation by means of a Paris law is intended for steady-state fatigue conditions. While rates of crack propagation are primarily driven by 𝐽𝐽𝑚𝑚𝑚𝑚𝑚𝑚, in the presence of evolving process zones and fiber bridging they also depend on the crack propagation length: cracks with fully formed process zones propagate more slowly. Analyses of the rates of crack propagation in laminates secondarily bonded with AF163-2K film adhesive were conducted to identify the rate transients that occur during the formation of the process zone in the adhesive and during the evolution of bridging of the adhesive’s carrier knit across the crack faces. Dávila's CF22 cohesive fatigue damage model was used in finite element analyses of the double cantilever beam test to explain discrepancies in experimental rates of propagation obtained from tests conducted with either increasing or decreasing 𝐽max. The results indicate that the analyses accurately reproduce transient and steady-state rates of propagation observed in the experimental results, as well as the threshold of propagation. Transient rates of propagation always exceed those in steady state, which underscores the need to understand transients and account for them in fatigue characterization. The results also highlight that a Paris curve for steady-state propagation can be obtained most quickly using a decreasing load function, but that this curve consists of an unconservative bound to the rates of propagation. A lateral shift in the Paris curve is proposed for a conservative characterization that removes the effect of bridging and other R-curve effects.

fracture mechanics; delamination; J-integral↗

NASA NTRSDate not supplied

Analysis of Crack Propagation Transients in Accelerated Fatigue Characterization of Adhesives

fracture mechanics; delamination; J-integral↗

NASA NTRS2012

Getting a Cohesive Answer from a Common Start: Scalable Multidisciplinary Analysis through Transformation of a System Model

One of the challenges of systems engineering is in working multidisciplinary problems in a cohesive manner. When planning analysis of these problems, system engineers must tradeoff time and cost for analysis quality and quantity. The quality is associated with the fidelity of the multidisciplinary models and the quantity is associated with the design space that can be analyzed. The tradeoff is due to the resource intensive process of creating a cohesive multidisciplinary system model and analysis. Furthermore, reuse or extension of the models used in one stage of a product life cycle for another is a major challenge. Recent developments have enabled a much less resource-intensive and more rigorous approach than handwritten translation scripts or codes of multidisciplinary models and their analyses. The key is to work from a core system model defined in a MOF-based language such as SysML and in leveraging the emerging tool ecosystem, such as Query-View- Transform (QVT), from the OMG community. SysML was designed to model multidisciplinary systems and analyses. The QVT standard was designed to transform SysML models. The Europa Hability Mission (EHM) team has begun to exploit these capabilities. In one case, a Matlab/Simulink model is generated on the fly from a system description for power analysis written in SysML. In a more general case, a symbolic mathematical framework (supported by Wolfram Mathematica) is coordinated by data objects transformed from the system model, enabling extremely flexible and powerful tradespace exploration and analytical investigations of expected system performance.

Cole, Bjorn↗

NASA NTRS2012

Getting a Cohesive Answer from a Common Start: Scalable Multidisciplinary Analysis through Transformation of a Systems Model

One of the challenges of systems engineering is in working multidisciplinary problems in a cohesive manner. When planning analysis of these problems, system engineers must trade between time and cost for analysis quality and quantity. The quality often correlates with greater run time in multidisciplinary models and the quantity is associated with the number of alternatives that can be analyzed. The trade-off is due to the resource intensive process of creating a cohesive multidisciplinary systems model and analysis. Furthermore, reuse or extension of the models used in one stage of a product life cycle for another is a major challenge. Recent developments have enabled a much less resource-intensive and more rigorous approach than hand-written translation scripts between multi-disciplinary models and their analyses. The key is to work from a core systems model defined in a MOF-based language such as SysML and in leveraging the emerging tool ecosystem, such as Query/View/Transformation (QVT), from the OMG community. SysML was designed to model multidisciplinary systems. The QVT standard was designed to transform SysML models into other models, including those leveraged by engineering analyses. The Europa Habitability Mission (EHM) team has begun to exploit these capabilities. In one case, a Matlab/Simulink model is generated on the fly from a system description for power analysis written in SysML. In a more general case, symbolic analysis (supported by Wolfram Mathematica) is coordinated by data objects transformed from the systems model, enabling extremely flexible and powerful design exploration and analytical investigations of expected system performance.

Operational QVT↗

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↗

NASA NTRS2007

A Selection of Composites Simulation Practices at NASA Langley Research Center

One of the major areas of study at NASA Langley Research Center is the development of technologies that support the use of advanced composite materials in aerospace applications. Amongst the supporting technologies are analysis tools used to simulate the behavior of these materials. This presentation will discuss a number of examples of analysis tools and simulation practices conducted at NASA Langley. The presentation will include examples of damage tolerance analyses for both interlaminar and intralaminar failure modes. Tools for modeling interlaminar failure modes include fracture mechanics and cohesive methods, whilst tools for modeling intralaminar failure involve the development of various progressive failure analyses. Other examples of analyses developed at NASA Langley include a thermo-mechanical model of an orthotropic material and the simulation of delamination growth in z-pin reinforced laminates.

Ratcliffe, James G.↗

NASA NTRS2017

Simulating Matrix Crack and Delamination Interaction in a Clamped Tapered Beam

Blind predictions were conducted to validate a discrete crack methodology based on the Floating Node Method to simulate matrix-crack/delamination interaction. The main novel aspects of the approach are: (1) the implementation of the floating node method via an 'extended interface element' to represent delaminations, matrix-cracks and their interaction, (2) application of directional cohesive elements to infer overall delamination direction, and (3) use of delamination direction and stress state at the delamination front to determine migration onset. Overall, good agreement was obtained between simulations and experiments. However, the validation exercise revealed the strong dependence of the simulation of matrix-crack/delamination interaction on the strength data (in this case transverse interlaminar strength, YT) used within the cohesive zone approach applied in this work. This strength value, YT, is itself dependent on the test geometry from which the strength measurement is taken. Thus, choosing an appropriate strength value becomes an ad-hoc step. As a consequence, further work is needed to adequately characterize and assess the accuracy and adequacy of cohesive zone approaches to model small crack growth and crack onset. Additionally, often when simulating damage progression with cohesive zone elements, the strength is lowered while keeping the fracture toughness constant to enable the use of coarser meshes. Results from the present study suggest that this approach is not recommended for any problem involving crack initiation, small crack growth or multiple crack interaction.

De Carvalho, N. V.↗

NASA NTRS1989

Nonlinear fracture of concrete and ceramics

The nonlinear fracture process zones in an impacted unnotched concrete bend specimen, a prenotched ceramic bend specimen, and an unnotched ceramic/ceramic composite bend specimen were estimated through hybrid experimental numerical analysis. Aggregate bridging in concrete, particulate bridging in ceramics, and fiber bridging in ceramic/ceramic composite are modeled by Barenblatt-type cohesive zones which are incorporated into the finite-element models of the bend specimens. Both generation and propagation analyses are used to estimate the distribution of crack closure stresses in the nonlinear fracture process zones. The finite-element models are then used to simulate fracture tests consisting of rapid crack propagation in an impacted concrete bend specimen, and stable crack growth and strain softening in a ceramic and ceramic/ceramic composite bend specimens.

Kobayashi, Albert S.↗

NASA NTRS2010

Finite Element Model for Failure Study of Two-Dimensional Triaxially Braided Composite

A new three-dimensional finite element model of two-dimensional triaxially braided composites is presented in this paper. This meso-scale modeling technique is used to examine and predict the deformation and damage observed in tests of straight sided specimens. A unit cell based approach is used to take into account the braiding architecture as well as the mechanical properties of the fiber tows, the matrix and the fiber tow-matrix interface. A 0 deg / plus or minus 60 deg. braiding configuration has been investigated by conducting static finite element analyses. Failure initiation and progressive degradation has been simulated in the fiber tows by use of the Hashin failure criteria and a damage evolution law. The fiber tow-matrix interface was modeled by using a cohesive zone approach to capture any fiber-matrix debonding. By comparing the analytical results to those obtained experimentally, the applicability of the developed model was assessed and the failure process was investigated.

Li, Xuetao↗