Tuesday, June 10, 2008

Predicting Post-Buckling Response and Ultimate Failure of Composite 2-Stringer Panels


Software Suite for Durability, Damage Tolerance, Reliability & Life Prediction
Enhances FEA Solvers MSC Nastran*, ABAQUS, ANSYS, RADIOSS & LS-DYNA

* Best Performance and Verified Solutions with MSC Nastran



This Week's Feature Composite Example

Predicting Post-Buckling Response and
Ultimate Failure of Composite 2-Stringer Panels


Figure 1 - Damage & Durability Analysis of a Composite 2-Stringer Panel Subjected to Post-Buckling Analysis

Predicting the failure load of composite panels with post-buckled regions is a complicated undertaking as it involves buckling and damage evolution (Figure 1). The difficulty is complicated further by the fact that the ultimate failure (structural collapse) is driven by localized failures produced by damage in fiber and matrix in the post-buckled area of the structure. Reliable assessment of these load limits requires the application of advanced computational simulation that integrates composite mechanics at the micro-scale of fiber and matrix with finite elements, buckling, damage tracking and fracture analysis. The technical approach must rely on physics-based composite failure criteria capable of detecting all types of damages including non-visible ones. If localized buckling takes place before damage, the buckling mode shape is superimposed on the structure's geometry. Inception of damage prior to buckling requires repeated buckling analysis with degraded stiffness to account for detected damage. This analysis capability, integrated in the GENOA Software, was applied successfully to a two-stiffener composite panel under shear loading. The test results were closely reproduced with the analysis simulation. Including post-buckling effect was instrumental to replicate the test, since otherwise we risk over-predicting the failure load. As shown next, this analytical approach can be very effective in guiding the design and in reducing the number of tests of post-buckled composite panels.
 
Problem Description

The two-stringer J-stiffened composite panel evaluated for durability with post-buckling consideration is presented in Figure 2. It was made of AS-4/3501-6 carbon/epoxy unidirectional fabric. The fabric consisted of a wide sheet of unidirectional tows of fibers based together with polyester thread to keep it from unraveling. Twelve plies of the skin were laid up to form a 1.83 mm (0.072 in.) thick quasi-isotropic [0,90,45,0,-45,90]s laminate. The "J" shaped stiffeners were constructed with the same lay-up except for the flanges that were only half the thickness of the basic laminate. The skin panel and the stiffeners were stitched together with Kevlar and fiberglass threads [1,2]. The experimental results were adopted from the literature [1].
 
Figure 2 - Geometry of Two-stiffener composite panel [1]


The panel was modeled with Mindlin-Reissner thick shell elements and loaded in tension in the diagonal direction at one corner of the panel while fixing the opposite corner [2]. Details of the technical approach are provided next followed by discussion of results.

Technical Approach

Using a representative finite element model with appropriate stiffness and strength limits of the various materials, the software assesses the performance of composite panels with post-buckling consideration by:
  1. Determining the buckling mode shape of the panel (in this case skin buckling mode). Typically selects the buckling mode shape with the lowest eigenvalue.
  2. Initiating durability analysis through the sequential application of incremented static loading till structural fracture takes place. During the load stepping process, the software checks and tracks damage resulting from: fiber failure in tension or compression, matrix cracking (transverse tension/compression and in-plane shear), and delamination (normal tension and transverse and longitudinal out-of-plane shear). The analysis strategy is described here: 
  1. GENOA's progressive Failure Analysis (PFA) module starts with the full-scale finite element model and reduces the material properties down to the micro-scale of fiber and matrix. With every load step material properties are updated, reflecting any changes resulting from damage or cracks due to the applied loading. 
  2. If the applied load reaches the skin buckling load, the buckled shape is superimposed onto the panel. In this way the panel is placed into its bifurcated and stable lower energy state. In a physical test it would naturally move to this state. 
  3. If material damage is detected prior to local buckling, the code updates the stiffness and repeats the buckling analysis with reduced stiffness to extract a new buckling load. 
  4. When material and structural equilibrium states are reached, the load is incremented to the next level. With every additional load step, buckling analysis is repeated to update the mode shape. As the load increases, damage in the panel will initiate, grow and accumulate leading to ultimate panel failure, as shown next in the simulation results. 
Results
For the considered two-stiffener composite panel, experimental data indicated that local buckling of the skin occurred prior to reaching the ultimate load. Figure 3 shows the analytical results illustrating damage in the panel at initiation stage in the post-buckled region and damage progression when the ultimate load is reached. The software was used first to perform a linear bifurcation buckling analysis. The Buckling analysis showed that the initial buckling load is 8.275 kips (36.81 kN), which is very close to the experimental results [2]. As the load was incremented during the progressive failure analysis, the software superimposed the buckling mode shape on the structure's geometry when the applied load reached the buckling load of 8.275 kips. From that forward, the buckling analysis was repeated every time the load was increased followed by update in the structure's geometry to include the mode shape effect. Structural damage was detected at a load higher than the buckling load. The PFA process allows for localized failures to occur in the panel that ultimately results in panel collapse. Delamination due to relative rotation of the plies is shown in Figure 4.
(a) Damage Initiation                                         (b) Damage Progression

