Shock and Vibration
Volume 2020 (2020), Article ID 8854076, 15 pages
https://doi.org/10.1155/2020/8854076

Research Article

Dynamic Behavior of a Bidirectional Functionally Graded Sandwich Beam under Nonuniform Motion of a Moving Load

Dinh Kien Nguyen,1,2 An Ninh Thi Vu,³ Ngoc Anh Thi Le,2,4 and Vu Nam Pham2,5

¹Institute of Mechanics, VAST, 18 Hoang Quoc Viet, Hanoi, Vietnam
²Graduate University of Science and Technology, VAST, 18 Hoang Quoc Viet, Hanoi, Vietnam
³University of Transport and Communications, 3 Cau Giay, Dong Da, Hanoi, Vietnam
⁴Institute of Applied Mechanics and Informatics, 291 Dien Bien Phu, Ho Chi Minh City, Vietnam
⁵Thuyloi University, 175 Tay Son, Dong Da, Hanoi, Vietnam

Correspondence should be addressed to Dinh Kien Nguyen; ndkien@imech.vast.vn

Received 13 April 2020; Revised 11 June 2020; Accepted 16 June 2020; Published 14 July 2020

Academic Editor: Marco Alfano

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.

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 first-order 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 Rayleigh-Ritz method was used in combination with the differential quadrature method by Khalili et al. [10] in examining dynamic behavior of an Euler-Bernoulli beam due to a moving mass. The material properties are considered to follow an exponential or a power-law function. The forced vibration of an Euler-Bernoulli beam under a moving oscillator was studied by Rajabi et al. [11] using the Petrov-Galerkin 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 higher-order 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 element-free Galerkin and mesh-free radial point interpolation methods. Vo et al. [25–27] and Nguyen et al. [28–30] presented various higher-order 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 Al-shujairi [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 Euler-Bernoulli 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 quasi-3D 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 first-order shear deformation beam theory, a finite beam element is formulated and employed in computing the dynamic response of the beam. The first-order 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₀, 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 x-axis is chosen on the midplane, and z₀ = −h/2, z₁, z₂ and z₃ = h/2 are, respectively, the vertical coordinates of the bottom layer and the interfaces of the layers and the top layer.

Figure 1: A simply supported BFGSW beam under a moving load.

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 z-directions according to [43, 44]:

where L is the beam length; V₁, V₂, and V₃ are, respectively, the volume fraction of M1, M2, and M3; nx and nz 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₁, V₂, and V₃ for nx = nz = 0.5 and z₂ = −z₁ = h/5.

(a)

(b)

(c)

Figure 2: Distribution of V₁, V₂, and V₃ for nx = nz = 0.5 and z₂ = −z₁ = h/5.

The effective material properties, Pf, such as the elastic moduli and mass density evaluated by Voigt’s model are of the following form [44]:

where P₁, P₂, and P₃ are the properties of the M1, M2, and M3, respectively.

Substituting (1) into (2), one gets

with

One can easily verify that if nx = 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 Ef and mass density ρf of the BFGSW beam made from a mixture of alumina (Al₂O₃) as M1, stainless steel (SUS304) as M2, and aluminum (Al) as M3 for nx = nz = 0.5 and nx = nz = 3. The properties of Al₂O₃, SUS304, and Al adopted from [32] are listed in Table 1.

(a)

(b)

(c)

(d)

Figure 3: Variation of Ef and ρf in the thickness and length directions. (a) Ef with nx = nz = 0.5. (b) Ef with nx = nz = 3. (c) ρf with nx = nz = 0.5. (d) ρf with nx = nz = 3.

Table 1: Properties of constituent materials of BFGSW beam.

Materials
Role
E (GPa)
ρ (kg/m³)
υ

Alumina (Al₂O₃)
M1
380
3960
0.3

Steel (SUS304)
M2
210
7800
0.3

Aluminum (Al)
M3
70
2702
0.3

3. Mathematical Formulation

Based on the first-order shear deformation theory, the displacements in the x- and z-directions, u₁(x, z, t), and u₃(x, z, t), respectively, are given by

where u(x, t) and (x, t) are, respectively, the axial and transverse displacements of the point on the x-axis, θ(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 form

where σxx and τxz are, respectively, the axial and shear stresses; Ef (x, z) and Gf (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 form

where A is the cross-sectional area; A₁₁, A₁₂, A₂₂, and A₃₃ are, respectively, the extensional, extensional-bending coupling, bending, and shear rigidities, defined as

Substituting Ef (x, z) and Gf (x, z) from (3) into (9), one can write the rigidities in the form

where , 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 form

where 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₁₁, I₁₂, and I₂₂ are the mass moments, defined as

As the rigidities, the above mass moments can also be written in the form

with , and being given by (A.5)–(A.7) in Appendix.

The beam is assumed under action of the load F₀, moving from the left end to the right end as shown in Figure 1. The potential of the load F₀ is simply given by

where 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₀ at the time t, which can be calculated as

where 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₀ 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 closed-form 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 two-node beam element contains six components as

where du and db denote the vectors of nodal axial and bending displacements, respectively; ui, , 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 to

where Nu = ,  = , 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 form

where 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 as

with kaa, kab, kbb, and kss being, respectively, the element stiffness matrices stemming from the axial stretching, axial-bending 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 form

where the element mass matrix m can be written in submatrices as follows:

with

being, respectively, the element mass matrices resulting from the translations in axial and transverse directions, axial translation-sectional rotation coupling, and cross-sectional 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; Fex is the vector of the nodal external force with the following form:

where fex is the element nodal force vector. Except for the element under the load F₀, the element nodal force vector fex is zero for all other elements, and the interpolation functions in (28) are evaluated at the current position of the force F₀. 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₂O₃, 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³F₀/48EsI is the static deflection of a fully SUS304 beam under the load F₀, acting at the midspan. Three numbers in brackets are used to denote the layer thickness ratio; for example, (2-1-2) 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₂O₃ 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 Al-shujairi 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.

Table 2: Comparison of frequency parameter of unidirectional FGM Sandwich beam.

nz
Source
L/h = 5
L/h = 20

(1-2-1)
(1-8-1)
(1-1-1)
(1-2-1)
(1-8-1)
(1-1-1)

0.5
Reference [32]
4.4711
4.8396
4.3192
4.6972
5.1064
4.5316

Present work
4.4715
4.84
4.3196
4.6973
5.1065
4.5318

1
Reference [32]
4.1002
4.6762
3.8627
4.2881
4.9230
4.0318

Present work
4.1006
4.6766
3.8632
4.2884
4.9231
4.0321

2
Reference [32]
3.7223
4.4510
3.4055
3.8760
4.7378
3.5379

Present work
3.7228
4.5107
3.4061
3.8764
4.738
3.5383

5
Reference [32]
3.3652
4.3456
3.0039
3.4912
4.5550
3.1100

Present work
3.3658
4.346
3.0045
3.4917
4.5552
3.1106

10
Reference [32]
3.2233
4.2727
2.8659
3.3397
4.4745
2.9651

Present work
3.2239
4.2731
2.8666
3.3402
4.4747
2.9656

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₂O₃ 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.

Figure 4: Comparison of time histories for midspan dimensionless deflection of unidirectional FGM Sandwich beam (L/h = 10, nz = 0.5,  = 50 m/s).

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.

Table 3: Convergence of the element in evaluating frequencies of BFGSW beam (L/h = 20).

nx
nz
NE = 16
NE = 18
NE = 20
NE = 22
NE = 24
NE = 26

(2-1-2)
0.5
0.3
4.7097
4.7096
4.7096
4.7096
4.7096
4.7096

1
3.8637
3.8637
3.8636
3.8635
3.8634
3.8634

3
3.2115
3.2113
3.2112
3.2111
3.2110
3.2110

1
0.3
4.8366
4.8366
4.8366
4.8366
4.8365
4.8365

1
3.8574
3.8573
3.8573
3.8572
3.8571
3.8571

3
3.2521
3.2520
3.2519
3.2518
3.2517
3.2517

(2-2-1)
0.5
0.3
4.8576
4.8575
4.8575
4.8575
4.8575
4.8575

1
4.1704
4.1703
4.1702
4.1702
4.1701
4.1701

3
3.6098
3.6097
3.6096
3.6095
3.6094
3.6094

1
0.3
3.6066
3.6066
3.6065
3.6095
3.6095
3.6095

1
4.1599
4.15998
4.1598
4.1597
4.1597
4.1597

3
3.6292
3.6290
3.6290
3.6289
3.6288
3.6288

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, nx = nz = 0.5 and nx = nz = 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 nx = nz = 3 is considerably higher than that of the beam with nx = nz = 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₂O₃. Since Young’s modulus of Al₂O₃ 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)

Figure 5: Time histories for midspan dimensionless deflection for various moving load speeds. (a) (4-1-4), nx = nz = 0.5. (b) (2-2-1), nx = nz = 0.5. (c) (4-1-4), nx = nz = 3. (d) (2-2-1), nx = nz = 3.

In Figure 6, the transverse grading index nz versus the DMF of the (2-2-1) beam is shown for different moving load speeds and two values of the axial index, nx = 0.5 and nx = 3. The influence of the transverse index nz 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 nz, irrespective of the moving load speed. The increase of the DMF with the increase of the nz 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 nz in range from 0 to 5. The axial index nx 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 nz. 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 nz, but it decreases with the increase of the axial index nx. The dependence of the DMF factor upon the indexes nx and nz 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)

Figure 6: Grading index nz versus DMF of (2-2-1) beam for various moving load speeds. (a) nx = 0.5. (b) nx = 3.

(a)

(b)

Figure 7: Variation of DMF with moving load speed of (2-2-1) beam. (a) nx = 0.5, nz is variable. (b) nz = 0.5, nx is variable.

(a)

(b)

Figure 8: Variation of DMF with material grading indexes of BFGSW beam for  = 50 m/s. (a) (2-1-2) beam. (b) (2-2-1) beam.

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 (2-1-2) and (2-2-1) beams for various values of the aspect ratio, L/h = 5, 10, 15, and 20, and for nx = nz = 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)

