P. G. Student, 2,3Faculty members

Department of Mechanical Engineering, Jadavpur University, Kolkata – 700032, India [email protected], [email protected], [email protected]

Abstract In the present work, the system modelled is a rotating disk of exponentially varying thickness, containing attached masses distributed radially at various locations. The limit angular speed for given geometry is obtained for various loading conditions. The problem is formulated through a variational method, where the radial displacement field is taken as unknown variable. Assuming a series solution and using Galerkin’s error minimization principle, the solution of the governing partial differential equation is obtained. Corresponding to various rotational speeds, the stress distribution in the disk is estimated and the limit angular speed is determined using Von Mises yield criterion. The effect of geometry and loading parameters are studied. The values of optimum geometry parameters for given loading conditions are reported. It is also observed that radial shifting of attached mass location changes the location of yield initiation. The relevant results are reported in dimensionless form, and hence one can readily obtain the corresponding dimensional value for any geometry through appropriate normalizing parameters. It is observed that the present methodology is quite robust, stable and realistic and the documented results may be used by the practicing engineers as design monograms. Keywords: Limit angular speed, Von Mises criterion, Annular Disk.

Nomenclature a, b

Inner and outer radii of the disk.

{c }

The vector of unknown coefficients.


Thickness of the disk at radius r.


Root thickness of the annular disk.


Thickness of equivalent uniform annular disk, β b.


Mass of the rotating disk.

n, k

Geometry parameters controlling the thickness variation of disk.


Radii of the attached mass.


Displacement field of the disk.


Ratio of attached mass to disk mass.


Thickness ratio, hc / b.


Aspect ratio, a /b.

Parameter ( b - a ).

µ, E

Poisson’s ratio and elasticity modulus.


Normalized coordinate in radial direction, (r - a ) / ∆ .


Density of the disk material.


Yield stress of the disk material.


Radial stress.


Tangential stress


The set of orthogonal polynomials used as coordinate functions.


Angular velocity of the disk.


Non-dimensional radii of the attached mass location.

Dimensionless angular velocity, ω b ρ / σ 0 .

1. Introduction High speed rotating disks find widespread application in mechanical engineering as a precision component, some common examples being high speed impellers, gears and flywheels. The stress distribution in a rotating disk depends on the rotational speed. At a particular speed the developed stresses exceed the yield stress of the material. This particular speed is termed as limit angular speed of the rotating disk. It is well understood that stresses in a rotating disk are maximum at its root. It is this theory that initiated the idea of radially varying the thickness of a disk without hampering its axisymmetric nature. Depending on the need and method of mounting, the rotating disks can be made solid or annular. One specific advantage of annular disks arises from the fact that, since usually disks are manufactured using forging processes, very large forgings are likely to have defects in the material at the center. Hence, to eliminate uncertainties it is a common practice to bore a central hole along the axis of the rotor [1]. Analysis of elastic stress state in rotating disks having variable thickness up to the point of yielding with the objective to determine limit angular speed has been extensively carried out by many researchers. Analysis based on dynamic relaxation technique [2] has been applied effectively to study elastic stress behaviour for isotropic material possessing non-linear stress-strain relationship. Exponential variation of thickness was first discussed in [3] assuming Tresca’s criterion and its associated flow rule along with linear strain hardening material behaviour, although solution satisfying all boundary and continuity conditions was not obtained. However in [4], a similar problem was studied and an analytical solution was presented that satisfy all boundary and continuity conditions. The analytical method was further applied in [5] to obtain solution for elastic deformation of solid and annular disk having parabolically varying thickness subjected to different boundary conditions. In

another paper, application of variational principle [6] involving Galerkin’s principle and assuming a series solution has provided an efficient numerical scheme for the estimation of limit angular speed of disks with and without attached masses. The analysis is based on Von Mises stress criterion and has been extended for design and optimization of various types of disk geometries [7, 8]. The application of variational principle provides an advantage over other methods in terms of simplicity and ease with which various complicating effects such as addition of attached masses at various radial locations; rigid inclusions, etc. are incorporated. The present work employs variational principle to study the stress distribution and estimate the limit angular speed of rotating annular disks with exponentially varying thickness. The variables related to the geometry of the disk are identified as geometry parameters and those controlling the external loading of the disk are termed as loading parameters. The effect of geometry and loading parameters is studied individually with an objective of maximizing the limit angular speed.