Figure 3 - Net damage due to progressive failure for composite two-stringer panel under diagonal tension. (Red indicates areas where failures have occurred)  a) at initiation of post-buckled region and b) near ultimate load. 

 
      
(a) Delamination Initiation                                         (b) Delamination Propagation
Figure 4 - Delamination location with respect to individual plies due to progressive failure analysis for composite two-stringer panel under diagonal tension. (Red indicates delamination locations)  a) at initiation of post-buckled region and b) near ultimate load. 

Load deflection results obtained from the simulation are compared to test and shown in Figure 5. The technical approach employed here yielded excellent agreement between the simulation and experimental results when progressive failure technique is used in a post-buckling analysis; i.e., when the panel is allowed to go into its lower energy buckled state prior to performing the progressive failure analysis. Figure5 also shows the diagonal deflection versus applied load obtained via simulating the progressive failure analysis without any buckling and post-buckling effects. The simulation results for the unbuckled panel clearly indicates that the ultimate load is over-predicted by about 35% and the stiffness by about 7.6% [2]. The results presented here demonstrate the effectiveness of the methodology integrated in GENOA. It can be reliably used to perform iterative designs prior to committing to a large test program. Testing can be done once satisfactory performance is derived from the software.
 
Figure 5 - Comparison of experimental and analytical load deflection curves for composite two-stringer panel [2]

 
References:
1. Yeh, H-Y and Chen, V., 1996. Experimental Study and Simple Failure Analysis of Stitched J-Stiffened Composite Shear Panels. Journal of Reinforced Plastics and Composites, Vol. 15, pp. 1070-1087.  Click here to email us for the technical publication.

2. Minnetyan, L. and Huang, D., 2001. Progressive Fracture of Stitched Stiffened Composite Shear Panels in the Postbuckling Range. Journal of Reinforced Plastics and Composites, Vol. 20, pp. 1617-1632. Click here to email us for the technical publication.


This issue was brought to you by Alpha STAR Corporation.

If you do not want to receive this newsletter, please click here.  If you still have trouble unsubscribing from our website, please email genoa@ascgenoa.com with a "UNSUBSCRIBE" subject header.

Monday, April 28, 2008

Composite Storage Module Joint Analysis and Test Verification


Software Suite for Durability, Damage Tolerance, Reliability & Life Prediction
Enhances FEA Solvers MSC Nastran*, ABAQUS, ANSYS, RADIOSS & LS-DYNA

* Best Performance and Verified Solutions with MSC Nastran



This Week's Feature Composite Example

Composite Storage Module Joint Analysis and Test Verification


Figure 1 -Multi-side Adhesive and Cohesive Failure of Joint and (Bottom) Composite Storage Module and Joint Geometry, Relatively Clean Ti Interface Indicates a Premature Adhesive Failure in this Location [1]

The premature failure of adhesively bonded composite joint tests, suggest adhesive failure as opposed to cohesive failure. This type of joint failure may be due to a) non-clean interface surfaces preparation, b) thick bond line, c) existence of voids in the adhesive bond line or a combination of all three

    In support of the ONR (Office of Naval Research) SBIR Phase II program, several Naval joint structures were studied. These joints included a General Dynamics designed joint for use on a composite storage module, and two test articles of a General Dynamics/Boeing designed joint attaching the walls of a director's room to the deck of a ship. During testing, the General Dynamics Composite Storage Module (CSM) joint failed prematurely and progressive failure analysis was employed to shed light on the possible causes of the premature failure.

    The General Dynamics Composite Storage Module (CSM) joint was a particularly useful example structure as a thorough analysis and testing program that accompanied the design and verification of the joint development. A significant portion of this program was developed into a paper [1]. The CSM joint version of interest was composed of a graphite/epoxy quasi-isotropic laminate adhesively bonded to a tapered titanium fitting. The CSM and joint details are shown in Figure 1. Details of the FE model and materials used in the analysis are shown in Figure 2. The tensile load on the joint bends the joint resulting in peel stresses in the adhesive EA9394 layer. The test indicated that the primary failure mode of the joint was adhesive failure within the EA9394 layer between the titanium and fiberglass.
 
Figure 2 -  FE Model Details and Materials for Progressive failure Analysis (PFA) [1]

    Several steps were taken before the final GENOA Progressive Failure Analysis (PFA) was performed in order to develop and demonstrate confidence in the software. For example, as a baseline check, linear analyses were performed using standard FEA software packages ANSYS and ABAQUS along with GENOA. General Dynamics performed the ABAQUS simulation. Stress results along the adhesive centerline of the EA9394 adhesive are shown in Figure 3 for all three codes. The validation of GENOA simulation to capture the nonlinear material behavior of the adhesive was achieved by duplicating the ABAQUS analysis performed by General Dynamics.
