**Materials Sciences and Applications**

Vol.09 No.04(2018), Article ID:83797,13 pages

10.4236/msa.2018.94025

A Method of Determining Realistic Stress S/N Curves by Interpolations and Extrapolations of Two Known Best-Fit S/N Curves for Fatigue Life Predictions

Claudiu Gudas^{ }

C & D Gudas, Sydney, Australia

Copyright © 2018 by author and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: February 10, 2018; Accepted: April 15, 2018; Published: April 18, 2018

ABSTRACT

The design and sizing of new mechanical components which are intended to operate under cyclic loads often require an acceptable level of confidence that the components will meet pre-defined fatigue strength objectives for crack initiation. For loads with multiple amplitudes and mean values, models based on Palmgren-Miner’s linear cumulative damage hypothesis and on multiple S/N curves (stress-no. cycles) are widely used in estimations for the duration of crack initiation. In this paper a procedure for generating S/N curves for multiple stress ratios by interpolation or extrapolation from the data available for two such curves is proposed. At any number of cycles, the stress for crack initiation is calculated from the far field macroscopic stresses of the known curves using Dang Van fatigue criterion and microscopic stresses evaluated at the grain level. An algorithm is presented for uniaxial loading and results verifications against curves established from laboratory tests with elastic and plastic stresses are shown and discussed for notched and un-notched specimens.

**Keywords:**

Crack Initiation, Extrapolation, Fatigue Life, Interpolation, Stress Ratio, The Dang Van Criterion

1. Introduction

With the rise of fracture mechanics, the certification against fatigue has largely switched towards crack-propagation models which correlate parameters like the crack tip stress intensity factor to its growth rate. At the same time new stronger materials have led to higher working stress and, in many situations, having higher crack growth rates, thus once initiated the duration of crack growth to a critical size is reduced [1] . Many commercial software packages dedicated to crack growth analysis incorporate modules for estimation of crack initiation. These modules are using strain-life initiation methodologies like Smith-Watson-Topper and require stabilized cyclic stress-strain data [2] . Crack growth analyses are significantly more laborious than traditional fatigue analyses. Often they can be done efficiently only after the design is finalized. In these situations also, the ability to estimate early from a traditional fatigue life approach the duration for crack initiation and growth to a predefined size is very beneficial during early design stages.

In the traditional approaches to fatigue strength evaluation for loading spectra containing cycles with multiple maximum and mean stress values, the degree of cumulative damage at each stress level is calculated from S/N curves. There are many forms in which these curves could be presented and most of them are obtained from uni-axial cupon tests after statistical post-processing of the results to account for scatter. For uni-axial tests, all these forms can be converted to express the cycles in terms of maximum stress and stress ratio (R). In case of multiaxial loading, more parameters are needed to specify the phase characteristics or the principal stresses and angles. The failure criteria used during testing can represent either the crack initiation, when the length of the crack is about 1.0 mm or 0.040 inch or complete fracture. In this approach only curves of first criteria are considered.

Due to its simplicity and larger availability of uni-axial S/N data, the evaluation of crack initiation using Miner’s rule is a pragmatic and quick alternative to other stress based methods like the critical plane approach, Crosslands, Sines.

In practice however, during and sizing, it is frequent that the data available is for a limited number of stress ratios which often do not cover the range of stresses required. The method presented below enables generation of S/N curves for new stress ratios from two known curves of stress ratios R_{1} and R_{2}.

2. Method Outline and Verification Details

2.1. Brief Overview of Dang Van Fatigue Criterion

During fatigue testing of materials, when the peak stresses exceed the yield limit they generate plastic flow and residual stresses [3] . At the same time, due to Bauschinger effects, after a certain amount of plastic deformation in one direction, the material yields at a lower level in the opposite direction. Within several cycles the stress-strain variations stabilize to a shake down state where the residual stresses alleviate the applied loads and induce a state of elastic response. In an elastic or shacked-down state no fatigue damage is supposed to occur because of no further strain accumulations. Tests however clearly demonstrate that fatigue occurs during cycles with peak stresses less than half of the material’s yield values.

