# Thermal Effect on Vibration of Tapered Rectangular Plate

Authors: Anupam Khanna; Ashish Singhal
Journal: Journal of Structures
 (2014)

---

## Abstract

A mathematical model is constructed to help the engineers in designing various mechanical structures mostly used in satellite and aeronautical engineering. In the present model, vibration of rectangular plate with nonuniform thickness is discussed. Temperature variations are considered biparabolic, that is, parabolic inx-direction and parabolic in y-direction. The fourth-order differential equation of the motion is solved by Rayleigh Ritz method for three different boundary conditions around the boundary of plate. Numerical values of frequencies for the first two modes of vibration are presented in tabular form for different values of thermal gradient, taper constants, and aspect ratio.


---

## Body

## 1. Introduction

The structures are designed to support the high speed engines and turbines subjected to the vibration. Due to faulty design, there is unbalance in the engines which causes excessive and unpleasant stresses because of vibration. In the field of science and technology, it is preferred to design large machines and structures for smooth operations with controlled vibration. Recent studies in the field of vibrational behavior create a huge interest for scientists and engineers in the designs and constructions of complex systems or structures such as ships, submarines, aircrafts, launch vehicles, missiles, and satellites.Tapered plates or plates of nonuniform thickness are commonly used in many engineering applications such as nuclear engineering, aeronautical engineering, and chemical plants, under the influence of elevated temperature to control high vibration.A collection of research papers on vibration of plates with different shapes and boundary conditions is given by Leissa [1] in his monograph. Leissa [2] discussed different models on free vibration of rectangular plates. Jain and Soni [3] analyzed the free vibrations of rectangular plates with parabolically varying thickness. Tomar and Gupta [4] studied the effect of thermal gradient on the vibration of a rectangular plate with bidirectional variation in thickness. Leissa [5] investigated the effect of thermal gradient on the vibration of parallelogram plate with bidirectional thickness variation in both directions. Singh and Chakraverty [6] studied the transverse vibration of circular and elliptical plates with variable thickness. Sharma and Chand [7] analyzed the vibrations in transversely isotropic plates due to suddenly punched hole. Leissa [8] discussed the historical bases of the Rayleigh and Ritz methods for the vibrations of plates. Chakraverty et al. [9] studied the flexural vibrations of nonhomogeneous elliptical plates. Gupta and Khanna [10] analyzed time period and deflection for the first two modes of vibrations of viscoelastic rectangular plate with linear thickness variations in both directions. Gupta and Khanna [11] had evaluated time period and deflection for the first two modes of vibration of viscoelastic rectangular plate for biparabolic thickness variation. Lal et al. [12] analyzed the transverse vibrations of nonhomogeneous rectangular plates of uniform thickness using boundary characteristic orthogonal polynomials. Khanna et al. [13] discussed the effect of nonhomogeneity on thermally induced vibrations of nonhomogeneous rectangular plate.In the present model, the authors analyzed vibration of rectangular plate with nonuniform thickness (Figure1). Nonuniformity in thickness is considered biparabolic. The authors calculated frequency for the first two modes of vibration at various values of thermal gradient (due to variation in temperature), taper constant (due to nonuniform thickness), and aspect ratio (ratio of length and breadth of the plate) for three different boundary conditions, that is, C-C-C-C, S-S-S-S, and C-S-C-S where C and S stand for clamped and simply supported, respectively. Results are given in the tabular form.Figure 1
Rectangular plate with 2-directional thickness variation.


## 2. Materials and Methods

### 2.1. Analysis of Motion