Figure 3 - Comparison of Centerline Peel Stress Predictions Within the (Linear) EA9394 [1]

    A full PFA with GENOA was then performed and compared to the test data. The predicted load deflection curve and associated damage mechanisms are shown in Figures 4a and b. While the failure mechanisms predicted by the analysis were very similar to those observed during the test, the analysis predicted a much higher strength for the joint as shown in Figure 4a. The tested joint may have failed at a much lower load (see blue curve in Figure 4a) because (1) the titanium surfaces were improperly prepared for the adhesive and/or (2) the maximum value of strain the adhesive could withstand was lower than that assumed from measurements of the adhesive done in 1995 and provided in reference 3. Without some relatively inexpensive coupon-level tests to establish a more relevant value of adhesive strength/strain limit, the failure strain criterion was not adjusted. However, improper surface preparation was investigated further.

    Subsequent review of the failure surfaces revealed a clean titanium surface indicating a premature interfacial failure of the adhesive to titanium bond (Figure 1). Analyses were performed to introduce the clean surface interfaces on the titanium fitting (test suggest adhesive failure as opposed to cohesive failure) using: 1) Virtual Crack Closure Technique (VCCT), and 2) using PFA and degrading the adhesive properties assuming using 20% void formation [2].

    VCCT required a predetermined crack path, and was modeled in two ways; a) multiple crack locations, and (b) single crack initiation point. For these simulations the crack path was determined from PFA (Figures 4b and 5a). The results obtained predicted higher strength value than predicted by PFA strain based analysis (blue curve in Figure 4a). Plane strain fracture toughness for peel and shear, (K1C and KIIC), values were required to run the analysis using the VCCT approach. The K1C and KIIC values for EA9394 adhesive were obtained from somewhere else [3]. Note that PFA predictions are usually close but are not considered very accurate when pre-cracked conditions exist in a simulation due to infinite stresses near a sharp crack tip.

 
Figure 4 - Load-deflection curve and associated damage events [2]

    Thereafter to capture the effect of improper bonding, the VCCT analysis was defined on limited node pairs as shown in Figure 5b. The trend of the results indicated that defining the partial bonding can simulate the effect of improper bonding, as shown in Figure 6.

    Similarly a void content of 20% was used to simulate the improper bonding in PFA. Again the results were found to be a close match with that of the test (see Figure 6).
Figure 5 - (a) Single VCCT approach. Note that the crack initiates from the right hand side and propagates to the left. (b) Partial Single VCCT approach to simulate partial (improper) bonding. The left hand side images are the initial setup and the right hand side is the final failure of the bonding [2]. 
 
Figure 6 - Comparison of the simulated load displacement curve with test data [2].

    In conclusion, if proper contact surface area information is available before the analysis, PFA and VCCT can be used judiciously to simulate the effect of improper bonding. Moreover, insight into potential improper bonding can can be achieved by varying the voids percentage in the adhesive to simulate the improper bonding scenario if surface area is not available.

References:

1. George F. Leon, Michael F. Trezza, Jeffrey C. Hall,1 and Kelli Bittick, "Evaluation of a Carbon Thermoplastic to Titanium Bonded Joint", ASTM-STP1455-11707.081503.  Click here to read technical publication.

2. Xie, D., Garg, M., Huang, D., and Abdi, F., "Cohesive Zone Model for Surface Cracks using Finite Element Analysis," AIAA-49SDM-106742-2008.  Click here to read technical publication.

3. T. R. Guess, E. D. Reedy, M. E. Stavig, 1995. "Mechanical Properties of Hysol EA-9394 Structural Adhesive," SANDIA REPORT, SAND95-0229. UC-704. Click here to read technical publication.

This issue was brought to you by Alpha STAR Corporation.

If you do not want to receive this newsletter, please click here.  If you still have trouble unsubscribing from our website, please email genoa@ascgenoa.com with a "UNSUBSCRIBE" subject header.

Monday, February 25, 2008

GENOA 4.3 Release with A- and B-Basis Allowables


Software Suite for Durability, Damage Tolerance, and Life Prediction
Augments FEA Solvers MSC Nastran*, ABAQUS, ANSYS & LS-DYNA

* Best Performance and Verified Solutions with MSC Nastran



This Week's Feature Composite Example

GENOA 4.3 Release Now Available 

Figure 1 - Progressive Failure Analysis of T-Joint

    The latest 4.3 release contains significant improvements in its popular Progressive Failure Analysis (PFA) module and translation support for NASTRAN, ABAQUS and ANSYS. Included in this release are 50 Step-by-Step Tutorial examples that fulfill a diverse range of applications and engineering needs.
    4.3 also features the new A- & B- Basis Allowable module which may be viewed through GENOA's Web Demo.  Users who are interested in applying their own data in this module may try GENOA through the on-line CVT web service for 3 free months.  For more information, please contactsales@ascgenoa.com.
    GENOA is an engineering software suite that augments commercial Finite Element Analysis (FEA) packages by providing Progressive Failure Analysis (PFA) capability. It can investigate structural responses to material degradation from damage induced by static, cyclic (low and high cycle fatigue), impact, Power Spectrum Density (PSD) and thermal loading.
 
Figure 2 - Over 100 materials in the new Archived Material Databank (Includes a large selection of commonly used fiber/resin and lamina properties) 

List of Significant Features and Updates

