Home Contact Sitemap

DTD Handbook

Handbook for Damage Tolerant Design

  • DTDHandbook
    • About
    • Contact
    • Contributors
    • PDF Versions
    • Related Links
    • Sections
    • Examples
      • 1. FAC Problems
      • 2. Merc Problems
      • 3. NRC Problems
      • 4. SIE Problems
        • 0. Crack Growth Analysis - Main Cargo Door (MSD Cracks)
        • 1. Damage Tolerance Analysis - Wing Main Spars
        • 2. Crack Growth Analysis - Main Cargo Door (Primary and Continuing Damage)
        • 3. Damage Tolerance Analysis - Rear Wing Spar
        • 4. Crack Growth and Residual Strength Analyses - Lap Joint
      • 5. UDRI Problems

Problem SIE-3

Title: Crack Growth Analysis of Main Cargo Door Surround Doubler Attachment to Fuselage Skin with Primary and Continuing Damage Cracks


To illustrate the process of estimating crack growth behavior to set inspection limits.

General Description:

This problem focuses on a damage tolerance assessment of a main cargo door surround doubler attachment to fuselage structure for the purpose of establishing inspection intervals for crack growth with primary and continuing damage cracks growing from opposite sides of a hole.  The critical area includes the main cargo door surround doubler and the existing fuselage skin.  The stresses acting at the doubler attachment at the upper edge are derived from a conservative loading spectrum based on pressure loading.  The critical area was modeled using a standard NASGRO 3.0 stress intensity factor solution and crack growth model. 

Topics Covered:          Damage tolerance assessment, crack growth analysis, inspection intervals

Type of Structure:       fuselage skin, main cargo door surround doubler

Relevant Sections of Handbook:  Sections 2, 5, 11

Author:                     Lesley Camblin

Company Name:   Structural Integrity Engineering
                                    9525 Vassar
                                    Chatsworth, CA  91311

Contact Point:       Matthew Creager

Phone:       818-718-2195

e-Mail:      mcreager@sieinc.com


Overview of Problem Description

This problem focuses on the main cargo surround doubler attachment to the existing fuselage skin at stringers 2R and 26L.  The skin is considered to be a single load path structure under the total hoop stress before the doubler attachment.  The critical location is in the skin at the first row of fasteners because the skin sees both bypass and bearing stresses at this row, where as, at the other fastener rows the load is in both the doubler and the skin with each row having lower load transfer. 

The fuselage skin was fabricated from 2024-T3 aluminum.  The fasteners are 0.188 in diameter, and join the skin and surround doubler. 

The specific area is shown in View A of Figure SIE-3.1, with the specific details and the primary and continuing damage cracks shown in Figure SIE-3.3.  Note that the skin at this first row of fasteners is a single load path as shown in Figure SIE-3.2.


Figure SIE-3.1.  Main Cargo Door Doubler Installation



Figure SIE-3.2.  Structural Detail for Critical Area


Figure SIE-3.3.  Detail Geometry of Critical Location, View A.


Model Geometry Description

The crack growth analysis is based on the Fatigue Crack Growth Computer Program NASGRO3.0.  This computer program calculates crack growth for a single crack for several standard crack cases.  Crack growth rate calculations use the “NASGRO” equation with elements developed by Forman, Newman, de Koning, and Henriksen (see NASGRO reference manual).  This is a modified Paris equation to account for fatigue crack closure, stress ratio effects, and upper and lower fatigue crack growth rate asymptotes for threshold and critical crack growth. 

The analysis uses the NASGRO3.0 material libraries for the crack growth rate equation constants.  Non-interaction of loads and constants for the Forman crack growth rate equation are used.

Since the standard crack models in NASGRO3.0 are for crack growth of single cracks, no influence of one crack upon another is calculated in NASGRO3.0 for these standard cases.  The analysis presented here includes crack interaction effects between the primary damage crack and the continuing damage crack.  This is accomplished by iterating though a series of NASGRO3.0 computer runs tracking the growth of both cracks and modifying the stress intensity factors appropriately.  The increased stress intensity factors are based on the crack sizes of the interacting cracks from the previous iteration and correction factors based on the compounding of analytical stress intensity solutions.

