Reliability Analysis of the Shaft Subjected to Twisting Moment and Bending Moment for Normally Distributed
Strength and Stress
Md. Yakoob Pasha1, M. Tirumala Devi2, T. Sumathi Uma Maheswari3.
1, 2, 3Department of Mathematics Kakatiya University, Warangal-506009, Telangana, India. 1mypasha139@gmail.com, 2oramdevi@yahoo.com, 3tsumathiuma@gmail.com
Abstract
Shaft is the rotating component that transmits power from one place to another. The shafts are commonly subject to torsional and bending moments and combinations of these moments. In general, shafts are subjected to a combination of torsional and bending stresses. The design of a shaft is essential, subject to its strength and stress. This paper presents the reliability analysis of the shaft subjected to (a) twisting moment, (b) bending moment and (c) combined twisting and bending moment for which stress and strength follow the normal distribution.
Keywords: Reliability, normal distribution, twisting moment, bending moment, round solid shaft, maximum shear stress theory, maximum normal stress theory.
1. Introduction
Shaft is the most vital component used in almost every mechanical system and machine. Out of all power transmission components, the shaft is the main component that must be designed carefully for the efficient working of the machine. The shafts are designed based on their strength and rigidity considerations. Based strength includes shafts subjected to bending moments, twisting moments, combined bending and twisting moments, and fluctuating loads. While considering rigidity, torsional and lateral rigidity are considered for design. Adekunle A. A. et al. [1] studied the development of CAD software for shafts under various loading conditions. Dr. Edward E. Osakue et al. [2] studied fatigue shaft design verification for bending and torsion. Dr. Edward E. Osakue et al. [3] studied the probabilistic fatigue design of shafts for bending and torsion. Frydrysek K. et al. [4] studied performance-based design applied to a shaft subjected to combined stress. Gowtham et al. [5] studied drive shaft design and analysis. Kamboh, M. S., et al. [7] discussed the design and analysis of the drive shaft with a critical review of advanced composite materials and the root causes of shaft failure. Kececioglu D. et al. [8] studied the reliability analysis of mechanical components and systems. Misra, A., et al. [9] studied the reliability analysis of drilled shaft behaviour using the finite difference method and Monte Carlo simulation. Munoz-Abella, B., et al. [10] discussed the determination of the critical speed of a cracked shaft from experimental data. Nayek et al. [11]
studied the reliability approximation for a solid shaft under a gamma setup. Patel, B., et al. [12] studied a critical review of the design of a shaft with multiple discontinuities and combined loadings. Villa-Covarrubias, B., et al. [14] discussed the probabilistic method to determine the shaft diameter and design reliability.
2. Statistical model
The normal distribution takes the well-known bell shape. This distribution is symmetrical about its mean value. The probability density function for a normally distributed stress x and normally distributed strength ^ is given by [4]
W) =
1
exp
1 X — ßx
ift-N
for — œ < x < œ for - œ < £ < œ
where = mean value of the stress ax = standard deviation of the stress ^ = mean value for the strength a^ = standard deviation of the strength
Let us define y = % — x- It is well known that the random variable y is normally distributed with mean of = ^ — nx and standard deviation of ay = J&2 + ax-The reliability R can be expressed in terms of y as
R = P(y> 0)=—^[
J22n.
c
If we let z = (y — ßy)/ay, then ay dz = dy. When y = 0, the lower limit of z is given by [4]
exp
1(y — Vy
dy
z =
0 — ßy_ fa- ßx)
aï + ai
and y ^ the upper limit of z ^ Then the reliability is given by [4]
=è [ exp{-r
dz
at+ax
(1)
2
2
2
a
f
2
a
y
ÜO
R
The random variable z = (y ^ is the standard normal variable.
ay
a. Shaft subjected to twisting moment only
If the shaft is subjected to a twisting moment only, then torsion equation is given by [6]
T = ai J r
where T = Twisting moment acting upon the shaft
Md. Yakoob Pasha, M. Tirumala Devi, T. Sumathi Uma Maheswari
RELIABILITY ANALYSIS OF THE SHAFT RT&A, No 4 (76)
FOR NORMALLY DISTRIBUTED STRENGTH AND STRESS_Volume 18, December 2023
J = Polar moment of inertia of the shaft about the axis of rotation at = Torsional shear stress and r = Distance from neutral axis to the outer-most fibre
= - where d is the diameter of the shaft
2
For round solid shaft, polar moment of inertia is given by [6]
J 32
From the torsion equation twisting moment is
n
16
The shear stress due to twisting moment is the normally distributed stress
16 XT
T = ^-rXatxd3 and
Then the reliability of the shaft subjected to twisting moment for normally distributed strength and stress is
Rt = — I exp I —) dz
=è / exp(=r)
( 16XT\ v 7
l^-^) (2)
2
Gr +G
f
Shaft subjected to bending moment only
If the shaft is subjected to a bending moment only, then bending equation is given by [6]
I y
where M = Bending moment
I = Moment of inertia of cross-sectional area of the shaft about the axis of rotation ob = Bending stress and
y = Distance from neutral axis to the outer-most fibre For solid shaft, moment of inertia is given by [6]
n d
I = — xd4 and v = ■ 64 2
From the bending equation bending moment is given by [8]
n
M = — xohxd3 32 b
The shear stress due to bending moment is the normally distributed stress
32 x M
Then the reliability of the shaft subjected to bending moment for normally distributed strength and stress is
Rb=2n / exp(-r) dz
( 32 XM\ y '
(3)
c. Shaft subjected to combined twisting moment and bending moment
When the shaft is subjected to combined twisting moment and bending moment, then the shaft must be designed on the basis of the two moments simultaneously.
According to maximum shear stress theory, the maximum shear stress in the shaft is
1
atmax = ~J{Ob)2 + 4<,OtY
where at = shear stress induced due to twisting moment ob = bending stress induced due to bending moment Therefore
16
o.
;^M2+T2
tmax nXd3
The expression lM2 + T2 is known as equivalent twisting moment and is denoted by Te. The equivalent twisting moment may be defined as that twisting moment, which when acting alone, produces same shear stress (at) as the actual twisting moment.
By limiting the maximum shear stress (atmax) equal to the allowable shear stress (ots) for the material is the normally distributed stress
%s
= Otr =■
16 x Te
n x d3
Then the reliability of the shaft subjected to combined twisting moment and bending moment for the normally distributed strength and stress is
72
R«=è S exp{^)dz
w-ozdf) (4)
According to maximum normal stress theory, the maximum normal stress in the shaft is
1 1
°bmax = 2°b + 2^((Jb)2 + 4(cJt)2 where at = shear stress induced due to twisting moment ob = bending stress induced due to bending moment Therefore
32
■\1(M + VM2 + T2)
abmax nxd3[2
The expression 1(M + lM2 + T2) is known as equivalent bending moment and is denoted by Me.
The equivalent bending moment may be defined as that moment which when acting alone produces the same tensile or compressive stress (ab) as the actual bending moment. By limiting the maximum normal stress (abmax) equal to the allowable bending stress (obn) for the material is the normally distributed stress
_ _ 32XMe
VXn = °bn = -¿Xi3
Then the reliability of the shaft subjected to combined twisting moment and bending moment for normally distributed strength and stress is
Rbn=2n S expl^jdz
t 32XMe\ V '
c
CO
3. Numerical results and discussion
a. For twisting moment only
Table-1: d = 110 mm, ^ = 119.6584 N/mm2.
T zt Rt
100000 -1.993596824 0.976901934
2500000 -1.834327549 0.966697306
5000000 -1.659291647 0.951471481
7500000 -1.478628009 0.930380119
10000000 -1.296096181 0.902528825
13000000 -1.079424515 0.859800736
18000000 -0.736112792 0.769168971
20000000 -0.607809459 0.728343073
25000000 -0.313849897 0.623182477
30000000 -0.059258971 0.523627080
Figure 1: Reliability and Twisting Moment
Table-2: T = 104 N-mm, ^ = 119.6584 N/mm2.
d Zt Rt
18 -0.437034097 0.668956690
20 -0.826725375 0.795803632
24 -1.323221625 0.907119157
30 -1.664313463 0.951975098
32 -1.725837528 0.957811677
40 -1.862933958 0.968764221
50 -1.930809365 0.973246684
58 -1.955923837 0.974762937
70 -1.975040282 0.975868212
82
-1.984508267
0.976400396
Diameter of the Shaft
Figure 2: Reliability and Diameter of the Shaft
b. For bending moment only
Table-3: d = 110 mm, ^ = 119.6584 N/mm2.
M zb Rb
100000 -1.98717341 0.976548408
1200000 -1.841158255 0.967200815
2000000 -1.730194966 0.958202276
2400000 -1.673552215 0.952890682
3800000 -1.471338225 0.929400165
4800000 -1.325278842 0.907460658
6500000 -1.079424515 0.859800736
7800000 -0.897319098 0.815225666
8600000 -0.789017221 0.784949029
12000000 -0.369480494 0.644115195
Figure 3: Reliability and Bending Moment
Table-4: M = 104 N-mm, ^ = 119.6584N/mm2.
d zb Rb
22 -0.313849897 0.623182477
26 -0.928635716 0.823461047
32 -1.433116930 0.924087788
36 -1.608706553 0.946159739
40 -1.718965160 0.957189642
52 -1.875536297 0.969640510
62 -1.927365578 0.973032957
82 -1.968900245 0.975517726
100 -1.982910074 0.976311262
110 -1.987173410 0.976548408
Diameter of the Shaft
Figure 4: Reliability and Diameter of the Shaft
c. For combined twisting moment and bending moment
Table-5: M = 104 N-mm, d = 110 mm, ^ = 119.6584 N/mm2.
T zts Rts Rbn
100000 -1.990938603 -1.984506864 0.976756181 0.976400318
2500000 -1.834190835 -1.827344797 0.966687164 0.966176028
5000000 -1.659220259 -1.652076253 0.951464292 0.950740496
7500000 -1.478579423 -1.471289619 0.930373622 0.929393595
10000000 -1.296059727 -1.288770095 0.902522546 0.901260987
13000000 -1.079397107 -1.072274996 0.859794630 0.858201733
18000000 -0.736094575 -0.729543208 0.769163428 0.767165276
20000000 -0.607793789 -0.601533305 0.728337875 0.726257582
25000000 -0.31383893 -0.308363433 0.623178312 0.621097098
30000000 -0.059251123 -0.054549707 0.523623955 0.521751396
1
0.95 0.9 0.85
it
U 0.7
Pi
0.65 0.6 0.55 0.5,
Twisting Moment
Figure 5: Reliability and Twisting Moment
Table-6: T = 104 N-mm, d = 110 mm, ^ = 119.6584 N/mm2.
0.8
"§ 07
0.5
1.5
2
2.5
3
x 10
M zts Rts Rbn
100000 -1.990938603 -1.984506864 0.976756181 0.976400318
1200000 -1.921595252 -1.840874437 0.972671647 0.967180019
2000000 -1.868157213 -1.730019145 0.969129920 0.958186573
2400000 -1.841016165 -1.673403924 0.967190405 0.952876097
3800000 -1.744153719 -1.471242305 0.959433855 0.929387200
4800000 -1.673478046 -1.325202811 0.952883388 0.907448053
6500000 -1.551268837 -1.079369702 0.939581364 0.859788524
7800000 -1.456703987 -0.897275084 0.927400946 0.815213926
8600000 -1.398287099 -0.788978462 0.918986565 0.784937703
12000000 -1.151035866 -0.369456983 0.875141260 0.644106434
1
0.95 0.9
is 0 8
0.75 0.7 0.65
Bending Moment
12
„6
0.85
0
2
4
6
8
10
Figure 6: Reliability and Bending Moment
Table-7: T = 104 N- mm, M = 15 x 104 N-mm, ^ = 119.6584 N/mm2.
d zts zbn Rts Rbn
24 -0.778536422 -0.208090258 0.781873578 0.582420753
28 -1.228383905 -0.802787694 0.890348557 0.788951272
32 -1.491547172 -1.195487826 0.932091052 0.884051755
36 -1.649076771 -1.442464991 0.950434046 0.925414380
40 -1.747806723 -1.599826782 0.959751250 0.945181492
52 -1.888095775 -1.823839541 0.970493453 0.965911833
62 -1.934633552 -1.897558264 0.973482360 0.971122852
82 -1.971987598 -1.956288280 0.975694489 0.974784398
100 -1.984601679 -1.976010500 0.976405597 0.975923206
110 -1.988441857 -1.982002340 0.976618578 0.976260511
Diameter of the Shaft
Figure 7: Reliability and Diameter of the Shaft
Table-8: T = 4x104 N - mm,M = 2 x 104 N-mm, d = 110 mm.
zts Rts Rbn
2.5 -0.521239161 -0.013920898 0.698899912 0.505553455
3.5 -0.918703097 -0.477954094 0.820874555 0.683658561
6.5 -1.425158981 -1.157220313 0.922944375 0.87640882
10.5 -1.652400548 -1.487687204 0.950773538 0.931583298
20.5 -1.826770123 -1.745803964 0.966132830 0.959577489
40.5 -1.913822863 -1.874254579 0.972178603 0.969552328
60.5 -1.942677032 -1.916561189 0.973972403 0.972353149
120.5 -1.971410740 -1.958499434 0.975661543 0.974914281
210.5 -1.983682381 -1.976343042 0.976354371 0.975942031
350.5 -1.990215750 -1.985824871 0.976716413 0.976473613
Mean Strength
Figure 8: Reliability and Diameter of the Shaft
4. Conclusion
The reliability of the shaft is derived from the twisting and bending moments and the combined twisting and bending moments. According to the computations, the reliability of the shaft decreases when the twisting and bending moments increase. Reliability increases when the diameter and strength of the shaft increase.
References
[1] Adekunle, A. A., Adejuyigbe, S. B., & Arulogun, O. T. (2012). Development of CAD software for shaft under various loading conditions. Procedia Engineering, 38, 1962-1983.
[2] Dr. Edward E. Osakue et al., "Fatigue Shaft Design Verification for Bending and Torsion" International Journal of Engineering Innovation and Research, pp.197-206, 2015.
[3] Dr. Edward E. Osakue, "Probabilistic Fatigue Design of Shaft for Bending and Torsion", International Journal of Research in Engineering and Technology, Volume: 03, Issue: 09, pp. 370-386, 2014, DOI: 10.15623/ijret.2014.0309059.
[4] Frydrysek, K. Performance-Based design applied for a Shaft subjected to combined stress.
[5] Gowtham, S., Raaghul, V., Sreeharan, B. N., Kumar, R. P., & Kasim, S. M. (2021, April). Drive Shaft Design and Analysis for Mini Baja. In IOP Conference Series: Materials Science and Engineering (Vol. 1123, No. 1, p. 012015). IOP Publishing.
[6] K. C. Kapur and L. R. Lamberson, "Reliability in Engineering Design", Wiley India Private Limited, New Delhi, 1997.
[7] Kamboh, M. S., Machhi, M. A., & Kamboh, M. F. (2020). Design and Analysis of Drive Shaft with a Critical Review of Advance Composite Materials and the Root Causes of Shaft Failure.
[8] Kececioglu, D. (1972)., Reliability analysis of mechanical components and systems. Nuclear Engineering and Design, 19(2), 259-290.
[9] Misra, A., Roberts, L.A. & Levorson, S.M. Reliability analysis of drilled shaft behaviour using finite difference method and Monte Carlo simulation. Geotech Geol Eng 25, 65-77 (2007). https://doi.org/10.1007/s10706-006-0007-2.
[10] Munoz-Abella, B., Montero, L., Rubio, P., & Rubio, L. (2022). Determination of the Critical Speed of a Cracked Shaft from Experimental Data. Sensors, 22(24), 9777.
[11]Nayek, S., Seal, B., & Roy, D. (2014). Reliability approximation for solid shaft under Gamma setup. Journal of Reliability and Statistical Studies, 11-17.
[12] Patel, B., Prajapati, H. R., & Thakar, D. B. (2014). Critical Review on design of shaft with multiple discontinuities and combined loadings ICCIET-2014.
[13] R. S. Khurmi and J. K. Gupta "A Textbook of Machine Design", Eurasia Publishing House (PVT.), Limited, Ram Nagar, New Delhi, 2005.
[14] Villa-Covarrubias, B., Baro-Tijerina, M., & Pina-Monarrez, M. R. Probabilistic Methodology to Determine the Shaft's Diameter and Designed Reliability.