Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interface properties”

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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

silicone intumescence↗
DOE OSTI2025

Entrapment Behavior of Solid Surrogate Fission Products at Engineered UN Nano‐Hetero‐Interfaces Within Metallic Nuclear Fuels

Nanometric hetero-interfaces provide a wealth of scientific and engineering opportunities due to their complex and often misunderstood properties that can differ from their respective bulk constituents. In this work, the ability for engineered nanostructures within a bulk U─Mo alloy to arrest simulant fission products is investigated experimentally and computationally. Nanostructured 90 wt% U/ 10 wt% Mo (U-10Mo) with 7.1 at% Nd is consolidated using spark-plasma- sintering (SPS) techniques and is heat-treated at 500 °C under vacuum for 24, 100, 500, and 1000 h. Analysis on the sintered and heat-treated U-10Mo reveals rapid kinetics in Nd diffusion to nanocluster sites, with evidence of Nd diffusion occurring during sintering and during the following heat-treatment. The segregation behavior of Nd at two different U─Mo/UN interfaces is computationally verified using density functional theory (DFT) to reinforce experimental data.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1987

Comparison of flexural properties of aramid-reinforced pultrusions having varied matrices, pretreatments and postcures

Aramid-reinforced composite materials of equal fiber volume and varied polymer thermoset matrices were pultruded and flexurally tested to failure. The objective was to improve the flexural properties of aramid-reinforced pultrusions. Pultrusions of both sized and unsized aramid fiber with four different resin systems were compared to determine the effects of sizing compounds and postcuring on flexural strength, fiber wettability, and fiber-to-resin interface bonding. Improvements in flexural strength resulting from pretreatments with the sizing solutions used were marginal. The most significant improvements in flexural properties resulted from postcuring. Flexural strengths ranged from a low of 39,647 psi (273MPa) to a high of 80,390 psi (554 MPa), an overall increase of 103 percent. The fact that postcuring improved the flexural properties of the pultrusions of the four resin systems indicates that a full cure did not occur in any of the resin systems during the pultrusion process. The increased flexural strengths of the polyester and vinyl ester pultrusions were the most surprising. The four resin systems examined were Interplastic Corporation VE 8300 vinyl ester, Ashland Chemical Company Aropol 7430 Polyester, and Shell Chemical Company Epon 9302 and Epon 9310 epoxides.

Wilson, Maywood L.↗
NASA NTRS1993

An interface configuration experiment on USML-1

Experiments were carried out for 'exotic' rotationally symmetric containers aboard the first NASA United States Microgravity Laboratory Space Shuttle flight (USML-1). The containers have the property that they admit an entire continuum of distinct equilibrium rotationally-symmetric capillary free-surfaces for a given liquid volume and contact angle. It was found, after the containers were filled in orbit, that an initial equilibrium interface from the symmetric continuum reoriented, when perturbed, to a stable interface that was not rotationally symmetric, as predicted by the mathematical theory.

Concus, Paul↗
DOE OSTI2022

Beyond Local Solvation Structure: Nanometric Aggregates in Battery Electrolytes and their Effect on Electrolyte Properties

Electrolytes are an essential component of all electrochemical storage and conversion devices, such as batteries. In the history of battery development, the complex nature of electrolytes has often been a bottleneck. Fundamental knowledge of electrolyte systems encompasses elucidation of structure-property relationships of the solution species. Recently, nanometric aggregates have been observed in several classes of electrolytes, including super-concentrated, redox-flow, multivalent, polymer, and ionic liquid-based electrolytes. Compared with the well-studied local solvation structures such as contact ion pairs and solvent-separated ions, these aggregates impose unique effects on the ion distribution and transport both within bulk electrolytes and at electrode/electrolyte interfaces. This Perspective highlights the discovery of the aggregates in various battery electrolytes and their impact on electrolyte properties. We also present an outlook for future studies of this emerging field of nanometric aggregates and the need for the development of new experimental and computational tools to study their properties.

Yu, Zhou↗
DOE OSTI2023

