SEARCH · Search NASA

Results for “mixed mode fracture”

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 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Sublaminate analysis of interlaminar fracture in composites. I - Analytical model

A simple analysis method based on a transverse shear-deformation theory and a sublaminate approach is utilized to analyze a mixed-mode edge delamination specimen. The analysis provides closed-form expressions for the interlaminar shear stresses ahead of the crack, the total strain-energy release rate, and the strain-energy release-rate components. The parameters controlling the behavior are identified. The effect of specimen stacking sequence and delamination interface on the strain energy release rate components is investigated. Results are compared with a finite-element simulation for reference. The simple nature of the method makes it suitable for preliminary design analyses which require a large number of configurations to be evaluated quickly and economically.

Armanios, Erian A.↗

NASA NTRS2001

Fatigue Life Methodology for Bonded Composite Skin/Stringer Configurations

A methodology is presented for determining the fatigue life of composite structures based on fatigue characterization data and geometric nonlinear finite element (FE) analyses. To demonstrate the approach, predicted results were compared to fatigue tests performed on specimens which represented a tapered composite flange bonded onto a composite skin. In a first step, tension tests were performed to evaluate the debonding mechanisms between the flange and the skin. In a second step, a 2D FE model was developed to analyze the tests. To predict matrix cracking onset, the relationship between the tension load and the maximum principal stresses transverse to the fiber direction was determined through FE analysis. Transverse tension fatigue life data were used to -enerate an onset fatigue life P-N curve for matrix cracking. The resulting prediction was in good agreement with data from the fatigue tests. In a third step, a fracture mechanics approach based on FE analysis was used to determine the relationship between the tension load and the critical energy release rate. Mixed mode energy release rate fatigue life data were used to create a fatigue life onset G-N curve for delamination. The resulting prediction was in good agreement with data from the fatigue tests. Further, the prediction curve for cumulative life to failure was generated from the previous onset fatigue life curves. The results showed that the methodology offers a significant potential to Predict cumulative fatigue life of composite structures.

Krueger, Ronald↗

NASA NTRS2024

Impact Response of Thermoplastic Composites – Experiments and Modeling

This study investigates the impact response of thermoplastic composite materials. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK reinforced with T700G (T700/LMPAEK ). Several High Energy Dynamic Impact (HEDI) tests were performed wherein an aluminum projectile impacted T700/LMPAEK panels at different velocities. In the impact experiment studied herein, the projectile impacted the T700/LMAPEK panel at 73.2 m/s and rebounded. The ensuing damage to the T700/LMPAEK panel was characterized via ultrasonic C-scans, visualizing damage modes including matrix cracking and delaminations. Since the projectile rebounded, the fiber breakage damage mode was not evident. The impact event was modeled using a combined continuum damage mechanics and cohesive zone analysis approach with the commercially available hydrocode LS-DYNA®. The dynamic deformation and damage of the continuum, i.e., the plies in the T700/LMPAEK laminate, was modeled using a Deformation Gradient Decomposition (DGD) material model, i.e., MAT299 available within LS-DYNA. The interlaminar regions were treated as zero-thickness cohesive zones, wherein the dynamic delamination mechanics was modeled using a traction-separation formulation. Mixed-mode delamination growth was captured using the Benzegaggh-Kenane (B-K) traction-separation interpolation law. Strain-rate sensitivity of interlaminar fracture toughness was included using approximations based on published literature. The outcomes of the modeling effort are discussed herein which demonstrate good correlations with test data regarding panel deflection, panel-projectile interactions, damage modes and damage extent.

Impact modeling↗

NASA NTRS1984

Application of Iosipescu specimen geometry to determination of shear strength in unidirectional composites

The Iosipescu (1967) test method, in conjunction with a finite element analysis of anisotropic, unidirectionally reinforced composite specimens, can be used to derive composite shear strength values despite the small size of the samples used. When such a specimen is loaded parallel to the fiber direction, fracture occurs by shear within a narrow region near the center of the specimen; this is consistent with the maximum shear area determined by finite element analysis. The resulting fracture surface, as observed, is characteristic of shear failure, occurring either within the matrix or at the fiber-matrix interface. Problems of compressive and mixed mode failure which accrue to short beam shear specimens are by the present method precluded.

Hurwitz, F. I.↗

NASA NTRS2012

Design and Characterization of Thin Stainless Steel Burst Disks for Increasing Two-Stage Light Gas Launcher Efficiency

