# Free Vibration Analysis of Thick Rectangular and Elliptical Plates with Concentric Cut-Out

Authors: Anjibabu Merneedi; Nalluri Mohan Rao; L. Natrayan; L. Yuvaraj; Prabhu Paramasivam
Journal: Advances in Materials Science and Engineering
 (2021)

---

## Abstract

This research paper deals with a numerical method which is modified and applied, by the authors to derive an eigenvalue of a thick plate having cut-out in which geometries of plate and cut-outs are different, through a deflection matching condition by including shear deformation and rotary inertia effects, with less computational efforts and high accuracy. The modified Independent Coordinate Coupling Method (ICCM) is validated with FEM package (ANSYS) and applied to know the change in eigenvalues for a plate with cut-out by varying various parameters like aspect ratios, cut-out size, and thickness ratios. Trigonometric functions considered at the boundary level conditions of a simply supported plate should be satisfied. Free vibrational exploration on a thick isotropic plate with various aspect ratios and an elliptical plate with various sizes is carried out through the modified ICCM. Independent coordinates are applied for a plate domain and for a hole domain individually followed by equating the deflection condition of hole and plate, a reduced mass to express with cut-out from which eigenvalues can be obtained. The deflection matching condition facilitates the analysis even though the geometries of plate and cut-outs are different.


---

## Body

## 1. Introduction

A plate with a cut-out is a commonly existing structure in the domain of aeronautics, civil, marine, and mechanical engineering to accommodate mountings or accessories under the action of various forces. In dealing with such cases, knowledge on dynamics of plates plays a vital role. Plates with cut-outs are subjected to free vibration. Eigenvalues of thin plates may be calculated using classical plate theory and easily obtained, whereas they cannot be applied directly in thick plate analysis due to the effect of rotary inertia. This leads to an overestimation of results. The expression of strain and potential energies of plate must contain the terms of shear and rotary inertia to improve the accuracy of eigenvalues while computing.From the literature available on the Mindlin plates, Senjanovic et al. [1, 2] derived exact formulas for a plate by reducing a three-equation system to one equation expressing a potential function in terms of bending deflection. Xiang et al. [3] analyzed Mindlin plates by implementing delta type of kernel wavelet in lateral vibration of plates. Merneedi et al. [4, 5] applied the ICCM method of elliptical and rectangular thin plate by changing the size, position, and number of holes at all possible ways. Aksu [6] applied an energy approach along with a finite difference in a Mindlin’s dynamic analysis approach of plate having cut-outs. Rui et al. [7] applied a rational method of superposition in the symplectic space for Mindlin plate’s vibration analysis. Xing and Liu [8, 9], implemented a new two-eigenfunction theory to obtain vibrational solutions of a rectangular Mindlin plate, in which amplitude and generalized curvatures are considered as functions.Liew et al. [10–13] minimized energy function of Mindlin’s theory by using the Rayleigh–Ritz method to formulate the governing eigenvalue equations. Dawe and Roufaeil [14] derived by assuming that plate deflection and rotations are a sequence. The Rayleigh–Ritz approach may be used to compute the natural frequencies of flexural vibrations of a square plate using products of appropriate Timoshenko beam functions. Leissa [15, 16] investigates the free vibrational exploration for twenty-one cases of various rectangular plates through precise analytical techniques. Sayyad and Ghugal [17], proposed a displacement-based theory in which thickness coordinates are used as an exponential function so that transverse shear stress can be obtained directly without any shear correction factor. Kanak and Salil [18] and Merneedi et al. [5] consider the nine-node isoparametric element. The natural frequencies of thin and thick plates with cut-outs were analysed using the first-order shear deformation theory. Prasad et al. [19] utilized the Rayleigh–Ritz approach to apply a three-term deflection function on clamped and simply supported elliptical plates. Leissa [20] calculated the precise fundamental natural frequencies of a simple sup. Maurizi et al. [21], derived fundamental natural frequencies for a vibration of elastic plate using simple static in the presence of a clamped boundary condition. Lam et al. [22] created a computational tool to determine the natural frequencies of circular and elliptical plates; the researchers used the Rayleigh–Ritz method using a set of orthogonal functions called acceptable functions. Shibaoka [23] creates a fundamental normal mode of flexural vibration of an elliptical plate under clamped-clamped boundary conditions. Mcnitt [24] utilized the Galerkin method and an ordinary product solution by ignoring the impact of rotational inertia in the investigation of free vibration of an elliptical plate under clamped boundary conditions using traditional small-deflection theory. Kwak and Han [25] used ICCM to do a free vibrational exploration on a rectangular plate with a hole obtaining energy linked to the hole with plate separately and equated by matching condition of deflection. Hasheminejad et al. [26] formulated the analytical expression for the elliptical plate with an eccentric cut-out. The solution is derived using Helmholtz’s decomposition theorem, and this technique is based on Navier’s displacement equation of motion. Mathieu functions of translational addition theorems are used to separate the variables.According to the literature, eigenvalues of a thick rectangular plate with a hole were estimated using several numerical tools such as the Rayleigh–Ritz methodology, FEM, and finite difference techniques and the plate is same. It is identified that if hole and plate are of different geometries, the RRM may lead to complex integrals with tedious numerical computations because of its integral limits [27]. Keeping in view of this, in the present work, the authors proposed a novel modified ICCM to determine the natural frequencies of a structure with cut-out of different shapes by implementing a deflection matching condition along with closed-form integrals to reduce the computational efforts. Shear deformation and rotational inertia effects are among them. Even if the shape of the plate and hole differs, the deflection matching condition simplifies the analysis, which is a key benefit of the ICCM. Using the ICCM, a Mindlin elliptical and rectangular plate’s first six natural frequencies with various sizes of cut-outs under simply supported boundary conditions were derived to show the variation of natural frequencies with respect to various hole sizes, and the same is validated through ANSYS software including mode shapes, which was not covered in the literature in view of the authors.


### 1.1. Free Vibrational Exploration of Rectangle Plate with Cut-Out

Separate energy expressions are formulated for a hole domain and for a plate domain exploration of a plate with hole, and hence deflection matching condition is essential to establish a kinematic relation between these independent coordinate domains. It will also enable us to unify these two independent energies.


#### 1.1.1. Expressions for Rectangular Plate