Fourth-order differential equation of motion for isotropic rectangular plate in Cartesian coordinate is [14]
(1)
[
D
1
(
∂
4
W
∂
x
4
+
2
∂
4
W
∂
x
2
∂
y
2
+
∂
4
W
∂
y
4
)
+
2
∂
D
1
∂
x
(
∂
3
W
∂
x
3
+
∂
3
W
∂
x
∂
y
2
)
+
2
∂
D
1
∂
y
(
∂
3
W
∂
y
3
+
∂
3
W
∂
y
∂
x
2
)
+
∂
2
D
1
∂
x
2
(
∂
2
W
∂
x
2
+
υ
∂
2
W
∂
y
2
)
+
∂
2
D
1
∂
y
2
(
∂
2
W
∂
y
2
+
υ
∂
2
W
∂
x
2
)
+
2
(
1
-
υ
)
∂
2
D
1
∂
x
∂
y
∂
2
W
∂
x
∂
y
]
-
ρ
p
2
d
W
=
0
,
where D
1
=
E
d
3
/
(
12
(
1
-
υ
2
)
) is the flexural rigidity of the plate’s material, W
=
W
(
x
,
y
) is the deflection function, υ is Poisson ratio, ρ is mass per unit volume, and d is thickness of the plate.It is assumed that thickness of the rectangular plate varies parabolically in both directions [11]; that is,
(2)
d
=
d
0
(
1
+
β
1
x
2
a
2
)
(
1
+
β
2
y
2
b
2
)
,
where a and b are length and breadth of rectangular plate, respectively, and β
1 and β
2 are taper parameters in x-direction and y-direction, respectively.The authors assumed biparabolic temperature variations as(3)
τ
=
τ
0
(
1
-
x
2
a
2
)
(
1
-
y
2
b
2
)
,
where τ denotes the temperature excess above the reference temperature at any point on the plate and τ
0 denotes the temperature excess above the reference temperature at x
=
y
=
0.The temperature dependence of the modulus of elasticity for most of engineering materials can be expressed as follows [15]:
(4)
E
=
E
0
(
1
-
γ
τ
)
,
where E
0 is value of Young’s modulus at reference temperature, that is, τ
=
0, and γ is slope of variation of E and τ. On using (3) in (4), one obtains
(5)
E
=
E
0
[
1
-
α
(
1
-
x
2
a
2
)
(
1
-
y
2
b
2
)
]
,
where α
=
γ
τ
0 (0
≤
α
<
1) is thermal gradient.On substituting the values ofd and E from (2) and (5), respectively, expression of flexural rigidity (D
1) becomes
(6)
D
1
=
(
(
1
+
β
1
(
x
2
/
a
2
)
)
3
E
0
[
1
-
α
(
1
-
(
x
2
a
2
)
)
(
1
-
(
y
2
b
2
)
)
]
s
×
d
0
3
(
1
+
β
1
(
x
2
a
2
)
)
3
(
1
+
β
2
(
y
2
b
2
)
)
3
)
×
(
12
(
1
-
υ
2
)
)
-
1
.


### 2.2. Boundary Conditions and Corresponding Deflection Function

In this model, the authors discussed three different boundary conditions, that is, C-C-C-C, S-S-S-S, and C-S-C-S around the boundary of plate.To satisfy all these boundary conditions, deflection function (W) is taken as shown in Table 1.Table 1
Boundary conditions
Deflection function (W)
C-C-C-C
W
=
[
(
x
a
)
2
(
y
b
)
2
(
1
-
x
a
)
2
(
1
-
y
b
)
2
]
×
[
A
1
+
A
2
(
x
a
)
(
y
b
)
(
1
-
x
a
)
(
1
-
y
b
)
]
S-S-S-S
W
=
[
(
x
a
)
(
y
b
)
(
1
-
x
a
)
(
1
-
y
b
)
]
×
[
A
1
+
A
2
(
x
a
)
(
y
b
)
(
1
-
x
a
)
(
1
-
y
b
)
]
C-S-C-S
W
=
[
(
x
a
)
2
(
y
b
)
(
1
-
x
a
)
2
(
1
-
y
b
)
]
×
[
A
1
+
A
2
(
x
a
)
(
y
b
)
(
1
-
x
a
)
(
1
-
y
b
)
]


### 2.3. Solution of Frequency Equation