Figure 9: Variation of DMF with moving load speed for nx = nz = 2 and various aspect ratios. (a) (2-1-2) beam. (b) (2-2-1) beam.

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 (2-1-2) and (2-2-1) beams with nx = nz = 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 (2-1-2) beam with nx = nz = 0.5 and nx = nz = 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)

Figure 10: Time histories for midspan dimensionless deflection of BFGSW beam under different motions for nx = nz = 0.5 and  = 100 m/s. (a) (2-1-2) beam. (b) (2-2-1) beam.

(a)

(b)

Figure 11: Moving load speed versus DMF of (2-1-2) beam under different motions of moving load. (a) nx = nz = 0.5; (b) nx = nz = 3.

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 first-order 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.02-2018.23.

References

M. Koizumi, “FGM activities in Japan,” Composites Part B: Engineering, vol. 28, no. 1-2, pp. 1–4, 1997.

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.

A. Mahi, E. A. Adda Bedia, A. Tounsi, and I. Mechab, “An analytical method for temperature-dependent free vibration analysis of functionally graded beams with general boundary conditions,” Composite Structures, vol. 92, no. 8, pp. 1877–1887, 2010.

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.

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 non-classical boundary conditions,” Composites Part B: Engineering, vol. 42, no. 4, pp. 801–808, 2011.

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.

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.

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.

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.

