Microsoft Word - ETASR_V13_N2_pp10310-10315 Engineering, Technology & Applied Science Research Vol. 13, No. 2, 2023, 10310-10315 10310 www.etasr.com Hien et al.: Finite Element Analysis of a Continuous Sandwich Beam resting on Elastic Support and … Finite Element Analysis of a Continuous Sandwich Beam resting on Elastic Support and Subjected to Two Degree of Freedom Sprung Vehicles Ta Duy Hien University of Transport and Communications, Vietnam tdhien@utc.edu.vn Nguyen Duy Hung Campus in Ho Chi Minh City, University of Transport and Communications, Vietnam hungnd_ph@utc.edu.vn Nguyen Trong Hiep University of Transport and Communications, Vietnam nguyentronghiep@utc.edu.vn Giap Van Tan Hanoi Architectural University, Vietnam tangv@hau.edu.vn Nguyen Van Thuan Nha Trang University, Vietnam thuannv@ntu.edu.vn (corresponding author) Received: 3 November 2022 | Revised: 4 December 2022 | Accepted: 12 December 2022 ABSTRACT This paper has developed a Finite Element Method (FEM) to calculate the dynamic response of a continuous sandwich beam resting on elastic support subjected to moving vehicles. The equation of motion is derived using the classical beam theory and FEM. The vehicle model is a two Degree of Freedom (2DOF) system that moves with a constant velocity. The governing equation of motion is integrated by applying the Wilson-θ time integration method to obtain the dynamic response in each time step. Numerical examples investigate the displacement of the sandwich beam with various values of the structure and vehicle velocity. The effects of the stiffness of elastic support and the vehicle velocity on displacement are studied. Keywords-FEM; forced vibration; continuous beam; elastic support I. INTRODUCTION Structures in the construction industry can be multiform and diverse. The most commonly used structures in civil engineering are columns [1, 2], piles [3], beams [4-8], frames [9, 10], and plates [11-19]. In general, a beam is a flexural member and the normal stress has maximum values at the top and the bottom of the beam. So, sandwich beams will have optimal bearing capacity if the outer layer is a high-strength material. In practice, the fibers at the top and bottom of the beam can be made of steel and the core layer from concrete or wood. If the sandwich beam is made of a suitable material, it will have a better load capacity, reduced height, and a larger span. Thus, the sandwich beam has many advantages if it is made of appropriate materials. A sandwich beam usually has two layers, i.e. the outer layer and the core layer, as shown in Figure 1. Structures are subjected to various dynamic loads, e.g. from wind, moving vehicles, impulsive loading, etc. Structure Engineering, Technology & Applied Science Research Vol. 13, No. 2, 2023, 10310-10315 10311 www.etasr.com Hien et al.: Finite Element Analysis of a Continuous Sandwich Beam resting on Elastic Support and … dynamics have many practical aspects, e.g. the dynamic behavior of beams [20-23] and the dynamics of plates [24, 25]. The dynamic response of an elastic plate in a viscoelastic medium, resting on a viscoelastic Winkler foundation is investigated in [26]. The forced vibration of functionally graded plates resting on a viscoelastic elastic foundation was solved analytically in [27]. The effect of randomness in the elastic modulus on the eigenvalue of free vibration of non- uniform beams is studied using stochastic finite elements in [28]. The variability dynamic response of a beam subjected to a moving load with various uncertain parameters is investigated by Monte Carlo simulation in [29]. The dynamic response of a beam resting on a viscoelastic foundation was analyzed with a novel state-space formulation considering the interaction of the oscillator-beam-foundation system in [30]. The dynamic response of sandwich beams with a viscoelastic core subjected to moving loads was studied by FEM in [31]. An analytical solution was preformed to investigate the transient response of sandwich beams in [32]. A new semi-analytical method was applied in [33] to study the dynamics of beams. The transient response of double beams with viscoelastic boundary conditions under a moving load was investigated by the analytical method in [34]. The nonlinear dynamic response of the coupled vehicle-pavement system was computed using the Galerkin truncation method in [35], in which the pavement was modeled as a Timoshenko beam on the six-parameter foundation. The dynamic response of bridges to moving vehicle loads was investigated using a semi-analytical approach in [36]. The arch bridge subjected to the moving loads was modeled as a continuous beam resting on interconnected springs to compute the dynamic response using the analytical method [37]. Author in [38] presented an analytical and numerical method to study the dynamics of a beam–mass system. Experiments and finite element simulations studied the dynamic response of metallic sandwich beams under impact loading in [39]. Fig. 1. Cross section of a sandwich beam. In this study, FEM is used to study the forced vibration of a continuous sandwich beam resting on an elastic support. II. FINITE ELEMENT FORMULATION Let us consider a continuous sandwich beam resting on elastic supports, subjected to moving vehicles, with the coordinate system placed at the mid-plane of the beam as shown in Figure 2. The deflection w of the beam can be written in a matrix form as [40]:                          1 1 1 2 3 4 2 2 , e w t t w x t N x N x N x N x w t t N x q t                   (1) where    e q t denotes the displacement vector of the finite element and   N x are Hermite interpolation functions, defined as follows:         2 3 2 1 22 3 2 2 3 2 3 42 3 2 1 3 2 ; 1 2 3 2 ; e e e e e e e e x x x x N x N x x L L L L x x x x N x N x x L L L L                          (2) Fig. 2. Continuous sandwich beam resting on elastic supports and subjected to a moving vehicle. The strain energy of the beam element is: 2 2 0 2 20 2 0 2 1 2 e Face A L Face e Core A Core w E y dA x U dx w E y dA x                                        ����������� ����������� (3) The potential energy of the elastic support is defined as:   2 0 1 2 s sU k w x    (4) where ks is the stiffness of the elastic support. The kinetic energy of the beam element is given as: 2 2 0 0 0 1 2 eL e Face Core A A Face Core w w dx t t                                              ��������� ��������� (5) The governing equation of the sandwich beam is derived as: Core Face Face Engineering, Technology & Applied Science Research Vol. 13, No. 2, 2023, 10310-10315 10312 www.etasr.com Hien et al.: Finite Element Analysis of a Continuous Sandwich Beam resting on Elastic Support and …       Tb b s wM w K K w N f            ɺɺ (6) where , , s b b K K M           are the stiffness matrix of the elastic support and the stiffness and mass matrices of the sandwich beam, respectively. The equation of motion of the moving mass m1 and the contact force fw are determined as follows:    1 0m m m mm y c y w k y w    ɺɺ ɺ ɺ (7)  1 2 2 1w mf m m g m w m y   ɺɺ ɺɺ (8) where wm is the vertical dynamic deflection of the beam at mass m2. Combining (6) to (8), the governing equation is obtained:                     1 1 1 2 * 0* 0 * 0 ( ) + 0 b T m m b s T m m mx cw wM m N m c N cy ym K K k w N m m g yc x N k N k                                            ɺɺ ɺ ɺɺ ɺ ɺɺ (9) where:                 2 2 2 2 2 * * 2 ( ) * ( ) ( ) T T x T T xx x m m N N c m x t N N k m x t N N m x t N N     ɺ ɺ ɺɺ (10) III. SOLVING THE EQUATION OF MOTION WITH THE WILSON- METHOD The equation of motion of the vehicle-bridge interaction can be written in the following general form:                  M u t C u t K u t F t  ɺɺ ɺ (11) where         , , ,M C K F t are the mass, damping, and stiffness matrices, and the force vector, respectively. In order to use the Wilson scheme, we introduce the effective stiffness matrix as form:      0 1K K a M a C     ɶ (12) where:   0 12 1 ; a a tt       (13) The vector of effective forces at time t+θ can be calculated as follows:                                      1 4 5 0 1 2 F t F t F t P t C a u t a u t a u t M a u t a u t a u t               ɶ ɺ ɺ ɺɺ ɺ ɺ ɺɺ (14) where 2 4 5 1 ; 1; 2 2 t a a a t                   . The values of the parameters are selected as 1 / 2, =1 / 6  , and 1.37  . The following linear system of equations computes the displacements at time after initialization of displacement, velocity, and acceleration vectors:      K̂ u t F t       ɶ (15) The vector of acceleration at the time t+θ is calculated by:                 0 2 3 u t a u t u t a u t a u t        ɺɺ ɺ ɺɺ (16) where 3 1 1 2 a    . Acceleration, velocity, and displacement at time t t are calculated from:                                         6 7 8 9 1 u t t u t u t u t u t t u t a u t a u t t u t t u t t u t a u t a u t t                          ɺɺ ɺɺ ɺɺ ɺɺ ɺ ɺ ɺɺ ɺɺ ɺ ɺɺ ɺɺ (17) where:   2 26 7 8 9 1 1 ; ; ; 2 a t a t a t a t                   . IV. NUMERICAL EXAMPLES A. Example 1 Considering the continuous sandwich beam shown in Figure 3, subjected to a moving vehicle with 2DOF. The rectangular sandwich beam has 3 uniform spans of 4m, the dimensions of cross-section are h=60cm, b=25cm, and the thickness of the face layer is 1cm. The parameters for the vehicle in Figure 3 are: m1=1680kg and m2=840kg. The structural properties of the sandwich beam are given in Table I. TABLE I. MATERIAL PROPERTIES OF THE SANDWICH BEAM Fiber Material properties Material Young's modulus (GPa) Mass density (kg/m 3 ) Face layer Steel 200 7800 Core Wood 10 700 Figure 4 shows the displacement at the center of the beam for 2 cases of stiffness of elastic support. The speed of the vehicle is set to 5, 10, and 15m/s. The Figure clearly shows that the displacements at the center of the beam are highly similar. So, the effect of speed of the vehicle on the displacement of the beam is small. Engineering, Technology & Applied Science Research Vol. 13, No. 2, 2023, 10310-10315 10313 www.etasr.com Hien et al.: Finite Element Analysis of a Continuous Sandwich Beam resting on Elastic Support and … Fig. 3. Three-span sandwich beam subjected to a moving vehicle. (a) (b) Fig. 4. Displacement at the center of the beam for Kgoi equal to (a) 10 7 , (b) 10 8 N/m. (a) (b) Fig. 5. Displacement at the contact point with the m2 mass of the vehicle for Kgoi equal to (a) 10 7 , (b) 10 8 N/m. The displacement of the mass of the vehicle m2 is investigated with vehicle speed values of 5, 10, and 15m/s and stiffness of the elastic support Kgoi equal to 10 7 and 10 8 N/m. When the speed is slow, the displacement is close to the static displacement. At high vehicle speed, the displacement is fluctuated around the displacement line at low speed. Figure 6 illustrates the displacement at the center of the beam and at the contact point m2 of the moving vehicle with speed of 10m/s for 3 stiffening cases of elastic support. As the stiffness of the support increases, the stiffness of the structure increases, leading to a decrease in displacement. When the stiffness of the support is quite large, the displacement at the center of the beam has an almost sinusoidal shape. It is clear that the stiffness of elastic support affects the displacement of the beam. (a) (b) Fig. 6. Displacements for 3 stiffening cases of elastic support: (a) at the center of the beam, (b) at the contact point. Fig. 7. Displacements of sandwich beam resting on elastic support. B. Example 2 In this example, the dynamics of a sandwich beam is investigated by replacing two rigid supports of the example 1 L k c m 1 m 2 y wm v m m w x Kgoi L/3 L/3L/3 Kgoi Engineering, Technology & Applied Science Research Vol. 13, No. 2, 2023, 10310-10315 10314 www.etasr.com Hien et al.: Finite Element Analysis of a Continuous Sandwich Beam resting on Elastic Support and … with elastic supports as described in Figure 7. The displacements at the center of the beam and at the contact point are shown in Figure 8 for 3 stiffening cases of elastic support. Figure 8 clearly shows that the displacements are greatly increased. It also shows the strong influence of the stiffness of the elastic support on the displacements. Increased stiffness causes the displacement to decrease drastically. (a) (b) Fig. 8. Displacements for 3 stiffening cases of elastic support. (a) at the center of the beam, (b) at contact point. V. CONCLUSIONS The forced vibration of continuous sandwich beams traversed by a moving vehicle has been studied in this paper. Finite element method based on the classical beam theory was used to derive the governing equations of motion of the sandwich beams. The governing equations were solved with the Wilson method. Numerical examples were employed to investigate the effects of vehicle velocity and the stiffness of elastic support on the displacement response of the beam. The result of the numerical examples showed that changes in the stiffness of the elastic support have a significant effect on the displacement response of the beam. ACKNOWLEDGMENT This research is funded by University of Transport and Communications (UTC) under grant number T2021-CT- 007TĐ. REFERENCES [1] P. C. Nguyen, D. D. Pham, T. T. Tran, and T. Nghia-Nguyen, "Modified Numerical Modeling of Axially Loaded Concrete-Filled Steel Circular- Tube Columns," Engineering, Technology & Applied Science Research, vol. 11, no. 3, pp. 7094–7099, Jun. 2021, https://doi.org/10.48084/ etasr.4157. [2] D. T. Hang, X. T. Nguyen, and D. N. Tien, "Stochastic Buckling Analysis of Non-Uniform Columns Using Stochastic Finite Elements with Discretization Random Field by the Point Method," Engineering, Technology & Applied Science Research, vol. 12, no. 2, pp. 8458–8462, Apr. 2022, https://doi.org/10.48084/etasr.4819. [3] P. H. V. Nguyen and P. C. Nguyen, "Effects of Shaft Grouting on the Bearing Behavior of Barrette Piles: A Case Study in Ho Chi Minh City," Engineering, Technology & Applied Science Research, vol. 11, no. 5, pp. 7653–7657, Oct. 2021, https://doi.org/10.48084/etasr.4389. [4] P. B. Thang and L. V. Anh, "Structural analysis of steel-concrete composite beam bridges utilizing the shear connection model," Transport and Communications Science Journal, vol. 72, no. 7, pp. 811– 823, 2021. [5] D. X. Quy and V. T. Nga, "Static analysis of beam resting on elastic foundation by anisotropic beam-foundation element taking into account non-contact between beam and foundation," Transport and Communications Science Journal, vol. 72, no. 5, pp. 552–564, 2021. [6] P.-C. Nguyen, T. N. Van, and H. T. Duy, "Stochastic Free Vibration Analysis of Beam on Elastic Foundation with the Random Field of Young’s Modulus Using Finite Element Method and Monte Carlo Simulation," in 6th International Conference on Geotechnics, Civil Engineering and Structures, Ha Long, Vietnam, Oct. 2021, pp. 499–506, https://doi.org/10.1007/978-981-16-7160-9_50. [7] G. Qiao and S. Rahmatalla, "Dynamics of Euler-Bernoulli beams with unknown viscoelastic boundary conditions under a moving load," Journal of Sound and Vibration, vol. 491, Jan. 2021, Art. no. 115771, https://doi.org/10.1016/j.jsv.2020.115771. [8] P. Hoa and P.-C. Nguyen, "Effects of size-dependence on static and free vibration of FGP nanobeams using finite element method based on nonlocal strain gradient theory," Steel and Composite Structures, vol. 45, pp. 331–348, Nov. 2022, https://doi.org/10.12989/scs.2022.45.3.331. [9] P. C. Nguyen, "Nonlinear Inelastic Earthquake Analysis of 2D Steel Frames," Engineering, Technology & Applied Science Research, vol. 10, no. 6, pp. 6393–6398, Dec. 2020, https://doi.org/10.48084/etasr.3855. [10] M. Usarov, G. Mamatisaev, and G. Ayubov, "Forced vibrations of a box element of a multi-story building under dynamic impact," Magazine of Civil Engineering, vol. 114, no. 6, pp. 11406–11406, 2022, https://doi.org/10.34910/MCE.114.6. [11] J. G. R. Iniguez, M. L. Daza-Torres, A. P. Gonzalez, and A. Cros, "Natural frequency of a heavy flexible plate: power law evolution as a function of length," Latin American Journal of Solids and Structures, vol. 18, no. 5, Jun. 2021, Art. no. e377, https://doi.org/10.1590/1679- 78256479. [12] A. P. Yankovskii, "Analysis of Thermal Response in Reinforced Plates under a Dynamic Explosion-Type Loading," Mechanics of Composite Materials, vol. 57, no. 4, pp. 439–448, Sep. 2021, https://doi.org/ 10.1007/s11029-021-09967-w. [13] L. Zhang et al., "High-hardness polyurea coated steel plates subjected to combined loadings of shock wave and fragments.," Latin American Journal of Solids and Structures, vol. 19, no. 2, Feb. 2022, Art. no. e433, https://doi.org/10.1590/1679-78256882. [14] A. Dasdemir, "A Modal Analysis of Forced Vibration of a Piezoelectric Plate with Initial Stress by the Finite-Element Simulation," Mechanics of Composite Materials, vol. 58, no. 1, pp. 69–80, Mar. 2022, https://doi.org/10.1007/s11029-022-10012-7. [15] Q.-H. Pham and P.-C. Nguyen, "Dynamic stability analysis of porous functionally graded microplates using a refined isogeometric approach," Composite Structures, vol. 284, Mar. 2022, Art. no. 115086, https://doi.org/10.1016/j.compstruct.2021.115086. [16] V. Karpilovsky, "Finite elements for the analysis of reissne-rmindlin plates with joint interpolation of displacements and rotations," International Journal for Computational Civil and Structural Engineering, vol. 17, no. 3, pp. 48–62, Sep. 2021, https://doi.org/ 10.22337/2587-9618-2021-17-3-48-62. [17] Q.-H. Pham, P.-C. Nguyen, and T. Thanh Tran, "Dynamic response of porous functionally graded sandwich nanoplates using nonlocal higher- order isogeometric analysis," Composite Structures, vol. 290, Jun. 2022, Art. no. 115565, https://doi.org/10.1016/j.compstruct.2022.115565. [18] B. Uymaz, "Buckling Characteristics of FGM Plates Subjected to Linearly Varying In-Plane Loads," Mechanics of Composite Materials, vol. 57, no. 1, pp. 69–80, Mar. 2021, https://doi.org/10.1007/s11029- 021-09934-5. Engineering, Technology & Applied Science Research Vol. 13, No. 2, 2023, 10310-10315 10315 www.etasr.com Hien et al.: Finite Element Analysis of a Continuous Sandwich Beam resting on Elastic Support and … [19] H. D. Ta and P.-C. Nguyen, "Perturbation based stochastic isogeometric analysis for bending of functionally graded plates with the randomness of elastic modulus," Latin American Journal of Solids and Structures, vol. 17, no. 7, Sep. 2020, Art. no. e306, https://doi.org/10.1590/1679- 78256066. [20] N. V. Thuan and T. D. Hien, "Stochastic Perturbation-Based Finite Element for Free Vibration of Functionally Graded Beams with an Uncertain Elastic Modulus," Mechanics of Composite Materials, vol. 56, no. 4, pp. 485–496, Sep. 2020, https://doi.org/10.1007/s11029-020- 09897-z. [21] D. K. Nguyen, T. T. Tran, V. N. Pham, and N. A. T. Le, "Dynamic analysis of an inclined sandwich beam with bidirectional functionally graded face sheets under a moving mass," European Journal of Mechanics - A/Solids, vol. 88, Jul. 2021, Art. no. 104276, https://doi.org/10.1016/j.euromechsol.2021.104276. [22] G. Qiao and S. Rahmatalla, "Dynamics of interaction between an Euler- Bernoulli beam and a moving damped sprung mass: Effect of beam surface roughness," Structures, vol. 32, pp. 2247–2265, Aug. 2021, https://doi.org/10.1016/j.istruc.2021.04.020. [23] Q.-H. Pham, V. K. Tran, and P.-C. Nguyen, "Hygro-thermal vibration of bidirectional functionally graded porous curved beams on variable elastic foundation using generalized finite element method," Case Studies in Thermal Engineering, vol. 40, Dec. 2022, Art. no. 102478, https://doi.org/10.1016/j.csite.2022.102478. [24] A. W. de Q. R. Reis, R. B. Burgos, and M. F. F. de Oliveira, "Nonlinear Dynamic Analysis of Plates Subjected to Explosive Loads," Latin American Journal of Solids and Structures, vol. 19, no. 1, Jan. 2022, Art. no. e422, https://doi.org/10.1590/1679-78256706. [25] Q.-H. Pham, P.-C. Nguyen, V.-K. Tran, and T. Nguyen-Thoi, "Finite element analysis for functionally graded porous nano-plates resting on elastic foundation," Steel and Composite Structures, vol. 41, no. 2, pp. 149–166, 2021, https://doi.org/10.12989/scs.2021.41.2.149. [26] M. V. Shitikova and A. I. Krusser, "Force driven vibrations of nonlinear plates on a viscoelastic winkler foundation under the harmonic moving load," International Journal for Computational Civil and Structural Engineering, vol. 17, no. 4, pp. 161–180, Dec. 2021, https://doi.org/10.22337/2587-9618-2021-17-4-161-180. [27] T. D. Hien and N. N. Lam, "Vibration of functionally graded plate resting on viscoelastic elastic foundation subjected to moving loads," IOP Conference Series: Earth and Environmental Science, vol. 143, no. 1, Dec. 2018, Art. no. 012024, https://doi.org/10.1088/1755-1315/ 143/1/012024. [28] N. T. Nguyen, H. D. Ta, T. N. Van, and T. N. Dao, "Stochastic finite element analysis of the free vibration of non-uniform beams with uncertain material," Journal of Materials and Engineering Structures, vol. 9, no. 1, pp. 29–37, Apr. 2022. [29] T. D. Hien, N. D. Hung, N. T. Kien, and H. C. Noh, "The variability of dynamic responses of beams resting on elastic foundation subjected to vehicle with random system parameters," Applied Mathematical Modelling, vol. 67, pp. 676–687, Mar. 2019, https://doi.org/10.1016/ j.apm.2018.11.018. [30] G. Muscolino and A. Palmeri, "Response of beams resting on viscoelastically damped foundation to moving oscillators," International Journal of Solids and Structures, vol. 44, no. 5, pp. 1317–1336, Mar. 2007, https://doi.org/10.1016/j.ijsolstr.2006.06.013. [31] Y. Karmi, Y. Khadri, S. Tekili, A. Daouadji, and E. M. Daya, "Dynamic Analysis of Composite Sandwich Beams with a Frequency-Dependent Viscoelastic Core under the Action of a Moving Load," Mechanics of Composite Materials, vol. 56, no. 6, pp. 755–768, Jan. 2021, https://doi.org/10.1007/s11029-021-09921-w. [32] L. Li, B. Han, Q.-C. Zhang, Z.-Y. Zhao, and T. J. Lu, "Dynamic response of clamped sandwich beams: analytical modeling," Theoretical and Applied Mechanics Letters, vol. 9, no. 6, pp. 391–396, Nov. 2019, https://doi.org/10.1016/j.taml.2019.06.002. [33] Z. Dimitrovova, "Dynamic interaction and instability of two moving proximate masses on a beam on a Pasternak viscoelastic foundation," Applied Mathematical Modelling, vol. 100, pp. 192–217, Dec. 2021, https://doi.org/10.1016/j.apm.2021.07.022. [34] R. Chawla and V. Pakrashi, "Dynamic responses of a damaged double Euler–Bernoulli beam traversed by a ‘phantom’ vehicle," Structural Control and Health Monitoring, vol. 29, no. 5, 2022, Art. no. e2933, https://doi.org/10.1002/stc.2933. [35] H. Ding, Y. Yang, L.-Q. Chen, and S.-P. Yang, "Vibration of vehicle– pavement coupled system based on a Timoshenko beam on a nonlinear foundation," Journal of Sound and Vibration, vol. 333, no. 24, pp. 6623– 6636, Dec. 2014, https://doi.org/10.1016/j.jsv.2014.07.016. [36] J. L. Humar and A. M. Kashif, "Dynamic response of bridges under travelling loads," Canadian Journal of Civil Engineering, vol. 20, no. 2, pp. 287–298, Apr. 1993, https://doi.org/10.1139/l93-033. [37] P. K. Chatterjee and T. K. Datta, "Dynamic analysis of arch bridges under travelling loads," International Journal of Solids and Structures, vol. 32, no. 11, pp. 1585–1594, Jun. 1995, https://doi.org/10.1016/0020- 7683(94)00193-Z. [38] Y.-M. Wang, "The dynamic analysis of a beam-mass system due to the occurrence of two-component parametric resonance," Journal of Sound and Vibration, vol. 258, no. 5, pp. 951–967, Dec. 2002, https://doi.org/10.1006/jsvi.2002.5183. [39] L. Jing, Z. Wang, J. Ning, and L. Zhao, "The mechanical response of metallic sandwich beams under foam projectile impact loading," Latin American Journal of Solids and Structures, vol. 8, pp. 107–120, 2011, https://doi.org/10.1590/S1679-78252011000100006. [40] J. Reddy, An Introduction to the Finite Element Method, 3rd ed. New York, NY, USA: McGraw-Hill Education, 2005.