Abstract
In practice, some components in large structures such as the connecting rods between the rotating parts in the engines, turbines, and so on, can model as beam structures rotating around the fixed axis and subject to the axial compression load; therefore, the study of mechanical behavior to these structures has a significant meaning in practice. This paper analyzes the vibration responses of rotating FGM beams subjected to axial compressive loads, in which the beam is resting on the two-parameter elastic foundation, taking into account the initial geometrical imperfection. Finite element formulations are established by using the new shear deformation theory type of hyperbolic sine functions and the finite element method. The materials are assumed to be varied smoothly in the thickness direction of the beam based on the power-law function with the porosity. Verification problems are conducted to evaluate the accuracy of the theory, proposed mechanical structures, and the calculation programs coded in the MATLAB environment. Then, a parameter study is carried to explore the effects of geometrical and material properties on the vibration behavior of FGM beams, especially the influences of the rotational speed and axial compressive load.
1. Introduction
Functionally graded material (FGM) is commonly made of ceramic and metal in a certain proportion of ingredients; this is a new material that fully carries the characteristics of the material components such as good heat resistance, corrosive resistance, and withstanding thermal stress. Therefore, studying the mechanical behavior of structures made of FG materials has great practical value, especially in designing, manufacturing, and using them in the most effective way. Also, large mechanical systems such as helicopter blades, aircraft, and ship propellers and even nanosize devices such as microturbines may include blade-like structures that occur in rotational movements, where they can be considered as rotating beam structures. Hence, investigations of rotating integrated beam structures play a very important role in science. There have been a number of research projects and exploration results on these issues. Saito and Otomi [1] researched the stability and vibration of a beam under axial and tangential loads, where the analysis is based on the Timoshenko beam theory, and the effects of the attached mass are expressed with Dirac delta functions. Esmailzadeh and Ohadi [2] introduced a solution for stability and vibration response analysis of a beam subjected to both tangential and axial loads using the Frobenius method and Timoshenko beam theory, in which the axial force is assumed to be tangential to the axis of the beam-column. Sloss et al. [3] used Green’s function to find the solution of the vibration control problem of piezoelectric beams subjected to axial load. Zhang et al. [4] studied the vibration problem of carbon nanotubes subjected to axial load using Bernoulli–Euler theory. Alkharabsheh and Younis [5] investigated the effect of axial forces on the static behavior and the fundamental natural frequency of electrostatically actuated MEMS arches using Euler–Bernoulli theory. Babilio [6] employed the Galerkin procedure to study the dynamics of an FGM beam under axial load, in which Taylor expansion was used to establish the partial differential equation. Barry et al. [7] gave out equations of motion of a beam under axial loads; then, the solution of the vibration response of this beam was introduced, in which this beam was attached to a mass-spring-mass system. Farokhi and Ghayesh [8] researched the nonlinear bending and resonant behavior of a microarch based on Euler–Bernoulli theory and the Galerkin method. Ghannadiasl and Mofid [9] used Green function to investigate the free vibration of a beam subjected to axial loads using Timoshenko beam theory and an analytical method. Nejati et al. [10] employed the generalized differential quadrature approach to explore the buckling and vibration behavior of the FGM beam reinforced with carbon nanotubes. A two-node C1 beam element was used by Thuc et al. [11] to analyze the free vibration of composite beams under axial load using a four-unknown shear and normal deformation theory. Waddar et al. [12] showed an experimental study of buckling and the dynamic response of the cenosphere reinforced epoxy composite core sandwich beam with sisal fabric/epoxy composite facings subjected to axial load. Zhao et al. [13] presented the exact solutions of the forced transverse vibration of a coupled Timoshenko double-beam system under compressive axial load based on Timoshenko beam theory and Green’s functions; the solution could be suitable for any boundary conditions. Szymczak and Kujawa [14] presented an exact formulation of the frequency for free torsional vibration analysis of laminated beams subjected to axial load. Xu et al. [15] used the energy principle in conjunction with the Rayleigh–Ritz procedure to investigate the flexural vibration and its power flow of an axially loaded beam with arbitrary boundary and nonuniform elastic foundation. Brahimi-Mamaghani et al. [16] employed an analytical approach to finding the solution of dynamic and vibration stability of a viscoelastic FGM beam under axial load, where the Laplace transform and a Galerkin discretization scheme were combined to capture the solution. Khosravi et al. [17] used the first-order shear deformation plate theory to investigate the nonaxisymmetric buckling behaviour of isotropic homogeneous annular plates subjected to simultaneous effects of uniform temperature rise and constant angular speed.
Recently, the explorations of mechanical responses of rotating beam structures have also been in interest and researched by scientists. Khosravi et al. [18] employed Timoshenko beam theory and von Kármán type of kinematic assumptions to examine the influence of uniform temperature elevation on the vibration of rotating composite beams reinforced with carbon nanotubes which may lead to instability. Timoshenko beam theory and von Kármán type of kinematic assumptions were also used by Khosravi and his coworkers [19] to study the thermal buckling of rotating carbon nanotube-reinforced composite beams. Kiani et al. [20] investigated the buckling behaviour of an isotropic homogeneous rotating annular plate subjected to uniform compression on both inner and outer edges. Bagheri et al. [21] researched the free vibration treatment of pre- and postbuckled rotating functionally graded beams using the Euler–Bernoulli beam theory. Hosseini et al. [22] analyzed the buckling and postbuckling of rotating clamped-clamped functionally graded beams due to the interaction between the thermal environment and the rotation based on the Euler–Bernoulli beam theory. From these works, one can see that the boundary condition has a significant influence on the mechanical responses of the rotating beam, especially to the natural frequencies as well as the vibration mode shapes of this structure. In addition, boundary details and material degradation also have a significant effect on these responses, which can be better understood in [23, 24].
In the fabrication process, there can be many types of defects, in which the porosities are easy to appear, dramatically altering the mechanical response of these structures. There have been a number of scientific studies on the free and forced vibration responses and static bending of beams and plates made of FG materials with pore defects [25–32].
From the published research studies, readers can find that previous studies only mentioned the mechanical behavior of porosity beam structures made of FG materials, whether or not resting on elastic foundations, while some works investigated the effect of initial compressive force with rotational movement around one fixed axis, where these studies have not fully presented beams with both pore defects and imperfections subjected to both axial compression and elastic foundation effect. Besides, this proposed mathematical model is close to the technical reality models, which can help to guide in designing, manufacturing, and using them effectively.
Therefore, this work is about to focus on these proposed problems by using the finite element method and a new type of shear deformation theory in the form of hyperbolic sine functions. This is a simple theory, which takes into account the effect of shear strain without any shear correction factors, and it is an advantage of computational theory in this work. This is the first work using this beam theory to solve the problem of free vibration of rotating FGM beams considering both shape and material defects, thereby giving the influence of different parameters on the free vibration response of this structure, especially the influence of the rotational movement. This is the new point of this work, in which the calculation results are of great value for the designers and users as mentioned above.
The rest of this work is structured as follows. Finite element formulations for free vibration analysis of rotating FGM beams resting on two-parameter elastic foundations are fully presented in Section 2, in which the initial geometrical imperfection and the porosity are considered. Verification examples are introduced in Section 3. Numerical results and discussion are abundantly presented in Section 4. Section 5 draws out some main important highlight explorations of this paper.
2. Finite Element Model of Rotating FGM Beams with Geometrical Imperfection and Porosity
This work is about to focus on the vibration analysis of the FGM beam rotating with the rotational speed , in which one side of the beam is at distance r from a fixed axis The beam has an initial geometrical imperfection and the porosity. The structure is subjected to an axial compressive load N0 and resting on the two-parameter elastic foundation with two parameters and as shown in Figure 1. The geometrical dimensions are the length L, width b, and thickness h.

