Abstract

Many numerical methods have been developed for in-plane vibration of orthotropic rectangular plates with various boundary conditions; however, the exact results for such structures with elastic boundary conditions are very scarce. Therefore, the object of this paper is to present an accurate solution for free in-plane vibration of orthotropic rectangular plates with various boundary conditions by the method of reverberation ray matrix (MRRM) and improved golden section search (IGSS) algorithm. The boundary condition studied in this paper is defined as that a set of opposite edges is with one kind of simply supported boundary conditions, while the other set is with any kind of classical and general elastic boundary conditions or their combination. Its accuracy, reliability, and efficiency are verified by some numerical examples where the results are compared with other exact solutions in the published literature and the FEA results based on the ABAQUS software. Finally, some new accurate results for free in-plane vibration of orthotropic rectangular plates with elastic boundary conditions are examined and further can be treated as the reference data for other approximate methods or accurate solutions.

1. Introduction

The orthotropic rectangular plates, as the fundamental components, have achieved a widespread use in aerospace, military, and marine industries and various engineering fields, and their structural stability and safety performance rely on the vibration characteristics. Thus, in the past decades, a lot of studies for free vibrations of rectangular plates have been carried out. Among the existing jobs, most of the researchers focus on the transverse vibrations of plates while the in-plane vibration research is very scarce. For example, Li et al. [14] proposed the analytical symplectic superposition method to solve free vibration solutions of rectangular plates. Zhang et al. [5, 6] presented a finite integral transform solution to determine the vibration behavior of rectangular orthotropic thin plates with different boundary conditions. Ullah et al. [7, 8] extended the integral transform method to study the buckling behavior of rectangular plates. Of course, there are many other excellent achievements, and there is no detailed discussion here. In theory, the aforementioned methods for studying transverse vibration characteristics of rectangular plates also have potential to conduct the analysis of in-plane vibration characteristics. According to the current research situation, there are relatively few related literatures on in-plane vibration characteristics of rectangular plates; the importance of which, however, has been shown in the complicated plate structures and sandwich plates [916]. As a matter of fact, the publications pertaining to in-plane vibrations have been increasing. In recent years, the study on in-plane vibrations has increased greatly and a lot of research results have been published.

Bardell et al. [9] studied in-plane vibrations of isotropic rectangular plates due to the lack of the related valid results at that time and provided some useful benchmark data of the in-plane frequencies of rectangular plates with simply supported, clamped, and free boundary conditions, where the Rayleigh–Ritz method is adopted. Du et al. [12] studied the in-plane vibration of isotropic rectangular plates with elastically restrained edges by the Rayleigh–Ritz method in conjunction with the Fourier series method. Gorman et al. [17] studied in-plane vibration of rectangular plates subject to completely free boundaries by extending the superposition method. Applying the variational approximation method, Seok et al. [18] studied free in-plane vibration of a cantilevered rectangular plate. Dozio [19] studied the accurate in-plane modal properties of rectangular plates using the Ritz method, where arbitrary nonuniform elastic edge restraints were considered. Singh and Muhammad [20] investigated free in-plane vibration of isotropic nonrectangular plates subjected to classical boundary conditions using the variational method. Du et al. [21] applied a generalized Fourier method to study the free vibration of the rectangular plate with different point-supported boundary conditions, which enriched the in-plane vibration of the plate. Based on the dynamic stiffness method, Nefovska-Danilovic and Petronijevic [22] carried out the in-plane free vibration and response analysis of isotropic rectangular plates, in which the numerical results were very consistent with FEM.

