Structural Fatigue Engineering, ASTM E1049 Spectrum Analysis & Miner Damage
ASTM E1049 Four-Point Rainflow & Equivalent Stress Formula
Evaluate ASTM E1049 Four-Point Rainflow Closed-Hysteresis Bounding Condition on Four Consecutive Stress Reversals (S1..S4), Extracted Inner Cycle Range ΔS_23 vs Outer Envelope ΔS_14 (MPa), and Basquin m-Power Equivalent Constant-Amplitude Stress Range ΔS_eq.
Extract closed stress-strain hysteresis loops with ASTM E1049’s 4-point bounding rule!
Real wind-turbine blades, bridge girders, and vehicle suspension arms experience irregular variable-amplitude stress reversals rather than clean sine waves! Under ASTM E1049’s Four-Point Rainflow Counting algorithm, four consecutive alternating peak/valley reversals (S1, S2, S3, S4) enclose a CLOSED inner stress-strain hysteresis loop (S2, S3) if and only if the inner pair lies completely inside the outer pair’s span: MIN(S2, S3) >= MIN(S1, S4) AND MAX(S2, S3) <= MAX(S1, S4)! Extracting that full inner cycle ΔS_inner = ABS(S3 - S2) leaves the outer path S1 -> S4 (ΔS_outer = ABS(S4 - S1)), and under a repeating periodic block with S-N fatigue slope m (e.g., m = 3 for welded steel), Miner’s rule combines them via the m-th power mean ΔS_eq = ((ΔS_inner^m + ΔS_outer^m) / 2)^(1/m)—NEVER a simple arithmetic average!
Verify on a Welded Bridge Girder Strain History—S1_Valley = 40.0 MPa (B2), S2_Peak = 260.0 MPa (B3), S3_Intermediate_Valley = 100.0 MPa (B4), S4_Major_Peak = 340.0 MPa (B5), and Eurocode_Weld_Slope_m = 3.0 (B6): because MIN(260,100) = 100 >= 40 and MAX(260,100) = 260 <= 340, E2 confirms "Closed Full Cycle (S2-S3 Extracted)", extracting Inner Range ΔS_23 = 160.00 MPa (E3) inside Outer Envelope ΔS_14 = 300.00 MPa (E4), with Periodic 2-Cycle Equivalent Stress Range ΔS_eq = ((160^3 + 300^3) / 2)^(1/3) = 249.59 MPa (E5—vs 230 MPa naive average that understates m=3 fatigue damage by 27.8%)!
| ASTM E1049 4-Point Metric (S = 40 -> 260 -> 100 -> 340 MPa, Slope m = 3) | Formula | Verified Output |
|---|---|---|
| 4-Point Bounding Gate | MIN(S2,S3)>=MIN(S1,S4) & MAX(S2,S3)<=MAX(S1,S4) | Closed Full Cycle (S2-S3 Extracted) |
| Inner Cycle Range ΔS_23 (MPa) | ABS(S3 - S2) | 160.00 MPa |
| Outer Span Range ΔS_14 (MPa) | ABS(S4 - S1) | 300.00 MPa |
| Equivalent Stress ΔS_eq (MPa) | (((ΔS_23^m) + (ΔS_14^m)) / 2)^(1/m) | 249.59 MPa (vs 230 MPa Naive Avg) |
ASTM E1049 Four-Point Rainflow Closure & Equivalent Stress Formula
Validate alternating peak/valley reversals ((S2 - S1)*(S3 - S2) < 0 and (S3 - S2)*(S4 - S3) < 0) and fatigue slope m > 0!
=LET(s_1, B2, s_2, B3, s_3, B4, s_4, B5, sl_m, B6, IF(OR((s_2 - s_1) * (s_3 - s_2) >= 0, (s_3 - s_2) * (s_4 - s_3) >= 0, sl_m <= 0), "Invalid Reversal Sequence", IF(AND(MIN(s_2, s_3) >= MIN(s_1, s_4), MAX(s_2, s_3) <= MAX(s_1, s_4)), "Closed Full Cycle (S2-S3 Extracted)", "Open: Wait for Next Reversal")))Unbounded Valley Excursion Check (Set S3 = 10.0 MPa Below S1 = 40.0 MPa)
Drop S3 B4 from 100.0 MPa down to 10.0 MPa (dipping below starting valley S1 = 40.0 MPa): MIN(260, 10) = 10 < MIN(40, 340) = 40 violates the lower rainflow bound, switching E2 cleanly to "Open: Wait for Next Reversal"!
=IF(E2 = "Invalid Reversal Sequence", E2, ROUND((((ABS(B3 - B4) ^ B6) + (ABS(B5 - B2) ^ B6)) / 2) ^ (1 / B6), 2))Try the example
Use the original values or change them to test your understanding.