(a)

(b)

(c)
The beam is made from two main components including ceramic (denoted as c) and metal (denoted as m) with the corresponding proportions Vc and Vm, which are calculated by the following function [33–38]:in which z is the thickness coordinate variable with and n is the volume fraction gradient index and its variation.
Then, the mechanical properties are changed in the z-direction based on the power-law distribution aswhere E is Young’s modulus, is the density, is Poisson’s ratio, and is a porosity volume fraction ( < 1).
There have been many shear deformation theories to investigate the mechanical responses of beam structures, where each one has its advantages. Recently, a new type of shear deformation theory has been developed and used widely, which is hyperbolic sine functions [39, 40]. This work also employs this new theory to establish the finite element equations of free vibration problems of rotating FGM beams. So, the displacement field at any point (x, z) within the beam has the following form:where
The expressions of strain components of the beam are calculated asin which
For the linear problem, the stress-strain relation according to Hooke’s law is expressed as follows:where Tx is the normal stress and is the shear stress.
The elastic strain generates the elastic energy of the FGM beam which is calculated as
The work done by the external compressive load is expressed as follows [41]:
The energy of the elastic foundation is calculated as
When the FGM beam is rotating around one fixed axis with the rotational speed , the potential energy of this beam generated by the rotational movement is calculated as [41, 42]with the centrifugal force [42].and is the density of the material.
The kinetic energy of the beam is expressed as
To determine the equilibrium equation of the FGM beam, this work uses Hamilton’s principle as follows:
Herein, this paper uses a two-node beam element, where each node contains five degrees of freedom:in which the displacement components at each point in the beam element are calculated through Lagrange and Hermite interpolation functions Ni and Hi:
Equation (15) can be rewritten in the matrix form as follows:
Strain components are expressed according to the nodal displacement as
And now, the elastic energy of the FGM beam element has the following expression:
Equation (18) can be rewritten in the matrix form as
The energies of the centrifugal inertia force and the elastic foundation are defined as follows:
The energy expression due to the axial compressive load has the following form:
The kinetic energy expression of the beam element is calculated as follows:in which is the element stiffness matrix.
Substituting equations (18)–(20) into equation (13), the equation to get the fundamental frequency, as well as the vibration mode shape of the rotating FGM beam, is obtained as follows:
Equation (23) demonstrates that all components relating to rotational speed, elastic foundation parameter, initial geometrical imperfection, and axial compressive loads are presented in the equation for finding the fundamental vibration response of the FGM beam, which is completely different from traditional problems; this makes the computation more difficult in comparison with previous publications.
Each side of the beam is under the following boundary conditions: Simply supported (denoted as S): Clamped (denoted as C): This paper employs three boundary conditions for calculations of FGM beams. One side is clamped, and the other side is free: C-F. Fully simply supported beam: S-S. Fully clamped beam: C-C.
3. Verification Examples
This section conducts five examples to evaluate the accuracy of the proposed theory and mathematical models, in which the numerical data of these examples are compared with those of other trustful publications.
Example 1. This example compares the frequency of the homogeneous beam under the axial compressive load. Geometrical and material properties are the length L, cross section bxh, L/h = 100, E = 380 GPa, = 0.3, and = 3960 kg/m3. The beam is subjected to an axial compressive load N0 with = 0.6. The nondimensional frequency is defined by the formula (with ). The results obtained from this work, exact expressions [43, 44], and FEM modeling [45] are listed in Table 1, in which they are calculated with different meshes. It can be seen that when the mesh size increases, the calculation results ensure the necessary convergence. The computational results in this work with the 8-element mesh are different from the results of the same mesh in reference [45]. The reason is that document [45] used the first-order shear deformation theory of Timoshenko. And for the 8-element mesh, the accuracy is acceptable; therefore, this mesh will be used for all the following related investigations.
Example 2. Next, this problem compares the results of porosity FGM beams from this work and Chebyshev collocation solution (CCS) [46]. The ratio L/h obtains the values of 5, 10, and 20; the porosity parameter is 0 and 0.2. Geometrical and material characteristics of the fully clamped FGM beam are Ec = 380 GPa, = 0.3, = 3960 kg/m3, Em = 70 GPa, = 0.3, and = 2702 kg/m3. The comparative nondimensional frequencies are presented in Table 2. It can be observed that the data have a very good similarity.
Example 3. Now, the frequency of the fully simply supported FGM (Al/Al2O3) beam resting on the two-parameter elastic foundation is considered. Geometrical and material properties are the length L, thickness h, L/h = 100, width b, Em = 70 GPa, = 2702 kg/m3, Ec = 380 GPa, and = 396 kg/m3. Two nondimensional elastic foundation parameters are defined as follows:withThe nondimensional frequency of the beam is calculated as follows:Table 3 presents the comparative first frequencies of the beam between this work and the analytical method [47]. One can see that, with the 10 elements, the accuracy is acceptable; therefore, this work will use this mesh for all the following related explorations.
Example 4. Next, this example compares the frequency of the S-S beam with an initial geometrical imperfection. The beam has the following geometrical and material parameters: L = 288.7h, h = 0.02 m, b = 0.04 m, E = 971 GPa, and mass density = 2300 kg/m3. The initial imperfection of the beam is , in which J0 is the amplitude of the imperfection. The nondimensional frequency and the initial imperfection coefficient are defined as follows:The numerical data obtained from this work and the pseudo-arc-length continuation technique [48] in the case of increasing the value of are shown in Figure 2. It can be seen that the data meet a good agreement.

