An approximate solution for interlaminar stresses in laminated composites: Applied mechanics program

An approximate solution for interlaminar stresses in finite width, laminated composites subjected to uniform extensional, and bending loads is presented. The solution is based upon the principle of minimum complementary energy and an assumed, statically admissible stress state, derived by considering local material mismatch effects and global equilibrium requirements. The stresses in each layer are approximated by polynomial functions of the thickness coordinate, multiplied by combinations of exponential functions of the in-plane coordinate, expressed in terms of fourteen unknown decay parameters. Imposing the stationary condition of the laminate complementary energy with respect to the unknown variables yields a system of fourteen non-linear algebraic equations for the parameters. Newton's method is implemented to solve this system. Once the parameters are known, the stresses can be easily determined at any point in the laminate. Results are presented for through-thickness and interlaminar stress distributions for angle-ply, cross-ply (symmetric and unsymmetric laminates), and quasi-isotropic laminates subjected to uniform extension and bending. It is shown that the solution compares well with existing finite element solutions and represents an improved approximate solution for interlaminar stresses, primarily at interfaces where global equilibrium is satisfied by the in-plane stresses, but large local mismatch in properties requires the presence of interlaminar stresses

Estimation of Delamination Crack Depth Using Passive Thermography

Passive thermography is used to monitor small increases in temperature resulting from delamination damage formation in a composite hat-stiffened panel during quasi-static loading. The heating is composed of two heat generation components. The first component is an instantaneous response due to a strain release during quasi-static loading. The second component is mechanical heating, at the interface of failure, due to fracture damage. This second component produces a transient rise in temperature that is a function of the damage depth and thermal diffusivity. The first component defines the thermal start time for the transient response. A one-dimensional thermal model is used to determine the damage depth. The results are compared to ultrasonic and X-ray CT data. The advantages and limitations of the thermal technique for damage depth detection are discussed

Continuum Damage Mechanics Models for the Analysis of Progressive Failure in Open-Hole Tension Laminates

The performance of a state-of-the-art continuum damage mechanics model for interlaminar damage, coupled with a cohesive zone model for delamination is examined for failure prediction of quasi-isotropic open-hole tension laminates. Limitations of continuum representations of intra-ply damage and the effect of mesh orientation on the analysis predictions are discussed. It is shown that accurate prediction of matrix crack paths and stress redistribution after cracking requires a mesh aligned with the fiber orientation. Based on these results, an aligned mesh is proposed for analysis of the open-hole tension specimens consisting of different meshes within the individual plies, such that the element edges are aligned with the ply fiber direction. The modeling approach is assessed by comparison of analysis predictions to experimental data for specimen configurations in which failure is dominated by complex interactions between matrix cracks and delaminations. It is shown that the different failure mechanisms observed in the tests are well predicted. In addition, the modeling approach is demonstrated to predict proper trends in the effect of scaling on strength and failure mechanisms of quasi-isotropic open-hole tension laminates

STAGS Developments for Residual Strength Analysis Methods for Metallic Fuselage Structures

A summary of advances in the Structural Analysis of General Shells (STAGS) finite element code for the residual strength analysis of metallic fuselage structures, that were realized through collaboration between the structures group at NASA Langley, and Dr. Charles Rankin is presented. The majority of the advancements described were made in the 1990's under the NASA Airframe Structural Integrity Program (NASIP). Example results from studies that were conducted using the STAGS code to develop improved understanding of the nonlinear response of cracked fuselage structures subjected to combined loads are presented. An integrated residual strength analysis methodology for metallic structure that models crack growth to predict the effect of cracks on structural integrity is demonstrate

Guidelines and Parameter Selection for the Simulation of Progressive Delamination

Turon s methodology for determining optimal analysis parameters for the simulation of progressive delamination is reviewed. Recommended procedures for determining analysis parameters for efficient delamination growth predictions using the Abaqus/Standard cohesive element and relatively coarse meshes are provided for single and mixed-mode loading. The Abaqus cohesive element, COH3D8, and a user-defined cohesive element are used to develop finite element models of the double cantilever beam specimen, the end-notched flexure specimen, and the mixed-mode bending specimen to simulate progressive delamination growth in Mode I, Mode II, and mixed-mode fracture, respectively. The predicted responses are compared with their analytical solutions. The results show that for single-mode fracture, the predicted responses obtained with the Abaqus cohesive element correlate well with the analytical solutions. For mixed-mode fracture, it was found that the response predicted using COH3D8 elements depends on the damage evolution criterion that is used. The energy-based criterion overpredicts the peak loads and load-deflection response. The results predicted using a tabulated form of the BK criterion correlate well with the analytical solution and with the results predicted with the user-written element