2. Mathematical Formulation An annular disk of exponentially varying thickness subjected to centrifugal loading is shown in fig.1. The mathematical model is framed based on the assumptions that the material of disk is homogeneous, isotropic and linear-elastic.

Fig. 1. An annular disk having exponentially varying thickness The thickness variation of the disk in fig.1 is given by h (ξ ) = h 0 e xp [ − n (ξ ) k ] where the parameters n and k control the disk geometry. In each case of thickness variation, weight of the disk is kept fixed by adjusting the value of h0, i.e. the disk thickness at r = a. The purpose of this study is to maximize the operating speed range of a disk having mass equivalent to that of a uniform disk by varying its thickness. For this purpose, a parameter β is introduced which defines the mass of the equivalent uniform disk having same inner and outer radii. As the analysis is carried out under plane stress assumption, thickness ratio, β is assumed to be 0.10 in the present study. The variation of profile with geometry parameters k and n is shown in fig. 2 (a, b) for disks with constant root thickness. With increase in rotational speed the stresses increase due to enhanced effect of centrifugal force, the magnitude of which is governed by the boundary conditions of the disk. The

solution for the displacement field is obtained from the minimum potential energy principle δ(π) =

δ(U+V) = 0, where, U and V are strain energy stored in the disk and potential energy arising out of centrifugal force, respectively. Substituting the expressions of U and V in the energy principle δ(π) = 0, the governing equilibrium equation is obtained [1]. ∆ π E δ 1 − µ 2 

2 1  u 2  2 µ u  d u  ( ∆ ξ + a)  d u   2 2   +   hd ξ 2 π ρ ω ( ξ a ) uhd ξ − ∆ ∆ + +  =0   ( ∆ ξ + a) ∆ d ξ  ∆ 2  d ξ    0  0   1

The displacement function u(ξ) in equation (1), is approximated by u (ξ ) ≅


∑ c φ , i=1, 2,…, n, i i

where φi is the set of orthogonal functions developed through Gram-Schmidt scheme. The necessary starting function to generate the higher order orthogonal functions is selected by satisfying the relevant boundary conditions;


(a )

= 0 and σ r

( b)

= 0.

Fig. 2 (a). Profile variation keeping n fixed at 1.5

Fig. 2 (b). Profile variation keeping k fixed at 0.5

while k varies from 0.0-1.0 (λ = 0.20)

while n varies from 0.0-2.0 (λ = 0.20)

The effect of attached mass is taken into consideration through a factor α, which indicates the ratio of the attached mass to the mass of the total disk m. This attached mass will contribute to the centrifugal force by an additional amount, α m ω rm , where rm indicates the radius of the location of the attached 2

mass. The contribution of this term is included in the expression of equation (1) and substituting the assumed series solution of u (ξ), followed by Galerkin’s error minimization principle, the following equation is obtained.

E 1− µ




∑∑ i =1

= ρω


j =1 2


φ φ

