Approximate natural vibration analysis of rectangular plates with openings using assumed mode method

International Journal of Naval Architecture and Ocean Engineering.
2013.
Sep,
5(3):
478-491

- Published : September 30, 2013

Download

PDF

e-PUB

PubReader

PPT

Export by style

Share

Article

Metrics

Cited by

TagCloud

Natural vibration analysis of plates with openings of different shape represents an important issue in naval architecture and ocean engineering applications. In this paper, a procedure for vibration analysis of plates with openings and arbitrary edge constraints is presented. It is based on the assumed mode method, where natural frequencies and modes are determined by solving an eigenvalue problem of a multi-degree-of-freedom system matrix equation derived by using Lagrange's equations of motion. The presented solution represents an extension of a procedure for natural vibration analysis of rectangular plates without openings, which has been recently presented in the literature. The effect of an opening is taken into account in an intuitive way, i.e. by subtracting its energy from the total plate energy without opening. Illustrative numerical examples include dynamic analysis of rectangular plates with rectangular, elliptic, circular as well as oval openings with various plate thicknesses and different combinations of boundary conditions. The results are compared with those obtained by the finite element method (FEM) as well as those available in the relevant literature, and very good agreement is achieved.
w
, and angles of cross-section rotation about
x
and
y
axis,
ψ_{x}
and
ψ_{y}
, respectively. From the equilibrium of sectional forces (bending moments, torsional moments and shear forces) and inertia forces, equations of motions are derived:
where
ρ
is plate density,
h
is plate thickness,
k
is shear coefficient, while
ν
is Poisson’s ratio. Further on,
D
represents plate flexural rigidity
D
=
Eh
^{3}
/ (12(1−
ν
^{2}
)) , while
E
and
G
=
E
/(2(1+
ν
)) are Young’s and shear modulus, respectively.
The basic idea to account for an opening effect is to subtract the energy of the opening part from the total plate energy without openings. In that sense, following expressions for the strain and kinetic energy of rectangular plate with openings, respectively, are valid:
where
V_{p}
and
V_{o}
are strain energies and
T_{p}
and
T_{o}
are kinetic energies of plate and opening, respectively.
ξ
=
x
/
a
,
η
=
y
/
b
,
α
=
a
/
b
and
S
=
kGh
/
D
for a rectangular plate of length
a
and breadth
b
, one can write for the strain and kinetic energy of rectangular plate without openings,
V
and
T
, respectively (Kim et al., 2012):
where (
K
_{Tx1}
) = (
K
_{Tx1}
a
/
kGh
), (
K
_{Tx2}
) = (
K
_{Tx2}
a
/
kGh
), (
K
_{Ty1}
) = (
K
_{Ty1}
b
/
kGh
) and (
K
_{Ty2}
) = (
K
_{Ty2}
b
/
kGh
) re non-dimensional stiffness at
x
= 0,
x
=
a
,
y
= 0 and
y
=
b
, respectively, and correspond to the translational spring constants per unit length
K
_{Tx1}
,
K
_{Tx2}
,
K
_{Ty1}
and
K
_{Ty2}
shown in
Fig. 1
. In the same manner, (
K
_{Rx1}
) = (
K
_{Rx1}
a
/
D
), (
K
_{Rx2}
) = (
K
_{Rx2}
a
/
D
), (
K
_{Ry1}
) = (
K
_{Ry1}
b
/
D
) and (
K
_{Ry2}
) = (
K
_{Ry2}
b
/
D
), are for the rotational spring constants per unit length
K
_{Rx1}
,
K
_{Rx2}
,
K
_{Ry1}
and
K
_{Ry2}
shown in
Fig. 1
, respectively. In
Fig. 1
the following quantities are related to
x
direction:
a
,
K
_{Tx1}
,
K
_{Tx2}
,
K
_{Rx1}
and
K
_{Rx2}
. Other quantities (
b
,
K
_{Ty1}
,
K
_{Ty2}
K
_{Ry1}
and
K
_{Ry2}
) are relevant for
y
direction.
Rectangular plate elastically supported in translational and rotational directions.
In
Fig. 2
x_{ro}
and
y_{ro}
denote longitudinal and transverse coordinate of rectangular opening centre of gravity and
a_{ro}
and
b_{ro}
are one half of its length and breadth.
Rectangular plate with rectangular opening.
Rectangular plate with elliptic opening.
If one introduces
x_{eo}
=
y_{eo}
in the Eq. (10) and (11), strain and kinetic energies of circular opening is obtained. Similarly as in previous case, in
Fig. 3
x_{eo}
and
y_{eo}
denote longitudinal and transverse coordinate of elliptic opening centre of gravity, while and
a_{eo}
and
b_{eo}
are opening major and minor semi-axis, respectively.
x_{oo}
and
yo_{eo}
are longitudinal and transverse coordinate of oval opening centre of gravity, while and
a_{oo}
and
b_{oo}
are equal to one half of its length and breadth, respectively.
In that way its strain and kinetic energies can be written in following form:
and can be calculated by above listed formulae.
Rectangular plate with oval opening.
If one wants to apply developed procedure for more than one opening, opening strain and kinetic energies have to be calculated separately, and subtracted from the corresponding energies of plate without openings.
where
X_{m}
(
ξ
),
Y_{n}
(
η
),
Θ_{m}
(
ξ
) and
Φ_{n}
(
η
) are the orthogonal polynomials satisfying the specified elastic edge constraints with respect to
ξ
and
η
. Furthermore,
a_{mn}
(
t
),
b_{mn}
(
t
) and
c_{mn}
(
t
) are the influence coefficients of orthogonal polynomials. Also,
M
and
N
are the number of orthogonal polynomials used for approximate solution in
ξ
and
η
direction, respectively. Eqs. (14), (15) and (16) can be alternatively written in matrix form:
where
By substituting the Eq. (4) and (5) into Lagrange's equation of motion below
the discrete matrix equation with 3x
M
x
N
degrees of freedom can be obtained as the following equation:
where [M] and [K] are the mass and the stiffness matrix, respectively, and their constitution is described in details in Kim et al. (2012).
If we assume harmonic vibrations, i.e.
Eq. (22) leads to an eigenvalue problem which gives natural frequencies and eigenvectors of the system, where
j
is the unit imaginary and
ω
represents the angular frequency. The mode shape corresponding to each natural frequency is obtained from the following equation:
where
l
represents the order of mode. It should be noted here that when using the assumed mode method, the orthogonal polynomials corresponding to the property of target model should be applied to achieve accurate analysis. For this purpose, the characteristic orthogonal polynomials having the property of Timoshenko beam functions which satisfies the specified edge constraints are used, and their complete derivation is presented by Kim et al. (2012).
In addition, it should be mentioned that developed procedure is not fully physical consistent, since in the Rayleigh-Ritz method all assumed shape function should satisfy at least geometric boundary conditions (Szilard, 2004), which is not unconditionally the case in the presented procedure at the opening boundaries. Namely, it is necessary to say that there exists differrence in real boundary conditions at the free opening edge and remained after excluding that part. In the former case boundary forces are zero, while in the latter case they are equal to the sectional forces. However, if opening dimensions are not too large comparing the plate dimensions, the procedure should normally give quite accurate and reliable results.
ξ
and
η
directions, respectively, after convergence test. The values of Young's modulus, material density and Poisson's ratio are the same in all numerical examples and their values are set to 2.1×10
^{11}
N
/
m^{2}
, 7850
kg
/
m^{3}
and 0.3, respectively. Furthermore, the value of shear factor
k
is adjusted to 0.86667 in all calculations. In all numerical examples, the results obtained by the proposed procedure are denoted with PS (Present Solution).
λ
is presented for two different relative thickness ratios
h
/
b
, which include thin and moderately thick plate. Non-dimensional frequency parameter is calculated according to following formula:
where the value of natural frequency
ω
is given in
rad/s
. In this case the obtained results are compared with recently presented analytical solutions (Xing and Liu, 2009; Senjanović et al., 2013), as well as with numerical results obtained by finite element method (MSC, 2010). Generally, it is obvious from the shown results,
Table 1
, that the developed procedure provides very good results comparing to the analytical solution.
Frequency parameter λ of rectangular plate without openings, SCSC edge constraints.
m
and 0.7
m
, respectively, is considered. Dimensions of the central opening are 0.2×0.1
m
, and different relative thickness ratios
h
/
a
are taken into account. Values of the first 3 natural frequencies, calculated according to the developed procedure together with finite element solutions obtained by MSC/NASTRAN (MSC, 2010), are listed in
Table 2
. Different combinations of edge constraints, i.e. SSSS (simple-simple-simple-simple), CCCC (clamped- clamped-clamped-clamped) and FCSC (free-clamped-simple-clamped), have been analysed. The agreement of the natural frequencies is quite good, especially bearing in mind the fact that dimensions of an opening are quite large comparing to plate dimensions. Natural modes for some selected cases are presented in
Figs. 5
and
6
, where also good agreement is obvious.
Natural frequency f (Hz ) of rectangular plate with central rectangular opening.
Mode shapes of rectangular plate with centralrectangular opening, h = 0.01m , SSSS.
Mode shapes of rectangular plate with centralrectangular opening, h = 0.01m , FCSC.
a_{eo}
=
m
and 0.1
b_{eo}
=
m
, while in the latter one 0.08
a_{eo}
=
m
and 0.15
b_{eo}
=
m
. Due to reason of simplicity, only FSCS case is analysed. In this case also reasonable agreement between both, natural frequencies listed in
Table 3
, and natural modes shown in
Figs. 7
and
8
, is achieved.
Natural frequencies f (Hz ) of rectangular plate with elliptic openings.
Mode shapes of rectangular plate with centralelliptic opening, h = 0.05m , FSCS.
Mode shapes of rectangular plate with twoelliptic openings, h = 0.05m , FSCS.
R
represents a special case of a plate with elliptic opening, and that means the above formulae for case of elliptic opening are valid if
a_{eo}
=
b_{bo}
=
R
. In this case the obtained results are compared to those presented by Grossi et al. (1997), where natural frequencies, obtained by optimized Rayleigh-Ritz method, are converted in dimensionless frequency parameters. The comparison with FEM results is also provided. Non-dimensional frequency parameters
Ω_{ij}
obtained for a rectangular plate with different length to breadth,
a
/
b
, and diameter to length ratios,
η
=
2R
/
a
, respectively, are calculated according to following expression:
and shown in
Tables 4
and
5
, for different combinations of classical boundary conditions, where differences in % between FEM and other results are given in parentheses.
Frequency parameter Ω_{ij} and its difference in % between FEM and other method of rectangular plate with circular opening, CCCC edge constraints.
Frequency parameter Ω_{ij} and its difference in % between FEM and other method of rectangular plate with circular opening, CSCS edge constraints.
Comparison of non-dimensional vibration parameters for rectangular plate with circular openings,
Tables 4
and
5
, actually shows that developed procedure as well as optimized Rayleigh-Ritz method (Grossi et al., 1997) overestimate natural frequencies. However, comparing to optimized Rayleigh-Ritz method which is more complex, developed procedure gives closer results to those obtained by finite element method, which may be considered as the most accurate ones. Also, comparison of natural modes in this case,
Fig. 9
, shows good agreement.
Mode shapes of rectangular plate with central circular opening, a / b =1, η = 0.2 , CCCC.
Natural frequencies f (Hz) of rectangular plate with oval opening.
Mode shapes of rectangular plate with oval opening, h = 0.05m , CFCC.
type opening is widely adopted in the ship and offshore structures, but at the same time its FE modelling is rather time-consuming task. In that sense, application of the presented method to this problem is advantageous.
Plate length and breadth yield 0.7
m
and 0.5
m
, respectively, while oval opening dimensions are 0.1
a_{oo}
=
m
and 0.05
b_{oo}
=
m
. First three natural frequencies of rectangular plate with oval opening, for different relative thickness and boundary conditions, respectively, are shown in
Table 6
. Corresponding mode shapes for a selected case are presented in
Fig. 10
. Similarly as in previous cases, acceptable agreement of the results is achieved.
m
and 1.0
m
, respectively, and plate thickness equal to 0.03
m
. Due to reason of simplicity, CCCC boundary conditions are applied, and ratios of opening and plate areas,
A_{o}
, and,
A_{p}
, respectively, has been varied from zero to 0.20. As expected, the best agreement of the result is achieved in case of plate without openings,
Fig. 11
. Although the difference increases with the opening increase, it still remains acceptable for larger opening size.
Natural frequencies of a plate with different sizes of rectangular opening, CCCC.
Concerning the opening size influence on natural frequency, it differs from one natural mode to another. One can notice that natural frequency of the fundamental mode slightly increases, while the second doesn’t changes too much. Third natural mode is particularly interesting when increasing
A_{o}
/
A_{p}
ratio from 0.1 to 0.15, where natural frequency doesn’t increase. This effect is actually a result of mode shape change as can be seen in
Fig. 12
.
Third natural mode of a plate with different sizes of rectangular opening, CCCC.

