Abstract

This paper presents the first known vibration characteristic of rectangular thick plates on Pasternak foundation with arbitrary boundary conditions on the basis of the three-dimensional elasticity theory. The arbitrary boundary conditions are obtained by laying out three types of linear springs on all edges. The modified Fourier series are chosen as the basis functions of the admissible function of the thick plates to eliminate all the relevant discontinuities of the displacements and their derivatives at the edges. The exact solution is obtained based on the Rayleigh–Ritz procedure by the energy functions of the thick plate. The excellent accuracy and reliability of current solutions are demonstrated by numerical examples and comparisons with the results available in the literature. In addition, the influence of the foundation coefficients as well as the boundary restraint parameters is also analyzed, which can serve as the benchmark data for the future research technique.

1. Introduction

Rectangular plates as the common structural components have been extensively used in various engineering fields such as aerospace, military, and marine industries. Thus, the knowledge of vibration characteristics of the rectangular plates is of particular importance for the predesign of the engineering structures. In general, the analysis for the vibration problem of the rectangular plates is based on the classical thin plate theory [1] and two-dimensional approximate theories [25], that is, first-order shear deformation theory and higher-order deformation theory. However, the present results indicate that the two-dimensional plate theories have shortcomings [68]. For the classical thin plate theory, it neglects the effect of the transverse shear deformation by the simplified assumption that the normal to the undeformed mid-plane remains normal after deformation. For case of the two-dimensional approximate theories, the shear factor strongly relies on the boundary condition and they often need more hardware resources to obtain proper accuracy solution. Based on the above analysis, the three-dimensional elastic theory is presented to overcome the weakness of the two-dimensional plate theories, which does not rely on any hypotheses or deserve any numerical precision and can be used to solve the vibration problem of thick rectangular plates.

Over the past several decades, the three-dimensional vibration analyses of thick rectangular plates have been investigated by many researchers. Srinivas et al. [9] analyzed the free vibration of thick rectangular plates by the direct method based on the three-dimensional linear, small deformation theory. In their study, the boundary condition is limited to the all simply supported case on four edges. Cheung and Chakrabarti [10] used the finite layer method which is an extension of the well-known finite element method to study the free vibration of rectangular plates with various classical boundary conditions on the basis of the three-dimensional linear, small deformation theory. Hutchinson and Zillmer [11] used the series solution method to analyze the free vibration of a completely free parallelepiped. Fromme and Leissa [12] extended the Fourier method to investigate the free vibration of the rectangular parallelepiped with simple classical boundary conditions based on the three-dimensional elasticity theory. On the basis of the three-dimensional elasticity theory and differential quadrature method (DQM), Malik and Bert [13], Liew and Teo [14], and Liew et al. [15] investigated the free vibration characteristics of rectangular plates with some selected classical boundary conditions. Leissa and Zhang [16] employed the Ritz method to study the three-dimensional problem of determining the free vibration frequencies and mode shapes for a rectangular parallelepiped which is completely fixed on one face and free on other five faces. In this research, the displacements of the rectangular parallelepiped are assumed in the form of algebraic polynomials. Also, Liew et al. [1719] applied the Ritz method and the three-dimensional elastic theory to establish the continuum three-dimensional Ritz formulation to study the vibration characteristic of the thick rectangular plates with classical boundary conditions. The orthogonally generated polynomial functions are used as admissible functions to make analysis in such plates. Itakura [20] adopted the three-dimensional theory to study the free vibration of rectangular thick plate with classical boundary conditions. Lim et al. [21, 22] presented a concise formulation and an efficient method of solution to study the free vibration of thick, shear deformable plates with classical boundary conditions. Zhou et al. [23] used the Chebyshev polynomials as admissible functions and applied Ritz method to study the three-dimensional vibration of rectangular plates with classical boundary conditions.

