Published: 30 June 2015

Rayleigh-Ritz method for analyzing free vibration of orthotropic rectangular plate with 2D thickness and temperature variation

Subodh Kumar Sharma1
Ashish Kumar Sharma2
1Government Post Graduate College, Ambala Cantt, Haryana, India
2Research Scholar Department of Mathematics, Pacific University, Udaipur, Rajasthan, India
Corresponding Author:
Subodh Kumar Sharma
Views 538
Reads 236
Downloads 1272

Abstract

A plate structure is of key interest to aerospace, mechanical and civil engineers. Vibration reduction is a major challenge pertaining to these fields; especially in aerospace applications, such a reduction must be achieved with a minimal increase in weight. The vibration of plates is a special case of the more general problem of mechanical vibrations. In this paper effect of bi-parabolic variation in temperature is premeditated on vibration of an orthotropic rectangular plate and whose thickness also varies bi-linear as in two dimensional. Frequency equation is derived by using Rayleigh-Ritz technique with a two-term deflection function. Time period, Deflection and Logarithmic decrement at different points for the first two modes of vibration are calculated for various values of thermal gradients, aspect ratio and taper constants.

1. Introduction

With the upsurge of interest by the aerospace industry in using advanced composites for primary load-bearing structures, accurate design and analysis methods are becoming more critical. Most designs and analysis techniques for composites assume orthotropic elastic behavior and for case where most of the load is carried by the fibers, the assumption is satisfactory. However, properties transverse to the fibers and shear properties are matrix controlled and exhibit strong elastic/visco-elastic behavior. There are many cases where the transverse and shear stiff nesses significantly affect the material response.

With the advancement of space technology, the need of the study of vibrations of plates of certain aspect ratios with some simple restraints on the boundaries has increased. The development of solid propellant rocket motors, increased use of soft filaments in aerospace structures and the building activities in the cold regions have intensified the need for solutions of various problems of plates continuously supported by elastic or visco-elastic media. The information for the first few modes of vibration is essential for a construction engineer before finalizing a design.

The term vibration describes repetitive motion that can be measured and observed in a structure. Unwanted vibration can cause fatigue or degrade the performance of the structure. Therefore it is desirable to eliminate or reduce the effects of vibration. In other cases, the goal may be to understand the effect on the structure or to control or modify the vibration, or to isolate it from the structure and minimize the structural responses.

Recently, Leissa [1, 2] has given the solution for rectangular plate of variable thickness. Gupta, Johri and Vats [3] have discussed the thermal effect on vibration of non-homogeneous orthotropic rectangular plate having bi-directional parabolically varying thickness. Gupta and Khanna [4] have solved the problem of free vibration of visco-elastic rectangular plate with linearly thickness variations in both directions. Singh and Saxena [5] have discussed the transverse vibration of rectangular plate with bi-directional thickness variation. Sobotka [6] has investigated the vibration of rectangular orthotropic visco-elastic plates. Lal [7] studied transverse vibrations of orthotropic non-uniform rectangular plates with continuously varying density. Warade and Deshmukh [8] discussed thermal deflection of a thin clamped circular plate due to partially distributive heat supply. Sobotka [9] discussed rheology of orthotropic visco-elastic plates. Gupta and Kumar [10] analyzed vibration of non-homogeneous visco-elastic rectangular plates with linearly varying thickness. Khanna A., Kaur N., and Sharma A. K. [11] have discussed effect of varying poisson ratio on thermally induced vibrations of non-homogeneous rectangular plate. Sharma S. K., and Sharma A. K. [12] have discussed the mechanical vibration of orthotropic rectangular plate with 2D linearly varying thickness and thermal Effect. Khanna A., and Sharma A. K. [13] have solved the problem on vibration analysis of visco-elastic square plate of variable thickness with thermal gradient. Kumar Sharma A., and Sharma S. K. [14] have discussed the vibration computational of visco-elastic plate with sinusoidal thickness variation and linearly thermal effect in 2d. Khanna A., Kumar A., and Bhatia M. [15] has investigated the computational prediction on two dimensional thermal effect on vibration of visco-elastic square plate of variable thickness. Khanna A., and Sharma A. K. [16] studied natural vibration of visco-elastic plate of varying thickness with thermal effect. Kumar Sharma A., and Sharma S. K. [17] discussed free vibration analysis of visco-elastic orthotropic rectangular plate with bi-parabolic thermal effect and bi-linear thickness variation. Sharma S. K. and Sharma A. K. [18] discussed effect of bi-parabolic thermal and thickness variation on vibration of visco-elastic orthotropic rectangular plate. Khanna A., and Sharma A. K. [19] Mechanical Vibration of Visco-Elastic Plate with Thickness Variation. Khanna A., Kaur N., and Sharma A. K. [20] Effect of varying poisson ratio on thermally induced vibrations of non-homogeneous rectangular plate. Khanna A., and Sharma A. K. [21] Analysis of free vibrations of visco elastic square plate of variable thickness with temperature effect.