In the aforementioned works, the solutions by the related solving methods like the Rayleigh–Ritz method are classified as approximate solutions. Although the methods for approximate solutions have the merits of the simple solving process and high precision, considering that the in-plane natural frequencies of plates are generally in high-frequency regions, the minor errors of results may also cause the poor structure stability and safety performance. Thus, to achieve higher precision of in-plane natural frequencies is the ultimate goal and the corresponding solution is also called the exact solution. Gorman [23] presented the exact solutions for free in-plane vibration of rectangular plates whose one set of opposite edges is simply supported and the other set both clamped or both free. Xing et al. [24, 25] obtained exact solutions for free in-plane vibrations of rectangular plates employing the direct separation of variables, where SS1-SS2-SS1-SS2, SS1-SS1-SS1-SS1, SS2-SS2-SS2-SS2, SS1-C-SS1-C, SS1-F-SS1-F, SS1-SS1-SS1-C, SS1-SS2-SS1-C, SS1-SS1-SS1-F, SS1-SS2-SS1-F, SS1-C-SS1-F, and other possible classical boundary conditions were considered. Later, Xing’s [26] group obtained the exact solutions of the free in-plane vibration of orthotropic rectangular plates with SS1-C-SS1-C, SS2-F-SS2-F, SS1-C-SS2-F, SS2-SS1-SS2-C, SS2-SS2-SS2-F, and SS2-SS1-SS1-F boundary conditions. The symbols S, C, F, and E represent simply supported boundary, clamped boundary, free boundary, and elastically restrained boundary conditions, respectively. Two forms of simply supported boundary conditions are as follows [23]: SS1—displacement parallel to the edge and normal stress (Nxx) perpendicular to the edge are forbidden; SS2—shear stress (Nxy) along the edge and displacement normal to the edge are forbidden. Take boundary x = 0 as an example, the symbol SS1 indicates that Nxx= 0 and and SS2 means that u= 0 and Nxy= 0.

Through the above review, it can be known that the exact solutions for in-plane vibration analysis of orthotropic rectangular plates are deficient. By investigating the existing literature, it is found that researchers have been pursuing more accurate in-plane vibration results because it can better promote the development of numerical algorithms. As far as authors know, only “comprehensive exact solutions for free in-plane vibrations of orthotropic rectangular plates” [26] presents the exact solutions for the titled problem. However, in this paper, only the classical boundary conditions are considered. In the existing engineering applications, there are still a large number of elastic boundaries except for the classical boundary conditions. At present, numerical methods can solve many problems well, but researchers have been pursuing the exact solution. Thus, it is necessary and urgent to establish an accurate analytical model to conduct in-plane vibration analysis of orthotropic rectangular plates with elastic boundary conditions. Before that, the members of the author’s team have done some innovative work around the method of reverberation ray matrix (MRRM) in solving the exact solutions of structural vibration, such as the exact solutions of isotropic rectangular plates with elastic constraints [27]. However, it should be pointed out that the mechanical properties of orthotropic rectangular plates are superior to those of isotropic rectangular plates due to the difference of material properties. Therefore, it is necessary to study the in-plane exact solutions of orthotropic rectangular plates. This paper can be regarded as an extension of previous work.

In this paper, the method of reverberation ray matrix (MRRM) and the improved golden section search (IGSS) algorithm are adopted to solve accurate solutions for free in-plane vibration of orthotropic rectangular plates. MRRM is used to obtain the natural frequency characteristic equation and IGSS algorithm to acquire in-plane natural frequencies and modal shapes. The results are compared with existing published results and FEA results, by which the accuracy, reliability, and efficiency can be validated. On this basis, some new results for the orthotropic rectangular plates with elastic boundary conditions are presented for the first time, which may be worked as the benchmark data.

2. Theoretical Formulations of Accurate Solutions

2.1. Differential Equations and Boundary Conditions

Figure 1 shows an orthotropic rectangular plate characterized by length Lx, width Ly, and thickness h in x-, y-, and z-directions, respectively. According to the small deformation, the stresses are written asin whichwhere Ex(Ey) and μx(μy) are Young’s modulus and Poisson’s ratio in x-direction (y-direction), respectively, and Gxy is the shear modulus.

The constitutive relations of the in-plane vibration of the plate, by Hooke’s law, can be expressed as

The strain energy and kinetic energy arewhere ρ is the mass density. As mentioned before, the aim of this study is to present an accurate solution for the free in-plane vibration of the orthotropic rectangular plate with various classical and elastic boundary conditions. Thus, to meet such requirements, the boundary simulation technique is introduced, where two opposite plate edges are implemented by the general elastic restraints (E) and other edges are simply supported restraints (S) containing the simply supported boundary 1 (SS1) and simply supported boundary 2 (SS2) to accord with the accurate solutions. The detailed information of the boundary conditions (two types: E-S-E-S and S-E-S-E) are shown as in Figure 2. Hence, the strain energy stored in these boundary springs during vibration is