Consider a thick rectangular plate with a cut-out inx−y plane, having dimensions of a×b and ac×bc, respectively. Expressions for rectangular plate strain energy and kinetic energy from Reference [11] are as follows:(1)Vp=12∫vEh21−μ2∂ψx∂x+∂ψy∂y2−21−μ∂ψx∂x∂ψy∂y−14∂ψx∂y+∂ψy∂x2+κGψx+∂w∂x2+ψy+∂w∂y2dv,Tp=12ω2∫Aρhw2+112ρh3ψx2+ψy2dA,where h = plate thickness, ρ = mass density, ϑ=E = modulus of elasticity, w=transverse displacement, G = shear modulus, μ=Poisson’s ratio, ψx=−∂w/∂x = rotation about y direction, ψy=−∂w/∂y = rotation about x direction, κ = shear correction factor = 5/6, and ω = angular frequency.Introducing a non-dimensional parameter, thickness to width ratioβ=h/b, aspect ratio α=a/b, and non-dimensional variables ξ=x/a, η=y/b, the energy expression of plate is a non-dimensional representation of plate which can be written as follows:(2a)Vp=∫01∫011α∂ψx∂ξ+∂ψy∂η2−21−μ1α∂ψx∂ξ∂ψy∂η−14∂ψx∂η+1α∂ψy∂ξ2+61−μκβ2ψx+1α∂w∂ξ2+ψy+∂w∂η2abdξdη,(2b)Tp=∫01∫01w∗2+112β2ψx2+ψy2abdξdη,where w∗=w/b.Admissible functions are considered to expand transverse deflection and the rotations in equations (2a) and (2b) as(3)wξ,η,t=A¯iξ,ηΩt,ψxξ,η,t=B¯iξ,ηΩt,ψyξ,η,t=C¯iξ,ηΩt,where A¯iξ,η=A1A2A3A4…Am is a 1×m matrix and Ωiξ,η=Ω1Ω2Ω3Ω4…ΩmT is m×1 generalized coordinate vector, where “m” denotes admissible functions required for the approximating the deflection curve.Substituting equation (3) in equations (2a) and (2b) yields(4a)Tp=12Ω˙TMΩ˙,(4b)Vp=12ΩTKΩ,in which(5a)M=ρhb2M¯,(5b)K=Db2K¯,where(6a)M¯=∫01∫01A¯T.A¯+β122B¯T.B¯+C¯T.C¯dξdη,(6b)K¯=∫01∫011α2∂B¯T∂ξ⋅∂B¯∂ξ+∂C¯T∂η⋅∂C¯∂η+1α∂B¯T∂ξ⋅∂C¯∂η+∂C¯T∂η⋅∂B¯∂ξ+1−μ⋅1α∂B¯T∂ξ⋅∂C¯∂η+∂C¯T∂η⋅∂B¯∂ξ−12∂B¯T∂η⋅∂B¯∂η+1α2∂C¯T∂ξ⋅∂C¯∂ξ+1α∂B¯T∂η⋅∂C¯∂ξ+∂C¯T∂ξ⋅∂C¯∂η+61−μκβ2B¯T⋅B¯+1α2∂A¯T∂ξ⋅∂A¯∂ξ+1αB¯T⋅∂A¯∂ξ+∂A¯T∂ξB¯+C¯T⋅C¯+∂A¯T∂η⋅∂A¯∂η+C¯T⋅∂A¯∂η+∂A¯T∂ηC¯dξdη.In equations (6a) and (6b), stiffness matrices are represented by M¯, K¯, respectively, along with non-dimensional variables. Equation of motion of plate is obtained by substituting equation (3) in Lagrange’s expression, and its eigenvalue can be computed from(7)K−ω2MΛ=0.Using non-dimensional matrices introduced in equations (4a) and (4b), the above eigenvalue problem in non-dimensional form(8)K¯−ω¯2M¯Λ=0will give non-dimensional eigenvalues of a plate, in which “Λ” is an eigenvector. In equation (7),(9)ω¯=ωb2π2ρhDis known as the frequency parameter which represents the flexural rigidity in D=Eh3/121−μ2, which is known.To obtain in terms of admissible function matrices in theX and Y directions, an admissible function matrix from equation (3) has to be changed.(10)A¯iξ,η=Aiξaiη,B¯iξ,η=Biξbiη,C¯iξ,η=Ciξciη,where,i=1,2,3…m.


#### 1.1.2. Boundary Conditions

A thick simply supported plate is considered and its eigenvalues are derived by considering an individual eigenfunction in each direction for each admissible function.(11)Aiξ=2siniπξ,aiη=2siniπη,Biξ=2cosiπξ,biη=2siniπη,Ciξ=2siniπξ,ciη=2cosiπη.Admissible functions in equation (11) are taken into account in such a way that they meet the assumed boundary conditions of the plate. All geometrical end points of ξ,η in equation (3) along with equation (11) will result in a zero deflection.


#### 1.1.3. Energy Expressions for Rectangular Plate with Cut-Out Using ICCM

Kinematic relations are formulated for a plate having cut-out using ICCM; by matching the plate’s deflection to the hole’s deflection, these relationships may be used to produce updated mass and stiffness matrices as well as distinct plate and hole domains having their own coordinates. Fixing local coordinates to the hole domain is presented for rectangular hole in rectangular plate as illustrated in Figure1.Figure 1
Rectangular plate opening with cut-out.Taking non-dimensional coordinates into account, displacement and rotations are expressed asξh=xc/ac and ηh=yc/bc.(12)whξh,ηh,t=A¯hiξh,ηhΩt,ψhxξh,ηh,t=B¯hiξh,ηhΩt,ψhxξh,ηh,t=B¯hiξh,ηhΩt,where A¯hiξh,ηh=Ah1Ah2Ah3…Ahm is an 1×mh function matrix and Ωht=Ωh1Ωh2Ωh3…ΩhmhT is a mh×1, here “mh” is the number of acceptable functions used to approximate deflection in a generalized coordinate vector.(13)A¯hiξh,ηh=Ahiξhahiηh,B¯hiξh,ηh=Bhiξhbhiηh,C¯hiξh,ηh=Chiξhchiηh,where,i=1,2,3,…mh.Equations (12) and (13) represented the potential and kinetic energies.(14a)TH=12Ω˙hTMhΩ˙h,(14b)VH=12ΩhTKhΩh.Hence, total energies of plate with hole are(15a)Ttotal=12Ω˙TMΩ˙−12Ω˙hTMhΩ˙h,(15b)Vtotal=12ΩTKΩ−12ΩhTKhΩh,where(16)Mh=ρhbc2M¯h,Kh=Dbc2K¯h,in which(17)M¯h=∫01∫01A¯hT⋅A¯h+βc122B¯hT⋅B¯h+C¯hT⋅C¯hdξhdηh,(18)K¯h=∫01∫011αc2∂B¯hT∂ξh⋅∂B¯h∂ξh+∂C¯hT∂ηh⋅∂C¯h∂ηh+1αc∂B¯hT∂ξh⋅∂C¯h∂ηh+∂C¯hT∂ηh⋅∂B¯h∂ξh+1−μ⋅1αc∂B¯hT∂ξh⋅∂C¯h∂ηh+∂C¯hT∂ηh⋅∂B¯h∂ξh−12∂B¯hT∂ηh⋅∂B¯h∂ηh+1αc2∂C¯hT∂ξh⋅∂C¯h∂ξh+1αc∂B¯hT∂ηh⋅∂C¯h∂ξh+∂C¯hT∂ξh⋅∂C¯h∂ηh+61−μκβc2B¯hT⋅B¯h+1αc2∂A¯hT∂ξh⋅∂A¯h∂ξh+1αcB¯hT⋅∂A¯h∂ξh+∂A¯hT∂ξhB¯h+C¯hT⋅C¯h+∂A¯hT∂ηh⋅∂A¯h∂ηh+C¯hT⋅∂A¯h∂ηh+∂A¯hT∂ηhC¯hdξhdηh,where thickness to width ratio βc=h/bc, aspect ratio αc=ac/bc, and non-dimensional variables ξh=xc/ac and ηh=yc/bc.Equations (17) and (18) have a distinctive advantage that the integral bounds range from 0 to 1, allowing for a closed-form solution. In terms of plate domain, the displacement matching requirement of the hole may be represented as follows:(19)whξh,ηh=wξ,η,ψhxξh,ηh=ψxξ,η,ψhyξh,ηh=ψyξ,η.The non-dimensional global and local coordinates have the following relationship:(20)ξ=rxa+acaξh,η=ryb+bcbηh.Considering equations (3), (10), (12), and (13) and by plugging them into equation (19), we can get(21)∑i=1mhA¯hiξh,ηhB¯hiξh,ηhC¯hiξh,ηhΩhit=∑i=1mhAhiξhahiηhBhiξhbhiηhChiξhchiηhΩhit=∑k=1mA¯kξ,ηB¯kξ,ηC¯kξ,ηΩkt=∑k=1mAkξakηBkξbkηCkξckηΩkt.Multiplying equation (21) by AhjξhahjηhBhjξhbhjηhChjξhchjηh and performing integration,(22)∑i=1mh∫01∫01AhiξhahiηhBhiξhbhiηhChiξhchiηhAhjξhahjηhBhjξhbhjηhChjξhchjηhdξhdηhΩhjt=∑k=1m∫01∫01AhiξhahiηhBhiξhbhiηhChiξhchiηhAkξakηBkξbkηCkξckηdξhdηhΩkt.Using orthogonal property, equation (22) may be expressed as eigenfunctions of beam(23)=∑k=1m∫01AhiξhAkξdξh∫01ahiηhakηdηh∫01BhiξhBkξdξh∫01bhiηhbkηdηh∫01ChiξhCkξdξh∫01chiηhckηdηhΩkt=∑k=1mThikΩkt,i=1,2,3,4…mh.If we write equation (23) as a matrix, we get(24)Ωh=ThiΩ,in which Th is the mh×m matrix for transforming coordinates. Inserting equation (24) in to equations (15a) and (15b), we can derive(25)Ttotal=12Ω˙TMΩ˙−12Ω˙TThTMhThΩ˙=12Ω˙TMrΩ˙,Vtotal=12ΩTKΩ−12ΩhTThTKhThΩh=12ΩTKrΩ,where Mr=M−ThTMhTh, Kr=K−ThTKhTh.In non-dimensional parameters, equations can be stated.(26a)Mr=ρhb2M¯r,Kr=Db2K¯r,whereM¯r=M¯−b¯c2ThTM¯hTh,(26b)K¯r=K¯−1b¯c2ThTK¯hTh,in which b¯ci=bci/bi=1,2,3,…n.Hence,(27)K¯r−ω¯2M¯rΛ=0.From equation (23), ‘Th’ is the transformation matrix because of the integral limits. In equation (25), the eigenvalue of the plate with cut-out can be derived. Equation (26) is used to find the natural frequencies.