Natural vibration analysis
;
Flexural vibrations
;
Assumed mode method
;
Lagrange's equations
;
Rectangular plate
;
Plate opening
;
Arbitrary boundary constraints.

INTRODUCTION

Dynamic analysis of plates with openings represents an important issue in naval architecture and civil, mechanical as well as ocean engineering. Problem of determining flexural natural frequencies of plates with central free openings has been studied by many researchers (Leissa, 1969; Monahan et al., 1970; Bathe et al., 1973; Paramasivam, 1973; Hegarty and Ariman, 1975; Ali and Atwal, 1980; Laura et al., 1987; Gutierrez et al., 1987). Nowadays, the finite element method gives complete solution to a problem of vibration of plates with openings, and although it applied for very long time (Monahan et al., 1970), it still suffers from the computational time consumption. Finite difference method has been applied to a plate with opening and simply supported and clamped edges by Paramasivam (1973). Reddy (1982) analysed large amplitude flexural vibrations of layered composite plates with openings. Grossi et al. (1997) applied optimized Rayleigh-Ritz method to generate values of the fundamental frequency coefficient corresponding to the first fully antisymmetric mode for rectangular plates with circular holes and elastically restrained against rotation. In all these references positions of the openings are limited to plate central part. However, finite difference method has been applied to a plate with more than two openings by Aksu and Ali (1976). Huang and Sakiyama (1999) proposed approximate method for free vibration analysis of rectangular plates with openings of different shapes (circular, semi-circular, elliptic, square, rectangular, triangular, etc.). They consider a plate with an opening as a body with non-uniform thickness, where the opening is treated as an extremely thin part of the plate. In that sense the problem is translated into the free vibration analysis of the equivalent rectangular plate with non-uniform thickness by applying the discrete solution (Huang and Sakiyama, 1998). It should be noted that the problem is solved for arbitrary position of the opening within the plate area. An extensive literature survey on vibration analysis of rectangular plates with openings, with particular emphasis to application of Rayleigh-Ritz method, has been presented by Kwak and Han (2007). They proposed a Rayleigh-Ritz based method with different coordinate systems for a plate and opening with particular aim to simplify integration process in total energy determination.
Recently, assumed mode method using characteristic orthogonal polynomials is applied to problem of free vibration of rectangular plates of arbitrary edge constraints (Chung et al., 1993; Kim et al., 2012). In this paper, assumed mode method is applied to natural vibration analysis of plates with openings of different shape and arbitrary boundary conditions. The opening effect is accounted for by subtracting its strain and kinetic energies from corresponding total plate energies without openings. Natural frequencies and corresponding mode shapes are obtained by solving an eigenvalue problem of a multi-degree-offreedom system matrix equation, which has been derived by Lagrange's equations of motion. It should be noted also that the proposed method can be applied for arbitrary number of openings and their arbitrary positions. Numerical examples include analysis of natural vibrations of rectangular plates with rectangular, elliptic, circular and oval openings. Both, thin and thick plates are considered. The results are compared with those obtained by the finite element method as well as some results available in the relevant literature. It is shown that the presented procedure can be used as a reliable alternative to widely used FEM, which may require much more time in model preparation, as for instance in case of plate with oval opening.
MATHEMATICAL MODEL