Dang Van took the above concept to the grain level and noted that plasticity at this level can exist even when the macroscopic stresses are in the elastic range. He noted that a crack will initiate when the microscopic shake-down limits are just exceeded (called it a pseudo shake-down state) [4] [5] [6] .

A brief outline of the mathematical formulation for estimating the microscopic stresses:

・ The macroscopic stress tensor Σ is expressed in terms of hydrostatic pressure PH and the deviatoric S tensor;

・ Similarly, the microscopic stress tensor σ, is formulated in terms of microscopic hydrostatic pressure and deviatoric tensor, ph and s, respectively;

・ The relations between the macro and microscopic stresses:

o the hydrostatic pressure is the same at macro and microscopic levels (PH = ph);

o the macroscopic deviatoric stresses are given by the sum of the microscopic deviatoric and stabilized residual stresses q;

It was postulated in [4] that the stabilized residual stresses are at the center of the smallest hypersphere in the six-dimensional deviatoric stress space that completely encloses the load path. Tensor q needs to include only the deviatoric stresses q*, as only these stresses affect sliding along intergranular preferential slip bands (PSB) and the initiation of shear induced cracks.

$\underset{\_}{S}=\underset{\_}{s}+{\underset{\_}{q}}^{\ast}$ (1)

Noting that both, the hydrostatic pressure ph and maximum shear σ_{max} of the microscopic stresses affect crack initiation, Dang Van et al. [2] formulated a fatigue crack criterion in the ph-τ space as being defined by the two lines D and D' shown in Figure 1. The lines are given by Equation (2), where the coefficients a and b are material dependent real constants.

${\tau}_{\mathrm{min}/\mathrm{max}}=\pm a\ast ph+b$ , (2)

Figure 1. Fatigue initiation lines in the microscopic ph - τ space.

and:

$\left|{\tau}_{\mathrm{min}}\right|=\left|{\tau}_{\mathrm{max}}\right|=Tresca\left(s\right)$ . (3)

Test verifications of this criterion show that coefficients a and b in Equation (2) remain constant with the change of stress ratios from −1.75 to +0.2 and with type of loading (shear, tension and biaxial). It is also shown that with fatigue life decrease there is a substantial decrease in |a| at the same time as an increase in b [5] .

2.2. Method Outline

For any of the known uni-axial stress S/N curves of stress ratio R_{k}, k = {1, 2}, the temporal variation of the stress is assumed to be of the form:

${\Sigma}_{k}=\left(\begin{array}{c}{\sigma}_{k\_x}\\ {\sigma}_{k\_y}\\ {\tau}_{k\_xy}\end{array}\right)=\left(\begin{array}{c}\frac{{A}_{k}}{2}\ast \left[\frac{1+{R}_{k}}{1-{R}_{k}}+\mathrm{sin}\left({\omega}_{k}t\right)\right]\\ 0\\ 0\end{array}\right)$ , (4)

where the amplitude of cyclic stresses A_{k} is defined by:

${A}_{k}=\frac{2\ast {\sigma}_{k\_\mathrm{max}}}{\frac{1+{R}_{k}}{1-{R}_{k}}+1}$ , if ${R}_{k}<1$ and ${A}_{k}=\frac{2\ast {\sigma}_{k\_\mathrm{max}}}{\frac{1+{R}_{k}}{1-{R}_{k}}-1}$ , if ${R}_{k}>1$ (5)

The macro and microscopic hydrostatic pressures at the measurement locations are:

$\begin{array}{c}P{H}_{k}\left(t\right)=p{h}_{k}\left(t\right)\\ =\frac{{A}_{k}}{6}\ast \left[\frac{1+{R}_{k}}{1-{R}_{k}}+\mathrm{sin}\left({\omega}_{k}t\right)\right]\end{array}$ , (6)

and the macroscopic deviatoric stresses,

$\left(\begin{array}{c}{S}_{k\_1}\\ {S}_{k\_2}\\ {S}_{k\_3}\end{array}\right)=\left[\frac{1+{R}_{k}}{1-{R}_{k}}+\mathrm{sin}\left({\omega}_{k}t\right)\right]\ast \left(\begin{array}{c}\frac{{A}_{k}}{3}\\ \frac{-{A}_{k}}{6}\\ \frac{-{A}_{k}}{6}\end{array}\right)$ (7)

