Published: 13 December 2018

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

Anup Kumar1
Neeraj Lather2
Reeta Bhardwaj3
Naveen Mani4
Amit Sharma5
1, 2, 3, 5Amity University Haryana, Gurgaon, India
4Sandip University, Nashik, Maharashtra, India
Corresponding Author:
Amit Sharma
Views 145
Reads 43
Downloads 1292

Abstract

In this paper, a theoretical analysis is carried out to investigate the effect of linear variation in density and circular variation in Poisson’s ratio on time period of frequency modes of rectangular plate with variable thickness under temperature field. The thickness variation is considered to be circular and temperature variation on the plate is assumed to be bi-linear. Rayleigh Ritz method is used to solve the differential equation. All the results (time period for first two modes of vibration) are presented with the help of tables.

1. Introduction

The study of vibration of non-homogeneous plate is essential in these days because non-homogeneous plate with variable thickness are used in almost all engineering structures such as power plants, wings of an aircrafts, machines, bridges etc. Plates with variable thickness along with non-homogeneity have great impact when compared to homogeneous plate with uniform thickness because of their efficiency and strength. The first few modes of vibration provide us good information about the behavior of systems/structures. Therefore, in order to design perfect structures/systems, it is essential to determine natural frequencies and mode shapes. A significant work has been reported in these directions.

The natural vibration of cantilever plates with variable thickness is analyzed by using mixed boundary grid method (FBGM) [1] and obtained characteristic equation and frequency parameters. A model [2] is presented to analyze the nonlinear vibrations of visco elastic thin rectangular plates by using von Kármán nonlinear strain–displacement relationships and obtained fundamental modes of a simply supported square plate with immovable edges. The nonlinear damping of visco elastic rectangular plate [3] is studied and equation of motion is derived using Lagrange equations. The results are also obtained theoretically as well as experimentally. Vibration analysis of rectangular plates with rectangular cutouts is investigated by using extended Hencky bar-net method (HBM) [4]. Method of reverberation ray matrix (MRRM) and golden section search (GSS) algorithm [5] is applied to obtain the exact solution of rectangular plates with arbitrary boundary conditions. Accurate analytic solutions for natural vibration of thick rectangular plates with a free edge is presented in [6]. Natural vibration of thick rectangular plate [7] without two parallel simply supported edges is studied and new analytic solutions are obtained. The effect of crack defects and temperature on vibration of thin isotropic and orthotropic rectangular plates is studied in [8]. The effect of circular variation in density and exponential variation in Poisson’s ratio on vibrational frequency of parallelogram plate under temperature field is examined in [9]. The effect of two-dimensional thickness and temperature effect on natural vibration of parallelogram plate on clamped edges is studied using Rayleigh Ritz method [10]. Rayleigh Ritz method [11] is used to analyze the natural vibration of isotropic rectangular plate with circular variation in thickness and exponential variation in Poisson’s ratio under temperature field. Free vibration of moderately thick laminated composite rectangular plate with non-uniform boundary conditions is presented by using an improved Fourier series method [12]. Closed-form solution [13] is presented to study the effects of rotary inertia and shear deformation on frequency of natural vibration problems. A mathematical model [14] is presented to analyze the free vibration of rectangular plates with various rectangular cutouts and variable thickness.

In this study, authors show the effect of linear variation in density and circular variation in Poisson’s ratio on time period of vibration of non homogeneous rectangular plate on clamped edges. Authors also calculate the time period of vibration corresponding to circular variation in thickness and bi linear variation in temperature. The results are presented with the help of tables.

2. Analysis and assumptions

The differential equation for transverse motion of the plate with variable D1 and ν i.e., flexural rigidity and Poisson’s ratio is given by:

1
D14Φζ4+24Φζ2ψ2+4Φψ4+2νζ22Φψ2+2D1ζ3Φζ3+3Φζψ2+νζ2Φψ2
+2D1ζ22Φζ2+ν2Φψ2+2D1ψ3Φψ3+3Φψζ2-νζ2Φζψ
+2D1ψ22Φψ2+ν2Φζ2+21-ν2D1ζψ2Φζψ-ρlω2Φ=0,

where Φis known as deflection function. The expression for flexural rigidity is D1=El3/121-ν2, where E, l are known as Young’s modulus and thickness of the plate.

In order to avoid complexity, the present study requires some assumptions as:

a) Since the plate has variable thickness l, therefore authors considered one dimensional circular variation in thickness as shown in Fig. 1 as:

2
l=l01+β1-1-ζ2a2,

where β, (0 ≤ β ≤ 1) is known as tapering parameter. Thickness of plate becomes constant at ζ= 0.

b) For non-homogeneity consideration, authors assumed one dimensional linear variation in density and one dimensional circular variation in Poisson’s ratio as:

3
ρ=ρ01+m1ζa,
4
ν=ν01-m21-1-ζ2a2,