- Equations of motion

The Mindlin thick plate theory, which takes shear influence and rotary inertia into account, is introduced (Mindlin et al., 1956). In that way the method can be applied not only to thin plates, but also to moderately thick and thick ones. The Mindlin theory deals with three general displacements, i.e. plate deflection
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

- Energy of a rectangular plate without openings

By introducing non-dimensional parameters
Lager Image

Lager Image

Lager Image

- Energy of rectangular opening

In similar way, as for a plate without openings, one can write for strain and kinetic energy of rectangular opening shown in
Fig. 2
, respectively:
Lager Image

Lager Image

Lager Image

Lager Image

- Energy of elliptic opening

Analogously to the rectangular opening, strain and kinetic energies of elliptical opening shown in
Fig. 3
yield:
Lager Image

Lager Image

- Energy of oval opening

Oval opening, which often appears in ships and offshore structures, can be treated as a combination of rectangular and circular openings,
Fig. 4
., where
Lager Image

Lager Image

Lager Image

- Outline of assumed mode method

In the assumed mode method, lateral displacement and the rotational angles are expressed by superposing products of orthogonal polynomials:
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

NUMERICAL EXAMPLES

In order to illustrate applicability of the developed procedure, natural vibration analyses of plates with different openings are carried out. Just for the illustration purposes, first, a rectangular plate without openings as a very simple case with available analytical solution is analysed. Then, rectangular plate with single central rectangular opening is analysed. Further on, the procedure is applied to rectangular plate with central elliptic opening. Also, plate with two elliptic openings is analysed as well as a rectangular plate with central circular opening. Finally, results of vibration analysis of a rectangular plate with oval opening, which is typical in naval architecture, are presented. The number of polynomials has been set to 13 in all examples, for both
- Rectangular plate without openings