Distance-dependent dielectric constant at the calcite/electrolyte interface: Implication for surface complexation modeling

Hypothesis: The electrical double layer formed at the mineral/electrolyte interface is often modeled using mean-field approaches based on a continuum description of the solvent whose dielectric constant is assumed to decrease monotonically with decreasing distance to the surface. In contrast, molecular simulations show that the solvent polarizability oscillates near the surface similar to the water density profile - as shown previously, for example, by Bonthuis et al. (D.J. Bonthuis, S. Gekle, R.R. Netz, Dielectric Profile of Interfacial Water and its Effect on Double-Layer Capacitance, Phys Rev Lett 107(16) (2011) 166102). Here, we showed that molecular and mesoscale pictures agree by spatially averaging the dielectric constant obtained from molecular dynamics simulations over the distances relevant to the mean-field representation. In addition, the values of capacitances used to describe the electrical double layer in Surface Complexation Models (SCMs) of the mineral/electrolyte interface can be estimated using molecularly informed spatially averaged dielectric constants and positions of hydration layers. Experiments: First, we used molecular dynamics simulations to model the calcite 101¯4/electrolyte interface. Next, by using atomistic trajectories, we calculated the distance-dependent static dielectric constant and water density in the direction normal to the. Finally, we applied spatial compartmentalization consistent with the model of parallel-plate capacitors connected in series to estimate SCM capacitances. Findings: Computationally expensive simulations are required to determine the dielectric constant profile of interfacial water near the mineral surface. On the other hand, water density profiles are readily assessable from much shorter simulation trajectories. Our simulations confirmed that dielectric and water density oscillations at the interface are correlated. Here, we parametrized linear regression models to estimate the dielectric constant directly from the local water density. This is a significant computational shortcut compared to slowly converging calculations relying on total dipole moment fluctuations. The amplitude of the interfacial dielectric constant oscillation can exceed the dielectric constant of the bulk water, suggesting an ice-like frozen state, but only if there are no electrolyte ions. The interfacial accumulation of electrolyte ions causes a decrease in the dielectric constant due to the reduction of water density and re-orientation of water dipoles in ion hydration shells. Finally, we show how to use the computed dielectric properties to estimate SCM's capacitances.

58 GEOSCIENCES↗
NASA NTRS2000

Solid Lubrication Fundamentals and Applications

This chapter focuses attention on the friction and wear properties of selected solid lubricating films to aid users in choosing the best lubricant, deposition conditions, and operational variables. For simplicity, discussion of the tribological properties of concern is separated into two parts. The first part of the chapter discusses the different solid lubricating films selected for study including commercially developed solid film lubricants: (1) bonded molybdenum disulfide (MoS2), (2) magnetron-sputtered MoS2, (3) ion-plated silver, (4) ion-plated lead, (5) magnetron-sputtered diamondlike carbon (MS DLC), and (6) plasma-assisted, chemical-vapor-deposited diamondlike carbon (PACVD DEC) films. Marked differences in the friction and wear properties of the different films resulted from the different environmental conditions (ultrahigh vacuum, humid air, and dry nitrogen) and the solid film lubricant materials. The second part of the chapter discusses the physical and chemical characteristics, friction behavior, and endurance life of the magnetron-sputtered MoS2 films. The role of interface species and the effects of applied load, film thickness, oxygen pressure, environment, and temperature on the friction and wear properties are considered.

Miyoshi, Kazuhisa↗
DOE OSTI2025

Unveiling the Reactivity of Fluoropolymers with Sodium Metal: Mechanistic Insights and Battery Implications

Poly(vinylidene fluoride-co-hexafluoropropylene) (PVDF-HFP) has broad applications across various metal-ion battery systems, such as a binder for electrodes, a supporting matrix for electrolytes, and a separator material. Due to its excellent mechanical properties, PVDF-HFP has become an excellent candidate for fabricating gel and solid-state electrolytes in sodium-based batteries. However, in this study, we noticed notable side reactions occurring at the interface of PVDF-HFP membranes and Na metal. These reactions not only alter chemical compositions but also further affect the surface morphology and adhesive properties of membranes. Similar phenomena are observed in other polyfluoroalkyl-based membranes (PVDF and PTFE). Therefore, we systematically studied the reaction mechanisms between the Na metal and these polymers. The influence of different functional groups (−F, −CF 3 , −H) and their arrangement on the reaction extent has also been discussed. Finally, we concluded with the key factors driving these side reactions and provided new perspectives for designing polymers tailored for sodium-based batteries.