Rayleigh Ritz technique is applied to solve the frequency equation. In this method, one requires that the maximum strain energy(
S
E
) must be equal to the maximum kinetic energy (K
E). So it is necessary for the problem under consideration that [16]
(7)
δ
(
S
E
-
K
E
)
=
0
.
Now assume the nondimensional variables as
(8)
X
=
x
a
,
Y
=
y
a
,
W
-
=
W
a
,
d
-
=
d
a
.
The expressions for kinetic energy (K
E) and strain energy (S
E) are
(9)
K
E
=
1
2
ρ
p
2
d
0
-
a
5
×
∫
0
1
∫
0
b
/
a
[
(
1
+
β
1
X
2
)
(
1
+
β
2
a
2
b
2
Y
2
)
W
-
2
]
d
Y
d
X
,
S
E
=
E
0
d
0
-
3
a
3
24
(
1
-
υ
2
)
×
∫
0
1
∫
0
b
/
a
{
1
-
α
(
1
-
X
2
)
(
1
-
a
2
b
2
Y
2
)
}
×
(
1
+
β
1
X
2
)
3
(
1
+
β
2
a
2
b
2
Y
2
)
3
×
[
(
∂
2
W
-
∂
X
2
)
2
+
(
∂
2
W
-
∂
Y
2
)
2
+
2
υ
∂
2
W
-
∂
X
2
∂
2
W
-
∂
Y
2
+
2
(
1
-
υ
)
(
∂
2
W
-
∂
X
∂
Y
)
2
]
d
Y
d
X
On using (9) in (7), one gets
(10)
S
E
*
-
Ω
2
K
E
*
=
0
,
where
(11)
K
E
*
=
∫
0
1
∫
0
b
/
a
[
(
1
+
β
1
X
2
)
(
1
+
β
2
a
2
b
2
Y
2
)
W
-
2
]
d
Y
d
X
,
S
E
*
=
∫
0
1
∫
0
b
/
a
{
1
-
α
(
1
-
X
2
)
(
1
-
a
2
b
2
Y
2
)
}
×
(
1
+
β
1
X
2
)
3
(
1
+
β
2
a
2
b
2
Y
2
)
3
×
[
(
∂
2
W
-
∂
X
2
)
2
+
(
∂
2
W
-
∂
Y
2
)
2
+
2
υ
∂
2
W
-
∂
X
2
∂
2
W
-
∂
Y
2
+
2
(
1
-
υ
)
(
∂
2
W
-
∂
X
∂
Y
)
2
]
d
Y
d
X
.
Here, Ω
2
=
(
12
ρ
p
2
a
2
(
1
-
υ
2
)
)
/
(
E
0
d
0
-
2
) is a frequency parameter.Equation (10) consists of two unknown constants, that is, A
1 and A
2, arising due to the substitution of W. These two constants are to be determined as follows:
(12)
∂
∂
A
n
(
S
E
*
-
Ω
2
K
E
*
)
=
0
,
n
=
1
,
2
.
On simplifying (12), one gets
(13)
b
n
1
A
1
+
b
n
2
A
2
=
0
,
n
=
1
,
2
,
where b
n
1 and b
n
2 include parametric constant, that is, taper constants, thermal gradient, aspect ratio, and frequency parameter.For a nontrivial solution, determinant of the coefficient must be zero. So, one gets the frequency equation as(14)
|
b
11
b
12
b
21
b
22
|
=
0
.
With the help of (14), one can obtain a quadratic equation in Ω
2 from which the two values of frequency parameter for both the modes of vibration can be evaluated easily.


## 2.1. Analysis of Motion

