Integral Analysis of Conveyor Pulley Using ®nite Element Method
Integral Analysis of Conveyor Pulley Using ®nite Element Method
Abstract
Pulleys are functionally critical items in belt conveyors. Earlier researchers have mostly analyzed the pulley in
parts rather than as an integral structure. Since the stresses and de¯ections of its various parts are inter-related, an
integral analysis is mandatory for reliable pulley design. A ®nite element software has been developed for the
analysis of pulleys utilizing a 3-noded thick conical shell element and a triangular ring element. This work considers
all the well-known assuctmptions on load variations along the pulley face-width and around its periphery. A semi-
analytical approach using Fourier analysis has been adopted. To validate the software, comparison has been carried
out with published experimental and theoretical investigations. A few parametric studies are also reported to
elucidate the mutual inter-dependence of stresses in the pulley components. # 1999 Elsevier Science Ltd. All rights
reserved.
0045-7949/99/$ - see front matter # 1999 Elsevier Science Ltd. All rights reserved.
PII: S 0 0 4 5 - 7 9 4 9 ( 9 8 ) 0 0 1 4 5 - X
304 M. Ravikumar, A. Chattopadhyay / Computers and Structures 71 (1999) 303±332
2. Loading details the belt wrap around its rim. Whereas the driven
pulley, having equal belt tension on either side
In a belt±pulley drive (with no-slip, and friction [Fig. 2(b)], can be considered to be subjected to only
coecient = m) the belt tensions (Fig. 2) are given by radial load.
the Euler's equation:
T1 T2 emyr :
1 2.1. Loading around the circumference
The varying belt tension around a drive pulley No well established experimental or real-life data is
[Fig. 2(a)] results in radial and tangential loads over available on the actual load variation around the
(drive) pulley circumference. From the Eulerian con- Both ®rst order (m = 1) and third order (m = 3) har-
cept, the variation is considered to be exponential; for monics are considered in the present work. Besides
practical purposes it is often assumed to be linear. these sinusoidal variations, a uniform load distribution
Both the types of variations are studied here. For any along the face-width has also been considered as fol-
angle, y (y1EyE y2), the tension is given by: lows:
exponential T
y T2 e m
y2 ÿy ;
2
T
z T
y=b:
7
T1 ÿ T2
linear T
y T2
y2 ÿ y;
3
y2 ÿ y1
for 0E y < y1 and y2 < yE2p, T(y) = 0.
The corresponding radial force per unit area on the
rim surface is obtained as 3. Finite element formulation
3.1. Element details the n harmonics) can be expressed as the sums of sym-
metric and anti-symmetric components:
Two types of axisymmetric elements have been
X
3 X
3
adopted for this work: u0 Ni u 0i cos ny Ni u0i sin ny;
i1 i1
1. conical shell element, since it is computationally e- X
3 X
3
cient; and u0 Ni u 0i cos ny Ni u 0i sin ny;
2. triangular ring element, as it can approximate com- i1 i1
plicated shapes. X
3 X
3
v0 Ni v0i sin ny ÿ Ni v0i cos ny;
3.1.1. Conical shell element i1 i1
Fig. 4. Details of elements with coordinate systems: (a) view from right end; (b) 3-noded thick conical shell; (c) triangular ring.
M. Ravikumar, A. Chattopadhyay / Computers and Structures 71 (1999) 303±332 307
x2 ÿ x x2 x
N1 ; N2 1 ÿ x2 ; N3 ; x s=l:
10
2 2
Substituting the above symmetric and anti-symmetric displacement components in the general strain±displacement
relations, the following matrices are obtained for the corresponding cases.
8 9
> X 3
> @Ni @Ni 0 >
>
>
> u0i z ui cos ny >
>
>
> @s @s >
>
>
> i1 >
>
>
> >
>
>
> X3 >
>
>
> 1 n u w n u0
w0
z2
w 00 >
>
>
> N v
0i
sin f
0i
z N v 0
i
sin f i
i
cos ny >
>
>
> 1 z r
i 0i
r l r
i i
r l 2 l >
>
>
> l i1 >
>
8 9 >
> >
>
> e s > >
> X3 >
>
>
> > > >
>
>
> e y >>
>
>
>
> 0 00
wi zwi cos ny
>
>
>
< e > = >
< i1
>
=
z
;
11
> g sy >
> > > >X 3 >
>
>
> >
> >
> @Ni @Ni 0 1 n v0i n 0 v0i >
>
>
> g > > v0i z v ÿ u0i ÿ sin f z ÿ ui ÿ sin f sin ny >
: sz > ; > >
> @s @s i 1 lz r r r r >
>
>
g yz >
> i1 >
>
>
> >
>
>
> X 3
@N @N z2
@N >
>
>
> i i 0 i 00 0 >
>
>
> w0i z wi wi ui cos ny >
>
>
> @s @s 2 @s >
>
>
> i1 >
>
>
> >
>
>
> X3 >
>
>
> 1 n n 0 2 n 00 v0i 0 >
>
>
: z
ÿ N w
i 0i ÿ z N w
i i ÿ z N w
i i ÿ v i sin ny >
;
1 l i1 r r 2r l
8 9
> X 3
>
>
> @Ni @Ni 0 >
>
>
> u0i z ui sin ny >
>
>
> @s @s >
>
>
> i1 >
>
>
> >
>
>
> 1 X3
n u w n u 0
w 0
z 2
w 00 >
>
>
> 0i 0i >
>
> ÿ Ni v0i sin f z ÿ Ni vi sin f 0 i i
i
sin ny >
>
>
>
> 1 z r r l r r l 2 l >
>
8 9 > >
>
l i1 >
>
>
> e > >
> X3 >
>
>
>
s >
> >
> >
>
>
> e y >
> >
> w 0
zw 00
sin ny >
>
>
< = < > > i i >
=
ez i1
:
12
>
> gsy >> >
> X 3
@N @N 1 n v n v 0 >
>
>
> >
> > >
> i i 0 0i 0 i >
>
>
> g > v0i z v u0i ÿ sin f z u ÿ sin f cos ny >
: sz > ; > >
> @s @s i 1 lz r r r i r >
>
>
gyz >
> i1 >
>
>
> >
>
>
> X3
@N @N z 2
@N >
>
>
> i i 0 i 00 0 >
>
>
> w 0i z w i w i u i sin ny >
>
>
> @s @s 2 @s >
>
>
> i1 >
>
>
> >
>
>
> X3 >
>
>
> 1 n n 0 2 n 00 v0i 0 >
>
>
: z
N i w 0i z N i wi z N i w i ÿ vi cos ny >
;
1 l i1 r r 2r l
Assuming linear displacement models for shape function and using isoparametric form, the strain±displacement
matrices for the symmetric and anti-symmetric cases are:
308 M. Ravikumar, A. Chattopadhyay / Computers and Structures 71 (1999) 303±332
8 9
> @u >
>
> cos ny >
> has been undertaken to select the number of Fourier
>
> @z >
>
>
> >
> harmonics n.
>
> w v >
>
>
> n cos ny >
>
>
> >
>
8 9 >
> r r >
>
>
e >>
> >
> 4.1. Rim
>
>
z >
>>
> @w >
>
>
>
e y >
>>
> cos ny >
>
>
< = >>
< >
=
e r @r
;
14 Fig. 5 describes the axial and circumferential stresses
>
> g zy >> >
> u @v >
> around the rim periphery as evaluated from the present
>
> > > ÿ n sin ny >
>
> g > > > > r @z >
> software. Here, all the above mentioned types of load
: zr > ; > >
>
>
>
>
g yr >
> @w @u >
> distributions along the face-width have been considered.
>
> >
>
>
> cos ny >
> Lange's experimental values are also shown therein. It
>
> @z @r >
>
>
> >
>
>
> >
>
>
> w @v v >
: ÿ n ÿ sin ny > ;
r @r r
8 9
> @u >
>
> sin ny >
>
>
> @z >
>
>
> >
>
>
> w v >
>
>
> ÿ n sin ny >
>
>
> >
>
8 9 >
> r r >
>
> e z >
>
> >
>
>
> >
> >
> @w >
>
>
> e >
> >
> sin ny >
>
>
< = y> >
< >
=
er @r
:
15
>
> gzy >
> >
> u @v >
>
> >
> >
> n cos ny >
>
>
>
> g > > r @z >
: zr >; > >
>
>
>
>
gyr >
> @w @u >
>
>
> >
>
>
> sin ny >
>
>
> @z @r >
>
>
> >
>
>
> >
>
>
> w @v v >
: n ÿ cos ny > ;
r @r r
4. Validatory studies
may be observed that the peak stresses calculated pulley with exponential and linear increase in tension
assuming uniform distribution are quite close to the ex- are furnished in Fig. 9. Both the cases provide almost
perimental values as against those obtained considering the same results throughout.
®rst and third order sinusoidal variations. This is A comparison of the stresses obtained by considering
because the tensioning arrangement in the experimental the conical shell and triangular ring elements is reported
set-up provided uniform load along the pulley face. in Fig. 10. Both match quite closely. But, it is noted
Fig. 6 compares Lange's experimental plots of stresses that computation with a conical shell element is about 3
along the face-width at dierent peripheral angles with times faster than with the triangular element and it
the prediction of this software. While the agreement is requires about 2 times less memory.
reasonable in most of the cases, for y = 958, the trend Fig. 11 provides a comparison of the stresses around
of experimental plot diers marginally. However, in the the rim periphery calculated by assuming the rim as
absence of data at all points, the experimental graph simply supported at the rim±web junction, and consider-
appears to be approximate. ing the rim as an integral part of the entire pulley.
The stress variations around the rim at the mid-sec- Lange's experimental values are also superimposed
tion due to various wrap angles (60, 120, 180, and 2408) thereon. It is interesting to observe that the results from
are furnished in Fig. 7. In line with the observations analysis with simply supported consideration is in very
reported earlier [7, 13], the stress magnitude tends to go good agreement with those from the integral analysis
up with reduced wrap angle due to localization of the and experiment.
applied load. The location of maximum stress also
depends signi®cantly on the wrap angle. 4.2. Web
Fig. 8 compares the rim stress variation along the
face-width as calculated for a drive pulley reported in The web stresses are closely dependent on the overall
Ref. [7]. The stresses around the periphery of the drive rigidity of the pulley structure [3]. Hence, unlike in the
Fig. 7. Rim stresses around the circumference at mid-section for various wrap angles. Yr = 608, - - - - - Yr = 1208,
± ± ± ± Yr = 1808,Ð Ð Yr = 2408.
M. Ravikumar, A. Chattopadhyay / Computers and Structures 71 (1999) 303±332 311
rim, the stresses in the web need to be analyzed consid- not reported in Lange's work, and hence are not com-
ering the entire pulley as an integral structure. For a pared. Fig. 14 shows the de¯ections at those points as
complex shape such as an entire pulley, the shell element compared with measured values.
can not be used. Therefore, the triangular ring element
has been adopted for the integral analysis. 4.3. In¯uence of aspect ratio
Fig. 12 indicates that the radial and circumferential
stresses in the web compare well with Lange's exper- For a reliable ®nite element analysis, usually the
imental data. The radial stresses at various radii on the aspect ratios should be within 1±3. Since the pulley is
web (presented in Fig. 13) also agree closely with the ex- constituted of parts mostly having thin cross sections, it
periments. Circumferential stresses at these points are is very convenient to have meshes with higher aspect
ratios. This reduces the total number of elements, mini-
mizing the required computational memory and time. A
comparative study with 3 dierent ratios, 1.25, 2 and
2.5 was made to ascertain a suitable value. Veri®cation
with experimental investigations on web stresses (Fig. 15)
shows that analysis with an aspect ratio of around 1.25
produces good results. Hence, aspect ratios of 1±1.5
have generally been adopted in the areas of interest.
Fig. 16 shows a typical mesh pattern used for estimating
the stresses at the taper-lock interfaces with shaft and
hub.
5. Design studies
Fig. 9. Rim stresses around the circumference at mid-section for linear and exponential variation in belt tensions. Linear,
Y Y Y Exponential.
M. Ravikumar, A. Chattopadhyay / Computers and Structures 71 (1999) 303±332 313
Fig. 10. Rim stresses around the circumference at mid-section Fig. 11. Rim stresses around the circumference at mid-section
with dierent elements; simply supported rim. . . . for dierent supporting conditions. . . . Experiments,
Experiments, - - - - Triangular ring element, Conical Y Y Y Rim simply supported at webs, Integral pulley.
shell element.
A comparison of the stresses in the rim mid-section
project beyond the rim. MHEA shows that when the around the periphery due to the three web locations is
distance between the webs equals one-half the bearing shown in Fig. 18. Similarly, Fig. 19 depicts the maxi-
centers, the web stress is highest. Since the location of mum stresses in the web. The maximum radial stresses
the web in¯uences the stresses in other parts of the in the web obtained from this study were found to com-
pulley as well, a rational approach is needed to ®nd out pare well with the magnitude calculated from the
the web location, minimizing the stresses in all related MHEA formula for the assumed loading on this model
parts. Hence, the study has been undertaken with the pulley. It is clear that the ®rst location causes the mini-
following three locations of the web, shown in Fig. 17. mum stresses. It is important to observe from these
plots that, though the web stresses at the outer radii are
1. Web placed along the ideal belt edge; very low, they increase considerably near the hub±web
2. web shifted outside, near to the pulley edge; and junction. As reported elsewhere [3, 4, 14], most pulley
3. web shifted inside, away from the pulley edge. failures are initiated from this region only.
314 M. Ravikumar, A. Chattopadhyay / Computers and Structures 71 (1999) 303±332
5.3. Shaft
Fig. 15. In¯uence of aspect ratio on web stresses. . . . Experiments, - - - - Aspect ratio = 1.25, Aspect ratio = 2, ±±±Aspect
ratio = 2.5.
318 M. Ravikumar, A. Chattopadhyay / Computers and Structures 71 (1999) 303±332
Fig. 16. Typical mesh considered for analyzing stresses at taper-lock interfaces.
M. Ravikumar, A. Chattopadhyay / Computers and Structures 71 (1999) 303±332 319
Fig. 18. Variation in rim stresses at mid-section around the circumference for dierent web locations. Web locationÐ1,
- - - Web locationÐ2, ±±±Web locationÐ3.
M. Ravikumar, A. Chattopadhyay / Computers and Structures 71 (1999) 303±332 321
Fig. 19. Variation in maximum web stresses for dierent web locations. - - - Web locationÐ1, Web locationÐ2, ±±±Web
locationÐ3.
322 M. Ravikumar, A. Chattopadhyay / Computers and Structures 71 (1999) 303±332
Fig. 20. Variation in stresses at rim±web junction for dierent web locations. Web locationÐ1, - - - Web locationÐ2,
±±±Web locationÐ3.
M. Ravikumar, A. Chattopadhyay / Computers and Structures 71 (1999) 303±332 323
Fig. 21. Variation in stresses at plain-hub/shaft interface for dierent web locations, at y = 1808. - - - Web locationÐ1, ±±±Web
locationÐ2, Web location 3.
324 M. Ravikumar, A. Chattopadhyay / Computers and Structures 71 (1999) 303±332
Fig. 22. Variation in stresses at hub/taper-lock interface for dierent web locations, at y = 1808. - - - Web locationÐ1, Web
locationÐ2, ±±±Web locationÐ3.
M. Ravikumar, A. Chattopadhyay / Computers and Structures 71 (1999) 303±332 325
Fig. 23. Variation in stresses at shaft/taper-lock interface for dierent web locations, at y = 1808. - - - Web locationÐ1, Web
locationÐ2, ±±±Web locationÐ3.
326 M. Ravikumar, A. Chattopadhyay / Computers and Structures 71 (1999) 303±332
Fig. 24. Variation in rim stresses around the circumference at mid-section due to dierent stieners. Stiener caseÐ1,
- - - Stiener caseÐ2, ±±±Without stiener.
M. Ravikumar, A. Chattopadhyay / Computers and Structures 71 (1999) 303±332 327
Fig. 25. Variation in web stresses due to dierent stieners. With stiener, - - - Without stiener.
328 M. Ravikumar, A. Chattopadhyay / Computers and Structures 71 (1999) 303±332
Fig. 26. Variation in web stresses due to dierent shaft diameters. 90 mm F, - - - 100 mm F, ±±±135 mm F.
M. Ravikumar, A. Chattopadhyay / Computers and Structures 71 (1999) 303±332 329
Fig. 27. Typical tapered webs; (a) type A; (b) type B; (c) type C.
330 M. Ravikumar, A. Chattopadhyay / Computers and Structures 71 (1999) 303±332
References
[11] Vodstrcil R. Analysis of a belt conveyor pulley using ®nite [14] Ravikumar M. Studies on conveyor pulley design using
element methods. In: 4th International Conference in ®nite element method. M.S. thesis. IIT Madras, India,
Australia on Finite Element Methods. Melbourne, 1995.
Australia, 1982:143±148. [15] Zienkiewicz O.C. The ®nite element method in engineering
[12] Qiu Xiangjun, Sethi V. A new pulley stress analysis science. London: McGraw Hill, 1971.
method based on modi®ed transfer matrix. Bulk Solids [16] Sivadas KR, Ganesan N. Vibration analysis of thick com-
Handling 1993;13:713±24. posite clamped conical shells of varying thicknesses. J
[13] Siva Prasad N, Sarma Radha. A ®nite element analysis for Sound Vibr 1992;152:27±37.
the design of a conveyor pulley shell. Comp Struct [17] BIKON-Technik GmbH Germany. Welle-Nabe-
1990;35:267±77. Verbindungen.