Shock and Vibration
Volume 2021 (2021), Article ID 5540123, 20 pages
https://doi.org/10.1155/2021/5540123

Research Article

In-Plane Vibration Analysis of Square Plate with Multiple Cutouts

Shuangxia Shi,¹ Yinlang Qi,¹ Changqing Zhang,¹ Qingtai Meng,¹ Yinbei Liu,¹ Guoyong Jin,² and Bin Xiao¹

¹School of Energy and Power Engineering, Northeast Electric Power University, Jilin 132012, China
²College of Power and Energy Engineering, Harbin Engineering University, Harbin 150001, China

Correspondence should be addressed to Bin Xiao; dr.binxiao@yahoo.com

Received 23 January 2021; Revised 16 March 2021; Accepted 26 April 2021; Published 5 May 2021

Academic Editor: Sumeet S. Aphale

Copyright © 2021 Shuangxia Shi 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.

Abstract.

A theoretical solution method for the in-plane vibration characteristics of the plate with cutouts is proposed in this paper. The energy principle in conjunction with the Rayleigh–Ritz solution technique is adopted for the theoretical modeling of the in-plane vibration of the structure. The energy functions for the plate with cutouts are established by subtracting the energies of the cutout domains from the total energies of the whole plate. To ensure continuity over the entire solution domain, the in-plane displacements are composed of two-dimensional standard Fourier series and supplementary functions. The in-plane eigenmodes of the plates with different square cutouts are compared with the results obtained from the finite element method (FEM), with good agreements. The influences of the cutouts on the in-plane vibration characteristics of the plates with cutouts are investigated by varying the number, size, and position of the cutouts.

1. Introduction

The plates with cutouts are a kind of structure widely used in engineering fields, such as in the field of aeronautics, marine, and civil engineering. The cutouts within the structures are usually used to facilitate the connection of substructures. The in-plane vibration of the plates could not be ignored in the theoretical modeling for the connected structures; thus, investigating the in-plane vibration characteristics of the plate with cutouts is of great significance in the field of vibration control in engineering practice. If the in-plane vibration characteristics of the plate with cutouts are calculated by the finite element method (FEM), the results are more effective for low frequencies but require high computational cost for high frequencies. To reduce the computation time and reveal the physical mechanism of the influence of the cutouts within the plate on the in-plane vibration, analytical methods are more applicable for predicting the in-plane vibration characteristics of the plate with cutouts.

The number and size of cutouts have a great influence on the transverse vibration characteristics of the plate structures [1–9]. The rigidity of the structure with cutouts is mostly dependent on the number of cutouts. The vibrations of rectangular plates with one and two cutouts are, respectively, investigated by Avalos [1, 2]. Lam [3] used the Rayleigh–Ritz method to resolve the vibration problems of single opening and double opening structures. Moon [4] is concerned with the vibration analysis of a rectangular plate with a rectangular or circular hole. Lu [5] proposed an analytical solution of an elastic infinite plate containing an ellipse hole and a circle hole subjected to internal tractions and shears at infinity. The investigations on the effects of the cutout size are focused on the cutouts located in the center of the flexible plate [6–9]. Takahashi [6] used the Ritz method to analyze the rectangular plate with a central circular hole. Paramasivam [7] studied the effects of rectangular cutout sizes on fundamental frequencies of the structure with different types of boundary conditions by extending the grid framework model.

By comparison, the theoretical investigations on the in-plane vibrations of the plates are less than the transverse vibrations of the flexible plates. The Rayleigh–Ritz method [10, 11] and dynamic stiffness method [12–14] are the common theoretical methods for solving the vibration problems of structures. The in-plane vibration of isotropic rectangular plates with various boundary conditions is studied by the Rayleigh–Ritz method [10]. Gorman [12, 13] used the superposition method to analyze the in-plane vibration characteristics of the plates with different types of boundary conditions. At present, the researchers pay more attention to the theoretical solutions for the in-plane vibration analysis of the plate with cutouts [15–20]. The natural frequencies and mode shapes of the square plate with an inner cutout are obtained by the discrete Ritz method [16]. Huang [15] used the receptance approach to investigate the in-plane vibration characteristics of plates with circular cutouts. Zhou [18] investigated the free in-plane vibration problem of rectangular plates with classical boundary conditions using various numerical methods. A modified Fourier–Ritz approach is developed to investigate the free in-plane vibration characteristics of the orthotropic circular, annular, and sector plates with general boundary conditions [19]. Shi [20] adopted an improved Fourier series method and the first-order beam theory to study the in-plane free vibration behaviors of functionally graded carbon nanotube-reinforced composite circular arches and curved beams under elastic boundary restraints.

