effect of twist and rotation on vibration of functionally...
TRANSCRIPT
Effect of twist and rotation on vibration of functionallygraded conical shells
Sudip Dey • Subhrajyoti Sarkar • Apurba Das •
Amit Karmakar • Sondipon Adhikari
Received: 15 February 2014 / Accepted: 10 June 2014 / Published online: 28 August 2014
� Springer Science+Business Media Dordrecht 2014
Abstract This paper presents the effect of twist and
rotational speeds on free vibration characteristics of
functionally graded conical shells employing finite
element method. The objective is to study the natural
frequencies, the influence of constituent volume frac-
tions and the effects of configuration of constituent
materials on the frequencies. The equation of dynamic
equilibrium is derived from Lagrange’s equation
neglecting the Coriolis effect for moderate rotational
speeds. The properties of the conical shell materials are
presumed to vary continuously through their thickness
with power-law distribution of the volume fractions of
their constituents. The QR iteration algorithm is used
for solution of standard eigenvalue problem. Computer
codes developed are employed to obtain the numerical
results concerning the combined effects of twist angle
and rotational speed on the natural frequencies of
functionally graded conical shells. The mode shapes
for typical shells are also depicted. Numerical results
obtained are the first known non-dimensional frequen-
cies for the type of analyses carried out here.
Keywords Functionally graded material �Vibration � Rotational speed � Finite element � Twistangle � Conical shell
Abbreviations
A0, A-1, A1, A2, A3 Temperature coefficients
T Temperature in Kelvin
Ai Material property
Vf Volume fraction
N Material property graded index
(or Power law exponent)
t Thickness
E Young’s modulus
G Shear modulus
t Poisson’s ratio
q Mass density
r1 , r2 Radius of curvature (r1[ r2)
rx Radius of curvature in x-direction
ry Radius of curvature in y-direction
rxy Radius of twist
Lo Length
bo Reference width
S. Dey (&) � S. AdhikariCollege of Engineering, Swansea University,
Swansea SA2 8PP, UK
e-mail: [email protected]; [email protected]
S. Adhikari
e-mail: [email protected];
S. Sarkar � A. Das � A. Karmakar
Mechanical Engineering Department, Jadavpur
University, Kolkata 700032, India
e-mail: [email protected]
A. Das
e-mail: [email protected]
A. Karmakar
e-mail: [email protected]
123
Int J Mech Mater Des (2015) 11:425–437
DOI 10.1007/s10999-014-9266-x
/ve Vertex angle
/o Base subtended angle of cone
W Twist angle
Lf Lagrangian function
X Non-dimensional speed of
rotation (X0/ xo)
X0 Actual angular speed of rotation
xo Fundamental natural frequency of
a non-rotating shell
x Non-dimensional frequency
param-
eter
k Non-dimensional frequency
[Me], [M] Element mass matrix and Global
mass matrix
[Ce], [C] Element Coriolis matrix and
Global Coriolis matrix
[Ke], [K] Element elastic stiffness matrix
and Global elastic stiffness matrix
[Kre], [Kr] Element geometric stiffnessmatrix,
Global geometric stiffness matrix
[KRe], [KR] Element rotational stiffness matrix
and Global rotational stiffness
matrix
{F(X2)} Global vector of nodal equivalent
centrifugal forces
{Fce}, {Fe}, {F} Element load vector due to
centrifugal force, Element load
vector due to externally applied
load and Global vector of
externally applied load
{de},{d} Element displacement vector and
Global displacement vector
1, v Local natural coordinates of the
element
x, y, z Local coordinate axes (plate
coordinate system)
Lo/s Aspect ratio
xn Natural frequency of rotating
shell
NDFF Non-dimensional fundamental
natural frequency
NDSF Non-dimensional second natural
frequency
1 Introduction
Functionally graded materials (FGM) have gained
immense popularity as advanced engineering materi-
als due to its high corrosion and heat resistance
properties along with high strength and stiffness. Such
materials are the novel inhomogeneous materials
consisting of a mixture of isotropic materials in which
the material properties are graded but continuous
particularly along a given direction to improve the
quality of performance. In general, these materials are
made from a mixture of ceramic and metal to take the
advantage of coupling the desirable properties such as
heat and corrosion resistance of ceramics and high
toughness, tensile strength and bonding capability of
metals. Functionally graded shallow pretwisted con-
ical shells could be employed as rotating turboma-
chinery blades in aerospace, nuclear and mechanical
industries. Hence, for the purpose of design and
manufacturing, the accurate evaluation of free vibra-
tion characteristics of such pretwisted shells under
rotation is essential to arrest the zone of resonance.
Therefore, a profound understanding of dynamic
behaviour of such functionally graded (FG) conical
shell is important for the designers to ensure the
operational safety. The investigations on FGMs for the
vibration analysis were carried out by Pradhan et al.
(2000), Loy et al. (1999), Ng et al. (2001). The free
vibration analysis of rotating truncated circular ortho-
tropic conical shells has been studied by Hua (2000).
However, to the best of authors’ knowledge, the only
available works regarding the free vibration analysis
of rotating FG shells of revolution are mostly limited
to FG cylindrical shells (Ahmad and Naeem 2009;
Malekzadeh and Heydarpour 2012). On the other
hand, the free vibration behaviours of stationary FG
truncated conical shells (Bhangale et al. 2006; Sofiyev
2009; Tornabene et al. 2009; Setoodeh et al. 2012)
have been investigated in the recent years. It is to be
noted that these investigations are concerned with the
stress and deformation analysis of the rotating FG
cylindrical shells (Ghafoori and Asghari 2012; Hey-
darpour et al. 2012). Other studies on thermo-elastic
buckling analysis of plates and shells are extensive.
Many of these studies are for isotropic plates and
shells which include Birman and Bert (1993), Tauc-
hert (1991), Thornton (1993), Gotsis and Guptill
(1994), Eslami et al. (1996) and Sofiyev (2007). Finite
element method employed for the analysis of shells
426 S. Dey et al.
123
subjected to vibration was also examined by He et al.
(2002), Ng et al. (2002), Liew et al. (2004), Shu (1996)
and Bardell et al. (1998). Much effort has been
devoted to various structural analyses of functionally
graded structures, such as thermal stresses, static,
buckling, and vibration analyses. Remarkable inves-
tigation is also carried out by Noda (1999), Fukui et al.
(1993) and Praveen and Reddy (1998).
In the present study, two types of functionally
graded materials, namely Stainless Steel/Ni FG and
Ti–6Al–4V/aluminium oxide FGM are investigated to
address the effects of the volume fraction and bound-
ary condition on their natural frequencies. To bridge
the research gap identified in the review of open
literature, the effect of rotational speeds on free
vibration analysis of functionally graded conical shells
is studied. To the best of authors’ knowledge and
considering the review of open literature, it is revealed
that the combined effect of rotation and twist on
functionally graded shallow conical shell structure has
not received due attention. To fill up this apparent
void, present study employed a finite element based
numerical approach to study the free vibration char-
acteristics of pretwisted functionally graded rotating
shallow conical shells which could be idealized as
turbomachinery blades. The shell elements are devel-
oped by superimposing the stiffness of the membrane
element and plate bending element. The finite element
analyses are carried out using an eight-noded isopara-
metric quadratic element and QR iteration algorithm
(Bathe 1990) is utilized to solve the standard eigen-
value problem. The non-dimensional fundamental
natural frequencies are obtained considering the
effects of triggering parameters like the twist angle
and rotational speed.
2 Theoretical formulation
Functionally graded materials (FGM) are defined as the
combination of two or more materials. Most of the
functionally graded materials are employed in high-
temperature environments and many of the constituent
materials may possess temperature-dependent properties.
Thematerial properties (A) can be expressed as a function
of temperature (Touloukian 1967) as,
A ¼ A0 þ A�1 T�1 þ 1 þ A1T þ A2T
2 þ A3T3
ð1Þ
where A0, A-1, A1, A2 and A3 are the coefficients of
temperature T(K) expressed in Kelvin and are unique
to the constituent materials and ‘A’ is a function of the
material properties and volume fractions of the
constituent materials, and is expressed as
A ¼Xk
i¼1
Ai Vf i ð2Þ
where Ai and Vf are material property and volume
fraction of the constituent materials i, respectively.
The volume fractions of all the constituent materials
can be related as
Xk
i¼1
Vf i ¼ 1 and Vfi ¼z þ t
2
t
� �Nð3Þ
where N is the material property graded index which is
a positive real number (0 B N B a). The material
properties of the FG conical shells is varied continu-
ously and smoothly in the thickness direction, i.e. z-
direction, of FG conical shell. Without loss of gener-
ality of the formulation and the method of solution in
this study, the effective material properties are obtained
using the power law distribution (Loy et al. 1999),
EðzÞ ¼ Em þ ðEc � EmÞ ½ð2z þ tÞ=2t�N ð4Þ
Fig. 1 a Turbo-machinery
blade idealized as pretwisted
conical shells, b Conical
shell model
Effect of twist and rotation on vibration 427
123
tðzÞ ¼ tm þ ðtc � tmÞ ½ð2z þ tÞ=2t�N ð5Þ
qðzÞ ¼ qm þ ðqc � qmÞ ½ð2z þ tÞ=2t�N ð6Þ
where E, t and q denote Young’s modulus, Poisson’s
ratio and mass density with suffix as ‘c’ and ‘m’
indicating the corresponding values at the outer
surface (ceramic rich) and inner surface (metal rich)
of the FG conical shell. A functionally graded conical
shell panel is shown in Fig. 1, wherein a coordinate
system (x, y, z) is established on the middle surface of
the conical shell panel. The geometric parameters of
FG conical shell are represented by length Lo, semi-
vertex angle/ve , base subtended angle/o , thickness t,
and the radii at the two ends r1 and r2 (r1[ r2). In the
spanwise direction of FG conical shell, no curvature is
considered (i.e., rx = ?). The cone radius at any point
along its length is expressed as (Zhao and Liew 2011)
ryðzÞ ¼ r2 � z Sinð/ve=2Þ ð7Þ
A shallow shell is characterized by its middle
surface which is defined by the equation (Leissa et al.
1984),
z ¼ �0:5½ðx2=rxÞ þ ð2xy=rxyÞ þ ðy2=ryÞ� ð8Þ
Considering non-dimensional coordinate (1, v) sys-
tem, the radius of twist (rxy), length (Lo) of the shell
and twist angle (W) are related as (Liew et al. 1994)
1 ¼ x
Lov ¼ y
bof ð1; vÞ ¼ s tanð/ve=2Þ
ryð1; vÞrxy ¼ �Lo= tanw
ð9Þ
Hamilton’s principle applicable to non-conservative
system can be expressed as,
dH ¼Zt2
t1
½dT � dU � dW �dt ¼ 0 ð10Þ
where T, U and W are total kinetic energy, total strain
energy and total potential of the applied load, respec-
tively wherein dW = 0 for free vibration analysis. The
Lagrange’s equation of motion is given by
d
dt
oLf
ode
� �� oLf
ode
� �fFeg ð11Þ
where {Fe} is the applied external force vector of an
element and Lf is the Lagrangian function. Substituting
Lf = T - U in Lagrange’s equation incorporating the
corresponding expressions for T and U, the dynamic
equilibrium equation for each element in the following
form (Karmakar and Sinha 2001)
½Me�f€deg þ ½Ce�f _deg þ ð½Ke� þ ½Kre� � ½KRe�Þfdeg¼ fFceg þ fFeg ð12Þ
where {de} is the element displacement vector, [Me] is
the element mass matrix, [Ce] is the element Coriolis
matrix (skew symmetric), [Ke] is the element stiffness
matrix, [Kre] is the element geometric stiffness matrix,
[KRe] is the element rotational stiffness matrix (sym-
metric and positive definite), {Fe} is the element force
vector, {Fce}is the element centrifugal force vector.
After assembling all the element matrices and the force
vectors with respect to common global coordinates, the
resulting equilibrium equation of the structure becomes
½M�f€dg þ ½C�f _dg þ ð½K� þ ½Kr� � ½KR�Þ fdg¼ fFðX2Þg þ fFg ð13Þ
where [M], [C], [K], [Kr] and [KR] are the respective
global matrices, {F(X2)} is the vector of nodal
equivalent centrifugal forces, {F} is the global vector
of externally applied load and {d} is the global
displacement vector. For moderate rotational speed,
Coriolis matrix and rotational stiffness matrix are
neglected and finally the dynamic equilibrium equa-
tion in global form in absence of externally applied
load ({F}) reduces to (Sreenivasamurthy and Ramam-
urti 1981)
½M�f€dg þ ð½K� þ ½Kr�Þfdg ¼ fFðX2Þg ð14Þ
where [Kr] depends on initial stress distribution and is
obtained by the iterative procedure upon solving
ð½K� þ ½Kr�Þfdg ¼ fFðX2Þg ð15Þ
The natural frequencies (xn) are determined from
the standard eigenvalue problem (Bathe 1990) which
is represented below and is solved by QR iteration
algorithm
Xneq
� �df g ¼ k df g ð16Þ
where; Xneq
� �¼ K½ � þ Kr½ �ð Þ�1
M½ � ð17Þ
k ¼ 1=x2n ð18Þ
428 S. Dey et al.
123
An eight noded isoparametric quadratic element
with five degrees of freedom at each node (three
translation and two rotations) is employed wherein the
shape functions (Sj) are considered as function of local
natural coordinates of the element (1, v). Here
vj = ? 1 for nodes 2, 3, 6 and vj = -1 for nodes 1,
4, 8 and 1j = ?1 for nodes 3, 4, 7 and 1j = -1 for
nodes 1, 2, 5 as furnished in Fig. 2a. The variation of
volume fraction through the thickness with different
volume fraction exponents for Al–ZrO2 is illustrated
in Fig. 2b, c which indicates the variation of ceramic
and metallic composition (Vc, Vm) across the thick-
ness parameter (z/t) as per power law of functionally
graded material as furnished in Eq. (3). By virtue of
FGM properties, linear variation of volume fractions
for both ceramic and metallic composition is only
observed for N = 1.
3 Results and discussion
Parametric studies are carried out with respect to
rotational speed, twist angle on natural frequencies of
functionally graded rotating shallow conical shells.
Non-dimensional frequencies for conical shells
(rx = a) having square plan-form (Lo/bo = 1), cur-
vature ratio (bo/ry) of 0.5 and thickness ratio (s/t) of
1000 are obtained corresponding to non-dimensional
speeds of rotation, X = (X0/xo) = 0.0, 0.5 and 1.0,
Fig. 2 a Finite element discretisation of 8 9 8 mesh on plan
area with elements, node numbers and the natural coordinates
considering isoparametric element. b Variation of volume
fraction of ceramic (Vc) with respect to the thickness parameter
(z/t) of Al–ZrO2 FG conical shell. cVariation of volume fraction
of metal (Vm) with respect to the thickness parameter (z/t) of
Al–ZrO2 FG conical shell
Effect of twist and rotation on vibration 429
123
considering three different angles of twist of conical
shells, namely w = 15�, 30� and 45�, in addition to theuntwisted one (w = 0�). The variation of volume
fractions (Vc, Vm) through the thickness for Al–ZrO2
are furnished in Fig. 2b, c. An eight noded isopara-
metric quadratic element with five degrees of freedom
at each node is considered for finite element formu-
lation. Material properties of functionally graded
materials (Loy et al. 1999; Zhao and Liew 2011) are
adopted as stipulated in Table 1.
Even though the main application of FG material is
found in the area of thermal field where the temper-
ature gradient is significantly high, the present study is
aimed to address the combined effect of twist and
rotational speeds on the free vibration characteristics
of the functionally graded conical shell which is
idealized as turbomachinery blades. The constant
uniform temperature (300 K) is considered as the
temperature transient of the blades during engine
startup (i.e., at the beginning from uniform room
temperature) and the blades are initially exposed to
operate in external environment. Computer codes are
developed based on present finite element method.
The numerical results obtained are compared and
validated with the results of published literature (Zhao
and Liew 2011; Sreenivasamurthy and Ramamurti
1981) as furnished in Tables 2 and 3. Table 2 presents
the validation of non-dimensional fundamental natural
frequencies of functionally graded Al/ZrO2 conical
shells (Zhao and Liew 2011), while the effect of
rotation on non-dimensional fundamental natural
frequencies of an isotropic rotating cantilever plate
(Sreenivasamurthy and Ramamurti 1981) is furnished
in Table 3. The comparative study shows an excellent
agreement with the previously published results and
hence it demonstrates the capability of the computer
codes developed and proves the accuracy of the
analyses. The predictive capability of the computer
programs in respect of FG conical shells with rotating
and twisted models is thus confirmed. Convergence
studies are also performed to determine the converged
mesh size as furnished in Table 2. It is also observed
from the convergence study that uniform mesh
divisions of (6 9 6), (8 9 8) and (10 9 10) give
closer values to the results obtained by Zhao and Liew
et al. (2011) considering the complete plan-form of the
shell and provide nearly equal results the difference
being around one percent (1 %). To save the computer
time, (6 9 6) mesh division was chosen which gave
closest results to the benchmark results (Zhao and
Liew 2011). It is to be noted that the numerical results
depicted in Table 2 indicates that the accuracy of
numerical results of increased mesh size (viz.
10 9 10) is adverse compared to that of the coarse
mesh size (viz. 6 9 6). In conformity of the same, the
percentage difference of the results obtained by
(8 9 8) mesh and (10 9 0) mesh with respect to the
benchmark values were found much more than that of
Table 1 Properties of the FGM components at 300 K
Material Properties
E (N/m2) t q (Kg/m3)
Ti–6Al–4V 105.700 9 109 0.298000 4,429
Aluminum oxide
(Al2O3)
320.200 9 109 0.260000 3,750
Stainless steel (SS) 207.788 9 109 0.317756 8,166
Nickel (Ni) 205.098 9 109 0.310000 8,900
Aluminum (Al) 70.000 9 109 0.300000 2,707
Zirconia (ZrO2) 151.000 9 109 0.300000 3,000
Table 2 Convergence study for NDFF [x = xn�a2 H(qc/Ect2)/2p,
for Al/ZrO2 functionally graded conical shells, considering
r1 = 0.2 m, t = 0.01 m, Lo = 0.8 m, /ve = 30�, /o = 120�
N Zhao and Liew (Zhao and
Liew 2011)
Present FEM
6 9 6 8 9 8 10 9 10
0.0 1.3666 1.3608 1.3449 1.2356
0.5 1.2486 1.2387 1.2271 1.0848
1.0 1.1893 1.1792 1.1641 1.1038
5.0 1.0737 1.0578 1.0492 0.9910
10.0 1.0404 1.0222 1.0102 0.9566
Table 3 Non-dimensional fundamental natural frequencies
[x = xnLo2H(qt/D)] of an isotropic rotating cantilever plate,
Lo/b = 1, t/Lo = 0.12, D = Et3/{12(1 - m2)}, m = 0.3
Non-dimensional
speed (X)Present
FEM
Sreenivasamurthy and
Ramamurti (Sreenivasamurthy
and Ramamurti 1981)
0.0 3.4174 3.4368
0.2 3.4933 3.5185
0.4 3.7110 3.7528
0.6 4.0458 4.1287
0.8 4.4690 4.5678
1.0 4.9549 5.0916
430 S. Dey et al.
123
(6 9 6) mesh. Hence, (8 9 8) mesh and (10 9 10)
mesh were excluded in order to achieve better
accuracy of results and to reduce computational cost.
The lower mesh size (6 9 6) consisting of 36 elements
and 133 nodes, has been used for the analysis due to
computational efficiency. The total number of degrees
of freedom involved in the computation is 665 as each
node of the isoparametric quadratic element is having
five degrees of freedom comprising of three transla-
tions and two rotations. Table 4 represents the non-
dimensional fundamental natural frequencies of sim-
ply supported aluminum-zirconia FG plate (Prady-
umna and Bandyopadhyay 2008; Matsunaga 2008;
Ferreira et al. 2006).
The parametric study is carried out for functionally
graded (stainless steel–nickel and Ti–6Al–4V and
Aluminum oxide) conical shells for different twist
angles at both stationary and rotating condition, as
furnished in Tables 5 and 6. At stationary condition,
non-dimensional fundamental natural frequencies are
identified to attain lowest value at w = 0� and
gradually are found to increase to maximum value at
w = 45�. For both stationary and rotational cases,
non-dimensional fundamental frequencies are found
to decrease with increase of material property graded
index (N), irrespective of twist angle. Due to decrease
in elastic stiffness at stationary condition, minimum
value of non-dimensional fundamental natural fre-
quencies of stainless steel–nickel FG conical shells is
Table 4 NDFF [x = xnt2H(qm/Em)] of simply supported aluminum–zirconia functionally graded plate a = b = 1 m, N = 1 (For
aluminum: Em = 70 GPa, mm = 0.3, and qm = 2,707 kg/m3 and for zirconia: Ec = 200 GPa, mc = 0.3 and qc = 2,702 kg/m3)
a/h Present FEM Pradyumna and Bandyopadhyay
(Pradyumna and Bandyopadhyay 2008)
Matsunaga
(Matsunaga 2008)
Ferreira et al.
(Ferreira et al. 2006)
5 0.2235 0.2257 0.2285 0.2188
10 0.0608 0.0613 0.0618 0.0592
20 0.0156 0.0157 0.0158 0.0147
Table 5 NDFF and NDSF [x = xnLo2H(q/E1t
2)] of rotating stainless steel–nickel functionally graded conical shells for various
twist angles considering Lo/s = 0.7, r1 = 0.2 m, Lo = 0.8 m, t = 0.01 m, /o = 45�, /ve = 20�
W N Non-dimensional natural frequencies
NDFF NDSF
X = 0 X = 0.5 X = 1.0 X = 0 X = 0.5 X = 1.0
0� 0 0.8121 1.0097 1.4504 2.0889 2.3701 3.0542
1 0.7916 0.9915 1.4326 2.0364 2.3212 3.0078
5 0.7792 0.9714 1.3982 2.0046 2.2782 2.9415
10 0.7763 0.9666 1.3899 1.9970 2.2679 2.9257
15� 0 0.8536 1.0626 1.5229 2.0891 2.3979 3.1347
1 0.8321 1.0434 1.5041 2.0366 2.3493 3.0883
5 0.8190 1.0222 1.4679 2.0047 2.3052 3.0194
10 0.8159 1.0171 1.4593 1.9971 2.2946 3.0029
30� 0 1.0525 1.3198 1.8759 2.0578 2.5110 3.4848
1 1.0260 1.2960 1.8521 2.0061 2.4642 3.4377
5 1.0097 1.2694 1.8076 1.9747 2.4151 3.3579
10 1.0059 1.2632 1.7972 1.9672 2.4035 3.3390
45� 0 1.4253 1.8202 2.5483 1.9931 2.7358 4.0867
1 1.3896 1.7873 2.5146 1.9431 2.6917 4.0364
5 1.3673 1.7504 2.4549 1.9126 2.6333 3.9389
10 1.3622 1.7419 2.4410 1.9054 2.6196 3.9161
Effect of twist and rotation on vibration 431
123
identified at N = 10 while due to increase in elastic
stiffness the maximum value is found at N = 0
irrespective of twist angle.
Considering N = 0, 1, 5 and 10, the differences
between maximum and minimum value of non-
dimensional fundamental natural frequencies at sta-
tionary condition for stainless steel–nickel FG conical
shells are found 4.41, 4.42, 4.43 and 4.44 % for
w = 0�, 15�, 30� and 45�, respectively. Here, the
maximum difference (in percentage) between maxi-
mum and minimum NDFF values at stationary con-
dition for stainless steel–nickel FG conical shell is
identified for w = 45� while the minimum difference
(in percentage) is observed at w = 0�.In contrast, the differences between maximum and
minimum values of non-dimensional fundamental
natural frequencies at lower rotational speeds
(X = 0.5) of stainless steel–nickel FG conical shells
are identified 4.27, 4.28, 4.29 and 4.30 % for w = 0�,15�, 30� and 45�, respectively, while at higher
rotational speeds the differences between maximum
and minimum NDFF values are found 4.17, 4.18, 4.20
and 4.21 % forw = 0�, 15�, 30� and 45�, respectively.
Therefore, it is noted that the percentage difference
between maximum and minimum value of fundamen-
tal natural frequencies obtained from the analyses of
four different material property graded index (N = 0,
1, 5 and 10) for stainless steel–nickel FG conical shell
decreases as the twist angle increases. Hence it leads to
the fact that for untwisted conical shells, material
property graded index has pronounced effect on non-
dimensional natural frequencies. At rotating condi-
tion, NDFF and NDSF are found to reduce with
increase of material property graded index. Table 6
stipulates the trends of NDFF and NDSF of stainless
steel–nickel functionally graded conical shells for
different twist angles at rotating condition. At both
lower rotational speeds (X = 0.5) and higher rota-
tional speeds (X = 1.0), non-dimensional fundamen-
tal natural frequencies and non-dimensional second
natural frequencies at twisted cases are found pre-
dominantly higher compared to respective NDFF and
NDSF values of untwisted cases. The trend of relative
frequencies (ratio of rotating natural frequency and
stationary natural frequency) for both SS–Ni and
Ti–6Al–4V and Aluminum oxide FG configuration at
Table 6 NDFF and NDSF [x = xnL2H(q/E1h
2)] of rotating Ti–6Al–4V and aluminum oxide (Al2O3) functionally graded conical
shells for various twist angles considering Lo/s = 0.7, r1 = 0.2 m, Lo = 0.8 m, t = 0.01 m, /o = 45�, /ve = 20�
W N Mode
NDFF NDSF
X = 0 X = 0.5 X = 1.0 X = 0 X = 0.5 X = 1.0
0� 0 2.6948 3.3537 4.8139 6.9362 7.8761 10.1525
1 2.0107 2.6506 3.6314 5.1641 6.2044 7.9374
5 1.6874 2.0767 3.0300 4.3399 4.9093 6.3692
10 1.5967 1.9830 2.8640 4.1102 4.6651 6.0284
15� 0 2.8328 3.5291 5.0547 6.9368 7.9585 10.4209
1 2.1254 2.9653 3.8425 5.1655 6.3297 8.0864
5 1.7708 2.1852 3.1796 4.3401 4.9648 6.5331
10 1.6720 2.0803 2.9963 4.1110 4.7167 6.1765
30� 0 3.4961 4.3854 6.2307 6.8332 8.3447 11.5930
1 2.6536 3.8049 4.8216 5.0899 6.6928 8.9356
5 2.1613 2.6968 3.8805 4.2736 5.1850 7.2339
10 2.0289 2.5223 3.6363 4.0463 4.7324 6.8114
45� 0 4.7372 5.9051 8.4696 6.6211 9.0916 13.5957
1 3.6219 5.2735 6.6290 4.9389 6.9492 10.6233
5 2.9028 3.7335 5.2286 4.1297 5.6734 8.4363
10 2.7133 3.4425 4.8727 3.9071 5.2637 7.8817
432 S. Dey et al.
123
X = 0.5, 1.0 are furnished in Figs. 3 and 4, respec-
tively. The centrifugal stiffening effect (i.e., increase
of structural stiffness with increase of rotational
speed) is predominantly observed for both NDFF
and NDSF values irrespective of twist angle for both
stainless steel–nickel and Ti–6Al–4V and Aluminum
oxide functionally graded conical shells. This phe-
nomenon can be attributed to the fact that the variation
in natural frequencies is occurred due to contribution
of geometric stiffness matrix which changes the
effective stiffness matrix of the functionally graded
material. The geometric stiffness matrix depends on
initial stress distribution and influences the elastic
stiffness differently, yielding to different frequency
parameters irrespective of twist angle.
For both FG conical shells, it is observed that
relative frequencies [ratio of natural frequency at
rotating condition (i.e., X = 0.5 or 1.0) and respective
natural frequency at stationary condition (i.e.,
X = 0.0)] corresponding to untwisted cases are
mostly in the lower range compared to those of
twisted ones in respect of lower (X = 0.5) and higher
(X = 1.0) rotational speeds. In other words, rotating
effect is more pronounced for twisted shell in com-
parison to untwisted one. The centrifugal stiffening
effect (i.e., increase of structural stiffness with
increase of rotational speeds) with respect to both
NDFF and NDSF are observed for Ti–6Al–4V and
Aluminum oxide FG configuration irrespective of
twist angle. Considering w = 0�, 15�, 30� and 45�, thedifferences between maximum and minimum value of
non-dimensional fundamental natural frequencies at
stationary condition for Ti–6Al–4V and Aluminum
oxide FG conical shells are found 43.1, 44.5, 41.9 and
Fig. 3 Variation of relative frequencies (NDFF and NDSF)
with material property graded index (N) at rotating conditions
(X = 0.5 and 1.0) for stainless steel–nickel functionally graded
conical shells with various twist angles considering Lo/s = 0.7,
r1 = 0.2 m, Lo = 0.8 m, t = 0.01 m, /o = 45�, /ve = 20�
Effect of twist and rotation on vibration 433
123
41.2 % for N = 0, 1, 5 and 10, respectively. It is also
observed that the maximum value of relative fre-
quency with respect to both NDFF and NDSF are
obtained at w = 15� for N = 1.
Functionally graded conical shells, because of the
possible different volume fractions for the constituent
materials, also exhibit interesting frequency charac-
teristics that are not seen in homogeneous isotropic
conical shells. The influence of the value of material
property graded index (N), which affects the constit-
uent volume fraction, can be seen from the above
tables (Tables 5, 6). As the value of ‘N’ increases, the
natural frequency decreases irrespective of twist
angle. This corroborates with the results obtained by
Loy et al. (1999). Due to variation in elastic stiffness of
FG conical shell panels with certain material graded
index or power law exponent, the variation in NDSF
has lesser pronounced effect compared to NDFF and
subsequently observed to decrease the effect as the
twist angle increases.
4 Mode shapes
The mode shapes corresponding to the natural fre-
quencies are furnished for various twist angles
(w = 0�, 15�, 30� and 45�) and rotational speeds
(X = 0.0, 0.5 and 1.0), considering stainless steel–
nickel and Ti–6Al–4V and Aluminum oxide FG
conical shells, as shown in Figs. 5 and 6, respectively.
The fundamental frequency corresponds to the first
torsional mode for all cases. It is identified that the
Fig. 4 Variation of relative frequencies (NDFF and NDSF)
with material property graded index (N) at rotating conditions
(X = 0.5 and 1.0) for Ti–6Al–4V and Aluminum oxide (Al2O3)
functionally graded conical shells for various twist angles
considering Lo/s = 0.7, r1 = 0.2 m, Lo = 0.8 m, t = 0.01 m,
/o = 45�, /ve = 20�
434 S. Dey et al.
123
symmetry modes are absent when twist angle is non-
zero. The first span wise bending is observed at higher
rotational speeds (X = 1.0) for stainless steel–nickel
FG conical shells corresponding to both NDFF and
NDSF, irrespective of twist angle. In contrast, for
Ti–6Al–4V and Aluminum oxide FG conical shells,
the first span wise bending is restricted to untwisted
cases of both lower rotational speeds (X = 0.5) and
higher rotational speeds (X = 1.0) corresponding to
NDFF. Due to variation of effective elastic stiffness of
FG conical shell, the centrifugal stiffening effect is
found to comply with the mode shapes of both twisted
and untwisted cases. At higher rotational speeds
(X = 1.0), the influence of geometric stiffness matrix
is influenced to transform the mode shape of both
fundamental and second natural frequency of SS–Ni FG
conical shell from torsional mode to first spanwise
bending mode. In contrast, for Ti–6Al–4V and Alumi-
num oxide FG conical shells, due to dominant effect of
twist angle say, at w = 30� and 45�, the mode shape of
fundamental natural frequency remain under torsional
mode throughout at higher rotational speeds (X = 1.0)
while due to predominance of geometric stiffness
matrix, first spanwise bendingmode is identified only at
twist anglew = 0� and 15� for higher rotational speeds.
5 Conclusions
This paper illustrated the combined effect of twist and
rotational speeds for functionally graded conical
shells. A validation of the analysis has been carried
out by comparing results with those in the literature
and has found to be accurate. Non-dimensional natural
frequencies of both stainless steel–nickel and Ti–6Al–4V
and Aluminum oxide FG conical shells are consis-
tently observed to increase with rise of rotational
speeds irrespective of twist angle. At both stationary
and rotating condition, the NDFF and NDSF are found
to decrease with increase of material property graded
index (viz. N = 0, 1, 5, 10) corresponding to both
twisted and untwisted cases. For both stationary and
rotating condition, it is observed that NDFF and NDSF
are obtained to increase with increase of twist and the
Fig. 5 Effect of twist and rotation on NDFF and NDSF for mode shapes for SS–Ni FG conical shells for different twist angles,
considering r1 = 0.2 m, t = 0.01 m, N = 1, s/t = 1000, Lo/s = 0.7, /o = 45�, /ve = 20�
Effect of twist and rotation on vibration 435
123
rotating effect is found to be more pronounced for
twisted FG conical shells compared to untwisted one.
The centrifugal stiffening effect is predominantly
observed for both NDFF and NDSF values irrespec-
tive of twist angle for both stainless steel–nickel and
Ti–6Al–4V and aluminium oxide functionally graded
conical shells. It is noted from the mode shapes of the
FG conical shells that the fundamental and second
natural frequencies at stationary condition corre-
sponds to the first torsional mode for all cases
irrespective of twist angle. It is identified that the
symmetry modes are absent when twist angle is non-
zero. The non-dimensional frequencies obtained are
the first known results of the type of analyses carried
out here and the results could serve as reference
solutions for future investigators.
References
Ahmad, M., Naeem, M.N.: Vibration characteristics of rotating
FGM circular cylindrical shells using wave propagation
method. Eur. J. Sci. Res. 36, 184–235 (2009)
Bardell, N.S., Dunsdon, J.M., Langley, R.S.: Free vibration of
thin, isotropic, open, conical panels. J. Sound and Vibration
217, 297–320 (1998)
Bathe, K.J.: Finite Element Procedures in Engineering Analysis.
PHI, New Delhi (1990)
Bhangale, R.K., Ganesan, K., Padmanabhan, C.: Linear thermo
elastic buckling and free vibration behaviour of function-
ally graded truncated conical shells. J. Sound and Vibration
292, 341–371 (2006)
Birman, V., Bert, C.W.: Buckling and postbuckling of com-
posite plates and shells subjected to elevated temperature.
ASME J. Appl. Mech. 60, 514–519 (1993)
Eslami, M.R., Ziaii, A.R., Chorbarpour, A.: Thermo-elastic
buckling of thin cylindrical shells based on improved sta-
bility equations. J. Therm. Stresses 19(4), 299–315 (1996)
Ferreira, A.J.M., Batra, R.C., Roque, C.M.C., Qian, L.F., Jorge,
R.M.N.: Natural frequencies of functionally graded plates
by meshless method. Composite Struct. 75, 593–600 (2006)Fukui, Y., Yamanaka, N., Wakashima, K.: The stresses and
strains in a thick-walled tube for functionally graded under
Fig. 6 Effect of twist and rotation on NDFF and NDSF for mode shapes of Ti–6Al–4V and Aluminum oxide (Al2O3) FG conical shells
for different twist angles, considering r1 = 0.2 m, t = 0.01 m, N = 1, s/t = 1000, Lo/s = 0.7, /o = 45�, /ve = 20�
436 S. Dey et al.
123
uniform thermal loading. Int. J. Jpn. Soc. Mech. Eng. Ser.
A 36, 156–162 (1993)
Ghafoori, E., Asghari, M.: Three-dimensional elasticity analysis
of functionally graded rotating cylinders with variable
thickness profile. Proc IME C J. Mech Eng Sci 226,585–593 (2012)
Gotsis, P.K., Guptill, J.D.: Fiber-composite thin shells subjected
to thermal buckling loads. Comput. Struct. 53(6),1263–1274 (1994)
He, X.Q., Liew, K.M., Ng, T.Y., Sivashanker, S.: A FEMmodel
for the active control of curved FGM shells using piezo-
electric sensor/actuator layers. Int J. Numer Methods Eng
54, 853–870 (2002)
Heydarpour, Y., Malekzadeh, P., Golbahar Haghighi, M.R.,
Vaghefi, M.: Thermoelastic analysis of rotating laminated
functionally graded cylindrical shells using layerwise dif-
ferential quadrature method. Acta Mech. 223, 81–93
(2012)
Hua, L.: Frequency analysis of rotating truncated circular
orthotropic conical shells with different boundary condi-
tions. Compos Sci Tech 60, 2945–2955 (2000)
Karmakar, A., Sinha, P.K.: Failure analysis of laminated com-
posite pretwisted rotating plates. J. Reinforced Plastics and
Composites 20, 1326–1357 (2001)
Leissa, A.W., Lee, J.K., Wang, A.J.: Vibrations of Twisted
Rotating Blades. J. Vibration, Acoustics, Stress and Reli-
ability in Design, Trans., ASME 106(2), 251–257 (1984)
Liew, K.M., Lim, C.M., Ong, L.S.: Vibration of pretwisted
cantilever shallow conical shells. I. J. Solids Structures 31,2463–2474 (1994)
Liew, K.M., He, X.Q., Kitipornchai, S.: Finite element method
for the feedback control of FGM shells in the frequency
domain via piezoelectric sensors and actuators. Comput.
Methods Appl. Mech. Eng. 193, 257–273 (2004)
Loy, C.T., Lam, J.N., Reddy, J.N.: Vibration of functionally
graded cylindrical shells. Int. J. Mech. Sci. 41, 309–324(1999)
Malekzadeh, P., Heydarpour, Y.: Free vibration analysis of
rotating functionally graded cylindrical shells in thermal
environment. Composite Structure 94, 2971–2981 (2012)
Matsunaga, H.: Free vibration and stability of functionally
graded plates according to a 2-D higher-order deformation
theory. Composite Struct. 82, 499–512 (2008)
Ng, T.Y., Lam, K.Y., Liew, K.M., Reddy, J.N.: Dynamic sta-
bility analysis of functionally graded cylindrical shells
under periodic axial loading. Int. J. Solids Struct. 38,1295–1309 (2001)
Ng, T.Y., He, X.Q., Liew, K.M.: Finite element modeling of
active control of functionally graded shells in frequency
domain via piezoelectric sensors and actuators. Comput.
Mech. 28, 1–9 (2002)
Noda, N.: Thermal stresses in functionally graded materials.
J. Thermal Stress 22, 477–512 (1999)
Pradhan, S.C., Loy, C.T., Lam, K.Y., Reddy, J.N.: Vibration
characteristics of functionally graded cylindrical shells
under various boundary conditions. Appl. Acoust. 61(1),119–129 (2000)
Pradyumna, S., Bandyopadhyay, J.N.: Free vibration analysis of
functionally graded curved panels using a higher-order
finite element formulation. J. Sound and Vibration 318,176–192 (2008)
Praveen, G.N., Reddy, J.N.: Nonlinear transient thermoelastic
analysis of functionally graded ceramic–metal plates. Int.
J. Solids Struct. 33, 4457–4476 (1998)
Setoodeh, A.R., Tahani, M., Selahi, E.: Transient dynamic and
free vibration analysis of functionally graded truncated
conical shells with non-uniform thickness subjected to
mechanical shock loading. Compos Part B: Eng. 43,2161–2171 (2012)
Shu, C.: Free vibration analysis of composite laminated conical
shells by generalized differential quadrature. J. Sound and
Vibration. 194, 587–604 (1996)
Sofiyev, A.H.: Thermoelastic stability of functionally graded
truncated conical shells. Composite Struct. 77, 56–65
(2007)
Sofiyev, A.H.: The vibration and stability behavior of freely
supported FGM conical shells subjected to external pres-
sure. Compos. Struct. 89, 356–366 (2009)
Sreenivasamurthy, S., Ramamurti, V.: Coriolis Effect on the
Vibration of Flat Rotating low Aspect Ratio Cantilever
Plates. Journal of Strain Analysis 16, 97–106 (1981)
Tauchert, T.R.: Thermally induced flexure, buckling, and
vibration of plates. Appl. Mech. Rev. 44, 347–360 (1991)
Thornton, E.A.: Thermal buckling of plates and shells. Appl.
Mech. Rev. 46(10), 485–506 (1993)
Tornabene, F., Viola, E., Inman, D.J.: 2-D differential quadra-
ture solution for vibration analysis of functionally graded
conical, cylindrical shell and annular plate structures.
J. Sound and Vibration 328, 259–290 (2009)
Touloukian, Y.S.: Thermophysical properties of high tempera-
ture solid materials. McMillan, New York (1967)
Zhao, X., Liew, K.M.: Free vibration analysis of functionally
graded conical shell panels by a meshless method. Com-
posite Struct. 93, 649–664 (2011)
Effect of twist and rotation on vibration 437
123