In response to the requests of our existing and potential customers, the following new features and improvements have been made in GENOA 4.3:

New Features
  • A- & B-Basis Allowables
  • Elliptic dome shapes
  • Strain Gage feature extended to both element- and node-based FE models
  • Thermal databank
  • Low and high fidelity filament winding optimization
  • Design envelope for Damage Plies % vs. Loading 
  • Archived databank of 100+ materials
Other Improvements
    GUI:
  • Results Model panel displays the FE model volume, weight and current iteration number next to GENOA logo
  • Importing of external models into Filament Winding module is made easier
  • Fixed concerns of numerical translation problems in Element Material Coordinate card due to Regional Settings Differences (GUI)
  • Windows VISTA support (Environment variables adjusted)
  • Improved table editing fields in all tables
  • Strain Gauge added for Stress and Strain panels (including Ply Stress and Strain panels)
  • Added the out-of-plane ply properties entry fields 
    PFA:
  • NASTRAN Element 35 (C-Gap element) added to GUI
  • Offset feature for shell elements added for both NASTRAN and ABAQUS
  • Extended PFA's capability to account for Virtual Crack Closure Technique (VCCT) analysis when mixed ply properties are used
  • Added *INCLUDE keyword to write user additional keywords for FE solvers used in GENOA
  • Added material definition for temperature variations
  • Evaluate inter-lamina strength for Ply Option
  • Added out-of-plane material properties entries for ply properties in material databank file
Figure 3 - New A&B Module (Effectively Reduces Coupon Testing for Material Qualification)

Did You Know?

Trying GENOA with Your Own Models

Alpha Star Corporation's Collaborative Virtual Testing (CVT) is a web service that runs GENOA through the Internet.  If you are interested in trying GENOA for free through the web with your own models, please contact Alpha Star Corporation for a free trial.  For more information, contact sales@ascgenoa.com.

This issue was brought to you by Alpha STAR Corporation.

If you do not want to receive this newsletter, please click here.  If you still have trouble unsubscribing from our website, please email genoa@ascgenoa.com with a "UNSUBSCRIBE" subject header.

Monday, January 28, 2008

A- and B-Basis Allowables for Composite Materials


Software Suite for Durability, Damage Tolerance, and Life Prediction
Augments FEA Solvers MSC Nastran*, ABAQUS, ANSYS & LS-DYNA

* Best Performance and Verified Solutions with MSC Nastran



This Week's Feature Composite Example

A Novel Approach to Determine A- and B-Basis Allowables for
Composite Materials 

Figure 1 - V-Notched Shear Test (Equivalent to ASTM D5379-Isopescu Test)

Many aircraft and spacecraft have advanced designs based on a well-matured and systematic testing technology; however, the manufacture and service life risks still remain. In production, testing to reduce risks due to manufacturing defects is performed at many steps up to and including the entire assembly / integration system levels. Even with all this testing, costly failures still occur and these are increasingly prone to public scrutiny. Table 1 shows the possible mission outcomes given a test prior to service. These tests could be at the part, sub-component, component, sub-assembly, assembly and product levels. As the system moves through production, the cost of executing the mitigation strategy increases due to the increased time and labor required for disassembly and reassembly. Composites are particularly vulnerable to increased risk from the execution of this mitigation strategy from both the additional testing that the adjacent structure will be subjected to and the risk of damage due to increased handling. Sometimes the additional damage done by repeated ground tests is not considered in the overall "remaining life" assessment of the system.
 
Table 1.  Possible Outcomes of the Item Tests [1]

A-Basis and B-Basis allowable values are of great importance in the engineering industries involving composite structures. In most cases, strength Allowables are defined using "design point strategy," meaning that Allowables will be obtained for certain predefined configurations, after which only these configurations are allowed for use in the design. In order to reduce cost, limited amount of configurations will be tested, reducing the choice in configurations and with that the optimization possibilities. A solution to this problem would be a "design space strategy" where few tests are needed to calibrate and validate a failure analysis method, after which analysis can be used to analyze a wide range of possible configurations.

The use of A-Basis and B-Basis allowable properties recognizes that material properties are statistical in nature. The two statistically based tolerance bounds are:
A-Basis or T99: At least 99% of the population of material values is expected to equal or exceed this tolerance bound [2-3] with 95% confidence (single point catastrophic failure with no-load redistribution)

B-Basis or T90: At least 90% of the population of material values is expected to equal or exceed this tolerance bound with 95% confidence (redundant load path with load redistribution)  
These allowable values are important for reducing risk in structural designs. Allowable determination is a time consuming and expensive process, since a large number of tests are required. In order to reduce cost and product lead-time, Virtual Testing is used to reduce necessary physical tests by replacing them with analysis. The objective is to reproduce scatter from multiple sources: manufacturing, material properties, and test. This is done analytically using combined multi-scale Progressive Failure Analysis and Probabilistic Analysis methods to generate reliably the scatter (variability) in material and coupon strength. To assure high accuracy, the multi-scale analysis is used based on a hierarchical analysis, where a combination of macro-mechanics and micro-mechanics are coupled with finite element analysis, damage tracking and fracture, and material degradation capability to analyze structures in great detail [4]: The procedure for predicting material Allowables is: 
1. Calibration of Constituents (fiber and matrix) to verify the lamina/laminate properties