where m1, (0 ≤ m1 ≤ 1) and m2, (0 ≤ m2 < 1) are known as non-homogeneity constants corresponding to density and Poisson’s ratio.

c) The plate is subjected to steady two-dimensional linear temperature distributions as:

5
τ=τ01-ζa1-ψ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:

6
E=E01-γτ,

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

Using Eq. (5), Eq. (6) becomes:

7
E=E01-α1-ζa1-ψb,

where α, (0 ≤ α < 1) is called temperature gradient, which is the product of temperature at origin and slope of variation i.e., α=γτ0.

Fig. 1Rectangular plate with one dimensional circular variation

Rectangular plate with one dimensional circular variation

3. Solution for frequency equation and time period

We are using Rayleigh Ritz technique (i.e., maximum strain energy Vs must equal to maximum kinetic energy Ts) in order to obtain frequency equation and time period for both modes of vibrations. Therefore, we must have:

8
δVs-Ts=0.

Here the expression for Vs and Ts are given by:

9
Vs=120a 0bD1×2Φζ22+2Φψ22+2ν2Φζ22Φψ2+21-ν2Φζψ2dψdζ,
10
Ts=12ω20a 0bρlΦ2dψdζ.

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

11
Φζ,ψ=Φζ,ψζ=0, ζ=0, a, Φζ,ψ=Φζ,ψψ=0, ψ=0, b.

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

12
Φζ,ψ=ζa2ψb21-ζa21-ψb2Ω1+Ω2ζaψb1-ζa1-ψb,

where Ω1 and Ω2 are arbitrary constants.

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

13
ζ1=ζa, ψ1=ψa.

Using Eqs. (2)-(4), (7), (9), (10) and (13), Eq. (8) becomes:

14
δVs*-λ2Ts*=0,

where:

Vs*=010ba1-α1-ζ11-aψ1b1+β1-1-ζ1231-ν021-m21-1-ζ122
2Φζ122+2Φψ122+2ν01-m21-1-ζ122Φζ122Φψ12+21-ν01-m21-1-ζ122Φζ1ψ12dψ1dζ1,
Ts*=01 0ba1+m1ζ11+β1-1-ζ12Φ2dψ1dζ1,

and λ2=12ρ0ω2a4/E0l02 is known as frequency parameter. Eq. (14) consists of two unknown constants Ω1 and Ω2 (because of substitution of deflection function Φζ,ψ). These two unknowns could be calculated as follows:

15
ΩnVs*-λ2Ts*=0, n=1,2.

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

16
c11Ω1+c12Ω2=0,c21Ω2+c22Ω2=0.

To obtain non-zero solution (frequency equation), the determinant of coefficient matrix obtain from Eq. (16) must be zero i.e.,

17
c11c12c21c22=0.

From Eq. (17), we get a quadratic equation from which we get frequency modes.

The time period of frequency modes is calculated as:

18
K=2πλ,

where λ is frequeny modes obtained from Eq. (17).

4. Results and discussion

The effect of variation of different plate parameters (non-homogeneity constants m1, m2, thermal gradient α and tapering parameter β) on time period (in seconds) of vibrations are calculated for fixed value of aspect ratio a/b= 1.5 and presented in tabular form.

Table 1 provides the time period K of vibration corresponding to non-homogeneity m1 for two cases i.e., β=m2=α= 0.2 and β=m2=α= 0.6. From Table 1, we can see that time period is increasing with the increasing value of non-homogeneity constant m1 for both the cases. While the time period is decreasing with the combined increasing value of taper constant β, non-homogeneity m2 and thermal gradient α.

Table 2 accommodates the time period K corresponding to non-homogeneity constant m2 for two different cases i.e., m1=β=α= 0.2 and m1=β=α= 0.6. Here, time period is increasing (almost constant) for first case. But for second case, time period is increasing for first mode and time period is decreasing for second mode. The time period is also decreasing when the combined value of thermal gradient α, taper constant β and non-homogeneity constant m1 varies from 0.2 to 0.6.

Table 1Non-homogeneity m1 vs. time period K[s] for a/b= 1.5

m1
α=β=m2= 0.2
α=β=m2= 0.6
λ1
λ2
λ1
λ2
0.0
0.02214
0.08729
0.01848
0.07226
0.2
0.02323
0.09157
0.01941
0.07584
0.4
0.02428
0.09566
0.02030
0.07925
0.6
0.02528
0.09959
0.02115
0.08253
0.8
0.02646
0.10336
0.02197
0.08568
1.0
0.02717
0.10701
0.02276
0.08872

Table 2Non-homogeneity m2 vs. time period K[s] for a/b= 1.5

m2
α=β=m1= 0.2
α=β=m1= 0.6
λ1
λ2
λ1
λ2
0.0
0.02322
0.09155
0.02113
0.08272
0.2
0.02323
0.09157
0.02114
0.08269
0.4
0.02324
0.09158
0.02114
0.08263
0.6
0.02325
0.09156
0.02115
0.08253
0.8
0.02326
0.09151
0.02115
0.08239