Laser etched 300 series Stainless Steel Burst Disks (SSBD) ranging between 0.178 mm (0.007-in.) and 0.508mm (0.020-in.) thick were designed for use in a 17-caliber two-stage light gas launcher. First, a disk manufacturing method was selected using a combination of wire electrical discharge machining (EDM) to form the blank disks and laser etching to define the pedaling fracture pattern. Second, a replaceable insert was designed to go between the SSDB and the barrel. This insert reduced the stress concentration between the SSBD and the barrel, providing a place for the petals of the SSDB to open, and protecting the rifling on the inside of the barrel. Thereafter, a design of experiments was implemented to test and characterize the burst characteristics of SSBDs. Extensive hydrostatic burst testing of the SSBDs was performed to complete the design of experiments study with one-hundred and seven burst tests. The experiment simultaneously tested the effects of the following: two SSBD material states (full hard, annealed); five SSBD thicknesses 0.178, 0.254, 0.305, 0.381 mm (0.007, 0.010, 0.012, 0.015, 0.020-in.); two grain directions relative); number of times the laser etch pattern was repeated (varies between 5-200 times); two heat sink configurations (with and without heat sink); and, two barrel configurations (with and without insert). These tests resulted in the quantification of the relationship between SSBD thickness, laser etch parameters, and desired burst pressure. Of the factors investigated only thickness and number of laser etches were needed to develop a mathematical relationship predicting hydrostatic burst pressure of disks using the same barrel configuration. The fracture surfaces of two representative SSBD bursts were then investigated with a scanning electron microscope, one burst hydrostatically in a fixture and another dynamically in the launcher. The fracture analysis verified that both burst conditions resulted in a ductile overload failure indicated by transgranular microvoid coalescence, non-fragmenting rupture and mixed tensile and shear failure modes, regardless of the material states tested. More testing is underway to determine the relationship between SSBD burst pressure and projectile velocity.

Tylka, Jonathan M.↗

NASA NTRS2006

Fracture Mechanics for Composites: State of the Art and Challenges

Interlaminar fracture mechanics has proven useful for characterizing the onset of delaminations in composites and has been used with limited success primarily to investigate onset in fracture toughness specimens and laboratory size coupon type specimens. Future acceptance of the methodology by industry and certification authorities however, requires the successful demonstration of the methodology on the structural level. In this paper, the state-of-the-art in fracture toughness characterization, and interlaminar fracture mechanics analysis tools are described. To demonstrate the application on the structural level, a panel was selected which is reinforced with stringers. Full implementation of interlaminar fracture mechanics in design however remains a challenge and requires a continuing development effort of codes to calculate energy release rates and advancements in delamination onset and growth criteria under mixed mode conditions.

Krueger, Ronald↗

NASA NTRS1984

On stress field near a stationary crack tip

It is well known that the stress and elastic-plastic deformation fields near a crack tip have important roles in the corresponding fracture process. For elastic-perfectly-plastic solids, different solutions are given in the literature. In this work several of these solutions are examined and compared for Mode I (tension), Mode II (shear), and mixed Modes I and II loading conditions in plane strain. By consideration of the dynamic solution, it is shown that the assumption that the material is yielding all around a crack tip may not be reasonable in all cases. By admitting the existence of some elastic sectors, continuous stress fields are obtained even for mixed Modes I and II.

Nemat-Nasser, S.↗

NASA NTRS2019

Benchmarking Mixed Mode Failure in Progressive Damage and Failure Analysis Methods

The verification and validation of progressive-damage-analysis finite element methods are difficult but critical tasks to undertake during their development. Verification exercises assess whether a predictive analysis tool produces results that are consistent with the fundamental concepts and assumptions of the tool under evaluation. Ideally, closed-form analytical solutions can be derived for which method verification results can be compared. Problems selected for computational tool verification are often simple and isolate individual features of the tool. In the case of progressive damage finite element methods, verifications should be performed to evaluate the ability of the model to predict the initiation of damage and its growth through the finite element mesh under a variety of conditions. Mabson et al. proposed a test case of a unidirectional, fiber-reinforced plate with a center crack subjected to tensile loads to evaluate matrix crack propagation predictions. The problem was modeled using the Abaqus Hashin continuum damage mechanics (CDM) model for fiber-reinforced composites. Different combinations of matrix strengths and element sizes were used in the simulations, and the results were compared to a closed-form solution based on linear elastic fracture mechanics (LEFM). It was determined that the Abaqus CDM model could predict the LEFM solution of Mode I cracks only when the finite element mesh density met specific requirements based on the material properties. This paper presents closed-form LEFM solutions for a center notch mixed mode (CNMM) verification problem. Parametric finite element analyses were developed using progressive damage analysis methods of both the Discrete Damage Mechanics (DDM) and CDM classes. The progressive damage analysis methods applied in the analyses of the CNMM problem include CompDam and the Floating Node Method. Analyses were conducted with various mode mixities and element sizes to verify that the damage models were working as intended and to identify any limits of applicability.