### 1.2. Free Vibrational Exploration of an Elliptical Plate Using ICCM

In the previous section, ICCM was used to evaluate a rectangular plate with a rectangular cut, where the plate and hole geometry were identical. In this section, the natural frequencies of an elliptical plate with a cut-out are computed using ICCM, with the hole and plate geometries modified. An elliptical plate with a hole is assumed in the plane. The semi-major axis “a” and semi-minor axis “b” are considered in direction. The border of an elliptical plate is written as ξ2+η2=1.The strain energy and kinetic energy of an elliptical plate is(28)Vp=∫ξ=01∫η=01−ξ21α∂ψx∂ξ+∂ψy∂η2−21−μ1α∂ψx∂ξ∂ψy∂η−14∂ψx∂η+1α∂ψy∂ξ2+61−μκβ2ψx+1α∂w∂ξ2+ψy+∂w∂η2abdξdη,Tp=∫ξ−01∫η−01−ξ2w∗2+112β2ψx2+ψy2abdξdη.Deflection and rotation terms in the energy expressions (26a) and (26b) of an elliptical plate, to get the non-dimensional parameters, may be written in terms of acceptable functions that are introduced in equation (3) of an elliptical plate. Using relationship between local and global coordinates of hole and plate in equation (20), transformation matrix can be derived by applying an orthogonal property on eigenfunctions. Reduced non-dimensional and elliptical plate with hole stiffness matrices with hole are obtained from equation (25). Equations (26a) and (26b) give the non-dimensional eigenvalues of the plate with cut-out. Figure 2 shows the elliptical plate with a rectangular cut-out.Figure 2
Elliptical plate with a rectangular cut-out.Table1 shows the frequency parameter for a supported elliptical plate a/b = 2 and h/b = 0.2 of various sizes. Tables 2 and 3 show various sizes of h/b = 0.4 and h/b = 0.6.Table 1
The frequency parameter,ω¯=ωb2/π2ρh/D, for a simply supported elliptical plate a/b=2, h/b = 0.2 of various sizes, ac/a=bc/b=0,0.1,0.2,0.3,0.4,0.5,0.6.
ω¯00.10.20.30.40.50.6ω¯10.32530.31290.29590.28890.29490.31450.3399ω¯20.56840.56620.55240.53120.52290.48190.4272ω¯30.89970.89270.85390.66020.54590.54190.5966ω¯41.07841.04270.90190.89990.83950.76530.7087ω¯51.31241.30031.26521.24871.14741.10671.1334ω¯61.42241.40261.35121.29591.36541.31011.2237Table 2
The frequency parameter,ω¯=ωb2/π2ρh/D, for a simply supported elliptical plate a/b=2,h/b = 0.4 of various sizes, ac/a=bc/b=0,0.1,0.2,0.3,0.4,0.5,0.6.
ω¯00.10.20.30.40.50.6ω¯10.30430.29260.27610.26740.26770.27310.2636ω¯20.50960.50620.49020.46790.43770.36990.2972ω¯30.77240.76830.68830.53390.45620.46540.5002ω¯40.90990.86340.77140.75050.67860.60360.5544ω¯51.07761.06431.03871.01380.92160.85610.7968ω¯61.15311.13791.09871.06471.08650.98210.8469Table 3
The frequency parameter,ω¯=ωb2/π2ρh/D, for a simply supported elliptical plate a/b=2,h/b = 0.6 of various sizes, ac/a=bc/b=0,0.1,0.2,0.3,0.4,0.5,0.6.
ω¯00.10.20.30.40.50.6ω¯10.27850.26850.25310.24270.23710.22960.2017ω¯20.44730.44310.42690.40580.34700.28230.2155ω¯30.65280.65090.55410.42960.39350.39640.4160ω¯40.75620.70660.64950.61710.54640.48310.4504ω¯50.88020.86730.84990.81990.73450.65490.5630ω¯60.93190.92120.89090.87150.85600.73950.5994


### 1.3. Numerical Study