The aim is to study two dimensional thermal effects on the vibration of visco-elastic rectangular plate whose thickness varies linearly in both directions and temperature varies bi-parabolically in another direction. All results presents in graphical form.

2. Methodology

Consider a plate which is rectangular in shape and assumed to be made up of orthotropic material. Let a and b be the length and breadth of the plate.

Let the plate under consideration is subjected to a steady two dimensional parabolic temperature distribution τ along x-axis and y-axis, then as [2]:

1
τ=τ01-x2a21-y2b2,

where τ the temperature is excess above the reference temperature at a distance x/a and τ0 is the temperature excess above the reference temperature at the end of the plate. For most orthotropic materials, moduli of elasticity are defined as [3]:

2
Exτ=E11-γτ,Eyτ=E21-γτ,Gxyτ=G01-γτ,

where Ex and Ey are Young’s moduli in x-and y-directions respectively, Gxy is shear modulus and γ is slope of variation of moduli with temperature. Using Eq. (1) in Eq. (2), one has:

3
Exx=E11-α1-x2a21-y2b2,Eyx=E21-α1-x2a21-y2b2,Gxyx=G01-α1-x2a21-y2b2,

where α=γT0(0α<1), a thermal gradient.

The governing differential equation of transverse motion of an orthotropic rectangular plate of variable thickness in Cartesian coordinate as [12]:

4
Dx4wx4+Dy4wy4+2H2wx22wy2+2Hxxwx2wy2+2Hyy2wx2wy+2Dxx3wx3
+2Dyy3wy3+2Dxx22wx2+2Dyy22wy2+2D1x22wy2+2D1y22wx2+42Dxyxy2wxy
+ρh2wt2=0,

where Dx and Dy are flexural rigidities in x-and y-directions respectively and Dxy is torsional rigidity. Here:

5
Dx=Exh3121-vxvy, Dy=Eyh312(1-vxvy), Dxy=Gxyh312,
6
D1=vxDy=vyDx,
7
H=D1+2Dxy,

where vx and vy are Poisson’s ratio.

For free transverse vibrations of the plate, wx, y, t can be defined as:

8
wx,y,t=Wx,yTt,

where p is radian frequency of vibration. Two terms deflection function for clamped rectangular plate is written as [5]:

9
Wx,y=A1xa2yb21-xa21-yb2+A2xa3yb31-xa31-yb3,

where A1 and A2 are constants to be evaluated. On using Eq. (3) in Eq. (5), we have:

10
Dx=E1h3121-νxνy1-α1-x2a21-y2b2,
Dy=E2h3121-νxνy1-α1-x2a21-y2b2,
Dxy=G0h3121-α1-x2a21-y2b2.

When the plate is executing transverse vibration of mode shape W(x, y) then Strain energy V and Kinetic energy T1 are respectively expressed as [16]:

11
S=120a0bDx2Wx22+Dy2Wy22+2D12Wx22Wy2+4Dxy2Wxy2dydx,
12
V=12p2ρ0a0bhW2dydx,

where ρ is the mass density. Using Eqs. (6), (9) and (10) in Eq. (11), we have:

13
S=12E112(1-νxνy)0a0bh31-α1-x2a21-y2b2
×2Wx22+E2E12Wy22+2vxE2E12Wx22Wy2
+4G0E1(1-vxvy)2Wxy2dydx.

In the present study, variation in thickness h of the plate is assumed to be varying parabolic in both directions, i.e.:

14
h=h01+β1xa1+β2yb,

where h0=h at x= 0 and y= 0.

Using Eq. (14) in Eqs. (13) and (12), one has:

15
S=12E1h0312(1-νxνy)0a0b1-α1-x2a21-y2b2
×1+β1xa31+β2yb3×2Wx22+E2E12Wy22+2vxE2E12Wx22Wy2
+4G0E1(1-vxvy)2Wxy2dydx,
16
V=12h0ρp20a0b1+β1xa1+β2ybW2dydx.

3. Method of solution

The solution to the current problem is given by the application of Rayleigh-Ritz method. In order to apply their procedure, maximum Strain energy must be equal to maximum Kinetic energy. Therefore it is desired that following equation must be satisfied:

17
δS-λ2V=0.

On substituting the values of ‘S’ and ‘V’ from Eqs. (15) and (16) in Eq. (17), we have:

18
δS1-λ2V1=0,

where:

19
λ2=12a4ρp2(1-vxvy)E1h02,

is a frequency parameter:

20
S1=0a0b1-α1-x2a21-y2b2×1+β1xa31+β2yb3
×2Wx22+E2E12Wy22+2vxE2E12Wx22Wy2
+4G0E1(1-vxvy)2Wxy2dydx,
21
V1=0a0b1+β1xa1+β2ybW2dydx.

4. Boundary condition and equation of frequency

For a clamped rectangular plate, boundary conditions are:

22
W=Wx=0,x=0,a,W=Wy=0,y=0,b.

Eq. (18) contains two unknown parameters A1 and A2 to be evaluated. Values of these constants may be evaluated by the following procedure:

23
AqS1-λ2V1=0, q=1, 2.

On simplifying Eq. (23), one gets the following equation:

24
cq1A1+cq2A2=0,

where cq1 and cq2 involves the parametric constants and the frequency parameter. For a non-zero solution, determinant of coefficients of Eq. (24) must vanish. In this way frequency equation comes out to be:

25
c11c21 c12c22=0.

On solving Eq. (25) one gets a quadratic equation in p2, so it will give two roots. On substituting the value of A1=1, by choice, in Eq. (9) one get A2=-c11/c12 and hence W becomes:

26
Wx,y=xa2yb21-xa21-yb2+-c11c12xa3yb31-xa31-yb3.

5. Time function

Time functions of free vibrations of visco-elastic plates are defined by the general ordinary differential equation:

27
T ¨+p2D~T=0.

Their form depends on visco-elastic operator and which for Kelvin’s model, one can be taken as:

28
D~=1+ηGddt,

where η is visco-elastic constant and G is shear modulus. Taking temperature dependence of visco-elastic constant η and shear modulus G is the same form as that of Young’s moduli, we have:

29
Gτ=G01-γ1τ,ητ=η01-γ2τ,

where G0 is shear modulus and η0 is visco-elastic constant at some reference temperature, i.e., T= 0, γ1 and γ2 are slope variation of τ from Eq. (1) in Eq. (29), one gets:

30
G=G01-α11-x2a21-y2b2, α1= γ1 τ0, 0α11,
η=η01-α21-x2a21-y2b2, α2= γ2 τ0, 0α21,

here α1 and α2 are thermal constants.

After using Eq. (28) in Eq. (27), one gets:

31
T¨+p2kT˙+p2T=0,

where:

32
k=ηG=η0α21-x2a21-y2b2G0α11-x2a21-y2b2.