Leone, Frank↗

NASA NTRS1985

On the effect of quarter-point element size on fracture criteria

The effect of the quarter-point element size on the solution parameters of fracture criteria is investigated in view of its effect on the computed stress intensity factors. An analytical relationship between the error in calculating crack propagation increment and stress intensity for single mode problems is obtained while, for mixed mode problems, the investigation is based on numerical experimentation. The singular element size is found to have a substantial effect on both crack propagation angle and increment for the geometries considered.

Yehia, N. A. B.↗

NASA NTRS1989

A simplified approach to strain energy release rate computations for interlaminar fracture of composites

A simple approach for the strain energy release rate computations based on the finite element method and a singular fitting model is presented. The model uses the stress and displacement distributions at the delamination front. The method is applied to a mixed-mode double cracked-lap-shear composite configuration. The strain energy release rate components predicted by the model are compared with the finite element crack-closure method. The effect of the mesh size on the stress and displacement distribution is isolated. The strain energy release rates predicted by relatively coarse mesh sizes are in good agreement with the finite element crack closure method.

Armanios, Erian A.↗

NASA NTRS2008

Assessment of Crack Path Prediction in Non-Proportional Mixed-Mode Fatigue

Non-proportional mixed-mode loading is present in many systems and a growing crack can experience any manner of mixed-mode loading. Prediction of the resulting crack path is important when assessing potential failure modes or when performing a failure investigation. Current crack path selection criteria are presented along with data for Inconel 718 under non-proportional mixed-mode loading. Mixed-mode crack growth can transition between path deflection mechanisms with very different orientations. Non-proportional fatigue loadings lack a single parameter for input to current crack path criteria. Crack growth transitions were observed in proportional and non-proportional FCG tests. Different paths displayed distinct fracture surface morphologies. New crack path drivers & transition criteria must be developed.

Highsmith, Shelby, Jr.↗

NASA NTRS2001

Numerical Simulation of Delamination Growth in Composite Materials

The use of decohesion elements for the simulation of delamination in composite materials is reviewed. The test methods available to measure the interfacial fracture toughness used in the formulation of decohesion elements are described initially. After a brief presentation of the virtual crack closure technique, the technique most widely used to simulate delamination growth, the formulation of interfacial decohesion elements is described. Problems related with decohesion element constitutive equations, mixed-mode crack growth, element numerical integration and solution procedures are discussed. Based on these investigations, it is concluded that the use of interfacial decohesion elements is a promising technique that avoids the need for a pre-existing crack and pre-defined crack paths, and that these elements can be used to simulate both delamination onset and growth.

Camanho, P. P.↗

NASA NTRS1994

Dynamic fracture mechanics analysis for an edge delamination crack

A global/local analysis is applied to the problem of a panel with an edge delamination crack subject to an impulse loading to ascertain the dynamic J integral. The approach uses the spectral element method to obtain the global dynamic response and local resultants to obtain the J integral. The variation of J integral along the crack front is shown. The crack behavior is mixed mode (Mode 2 and Mode 3), but is dominated by the Mode 2 behavior.

Rizzi, Stephen A.↗

NASA NTRS2007

Analysis of Composite Panel-Stiffener Debonding Using a Shell/3D Modeling Technique

Interlaminar fracture mechanics has proven useful for characterizing the onset of delaminations in composites and has been used successfully primarily to investigate onset in fracture toughness specimens and laboratory size coupon type specimens. Future acceptance of the methodology by industry and certification authorities, however, requires the successful demonstration of the methodology on the structural level. For this purpose, a panel was selected that is reinforced with stiffeners. Shear loading causes the panel to buckle, and the resulting out-of-plane deformations initiate skin/stiffener separation at the location of an embedded defect. A small section of the stiffener foot, web and noodle as well as the panel skin in the vicinity of the delamination front were modeled with a local 3D solid model. Across the width of the stiffener foot, the mixedmode strain energy release rates were calculated using the virtual crack closure technique. A failure index was calculated by correlating the results with a mixed-mode failure criterion of the graphite/epoxy material. Computed failure indices were compared to corresponding results where the entire web was modeled with shell elements and only a small section of the stiffener foot and panel were modeled locally with solid elements. Including the stiffener web in the local 3D solid model increased the computed failure index. Further including the noodle and transition radius in the local 3D solid model changed the local distribution across the width. The magnitude of the failure index decreased with increasing transition radius and noodle area. For the transition radii modeled, the material properties used for the noodle area had a negligible effect on the results. The results of this study are intended to be used as a guide for conducting finite element and fracture mechanics analyses of delamination and debonding in complex structures such as integrally stiffened panels.

