Natural vibration of tapered rectangular plate with exponential variation in non homogeneity

Sharma, Amit; Mani, Naveen; Bhardwaj, Reeta

doi:10.21595/jve.2018.19589

Journal of Vibroengineering

Browse Journal

Submit article

Published: 15 February 2019

Check for updates

Natural vibration of tapered rectangular plate with exponential variation in non homogeneity

Amit Sharma¹

Naveen Mani²

Reeta Bhardwaj³

^{1, 3}Amity University Haryana, Gurgaon, India

²Sandip University, Nasik, Maharashtra, India

Corresponding Author:

Amit Sharma

Cite the article Download PDF

Downloads 1497

WoS Core Citations 6

CrossRef Citations 3

Abstract

In this paper, authors studied the natural vibration of tapered non homogeneous rectangular plate on clamped edges. For tapering in plate, authors considered circular variation in thickness and for non-homogeneity (in plate’s material) Poisson’s ratio varies exponentially. Bilinear temperature (linear along both the axes) variation on the plate is being viewed. Rayleigh Ritz method is used to solve differential equation of motion. All the results are presented with the help of tables and graphs. A comparison of results is also given to support the present study.

1. Introduction

The non uniform and non homogeneous rectangular plate are used in various engineering structures such as aerospace engineering, marine engineering, ocean engineering, optical instruments and mechanical engineering. Non uniform and non homogeneous plates have greater efficiency when compared to uniform and homogeneous plates because of reduction in weight, high tensile strength, durability, elastic behavior and size. The study of plate vibration without consideration of temperature is nothing because almost all engineering structure worked under great influence of temperature. Thus, to design an accurate structure, determination of natural frequencies and mode shape are essential. A significant work has been provided in these directions.

Timoshenko and Woinowsky-Krieger [1] discussed comprehensive background on plate and shells theory in their book and explained bending (pure and symmetrical) of long, circular, anisotropic and plates having lateral loads, forces in the middle. They also studied rectangular plates (on simply supported and on various edge conditions), plates on elastic foundation and plates of various shapes. Bhat [2] provided natural frequencies of rectangular plates on clamped, simply supported and combination of simply supported and free edges by applying beam characteristic orthogonal polynomials in the Rayleigh Ritz method and showed that when some edges of plate are free, the present method acquires superior result for lower modes. Chakraverty and Petyt [3] discussed transverse vibration of elliptic and circular plates on clamped, simply supported and free boundary using two-dimensional boundary characteristic orthogonal polynomials in Rayleigh Ritz method and evaluated first five natural frequencies for various combinations of non homogeneity and aspect ratios of the ellipse. Li [4] applied Rayleigh Ritz method to analyze the modal characteristics of a rectangular plate with general elastic supports alone its edges and provided numerical examples to demonstrate the accuracy and convergence of the current solution. Ranji and Shahbaztab [5] analyzed free vibration of non homogeneous and orthotropic plate resting on a Pasternak type of elastic foundation by using Rayleigh Ritz method and provided the effect of foundation, density and non homogeneity parameters on natural frequency. Singh and Sexena [6] applied Rayleigh Ritz method to study the transverse vibrations of tapered skew plates with different combinations of boundary conditions and obtained first three modes of frequency and mode shapes. Sharma et al. [7, 8] studied the natural vibration of tapered (exponential and circular variation in thickness) square plate with non homogeneity (variation in density and Poisson’s ratio) under temperature variation (bi linear and bi parabolic). An effect of circular variation in thickness on vibrational frequencies of parallelogram plate has been studied by Sharma [9] and obtained first two modes of vibrations for different combinations of plate parameters on clamped edges. Sharma et al. [10] also studied the vibration of tapered rectangular plate with linear variation in density and evaluated first two modes of vibration on clamped edges. Khanna and Kaur [11, 12] analyzed the free vibration of non uniform (linear and exponential variation in thickness) and non homogeneous (exponential variation in Poisson’s ratio) of rectangular plate with temperature (linear and exponential) variation by using Rayleigh Ritz method on clamped boundary condition. Hwu [13] presented closed-form solution for the free vibration problems, to study the effects of rotary inertia and shear deformation on frequency and provided numerical examples which shows the importance of considering the effects of rotary inertia and shear deformation. Avalos and Laura [14] studied the transverse vibrations of isotropic, orthotropic and anisotropic rectangular plates with two rectangular cutouts on simply supported edge using Classical plate theory. Shufrin and Eisenberger [15] presented a new semi-analytical method for modeling rectangular plates (with one dimensional variable thickness and cutouts) to analyze the free vibration of rectangular plates with various rectangular cutouts and variable thickness. In order to show the accuracy and convergence of the solution, they compared their results with other semi-analytical methods. Sakiyama and Huang [16] presented an approximate method to analyze the free vibration of thin and moderately thick rectangular plates with arbitrary variable thickness using Green function and showed that the numerical solution for the Green function has good convergence and accuracy. Wang and Zu [17] studied the vibrations of functionally graded material (FGM) rectangular plates with porosities and moving in thermal environment by using von Kármán nonlinear plate theory and obtained vibration characteristics such as natural frequency and nonlinear frequency response. Huang et al. [18] developed a discrete method to analyze the free vibration problem of orthotropic rectangular plates with variable thickness using Green function and obtained the effects of the aspect ratios, boundary conditions and the variation of the thickness on the frequencies. Huang et al. [19] proposed new method (FBGM) to analyze the free vibration of cantilever plates with variable thickness and obtained the solution of characteristic equation of free vibration. Taher et al. [20] evaluated first nine frequency parameters of circular and annular plates with variable thickness and combined boundary conditions, for different thickness to radius ratios by using three-dimensional elasticity theory. Zur [21] discussed the free vibration of homogeneous and isotropic annular thin plates with variable distributions of parameters by using Green's function and Neumann series, and obtained numerical solutions of the characteristic equations for constant and hyperbolic varying thickness on different boundary conditions. Zur [22] discussed the free vibration analysis of homogeneous and isotropic circular thin plates with variable distribution of parameters by using Green's functions and evaluated six lower natural dimensionless frequencies of axisymmetric vibration of circular plates of constant and variable thickness.