Example 5. Finally, this last example carries the verification of the nondimensional frequency of the cantilever rotating beam with the speed . Consider a beam with the following geometrical and material parameters: length L, thickness h = b = L/100, r = 0, E = 70 GPa, and mass density = 2700 kg/m3. The nondimensional frequency is normalized as follows:The first three nondimensional frequencies of the rotating beam with different values of the speed ratio obtained from this work, an exact solution [49], a new dynamic modeling method (DMM) [50], and an isogeometric analysis [42] are presented in Table 4.
4. Numerical Results
This section carries out the vibration analysis of the rotating FGM porosity beam resting on the two-parameter elastic foundation subjected to the preaxial compressive load, where the initial geometrical imperfection is taken into calculations. Consider a beam with the following geometrical and material properties: length L, cross section bxh, Ec = 380 GPa, = 0.3, = 3960 kg/m3, Em = 70 GPa, = 0.3, and = 2702 kg/m3. The imperfection of the beam is , in which G0 is the amplitude of the imperfection, and the imperfection ratio . The nondimensional frequency and other parameters are normalized as follows:
4.1. Effect of the Rotational Speed
Firstly, the influence of the rotational speed on the free vibration response of the rotating FGM beam is investigated. The rotational speed is increased gradually so that the coefficient changes in a range of 0 to 10; the volume fraction index n obtains the values from 0 to 10; other parameters are L/h = 20, imperfection ratio = 0.001, axial compressive load = 0.5, = 5, = 1, distance ratio r/L = 1, and = 0.2. Figures 3–5 present the dependence of the first frequency of the beam on the rotational speed ratio and the volume fraction index n. The first four vibration mode shapes of the S-S and C-C rotating FGM beams with the speed ratio = 0 and 8 are presented in Figures 6 and 7. These figures point out the following.

