Modeling Progressive Damage In Laminated Composites
Report describes methodology for predicting strengths and service lifetimes of structural components made of laminate composite materials.
SEARCH · Engineering Papers
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.
Report describes methodology for predicting strengths and service lifetimes of structural components made of laminate composite materials.
The experiments described herein were performed to determine whether damage imposed by axial loading interacts with damage imposed by torsional loading. This paper is a follow on to a study that investigated effects of load-type sequencing on the cumulative fatigue behavior of a cobalt base superalloy, Haynes 188 at 538 C Both the current and the previous study were used to test the applicability of cumulative fatigue damage models to conditions where damage is imposed by different loading modes. In the previous study, axial and torsional two load level cumulative fatigue experiments were conducted, in varied combinations, with the low-cycle fatigue (high amplitude loading) applied first. In present study, the high-cycle fatigue (low amplitude loading) is applied initially. As in the previous study, four sequences (axial/axial, torsion/torsion, axial/torsion, and torsion/axial) of two load level cumulative fatigue experiments were performed. The amount of fatigue damage contributed by each of the imposed loads was estimated by both the Palmgren-Miner linear damage rule (LDR) and the non-linear damage curve approach (DCA). Life predictions for the various cumulative loading combinations are compared with experimental results.
In the present study, the assessment and evaluation of various acoustic tile designs were conducted using three-dimensional finite element analysis, which included static analysis, thermal analysis and modal analysis of integral and non-integral tile design options. Various benchmark specimens for acoustic tile designs, including CMC integral T-joint and notched CMC plate, were tested in both room and elevated temperature environment. Various candidate ceramic matrix composite materials were used in the numerical modeling and experimental study. The research effort in this program evolved from numerical modeling and concept design to a combined numerical analysis and experimental study. Many subjects associated with the design and performance of the acoustic tile in jet engine exhaust nozzle have been investigated.
An automated procedure has been developed to derive relevant factors, which can increase the ability to produce objective, repeatable methods for determining aerodynamic surface roughness. Aerodynamic surface roughness is used for many applications, like atmospheric dispersive models and wind-damage models. For this technique, existing lidar data was used that was originally collected for terrain analysis, and demonstrated that surface roughness values can be automatically derived, and then subsequently utilized in disaster-management and homeland security models. The developed lidar-processing algorithm effectively distinguishes buildings from trees and characterizes their size, density, orientation, and spacing (see figure); all of these variables are parameters that are required to calculate the estimated surface roughness for a specified area. By using this algorithm, aerodynamic surface roughness values in urban areas can then be extracted automatically. The user can also adjust the algorithm for local conditions and lidar characteristics, like summer/winter vegetation and dense/sparse lidar point spacing. Additionally, the user can also survey variations in surface roughness that occurs due to wind direction; for example, during a hurricane, when wind direction can change dramatically, this variable can be extremely significant. In its current state, the algorithm calculates an estimated surface roughness for a square kilometer area; techniques using the lidar data to calculate the surface roughness for a point, whereby only roughness elements that are upstream from the point of interest are used and the wind direction is a vital concern, are being investigated. This technological advancement will improve the reliability and accuracy of models that use and incorporate surface roughness.
The increased use of high-temperature composite materials in modern and next generation aircraft and spacecraft have led to the need for improved nondestructive evaluation and health monitoring techniques. Such technologies are desirable to improve quality control, damage detection, stress evaluation and temperature measurement capabilities. Novel eddy current sensors and sensor arrays, such as Meandering Winding Magnetometers (MWMs) have provided alternate or complimentary techniques to ultrasound and thermography for both nondestructive evaluation (NDE) and structural health monitoring (SHM). This includes imaging of composite material quality, damage detection and .the monitoring of fiber temperatures and multidirectional stresses. Historically, implementation of MWM technology for the inspection of the Space Shuttle Orbiter Reinforced Carbon-Carbon Composite (RCC) leading edge panels was developed by JENTEK Sensors and was subsequently transitioned by NASA as an operational pre and post flight in-situ inspection at the Kennedy Space Center. A manual scanner, which conformed'automatically to the curvature of the RCC panels was developed and used as a secondary technique if a defect was found during an infrared thermography screening, During a recent proof of concept study on composite overwrapped pressure vessels (COPV's), three different MWM sensors were tested at three orientations to demonstrate the ability of the technology to measure stresses at various fiber orientations and depths. These results showed excellent correlation with actual surface strain gage measurements. Recent advancements of this technology have been made applying MWM sensor technology for scanning COPVs for mechanical damage. This presentation will outline the recent advance~ in the MWM.technology and the development of MWM techniques for NDE and SHM of carbon wra~ped composite overwrapped pressure vessels (COPVs) including the measurement of internal stresses via a surface mounted sensor array. In addition, this paper will outline recent efforts to produce sensors capable of making real-time measurements at temperatures up to 850 C, and discuss previous results demonstrating capability to monitor carbon fiber temperature changes within a composite material.
A constitutive model has been developed for predicting the dynamic properties of composite beams which accounts for the influence of microstructural damage. The model assumes the damage state to be fixed. The model can account for both the damping increase and stiffness loss due to microstructural damage. Damage dependent material constants obtained from a baseline data set can be used to predict damping changes for several composite layups. The validity of the model has been verified by comparing the experimental results to model predictions for crossply laminates with a variety of stacking sequences.
Described here is the procedure for utilizing a damage constitutive model to predict progressive damage growth in laminated composites. In this model, the effects of the internal damage are represented by strain-like second order tensorial damage variables and enter the analysis through damage dependent ply level and laminate level constitutive equations. The growth of matrix cracks due to fatigue loading is predicted by an experimentally based damage evolutionary relationship. This model is incorporated into a computer code called FLAMSTR. This code is capable of predicting the constitutive response and matrix crack damage accumulation in fatigue loaded laminated composites. The structure and usage of FLAMSTR are presented along with sample input and output files to assist the code user. As an example problem, an analysis of crossply laminates subjected to two stage fatigue loading was conducted and the resulting damage accumulation and stress redistribution were examined to determine the effect of variations in fatigue load amplitude applied during the first stage of the load history. It was found that the model predicts a significant loading history effect on damage evolution.
A constitutive model was developed for predicting damping as a function of damage in continuous fiber reinforced laminated composites. The damage model is a continuum formulation, and uses internal state variables to quantify damage and its subsequent effect on material response. The model is sensitive to the stacking sequence of the laminate. Given appropriate baseline data from unidirectional material, and damping as a function of damage in one crossply laminate, damage can be predicted as a function of damage in other crossply laminates. Agreement between theory and experiment was quite good. A micromechanics model was also developed for examining the influence of damage on damping. This model explicitly includes crack surfaces. The model provides reasonable predictions of bending stiffness as a function of damage. Damping predictions are not in agreement with the experiment. This is thought to be a result of dissipation mechanisms such as friction, which are not presently included in the analysis.
The ability of a material model to capture in-plane matrix mode I and mode II crack growth is an essential component for modeling ply level damage evolution in composite structures. Previous studies using a continuum damage mechanics (CDM) approach have shown success in satisfying benchmark solutions for mode I and II crack growth. However, success was shown using a fiber-aligned meshing strategy, which encourages matrix cracks to propagate in a single band of elements, along the fiber direction. Generating a fiber-aligned mesh becomes a highly involved process for laminates including off-axis (non 0° or 90°) plies. The objective of this study is to quantify the effect of non-fiber aligned mesh discretization on predictions of in-plane matrix crack propagation. The approach taken incrementally varies the mesh orientation angle relative to the fiber orientation; more specifically, misaligned meshes are used to quantify the effect of element angle orientation relative to the initial crack orientation on the energy released during matrix crack propagation simulations using a CDM method. CDM solutions obtained with the misaligned meshes are evaluated against known benchmarks for mode I and II matrix crack growth. The CDM solutions reveal a near-polynomial trend of increased predicted failure stress with increased mesh misalignment angle; hence implying a potential relationship between element orientation angle and apparent fracture toughness.
The overall objective of this program is to demonstrate the applicability of NASA-developed advanced constitutive and life damage models for calculating cyclic structural response and crack initiation in selected components of reusable space propulsion systems. The computer model resulting from this program will enable the user to produce an accurate life prediction of hot gas path, life limiting components of propulsion systems such as the space shuttle main engine (SSME). Previously developed computer models addressing constitutive modeling and life damage will be combined in an advanced finite element analysis to generate a sophisticated baseline life prediction program. A material data base will be established for the constitutive and life models parametrically involving temperature, strain range, strain rate, mean strain/stress, and dwell time. The verified computer program will be used to accomplish the life predictions of three SSME critical components as evidence of the model functionality.
A dual spring model is developed for the static and fatigue damage predictions of laminates interface in composite structures. Stress concentrations can be induced by the defects formed in the fabrication or service process. A conventional S-N based fatigue damage model may not be accurate to predict the fatigue life of a structure with high stress concentration. With the dual spring model, static delamination failure can be simulated using springs of cohesive type material model while fatigue delamination development can be predicted using linear springs, where the crack driving force is computed based on virtual crack closure technique (VCCT). A Paris law type fatigue growth law with its mode mixity is applied for fatigue crack growth prediction. After verified using benchmark examples, including Double Cantilever Beam (DCB), End-Notched Flexure (ENF) and Mix-Mode Bending (MMB), the proposed dual spring model is applied in the static and fatigue damage prediction of NASA/Boeing sub-elements and UTC sub-elements.
The relationship between fatigue life and fatigue crack propagation rate is explored with a new cohesive damage model. The parameters of the model are obtained from idealizations of S-N diagrams used in engineering design. The model is based on the hypothesis that both stable tearing damage and damage due to cyclic loading are representations of a density of microcracks and, therefore, a single damage variable can describe the state of damage. This assumption implies that the quasi-static cohesive law that describes tearing is also the envelope of the fatigue damage. Fatigue damage within the cohesive envelope is assumed to accumulate at a rate that depends on the displacement jumps. The fatigue model was implemented as a UMAT subroutine for Abaqus cohesive elements by adding fatigue damage accumulation to a cohesive model based on the Turon mixed-mode cohesive laws. The analyses were conducted using a simplified cyclic loading procedure in which the maximum applied load is kept constant and the computational expense of cycling the load is avoided. The predicted propagation rates in mode I and mixed mode were compared to experimental results for IM7/8552 graphite/epoxy tape. Several aspects of the results were investigated, including the effect of R-curves, the stress ratio R, and the difference between displacement control and force control.
The effects of nuclear reactions on integral low-linear-energy-transfer (LET) protons spectra are studied, behind typical levels of spacecraft and body shielding, for the historically largest flares using the high-energy transport code BRYNTRN in conjunction with several biological damage models. The cellular track model of Katz provides an accurate description of cellular damage from heavy ion exposure. The track model is applied with BRYNTRN to provide a LET decomposition of survival and transformation rates for solar proton events. In addition, a fluence-based risk coefficient formalism is used to estimate Harderian gland-tumor induction in rodents and cataractogenesis in rabbits from solar flares, and a LET analysis is used to assess the relative contribution from target fragments on these biological endpoints.
The short-circuit current reduction in GaAs shallow junction heteroface solar cells was calculated according to a simplified solar cell damage model in which the nonuniformity of the damage as a function of penetration depth is treated explicitly. Although the equivalent electron fluence was not uniquely defined for low-energy monoenergetic proton exposure, an equivalent electron fluence is found for proton spectra characteristic of the space environment. The equivalent electron fluence ratio was calculated for a typical large solar flare event for which the proton spectrum is PHI(sub p)(E) = A/E(p/sq. cm) where E is in MeV. The equivalent fluence ratio is a function of the cover glass shield thickness or the corresponding cutoff energy E(sub c). In terms of the cutoff energy, the equivalent 1 MeV electron fluence ratio is r(sub p)(E sub c) = 10(9)/E(sub c)(1.8) where E(sub c) is in units of KeV.
The objective of this work was to understand matrix cracking mechanisms in a unidirectional alpha(sub 2) TMC in possible hypersonic applications. A (0)(sub 8) SCS-6/Ti-24Al-11Nb (at. percent) TMC was first subjected to a variety of simple isothermal and nonisothermal fatigue cycles to evaluate the damage mechanisms in simple conditions. A modified ascent mission cycle test was then performed to evaluate the combined effects of loading modes. This cycle mixes mechanical cycling at 150 and 483 C, sustained loads, and a slow thermal cycle to 815 C. At low cyclic stresses and strains more common in hypersonic applications, environment-assisted surface cracking limited fatigue resistance. This damage mechanism was most acute for out-of-phase nonisothermal cycles having extended cycle periods and the ascent mission cycle. A simple linear fraction damage model was employed to help understand this damage mechanism. Time-dependent environmental damage was found to strongly influence out-of-phase and mission life, with mechanical cycling damage due to the combination of external loading and CTE mismatch stresses playing a smaller role. The mechanical cycling and sustained loads in the mission cycle also had a smaller role.
This report provides a brief review of experimentally observed hereditary and nonhereditary material behavior along with background information on standard as well as advanced internal state variable constitutive modeling at elevated temperature. A description of exploratory, characterization, and validation testing is presented along with a detailed outline of what constitutes “sufficient” data content (i.e., quality and quantity) for developing or enhancing, characterizing, and validating a particular sophisticated nonlinear time- and history-dependent (hereditary) class of constitutive models known as GVIPS (generalized viscoplasticity with potential structure). The tests described are necessary to reveal a material’s behavior in the reversible (or viscoelastic) and irreversible (or viscoplastic) regimes, for the identification of both deformation and damage model parameters. Results presented are primarily for metallic materials. In all cases, both uniaxial and multiaxial tests are described, and the linkage between the specific tests and the parameters within the model that can be characterized from the results of these tests are also discussed. Discussion is also provided relative to the role information management must play relative to material data collection, analysis, maintenance, and dissemination. The need for such an information system is particularly important as both analyst and designer move toward utilizations of sophisticated, nonlinear time- and history dependent (hereditary) constitutive models. Lastly, the concept of state space and its utility in understanding and establishing constitutive models is addressed throughout. The intent behind this document is to help both the modeler and experimentalist understand each other’s specific points of view and provide guidance for both model development and characterization. Emphasis has been placed on providing the mechanician with information regarding how tests are performed and what issues to be aware of when interpreting results. It is hoped that experimentalists will take away a new perspective on the types of information that modelers and analysts are looking for from them.
Multiscale failure simulations have been performed for a Three-dimensional woven composite unit cell considering five, or more, length scales spanning the woven composite mesoscale to the sub-microscale voids. The multiscale recursive micromechanics approach, which enables recursive integration of general micromechanics theories over an arbitrary number of length scales, has been employed within the NASA Multiscale Analysis Tool. The multiscale model uses both the generalized method of cells and Mori-Tanaka micromechanics theories, and considers failure in the constituent materials using a simple damage model. Baseline results, containing distributed voids, are compared to uniaxial experimental data for an AS4 carbon fiber/ RTM6 epoxy matrix 3D orthogonal woven composite with good agreement in terms of global stiffness and global failure stress. The simulations demonstrate that the 3D woven composite exhibits damage tolerance through sustaining increasing axial load far beyond the first initiation of damage. The multiscale model is used to examine the nonlinear response of the material to other loading conditions. Case studies, motivated by X-ray computed tomography data, are presented on the effects of manufacturing induced voids and cracks.