Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2020 / Article
Special Issue

Structure, Dynamics, and Applications of Complex Networks in Software Engineering

View this Special Issue

Research Article | Open Access

Volume 2020 |Article ID 8278743 | https://doi.org/10.1155/2020/8278743

Trung Thanh Tran, Quoc-Hoa Pham, Trung Nguyen-Thoi, "An Edge-Based Smoothed Finite Element for Free Vibration Analysis of Functionally Graded Porous (FGP) Plates on Elastic Foundation Taking into Mass (EFTIM)", Mathematical Problems in Engineering, vol. 2020, Article ID 8278743, 17 pages, 2020. https://doi.org/10.1155/2020/8278743

An Edge-Based Smoothed Finite Element for Free Vibration Analysis of Functionally Graded Porous (FGP) Plates on Elastic Foundation Taking into Mass (EFTIM)

Academic Editor: Xesús Nogueira
Received31 Dec 2019
Revised10 Mar 2020
Accepted17 Mar 2020
Published25 Apr 2020

Abstract

In this paper, free vibration analysis of the functionally graded porous (FGP) plates on the elastic foundation taking into mass (EFTIM) is presented. The fundamental equations of the FGP plate are derived using Hamilton’s principle. The mixed interpolation of the tensorial components (MITC) approach and the edge-based smoothed finite element method (ES-FEM) is employed to avoid the shear locking as well as to improve the accuracy for the triangular element. The EFTIM is a foundation model based on the two-parameter Winkler–Pasternak model but added a mass parameter of foundation. Materials of the plate are FGP with a power-law distribution and maximum porosity distributions in the forms of cosine functions. Some numerical examples are examined to demonstrate the accuracy and reliability of the proposed method in comparison with those available in the literature.

1. Introduction

The plate resting on the elastic foundation (EF) is one of the most common types of structures which have practical applications in civil and industrial constructions, especially in transportation and irrigation. In particular, the structures of beams and plates on the EF are subjected to moving loads of means of transport such as roadbeds affected by vehicles, railway tracks, and aircraft runways. In most publications, when investigating the mechanical behavior of structures on the EF, the researchers mainly use the one-parameter Winkler foundation model [1] or two-parameter Winkler–Pasternak foundation model [2, 3]. The analysis of plates resting on the Winkler–Pasternak foundation was previously addressed by several authors. For instance, Fazzolari [3] used an analytical method to consider free vibration and buckling of porous FG Sandwich beams resting on the EF with the Winkler–Pasternak foundation model. Leissa [4] presented results for the free vibration of rectangular plates. Xiang et al. [5] studied free vibration for Mindlin plates on the Winkler–Pasternak foundation using an analytical method. Omurtag et al. [6] used the finite element method (FEM) for the free vibration analysis of the Kirchhoff plates resting on the EF. Özçelikörs et al. [7] analyzed the exact solutions of bending, buckling, and vibration problems of a Levy-plate on the two-parameter foundation. Matsunaga [8] used a special higher-order plate theory (HSDT) to analyze vibration and buckling of thick plates on the EF. Ayvaz et al. [9] developed a modified Vlasov model to consider the earthquake response of thin plates on the EF. Shen et al. [10] based on the Rayleigh–Ritz method to study free and forced vibration of the Reissner–Mindlin plates resting on the EF. Liew et al. [11, 12] and Han and Liew [13] analyzed the free vibration of rectangular plates resting on the EF using a differential quadrature method. Zhou et al. [14] considered the vibration of rectangular plates on the EF using the Ritz method. Chucheepsakul and Chinnaboon [15] investigated plates by a two-parameter foundation model using a boundary element method. Civalek and Acar [16] investigated the bending of Mindlin plates on the EF by developing a singular convolution method. Ferreira et al. [17] presented bending and free vibration analyses of the FGP plates on the Winkler–Pasternak foundation by using radial basis functions. Shahsavari et al. [18] used a new quasi-3D hyperbolic theory for the free vibration analysis of the FGP plates resting on the EF. Zenkour and Radwan [19] proposed an exact analytical approach for free vibration analysis of laminated composite and Sandwich plates resting on the EF using a four-unknown plate theory. Duc et al. [20] presented the analysis of nonlinear thermal dynamic response of shear deformable functionally graded plates on the EF. Mahmoudi et al. [21] developed a refined quasi-three-dimensional shear deformation theory to analyze the functionally graded Sandwich plates resting on the two-parameter EF under thermomechanical loading. Duc et al. [22] used an analytical method to calculate static bending and free vibration of FG carbon nanotube-reinforced composite plate resting on Winkler–Pasternak foundations. Thang et al. [23] considered the effects of variable thickness on buckling and postbuckling behavior of imperfect sigmoid FGM plates on elastic medium subjected to compressive loading. Banh-Thien et al. [24] developed the isogeometric analysis for buckling analysis of nonuniform thickness nanoplates in an elastic medium.

