Chinese Journal of Engineering

Volume 2014 (2014), Article ID 894902, 13 pages

http://dx.doi.org/10.1155/2014/894902

Research Article

Transient Thermoelastic Analysis of Pressurized Rotating Disks Subjected to Arbitrary Boundary and Initial Conditions

Mohammad Zamani Nejad and Azam Afshin

Mechanical Engineering Department, Yasouj University, P.O. Box 75914-353, Yasouj, Iran

Received 6 December 2013; Accepted 9 February 2014; Published 17 April 2014

Academic Editors: X. Li and B. Sun

Copyright © 2014 Mohammad Zamani Nejad and Azam Afshin. 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.

This paper focuses on exact analytical solution of transient thermoelastic behaviors of rotating pressurized disks subjected to arbitrary boundary and initial conditions. The pressure, inner radius, and outer radius are considered constant. The basic thermoelasticity theory under generalized assumptions is used to solve the thermoelastic problem. Using the method of the separation of variables, the relations of temperature and transient thermal stresses in the radial direction are obtained. In the case study, the disk is considered under heat flux. Some useful discussions and numerical examples are presented. The analytical results were compared with those of the finite element method and good agreement was found. The relations obtained in this paper can be applied to any arbitrary boundary and initial conditions.

1. Introduction

The rotating disks subjected to mechanical and transient thermal loads are widely used in many engineering fields such as aerospace, mechanical, naval, nuclear energy, chemical plant, electronics, and biomaterials. Therefore, much research has been conducted in this field.

Cheung et al. [1] investigated the problem of transient thermal stresses in a solid elastic homogeneous and isotropic sphere for uniform and nonuniform local surface heating. Sugano [2] analyzed the problem of transient thermal stresses in a homogeneous transversely isotropic finite cylinder due to an arbitrary internal heat generation. A transient response of one-dimensional axisymmetric quasistatic-coupled thermoelastic problems was studied by Yang et al. [3]. Transient thermal stresses in cylindrically orthotropic tubes were obtained by Kardomateas [4]. The transient thermal stresses in a homogeneous cylindrically orthotropic hollow cylinder due to a constant temperature were imposed on one surface and heat convection into a medium at the other surface was studied by Kardomateas [5]. Ashida et al. [6] presented a general solution technique for transient thermoelasticity of transversely isotropic solids in cylindrical coordinates. For quasistatic problems, and without considering the effect of inertial term, Yee and Moon [7] obtained the closed-form solutions of the orthotropic hollow cylinder subjected to an arbitrary axisymmetric temperature distribution. Assuming that the material properties vary nonlinearly in the radial direction and that Poisson’s ratio is constant, Zamani Nejad and Rahimi [8] obtained closed-form solutions for one-dimensional steady-state thermal stresses in a rotating functionally graded pressurized thick-walled hollow circular cylinder. In another study, they obtained a complete and consistent three-dimensional set of field equations by tensor analysis to characterize the behavior of FGM thick shells of revolution with arbitrary curvature and variable thickness along the meridional direction [9]. In this work, transient thermoelastic analysis of homogeneous and isotropic rotating disks under general boundary conditions was obtained.

2. Temperature Analysis

Consider a homogeneous isotropic hollow disk of inner and outer radii and , respectively. One-dimensional transient heat conduction equation of disk, in the absence of heat source, for axisymmetric case is

where is the temperature distribution, is thermal diffusivity, and , , are the thermal conductivity, mass density, and specific heat capacity, respectively.

The initial and boundary conditions are as follows:

where    and    are the constants dependent on the thermal boundary condition, and is the given initial condition.

The solution of (1) may be obtained by the method of separation of variables, generalized Bessel function, and eigenfunction method as

The boundary conditions are defined as follows:

Integrating (4) twice yields

Using the boundary conditions (see (6)) to determine the constants and yields

or

For , boundary and initial conditions are as follows:

Thus, the solution for is expressed as

where

Here are eigenvalue and are positive roots of the following equation:

Here and are Bessel functions of the first and the second kinds.

Eigenfunction is as follows:

and parameters are constant and are defined as

or

is the norm of eigenfunction and are obtained as

3. Transient Thermoelastic Analysis

Consider thick hollow disk in the previous section at a constant angular velocity about its central axis. In the cylindrical coordinate system , for the axisymmetric problem, radial and circumferential strains are as follows:

Here, radial displacement is the only nonzero component of the displacement and is only the function of radial coordinate .

For homogenous and isotropic disk, the linear thermoelastic constitutive relations can be expressed as follows:

where is modulus of elasticity, is the coefficient of linear thermal expansion, is Poisson’s ratio and and are components of stresses in radial and hoop directions, respectively, and is the reference temperature. Although, in this paper the reference temperature is taken as zero, the analysis would be valid for any nonzero value.

The equilibrium equation of the rotating hollow disk, in the absence of body forces, is expressed as

where is angular velocity.

Substituting (18) into (19)-(20) and utilizing (22) yield

Integrating (23) twice yields

Substituting (24) into (19)-(20), the radial and circumferential stresses of the rotating thick hollow disk are obtained as

where

The mechanical boundary conditions are expressed as

Substituting boundary conditions into (25) gives

Therefore,

where

4. Results and Discussion

Consider a disk with  m and  m. The modulus of elasticity and Poisson’s ratio, respectively, have values of  GPa and . The applied internal pressure is  MPa. The mass density is  Kg/m³. The thermal diffusivity, thermal conductivity, and specific heat capacity are ,   W/mK, and  J/KgK, respectively. In the following, the analytical solution described in the preceding section may be checked at two examples.

Example 1. Assume that the inner surface of disk is under heat flux and the outer surface of disc is under convection. The initial condition is a linear function in terms of the radius. Thus, thermal boundary condition is expressed as follows:

Temperature distribution for this thermal boundary condition is

where

Figure 1 indicates temperature distribution for the course of 20 seconds for  m. Figure 2 illustrates temperature distribution in the radial direction at  sec.

Figure 1: Temperature distribution ( m).

Figure 2: Temperature distribution in the radial direction ( sec).

Figure 3 shows the comparison of the temperature distribution for different radii. As can be seen, temperature increases as radius increases. Figure 4 compares the temperature distribution for different times at point  m. When , the variation of temperature versus dimensionless radius becomes linear.

Figure 3: Comparison of the temperature distribution for different radius.

Figure 4: Comparison of the temperature distribution for different times.

Figures 5, 6, and 7 show the distribution of radial displacement, radial stress, and hoop stress versus radial direction. Meanwhile, this result is compared with FEM.

Figure 5: Radial displacement distribution ( sec).

Figure 6: Radial stress distribution ( sec).

Figure 7: Hoop stress distribution ( sec).

Figures 8, 9, and 10 show the distribution of radial displacement, radial stress, and hoop stress versus time.

Figure 8: Radial displacement distribution ( m).

Figure 9: Radial stress distribution ( m).

Figure 10: Hoop stress distribution ( m).

Figures 11, 12, 13, 14, 15, and 16 show comparison of the displacement and stresses distribution for different time and different radius. Figures 14 and 16 show that variation of radial and hoop stresses into variation of time is very small. Figures 17, 18, 19, 20, 21, and 22 show comparison of the displacement and stresses distribution for different angular velocity. (In these figures, is angular velocity and its dimension is radian per second ().) In Figure 20 it is seen that angular velocity increases as radial stress decreases and if angular velocity increases more, radial stress changes from pressure to tension.

Figure 11: Comparison of the radial displacement distribution for different radius.

Figure 12: Comparison of the radial displacement distribution for different times.

Figure 13: Comparison of the radial stress distribution for different radius.

Figure 14: Comparison of the radial stress distribution for different times.

Figure 15: Comparison of the hoop stress distribution for different radius.

Figure 16: Comparison of the hoop stress distribution for different times.

Figure 17: Comparison of the radial displacement distribution for different angular velocity ( m).

Figure 18: Comparison of the radial displacement distribution for different angular velocity ( sec).

Figure 19: Comparison of the radial stress distribution for different angular velocity ( m).

Figure 20: Comparison of the radial stress distribution for different angular velocity ( sec).

Figure 21: Comparison of the hoop stress distribution for different angular velocity ( m).

Figure 22: Comparison of the hoop stress distribution for different angular velocity ( sec).

Example 2. Thermal boundary and initial conditions for the second example are

Temperature distribution for this thermal boundary condition is

Figure 23 indicates temperature distribution for the course of 20 seconds and Figure 24 illustrates temperature distribution in thickness; meanwhile the result of this state is compared with FEM.

Figure 23: Temperature distribution ( m).

Figure 24: Temperature distribution in the radial direction ( sec).

Figure 25 shows comparison of the temperature distribution for different radius. As can be seen, temperature increases as radius increases. Figure 26 shows comparison of the temperature distribution for different times at point   m. In the initial seconds and toward the inner radius the variation of temperature is intensive.