In studying the effects of cutouts within the structure, the cutouts are mostly located in the center of the structure, and the number of cutouts is less than two. Comparing with the investigations above, the theory proposed in this paper mainly could be used to predict the in-plane modal characteristics of the plate having arbitrary number and position of cutouts. In the existing methods, the cutouts within the structure are considered as a kind of structure with zero thickness. To behave the physical characteristics of the cutouts within the structure, the influence of cutouts on the vibration characteristic of the structure is considered in terms of energy [21]. The energy principle combined with the Rayleigh–Ritz solution is adopted to determine the vibration characteristics of the structures with cutouts. The effect of the cutouts is taken into account by subtracting the energies of the cutout domains from the total energies of the whole plate. The modal characteristics of square plates with different cutouts are predicted by the present method and verified by FEM. In the Results and Discussion, the influences of the number, size, and position of the cutouts on the natural frequencies and mode shapes of the plate with cutouts are investigated.

While previous researches focused mainly on relatively simple sandwich constructions and approximated the rib-stiffeners as an Euler-beam or a combination of translational spring and rotational spring, an exact theoretical model concerning the vibroacoustic response of more complex structures (e.g., two-dimensional sandwich structures orthogonally reinforced by periodic rib-stiffeners) is desirable. In addition to helping to explore the underlying physical subtleties, the model should also serve as benchmark checking for approximate analytical approaches, with a small computational expense afforded compared to numerical methods such as the FEM and the boundary element method (BEM).

2. Theoretical Formulations

2.1. Description of the Numerical Model

The geometric sketch of the plate with cutouts investigated in this paper is given in Figure 1. The boundary conditions of the rectangular plate with cutouts shown in Figure 1 can be represented by two kinds of springs with determined stiffness constants. The linear springs are attached to the edges in the normal or tangential direction. The dimensions in x and y directions for the rectangular plate with cutouts are represented as a and b, respectively. The rectangular cutouts positions are unrestricted in Figure 1. The length and width of the ith cutout are, respectively, (Xdi-Xci) and (Ydi-Yci), in which the size of the rectangular cutout is variable.

Figure 1: The geometric sketch of the rectangular plate with cutouts.

When the in-plane displacements of the rectangular plate are periodically extended as a standard Fourier series on the x-y surface, the discontinuities problems for the original displacement functions and their derivatives will encounter along the plate structure boundary. Therefore, the in-plane displacement functions of the plate structure are expressed as the modified Fourier series [22] to overcome the potential discontinuity of the original displacement functions and their derivatives. The in-plane displacement functions u1 in the x direction and v1 in the y direction for the rectangular plate are, respectively, expressed as

where and the supplementary functions in the in-plane displacement functions are expressed as follows:

It can be inferred that the supplementary functions are introduced to remove all the discontinuities potentially associated with the original displacement and their derivatives throughout the entire solution domain including the boundaries and then to effectively enhance the convergence of the results.

2.2. Solution Procedure of the Plate with Cutouts

Among the solutions of vibration analysis for the continuous systems, the energy principle combined with the Rayleigh–Ritz method is often adopted to solve the vibration characteristics of the rectangular plate due to the reliability and efficiency of the results [11, 22]. The energy function of the plate with cutouts is obtained by deleting the potential and kinetic energies on the cutout domains from the energy functions of the whole plate. Therefore, the Lagrange function for the rectangular plate with n cutouts is given as

where and T denote the potential and kinetic energy of the rectangular plate without cutouts, respectively. and Ti denote the potential and kinetic energy of the rectangular plate occupying the ith cutout, respectively.

The potential energy of the rectangular plate without cutouts can be expressed as

The kinetic energy of the rectangular plate without cutouts can be expressed as