In recent years, the FGP materials have attracted great interest from many researchers over the world due to their lightness and high strength. As a result, they are widely applied for civil engineering, aerospace structures, nuclear plants, and other applications. Kim et al. [25] investigated bending, vibration, and buckling of the FGP microplates using a modified couples stress based on the analytical method. Coskun et al. [26] presented analytical solutions to analyze the bending, vibration, and buckling of the FG microplates based on the third-order shear deformation theory (TSDT). Chen et al. [27] investigated the static bending and buckling of the FGP beams by using the Timoshenko beam theory. Rezaei and Saidi [28, 29] studied the vibration of rectangular and porous-cellular plates based on an analytical method. The vibration of the FGP shallow shells using an improve Fourier method was examined by Zhao et al. [30]. Moreover, the dynamics of the FGP doubly-curved panels and shells were also investigated in [31]. Li et al. [32] analyzed the nonlinear vibration and dynamic buckling of the Sandwich FGP plates with graphene platelet reinforcement (GPL) on the EF. For nonlinear problems, Sahmani et al. [33] used the nonlocal method to analyze nonlinear large-amplitude vibrations of the FGP micro-/nanoplates with GPL reinforce. Wu et al. [34] studied the dynamic of the FGP structures by using FEM. Thang et al. [35] investigated the elastic buckling and free vibration of porous cellular plates based on the first-order shear deformation theory (FSDT). Although the FGP materials have many different types, in this paper, the authors only use the distribution of porosity as presented in [25, 26].

In the other front of the development of numerical methods for computational mechanics, Liu et al. [36] have recently proposed an edge-based smoothed FEM (ES-FEM) using triangular elements which show some following excellent properties for the 2D solid mechanics analyses such as (1) the numerical results are often found superconvergent and very accurate; (2) the method is stable and works well for dynamic analysis; and (3) the implementation of the method is straightforward and no penalty parameter is used. The ES-FEM has been developed for n-sided polygonal elements [37], viscoelastoplastic analyses [38], 2D piezoelectric [39], primal-dual shakedown analyses [40], fluid structure interaction [41, 42], and various applications [4345]. Recently, in an effort to improve the accuracy of the plate and shell structural analyses, the classical MITC3 element [46] incorporated with the ES-FEM [36], has been proposed to give the so-called ES-MITC3 element [4751]. In the formulation of the ES-MITC3, the system stiffness matrix is employed using strains smoothed over the smoothing domains associated with the edges of the MITC3 elements. The numerical results demonstrated that the ES-MITC3 has the following great properties [47]: (1) the ES-MITC3 can eliminate transverse shear locking even with the ratio of the thickness to the length of the structures reach and (2) the ES-MITC3 has better accuracy than the existing triangular elements such as MITC3 [46], DSG3 [52], and CS-DSG3 [53] and is a good competitor with the quadrilateral element MITC4 element [54].

The objective of this research now is to further extend the ES-MITC3 method for free vibration analyses of the FGP plates resting on the EFTIM. The governing equations are derived from the FSDT and the Reissner–Mindlin plate theory due to simplicity and computational efficiency. Besides, the EFTIM is modelled based on a two-parameter Winkler–Pasternak foundation model and added in a mass parameter of foundation. The plate is made from the FGP materials with a power-law distribution () and maximum porosity distributions () in the forms of cosine functions. The accuracy and reliability of the present formulation are verified by comparing with those of other available numerical results. Moreover, the effects of some geometric parameters and material properties on the free vibration of the FGP plates are also examined in detail.

2. Functionally Graded Porous Material Plates on Elastic Foundation

Let us consider an FGP plate resting on EFTIM, as shown in Figure 1. The FGP materials with a variation of two constituents and three different distributions of porosity through-thickness are presented as follows [25, 26]:where is the maximum porosity value. A typical material property of the FGP materials can be considered as in the following power-law relations:where and are the typical material properties at the top and the bottom surfaces, respectively; and is the power-law index. The normalized distributions of porosity through the thickness are shown in Figure 2(a). As shown in Figure 2(a), the porosity distribution of Case 1 is symmetric with respect to the midplane of plates. Case 2 and Case 3 are bottom and top surface-enhanced distributions, respectively. Besides, Figures 2(b)2(d) show the distributions of a normalized typical property associated with three different cases of porosity distributions with parameters: , k = 1, 5, 10, and Pt/Pb= 10.

The EFTIM is built based on the Winkler–Pasternak foundation by adding a mass parameter of foundation:where is the displacement of FGP plate; and are, respectively, Winkler foundation stiffness and shear layer stiffness of the Pasternak foundation. In order to mention the effectiveness of the foundation mass involved in the oscillation as well as the continuous interaction of the spring with the plate, the parameter β with unit kg−1 is added. It characterizes the effective level of the foundation mass involved in vibration, which is determined based on an experimental basis and the ratio of the density of the foundation to the density of plate material which is defined as . Thus, the density of mass involved vibration with the foundation is determined . From equation (3) we see that, for the static problem, the EFTIM model and Winkler–Pasternak foundation model are the same. But, for the dynamic problems, these two models have differences, and when omitting the influence parameters of the foundation mass, the EFTIM model is equivalent to the Winkler–Pasternak foundation model. In addition, this foundation model also covers the Winkler foundation model when the influence of shear parameters and foundation mass parameters ignored. It was found that the EFTIM model closely resembles the true feature of the foundation, including the Pasternak and Winkler foundation models.

3. The First-Order Shear Deformation Theory and Weak Form of the FGP Plates

3.1. First-Order Shear Deformation Theory for FGP Plates

The displacement of the FGP plates in the present work based on the FSDT model can be expressed aswhere , , , , and are five unknown displacements of the midsurface of the plate. For a bending plate, the strain field can be expressed as follows:where the membrane strain is given as

The bending and transverse shear strains are written as