In this study, authors show the effect of circular variation in thickness on free vibration of non homogeneous rectangular plate with linear (along both the axes) variation in temperature and compute the vibrational frequency modes (first two modes). The results are presented with the help of tables and figures.

2. Differential equation of motion

The differential equation for transverse motion of the plate is:

1

\frac{\partial^{2} M_{ζ}}{\partial ζ^{2}} + 2 \frac{\partial^{2} M_{ζ ψ}}{\partial ζ \partial ψ} + \frac{\partial^{2} M_{ψ}}{\partial ψ^{2}} = ρ l \frac{\partial^{2} ϕ}{\partial t^{2}},

where:

2

M_{ζ} = - D_{1} [\frac{\partial^{2} ϕ}{\partial ζ^{2}} + ν \frac{\partial^{2} ϕ}{\partial ψ^{2}}], M_{ψ} = - D_{1} [\frac{\partial^{2} ϕ}{\partial ψ^{2}} + ν \frac{\partial^{2} ϕ}{\partial ζ^{2}}], M_{ζ ψ} = - D_{1} (1 - ν) \frac{\partial^{2} ϕ}{\partial ζ \partial ψ} .

Using Eq. (2), Eq. (1) becomes:

3

[D_{1} (\frac{\partial^{4} ϕ}{\partial ζ^{4}} + 2 \frac{\partial^{4} ϕ}{\partial ζ^{2} \partial ψ^{2}} + \frac{\partial^{4} ϕ}{\partial ψ^{4}} + \frac{\partial^{2} ν}{\partial ζ^{2}} \frac{\partial^{2} ϕ}{\partial ψ^{2}}) + 2 \frac{\partial D_{1}}{\partial ζ} (\frac{\partial^{3} ϕ}{\partial ζ^{3}} + \frac{\partial^{3} ϕ}{\partial ζ \partial ψ^{2}} + \frac{\partial ν}{\partial ζ} \frac{\partial^{2} ϕ}{\partial ψ^{2}})

+ 2 \frac{\partial D_{1}}{\partial ψ} (\frac{\partial^{3} ϕ}{\partial ψ^{3}} + \frac{\partial^{3} ϕ}{\partial ψ \partial ζ^{2}} - \frac{\partial ν}{\partial ζ} \frac{\partial^{2} ϕ}{\partial ζ \partial ψ}) + \frac{\partial^{2} D_{1}}{\partial ζ^{2}} (\frac{\partial^{2} ϕ}{\partial ζ^{2}} + ν \frac{\partial^{2} ϕ}{\partial ψ^{2}})

+ \frac{\partial^{2} D_{1}}{\partial ψ^{2}} (\frac{\partial^{2} ϕ}{\partial ψ^{2}} + ν \frac{\partial^{2} ϕ}{\partial ζ^{2}}) + 2 (1 - ν) \frac{\partial^{2} D_{1}}{\partial ζ \partial ψ} \frac{\partial^{2} ϕ}{\partial ζ \partial ψ}] + ρ l \frac{\partial^{2} ϕ}{\partial t^{2}} = 0 .

For solution of Eq. (3), we can take deflection function as:

4

ϕ (ζ, ψ, t) = Φ (ζ, ψ) T (t) .

Using Eq. (4) in Eq. (3), we have:

5