This iteration procedure is accomplished in an Excelâ Spreadsheet utilizing Visual Basic Programming to submit a NASGRO3.0 computer run for each crack at each iteration.  The spreadsheet reads the NASGRO3.0 output files for cycles and current crack lengths.  Based on these crack lengths, correction factors are calculated and input into the NASGRO3.0 input file for the next iteration, which is automatically submitted by the spreadsheet. 

The correction factors are accounted for by increasing the stress scaling factors input into NASGRO3.0.  These increased stress scaling factors can be input based on the following:

These correction factors for crack interaction account for interactions between the primary damage crack (rogue flaw) and the continuing damage crack.  This is done assuming both cracks are in the same part.

Note that interactions between the rogue flaw and the continuing damage crack have historically not been done in crack growth analysis.  This method of including these interactions from the onset of the crack growth of the rogue flaw is conservative since it does not account for any fatigue life due to the nucleation of the continuing damage crack.

As previously discussed, the correction factors for crack interaction are based on comparison of analytical stress intensity solutions.  The correction factor for this analysis is termed “INT”, and is used for unequal length cracks growing from opposite sides of a hole.

Figure SIE-3.4. INT correction factor


These correction factors are based on through the thickness cracks.  They are used for part through cracks when defined with an equivalent crack length.  The equivalent crack length is based on equating the area of a part through crack as a quarter ellipse to that of an equivalent through crack as a rectangular area with thickness, t:

The “INT” correction factor is derived based on comparing the stress intensity solution of a center cracked panel for two different crack lengths, a1 and a2.  Including the diameter of the hole, D, in the total crack lengths, yields:

The crack growth model for the main cargo door surround doubler attachment to the fuselage skin at stringer 2R employs the NASGRO3.0 corner crack from a hole centered in a plate, CC02, with the correction factors for the influence of unequal length cracks growing from opposites sides of a hole.  When the initial crack reaches the edge of the plate, the crack growth is continued in the opposite direction as a through crack from the edge of a plate, TC02.

The crack growth model CC02 was used with the following dimensional values. 


Figure SIE-3.5.  NASGRO3.0 Crack Model, CC02.

Two NASGRO files are created for the primary damage crack and the continuing damage crack and submitted to the Excel interaction spreadsheet.  The spreadsheet accesses NASGRO and grows both cracks for 100 flights.  The b correction factors are calculated for the crack lengths at that time and the resulting increased stress scaling factors are plugged back into the NASGRO files.  The interaction spreadsheet grows the two cracks until there is a 10% change in the primary damage crack length (this could also be done in increments of flights), recalculates the b correction and stress scaling factors, and continues to grow the cracks until the primary damage crack reaches the edge of the plate.

Once the primary damage crack reaches the edge and transitions into an edge crack growing in the opposite direction, the crack growth model TC02 was used with the following dimensional values.

Figure SIE-3.6. NASGRO3.0 Crack Model, TC02

Inspection Capabilities and Crack Limits

The holes in the fuselage skin at the attachment of the first row of the surround doubler attachment at 2R (and 26L) are directly accessible from the inside.  Therefore, these areas are inspected by HFEC surface probe.  With a HFEC inspection, the minimum detectable crack size in the field is assumed to be a 0.0625 inch crack past the fastener head.

Structural Loading and Stress History Description

The stress spectrum is considered to have a remote stress due to cabin pressurization.  Cabin pressurization primarily causes hoop tension in the fuselage.  The GAG pressurization load is based on FAR25.571. The pressure condition is comprised of a 7.8 psi normal operating differential pressure and an additional 0.5 psi external aerodynamic pressure.  A factor of 1.1 is only applied to the normal operating pressure for residual strength.

The bypass and bearing load at the critical fastener row is calculated using a displacement compatibility analysis as described by Swift (“Repairs to Damage Tolerant Aircraft,” presented to the International Symposium on Structural Integrity of Aging Airplanes, Atlanta, Georgia, USA, 1990).  Layer “a” is the fuselage skin and an existing bonded doubler.  Layer “b” is the main cargo surround doubler.  The surround doubler becomes fully effective after the first three rows.  This analysis shows the most critical fastener location is the first row of fasteners.

Table SIE-3.1.  Fastener Transfer Calculations.