By means of the review of the above references, most of the studies for the vibration problem of the thick rectangular plates do not consider the elastic foundations. However, a lot of engineering problems boil down to a rectangular plate on elastic foundations, such as footing of buildings, pavement of roads, and base of heavy machines. In the practical application, the Pasternak model (also referred to as the two-parameter model) [2426] is widely used to describe the mechanical behavior of the foundation, in which the well-known Winkler model [27] is a special case. In addition, the boundary condition may not always be classical case in reality, and a variety of possible boundary conditions including classical boundary conditions, elastic boundary restraints, and the combinations of two or more conditions may be encountered [68, 2847]. Based on the open published paper, we can know that merely Zhou et al. [48, 49] extended the Chebyshev-Ritz method to study the free vibration characteristics of rectangular thick plates and thick circular plates resting on elastic foundations, respectively. However, it should be noted that the boundary conditions are limited to the FFFF, SSSS, and CCCC classical boundary conditions. There is no work reported on the free vibration behavior of rectangular thick plates on Pasternak foundation with general boundary conditions. Thus, to present a new three-dimensional exact solution for the free vibrations of rectangular thick plates with arbitrary boundary conditions is of crucial importance.

In this paper, the first known three-dimensional vibration characteristic of rectangular thick plates on Pasternak foundation with arbitrary boundary conditions is studied. The theoretical formulations are based on the three-dimensional linear, small-strain elasticity theory. The modified Fourier series are chosen as the basis functions of the admissible function of the thick plates in the Ritz method. The arbitrary boundary conditions are obtained by laying out three types of liner springs on all edges. The excellent accuracy and reliability of current solutions are demonstrated by numerical examples and comparisons with the results available in the literature. In addition, the influence of the foundation coefficients as well as the boundary restraint parameters are also analyzed, which can serve as the benchmark data for the future research technique.

2. Theoretical Formulations

Consider a rectangular thick plate on Pasternak foundation with arbitrary boundary conditions as shown in Figure 1. A Cartesian coordinate is also shown in the Figure 1, which will be used in the analysis. , , and denote the displacement components in , , and directions. Length, width, and thickness of the rectangular thick plates are assumed as , , and , respectively. As shown in Figure 1, three groups of linear restraint springs are arranged at all sides of the rectangular thick plates to separately simulate the arbitrary boundary conditions. The undersurface of the plate is continuously rested on an elastic foundation represented by the Winkler/Pasternak model, in which the Winkler and shear stiffness are denoted by and , respectively, seen from Figure 1.

According to the three-dimensional elasticity theory, the strain components of rectangular thick plates can be expressed as

Based on Hooke’s law, the corresponding stress-strain relations of the rectangular thick plates can be written aswhere and are as follows:in which and are Young’s modulus and Poisson’s ratio of the rectangular thick plates. The strain energy for rectangular thick plates is given in integral form by

The kinetic energy of the rectangular thick plates is depicted as

As mentioned before, in this paper, the boundary conditions of the rectangular thick plates are the arbitrary boundary condition and obtained by laying out three types of linear springs on all the edges. Thus, the corresponding boundary conditions can be described asand the potential energy P stored in the boundary springs can be expressed as

The potential energy of the elastic foundation () is given aswhere Kw and Ks are the stiffness of the Winkler layer and shear layer, respectively.

The total energy functional for the rectangular thick plates can be expressed as

For the free vibration of the plate, the displacement components can be expressed in terms of the displacement amplitude functions:where denotes the eigenfrequency of the plate and .

Each of the displacement amplitude functions , , and is expressed, respectively, in the form of a three-dimensional (3D) Fourier cosine series supplemented with closed-form auxiliary functions introduced to eliminate all the relevant discontinuities of the displacement function and its derivatives at the edges; that is,where , , , and , , and are the Fourier coefficients of three-dimensional Fourier series expansions for the displacements functions, respectively., , , , , ,, , , , , ,, , , , , and are the supplemented coefficients of the auxiliary functions. The closed-form auxiliary functions are given as follows:

It is easy to verify that

Substituting (11) into (9) and then minimizing the functional with respect to the coefficients of the admissible functions, that isthe equations of motion for the rectangular plates can be yielded and are given in the matrix form:where and , respectively, represent the stiffness submatrix and mass submatrix. The column vectors of the unknown coefficients , , and are shown as follows:

The exact natural frequencies and mode shapes of the rectangular plates can be easily obtained by solving (15).

3. Numerical Results and Discussion