S. M. R. Khalili, A. A. Jafari, and S. A. Eftekhari, “A mixed Ritz-DQ method for forced vibration of functionally graded beams carrying moving loads,” Composite Structures, vol. 92, no. 10, pp. 2497–2511, 2010.

K. Rajabi, M. H. Kargarnovin, and M. Gharini, “Dynamic analysis of a functionally graded simply supported Euler-Bernoulli beam subjected to a moving oscillator,” Acta Mechanica, vol. 224, no. 2, pp. 425–446, 2013.

B. S. Gan, T.-H. Trinh, T.-H. Le, and D.-K. Nguyen, “Dynamic response of non-uniform Timoshenko beams made of axially FGM subjected to multiple moving point loads,” Structural Engineering and Mechanics, vol. 53, no. 5, pp. 981–995, 2015.

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.

D. K. Nguyen and V. T. Bui, “Dynamic analysis of functionally graded Timoshenko beams in thermal environment using a higher-order hierarchical beam element,” Mathematical Problems in Engineering, vol. 2017, Article ID 7025750, 12 pages, 2017.

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.

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.

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.

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.

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.

S. C. Pradhan and T. Murmu, “Thermo-mechanical vibration of an FGM sandwich beam under variable elastic foundations using differential quadrature method,” Journal of Sound and Vibration, vol. 321, no. 1-2, pp. 342–362, 2009.

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. 3-4, pp. 233–252, 2010.

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.

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.