The ICCM is validated with the available literature for a rectangular plate with various aspect ratios and a thickness ofh/b = 0.2 in Table 4. The current ICCM results seem to match the available data with a high degree of precision.Table 4
Frequency parameterω¯=ωb2/π2ρh/D of various aspect ratios having a plate thickness of h/b = 0.2.
Aspect ratio (a/b)0.40.811.52Present method5.183182.201141.7681.31641.15218Ref. [1, 11]5.18312.20111.76791.31641.1521% error—————The results of a square plate with a cut-out of various sizes are validated with literature in Table5.Table 5
Frequency parameterω¯=ωb2/π2ρh/D of a square plate with a cut-out thickness of h = 0.2.
ac/a=bc/b=0.2ac/a=bc/b=0.4ac/a=bc/b=0.6a/3×a/3Present method1.7181.83552.3245137.428 (rad/sec)Ref. [18]1.7161.83452.323—Ref. [6]———5136.5 (rad/sec)% error0.11%0.054%0.043%0.02%From the Tables4 and 5, it is clear that the improved ICCM produces correct results when compared to techniques such as RRM, FEM, and FDM, all of which have a faster convergence time. As a result, the ICCM is used to find the natural frequencies of a thick rectangular plate with different aspect ratios and a cut-out elliptical plate by increasing the size of hole, in foregoing tables. Graphs are plotted in Figures 3 and 5 to show how these natural frequencies are varying with increase in size of cut-out for a rectangular plate and for an elliptical plate, respectively. Eigenvectors were also plotted to observe the deflection in Z-direction for each mode of vibration, and it is validated with the mode shape obtained from ANSYS in Figures 4 and 6.Figure 3
Aspect ratio variation of nondimensional eigenvalues. (a) Aspect ratioa/b = 0.4 and h/b = 0.2. (b) Aspect ratio a/b = 0.8 and h/b = 0.2. (c) Aspect ratio a/b = 1 and h/b = 0.2. (d) Aspect ratio a/b = 1.5 and h/b = 0.2. (e) Aspect ratio a/b = 2 and h/b = 0.2.
(a)(b)(c)(d)(e)Figure 4
Comparison of transverse displacement of various mode shapes of a rectangular plate (a = 2, b = 1, and h/b = 0.2) with a concentric cut-out (ac/a=bc/b=0.1) using MATLAB and ANSYS.
(a)(b)Figure 5
Variation of first six natural frequencies of an elliptical plate (a/b = 2) with respect to the size of cut-out. (a) Thickness (h/b = 0.2). (b) Thickness (h/b = 0.4). (c) Thickness (h/b = 0.6).
(a)(b)(c)Figure 6
Comparison of transverse displacement of various mode shapes of an elliptical plate (a = 2, b = 1, and h/b = 0.2) with a concentric cut-out (ac/a=bc/b=0.1) using MATLAB and ANSYS.Graphs were plotted for a numerical value in Tables6–10, as a function of cut-out size.Table 6
The frequency parameter,ω¯=ωb2/π2ρh/D, for a simply supported rectangular plate with a/b=0.4 and h/b = 0.2 with a concentric cut-out of various sizes, ac/a=bc/b=0,0.1,0.2,0.3,0.4,0.5,0.6.
ω¯00.10.20.30.40..50.6ω¯15.18315.02564.76754.57614.45984.35514.1113ω¯26.72126.67996.54286.36995.71654.96714.3116ω¯38.91388.89138.62996.79946.27056.32486.5839ω¯411.487511.33668.88088.64758.02467.43107.2052ω¯512.704011.586411.154511.410310.89549.67138.5177ω¯613.614513.546913.280512.305411.830910.7999.2534Table 7
The frequency parameter,ω¯=ωb2/π2ρh/D, for a simply supported rectangular plate with a/b=0.8 and h/b = 0.2 with a concentric cut-out of various sizes, ac/a=bc/b=0,0.1,0.2,0.3,0.4,0.5,0.6.
ω¯00.10.20.30.40.50.6ω¯12.20112.16882.13142.15172.25742.46692.7971ω¯24.21134.16003.96253.68823.49643.46143.4616ω¯35.18315.06624.60874.04333.65483.47493.6022ω¯46.72126.65386.53986.36305.72705.04724.6080ω¯56.87036.85046.87946.53536.15916.06696.2075ω¯68.83818.77158.39988.22818.03597.76687.4985Table 8
Frequency parameter for supported elliptical plate witha/b=1 and h/b = 0.2.
ω¯00.10.20.30.40..50.6ω¯11.76801.74461.71801.74021.83552.02222.3245ω¯23.86673.80893.56683.23362.99382.91182.9903ω¯33.86673.80933.56743.23422.99412.91212.9905ω¯45.58925.53175.43085.28814.98464.36693.9942ω¯56.60466.54186.37725.80405.13195.07825.2457ω¯66.60466.61076.92627.44887.16996.75936.5139Table 9
Frequency parameter for supported elliptical plate witha/b=1.5 and h/b = 0.2.
ω¯00.10.20.30.40.50.6ω¯11.31641.28921.25491.25771.31461.43881.6507ω¯22.36122.34412.26732.14962.06862.07642.1431ω¯33.51193.42853.05052.58842.28522.14702.2043ω¯43.86593.84653.88233.83633.51043.16202.9547ω¯54.34064.29914.22034.08073.90963.81683.8947ω¯65.58815.48955.22565.11615.02624.91384.8277Table 10
Frequency parameter for supported elliptical plate witha/b=2 and h/b = 0.2.
ω¯00.10.20.30.40.50.6ω¯11.15211.11671.06881.05281.08101.16221.3106ω¯21.76801.76061.72441.66741.63381.66041.6913ω¯32.70272.68482.70792.17361.86321.72081.7764ω¯43.38583.27002.71812.72952.61712.44522.3393ω¯53.86673.82223.71563.58973.35123.19343.1870ω¯63.86673.83783.76523.78434.04804.02533.8016Case 1.
Rectangular plate of different aspect ratios with concentric rectangular hole with increase in its size.Case 2.
An elliptical plate of aspect ratio (a/b = 2) of various thickness having a concentric rectangular hole with increase in its size.


## 1.1. Free Vibrational Exploration of Rectangle Plate with Cut-Out

Separate energy expressions are formulated for a hole domain and for a plate domain exploration of a plate with hole, and hence deflection matching condition is essential to establish a kinematic relation between these independent coordinate domains. It will also enable us to unify these two independent energies.


### 1.1.1. Expressions for Rectangular Plate

Consider a thick rectangular plate with a cut-out inx−y plane, having dimensions of a×b and ac×bc, respectively. Expressions for rectangular plate strain energy and kinetic energy from Reference [11] are as follows:(1)Vp=12∫vEh21−μ2∂ψx∂x+∂ψy∂y2−21−μ∂ψx∂x∂ψy∂y−14∂ψx∂y+∂ψy∂x2+κGψx+∂w∂x2+ψy+∂w∂y2dv,Tp=12ω2∫Aρhw2+112ρh3ψx2+ψy2dA,where h = plate thickness, ρ = mass density, ϑ=E = modulus of elasticity, w=transverse displacement, G = shear modulus, μ=Poisson’s ratio, ψx=−∂w/∂x = rotation about y direction, ψy=−∂w/∂y = rotation about x direction, κ = shear correction factor = 5/6, and ω = angular frequency.Introducing a non-dimensional parameter, thickness to width ratioβ=h/b, aspect ratio α=a/b, and non-dimensional variables ξ=x/a, η=y/b, the energy expression of plate is a non-dimensional representation of plate which can be written as follows:(2a)Vp=∫01∫011α∂ψx∂ξ+∂ψy∂η2−21−μ1α∂ψx∂ξ∂ψy∂η−14∂ψx∂η+1α∂ψy∂ξ2+61−μκβ2ψx+1α∂w∂ξ2+ψy+∂w∂η2abdξdη,(2b)Tp=∫01∫01w∗2+112β2ψx2+ψy2abdξdη,where w∗=w/b.Admissible functions are considered to expand transverse deflection and the rotations in equations (2a) and (2b) as(3)wξ,η,t=A¯iξ,ηΩt,ψxξ,η,t=B¯iξ,ηΩt,ψyξ,η,t=C¯iξ,ηΩt,where A¯iξ,η=A1A2A3A4…Am is a 1×m matrix and Ωiξ,η=Ω1Ω2Ω3Ω4…ΩmT is m×1 generalized coordinate vector, where “m” denotes admissible functions required for the approximating the deflection curve.Substituting equation (3) in equations (2a) and (2b) yields(4a)Tp=12Ω˙TMΩ˙,(4b)Vp=12ΩTKΩ,in which(5a)M=ρhb2M¯,(5b)K=Db2K¯,where(6a)M¯=∫01∫01A¯T.A¯+β122B¯T.B¯+C¯T.C¯dξdη,(6b)K¯=∫01∫011α2∂B¯T∂ξ⋅∂B¯∂ξ+∂C¯T∂η⋅∂C¯∂η+1α∂B¯T∂ξ⋅∂C¯∂η+∂C¯T∂η⋅∂B¯∂ξ+1−μ⋅1α∂B¯T∂ξ⋅∂C¯∂η+∂C¯T∂η⋅∂B¯∂ξ−12∂B¯T∂η⋅∂B¯∂η+1α2∂C¯T∂ξ⋅∂C¯∂ξ+1α∂B¯T∂η⋅∂C¯∂ξ+∂C¯T∂ξ⋅∂C¯∂η+61−μκβ2B¯T⋅B¯+1α2∂A¯T∂ξ⋅∂A¯∂ξ+1αB¯T⋅∂A¯∂ξ+∂A¯T∂ξB¯+C¯T⋅C¯+∂A¯T∂η⋅∂A¯∂η+C¯T⋅∂A¯∂η+∂A¯T∂ηC¯dξdη.In equations (6a) and (6b), stiffness matrices are represented by M¯, K¯, respectively, along with non-dimensional variables. Equation of motion of plate is obtained by substituting equation (3) in Lagrange’s expression, and its eigenvalue can be computed from(7)K−ω2MΛ=0.Using non-dimensional matrices introduced in equations (4a) and (4b), the above eigenvalue problem in non-dimensional form(8)K¯−ω¯2M¯Λ=0will give non-dimensional eigenvalues of a plate, in which “Λ” is an eigenvector. In equation (7),(9)ω¯=ωb2π2ρhDis known as the frequency parameter which represents the flexural rigidity in D=Eh3/121−μ2, which is known.To obtain in terms of admissible function matrices in theX and Y directions, an admissible function matrix from equation (3) has to be changed.(10)A¯iξ,η=Aiξaiη,B¯iξ,η=Biξbiη,C¯iξ,η=Ciξciη,where,i=1,2,3…m.


### 1.1.2. Boundary Conditions