T [D_{1} (\frac{\partial^{4} Φ}{\partial ζ^{4}} + 2 \frac{\partial^{4} Φ}{\partial ζ^{2} \partial ψ^{2}} + \frac{\partial^{4} Φ}{\partial ψ^{4}} + \frac{\partial^{2} ν}{\partial ζ^{2}} \frac{\partial^{2} Φ}{\partial ψ^{2}}) + 2 \frac{\partial D_{1}}{\partial ζ} (\frac{\partial^{3} Φ}{\partial ζ^{3}} + \frac{\partial^{3} Φ}{\partial ζ \partial ψ^{2}} + \frac{\partial ν}{\partial ζ} \frac{\partial^{2} Φ}{\partial ψ^{2}})

+ 2 \frac{\partial D_{1}}{\partial ψ} (\frac{\partial^{3} Φ}{\partial ψ^{3}} + \frac{\partial^{3} Φ}{\partial ψ \partial ζ^{2}} - \frac{\partial ν}{\partial ζ} \frac{\partial^{2} Φ}{\partial ζ \partial ψ}) + \frac{\partial^{2} D_{1}}{\partial ζ^{2}} (\frac{\partial^{2} Φ}{\partial ζ^{2}} + ν \frac{\partial^{2} Φ}{\partial ψ^{2}})

+ \frac{\partial^{2} D_{1}}{\partial ψ^{2}} (\frac{\partial^{2} Φ}{\partial ψ^{2}} + ν \frac{\partial^{2} Φ}{\partial ζ^{2}}) + 2 (1 - ν) \frac{\partial^{2} D_{1}}{\partial ζ \partial ψ} \frac{\partial^{2} Φ}{\partial ζ \partial ψ}] + ρ l Φ \frac{\partial^{2} T}{\partial t^{2}} = 0 .

Now by using variable separable technique Eq. (5) becomes:

6

\frac{[\begin{matrix} D_{1} (\frac{\partial^{4} Φ}{\partial ζ^{4}} + 2 \frac{\partial^{4} Φ}{\partial ζ^{2} \partial ψ^{2}} + \frac{\partial^{4} Φ}{\partial ψ^{4}} + \frac{\partial^{2} ν}{\partial ζ^{2}} \frac{\partial^{2} Φ}{\partial ψ^{2}}) + 2 \frac{\partial D_{1}}{\partial ζ} (\frac{\partial^{3} Φ}{\partial ζ^{3}} + \frac{\partial^{3} Φ}{\partial ζ \partial ψ^{2}} + \frac{\partial ν}{\partial ζ} \frac{\partial^{2} Φ}{\partial ψ^{2}}) \\ + 2 \frac{\partial D_{1}}{\partial ψ} (\frac{\partial^{3} Φ}{\partial ψ^{3}} + \frac{\partial^{3} Φ}{\partial ψ \partial ζ^{2}} - \frac{\partial ν}{\partial ζ} \frac{\partial^{2} Φ}{\partial ζ \partial ψ}) + \frac{\partial^{2} D_{1}}{\partial ζ^{2}} (\frac{\partial^{2} Φ}{\partial ζ^{2}} + ν \frac{\partial^{2} Φ}{\partial ψ^{2}}) \\ + \frac{\partial^{2} D_{1}}{\partial ψ^{2}} (\frac{\partial^{2} Φ}{\partial ψ^{2}} + ν \frac{\partial^{2} Φ}{\partial ζ^{2}}) + 2 (1 - ν) \frac{\partial^{2} D_{1}}{\partial ζ \partial ψ} \frac{\partial^{2} Φ}{\partial ζ \partial ψ} \end{matrix}]}{ρ l Φ} = - \frac{1}{T} \frac{\partial^{2} T}{\partial t^{2}} = ω^{2} .

Taking first and last expression of Eq. (6), we have:

7

[D_{1} (\frac{\partial^{4} Φ}{\partial ζ^{4}} + 2 \frac{\partial^{4} Φ}{\partial ζ^{2} \partial ψ^{2}} + \frac{\partial^{4} Φ}{\partial ψ^{4}} + \frac{\partial^{2} ν}{\partial ζ^{2}} \frac{\partial^{2} Φ}{\partial ψ^{2}}) + 2 \frac{\partial D_{1}}{\partial ζ} (\frac{\partial^{3} Φ}{\partial ζ^{3}} + \frac{\partial^{3} Φ}{\partial ζ \partial ψ^{2}} + \frac{\partial ν}{\partial ζ} \frac{\partial^{2} Φ}{\partial ψ^{2}})

+ 2 \frac{\partial D_{1}}{\partial ψ} (\frac{\partial^{3} Φ}{\partial ψ^{3}} + \frac{\partial^{3} Φ}{\partial ψ \partial ζ^{2}} - \frac{\partial ν}{\partial ζ} \frac{\partial^{2} Φ}{\partial ζ \partial ψ}) + \frac{\partial^{2} D_{1}}{\partial ζ^{2}} (\frac{\partial^{2} Φ}{\partial ζ^{2}} + ν \frac{\partial^{2} Φ}{\partial ψ^{2}})

+ \frac{\partial^{2} D_{1}}{\partial ψ^{2}} (\frac{\partial^{2} Φ}{\partial ψ^{2}} + ν \frac{\partial^{2} Φ}{\partial ζ^{2}}) + 2 (1 - ν) \frac{\partial^{2} D_{1}}{\partial ζ \partial ψ} \frac{\partial^{2} Φ}{\partial ζ \partial ψ}] - ρ l Φ ω^{2} = 0 .

Eq. (7) is differential equation for transverse motion of plate with variable $D_{1}$ and $ν$ i.e., flexural rigidity and Poisson’s ratio. The expression for flexural rigidity is given by $D_{1} = E l^{3} / 12 (1 - ν^{2})$ .

3. Construction of problem

Consider a non homogeneous rectangular plate of length $a$ and breadth $b$ with one dimensional circular variation in thickness $l$ as shown in Fig. 1 as:

8

l = l_{0} [1 + β (1 - \sqrt{1 - \frac{ζ^{2}}{a^{2}}})],

where $β$ , $(0 \leq β \leq 1)$ is known as tapering parameter. Thickness of plate become constant at $ζ =$ 0. The plate is subjected to steady two dimensional linear temperature distributions as:

9

τ = τ_{0} (1 - \frac{ζ}{a}) (1 - \frac{ψ}{b}),

where $τ$ and $τ_{0}$ denotes the temperature excess above the reference temperature on the plate at any point and at the origin respectively. The temperature dependence modulus of elasticity for engineering structures is given by:

10

E = E_{0} (1 - γ τ),

where $E_{0}$ is the Young’s modulus at mentioned temperature (i.e., $τ =$ 0) and $γ$ is called slope of variation.

Fig. 1Rectangular plate with one dimensional circular variation

Using Eq. (9), Eq. (10) becomes:

11

E = E_{0} [1 - α \{1 - \frac{ζ}{a}\} \{1 - \frac{ψ}{b}\}],

where $α$ , (0 ≤ $α$ < 1) is called temperature gradient, which is the product of temperature at origin and slope of variation i.e., $α = γ τ_{0}$ . For non homogeneous consideration in plate’s material, Poisson’s ratio of the plate varies exponentially in one direction as:

12

ν = ν_{0} [e^{m \frac{ζ}{a}}],

where $m$ , (0 ≤ $m$ ≤ 1) are known as non homogeneity constant.

Using Eqs. (8), (11) and (12), flexural rigidity of the plate becomes:

13

D_{1} = E_{0} {l_{0}}^{3} [1 - α \{1 - \frac{ζ}{a}\} \{1 - \frac{ψ}{b}\}] {[1 + β (1 - \sqrt{1 - \frac{ζ^{2}}{a^{2}}})]}^{3} / 12 (1 - {ν_{0}}^{2} e^{2 m \frac{ζ}{a}}) .

4. Solution of the problem

We are using Rayleigh Ritz technique (i.e., maximum strain energy $V_{s}$ must equal to maximum kinetic energy $T_{s}$ ) in order to obtain frequency equation and frequencies for both modes of vibrations. Therefore, we must have:

14

δ (V_{s} - T_{s}) = 0 .

Here the expression for $V_{s}$ and $T_{s}$ are given by:

15

V_{s} = \frac{1}{2} \int_{0}^{a} ‍ \int_{0}^{b} ‍ D_{1} \times [{(\frac{\partial^{2} Φ}{\partial ζ^{2}})}^{2} + {(\frac{\partial^{2} Φ}{\partial ψ^{2}})}^{2} + 2 ν \frac{\partial^{2} Φ}{\partial ζ^{2}} \frac{\partial^{2} Φ}{\partial ψ^{2}} + 2 (1 - ν) {(\frac{\partial^{2} Φ}{\partial ζ \partial ψ})}^{2}] d ψ d ζ,

16

T_{s} = \frac{1}{2} ω^{2} ρ \int_{0}^{a} \int_{0}^{b} l Φ^{2} d ψ d ζ .

Here, we are computing frequency on C-C-C-C condition (i.e., all the four edges are clamped), therefore the boundary conditions are:

17

Φ (ζ, ψ) = \frac{\partial Φ (ζ, ψ)}{\partial ζ} = 0, ζ = 0, a,

Φ (ζ, ψ) = \frac{\partial Φ (ζ, ψ)}{\partial ψ} = 0, ψ = 0, b .

Therefore, deflection function (i.e, maximum displacement) which satisfy boundary condition given in Eq. (17) is taken as [11]:

18

Φ (ζ, ψ) = {(\frac{ζ}{a})}^{2} {(\frac{ψ}{b})}^{2} {(1 - \frac{ζ}{a})}^{2} {(1 - \frac{ψ}{b})}^{2} [Ω_{1} + Ω_{2} (\frac{ζ}{a}) (\frac{ψ}{b}) (1 - \frac{ζ}{a}) (1 - \frac{ψ}{b})],

where $Ω_{1}$ and $Ω_{2}$ are arbitrary constants. On substituting Eqs. (8), (12) and (13) in Eqs. (15) and (16) we have:

19

V_{s} = \frac{E_{0} {l_{0}}^{3}}{24} \int_{0}^{a} \int_{0}^{b} ‍ [\begin{array}{l} \{\frac{[1 - α \{1 - \frac{ζ}{a}\} \{1 - \frac{ψ}{b}\}] {[1 + β (1 - \sqrt{1 - \frac{ζ^{2}}{a^{2}}})]}^{3}}{1 - {ν_{0}}^{2} e^{2 m \frac{ζ}{a}}}\} \\ ∙ [{(\frac{\partial^{2} Φ}{\partial ζ^{2}})}^{2} + {(\frac{\partial^{2} Φ}{\partial ψ^{2}})}^{2} + 2 ν_{0} e^{m \frac{ζ}{a}} \frac{\partial^{2} Φ}{\partial ζ^{2}} \frac{\partial^{2} Φ}{\partial ψ^{2}} + 2 (1 - ν_{0} e^{m \frac{ζ}{a}}) {(\frac{\partial^{2} Φ}{\partial ζ \partial ψ})}^{2}] \end{array}] d ψ d ζ,

