CS 2006
CS 2006
CS 2006
www.elsevier.com/locate/compstruct
Department of Mechanical Engineering and Applied Mechanics, North Dakota State University, Fargo, ND 58105, USA
Available online 19 December 2004
Abstract
In this paper, natural frequencies of square laminated composite plates for dierent supports at edges are presented. Using a third
order shear deformation theory of plates (TSDT), which is categorized in equivalent single layer theories (ESL), a new set of linear
equations of motion for square multi-layered composite plates has been derived. Laminated plates are supposed to be either angleply or cross-ply. Moreover, FEM is used to solve the equations and nd the fundamental natural frequencies. Finally some results
for plates with dierent combination of layers and supports are reported. The results are compared to the results of other ESL.
2004 Elsevier Ltd. All rights reserved.
Keywords: Natural frequency; Third order shear deformation (TSDT); Laminated composite plate; Finite element method (FEM)
1. Introduction
Laminated composite plates are widely used in industry and new elds of technology. Due to high degrees of
anisotropy and low rigidity in transverse shear, Kirchho hypothesis as a classical theory is no longer adequate. The hypothesis states that transverse normal to
the mid-plane of a plate remains straight and normal
after deformation because of the negligible transverse
shear eects. Rened theories based on removing those
restrictions of transverse normal have been recently used.
As a result, the free vibration frequencies calculated by
using the classical thin plate theory are higher than those
obtained by Mindlin plate theory in which transverse
shear and rotary inertia eects are included [1].
Most of the structural theories used up to now, characterize the behavior of composite laminates fall into the
category of equivalent single layer (ESL) theories.
Dierent techniques for analysis of free vibration of
*
274
Nomenclature
C1
o
ei
e0i
Ei
{F}
h
I1. . .7
[K]
k ji
[M]
Mi
Ni
Ri
q
Qij
Qij
Qi
Pi
u
u0
u
v
v0
v
w
w0
w
x
xn
{X}
fX g
y
yn
z
zn
d
/x, /y
,/
/
x y
/z
/
z
uz
z
u
ri
mij
h
xm
x
X
Xn
y
z
275
6a
r5 Q45 e4 Q55 e5 0
6b
r4 Q44 e4 Q45 e5 0
7a
3C 1 h2
ow0
ow0
r5 1
/y Q55
/x
Q45
4
oy
ox
The following displacement eld is a third order displacement eld and has seven independent variables u, v,
w, /x, /y, /z and uz.
ow0
1 ouz
2 1 o/z
3
u u0 z/x z
/x
z C1
2 ox
3 ox
ox
1
2
v v0 z/y z
1 o/z
2 oy
ow0
1 ouz
/y
z C1
3 oy
oy
3
2
w w0 z/z z2 uz
where
u0 ux; y; 0; t; v0 vx; y; 0; t and w0 wx; y; 0; t
4
are the displacements of transverse normal on plane
z = 0. /x and /y are rotations of transverse normal on
plane z = 0, /z is extension of transverse normal, and
uz is interpreted as a higher order rotation of transverse
normal. According to displacement eld Eqs. (1)(3) the
linear strains are
e1
e2
e3
e4
e5
e6
e01
e02
e03
e04
e05
e06
zk 01 zk 11 z2 k 21
zk 02 zk 12 z2 k 22
zk 03
zk 14 2e23
zk 15 2e13
zk 06 k 16 k 26
In Appendix A, the relationships between strain components (5) and displacement eld Eqs. (1)(3) are
presented.
7b
w0
1 I 3 o/z C 1 I 4 o
N 1;x N 6;y I 1 u0 I 2 /
x
2 ox
ox
1
o
u
z
I4
C1I 4/
x
3 ox
dv0 :
w0
1 I 3 o/z C 1 I 4 o
N 2;y N 6;x I 1v0 I 2 /
y
2 oy
oy
uz
1 I 4 o
C1I 4/
y
3 oy
10
276
dw0 :
C 1 P 1;xx C 1 P 2;yy Q4;y Q5;y 3C 1 R4;y
3C 1 R5;x 2C 1 P 6;xy q
o
u0
ov0
1u
C1I 4
0 I 2/
C1I 4
I 1w
z
3 z
ox
oy
2
o/y 1
o/
o/
1
o2 /
C1I 5 x C1I 5
C 1 I 6 2z C 1 I 6 2z
2
2
ox
oy
ox
oy
2
2
o/y
0
0
ow
ow
o/
C 21 I 7 2 C 21 I 7 2 C 21 I 7 x C 21 I 7
ox
oy
ox
oy
2
2
1 ou
1 ou
I 7 2z I 7 2z
11
3 ox
3 oy
vi hvi ifN g
MfX g KfX g fF g
16
I 2 C 1 I 4 v0 I 3 2C 1 I 5 C 2 I 7 /
y
1
o/
o
w0
I 4 C 1 I 6 z C 1 I 5 C 21 I 7
2
oy
oy
1
o
u
I 5 C 1 I 7 z
3
oy
17
12
d/y :
13
d/z :
1
1
h
R1;xx R2;yy R6;xy N 3 q
2
2
2
o
u0
ov0
I 3u
z C1I 4
0 I 2/
I 1w
C1I 4
z
ox
oy
o/
o/
y
C 1 I 5 C 21 I 7 x C 1 I 5 C 21 I 7
ox
oy
0
0 1
o2 w
o2 w
o2 /
1
o2 /
C 21 I 7 2 C 21 I 7 2 C 1 I 6 2z C 1 I 6 2z
2
2
ox
oy
ox
oy
1 o2 u
1 o2 u
I 7 2z I 7 2z
3 ox
3 oy
1
o/
o
w0
I 2 C 1 I 4 z C 1 I 5 C 21 I 7
2
ox
ox
1
o
uz
I 5 C 1 I 7
3
ox
ui hui ifN g
d/x :
M 1;x C 1 P 1;x Q5 Q5;y 3C 1 R5 M 6;y C 1 P 6;y
I 2 C 1 I 4
u0 I 3 2C 1 I 5 C 2 I 7 /
1
14
d/z :
1
1
2
h2
P 1;xx R2;yy R6;xy 2M 3 q
3
3
3
4
1
o
u
1 ov0
ov0
0
I 5u
z I4
0 I 4/
I4
C1I 4
I 3w
z
3 ox 3 oy
oy
o
/
0
1
o/
1
1
o2 w
y
C1I 7 2
I 5 C 1 I 7 x I 5 C 1 I 7
3
3
ox 3
oy
oy
1 o2 /
1 o2 /
1 o2 u
1 o2 u
I 6 2z I 6 2z I 7 2z I 7 2z
15
6 ox
6 oy
9 ox
9 oy
18
19
v0 0; y v0 a; y 0
w0 0; y w0 a; y 0
277
Table 2
Nondimensional natural frequencies of a square angle-ply (45/45/45/
45/45) laminated composite plate with dierent support conditions
/y 0; y /y a; y 0
/z 0; y /z a; y 0
uz 0; y uz a; y 0
20
N 1 0; y N 1 a; y 0
M 1 0; y M 1 a; y 0
P 1 0; y P 1 a; y 0
a/h
10
20
50
100
SSSS
SSSC
SCSC
SFSF
SFSS
SFSC
11.196
11.752
11.595
4.413
6.953
6.988
19.059
20.501
19.333
6.123
9.593
9.788
21.371
27.466
26.303
7.161
11.790
11.847
24.302
29.465
29.581
7.873
12.480
13.477
25.949
30.830
30.018
8.163
13.001
13.995
S 1 0; y S 1 a; y 0
For the other two edges (y = 0 and y = b), boundary
conditions are as follow:
For S (simply supported)
Table 3
Nondimensional natural frequencies of a square angle-ply (45/45/45/
45) laminated composite plate with dierent support conditions
a/h
10
20
50
100
u0 w0 /x /z uz 0
SSSS
SSSC
SCSC
SFSF
SFSS
SFSC
11.151
11.705
11.549
4.395
6.925
6.960
18.964
20.399
19.237
6.093
9.545
9.739
21.296
27.370
26.211
7.136
11.749
11.806
24.133
29.260
29.375
7.818
12.393
13.383
25.897
30.768
29.958
8.147
12.975
13.967
21
N 2 M 2 P 2 S2 0
For C (clamped)
u0 v0 w0 /x /y /z uz 0
N 2 N 6 N 4 C1N 4 0
22
and
M 2 M 6 S2 S6 P 2 P 6 0
Table 4
Nondimensional natural frequencies of a square angle-ply (45/45/45)
laminated composite plate with dierent support conditions
4. Numerical solutions
In Table 1, natural frequencies of dierent approaches that had been reported in other papers have
been compared to the results of the present work.
Results of present work (TSDT) for dierent ratio of
longitudinal modulus to transverse modulus are presented. TSDT results show a small dierence to the analytical Levy solution for other third order displacement
elds. A signicant improvement in the result in comparison to the FSDT and CLPT result is shown. Also
it is seen that by increasing E1/E2, the dierence between
results of methods are decreasing.
Tables 29 present nondimensional natural frequencies of square laminated plates with dierent layers
and boundary conditions. Plates are squared with sides
(a = b).
Table 1
Nondimensional natural frequencies of a SSSS square laminated
composite cross-ply (0/90/90/0) with dierent E1/E2 ratios and a/h = 5
E1/E2
TSDTa
HSDT [13]b
FSDT [13]
CLPT [13]
10
20
30
40
8.2741
9.5312
10.2651
10.7912
8.2940
9.5439
10.2840
10.7940
8.2982
9.5671
10.3260
10.8540
10.650
13.948
16.605
18.891
a/h
10
20
50
100
SSSS
SSSC
SCSC
SFSF
SFSS
SFSC
10.560
11.085
10.937
4.162
6.558
6.591
17.959
19.318
18.217
5.770
9.039
9.223
20.167
25.919
24.822
6.758
11.126
11.180
22.854
27.709
27.818
7.404
11.736
12.674
24.524
29.137
28.370
7.715
12.287
13.227
Table 5
Nondimensional natural frequencies of a square angle-ply (45/45)
laminated composite plate with dierent support conditions
a/h
10
20
50
100
SSSS
SSSC
SCSC
SFSF
SFSS
SFSC
9.701
10.183
10.048
3.824
6.025
6.055
16.499
17.747
16.736
5.301
8.304
8.473
18.528
23.812
22.804
6.208
10.222
10.271
20.996
25.456
25.556
6.802
10.782
11.643
22.530
26.768
26.063
7.088
11.288
12.151
Table 6
Nondimensional natural frequencies of a square cross-ply (0/90/0/90/0/
90) laminated composite plate with dierent support conditions
a/h
10
20
50
100
SSSS
SSSC
SCSC
SFSF
SFSS
SFSC
9.016
9.256
10.287
3.450
3.818
5.189
12.300
14.974
17.522
3.782
4.283
6.243
14.130
17.343
23.448
3.828
4.406
6.516
15.040
18.781
30.392
3.909
4.461
6.809
15.177
18.960
34.179
3.936
4.617
6.830
278
Table 7
Nondimensional natural frequencies of a square cross-ply (0/90/0/90)
laminated composite plate with dierent support conditions
a/h
10
20
50
100
SSSS
SSSC
SCSC
SFSF
SFSS
SFSC
8.953
9.191
10.215
3.426
3.792
5.153
12.214
14.870
17.400
3.755
4.253
6.199
14.032
17.222
23.285
3.802
4.376
6.471
14.936
18.650
30.180
3.882
4.430
6.761
15.071
18.829
33.942
3.909
4.585
6.783
Table 8
Nondimensional natural frequencies of a square cross-ply (0/90/0)
laminated composite plate with dierent support conditions
a/h
10
20
50
100
SSSS
SSSC
SCSC
SFSF
SFSS
SFSC
8.935
9.173
10.195
3.419
3.784
5.143
12.19
14.84
17.365
3.748
4.245
6.187
14.004
17.188
23.239
3.794
4.367
6.458
14.906
18.613
30.12
3.874
4.421
6.748
15.041
18.791
33.874
3.901
4.576
6.769
Table 9
Nondimensional natural frequencies of a square cross-ply (0/90)
laminated composite plate with dierent support conditions
a/h
10
20
50
100
SSSS
SSSC
SCSC
SFSF
SFSS
SFSC
8.220
8.439
9.379
3.145
3.481
4.732
11.215
13.653
15.976
3.448
3.905
5.692
12.884
15.813
21.380
3.490
4.018
5.941
13.714
17.124
27.710
3.564
4.067
6.208
13.838
17.288
31.164
3.589
4.210
6.227
(45/45/45/45/45) angle-ply. Table 3 is the natural frequencies of a (45/45/45/45) angle-ply, Table 4 shows
the natural frequencies of a (45/45/45) and Table 5
Shows the results for a (45/45) angle-plies. The material properties for these plates are those of material 1.
It is seen that growth of fundamental frequencies with
respect to a/h are decreased when a/h is increased. Furthermore, with an increase in the number of layers, the
rate of increase in natural frequencies will be reduced.
The dierence between natural frequencies of a two
layer angle-ply and a three layer angle-ply is much more
than the dierence between a four layer and a ve layer
plate respectively. By considering the natural frequencies at a/h = 10, the relative dierence between a three
layer and a two layer angle-ply is 8.84%, between a four
layer and a three layer is 5.59%, and between a ve layer
and a four layer is 0.5%. Tables 58 are natural frequencies of dierent cross-ply laminates. Table 6 is the results
of a six layered (0/90/0/90/0/90) cross-ply. Tables 79 are
natural frequencies of a (0/90/0/90), (0/90/0) and (0/90)
cross-ply respectively. The material of the cross-ply laminates is material 2. The same pattern and a close range
of relative dierence can be seen in natural frequencies
of plates when other support conditions presents. Same
to an angle-ply plate; the rate of change in natural frequencies reduces by increasing a/h. Taking nondimensional fundamental frequencies at a/h = 10 shows a
relative dierence of 8.6% between a three layer and a
two layer cross-plies. The relative dierence is 1.9% between a four layer and a three layer and 0.7% between
a six layer and a four layer respectively.
5. Conclusions
When comparing the same number of layers, angleplies have higher natural frequencies than cross-plies.
Increasing the number of layers causes a larger dierence between their natural frequencies. It is seen that a
cross-ply plate with SCSC boundary conditions, has a
higher natural frequency than other types of boundary
conditions, but for an angle-ply, the SSSC boundary
conditions show a higher natural frequency than other
boundary conditions. Also for a cross-ply, the natural
frequencies of SFSF and SFSS boundary conditions
are closer to each other and smaller than SFSC type.
On the contrary, natural frequencies for an angle-ply,
SFSS and SFSC boundary conditions are closer to each
other and higher than the natural frequency of a plate
with SFSF boundary conditions. Comparison between
the natural frequencies of SSSS plates (a/h = 5) as a
function of the orthotropicity ratio with the results of
analytical levy solution for other displacement elds reported in [12,13] shows compatibility of the results of
this new set of equations to the other shear deformation
models.
e01
k 01
k 21
o/x
ox"
C1
e02
k 12
ov0
oy
1 o2 /z
2 oy 2
k 02
k 22
C1
k 03 2uz
k 13 0
ow0
/
e04
oy y
ow0
1
k 4 3C 1
/y
oy
ow0
/x
e05
ox
ow0
/x
k 15 3C 1
ox
k 23 0
ou0 ov0
oy
ox
k 16
o2 u0
oxoy
o2 w0 o/y
oy 2
oy
#
1 o2 wz
3 oy 2
e03 /z
e06
References
1 o2 wz
3 ox2
o/y
oy
"
o2 w0 o/x
ox2
ox
#
k 04 0
k 24 0
k 05 0
k 25 0
o/x o/y
oy
ox
"
o2 w0 o/x
k 26 C 1 2
oxoy
oy
#
o/y
2 o2 uz
3 oxoy
ox
k 06
279