Eq. (31) is a differential equation of second order for time function T.

On solving Eq. (31), one gets:

33
Tt=ea1tC1cosb1t+C2sinb1t,

where:

34
a1=-p2k2,b1=p1-pk22,

and C1, C2 are constants which can be determined easily from initial conditions of the plate.

Let us take initial conditions as:

35
T=1, T˙=0, at t=0.

Using Eq. (35) in Eq. (33), one gets:

36
C1=1, C2=p2ηG2b1=-a1b1.

Using Eq. (36) in Eq. (33), one gets:

37
Tt=ea1tcosb1t+-a1b1sinb1t.

On using Eqs. (26) and (37) in Eq. (8), one gets:

38
w=xa2yb21-xa21-yb2+-c11c12xa3yb31-xa31-yb3
×xea1tcosb1t+-a1b1sinb1t.

Time period of the vibration of the plate is given by:

39
K=2πp,

where p is the frequency given by Eq. (19).

Logarithmic decrement of the vibrations given by the standard formula:

40
Λ=logew2w1,

where w1 is the deflection at any point on the plate at time period K=K1 and w2 is the deflection at same point at the time period succeeding K1.

6. Numerical evaluations

For calculations, the material of ‘Duralium’ which is an alloy of Aluminium, Copper, Magnesium and Manganese have been taken. Computations have been made for calculating the values of logarithmic decrement (Λ), time period (K) and deflection (w) for a isotropic visco-elastic rectangular plate for different values of taper constants β1 and β2 and aspect ratio a/b at different points for first two modes of vibrations. In calculations, the following material parameters are used: E= 7.08×1010 N/M2, G= 2:632×1010 N/M2, η= 14:612×105 Ns/M2, ρ= 2:80×103 kg/M3, v= 0:345.

The thickness of the plate at the center is taken as h0= 0.01 m.

7. Results and discussion

For calculating the values of time period (K) for a orthotropic rectangular plate with different values of aspect ratio (a/b), thermal gradient (α) and taper constant (β) for first two modes of vibrations. In the present problem, latest software technology “MAT LAB” is used to get the numerical results with great accuracy and concentration. We had considered the various cases of time period against taper constant, aspect ratio and thermal gradient which are stated as below:

Fig. 1 illustrates the result of time period K with different values of thermal gradient α for first two modes of vibration. It can be seen from figures that as thermal gradient α increases time period increases time period decreases for both modes of vibration.

Fig. 2 illustrates the result of time period K with different values of aspect ratio a/b=1.5 for first two modes of vibration. It can be seen from figures that as aspect ratio increases then time period decreases for both modes of vibration.

Fig. 1Vibration of time period K∙10-5 with different values of thermal gradient α and aspect ratio a/b=1.5

Vibration of time period K∙10-5 with different values of thermal gradient α  and aspect ratio a/b=1.5

Fig. 2Vibration of time period K with different values of aspect ratio (a/b)

Vibration of time period K with different values of aspect ratio (a/b)

Fig. 3Time period of variation K with different values of taper constant β1 and constant aspect ratio (a/b= 1.5)

Time period of variation K with different values of taper constant β1  and constant aspect ratio (a/b= 1.5)

Figs. 3-4 illustrates the results of time period K with different values of taper constants β1 and β2 for first two modes of vibration. It can be seen from figures that as β1 and β2 increases then time period decreases for both modes of vibration.

Fig. 5 illustrate the result of deflection w for first two modes of vibration for aspect ratio a/b=1.5 with other different values: β1=β2= 0.0, α= 0.0, α1= 0.2, α2= 0.3, Y= 0.6 and time = 0K and 5K. It is interesting to see that as aspect ratio increases from 0.1 to 0.5 then deflection w increases but as the value of aspect ratio increases from 0.5 to 1.0 then we clearly see that deflection w decreases for both modes of vibration.

Fig. 4Time period of variation K with different values of β2 and constant aspect ratio (a/b= 1.5)