Since the natural frequencies and mode shapes are obtained from the Ritz method and the number of terms of the admissible functions being chosen as infinity is unrealistic in the actual calculations, thus, it is very important to make the convergence study of the present method to understand the convergence rate and the accuracy of the method. For general purposes in the future, the nondimensional foundation parameters and frequency parameters are used for numerical results as follows: , , and , where . Figure 2 shows the convergence studies of the first four frequency parameters for a rectangular plate with CFCF boundary condition and different truncated numbers , , and . The geometrical dimensions of the plates are used:  m,  m, and  m. The foundation coefficients are and . It is observed that the proposed method has fast convergence and good stability. In view of the excellent numerical behavior of the current solution, the truncation numbers will be simply set as in Figure 2.

Next, the comparison study will be carried out by the present method and other methods presented in the existing publications. Table 1 shows the comparison of the first seven frequency parameters for a square thick plate with SSSS, CCCC, and FFFF boundary conditions. In this example, the thick plate is without considering the foundation parameters. In order to check the present method, the results are compared with other published solutions by using the 3D Ritz method with general orthogonal polynomials [17], the 3D Ritz method with general orthogonal polynomials using the Gram–Schmidt process [18], and the 3D Ritz method with simple algebraic polynomials [20]. From the comparison, we can see a consistent agreement of the results taken from the current method and referential data.

Furthermore, the rectangular thick plate with various boundary conditions and elastic foundations is examined. Comparison is presented in Table 2. In this table, first four frequency parameters are obtained for various boundary conditions and foundation coefficients according to the present formulation and compared with those given by Zhou et al. [48] and Omurtag et al. [50] based on different numerical methods and elastic theories. It is observed that the frequencies are in excellent agreement with those obtained from FEM, which verifies the accuracy and efficiency of the proposed model.

From Figure 2 and Tables 1 and 2, we can draw the conclusion that the present approach has good convergence and excellent accuracy and reliability to study the vibration analysis of rectangular thick plates on Pasternak foundation with arbitrary boundary conditions. Thus, based on the comparison studies, the first known results of rectangular thick plates rested on Pasternak foundation with various elastic boundary conditions as well as foundation coefficients are presented in Table 3. The geometrical dimensions of the plates are the same as those in Figure 1. The boundary condition is clamped at = constants and ; and the elastic restraints on the are uniformly () changed from 104 to 1020. From the table, first we can see that the frequency parameters gradually increase as the elastic restraints increase and then they almost keep unchanged while the stiffness is larger than 1018. In order to enhance the understanding of the phenomenon, the corresponding fundamental mode shapes of Table 3 are depicted in Figure 3. From Figure 3, it is directly seen that the boundary elastic restraints have a significant effect on the vibration behavior of the rectangular thick plate.

Lastly, the influence of the foundation coefficients of the Pasternak foundation on the rectangular thick plates subjected to various boundary conditions is reported for the first time, as shown in Figure 4. The elastic boundary condition is defined as . The geometric parameters are the same as Figure 2. It can be easily obtained that, regardless of the boundary conditions of the rectangular thick plates, there exists a certain range of the elastic foundation coefficients during which the frequency parameter Ω increases and out of which the influence on frequency parameter is negligible.

4. Conclusions

In this paper, a unified method is developed for the vibration analysis of the rectangular thick plates on Pasternak foundation with arbitrary boundary conditions based on the linear, small-strain 3D elasticity theory. The method combines the advantages of the Ritz method and the modified Fourier expansion. Under the current framework, regardless of boundary conditions, each displacement function of the plate is expanded as the superposition of a standard three-dimensional cosine Fourier series and several auxiliary functions are introduced to remove any potential discontinuousness of the original rectangular plates unknown and its derivatives at the edges. The good convergence and excellent accuracy and reliability are checked and validated by the comparison with the results presented by other contributors. New results for free vibration of rectangular thick plates with elastic boundary conditions are presented, which may be used for benchmarking of researchers in the field. In addition, by means of the influences of elastic restrain parameters and foundation coefficients on free vibration characteristic of the rectangular thick plates, we can know that the frequency parameters increase rapidly as the restraint parameters and foundation coefficients increase in the certain range.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors gratefully acknowledge the financial support from the National Natural Science Foundation of China (no. 61633008 and no. 61503073) and the Natural Research Fund of Jilin Province (no. 20170101125JC).