A numerical calculation of the natural frequencies has been carried out for a plate without openings, simply supported on the longitudinal edges and clamped on the transversal ones SCSC (simple-clamped-simple-clamped), according to the procedure described in previous section. In
Table 1
, the frequency parameter
Lager Image

Frequency parameterλof rectangular plate without openings, SCSC edge constraints.

Lager Image

- Rectangular plate with rectangular opening

Rectangular plate with length and breadth 1.0
Natural frequencyf(Hz) of rectangular plate with central rectangular opening.

Lager Image

Lager Image

Lager Image

- Rectangular plate with elliptic openings

Natural vibration analysis of rectangular plate of same dimensions as in previous case, with central elliptic opening and two elliptic openings at arbitrary defined locations has been also performed. In the former case, the dimensions of the opening are 0.2
Natural frequenciesf(Hz) of rectangular plate with elliptic openings.

Lager Image

Lager Image

Lager Image

- Rectangular plate with circular opening

Mathematical model for vibration analysis of rectangular plates with circular openings of radius
Lager Image

Frequency parameterΩijand its difference in % between FEM and other method of rectangular plate with circular opening, CCCC edge constraints.

Lager Image

Frequency parameterΩijand its difference in % between FEM and other method of rectangular plate with circular opening, CSCS edge constraints.