Fourth-order differential equation of motion for isotropic rectangular plate in Cartesian coordinate is [14]
(1)
[
D
1
(
∂
4
W
∂
x
4
+
2
∂
4
W
∂
x
2
∂
y
2
+
∂
4
W
∂
y
4
)
+
2
∂
D
1
∂
x
(
∂
3
W
∂
x
3
+
∂
3
W
∂
x
∂
y
2
)
+
2
∂
D
1
∂
y
(
∂
3
W
∂
y
3
+
∂
3
W
∂
y
∂
x
2
)
+
∂
2
D
1
∂
x
2
(
∂
2
W
∂
x
2
+
υ
∂
2
W
∂
y
2
)
+
∂
2
D
1
∂
y
2
(
∂
2
W
∂
y
2
+
υ
∂
2
W
∂
x
2
)
+
2
(
1
-
υ
)
∂
2
D
1
∂
x
∂
y
∂
2
W
∂
x
∂
y
]
-
ρ
p
2
d
W
=
0
,
where D
1
=
E
d
3
/
(
12
(
1
-
υ
2
)
) is the flexural rigidity of the plate’s material, W
=
W
(
x
,
y
) is the deflection function, υ is Poisson ratio, ρ is mass per unit volume, and d is thickness of the plate.It is assumed that thickness of the rectangular plate varies parabolically in both directions [11]; that is,
(2)
d
=
d
0
(
1
+
β
1
x
2
a
2
)
(
1
+
β
2
y
2
b
2
)
,
where a and b are length and breadth of rectangular plate, respectively, and β
1 and β
2 are taper parameters in x-direction and y-direction, respectively.The authors assumed biparabolic temperature variations as(3)
τ
=
τ
0
(
1
-
x
2
a
2
)
(
1
-
y
2
b
2
)
,
where τ denotes the temperature excess above the reference temperature at any point on the plate and τ
0 denotes the temperature excess above the reference temperature at x
=
y
=
0.The temperature dependence of the modulus of elasticity for most of engineering materials can be expressed as follows [15]:
(4)
E
=
E
0
(
1
-
γ
τ
)
,
where E
0 is value of Young’s modulus at reference temperature, that is, τ
=
0, and γ is slope of variation of E and τ. On using (3) in (4), one obtains
(5)
E
=
E
0
[
1
-
α
(
1
-
x
2
a
2
)
(
1
-
y
2
b
2
)
]
,
where α
=
γ
τ
0 (0
≤
α
<
1) is thermal gradient.On substituting the values ofd and E from (2) and (5), respectively, expression of flexural rigidity (D
1) becomes
(6)
D
1
=
(
(
1
+
β
1
(
x
2
/
a
2
)
)
3
E
0
[
1
-
α
(
1
-
(
x
2
a
2
)
)
(
1
-
(
y
2
b
2
)
)
]
s
×
d
0
3
(
1
+
β
1
(
x
2
a
2
)
)
3
(
1
+
β
2
(
y
2
b
2
)
)
3
)
×
(
12
(
1
-
υ
2
)
)
-
1
.


## 2.2. Boundary Conditions and Corresponding Deflection Function

In this model, the authors discussed three different boundary conditions, that is, C-C-C-C, S-S-S-S, and C-S-C-S around the boundary of plate.To satisfy all these boundary conditions, deflection function (W) is taken as shown in Table 1.Table 1
Boundary conditions
Deflection function (W)
C-C-C-C
W
=
[
(
x
a
)
2
(
y
b
)
2
(
1
-
x
a
)
2
(
1
-
y
b
)
2
]
×
[
A
1
+
A
2
(
x
a
)
(
y
b
)
(
1
-
x
a
)
(
1
-
y
b
)
]
S-S-S-S
W
=
[
(
x
a
)
(
y
b
)
(
1
-
x
a
)
(
1
-
y
b
)
]
×
[
A
1
+
A
2
(
x
a
)
(
y
b
)
(
1
-
x
a
)
(
1
-
y
b
)
]
C-S-C-S
W
=
[
(
x
a
)
2
(
y
b
)
(
1
-
x
a
)
2
(
1
-
y
b
)
]
×
[
A
1
+
A
2
(
x
a
)
(
y
b
)
(
1
-
x
a
)
(
1
-
y
b
)
]


## 2.3. Solution of Frequency Equation