For uniaxial loading the load path is a straight line and the deviatoric of the microscopic residual stress tensor is:

$\left(\begin{array}{c}{q}_{k\_1}^{*}\\ {q}_{k\_2}^{*}\\ {q}_{k\_3}^{*}\end{array}\right)=\frac{1+{R}_{k}}{1-{R}_{k}}\ast \left(\begin{array}{c}\frac{-{A}_{k}}{3}\\ \frac{{A}_{k}}{6}\\ \frac{{A}_{k}}{6}\end{array}\right)$ (8)

From Equation (1) the microscopic deviatoric stresses are:

$\left(\begin{array}{c}{s}_{k\_1}\\ {s}_{k\_2}\\ {s}_{k\_3}\end{array}\right)=\mathrm{sin}\left({\omega}_{k}t\right)\ast \left(\begin{array}{c}\frac{{A}_{k}}{3}\\ \frac{-{A}_{k}}{6}\\ \frac{-{A}_{k}}{6}\end{array}\right)$ (9)

The resulting maximum microscopic shear is

${\tau}_{k\_\mathrm{max}}=Tresca\left({s}_{k}\right)=\left|\frac{{A}_{k}}{4}\ast \mathrm{sin}\left({\omega}_{k}t\right)\right|,\text{\hspace{0.17em}}\text{\hspace{0.17em}}k=\left\{1,2\right\}$ (10)

With
$\left(p{h}_{1},{\tau}_{1\_\mathrm{max}}\right)$
and
$\left(p{h}_{2},{\tau}_{2\_\mathrm{max}}\right)$
calculated for the same life N_{f}, the values of a and b defining the lines for the Dang Van fatigue criterion at N_{f} are given by Equation (2).

From Equation (5), for a life of N_{f} cycles, the amplitude of cycles of other stress ratio R is given by:

$A=\frac{b}{\frac{1}{4}-\frac{a}{6}\ast \left(\frac{1+R}{1-R}+1\right)}$ , if ${R}_{k}<1$ and $A=\frac{b}{\frac{1}{4}-\frac{a}{6}\ast \left(\frac{1+R}{1-R}-1\right)}$ , if ${R}_{k}>1$ (11)

Based on the study in reference [7] the equations above were assessed as suitable for finite lifetimes and for plastic macroscopic stresses.

Referring to Figure 1, for any cycle of stress ratios R, the points 1, 2 and 3 defined by Equation 6 and 10 when sin(ωt) takes the values of 1, 0 and −1, respectively are:

・ Point 1 of coordinates $p{h}_{1}=\frac{A}{6}\ast \left[\frac{1+R}{1-R}+1\right]$ and ${\tau}_{\mathrm{max}1}=\frac{A}{4}$ ;

・ Point 2 of coordinates $p{h}_{2}=\frac{A}{6}\ast \left[\frac{1+R}{1-R}\right]$ and ${\tau}_{\mathrm{max}2}=0$ ;

・ Point 3 of coordinates $p{h}_{3}=\frac{A}{6}\ast \left[\frac{1+R}{1-R}-1\right]$ and ${\tau}_{\mathrm{max}3}=-\frac{A}{4}$ .

The slopes m_{12} and m_{23} for segments 12 and 23 show that the three points are on the same line independent of R:

${m}_{12}=\frac{{\tau}_{\mathrm{max}1}}{p{h}_{1}-p{h}_{2}}=\frac{\frac{A}{4}}{\frac{A}{6}\ast \left[\frac{1+R}{1-R}+1\right]-\frac{A}{6}\ast \left[\frac{1+R}{1-R}\right]}=\frac{6}{4}=1.5$

${m}_{23}=\frac{-{\tau}_{\mathrm{max}3}}{p{h}_{2}-p{h}_{3}}=\frac{\frac{A}{4}}{\frac{A}{6}\ast \left[\frac{1+R}{1-R}\right]-\frac{A}{6}\ast \left[\frac{1+R}{1-R}-1\right]}=\frac{6}{4}=1.5$

For the particular case when σ_{max} = 0 (R = ∞), the microscopic hydrostatic pressure and shear stress for σ_{min} are given by the intersection of line D' with the line 1’2’3’ passing though (0, b):

