Fastran

From HandWiki
Short description: Crack growth calculation program

Fastran is a computer program for calculating the rate of fatigue crack growth by combining crack growth equations and a simulation of the plasticity at the crack tip.

Fastran models accelerations and retardation and other variable amplitude loading effects in crack growth using a crack closure model. The program uses a strip yield model of the crack tip that was first proposed by Dugdale to calculate the size of the plastic zone ahead of a crack tip. A series of elastic-perfectly plastic strips (originally 30 strips were used) that model the region both ahead and behind the crack tip is used to keep track of the plasticity produced at the crack tip.[1] As the crack grows, the strips are cut and leave a region of raised plastic material in the crack wake that prevents the complete closure of a crack. This profile of the crack is used to calculate the stress intensity factor level [math]\displaystyle{ K_\text{op} }[/math] at which the crack tip is fully open. The effective stress intensity factor range is then

[math]\displaystyle{ \Delta K_\text{eff} = K_\text{max} - \text{max}(K_\text{op}, K_\text{min}) }[/math]

which allows the rate of growth for the loading cycle to be obtained from the crack growth equation. The rate of crack growth is then calculated from

[math]\displaystyle{ {da \over dN} = f(\Delta K_\text{eff}) }[/math]

History

Fastran was written in the 1980's by James C. Newman while at NASA and is an acronym derived from NASA FATIGUE CRACK GROWTH STRUCTURAL ANALYSIS.[2] Crack closure was first observed by Wolf Elber as propping open a crack tip resulting in a reduction of the full stress intensity range or crack tip driving force.[3] It was assumed this was due to plasticity at the crack tip preventing the fracture surfaces from fully closing.

A similar program CORPUS was also developed around the same time by de Koning.[4]

FASTRAN is written in the Fortran programming language.

Features

Geometry Factors

The geometry factor [math]\displaystyle{ \beta }[/math] relates the far-field stresses to the region near the crack tip. Many standard geometry factors are supplied in the program. These scaling factors allow the calculation of the stress intensity factor from the applied loading sequence using

[math]\displaystyle{ K = \beta \sigma \sqrt{\pi a} }[/math]

where [math]\displaystyle{ \sigma }[/math] is the applied far field stress and [math]\displaystyle{ a }[/math] is the crack length.

The loading sequence is given as a file of sequential turning points that represent the loading sequence. This in combination with a load factor is used to supply the far-field stress of the given geometry. The load sequence is converted into a series of individual load cycles by a method known as rainflow on the fly which is a modified form of the standard rainflow-counting algorithm.

The closure model has also been used to explain the increase rate of growth seen with small cracks known as the small crack effect.

Crack growth equations

Fastran has a variety of crack growth equations built in along with piece wise linear equations that can be read from file.

Theory

This model allows the calculation of the stress ratio [math]\displaystyle{ R }[/math] or mean stress effect that gives rise to the increased rate of crack growth at higher stress ratios.[5] Experiments have shown the crack is typically open at [math]\displaystyle{ R\gt 0.7 }[/math]. In addition the model is able to predict retardation due to overloads which increase the plastic material in the wake of the crack. It also explains the acceleration due to underloads where the crack growth rate increases following an underload which compresses the crack faces together and reduced the degree of interference lowering [math]\displaystyle{ K_\text{op} }[/math].

The onset of plasticity is given by the flow stress whose value typically lies mid-way between the yield and ultimate stresses. The flow stress scaling parameter [math]\displaystyle{ \alpha }[/math] is used to adjust the flow stress to the degree of restraint experienced at the crack tip. This value reflects the stress state at the crack tip and typically lies between a value of [math]\displaystyle{ \pi }[/math] for plane stress and [math]\displaystyle{ 2\pi }[/math] for plane strain. The parameter is also used as an adjustment variable to correct the rate of crack to match test data.

Limitations

Plasticity will be greater in regions of plane stress but Fastran only models the crack as a 2d cross section.

Usage

Fastran has been used in the research community and for maintaining the safe life of aircraft as the C-130 used by the USAF, RAF and RAAF. If forms a component of the crack growth program Nasgro.[6]

References

  1. Dugdale, D. S. (1960). "Yielding of steel sheets containing slits". Journal of the Mechanics and Physics of Solids 8 (2): 100-104. 
  2. Newman, Jr., J. C. (1992). FASTRAN II - A fatigue crack growth structural analysis program. NASA. https://docs.google.com/viewer?a=v&pid=sites&srcid=ZnJhY3R1cmVsYWIuY29tfGhvbWV8Z3g6NjBiZTcwOGNkNDc1MmM2NA. Retrieved 6 Jan 2020. 
  3. Elber, Wolf (1971). "The significance of fatigue crack closure". Damage tolerance in aircraft structures, ASTM International: 230–242. 
  4. de Koning, A. U. (1981). "A simple crack closure model for prediction of fatigue crack growth rates under variable-amplitude loading". Fracture Mechanics (ASTM) (STP 743): 63-85. 
  5. Maddox, S. J. (1975). "The effect of mean stress on fatigue crack propagation—A literature review". International Journal of Fracture 1 (3). 
  6. NASGRO Fracture Mechanics and Fatigue Crack Growth Analysis Software, Version 4.02. SwRI. 2002.