Rainflow-counting algorithm

Last updated
Rainflow counting identifies the closed cycles in a stress-strain curve Rainflow counting vs stress-strain curve.svg
Rainflow counting identifies the closed cycles in a stress-strain curve

The rainflow-counting algorithm is used in calculating the fatigue life of a component in order to convert a loading sequence of varying stress into a set of constant amplitude stress reversals with equivalent fatigue damage. The method successively extracts the smaller interruption cycles from a sequence, which models the material memory effect seen with stress-strain hysteresis cycles. [1] This simplification allows the number of cycles until failure of a component to be determined for each rainflow cycle using either Miner's rule to calculate the fatigue damage, or in a crack growth equation to calculate the crack increments. [2] Both methods give an estimate of the fatigue life of a component. In cases of multiaxial loading, critical plane analysis can be used together with rainflow counting to identify the uniaxial history associated with the plane that maximizes damage. The algorithm was developed by Tatsuo Endo and M. Matsuishi in 1968. [3]

Contents

The rainflow method is compatible with the cycles obtained from examination of the stress-strain hysteresis cycles. When a material is cyclically strained, a plot of stress against strain shows loops forming from the smaller interruption cycles. At the end of the smaller cycle, the material resumes the stress-strain path of the original cycle, as if the interruption had not occurred. The closed loops represent the energy dissipated by the material. [1]

Figure 1: Uniform alternating loading Rainflow fig1.PNG
Figure 1: Uniform alternating loading
Figure 2: Spectrum loading Rainflow fig2.PNG
Figure 2: Spectrum loading

History

The rainflow algorithm was developed by T. Endo and M. Matsuishi (an M.S. student at the time) in 1968 and presented in a Japanese paper. The first english presentation by the authors was in 1974. They communicated the technique to N. E. Dowling and J. Morrow in the U.S. who verified the technique and further popularised its use. [1]

Downing and Socie created one of the more widely referenced and utilized rainflow cycle-counting algorithms in 1982, [4] which was included as one of many cycle-counting algorithms in ASTM E1049-85. [5]

Igor Rychlik gave a mathematical definition for the rainflow counting method, [6] thus enabling closed-form computations from the statistical properties of the load signal.

Algorithms

There are a number of different algorithms for identifying the rainflow cycles within a sequence. They all find the closed cycles and may be left with half closed residual cycles at the end. All methods start with the process of eliminating non turning points from the sequence. A completely closed set of rainflow cycles can be obtained for a repeated load sequence such as used in fatigue testing by starting at the largest peak and continue to the end and wrapping around to the beginning.

Four point method

Rainflow counting using the four point method. Any pair of turning points B,C that lie between adjacent points A and D is a rainflow cycle. Count and eliminate the pair B,C and continue processing the sequence until no more cycles can be extracted. Rainflow-4point-method.png
Rainflow counting using the four point method. Any pair of turning points B,C that lie between adjacent points A and D is a rainflow cycle. Count and eliminate the pair B,C and continue processing the sequence until no more cycles can be extracted.

This method evaluates each set of 4 adjacent turning points A-B-C-D in turn: [7]

  1. Any pair of points B-C that lies within or equal to A-D is a rainflow cycle.
  2. Remove the pair B-C and re-evaluate the sequence from the beginning.
  3. Continue until no further pairs can be identified.

Pagoda roof method

This method considers the flow of water down of a series of pagoda roofs. Regions where the water will not flow identify the rainflow cycles which are seen as an interruption to the main cycle.

  1. Reduce the time history to a sequence of (tensile) peaks and (compressive) valleys.
  2. Imagine that the time history is a template for a rigid sheet (pagoda roof).
  3. Turn the sheet clockwise 90° (earliest time to the top).
  4. Each "tensile peak" is imagined as a source of water that "drips" down the pagoda.
  5. Count the number of half-cycles by looking for terminations in the flow occurring when either:
    • case (a) It reaches the end of the time history;
    • case (b) It merges with a flow that started at an earlier tensile peak; or
    • case (c) An opposite tensile peak has greater or equal magnitude.
  6. Repeat step 5 for compressive valleys.
  7. Assign a magnitude to each half-cycle equal to the stress difference between its start and termination.
  8. Pair up half-cycles of identical magnitude (but opposite sense) to count the number of complete cycles. Typically, there are some residual half-cycles.