Interfacial side reactions↗
DOE OSTI2025

Electrical Transport Interplay with Charge Density Waves, Magnetization, and Disorder Tuned by 2D van der Waals Interface Modification via Elemental Intercalation and Substitution in ZrTe 3 , 2H-TaS 2 , and Cr 2 Si 2 Te 6 Crystals

Electrical transport in 2D materials exhibits unique behaviors due to reduced dimensionality, broken symmetries, and quantum confinement. It serves as both a sensitive probe for the emergence of coherent electronic phases and a tool to actively manipulate many-body correlated states. Exploring their interplay and interdependence is crucial but remains underexplored. This review integratively cross-examines the atomic and electronic structures and transport properties of van der Waals-layered crystals ZrTe 3 , 2H-TaS 2 , and Cr 2 Si 2 Te 6 , providing a comprehensive understanding and uncovering new discoveries and insights. A common observation from these crystals is that modifying the atomic and electronic interface structures of 2D van der Waals interfaces using heteroatoms significantly influences the emergence and stability of coherent phases, as well as phase-sensitive transport responses. In ZrTe 3 , substitution and intercalation with Se, Hf, Cu, or Ni at the 2D vdW interface alter phonon–electron coupling, valence states, and the quasi-1D interface Fermi band, affecting the onset of CDW and SC, manifested as resistance upturns and zero-resistance states. We conclude here that these phenomena originate from dopant-induced variations in the lattice spacing of the quasi-1D Te chains of the 2D vdW interface, and propose an unconventional superconducting mechanism driven by valence fluctuations at the van Hove singularity, arising from quasi-1D lattice vibrations. Short-range in-plane electronic heterostructures at the vdW interface of Cr 2 Si 2 Te 6 result in a narrowed band gap. The sharp increase in in-plane resistance is found to be linked to the emergence and development of out-of-plane ferromagnetism. The insertion of 2D magnetic layers such as Mn, Fe, and Co into the vdW gap of 2H-TaS 2 induces anisotropic magnetism and associated transport responses to magnetic transitions. Overall, 2D vdW interface modification offers control over collective electronic behavior, transport properties, and their interplays, advancing fundamental science and nanoelectronic devices.

2D van der Waals interfaces↗
DOE OSTI2021

Characterize traction–separation relation and interfacial imperfections by data-driven machine learning models

Abstract Interfacial mechanical properties are important in composite materials and their applications, including vehicle structures, soft robotics, and aerospace. Determination of traction–separation (T–S) relations at interfaces in composites can lead to evaluations of structural reliability, mechanical robustness, and failures criteria. Accurate measurements on T–S relations remain challenging, since the interface interaction generally happens at microscale. With the emergence of machine learning (ML), data-driven model becomes an efficient method to predict the interfacial behaviors of composite materials and establish their mechanical models. Here, we combine ML, finite element analysis (FEA), and empirical experiments to develop data-driven models that characterize interfacial mechanical properties precisely. Specifically, eXtreme Gradient Boosting (XGBoost) multi-output regressions and classifier models are harnessed to investigate T–S relations and identify the imperfection locations at interface, respectively. The ML models are trained by macroscale force–displacement curves, which can be obtained from FEA and standard mechanical tests. The results show accurate predictions of T–S relations ( R 2 = 0.988) and identification of imperfection locations with 81% accuracy. Our models are experimentally validated by 3D printed double cantilever beam specimens from different materials. Furthermore, we provide a code package containing trained ML models, allowing other researchers to establish T–S relations for different material interfaces.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Effects of Cooling Rates on Self-Reinforced Polyethylene Composites via Multiscale Experimental Characterizations

