Research Article  Open Access
Dinh Kien Nguyen, An Ninh Thi Vu, Ngoc Anh Thi Le, Vu Nam Pham, "Dynamic Behavior of a Bidirectional Functionally Graded Sandwich Beam under Nonuniform Motion of a Moving Load", Shock and Vibration, vol. 2020, Article ID 8854076, 15 pages, 2020. https://doi.org/10.1155/2020/8854076
Dynamic Behavior of a Bidirectional Functionally Graded Sandwich Beam under Nonuniform Motion of a Moving Load
Abstract
A bidirectional functionally graded Sandwich (BFGSW) beam model made from three distinct materials is proposed and its dynamic behavior due to nonuniform motion of a moving point load is investigated for the first time. The beam consists of three layers, a homogeneous core, and two functionally graded face sheets with material properties varying in both the thickness and longitudinal directions by power gradation laws. Based on the firstorder shear deformation beam theory, a finite beam element is derived and employed in computing dynamic response of the beam. The element which used the shear correction factor is simple with the stiffness and mass matrices evaluated analytically. The numerical result reveals that the material distribution plays an important role in the dynamic response of the beam, and the beam can be designed to meet the desired dynamic magnification factor by appropriately choosing the material grading indexes. A parametric study is carried out to highlight the effects of the material distribution, the beam layer thickness and aspect ratios, and the moving load speed on the dynamic characteristics. The influence of acceleration and deceleration of the moving load on the dynamic behavior of the beam is also examined and highlighted.
1. Introduction
Investigations on vibration of beams made of functionally graded materials (FGMs), a new type of composite materials initiated by Japanese researchers in 1984 [1], have been extensively carried out in the last two decades. It has been shown that the natural frequencies of FGM beams are governed by the variation of the material properties in the beam thickness or length direction [2–6]. Dynamic behavior of FGM beams under moving loads, the topic discussed herein, has also been reported in recent years. In this line of works, Şimşek and Kocatürk [7] and Şimşek [8] employed polynomials to approximate the displacement field in computing the dynamic response of FGM beams to moving loads. The authors showed that the dynamic deflection of the beams is significantly influenced by the material gradation in the beam thickness. The dynamic behavior of beams due to a moving harmonic load is also significantly influenced by the variation of the material properties in the beam length [9]. The RayleighRitz method was used in combination with the differential quadrature method by Khalili et al. [10] in examining dynamic behavior of an EulerBernoulli beam due to a moving mass. The material properties are considered to follow an exponential or a powerlaw function. The forced vibration of an EulerBernoulli beam under a moving oscillator was studied by Rajabi et al. [11] using the PetrovGalerkin method. Gan et al. [12] formulated a finite element formulation for computing dynamic response of a Timoshenko beam with material properties varying along the beam length. Linear variation of the beam width was considered by the authors. The dynamic response of an axially FGM beam to a moving harmonic force was examined by Wang and Wu [13] using the Lagrange method. The effect of axial compressive force and the temperature rise were taken into consideration by the authors. Nguyen and Bui [14] studied the effect of temperature rise on dynamic response of a FGM Timoshenko beam under a moving load using a hierarchical beam element. Finite element method was also used by Esen et al. [15] in computing the dynamic response of a FGM Timoshenko to an accelerating moving mass.
With the development of the advanced manufacturing methods [16], FGMs can be employed to fabricate structural Sandwich elements to improve performance of structures. FGM Sandwich structures with smooth variation of material properties can eliminate the interface separation problems which are often seen in the conventional Sandwich structures. Analysis of FGM Sandwich beams has drawn much attention from researchers recently. Chakraborty et al. [17] presented a finite element procedure for thermoelastic analysis of FGM and FGM Sandwich Timoshenko beams. Based on different beam models, Apetre et al. [18] examined the bending behavior of Sandwich beams with a FGM core. Adopting a higherorder sandwich panel theory, Rahmani et al. [19] studied free vibration of Sandwich beams with a syntactic core. Pradhan and Murmu [20] used the modified differential quadrature method to study thermomechanical vibration of FGM Sandwich beams resting on an elastic foundation. A refined sinusoidal shear deformation beam theory was used by Zenkour et al. [21] to study static bending of a FGM Sandwich beam resting on a Pasternak foundation. Su et al. [22] employed the general Fourier formulation to calculate frequencies of FGM Sandwich beams on an elastic foundation. Free vibration analysis of FGM Sandwich beams was carried out by Amirani et al. [23] and Yang et al. [24] using the elementfree Galerkin and meshfree radial point interpolation methods. Vo et al. [25–27] and Nguyen et al. [28–30] presented various higherorder shear deformation theories for buckling, free vibration, and bending analyses of FGSW beams. The transverse displacement in the theories is split into bending and shear parts, and the effect of the thickness stretching is taken into account. The bucking and free vibration of FGM Sandwich beams were also considered in [31] using a refined hyperbolic shear and normal deformation beam theory. Lagrange multiplier method was used in combination with Newmark method by Şimşek and Alshujairi [32] to evaluate dynamic response of FGM Sandwich beams under two moving forces. Songsuwan et al. [33] employed the Ritz and Newmark methods to study dynamic behavior of FGM Sandwich beams under a moving harmonic load.
Development of structural elements with material properties varying in two or more directions for withstanding severe general loadings is of great important in practice. Several bidirectional FGM beam models and their mechanical behavior have been considered in recent years. Şimşek [34] and Hao and Wei [35] assumed an exponential variation for the material properties in both the thickness and length directions in their vibration analysis of FGM beams. The authors showed that the dynamic characteristics of the beams are greatly influenced by the material gradient distribution, and the material properties of the bidirectional FGM beams can be designed to meet the goals of optimization. Wang et al. [36] considered the material properties varying by an exponential function along the beam length and a power law though the thickness in their free vibration study of bidirectional FGM EulerBernoulli beams. Free vibration analysis of FGM beams with material properties varying in the length and thickness directions by various gradation laws was carried out by Huynh et al. [37] using the NURBS isogeometric finite element approach. Based on a quasi3D beam theory, Karamanli [38] studied static bending of a FGM Sandwich beam with material properties varying in the length and thickness direction by the power laws. The symmetric smoothed particle hydrodynamics (SSPH) method was employed to evaluate the response of the beam to a distributed load. Finite element method was used by Nguyen et al. [39] and Nguyen and Tran [40] to study vibration of bidirectional FGM Timoshenko beams. The beams were considered to be made from four materials with volume fraction varying in both the thickness and length directions by the power gradation laws.
It is clear from the above literature review that, except for the work in [38], the mechanical behavior of bidirectional functionally graded Sandwich (BFGSW) beams has not been studied so far. This paper tries to make an effort to address this shortage. To this end, a BFGSW beam model is proposed and its dynamic behavior due to nonuniform motion of a moving point load is studied. The beam made from three distinct materials has three layers, a homogeneous core, and two FGM face sheets with material properties varying in both the thickness and length directions by power gradation laws. The conventional unidirectional FGM Sandwich beam, for example, the Sandwich beam in [32, 33], is a special case of the present beam model. In addition to the BFGSW beam model, the main novelty of this paper is that the effect of variation of the material properties in both the thickness and longitudinal directions on the dynamic behavior of FGM Sandwich beams under moving loads is taken into consideration herein for the first time. Based on the firstorder shear deformation beam theory, a finite beam element is formulated and employed in computing the dynamic response of the beam. The firstorder shear deformation beam theory requires a shear correction factor [41, 42], but it leads to a simple element. The element with stiffness and mass matrices evaluated analytically is derived using quadratic and cubic polynomials to interpolate the rotation and transverse displacement, respectively. The effects of the material distribution, the layer thickness and aspect ratios, and the moving load parameters on the dynamic behavior of the beam are investigated in detail. The influence of the acceleration and deceleration of the moving load on the dynamic response of the beam is also examined and discussed.
2. BFGSW Beam
A simply supported BFGSW beam with rectangular cross section (b × h) subjected to a load F_{0}, moving from left to right as depicted in Figure 1, is considered. The beam consists of three layers, a homogeneous core, and two FGM face layers with material properties varying in both the length and thickness directions. In the figure, the xaxis is chosen on the midplane, and z_{0} = −h/2, z_{1}, z_{2} and z_{3} = h/2 are, respectively, the vertical coordinates of the bottom layer and the interfaces of the layers and the top layer.
The beam is made from three materials, referred to as M1, M2, and M3. The core of the beam is pure M1, while the two face layers are bidirectional FGM with volume fraction of the constituents smoothly varying in both the x and zdirections according to [43, 44]:where L is the beam length; V_{1}, V_{2}, and V_{3} are, respectively, the volume fraction of M1, M2, and M3; n_{x} and n_{z} are the axial and transverse grading indexes. It can be verified that (1) defines the volume fraction of the conventional unidirectional FGM Sandwich beams in [26, 27] if nx = 0 or M2 is identical to M3. Figure 2 shows the distribution in the thickness and longitudinal directions of V_{1}, V_{2}, and V_{3} for n_{x} = n_{z} = 0.5 and z_{2} = −z_{1} = h/5.
(a)
(b)
(c)
The effective material properties, P_{f}, such as the elastic moduli and mass density evaluated by Voigt’s model are of the following form [44]:where P_{1}, P_{2}, and P_{3} are the properties of the M1, M2, and M3, respectively.
Substituting (1) into (2), one getswith
One can easily verify that if n_{x} = 0 or M2 is identical to M3, (3) returns to the effective properties of the unidirectional transverse FGM Sandwich beams made from M1 and M3. Figure 3 shows the variation of the effective Young’s modulus E_{f} and mass density ρ_{f} of the BFGSW beam made from a mixture of alumina (Al_{2}O_{3}) as M1, stainless steel (SUS304) as M2, and aluminum (Al) as M3 for n_{x} = n_{z} = 0.5 and n_{x} = n_{z} = 3. The properties of Al_{2}O_{3}, SUS304, and Al adopted from [32] are listed in Table 1.
(a)
(b)
(c)
(d)