In the next theoretical formulations, the boundary condition is the E-S-E-S case. Within any time interval, 0 to t0, Hamilton’s principle iswhere δUV, δT, and δUS are the virtual strain energy, the virtual kinetic energy, and the virtual spring potential energy, respectively, and δ is the first variation operator. Particularly, it should be pointed out here that the Hamilton's principle is a widely used principle in structural dynamics, such as linear and nonlinear vibration analysis of orthotropic or functionally graded plates.

Substituting for UV, T, and US from equations (4)–(6) into equation (7) yields

Relieving the virtual displacements, and , in equation (8) via the integration by parts and then applying the calculus of variations yields

The necessary and sufficient condition for the above equation is that the virtual displacements equal to zeros because of the arbitrariness of virtual displacements. So, from the first part, the governing equations of the plate can be reduced as

The elastically supported boundary conditions are expressed as

In the same way, the elastically supported condition for S-E-S-E case (Figure 2(b)) can also be obtained as follows:

2.2. Wave Solutions

In Section 2.1, the derivations are in the time domain. But, the free in-plane vibration analysis is in the frequency domain. Therefore, the Laplace transform is primarily applied to obtain the corresponding differential equations and boundary equations in the frequency domain. The expressions of Laplace transform and inverse Laplace transform arewhere β is the intersection point between the Bromwich contour of integration and the real axis on the complex s-plane.

The in-plane governing equations in equation (10) can be rewritten as the following form by the Laplace transform:where the coefficients of the linear operator L (Lij = Lji) are given as follows:

The accurate forms of generalized displacement functions for in-plane vibration of the orthotropic rectangular plate whose edges of y = 0 and y = Ly are with classical boundary conditions arewhere An and Bn, the arbitrary coefficients, correspond to the generalized displacement amplitudes for the nth mode of the plate, λn is the characteristic wave number in the x-direction, and Uyn(y) and Vyn(y) are in-plane modal functions for the plate whose y-direction is with different boundaries; and for the case of two parallel edges y = 0 and y = Ly with the simply supported boundary, their expressions are written as

Then, substituting the above wave solutions into equation (12), two homogeneous linear equations are obtained, which can be also written in the following matrix form:where the elements of Tij (i, j = 1,2) are

For equation (17), if and only if the determinant of the coefficient matrix equals to zero, its nontrivial solutions can be obtained and the characteristic equation of λn iswhere

For equation (19), when the wave number is zero, the rigid-body displacements occur, but in free vibration, such displacements will be eliminated. By solving equation (19), two pairs of characteristic wave numbers ±λ1n and ±λ2n can be obtained, and the corresponding basic solution vectors of equation (17), {αa,din, 1}T (i = 1, 2), are defined by

Then, the general solutions for the orthotropic rectangular plate displacements are obtained by the superposition principle:

To apply the MRRM further, a generalized displacement vector is introduced:wherewhere is the phase propagation matrix and and are the displacement coefficient matrixes. The boundary forces, similarly, can be written as a generalized displacement vector fn:where

2.3. Accurate Solutions for Free In-Plane Vibration of Rectangular Plates

The artificial spring boundary technique and the method of reverberation ray matrix are used together to derive the reverberation ray matrix approach in this section and then to conduct the analysis of in-plane vibration for orthotropic rectangular plates with general boundary conditions [2833].

In the orthotropic rectangular plates, the dual local coordinates (o-x12y12z12) and (o-x21y21z21) should be established to determine the phase relation of the reverberation ray groups. For information on double local coordinates, please refer to reference [27]. The origins of the two dual local coordinates, nodes “1” and “2,” are set at the right and left edges of the plate. For brevity, it is assumed that (o-x12y12z12) and (o-x21y21z21) are both paralleled with the global coordinate (o-xyz) but their x- and z-coordinates are in the opposite directions, and dual local coordinates satisfy the following relation:

In the dual local coordinate systems, the generalized state vectors and fn in equations (23) and (25) are written as