$\tau =1.5\ast p{h}_{\mathrm{min}}+b,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\tau =a\ast p{h}_{\mathrm{min}}+b$ (12)

2.3. Verification Details

The algorithm outlined above was verified using Mathcad v.14, for the materials and specimen forms indicated in Table 1. The data used for these verifications contains S/N curves for multiple stress ratios, explicit equations for the fatigue models used and ranges of applicability. Many of the curves have stresses in plastic range enabling verifications of the method at stress levels typical for fail-safe designs.

For each test, two curves R_{1} < R_{2}, were selected as “known data” and discretized in three columns tables containing stress, life and log_{10}(life). For each of the “new” stress ratio R, relative differences were calculated in terms of log_{10}(life) and in terms of maximum stress, using:

${\delta}_{\text{LogLife}}\left(\sigma \right)=\frac{{\mathrm{log}}_{10}\left({n}_{0}\right)-{\mathrm{log}}_{10}\left(n\right)}{{\mathrm{log}}_{10}\left({n}_{0}\right)}$ and ${\delta}_{\sigma}\left(\text{life}\right)=\frac{{\sigma}_{0}-\sigma}{{\sigma}_{0}}$ (13)

n_{0}―life at σ, for the curve generated from published equation

n―life at σ, for the curve generated from two known S/N curves

σ_{0}―stress at given life, for the curve generated from published equation

σ―stress at given life, for the curve generated from two known S/N curves

During the verifications, all tables were limited to 15 lines with constant increments in stress or in log_{10}(life). The relative differences were calculated at each entry in the output tables and checked for direction and trends. Cubic splines in stress and logarithm of life were used when additional points were required to extend the curves.

The testing process included extrapolations and interpolations. For extrapolations R_{1} and R_{2} were selected for two enveloping cases and for several intermediate ones.

・ In the first enveloping case R_{1} was selected at the lowest available stress ratio and R_{2} was 0.1 higher than R_{1}. New curves, in increments of 0.1, were extrapolated to the highest stress ratio in the range of applicability and compared

Table 1. Materials and configurations used for verifications.

^{[a]}kilo pound /in^{2}, (1 ksi = 6.895 MPa); ^{[b]}Range of stress ratios; ^{[c]}From reference [8] .

with the curves generated from equations;

・ In the reverse case, the two R_{1} and R_{2} curves were selected at the highest end and the new curves were extrapolated in a similar manner to the lowest available stress ratio;

・ In the intermediate cases, like that exemplified in Figure 2, the two known curves were aleatory placed and extrapolations were conducted upwards and/or downwards.

For interpolations, the extreme envelope considered was when the two sets of known data were at the lowest and respectively the highest values of the published range of applicability. Additional tests were conducted between stress ratio of zero and the minimum or maximum available.

3. Results

The results for all S/N curves generated by extrapolation were conservative (safer than the results of the curves constructed from the published equations). It was observed that as the distance of extrapolation increases so does the level of conservatism. For interpolations the resulting curves are optimistic with the relative differences increasing towards the middle of interval between R_{1} and R_{2}.

Focusing on the enveloping tests for which the relative differences have the

Figure 2. Graphical sample of S/N curves in verifications. (a) S/N curves generated by interpolation and extrapolation from R_{1} and R_{2}; (b) Reference data for R_{1}, R_{2} and for verifications of relative differences (extract from [8] )

largest values and using as parameter δR, defined by the distance between the new stress ratio and the closest of R_{1} or R_{2}:

・ For extrapolations: At the left end of the curves representing the high stresses and shorter lives, both δ_{σ} (life) and δ_{LogLife}(σ) are less than 10% up to a δR of around 0.4 - 0.6 and they increase afterwards to 25% as δR approaches 1.5. At the other end (of low stresses and long lives) the estimates are smaller than 15% for δR less or equal to 0.5 and increase to 30% afterwards. At this end however, the runout or the endurance limits will frequently reduce these differences.

・ For interpolations: When the distance between R_{1} and R_{2} is less than 1.2 - 1.3 the relative differences are below 8% in both measures. They increase to 15% for larger distances between R_{1} and R_{2}.