2. Selection of Random Variables: appropriate variables are selected and reasonable assumptions are made to reproduce natural scatter in the material strength (i.e., void and fiber volume fractions). (See Reference 4 for further details.)

3. Use Virtual Simulation capability in GENOA Software to Predict Probabilistic Scatter in the Failure Load for Each Type of Test: combine progressive failure, as shown in Figure 1, analysis with probabilistic methods

4. Plot and Compare Generated Cumulative Distribution and Probability Density Functions for each test category, as shown in Figure 2.

5. Obtain A and B Basis Allowables from 0.01 and 0.1 Probabilities with 95% confidence, as shown in Figure 2. 
Figure 2. GENOA Results: Cumulative Distribution Function Vs. Test Data [3]

There are three distinct capabilities (options) for generating the A- and B-Basis Allowables using GENOA Software using:
1. User Provided Test Data - Approach combines deterministic and statistical methods as specified in: 
- MIL HDBK-17E: Military Handbook for Polymer Matrix Composites
- FAA CFR 14: Aeronautics and Space

2. Limited Test Data (1 coupon from each batch) - Approach combines progressive failure analysis with probabilistic methods. Reduce the number of tests for material qualification. 

3. Without Test Data - For trade study approach combines progressive failure analysis with probabilistic methods (assume scatter based on standard practices).
Figure 3. GENOA Probabilistic Sensitivities of Random Variables

Results & Discussions

Example of generating the A- and B-Basis Allowables using GENOA Virtual Simulation, based on limited test data [option 2] is shown in Figure 2. The material Allowables are generated and compared with published test data using 3 compressive fiber-resin woven composite coupon tests (1st coupon of the 1st panel of the 1st batch). In addition, the software generates the sensitivity of design parameters to the Allowables (Figure 3). The selected variables in simulation were: fiber longitudinal modulus, fiber compressive strength, matrix modulus, matrix compressive strength, fiber volume ratio and void volume ratio. The computational capability presented here results in three major benefits:
1. Reduction in the number of coupon tests for material qualification

2. Accelerate the trade-off studies in the down-selection process for certifying new materials

3. Identification of critical material and manufacturing variables to maximize material performance for a given application
References:

1. Frank Abdi, Tina Castillo, Edward Shroyer "Risk Management of Composite Structure" Book Chapter 45, CRC Handbook, January 2005. Click here to read technical publication. 

2. R. Rice, R. Randall, J. Bakuckas, S. Thompson., "Development of MMPDS Handbook Aircraft Design Allowables". Prepared for the 7th Joint DOD/FAA/NASA Conference on Aging Aircraft, September 8-11, 2003, New Orleans, LA. Click here to read technical publication.

3. DOT/FAA/AR-03/19, Final Report, "Material Qualification and Equivalency for Polymer Matrix Composite Material System: Updated Procedure" Office of Aviation Research, Washington, D.C. 20591, U.S. Department of Transportation Federal Aviation Administration, September, 2003. Click here to read technical publication.

4. M. R. Talagani, Z. Gurdal, and F. Abdi, S. Verhoef "Obtaining A-Basis and B-Basis Allowable Values for Open-Hole Specimens Using Virtual Testing" AIAAC-2007-127, 4. Ankara International Aerospace Conference, 10-12 September, 2007 - METU, Ankara. 
Click here to read technical publication.

This software feature will be available in our upcoming GENOA 4.3 release and utilized through our on-line web service.
This issue was brought to you by Alpha STAR Corporation.

If you do not want to receive this newsletter, please email genoa@ascgenoa.com with a "UNSUBSCRIBE" subject header.

Thursday, December 13, 2007

Composite T-Joint Design Analysis


Software Suite for Durability, Damage Tolerance, and Life Prediction
Augments FEA Solvers MSC Nastran*, ABAQUS, ANSYS & LS-DYNA

* Best Performance and Verified Solutions with MSC Nastran



This Week's Feature Composite Example

Composite T-Joint Design Analysis

Figure 1 - Configuration and laminate layups of the British Naval composite joint [1]. 

GENOA-PFA is virtual testing software tool simulating both the detailed micro and macro failures occurred in a composite structure throughout the entire loading process. It can greatly reduce the experimental investigation effort and cost involved in structural design. A British Naval composite joint [1, 2] was employed to demonstrate the application of GENOA-PFA in composite structure durability analysis and design. 
Figure 2 - Configuration and laminate layups of the British Naval composite joint [1]. 

Configuration of the T-Joint and the FEM -
 The joint is comprised of two FRP overlaminates bonded to either side of a web plate and then bonded to the base plate or flange (Figure 1). These overlaminates form a boundary angle connection and are comprised of alternating polyester/E-glass woven roving (WR) and chopped strand mat (CSM) layers. The gap within the boundary angle is filled with a compatible resin. The members being joined are comprised solely of polyester/E-glass woven roving. There are 14616 solid elements in the British joint model (Figure 2). The loading condition is three point bending. The load type is forced displacement. 