Furthermore, because of the uniqueness of the physical quantity, the compatibility condition of generalized state variables can be determined on the basis of the sign convention (seen in Figure 3):where Tδ = diag{−1, 1} is the dual transformation matrix between and and Tf = diag{1, −1} that between and . Substituting equations (27)–(29) into equations (30a) and (30b) and considering the wave numbers are independent of two dual local coordinates, i.e., (i = 1, 2), we havewhich can be regrouped as

Also, equation (32) can be further simplified aswhere and are the wave amplitude vectors andis the phase matrix expressing the phase relations of wave amplitude vectors. To ensure numerical stability of calculation, the exponentially growing function should not occur in the phase matrix, which indicates that the wave numbers λi (i = 1, 2) should have the positive real part.

Next, the system wave amplitude vectors are completely determined by the scattering relations at each of the nodes. According to the boundary conditions given in equations (11a) and (11b), the displacement continuity and force equilibrium conditions of the node are written aswhere and are stiffness matrices of the two sets of support springs. By changing the stiffness values, all possible boundary conditions can be simulated. Here, the boundary stiffness matrixes shown in Table 1 correspond to several boundary conditions often encountered in practice. Then, by classifications and reorganizations further, equations (33) and (34) are expressed in the following form of matrix:whereis the scattering matrix. Further combining equations (33) and (37) yieldswhere R = SPh is the reverberation ray matrix.

The equation of natural frequencies for free vibration is expressed as

For equation (40), it can be known that it is transcendental, which brings difficulty in solving analytically. Thus, an effective search technique is established by the author’s group based on the extrapolation method and golden section search (GSS) algorithm to obtain nontrivial solutions of equation (40). The detailed principle of the effective search technique can be seen in reference [27, 34].

3. Numerical Results and Discussion

Section 2 shows the theoretical formulations of the present accurate solutions. Thus, in this section, the primary purpose is to verify its accuracy and reliability via several numerical examples including the comparison of in-plane frequency parameters and modal shapes. On the basis of that some new in-plane vibration results of orthotropic rectangular plates with elastic boundary condition are presented.

Tables 25 show the first ten in-plane frequency parameters of orthotropic rectangular plates with different classical boundary conditions. In Tables 24, the finite element method (FEM) results taken from the ABAQUS software are also given as comparison due to the lack of the accurate benchmark data. In the FEM results, two types of the modeling methods are considered: FEMI: mesh size 0.025  m × 0.025  m, mesh type: S4R, and mesh number: 1600; FEMII: mesh size 0.01  m × 0.01  m, mesh type: S4R, and mesh number: 10000. In addition, the corresponding in-plane modal shapes of Tables 24 are presented in Figures 35, respectively. In Table 5, the results by Du et al. [12] using a semianalytical method and Liu and Xing [26] using the exact solutions are also considered here. In Tables 25, the dimensionless formula Ω is defined as ωLx(ρ/Gxy)1/2/π, and the geometric parameters and material parameters are given as follows: Lx/Ly = 1, μx = 0.3, Ey = 70 GPa, Ex = 2Ey, and Gxy = Ey/2/(1 − μxμy). From Tables 25, it is obvious that the results of present accurate solutions are close to the referential solutions. In addition, it is also clear from Tables 24 that the accuracy of the finite element method depends strongly on the mesh type and size of the finite element model. In theory, the smaller the size of the element, the more the number of the element, and the higher the accuracy of the calculation at this time; however, more computing resources need to be paid. This is why the purpose of developing exact solutions is to provide reference for numerical solutions [3, 4].