Time period of variation K with different values of β2 and constant aspect ratio (a/b= 1.5)

Fig. 5Deflection (w∙10-5) of a clamped visco-elastic rectangular plate for different values of X and Y, a constant aspect ratio (a/b= 1.5) and β1=β2= 0.0, α= 0.0, α1= 0.2, α2= 0.3 and time = 0K and 5K

Deflection (w∙10-5) of a clamped visco-elastic rectangular plate for different values of X and Y, a constant aspect ratio (a/b= 1.5) and β1=β2= 0.0, α= 0.0, α1= 0.2, α2= 0.3 and time = 0K and 5K

Fig. 6 illustrate the result of deflection w for first two modes of vibration for aspect ratio a/b= 1.5 with other different values: β1=β2=α=α1=α2= 0.0 and time = 0K and 5K, X=Y= 0.6. It is interesting to see that as aspect ratio increases then deflection increases for both modes of vibration.

Fig. 7 illustrate results for logarithmic decrement Λ for aspect ratio a/b= 1.5 for first two modes of vibration for different values of taper constant β1. It is interesting to note that as taper constant increases then logarithmic decrement decreases.

The accuracy of the present computations as shown in above figures is compared with the published results [18]. Figs. 1-7 shows a comparison of the values of frequency parameter obtained in the present problem and published paper of the authors [18]. A very close agreement is seen between the present results and of the published paper in which effect of bi-parabolic thermal and thickness variation on vibration of visco-elastic orthotropic rectangular plate has been studied.

Fig. 6Deflection (w∙10-5) of a clamped visco-elastic rectangular plate for different values of aspect ratio (a/b) and β1=β2=α=α1=α2= 0.0

Deflection (w∙10-5) of a clamped visco-elastic rectangular plate for different values of aspect ratio (a/b) and β1=β2=α=α1=α2= 0.0

Fig. 7Logarithmic decrement (Λ) of a clamped visco-elastic rectangular plate for different values of taper constant β1 and β2 and α=α1=α2= 0.0 and aspect ratio a/b= 1.5

Logarithmic decrement (Λ) of a clamped visco-elastic rectangular plate for different values of taper constant β1 and β2 and α=α1=α2= 0.0 and aspect ratio a/b= 1.5

8. Conclusions

Our main aim is to provide such kind of a mathematical design so that scientist can perceive their potential in mechanical engineering field and increase strength, durability and efficiency of mechanical design and structuring with a practical approach .Actually this is the need of the hour to develop more but authentic mathematical model for the help of mechanical engineers. Rectangular plates with variable thickness are of great importance in a wide variety of engineering applications i.e. construction of wings, fins of rockets, missiles. As a result, the analysis of plate’s vibrations has attracted many research works, and has been considerably improved to achieve realistic results. As space technology has advanced, the need of the study of vibration of plates of certain aspect ratios with some simple restraints on the boundaries has also increased. The information for the first few modes of vibrations is essential for a constitution engineer before finalizing a design. Therefore mechanical engineers and technocrats are advised to study and get the practical importance of the present paper and to provide much better structure and machines with more safety and economy.