This study investigates cooling rate effects on self-reinforced polyethylene composites (PE-SRCs), focusing on macroscopic thermal and mechanical properties, and micro- to sub-micron morphologies. Unlike conventional fiber-reinforced composites (e.g., carbon and glass fiber based), cooling rate effects on SRCs remain less studied. PE-SRCs were fabricated from ultrahigh molecular weight polyethylene (UHMWPE) fabric and high-density polyethylene (HDPE) matrix with cooling rates from 0.76°C/min to 425°C/min. Tensile properties, including modulus (~9.5 GPa), strength (~360 MPa), and ultimate elongation (~6%), showed no significant cooling rate dependence. However, thermal analysis revealed melting temperature of HDPE matrix decreased from ~130°C to ~126.5°C with total crystallinity reducing from 84% to 75%. Wide-angle X-ray scattering confirmed decreasing crystallinity but at lower values (67%–59%), as thermal analysis might overestimate crystallinity. Small-angle X-ray scattering showed long period decreased from 24.3 to 17.4 nm, attributed to thinner lamellae and closer packing, suggesting limited crystal growth and faster nucleation. While cooling rates alter melting point and crystalline structures of HDPE matrix, these changes do not translate to substantial impact on macroscale mechanical properties of PE-SRCs. This insensitivity might be related to low interfacial shear strength at PE/PE interfaces, enabling process optimization while maintaining mechanical performance.

Cooling rate↗
DOE OSTI2023

High-temperature thermal conductivity measurements of macro-porous graphite

Graphite is a unique material for high temperature applications and will likely become increasingly important as we attempt to electrify industrial applications. Here, we investigate the thermal properties of low-quality, macro-porous graphite to determine the tradeoff between quality and thermal performance. We use laser flash analysis (LFA) to measure the thermal diffusivity of graphite at high temperatures. Due to the large pores in the graphite samples preventing uniform laser flash heating, we first apply a thick coating to achieve the required flat, parallel surfaces for LFA measurements. We then develop a methodology to extract properties of the sample from the diffusivity measurements, based on finite element modeling of a variety of sample/coating interface profiles. Validating the methodology against a reference sample demonstrates a mean absolute percentage error of 8.5%, with potential improvement with better sample characterization. We show low-quality graphite has a thermal conductivity of ~10 W/m/K up to 1000 °C, which is an order of magnitude lower than high-quality graphite, but contributions from photon conductivity may result in higher conductivities at higher temperatures. Overall, we demonstrate an approach for measuring thermal properties of macro-porous materials at high temperatures, and apply the approach to measuring thermal conductivity of porous graphite, which will aid in the design of hightemperature systems for cost-competitive decarbonization.

Verma, Shomik↗
DOE OSTI2020

Emerging Potassium Metal Anodes: Perspectives on Control of the Electrochemical Interfaces