Experimental and Numerical Analysis of Skin-Stiffener Separation Using a Seven-Point Bend Configuration

Skin-stiffener separation in stiffened composite panels consists of a complex interaction between multiple scales of progressive damage and failure mechanisms. This work used a superposed cohesive element method of modeling resistance curve effects that represents the structural interface of a unidirectional tape skin and a fabric stiffener. Finite element models using cohesive elements with input properties obtained from tape-to-fabric interface characterization tests were compared to experimental results of a stable skin-stiffener separation characterization test using a seven-point bend fixture. This fixture deformed the stiffened panel specimens into a buckled configuration which induced mixed-mode interlaminar stress states at the skin-stiffener interface. The advantage of this fixture was the potential for more stable damage initiation and delamination growth compared to a stringer-stiffened panel under axial compression. Use of the superposed cohesive elements showed promise, but the characterization of interface material properties as input to the cohesive elements remains a crucial component to be investigated to ensure accurate prediction of skin-stiffener separation

Stable Tearing and Buckling Responses of Unstiffened Aluminum Shells with Long Cracks

The results of an analytical and experimental study of the nonlinear response of thin, unstiffened, aluminum cylindrical shells with a long longitudinal crack are presented. The shells are analyzed with a nonlinear shell analysis code that accurately accounts for global and local structural response phenomena. Results are presented for internal pressure and for axial compression loads. The effect of initial crack length on the initiation of stable crack growth and unstable crack growth in typical shells subjected to internal pressure loads is predicted using geometrically nonlinear elastic-plastic finite element analyses and the crack-tip-opening angle (CTOA) fracture criterion. The results of these analyses and of the experiments indicate that the pressure required to initiate stable crack growth and unstable crack growth in a shell subjected to internal pressure loads decreases as the initial crack length increases. The effects of crack length on the prebuckling, buckling and postbuckling responses of typical shells subjected to axial compression loads are also described. For this loading condition, the crack length was not allowed to increase as the load was increased. The results of the analyses and of the experiments indicate that the initial buckling load and collapse load for a shell subjected to axial compression loads decrease as the initial crack length increases. Initial buckling causes general instability or collapse of a shell for shorter initial crack lengths. Initial buckling is a stable local response mode for longer initial crack lengths. This stable local buckling response is followed by a stable postbuckling response, which is followed by general or overall instability of the shell

Effect of Buckling Modes on the Fatigue Life and Damage Tolerance of Stiffened Structures

The postbuckling response and the collapse of composite specimens with a co-cured hat stringer are investigated experimentally and numerically. These specimens are designed to evaluate the postbuckling response and the effect of an embedded defect on the collapse load and the mode of failure. Tests performed using controlled conditions and detailed instrumentation demonstrate that the damage tolerance, fatigue life, and collapse loads are closely tied with the mode of the postbuckling deformation, which can be different between two nominally identical specimens. Modes that tend to open skin/stringer defects are the most damaging to the structure. However, skin/stringer bond defects can also propagate under shearing modes. In the proposed paper, the effects of initial shape imperfections on the postbuckling modes and the interaction between different postbuckling deformations and the propagation of skin/stringer bond defects under quasi-static or fatigue loads will be examined

Analysis Methods for Progressive Damage of Composite Structures

This document provides an overview of recent accomplishments and lessons learned in the development of general progressive damage analysis methods for predicting the residual strength and life of composite structures. These developments are described within their State-of-the-Art (SoA) context and the associated technology barriers. The emphasis of the authors is on developing these analysis tools for application at the structural level. Hence, modeling of damage progression is undertaken at the mesoscale, where the plies of a laminate are represented as a homogenous orthotropic continuum. The aim of the present effort is establish the ranges of validity of available models, to identify technology barriers, and to establish the foundations of the future investigation efforts. Such are the necessary steps towards accurate and robust simulations that can replace some of the expensive and time-consuming "building block" tests that are currently required for the design and certification of aerospace structures

Comparison of Damage Models for Predicting the Non-Linear Response of Laminates Under Matrix Dominated Loading Conditions

Five models for matrix damage in fiber reinforced laminates are evaluated for matrix-dominated loading conditions under plane stress and are compared both qualitatively and quantitatively. The emphasis of this study is on a comparison of the response of embedded plies subjected to a homogeneous stress state. Three of the models are specifically designed for modeling the non-linear response due to distributed matrix cracking under homogeneous loading, and also account for non-linear (shear) behavior prior to the onset of cracking. The remaining two models are localized damage models intended for predicting local failure at stress concentrations. The modeling approaches of distributed vs. localized cracking as well as the different formulations of damage initiation and damage progression are compared and discussed