20

T_{s} = \frac{1}{2} ω^{2} ρ l_{0} \int_{0}^{a} \int_{0}^{b} [1 + β (1 - \sqrt{1 - \frac{ζ^{2}}{a^{2}}})] Φ^{2} d ψ d ζ .

Now converting $ζ$ and $ψ$ into non dimensional variable as:

21

ζ_{1} = \frac{ζ}{a}, ψ_{1} = \frac{ψ}{a} .

Using Eq. (21), Eqs. (19) and (20) converted into:

22

{V_{s}}^{*} = \frac{E_{0} {l_{0}}^{3}}{24 a^{2}} \int_{0}^{1} ‍ \int_{0}^{\frac{b}{a}} ‍ [\begin{matrix} \{\frac{[1 - α \{1 - ζ_{1}\} \{1 - \frac{a ψ_{1}}{b}\}] {[1 + β (1 - \sqrt{1 - {ζ_{1}}^{2}})]}^{3}}{1 - {ν_{0}}^{2} e^{2 m ζ_{1}}}\} \\ ∙ [\begin{matrix} {(\frac{\partial^{2} Φ}{\partial {ζ_{1}}^{2}})}^{2} + {(\frac{\partial^{2} Φ}{\partial {ψ_{1}}^{2}})}^{2} + 2 ν_{0} e^{m ζ_{1}} \frac{\partial^{2} Φ}{\partial {ζ_{1}}^{2}} \frac{\partial^{2} Φ}{\partial {ψ_{1}}^{2}} \\ + 2 (1 - ν_{0} e^{m ζ_{1}}) {(\frac{\partial^{2} Φ}{\partial ζ_{1} \partial ψ_{1}})}^{2} \end{matrix}] \end{matrix}] d ψ_{1} d ζ_{1},

23

{T_{s}}^{*} = \frac{1}{2} ω^{2} ρ a^{2} l_{0} \int_{0}^{1} \int_{0}^{\frac{b}{a}} [1 + β (1 - \sqrt{1 - {ζ_{1}}^{2}})] Φ^{2} d ψ_{1} d ζ_{1} .

Using Eqs. (22) and (23), Eq. (14) becomes:

24

δ ({V_{s}}^{*} - λ^{2} {T_{s}}^{*}) = 0,

where $λ^{2} = 12 ρ ω^{2} a^{4} / E_{0} {l_{0}}^{2}$ is known as frequency parameter. Eq. (24) consists of two unknown constants $Ω_{1}$ and $Ω_{2}$ (because of substitution of deflection function $Φ (ζ, ψ)$ ). These two unknowns could be calculated as follows:

25

\frac{\partial}{\partial Ω_{n}} ({V_{s}}^{*} - λ^{2} {T_{s}}^{*}) = 0, n = 1,2 .

After simplifying Eq. (25), we get system of homogeneous equations as:

26

\begin{array}{l} c_{11} Ω_{1} + c_{12} Ω_{2} = 0, \\ c_{21} Ω_{2} + c_{22} Ω_{2} = 0 . \end{array}

To obtain non zero solution (frequency equation), the determinant of coffiecient matrix (symmetric matrix) of Eq. (26) must zero i.e.:

27

|\begin{array}{l} c_{11} & c_{12} \\ c_{21} & c_{22} \end{array}| = 0 .

Eq. (27) is quadratic equation from which we get two modes as $λ_{1}$ (first mode) and $λ_{2}$ (second mode).

5. Results and discussion

The first two modes of vibrations are calculated and presented tabularly as well as graphically for different values of thermal gradient $α$ , non homogeneity constant $m$ and tapering parameters $β$ . The parameters used in numerical calculation are $E_{0} =$ 7.08×10¹⁰ nm^-2, $ν_{0} =$ 0.345, $l_{0} =$ 0.01 m, $ρ =$ 2.80×10³ kgm^-3.

Table 1 provides the natural vibration (first two modes) of plate corresponding to thermal gradient $α$ for three different combination of non homogeneity $m$ and taper constant $β$ i.e., $m = β =$ 0.2, $m = β =$ 0.4 and $m = β =$ 0.8. From Table 1, we can see that frequency of both modes decreases with the increasing value of thermal gradient $α$ for all the combination of non homogeneity $m$ and taper constant $β$ . While the frequency modes of vibration increases with the increasing value of thermal gradient $α$ when the combined value of non homogeneity $m$ and taper constant $β$ varies from 0.2 to 0.8.

Table 2 shows the frequency modes of vibration corresponding to non homogeneity constant $m$ for the following values of taper constant $β$ and thermal gradient $α$ . $α = β =$ 0.2, $α = β =$ 0.4 and $α = β =$ 0.8.

Table 2 provides the fact that vibrational frequency increases rapidly (not in linear fashion) when non homogeneity $m$ in plate’s material increases 0 to 1 for all the above mentioned values of thermal gradient $α$ and tapering parameter $β$ . The vibrational frequency also increases when the combined value of thermal gradient $α$ and tapering parameter $β$ varies from 0.2 to 0.8 with the increasing value of non homogeneity constant $m$ .