where and represent the spring stiffness acting in the normal direction along the edge of y = 0 and the spring stiffness acting in the tangential direction along the edge of y = b, respectively. and represent the spring stiffness in the normal direction along the edge of x = 0 and the spring stiffness acting in the tangential direction along the edge of x = a, respectively. ω is the circular frequency. μ and ρ denote Poisson’s ratio and mass density of the rectangular plate, respectively. h denotes the thickness of the rectangular plate.

The potential energy of the ith cutout within the plate is expressed as

The kinetic energy of the ith cutout within the plate is expressed as

The boundary condition for the cutout in the rectangular plate is sensible with its geometric positions. If the cutout is located in the rectangular plate, then the boundary condition of the cutout is free. When the cutout is located on the edge of the rectangular plate, the edge of the boundary condition of the cutout will be the same as that of the rectangular plate.

Substituting equations (5)–(8) into equation (4) and applying the Rayleigh–Ritz procedure against each of the unknown Fourier series coefficients, the linear system in the matrix form can be described as

where

where K is the stiffness matrix of the structure studied, M is the mass matrix of the structure studied, Ki is the stiffness matrix of the ith cutout within the plate, and Mi is the mass matrix of the ith cutout within the plate. In the numerical calculations, the Fourier series is truncated to m = M and n = N. Once the coefficient vector E is obtained from the eigenvalue equation, the in-plane displacement of the rectangular plate with cutouts can be easily derived by substituting the coefficient vector into equations (1) and (2).

3. Results and Discussion

3.1. Convergence and Validation

The convergence of the theoretical modeling is studied by taking a square plate with an inner cutout as an example. The material properties of the square plate with clamped boundary conditions are taken as follows: length (a) = 3 m, thickness (h) = 0.0025 m, Poisson’s ratio (υ) = 0.33, density (ρ) = 2700 kg/m³, and Young’s modulus (E) = 7 × 10¹⁰ N/m². As shown in Figure 1, the size and original point of the inner cutout are (0.4 m × 0.4 m) and (2.2 m, 2.2 m), respectively. The modal characteristics for the plate with an inner cutout are given in Table 1 and Figure 2. The simulation results of the rectangular plate with an inner cutout are achieved from the software ANSYS, and the model is constructed by setting the unit size of PLANE42 as 0.05 m× 0.05 m. It can be seen from Table 1 that the natural frequencies for the rectangular plate with an inner cutout are verified by the comparisons with the results obtained from FEM and converge at M = N = 40. The phenomenon could be seen from the comparisons of mode shapes in Figure 2 that the mode shapes calculated by the current method are well-validated.

Table 1: The natural frequencies of the rectangular plate with an inner cutout.

Mode
M = N = 30
M = N = 32
M = N = 34
M = N = 36
M = N = 40
FEM

1
978.6
978.3
978.0
977.7
977.3
973.98

2
1003.4
1003.4
1003.4
1003.4
1003.4
1003.3

3
1162.6
1162.5
1162.5
1162.4
1162.3
1161.9

4
1453.6
1453.3
1453.0
1452.8
1452.5
1449.8

5
1634.8
1634.4
1634.2
1634.0
1633.6
1631.1

6
1637.7
1637.5
1637.5
1637.4
1637.2
1636.5

7
1660.3
1660.3
1660.3
1660.3
1660.3
1660.2

8
1921.5
1921.3
1921.1
1920.9
1920.6
1918.7

9
2013.1
2013.0
2013.0
2012.9
2012.8
2012.3

10
2023.6
2023.2
2023.0
2022.6
2022.1
2018.6

20
2737.2
2737.0
2736.9
2736.8
2736.6
2735.6

50
3939.1
3938.7
3938.5
3938.2
3937.9
3937.0

100
5489.0
5488.3
5488.1
5488.0
5487.8
5488.2

(a)

(b)

(c)

(d)

(e)

(f)

Figure 2: The mode shapes of the square plate with an inner cutout. (a) The first mode shape; (b) the second mode shape; (c) the third mode shape; (d) the fourth mode shape; (e) the fifth mode shape; (f) the sixth mode shape. Left: the ANSYS results; right: the present results.

3.2. The Influence of the Cutout Parameters