3. Mathematical Formulation
Based on the firstorder shear deformation theory, the displacements in the x and zdirections, u_{1}(x, z, t), and u_{3}(x, z, t), respectively, are given bywhere u(x, t) and (x, t) are, respectively, the axial and transverse displacements of the point on the xaxis, θ(x, t) is the rotation of the cross section, and t is the time variable. Equation (5) leads to the axial train ε_{xx} and the shear strain γ_{xz} in the following forms:
In (6) and hereafter, a subscript comma is used to denote the derivative with respect to the variable which follows.
Constitutive equation based on linear behavior of the beam material is of the formwhere σ_{xx} and τ_{xz} are, respectively, the axial and shear stresses; E_{f} (x, z) and G_{f} (x, z) are the effective Young and shear moduli, defined by (3); is the shear correction factor, chosen by 5/6 for the beam with rectangular cross section herein.
From (6) and (7), one can write the strain energy of the beam in the formwhere A is the crosssectional area; A_{11}, A_{12}, A_{22}, and A_{33} are, respectively, the extensional, extensionalbending coupling, bending, and shear rigidities, defined as
Substituting E_{f} (x, z) and G_{f} (x, z) from (3) into (9), one can write the rigidities in the formwhere , and are, respectively, the rigidities of homogeneous and unidirectional transverse FGM beams made from M1, M2, and M3. These terms can be explicitly evaluated and their expressions are given by (A.1)–(A.4) in Appendix.
The kinetic energy resulted from (5) is of the formwhere an over dot is used to denote the derivative with respect to the time variable t; ρ_{f} (x,z) is the effective mass density; I_{11}, I_{12}, and I_{22} are the mass moments, defined as
As the rigidities, the above mass moments can also be written in the formwith , and being given by (A.5)–(A.7) in Appendix.
The beam is assumed under action of the load F_{0}, moving from the left end to the right end as shown in Figure 1. The potential of the load F_{0} is simply given bywhere is Dirac delta function; x is the abscissa, measured from the left end of the beam; s(t) is the function describing motion of the force F_{0} at the time t, which can be calculated aswhere is the initial speed of the moving load and a is the acceleration of the load which is assumed to be constant herein. The motion of the load F_{0} is defined by the sign of the acceleration a; namely, if a = 0, the load is in a uniform motion (constant moving speed), a > 0 defines an acceleration motion, and a < 0 defines a deceleration motion.
4. Finite Element Formulation
The differential equation of motion for the beam can be obtained by applying Hamilton’s principle to (8), (11), and (14). However, due to the rigidities and mass moments, as seen from (10) and (13), functions of x, a closedform solution for such equation is difficult to obtain. A finite element formulation is derived in this section for computing the dynamic response of the beam. To this end, the beam is assumed to be divided into a number of elements with length l. The vector of nodal displacements for a twonode beam element contains six components aswhere d_{u} and d_{b} denote the vectors of nodal axial and bending displacements, respectively; u_{i}, , and θ_{i} (i = 1, 2) are, respectively, the axial displacement, transverse displacement, and rotation at the node i. The order of the nodal values is not necessary as in (16), but it is convenient to separate the nodal axial displacements from the bending ones. A superscript “T” in (16) and hereafter is used to denote the transpose of a vector or a matrix.
The displacements and rotation are interpolated from their nodal values according towhere N_{u} = , = , and N_{θ} = are the matrices of interpolation functions for u, , and θ, respectively. The following polynomials are employed as the interpolation functions in the present work.
For the axial displacement u,
For the transverse displacement ,
For the rotation θ,
In (19) and (20), is the shear deformation parameter. The cubic and quadratic polynomials in (19) and (20) were firstly derived by Kosmatka in [45] for homogeneous beams, and they have been used by several authors [5, 39] to formulate the finite element formulations for analyzing FGM beams.
Using the interpolation functions (18)–(20), one can write the strain energy (8) in the formwhere NE is the total number of elements used to discretize the beam and k is the element stiffness matrix, which can be written in submatrices aswith k_{aa}, k_{ab}, k_{bb}, and k_{ss} being, respectively, the element stiffness matrices stemming from the axial stretching, axialbending coupling, bending, and shear deformation. The expressions for these submatrices are as follows:
Making use of (10), the integrals in (23) can be evaluated explicitly, and thus explicit expressions for the submatrices are obtained.
The kinetic energy (11) can also be written in the formwhere the element mass matrix m can be written in submatrices as follows:withbeing, respectively, the element mass matrices resulting from the translations in axial and transverse directions, axial translationsectional rotation coupling, and crosssectional rotation. Making use of (13), explicit expressions for the submatrices in (26) can also be obtained.
The discrete equation of motion for the dynamic analysis of the beam can be written in the form [46]where M and K are, respectively, the global mass and stiffness matrices, obtained by assembling the matrices m and k over the elements; D and are, respectively, the vectors of nodal displacements and accelerations; F_{ex} is the vector of the nodal external force with the following form:where f_{ex} is the element nodal force vector. Except for the element under the load F_{0}, the element nodal force vector f_{ex} is zero for all other elements, and the interpolation functions in (28) are evaluated at the current position of the force F_{0}. Equation (27) can be solved by the Newmark method. The average acceleration method which ensures the numerically unconditional stability [46] is adopted herein.
5. Numerical Investigation
Numerical investigation is carried out in this section to examine the effects of various parameters such as the material grading indexes, layer thickness ratio, and moving load speed on the dynamic behavior of the beam. To this end, a simply supported beam with b = 0.5 m and h = 1 m, made from Al_{2}O_{3}, SUS304, and Al with the material data in Table 1, is employed herewith. Otherwise stated, an aspect ratio L/h = 20 is assumed. A uniform increment time step ∆t = ∆T/200 with ∆T being the total time, necessary for the load completely crossing the beam, is used for the Newmark procedure. For the convenience of discussion, the following dynamic magnification factor (DMF) is introduced:where = L^{3}F_{0}/48E_{s}I is the static deflection of a fully SUS304 beam under the load F_{0}, acting at the midspan. Three numbers in brackets are used to denote the layer thickness ratio; for example, (212) means that the thickness ratio of the layers from bottom to top surfaces is 2 : 1 : 2.
5.1. Accuracy and Convergence Studies
Before computing the dynamic response of the beam, the accuracy and convergence of the derived element are firstly verified. Since there is no data on the present beam model in the literature, the verification is carried out for a unidirectional transverse FGM Sandwich beam made from Al_{2}O_{3} and Al in [32, 33], a special case of the present BFGSW beam when M2 is identical to M3.
In Table 2, the frequency parameter, with ω being the fundamental frequency, obtained in the present work is compared to the result of Şimşek and Alshujairi in [32], where a semianalytical method has been used. Very good agreement between the frequency parameter of the present work with that of [32] is noted from Table 2, regardless of the material index and the layer thickness ratio.