Table 1Thermal gradient α vs vibrational frequency λ for a/b= 1.5

$α$	$β = m =$ 0.2		$β = m =$ 0.4		$β = m =$ 0.8
$α$	$λ_{1}$	$λ_{2}$	$λ_{1}$	$λ_{2}$	$λ_{1}$	$λ_{2}$
0.0	68.16	268.66	72.15	283.85	82.92	327.10
0.2	66.50	262.09	70.45	277.19	81.14	320.21
0.4	64.79	255.36	68.72	270.37	79.31	313.16
0.6	63.04	248.45	66.94	263.37	77.44	305.96
0.8	61.24	241.45	65.11	256.18	75.52	298.58

Table 2Non homogeneity m vs vibrational frequency λ for a/b= 1.5

$m$	$α = β =$ 0.2		$α = β =$ 0.4		$α = β =$ 0.8
$m$	$λ_{1}$	$λ_{2}$	$λ_{1}$	$λ_{2}$	$λ_{1}$	$λ_{2}$
0.0	65.51	258.12	66.35	260.63	68.24	266.63
0.2	66.50	262.09	67.36	264.79	69.30	271.18
0.4	67.81	267.38	68.72	270.37	70.75	277.43
0.6	69.55	274.51	70.55	278.00	72.75	286.16
0.8	71.91	284.31	73.04	288.60	75.52	298.58
1.0	75.10	297.96	76.45	303.56	79.34	316.50

Table 3 gives the frequency modes corresponding to thickness variation (tapering parameter $β$ ) in plate for the following values of non homogeneity constant $m$ and thermal gradient $α$ . $m = α =$ 0.2, $m = α =$ 0.4 and $m = α =$ 0.8.

Form Table 3, one can easily concludes that when tapering parameter $β$ increases from 0 to 1, frequency modes increases for all the above mentioned values of non homogeneity $m$ and thermal gradient $α$ . The frequency modes decrease when the combined value of non homogeneity $m$ and thermal gradient $α$ varies from 0.2 to 0.4. But when the combined value of non homogeneity $m$ and thermal gradient $α$ varies from 0.4 to 0.8, the frequency modes again increases.

In order to get good understanding of results and discussion (variation of plate parameter), graphical representation of Tables 1-3 are presented in the form of Figs. 2-4.

Table 3Taper constant β vs vibrational frequency λ for a/b= 1.5

$β$	$m = α =$ 0.2		$m = α =$ 0.4		$m = α =$ 0.8
$β$	$λ_{1}$	$λ_{2}$	$λ_{1}$	$λ_{2}$	$λ_{1}$	$λ_{2}$
0.0	64.07	253.24	63.58	251.36	63.61	252.08
0.2	66.50	262.09	66.07	260.57	66.38	262.69
0.4	69.08	271.53	68.72	270.37	69.30	274.01
0.6	71.78	281.50	71.49	280.73	72.35	285.99
0.8	74.61	291.98	74.39	291.63	75.52	298.58
1.0	77.55	302.94	77.39	303.01	78.79	311.73

6. Results comparison

A comparison of present analysis with [11, 12] are presented in tabular (Table 4 and Table 5) as well as in graphical form (Fig. 5 and Fig. 6).

Fig. 2Thermal gradient (α) vs. frequency (λ) for fixed a/b= 1.5

a)

b)

Fig. 3Non homogeneity constant (m) vs. frequency (λ) for fixed a/b= 1.5

a)

b)

Fig. 4Taper constant (β) vs. frequency (λ) for fixed a/b= 1.5

a)

b)

Table 4 shows the comparison of frequency modes of present paper with [11, 12] corresponding to tapering parameter $β$ for fixed value of thermal gradient $α$ and non homogeneity constant $m$ i.e., $α = m =$ 0. Table 4 provides the information that frequency modes in present paper are less when compared with [11, 12]. The frequency modes of present paper, [11, 12] are coincides at $α = m = β =$ 0.

Fig. 5Comparison of frequency modes with [11, 12] corresponding to taper constant (β)

a)

b)

Fig. 6Comparison of frequency modes with [12] corresponding to thermal gradient (α)

Table 4Comparison of frequency modes of presented study with [11, 12] corresponding to taper constant β

$β$	$m = α =$ 0.0
$β$	$λ_{1}$	$λ_{2}$
0.0	64.77	255.98
	64.77	255.98
	64.77	255.98
0.2	71.40	282.12
	67.16	264.61
	71.84	283.76
0.4	78.27	309.17
	69.69	273.81
	80.25	316.46
0.6	85.34	336.90
	72.36	283.53
	90.21	354.87
Note: The value written in bold and italic are from [11, 12] respectively

Table 5 portrays the comparison of frequency modes of present paper with [12] corresponding to thermal gradient $α$ for fixed value of tapering parameter $β$ and non homogeneity constant $m$ i.e., $m = β =$ 0. It is noticed from Table 5 that frequency modes of present paper are less when compared with [12]. Frequency modes of present paper and [12] are coincides at $α = m = β =$ 0.