Lager Image

Lager Image

- Rectangular plate with oval opening

Finally, natural vibration analysis of rectangular plate with oval opening has been performed. It should be noted that oval
Natural frequenciesf(Hz) of rectangular plate with oval opening.

Lager Image

Lager Image

- Influence of opening size on plate dynamic response

Influence of the opening size on plate natural frequencies is considered analyzing natural vibrations of rectangular plate with length and breadth 2.0
Lager Image

Lager Image

CONCLUSION

An approximate procedure for natural vibration analysis of rectangular plates with openings and arbitrary edge constraints is presented. The opening effect is accounted for in a very simple and intuitive way, i.e. by subtracting its strain and kinetic energies from the corresponding total plate energies. Numerical examples, which include analysis of natural vibrations of rectangular plates with variable thicknesses and with different opening shapes (rectangular, elliptic, circular, oval) and sizes show very good agreement with the results obtained by FEM as well as to some results that are presented in the relevant literature. Moreover, despite its simplicity, the presented procedure can give even closer results to FEM comparing to some more complex methods, as shown in the case of rectangular plate with circular opening. Anyhow, it should be mentioned that the procedure should be applied to plates with smaller openings comparing to plate dimensions, as for instance to vibration analysis of swash bulkheads with small openings etc. Thus, it can be used as a reliable alternative to widely used FEM, which may require much more time in model preparation, which is especially pronounced in case of rectangular plate with oval opening. Future investigations should be focused on the application of assumed mode method to vibration analysis of more complex structural members of ships and offshore structures, as for instance stiffened panels, stiffened panels with openings of different shape, rectangular plates with openings or stiffened panels partially or fully immersed into the water, etc.
Acknowledgements

