Research Article | Open Access
Transient Thermoelastic Analysis of Pressurized Rotating Disks Subjected to Arbitrary Boundary and Initial Conditions
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.
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.  investigated the problem of transient thermal stresses in a solid elastic homogeneous and isotropic sphere for uniform and nonuniform local surface heating. Sugano  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. . Transient thermal stresses in cylindrically orthotropic tubes were obtained by Kardomateas . 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 . Ashida et al.  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  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  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 . 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/m3. 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 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.
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.
Example 2. Thermal boundary and initial conditions for the second example are Temperature distribution for this thermal boundary condition is
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.
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.
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.
- 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.
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.