Based on these results, 30% of the load is taken through bearing in the first row of fasteners.  This first row of fasteners therefore has 30% as a bearing load and the remaining 70% as a bypass load.

The axial stress and bearing stress acting on this section are:

The limit stress used for residual strength purposes in this scenario is calculated, as stated earlier, according to FAR25.571.

The residual strength axial stress and bearing stress acting on this section are:

Material Property Description

The outer skin and doubler are made from 2024-T3 IAW QQ-A-250/5.  The material properties from the NASGRO3.0 libraries are used for the fracture toughness and the crack growth rate properties.  The material properties used are for 2024-T3; Clad, Plate and Sheet; T-L; LA & HHA NASGRO material code M2EA12AB1.

Table SIE-3.2.  Material Properties and Growth Rate Data.


          MATL 1: 2024-T3                                                      

                  Clad Plt & Sht; L-T; LA & HHA                                


          Material Properties:                                                 


          :Matl:  UTS :  YS  :  K1e :  K1c :  Ak :  Bk :  Thk  :  Kc  : Keac : 

          : No.:      :      :      :      :     :     :       :      :      : 


          :  1 :  66.0:  53.0:  46.0:  33.0: 1.00: 1.00:  0.036:  66.0:      : 


          :Matl:--------------- Crack Growth Eqn Constants -------------------:

          : No.:    C    :  n  :  p :  q :  DKo : Cth+ :Cth- : Rcl:Alpha:Smax/:

          :    :         :     :    :    :      :      :     :    :     :SIGo :


          :  1 :0.829D-08:3.284:0.50:1.00:  2.90:  1.50: 0.10:0.70: 1.50: 0.30:


Solution Technique

This type of problem is conveniently solved using NASGRO3.0 with the crack growth interactions previously discussed.  The input files for the equal length cracks growing from opposites sides of a hole are identical for the NASGRO3.0 analysis shown in Table SIE-3.3.  The spectrum is included as a constant amplitude GAG cycle with 100 flights per block, with a single block applied per schedule.


Table SIE-3.3.  NASGRO Input File for Problem SIE-3.




Output file name


1=US units



71fc1-2 skin at upper and lower doubler edges

Problem name


Crack model type


Crack model no.


Thk, t






Hole center to edge


Poisson's Ratio


U=User defined crack


Initial a


Initial a/c


Number of materials


Non Interaction


Matl input choice


File input choice


Material Category


Material type


Material alloy


Material heat treat information

Stress on skin at upper/lower edged

Spectrum name


Flag for identifying steps


No. times to apply schedule


No. distinct blocks


Don't display spec blocks


Num steps/block


Schedule option


Load step number


Number of cycles


FMIN(1) t1 S0


FMAX(1) t2 S0


FMIN(2) t1 S1


FMAX(2) t2 S1


FMIN(3) t1 S3


FMAX(3) t2 S3


End manual input


Scaling Factor S0


Scaling Factor S1


Scaling Factor S3


Reference stress input


REFACT(1,1,1) S0


Ref Stress at t2


REFACT(2,1,1) S1


Ref Stress at t2


REFACT(4,1,1) S3


Ref Stress at t2


Do not enter schedule from file


Sblock case


Number of times to apply


End Spectrum input




Critical crack size/Residual Strength

The primary damage crack, crack A, is assumed to grow from a hole to the edge, during which the continuing damage crack, crack B, is growing from the opposite side of the hole towards an adjacent hole.  Once crack A reaches the edge it transitions into an edge crack with the crack tip at the tip of crack B.



Based on the calculations for growing the crack in NASGRO and the crack growth interactions, the life from initial crack size to failure is determined to be 41,412 flights.  The results of crack length versus life are shown in Figure SIE-3.7.  The life is given in numbers of flights.

Figure SIE-3.7.  Crack Growth Life for Problem SIE-3.


Inspection Intervals

The threshold and repeat intervals are calculated using the life reduction factors shown below. 

Life Reduction Factors:

K1 = 2.0

K2 = 3.0


Detectable crack length (HFEC around fastener head):

Number of flights @ detectable crack length, Ndet = 13,758 flights

Critical crack length (distance to next adjacent hole):

Number of flights @ critical crack length, Ncrit = 41,412 flights