This work was supported by the National Research Foundation of Korea (NRF), grant funded by the Korean Government(MSIP) through GCRC-SOP (Grant No. 2011-0030013).

Aksu G.
,
Ali R.
1976
Determination of dynamic characteristics of rectangular plates with cut-outs using a finitedifference formulation.
Journal of Sound and Vibration
44
147 -
158

Ali R.
,
Atwal S.J.
1980
Prediction of natural frequencies of rectangular plates with circular cutouts.
Computers and Structures
12
(6)
819 -
823

Bathe K.J.
,
Wilson E.L.
,
Peterson F.E.
1973
A structural analysis program for static and dynamic response of linearsystems.
University of California
Berkeley, California
Report No. EERC 73-11

Chung J.H.
,
Chung T.Y.
,
Kim K.C.
1993
Vibration analysis of orthotropic Mindlin plates with edges elastically restrainedagainst rotation.
Journal of Sound and Vibration
163
(1)
151 -
163

Gutierrez R.H.
,
Laura P.A.A.
,
Pombo J.L.
1987
Higher frequencies of transverse vibration of rectangular plates elasticallyrestrained against rotation at the edges with a central free hole.
Journal of Sound and Vibration
117
(1)
202 -
206

Huang M.
,
Sakiyama T.
1998
Free vibration analysis of rectangular plates with variable thickness.
Journal of Sound and Vibration
216
(3)
379 -
397

Huang M.
,
Sakiyama T.
1999
Free vibration analysis of rectangular plates with variously-shaped holes.
Journal of Sound and Vibration
226
(4)
769 -
786

Kwak M.K.
,
Han S.
2007
Free vibration analysis of rectangular plate with a hole by means of independent coordinatecoupling method.
Journal of Sound and Vibration
306
(1-2)
12 -
30

Leissa A.W.
1969
Vibration of plates, NASA Report No. SP-160.
Scientific and Technical Information Division, National Aeronautics and Space Administration

Mindlin R.D.
,
Schacknow A.
,
Deresiewicz H.
1956
Flexural vibrations of rectangular plates.
Journal of Applied Mechanics
23
430 -
436

Monahan L.J.
,
Nemergut P.J.
,
Maddux G.E.
1970
Natural frequencies and mode shapes of plates with interior cutouts.
The Shock and Vibration Bulletin
41
(3)
37 -
49

MSC
2010
MD Nastran 2010 Dynamic analysis user's guide.
MSC Software

Paramasivam P.
1973
Free vibration of square plates with square openings.
Journal of Sound and Vibration
30
(2)
173 -
178

Reddy J.N.
1982
Large amplitude flexural vibration of layered composite plates with cutouts.
Journal of Sound and Vibration
83
(1)
1 -
10

Senjanović I.
,
Vladimir N.
,
Tomić M.
2013
An advanced theory of moderately thick plate vibrations.
Journal of Sound and Vibration
332
(7)
1868 -
1880

Szilard R.
2004
Theories and applications of plate analysis.
John Wiley & Sons

Citing 'Approximate natural vibration analysis of rectangular plates with openings using assumed mode method
'

@article{ E1JSE6_2013_v5n3_478}
,title={Approximate natural vibration analysis of rectangular plates with openings using assumed mode method}
,volume={3}
, url={http://dx.doi.org/10.3744/JNAOE.2013.5.3.478}, DOI={10.3744/JNAOE.2013.5.3.478}
, number= {3}
, journal={International Journal of Naval Architecture and Ocean Engineering}
, publisher={The Society of Naval Architects of Korea}
, author={Cho, Dae Seung
and
Vladimir, Nikola
and
Choi, Tae Muk}
, year={2013}
, month={Sep}