Krueger, Ronald↗

NASA NTRS1992

Determination of stress intensity factors for interface cracks under mixed-mode loading

A simple technique was developed using conventional finite element analysis to determine stress intensity factors, K1 and K2, for interface cracks under mixed-mode loading. This technique involves the calculation of crack tip stresses using non-singular finite elements. These stresses are then combined and used in a linear regression procedure to calculate K1 and K2. The technique was demonstrated by calculating three different bimaterial combinations. For the normal loading case, the K's were within 2.6 percent of an exact solution. The normalized K's under shear loading were shown to be related to the normalized K's under normal loading. Based on these relations, a simple equation was derived for calculating K1 and K2 for mixed-mode loading from knowledge of the K's under normal loading. The equation was verified by computing the K's for a mixed-mode case with equal and normal shear loading. The correlation between exact and finite element solutions is within 3.7 percent. This study provides a simple procedure to compute K2/K1 ratio which has been used to characterize the stress state at the crack tip for various combinations of materials and loadings. Tests conducted over a range of K2/K1 ratios could be used to fully characterize interface fracture toughness.

Naik, Rajiv A.↗

NASA NTRS1998

Three-Dimensional Gear Crack Propagation Studies

Three-dimensional crack growth simulation was performed on a split-tooth gear design using boundary element modeling and linear elastic fracture mechanics. Initial cracks in the fillet of the teeth produced stress intensity factors of greater magnitude (and thus, greater crack growth rates) than those in the root or groove areas of the teeth. Crack growth simulation was performed on a case study to evaluate crack propagation paths. Tooth fracture was predicted from the crack growth simulation for an initial crack in the tooth fillet region. Tooth loads on the uncracked mesh of the split-tooth design were up to five times greater than those on the cracked mesh if equal deflections of the cracked and uncracked teeth were considered. Predicted crack shapes as well as crack propagation life are presented based on calculated stress intensity factors, mixed-mode crack propagation trajectory theories, and fatigue crack growth theories.

Lewicki, David G.↗

NASA NTRS1997

An Insidious Mode of Oxidative Degradation in a SiC-SiC Composite

The oxidative durability of a SiC-SiC composite with Hi-Nicalon fiber and BN interphase was investigated at 800 C (where pesting is known to occur in SiC-SiC composites) for exposure durations of up to 500 hours and in a variety of oxidant mixes and flow rates, ranging from quasi-stagnant room air, through slow flowing O2 containing 30-90% H2O, to the high-velocity flame of a burner rig. Degradation of the composite was determined from residual strength and fracture strain in post-exposure mechanical tests and correlated with microstructural evidence of damage to fiber and interphase. The severest degradation of composite behavior was found to occur in the bumer rig, and is shown to be connected with the high oxidant velocity and substantial moisture content, as well as a thin sublayer of carbon indicated to form between fiber and interphase during composite processing.

Ogbuji, Linus U. J. T.↗

NASA NTRS1998

Fracture Mechanics Analysis of Stitched Stiffener-Skin Debonding

An analysis based on plate finite elements and the virtual crack closure technique has been implemented to study the effect of stitching on mode I and mode II strain energy release rates for debond configurations. The stitches were modeled as discrete nonlinear fastener elements with a compliance determined by experiment. The axial and shear behavior of the stitches was considered, however, the two compliances and failure loads were assumed to be independent. Both a double cantilever beam (mode I) and a mixed mode skin-stiffener debond configuration were studied. In the double cantilever beam configurations, G(sub I) began to decrease once the debond had grown beyond the first row of stitches and was reduced to zero for long debonds. In the mixed-mode skin-stiffener configurations, G(sub I) showed a similar behavior as in the double cantilever beam configurations, however, G(sub u), continued to increase with increasing debond length.

Glaessgen, E. H.↗