From Tables 25, it can also be concluded that the present approach has good accuracy and reliability in the analysis of in-plane vibration of orthotropic rectangular plates. However, in Section 1, it can be known that although relevant research has been carried out in the study of isotropic rectangular plates [27], no parametric study of the exact elastic solutions for the orthogonal materials has been found. Therefore, in order to enrich the existing research system and results, some parametric research work on elastic and exact solutions will be implemented based on the established model in Section 2. Firstly, the effects of the boundary springs on fundamental frequency parameters of E-SS1-E-SS1, E-SS2-E-SS2, SS1-E-SS1-E, and SS2-E-SS2-E orthotropic rectangular plates are presented in Figures 69, respectively. For simplicity, every time only one type of boundary springs changes from 10−2 to 1016 uniformly and the resting springs keep infinity. The plate is characterized by Lx = 1 m and Ly = 1.2 m. From Figures 69, it can be found that different spring parameters have different effects on the vibration characteristics of the orthotropic rectangular plates. For the convenience of expression, two different types of elastic stiffness coefficients are represented by unified symbols kn and kp. It must be noted here that the symbols kn and kp represent the spring stiffness in the normal and parallel directions, respectively. For examples of E-SS1-E-SS1 and E-SS2-E-SS2, kn is equivalent to and and kp is equivalent to and . Similarly, for SS1-E-SS1-E and SS2-E-SS2-E, kn is equivalent to and and kp is equivalent to and . For the elastic parameters kn and kp, the influence of spring kp on in-plane vibration characteristics is more severe. In addition, it can be found that regardless of the type of spring, when the spring stiffness coefficient is less than 105, the change of spring stiffness coefficient has little effect on the in-plane vibration characteristics of orthotropic rectangular plates in this region, while when the stiffness coefficient is between 105 and 1011, the influence of spring stiffness coefficient variation on in-plane vibration characteristics of orthotropic rectangular plates is very severe in this region. When the stiffness coefficient exceeds 1011, the change of spring stiffness coefficient has little effect on the in-plane vibration characteristics of orthotropic rectangular plates, so it can be considered as a clamped boundary condition. Based on this analysis, three types of elastic boundary restraints are defined and used in next calculation: E1: kn = kp = 106 N/m2, E2: kn = kp = 108 N/m2, and E3: kn = kp = 1010 N/m2.

Some new results of orthotropic rectangular plates with different elastically restrained edges are also performed in Tables 69. The geometrical dimensions of the plate are the same as those in Tables 69. The orthotropic ratios of the plate are Ex/Ey = 2, 10, 20, and 40. From Tables 69, it can be found that the frequency parameters of orthotropic rectangular plates increase gradually with the increase of orthotropic stiffness Ex/Ey, regardless of the boundary conditions. However, it is also obvious that when the orthotropic stiffness Ex/Ey value is small, its change has a great influence on the in-plane vibration characteristics of the orthotropic rectangular plates, such as the orthotropic stiffness Ex/Ey value changing from 2 to 10. When the stiffness value is large, its change has little effect on the in-plane vibration characteristics, such as the orthotropic stiffness Ex/Ey value changing from 20 to 40. Since the in-plane vibration analysis of orthotropic rectangular plates with elastically restrained edges is conducted for the first time, thus, the results in Tables 69 can provide the reference for future research on the titled problem.

4. Conclusions

Based on the method of reverberation ray matrix (MRRM) and improved golden section search (IGSS) algorithm, an accurate solution is proposed and free in-plane vibration of orthotropic rectangular plates subjected to various boundary conditions is studied. The boundary condition is defined as that a set of opposite edges are with one kind of simply supported boundary conditions, where the rest two edges are arbitrary boundary conditions. In the current framework, the boundary simulation technique is introduced, where two opposite plate edges are implemented by the general elastic restraints (E) and the other two edges are simply supported restraints (S) containing the simply supported boundary 1 (SS1) and simply supported boundary 2 (SS2) to accord with the accurate solutions. Thus, the accurate solution with arbitrary specified boundary conditions can be easily obtained. The present accurate solutions are validated by both mathematical proof and comparisons with FEM solutions. By comparing with the existing literature, the main academic contributions of this paper can be summarized as the following two points: (a) By contrast with the existing accurate solution methods like direct separation of variables, the present accurate method breaks the barrier that the existing accurate solutions are always restricted to classical boundary conditions and adapts to more widespread boundary conditions on the premise of guaranteeing the accuracy of calculation. (b) Based on the existing in-plane accurate calculation results of orthotropic rectangular plates, some new in-plane accurate results are presented in this paper, which further enriches the solution system and calculation results of in-plane characteristics of orthotropic rectangular plates.

Data Availability

All the underlying data related to this article are available upon request.

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 (grant no. 51905511).