Abstract

A modified Fourier–Ritz method is developed for the flexural and in-plane vibration analysis of plates with two rectangular cutouts with arbitrary boundary conditions, aiming to provide a unified solving process for cases that the plate has various locations or sizes of cutout, and different kinds of boundary conditions. Under the current framework, modifying the position of the cutout or the boundary conditions of the plate is just as changing the geometric parameters of the plate, and there is no need to change the solution procedures. The arbitrary boundary conditions can be obtained by setting the stiffness constant of the boundary springs which are fixed uniformly along the edges of the plate at proper values. The strain and kinetic energy functions of a plate with rectangular cutout are derived in detail. The convergence and accuracy of the present method are demonstrated by comparing the present results with those obtained from the FEM software. In this paper, free in-plane and flexural vibration characteristics of the plate with rectangular cutout under general boundary conditions are studied. From the results, it can be found that the geometric parameters and positions of the cutout and the boundary conditions of the plate will obviously influence the natural vibration characteristics of the structures.

1. Introduction

As one of the most commonly used structures, the rectangular plate has been in study for a long time. It is well known that, in a vibrating thin plate, there exist three types of modes: bending mode, longitudinal mode, and shear mode. The bending mode is referred to the out-of-plane vibration, while the longitudinal mode and shear mode are referred to the in-plane vibration [1]. The flexural vibration of thin plates with various boundary conditions has been attracting researchers for a long time. And a large number of papers on the flexural vibration have been published, such as the highly referenced [212]. In these papers, various methods including numerical, analytical, and semianalytical methods have been developed to analyze the vibration characteristics of plates under various boundary conditions. Maury et al. [3] studied the response of excited panels by the wave propagation method. Gorman [4, 5] and Mochida [8] used the superposition-Galerkin method to study the vibration characteristics of rectangular plates under free boundary condition. Li [9] proposed an exact series solution for the transverse vibration of beams under general elastic boundary restraints. Besides, flexural vibration of thin plates is also studied by other researchers with FE method, DQ method DSC method [1012], and so on.

Probably since the first nonzero natural frequency of in-plane vibration of a thin plate is much higher than the flexural vibration of the plate, compared with the transverse vibration problem, the in-plane vibration problem of plates is far less investigated in the past. As the in-plane vibration plays an important role in structural vibration, the in-plane vibration analysis has attracted attention of researchers in recent years. Many researchers have investigated the in-plane vibrations of plates with various methods [1320]. Gorman [13, 14] analyzed the free in-plane vibration characteristics of plates with different kinds of boundary conditions by the superposition method. Farag and Pan [17, 18] studied the modal characteristics of the plate which has two clamped parallel edges and two completely free parallel edges. It should be mentioned that Xing and Liu [19] solved all possible exact solutions about in-plane vibrations of rectangular plates.

Due to practical engineering requirements, it is often necessary to open a hole on a plate, for example, opening a hole on a steel plate to reduce the weight of the structure. So, it is necessary to study the in-plane and flexural vibration of rectangular plates with cutouts. As the transverse vibration of the plate has been studied for a long time, the transverse vibration of the rectangular plate with all kinds of cutouts also has been extensively studied. Several ways have been applied to handle the vibration problem of plates with cutouts. As illustrated in [21], the cutout on the plate is considered as an extremely thin plate. Kwaka and Han [22] used the coordinate coupling method to study the vibration of the plate with cutout under free boundary condition. Laura et al. [23] studied the out-of-plane vibration of the plate with rectangular cutout with simple supported boundary condition. And the out-of-plane vibration of plates with two rectangular cutouts is studied in Reference [24]. In this paper, plates with different kinds of rectangular cutouts are studied.

Just as illustrated previously, neither the in-plane vibrations of the plate nor the plate with different kinds of cutouts is seldom studied. Among the available papers, Shufrin and Esienberger [25] used the multiterm extended Kantorovich method to study the in-plane vibration of the plate with rectangular cutouts. In their paper, the in-plane vibration of the rectangular plate with inner or corner cutout is studied. Chen et al. [26] used the Chebyshev–Lagrangian method to study the out-of-plane and in-plane vibration of the rectangular plate with rectangular cutout.