From Hooke’s law, the linear stress-strain relations of the FGP plates can be expressed aswherewhere presents for effective Young’s modulus and represents Poisson’s ratio.

3.2. Weak Form Equations

To obtain the motion equations of the FGP plates for the free vibration analysis, Hamilton’s principle is applied with the following form:where and are the strain and kinetic energies, respectively. The strain energy is expressed aswhere is the strain energyand is the strain energywhere andand , , , and can be given by

The kinetic energy in equation (11) is given bywhere is the kinetic energywhere is the displacement field and is the mass matrix defined bywhere

In equation (11), the kinetic energy of the mass of foundation is defined as

Substituting equations (12) and (18) into equation (11), the weak formulation for the free vibration of the FGP plate is finally obtained as

4. Formulation of an ES-MITC3 Method for FGP Plates

4.1. Formulation of the Finite Element Using the MITC3 Element

The middle surface of plate is discretized into finite three-node triangular elements with nodes such that and . Then, the generalized displacements at any point of element can be approximated aswhere is the number of nodes of ; is the shape function matrix; and is the nodal degrees of freedom (d.o.f) associated with the jth node of .

The membrane bending strains of MITC3 element can be expressed in matrix forms as follows:where

To eliminate the shear locking phenomenon as the thickness of the plate becomes small, the formulation of the transverse shear strains of the MITC3 element based on FSDT [36] in this study can be written as follows:wherewithwhere

Here, , , , and are pointed out in and is the area of the three-node triangular element as shown in Figure 3.

Substituting the discrete displacement field into equation (22), we obtained the discrete system equations for free vibration analysis of FGP plate resting on the EF, respectively, aswhere and are the stiffness and mass matrices, respectively.

The stiffness matrix in equation (38) can be written aswherewhere

Next, the mass matrix in equation (38) can be defined aswhere

4.2. Formulation of an ES-MITC3 Method for FGP Plates

The smoothing domains is constructed based on edges of the triangular elements such that and . An edge-based smoothing domain for the inner edge is formed by connecting two end-nodes of the edge to the centroids of adjacent triangular MITC3 elements, as shown in Figure 4.

Applying the edge-based smooth technique [36], the smoothed membrane, bending, and shear strain , over the smoothing domain can be created bywhere , , and the compatible membrane, bending, and the shear strains, respectively; is a given smoothing function that satisfies at least unity property .

In this study, we use the constant smoothing functionwhere is the area of the smoothing domain and is given bywhere is the number of the adjacent triangular elements in the smoothing domain ; and is the area of the triangular element attached to the edge