Rayleigh Ritz technique is applied to solve the frequency equation. In this method, one requires that the maximum strain energy(
S
E
) must be equal to the maximum kinetic energy (K
E). So it is necessary for the problem under consideration that [16]
(7)
δ
(
S
E
-
K
E
)
=
0
.
Now assume the nondimensional variables as
(8)
X
=
x
a
,
Y
=
y
a
,
W
-
=
W
a
,
d
-
=
d
a
.
The expressions for kinetic energy (K
E) and strain energy (S
E) are
(9)
K
E
=
1
2
ρ
p
2
d
0
-
a
5
×
∫
0
1
∫
0
b
/
a
[
(
1
+
β
1
X
2
)
(
1
+
β
2
a
2
b
2
Y
2
)
W
-
2
]
d
Y
d
X
,
S
E
=
E
0
d
0
-
3
a
3
24
(
1
-
υ
2
)
×
∫
0
1
∫
0
b
/
a
{
1
-
α
(
1
-
X
2
)
(
1
-
a
2
b
2
Y
2
)
}
×
(
1
+
β
1
X
2
)
3
(
1
+
β
2
a
2
b
2
Y
2
)
3
×
[
(
∂
2
W
-
∂
X
2
)
2
+
(
∂
2
W
-
∂
Y
2
)
2
+
2
υ
∂
2
W
-
∂
X
2
∂
2
W
-
∂
Y
2
+
2
(
1
-
υ
)
(
∂
2
W
-
∂
X
∂
Y
)
2
]
d
Y
d
X
On using (9) in (7), one gets
(10)
S
E
*
-
Ω
2
K
E
*
=
0
,
where
(11)
K
E
*
=
∫
0
1
∫
0
b
/
a
[
(
1
+
β
1
X
2
)
(
1
+
β
2
a
2
b
2
Y
2
)
W
-
2
]
d
Y
d
X
,
S
E
*
=
∫
0
1
∫
0
b
/
a
{
1
-
α
(
1
-
X
2
)
(
1
-
a
2
b
2
Y
2
)
}
×
(
1
+
β
1
X
2
)
3
(
1
+
β
2
a
2
b
2
Y
2
)
3
×
[
(
∂
2
W
-
∂
X
2
)
2
+
(
∂
2
W
-
∂
Y
2
)
2
+
2
υ
∂
2
W
-
∂
X
2
∂
2
W
-
∂
Y
2
+
2
(
1
-
υ
)
(
∂
2
W
-
∂
X
∂
Y
)
2
]
d
Y
d
X
.
Here, Ω
2
=
(
12
ρ
p
2
a
2
(
1
-
υ
2
)
)
/
(
E
0
d
0
-
2
) is a frequency parameter.Equation (10) consists of two unknown constants, that is, A
1 and A
2, arising due to the substitution of W. These two constants are to be determined as follows:
(12)
∂
∂
A
n
(
S
E
*
-
Ω
2
K
E
*
)
=
0
,
n
=
1
,
2
.
On simplifying (12), one gets
(13)
b
n
1
A
1
+
b
n
2
A
2
=
0
,
n
=
1
,
2
,
where b
n
1 and b
n
2 include parametric constant, that is, taper constants, thermal gradient, aspect ratio, and frequency parameter.For a nontrivial solution, determinant of the coefficient must be zero. So, one gets the frequency equation as(14)
|
b
11
b
12
b
21
b
22
|
=
0
.
With the help of (14), one can obtain a quadratic equation in Ω
2 from which the two values of frequency parameter for both the modes of vibration can be evaluated easily.


## 3. Results and Discussion