A thick simply supported plate is considered and its eigenvalues are derived by considering an individual eigenfunction in each direction for each admissible function.(11)Aiξ=2siniπξ,aiη=2siniπη,Biξ=2cosiπξ,biη=2siniπη,Ciξ=2siniπξ,ciη=2cosiπη.Admissible functions in equation (11) are taken into account in such a way that they meet the assumed boundary conditions of the plate. All geometrical end points of ξ,η in equation (3) along with equation (11) will result in a zero deflection.


### 1.1.3. Energy Expressions for Rectangular Plate with Cut-Out Using ICCM

Kinematic relations are formulated for a plate having cut-out using ICCM; by matching the plate’s deflection to the hole’s deflection, these relationships may be used to produce updated mass and stiffness matrices as well as distinct plate and hole domains having their own coordinates. Fixing local coordinates to the hole domain is presented for rectangular hole in rectangular plate as illustrated in Figure1.Figure 1
Rectangular plate opening with cut-out.Taking non-dimensional coordinates into account, displacement and rotations are expressed asξh=xc/ac and ηh=yc/bc.(12)whξh,ηh,t=A¯hiξh,ηhΩt,ψhxξh,ηh,t=B¯hiξh,ηhΩt,ψhxξh,ηh,t=B¯hiξh,ηhΩt,where A¯hiξh,ηh=Ah1Ah2Ah3…Ahm is an 1×mh function matrix and Ωht=Ωh1Ωh2Ωh3…ΩhmhT is a mh×1, here “mh” is the number of acceptable functions used to approximate deflection in a generalized coordinate vector.(13)A¯hiξh,ηh=Ahiξhahiηh,B¯hiξh,ηh=Bhiξhbhiηh,C¯hiξh,ηh=Chiξhchiηh,where,i=1,2,3,…mh.Equations (12) and (13) represented the potential and kinetic energies.(14a)TH=12Ω˙hTMhΩ˙h,(14b)VH=12ΩhTKhΩh.Hence, total energies of plate with hole are(15a)Ttotal=12Ω˙TMΩ˙−12Ω˙hTMhΩ˙h,(15b)Vtotal=12ΩTKΩ−12ΩhTKhΩh,where(16)Mh=ρhbc2M¯h,Kh=Dbc2K¯h,in which(17)M¯h=∫01∫01A¯hT⋅A¯h+βc122B¯hT⋅B¯h+C¯hT⋅C¯hdξhdηh,(18)K¯h=∫01∫011αc2∂B¯hT∂ξh⋅∂B¯h∂ξh+∂C¯hT∂ηh⋅∂C¯h∂ηh+1αc∂B¯hT∂ξh⋅∂C¯h∂ηh+∂C¯hT∂ηh⋅∂B¯h∂ξh+1−μ⋅1αc∂B¯hT∂ξh⋅∂C¯h∂ηh+∂C¯hT∂ηh⋅∂B¯h∂ξh−12∂B¯hT∂ηh⋅∂B¯h∂ηh+1αc2∂C¯hT∂ξh⋅∂C¯h∂ξh+1αc∂B¯hT∂ηh⋅∂C¯h∂ξh+∂C¯hT∂ξh⋅∂C¯h∂ηh+61−μκβc2B¯hT⋅B¯h+1αc2∂A¯hT∂ξh⋅∂A¯h∂ξh+1αcB¯hT⋅∂A¯h∂ξh+∂A¯hT∂ξhB¯h+C¯hT⋅C¯h+∂A¯hT∂ηh⋅∂A¯h∂ηh+C¯hT⋅∂A¯h∂ηh+∂A¯hT∂ηhC¯hdξhdηh,where thickness to width ratio βc=h/bc, aspect ratio αc=ac/bc, and non-dimensional variables ξh=xc/ac and ηh=yc/bc.Equations (17) and (18) have a distinctive advantage that the integral bounds range from 0 to 1, allowing for a closed-form solution. In terms of plate domain, the displacement matching requirement of the hole may be represented as follows:(19)whξh,ηh=wξ,η,ψhxξh,ηh=ψxξ,η,ψhyξh,ηh=ψyξ,η.The non-dimensional global and local coordinates have the following relationship:(20)ξ=rxa+acaξh,η=ryb+bcbηh.Considering equations (3), (10), (12), and (13) and by plugging them into equation (19), we can get(21)∑i=1mhA¯hiξh,ηhB¯hiξh,ηhC¯hiξh,ηhΩhit=∑i=1mhAhiξhahiηhBhiξhbhiηhChiξhchiηhΩhit=∑k=1mA¯kξ,ηB¯kξ,ηC¯kξ,ηΩkt=∑k=1mAkξakηBkξbkηCkξckηΩkt.Multiplying equation (21) by AhjξhahjηhBhjξhbhjηhChjξhchjηh and performing integration,(22)∑i=1mh∫01∫01AhiξhahiηhBhiξhbhiηhChiξhchiηhAhjξhahjηhBhjξhbhjηhChjξhchjηhdξhdηhΩhjt=∑k=1m∫01∫01AhiξhahiηhBhiξhbhiηhChiξhchiηhAkξakηBkξbkηCkξckηdξhdηhΩkt.Using orthogonal property, equation (22) may be expressed as eigenfunctions of beam(23)=∑k=1m∫01AhiξhAkξdξh∫01ahiηhakηdηh∫01BhiξhBkξdξh∫01bhiηhbkηdηh∫01ChiξhCkξdξh∫01chiηhckηdηhΩkt=∑k=1mThikΩkt,i=1,2,3,4…mh.If we write equation (23) as a matrix, we get(24)Ωh=ThiΩ,in which Th is the mh×m matrix for transforming coordinates. Inserting equation (24) in to equations (15a) and (15b), we can derive(25)Ttotal=12Ω˙TMΩ˙−12Ω˙TThTMhThΩ˙=12Ω˙TMrΩ˙,Vtotal=12ΩTKΩ−12ΩhTThTKhThΩh=12ΩTKrΩ,where Mr=M−ThTMhTh, Kr=K−ThTKhTh.In non-dimensional parameters, equations can be stated.(26a)Mr=ρhb2M¯r,Kr=Db2K¯r,whereM¯r=M¯−b¯c2ThTM¯hTh,(26b)K¯r=K¯−1b¯c2ThTK¯hTh,in which b¯ci=bci/bi=1,2,3,…n.Hence,(27)K¯r−ω¯2M¯rΛ=0.From equation (23), ‘Th’ is the transformation matrix because of the integral limits. In equation (25), the eigenvalue of the plate with cut-out can be derived. Equation (26) is used to find the natural frequencies.


## 1.1.1. Expressions for Rectangular Plate