Potassium metal acts as the anode in emerging potassium metal batteries (KMBs). It further serves as the counter-electrode for potassium ion battery (KIB) half-cells, with its reliable performance being critical for assessing the working electrode material. This first-of-its-kind critical review focuses on the dual challenge of controlling the potassium metal-substrate and the potassium metal-electrolyte interface so as to prevent dendrites. The discussion begins with a comparison of the physical and chemical properties of K metal anodes versus the much oft studied Li and Na metal anodes. Based on established descriptions for root causes of dendrites, the problem should be less severe for K than for Li or Na, while in fact the opposite is observed. The key reason that the K metal surface rapidly becomes dendritic in common electrolytes is its unstable solid electrolyte interphase (SEI). An unstable SEI layer is defined as being non self-passivating. No SEI is perfectly stable during cycling, and all SEI structures are heterogenous both vertically and horizontally relative to the electrolyte interface. The difference between a "stable" and an "unstable" SEI may be viewed as the relative degree to which during cycling it thickens and becomes further heterogeneous. The unstable SEI on K metal leads to a number of interrelated problems, such as low cycling Coulombic efficiency (CE), a severe impedance rise, large overpotentials, and possibly electrical shorting; all of which have been reported to occur as early as in the first ten plating/stripping cycles. Many of the traditional "interface fixes" employed for Li and Na metal anodes, such as various artificial SEIs, surface membranes, barrier layers, the secondary separators, etc. have not been attempted or optimized for the case of K. This is an important area for further exploration, with an understanding that success may come harder than with Li due to K-based SEI reactivity with both ether and ester solvents. The second critical problem with K metal anodes is that they do not thermally or electrochemically wet a standard (untreated) Cu foil current collector. Published experimental and modeling research directly highlights the weak bonding between the K atoms and a Cu surface. Existing surface treatment approaches that achieve improved K wetting are discussed, along with the general design rules for future studies. Also discussed are geometry-based methods to tune nucleation, as well dual approaches where nucleation and SEI structure are tuned through complementary schemes to achieve extended half-cell and full battery stability. We hypothesize that K metal never achieves a planar wetting morphology even at cycle one making the dendrites "baked-in". We propose that classical thin films growth models, Frank van der Merwe (F-M), Volmer Weber (V-W), and Stranski-Krastanov (S-K) can be employed to describe early stage plating behavior. It is demonstrated that island-like V-W growth is the applicable description for the natural plating behavior of K on pristine Cu. Moving forward, there are three inter-related thrusts to be pursued: First, K salt - based electrolyte formulations have to mature and become further tailored to handle the increased reactivity of a metal rather than an ion anode. Second, the K-based SEI structure needs to be further understood and ultimately tuned to be less reactive. Third, the energetics of the K metal - current collector interface must be controlled to promote planar wetting/dewetting throughout cycling.

25 ENERGY STORAGE↗
DOE OSTI2021

Process Robustness of Friction Stir Dovetailing of AA7099 to Steel with In Situ AA6061 Interlayer Linking

Rolled homogeneous armor (RHA) steel was joined to aluminum alloy AA7099-T7451 with an interlayer of AA6061-T651 by single-pass friction stir dovetailing (FSD). The evolution of the AA7099/AA6061 and AA6061/steel interfacial microstructures and mechanical properties of the joints were examined. The AA6061 interlayer prevented Zn migration from the AA7099 to the steel during FSD and inhibited formation of brittle, Zn-rich, Fe-Al intermetallic compounds within the joint. The flow of the dissimilar aluminum alloys at the interface was controlled with a partially threaded friction stir welding pin and a dovetail geometric configuration optimized in previous work. The process robustness was explored to establish parameter windows by elucidating parametric effects on the relationship between process conditions and the joint’s microstructure and mechanical properties. Ranges of process forge force (50–70 kN) and interface temperature (475°C–515°C) consistently yielded joints with adequate performance. The load-carrying capacity in lap shear tensile tests was 2025 ? 100 N/mm.

Friction Stir Dovetailing, FSW, AA7099, AA6061, Ro↗
DOE OSTI2026

The Impact of Catalyst Layer Properties and Transport Layer Interactions on Rutile Iridium Oxide Anode Durability in Proton Exchange Membrane Electrolysis

The performance and durability of proton exchange membrane electrolyzers are governed by the interplay between catalyst inks and integration, the resulting catalyst layer (CL) structure, and the interface between the CL and porous transport layer (PTL). In this study, we examine the impact of these factors on cell performance and degradation using rutile IrO2 as the anode catalyst. The two commercial rutile IrO2 catalysts differ from each other and from the amorphous benchmark in agglomerate size and ink stability. This has implications for CL morphology and electronic connectivity. In the context of long-term pressurized operation for up to 1000 h, the selection of the PTL has a marginal effect on the durability of CLs with good structural and compositional properties, whereas poorly structured CLs demonstrate evident degradation when used with certain PTLs. Poor CL properties amplify the negative impact of non-ideal CL/PTL interfaces, demonstrating the importance of CL/PTL interface engineering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Davis Computational Spectroscopy Workflow—From Structure to Spectra

