Stress-Life (S-N) Approach
S-N Curve
The S-N curve, first developed by Wöhler, defines a relationship between stress and number of cycles to failure.
for segment 1
Where, is the nominal stress range, are the fatigue cycles to failure, is the first fatigue strength exponent, and is the fatigue strength coefficient.
The S-N approach is based on elastic cyclic loading, inferring that the S-N curve should be confined, on the life axis, to numbers greater than 1000 cycles. This ensures that no significant plasticity is occurring. This is commonly referred to as high-cycle fatigue.
S-N curve data is provided for a given material using the Materials module.
Multiple SN Curves
- Multi-mean S-N curves: group of S-N curves defined at different mean stress.
- Multi-ratio S-N curves: group of S-N curves defined at different stress ratio R.
- Multi-Haigh Diagram: group of Haigh curves defined at different Number of Cycles.
Rainflow Cycle Counting
Cycle counting is used to extract discrete simple "equivalent" constant amplitude cycles from a random loading sequence.
One way to understand "cycle counting" is as a changing stress-strain versus time signal. Cycle counting will count the number of stress-strain hysteresis loops and keep track of their range/mean or maximum/minimum values.
- Simple Load History:
Since this load history is continuous, it is converted into a load history consisting of peaks and valleys only.It is clear that point 4 is the peak stress in the load history, and it will be moved to the front during rearrangement (Figure 8). After rearrangement, the peaks and valleys are renumbered for convenience.
Next, pick the first three stress values (1, 2, and 3) and determine if a cycle is present.
If represents the stress value, point then:(2) (3) As you can see from Figure 8 , ; therefore, no cycle is extracted from point 1 to 2. Now consider the next three points (2, 3, and 4).(4) (5) , hence a cycle is extracted from point 2 to 3. Now that a cycle has been extracted, the two points are deleted from the graph.The same process is applied to the remaining points:(6) (7) In this case, , so another cycle is extracted from point 1 to 4. After these two points are also discarded, only point 5 remains; therefore, the rainflow counting process is completed.
Two cycles (2→3 and 1→4) have been extracted from this load history. One of the main reasons for choosing the highest peak/valley and rearranging the load history is to guarantee that the largest cycle is always extracted (in this case, it is 1→4). If you observe the load history prior to rearrangement, and conduct the same rainflow counting process on it, then clearly, the 1→4 cycle is not extracted.
- Complex Load HistoryThe rainflow counting process is the same regardless of the number of load history points. However, depending on the location of the highest peak/valley used for rearrangement, it may not be obvious how the rearrangement process is conducted.Figure 10 shows just the rearrangement process for a more complex load history. The subsequent rainflow counting is just an extrapolation of the process mentioned in the simple example above, and is not repeated here.Since this load history is continuous, it is converted into a load history consisting of peaks and valleys only:Clearly, load point 11 is the highest valued load and therefore, the load history is now rearranged and renumbered.
The load history is rearranged such that all points including and after the highest load are moved to the beginning of the load history and are removed from the end of the load history.
Equivalent Nominal Stress
Since S-N theory deals with uniaxial stress, the stress components need to be resolved into one combined value for each calculation point, at each time step, and then used as equivalent nominal stress applied on the S-N curve.
Various stress combination types are available with the default being "Absolute maximum principle stress". "Absolute maximum principle stress" is recommended for brittle materials, while "Signed von Mises stress" is recommended for ductile material. The sign on the signed parameters is taken from the sign of the Maximum Absolute Principal value.
"Critical plane stress" is also available as a stress combination for uniaxial calculations (stress life and strain life ).
For example, if number of planes requested is 20, then stress is calculated every 10 degrees.
By default, HyperLife also calculates at 𝜃 = 45 and 135-degree planes in addition to the requested number of planes. This is to include the worst possible damage if occurring on these planes.
Mean Stress Correction
Generally, S-N curves are obtained from standard experiments with fully reversed cyclic loading. However, the real fatigue loading could not be fully-reversed, and the normal mean stresses have significant effect on fatigue performance of components. Tensile normal mean stresses are detrimental and compressive normal mean stresses are beneficial, in terms of fatigue strength. Mean stress correction is used to take into account the effect of non-zero mean stresses.
The Gerber parabola and the Goodman line in Haigh's coordinates are widely used when considering mean stress influence, and can be expressed as:
Multi-mean SN curves
There are 2 available options for FKM correction in HyperLife. They are activated by setting FKM MSS to 1 slope/4 slopes in the Assign Material dialog.
- Regime 1 (R > 1.0)
- Regime 2 (-∞ ≤ R ≤ 0.0)
- Regime 3 (0.0 < R < 0.5)
- Regime 4 (R ≥ 0.5)
- Stress amplitude after mean stress correction (Endurance stress)
- Mean stress
- Stress amplitude
- Slope entered for region 2
If all four slopes are specified for mean stress correction, the corresponding Mean Stress Sensitivity values are slopes for controlling all four regimes. Based on FKM-Guidelines, the Haigh diagram is divided into four regimes based on the Stress ratio () values. The Corrected value is then used to choose the S-N curve for the damage and life calculation stage.
- Regime 1 (R > 1.0)
- Regime 2 (-∞ ≤ R ≤ 0.0)
- Regime 3 (0.0 < R < 0.5)
- Regime 4 (R ≥ 0.5)
- Fully reversed fatigue strength (Endurance stress)
- Mean stress
- Stress amplitude
- Slopes at each region
Gerber
Goodman
When SN curve is of the Stress Ratio R = -1
Gerber2
Improves the Gerber method by ignoring the effect of negative mean stress.
When SN curve is of the Stress Ratio R != -1
If , Same as Gerber
If ,
SODERBERG
Is slightly different from GOODMAN; the mean stress is normalized by yield stress instead of ultimate tensile stress.
- Equivalent stress amplitude
- Stress amplitude
- Mean stress
- Yield stress
FKM
Interpolate
Multi-Mean SN Curves
Case A
If a cycle has a mean stress of 150MPa at point A, HyperLife locates point 1 and point 2 in Figure 16. Then HyperLife linearly interpolates logN1 and logN2 with respect to mean stress in order to determine logN_A at mean stress 150MPa. Once logN_A is determined, life (N_A) and corresponding damage can be determined.
Case B
- Option 1 , Curve Extrapolation = False
- Use an SN curve of the maximum mean stress (the SN curve of mean stress 180 MPa in this case). In the example in HyperLife, N1 is the life HyperLife will report.
- Option 2 , Curve Extrapolation = True
- Extrapolate log(N) of the two SN curves with the highest mean stress values. In the example in Figure 16, log(N) will be extrapolated from log(N1) and log(N2) with respect to mean stress.
Case C
If the cycle has a mean stress less than the minimum mean stress of the curve set (90MPa in this case), HyperLife will use the SN curve of the minimum mean stress to determine life. In the example in Figure 16, life will be N2.
Multi-Stress Ratio SN Curve
Case A
If a cycle has R = -0.2 at point A, HyperLife locates point 1 and point 2 in Figure 17. Then HyperLife linearly interpolates logN1 and logN2 with respect to mean stress in order to determine logN_A at R = -0.2. Once R value and stress amplitude of the cycle are given, we can always calculate mean stress of the cycle. Once logN_A is determined, life (N_A) and corresponding damage can be determined. It is worthwhile to mention that HyperLife does not use stress ratio for interpolation because R can be an infinite value when maximum stress is zero.
Case B
- Option 1, Curve Extrapolation = False
- Use an SN curve of the maximum R (the SN curve of R= 0 in this case). In the example in Figure 17, N1 is the life HyperLife will report.
- Option 2, Curve Extrapolation = True
- Extrapolate log(N) of the two SN curves with the highest R values. In the example in Figure 17, log(N) will be extrapolated from log(N1) and log(N2) with respect to mean stress.
Case C
If the cycle has R less than the minimum R of the curve set (R= -1 in this case), HyperLife will use the SN curve of the minimum R to determine life. In the example in Figure 17, life will be N2.
Constant Life Haigh Diagram
Case A
If a cycle has a mean stress and stress amplitude at point A, HyperLife locates point 1 and point 2 in Figure 18. Life of point A should be between 1000 and 100000. HyperLife linearly interpolates log(1000) and log(100000) with respect to stress amplitude along Sm_A constant mean stress line in order to determine logN_A at point A. Once logN_A is determined, life (N_A) and corresponding damage can be determined.
Case B
If a point (mean stress, stress amplitude) is located above or below all the Haigh diagrams, life of the point is calculated by extrapolation of 2 highest or 2 lowest curves. In the example in Figure 18, log(1000) and log(100000) will be extrapolated with respect to stress amplitude along Sm_B constant mean stress line.
Case C
In this case, stress amplitude at point 5 and point 6 may be calculated from extrapolation. Once stress amplitudes become available at the 2 points, a procedure described in case A is applied.
Damage Accumulation Model
- Materials fatigue life (number of cycles to failure) from its S-N curve at a combination of stress amplitude and means stress level .
- Number of stress cycles at load level .
- Cumulative damage under load cycle.
The linear damage summation rule does not take into account the effect of the load sequence on the accumulation of damage, due to cyclic fatigue loading. However, it has been proved to work well for many applications.
Safety Factor
Safety factor is calculated based on the endurance limit or target stress (at target life) against the stress amplitude from the working stress history.
- Mean Stress = Constant
- Stress Ratio = Constant
The safety factor (SF) based on the mean stress correction applied is given by the following equations.
- Mean Stress = Constant
-
- Goodman or Soderberg
When SN curve is of the Stress Ratio R = -1
(20) = Target stress amplitude against the target life from the modified SN curve
= Stress amplitude after mean stress correction
When SN curve is of the Stress Ratio R != -1= Stress Amplitude
= Mean Stress
= Endurance limit obtained from SN curve with R ratio
= Mean Stress corresponding to
If ,
(21) If ,
If ,
If ,(22) - Gerber
(23) When SN curve is of the Stress Ratio R != -1(24) (25) (26) (27) - Gerber2
-
(28) -
(29)
When SN curve is of the Stress Ratio R != -1
If(30) (31) If ,
If ,
If ,(32) -
- FKM
(33) -
(34) -
(35) -
(36) -
(37)
-
- No Mean Stress Correction
(38) - Interpolate
- Goodman or Soderberg
Safety Factor with Multi-Mean
Safety Factor with Multi-Ratio
Safety Factor with Haigh
- Constant R : SF = OB/OA
- Constant mean : SF = OD/OC
If Haigh diagram for a target life is not defined by user, OptiStruct creates Haigh diagram for the target life. In Figure 26, if target life is 10000, and Haigh diagram for N=10000 is not defined, OptiStruct will created dashed curve to calculate Safety factor.
- Stress Ratio = Constant
-
- Goodman
When SN curve is of the Stress Ratio R = -1
(39) When SN curve is of the Stress Ratio R != -1
If ,(40) If ,
If ,
If ,
- Gerber
When SN curve is of the Stress Ratio R = -1
-
(41) -
(42)
When SN curve is of the Stress Ratio R != -1(43) (44) If ,
If ,
-
- Gerber2
-
(45) -
(46)
When SN curve is of the Stress Ratio R != -1
If(47) (48) If ,
If ,
If ,
-
- FKM
(49) = Corrected Stress Amplitude in Constant R mean stress correction
- No Mean Stress Correction
(50)
- Goodman
Safety Factor with Multi-Mean
Safety Factor with Multi-Ratio
Safety Factor with Haigh
- Constant R : SF = OB/OA
- Constant mean : SF = OD/OC
If Haigh diagram for a target life is not defined by user, OptiStruct creates Haigh diagram for the target life. In Figure 30, if target life is 10000, and Haigh diagram for N=10000 is not defined, OptiStruct will created dashed curve to calculate Safety factor.