Example

Figure 3: Rainflow analysis for tensile peaks Rainflow analysis for tensile peaks.svg
Figure 3: Rainflow analysis for tensile peaks

The stress history in Figure 2 is reduced to tensile peaks in Figure 3 and compressive valleys in Figure 4. From the tensile peaks in Figure 3:

  • The first half-cycle starts at tensile peak 1 and terminates opposite a greater tensile stress, peak 3 (case c); its magnitude is 16 MPa (2 - (-14) = 16).
  • The half-cycle starting at peak 9 terminates where it is interrupted by a flow from earlier peak 8 (case b); its magnitude is 16 MPa (8 - (-8) = 16).
  • The half-cycle starting at peak 11 terminates at the end of the time history (case a); its magnitude is 19 MPa (15 - (-4) = 19).

Similar half-cycles are calculated for compressive stresses (Figure 4) and the half-cycles are then matched.

Figure 4: Rainflow analysis for compressive valleys Rainflow analysis for compressive valleys.svg
Figure 4: Rainflow analysis for compressive valleys
Stress (MPa)Whole cyclesHalf cycles
1020
1301
1611
1701
1901
2010
2210
2901

Related Research Articles

<span class="mw-page-title-main">Stress (mechanics)</span> Physical quantity that expresses internal forces in a continuous material

In continuum mechanics, stress is a physical quantity that describes forces present during deformation. For example, an object being pulled apart, such as a stretched elastic band, is subject to tensile stress and may undergo elongation. An object being pushed together, such as a crumpled sponge, is subject to compressive stress and may undergo shortening. The greater the force and the smaller the cross-sectional area of the body on which it acts, the greater the stress. Stress has dimension of force per area, with SI units of newtons per square meter (N/m2) or pascal (Pa).

The field of strength of materials typically refers to various methods of calculating the stresses and strains in structural members, such as beams, columns, and shafts. The methods employed to predict the response of a structure under loading and its susceptibility to various failure modes takes into account the properties of the materials such as its yield strength, ultimate strength, Young's modulus, and Poisson's ratio. In addition, the mechanical element's macroscopic properties such as its length, width, thickness, boundary constraints and abrupt changes in geometry such as holes are considered.

<span class="mw-page-title-main">Fatigue (material)</span> Initiation and propagation of cracks in a material due to cyclic loading

In materials science, fatigue is the initiation and propagation of cracks in a material due to cyclic loading. Once a fatigue crack has initiated, it grows a small amount with each loading cycle, typically producing striations on some parts of the fracture surface. The crack will continue to grow until it reaches a critical size, which occurs when the stress intensity factor of the crack exceeds the fracture toughness of the material, producing rapid propagation and typically complete fracture of the structure.

Stress–strain analysis is an engineering discipline that uses many methods to determine the stresses and strains in materials and structures subjected to forces. In continuum mechanics, stress is a physical quantity that expresses the internal forces that neighboring particles of a continuous material exert on each other, while strain is the measure of the deformation of the material.

<span class="mw-page-title-main">Hydrogen embrittlement</span> Reduction in ductility of a metal exposed to hydrogen

Hydrogen embrittlement (HE), also known as hydrogen-assisted cracking or hydrogen-induced cracking (HIC), is a reduction in the ductility of a metal due to absorbed hydrogen. Hydrogen atoms are small and can permeate solid metals. Once absorbed, hydrogen lowers the stress required for cracks in the metal to initiate and propagate, resulting in embrittlement. Hydrogen embrittlement occurs most notably in steels, as well as in iron, nickel, titanium, cobalt, and their alloys. Copper, aluminium, and stainless steels are less susceptible to hydrogen embrittlement.