(a)

(b)

(a)

(b)

(a)

(b)

(a)

(b)

(a)

(b)
When the rotational speed increases, nondimensional frequencies of the beam in the case of both boundary conditions (S-S and C-C) increase. This is because as the rotational speed of the beam increases, the inertial force component increases; therefore, the stiffness matrix of the FGM beam is increased by the addition of this inertial component. The change of the first nondimensional frequency of the structure depends strongly on the boundary conditions and the volume fraction index n. For the S-S boundary, when n increases, the first frequency of the beam decreases. In contrast, for the C-C beam, when n increases, the first frequency of the beam decreases sharply at first, then its value reaches the minimum point, and it increases gradually.
Besides, the rotational speed has a strong effect on the vibration mode shape of the beam, especially the first vibration mode shape. Also, due to the appearance of the inertial force, the first vibration mode shape is no longer symmetrical through the position x = L/2. The influence of the rotational speed on the vibration mode shapes depends strongly on the boundary conditions; for the S-S boundary condition, the influence level is more pronounced than the beam under the C-C boundary condition.
4.2. Influence of the Axial Load
This section is about to capture the effect of the axial compressive load on the vibration behavior of the FGM beam. The axial compressive load is changed so that the nondimensional load obtains the values from 0 to 0.8. The dependence of the first frequency of the FGM beam on the axial load and the volume fraction index n is presented in Table 5, and the dependence of the first frequency on the axial load and the rotational speed is presented in Table 6. As the axial compressive load increases, the frequency of the FGM beam decreases; this is because as the compressive load increases, the stiffness matrix of this beam decreases. Table 7 shows that when the value of the compression load is increased, the first nondimensional frequency changes the most, while the higher-order frequencies change slightly.
4.3. Influence of the Thickness of the FGM Beam
Consider an FGM beam with L/h varying from 5 to 80, imperfect parameter = 0.001, = 0.5, = 5, = 1, distance ratio r/L = 1, and volume fraction index n = 0.2. The dependence of the first nondimensional frequency of the FGM beam on L/h and is presented in Table 8, and the dependence of the first nondimensional frequency on L/h and the rotational speed parameter is presented in Table 9. Table 10 shows the dependence of the first six nondimensional frequencies of the FGM beam on L/h. These numerical results point out the following: When the thickness of the FGM beam is reduced, the first nondimensional frequency of the beam decreases; however, the strongest decrease of the first frequency is when the thickness decreases from L/h = 5 to L/h = 10. The frequency of the C-C beam is always greater than that of the S-S beam in the same other conditions. When increasing parameter , the first nondimensional frequency of the beam increases; this is because when increasing the porosity ratio of the material, both mass and stiffness of the structure decrease; however, the stiffness of the beam decreases sharply than the decrease of the mass. This also means that the porosity ratio is a parameter to control the vibration of the FGM beam; in other words, changing the porosity ratio will change the working performance of the structure.
5. Conclusions
By using the new shear deformation theory type of hyperbolic sine functions, finite element formulations for free vibration analysis of rotating FGM beams resting on two-parameter elastic foundations are established in this paper, in which the initial geometrical imperfection and the porosity of the material are considered. Based on the numerical results, some important novel explorations are drawn as follows: The rotational speed has a significant effect on the free vibration response as well as the vibration mode shape of the structures. When the rotational speed increases, the frequencies of the beam increase, especially the first frequency, and the vibration mode shape is no longer symmetrical through position x = L/2. As the axial compressive load value increases, the frequency decreases, and the strongest change manifests itself at the first frequency as the axial load increases. When reducing the beam thickness and the porosity ratio, the frequencies of the FGM beam are reduced.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The author declares that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the University of Transport Technology Foundation for Science and Technology Development (Grant no. 1139/QD-DHCNGTVT).