Here, we describe an automated workflow that connects a series of atomic simulation tools to investigate the relationship between atomic structure, lattice dynamics, materials properties, and inelastic neutron scattering (INS) spectra. Starting from the atomic simulation environment (ASE) as an interface, we demonstrate the use of a selection of calculators, including density functional theory (DFT) and density functional tight binding (DFTB), to optimize the structures and calculate interatomic force constants. We present the use of our workflow to compute the phonon frequencies and eigenvectors, which are required to accurately simulate the INS spectra in crystalline solids like diamond and graphite as well as molecular solids like rubrene. We have also implemented a machine-learning force field based on Chebyshev polynomials called the Chebyshev interaction model for efficient simulation (ChIMES) to improve the accuracy of the DFTB simulations. We then explore the transferability of our DFTB/ChIMES models by comparing simulations derived from different training sets. We show that DFTB/ChIMES demonstrates ~100× reduction in computational expense while retaining most of the accuracy of DFT as well as yielding high accuracy for different materials outside of our training sets. The DFTB/ChIMES method within the workflow expands the possibilities to use simulations to accurately predict materials properties of increasingly complex structures that would be unfeasible with ab initio methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1986

Analysis of Delamination Growth from Matrix Cracks in Laminates Subjected to Bending Loads

A major source of delamination damage in laminated composite materials is from low-velocity impact. In thin composite laminates under point loads, matrix cracks develop first in the plies, and delaminations then grow from these cracks at the ply interfaces. The purpose of this study was to quantify the combined effects of bending and transverse shear loads on delamination initiation from matrix cracks. Graphite-epoxy laminates with 90 deg. plies on the outside were used to provide a two-dimensional simulation of the damage due to low-velocity impact. Three plate bending problems were considered: a 4-point bending, 3-point bending, and an end-clamped center-loaded plate. Under bending, a matrix crack will form on the tension side of the laminate, through the outer 90 deg. plies and parallel to the fibers. Delaminations will then grow in the interface between the cracked 90 deg. ply and the next adjacent ply. Laminate plate theory was used to derive simple equations relating the total strain energy release rate, G, associated with the delamination growth from a 90 deg. ply crack to the applied bending load and laminate stiffness properties. Three different lay-ups were tested and results compared. Test results verified that the delamination always formed at the interface between the cracked 90 deg. ply and the next adjacent ply. Calculated values for total G sub c from the analysis showed good agreement for all configurations. The analysis was able to predict the delamination onset load for the cases considered. The result indicated that the opening mode component (Mode I) for delamination growth from a matrix crack may be much larger than the component due to interlaminar shear (Mode II).

Murri, G. B.↗
NASA NTRS1988

Analysis of delamination growth from matrix cracks in laminates subjected to bending loads

A major source of delamination damage in laminated composite materials is from low-velocity impact. In thin composite laminates under point loads, matrix cracks develop first in the plies, and delaminations then grow from these cracks at the ply interfaces. The purpose of this study was to quantify the combined effects of bending and transverse shear loads on delamination initiation from matrix cracks. Graphite-epoxy laminates with 90 deg plies on the outside were used to provide a two-dimensional simulation of the damage due to low-velocity impact. Three plate bending problems were considered: a 4-point bending, 3-point bending, and an end-clamped center-loaded plate. Under bending, a matrix crack will form on the tension side of the laminate, through the outer 90 deg plies and parallel to the fibers. Delaminations will then grow in the interface between the cracked 90 deg ply and the next adjacent ply. Laminate plate theory was used to derive simple equations relating the total strain energy release rate, G, associated with the delamination growth from a 90 deg ply crack to the applied bending load and laminate stiffness properties. Three different lay-ups were tested and results compared. Test results verified that the delamination always formed at the interface between the cracked 90 deg ply and the next adjacent ply. Calculated values for total G sub c from the analysis showed good agreement for all configurations. The analysis was able to predict the delamination onset load for the cases considered. The result indicated that the opening mode component (Mode I) for delamination growth from a matrix crack may be much larger than the component due to interlaminar shear (Mode II).

Murri, Gretchen Bostaph↗