Consider a thick rectangular plate with a cut-out inx−y plane, having dimensions of a×b and ac×bc, respectively. Expressions for rectangular plate strain energy and kinetic energy from Reference [11] are as follows:(1)Vp=12∫vEh21−μ2∂ψx∂x+∂ψy∂y2−21−μ∂ψx∂x∂ψy∂y−14∂ψx∂y+∂ψy∂x2+κGψx+∂w∂x2+ψy+∂w∂y2dv,Tp=12ω2∫Aρhw2+112ρh3ψx2+ψy2dA,where h = plate thickness, ρ = mass density, ϑ=E = modulus of elasticity, w=transverse displacement, G = shear modulus, μ=Poisson’s ratio, ψx=−∂w/∂x = rotation about y direction, ψy=−∂w/∂y = rotation about x direction, κ = shear correction factor = 5/6, and ω = angular frequency.Introducing a non-dimensional parameter, thickness to width ratioβ=h/b, aspect ratio α=a/b, and non-dimensional variables ξ=x/a, η=y/b, the energy expression of plate is a non-dimensional representation of plate which can be written as follows:(2a)Vp=∫01∫011α∂ψx∂ξ+∂ψy∂η2−21−μ1α∂ψx∂ξ∂ψy∂η−14∂ψx∂η+1α∂ψy∂ξ2+61−μκβ2ψx+1α∂w∂ξ2+ψy+∂w∂η2abdξdη,(2b)Tp=∫01∫01w∗2+112β2ψx2+ψy2abdξdη,where w∗=w/b.Admissible functions are considered to expand transverse deflection and the rotations in equations (2a) and (2b) as(3)wξ,η,t=A¯iξ,ηΩt,ψxξ,η,t=B¯iξ,ηΩt,ψyξ,η,t=C¯iξ,ηΩt,where A¯iξ,η=A1A2A3A4…Am is a 1×m matrix and Ωiξ,η=Ω1Ω2Ω3Ω4…ΩmT is m×1 generalized coordinate vector, where “m” denotes admissible functions required for the approximating the deflection curve.Substituting equation (3) in equations (2a) and (2b) yields(4a)Tp=12Ω˙TMΩ˙,(4b)Vp=12ΩTKΩ,in which(5a)M=ρhb2M¯,(5b)K=Db2K¯,where(6a)M¯=∫01∫01A¯T.A¯+β122B¯T.B¯+C¯T.C¯dξdη,(6b)K¯=∫01∫011α2∂B¯T∂ξ⋅∂B¯∂ξ+∂C¯T∂η⋅∂C¯∂η+1α∂B¯T∂ξ⋅∂C¯∂η+∂C¯T∂η⋅∂B¯∂ξ+1−μ⋅1α∂B¯T∂ξ⋅∂C¯∂η+∂C¯T∂η⋅∂B¯∂ξ−12∂B¯T∂η⋅∂B¯∂η+1α2∂C¯T∂ξ⋅∂C¯∂ξ+1α∂B¯T∂η⋅∂C¯∂ξ+∂C¯T∂ξ⋅∂C¯∂η+61−μκβ2B¯T⋅B¯+1α2∂A¯T∂ξ⋅∂A¯∂ξ+1αB¯T⋅∂A¯∂ξ+∂A¯T∂ξB¯+C¯T⋅C¯+∂A¯T∂η⋅∂A¯∂η+C¯T⋅∂A¯∂η+∂A¯T∂ηC¯dξdη.In equations (6a) and (6b), stiffness matrices are represented by M¯, K¯, respectively, along with non-dimensional variables. Equation of motion of plate is obtained by substituting equation (3) in Lagrange’s expression, and its eigenvalue can be computed from(7)K−ω2MΛ=0.Using non-dimensional matrices introduced in equations (4a) and (4b), the above eigenvalue problem in non-dimensional form(8)K¯−ω¯2M¯Λ=0will give non-dimensional eigenvalues of a plate, in which “Λ” is an eigenvector. In equation (7),(9)ω¯=ωb2π2ρhDis known as the frequency parameter which represents the flexural rigidity in D=Eh3/121−μ2, which is known.To obtain in terms of admissible function matrices in theX and Y directions, an admissible function matrix from equation (3) has to be changed.(10)A¯iξ,η=Aiξaiη,B¯iξ,η=Biξbiη,C¯iξ,η=Ciξciη,where,i=1,2,3…m.


## 1.1.2. Boundary Conditions

A thick simply supported plate is considered and its eigenvalues are derived by considering an individual eigenfunction in each direction for each admissible function.(11)Aiξ=2siniπξ,aiη=2siniπη,Biξ=2cosiπξ,biη=2siniπη,Ciξ=2siniπξ,ciη=2cosiπη.Admissible functions in equation (11) are taken into account in such a way that they meet the assumed boundary conditions of the plate. All geometrical end points of ξ,η in equation (3) along with equation (11) will result in a zero deflection.


## 1.1.3. Energy Expressions for Rectangular Plate with Cut-Out Using ICCM

Kinematic relations are formulated for a plate having cut-out using ICCM; by matching the plate’s deflection to the hole’s deflection, these relationships may be used to produce updated mass and stiffness matrices as well as distinct plate and hole domains having their own coordinates. Fixing local coordinates to the hole domain is presented for rectangular hole in rectangular plate as illustrated in Figure1.Figure 1
Rectangular plate opening with cut-out.Taking non-dimensional coordinates into account, displacement and rotations are expressed asξh=xc/ac and ηh=yc/bc.(12)whξh,ηh,t=A¯hiξh,ηhΩt,ψhxξh,ηh,t=B¯hiξh,ηhΩt,ψhxξh,ηh,t=B¯hiξh,ηhΩt,where A¯hiξh,ηh=Ah1Ah2Ah3…Ahm is an 1×mh function matrix and Ωht=Ωh1Ωh2Ωh3…ΩhmhT is a mh×1, here “mh” is the number of acceptable functions used to approximate deflection in a generalized coordinate vector.(13)A¯hiξh,ηh=Ahiξhahiηh,B¯hiξh,ηh=Bhiξhbhiηh,C¯hiξh,ηh=Chiξhchiηh,where,i=1,2,3,…mh.Equations (12) and (13) represented the potential and kinetic energies.(14a)TH=12Ω˙hTMhΩ˙h,(14b)VH=12ΩhTKhΩh.Hence, total energies of plate with hole are(15a)Ttotal=12Ω˙TMΩ˙−12Ω˙hTMhΩ˙h,(15b)Vtotal=12ΩTKΩ−12ΩhTKhΩh,where(16)Mh=ρhbc2M¯h,Kh=Dbc2K¯h,in which(17)M¯h=∫01∫01A¯hT⋅A¯h+βc122B¯hT⋅B¯h+C¯hT⋅C¯hdξhdηh,(18)K¯h=∫01∫011αc2∂B¯hT∂ξh⋅∂B¯h∂ξh+∂C¯hT∂ηh⋅∂C¯h∂ηh+1αc∂B¯hT∂ξh⋅∂C¯h∂ηh+∂C¯hT∂ηh⋅∂B¯h∂ξh+1−μ⋅1αc∂B¯hT∂ξh⋅∂C¯h∂ηh+∂C¯hT∂ηh⋅∂B¯h∂ξh−12∂B¯hT∂ηh⋅∂B¯h∂ηh+1αc2∂C¯hT∂ξh⋅∂C¯h∂ξh+1αc∂B¯hT∂ηh⋅∂C¯h∂ξh+∂C¯hT∂ξh⋅∂C¯h∂ηh+61−μκβc2B¯hT⋅B¯h+1αc2∂A¯hT∂ξh⋅∂A¯h∂ξh+1αcB¯hT⋅∂A¯h∂ξh+∂A¯hT∂ξhB¯h+C¯hT⋅C¯h+∂A¯hT∂ηh⋅∂A¯h∂ηh+C¯hT⋅∂A¯h∂ηh+∂A¯hT∂ηhC¯hdξhdηh,where thickness to width ratio βc=h/bc, aspect ratio αc=ac/bc, and non-dimensional variables ξh=xc/ac and ηh=yc/bc.Equations (17) and (18) have a distinctive advantage that the integral bounds range from 0 to 1, allowing for a closed-form solution. In terms of plate domain, the displacement matching requirement of the hole may be represented as follows:(19)whξh,ηh=wξ,η,ψhxξh,ηh=ψxξ,η,ψhyξh,ηh=ψyξ,η.The non-dimensional global and local coordinates have the following relationship:(20)ξ=rxa+acaξh,η=ryb+bcbηh.Considering equations (3), (10), (12), and (13) and by plugging them into equation (19), we can get(21)∑i=1mhA¯hiξh,ηhB¯hiξh,ηhC¯hiξh,ηhΩhit=∑i=1mhAhiξhahiηhBhiξhbhiηhChiξhchiηhΩhit=∑k=1mA¯kξ,ηB¯kξ,ηC¯kξ,ηΩkt=∑k=1mAkξakηBkξbkηCkξckηΩkt.Multiplying equation (21) by AhjξhahjηhBhjξhbhjηhChjξhchjηh and performing integration,(22)∑i=1mh∫01∫01AhiξhahiηhBhiξhbhiηhChiξhchiηhAhjξhahjηhBhjξhbhjηhChjξhchjηhdξhdηhΩhjt=∑k=1m∫01∫01AhiξhahiηhBhiξhbhiηhChiξhchiηhAkξakηBkξbkηCkξckηdξhdηhΩkt.Using orthogonal property, equation (22) may be expressed as eigenfunctions of beam(23)=∑k=1m∫01AhiξhAkξdξh∫01ahiηhakηdηh∫01BhiξhBkξdξh∫01bhiηhbkηdηh∫01ChiξhCkξdξh∫01chiηhckηdηhΩkt=∑k=1mThikΩkt,i=1,2,3,4…mh.If we write equation (23) as a matrix, we get(24)Ωh=ThiΩ,in which Th is the mh×m matrix for transforming coordinates. Inserting equation (24) in to equations (15a) and (15b), we can derive(25)Ttotal=12Ω˙TMΩ˙−12Ω˙TThTMhThΩ˙=12Ω˙TMrΩ˙,Vtotal=12ΩTKΩ−12ΩhTThTKhThΩh=12ΩTKrΩ,where Mr=M−ThTMhTh, Kr=K−ThTKhTh.In non-dimensional parameters, equations can be stated.(26a)Mr=ρhb2M¯r,Kr=Db2K¯r,whereM¯r=M¯−b¯c2ThTM¯hTh,(26b)K¯r=K¯−1b¯c2ThTK¯hTh,in which b¯ci=bci/bi=1,2,3,…n.Hence,(27)K¯r−ω¯2M¯rΛ=0.From equation (23), ‘Th’ is the transformation matrix because of the integral limits. In equation (25), the eigenvalue of the plate with cut-out can be derived. Equation (26) is used to find the natural frequencies.