References

  • Leissa A. W. Vibration of Plate. NASA SP-160, 1969.
  • Leissa A. W. Recent studies in plate vibration 1981-1985 part II, complicating effects. Shock and Vibration Digest, Vol. 19, 1987, p. 10-24.
  • Gupta A. K., Johri T., Vats R. P. Thermal effect on vibration of non-homogeneous orthotropic rectangular plate having bi-directional parabolically varying thickness. Proceeding of International Conference in World Congress on Engineering and Computer Science, San Francisco, USA, 2007, p. 784-787.
  • Gupta A. K., Khanna A. Vibration of visco-elastic rectangular plate with linearly thickness variations in both directions. Journal of Sound and Vibration, Vol. 301, 2007, p. 450-457.
  • Singh B., Saxena V. Transverse vibration of rectangular plate with bi-directional thickness variation. Journal of Sound and Vibration, Vol. 198, 1996, p. 51-65.
  • Sobotka Z. Vibration of rectangular orthotropic visco-elastic plate. Acta Techanica CSAV, 1969, p. 511-514.
  • Lal R. Transverse vibrations of orthotropic non-uniform rectangular plates with continuously varying density. Indian Journal of Pure and Applied Mathematics, Vol. 34, 2003, p. 587-606.
  • Warade R. W., Deshmukh K. C. Thermal deflection of a thin clamped circular plate due to a partially distributive heat supply. Ganita, Vol. 55, 2004, p. 179-186.
  • Sobotka Z. Free vibration of visco-elastic orthotropic rectangular plates. Acta Techanica CSAV, 1978, p. 678-705.
  • Gupta A. K., Kumar L. Thermal effect on vibration of non-homogenous visco-elastic rectangular plate of linear varying thickness. Meccanica, Vol. 43, 2008, p. 47-54.
  • Khanna A., Kaur N., Sharma A. K. Effect of varying poisson ratio on thermally induced vibrations of non-homogeneous rectangular plate. Indian Journal of Science and Technology, Vol. 5, Issue 9, 2012, p. 3263-3267.
  • Sharma S. K., Sharma A. K. Mechanical vibration of orthotropic rectangular plate with 2d linearly varying thickness and thermal effect. International Journal of Research in Advent Technology, Vol. 2, Issue 6, 2014, p. 184-190.
  • Khanna A., Sharma A. K. Vibration analysis of visco-elastic square plate of variable thickness with thermal gradient. International Journal of Engineering and Applied Sciences, Vol. 3, Issue 4, 2011, p. 1-6.
  • Kumar Sharma A., Sharma S. K. Vibration computational of visco-elastic plate with sinusoidal thickness variation and linearly thermal effect in 2D. Journal of Advanced Research in Applied Mechanics and Computational Fluid Dynamics, Vol. 1, Issue 1, 2014, p. 46-54.
  • Khanna A., Kumar A., Bhatia M. A computational prediction on two dimensional thermal effect on vibration of visco-elastic square plate of variable thickness. Proceeding of CONIAPS XIII, UPES, Deharadun, 2011.
  • Khanna A., Sharma A. K. Natural vibration of visco-elastic plate of varying thickness with thermal effect. Journal of Applied Science and Engineering, Vol. 16, Issue 2, 2013, p. 135-140.
  • Kumar Sharma A., Sharma S. K. Free vibration analysis of visco-elastic orthotropic rectangular plate with bi-parabolic thermal effect and bi-linear thickness variation. Journal of Advanced Research in Applied Mechanics and Computational Fluid Dynamics, Vol. 1, Issue 1, 2014, p. 10-23.
  • Sharma S. K., Sharma A. K. Effect of bi-parabolic thermal and thickness variation on vibration of visco-elastic orthotropic rectangular plate. Journal of Advanced Research in Manufacturing, Material Science and Metallurgical Engineering, Vol. 1, Issue 2, 2014, p. 26-38.
  • Khanna A., Sharma A. K. Mechanical vibration of visco-elastic plate with thickness variation. International Journal of Applied Mathematical Research, Vol. 1, Issue 2, 2012, p. 150-158.
  • Khanna A., Kaur N., Sharma A. K. Effect of varying poisson ratio on thermally induced vibrations of non-homogeneous rectangular plate. Indian Journal of Science and Technology, Vol. 5, Issue 9, 2012, p. 3263-3267.
  • Khanna A., Sharma A. K. Analysis of free vibrations of visco elastic square plate of variable thickness with temperature effect. International Journal of Applied Engineering Research, Vol. 2, Issue 2, 2011, p. 312-317.

About this article

Received
01 September 2014
Accepted
19 November 2014
Published
30 June 2015
SUBJECTS
Modal analysis
Keywords
visco-elastic
rectangular plate
parabolically
thermal gradient
taper constant
aspect ratio