Computations have been made for calculating frequencies for different values of thermal gradient (α), taper constants (β
1 and β
2), and aspect ratio (a
/
b) for the first two modes of vibration. In calculation, Poisson ratio is considered as constant, that is, 0.345, and thickness of the plate at X
=
Y
=
0 is taken as 0.01 m. In Tables 2, 3, and 4, frequency for the first two modes of vibration is reported for various values of taper constant (β
1 and β
2) and thermal gradient (α) at fixed aspect ratio (
a
/
b
=
1.5
). Also, frequency is tabulated for various values of aspect ratio in Table 5 for the following cases: (i) α = β
1 = β
2 = 0.0, (ii) α = β
1= β
2 = 0.2, (iii) α = β
1 = β
2 = 0.4, (iv) α = β
1 = β
2 = 0.6, and (v) α = β
1 = β
2= 0.8.Table 2
Frequency versus taper constants at fixeda
/
b
=
1.5
,
α
=
0.0.
Boundary conditions
β
2
β
1
0.0
0.2
0.4
0.6
0.8
Mode I
Mode II
Mode I
Mode II
Mode I
Mode II
Mode I
Mode II
Mode I
Mode II
C-C-C-C
0.0
240.47
60.85
271.43
68.92
287.81
73.31
305.20
77.97
323.46
82.85
S-S-S-S
256.91
34.25
273.20
36.62
291.08
39.22
310.34
42.01
330.77
44.97
C-S-C-S
288.82
41.94
305.18
44.90
322.47
48.16
340.57
51.66
359.37
55.39
C-C-C-C
0.2
277.75
70.26
294.19
74.65
311.85
79.36
330.58
84.35
350.23
89.56
S-S-S-S
281.89
36.63
299.57
39.17
318.95
41.97
339.79
44.98
361.87
48.17
C-S-C-S
317.84
44.79
335.78
47.94
354.71
51.41
374.50
55.15
395.06
59.13
C-C-C-C
0.4
302.35
76.32
320.15
81.04
339.26
86.10
359.50
91.44
380.72
97.02
S-S-S-S
310.85
39.22
330.13
41.97
351.22
44.99
373.86
48.24
397.83
51.69
C-S-C-S
351.62
47.88
371.38
51.25
392.22
54.96
413.99
58.96
436.57
63.20
C-C-C-C
0.6
329.48
82.85
348.78
87.94
369.47
93.37
391.37
99.09
414.32
105.07
S-S-S-S
342.97
42.01
364.04
44.97
387.02
48.24
411.65
51.76
437.70
55.49
C-S-C-S
389.13
51.19
410.93
54.80
433.88
58.76
457.84
63.04
482.67
67.57
C-C-C-C
0.8
358.70
89.75
379.61
95.21
402.00
101.04
425.68
107.17
450.47
113.58
S-S-S-S
377.61
44.96
400.59
48.15
425.61
51.68
452.40
55.48
480.69
59.51
C-S-C-S
429.56
54.68
453.54
58.53
478.77
62.77
505.09
67.33
532.36
72.18Table 3
Frequency versus taper constants at fixeda
/
b
=
1.5
,
α
=
0.2.
Boundary conditions
β
2
β
1
0
0.2
0.4
0.6
0.8
Mode I
Mode II
Mode I
Mode II
Mode I
Mode II
Mode I
Mode II
Mode I
Mode II
C-C-C-C
0
244.24
61.78
259.23
65.80
275.35
70.12
292.47
74.69
310.43
79.47
S-S-S-S
246.02
32.58
262.23
34.93
280.01
37.51
299.14
40.27
319.4
43.19
C-S-C-S
276.34
39.90
292.34
42.84
309.24
46.07
326.93
49.55
345.32
53.24
C-C-C-C
0.2
266.18
67.25
282.38
71.57
299.79
76.20
318.24
81.08
337.61
86.19
S-S-S-S
271.74
34.94
289.35
37.47
308.63
40.24
329.34
43.22
351.27
46.37
C-S-C-S
306.35
42.71
323.95
45.85
342.53
49.30
361.95
53.01
382.12
56.94
C-C-C-C
0.4
291.10
73.32
308.67
77.96
327.53
82.92
347.50
88.16
368.44
93.63
S-S-S-S
301.30
37.51
320.52
40.24
341.51
43.23
364.03
46.46
387.85
49.87
C-S-C-S
340.94
45.77
360.39
49.13
380.89
52.81
402.31
56.78
424.53
60.98
C-C-C-C
0.6
318.49
79.81
337.56
84.81
358.00
90.14
379.63
95.76
402.30
101.62
S-S-S-S
333.89
40.27
354.89
43.21
377.77
46.45
402.29
49.94
428.19
53.63
C-S-C-S
379.07
49.04
400.57
52.63
423.20
56.56
446.83
60.80
471.322
65.30
C-C-C-C
0.8
347.88
86.65
368.57
92.02
390.71
97.74
414.12
103.77
438.64
110.04
S-S-S-S
368.87
43.18
391.78
46.35
416.71
49.85
443.36
53.62
471.51
57.62
C-S-C-S
419.95
52.47
443.65
56.30
468.57
60.51
494.57
65.05
521.51
69.85Table 4
Frequency versus taper constants at fixeda
/
b
=
1.5
,
α
=
0.6.
Boundary conditions
β
2
β
1
0
0.2
0.4
0.6
0.8
Mode I
Mode II
Mode I
Mode II
Mode I
Mode II
Mode I
Mode II
Mode I
Mode II
C-C-C-C
0
218.37
55.16
232.93
59.04
248.58
63.19
265.19
67.57
282.61
72.14
S-S-S-S
222.63
28.95
238.78
31.28
256.45
33.82
275.40
36.53
295.43
39.39
C-S-C-S
249.51
35.45
264.78
38.37
280.92
41.58
297.81
45.01
315.35
48.64
C-C-C-C
0.2
241.39
60.73
257.16
64.89
274.10
69.33
292.05
74.02
310.88
78.90
S-S-S-S
250.19
31.28
267.73
33.78
286.88
36.52
307.39
39.45
329.06
42.53
C-S-C-S
281.98
38.22
298.91
41.34
316.77
44.75
335.44
48.42
354.82
52.28
C-C-C-C
0.4
267.22
66.80
284.36
71.27
302.74
76.04
322.21
81.06
342.61
86.29
S-S-S-S
281.23
33.81
300.37
36.51
321.21
39.47
343.51
42.64
367.06
45.98
C-S-C-S
318.52
41.21
337.33
44.54
357.16
48.19
377.88
52.11
399.36
56.24
C-C-C-C
0.6
295.31
73.22
313.95
78.03
333.92
83.16
355.05
88.55
377.17
94.16
S-S-S-S
314.95
36.52
335.85
39.43
358.57
42.63
382.86
46.06
408.48
49.69
C-S-C-S
358.12
44.40
379.01
47.96
401.01
51.86
423.96
56.04
447.75
60.47
C-C-C-C
0.8
325.22
79.92
345.47
85.10
367.14
90.60
390.04
96.38
414.01
102.40
S-S-S-S
350.74
39.36
373.54
42.51
398.30
45.97
424.72
49.68
452.58
53.62
C-S-C-S
400.05
47.74
423.17
51.54
447.48
55.71
472.83
60.18
499.09
64.92Table 5
Frequency versus aspect ratio.
Boundary conditions
a
/
b
α
=
β
1
=
β
2
=
0.0
α
=
β
1
=
β
2
=
0.2
α
=
β
1
=
β
2
=
0.4
α
=
β
1
=
β
2
=
0.6
α
=
β
1
=
β
2
=
0.8
Mode I
Mode II
Mode I
Mode II
Mode I
Mode II
Mode I
Mode II
Mode I
Mode II
C-C-C-C
0.5
104.94
26.24
115.80
28.96
129.50
32.20
145.99
35.82
165.16
39.72
S-S-S-S
105.91
13.17
119.29
14.32
136.67
15.69
157.82
17.29
182.48
19.13
C-S-C-S
86.91
25.38
96.36
28.13
108.49
31.44
123.29
35.21
140.66
39.34
C-C-C-C
1
150.10
38.35
165.27
42.35
184.37
47.10
207.31
52.41
233.97
58.14
S-S-S-S
149.67
21.08
168.55
23.15
193.11
25.72
223.06
28.80
258.06
32.39
C-S-C-S
151.36
30.95
169.11
34.14
191.95
38.04
219.64
42.58
251.88
47.74
C-C-C-C
1.5
256.20
60.85
282.38
71.57
315.37
79.58
355.05
88.55
401.16
98.21
S-S-S-S
256.91
34.25
289.35
37.47
331.51
41.40
382.86
46.06
442.82
51.49
C-S-C-S
288.82
41.94
323.95
45.85
369.22
50.56
423.96
56.04
487.50
62.28
C-C-C-C
2
419.78
104.96
463.20
115.86
518.02
128.80
583.98
143.29
660.66
158.90
S-S-S-S
423.64
52.71
477.17
57.29
546.71
62.78
631.29
69.19
729.93
76.55
C-S-C-S
491.36
58.84
551.87
63.78
629.85
69.58
724.11
76.22
833.37
83.67
C-C-C-C
2.5
636.02
157.78
702.29
174.16
786.05
193.58
886.90
215.34
1004.13
238.79
S-S-S-S
644.38
76.44
725.84
82.66
831.61
89.98
960.19
98.38
1110.02
107.88
C-S-C-S
755.16
81.55
848.63
87.86
969.10
95.07
1114.67
103.146
1283.32
112.03In Tables2–4, the authors found that frequency for both the modes of vibration continuously increases as taper constants (β
1 and β
2) increase from 0.0 to 0.8 for increasing values of thermal gradient (α).It is interesting to notice that the first mode of frequency is maximum at C-S-C-S boundary condition while it is minimum at C-C-C-C boundary condition. But for the second mode of frequency, it is vice versa; that is, the second mode is maximum at C-C-C-C boundary condition while it is minimum at S-S-S-S boundary condition.Here, the authors also noticed that as thermal gradient (α) increases from 0.0 (Table 2) to 0.6 (Table 4), frequency for both the modes decreases at all values of taper parameters for these boundary conditions.From Table5, the authors found that frequency increases for both the modes of vibration for these boundary conditions as aspect ratio increases from 0.5 to 2.5. Frequency also increases as combined values of α, β
1 and β
2 increase from 0.0 to 0.8.