Figure 3 shows the maximum values of δ_{LogLife} (σ) with the increase in δR for the extrapolation envelopes while Figure 4 shows maximum values of δ_{σ} (life). When the distance between R_{1} and R_{2} is larger than 0.1 the maximum relative differences vs. δR are smaller than those shown in two figures referenced above.

For interpolations the results of the enveloping cases are shown in Figure 5. Table 2 and Table 3 summarize the differences when one of the known stress ratios is zero.

Figure 3. Extrapolations―Maximum relative differences in δ_{LogLife}(σ). (left) for upward extrapolations; (right) for downward extrapolations (top) for high stresses; (bottom) for low stresses.

Figure 4. Extrapolations―Maximum relative differences in δ_{σ} (life). (left) for upward extrapolations; (right) for downward extrapolations (top) for shot lives; (bottom) for long lives.

4. Discussion

The tests show that for interpolations and extrapolations with moderate δR (≤0.5) the results are comparable with test data. For very large extrapolations, the curves might be too conservative from a weight or economic perspective and another curve reducing the span between the known stress ratios is required. When the known data was for high tension stress ratios, cubic splines proved to be effective and accurate in extending the new curves of lower stress ratios.

These tests also showed realistic results for notched specimens with moderate stress concentration factors (orange and black entries in Figure 3 to Figure 5). Note that for these situations the cyclic stress are specified in terms of net or gross area averages and the lines D and D’ in Figure 1, are only far-away indicators for the crack initiation at the notch.

It is therefore noted that in current form the algorithm presented in Equations (1) to (12) be limited to notch specimens with a stress concentration factor not larger than 2.4 when the “known” data represents the average area stresses.

For larger concentration factors, the macroscopic stresses at the notch enter quickly in plastic range and the changes of the local values from the net or gross area far field average become significant. In these situations, for part sizing

Figure 5. Interpolations―Minimum relative differences. (left) measured as δ_{LogLife} (σ); (right) measured as δ_{σ} (life).

Table 2. Relative differences for interpolations between R_{1} & R_{2} = 0.00.

Table 3. Relative differences for interpolations between R_{1} = 0.00 & R_{2}.

estimations, it is necessary to work with the stresses at the notch established by using non-linear finite element results calculated with material’s stress-strain curve or a Ramberg-Osgood type idealization and generate new input S/N curves from the known data. After computation, the curves for the new stress ratios may converted back to net or gross area average stress.

Often for high concentration factors in sheet/plate specimens, the S/N data available is for curves of constant mean stress not for constant R. In these situations the term R for the stress ratio in Equation (4) to (11) has to be replaced with:

${\sigma}_{\text{mean}}=\frac{{\sigma}_{\mathrm{max}}}{R-{\sigma}_{\mathrm{max}}}$ , (14)

and the curves for the notch area derived as described above. In these cases, for uni-axial estimates, principal or von Mises stresses might be more representative.

While the estimates by interpolation of S/N curves can be done quite easy graphically or algebraically or, they can be simply avoided at the design stage by using the conservative curve, extrapolations of the same nature are rarely advisable. The algorithm proposed provides a mean for generating conservative S/N curves by extrapolation and the verifications conducted enable a quantifiable estimate of maximum errors. For predictions of fatigue performance the normal safety factors are required to address the scatter in fatigue test results.

5. Summary

The aim of this study was to develop a method that enables quick evaluation of the duration for crack initiation during the design stage of parts expected to see loading spectra with multiple amplitude and mean values. For this purpose a method for generating S/N curves for multiple stress ratios by interpolation or extrapolation from the data available in two such curves was developed and tested. Generation of similar curves for multiple mean stresses was discussed in the context of notches with high concentration factors. Combining this procedure with Miner’s rule and the knock down factors specific to particular applications, durations for crack initiation can be estimated for spectra with many loads. The procedure proposed has the following characteristics:

・ Is based entirely on net or gross area stresses which can be measured or derived conveniently;

・ Is applicable to cycles with bulk stresses in elastic and/or plastic regime;

・ Is working with digitized curves;

・ Provides results of a quantifiable precision, conservative in all cases of extrapolations and slightly optimistic for interpolations;