Table 5Comparison of frequency modes of present study with [12] corresponding to thermal gradient α

$α$	$β = m =$ 0.0
$α$	$λ_{1}$	$λ_{2}$
0.0	64.77	255.98
0.0	64.77	255.98
0.2	63.13	249.50
0.2	69.09	272.86
0.4	61.45	242.84
0.4	73.10	288.75
0.6	59.71	236.00
0.6	76.92	303.82
Note: The value written in bold are from [12]

7. Conclusions

From the results discussion and comparison authors would like to record the following points.

1) The frequency modes are less in case of circular variation in thickness (present paper) when compared to linear variation in thickness [11] and exponential variation in thickness [12] as shown in Table 4 and Fig. 5.

2) The frequency modes are less in case of bi linear temperature variation on plate (present paper) when compared to exponential variation in temperature on plate [12] as shown in Table 5 and Fig. 6.

3) The frequency modes decrease when the temperature increases on the plate as shown in Table 1 and Fig. 2.

4) When the non homogeneity increases in plate materials, frequency modes increases as shown in Table 2 and Fig. 3.

5) The frequency modes also increase when the thickness of plate increases as shown in Table 3 and Fig. 4.

References

Timoshenko S., Woinowsky Krieger S. Theory of plate and shells. McGraw-Hill, 1959.

Search CrossRef
Bhat R. B. Natural frequencies of rectangular plates using characteristic orthogonal polynomials in Rayleigh Ritz method. Journal of Sound and Vibration, Vol. 102, Issue 4, 1985, p. 493-499.

Publisher
Chakraverty S., Petyt M. Natural frequencies for free vibration of non homogeneous elliptic and circular plates using two-dimensional orthogonal polynomials. Applied Mathematical Modelling, Vol. 21, 1997, p. 399-417.

Publisher
Li W. L. Vibration analysis of rectangular plate with general elastic boundary support. Journal of Vibration and Control, Vol. 273, 2004, p. 619-635.

Publisher
Rahbar Ranji A., Shahbaztabar A. Free vibration analysis of non-homogeneous orthotropic plates resting on Pasternak elastic foundation by Rayleigh Ritz method. Journal of Central South University of Technology, Vol. 22, 2015, p. 413-420.

Publisher
Singh B., Sexena V. Transverse vibration of skew plates with variable thickness. Journal of Sound and Vibration, Vol. 206, Issue 1, 1997, p. 1-13.

Publisher
Sharma A., Sharma A. K., Raghav A. K., Kumar V. Vibrational study of square plate with thermal effect and circular variation in density. Romanian Journal of Acoustics and Vibration, Vol. 13, Issue 2, 2016, p. 146-152.

Search CrossRef
Sharma A., Kumar V., Raghav A. K. Vibrational frequency of circular tapered square plate. Romanian Journal of Acoustics and Vibration, Vol. 14, Issue 1, 2017, p. 21-27.

Search CrossRef
Sharma A. Vibrational frequencies of parallelogram plate with circular variations in thickness. Soft Computing: Theories and Application, Advances in Intelligent System and Computing, Vol. 583, 2018, p. 317-326.

Publisher
Sharma A., Raghav A. K., Kumar V., Sharma A. K. Effect of circular variation in thickness and linear variation in density on vibrational frequencies. Progress in Advanced Computing and Intelligent Engineering, Vol. 564, 2018, p. 501-509.

Publisher
Khanna A. K., Kaur N. A study on vibration of tapered rectangular plate under non-uniform temperature field. Mechanika, Vol. 20, Issue 4, 2014, p. 376-381.

Publisher
Khanna A. K., Kaur N. Effect of non-homogeneity on free vibration of visco-elastic rectangular plate with varying structural parameters. Journal of Vibroengineering, Vol. 15, Issue 4, 2013, p. 2146-2155.

Search CrossRef
Hwu C., Chang W. C., Gai H. S. Vibration suppression of composite sandwich beams. Journal of Sound and Vibration, Vol. 272, Issues 1-2, 2004, p. 1-20.

Publisher
Avalos D. R., Laura P. A. A. Transverse vibration of a simply supported plate of generalized anisotropy with an oblique cut-out. Journal of Sound and Vibration, Vol. 258, Issue 4, 2002, p. 773-776.

Publisher
Shufrin I., Eisenberger M. Semi-analytical modeling of cutouts in rectangular plates with variable thickness–Free vibration analysis. Applied Mathematical Modelling, Vol. 40, Issues 15-16, 2016, p. 6983-7000.

Publisher
Sakiyama T., Huang M. Free vibrational analysis of rectangular plate with variable thickness. Journal of Sound and Vibration, Vol. 216, Issue 3, 1998, p. 379-397.

Publisher
Wang Y. Q., Zu J. W. Vibration behavior of functionally graded rectangular plate with porosities and moving in thermal environment. Aerospace Science and Technology, Vol. 69, 2017, p. 550-562.

Publisher
Huang M., Ma X. Q., Sakiyama T., Matuda H., Morita C. Free vibration of orthotropic rectangular plates with variable thickness and general boundary conditions. Journal of Sound and Vibration, Vol. 288, Issues 4-5, 2005, p. 931-955.