## 4. Comparison

In Tables6 and 7, frequency for both the modes of the present paper is compared with [17] for various values of taper constants, thermal gradient, and aspect ratio.Table 6
Frequency versus thermal gradient*.
α
β
1
=
0.0
,
β
2
=
0.0
β
1
=
0.2
,
β
2
=
0.6
Mode I
Mode II
Mode I
Mode II
0.0
240.47
60.84
340.87
85.87
240.47
64.84
348.78
87.94
0.2
228.13
57.72
322.11
81.42
244.24
61.78
337.56
84.81
0.4
215.08
54.42
305.99
76.88
231.67
58.57
325.96
81.52
0.6
201.97
51.06
284.00
71.90
218.37
55.16
313.95
78.03
0.8
186.26
47.13
265.05
66.19
204.22
51.50
301.74
74.30
*Results of the present paper are shown in bold text.Table 7
Frequency versus taper constant (β
1)*.
β
1
α
=
0.0
,
β
2
=
0.0
Mode I
Mode II
0.0
240.47
60.84
240.47
60.84
0.2
261.60
66.34
271.43
68.92
0.4
282.61
72.59
287.81
73.31
0.6
301.10
78.00
305.20
77.97
0.8
319.89
84.00
323.46
82.85
*Results of the present paper are shown in bold text.Table6 shows values of frequency for the first two modes of vibration with increasing values of thermal gradient for fixed values of taper constants. Here, the authors noticed that both the modes of frequency in the present paper are greater than [17] at corresponding values of thermal gradient and taper parameters.Thickness variation in present paper is taken as biparabolic while, in [17], it is linear in x-direction and parabolic in y-direction. Therefore, the authors compared both the modes of frequency of present paper with [17] for increasing values of β
1 only at α
=
0.0, β
2
=
0.0 as shown in Table 7.Here, the authors found that the first mode of frequency is greater than [17] and it is continuously increasing with increasing values of β
1 but the second mode of frequency is greater than [17] for β
1
=
0.0
,
0.2
,
0.4, and at β
1
=
0.6
,
0.8, it is smaller than [17].


## 5. Conclusions

On behalf of the above comparisons, the authors conclude the following.(i)
Frequency in the present paper (biparabolic temperature variation) is greater than [17] (linear temperature variation).
(ii)
The first mode of frequency in the present paper (biparabolic thickness variation) is greater than [17] (bidimensional thickness variation, that is, linear in x-direction and parabolic in y-direction).
(iii)
The main emphasis of the present study is to provide some numerical data about the first few modes of vibration for design engineers and researchers to fulfill their requirements. The authors suggest that the readers and engineers can obtain the desired values of frequency by appropriate tapering of the plate and can provide more realistic and reliable structures and designs.



---

*Source: 642926-2014-07-08.xml*