In this section, the influence of the number, size, and location of the cutouts on the in-plane vibration characteristics is investigated. In the parametric studies, the material properties and sizes of the square plates with clamped boundary conditions remain unchanged. The length, width, and thickness of the selected plate are a = 40 m, b = 40 m, and h = 0.0025 m. The material properties of the selected plate are considered as Poisson’s ratio (υ) = 0.33, ρ = 2700 kg/m³, and Young’s modulus (E) = 7 × 10¹⁰ N/m². Three plate structures with different numbers of cutouts are selected for the investigation on the influence of the number of cutouts. For the rectangular plate with four same cutouts having free boundary conditions, the coordinate positions of the origins for the square cutouts whose size is 5 m × 5 m are (12.5 m, 5 m), (22.5 m, 5 m), (12.5 m, 30 m), and (22.5 m, 30 m), respectively. For the rectangular plate with six same cutouts having free boundary conditions, the original points of the six cutouts whose size is 5 m × 5 m are (5 m, 5 m), (17.5 m, 5 m), (30 m, 5 m), (5 m, 30 m), (17.5 m, 30 m), and (30 m, 30 m), respectively. For the rectangular plate with eight same cutouts having free boundary conditions, the coordinate positions of the origins for the cutouts with 5 m × 5 m size are (4 m, 5 m), (13 m, 5 m), (22 m, 5 m), (31 m, 5 m), (4 m, 30 m), (13 m, 30 m), (22 m, 30 m), and (31 m, 30 m), respectively. The natural frequencies and mode shapes for the rectangular plate with different numbers of cutouts are shown in Table 2 and Figures 3–5, respectively. It can be seen from Table 2 and Figures 3–5 that the natural frequencies and mode shapes for the three cases obtained by the present method are in good agreement with those obtained by FEM. It can be found by the comparison in Table 3 that the more the number of cutouts, the lower the natural frequencies. The reason lies in that the increase of the cutouts makes the stiffness of the plate structure and the natural frequencies decrease. The comparison results between three cases in Figures 3–5 show that the mode shapes are influenced by the number of cutouts within the square plate.

Table 2: The natural frequencies of rectangular plates with different numbers of cutouts.

Mode number
Four cutouts
Six cutouts
Eight cutouts

Present
FEM
Present
FEM
Present
FEM

1
66.72
65.89
65.83
64.82
61.98
60.22

2
73.12
72.44
70.20
69.26
68.06
67.46

3
82.20
81.17
80.98
80.09
77.34
76.69

4
112.67
112.27
106.26
104.93
105.11
104.26

5
115.98
114.90
110.62
109.76
108.71
108.02

6
124.18
123.65
124.10
123.45
123.54
123.13

7
129.35
128.39
124.90
123.95
127.59
126.77

8
142.90
142.33
131.99
130.20
134.41
132.93

9
145.65
147.60
146.08
144.14
140.94
139.80

10
148.51
148.65
148.39
145.82
144.77
143.28

20
193.18
192.46
198.62
196.93
190.12
188.69

50
299.04
296.67
290.72
289.92
290.88
289.78

100
410.74
408.81
404.63
403.36
401.95
400.77

(a)

(b)

(c)

(d)

(e)

(f)

Figure 3: The first six mode shapes for the square plate with four cutouts. (a) The first mode; (b) the second mode; (c) the third mode; (d) the fourth mode; (e) the fifth mode; (f) the sixth mode. Upper: the present results; lower: the ANSYS results.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 4: The first six mode shapes for the square plate with six cutouts. (a) The first mode; (b) the second mode; (c) the third mode; (d) the fourth mode; (e) the fifth mode; (f) the sixth mode. Upper: the present results; lower: the ANSYS results.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 5: The first six mode shapes for the square plate with eight cutouts. (a) The first mode; (b) the second mode; (c) the third mode; (d) the fourth mode; (e) the fifth mode; (f) the sixth mode. Upper: the present results; lower: the ANSYS results.

Table 3: The natural frequencies of the rectangular plate with different cutout sizes.

Mode number
10 m ×  10 m
20 m ×  20 m
25 m ×  25 m

FEM
Present
FEM
Present
FEM
Present

1
70.61
70.95
83.29
83.54
103.38
103.82

2
74.30
74.32
86.94
87.04
104.36
104.52

3
85.96
86.13
97.35
97.61
136.02
136.08

4
107.49
108.27
102.06
102.63
137.12
137.29