Simulation Results - The joint performance was simulated and the failure mechanism, which occurred in the joint, was identified. Figure 3 illustrates the comparison between the simulated and tested load-deflection relationship of the joint. The two results agree reasonably well for both damage initiation and final fracture loads.
Figure 3 - Comparison between simulated and tested load-deflection relationship.

The damage initiated as tensile- driven delamination in the fillets at the juncture area, where the interlamina tensile stress was the largest under the bending condition. Then the delamination failure propagated through the fillet thickness due to stress redistribution to the undamaged layers from the failed layers. The delamination in the fillets reduced their contribution to bending resistance, which finally resulted in the fracture of the flange at the end of the loading process. The entire simulated damage and fracture process of the joint is illustrated in Figure 4 where red areas represent the damage.

Figure 4 - Simulated failure process of the British Naval joint. Red areas represent the failure which was caused by interlamina tensile stress (delamination).

Conclusions
GENOA-PFA computed detailed laminate failure in the composite joint throughout the entire loading process. The load capacity of the British Naval joint was accurately predicted and its underlying failure mechanism was clearly identified, namely, delamination due to interlamina tensile stresses. Hence, GENOA-PFA is a useful virtual testing tool for optimal design of composite structures, e.g. for a composite joints, the fillet radii and thickness can be optimized to reduce the delamination failure and thus maximize the joint durability. 
References:
1. Cody Godines, Frank Abdi, Steven Kiefer and Keith Kedward, "Simplified Analytical Procedure for Prediction of Fracture Damage in Composite Structures", ASTM COMMITTEE-D30 Symposium on Joining and Repair of Composite Structure March 17-18, 2003 Kansas City, MO. Click here to read technical publication.
2. Phillips, H.J., and Shenoi, R.A., "Damage Tolerance of Laminated Tee Joints in FRP Structures", Composites Part A - Applied Science and Manufacturing, Vol. 29, No. 4, pp. 465, 1998. Click here to read technical publication.
 

Did You Know?

Probabilistic Progressive Failure Analysis 

imageGENOA's Probabilistic Progressive Failure Analysis capability enables the prediction of structural reliability in presence of uncertainties in fabrication parameters, cure, material, geometry, and loading. First, perform low fidelity probabilistic analysis to identify influential random design variables. Second, reduce the list of variables to include critical ones and perform high fidelity simulation (e.g. Monte Carlo) to obtain a measure of reliability. As additional benefits, you will obtain a database of competing designs to improve the product performance and reduce the number of unnecessary tests.  For more information on this feature and trying out GENOA through our demos, please contact info@ascgenoa.com.
This issue was brought to you by Alpha STAR Corporation.

If you do not want to receive this newsletter, please email genoa@ascgenoa.com with a "UNSUBSCRIBE" subject header.

Tuesday, November 13, 2007

Delamination Initiation/Propagation Failure Analysis of Reinforced Carbon-Carbon Woven Composite Specimens


Software Suite for Durability, Damage Tolerance, and Life Prediction
Augments FEA Solvers MSC Nastran*, ABAQUS, ANSYS & LS-DYNA

* Best Performance and Verified Solutions with MSC Nastran



This Week's Feature Composite Example

Delamination Initiation/Propagation Failure Analysis of Reinforced Carbon-Carbon Woven Composite Specimens

Figure 1 - Space Shuttle and Close-Up View of the RCC Panel near the Leading Edge [1]. 

The Shuttle Reinforced Carbon-Carbon (RCC) leading edge (Figure 1) is a brittle composite system which is subject to cracking and delamination from foreign debris and other impacts. The combination of cracks and delaminations can provide pathways for hot gases to enter the interior of the Shuttle wing during reentry, leading to serious consequences. If an impact occurred during launch, delaminations would be unseen by surface inspection on orbit. An investigation was carried out to demonstrate delamination prediction in RCC when impacts are known to have occurred. 
Figure 2 - Three-Point Bending (a) Test Setup and (b) Schematic and Dimension of the Short Specimen [1].

A detailed failure analysis of two RCC specimens in a three-point bending configuration (Figure 2) was carried out by Alpha STAR Corporation (ASC) using GENOA and compared to NASA tests. GENOA progressive failure analysis (PFA) for 1.0 inch and 1.5 inch three-point bending specimens revealed that the twospecimens exhibit different failure behavior prior to final failure (different delamination initiation and damage/fracture progression). 
Figure 3 - Test (a) versus Predicted [1] (b): Major and Minor Delamination (Initiation & Propagation) Location. Red Indicates Damaged Locations Just Prior to Final Failure. Simulation was based on a Single Layer of Shell Elements as the Finite Element Model.

The finite element models of each specimen was constructed with 200 Mindlin-Reissner shell elements, simulating woven Reinforced Carbon-Carbon (RCC) material properties. Shell thickness was 0.229 inch. The loading was applied using displacement control. Each shell element consisted of 19 plies, significantly reducing the computational time without loss of analysis details such as stresses, strains, and damage information for the individual plies. The material properties were initially calibrated with the 1.0 inch specimen test data.  The material calibration process was presented in one of our previous issues.