To verify the accuracy of the derived beam element in evaluating the dynamic response, Figure 4 compares the time histories for midspan dimensionless deflection of a transverse FGM Sandwich beam made from Al_{2}O_{3} and Al obtained in the present paper with that of Songsuwan et al. [33] by using the Ritz method. Very good agreement between the present result with that of [33] is seen from Figure 4. Note that in Figure 4 is the maximum static deflection of the simply supported beam made from pure Al.
Table 3 shows the convergence of the derived beam element in evaluating the fundamental frequency parameter of the BFGSW beam. As seen from the table, the convergence is achieved by using 26 elements, regardless of the material grading indexes and the layer thickness ratio. Note that the convergence of the results in Table 2 and also in Figure 4 has been achieved by using only 16 elements. Thus, the longitudinal variation of the material properties of the BFGSW beam makes the convergence of the element considerably slower. Because of this convergence result, 26 elements are used for all the computations reported in the following.

5.2. Uniform Motion
The dynamic response of the BFGSW beam under uniform motion of the moving load (a = 0) is considered in this subsection. In Figure 5, the time histories for midspan dimensionless deflection of the BFGSW beam are depicted for various moving load speeds, different layer thickness ratios, and two pairs of the grading indexes, n_{x} = n_{z} = 0.5 and n_{x} = n_{z} = 3. In addition to the moving load speed, the figure also shows a significant influence of the grading indexes and the layer thickness ratio on the time histories of the beam. At a given moving load speed, the maximum deflection of the beam with grading indexes n_{x} = n_{z} = 3 is considerably higher than that of the beam with n_{x} = n_{z} = 0.5. The time at which the deflection attains the maximum value is also altered by the change of the grading indexes. The influence of the layer thickness ratio is also seen from Figure 5, where the maximum deflection is seen to be lower for the beam associated with a larger core thickness. This due to the fact that the beam with a larger core thickness has a higher percentage of Al_{2}O_{3}. Since Young’s modulus of Al_{2}O_{3} is higher than that of the other constituents, the rigidities of the beam with a larger core thickness are higher, and this results in the lower deflection.
(a)
(b)
(c)
(d)
In Figure 6, the transverse grading index n_{z} versus the DMF of the (221) beam is shown for different moving load speeds and two values of the axial index, n_{x} = 0.5 and n_{x} = 3. The influence of the transverse index n_{z} on the DMF of the BFGSW beam is similar to that of the transverse FGM sandwich beam, previously investigated in [33], and the DMF increases by the increase of the transverse index n_{z}, irrespective of the moving load speed. The increase of the DMF with the increase of the n_{z} can be explained by the decrease of the effective moduli, as can be seen from (3), and this leads to the decrease of the beam rigidities. The effect of the transverse index on the DMF, as seen from Figure 6, is more significant for n_{z} in range from 0 to 5. The axial index n_{x} and the moving load speed can change the amplitude of the DMF, but they hardly alter the variation of the factor with the transverse index n_{z}. Figure 6 also shows an important role of the moving load speed on the dynamic response, and the DMF of the beam increases with the increase of the moving load speed, regardless of the material grading indexes. The effects of the moving load speed and the material grading indexes on the DMF can be seen more clearly from Figures 7 and 8, where the variation of the DMF with the moving load speed and the material grading indexes are, respectively, depicted. As seen from the figures, the DMF increases with the increase of the transverse index n_{z}, but it decreases with the increase of the axial index n_{x}. The dependence of the DMF factor upon the indexes n_{x} and n_{z} of the BFGSW beam is similar to that of the bidirectional FGM beam as reported in [39]. The core thickness ratio, as seen from Figure 8, alters the amplitude of the DMF, but it hardly changes the dependence of this factor on the grading indexes. Based on the variation of the DMF with the grading indexes depicted in Figure 8, one can design the beam with a desired DMF by approximately choosing of the material grading indexes.
(a)
(b)
(a)
(b)
(a)
(b)
In order to examine the influence of the aspect ratio on the dynamic response of the beam, Figure 9 shows the variation of the DMF with the moving load speed of the (212) and (221) beams for various values of the aspect ratio, L/h = 5, 10, 15, and 20, and for n_{x} = n_{z} = 2. The aspect ratio, as seen from the figure, has significant influence on the relation between the DMF and the moving load speed, and the DMF attains the maximum value at a higher moving load speed for the beam having a lower aspect ratio.
(a)
(b)
5.3. Acceleration and Deceleration
The influence of the acceleration and deceleration of the moving load on the dynamic behavior of the BFGSW beam is examined herewith. It is assumed that for the deceleration motion the load enters the beam with a velocity , and it exits the beam with a velocity = 0. On the other hand, for the acceleration motion, the load is assumed to enter the beam with a velocity = 0 and it exits the beam with a velocity .
In Figure 10, the time histories for the midspan dimensionless deflection of the BFGSW beam under different motions of the moving load are shown for the (212) and (221) beams with n_{x} = n_{z} = 0.5 and = 100 m/s. The deceleration (a < 0) and acceleration (a > 0) motions, as seen from the figure, have a significant influence on the time history of the beam. As seen from Figure 10, the beam under the deceleration motion executives more vibration cycles than it does under the uniform and acceleration motions of the moving load. The maximum midspan deflection is also affected by the motion type, and the maximum deflection of the beam under the acceleration motion is considerably lower than that of the beam under the uniform and deceleration motions of the moving load. The effect of the deceleration and acceleration on the dynamic response of the BFGSW beam can also be seen from Figure 11, where the moving load speed versus the DMF is shown for the (212) beam with n_{x} = n_{z} = 0.5 and n_{x} = n_{z} = 3. For the most values of the moving load speed, the DMF of the beam under the acceleration motion of the moving load is considerably lower that than of the beam under the uniform and deceleration motions of the load.
(a)
(b)
(a)
(b)
6. Conclusions
A BFGSW beam model made from three distinct materials has been proposed and its dynamic response to nonuniform motion of a moving load was studied. The beam consists of three layers, a homogeneous core, and two bidirectional FGM face sheets with material properties varying in both the thickness and length directions by the power gradation laws. The conventional transverse FGM Sandwich beam is a special case of the present BFGSW beam. Based on the firstorder shear deformation beam theory, a finite beam element was formulated and employed in combination with Newmark method to compute the dynamic response of the beam. The element based on the quadratic and cubic interpolation for the rotation and transverse displacement is simple with the stiffness and mass matrices evaluated analytically. The accuracy of the formulated element has been confirmed through a comparison study. The obtained numerical results reveal that the material distribution and the layer thickness ratio play an important role on the dynamic behavior of the beam, and the beam can be designed to meet the desired DMF by choosing suitable grading indexes and a thickness ratio. It has been shown that the DMF increases with the increase of the transverse grading index, but it decreases with increasing the axial grading index. It was also shown that the dynamic characteristics, including the time histories for the midspan deflection and the DMF, are significantly influenced by the acceleration and deceleration of the moving load. The effects of the moving load speed, the material grading indexes, the layer thickness, and aspect ratios on the dynamic behavior of the beam have been examined in detail and highlighted. It is necessary to note that though the numerical investigation has been demonstrated for the simply supported beam only, the finite element formulation derived in the present work can be used to study dynamic response of BFGSW beam with other boundary conditions as well.
Appendix
(1)The rigidities in (10) are as follows:with .(2)Mass moments in (13) are as follows:with .
Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the Vietnam National Foundation for Science and Technology Development (NAFOSTED) under Grant no. 107.022018.23.
References
 M. Koizumi, “FGM activities in Japan,” Composites Part B: Engineering, vol. 28, no. 12, pp. 1–4, 1997. View at: Publisher Site  Google Scholar
 S. A. Sina, H. M. Navazi, and H. Haddadpour, “An analytical method for free vibration analysis of functionally graded beams,” Materials and Design, vol. 30, no. 3, pp. 741–747, 2009. View at: Publisher Site  Google Scholar
 A. Mahi, E. A. Adda Bedia, A. Tounsi, and I. Mechab, “An analytical method for temperaturedependent free vibration analysis of functionally graded beams with general boundary conditions,” Composite Structures, vol. 92, no. 8, pp. 1877–1887, 2010. View at: Publisher Site  Google Scholar
 A. E. Alshorbagy, M. A. Eltaher, and F. F. Mahmoud, “Free vibration characteristics of a functionally graded beam by finite element method,” Applied Mathematical Modelling, vol. 35, no. 1, pp. 412–425, 2011. View at: Publisher Site  Google Scholar
 A. Shahba, R. Attarnejad, M. T. Marvi, and S. Hajilar, “Free vibration and stability analysis of axially functionally graded tapered Timoshenko beams with classical and nonclassical boundary conditions,” Composites Part B: Engineering, vol. 42, no. 4, pp. 801–808, 2011. View at: Publisher Site  Google Scholar
 N. Wattanasakulpong, B. Gangadhara Prusty, D. W. Kelly, and M. Hoffman, “Free vibration analysis of layered functionally graded beams with experimental validation,” Materials & Design (1980–2015), vol. 36, pp. 182–190, 2012. View at: Publisher Site  Google Scholar
 M. Şimşek and T. Kocatürk, “Free and forced vibration of a functionally graded beam subjected to a concentrated moving harmonic load,” Composite Structures, vol. 90, no. 4, pp. 465–473, 2009. View at: Publisher Site  Google Scholar
 M. Şimşek, “Vibration analysis of a functionally graded beam under a moving mass by using different beam theories,” Composite Structures, vol. 92, no. 4, pp. 904–917, 2010. View at: Publisher Site  Google Scholar
 M. Şimşek, T. Kocatürk, and Ş. D. Akbas, “Dynamic behavior of an axially functionally graded beam under action of a moving harmonic load,” Composite Structures, vol. 94, no. 8, pp. 2358–2364, 2012. View at: Publisher Site  Google Scholar
 S. M. R. Khalili, A. A. Jafari, and S. A. Eftekhari, “A mixed RitzDQ method for forced vibration of functionally graded beams carrying moving loads,” Composite Structures, vol. 92, no. 10, pp. 2497–2511, 2010. View at: Publisher Site  Google Scholar
 K. Rajabi, M. H. Kargarnovin, and M. Gharini, “Dynamic analysis of a functionally graded simply supported EulerBernoulli beam subjected to a moving oscillator,” Acta Mechanica, vol. 224, no. 2, pp. 425–446, 2013. View at: Publisher Site  Google Scholar
 B. S. Gan, T.H. Trinh, T.H. Le, and D.K. Nguyen, “Dynamic response of nonuniform Timoshenko beams made of axially FGM subjected to multiple moving point loads,” Structural Engineering and Mechanics, vol. 53, no. 5, pp. 981–995, 2015. View at: Publisher Site  Google Scholar
 Y. Wang and D. Wu, “Thermal effect on the dynamic response of axially functionally graded beam subjected to a moving harmonic load,” Acta Astronautica, vol. 127, pp. 171–181, 2016. View at: Publisher Site  Google Scholar
 D. K. Nguyen and V. T. Bui, “Dynamic analysis of functionally graded Timoshenko beams in thermal environment using a higherorder hierarchical beam element,” Mathematical Problems in Engineering, vol. 2017, Article ID 7025750, 12 pages, 2017. View at: Publisher Site  Google Scholar
 I. Esen, M. A. Koc, and Y. C ay, “Finite element formulation and analysis of a functionally graded Timoshenko beam subjected to an accelerating mass including inertial effects of the mass,” Latin American Journal of Solids and Structures, vol. 15, no. 10, p. e119, 2018. View at: Publisher Site  Google Scholar
 Y. Fukui, “Fundamental investigation of functionally gradient material manufacturing system using centrifugal force,” JSME International Journal. Ser. 3, Vibration, Control Engineering, Engineering for Industry, vol. 34, no. 1, pp. 144–148, 1991. View at: Publisher Site  Google Scholar
 A. Chakraborty, S. Gopalakrishnan, and J. N. Reddy, “A new beam finite element for the analysis of functionally graded materials,” International Journal of Mechanical Sciences, vol. 45, no. 3, pp. 519–539, 2003. View at: Publisher Site  Google Scholar
 N. A. Apetre, B. V. Sankar, and D. R. Ambur, “Analytical modeling of sandwich beams with functionally graded core,” Journal of Sandwich Structures & Materials, vol. 10, no. 1, pp. 53–74, 2008. View at: Publisher Site  Google Scholar
 O. Rahmani, S. M. R. Khalili, K. Malekzadeh, and H. Hadavinia, “Free vibration analysis of sandwich structures with a flexible functionally graded syntactic core,” Composite Structures, vol. 91, no. 2, pp. 229–235, 2009. View at: Publisher Site  Google Scholar
 S. C. Pradhan and T. Murmu, “Thermomechanical vibration of an FGM sandwich beam under variable elastic foundations using differential quadrature method,” Journal of Sound and Vibration, vol. 321, no. 12, pp. 342–362, 2009. View at: Publisher Site  Google Scholar
 A. M. Zenkour, M. N. M. Allam, and M. Sobhy, “Bending analysis of FG viscoelastic sandwich beams with elastic cores resting on Pasternak’s elastic foundations,” Acta Mechanica, vol. 212, no. 34, pp. 233–252, 2010. View at: Publisher Site  Google Scholar
 Z. Su, G. Jin, Y. Wang, and X. Ye, “A general Fourier formulation for vibration analysis of functionally graded sandwich beams with arbitrary boundary condition and resting on elastic foundations,” Acta Mechanica, vol. 227, no. 5, pp. 1493–1514, 2016. View at: Publisher Site  Google Scholar
 M. C. Amirani, S. M. R Khalili, and N. Nemati, “Free vibration analysis of sandwich beam with FG core using the element free Galerkin method,” Composite Structures, vol. 90, no. 3, pp. 373–379, 2009. View at: Publisher Site  Google Scholar
 Y. Yang, C. C. Lam, K. P. Kou, and V. P. Iu, “Free vibration analysis of the functionally graded sandwich beams by a meshfree boundarydomain integral equation method,” Composite Structures, vol. 117, pp. 32–39, 2014. View at: Publisher Site  Google Scholar
 T. P. Vo, H.T. Thai, T.K. Nguyen, A. Maheri, and J. Lee, “Finite element model for vibration and buckling of functionally graded sandwich beams based on a refined shear deformation theory,” Engineering Structures, vol. 64, pp. 12–22, 2014. View at: Publisher Site  Google Scholar
 T. P. Vo, H.T. Thai, T.K. Nguyen, F. Inam, and J. Lee, “A quasi3D theory for vibration and buckling of functionally graded sandwich beams,” Composite Structures, vol. 119, pp. 1–12, 2015. View at: Publisher Site  Google Scholar
 T. P. Vo, H.T. Thai, T.K. Nguyen, F. Inam, and J. Lee, “Static behaviour of functionally graded sandwich beams using a quasi3D theory,” Composites Part B: Engineering, vol. 68, pp. 59–74, 2015. View at: Publisher Site  Google Scholar
 T.K. Nguyen, T. TruongPhong Nguyen, T. P. Vo, and H.T. Thai, “Vibration and buckling analysis of functionally graded sandwich beams by a new higherorder shear deformation theory,” Composites Part B: Engineering, vol. 76, pp. 273–285, 2015. View at: Publisher Site  Google Scholar
 T.K. Nguyen and B.D. Nguyen, “A new higherorder shear deformation theory for static, buckling and free vibration analysis of functionally graded sandwich beams,” Journal of Sandwich Structures & Materials, vol. 17, no. 6, pp. 613–631, 2015. View at: Publisher Site  Google Scholar
 T.K. Nguyen, T. P. Vo, B.D. Nguyen, and J. Lee, “An analytical solution for buckling and vibration analysis of functionally graded sandwich beams using a quasi3D shear deformation theory,” Composite Structures, vol. 156, pp. 238–252, 2016. View at: Publisher Site  Google Scholar
 R. Bennai, H. A. Atmane, and A. Tounsi, “A new higherorder shear and normal deformation theory for functionally graded sandwich beams,” Steel and Composite Structures, vol. 19, no. 3, pp. 521–546, 2015. View at: Publisher Site  Google Scholar
 M. Şimşek and M. Alshujairi, “Static, free and forced vibration of functionally graded (FG) sandwich beams excited by two successive moving harmonic loads,” Composites Part B, vol. 108, pp. 18–34, 2017. View at: Google Scholar
 W. Songsuwan, M. Pimsarn, and N. Wattanasakulpong, “Dynamic responses of functionally graded sandwich beams resting on elastic foundation under harmonic moving loads,” International Journal of Structural Stability and Dynamics, vol. 18, no. 9, Article ID 1850112, 2018. View at: Publisher Site  Google Scholar
 M. Şimşek, “Bidirectional functionally graded materials (BFGMs) for free and forced vibration of Timoshenko beams with various boundary conditions,” Composite Structures, vol. 133, pp. 969–997, 2015. View at: Publisher Site  Google Scholar
 D. Hao and C. Wei, “Dynamic characteristics analysis ofbidirectional functionally graded Timoshenko beams,” Composite Structures, vol. 141, pp. 253–263, 2016. View at: Publisher Site  Google Scholar
 Z.h. Wang, X.h. Wang, G.d. Xu, S. Cheng, and T. Zeng, “Free vibration of twodirectional functionally graded beams,” Composite Structures, vol. 135, pp. 191–198, 2016. View at: Publisher Site  Google Scholar
 T. A. Huynh, X. Q. Lieu, and J. Lee, “NURBSbased modeling of bidirectional functionally graded Timoshenko beams for free vibration problem,” Composite Structures, vol. 160, pp. 1178–1190, 2017. View at: Publisher Site  Google Scholar
 A. Karamanli, “Bending behaviour of two directional functionally graded sandwich beams by using a quasi3d shear deformation theory,” Composite Structures, vol. 174, pp. 70–86, 2017. View at: Google Scholar
 D. K. Nguyen, Q. H. Nguyen, T. T. Tran, and V. T. Bui, “Vibration of bidimensional functionally graded Timoshenko beams excited by a moving load,” Acta Mechanica, vol. 228, no. 1, pp. 141–155, 2017. View at: Publisher Site  Google Scholar
 D. K. Nguyen and T. T. Tran, “Free vibration of tapered BFGM beams using an efficient shear deformable finite element model,” Steel and Composite Structures, vol. 29, no. 3, pp. 363–377, 2018. View at: Google Scholar
 A. A. Bousahla, F. Bourada, S. R. Mahmoud et al., “Buckling and dynamic behavior of the simply supported CNTRC beams using an integralfirst shear deformation theory,” Computers and Concrete, vol. 25, no. 2, pp. 155–166, 2020. View at: Google Scholar
 H. Matouk, A. A. Bousahla, H. Heireche et al., “Investigation on hygrothermal vibration of PFG and symmetric SFG nanobeam using integral Timoshenko beam theory,” Advances in Nano Research, vol. 8, no. 4, pp. 293–305, 2020. View at: Google Scholar
 M. NematAlla, K. I. E. Ahmed, and I. HassabAllah, “Elasticplastic analysis of twodimensional functionally graded materials under thermal loading,” International Journal of Solids and Structures, vol. 46, no. 1415, pp. 2774–2786, 2009. View at: Publisher Site  Google Scholar
 T. V. Do, D. K. Nguyen, N. D. Duc, D. H. Doan, and T. Q. Bui, “Analysis of bidirectional functionally graded plates by FEM and a new thirdorder shear deformation plate theory,” ThinWalled Structures, vol. 119, pp. 687–699, 2017. View at: Google Scholar
 J. B. Kosmatka, “An improved twonode finite element for stability and natural frequencies of axialloaded Timoshenko beams,” Computers & Structures, vol. 57, no. 1, pp. 141–149, 1995. View at: Publisher Site  Google Scholar
 R. D. Cook, D. S. Malkus, M. E. Plesha, and R. J. Witt, Concepts and Applications of Finite Element Analysis, John Wiley & Sons, Inc., New York, NY, USA, 4th edition, 2002.
Copyright
Copyright © 2020 Dinh Kien Nguyen et al. 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.