・ Together with statistically based methods (like the staircase method), this procedure is able to provide testing laboratories with alternatives to physical testing when developing stress based life curves;

Acknowledgements

Many thanks to my son Michael for his assistance in with the research work and consultations and to Steve Swift for advice and telephone consultations during this project.

Cite this paper

Gudas, C. (2018) A Method of Determining Realistic Stress S/N Curves by Interpolations and Extrapolations of Two Known Best-Fit S/N Curves for Fatigue Life Predictions. Materials Sciences and Applications, 9, 368-380. https://doi.org/10.4236/msa.2018.94025

References

- 1. Peridas, G., Korsunsky, A.M. and Hills, D.A. (2003) The Relationship between the Dang Van Criterion and the Traditional Bulk Fatigue. http://journals.sagepub.com/doi/abs/10.1243/030932403765310536?journalCode=sdjb
- 2. LexTech Inc. (2018) AFGROW Verion 5.03.01.21—Users Guide and Technical Manual.
- 3. Suresh, S. (1998) Fatigue of Materials. 2nd Edition, Cambridge University Press.
- 4. Dang Van, K., Griveau, B. and Message, O. On a New Multiaxial Fatigue Limit Criterion: Theory and Application. Biaxial and Multiaxial Fatigue EGF3. In: Brown, M.W. and Miller, K.J., Eds., International Conference on Biaxial/Multiaxial Fatigue (2nd: 1985: University of Sheffield) Mechanical Engineering Publications, London, 459-478.
- 5. Dang Van, K., Cailletaud, G., Flavenot, J.F., Le Douaron, A. and Lieurade, H.P. Criterion for High Cycle Fatigue Failure under Multiaxial Loading. Biaxial and Multiaxial Fatigue EGF3. In: Brown, M.W. and Miller, K.J., Eds., International Conference on Biaxial/Multiaxial Fatigue (2nd: 1985: University of Sheffield) Mechanical Engineering Publications, London, 479-496.
- 6. Dang Van, K., Le Douaron, A. and Lieurade H.P. Criterion for High Multiaxial Fatigue Limit: A New Approach. https://www.sciencedirect.com/science/article/pii/B978148328440850185Xhttp://www.gruppofrattura.it/ocs/index.php/ICF/ICF6/paper/viewFile/3155/8217
- 7. Charkaluk, E., Constantinescu, A., Maitournam, H. and Dang Van, K. Revisiting the Dang Van Criterion. https://www.sciencedirect.com/science/article/pii/S1877705809000344
- 8. Federal Aviation Administration—Metallic Materials Properties Development and Standardization, Version 7 (MMPDS-7). http://www.worldcat.org/title/mmpds-07-metallic-materials-properties-development-and-standardization-mmpds/oclc/822046630

Notations

a real constant

b real constant

ksi kilo pound/in^{2}, (1 ksi = 6.895 MPa)

t time

A amplitude of cyclic stresses

D, D' lines representing crack initiation domain for unnotched specimens/

faraway indicators for crack initiation in notched specimens

F_{tu}/F_{cu} ultimate tensile/compressive stress

F_{ty}/F_{cy} yield tensile/compressive stress

L material longitudinal direction

LT material long transversal direction

N_{f}_{ }number of cycles to fatigue

R cycle stress ratio R = σ_{min}/σ_{max}_{ }

R_{1}, R_{2} stress ratios of the known best-fit curves

δR R_{2} − R_{1} (R_{2} > R_{1}), distance from R to the closest known curve

τ_{ }microscopic shear stress

σ microscopic stress tensor

σ_{min} minimum cycle stress (macroscopic)

σ_{max}_{ }maximum cycle stress (macroscopic)

σ_{mean}_{ }mean cycle stress (σ_{max} + σ_{min})/2

Σ macroscopic stress tensor

ω pulsation

ph microscopic hydrostatic pressure

PH macroscopic hydrostatic pressure

s microscopic deviatoric stress tensor

S macroscopic deviatoric stress tensor

q^{*} deviatoric component of the stabilized microscopic residual stress tensor

Anld. annealed form

Norm. normalized form

Wrght. Wrought form (drawn, rolled or extruded)