In this paper, a modified Fourier–Ritz method is employed to analyze the vibration of the rectangular plate with cutout under arbitrary boundary conditions, aiming to provide a unified procedure for many cases, such as plate with various locations (inner corner or edge) or sizes of rectangular cutout, or different kinds of boundary conditions. The arbitrary boundary conditions can be obtained by setting the stiffness constant of the boundary springs at proper values. Under the current framework, the in-plane and flexural displacements of the rectangular plate with cutout are expanded in the form of standard Fourier series and eight auxiliary polynomials. The use of the auxiliary polynomial functions can eliminate potential discontinuities and accelerate convergence of the expanded Fourier series. The natural frequencies and mode results of plates with rectangular cutouts are obtained by using the Rayleigh–Ritz procedure.

2. Theoretical Formulations

2.1. The Model and Stress-Strain Relations of the Restrained Plate

Consider an elastically restrained plate with length , width , and thickness , as depicted in Figure 1. A Descartes coordinate is fixed at the middle surface of the rectangular plate. x, y, and z represent the directions of length, width, and thickness, respectively. , , and are employed to denote the displacements at the midplane of the rectangular plate in directions of , , and , separately. Under the current framework, the arbitrary boundary conditions of a rectangular plate can be simulated by setting three groups of linear springs () and a group of rotational spring () which are fixed uniformly along the boundaries of the rectangular plate at proper values. And arbitrary boundary conditions can be achieved by setting the stiffness constants of the four sets of the springs at proper values. According to Kirchhoff assumptions, the displacements of a thin plate in , , and directions can be written as follows:where , , and represent the displacements of the plate in directions of , , and , respectively. Based on the thin plate theory, the relationship between the strains and displacements of thin plates can be written as follows:where and are the normal strains in the midplane of the rectangular plate and is the shear strain in the midplane of the rectangular plate. and represent curvature changes and represents twist change, separately. They can be expressed as

Equations (3)–(8) represent the relationship between the strain and displacement of a thin plate in Descartes coordinates. Based on Hook’s law, the stress-strain equations are described aswherewhere , and represent the normal stress of a rectangular plate in the direction of , , and , respectively. and are Young’s modulus and Poisson’s ratio of the plate, respectively.

Performing the integration of stresses () over the thickness of the plate, we get

Performing Equations (11) and (12) yields, we getwherewhere and are defined as the force resultants and is defined as the shear force resultant. are defined as bending resultants and is defined as the twist moment resultant, respectively.

By now, all the stress-strain relations and other needed equations are obtained, and they will be used to derive the energy functions of the rectangular plate with small deformation.

2.2. Energy Functions

In this paper, the Rayleigh–Ritz procedure is employed to solve the vibration problems of a rectangular plate. Thus, the energy functions of the plate should be obtained first. For vibration analysis, the Lagrange energy function () of the rectangular plate is defined aswhere , and are kinetic energy of the rectangular plate, potential energy of the rectangular plate, and the potential energy stored in the boundary springs, respectively.

Based on the classical thin plate theory (CPT), the potential energy () of a rectangular plate in Equation (15) can be rewritten as

Inserting Equation (13) into Equation (16), the potential energy equation of the rectangular plate can be described as

And the kinetic energy () of the rectangular plates is written aswhere and are the density and the thickness of the rectangular plate, respectively.

As Laura et al. illustrated in [23], a rectangular plate with a hole can be regarded as a structure combined by two kinds of different materials. And the material properties of the hole on the plate can be considered as , and . Thus, the Lagrange energy equation of a rectangular plate with rectangular cutout can be obtained in this method. In order to describe this problem clearly, one specific example is presented here. Consider a rectangular plate with one inner cutout, as shown in Figure 2. The length and width of the plate are and . The location of the cutout can be described by the coordinates of the corners of the cutout, and the coordinates of the corners are and . Thus, the potential energy equations of the structure depicted in Figure 2 can be written aswhere