## 1.2. Free Vibrational Exploration of an Elliptical Plate Using ICCM

In the previous section, ICCM was used to evaluate a rectangular plate with a rectangular cut, where the plate and hole geometry were identical. In this section, the natural frequencies of an elliptical plate with a cut-out are computed using ICCM, with the hole and plate geometries modified. An elliptical plate with a hole is assumed in the plane. The semi-major axis “a” and semi-minor axis “b” are considered in direction. The border of an elliptical plate is written as ξ2+η2=1.The strain energy and kinetic energy of an elliptical plate is(28)Vp=∫ξ=01∫η=01−ξ21α∂ψx∂ξ+∂ψy∂η2−21−μ1α∂ψx∂ξ∂ψy∂η−14∂ψx∂η+1α∂ψy∂ξ2+61−μκβ2ψx+1α∂w∂ξ2+ψy+∂w∂η2abdξdη,Tp=∫ξ−01∫η−01−ξ2w∗2+112β2ψx2+ψy2abdξdη.Deflection and rotation terms in the energy expressions (26a) and (26b) of an elliptical plate, to get the non-dimensional parameters, may be written in terms of acceptable functions that are introduced in equation (3) of an elliptical plate. Using relationship between local and global coordinates of hole and plate in equation (20), transformation matrix can be derived by applying an orthogonal property on eigenfunctions. Reduced non-dimensional and elliptical plate with hole stiffness matrices with hole are obtained from equation (25). Equations (26a) and (26b) give the non-dimensional eigenvalues of the plate with cut-out. Figure 2 shows the elliptical plate with a rectangular cut-out.Figure 2
Elliptical plate with a rectangular cut-out.Table1 shows the frequency parameter for a supported elliptical plate a/b = 2 and h/b = 0.2 of various sizes. Tables 2 and 3 show various sizes of h/b = 0.4 and h/b = 0.6.Table 1
The frequency parameter,ω¯=ωb2/π2ρh/D, for a simply supported elliptical plate a/b=2, h/b = 0.2 of various sizes, ac/a=bc/b=0,0.1,0.2,0.3,0.4,0.5,0.6.
ω¯00.10.20.30.40.50.6ω¯10.32530.31290.29590.28890.29490.31450.3399ω¯20.56840.56620.55240.53120.52290.48190.4272ω¯30.89970.89270.85390.66020.54590.54190.5966ω¯41.07841.04270.90190.89990.83950.76530.7087ω¯51.31241.30031.26521.24871.14741.10671.1334ω¯61.42241.40261.35121.29591.36541.31011.2237Table 2
The frequency parameter,ω¯=ωb2/π2ρh/D, for a simply supported elliptical plate a/b=2,h/b = 0.4 of various sizes, ac/a=bc/b=0,0.1,0.2,0.3,0.4,0.5,0.6.
ω¯00.10.20.30.40.50.6ω¯10.30430.29260.27610.26740.26770.27310.2636ω¯20.50960.50620.49020.46790.43770.36990.2972ω¯30.77240.76830.68830.53390.45620.46540.5002ω¯40.90990.86340.77140.75050.67860.60360.5544ω¯51.07761.06431.03871.01380.92160.85610.7968ω¯61.15311.13791.09871.06471.08650.98210.8469Table 3
The frequency parameter,ω¯=ωb2/π2ρh/D, for a simply supported elliptical plate a/b=2,h/b = 0.6 of various sizes, ac/a=bc/b=0,0.1,0.2,0.3,0.4,0.5,0.6.
ω¯00.10.20.30.40.50.6ω¯10.27850.26850.25310.24270.23710.22960.2017ω¯20.44730.44310.42690.40580.34700.28230.2155ω¯30.65280.65090.55410.42960.39350.39640.4160ω¯40.75620.70660.64950.61710.54640.48310.4504ω¯50.88020.86730.84990.81990.73450.65490.5630ω¯60.93190.92120.89090.87150.85600.73950.5994


## 1.3. Numerical Study

