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.


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.

Temperature Analysis
Consider a homogeneous isotropic hollow disk of inner and outer radii   and   , respectively.One-dimensional transient 2 Chinese Journal of Engineering 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   (,  = 1, 2) and   ( = 1, 2) 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  1 and (17)

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  0 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 Chinese Journal of Engineering Therefore, where (32)

Results and Discussion
Consider a disk with   = 0.4 m and   = 0.6 m.The modulus of elasticity and Poisson's ratio, respectively, have values of  = 200 GPa and ] = 0.3.The applied internal pressure is 70 MPa.The mass density is  = 7854 Kg/m 3 .The thermal diffusivity, thermal conductivity, and specific heat capacity are  = 1.17(10 −5 )1/ ∘ C,  = 60.5 W/m⋅K, and  = 434 J/Kg⋅K, 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: = 25 (26.85 −  (  , )) Temperature distribution for this thermal boundary condition is where Figure 1 indicates temperature distribution for the course of 20 seconds for  = 0.45 m. Figure 2 illustrates temperature distribution in the radial direction at  = 3 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  = 0.45 m.When  → ∞, the variation of temperature versus dimensionless radius becomes linear.
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.
Figures 8, 9, and 10 show the distribution of radial displacement, radial stress, and hoop stress versus time.
Figures 11, 12, 13, 14, 15, and 16 show comparison of the displacement and stresses distribution for different time and  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 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  = 0.45 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 (rad/sec).)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.

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 Chinese Journal of Engineering      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.

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

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