As illustrated previously, the general boundary conditions can be achieved by setting the stiffness constant of boundary springs at proper values. In this paper, , , and are adopted to represent the stiffness constants of the boundary springs which are fixed along the four edges of the rectangular plate x = 0, y = 0, x = a, and y = b, respectively. Thus, the potential energy stored in the boundary springs can be written as

According to the Rayleigh–Ritz procedure, choosing suitable displacement functions is of great importance to make the approximate solutions converge to the exact solutions.

The displacement of a rectangular plate can be theoretically expressed in any group of orthogonal polynomials. However, since computers have the round-off errors, the high-order polynomials generally become numerically instable. For the plates under arbitrary boundary conditions, the displacement functions expressed in the form of standard Fourier cosine or sine series can avoid this problem, but discontinuities of the displacement functions and their derivatives will potentially exist along the boundaries of the rectangular plates. This kind of problem is demonstrated in detail [2733]. In order to overcome this difficulty, the modified Fourier series is employed, and the displacement of the plate is expressed as follows:where , and are the coefficients of expansions in Equations (23)–(25); are the coefficients of eight auxiliary polynomials in Equations (26)–(28); and and are eight auxiliary polynomial functions. Displacements in the midplane expressed in the form of standard cosine Fourier series plus eight auxiliary polynomials can eliminate the potential discontinuity along the boundary of the plate. The eight auxiliary polynomial functions are

Specifically, symbols and (, ) represent independent variables and length, respectively. And it is obvious that

As the Lagrange energy function and the displacement functions of the rectangular plate are established, the next task is to obtain the value of the coefficients in the displacement functions. Substituting Equations (16)–(20) into the Lagrange energy function and minimizing the Lagrange energy functions by making the derivatives of the unknown variables of expansions equal to zero as

Finally, totaled equations are obtained and can be written in the form of matrix aswhere , and are stiffness matrix, mass matrix, and coefficient matrix, respectively. They are described as

The element in stiffness and mass matrix is given in Appendix. The vector can be written as follows:

Based on the classical thin plate theory (CPT), the in-plane and out-of-plane vibrations are decoupled, and thus the out-of-plane and in-plane vibration of a plate structure can be studied separately. When studying the flexural vibrations of the plate structure, the stiffness mass and coefficients matrix of the structure can be defined as

When studying the in-plane vibrations of a rectangular plate, the stiffness mass and coefficients matrix can be written as

3. Validation and Discussion

With the formulations derived in Section 2, several numerical results together with those available in the literature will be presented in this section to validate the accuracy and convergence of the current method. Besides, the vibration characteristics of the structure will also be studied.

3.1. Flexural Vibration Analysis of Rectangular Plate with Rectangular Cutouts

First, the flexural vibration of rectangular plate with one inner, edge, or corner cutout will be studied. These three kinds of models are depicted in Figure 2. The length and width of the considered rectangular plate are , and the thickness of the plate is . And the material parameters of the plate structure are defined as , , and . The coordinate of the points shown in Figure 2 are (a) , (b) , and (c) . For simplicity, a string consisting of several letters is adopted to represent the boundary conditions of the structure depicted in Figure 2. For example, “CCCC-FFFF” denotes the boundary conditions of the edges corresponding to number one to eight, and “C” and “F” represent clamped and free boundary conditions, respectively.

The first six nonzero natural frequencies of flexural vibration of plate structures with one cutout are listed in Table 1. The results obtained by software ANSYS are also presented here to verify the correctness of the present method. Comparing the results listed in Table 1, it can be observed that these two sets of results are very close. Besides, it can be also found that the position of the cutouts has an effect on the vibration characteristics of the structure. For example, frequency parameters in Table 1 show that the plate with one corner cutout has the largest nonzero first frequency while the plate with one edge cutout has the smallest nonzero first frequency; structures considered here are under completely free boundary conditions, and all the cutouts are also kept the same size.