<span class="mw-page-title-main">Residual stress</span> Stresses which remain in a solid material after the original cause is removed

In materials science and solid mechanics, residual stresses are stresses that remain in a solid material after the original cause of the stresses has been removed. Residual stress may be desirable or undesirable. For example, laser peening imparts deep beneficial compressive residual stresses into metal components such as turbine engine fan blades, and it is used in toughened glass to allow for large, thin, crack- and scratch-resistant glass displays on smartphones. However, unintended residual stress in a designed structure may cause it to fail prematurely.

<span class="mw-page-title-main">Fatigue limit</span> Maximum stress that wont cause fatigue failure

The fatigue limit or endurance limit is the stress level below which an infinite number of loading cycles can be applied to a material without causing fatigue failure. Some metals such as ferrous alloys and titanium alloys have a distinct limit, whereas others such as aluminium and copper do not and will eventually fail even from small stress amplitudes. Where materials do not have a distinct limit the term fatigue strength or endurance strength is used and is defined as the maximum value of completely reversed bending stress that a material can withstand for a specified number of cycles without a fatigue failure.

Mining accidents at the Markham Colliery at Staveley near Chesterfield, Derbyshire, England.

<span class="mw-page-title-main">Three-point flexural test</span> Standard procedure for measuring modulus of elasticity in bending

The three-point bending flexural test provides values for the modulus of elasticity in bending , flexural stress , flexural strain and the flexural stress–strain response of the material. This test is performed on a universal testing machine with a three-point or four-point bend fixture. The main advantage of a three-point flexural test is the ease of the specimen preparation and testing. However, this method has also some disadvantages: the results of the testing method are sensitive to specimen and loading geometry and strain rate.

<span class="mw-page-title-main">Paris' law</span> Formula in materials science

Paris' law is a crack growth equation that gives the rate of growth of a fatigue crack. The stress intensity factor characterises the load around a crack tip and the rate of crack growth is experimentally shown to be a function of the range of stress intensity seen in a loading cycle. The Paris equation is

<span class="mw-page-title-main">Universal testing machine</span> Type of equipment for determining tensile or compressive strength of a material

A universal testing machine (UTM), also known as a universal tester, universal tensile machine, materials testing machine, materials test frame, is used to test the tensile strength (pulling) and compressive strength (pushing), flexural strength, bending, shear, hardness, and torsion testing, providing valuable data for designing and ensuring the quality of materials. An earlier name for a tensile testing machine is a tensometer. The "universal" part of the name reflects that it can perform many standard tests application on materials, components, and structures.

6061 aluminium alloy is a precipitation-hardened aluminium alloy, containing magnesium and silicon as its major alloying elements. Originally called "Alloy 61S", it was developed in 1935. It has good mechanical properties, exhibits good weldability, and is very commonly extruded. It is one of the most common alloys of aluminium for general-purpose use.

<span class="mw-page-title-main">Flexural modulus</span> Intensive property in mechanics

In mechanics, the flexural modulus or bending modulus is an intensive property that is computed as the ratio of stress to strain in flexural deformation, or the tendency for a material to resist bending. It is determined from the slope of a stress-strain curve produced by a flexural test, and uses units of force per area. The flexural modulus defined using the 2-point (cantilever) and 3-point bend tests assumes a linear stress strain response.

<span class="mw-page-title-main">Vibration fatigue</span>

Vibration fatigue is a mechanical engineering term describing material fatigue, caused by forced vibration of random nature. An excited structure responds according to its natural-dynamics modes, which results in a dynamic stress load in the material points. The process of material fatigue is thus governed largely by the shape of the excitation profile and the response it produces. As the profiles of excitation and response are preferably analyzed in the frequency domain it is practical to use fatigue life evaluation methods, that can operate on the data in frequency-domain, s power spectral density (PSD).