By substituting equations (47)–(49) into equations (24), (25), and (32), the approximation of the smoothed strains on the smoothing domain can be expressed as follows:where is the total number of nodes of the triangular MITC3 elements attached to edge k ( for boundary edges and for inner edges as given in Figure 4; is the nodal d.o.f associated with the smoothing domain ; and , , and are the smoothed membrane, bending, and shear strain gradient matrices, respectively, at the jth node of the elements attached to edge k computed by

The stiffness matrix of the FGP plate using the ES-MITC3 is assembled bywhere is the ES-MITC3 stiffness matrix of the smoothing domain and given bywhere

5. Accuracy of the Proposed Method

In this section, the various numerical examples are solved to verify the reliability and accuracy of the proposed method. For convenience, the stiffness factors and nondimensional frequencies of the plates are defined as the following equations:

To demonstrate the performance of numerical results, the relative frequency error is defined bywhere and are nondimensional frequencies of present method and nondimensional frequencies in [17, 18], respectively.

The results of the convergence of the first two nondimensional frequencies of the plate in the case of fully simple support (SSSS) plate and a fully clamped (CCCC) plate with h/a = 0.1, K1 = 100, K2 = 10, respectively, are shown in Figure 5. From these results, it can be seen that almost all frequencies corresponding to different cases of boundary conditions (BC) converge with 18 × 18 element mesh. For 18 × 18 mesh, we compare the first three nondimensional frequencies of a plate resting on the Winkler–Pasternak foundation with the published results as shown in Table 1. It can be seen that the present results agree well with the results of the authors using analytical methods [5, 14, 17] and are more accurate than those using the original MITC3 element and FEM [6]. In addition, from Table 2, it is obvious that the relative error of the present results compared to [18] is less than 2%. In [18], they used a new quasi-3D hyperbolic theory to investigate the free vibration of the FGP plate resting on the EF. These results are the basis to analyze the free vibration of FGP plates on the EFTIM.


PlatesK1K2Author

SSSS

h/a= 0.01
10010Ferreira et al. [17]2.65595.57188.5384
Zhou et al. [14]2.65510.035.57170.008.54060.03
Xiang et al. [5]2.65510.035.57180.008.54050.02
MITC32.66040.175.61030.708.62961.07
Present2.65900.125.59200.378.60170.74
50010Ferreira [17]3.34065.92858.7754
Zhou et al. [14]3.33980.025.92850.008.77750.02
Xiang et al. [5]3.34000.025.92870.008.77750.02
MITC33.34410.105.96490.618.86421.01
Present3.34300.075.94770.328.83700.70

SSSS

h/a= 0.1
20010Ferreira et al. [17]2.79025.34527.8255
Zhou et al. [14]2.77560.525.29540.937.72791.25
Xiang et al. [5]2.78420.225.30430.777.72871.24
MITC32.78740.105.32580.367.77190.68
Present2.78870.055.33620.177.79710.36
100010Ferreira et al. [17]3.98446.04308.3112
Zhou et al. [14]3.95660.705.97571.118.19541.39
Xiang et al. [5]3.98050.106.00780.588.22141.08
MITC33.98270.046.02660.278.26190.59
Present3.98360.026.03580.128.28560.31

CCCC

h/a= 0.015
1390.2166.83Ferreira et al. [17]8.166912.82116.842
Zhou et al. [14]8.16750.0112.8230.0216.8330.05
Omurtag et al. [6]8.13750.3612.8980.6016.9320.53
MITC38.18420.2112.9090.6917.0101.00
Present8.17290.0712.8720.4016.9390.58


(K1, K2)h/a
Present[18]Present[18]

(0, 0)0.059.0109.0200.118.4858.3701.37
0.108.8238.8180.068.3198.2031.41
0.158.5418.5160.298.0697.9501.50
0.208.1968.1510.557.7627.6411.58

(100, 0)0.059.3899.4300.439.0208.9171.16
0.109.2079.2310.268.8588.7531.20
0.158.9338.9340.018.6148.5051.28
0.208.5998.5770.268.3158.2031.37

(100, 100)0.0515.38315.4390.3616.33816.3200.11
0.1015.21315.2450.2116.17516.1480.17
0.1514.96214.9660.0315.93215.8950.23
0.2014.66414.6400.1615.63915.5950.28

Next, we consider an SSSS FGP plate (Al/Al2O3) with its material properties as follows: metal (Al) , and ceramic (Al2O3) , . Poisson’s ratio is fixed at . The FGP plate with even porosities is expressed as in [18]:where presents the porosity volume fraction. The stiffness factor and nondimensional frequencies of the plates are shown in equation (58) with and . The first nondimensional frequencies of present work compared with [18] are shown in Table 2.

6. Numerical Results and Discussions

For free vibration problems, a fully simple support (SSSS) FGP plate is considered, wherein , , , , , and . The first six nondimensional frequencies of the FGP plate with porosity distribution of Case 1 and stiffener of foundation are shown in Table 3; and the first six mode shapes are presented in Figure 6. The stiffness factors and nondimensional frequencies of FGP plate are shown in equation (58) with .



0.85831.81181.81572.80153.38983.4248

6.1. Influence of the Parameters of the EFTIM to the Free Vibration for the FGP Plate

Firstly, in order to investigate the effect of the featured-index of the mass of foundation β to free vibration of the FGP plate, the featured-index of mass is changed from 0 to 1. In Figure 7 and Table 4, it is seen that all cases of porosity distribution featured-index of the mass of foundation β significantly influence to free vibration of the FGP plate. As β increases, the mass of the plate increases, and the frequencies of the plate decrease. For all cases of porosity distribution of the FGP plate, the porosity distribution of Case 3 leads to the maximum values of frequencies of the plate, while the porosity distribution of Case 2 leads to the minimum values. It can be observed that the FGP plate with porosity distribution Case 3 is stiffer than plates with other porosity distributions.


Parameter of plateβ00.250.50.751

Case 11.10150.95750.85830.78470.7273
(K1 = 100, K2 = 10)2.30782.01561.81181.65951.5400
(SSSS)3.56513.11542.80152.56652.3822
Case 21.02250.89880.81130.74530.6931
(K1 = 100, K2 = 10)2.07741.83501.66121.52891.4239
(SSSS)3.24102.86302.59212.38592.2221
Case 31.25371.07350.95380.86690.8001
(K1 = 100, K2 = 10)2.73992.35312.09421.90541.7599
(SSSS)4.19723.60973.21522.92692.7044

Next, the influence of nondimensional parameters of foundation stiffness K1 and K2 is investigated. We change from 100 to 1000 and from 10 to 100 with respect to and . The first nondimensional frequency of the FGP plate with three cases of porosity distribution is presented in Table 5 and shown in Figure 8. As shown in these figure and table, when K1 and K2 increase, the nondimensional frequency of plate also increases. We also examine the effect of (β, K1) and (β, K2) parameters to nondimensional frequency. The numerical result is presented in Table 6 and Figure 9. Consequently, Winkler foundation stiffness K1 and shear layer stiffness of Pasternak foundation K2 make stiffness of plate become greater and the mass of the EF involved in the plate’s vibration makes reduce frequencies.


Case of Porosity distributionK2K1
1002505007501000

Case 1
(SSSS)
100.85830.89290.94770.99951.0488
250.92580.95791.00921.05801.1047
501.02851.05751.10421.14901.1921
751.12181.14851.19161.23321.2734
1001.20791.23271.27301.31201.3499

Case 2
(SSSS)
100.81130.84520.89880.94930.9973
250.87750.90890.95891.00641.0518
500.97781.00611.05151.09501.1369
751.06871.09471.13651.17691.2160
1001.15251.17651.21561.25341.2902

Case 3
(SSSS)
100.95381.00091.07471.14381.2090
251.04521.08831.15661.22111.2823
501.18191.22021.28151.34001.3960
751.30441.33921.39531.44921.5011
1001.41631.44841.50041.55071.5993


βK1 (K2 = 10)K2 (K1 = 100)

Case 1
100250500750100010255075100
01.10151.14591.21631.28281.34601.10151.18821.31991.43971.5502
0.250.95750.99611.05721.11501.17000.95751.03281.14731.25141.3475
0.50.85830.89290.94770.99951.04880.85830.92581.02851.12181.2079
0.750.78470.81630.86640.91380.95880.78470.84640.94031.02561.1043
10.72730.75660.80300.84690.88870.72730.78440.87150.95051.0235

Case 2
01.02251.06521.13271.19641.25691.02251.10591.23231.34691.4525
0.250.89880.93630.99571.05171.10490.89880.97211.08331.18401.2768
0.50.81130.84520.89880.94930.99730.81130.87750.97781.06871.1525
0.750.74530.77630.82560.87200.91610.74530.80600.89820.98171.0586
10.69310.72200.76770.81090.85190.69310.74960.83530.91290.9845

Case 3
01.25371.31561.41271.50351.58911.25371.37391.55361.71461.8617
0.251.07351.12651.20961.28741.36071.07351.17641.33031.46811.5941
0.50.95381.00091.07471.14381.20900.95381.04521.18191.30441.4163
0.750.86690.90970.97681.03961.09880.86690.95001.07421.18551.2872
10.80010.83950.90150.95951.01410.80010.87670.99141.09421.1880

6.2. Influence of the Parameters-FGP to Free Vibration of the Plate on EFTIM

Let us consider the effect of materials property to free vibration of the FGP plate. The power-law index k is changed from 0 to 10, and maximum porosity distributions Ω has the value from 0 to 1. We examine the SSSS FGP plate and fully clamped (CCCC) plate resting on EFTIM. The parameters of EFTIM are given by , , , and . The first nondimensional frequencies of plate with three cases of porosity distribution is shown in Table 7 and Figure 10.


SSSSCCCC
kk
02.557.51002.557.510

Case 1
00.61460.84600.94681.01841.07061.02751.21531.34321.45711.5444
0.250.65840.88300.97771.04691.09841.09351.25181.36941.47831.5633
0.50.71260.92381.01051.07671.12731.17321.28691.39301.49581.5780
0.750.78200.96771.04461.10701.15641.27111.31541.41021.50611.5848
10.87451.01051.07761.13571.18361.39241.32261.41301.50211.5766

Case 2
00.61460.84600.94681.01841.07061.02751.21531.34321.45711.5444
0.250.66090.86040.95631.02881.08331.09771.20111.31561.42731.5162
0.50.71620.86850.95691.02991.08711.17071.16511.26111.36751.4564
0.750.77770.86000.93571.00731.06711.21791.08681.15561.24961.3342
10.78520.79720.84200.90090.95801.03470.91060.93100.99041.0547

Case 3
00.61460.84600.94681.01841.07061.02751.21531.34321.45711.5444
0.250.66090.94131.04751.11961.17141.09771.36611.50121.61471.6994
0.50.71621.07851.19081.26211.31221.17071.58361.72431.83171.9084
0.750.77771.30661.42561.49171.53631.21791.92972.06562.14842.2018
10.78521.85891.96311.97721.98071.03472.40832.41382.35882.3246

As shown in these figures and tables, when k and change, the values of nondimensional frequency change with no rule. It is understandable because with each case of change in porosity distributions k and , the stiffness and the weight of the plate changes. From Figure 10, in the case of the CCCC plate, nondimensional frequency depending on k and value varies more complex than the case of the SSSS plate. If k and values are same, the frequency of the CCCC plate is larger than that of the SSSS plate. The results are quite reasonable because the SSSS boundary condition inherently offers more flexible boundary conditions than the CCCC boundary condition.

7. Conclusions

In this paper, new numerical results of free vibration of the FGP plate resting on EFTIM are studied. We used the ES-MITC3 to establish the fundamental equation of the FGP plate. The computed results obtained by this approach are in excellent agreement with others published. Our work has the following advantages.

The novel ES-MITC3 which computes the free vibration of the plate on EF takes into account the mass of foundation.

The numerical results obtained by ES-MITC3 show good agreement with the reference solutions and are more accurate than those obtained by the original MITC3.

The elastic foundation of Pasternak with three-parameters is developed by adding the featured-index of mass β to accurately describe the actual elastic foundation.

The mass of the elastic foundation involved in the vibration of the plate reduces the frequency of vibration, while two parameters K1 and K2 effect the stiffness of the plate.

The material parameters k, and the case of porosity distribution effect of vibration of the plate. Numerical results are useful for calculation, design, and testing of material parameters in engineering and technologies.

This study suggests some further works on the dynamic response and heat transfer problems of the FGP plate resting on EF using different plate theories.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was funded by the Vietnam National Foundation for Science and Technology Development (NAFOSTED) under Grant no. 107.02-2019.330.

References

  1. E. Winkler, Die Lehre von der Elasticitaet und Festigkeit, Prag, Dominicus, 1867.
  2. P. L. Pasternak, On a New Method of Analysis of an Elastic Foundation by Means of Two Foundation Constants, Gosudarstvennoe Izdatelstvo Literaturi po Stroitelstvu i Arkhitekture, Moscow, Russia, 1954, in Russian.
  3. F. A. Fazzolari, “Generalized exponential, polynomial and trigonometric theories for vibration and stability analysis of porous FG sandwich beams resting on elastic foundations,” Composites Part B: Engineering, vol. 136, pp. 254–271, 2018. View at: Publisher Site | Google Scholar
  4. A. W. Leissa, “The free vibration of rectangular plates,” Journal of Sound and Vibration, vol. 31, no. 3, pp. 257–293, 1973. View at: Publisher Site | Google Scholar
  5. Y. Xiang, C. M. Wang, and S. Kitipornchai, “Exact vibration solution for initially stressed mindlin plates on Pasternak foundations,” International Journal of Mechanical Sciences, vol. 36, no. 4, pp. 311–316, 1994. View at: Publisher Site | Google Scholar
  6. M. H. Omurtag, A. Özütok, A. Y. Aköz, and Y. Özçelikörs, “Free vibration analysis of Kirchhoff plates resting on elastic foundation by mixed finite element formulation based on Gâteaux differential,” International Journal for Numerical Methods in Engineering, vol. 40, no. 2, pp. 295–317, 1997. View at: Publisher Site | Google Scholar
  7. K. Y. Özçelikörs, C. M. Wang, and X. Q. He, “Canonical exact solutions for Levy-plates on two-parameter foundation using Green’s functions,” Engineering Structures, vol. 22, no. 4, pp. 364–378, 2000. View at: Publisher Site | Google Scholar
  8. H. Matsunaga, “Vibration and stability of thick plates on elastic foundations,” Journal of Engineering Mechanics, vol. 126, no. 1, pp. 27–34, 2000. View at: Publisher Site | Google Scholar
  9. Y. Ayvaz, A. Daloglu, and A. Dogangün, “Application of a modified Vlasov model to earthquake analysis of plates resting on elastic foundations,” Journal of Sound and Vibration, vol. 212, no. 3, pp. 499–509, 1998. View at: Publisher Site | Google Scholar
  10. H.-S. Shen, J. Yang, and L. Zhang, “Free and forced vibration of Reissner-mindlin plates with free edges resting on elastic foundations,” Journal of Sound and Vibration, vol. 244, no. 2, pp. 299–320, 2001. View at: Publisher Site | Google Scholar
  11. K. M. Liew, J.-B. Han, Z. M. Xiao, and H. Du, “Differential quadrature method for mindlin plates on Winkler foundations,” International Journal of Mechanical Sciences, vol. 38, no. 4, pp. 405–421, 1996. View at: Publisher Site | Google Scholar
  12. K. M. Liew and T. M. Teo, “Differential cubature method for analysis of shear deformable rectangular plates on Pasternak foundations,” International Journal of Mechanical Sciences, vol. 44, no. 6, pp. 1179–1194, 2002. View at: Publisher Site | Google Scholar
  13. J.-B. Han and K. M. Liew, “Numerical differential quadrature method for Reissner/mindlin plates on two-parameter foundations,” International Journal of Mechanical Sciences, vol. 39, no. 9, pp. 977–989, 1997. View at: Publisher Site | Google Scholar
  14. D. Zhou, Y. K. Cheung, S. H. Lo, and F. T. K. Au, “Three-dimensional vibration analysis of rectangular thick plates on pasternak foundation,” International Journal for Numerical Methods in Engineering, vol. 59, no. 10, pp. 1313–1334, 2004. View at: Publisher Site | Google Scholar
  15. S. Chucheepsakul and B. Chinnaboon, “An alternative domain/boundary element technique for analyzing plates on two-parameter elastic foundations,” Engineering Analysis with Boundary Elements, vol. 26, no. 6, pp. 547–555, 2002. View at: Publisher Site | Google Scholar
  16. Ö. Civalek and M. H. Acar, “Discrete singular convolution method for the analysis of mindlin plates on elastic foundations,” International Journal of Pressure Vessels and Piping, vol. 84, no. 9, pp. 527–535, 2007. View at: Publisher Site | Google Scholar
  17. A. J. M. Ferreira, C. M. C. Roque, A. M. A. Neves, R. M. N. Jorge, and C. M. M. Soares, “Analysis of plates on Pasternak foundations by radial basis functions,” Computational Mechanics, vol. 46, no. 6, pp. 791–803, 2010. View at: Publisher Site | Google Scholar
  18. D. Shahsavari, M. Shahsavari, Li Li, and B. Karami, “A novel quasi-3D hyperbolic theory for free vibration of FG plates with porosities resting on Winkler/Pasternak/Kerr foundation,” Aerospace Science and Technology, vol. 72, pp. 134–149, 2018. View at: Publisher Site | Google Scholar
  19. A. Zenkour and A. Radwan, “Free vibration analysis of multilayered composite and soft core sandwich plates resting on Winkler-Pasternak foundations,” Journal of Sandwich Structures & Materials, vol. 20, no. 2, pp. 169–190, 2018. View at: Publisher Site | Google Scholar
  20. N. D. Duc, D. H. Bich, and P. H. Cong, “Nonlinear thermal dynamic response of shear deformable FGM plates on elastic foundations,” Journal of Thermal Stresses, vol. 39, no. 3, pp. 278–297, 2016. View at: Publisher Site | Google Scholar
  21. A. Mahmoudi, S. Benyoucef, A. Tounsi, A. Benachour, E. A. Adda Bedia, and S. Mahmoud, “A refined quasi-3D shear deformation theory for thermo-mechanical behavior of functionally graded sandwich plates on elastic foundations,” Journal of Sandwich Structures & Materials, vol. 21, no. 6, pp. 1906–1929, 2019. View at: Publisher Site | Google Scholar
  22. N. D. Duc, J. Lee, T. Nguyen-Thoi, and P. T. Thang, “Static response and free vibration of functionally graded carbon nanotube-reinforced composite rectangular plates resting on Winkler-Pasternak elastic foundations,” Aerospace Science and Technology, vol. 68, pp. 391–402, 2017. View at: Publisher Site | Google Scholar
  23. P.-T. Thang, T. Nguyen-Thoi, and J. Lee, “Closed-form expression for nonlinear analysis of imperfect sigmoid-FGM plates with variable thickness resting on elastic medium,” Composite Structures, vol. 143, pp. 143–150, 2016. View at: Publisher Site | Google Scholar
  24. T. Banh-Thien, H. Dang-Trung, L. Le-Anh, V. Ho-Huu, and T. Nguyen-Thoi, “Buckling analysis of non-uniform thickness nanoplates in an elastic medium using the isogeometric analysis,” Composite Structures, vol. 162, pp. 182–193, 2017. View at: Publisher Site | Google Scholar
  25. J. Kim, K. K. Żur, and J. N. Reddy, “Bending, free vibration, and buckling of modified couples stress-based functionally graded porous micro-plates,” Composite Structures, vol. 209, pp. 879–888, 2019. View at: Publisher Site | Google Scholar
  26. S. Coskun, J. Kim, and H. Toutanji, “Bending, free vibration, and buckling analysis of functionally graded porous micro-plates using a general third-order plate theory,” Journal of Composites Science, vol. 3, no. 1, p. 15, 2019. View at: Publisher Site | Google Scholar
  27. D. Chen, J. Yang, and S. Kitipornchai, “Elastic buckling and static bending of shear deformable functionally graded porous beam,” Composite Structures, vol. 133, pp. 54–61, 2015. View at: Publisher Site | Google Scholar
  28. A. S. Rezaei and A. R. Saidi, “Application of Carrera Unified Formulation to study the effect of porosity on natural frequencies of thick porous-cellular plates,” Composites Part B: Engineering, vol. 91, pp. 361–370, 2016. View at: Publisher Site | Google Scholar
  29. A. S. Rezaei and A. R. Saidi, “Exact solution for free vibration of thick rectangular plates made of porous materials,” Composite Structures, vol. 134, pp. 1051–1060, 2015. View at: Publisher Site | Google Scholar
  30. J. Zhao, F. Xie, A. Wang, C. Shuai, J. Tang, and Q. Wang, “A unified solution for the vibration analysis of functionally graded porous (FGP) shallow shells with general boundary conditions,” Composites Part B: Engineering, vol. 156, pp. 406–424, 2019. View at: Publisher Site | Google Scholar
  31. J. Zhao, F. Xie, A. Wang, C. Shuai, J. Tang, and Q. Wang, “Vibration behavior of the functionally graded porous (FGP) doubly-curved panels and shells of revolution by using a semi-analytical method,” Composites Part B: Engineering, vol. 157, pp. 219–238, 2019. View at: Publisher Site | Google Scholar
  32. Q. Li, D. Wu, X. Chen, L. Liu, Y. Yu, and W. Gao, “Nonlinear vibration and dynamic buckling analyses of sandwich functionally graded porous plate with graphene platelet reinforcement resting on Winkler-Pasternak elastic foundation,” International Journal of Mechanical Sciences, vol. 148, pp. 596–610, 2018. View at: Publisher Site | Google Scholar
  33. S. Sahmani, M. M. Aghdam, and T. Rabczuk, “Nonlocal strain gradient plate model for nonlinear large-amplitude vibrations of functionally graded porous micro/nano-plates reinforced with GPLs,” Composite Structures, vol. 198, pp. 51–62, 2018. View at: Publisher Site | Google Scholar
  34. D. Wu, A. Liu, Y. Huang, Y. Huang, Y. Pi, and W. Gao, “Dynamic analysis of functionally graded porous structures through finite element analysis,” Engineering Structures, vol. 165, pp. 287–301, 2018. View at: Publisher Site | Google Scholar
  35. P. T. Thang, T. Nguyen-Thoi, D. Lee, J. Kang, and J. Lee, “Elastic buckling and free vibration analyses of porous-cellular plates with uniform and non-uniform porosity distributions,” Aerospace Science and Technology, vol. 79, pp. 278–287, 2018. View at: Publisher Site | Google Scholar
  36. G. R. Liu, T. Nguyen-Thoi, and K. Y. Lam, “An edge-based smoothed finite element method (ES-FEM) for static, free and forced vibration analyses of solids,” Journal of Sound and Vibration, vol. 320, no. 4-5, pp. 1100–1130, 2009. View at: Publisher Site | Google Scholar
  37. T. Nguyen-Thoi, G. R. Liu, and H. Nguyen-Xuan, “An n-sided polygonal edge-based smoothed finite element method (nES-FEM) for solid mechanics,” International Journal for Numerical Methods in Biomedical Engineering, vol. 27, no. 9, pp. 1446–1472, 2011. View at: Publisher Site | Google Scholar
  38. T. Nguyen-Thoi, G. R. Liu, H. C. Vu-Do, and H. Nguyen-Xuan, “An edge-based smoothed finite element method for visco-elastoplastic analyses of 2D solids using triangular mesh,” Computational Mechanics, vol. 45, no. 1, pp. 23–44, 2009. View at: Publisher Site | Google Scholar
  39. H. Nguyen-Xuan, G. R. Liu, T. Nguyen-Thoi, and C. Nguyen-Tran, “An edge-based smoothed finite element method (ES-FEM) for analysis of two-dimensional piezoelectric structures,” Smart Materials and Structures, vol. 18, Article ID 065015, p. 12, 2009. View at: Publisher Site | Google Scholar
  40. T. N. Thanh, G. R. Liu, H. Nguyen-Xuan, and T. Nguyen-Thoi, “An edge-based smoothed finite element method for primal-dual shakedown analysis of structures,” International Journal for Numerical Methods in Engineering, vol. 82, no. 7, pp. 917–938, 2010. View at: Google Scholar
  41. T. Nguyen-Thoi, P. Phung-Van, T. Rabczuk, H. Nguyen-Xuan, and C. Le-Van, “An application of the ES-FEM in solid domain for dynamic analysis of 2D fluid-solid interaction problems,” International Journal of Computational Methods, vol. 10, no. 1, Article ID 1340003, 2013. View at: Publisher Site | Google Scholar
  42. T. Nguyen-Thoi, P. Phung-Van, V. Ho-Huu, and L. Le-Anh, “An edge-based smoothed finite element method (ES-FEM) for dynamic analysis of 2D Fluid-Solid interaction problems,” KSCE Journal of Civil Engineering, vol. 19, no. 3, pp. 641–650, 2015. View at: Publisher Site | Google Scholar
  43. C. V. Le, H. Nguyen-Xuan, H. Askes, T. Rabczuk, and T. Nguyen-Thoi, “Computation of limit load using edge-based smoothed finite element method and second-order cone programming,” International Journal of Computational Methods, vol. 10, no. 1, Article ID 1340004, 2013. View at: Publisher Site | Google Scholar
  44. H. Nguyen-Xuan, L. V. Tran, T. Nguyen-Thoi, and H. C. Vu-Do, “Analysis of functionally graded plates using an edge-based smoothed finite element method,” Composite Structures, vol. 93, no. 11, pp. 3019–3039, 2011. View at: Publisher Site | Google Scholar
  45. H. H. Phan-Dao, H. Nguyen-Xuan, C. Thai-Hoang, T. Nguyen-Thoi, and T. Rabczuk, “An edge-based smoothed finite element method for analysis of laminated composite plates,” International Journal of Computational Methods, vol. 10, no. 1, Article ID 1340005, 2013. View at: Publisher Site | Google Scholar
  46. P.-S. Lee and K.-J. Bathe, “Development of MITC isotropic triangular shell finite elements,” Computers & Structures, vol. 82, no. 11-12, pp. 945–962, 2004. View at: Publisher Site | Google Scholar
  47. T. Chau-Dinh, Q. Nguyen-Duy, and H. Nguyen-Xuan, “Improvement on MITC3 plate finite element using edge-based strain smoothing enhancement for plate analysis,” Acta Mechanica, vol. 228, no. 6, pp. 2141–2163, 2017. View at: Publisher Site | Google Scholar
  48. T.-K. Nguyen, V.-H. Nguyen, T. Chau-Dinh, T. P. Vo, and H. Nguyen-Xuan, “Static and vibration analysis of isotropic and functionally graded sandwich plates using an edge-based MITC3 finite elements,” Composites Part B: Engineering, vol. 107, pp. 162–173, 2016. View at: Publisher Site | Google Scholar
  49. Q.-H. Pham, T.-V. Tran, T.-D. Pham, and D.-H. Phan, “An edge-based smoothed MITC3 (ES-MITC3) shell finite element in laminated composite shell structures analysis,” International Journal of Computational Methods, vol. 15, no. 7, Article ID 1850060, 2018. View at: Publisher Site | Google Scholar
  50. Q.-H. Pham, T.-D. Pham, V. T. Quoc, and D.-H. Phan, “Geometrically nonlinear analysis of functionally graded shells using an edge-based smoothed MITC3 (ES-MITC3) finite elements,” Engineering with Computers, vol. 33, pp. 1–14, 2019. View at: Publisher Site | Google Scholar
  51. D. Pham-Tien, H. Pham-Quoc, T. Vu-Khac, and N. Nguyen-Van, “Transient analysis of laminated composite shells using an edge-based smoothed finite element method,” in Proceedings of the International Conference on Advances in Computational Mechanics 2017, pp. 1075–1094, Springer, Berlin, Germany, 2018. View at: Publisher Site | Google Scholar
  52. K.-U. Bletzinger, M. Bischoff, and E. Ramm, “A unified approach for shear-locking-free triangular and rectangular shell finite elements,” Computers & Structures, vol. 75, no. 3, pp. 321–334, 2000. View at: Publisher Site | Google Scholar
  53. T. Nguyen-Thoi, P. Phung-Van, H. Nguyen-Xuan, and C. Thai-Hoang, “A cell-based smoothed discrete shear gap method using triangular elements for static and free vibration analyses of Reissner-Mindlin plates,” International Journal for Numerical Methods in Engineering, vol. 91, no. 7, pp. 705–741, 2012. View at: Publisher Site | Google Scholar
  54. K.-J. Bathe and E. N. Dvorkin, “A formulation of general shell elements-the use of mixed interpolation of tensorial components,” International Journal for Numerical Methods in Engineering, vol. 22, no. 3, pp. 697–722, 1986. View at: Publisher Site | Google Scholar

Copyright © 2020 Trung Thanh Tran 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.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views278
Downloads126
Citations

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.