By solving Equation (32), the natural frequencies and eigenvectors () can be obtained simultaneously. By inserting the eigenvector into the displacement expressions (Equations (23)–(25)), the mode shapes of the plate with one inner cutout can be acquired. The mode shapes of CCCC-FFFF plate with one inner cutout are compared with those obtained by ANSYS in Figure 3 (the right one). It can be seen that the mode shapes obtained by the present method match with those acquired by software ANSYS very well. Besides, it can also be found that with the increase of the frequency, the mode shapes of the structure become more and more complex.

The nonzero frequency parameters of the rectangular plate with one rectangular cutout in Table 1 are calculated by selecting the truncate number as . In order to check the convergence of the current method, the natural frequencies of the plate with two inner cutouts are calculated with different truncation numbers (i.e., ). The models of the rectangular plate with two cutouts are depicted in Figure 4. In the following discussion, the geometric parameters of the plate with two cutouts are considered as 2 m in length, 2 m in width, and 0.003 m in thickness, and the material properties of the plate are defined the same as previous. The coordinates of the points shown in Figure 4 are (a) , (b) , and (c) . The frequency results of the plate with two inner cutouts are listed in Table 2. It shows that the natural frequency converges at , and the convergence and numerical stability of the current method are validated.

As mentioned previously, the current method can be used to study the vibration characteristics of the rectangular plate with rectangular cutout under general elastic-restrained boundary. And the general elastic restrained boundary can be obtained by setting the stiffness constants of the boundary springs at proper values. For example, the clamped boundary conditions can be achieved by setting the stiffness constants of the linear and rotational springs as infinite. In practical calculations, a sufficiently large number is always employed to represent the “infinite.” In the first example, the clamped boundary conditions are obtained by setting the stiffness constant at 1012. So, effects of the stiffness constant of boundary springs on the solution results ought to be studied. The frequency parameters for EEEE-FFFF-FFFF plate with two inner rectangular cutouts are listed in Table 3, and the alphabet “E” represents elastic boundary condition. It can be found that as the stiffness constants of the springs increase, the natural frequencies of the structure increase gradually, and when the stiffness constant is beyond 1012, the natural frequency remain unchanged. Thus, in the following discussion, 1012 is selected to represent the value of “infinite.”

The mode shapes of the structure that rectangular cutouts locate inside of the plate or on the boundary of the plate or on the corner of the plate are shown in Figures 57, and the geometric and material parameters of the plate are kept the same as mentioned previously. Figures 57 show that structures shown in Figure 4 has similar mode shapes in low frequency. Besides, the mode shapes of these structures have symmetric characteristics.

As illustrated previously, the current method is able to handle the cases with arbitrary boundary conditions. The first six nonzero frequency parameters of the structure with one rectangular cutout under four kinds of boundary conditions are listed in Table 4. The rectangular plate considered is in the length of 2 m, width of 1 m, and thickness of 0.003 m, and the corner points of the cutout are located at A (0.8 m, 0.4 m) and B (1.2 m, 0.6 m). The material properties of the structure are , , and . It should be mentioned that the alphabet “S” and “C” are used to represent the simply supported and clamped boundary conditions, respectively. From Table 4, it can be found that the present results are very close to the results obtained by software ANSYS.

3.2. In-Plane Vibration Analysis of Rectangular Plate with Cutouts

The flexural vibration of the rectangular plate with cutouts has been studied in the previous section. It proves that the current method is capable of investigating the flexural vibration characteristics of the plate with rectangular cutout. Next, the in-plane vibration of the structure will be investigated. The natural frequencies and mode shapes for in-plane vibration of rectangular plate with cutouts will be presented in this section.

In order to verify the correctness of the present method, plate with one cutout is studied firstly. The first example is a classical case: plate with one cutout under the clamped boundary condition. The clamped boundary conditions can be obtained by setting the stiffness constants of the springs at 1012. For this example, the geometric parameters of the rectangular plate considered here are and , the material properties of the plate are defined as , , and , and the truncated numbers are selected as . The coordinates of the opposite corners of the cutouts shown in Figure 2 are (a) , (b) , (c) . The first six frequency parameters for in-plane vibration of structures with one cutout are presented in Table 5. And the mode shapes of the CCCC-FFFF plate with one inner cutout are also presented here. From Table 5 and Figure 8, it can be seen that the present results show an excellent agreement with results acquired by software ANSYS.