Polymer fracture is the study of the fracture surface of an already failed material to determine the method of crack formation and extension in polymers both fiber reinforced and otherwise. Failure in polymer components can occur at relatively low stress levels, far below the tensile strength because of four major reasons: long term stress or creep rupture, cyclic stresses or fatigue, the presence of structural flaws and stress-cracking agents. Formations of submicroscopic cracks in polymers under load have been studied by x ray scattering techniques and the main regularities of crack formation under different loading conditions have been analyzed. The low strength of polymers compared to theoretically predicted values are mainly due to the many microscopic imperfections found in the material. These defects namely dislocations, crystalline boundaries, amorphous interlayers and block structure can all lead to the non-uniform distribution of mechanical stress.

Digital image correlation analyses have applications in material property characterization, displacement measurement, and strain mapping. As such, DIC is becoming an increasingly popular tool when evaluating the thermo-mechanical behavior of electronic components and systems.

<span class="mw-page-title-main">Striation (fatigue)</span>

Striations are marks produced on the fracture surface that show the incremental growth of a fatigue crack. A striation marks the position of the crack tip at the time it was made. The term striation generally refers to ductile striations which are rounded bands on the fracture surface separated by depressions or fissures and can have the same appearance on both sides of the mating surfaces of the fatigue crack. Although some research has suggested that many loading cycles are required to form a single striation, it is now generally thought that each striation is the result of a single loading cycle.

<span class="mw-page-title-main">Crack growth equation</span>

A crack growth equation is used for calculating the size of a fatigue crack growing from cyclic loads. The growth of a fatigue crack can result in catastrophic failure, particularly in the case of aircraft. When many growing fatigue cracks interact with one another it is known as widespread fatigue damage. A crack growth equation can be used to ensure safety, both in the design phase and during operation, by predicting the size of cracks. In critical structure, loads can be recorded and used to predict the size of cracks to ensure maintenance or retirement occurs prior to any of the cracks failing. Safety factors are used to reduce the predicted fatigue life to a service fatigue life because of the sensitivity of the fatigue life to the size and shape of crack initiating defects and the variability between assumed loading and actual loading experienced by a component.

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.

References

  1. 1 2 3 Endo, Tatsuo; Mitsunaga, Koichi; Takahashi, Kiyohum; Kobayashi, Kakuichi; Matsuishi, Masanori (1974). "Damage evaluation of metals for random or varying loading—three aspects of rain flow method". Mechanical Behavior of Materials. 1: 371–380.
  2. ↑ Sunder, R.; Seetharam, S. A.; Bhaskaran, T. A. (1984). "Cycle counting for fatigue crack growth analysis". International Journal of Fatigue. 6 (3): 147–156. doi:10.1016/0142-1123(84)90032-X.
  3. ↑ Matsuishi, M.; Endo, T. (1968). "Fatigue of metals subjected to varying stress". Japan Society of Mechanical Engineering.
  4. ↑ Downing, S.D.; Socie, D.F. (1982). "Simple rainflow counting algorithms". International Journal of Fatigue. 4 (1): 31–40. doi:10.1016/0142-1123(82)90018-4.
  5. ↑ Standard practices for cycle counting in fatigue analysis. ASTM E 1049-85. ASTM International. 2005.
  6. ↑ Rychlik, I. (1987). "A New Definition of the Rainflow Cycle Counting Method". International Journal of Fatigue. 9 (2): 119–121. doi:10.1016/0142-1123(87)90054-5.
  7. ↑ Lee, Yung-Li; Tjhung, Tana (2012). "Rainflow Cycle Counting Techniques". Metal Fatigue Analysis Handbook. pp. 89–114. doi:10.1016/B978-0-12-385204-5.00003-3. ISBN   9780123852045.
This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.