The ICCM is validated with the available literature for a rectangular plate with various aspect ratios and a thickness ofh/b = 0.2 in Table 4. The current ICCM results seem to match the available data with a high degree of precision.Table 4
Frequency parameterω¯=ωb2/π2ρh/D of various aspect ratios having a plate thickness of h/b = 0.2.
Aspect ratio (a/b)0.40.811.52Present method5.183182.201141.7681.31641.15218Ref. [1, 11]5.18312.20111.76791.31641.1521% error—————The results of a square plate with a cut-out of various sizes are validated with literature in Table5.Table 5
Frequency parameterω¯=ωb2/π2ρh/D of a square plate with a cut-out thickness of h = 0.2.
ac/a=bc/b=0.2ac/a=bc/b=0.4ac/a=bc/b=0.6a/3×a/3Present method1.7181.83552.3245137.428 (rad/sec)Ref. [18]1.7161.83452.323—Ref. [6]———5136.5 (rad/sec)% error0.11%0.054%0.043%0.02%From the Tables4 and 5, it is clear that the improved ICCM produces correct results when compared to techniques such as RRM, FEM, and FDM, all of which have a faster convergence time. As a result, the ICCM is used to find the natural frequencies of a thick rectangular plate with different aspect ratios and a cut-out elliptical plate by increasing the size of hole, in foregoing tables. Graphs are plotted in Figures 3 and 5 to show how these natural frequencies are varying with increase in size of cut-out for a rectangular plate and for an elliptical plate, respectively. Eigenvectors were also plotted to observe the deflection in Z-direction for each mode of vibration, and it is validated with the mode shape obtained from ANSYS in Figures 4 and 6.Figure 3
Aspect ratio variation of nondimensional eigenvalues. (a) Aspect ratioa/b = 0.4 and h/b = 0.2. (b) Aspect ratio a/b = 0.8 and h/b = 0.2. (c) Aspect ratio a/b = 1 and h/b = 0.2. (d) Aspect ratio a/b = 1.5 and h/b = 0.2. (e) Aspect ratio a/b = 2 and h/b = 0.2.
(a)(b)(c)(d)(e)Figure 4
Comparison of transverse displacement of various mode shapes of a rectangular plate (a = 2, b = 1, and h/b = 0.2) with a concentric cut-out (ac/a=bc/b=0.1) using MATLAB and ANSYS.
(a)(b)Figure 5
Variation of first six natural frequencies of an elliptical plate (a/b = 2) with respect to the size of cut-out. (a) Thickness (h/b = 0.2). (b) Thickness (h/b = 0.4). (c) Thickness (h/b = 0.6).
(a)(b)(c)Figure 6
Comparison of transverse displacement of various mode shapes of an elliptical plate (a = 2, b = 1, and h/b = 0.2) with a concentric cut-out (ac/a=bc/b=0.1) using MATLAB and ANSYS.Graphs were plotted for a numerical value in Tables6–10, as a function of cut-out size.Table 6
The frequency parameter,ω¯=ωb2/π2ρh/D, for a simply supported rectangular plate with a/b=0.4 and h/b = 0.2 with a concentric cut-out of various sizes, ac/a=bc/b=0,0.1,0.2,0.3,0.4,0.5,0.6.
ω¯00.10.20.30.40..50.6ω¯15.18315.02564.76754.57614.45984.35514.1113ω¯26.72126.67996.54286.36995.71654.96714.3116ω¯38.91388.89138.62996.79946.27056.32486.5839ω¯411.487511.33668.88088.64758.02467.43107.2052ω¯512.704011.586411.154511.410310.89549.67138.5177ω¯613.614513.546913.280512.305411.830910.7999.2534Table 7
The frequency parameter,ω¯=ωb2/π2ρh/D, for a simply supported rectangular plate with a/b=0.8 and h/b = 0.2 with a concentric cut-out of various sizes, ac/a=bc/b=0,0.1,0.2,0.3,0.4,0.5,0.6.
ω¯00.10.20.30.40.50.6ω¯12.20112.16882.13142.15172.25742.46692.7971ω¯24.21134.16003.96253.68823.49643.46143.4616ω¯35.18315.06624.60874.04333.65483.47493.6022ω¯46.72126.65386.53986.36305.72705.04724.6080ω¯56.87036.85046.87946.53536.15916.06696.2075ω¯68.83818.77158.39988.22818.03597.76687.4985Table 8
Frequency parameter for supported elliptical plate witha/b=1 and h/b = 0.2.
ω¯00.10.20.30.40..50.6ω¯11.76801.74461.71801.74021.83552.02222.3245ω¯23.86673.80893.56683.23362.99382.91182.9903ω¯33.86673.80933.56743.23422.99412.91212.9905ω¯45.58925.53175.43085.28814.98464.36693.9942ω¯56.60466.54186.37725.80405.13195.07825.2457ω¯66.60466.61076.92627.44887.16996.75936.5139Table 9
Frequency parameter for supported elliptical plate witha/b=1.5 and h/b = 0.2.
ω¯00.10.20.30.40.50.6ω¯11.31641.28921.25491.25771.31461.43881.6507ω¯22.36122.34412.26732.14962.06862.07642.1431ω¯33.51193.42853.05052.58842.28522.14702.2043ω¯43.86593.84653.88233.83633.51043.16202.9547ω¯54.34064.29914.22034.08073.90963.81683.8947ω¯65.58815.48955.22565.11615.02624.91384.8277Table 10
Frequency parameter for supported elliptical plate witha/b=2 and h/b = 0.2.
ω¯00.10.20.30.40.50.6ω¯11.15211.11671.06881.05281.08101.16221.3106ω¯21.76801.76061.72441.66741.63381.66041.6913ω¯32.70272.68482.70792.17361.86321.72081.7764ω¯43.38583.27002.71812.72952.61712.44522.3393ω¯53.86673.82223.71563.58973.35123.19343.1870ω¯63.86673.83783.76523.78434.04804.02533.8016Case 1.
Rectangular plate of different aspect ratios with concentric rectangular hole with increase in its size.Case 2.
An elliptical plate of aspect ratio (a/b = 2) of various thickness having a concentric rectangular hole with increase in its size.


## 2. Results and Discussion

This paper delves at the first six natural frequencies of a thick, rectangular, and elliptical plate. A cut-out is derived using ICCM. Based on the various combinations of aspect ratios and cut-out sizes, the derived numerical results were plotted in Figures3 and 5. The frequency occurs due to a large hole, and maximum variation was observed compared to a plate without a cut. These observations were tabulated as follows so that the designers can change the natural frequency of a structure up to maximum limit of scope simply by introducing a cut-out. Figure 6 shows the comparison of transverse displacement of various mode shapes of an elliptical plate.


### 2.1. Rectangular Plate with Cut-Out

The eigenvalue of a rectangular plate with a hole is lower when compared to a rectangular plate without a cut-out, as seen in the images. For an aspect ratio ofa/b = 0.4 and h/b = 0.2, the lowest natural frequency is obtained with a cut-out of size ac/a=bc/b=0.6. The cut-out frequency is reduced by 20.6%. Similarly, rectangular plates with cut-out frequency are increased. The combination of aspect ratios and cut-outs are as shown in Table 11.Table 11
Various aspect ratios of plate.
Aspect ratioa/bh/bCut-out size at which lowestω is obtained ac/a=bc/b=% ofω reduced, when compared with a plate without cut-out0.40.20.620.60.80.20.23.1610.20.22.821.50.20.24.6720.20.38.61


### 2.2. Elliptical Plate with Cut-Out

For an aspect ratio ofa/b = 2 and h/b = 0.2, the lowest natural frequency is obtained with a cut-out of size ac/a=bc/b=0.3. The frequency is reduced by 11.18% when compared to other combinations of thickness ratios and cut-outs. The lowest natural frequency is obtained as shown in Table 12.Table 12
Various thickness ratios of plate.
Aspect ratioa/bh/bCut-out size at which lowestω is obtained ac/a=bc/b% ofω reduced, when compared with a plate without cut-out20.20.311.1820.40.613.3720.60.627.57This information is useful to the designers to vary the natural frequency of a structure, up to a maximum limit of scope simply by introducing a cut-out.


## 2.1. Rectangular Plate with Cut-Out

The eigenvalue of a rectangular plate with a hole is lower when compared to a rectangular plate without a cut-out, as seen in the images. For an aspect ratio ofa/b = 0.4 and h/b = 0.2, the lowest natural frequency is obtained with a cut-out of size ac/a=bc/b=0.6. The cut-out frequency is reduced by 20.6%. Similarly, rectangular plates with cut-out frequency are increased. The combination of aspect ratios and cut-outs are as shown in Table 11.Table 11
Various aspect ratios of plate.
Aspect ratioa/bh/bCut-out size at which lowestω is obtained ac/a=bc/b=% ofω reduced, when compared with a plate without cut-out0.40.20.620.60.80.20.23.1610.20.22.821.50.20.24.6720.20.38.61


## 2.2. Elliptical Plate with Cut-Out

For an aspect ratio ofa/b = 2 and h/b = 0.2, the lowest natural frequency is obtained with a cut-out of size ac/a=bc/b=0.3. The frequency is reduced by 11.18% when compared to other combinations of thickness ratios and cut-outs. The lowest natural frequency is obtained as shown in Table 12.Table 12
Various thickness ratios of plate.
Aspect ratioa/bh/bCut-out size at which lowestω is obtained ac/a=bc/b% ofω reduced, when compared with a plate without cut-out20.20.311.1820.40.613.3720.60.627.57This information is useful to the designers to vary the natural frequency of a structure, up to a maximum limit of scope simply by introducing a cut-out.


## 3. Conclusions

Numerical results, plotted as curves in Figure3, reveal that as aspect ratio of the plate is increased, natural frequencies were decreased due to increase in mass without a significant increase in stiffness and the behavior of fundamental natural frequency is consistent with the variation of size of cut-out. It is also observed that in every aspect ratio, 2nd and 3rd mode of frequencies can be made same by including an appropriate size of cut-out. In the same manner, 4th and 5th mode of frequencies can also be made same in every aspect ratio of plate. Aspect ratio a/b = 0.8 with h/b = 0.2 has a high-frequency range that significantly increased the mass compared to aspect ratio 1. Due to the change in stiffness, the variation of frequencies exists that is smooth in manner. From Figure 5, it is worth noting that increasing the thickness ratio has no discernible influence on the plate’s basic natural frequency. In every thickness ratio, mode 3 natural frequency can be made equal to mode 2 natural frequency by introducing a cut-out of size nearer to 0.4, and hence by introducing a cut-out, structure can be made to vibrate in different modes for the same value of natural frequency. From the analysis, for the study of a thick plate with cut-outs, modified ICCM is an excellent numerical tool that converges to the solution faster than other traditional approaches.



---

*Source: 7212075-2021-10-16.xml*