Journal of MaterialsVolume 2013 (2013), Article ID 809205, 9 pageshttp://dx.doi.org/10.1155/2013/809205

Research Article

Finite Difference Solution of Elastic-Plastic Thin Rotating Annular Disk with Exponentially Variable Thickness and Exponentially Variable Density

Sanjeev Sharma and Yadav Sanehlata

Department of Mathematics, Jaypee Institute of Information Technology, A-10, Sector 62, Noida-201307, Uttar Pradesh, India

Received 10 December 2012; Accepted 14 February 2013

Academic Editor: Francois Peeters

Copyright © 2013 Sanjeev Sharma and Yadav Sanehlata. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Elastic-plastic stresses, strains, and displacements have been obtained for a thin rotating annular disk with exponentially variable thickness and exponentially variable density with nonlinear strain hardening material by finite difference method using Von-Mises' yield criterion. Results have been computed numerically and depicted graphically. From the numerical results, it can be concluded that disk whose thickness decreases radially and density increases radially is on the safer side of design as compared to the disk with exponentially varying thickness and exponentially varying density as well as to flat disk.

1. Introduction

Due to wide applications of rotating disk, circular disk, spherical shells, cylinders, and shafts in engineering, elastic-plastic analysis of rotating disk, becoming more and more active topic in the field of solid mechanics. The research on them is always an important topic, and their benefits have been included in some books [1–3]. Sharma et al. [4] uses the concept of transition theory to evaluate the stresses for disk with different parameters like variable thickness and variable densities. The elastic-plastic strain hardening problems of annular disks with constant thickness under external pressure were firstly studied by Gamer [5] using linear stress-plastic strain relation. Later on, the work was extended to annular disk with some especially thickness functions by Güven [6]. The elastic-plastic linear strain hardening problems of rotating annular disks subjected to angular velocity are easier to deal with as compared to nonlinear strain hardening problems of annular disks. Therefore, main focus of research now these days is on rotating annular disk made of strain hardening material. In contrast, very few researchers are there who investigated elastic-plastic deformation and stresses for rotating disks with nonlinear strain hardening.

The obvious advantage using a linear strain hardening stress-plastic strain is that a closed form solution can be obtained for annular disks with constant thickness and some especially variable thickness functions. However, most of the materials exhibit nonlinear strain hardening behavior; thus, this nonlinearity is obvious in the transition region from elastic to plastic parts of stress-strain curve. Due to the previously reason, a polynomial stress-strain relation of nonlinear strain hardening material is proposed in the papers of You et al. [7–9]. Numerical method such as finite difference method is an effective technique to the stresses and strains for these rotating disks. However, for scientific research and engineering analysis, analytical methods and numerical methods are still very active. Therefore, Sterner et al. [10] proposed truncated Taylor’s series numerical method to solve elastic problems of rotating disk with arbitrary variable thickness. You and Zhang [8] examine elastic-plastic stresses in rotating disks, using Runge-Kutta’s method.

In this paper, we proposed a more straightforward and most effective numerical method such as finite difference method which is most celebrated method to solve boundary value problems. The proposed method is used to analyze the stresses, strains, and displacements for annular disk having exponential variable thickness and exponential variable density with nonlinear strain hardening material behavior.

2. Mathematical Formulation

Assuming that the stresses vary over the thickness of the disk, the theory of the disks of variable thickness can give good result as that of the disks of constant thickness as long as they meet the assumption of plane stress. For disk profile, it is assumed that disk is symmetric with respect to the mid plane. This profile is defined by the thickness function and density function where , , and are geometric parameters, is the radius of the disk, is the thickness at the axis of the disk, and is the density of the material.

2.1. Basic Equations

For rotating disk with variable thickness and variable density, the governing equilibrium equation is where is the radial coordinate, and are the radial and circumferential stresses, and is the angular velocity of the disk.

The strains and radial displacement are where is the displacement component in the radial direction, is the radial strain, and is the circumferential strain.

The equation of compatibility can be derived from (3) as follows: The relations between stresses and elastic strains for plane stress problems can be determined according to theory of elasticity where , are the elastic radial and circumferential strains, and and are Lame’s constants.

For plastic deformation, the relation between the stresses and plastic strains can be determined according to the deformation theory of plasticity [3] where and are the plastic radial and circumferential strains, is the equivalent plastic strain, and is the equivalent stress.

The Von-Mises yield criterion is given by .