The predicted simulation results for the 1 and 1.5 inch specimens compared well with the reported test data (Figures 3 to 5). Alpha Star Corporation (ASC) did not have access to the 1.5 inch long specimen test results prior to submission of the analytic predictions to NASA (Figure 5).

Figure 4 - Calibrated: Load-Displacement Response of 1.0 Inch Three-Point Bending Specimen [1].

Figure 4 shows the test load-displacement behaviour and the corresponding progrssive failure analysis predictions for the three-point RCC bending specimens. Points A through F indicate the onset of various types of damage.


Figure 5 - Measured and Predicted Load-Displacement Behavior for 1.5 Inch Specimen [1].

Progressive Failure Analysis predicted failure loads of 350 lbf and 231 lbf for the 1 inch and the 1.5 inch specimens, respectively (Figures 4 and 5). Measured and predicted failure loads were in excellent agreement, with errors of less than 1% (1 inch specimen) and 5.7% (1.5 inch specimen). Both the shorter and longer specimens were predicted to fail at the specimen mid-section.This was confirmed by the test results.The analyses indicated that the failure mode of the shorter specimen was very different from that of the longer specimen. PFA of the shorter specimen predicted considerable delamination, but longer specimen PFA revealed virtually no delamination. In the shorter specimen, two delaminations were predicted (Points C, D, and E in Figure 4); a primary one and a secondary one.

Figure 3 shows the RCC damage when the load reaches point F (348 lbf) inFigure 4. The delamination, which initiated at load point C and is tracked as excessive relative rotation at ply 3, that is, between plies 2 and 3, and 3 and 4, grows toward the Y-axis (green arrow) as the applied load is further increased (Figure 3). It grows to become a secondary delamination of the specimen.

Also at an applied load of 348 lbf additional delamination was observed near the center of the ply stack, ply 9 (Figure 3). This new delamination occurs between plies 8 and 9, and 9 and 10 and is due to excessive transverse normal shear. The delamination originates where the shear is greatest, namely under the applied load and grows to become the primary delamination in the specimen. It reduces the bending stiffness of the specimen in the vicinity of the delamination.

Figure 6 - Photograph of damaged 1.5 inch specimen: (a) Test (b) Prediction at 231 lbf [1].

Photos of the interior of the failed test specimens confirmed the predicted failure behavior. No delamination occurred in the longer specimen (Figure 6), but primary and secondary delaminations occurred in the shorter specimen (Figure 3). Figure 6 shows the similarity between the predicted and the snaps taken at the failure. The simulation results are at 231 lbf.

Conclusions
GENOA PFA predictions confirmed that the 1-inch three-point bending specimen will exhibit two delaminations; a primary one and a secondary one. The primary one results in considerable loss of stiffness prior to failure, whereas the secondary one does not. Photographs of the interior damage in the tested specimens indicated that the two delaminations occur near their predicted locations in the 1-inch specimen.

Confidence in the analytical predictions was reinforced since the predicted load-displacement behavior of both specimens were in excellent agreement with the measured load-displacement behavior. 

References:
[1] Sokolinsky, V. S., Housner, J., Surdenas, J., and Abdi, F., 2006. "Progressive Failure Analysis of Shuttle Reinforced Carbon-Carbon Plate Specimens." AIAA-2006-1789, RI, May 1-5. Click here to read technical publication.
 

Did You Know?

GENOA Custom Failure Criteria

imageGENOA has a recommended 'Custom Failure Criteria' that consists of a set of well known failure criteria, such as Tsai Hill, Puck, Strain Invariant and many more that help predict the strength of composites made up of wide variety of materials and several configuration (2D, 3D braid and woven). The approach is simple to use and allows implication of user defined failure criteria subroutine.  For more information on this feature and trying out GENOA through our demos, please contact info@ascgenoa.com.
This issue was brought to you by Alpha STAR Corporation.

If you do not want to receive this newsletter, please email genoa@ascgenoa.com with a "UNSUBSCRIBE" subject header.

Thursday, October 18, 2007

Prediction and Verification of Metal Fracture Toughness Tests Using Non-linear Static Stress-Strain Curve


Software Suite for Durability, Damage Tolerance, and Life Prediction
Augments FEA Solvers MSC Nastran*, ABAQUS, ANSYS & LS-DYNA

* Best Performance and Verified Solutions with MSC Nastran



This Week's Feature Metal Example

Prediction & Verification of Metal Fracture Toughness Tests Using Non-linear Static Stress-Strain Curve

Figure 1 -  Full Stress-Strain Curve and Crack Tip Deformation [2] (a) Areas Associated With The Uniform (b) A Center Crack in Wide Panel and Non-Uniform Straining 
Safe life prediction of components must be conducted to ensure the life adequacy of parts during service usage. In many cases fracture properties of material are not available because of 1) cost associated with generating fatigue and fracture data, 2) Inability to conduct tests because of time limitation and deadline set forth by the customers, and 3) lack of analytical tools to conduct a comprehensive crack tip stress analysis. 

