In this paper, the exact analytical and numerical solutions for rotating variable-thickness annular disk are presented. The inner and outer edges of the rotating variable-thickness annular disk are considered to have free boundary conditions. Two different annular disks for the radially varying thickness are given. The numerical Runge-Kutta solution as well as the exact analytical solution is available for the first disk while the exact analytical solution is not available for the second annular disk. Both exact and numerical results for stress function, stresses, strains and radial displacement will be investigated for the first annular disk of variable thickness. The accuracy of the present numerical solution is discussed and its ability of use for the second rotating variable-thickness annular disk is investigated. Finally, the distributions of stress function, displacement, strains, and stresses will be presented. The appropriate comparisons and discussions are made at the same angular velocity.
The annular disks have received a great deal of attention because of their widely used in many mechanical and electronic devices, such as computer disk memory units, circular saws, and turbine rotors. The problems of rotating annular disks can be easily found in most of the standard elasticity books [1,2]. Most of the research works are concentrated on the analytical solutions of rotating disks with simple cross-section geometries of uniform thickness and especially variable thickness [3-8]. In [
As many rotating components in use have complex cross-sectional geometries, they cannot be dealt with using the existing analytical methods. Numerical methods, such as the finite element method [
This work is a predictive assessment of the stresses in and deformation of a rotating annular disk with variable thickness variation. Two methods of analysis, namely, the analytical method and Runge-Kutta numerical method of governing differential equations were used. This study concentrated first on the analytical solution for the rotating annular disk with arbitrary cross-section of continuously variable thickness. A unified governing equation will be first derived from the basic equations of the rotating disks and the proposed stress-strain relationship. Next, Runge-Kutta method is introduced to solve the governing equation. A comparison between both analytical and numerical solutions is made. The accuracy of the numerical solution is used to find the stress function, stresses, strains and displacement of rotating variable-thickness annular disk whose analytical solution is not available. Finally, a number of numerical examples are given to demonstrate the validity of the proposed method.
As the effect of thickness variation of rotating annular disks can be taken into account in their equation of motion, the theory of the variable-thickness annular disks can give good results as that of the uniform-thickness annular disks as long as they meet the assumption of plane stress. After considering this effect, the equation of motion of rotating disks with variable thickness can be written as
where and are the radial and circumferential stresses, r is the radial coordinate, is the density of the rotating disk, is the constant angular velocity, and h is the thickness which is function of the radial coordinate r.
The relations between the radial displacement u and the strains are irrespective of the thickness of the rotating disk. They can be written as
where and are the radial and circumferential strains, respectively. The above geometric relations lead to the following condition of deformation harmony:
For the elastic deformation, the constitutive equations for the rotating disk can be described with Hooke’s law
where E is Young’s modulus and is Poisson’s ratio. Introducing the stress function and assuming that the following relations hold between the stresses and the stress function
Substituting Equation (5) into Equation (4), one obtains
The substitution of Equation (6) into Equation (3) produces the following confluent hypergeometric differential equation for the stress function
The thickness of the annular disk is assumed to vary nonlinearly through the radial direction. It is assumed to be in terms of a simple exponential power law distribution according to the following two types:
Disk 1:
Disk 2:
For both disks, is the thickness at the inner edge of the disk, n and k are geometric parameters, a is the inner radius of the disk and b is the outer radius of the disk (see Figures 1 and 2). The value of n equal to zero represents a uniform-thickness annular disk. It is to be noted that the parameter n determines the thickness at the outer edge of the annular disk relative to while the parameter k determine the shape of the profile. The thickness of the annular disk 1 may be decreases at the outer edge (see
For small k and large n (case 3) the profile of the annular disk 1 is concave while it is convex for large k and small n (case 2). For large k and n (case 1) the profile of the annular disk 1 may be changing its shape from convexity to concavity. In addition, the profile of the annular disk 2 is convex for small k and large n (case 3) while it is concave for large k and small n (case 2). Also, for large k and n (case 3) the profile of the annular disk 2 may change its shape from concavity to convexity.
To simplify the solving process, we introduce the following dimensionless variables:
Then, Equation (7) may be written in the following simple forms according to the two cases:
for annular disk 1, and
for annular disk 2.
The exact general solution for Equation (11) may be available while it is not available for Equation (12). The modified Runge-Kutta numerical solution may be available for both cases. Firstly, we will get both the exact and numerical solutions for Equation (11). If the numerical solution is compared well with the exact one, we can use it to solve Equation (12) for the second case.
The exact general solution for Equation (11) can be written as
,(13)
where and are arbitrary constants and and are Whittaker’s functions (see Abramowitz and Stegun [
in which
It is to be noted that Whittaker’s functions may be given by
where H and K are the hypergeometric and Kummer’s functions, respectively, with n > 0. It is to be noted that for real values of i and j, Whittaker’s functions and converge for.
The term in Equation (13) can be written as
where
The substitution of Equation (13) into Equation (5) with the aid of the dimensionless forms given in Equation (10) gives the radial and circumferential stresses in the following forms:
Note that the first differentiations of Whittaker’s functions and are given by
Finally, the dimensionless strains and the corresponding radial displacement may be obtained easily using Equation (19) as well as the dimensionless forms given in Equation (10). So, all of stress function, stresses, strains, and radial displacement can be completely determined under the traction conditions on the inner and outer surfaces of the annular disk. They can be expressed as
The modified Runge-Kutta algorithm may be used here to solve both the differential Equations (11) and (12) with the aid of the boundary conditions given in Equation (21). Equations (11) or (12) can be written in the following general form:
where the prime denotes differentiation with respect to R and s denotes the disk number. The functions may be given by
Runge-Kutta’s iterative formulae for the second-order differential equation are [11-14]
where is the increment of the distance along the radial direction of the rotating annular disks, and can be determined with the following relations
The numerical simulation starts from the inner boundary, where a trial value of the first derivative of the stress function is assumed. With a small distance, the stress function and its first derivative at the new position can be obtained using Equation (25) and the radial stress is calculated. According to the difference between the computed radial stress and the known radial stress at the outer boundary, the initial trial value of the first derivative of the stress function at the inner boundary is modified and the next iteration is carried out in the same way. This iterative process is performed until both the boundary conditions are simultaneously satisfied. Once the stress function is obtained, the stresses, strains, and displacement in the rotating annular disks can be obtained using Equations (2), (5) and (6) with the aid of the dimensionless given in Equation (10). This process may be easily done after getting the continuous form of the stress function using the curve fitting and least square method. So, all other quantities can be easily determined.
Some numerical examples for the rotating variablethickness annular disk will be given according the exact and numerical solutions. Results determined as per the exact analytical solution are compared with those obtained by the numerical Runge-Kutta (R-K) solution in Figures 3-10 for disk 1. Additional results for stress function, radial displacement, strains and stresses of the rotating variable-thickness annular disk 2 using the numerical R-K solution are also presented in Figures 11-16. The inner and outer radii of the disk are taken to be a = 0.1 b (R = A = 0.1) and b (R = 1), and the results are given in terms of the rotating angular velocity.
The numerical applications will be carried out for the radial displacement, strains, stress function and stresses
that being reported herein, with help of Equation (10), are in the following dimensionless forms:
The distribution of the stress function and stresses are presented in
The distributions of the stress functions, radial stress, circumferential stress, radial displacement, radial strain as well as circumferential strain through the radial direction of the annular disk 1 are all presented in Figures 5-10, respectively, for different cases of the geometric parameters k and n.
The radial stress and circumferential stress in the rotating variable-thickness annular disk 1 are plotted in Figures 7 and 8. The maximum value of does not occur at the mid-plane of the disk 1. It occurs at R = 0.33, R = 0.43, R = 0.35 for cases 2, 3 and 1, respectively. The maximum value of occurs at the inner surface of the disk 1 while the smallest value occurs at outer surface. Case 1 gives the smallest stresses. The same conclusion may be applied on the plots of strains in Figures 9 and 10.
It can be seen from Figures 3-10 that the R-K method can describe all results through-the-thickness of the rotating annular disk 1 very well enough. This puts into evidence the great role played by R-K method in the modeling of rotating variable-thickness annular disks. So, it will trustily used to find the solution for Equation (12) of the rotating variable-thickness annular disk 2. As a result, the stresses and stress function in the rotating variable-thickness disk 2 are plotted in Figures 11-13 according to all cases. However, the radial displacement and strains are plotted in Figures 14-16.
This paper presents an exact analytical solution for a rotating variable-thickness annular disk. Also, it presents a unified numerical method based on Runge-Kutta forth-order method for the elastic calculation of different rotating annular disks with a general, arbitrary configuration. The governing equation was derived from the equilibrium equation and the stress-strain relationship. The calculation of the rotating annular disk was turned into finding the solution of second-order differential equations under the given conditions at two boundary fixed points. Runge-Kutta method algorithm was introduced to solve the governing equation for two types of annular disks, in which the exact solution of one of them only is available. A number of numerical examples were studied. The results from the exact analytical and R-K solutions were compared. The proposed R-K approach gives very agreeable results to the exact analytical solution. The modified R-K method presented here may be trustily used for the boundary value problems that their exact solutions are not available.