The total strains are the sum of elastic and plastic strain Let us now introduce a stress function and assume that the relations between stresses and the stress function are Substituting (8) into (5), and further substitution into (7), one obtains where .

By the substitution of (9) into (1), we have where , .

By the nonlinear strain-hardening material model proposed by You and Zhang [8], the stress-strain relationship can be written as where is the equivalent total strain and the yield strain and , are material constants.

Substitution of the second equation of (11) into (6), the plastic strains can be written as

The governing equation in the plastic region of the rotating disks in terms of stresses and stress function can be obtained by substituting (12) into (10) as The values of the stress function at the elastic-plastic interface radius are the same, and, therefore, the stress function is continuous at the interface radius [9]. It can be seen from the continuity of the stress function and (3), (5)–(7), and (8) that the continuity conditions of the stresses and displacements at the elastic-plastic interface radius are satisfied.

The boundary conditions for the rotating annular disks are where and are the inner and outer radii of the rotating disks, respectively.

3. Finite Difference Algorithm

To determine the elastic-plastic stresses, strains, and displacement in thin rotating disks with a nonlinear strain hardening material, we have to solve the second-order nonlinear differential equation (13) under the given boundary condition (14). The general form of (13) can be written as

3.1. Finite Difference Algorithm Steps

The second-order differential equation (13) with the given boundary conditions can be solved by using finite difference method.(i) First, partition the domain into subintervals of length . (ii) To express the differential operators and in a discrete form, we use the finite difference approximations as (iii) With , we have nodal points . The values at the end points are given by the boundary conditions; that is, we are given that , . Using the finite difference approximation, we get the following system of equations: (iv) After simplifying and collecting coefficients of , , and , in Step 3, the boundary value problem results in system of nonlinear equations.(v) The solution of the system of nonlinear equations can be obtained using Newton-Raphson’s method where the unknowns are . By this process, we will get the stress function and then stresses, strains, and displacement.

4. Numerical Illustration and Discussion

A two-dimensional plane stress analysis of rotating disk with nonuniform thickness and nonuniform density are carried out using finite difference method. The radius of the rotating disk is taken to be , and its materials properties are material density Kg/m³, Lame’s constant are GPa, GPa, and Poisson’s ratio .

4.1. Rotating Annular Disk with Exponentially Variable Thickness and Exponentially Variable Density

It has been observed from Figures 1(a) and 1(b) that as thickness decreases and density increases radially, circumferential stress is going on decreasing. Also, circumferential stress is maximum at the internal surface. With the increase in angular velocity, circumferential stress is going on increasing. Figures 2(a), 2(b) and 2(c) are sketched for uniform density and we observed that circumferential stresses are maximum at internal surface and these stresses decreases with decrease in thickness radially. As with the increase in density radially, circumferential goes on decreasing which gives appropriate design of the disk. It has been observed from Figures 3(a), 3(b), and 3(c) that for uniform thickness, circumferential stress is maximum at internal surface. Also, circumferential stress for the disk whose density decreases radially is maximum as compared to the disk with high density. Also, decrease in thickness radially yields decrease in circumferential stress, which gives appropriate design of disk.

(a)

(b)

Figure 1: Radial stress () and circumferential stress () for annular disk with exponentially variable thickness and exponentially variable density with and under (a) 500 rad/s and (b) 780 rad/s for both approaches.

(a)

(b)

(c)

Figure 2: Radial stress () and circumferential stress () for annular disk with exponentially variable thickness and exponentially variable density with (a) and , (b) , and (c) and , under 500 rad/s for both approaches.

(a)

(b)

(c)

Figure 3: Radial stress () and circumferential stress () for annular disk with exponentially variable thickness and exponentially variable density with (a) and , (b) and , and (c) and  rad/s for both approaches.

From all previous analysis, we can conclude that disk whose thickness decreases radially and density increases radially is on the safer side of the design as compared to other thickness and density parameters, because circumferential stress is less for the disk whose thickness decreases radially and density increases radially as compared to other thickness and density parameters.

4.2. Rotating Annular Disk with Variable Thickness and Variable Density Using Power Law

