Fatigue Design Criteria for Machine Elements
Fatigue Design Criteria for Machine Elements
For questions : Recep Durna Problem 2: A flat steel part machined from a plate, having a rectangular cross-
Room : D-112/D-103 Phone: 5216 section and made of AISI 1050 CD steel will be subjected to load F, which
E-mail : rdurna@[Link] Date: 25/11/2025 fluctuates between 4kN and 10kN, in a working environment of 100°C. Determine
the minimum thickness of the part for an infinite operation life taking factor of
ME 307 – MACHINE ELEMENTS I
safety as 2 and reliability as 0.95 according to the modified Goodman line. All
TUTORIAL 8
“FATIGUE DESIGN CRITERIA” dimensions are in mm.
R25
A
A
Figure 2
Figure 1
Figure 4
Stress Concentration Factor where Kf is fatigue stress concentration factor and To remember values
this factor is critical here because of obvious radius given to the work piece.
𝑆𝑦 = 400𝑀𝑃𝑎 𝑆𝑒 = 184.5𝑀𝑃𝑎 𝑆𝑢𝑡 = 550𝑀𝑃𝑎 𝑛=2
𝑟/𝑑 = 0.15
} ⇒ 𝐾𝑡 ≅ 1.5 (Figure A-15-9, the formula sheet)
𝐷/𝑑 = 1.5 𝜎𝑎′′ = 37.9 𝑀𝑃𝑎 and ′′
𝜎𝑚 = 162.1 𝑀𝑃𝑎
q 0.8 (Figure 6.20, the formula sheet) 𝛼 = atan(𝜎𝑎′′ /𝜎𝑚
′′ )
= 13.2°
This leads to 𝐾𝑓 = 1.4
Since 90° > 𝛽 > 𝛼; failure will be due to fatigue. Therefore, fatigue failure criteria
𝜎𝑚𝑎𝑥 + 𝜎𝑚𝑖𝑛 should be used.
′
𝜎𝑚 = 𝐾𝑓 = 1.778 ⋅ 106 ⋅ 𝑀
2
𝜎𝑚𝑎𝑥 − 𝜎𝑚𝑖𝑛 𝜎𝑎′ 𝜎𝑚
′
1
𝜎𝑎′ = 𝐾𝑓 = 5.32 ⋅ 106 ⋅ 𝑀 + =
2 𝑆𝑒 𝑆𝑢𝑡 𝑛
𝑀 ≅ 15.6 𝑁 ⋅ 𝑚
𝑑ℎ𝑜𝑙𝑒 3 Solving:
= = 0.1 ⇒ 𝐾𝑡 ≅ 2.25 (Figure A − 15 − 11)
𝐷 30 𝑆𝑦 𝜎𝑎′′ 𝜎𝑚
′′
1
𝜎𝑎′′ + 𝜎𝑚
′′
= and + = for 𝜎𝑎′′ , 𝜎𝑚
′′
𝑛 𝑆𝑒 𝑆𝑢𝑡 𝑛
𝐾𝑓 = 1 + 𝑞 ⋅ (𝐾𝑡 − 1)
To remember values
𝐾𝑓 = 1 + 0.8 ⋅ (2.25 − 1) = 2
𝑆𝑦 = 400𝑀𝑃𝑎 𝑆𝑒 = 150.5𝑀𝑃𝑎 𝑆𝑢𝑡 = 550𝑀𝑃𝑎 𝑛=2
′
𝜎𝑚𝑎𝑥 + 𝜎𝑚𝑖𝑛 2𝑀 𝜎𝑎′′ = 28.3 𝑀𝑃𝑎 and ′′
𝜎𝑚 = 171.7 𝑀𝑃𝑎
𝜎𝑚 = 𝐾𝑓 =
2 2.2 ⋅ 10−6
𝛼 = atan(𝜎𝑎′′ /𝜎𝑚
′′ )
= 9.36°
𝜎𝑚𝑎𝑥 − 𝜎𝑚𝑖𝑛 6𝑀
𝜎𝑎′ = 𝐾𝑓 =
2 2.2 ⋅ 10−6 Since 90° > 𝛽 > 𝛼 ; failure will be due to fatigue. Therefore fatigue failure criteria
shuld be used.
Slope of the load line is the same as in previous section. 𝜎𝑎′ 𝜎𝑚 ′
1
+ =
𝑆𝑒 𝑆𝑢𝑡 𝑛
𝛽 = atan (𝜎𝑎′ /𝜎𝑚
′
) = 71.6 °
𝑃𝑚𝑖𝑛 = 4 ⋅ 103 𝑁, 𝑃𝑚𝑎𝑥 = 10 ⋅ 103 𝑁 𝑀𝑚𝑎𝑥 ⋅ (𝑟𝑛 − 𝑟𝑖 ) 1500𝑁𝑚 ⋅ (0.0455 − 0.025) 5467 ⋅ 103
𝜎𝑏,𝑚𝑎𝑥 = = =
𝑀𝑚𝑖𝑛 = 4 ⋅ 103 ⋅ 0.15 = 600𝑁𝑚 (ℎ ⋅ 𝑡) ⋅ 𝑒 ⋅ 𝑟𝑖 (0.05 ⋅ 𝑡) ⋅ 0.0045 ⋅ 0.025 𝑡
𝑀𝑚𝑎𝑥 = 10 ⋅ 103 ⋅ 0.15 = 1500𝑁𝑚 For AISI 1050 CD steel: 𝑆𝑢𝑡 = 590𝑀𝑃𝑎 and 𝑆𝑦 = 490𝑀𝑃𝑎
At section A-A, the part is subjected to both axial and bending stresses. The Then 𝑆𝑒′ = 0.5 × 𝑆𝑢𝑡 → 𝑆𝑒′ = 0.5 × 590 = 295 𝑀𝑃𝑎
maximum and the minimum values of these stresses can be obtained in terms
of thickness t. The modifying factors for the endurance limit of the part can be determined. The
procedure follows the same logic of the previous question.
This gives value around 0.81 to 1. Let’s assume 0.85 as it is expected to be close 1640 ⋅ 103 60000 1700 ⋅ 103
𝜎𝑎′ = 1 . 𝜎𝑏−𝑎 + 1 . 𝜎𝑎𝑥−𝑎 = + =
to maximum limit (51 mm) and check here later. For axial stress kb is 1. 𝑡 𝑡 𝑡
140000 3830 ⋅ 103 3970 ⋅ 103
𝜎𝑚 ′ = 𝜎𝑎𝑥−𝑚 + 𝜎𝑏−𝑚 = + =
Loading factor: For bending: 𝑘𝑐 = 1. For axial load 𝑘𝑐 = 0.85 𝑡 𝑡 𝑡
Temperature Factor:𝑘𝑑 = 1.03 (for T=100°C, from figure 2.17 or Equation 6-
26) Note that the load factors are taken into account during the determination of
Reliability Factor: 𝑘𝑒 = 0.868 (for a reliability of 0.95; Table 6-5) alternating stresses since the values for axial and bending cases are different.
Then the total bending and alternating stresses are
Then, the endurance limit of the machine member can be calculated.
𝑆𝑒 = 𝑘𝑎 ⋅ 𝑘𝑏 ⋅ 𝑘𝑐 ⋅ 𝑘𝑑 ⋅ 𝑘𝑒 ⋅ 𝑆𝑒′
Then 𝑆𝑒 = 0.773 ⋅ 0.85 ⋅ 0.868 ⋅ 1.03 ⋅ 1 ⋅ 295 = 173.9 𝑀𝑃𝑎
The mean and alternating stresses for axial and bending loading are (for ductile
material):
𝜎𝑎𝑥−𝑚𝑎𝑥 + 𝜎𝑎𝑥−𝑚𝑖𝑛 1 200000 80000 140000
𝜎𝑎𝑥−𝑚 = = ( + )=
2 2 𝑡 𝑡 𝑡
𝜎𝑎𝑥−𝑚𝑎𝑥 − 𝜎𝑎𝑥−𝑚𝑖𝑛 1 200000 80000 60000
𝜎𝑎𝑥−𝑎 = = ( − )=
2 2 𝑡 𝑡 𝑡
𝑘𝑐 = 1.24𝑑−0.107 = 0.854
𝑆𝑒
𝜎𝑎′ Since 7.62 𝑚𝑚 < 32𝑚𝑚 < 51𝑚𝑚 and 𝑘𝑐 assumption is close to obtained value
of 𝑘𝑐 , the assumption is correct.
𝜎𝑎′′ 𝛽
𝛼
Solution 3:
0 ′
′′
𝜎𝑚 𝑆𝑦 𝑆𝑢𝑡 𝜎𝑚
𝜎𝑚 At first, one should perform the force analysis of the system. During operation,
torque on the shaft is steady, so does the torsional stress on the shaft. Also the forces
in y and z directions are steady. However, because of the rotation of the shaft
Slope of the load line is bending stresses on the shaft do change periodically.
𝜎𝑎′
𝛽 = tan−1 ( ′ ) = 23.18° T1 and T2 is found from the torque balance on the shaft ( T = 0 ) by using the
𝜎𝑚
relation between them:
The intersection point of the modified Goodman and the Yield lines is:
𝑑𝐴 𝑑𝐵
∑ 𝑇 = 0 : (𝑇2 − 𝑇1 ) = (270 − 50) = 220(0.15) = 33 𝑁. 𝑚
2 2
𝑆𝑦 𝜎𝑎 ′′ 𝜎 ′′ 1
𝜎𝑎 ′′ + 𝜎𝑚 ′′ = 𝑛
and 𝑆𝑒
+ 𝑆𝑚 =𝑛 It is given that force on the loose side of the pulley is 15% of the force on the
𝑢𝑡
tight side. Knowing T2>T1:
To remember values
𝑆𝑦 = 490𝑀𝑃𝑎 𝑆𝑒 = 173.9𝑀𝑃𝑎 𝑆𝑢𝑡 = 590𝑀𝑃𝑎 𝑛=2 dA
(T2 − 0.15T2 ) = 33 N ⋅ m
2
33 N. m
𝜎𝑎′′ = 20.9 𝑀𝑃𝑎 and 𝜎𝑚
′′
= 224 𝑀𝑃𝑎 0.85T2 = → T2 = 311 N and T1 = 46.6 N
0.25 𝑚
Since 90° > 𝛽 > 𝛼; the failure will be due to fatigue. Therefore fatigue failure 2
criteria should be used. The next step is to obtain the moment diagram. The forces act in y and z directions
on the shaft. The reactions at bearings also act in these directions. It is needed to
𝑆𝑎′ 𝑆′ 1 1700⋅103 /𝑡 3970⋅103 /𝑡 1
+ 𝑆𝑚 = 𝑛 ➔ 173.9 .103 𝑃𝑎 + =2 analyze the forces on the shaft in x-z and x-y planes separately.
𝑆𝑒 𝑢𝑡 590.103 𝑃𝑎
-107N -48.9N.m
-174N
Torque on the shaft is:
Moment 26.3N.m Torque
33N.m
x
-32.1N.m
Endurance limit of shafts used in machine elements is found by finding the where, 𝜏𝑎 = 0 𝑎𝑛𝑑 𝜎𝑎,𝑎𝑥𝑖𝑎𝑙 = 0 → 𝜎𝑎′ = 𝐾𝑓,𝑏𝑒𝑛𝑑𝑖𝑛𝑔 𝜎𝑎
modification factors:
′ 2 2
𝜎𝑚 = √(𝐾𝑓,𝑏𝑒𝑛𝑑𝑖𝑛𝑔 ⋅ 𝜎𝑚,𝑏𝑒𝑛𝑑𝑖𝑛𝑔 + 𝐾𝑓,𝑎𝑥𝑖𝑎𝑙 ⋅ 𝜎𝑚,𝑎𝑥𝑖𝑎𝑙 ) + 3(𝐾𝑓𝑠,𝑡𝑜𝑟𝑠𝑖𝑜𝑛 𝜏𝑚 )
Table 6-2 gives 𝛼 = 38.6 and 𝑏 = −0.650. Using Equation (6-18)
′
where, 𝜎𝑚,𝑏𝑒𝑛𝑑𝑖𝑛𝑔 = 0, 𝜎𝑚,𝑎𝑥𝑖𝑎𝑙 = 0 → 𝜎𝑚 = √3 ⋅ 𝐾𝑓𝑠 ⋅ 𝜏𝑚
𝑘𝑎 = 𝑎(𝑆𝑢𝑡 )𝑏 = 38.6 ⋅ 470−0.650 (6-18)
= 0.707
For the critical differential element at A:
Since the diameter is not known, assume that it is in the range 7.62<d<51mm, 𝑑
𝑀𝐴 𝑐 1.6 𝑀𝐴 953.4𝑁𝑚
then let’s assume kb=0.85 as we did in previous question. 𝜎𝑎′= 𝐾𝑓,𝑏𝑒𝑛𝑑𝑖𝑛𝑔 𝜎𝑎 = 1.6 ⋅ ( ) = 1.6 ⋅ 2 4 = ⋅ 32 ⋅ 3 =
𝑘𝑐 = 1 as cycling stress is bending. (6-25) 𝐼 𝜋𝑑 𝜋 𝑑 𝑑3
𝑘𝑑 = 1 as no temperature is given 64
′
𝑇𝑐
𝑘𝑒 = 0.753 for 99.9% reliability 𝜎𝑚 = √3 ⋅ 1.3 ⋅ 𝜏𝑚 = √3 ⋅ 1.3 ⋅
𝐽
𝑑
√3 𝑇 378.4𝑁𝑚
𝑆𝑒 = 𝑘𝑎 𝑘𝑏 𝑘𝑐 𝑘𝑑 𝑘𝑒 𝑆𝑒′ = 0.707 ⋅ 0.85 ⋅ 1 ⋅ 1 ⋅ 0.753 ⋅ 0.5 ⋅ 470𝑀𝑃𝑎 𝜎𝑚′
= √3 ⋅ 1.3 ⋅ 𝑇 ⋅ 2 4 = ⋅ 1.3 ⋅ 16 ⋅ 3 =
𝜋𝑑 𝜋 𝑑 𝑑3
= 106.34𝑀𝑃𝑎 32
Note that √3 and 1.6 constants in their respective equations are not of equivalent
purpose. Former is from Von mises while latter is from stress modifiers.
Although the slope of the load line is high enough, it is better to check the 𝑆𝑒 ′ = 0.5(620) = 310𝑀𝑃𝑎
intersection point of the Yield line and the modified Goodman line: For the modifying endurance limit, modifying factors are to be determined.
Equation of the Yield line: 'a + 'm = Sy
From Equation (6-18) and Table 6-2, for forged steel
' ' 𝑘𝑎 = 54.9 ⋅ 620−0.758 = 0.454
Equation of the Goodman line: a + m = 1
Se Sut From (6-20), since the bar in question is bending while rotating, there is no need to
Intersection point is: calculate equivalent diameter. (it writes in pages 312-313);
𝑆𝑦 𝑘𝑏 = 1.51 ⋅ 160−0.157 = 0.68 d > 51 mm.
𝜎𝑎′ 𝑆𝑦 − 𝜎′𝑎 1− Note that free body diagram has not been drawn yet so estimated weak spots are
𝑆𝑢𝑡
+ = 1 → 𝜎𝑎′ = = 61.41𝑀𝑃𝑎 not given yet. However, we have no value other than 160mm and even for let’s say
𝑆𝑒 𝑆𝑢𝑡 1 1
−
𝑆𝑒 𝑆𝑢𝑡 200mm, this constant is still over 0.65 so this is reasonable.
𝑘𝑐 = 1 as the bar is bending while 𝑘𝑑 = 𝑘𝑒 = 𝑘𝑓 = 1 as nothing mentioned in the
′
𝜎𝑚 = 𝑆𝑦 − 𝜎𝑎′ = 198.59𝑀𝑃𝑎 question.
𝑆𝑒 = 𝑘𝑎 ⋅ 𝑘𝑏 ⋅ 𝑘𝑐 ⋅ 𝑘𝑑 ⋅ 𝑘𝑒 ⋅ 𝑆𝑒′ = 0.453 ⋅ 0.68 ⋅ 1 ⋅ 1 ⋅ 1 ⋅ 1 ⋅ 0.5 ⋅ 620
Slope of the intersection: = 95.70𝑀𝑃𝑎
𝜎𝑎′ 61.41
tan−1 ( ′ ) = tan
−1
( ) = 17.18° The next step is to estimate the bending stress. Note that, the mean stress due to
𝜎𝑚 198.59
fully reversed bending moment will be zero. The alternating stress due to bending
moment:
Since 17.18°<68.35°, failure is expected to be from fatigue. Thus, the Goodman
line determines the factor of safety.
𝜎𝑎′ 𝜎𝑚′
1
+ = 𝜋𝑑4 𝜋(160)4
𝑆𝑒 𝑆𝑢𝑡 𝑛 𝐼= = = (32.17)106 𝑚𝑚4
64 64
953.4𝑁𝑚 378.4𝑁𝑚 1 1
+ = =
𝑑3 ⋅ 106.34𝑀𝑃𝑎 𝑑3 ⋅ 470𝑀𝑃𝑎 𝑛 3
𝑑 ≥ 30.8𝑚𝑚
𝑘𝑏 = 1.24𝑑−0.107 = 0.859
The guessed 𝑘𝑏 value of 0.85 is not so far away from the obtained value of
0.859.
𝑆𝑓 = 𝑎𝑁 𝑏
𝐹 𝐹
(𝑓𝑆𝑢𝑡 )2 (0.86 ⋅ 620)2
2 2 𝑎= = = 2970
𝑆𝑒 95.70
Figure 5
200 = 2970 ⋅ 𝑁 −0.2486 → 𝑁200 = 51678 𝑐𝑦𝑐𝑙𝑒𝑠
Lifetime, N1, corresponding to 138 MPa in the initial S-N diagram (Figure 7):
d)
The Miner’s rule does not consider the order of cumulative damage effect
introduced by different loadings. So, the order of loading is not important for The
Figure 8 Life reduction after second loading Miner’s rule. However, for The Manson’s rule, the order of application of the
loadings is taken into account. This is because The Manson’s rule takes the damage
of the part after each loading into account by considering the new damaged material
𝑁2 ≅ 173000 𝑐𝑦𝑐𝑙𝑒𝑠
𝑁1 ≅ 98000 𝑐𝑦𝑐𝑙𝑒𝑠