5
115.16
115.42
120.23
120.34
144.29
144.35

6
123.59
123.85
126.40
127.08
145.73
146.13

7
128.91
128.94
141.30
141.50
152.01
152.43

8
144.00
144.28
144.73
144.99
157.66
158.06

9
149.48
149.73
158.55
158.72
164.98
165.07

10
155.57
155.68
167.84
167.92
169.78
170.33

20
209.85
210.01
226.07
226.36
241.77
241.95

50
298.27
298.37
322.88
323.77
358.95
359.15

100
416.65
416.77
460.46
460.95
499.36
500.90

As the in-plane modal characteristics of the structure are influenced by the cutout sizes, three cases of plates with an inner cutout are considered in the numerical discussions. The physical parameters for the plate structure discussed in this section remain unchanged. The sizes of the inner cutouts whose original points are (5 m, 25 m), (5 m, 15 m), and (5 m, 10 m) are (10 m × 10 m), (20 m × 20 m), and (25 m × 25 m), respectively. The natural frequencies and mode shapes for the selected three cases are shown in Table 3 and Figures 6–8, respectively. It can be seen from Table 3 that the comparisons of the natural frequencies for the three cases obtained by the present method and FEM are satisfactory. The phenomenon can be seen from Table 3 that the natural frequencies increase with the cutout size, which is attributed to the fact that the mass of the plate structure is decreased with the increase of cutout sizes. As shown in Figures 6–8, the mode shapes obtained by the present method are in good agreement with those obtained by FEM. The comparison of the mode shapes between three cases indicates that the mode shapes are influenced significantly by the cutout size.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 6: The first six mode shapes for the square plate with a cutout whose size is 10 m × 10 m. (a) The first mode; (b) the second mode; (c) the third mode; (d) the fourth mode; (e) the fifth mode; (f) the sixth mode. Upper: the present results; lower: the ANSYS results.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 7: The first six mode shapes for the square plate with a cutout whose size is 20 m × 20 m. (a) The first mode; (b) the second mode; (c) the third mode; (d) the fourth mode; (e) the fifth mode; (f) the sixth mode. Upper: the present results; lower: the ANSYS results.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 8: The first six mode shapes for the square plate with a cutout whose size is 25 m × 25 m. (a) The first mode; (b) the second mode; (c) the third mode; (d) the fourth mode; (e) the fifth mode; (f) the sixth mode. Upper: the present results; lower: the ANSYS results.

To investigate the influence of the cutout positions on the modal characteristics, the natural frequencies and mode shapes for three plates with double cutouts having the same size are calculated. Three kinds of plates with cutouts are selected to investigate the influence of cutout positions. As shown in Figure 9, the coordinate positions of the origins for the double inner cutouts with 10 m × 10 m size are (5 m, 15 m) and (25 m, 15 m), respectively. As shown in Figure 10, the coordinate positions of the origins for the double edge cutouts with 10 m × 10 m size are (15 m, 0 m) and (15 m, 30 m), respectively. As shown in Figure 11, the coordinate positions of the origins for the double corner cutouts with 10 m × 10 m size are (0 m, 30 m) and (30 m, 30 m), respectively. As shown in Table 4, the natural frequencies for the three cases from the present method and FEM are listed. We can see from Table 4 that the natural frequencies are in good agreement with those from the FEM and are greatly affected by the cutout positions. The first six mode shapes for the three cases from the present method and FEM are shown in Figures 9–11. An interesting phenomenon can be found that the distributions of the mode shapes are dependent on the cutout positions by comparing Figures 9–11.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 9: The first six mode shapes for the square plate with double inner cutouts. (a) The first mode; (b) the second mode; (c) the third mode; (d) the fourth mode; (e) the fifth mode; (f) the sixth mode. Upper: the present results; lower: the ANSYS results.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 10: The first six mode shapes for the square plate with double edge cutouts. (a) The first mode; (b) the second mode; (c) the third mode; (d) the fourth mode; (e) the fifth mode; (f) the sixth mode. Upper: the present results; lower: the ANSYS results.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 11: The first six mode shapes for the square plate with double corner cutouts. (a) The first mode; (b) the second mode; (c) the third mode; (d) the fourth mode; (e) the fifth mode; (f) the sixth mode. Upper: the present results; lower: the ANSYS results.