Figure 25: Comparison of the temperature distribution for different radius.

Figure 26: Comparison of the temperature distribution for different times.

When the variation of temperature versus dimensionless radius becomes linear.

Figures 27, 28, 29, 30, 31, and 32 show radial displacement and radial stress and hoop stress distribution versus time and thickness; meanwhile this result is compared with FEM. Figure 32 shows that in the course of 20 seconds hoop stress graphs have a minimum point. Figures 33, 34, 35, 36, 37, and 38 show comparison of the displacement and stresses distribution for different time and different radius. Figures 36 and 38 show that variation of radial and hoop stresses into variation of time is very small. Figures 39, 40, 41, 42, 43, and 44 show comparison of the displacement and stresses distribution for different angular velocity. (In these figures, is angular velocity and its dimension is radian per second ().) In Figure 42 it is seen that angular velocity increases as radial stress decreases and if angular velocity increases more, radial stress changes from pressure to tension.

Figure 27: Radial displacement distribution ( sec).

Figure 28: Radial displacement distribution ( m).

Figure 29: Radial stress distribution ( sec).

Figure 30: Radial stress distribution ( m).

Figure 31: Hoop stress distribution ( sec).

Figure 32: Hoop stress distribution ( m).

Figure 33: Comparison of the radial displacement distribution for different radius.

Figure 34: Comparison of the radial displacement distribution for different times.

Figure 35: Comparison of the radial stress distribution for different radius.

Figure 36: Comparison of the radial stress distribution for different times.

Figure 37: Comparison of the hoop stress distribution for different radius.

Figure 38: Comparison of the hoop stress distribution for different times.

Figure 39: Comparison of the radial displacement distribution for different angular velocity ( m).

Figure 40: Comparison of the radial displacement distribution for different angular velocity ( sec).

Figure 41: Comparison of the radial stress distribution for different angular velocity ( m).

Figure 42: Comparison of the radial stress distribution for different angular velocity ( sec).

Figure 43: Comparison of the hoop stress distribution for different angular velocity ( m).

Figure 44: Comparison of the hoop stress distribution for different angular velocity ( sec).

5. Conclusions

In the present study, an exact solution procedure has been presented for transient thermoelastic analysis of isotropic and homogeneous hollow disks under general initial and boundary conditions. Some numerical results are shown in the figures. The distribution of the transient temperature, radial displacement, and thermal stresses through the radial direction of the disk is plotted and the results are compared with the solution using finite element method (FEM), which showed good agreement.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

J. B. Cheung, T. S. Chen, and K. Thirumalai, “Transient thermal stresses in a sphere by local heating,” Journal of Applied Mechanics, vol. 41, no. 4, pp. 930–934, 1974.

Y. Sugano, “Transient thermal stresses in a transversely isotropic finite circular cylinder due to an arbitrary internal heat-generation,” International Journal of Engineering Science, vol. 17, no. 8, pp. 927–939, 1979.

Y.-C. Yang, T.-S. Wang, and C.-K. Chen, “Thermoelastic transient response of an infinitely long annular cylinder,” Journal of Thermal Stresses, vol. 9, no. 1, pp. 19–30, 1986.

G. A. Kardomateas, “Transient thermal stresses in cylindrically orthotropic composite tubes,” Journal of Applied Mechanics, vol. 56, no. 2, pp. 411–417, 1989.

G. A. Kardomateas, “Initial phase of transient thermal stresses due to general boundary thermal loads in orthotropic hollow cylinders,” Journal of Applied Mechanics, vol. 57, no. 3, pp. 719–724, 1990.

F. Ashida, N. Noda, and I. A. Okumura, “General solution technique for transient thermo elasticity of transversely isotropic solids in cylindrical coordinates,” Acta Mechanica, vol. 101, no. 1–4, pp. 215–230, 1993.

K.-C. Yee and T. J. Moon, “Plane thermal stress analysis of an orthotropic cylinder subjected to an arbitrary, transient, asymmetric temperature distribution,” Journal of Applied Mechanics, vol. 69, no. 5, pp. 632–640, 2002.

M. Zamani Nejad and G. H. Rahimi, “Deformations and stresses in rotating FGM pressurized thick hollow cylinder under thermal load,” Scientific Research and Essays, vol. 4, no. 3, pp. 131–140, 2008.

M. Zamani Nejad, G. H. Rahimi, and M. Ghannad, “Set of field equations for thick shell of revolution made of functionally graded materials in curvilinear coordinate system,” Mechanika, vol. 77, no. 3, pp. 18–26, 2009.