∑ ∫ {( ∆ ξ n

i =1


i j + ( φ i′φ  0 ξ a) ( ∆ + ∆  1


+ a )2 φ i



+ φ i φ ′j ) +

(∆ξ + a ) φ i′φ ∆2

α m  2 hd ξ +  ω ξ m φ i  2π ∆ 


 ′j  hd ξ 


The above equations are the set of system governing equations which are expressed in matrix form as [K] {c} = {R}. Finally the unknown coefficients are obtained from {c} = [K]-1 {R}.

3. Results and discussion The numerical values of the different system parameters, considered are, E = 207 GPa, µ = 0.3, ρ = 7850 Kg/m3 and σ 0 = 350 MPa. The dimensionless angular speed corresponding to the onset of yielding is denoted by Ω1. Various parameters considered in the present study are classified into two categories. Thickness control parameters k and n, and aspect ratio λ are considered as geometry parameters and loading effect of attached mass α and its location ξ


are considered as loading

parameters. Extensive numerical simulation is carried out to observe the effect of different system parameters on Ω1. In each case, the analysis is carried out for disks of identical mass. Von Mises yield criterion is considered here to establish the onset of yielding, in which the yield condition is given by σ r




− σ rσθ + σθ ≥ σ 0 .

3.1 Disks without attached masses The effect of geometry parameters on limit angular speed at two different aspect ratios is investigated and plotted in fig.3-4. It is observed that for given λ, with increase in n, Ω1 also increases; however for each value of n, there exists an optimum value of k (kopt) that maximizes Ω1. The variation of Von Mises stress for different profiles and the representation of the profiles are provided in fig. 5 (a, b). The figure indicates that for each case maximum stresses occur at the root of the disk. For values of k less than kopt, root thickness changes abruptly giving rise to stress concentration effects, while for larger values of k, root thickness further reduces thus lowering the strength at root. This probably explains the existence of kopt for each value of n.

Fig. 3. Contour and 3-D plot of variation of Ω1 with geometry parameters (λ = 0.20) A similar observation is made in fig. 6 (a, b) for a different value of λ. It is also observed that increase in λ is accompanied by reduction in Ω1 and higher values of kopt.

Fig.4. Contour and 3-D plot of variation of Ω1 with geometry parameters (λ = 0.40)

Fig.5 (a). Radial distribution of Von-Mises stress

Fig. 5 (b). Disk profile for different k values

Fig. 6 (a). Radial distribution of Von-Mises stress

Fig. 6 (b). Disk profile for different k values

3.2 Disks with attached masses In the present study, attached mass ratio and its location is varied from 0.00 to 1.00 and 0.25 to 1.00 respectively. Detailed investigation reveals reduced values of Ω1 as well as kopt for increasing attached mass ratio at a given location. In fig. 7 (a. b), the effect of parameter k on limit angular speed is studied under different loading conditions for two different aspect ratios. In these figures, n is kept fixed at 1.0 but its effect for a given value of aspect ratio and loading parameters is presented later. It is evident from the figures that increase in λ is accompanied by fall in limit angular speed as well as kopt.

Fig. 7. Effect of k on Ω1 under different loading parameters keeping n constant at 1.0

Fig.8. Radial distribution of Von-Mises stress at different aspect ratios (n = 1.0, k = 0.5) It is further observed that for a given attached mass ratio, an increase in the value of ξ


also reduces

Ω1 and kopt (fig. 7 a-b). It is due to the fact that attached masses add to the centrifugal effect thereby increasing the stresses leading to reduced values of Ω1. Reduced kopt value results in larger root thickness, thereby increasing strength at the root which takes up the additional centrifugal effect. The Von-Mises stress distribution at given values of n and k for different aspect ratios under various

combination of loading parameters is plotted in fig. 8 (a, b). It is observed that the maximum stresses occur at the root of the disk irrespective of the loading conditions. The effects of geometry parameters on limit angular speed at λ = 0.20 for given values of attached mass ratio and location is studied in detail and some sample plots are shown in figs. 9-10. A similar study is carried out for λ = 0.40 and the results are plotted in figs. 11-12. These figures indicate once again, that increase in λ and α cause a reduction in the value of limit angular speed and kopt. Thus a detail study on the interactions of λ and α may yield come interesting findings towards an optimized disk design. In figs. 13-14, the effect of loading parameters on limit angular speed is plotted for a given values of geometry parameters and aspect ratios. These variations are almost identical to that of the case of a solid disk, as reported in an earlier study [8].

Fig. 9. Contour and 3-D plot of variation of Ω1 with geometry parameters (α = 0.50, ξ


0.75, λ = 0.20)

Fig. 10. Contour and 3-D plot of variation of Ω1 with geometry parameters (α = 1.00, ξ


0.75, λ = 0.20)

4. Conclusion Modeling and analysis of exponentially varying high speed annular disks is carried out using variational formulation and assuming a series solution of unknown displacement field. Assuming Von -

Fig. 11. Contour and 3-D plot of variation of Ω1 with geometry parameters (α = 0.50, ξ


0.75, λ = 0.40)

Fig. 12. Contour and 3-D plot of variation of Ω1 with geometry parameters (α = 1.00, ξ


0.75, λ = 0.40)

Fig. 13. Contour and 3-D plot of variation of Ω1 with loading parameters (n = 1.50, k = 0.50, λ = 0.20)

Fig. 14. Contour and 3-D plot of variation of Ω1 with loading parameters (n = 1.50, k = 0.50, λ = 0.40) Mises failure criterion the effect of various geometry and loading parameters on limit angular speed is investigated and corresponding parameter values leading to the optimization of limit speed are reported. The results obtained are in dimensionless form and presented graphically. The results provide a substantial insight in understanding the behavior of rotating disks with and without attached masses. The method of formulation gives a kernel for subsequent dynamic analysis and study of many other complicating effects. Further the method developed has application potential in various other problems, e.g., shrink fitted disks, disks made of anisotropic material, etc.

References 1. Timoshenko S., Strength of materials, advanced theory and problems Part-2, CBS, 2002. 2. Sherbourne A.N. and Murthy D.N.S., Stresses in discs with variable profile, Int. J. Mech. Sci., Vol. 16, pp. 449-459, 1974.

3. Güven U., On the applicability of Tresca’s yield condition to the linear hardening rotating solid disk of variable thickness, ZAMM, Vol. 75, pp. 397–398, 1995.

4. Eraslan A.N. and Orcan Y., Elastic–plastic deformation of a rotating solid disk of exponentially varying thickness, Mech. Mat., Vol. 34, pp. 423–432, 2002.

5. Apatay T. and Eraslan A.N., Elastic deformations of rotating parabolic disks: an analytical solution, Fac. Eng. Arch. Gazi Univ., Vol. 18, pp. 115-135, 2003.

6. Bhowmick S., Pohit G., Misra D. and Saha K.N., Design of high speed impellers, Proc. Int. Conf. Hyd. Engg. Res. Prac., IIT, Roorkee, pp. 229-241, 2004.

7. Bhowmick S., Misra D. and Saha K.N., On optimum disk profile for high speed impellers having parabolic thickness variation, Proc. 31st Nat. Conf. Fluid Mech. and Fluid Power, Jadavpur University, Kolkata, pp. 417-424, 2004.

8. Bhowmick S., Misra D. and Saha K.N., Modeling of high speed impellers under various loading parameters and optimization for exponentially varying disk profile, Proc. Nat. Conf. Front. App. Comp. Math., TIET, Patiala, pp. 449-462, 2005.

a study on the effect of geometry and loading ...

In the present work, the system modelled is a rotating disk of exponentially ..... Mises failure criterion the effect of various geometry and loading parameters on ...

505KB Sizes 1 Downloads 180 Views

Recommend Documents

Study and Investigate Effect of Input Parameters on Temperature and ...
Equipment Components AB SE–631 85 Eskilstuna, Sweden 2007. [12]. www.fuchs-europe.de. [13]. Industrial Gearbox Service Manual of LOCTITE. [14]. Industry Catalogue of HI-TECH DRIVES Pvt. Ltd. [15].Carl Byington, Ryan Brawrs, Sanket Amin, James Hopki

Study and Investigate Effect of Input Parameters on ... - IJRIT
to apply DOE techniques to achieve desired design of gearbox for control the temperature and noise .... replication total 32 experiments will be performed as shown in table ΙΙ. ... will be carried out using dB meter or by using ultrasonic sensor.

Effect Of Ecological Factors On The Growth And Chlorophyll A ...
Effect Of Ecological Factors On The Growth And Chlor ... ed Kappaphycus alvarezii In Coral Reef Ecosystem.pdf. Effect Of Ecological Factors On The Growth And ...

The effect of culture on perception and cognition - A conceptual ...
self-perceptions, with English-speaking bi-cultural people reporting a. perception of the self as independent of others and Chinese-speaking. bi-cultural people ...

A photoacoustic study of the effect of doping ...
Cochin University of Science and. Technology ... Abstract. We report a photoacoustic (PA) study of the thermal and ... The ex- perimental data obtained from the measurement of the PA signal as a ..... working towards his PhD degree at the In-.

effect of crop geometry, intercropping system and ...
MANAGEMENT ON COB YIELD AND NUTRIENT UPTAKE OF BABY CORN .... Crop geometry levels significantly altered the .... Effect of different levels of.

(d) via the gradient flow of the Frobenius condition number ˜µ(f,ζ) defined on this variety. We recall and introduce a few definitions before we state our principal ...

On the geometry of a generalized cross-correlation ...
defined for analyzing functional connectivity in brain images, where the main idea ... by a postdoctoral fellowship from the ISM-CRM, Montreal, Quebec, Canada.

On the geometry of a generalized cross-correlation ...
Random fields of multivariate test statistics, with an application to shape analysis. Annals of Statistics 36, 1—27. van Laarhoven, P., Kalker, T., 1988. On the computational of Lauricella functions of the fourth kind. Journal of. Computational and

An Experimental Study on the Capture Effect in ... - Semantic Scholar
A recent measurement work on the cap- ture effect in 802.11 networks [12] argues that the stronger frame can be successfully decoded only in two cases: (1) The.

Doing and Learning: The Effect of One on the Other
Technion, Israel Institute of Technology, Department of Education in Technology and Science. Introduction. During the last ... constructing, and debugging some amazing machinery, their social and communication skills improved. Another .... o A change

Doing and Learning: The Effect of One on the Other
design, to the amazement of his care-takers. Thinking and. Multi-Media. Doppelt, MA. Thesis, in press. 10 th graders, judged and classified by school system to low- level classes. Tagged as un-able to pass matriculation. Considered at school as troub

The effect of coherence and noise on the ...
LFMs, is shown to create large side lobes in the time domain. Alternative ..... free numerical simulations produce very similar focal patterns as shown in Fig.

The Effect of a Prerecorded Conversation on ...
speakers on a laptop computer at the maximum volume level. Participants used ... effect of auditory distraction on reading comprehension, F(1,9). 0.630, p 0.448.

the study of virechana karma and it‟s effect on body ... - PDFKUL.COM
the interstitial fluid through the pores of the capillary membranes. Important Constituents of Intracellular Fluid: The ICF is separated from the extra cellular fluid by a ...... 4. Complete evaluation of loose stools produced during. Virechana karma

The Effect of Crossflow on Vortex Rings
The trailing column enhances the entrainment significantly because of the high pressure gradient created by deformation of the column upon interacting with crossflow. It is shown that the crossflow reduces the stroke ratio beyond which the trailing c

The effect of Medicaid on Children's Health: A ...
design. I exploit the discontinuity generated by Medicaid's eligibility rule, based on family income, on program participation rates. In contrast with a standard regression discontinuity approach, here there are ... reau's March 2009 and 2010 Current

The Effect of Medicaid Premiums on Enrollment: A ...
However, this result is not robust to other choices of bandwidth. These patterns in the data are consistent with the conclusion that treatment status is unrelated to sample composition and supports identification. As a final check, I perform a set of

Hydrodynamic effect on the behaviour of a corrosion ...
Available online 27 December 2005. Abstract. Influence of hydrodynamic ... +33 4 72 43 81 87; fax: +33 4 72 43 87 15. E-mail address: [email protected] ..... trolyte conductivity evolution due to the greater account of iron.

The Effect of Crossflow on Vortex Rings
University of Minnesota, Minneapolis, MN, 55414, USA. DNS is performed to study passive scalar mixing in vortex rings in the presence, and ... crossflow x y z wall. Square wave excitation. Figure 1. A Schematic of the problem along with the time hist


The effect of production system and age on ...
(P < 0.05). Aspects of the fatty-acid patterns that are of relevance to human nutrition tended to favour the .... Data analysis employed a block design within the.

Selection and the effect of smoking on mortality
death T, and smoking S, controlling for observed individual characteristics X and .... Within this large data set, Statistic Sweden has constructed a smaller panel.

Effect of voriconazole on the pharmacokinetics and ...
of voriconazole (Vfend tablet; Pfizer, New York, NY) .... SD, except for tmax data, which are given as median and range. ..... Measurements of recovery from.