Y. Yang, C. C. Lam, K. P. Kou, and V. P. Iu, “Free vibration analysis of the functionally graded sandwich beams by a meshfree boundary-domain integral equation method,” Composite Structures, vol. 117, pp. 32–39, 2014.

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.

T. P. Vo, H.-T. Thai, T.-K. Nguyen, F. Inam, and J. Lee, “A quasi-3D theory for vibration and buckling of functionally graded sandwich beams,” Composite Structures, vol. 119, pp. 1–12, 2015.

T. P. Vo, H.-T. Thai, T.-K. Nguyen, F. Inam, and J. Lee, “Static behaviour of functionally graded sandwich beams using a quasi-3D theory,” Composites Part B: Engineering, vol. 68, pp. 59–74, 2015.

T.-K. Nguyen, T. Truong-Phong Nguyen, T. P. Vo, and H.-T. Thai, “Vibration and buckling analysis of functionally graded sandwich beams by a new higher-order shear deformation theory,” Composites Part B: Engineering, vol. 76, pp. 273–285, 2015.

T.-K. Nguyen and B.-D. Nguyen, “A new higher-order 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.

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 quasi-3D shear deformation theory,” Composite Structures, vol. 156, pp. 238–252, 2016.

R. Bennai, H. A. Atmane, and A. Tounsi, “A new higher-order shear and normal deformation theory for functionally graded sandwich beams,” Steel and Composite Structures, vol. 19, no. 3, pp. 521–546, 2015.

M. Şimşek and M. Al-shujairi, “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.

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.

M. Şimşek, “Bi-directional functionally graded materials (BFGMs) for free and forced vibration of Timoshenko beams with various boundary conditions,” Composite Structures, vol. 133, pp. 969–997, 2015.

D. Hao and C. Wei, “Dynamic characteristics analysis ofbi-directional functionally graded Timoshenko beams,” Composite Structures, vol. 141, pp. 253–263, 2016.

Z.-h. Wang, X.-h. Wang, G.-d. Xu, S. Cheng, and T. Zeng, “Free vibration of two-directional functionally graded beams,” Composite Structures, vol. 135, pp. 191–198, 2016.

T. A. Huynh, X. Q. Lieu, and J. Lee, “NURBS-based modeling of bidirectional functionally graded Timoshenko beams for free vibration problem,” Composite Structures, vol. 160, pp. 1178–1190, 2017.

A. Karamanli, “Bending behaviour of two directional functionally graded sandwich beams by using a quasi-3d shear deformation theory,” Composite Structures, vol. 174, pp. 70–86, 2017.

D. K. Nguyen, Q. H. Nguyen, T. T. Tran, and V. T. Bui, “Vibration of bi-dimensional functionally graded Timoshenko beams excited by a moving load,” Acta Mechanica, vol. 228, no. 1, pp. 141–155, 2017.

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.

A. A. Bousahla, F. Bourada, S. R. Mahmoud et al., “Buckling and dynamic behavior of the simply supported CNT-RC beams using an integral-first shear deformation theory,” Computers and Concrete, vol. 25, no. 2, pp. 155–166, 2020.

H. Matouk, A. A. Bousahla, H. Heireche et al., “Investigation on hygro-thermal vibration of P-FG and symmetric S-FG nanobeam using integral Timoshenko beam theory,” Advances in Nano Research, vol. 8, no. 4, pp. 293–305, 2020.

M. Nemat-Alla, K. I. E. Ahmed, and I. Hassab-Allah, “Elastic-plastic analysis of twodimensional functionally graded materials under thermal loading,” International Journal of Solids and Structures, vol. 46, no. 14-15, pp. 2774–2786, 2009.

T. V. Do, D. K. Nguyen, N. D. Duc, D. H. Doan, and T. Q. Bui, “Analysis of bi-directional functionally graded plates by FEM and a new third-order shear deformation plate theory,” Thin-Walled Structures, vol. 119, pp. 687–699, 2017.

J. B. Kosmatka, “An improved two-node finite element for stability and natural frequencies of axial-loaded Timoshenko beams,” Computers & Structures, vol. 57, no. 1, pp. 141–149, 1995.

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.