Microsoft Word - ETASR_V13_N1_pp9985-9990 Engineering, Technology & Applied Science Research Vol. 13, No. 1, 2023, 9985-9990 9985 www.etasr.com Hien et al.: Stochastic Higher-order Finite Element Model for the Free Vibration of a Continuous Beam … Stochastic Higher-order Finite Element Model for the Free Vibration of a Continuous Beam resting on Elastic Support with Uncertain Elastic Modulus Ta Duy Hien University of Transport and Communications, Vietnam tdhien@utc.edu.vn Nguyen Dang Diem Campus in Ho Chi Minh City, University of Transport and Communications, Vietnam diemnd_ph@utc.edu.vn Giap Van Tan Hanoi Architectural University, Vietnam tangv@hau.edu.vn Vu Van Hiep University of Transport and Communications, Vietnam vvhiep@utc.edu.vn Nguyen Van Thuan Nha Trang University, Vietnam thuannv@ntu.edu.vn (corresponding author) Received: 28 October 2022 | Revised: 3 December 2022 | Accepted: 4 December 2022 ABSTRACT This paper deals with a continuous beam resting on elastic support with elastic modulus derived from a random process. Governing equations of the stochastic higher-order finite element method of the free vibration of the continuous beam were derived from Hamilton's principle. The random process of elastic modulus was discretized by averaging random variables in each element. A solution for the stochastic eigenvalue problem for the free vibration of the continuous beam was obtained by using the perturbation technique, in conjunction with the finite element method. Spectral representation was used to generate a random process and employ the Monte Carlo simulation. A good agreement was obtained between the results of the first-order perturbation technique and the Monte Carlo simulation. Keywords-SFEM; free vibration; continuous beam; random field I. INTRODUCTION Beams and frames are used in many engineering applications, such as civil engineering [1-3], bridge engineering, and aerospace engineering. The importance of dynamic problems has been the subject of many research studies [4-9]. In civil engineering, the structures resting on foundations are very important [10-15]. Considering continuous beams resting on elastic foundations, elastic support is a common problem. Many practical problems are idealized to continuous beams resting on elastic foundations, elastic support such as rails in the railway, and strip foundations on the soil. The nonlinear dynamic of plates on a viscoelastic Winkler foundation under harmonic moving load is solved using the multiple time scales method [16]. Subway dynamic problem under loads from the ground and underground transport is modeled using finite element software in [17]. The dynamic of sandwich beams with a viscoelastic core subjected to a moving Engineering, Technology & Applied Science Research Vol. 13, No. 1, 2023, 9985-9990 9986 www.etasr.com Hien et al.: Stochastic Higher-order Finite Element Model for the Free Vibration of a Continuous Beam … load is investigated using the finite element method in [18]. The natural frequencies of non-uniform axially functionally graded beams were investigated using the Chebyshev collocation method in [19]. For stochastic problems, material properties, loads, and geometrical dimensions are assumed to be stochastic. Stochastic dynamics problems have been the subject of many studies [20-24]. There are several types of Stochastic Finite Element Methods (SFEM), e.g. the probabilistic finite element method [25, 26], the Spectral SFEMs (SSFEMs) using Karhunen–Loève expansion series [27], the homogeneous chaos expansion method [28] for the representation of random fields, and the SFEMs that use weighed integration techniques [29, 30]. The effect of randomness in elastic modulus on the stochastic free vibration of non-uniform beams is investigated by using an SFEM in [31]. Authors in [32] dealt with the stochastic dynamics of an infinite double beam resting on a random elastic foundation subjected to a moving load. Authors in [33] investigated a beam on a random foundation using the SFEM. The stochastic buckling of non-uniform columns was analyzed using stochastic finite elements in [34]. In the present work, the SFEM is used in the study of the free vibration of a continuous beam resting on an elastic support. The results are compared with the ones of the Monte Carlo simulation. II. STOCHASTIC FINITE ELEMENT FORMULATION FOR CONTINUOUS BEAM RESTING ON ELASTIC SUPPORT Consider a continuous beam resting on elastic support with length L as shown in Figure 1. Fig. 1. Continuous beam resting on elastic support. The displacement fields using Reddy's beam theory [35] present as formulation are:           3 0 2 0 4 , 3 , w u x z z x z x xh w x z w x                (1) The strain formulation is calculated from (1) as follows: 2 23 3 2 2 0 2 2 2 0 2 4 4 3 3 4 4 1 1 x xz w x x w h x z z z h h z z h                                 (2) The displacement fields are approximated by an interpolation functions as follows:     0 1 1 2 1 3 2 4 2 1 1 2 2 L L w x N w N N w N x N N              (3) The interpolation functions are defined as follows: 2 3 2 3 1 32 3 2 3 2 2 2 42 2 1 2 1 3 2 3 2 1 2 1 1 L L x x x x N N L L L L x x x x N x N x L L L L x x N N L L                                  (4) The displacement vector of finite element is:    1 1 1 2 2 2 T e q w w    (5) The strain energy of beam element is given as: 2 2 3 0 2 2 20 0 2 1 4 2 3 1 4 1 2 e x L e A w E y y x h x x U dx w G h x                                              (6) The potential energy of elastic support is defined as:   2 0 1 2 s sU k w x    (7) One dimensional homogeneous random field of elastic modulus is assumed:    0 1 EE x E r x    (8) where  Er x is one dimensional homogeneous random field with zero mean. The autocorrelation function of the random field is:   2 expR d           (9) where σ, d, are the coefficient of variation and the correlation distance for the random field of an elastic modulus, and the relative distance vector ξ is defined as j i  x x . The randomness in the elastic modulus of the beam element is illustrated in Figure 2.  0 1i iE E r  (10) The random field of elastic modulus in the element is approximated by averaging random variables ri at n points in the element: 1 2 0 ... 1 e nr r rE E n          (11) Engineering, Technology & Applied Science Research Vol. 13, No. 1, 2023, 9985-9990 9987 www.etasr.com Hien et al.: Stochastic Higher-order Finite Element Model for the Free Vibration of a Continuous Beam … Fig. 2. Beam element with randomness in elastic modulus. The stochastic stiffness matrix of the element is:           1 2 0 ... 1 e n K B D B d r r r B D B d n                     (12) The stiffness matrix of the structures consists of the stiffness matrix of the beam element and the stiffness of the elastic foundation:       supbeam port K K K  (13) The free vibration equation is:        2 0K M X  (14) The eigenvalue of the free vibration is 2  . The stiffness matrix, the eigenvector, and the eigenvalue are expanded to Taylor series as:                 2 0 1 1 1 2 0 1 1 1 2 0 1 1 1 1 2 1 2 1 2 Nr Nr Nr ei ei ej i i jei ei ej Nr Nr Nr ei ei ej i i jei ei ej Nr Nr Nr ei ei ej i i jei ei ej K K K K r r r r r r r r r r r r X X X X r r r r r r                                               ⋯ ⋯ ⋯ (15) and by substituting the series expressions in (6) into (5):                 2 0 1 1 1 2 0 1 1 1 2 0 1 1 1 1 2 1 2 1 0 2 Nr Nr Nr ei ei ej i i jei ei ej Nr Nr Nr ei ei ej i i jei ei ej Nr Nr Nr ei ei ej i i jei ei ej K K K r r r r r r r r r M r r r X X X r r r r r r                                                             ⋯ ⋯ ⋯ (16) Equation (16) is solved by comparing the terms of the random variables to obtain the solution:      0 00 0K - M X    (17) and    0 0 T ei ei K X X r r     (18) The first-order approximation of the variance of the eigenvalues is:           1 2 1 1 , , , Nr Nr Nr E ij i j Ei Ej Var r r r E E E R r r                         ⋯ (19) The variability of response COV of the eigenvalue is defined as:     Var COV E    (20) III. NUMERICAL EXAMPLES In this section, a continuous beam on 5 elastic supports, as shown in Figure 1, is investigated. The geometrical and material parameters of the beam are: length L=10m, height h=0.6m, mean of Young’s modulus E=30×10 5 MPa, coefficient of variation of a random field of elastic modulus σ=0.1, and mass density ρ=2400kg/m 3 . To valid the proposed approach, it is compared with Monte Carlo simulation. The random field of stiffness of the elastic modulus is generated using the spectral representation method [36, 37]:   1 0 2 cos( ) 2 ( ) , 0,1, 2,..., 1 N E n n n n n ff n u n r x A x A S N n n N                      (21) where u denotes the upper cut-off frequency beyond with the power spectral density function ( )ff nS  of the random field of the stiffness of the elastic modulus. The sample value of random field  Er x is substituted to (11) and the stiffness matrix in (12) is deterministic. The free vibration equation (14) is repeatedly solved 10000 times to get the eigenvalue. Figure 3 shows the first and second mode shapes with samples of 3 cases of stiffness of elastic support 6 7 8 10 , 10 , 10 N/msK  . It clearly shows that the first mode represents a symmetric mode shape and the second mode represents an antisymmetric mode shape. The mode shapes on Figure 3 clearly show the effect of the stiffness of support on mode shapes. The effect of the correlation distance d of the random field on the variability of the eigenvalue is shown in Figure 4, where the results of the proposed formulation were compared with those of the Monte Carlo simulation using 10,000 samples with two cases of stiffness of elastic support, i.e. 10 7 and 10 8 N/m. As seen in Figure 4, the COV varies depending on the correlation distance and the coefficient of stiffness of the elastic support. In all other cases, the obtained COV has a small value E(z) E1 E2 Ei En E Engineering, Technology & Applied Science Research Vol. 13, No. 1, 2023, 9985-9990 9988 www.etasr.com Hien et al.: Stochastic Higher-order Finite Element Model for the Free Vibration of a Continuous Beam … of the assumed standard deviation of the random field of elastic modulus. Also, the result on Figure 4 shows the good agreement between the proposed approach and Monte Carlo simulation. (a) (b) Fig. 3. Mode shapes. (a) First mode, (b) second mode. (a) (b) Fig. 4. Coefficient of variation of the eigenvalue. (a) First mode, (b) second mode. Figure 5 illustrates the statistical probabilistic distribution of the eigenvalues of the first 2 modes obtained from Monte Carlo simulation corresponding with the stiffness of elastic support Ks =10 7 N/m and two cases of correlation distance d=0.01 and d=0.1. On Figure 5, it is shown clearly that with smaller correlation distance, the histogram is quite scattered, whereas when the correlation distance is larger, the histogram is close to the normal distribution. (a) (b) (c) (d) Fig. 5. Probabilistic distribution of the eigenvalue. (a) first mode with d=0.01, (b) first mode with d=0.1, (c) second mode with d=0.01, (d) second mode with d=0.1. IV. CONCLUSIONS In this paper, the stochastic finite element method was employed to investigate the response variability of the eigenvalue of the free vibration of a continuous beam on elastic support. The free vibration equation used higher order beam theory to establish finite element formulas. The first-order perturbation solution exhibits a good agreement with Monte Carlo simulation. The response variability of eigenvalues, given in terms of the coefficient of variation, is smaller of the input standard deviation of the random process of the elastic modulus of the beam. Engineering, Technology & Applied Science Research Vol. 13, No. 1, 2023, 9985-9990 9989 www.etasr.com Hien et al.: Stochastic Higher-order Finite Element Model for the Free Vibration of a Continuous Beam … ACKNOWLEDGMENT This research is funded by University of Transport and Communications (UTC) under the grant number T2021-CT- 007TĐ. REFERENCES [1] 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. [2] 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. [3] P. C. Nguyen, B. Le-Van, and S. D. T. V. Thanh, "Nonlinear Inelastic Analysis of 2D Steel Frames : An Improvement of the Plastic Hinge Method," Engineering, Technology & Applied Science Research, vol. 10, no. 4, pp. 5974–5978, Aug. 2020, https://doi.org/10.48084/etasr. 3600. [4] D. X. Quy and N. V. Thi, "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, Jun. 2021, https://doi.org/10.47869/tcsj.72.5.4. [5] 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. [6] 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. [7] V. T. A. Ninh, "Fundamental frequencies of bidirectional functionally graded sandwich beams partially supported by foundation using different beam theories," Transport and Communications Science Journal, vol. 72, no. 4, pp. 452–467, 2021. [8] P. M. Phuc, "Using phase field and third-order shear deformation theory to study the effect of cracks on free vibration of rectangular plates with varying thickness," Transport and Communications Science Journal, vol. 71, no. 7, pp. 853–867, 2020. [9] Q.-H. Pham, V. K. Tran, T. T. Tran, P.-C. Nguyen, and P. Malekzadeh, "Dynamic instability of magnetically embedded functionally graded porous nanobeams using the strain gradient theory," Alexandria Engineering Journal, vol. 61, no. 12, pp. 10025–10044, Dec. 2022, https://doi.org/10.1016/j.aej.2022.03.007. [10] Q.-H. Pham, T. Thanh Tran, V. Ke Tran, P.-C. Nguyen, and T. Nguyen- Thoi, "Free vibration of functionally graded porous non-uniform thickness annular-nanoplates resting on elastic foundation using ES- MITC3 element," Alexandria Engineering Journal, vol. 61, no. 3, pp. 1788–1802, Mar. 2022, https://doi.org/10.1016/j.aej.2021.06.082. [11] Q.-H. Pham, V. K. Tran, T. T. Tran, T. Nguyen-Thoi, P.-C. Nguyen, and V. D. Pham, "A nonlocal quasi-3D theory for thermal free vibration analysis of functionally graded material nanoplates resting on elastic foundation," Case Studies in Thermal Engineering, vol. 26, Aug. 2021, Art. no. 101170, https://doi.org/10.1016/j.csite.2021.101170. [12] T. T. Tran, P.-C. Nguyen, and Q.-H. Pham, "Vibration analysis of FGM plates in thermal environment resting on elastic foundation using ES- MITC3 element and prediction of ANN," Case Studies in Thermal Engineering, vol. 24, Apr. 2021, Art. no. 100852, https://doi.org/ 10.1016/j.csite.2021.100852. [13] 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. [14] P.-C. Nguyen, Q. H. Pham, T. T. Tran, and T. Nguyen-Thoi, "Effects of partially supported elastic foundation on free vibration of FGP plates using ES-MITC3 elements," Ain Shams Engineering Journal, vol. 13, no. 3, May 2022, Art. no. 101615, https://doi.org/10.1016/ j.asej.2021.10.010. [15] 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, p. 115565, Jun. 2022, https://doi.org/10.1016/j.compstruct.2022.115565. [16] 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. [17] M. Barabash and P. Bogdan, "Modeling of the subway dynamic influance on the ground structure," International Journal for Computational Civil and Structural Engineering, vol. 17, no. 3, pp. 14– 23, Sep. 2021, https://doi.org/10.22337/2587-9618-2021-17-3-14-23. [18] 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. [19] W.-R. Chen, "Vibration Analysis of Axially Functionally Graded Timoshenko Beams with Non-uniform Cross-section," Latin American Journal of Solids and Structures, vol. 18, no. 7, Oct. 2021, Art. no. e397, https://doi.org/10.1590/1679-78256434. [20] T. D. Hien and P.-C. Nguyen, "Evaluation of Response Variability of Euler–Bernoulli Beam Resting on Foundation Due to Randomness in Elastic Modulus," in International Conference on Sustainable Civil Engineering and Architecture, Ho Chi Minh, Vietnam, Oct. 2019, pp. 1087–1092, https://doi.org/10.1007/978-981-15-5144-4_105. [21] 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. [22] N. V. Thuan and T. D. Hien, "Variablitity in frequencies of vehicle vibration anlysis with muiltiple random variables," Transport and Communications Science Journal, vol. 72, no. 2, pp. 215–226, 2021. [23] T.-T. Tran and D. Kim, "Uncertainty quantification for nonlinear seismic analysis of cabinet facility in nuclear power plants," Nuclear Engineering and Design, vol. 355, Dec. 2019, Art. no. 110309, https://doi.org/10.1016/j.nucengdes.2019.110309. [24] T.-T. Tran, K. Salman, S.-R. Han, and D. Kim, "Probabilistic Models for Uncertainty Quantification of Soil Properties on Site Response Analysis," ASCE-ASME Journal of Risk and Uncertainty in Engineering Systems, Part A: Civil Engineering, vol. 6, no. 3, Sep. 2020, Art. no. 04020030, https://doi.org/10.1061/AJRUA6.0001079. [25] W. K. Liu, T. Belytschko, and A. Mani, "Random field finite elements," International Journal for Numerical Methods in Engineering, vol. 23, no. 10, pp. 1831–1845, 1986, https://doi.org/10.1002/nme.1620231004. [26] 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. [27] R. G. Ghanem and P. D. Spanos, "Spectral Stochastic Finite‐Element Formulation for Reliability Analysis," Journal of Engineering Mechanics, vol. 117, no. 10, pp. 2351–2372, Oct. 1991, https://doi.org/ 10.1061/(ASCE)0733-9399(1991)117:10(2351). [28] O. H. Galal, W. El-Tahan, M. A. El-Tawil, and A. A. Mahmoud, "Spectral SFEM analysis of structures with stochastic parameters under stochastic excitation," Structural engineering and mechanics : An international journal, vol. 28, no. 3, pp. 281–294, 2008. [29] T. D. Hien, B. T. Thanh, N. N. Long, N. Van Thuan, and D. T. Hang, "Investigation Into The Response Variability of A Higher-Order Beam Resting on A Foundation Using A Stochastic Finite Element Method," in 5th International Conference on Geotechnics, Civil Engineering Works and Structures, 2020, pp. 117–122, https://doi.org/10.1007/978-981-15- 0802-8_15. Engineering, Technology & Applied Science Research Vol. 13, No. 1, 2023, 9985-9990 9990 www.etasr.com Hien et al.: Stochastic Higher-order Finite Element Model for the Free Vibration of a Continuous Beam … [30] C.-K. Choi and H.-C. Noh, "Weighted Integral SFEM Including Higher Order Terms," Journal of Engineering Mechanics, vol. 126, no. 8, pp. 859–866, Aug. 2000, https://doi.org/10.1061/(ASCE)0733-9399(2000) 126:8(859). [31] 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. [32] S. Mohammadzadeh, M. Esmaeili, and M. Mehrali, "Dynamic response of double beam rested on stochastic foundation under harmonic moving load," International Journal for Numerical and Analytical Methods in Geomechanics, vol. 38, no. 6, pp. 572–592, 2014, https://doi.org/ 10.1002/nag.2227. [33] T.-P. Chang, "Dynamic finite element analysis of a beam on random foundation," Computers & Structures, vol. 48, no. 4, pp. 583–589, Aug. 1993, https://doi.org/10.1016/0045-7949(93)90251-8. [34] 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. [35] J. Reddy, An Introduction to the Finite Element Method, 3rd ed. New York, NY, USA: McGraw-Hill Education, 2005. [36] M. Shinozuka and G. Deodatis, "Simulation of Stochastic Processes by Spectral Representation," Applied Mechanics Reviews, vol. 44, no. 4, pp. 191–204, Apr. 1991, https://doi.org/10.1115/1.3119501. [37] 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.