Multiaxial HCF Models
Multiaxial HCF Models
Fatigue
71
72 CHAPTER 6. FATIGUE
max(a,n + kn ) = f, (6.1)
where a,n is the shear amplitude on a plane with a unit normal n, i.e. a,n = n /2, and
n is the normal stress on the same plane. The two material parameters k and f can be
determined from two fatigue loading tests, which can be for example
alternating normal stress, m = 0, af = a,R=1 ,
pulsating normal stress m = a,R=0 .
Other test combinations are also possible, like loading under alternating normal stress and
alternating shear stress.
The idea of Findleys fatigue criterion on the critical plane is illustrated in Fig. 6.1.
Alternating normal stress. In this case the uniaxial normal stress in the direction of the
loading can be expressed as (t) = a sin t and thus in an arbitrary plane inclined by an
angle from the loading direction
Both the normal stress and shear stress attain their maximum at the same time. It is also
observed that in this alternating loading case a,n = n , thus
The subscript 1 indicates that the function refers to the alternating loading case. For the
existence of an extreme value of a function, the derivative has to vanish at the extremum
point, therefore
dg1
= cos2 sin2 2k sin cos = cos 2 k sin 2 = 0, (6.8)
d
giving the solution of the critical plane angle for the alternating loading case
1
tan 21 = . (6.9)
k
If the stress a is at the fatigue limit, i.e. a = 1 , the value for the f parameter is
obtained
Pulsating normal stress. In the case of pulsating normal stress m = a and the stress
history can be written as (t) = a (1cos t). As in the case of atrenating normal stress,
the shear stress and normal stress are in the same phase, thus attaing their maximum at
the same time instant. Now the maximum normal stress is n = 2a and a,n = a , thus
Denoting the angle dependent term as g0 () = sin cos + 2k cos2 , the condition for
the extremum is
dg0
= cos2 sin2 4k sin cos = cos 2 2k sin 2 = 0, (6.12)
d
which gives the solution of the critical plane angle for the pulsating normal stress case
1
tan 20 = . (6.13)
2k
In both cases (6.6) and (6.11), the parameter f is the same on the critical planes, which
gives the equality
Terms depending on angles 1 and 0 can be eliminated using the equations (6.9) and
(6.13).
For the elimination, the following trigonometric identities are usefull:
tan 2 1
tan = , and cos2 = . (6.16)
1 + 1 + tan2 2 1 + tan2
For the alternating stress case
tan 21 1/k 1
tan 1 = p = p = , (6.17)
1 + 1 + tan2 21 1 + 1 + 1/k 2 k + 1 + k2
1 (k + 1 + k 2 )2
cos2 1 = = . (6.18)
1 + tan2 1 1 + (k + 1 + k 2 )2
Thus
f = 1 cos2 1 (tan 1 + k) = 21 1 (k + 1 + k 2 ). (6.19)
Alternatively for the pulsating case, it is obtained
tan 20 1/2k 1
tan 0 = p = p = p , (6.20)
1 + 1 + tan2 20 1 + 1 + 1/(2k)2 2k + 1 + (2k)2
p
2 1 (2k + 1 + (2k)2 )2
cos 0 = = p . (6.21)
1 + tan2 0 1 + (2k + 1 + (2k)2 )2
Therefore p
f = 0 cos2 0 (tan 0 + 2k) = 21 0 (2k + 1 + (2k)2 ). (6.22)
Finally from (6.15), an equation
a,R=0 0 k + 1 + k2
= = p (6.23)
a,R=1 1 2k + 1 + (2k)2
is obtained. Denoting the fatigue stress ratio 0 /1 = , the parameter k can be found
as a function of the ratio as
1
(1 2 )
k=p 2 . (6.24)
(5 2 2 2 )
A typical range of value for metals is 0.2 k 0.3. Notice that when 1/2 then
k . When = 0.5 the pulsating and that the mean normal stress has no effect on
6.2. MULTIAXIAL HCF MODELS 75
1 45
40 alternating
0.8 35
30
0.6
25
k
20
0.4
15
0.2 10 tensile pulsating
5
0 0
0.5 0.6 0.7 0.8 0.9 1 0.5 0.6 0.7 0.8 0.9 1
0 /1 0 /1
(a) (b)
Figure 6.2: Findleys fatigue criterion. (a) Material parameter k as a function of the
fatigue stress ratio and (b) the angle of the critical plane in alternating and in tensile
pulsating normal stress states.
the fatigue limit, which is unrealistic. In the Haigh diagram this means that the fatigue
limit is a horizontal line. The value of k as a function of the fatigue stress ratio 0 /1 is
shown in Fig.6.2a. In addition the critical plane angles in both alternating and pulsating
tensile normal stresses as a function of the fatigue stress ratio are shown in Fig. 6.2b.
The parameter f can be obtained either from (6.19) or (6.22), i.e.
p
f = 21 1 (k + 1 + k 2 ) = 21 0 (2k + +(2k)2 ). (6.25)
It can be shown that the f -parameter can be given also using the fatigue stress amplitude
in pure shear stress as
f = 1 + k 2 1 , (6.26)
where 1 is the shear fatigue stress in alternating shear test. Derivation of this result is
left as an exercise for an interested reader. The ratio between alternating shear and normal
stress fatigue limits predicted by the Findleys model is therefore
1 k + 1 + k2
= . (6.27)
1 2 1 + k2
This relationship is shown in Fig. 6.3 as a function f the 0 /1 ratio.
The Findley model can also be used for finite life analysis. In this case the f -parameter
is replaced by
f = 1 + k 2 1 N b , (6.28)
where N is the number of cycles and b is an additional material parameter.
76 CHAPTER 6. FATIGUE
0.9
0.8
1 /1
0.7
0.6
0.5
0.5 0.6 0.7 0.8 0.9 1
0 /1
Figure 6.3: Alternating shear stress fatigue limit as a function of the fatigue stress ratio
0 /1 according to the Findley fatigue criterion.
where (t) and h (t) are the instantaneous shear and hydrostatic stresses on the slip band.
If (t) + ah (t) < b, t the structure will elastically shake-down to the applied loading.
The macro- and microscale stresses are related as
where is the stabilised residual stress tensor. During the elastic shakedown the mi-
croscale yield surface will expand and move according to isotropic and kinematic harden-
ing laws and the stabilized residual stress tensor designates the distance from the center
of the yield surcafe from the origin. The state of shakedown is found if the stress path
is completely inside the expanded and displaced yield surface. On the meridia plane the
situation is illustrated in Fig. 6.4.
6.2. MULTIAXIAL HCF MODELS 77
Figure 6.4: Dang Van fatigue failure criterion in the meridian plane of the micro-stress.
Safe stress path shown by solid green line and stress path resulting to fatigue fracture by
a dashed blue/red line.
In proportional loading the residual stress tensor can be computed from the opposite
value of the average of two extreme macroscopic stress tensors
1
= (S (t1 ) + S (t2 )) (6.31)
2
where t1 and t2 designates the two time instants at which the extreme values of the macro-
scopic stress S is obtained.
S2 S3 S2 S3
(a) (b)
Figure 6.5: Cyclic (a) alternating stress state and (b) tensile pulsating stress state in the
macroscopic stress principal stress space.
Both of these cyclic stress paths are illustarated in the principal stress space in Fig. 6.5.
It is clear that the stabilized residual stress tensor vanish in the alternating uniaxial
notrmal stress state and in the uniaxial pulsating tensile stress state it has only one nonzero
element 11 .
In the alternating uniaxial normal stress state the maximum shear stress is max =
1
2 1
and the hydrostatic stress at the same time instant is h = 13 1 . Thus, at the fatigue
limit the following equation has to be satisfied
1
2 1
+ a 31 1 = b. (6.35)
In the tensile pulsating uniaxial stress state the component 11 is clearly the mean
macroscopic stress S11,m . It can also be deduced from (6.31) where
20 0 0 0 0 0
S (t1 ) = 0 , and S (t2 ) = 0 0 0 , thus = 0 0 0 .
0 0 0 0 0 0
(6.36)
Deviatoric part of the stabilised residual stress tensor is
2
3 0 0 0
1
dev = 0
3 0
0 , (6.37)
1
0 0 3 0
and the micro-stress tensor is
1
3 0
0 cos t 0 0
1
(t) = S (t) + dev = 0
3 0
0 . (6.38)
1
0 0
3 0
6.2. MULTIAXIAL HCF MODELS 79
0.9
0.8
1 /1
0.7
0.6
0.5
0.5 0.6 0.7 0.8 0.9 1
0 /1
Figure 6.6: Alternating shear stress fatigue limit as a function of the fatigue stress ratio
0 /1 according to the Dang Van fatigue criterion.
and the most dangerous maximum values occur at the same time instant and are max max =
1
and max h = 32 0 . Therefore inserting these values to the fatigue criterion (6.29),
2 0
gives
1
+ a 23 1 = b.
2 0
(6.40)
Putting together (6.35) and (6.40), the a coefficient can be solved, resulting in expression
3(1 0 /1 ) 3(1 )
a= = . (6.41)
2(20 /1 1) 2(2 1)
Then, the b parameter can be solved either from (6.35) or (6.40), giving
0 0 1
b= = = . (6.42)
2(20 /1 1) 2(2 1) 2(2 1)
It can be noticed that the b parameter equals with the alternating shear fatigue limit 1 .
Dependency of alternating shear fatigue stress on the fatigue stress ratio 0 /1 is shown
in Fig. 6.6.
The cyclic uniaxial stress states in the medidian plane of the microtress is illustrated
in Fig. 6.7.
80 CHAPTER 6. FATIGUE
max
A ( 31 1 , 12 1 )
B ( 23 0 , 12 0 )
h
C ( 31 0 , 0)
Figure 6.7: Determining the material parameters using the Dang Van criterion. Fatigue
limit is shown by a red line and the uniaxial loading cases by blue lines. To fix the
failure line two tests have to be conducted. In the text the failure line is determined
from alternating and tensile pulsating normal stress tests, which are shown by blue solid
lines and the fatigue stresses from tests as solid black circles, points A and B. Alternative
uniaxial tests are alternating shear and compressive pulsating tests shown by blue dashed
lines. Notice that the meridian plane is for the microstress.
.
A
1
where s = 31 tr()I is the deviatoric stress tensor. Shape of the endurance surface in
the deviatoric plane is circular and the meridian planes are straight lines as in the case of
the Drucker-Prager model in plasticity. The center point in the deviatoric plane is defined
by the -tensor, which memorizes the load history.
The slope in the Haigh diagram A can be determined if the fatigue stress in tensile
pulsating loading is known as
1 0
A= . (6.45)
0
For the evolution of the deviatoric back-stress tensor , a hardening rule similar to
Zieglers kinematic hardening rule in plasticity theory is adopted
= C(s ), (6.46)
where C is a non-dimensional material parameter and the dot denotes the time rate. Evo-
lution of the back-stress (6.46) and damage (6.48) takes place only if the following con-
ditions are satisfied
0, and 0, (6.47)
which mean that damage accumulation can only occur when the stress is moving away
from the endurance surface. In contrast to plasticity, the stress state can lie outside the
endurance surface. The idea for the endurance surface movement is illustrated in Fig. 6.9.
Material damage is described with a scalar variable D [0, 1], for which the evolution
is governed by equation of the form
1 1
ds
>0 >0
0 ds < 0
s 6= 0 s = 0
d D 0 D = 0
2 (a) 3 2 (b) 3
Figure 6.9: Ottosens HCF model. (a) Movement of the endurance surface and damage
growth when the stress is outside the endurance surface and moving away from it. (b)
When the stress is outside the endurance surface, damage and back stress does not evolve.
Ottosen et al. [17] choose a simple exponential form for the damage evolution law
D = K exp(L). (6.49)
[1] V.V. Bolotin. Mechanics of Fatigue. CRC Mechanical Engineering Series. CRC Press, Boca
Raton, 1999.
[2] K. Dang Van. Macro-micro approach in high-cycle multiaxial fatigue. In D.L. McDow-
ell, editor, Advances in Multiaxial Fatigue, number 1191 in ASTM STP, pages 120130,
Philadelphia, 1993. Americal Society for Testing and Materials.
[3] K. Dang Van, G. Gailletaud, G. Flavenot, A. Le Douaron, and H.P. Lieurade. Criterion
for high cycle fatigue failure under multiaxial loading. In M.W. Brown and K.J. Miller,
editors, Biaxial and Multiaxial Fatigue, number 3 in EGF, pages 459478, London, 1989.
Mechanical Engineering Publications.
[4] W.N. Findley. A theory for the effect of mean stress on fatigue of metals under combined
torsion and axial load or bending. Journal of Engineering for Industry, pages 301306,
November 1959.
[5] W. Flgge. Viscoelasticity. Blaisdell Publishing Company, 1967.
[6] G.A. Holzapfel. Nonlinear Solid Mechanics - A Continuum Approach for Engineering. John
Wiley & Sons, 2000.
[7] J. Lemaitre. How to use damage mechanics. Nuclear Engineering and Design, 80(2):233
245, 1984.
[8] J. Lemaitre and J.-L. Chaboche. Mechanics of Solid Materials. Cambridge University Press,
1990.
[9] J. Lemaitre and R. Desmorat. Engineering Damage Mechanics, Ductile, Creep, Fatigue and
Brittle Failures. Springer-Verlag, Berlin, Heidelberg, 2005.
[10] J. Liu and H. Zenner. Fatigue limut of ductile metals under multiadial loading. In A Carpin-
teri, M. de Freitas, and A. Spagnoli, editors, Biaxial/multiaxial fatigue and fracture, pages
147163. Elsevier Science Ltd. and ESIS, 2003.
[11] J. Lubliner. Plasticity Theory. Pearson Education, Inc., 1990.
[12] L.E. Malvern. Introduction to the Mechanics of a Continuous Medium. Prentice Hall, En-
glewood Cliffs, New Jersey, 1969.
83
84 BIBLIOGRAPHY
[13] S. Murakami. Continuum Damage Mechanics, volume 185 of Solid Mechanics and Its Ap-
plications. Springer Netherlands, 2012.
[14] Y. Murakami. Metal Fatigue, Effects of Small defects and Nonmetallic Inclusions. Elsevier
Science, 2002.
[16] N.S. Ottosen and M. Ristinmaa. The Mechanics of Constitutive Modeling. Elsevier, 2005.
[17] N.S. Ottosen, R. Stenstrm, and M. Ristinmaa. Continuum approach to high-cycle fatigue
modeling. International Journal of Fatigue, 30(6):9961006, June 2008.
[18] I. V. Papadopoulos. Long life fatigue under multiaxial loading. International Journal of
Fatigue, 23(10):839 849, 2001.
[20] D.F. Socie and G.B. Marquis. Multiaxial Fatigue. Society of Automotive Engineers, Inc.,
Warrendale, Pa, 2000.
[21] A.J.M. Spencer. Continuum Mechanics. Dover Publications, Inc., 2004. First published by
Longman Group UK Limited in 1980.
[22] S. Suresh. Fatigue of Materials. Cambridge University Press, 2nd edition, 1998.