How to Obtain Fatigue Model Parameters

December 26, 2014

When simulating fatigue, you are faced with two main challenges. The first is to select a suitable fatigue model for your application and the second is to obtain the material data for the selected model. I recently addressed the first challenge in the blog post “Which Fatigue Model Should I Choose?“. Today, I will address the second challenge and discuss how you can obtain fatigue model parameters.

Predict Fatigue Using Many Different Models

Fatigue models are based on physical assumptions and are therefore said to be phenomenological. Since different micromechanical mechanisms govern fatigue under various conditions, many analytical and numerical relations are needed to cover the full spectrum of fatigue. These models, in turn, require dedicated material parameters.

It is well known that fatigue testing is expensive. Many test specimens are necessary since the impurities responsible for fatigue initiation are randomly distributed in the material. The difference in the fatigue life is clearly visible when you visualize all the test results in an S-N curve.

A picture of an S-N curve.
An S-N curve. The black squares represent individual fatigue tests.

Advice for Obtaining Model Parameters via the S-N Curve

Since the S-N curve — also called the Wöhler curve — is one of the oldest tools for fatigue prediction, there is a good chance that the material data is already available in this form. Many times, the data is given for a 50% failure risk. If you do not have access to the material data, you are faced with a testing campaign.

When you are done, pay attention to the statistical aspect and, at each load level, select the same reliability when constructing an S-N curve. This is important since the S-N curve is expressed in a logarithmic scale where a small difference in input has a large influence on the output. Then, S-N curves for different reliability levels fall under each other and you should select an appropriate level for your application. For noncritical structures, a failure rate of 50% might be acceptable. However, for critical structures, a significantly lower failure rate should be chosen.

Always pay attention when you combine fatigue data from different sources. Make sure that the testing conditions and the operating conditions are the same.

Advice for Running Fatigue Tests that Consider Mean Stress

Another aspect of fatigue testing considers the mean stress that has a substantial influence on the fatigue life. In general, fatigue tests performed at tensile mean stress will give a shorter life than tests performed at a compressive mean stress. This effect is also frequently expressed using the R-value (the ratio between the minimum and maximum stress in the load cycle). Thus, with decreasing mean stress (or R-value), the fatigue life increases.

In the Fatigue Module, the Stress-Life models do not take into account this effect. When using these models, you need to choose material data obtained under the same testing conditions as the operating one.

In the cumulative damage model, the Palmgren-Miner linear damage summation uses an S-N curve. However, in this model, the S-N curve is specified with the R-value dependence and the mean stress effect is accounted for.

A schematic highlighting the mean stress effect.
The mean stress effect.

In case you use a material library and the fatigue data is specified using the maximum stress, you can easily convert it to the stress amplitude using

\sigma_a=\frac{\sigma_{\textrm{max}}(1-R)}{2}

where \sigma_a is the stress amplitude, \sigma_{max} is the maximum stress, and R is the R-value.

Advice for Obtaining Parameters for Findley and Matake Critical Plane Models

The stress-based models seem to be fairly simple. For example, the Findley and the Matake models use the expressions

\left(\frac{\Delta\tau}{2}+k\sigma_{\textrm{n}}\right)_{\textrm{max}} =f

and

\left(\frac{\Delta\tau}{2}\right)_{\textrm{max}}+k\sigma_{\textrm{n}} =f

respectively. They depend on only two material constants: f and k. These material parameters are, however, nonstandard material data that can be related to the endurance limit of the material.

Note that the actual values of f and k differ between the two models. The analytical relation is somewhat cumbersome to obtain since the stress-based models are based on the critical plane approach and you need to find a plane where the left-hand sides of the above relations are maximized. This is basically done by expressing the shear and the normal stress as a function of the orientation using the Mohr’s stress circle, maximizing by setting the derivative to zero, and simplifying the resulting relation.

The different steps of the data manipulation will not be shown here. For the Findley model, the material parameters are related to the standard fatigue data using

\frac{f}{\sigma_U(R)} = \frac{(1-R)^2+2k\beta+4k^2}{2\beta(1-R)},\ \ \ \beta=\sqrt{(2k)^2+(1-R)^2}

Here, R is the R-value and \sigma_U(R) is the endurance limit. The argument of the endurance limit indicates that the stress is R-value dependent. For the Matake model, the relation is somewhat simpler and given by

\frac{f}{\sigma_U(R)}=0.5+\frac{k}{1-R}

Since both relations have two unknown material parameters, you need endurance limits from two different types of fatigue tests. To illustrate this, consider a case where one endurance limit is obtained by alternating the load between a tensile and a compressive value, R=-1. In the second case, the load is cycled between a zero load and a maximum load, R=0. For the Findley model, this leads to

\left\{
\begin{array}{lr}
\frac{f}{\sigma_U(-1)}=\frac{1}{2}\left(k+\sqrt{1+k^2}\right)\\
\frac{f}{\sigma_U(0)}=\frac{1}{2}\left(2k+\sqrt{1+4k^2}\right)
\end{array}
\right.

The pair of equations must be solved numerically. Here is the strategy:

  • Eliminate f between the two equations. This is trivial since it always appears as a linear term.
  • Now, you have a nonlinear equation for k only. Since k has a rather small variation (usually between 0.2 and 0.3), it is easy to solve even by pure trial and error.
  • Given the computed k, evaluate f using either of the original equations.

For the Matake model, the two fatigue tests lead to

\left\{
\begin{array}{lr}
\frac{f}{\sigma_U(-1)}=\frac{1}{2}+\frac{k}{2}\\
\frac{f}{\sigma_U(0)}=\frac{1}{2}+k
\end{array}
\right.

which you can solve analytically.

Fatigue Model Examples

I would like to share a few examples where the discussed fatigue models are used:


Comments (0)

Leave a Comment
Log In | Registration
Loading...
EXPLORE COMSOL BLOG