Rainflow method (cycle counting)
Rainflow method (cycle counting): Theoretical Foundations
What is the Rainflow Method?
Professor, what is the Rainflow Method?
Rainflow Counting is a technique to extract closed hysteresis loops (cycles) from a fluctuating load time history. It was named by Endo and Matsuishi (1968, Japan) by analogy to the flow of rainwater.
Why is Cycle Counting Necessary?
S-N curves and the Coffin-Manson equation give fatigue life for constant amplitude loading. Real loads are variable amplitude. The Rainflow Method decomposes variable loads into "a collection of many constant amplitude cycles," and cumulative damage is calculated using Miner's rule.
Algorithm
1. Extract peaks and valleys from the stress time history
2. Identify closed cycles as if "rainwater flows down a pagoda roof"
3. Record the stress range $\Delta\sigma$ and mean stress $\sigma_m$ for each cycle
4. Construct a histogram (Rainflow Matrix)
Summary
The "Rain Flowing Down a Roof" Algorithm
The name "Rainflow Method" comes from the fact that when the stress time history graph is rotated 90°, it resembles rainwater flowing down a roof. Proposed by Noboru Endo and Koji Matsuishi in 1968, it became a world standard as a method for accurately extracting fatigue cycles from irregular variable loads. It was standardized as ASTM E1049 in 1985.
Computational Methods for Rainflow method (cycle counting)
Rainflow Method Calculation
The Rainflow Method is built into all fatigue software. Input FEM stress time history → Rainflow → S-N/ε-N → Miner's rule automatic flow.
Also easily calculable in MATLAB/Python using the rainflow library.
Summary
How to Read a Rainflow Matrix
The result of the Rainflow Method is expressed as a Rainflow Matrix (RFC matrix), showing a frequency distribution with minimum stress on the vertical axis and maximum stress on the horizontal axis. Diagonal elements represent smaller amplitude, low cycle counts, while edge elements represent large amplitude, high damage cycles. Fatigue damage is dominated by elements far from the diagonal, and in actual evaluation, the accuracy of this edge data determines the accuracy of Miner's rule calculations.
Rainflow method (cycle counting) in Practice
Rainflow Method in Practice
Evaluating fatigue life from automotive road data (real driving stress time history). Aircraft flight spectra. Pressure cycles in piping.
Practical Checklist
Fatigue Evaluation from Real Road Data
In automotive durability testing, real road driving data is measured with strain gauges and analyzed using the Rainflow Method to evaluate fatigue damage. Ford standardized this approach with EDAP (Engineering Durability Analysis Package) in the 1980s, optimizing an algorithm that could process data from 10km of Belgian road driving in just seconds on 1990s PCs.
Rainflow method (cycle counting): Software & Solver Comparison
Tools
Built into all fatigue software. Also possible in MATLAB/Python.
Standard Differences and Implementation of Rainflow Counting
Cycle counting results for the same load history can differ by up to 15% between ASTM E1049-85 and DIN 45667. MSC Fatigue defaults to the 4-point method (ASTM compliant), nCode defaults to the HCM method (DIN compliant), and the automotive industry standard at ZF is the HCM method. Standard selection directly impacts fatigue life prediction accuracy.
Advanced Technology
Advanced Rainflow Method
Extension to Gaussian Processes and Non-stationary Loads
Traditional Rainflow Method assumes stationary random loads, but actual loads are often non-stationary processes where statistical properties change over time. The "windowed Rainflow Method" established in the 2010s performs cycle extraction within time windows, enabling handling of non-stationary loads. Applied to fatigue evaluation of wind turbines due to changing operating conditions.
Rainflow method (cycle counting): Common Issues & Debugging
Rainflow Method Troubles
Residual Cycle Problem with Truncated Data
When the beginning and end of measured data are incomplete, "residual" cycles occur in the Rainflow Method. Ignoring these leads to 10-20% underestimation of damage. ASTM E1049 recommends complementing residual cycles using the range-pair method, but in practice, designing the measurement protocol so the entire measurement interval is complete is the fundamental solution.