Free Vib 1
Free Vib 1
Free Vib 1
a
y
2
b
2
y
b
}
Figure 1. Conoid (CON).
168 S. PRADYUMNA AND J.N. BANDYOPADHYAY
higher-order facet quadrilateral composite shell element. Bhimaraddi [8], Mallikarjuna
and Kant [9], Cho et al. [10] are among the others to develop higher-order shear
deformable shell theory. It is observed that except for the theory of Yang [3], remaining
higher-order theories do not account for twist curvature (1/R
xy
), which is essential while
analyzing shell forms like hypar and conoid shells.
Static and dynamic analyses of shell panels were carried out by many researchers in the
past few decades. Choi [11], with modified isoparametric element, analyzed the
conoidal shells by adding four extra nonconforming displacement modes to transverse
displacement. Ghosh and Bandyopadhyay [12,13] studied bending behavior of conoidal
shells using the finite element method and Galerkin method. Simplified bending analysis of
doubly curved shells was done by Aditya and Bandyopadhyay [14] with modified
boundary conditions. Chakravorty et al. [15,16] carried out finite element dynamic
analysis of doubly curved shells which include hyperbolic paraboloids (Figure 2), hypars
(Figure 3) and conoids. Stavridis [17], by analytical treatment of Marguerre equations with
inertia terms added, studied free vibration of elliptical, hyperbolic paraboloid, hypar,
conoid and soap bubble isotropic shells. The effects of the presence of the stiffeners and
the cutouts on these types of shells were studied by Nayak and Bandyopadhyay [18,19].
Higher-order theories were widely used for free vibration analysis of laminated composite
shells by many researchers [4,6,8,20,21].
z
A
B
a
b
x
y
R
y
R
x
Figure 2. Hyparbolic paraboloid (HPR).
z
z = 4 +cx/a+cy/b
c
xy
ab
c
y
a
b
x
Figure 3. Hypar (HYP).
Static and Free Vibration Analyses of Laminated Shells using a Higher-order Theory 169
To the best of the authors knowledge, however, literature is not available related to the
application of higher-order theory for studying the static and dynamic behavior of
laminated composite shells with the combination of all three radii of curvature. Therefore,
in the present analysis, static and free vibration behavior of laminated composite
anticlastic shells are studied. A Higher-order shear deformation theory (HSDT), developed
by Kant and Khare [7], is used by extending it to the shells with all three radii of curvature.
THEORY AND FORMULATION
Let us consider a laminated shell element made of a finite number of uniformly thick
orthotropic layers (Figures 4 and 5), oriented arbitrarily with respect to the shell
co-ordinates (x, y, z). The co-ordinate system (x, y, z) is chosen such that the plane x-y at
z 0 coincides with the mid-plane. In order to approximate the three-dimensional
elasticity problem to a two-dimensional one, the displacement components u(x, y, z),
v(x, y, z) and w(x, y, z) at any point in the shell space are expanded in Taylors series in
terms of the thickness co-ordinates. The elasticity solution indicates that the transverse
shear stresses vary parabolically through the element thickness. This requires the use of a
displacement field in which the in-plane displacements are expanded as cubic functions of
the thickness co-ordinate. The displacement fields, which satisfy the above criteria are
assumed in the form as given by Kant and Khare [7]
ux, y, z u
0
x, y z
y
z
2
u
0
x, y z
3
y
x, y
vx, y, z v
0
x, y z
x
z
2
v
0
x, y z
3
x
x, y
wx, y, z w
0
1
a
R
x
R
y
x
z
y
h
b
Figure 4. Laminated composite doubly curved shell element.
2
1
x
z
Typical lamina
y
Figure 5. Lamina reference axis and fiber orientation.
170 S. PRADYUMNA AND J.N. BANDYOPADHYAY
where u, v and w are the displacements of a general point (x, y, z) in an element of the
laminate along x, y and z directions, respectively. The parameters u
0
, v
0
, w
0
,
x
and
y
are the displacements and rotations of the middle plane, while u
0
, v
0
,
x
and
y
are the
higher-order displacement parameters defined at the mid-plane.
The linear strain-displacement relations according to Sanders approximation are:
c
x
ou
ox
w
R
x
, c
y
ov
oy
w
R
y
, ,
xy
ov
ox
ou
oy
2w
R
xy
,
,
xz
ou
oz
ow
ox
C
1
u
R
x
C
1
v
R
xy
, ,
yz
ov
oz
ow
oy
C
1
v
R
y
C
1
u
R
xy
2
Substituting Equation 1 in Equation 2:
c
x
c
x0
zk
x
z
2
c
x0
z
3
k
x
c
y
c
y0
zk
y
z
2
c
y0
z
3
k
y
,
xy
,
xy0
zk
xy
z
2
,
xy0
z
3
k
xy
,
xz
x
zk
xz
z
2
x
z
3
k
xz
,
yz
y
zk
yz
z
2
y
z
3
k
yz
3
where:
c
x0
, c
y0
, ,
xy0
, c
x0
, c
y0
, ,
xy0
ou
0
ox
w
0
R
x
,
ov
0
oy
w
0
R
y
,
ou
0
oy
ov
0
ox
2w
0
R
xy
,
ou
0
ox
,
ov
0
oy
,
ov
0
ox
ou
0
oy
k
x
, k
y
, k
xy
, k
x
, k
y
, k
xy
o
y
ox
,
o
x
oy
,
o
y
oy
o
x
ox
C
0
ou
0
oy
C
0
ov
0
ox
,
o
y
ox
,
o
x
oy
,
o
y
oy
o
x
ox
x
,
y
ow
0
ox
y
C
1
u
0
R
x
C
1
v
0
R
xy
,
ow
oy
x
C
1
v
0
R
y
C
1
u
0
R
xy
x
,
y
3
y
C
1
u
0
R
x
C
1
v
0
R
xy
, 3
x
C
1
v
0
R
y
C
1
u
0
R
xy
k
xz
, k
yz
, k
xz
, k
yz
2u
o
C
1
y
R
x
C
1
x
R
xy
, 2v
o
C
1
x
R
y
C
1
y
R
xy
, C
1
y
R
x
C
1
x
R
xy
, C
1
x
R
x
C
1
y
R
xy
4
C
1
is a tracer by which the analysis can be reduced to that of shear deformable Loves first
approximation and C
0
0.51R
x
1R
y
is the result of Sanders theory which accounts
for the condition of zero strain for rigid body motion.
Static and Free Vibration Analyses of Laminated Shells using a Higher-order Theory 171
The constitutive relations for a typical lamina k with reference to the fibre-matrix
co-ordinate axis 12 (Figure 5) are written as:
o
1
o
2
t
12
t
13
t
23
Q
11
Q
12
0 0 0
Q
12
Q
22
0 0 0
0 0 Q
66
0 0
0 0 0 Q
44
0
0 0 0 0 Q
55
k
c
1
c
2
,
12
,
13
,
23
k
5
or, in matrix form:
o f g
k
Q
k
c f g
k
where, Q
11
E
1
1v
12
v
21
, Q
12
v
12
E
2
1v
12
v
21
, Q
22
E
2
1 v
12
v
21
, Q
66
G
12
, Q
44
G
13
, Q
55
G
23
and v
12
E
1
v
21
E
2
For the elastic constant matrix to any arbitrary axis with which the material principal
axes make an angle , standard co-ordinate transformation is required. Thus, the off-axis
elastic constant matrix is obtained from the on-axis elastic constant matrix using the
relation Q
ij
T
T
Q
ij
T, where [T] is the transformation matrix. Therefore, the stress-
strain relations for a lamina about any axis are given by:
o
x
o
y
t
xy
t
xz
t
yz
Q
11
Q
12
Q
16
0 0
Q
12
Q
22
Q
26
0 0
Q
16
Q
26
Q
66
0 0
0 0 0 Q
44
Q
45
0 0 0 Q
45
Q
55
k
c
x
c
y
,
xy
,
xz
,
yz
k
. 6
Integrating the stresses through the laminate thickness, the resultant forces and
moments acting on the laminate are obtained.
N, N
N
x
N
x
N
y
N
y
N
xy
N
xy
N
L
k1
z
k1
z
k
o
x
o
y
t
xy
1, z
2
dz
M, M
M
x
M
x
M
y
M
y
M
xy
M
xy
N
L
k1
z
k1
z
k
o
x
o
y
t
xy
z, z
3
dz
Q, S, Q
, S
Q
x
S
x
Q
x
S
x
Q
y
S
y
Q
y
S
N
L
k1
z
k1
z
k
t
xz
t
yz
1, z, z
2
, z
3
7
172 S. PRADYUMNA AND J.N. BANDYOPADHYAY
or o Dc, where:
o N
x
, N
y
, N
xy
, N
x
, N
y
, N
xy
, M
x
, M
y
, M
xy
, M
x
, M
y
, M
xy
, Q
x
, Q
y
, Q
x
, Q
y
, S
x
, S
y
, S
x
, S
y
T
c c
x0
, c
y0
, ,
xy0
, c
x0
, c
y0
, ,
xy0
, k
x
, k
y
, k
xy
, k
x
, k
y
, k
xy
,
x
,
y
,
x
,
y
, k
xz
, k
yz
, k
xz
, k
yz
T
D
D
m
D
c
0
D
T
c
D
b
0
0 0 D
s
where D
m
, D
b
, D
c
and D
s
are given in Appendix A.
FINITE ELEMENT FORMULATION
An eight-noded isoparametric C8 element with nine degrees of freedom per node is used.
The displacement vector d at any point on the mid-surface is given by:
d
8
i1
N
i
x, yd
i
where d
i
is the displacement vector corresponding to node i and N
i
is the interpolating
function associated with the node i.
Knowing the generalized displacement vector {d} at all points within the element, the
generalized mid-surface strains at any point given by Equation 3, are expressed in terms of
global displacements in the matrix form as:
c
8
i1
B
i
d
i
8
where B
i
is a differential operator matrix of interpolation functions and obtained from
Equation 4. The element stiffness matrix for the element e, which includes membrane,
flexure and the transverse shear effects, and the element mass matrix are given by the
following equations:
K
e
1
1
1
1
B
T
DB J j jdrds
M
e
1
1
1
1
N
T
mN J j jdrds
where [N] is the shape function matrix and [m] is the inertia matrix, as given in Appendix
B. In all the numerical computations, the selective integration rule is employed. A 3 3
Gaussian rule is used to compute in-plane, coupling between in-plane and bending
deformations, while a 2 2 rule is used to evaluate the terms associated with transverse
shear deformation. The element mass matrix is evaluated using a 3 3 Gaussian rule.
Static and Free Vibration Analyses of Laminated Shells using a Higher-order Theory 173
The element matrices are then assembled to obtain the global [K] and [M] matrices.
The free vibration analysis involves determination of natural frequencies from the
condition:
K o
2
n
M 0
This is a generalized eigenvalue problem and is solved by using the subspace iteration
algorithm.
RESULTS AND DISCUSSIONS
A computer program is developed based on the above formulation. A parallel program
is developed based on the first order shear deformation theory (FSDT) in order to
compare the results with those of HSDT. The shell-forms mainly considered here are
hyperbolic paraboloid, hypar and conoid shells.
The following two boundary conditions are used in the present analysis:
(1) Simply supported boundary (S-S-S-S) having SS1 for isotropic and cross-ply
laminates and SS2 for angle-ply laminates:
SS1: v
0
w
0
y
v
y
0, at x 0;
u
0
w
0
x
u
x
0, at y 0, b.
SS2: u
0
w
0
y
u
y
0, at x 0, a; and
v
0
w
0
x
v
x
0, at y 0, b.
(2) Clamped boundary (C-C-C-C): u
0
v
0
w
0
x
y
u
0
v
y
0, at
x0, a and y 0, b.
Non-dimensional length parameters are designated by x
x
a
and y
y
b
.
Non-dimensional center deflection parameter o
wh
3
E
2
qa
4
1000.
Non-dimensional frequency parameter o
oa
2
h
aE
2
N
L
k1
Q
11
H
1
Q
12
H
1
Q
16
H
1
Q
11
H
3
Q
12
H
3
Q
16
H
3
Q
22
H
1
Q
26
H
1
Q
12
H
3
Q
22
H
3
Q
26
H
3
Q
66
H
1
Q
16
H
3
Q
26
H
3
Q
66
H
3
Q
11
H
5
Q
12
H
5
Q
16
H
5
Symm. Q
22
H
5
Q
26
H
5
Q
66
H
5
N
L
k1
Q
44
H
1
Q
45
H
1
Q
44
H
3
Q
45
H
3
Q
44
H
2
Q
45
H
2
Q
44
H
4
Q
45
H
4
Q
55
H
1
Q
45
H
3
Q
55
H
3
Q
45
H
2
Q
55
H
2
Q
45
H
4
Q
55
H
4
Q
44
H
5
Q
45
H
5
Q
44
H
4
Q
45
H
4
Q
44
H
6
Q
45
H
6
Q
55
H
5
Q
45
H
4
Q
55
H
4
Q
45
H
6
Q
55
H
6
Q
44
H
3
Q
45
H
3
Q
44
H
5
Q
45
H
5
Symm. Q
55
H
3
Q
45
H
5
Q
55
H
5
Q
44
H
7
Q
45
H
7
Q
55
H
7
The parameters I
1
, I
3
and (I
5
, I
7
) are linear inertia, rotary inertia and higher-order inertia
terms, respectively. The parameters I
2
, I
4
and I
6
are the coupling inertia terms and are
expressed as follows:
I
1
, I
2
, I
3
, I
4
, I
5
, I
6
, I
7
N
L
k1
z
k1
z
k
1, z, z
2
, z
3
, z
4
, z
5
, z
6
a
k
dz,
where a
k
is the material density of the kth layer and the shape function matrix [N] is
given by:
N
N
i
N
i
N
i
N
i
N
i
N
i
N
i
N
i
N
i
i 1, 8 and N
i
is the shape function for the node i.
REFERENCES
1. Reissner, E. (1945). The Effect of Transverse Shear Deformation on the Bending of Elastic Plates, ASME
Journal of Applied Mechanics, 12(2): 6977.
2. Mindlin, R. D. (1951). Influence of Rotary Inertia and Shear in Flexural Motions of Isotropic Elastic Plates,
ASME Journal of Applied Mechanics, 18(2): 519540.
Static and Free Vibration Analyses of Laminated Shells using a Higher-order Theory 185
3. Yang, T. Y. (1973). High Order Rectangular Shallow Shell Finite Element, Journal of the Engineering
Mechanics Division, 99(EM1): 157181.
4. Reddy, J. N. and Liu, C. F. (1985). A Higher-order Shear Deformation Theory of Laminated Elastic Shells,
Int. J. Engng Sci., 23(3): 319330.
5. Shu, X. and Sun, L. (1994). An Improved Simple Higher-order Theory for Laminated Composite Plates,
Computers & Structures, 50(2): 231236.
6. Liew, K. M. and Lim, C. W. (1996). A Higher-order Theory for Vibration of Doubly Curved Shallow Shells,
ASME, Journal of Applied Mechanics, 63: 587593.
7. Kant, T. and Khare, R. K. (1997). A Higher-order Facet Quadrilateral Composite Shell Element, Int. J. for
Numerical Methods in Engng, 40: 44774499.
8. Bhimaraddi, A. (1984). A Higher Order Theory for Free Vibration Analysis of Circular Cylindrical Shells,
International Journal of Solids Structures, 20(7): 623630.
9. Mallikarjuna and Kant, T. (1992). A General Fibre-reinforced Composite Shell Element Based on a Refined
Shear Deformation Theory, Computers & Structures, 42(3): 381388.
10. Cho, M., Kim, K. and Kim, M. (1996). Efficient Higher-order Shell Theory for Laminated Composites,
Composite Structures, 34(2): 197212.
11. Choi, C. K. (1984). A Conoidal Shell Analysis by Modified Isoparametric Element, Computers & Structures,
18(5): 921924.
12. Ghosh, B. and Bandyopadhyay, J. N. (1989). Bending Analysis of Conoidal Shells using Curved Quadratic
Isoparametric Element, Computers & Structures, 33(4): 717728.
13. Ghosh, B. and Bandyopadhyay, J. N. (1990). Approximate Bending Analysis of Conoidal Shells Using the
Galerkin Method, Computers & Structures, 36(5): 801805.
14. Aditya, A. K. and Bandyopadhyay, J. N. (1989). Simplified Bending Analysis of Doubly Curved Shells,
Computers & Structures, 33(3): 781784.
15. Chakravorty, D., Bandyopadhyay, J. N. and Sinha, P. K. (1996). Finite Element Free Vibration Analysis of
Doubly Curved Laminated Composite Shells, Journal of Sound and Vibration, 191(4): 491504.
16. Chakravorty, D., Bandyopadhyay, J. N. and Sinha, P. K. (1998). Application of FEM on Free and Forced
Vibration of Laminated Shells, ASCE Journal of Engineering Mechanics, 124(1): 18.
17. Stavridis, L. T. (1998). Dynamic Analysis of Shallow Shells of Rectangular Base, Journal of Sound and
Vibration, 218(5): 861882.
18. Nayak, A. N. and Bandyopadhyay, J. N. (2002). Free Vibration Analysis and Design Aids of Stiffened
Conoid Shells, ASCE J. Eng. Mech., 124(4): 419427.
19. Nayak, A. N. and Bandyopadhyay, J. N. (2003). Free Vibration Analysis of Laminated Stiffened Shells,
ASCE J. Eng. Mech., 131(1): 100105.
20. Moita, J. S., Soares, C. M. M. and Soares, C. A. M. (1999). Buckling and Dynamic Behavior of
Laminated Composite Structures Using a Discrete Higher-order Displacement Model, Computers
& Structures, 73(15): 407423.
21. Khare, R. K., Kant, T. and Garg, A. K. (2004). Free Vibration of Composite and Sandwich Laminates with
a Higher-order Facet Shell Element, Composite Structures, 65: 405418.
22. Bhimaraddi, A. (1993). Three-dimensional Elasticity Solution for Static Response of Orthotropic Doubly
Curved Shallow Shells on Rectangular Planform, Composite Structures, 24(1): 6777.
23. Pagano, N. J. (1970). Exact Solutions for Rectangular Bidirectional Composite and Sandwich Plates,
Journal of Composite Materials, 4: 2034.
24. Reddy, J. N. (1984). A Simple Higher-order Theory for Laminated Composite Plates, ASME Journal of
Applied Mechanics, 51: 745752.
186 S. PRADYUMNA AND J.N. BANDYOPADHYAY