Table 3Taper constant β vs. time period K[s] for a/b= 1.5

β
α=m1=m2= 0.4
α=m1=m2= 0.8
λ1
λ2
λ1
λ2
0.0
0.02837
0.11209
0.03252
0.12384
0.2
0.02492
0.09815
0.02851
0.11193
0.4
0.02211
0.08677
0.02524
0.09864
0.6
0.01979
0.07741
0.02256
0.08776
0.8
0.01786
0.06963
0.02033
0.07877
1.0
0.01623
0.06310
0.01846
0.07125

Table 4Thermal gradient α vs. time period K[s] for a/b= 1.5

α
β=m1=m2= 0.4
β=m1=m2= 0.8
λ1
λ2
λ1
λ2
0.0
0.02107
0.08273
0.01848
0.07183
0.2
0.02157
0.08468
0.01890
0.07339
0.4
0.02211
0.08677
0.01934
0.07505
0.6
0.02269
0.08903
0.01982
0.07684
0.8
0.02333
0.09148
0.02033
0.07877

Time period K corresponding to taper constant β is tabulated in Table 3 for the two different cases i.e., m1=m2=α= 0.4 and m1=m2=α= 0.8. Time period is decreasing with the increasing value of taper constant β for both the cases. When we move from first case to second case, the time period is increasing with the increasing value of taper constant β.

Table 4 presents the time period K corresponding to thermal gradient α for two different cases i.e., m1=m2=β= 0.4 and m1=m2=β= 0.8. From Table 4, we conclude that time period is increasing with the increasing value of thermal gradient α for both the cases and decreasing with the combined increasing value of taper constant β and non-homogeneity constants m1, m2.

5. Conclusions

The present study reveals the effect of plate parameters especially density and Poisson’s ratio on time period of vibration of clamped rectangular plate with variable thickness under temperature field. The rate of increment in time period corresponding to non-homogeneity m1 (linear variation) is higher than the rate of increment in time period corresponding to non-homogeneity m2 (circular variation). Initially, time period (at m2= 0) is high in case of circular variation in non-homogeneity when compared to time period (at m1= 0) in case of linear variation in non-homogeneity. The time period is decreasing corresponding to circular variation in thickness and increasing corresponding to linear temperature variation on plate.

References

  • 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.
  • Amabili M. Nonlinear vibrations of viscoelastic rectangular plates. Journal of Sound and Vibration, Vol. 362, 2016, p. 142-156.
  • Amabili M. Nonlinear damping in nonlinear vibrations of rectangular plates: derivation from viscoelasticity and experimental validation. Journal of the Mechanics and Physics of Solids, Vol. 118, 2018, p. 275-292.
  • Zhang Y. P., Wang C. M., Pedroso D. M., Zhang H. Extension of Hencky bar-net model for vibration analysis of rectangular plates with rectangular cutouts. Journal of Sound and Vibration, Vol. 432, 2018, p. 65-87.
  • Yuansheng Z., Qingshan W., Dongyan S., Qian L., Zhongyu Z. Exact solutions for the free in-plane vibrations of rectangular plates with arbitrary boundary conditions. International Journal of Mechanical Sciences, Vol. 130, 2017, p. 1-10.
  • Li R., Wang P., Riye X., Xu G. New analytic solutions for free vibration of rectangular thick plates with an edge free. International Journal of Mechanical Sciences, Vol. 131, Issue 132, 2017, p. 179-190.
  • Li R., Wang P., Zheng X., Wang B. New benchmark solutions for free vibration of clamped rectangular thick plates and their variants. Applied Mathematics Letters, Vol. 78, 2018, p. 88-94.
  • Lai S. K., Zhang L. H. Thermal effect on vibration and buckling analysis of thin isotropic/orthotropic rectangular plates with crack defects. Engineering Structures, Vol. 177, 2018, p. 444-458.
  • Sharma A., Sharma A. K., Kumar V. Effect of density and Poisson’s ratio on thermal induced vibration of parallelogram plate. Journal of Vibroengineering, Vol. 20, Issue 3, 2018, p. 1288-1298.
  • SharmaA. 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.
  • Sharma A., Mani N., Bhardwaj R. Natural vibration of tapered rectangular plate with exponential variation in non homogeneity. Journal of Vibroengineering, 2018, (in Press).
  • Zhang H., Shi D., Wang Q. An improved Fourier series solution for free vibration analysis of the moderately thick laminated composite rectangular plate with non-uniform boundary conditions. International Journal of Mechanical Sciences, Vol. 121, 2017, p. 1-20.
  • 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.
  • 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.

About this article

Received
02 November 2018
Accepted
14 November 2018
Published
13 December 2018
SUBJECTS
Mechanical vibrations and applications
Keywords
density
Poisson’s ratio
rectangular plate
time period