It has been observed from Figures 4(a) and 4(b) that circumferential stress is maximum at internal surface. Circumferential stress is maximum for the disk whose density is less and thickness is high as compared to the disk whose thickness is less and density is high. Also, it has been observed from the figure that with the increase in angular speed, circumferential stress goes on increasing, but again circumferential stress is less as compared to other cases. It has been observed from Figures 5(a), 5(b), and 5(c) that for uniform density disk, circumferential stress is maximum at internal surface. With the decrease in thickness radially, circumferential stress goes on decreasing. Also with the increase in density, circumferential stress further decreases which provides an appropriate design of the disk. It has been observed from Figures 6(a), 6(b), and 6(c) that for uniform thickness, circumferential stresses are maximum at internal surface. Also, it has been noted that circumferential stress is less for disk profile with high density as compared to disk profiles with less density. With the decrease in thickness radially, circumferential stress decreases which yields that disk, whose thickness is less and density is high, is appropriate for the design as compared to other disk profiles.

(a)

(b)

Figure 4: Radial stress () and circumferential stress () for annular disk with variable thickness and variable density with and under (a) 500 rad/s and (b) 780 rad/s for both approaches.

(a)

(b)

(c)

Figure 5: Radial stress () and circumferential stress () for annular disk with variable thickness and variable density with (a) and , (b) and , and (c) and , under 500 rad/s for both approaches.

(a)

(b)

(c)

Figure 6: Radial stress () and circumferential stress () for annular disk with variable thickness and variable density with (a) and , (b) and , and (c) and , under 500 rad/s for both approaches.

From the previous analysis, we observed that disk with high density and less thickness is on the safer side of the design as compared to other parameters because circumferential stress is less for previous case as compared to other disk profiles. The results calculated for variable thickness and variable density using finite difference method are in very good agreement with Runge-Kutta’s method studied by You et. al [9] for all the rotating disks.

4.3. Annular Disk with Constant Thickness and Constant Density

It has been observed from Figures 7(a) and 7(b) that for constant thickness and constant density, circumferential stress is maximum at internal surface. Also, with the increase in angular speed, circumferential stress increases significantly.

(a)

(b)

Figure 7: Radial stress () and circumferential stress () for annular disk with constant thickness and constant density under (a) 500 rad/s and (b) 780 rad/s for both approaches.

5. Conclusion

After analyzing all the three disk profiles, it can be concluded that circumferential stresses are maximum at internal surface. It is also concluded that disk whose thickness decreases radially and density increases radially is on the safer side of design as compared to the disk with exponentially varying thickness and exponentially varying density as well as to flat disk. This is because circumferential stresses for the disks whose thickness decreases and density increases radially are less as compared to the disk with exponentially varying thickness and exponentially varying density as well as to flat disk.

References

S. P. Timoshenko and J. N. Goodier, Theory of Elasticity, McGraw-Hill, New York, NY, USA, 3rd edition, 1970.

S. C. Ugural and S. K. Fenster, Advanced Strength and Applied Elasticity, Elsevier, New York, NY, USA, 1987.

R. Hill, The Mathematical Theory of Plasticity, Oxford University Press, Oxford, UK, 1998.

S. Sharama, M. Sahni, and R. Kumar, “Elastic-plastic analysis of a thin rotating disk of exponentially variable thickness with inclusion,” WSEAS Transactions on Mathematics, vol. 9, no. 5, pp. 314–323, 2010.

U. Gamer, “Elastic-plastic deformation of the rotating solid disk,” Ingenieur-Archiv, vol. 54, no. 5, pp. 345–354, 1984.

U. Güven, “Elastic-plastic stresses in a rotating annular disk of variable thickness and variable density,” International Journal of Mechanical Sciences, vol. 34, no. 2, pp. 133–138, 1992.

L. H. You, S. Y. Long, and J. J. Zhang, “Perturbation solution of rotating solid disks with nonlinear strain-hardening,” Mechanics Research Communications, vol. 24, no. 6, pp. 649–658, 1997.

L. H. You and J. J. Zhang, “Elastic-plastic stresses in a rotating solid disk,” International Journal of Mechanical Sciences, vol. 41, no. 3, pp. 269–282, 1999.

L. H. You, Y. Y. Tang, J. J. Zhang, and C. Y. Zheng, “Numerical analysis of elastic-plastic rotating disks with arbitrary variable thickness and density,” International Journal of Solids and Structures, vol. 37, no. 52, pp. 7809–7820, 2000.

S. C. Sterner, S. Saigal, W. Kistler, and D. E. Dietrich, “A unified numerical approach for the analysis of rotating disks including turbine rotors,” International Journal of Solids and Structures, vol. 31, no. 2, pp. 269–277, 1994.