#### 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 in *x*-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.

#### 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 (Figure 1). 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.

#### 2. Materials and Methods

##### 2.1. Analysis of Motion

Fourth-order differential equation of motion for isotropic rectangular plate in Cartesian coordinate is [14] where is the flexural rigidity of the plate’s material, is the deflection function, is Poisson ratio, is mass per unit volume, and is thickness of the plate.

It is assumed that thickness of the rectangular plate varies parabolically in both directions [11]; that is, where and are length and breadth of rectangular plate, respectively, and and are taper parameters in -direction and -direction, respectively.

The authors assumed biparabolic temperature variations as where denotes the temperature excess above the reference temperature at any point on the plate and denotes the temperature excess above the reference temperature at .

The temperature dependence of the modulus of elasticity for most of engineering materials can be expressed as follows [15]: where is value of Young’s modulus at reference temperature, that is, , and is slope of variation of and . On using (3) in (4), one obtains where () is thermal gradient.

On substituting the values of and from (2) and (5), respectively, expression of flexural rigidity () becomes

##### 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 () is taken as shown in Table 1.

##### 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 must be equal to the maximum kinetic energy (). So it is necessary for the problem under consideration that [16] Now assume the nondimensional variables as The expressions for kinetic energy () and strain energy () are On using (9) in (7), one gets where Here, is a frequency parameter.

Equation (10) consists of two unknown constants, that is, and , arising due to the substitution of . These two constants are to be determined as follows: On simplifying (12), one gets where and 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 With the help of (14), one can obtain a quadratic equation in 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 ( and ), and aspect ratio () 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 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 ( and ) and thermal gradient () at fixed aspect ratio . Also, frequency is tabulated for various values of aspect ratio in Table 5 for the following cases: (i) = = = 0.0, (ii) = = = 0.2, (iii) = = = 0.4, (iv) = = = 0.6, and (v) = = = 0.8.

In Tables 2–4, the authors found that frequency for both the modes of vibration continuously increases as taper constants ( and ) 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 Table 5, 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 , and increase from 0.0 to 0.8.

#### 4. Comparison

In Tables 6 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 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 -direction and parabolic in -direction. Therefore, the authors compared both the modes of frequency of present paper with [17] for increasing values of only at , 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 but the second mode of frequency is greater than [17] for , and at , 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 -direction and parabolic in -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.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.