The effects of the boundary conditions on the vibration of the plate with rectangular cutout are considered here. According to [17], there exist two types of simply supported boundary conditions for in-plane vibration. The first type can be considered as stiffness constants of the tangential springs and normal springs are infinite and zero, respectively, while the second one can be seen as stiffness constants of the tangential and normal springs are zero and infinite, separately. For simplicity, these two kinds of boundary conditions are denoted by “”and “,” respectively. The geometric and material parameters of the plate are same as the ones defined in Table 5. Table 6 shows the natural frequencies of the plate with rectangular cutout located inside of the plate under different kinds of boundary conditions. It can be clearly found that boundary conditions influence the vibration characteristics of the plate structure obviously. The first natural frequency of the same structure varies with different boundary conditions.

Next, the convergence of the current method for in-plane vibration will be studied. Structures considered here have the same geometric and material parameters as the plate with two cutouts defined in Section 3.1. Table 7 shows the frequency parameters of the FFFFFF-FFF-FFF plate with two edge rectangular cutouts determined by different truncated numbers. It can be found that as the truncated number increases the natural frequencies for in-plane vibrations change little at first, and when the truncated number is beyond 28, the natural frequencies remain unchanged. The desired convergence and numerical stability of the current method for in-plane vibration are validated.

As mentioned previously, the mode shapes of the plate with cutouts can be obtained by inserting the eigenvector into the displacement expressions. The first six mode shapes for in-plane vibration of the plate with two cutouts are shown in Figures 911, respectively. Compared with the flexural vibration of the plate with two cutouts, it can be found that not all mode shapes for in-plane vibration of the structure have symmetric characteristic.

4. Conclusions

A modified Fourier–Ritz formulation is presented for the flexural and in-plane vibration of plates with different kinds of cutouts (inner, edge, or corner). Under the current framework, the in-plane and flexural displacements of rectangular plates are expanded in the form of standard Fourier cosine series plus eight auxiliary polynomial functions. The current method is capable of universally dealing with cases with various positions of cutouts or different kinds of boundary restraint. The modification of the position or size of the cutouts is just like modifying the geometric parameters of the plate, and there is no need to change the solution procedure. The arbitrary elastic boundary conditions can be achieved by setting the stiffness constant of the boundary springs at proper values, for example, clamped boundary condition can be achieved by setting the stiffness constant of the boundary springs to infinite. Although, only uniform elastic restraints along the edge of the rectangular plate are studied in this paper, the modified Fourier–Ritz method can also be used for nonuniform elastic supports.

Several examples are presented in this paper to validate the desired accuracy and convergence of the current method. The effects of stiffness constants of the boundary springs on the solution results are also studied. It should be mentioned that with current method, the vibration of the plates with different kinds of cutouts or boundary conditions can be studied in a unified solution procedure. The vibration characteristics of the plate structure with rectangular cutouts are studied. The results show that both the position of cutouts and boundary conditions of the plate will influence the vibration characteristics of the structure obviously. The mode shapes for flexural and in-plane vibration of the plates with cutouts located inside of the plate or on the edge of the plate or on the corner of the plate are plotted in this paper.

Appendix

Detailed Expressions for the Elements in Mass Matrix and Stiffness Matrix

The detailed expressions of the stiffness and mass matrices in Equations (33) and (34) are given as follows: to make the expressions simple and clear, six indexes are predefined:

The other elements in stiffness matrix are similar to those derived previously.

The other elements in mass matrix are similar to those mentioned previously.

Data Availability

All relevant data are within the paper, and data sources are public.

Conflicts of Interest

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

Acknowledgments

This paper was supported by the National Natural Science Foundation of China (Nos. 51822902, 51775125, and 51279035) and the Fundamental Research Funds for the Central Universities of China (No. HEUCFG201713).