Response of Multilayered Composite Laminates to Dynamic Surface
A theoretical investigation was performed of the response of.
SEARCH · Search NASA
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.
A theoretical investigation was performed of the response of.
To predict the mixed damage modes of unidirectional fiber-reinforced polymer (FRP) laminates under dynamic loading, an FEM-based peridynamic model is introduced in this paper. Based on its geometric structure and material composition, a long fiber lamina is considered a transversely isotropic medium as a result of homogenization at the meso-scale. The laminated structure is modeled by stacking surface mesh layers with arbitrary fiber angles along the thickness direction. The peridynamic bonds between Gauss points connect the separated elements. These bonds are classified as inner-layer bonds and inter-layer bonds. To represent the anisotropy of a laminate, the micro-elastic modulus of the inner-layer and inter-layer bonds is calculated from the anisotropic engineering material constants separately. To capture complex failure behaviors of laminate structures, an empirical damage model is proposed for the tension/ compression breakage of peridynamic bonds. This damage model can control the in-plane and delamination failure process. Finally, benchmark tests are conducted to validate the elastic response of laminates under dynamic loading. In terms of damage analysis, the proposed model can capture the complex damage modes and resistive force of laminate structures.
A unique dynamic delamination buckling and delamination propagation analysis capability has been developed and incorporated into a finite element computer program. This capability consists of the following: (1) a modification of the direct time integration solution sequence which provides a new analysis algorithm that can be used to predict delamination buckling in a laminate subjected to dynamic loading, and (2) a new method of modeling the composite laminate using plate bending elements and multipoint constraints. This computer program is used to predict both impact induced buckling in composite laminates with initial delaminations and the strain energy release rate due to extension of the delamination. It is shown that delaminations near the outer surface of a laminate are susceptible to local buckling and buckling-induced delamination propagation when the laminate is subjected to transverse impact loading. The capability now exists to predict the time at which the onset of dynamic delamination buckling occurs, the dynamic buckling mode shape, and the dynamic delamination strain energy release rate.
A unique dynamic delamination buckling and delamination propagation analysis capability has been developed and incorporated into a finite element computer program. This capability consists of the following: (1) a modification of the direct time integration solution sequence which provides a new analysis algorithm that can be used to predict delamination buckling in a laminate subjected to dynamic loading, and (2) a new method of modeling the composite laminate using plate bending elements and multipoint constraints. This computer program is used to predict both impact induced buckling in composite laminates with initial delaminations and the strain energy release rate due to extension of the delamination. It is shown that delaminations near the outer surface of a laminate are susceptible to local buckling and buckling-induced delamination propagation when the laminate is subjected to transverse impact loading. The capability now exists to predict the time at which the onset of dynamic delamination buckling occurs, the dynamic buckling mode shape, and the dynamic delamination strain energy release rate.
A material model has been developed which incorporates several key capabilities which have been identified as lacking in currently available composite impact models. The material model utilizes experimentally based tabulated input to define the evolution of plasticity and damage as opposed to specifying discrete input parameters (such as modulus and strength). The material model has been implemented into LS-DYNA® as MAT 213. The model can simulate the nonlinear deformation, damage and failure that takes place in a composite under dynamic loading conditions. The deformation model utilizes an orthotropic plasticity formulation. For the damage model, the nonlinear unloading response that is observed prior to the point where the peak stress is reached can be simulated, as well as the stress degradation response that occurs after the peak stress is reached. A variety of failure models, including a generalized tabulated failure model which facilitates the utilization of general failure surfaces, have been implemented into MAT 213. Recent studies at NASA have concentrated on using MAT 213 to analyze both the impact and crush response of a variety of laminated and textile architectures. Several of these studies will be discussed in the presented paper. For example, a woven carbon/Kevlar composite is currently being examined for use in an energy absorber system for rotorcraft structures. MAT 213 analyses have been conducted to examine the ability of the model to accurately simulate the dynamic crush response of this material. Studies are also being conducted to examine the ability of MAT 213 to simulate the ballistic impact response of representative laminated and woven thermoplastic and thermoset matrix composites. NASA efforts are also concentrated on developing methods and processes for improving characterization methods and developing “best practices” for using MAT 213, a summary of which will be discussed in the paper.
Predictive methods/computer codes for the computational simulation of acoustic fatigue resistance of hot composite structures subjected to acoustic excitation emanating from an adjacent vibrating component are discussed. Select codes developed over the past two decades at the NASA Lewis Research Center are used. The codes include computation of acoustic noise generated from a vibrating component, degradation in material properties of a composite laminate at use temperature, dynamic response of acoustically excited hot multilayered composite structure, degradation in the first ply strength of the excited structure due to acoustic loading, and acoustic fatigue resistance of the excited structure, including the propulsion environment. Effects of the laminate lay-up and environment on the acoustic fatigue life are evaluated. The results show that, by keeping the angled plies on the outer surface of the laminate, a substantial increase in the acoustic fatigue life is obtained. The effect of environment (temperature and moisture) is to relieve the residual stresses leading to an increase in the acoustic fatigue life of the excited panel.
This paper presents predictive methods/codes for computational simulation of acoustic fatigue resistance of hot composite structures subjected to acoustic excitation emanating from an adjacent vibrating component. Select codes developed over the past two decades at the NASA Lewis Research Center are used. The codes include computation of (1) acoustic noise generated from a vibrating component, (2) degradation in material properties of the composite laminate at use temperature, (3) dynamic response of acoustically excited hot multilayered composite structure, (4) degradation in the first-ply strength of the excited structure due to acoustic loading, and (5) acoustic fatigue resistance of the excited structure, including propulsion environment. Effects of the laminate lay-up and environment on the acoustic fatigue life are evaluated. The results show that, by keeping the angled plies on the outer surface of the laminate, a substantial increase in the acoustic fatigue life is obtained. The effect of environment (temperature and moisure) is to relieve the residual stresses leading to an increase in the acoustic fatigue life of the excited panel.
High-velocity impact tests were performed on graphite/epoxy laminates with embedded delaminations. High-speed photography of the transversely impacted laminates shows that dynamic local buckling of the sublaminate region formed by the delamination can cause a mode-I-dominated propagation of that delamination. Simultaneous recordings of the dynamic strain histories at several locations in the buckled sublaminate region are used to measure the laminate and sublaminate responses to the impact loading. Flexural motion resulting from transverse impact is shown to reduce local buckling of the sublaminate. This buckling causes a progressive separation of the delaminated surfaces and subsequent extension of the delamination in a mode-I dominated fashion.
Owing to the opaque nature of the laminated structures, traditional high-speed optical camera cannot be used to detect the dynamic process of sub-surface deformation. In this article, we report a study of using high speed X-ray imaging to study the high strain rate deformation in laminated Al structures. We used a Kolsky bar apparatus to apply dynamic compression and a high-speed synchrotron X-ray phase contrast imaging (PCI) setup to conduct the in situ X-ray imaging study. The in situ X-ray imaging captures the shock wave propagation in the laminated structures. After shock compression, we characterized the microstructures by using transmission electron microscopy (TEM), which demonstrates an increase of dislocation density. In conclusion, the micro-pillar compression tests show that the yield strength at 0.2% offset of laminated Al-graphene composite has a significant increase of 67%, from 30 to 50 MPa, compared to laminate Al after shock loading.
The damage response of resin-infused composites was investigated using quasi-static indentation (QSI) and instrumented drop-weight impact testing. The test specimens were composed of six layers of non-crimp fabrics where each layer contained +45°, 90°, −45°, and 0° plies. The goals of the testing were to identify damage responses that may differ from traditional thermoset materials and provide damage characterization data to evaluate current damage prediction analyses and aid improvements to current damage models. The damage was characterized using X-ray computed tomography (CT) and ultrasonic testing (UT). Incremental QSI loading was used on a single specimen to investigate damage growth in nine steps at progressively higher loads. Realtime acoustic emission data were used to determine the unloading points. Unlike typical thermoset composites, the force-displacement was smooth and showed no obvious indications of damage. The primary damage consisted of large lines of fiber fracture on the plies near the impact surface and delaminations through the thickness. Nine specimens were impacted at nominal impact energies of 25 ft-lbs, 30 ft-lbs, and 35 ft-lbs. The force histories at each energy followed nearly identical paths. Again, no obvious signs of damage were observed in any of the force histories. The general damage response in the dynamic tests was similar to the QSI-loaded specimens. However, repeat specimens impacted with a given impact energy showed variations in the damage. Some characteristics of the damage response were different than those typically observed in thermoset composites and will need to be evaluated using current analysis methods.
The damage response of resin-infused composites was investigated using quasi-static indentation (QSI) and instrumented drop-weight impact testing. The test specimens were composed of six layers of non-crimp fabrics where each layer contained +45°, 90°, −45°, and 0° plies. The goals of the testing were to identify damage responses that may differ from traditional thermoset materials and provide damage characterization data to evaluate current damage prediction analyses and aid improvements to current damage models. The damage was characterized using X-ray computed tomography (CT) and ultrasonic testing (UT). Incremental QSI loading was used on a single specimen to investigate damage growth in nine steps at progressively higher loads. Realtime acoustic emission data were used to determine the unloading points. Unlike typical thermoset composites, the force-displacement was smooth and showed no obvious indications of damage. The primary damage consisted of large lines of fiber fracture on the plies near the impact surface and delaminations through the thickness. Nine specimens were impacted at nominal impact energies of 25 ft-lbs, 30 ft-lbs, and 35 ft-lbs. The force histories at each energy followed nearly identical paths. Again, no obvious signs of damage were observed in any of the force histories. The general damage response in the dynamic tests was similar to the QSI-loaded specimens. However, repeat specimens impacted with a given impact energy showed variations in the damage. Some characteristics of the damage response were different than those typically observed in thermoset composites and will need to be evaluated using current analysis methods.
The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within commercial transient dynamic finite element codes, several features have been identified as being lacking in the currently available material models that could substantially enhance the predictive capability of the impact simulations. A specific desired feature pertains to the incorporation of both plasticity and damage within the material model. Another desired feature relates to using experimentally based tabulated stress-strain input to define the evolution of plasticity and damage as opposed to specifying discrete input properties (such as modulus and strength) and employing analytical functions to track the response of the material. To begin to address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed for implementation within the commercial code LS-DYNA. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. The effective plastic strain is computed by using the non-associative flow rule in combination with appropriate numerical methods. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used, in which a load in one direction results in a stiffness reduction in multiple coordinate directions. A specific laminated composite is examined to demonstrate the process of characterizing and analyzing the response of a composite using the developed model.
The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within commercial transient dynamic finite element codes, several features have been identified as being lacking in the currently available material models that could substantially enhance the predictive capability of the impact simulations. A specific desired feature pertains to the incorporation of both plasticity and damage within the material model. Another desired feature relates to using experimentally based tabulated stress-strain input to define the evolution of plasticity and damage as opposed to specifying discrete input properties (such as modulus and strength) and employing analytical functions to track the response of the material. To begin to address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed for implementation within the commercial code LS-DYNA. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. The effective plastic strain is computed by using the non-associative flow rule in combination with appropriate numerical methods. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used, in which a load in one direction results in a stiffness reduction in multiple coordinate directions. A specific laminated composite is examined to demonstrate the process of characterizing and analyzing the response of a composite using the developed model.
As the outlook for space exploration becomes more ambitious and spacecraft travel deeper into space than ever before, it is increasingly important that propulsion systems perform reliably within the space environment. The increased reliability compels designers to increase design margin at the expense of system mass, which contrasts with the need to limit vehicle mass to maximize payload. Such are the factors that motivate the integration of high specific strength composite materials in the construction of pressure vessels commonly referred to as composite overwrapped pressure vessels (COPV). The COPV consists of a metallic liner for the inner shell of the COPV that is stiff, negates fluid permeation and serves as the anchor for composite laminates or filaments, but the liner itself cannot contain the stresses from the pressurant it contains. The compo-site-fiber reinforced polymer (CFRP) is wound around the liner using a combination of hoop (circumferential) and helical orientations. Careful consideration of wrap orientation allows the composite to evenly bear structural loading and creates the COPV's characteristic high strength to weight ratio. As the CFRP overwrap carries most of the stresses induced by pressurization, damage to the overwrap can affect mission duration, mission success and potentially cause loss-of-vehicle/loss-of-crew. For this reason, it is critical to establish a fundamental understanding of the mechanisms involved in the failure of a stressed composite such as that of the COPV. One of the greatest external threats to the integrity of a spacecraft's COPV is an impact from the meteoroid and orbital debris environments (MMOD). These impacts, even from submillimeter particles, generate extremely high stress states in the CFRP that can damage numerous fibers. As a result of this possibility, initial assumptions in survivability analysis for some human-rated NASA space-craft have assumed that any alteration of the vessel due to impact is considered a catastrophic failure. This assumption is conservative and made due to lack of knowledge on the level of allow-able damage to the composite overwrap that can be sustained and still allow successful completion of the mission. To quantify the allowable damage level to the composite overwrap involves assessing stress redistribution following damage as well as evaluating possible time-dependent mechanisms involved in the COPV response to an impact event. Limited published work in this subject has shown that COPV can withstand at least some level of damage due to high energy impacts. These observations have been confirmed and expanded upon in recent experimental research performed by NASA. This research has demonstrated that there is not only robustness in a COPV to compensate for CFRP damage, but has also identified two significant failure modes for pressurized COPV. The lowest threshold failure mode involves the perforation of the vessel, and the highest threshold failure mode is the catastrophic rupture. While both of these failure modes mean a loss of the COPV, system robustness affords some tolerance to the venting as opposed to the more catastrophic rupture. As a consequence, it is necessary to understand the conditions that result in the transition between these failure modes. The aforementioned experimental research has been performed in both the unpressurized and pressurized condition to identify the damage level that triggered the failure thresh-old. This COPV test program was sponsored by the NASA Engineering and Safety Center (NESC), and tests were performed at NASA White Sands Test Facility (WSTF). Planning and coordination were provided by NASA JSC Hypervelocity Impact Technology (HVIT) group, and the COPVs were provided by the ISS Program. Unpressurized testing has been conducted at the pressure of the vacuum test chamber, while, the pressurized testing has been conducted at 290 +/- 10 bar (4,200 100 psi) using nitrogen as the pressurizing gas, which corresponds to the design pressure for the target COPV. In this research, spherical aluminum projectiles with varying diameter has been chosen as the impactor. For the unpressurized COPV, the dependence of penetration up to the dependence of hole size in the liner has been obtained as a function of impact conditions. For the pressurized research, the dependence of penetration up to rupture has been obtained as a function of im-pact conditions. Two representative post-test photographs of the failed COPV's from a nor-mal impact into the COPV surface are shown in Fig. 1. These images display the dramatic difference between failure modes, venting (Fig. 1a) and rupture (Fig. 1b). For venting, liner perforation, severed composite fibers/tows and ply delamination are commonly observed damage characteristics of this COPV failure mode. In the case of rupture, the COPV typically experienced a separation of its domed regions and severe break-up of the cylindrical region. Fully understanding the transition from venting to rupture experimentally is costly and potentially unachievable for conditions that cannot be generated in the laboratory. These shortcomings have motivated the performance of three-dimensional numerical simulations to expand the existing experimental database. These simulations have been carried out with the nonlinear-structural-dynamics, analysis-tool, CTH. A typical pressure contour plot from an impact simulation of an entire COPV is shown in Fig. 2. To generate the COPV stress state without initiating a shock wave, the pressure in the simulated COPV is ramped up to the final pressure over a millisecond prior to impact of the projectile with nitrogen gas. Figure 2a shows the system in this initial condition. After one millisecond, a projectile is initiated into the simulation and impacts the COPV. Figure 2b shows the system after this impact. In the figure, the onset of venting is represented as the change in pressure (μbar), red to green, at the perforation site. Also seen in the figure is the eroded projectile that had passed into the COPV vessel with the generated shock wave in the pressurant propagating just ahead of the material. In this paper, pertinent experimental details and the development of the material constitutive models necessary for this work along with the efforts to validate their use are dis-cussed. The simulation results are presented and compared with the NASA experimental observations. While work is on-going from this effort, early observations pertinent to the failure threshold are presented.