Publisher
Huang M., Xu Y., Cao B. Free Vibration analysis of cantilever rectangular plates with variable thickness. Applied Mechanics and Materials, Vol. 130, Issue 134, 2012, p. 2774-2777.

Publisher
Taher H. R. D., Omidi M., Zadpoor A. A., Nikoyaan A. A. Free vibration of circular and annular plates with variable thickness and different combination of boundary conditions. Journal of Sound and Vibration, Vol. 296, Issues 4-5, 2006, p. 1084-1092.

Publisher
Zur K. K. Green’s function for frequency analysis of thin annular plates with nonlinear variable thickness. Applied Mathematical Modelling, Vol. 40, Issues 5-6, 2016, p. 3601-3619.

Publisher
Zur K. K. Green’s function in frequency analysis of circular thin plates with variable thickness. Journal of Theoretical and Applied Mechanics, Vol. 53, Issue 4, 2015, p. 873-884.

Publisher

Cited by

Development of an analytical model for the flexural vibration of fish bone active camber structures with truncated, variable thickness partitions

Mahdi Nejati | Saeed Shokrollahi | Masoud Cheraghi

(2024)

Thermal Vibration of Rectangle Plate with Two-Dimensional Circular Thickness Effect

Amit Sharma | Naveen Mani | Reeta Bhardwaj | Kamal Kumar | Francesco Pellicano

(2022)

Effect of linear variation in density and circular variation in Poisson’s ratio on time period of vibration of rectangular plate

Anup Kumar | Neeraj Lather | Reeta Bhardwaj | Naveen Mani | Amit Sharma

(2018)

About this article

Received

30 December 2017

Accepted

26 August 2018

Published

15 February 2019

SUBJECTS

Modal analysis and applications

DOI

https://doi.org/10.21595/jve.2018.19589

Keywords

rectangular plate

exponential variation

natural vibration

non homogeneity

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Previous article in issue Previous Next article in issue Next

Research article

2023 11 27

Vibrational analysis of orthotropic rectangular plate having combination of circular thickness and parabolic temperature

Neeraj Lather, Ankit Kumar, Madhu Gupta, Pawan Joshi, Amit Sharma

Research article

2019 09 30

Analytical solution of a vibrational problem for visco-elastic plate with Kelvin type boundary conditions

Ashish Kumar Sharma, Manoj Kumar Dhiman, Silky Bensal

Research article

2019 03 15

Natural vibration of skew plate on different set of boundary conditions with temperature gradient

Neeraj Lather, Amit Sharma

Research article

2019 03 15

Modelling on vibration of skew plate with thickness and temperature variation

Reeta Bhardwaj, Naveen Mani

A. Sharma, N. Mani, and R. Bhardwaj, “Natural vibration of tapered rectangular plate with exponential variation in non homogeneity,” Journal of Vibroengineering, Vol. 21, No. 1, pp. 187–197, Feb. 2019, https://doi.org/10.21595/jve.2018.19589

Copy Extrica

Copied to clipboard!

TY  - JOUR
DO  - 10.21595/jve.2018.19589
UR  - https://doi.org/10.21595/jve.2018.19589
TI  - Natural vibration of tapered rectangular plate with exponential variation in non homogeneity
T2  - Journal of Vibroengineering
AU  - Sharma, Amit
AU  - Mani, Naveen
AU  - Bhardwaj, Reeta
PY  - 2019
DA  - 2019/02/15
PB  - JVE International Ltd.
SP  - 187-197
IS  - 1
VL  - 21
SN  - 1392-8716
SN  - 2538-8460
ER  - 

Copy Ris

Copied to clipboard!

@article{Sharma_2019,
	doi = {10.21595/jve.2018.19589},
	url = {https://doi.org/10.21595/jve.2018.19589},
	year = 2019,
	month = {feb},
	publisher = {{JVE} International Ltd.},
	volume = {21},
	number = {1},
	pages = {187--197},
	author = {Amit Sharma and Naveen Mani and Reeta Bhardwaj},
	title = {Natural vibration of tapered rectangular plate with exponential variation in non homogeneity},
	journal = {Journal of Vibroengineering}
}

Copy Bibtex

Copied to clipboard!

[1]A. Sharma, N. Mani, and R. Bhardwaj, “Natural vibration of tapered rectangular plate with exponential variation in non homogeneity,” Journal of Vibroengineering, vol. 21, no. 1, pp. 187–197, Feb. 2019, doi: 10.21595/jve.2018.19589.

Copy IEEE

Copied to clipboard!

Sharma, Amit, Naveen Mani, and Reeta Bhardwaj. “Natural Vibration of Tapered Rectangular Plate with Exponential Variation in Non Homogeneity.” Journal of Vibroengineering 21, no. 1 (February 15, 2019): 187–97. https://doi.org/10.21595/jve.2018.19589.

Copy Chicago

Copied to clipboard!

Natural vibration of tapered rectangular plate with exponential variation in non homogeneity

Abstract

1. Introduction

2. Differential equation of motion

3. Construction of problem

4. Solution of the problem

5. Results and discussion

6. Results comparison

7. Conclusions

References

Cited by

About this article

Related Articles