New material-physics based computational methodologies for assessing the plane-stress and the plane-strain fracture toughness (KC, KIC) and (da/dN versus Delta K curve) of the material have been used successfully in predicting these fracture allowables, knowing only the complete stress-strain behavior of the material of interest. 

Figure 2 - Typical Crack Growth Rate Versus Stress Intensity Range
Farahmand extended the Griffith theory to estimate fracture toughness value of metals from simple uniaxial tensile tests. Figure 1 illustrates the extended Griffith theory and regions of crack tip plastic deformation. Accounting for the energy absorption rate for plastic deformation at the crack tip is calculated and used to establish a relationship between fracture stress and half critical crack length [1].
On the other hand, fatigue crack growth properties of the material are determined using the well known Newman Forman and Koening (FNK) equation requiring input from fracture toughness theoretical model. The analytical procedure relies on both implicit and explicit computational schemes and evaluates points in the threshold, Paris, and accelerated regions (Figure 2).

Once the fracture allowables are determined, the values can be further used to predict the S-N curve of the component using finite element method approach with Virtual Crack Closure Technique (VCCT). The formulation extends to fatigue crack growth and strength life prediction of notched and unnotched components. Scientists at Alpha STAR and TU Delft are extending the algorithm to composites.
Figure 3 - Process and Comparison of Fracture Toughness Versus Thickness [FTD] and da/dN Versus Delta K Curve [FCG] With The Test Data Provided in NASGRO for Ti-6Al-4V (Mill Annealed) [2].

Figure 3 shows the process and comparison of the results obtained using the two methodologies (Fracture Toughness Determination [FTD] and Fatigue Crack Growth [FCG]) for a Titanium alloy. Similar verification has been done with several other pure and alloyed materials, such as Aluminum alloy, Inconel, Steel, and many more.

Figure 4 - Metallic Center Cracked Panel Subjected to Quasi-Static Fatigue Loading [2]
The virtual testing technique was later used to generate the high cycle fatigue data (the S-N curve). The fracture allowables predicted from the FTD/FCG modules in GENOA for a center-cracked specimen were used to assess the total life of uncracked specimen made of Ti6-4MA and 7075-T6 Titanium and Aluminum alloys, respectively (Figure 4). The panels were subjected to quasi-static fatigue loading using progressive failure analysis in conjunction with Virtual Crack Closure Technique (VCCT) (A link to the past news letter). The comparison of the predicted and test results is tabulated in Table 1 for four tests.

The capability of the three step approach: 1) Fracture Toughness Determination, 2) Fatigue Crack Growth Determination, and S-N curve Determination was demonstrated Life Assessment of Boeing 747 crown Panel Fuselage Section (Figure 5).

Table 1 - Comparison between Test and Simulation Results for Life Assessment of Metallic Center Cracked Panel, As Shown In Figure 3 [2]
The results validated that the novel Fracture Toughness Determination, Fatigue Crack Growth Determination and Life Assessment Methodologies in GENOA indicating further that GENOA can be reliably used to assess fracture allowables for a materials extremely quickly and with reasonable accuracies.
Figure 5. Three Step Approach to Assess Life of the Stiffened Curved Panel Made of Aluminum Alloy [2]
In addition, a probabilistic analysis of the fracture toughness and fatigue crack growth can also be performed using the Probabilistic Fracture Toughness (PFTD) and Probabilistic Fatigue Crack Growth (PFCG) modules in GENOA (Figures 6 & 7). The probabilistic capability allows monitoring the sensitivity of the response (KC, KIC, and da/dN versus DK curve) to different input variables [3].
Figure 6 - Variation of Fracture Toughness with Variable Thickness [3]
Figure 7 - Variation of Plane Stress (KC) and Threshold (Kth) Fracture Toughness Due to Variations in Material Properties. [3]

Click here to receive demo and presentation of FTD & FCG.

References:
1. Farahmand, B., Fatigue and Fracture Mechanics of High Risk Parts, Chapman and Hall, 1997.
2. Farahmand, B., Saff, C., Xie, D., and Abdi, F., 2007. Estimation of Fatigue and Fracture Allowables for Metallic Materials Under Cyclic Loading. AIAA-2007-2381.Click here to read technical publication.
3. Farahmand, B., and Abdi, F., 2002. Probabilistic Fracture Toughness, Fatigue Crack Growth Estimation Resulting From Material Uncertainties. ASTM International Paper. Click here to read technical publication.
 

Did You Know?

Benefits from Virtual Simulation of ASTM Tests

imageVirtual simulation of ASTM tests using GENOA enables the qualification and characterization of aerospace materials. The successful replication of these tests provides the designer and analyst with a reliable tool to evaluate the component performance. Iterative designs can be made with GENOA until satisfactory performance is achieved. This capability eliminates redundant tests thereby expediting the component certification and delivery to market.  For more information on this feature and trying out GENOA through our demos, please contact our sales at sales@ascgenoa.com.
This issue was brought to you by Alpha STAR Corporation.

If you do not want to receive this newsletter, please email genoa@ascgenoa.com with a "UNSUBSCRIBE" subject header.