Table 4: The natural frequencies of the rectangular plates for different cutout positions.

Mode number
Inner cutouts
Edge cutouts
Corner cutouts

ANSYS
Present
ANSYS
Present
ANSYS
Present

1
65.55
66.18
60.83
61.26
69.10
69.12

2
77.36
77.78
73.13
73.39
70.24
70.25

3
83.79
84.31
80.68
80.94
82.97
83.00

4
103.43
104.18
109.81
109.89
102.56
102.60

5
113.76
114.51
113.00
113.01
117.27
117.42

6
114.00
114.60
119.77
119.97
118.67
118.87

7
127.32
128.43
127.25
127.27
123.85
123.87

8
130.58
133.21
129.20
129.25
138.92
138.93

9
138.36
139.06
141.20
141.96
148.05
148.16

10
140.72
141.39
142.41
143.49
148.29
148.30

20
197.79
198.16
199.63
199.66
204.38
204.51

50
297.65
298.65
300.66
301.29
307.39
308.11

100
424.80
425.83
426.14
426.87
423.74
424.11

4. Conclusion

The energy principle combined with the Rayleigh–Ritz method has been adopted to predict the in-plane modal characteristics of the plate with cutouts, in which the number, size, and position of cutouts are variable. In the theoretical formulation, the cutouts within the structure are considered by subtracting the energy occupied by the cutouts area from the total energy of the whole plate. The in-plane displacements of the plate with cutouts are expressed as improved Fourier series to ensure the smoothness and continuity over the entire solution domain. The convergence of the present method is checked by calculating the natural frequencies for the square plate with an inner cutout under different truncation series. The accuracies of the modal results of all the numerical examples considered in the Results and Discussion have been validated by the comparisons with those obtained from FEM. By investigating the influence of the number of the cutouts on the in-plane modal characteristics of the square plate, we can find that the increase of the number of cutouts will decrease the structural stiffness and the natural frequencies. The reason why the square plate has a larger cutout with higher natural frequencies is attributed to the fact that the structure mass decreased with the increase of the cutout sizes. It can be found that the distributions of mode shapes are influenced by the cutout positions by investigating the positions of the multiple cutouts.

Appendix

The representative expressions of the first row of submatrices in the matrices Ki, Mi, K, and M are given in this section. For the convenience of description, we will introduce some indices listed in the following:

Taking the derivatives of the Lagrangian energy functions with respect to the undetermined coefficients Bmn and making them equal to zero, the submatrices in the first row of Ki can be obtained as

The submatrices can be derived by replacing in submatrices with .

The submatrices can be derived by replacing in submatrices with .

The submatrices can be derived by replacing in submatrices with .

The submatrices can be derived by replacing in submatrices with .

Taking the derivatives of the Lagrangian energy functions with respect to the undetermined coefficients Bmn and making them equal to zero, the submatrices in the first row of Mi can be obtained as

The submatrices can be derived by replacing in submatrices with .

The submatrices can be derived by replacing in submatrices with .

The stiffness and mass matrix of the whole plate K and M can be obtained by replacing the upper limits (Xdi, Ydi) and the lower limits (Xci, Yci) in the above equations with (a, b) and (0, 0), respectively.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors gratefully acknowledge the financial support from the National Natural Science Foundation of China (nos. 51705071, 51775125, and 51475088), the Natural Science Foundation of Jilin Province of China (Grants nos. 20190103058JH, 20190201109JC, and 20140101065JC), the Fundamental Research Funds for the Central Universities of China (no. HEUCFG201713), and the Research Foundation of Education Department of Jilin Province, China (no. 2015-244).

References

D. R. Avalos, H. A. Larrondo, P. A. A. Laura et al., “Transverse vibrations of simply supported rectangular plates with rectangular cutouts carrying an elastically mounted concentrated mass,” Journal of Sound and Vibration, vol. 202, no. 4, pp. 585–592, 1997.

D. R. Avalos and P. A. A. Laura, “Transverse vibrations of simply supported rectangular plates with two rectangular cutouts,” Journal of Sound and Vibration, vol. 267, no. 4, pp. 967–977, 2003.

K. Y. Lam, K. C. Hung, and S. T. Chow, “Vibration analysis of plates with cutouts by the modified Rayleigh-Ritz method,” Applied Acoustics, vol. 28, no. 1, pp. 49–60, 1989.

M. K. Kwak and S. Han, “Free vibration analysis of rectangular plate with a hole by means of independent coordinate coupling method,” Journal of Sound and Vibration, vol. 306, pp. 12–30, 2007.

A. Z. Lu, Z. Xu, and N. Zhang, “Stress analytical solution for an infinite plane containing two holes,” International Journal of Mechanical Sciences, vol. 128-129, pp. 224–234, 2017.

S. Takahashi, “Vibration of rectangular plates with circular holes,” Journal of the Japan Society of Mechanical Engineers, vol. 4, no. 1, pp. 380–385, 1958.

P. Paramasivam, “Free vibration of square plates with square openings,” Journal of Sound and Vibration, vol. 30, no. 2, pp. 173–178, 1973.

R. Ali and S. J. Atwal, “Prediction of natural frequencies of vibration of rectangular plates with rectangular cutouts,” Computers & Structures, vol. 12, no. 6, pp. 819–823, 1980.

J. Zhang, T. Li, and X. Zhu, “Free vibration analysis of rectangular Fgm plates with a cutout,” Earth and Environmental Science, vol. 283, Article ID 012037, 2019.

N. S. Bardell, R. S. Langley, and J. M. Dunsdon, “On the free in-plane vibration of isotropic rectangular plates,” Journal of Sound and Vibration, vol. 191, pp. 459–467, 1996.

Y. Chen, G. Jin, and Z. Liu, “Flexural and in-plane vibration analysis of elastically restrained thin rectangular plate with cutout using Chebyshev-Lagrangian method,” International Journal of Mechanical Sciences, vol. 89, pp. 264–278, 2014.

D. J. Gorman, “Free in-plane vibration analysis of rectangular plates with elastic support normal to the boundaries,” Journal of Sound and Vibration, vol. 285, pp. 941–966, 2005.

D. J. Gorman, “Free in-plane vibration analysis of rectangular plates by the method of superposition,” Journal of Sound and Vibration, vol. 272, pp. 831–851, 2004.

C. Zhang, G. Jin, Z. Wang et al., “Dynamic stiffness formulation for free vibration of truncated conical shell and its combinations with uniform boundary restraints,” Shock and Vibration, vol. 2021, Article ID 6655035, 20 pages, 2021.

D. T. Huang, “Effects of constraint, circular cutout and in-plane loading on vibration of rectangular plates,” International Journal of Mechanical Sciences, vol. 68, pp. 114–124, 2013.

K. M. Liew, S. Kitipornchai, and A. Y. T. Leung, “Analysis of the free vibration of rectangular plates with central cut-outs using the discrete Ritz method,” International Journal of Mechanical Sciences, vol. 45, pp. 941–959, 2003.

I. Shufrin and M. Eisenberger, “Semi-analytical modeling of cutouts in rectangular plates with variable thickness free vibration analysis,” Applied Mathematical Modeling, vol. 40, pp. 6983–7000, 2016.

Y. Zhou, Q. Wang, D. Shi et al., “Exact solutions for the free in-plane vibrations of rectangular plates with arbitrary boundary conditions,” International Journal of Mechanical Sciences, vol. 130, pp. 1–10, 2017.

E. Ahad Fazl, Q. Wang, Q. Liang, and D. Shi, “A unified solution for free in-plane vibration of orthotropic circular, annular and sector plates with general boundary conditions,” Applied Mathematical Modelling, vol. 40, no. 21-22, pp. 9228–9253, 2016.

Z. Shi, X. Yao, F. Pang et al., “A semi-analytical solution for in-plane free vibration analysis of functionally graded carbon nanotube reinforced composite circular arches with elastic restraints,” Composite Structures, vol. 182, pp. 420–434, 2017.

S. Shi, B. Xiao, G. Jin et al., “Modeling and simulation of transverse free vibration analysis of a rectangular plate with cutouts using energy principles,” Shock and Vibration, vol. 2018, Article ID 9609745, 16 pages, 2018.

J. Du, “An analytical method for the in-plane vibration analysis of rectangular plates with elastically restrained edges,” Journal of Sound and Vibration, vol. 306, pp. 908–927, 2007.