Published: 27 November 2023

Analyzing vibration modes in non-homogeneous parallelogram plates

Sapna1
Amit Sharma2
1, 2Amity University Haryana, Gurugram, India
Corresponding Author:
Amit Sharma
Views 28
Reads 6
Downloads 119

Abstract

This study analyzes vibration modes in non-homogeneous orthotropic parallelogram plate with a one-dimensional circular thickness, focusing on SCCC edge condition, where C and S represent the clamped and simply supported edges of the plate, respectively. Circular Poisson’s ratio variation is considered, along with linear temperature changes. The study demonstrates the advantages of variable Poisson’s ratio over density parameter variation in obtaining shorter vibration time periods. Orthotropic parallelogram plate, thermal gradient, circular tapering, nonhomogeneity.

1. Introduction

The analysis of vibrational modes in various plate configurations, including tapered ones, is vital in engineering, given their wide-ranging applications. Numerous studies have contributed to this area.

Di et al. [2] introduced a theory enhancing the accuracy of predicting the behavior of thick orthotropic plates compared to classical lamination and shear deformation theories. Gupta et al. [3] investigated transverse motion in an elastic plate with non-linear thickness and thermal gradient, utilizing the Rayleigh-Ritz technique. Gupta et al. [4] introduced a model for analyzing vibrations in parallelogram-shaped, viscoelastic, orthotropic plates with linear thickness variations in both directions. Khanna [6] used the Rayleigh Ritz method to study frequency modes of a viscoelastic isotropic rectangular plate, considering the impact of thermal gradient. The natural vibration of non-homogeneous tapered parallelogram plates with temperature variation was examined in [8]. Vibrational frequencies in a non-uniformly thick parallelogram plate were mathematically analyzed, considering linear temperature effects by authors in [9]. The influence of thermal effects on non-homogeneous parallelogram plates with two-dimensional circular variations in thickness was explored in a study referenced as [10]. In [11], the authors focused on the vibration of orthotropic square plates with varying thickness, using the Rayleigh-Ritz method and MATLAB software. Sharma et al. [12] examined an orthotropic parallelogram plate with bi-linear thickness variation and a parabolic temperature distribution, specifically focusing on the SSSS edge condition.

This paper investigates the impact of one-dimensional circular tapering, Poisson’s ratio, and linear temperature on vibrational modes in a non-homogeneous orthotropic parallelogram plate under SCCC boundary conditions.

2. Geometry and analytical approach

Consider a nonhomogeneous parallelogram plate depicted in Fig. 1 with dimensions a and b, thickness l, and density ρ.

The skew plate, represented in Fig. 1, has a circular thickness denoted as l in one dimension. Additionally, Poisson’s ratio ν is assumed to be circular in one dimension:

1
l=l01+β1-1-ξ2a2, νξ=νξ01-m1-1-ξ2a2,

where, l0 and νξ0 stand for the initial thickness and Poisson’s ratio of the plate at the origin. Furthermore, β (ranging from 0 to 1) and m (ranging from 0 to 1) represent the taper and non-homogeneity parameters, respectively. We assume a two-dimensional steady-state temperature distribution on the plate, following the parabolic model introduced in [7]:

2
τ=τ01-ξa1-ψb,

where, τ and τ0 denote the temperature excess above the reference temperature at any point on the plate and at the origin, respectively. The temperature dependence modulus of elasticity for engineering structures follows [7]:

3
Eξτ=E11-γτ, Eψτ=E21-γτ, Gξψ(τ)=G0(1-γτ),

where, Eξ and Eψ denote Young’s moduli in the ξ and ψ directions, while Gξψ represents the shear modulus. The parameter γ accounts for the slope variation of moduli with temperature. By substituting Eq. (2) in Eq. (3), we obtain the following expressions:

4
Eξτ=E11-α1-ξa1-ψb, Eψ(τ)=E21-α1-ξa1-ψb,
Gξψ(τ)=G01-α1-ξa1-ψb,

where α=γτ0, 0α<1 is called temperature gradient.

Fig. 1Parallelogram plate having 1-D circular thickness

Parallelogram plate having 1-D circular thickness

The flexural rigidities Dξ, Dψ and torsional rigidity Dξψ of the plate are taken as in [12]:

5
Dξ=Eξl3121-νξνψ, Dψ=Eηl3121-νξνψ, Dξψ=Gξψl312, D1=νξDψ=νψDξ,

where νξ and νη are Poisson’s ratios. Using Eqs. (1), (3) and (4) in Eq. (5), we get:

6
Dξ=E1h03121-νξ01-mΥ1vψ1-α1-ξa1-ψb1+β1-1-ξ2a23,
7
Dψ=E2h03121-νξ01-mΥ1vψ1-α1-ξa1-ψb1+β1-1-ξa3,
8
Dξψ=G0h03121-α1-ξa1-ψb1+β1-1-ξ2a23,
9
D1=E1h03νψ121-νξ01-mΥ1vψ1-α1-ξa1-ψb1+β1-1-ξ2a23,

Now, introducing non-dimensional variable as:

10
E1*=E11-νξ01-mΥ1νψ, E2*=E21-νξ01-mΥ1νψ, E*=νξ01-mΥ1E2*=νψE1*.

The equation for kinetic energy Ts and strain energy Vs for natural transverse vibration of non-uniform orthotropic parallelogram is taken as in [1]:

11
Ts=12ω20a0bρl Φ2cosθdξdψ,
12
Vs=120a0bDξ2Φξ22+Dψ2Φξ2tan2θ-22Φξψtanθsecθ+2Φψ2sec2θ2
+2D12Φξ22Φξ2tan2θ+22Φξψtanθsecθ+2Φψ2sec2θ
+4Dξψ-2Φξ2tanθ+2Φξψsecθ2cosθdξdψ.

Rayleigh-Ritz method requires that maximum strain energy must be equal to maximum kinetic energy i.e.,:

13
J=δVs-Ts=0.

From Eqs. (1), (6-9), we have:

14
J=0a0b1-α1-ξa1-ψ b1+βΥ13cos4θ+E2E1sin4θ+2νψsin2θcos2θ2Φξ22+4G0E11-νξ01-mΥ1νψsin2θcos2θ2Φξ22+E2E12Φψ22+4E2E1sin2θ+G0E11-νξ01-mΥ1νψcos2θ2Φξψ2+2E2E1sin2θ+νψcos2θ2Φξ22Φψ2-4E2E1sin3θ+2νψsinθcos2θ2Φξ22Φξψ+2G0E11-νξ01-mΥ1νψsinθcos2θ2Φξ22Φξψ
-4E2E1sinθ2Φψ22Φξψdξdψ-λ20a0b1-mΥ11+βΥ1Φ2dξdψ,

where λ2=12ρ0a2cos5θE1*h02, Υ1=1-1-ξ2a2.

The two term deflection function which satisfy all the edge conditions can be taken as in [7] The two term deflection function which satisfy all the edge conditions can be taken as in [7]:

15
Φξ,ψ=ξaeψbf1-ξag1-ψbh×i=0NΩiξaψb1-ξa1-ψbi.

This expression results from two components: one defining boundary conditions (0 for free edge, 1 for simply supported, 2 for clamped edge), and the other representing mode frequencies with constants Ωi for i=0,1,2,,N. To minimize the functional in Eq. (11), the following condition is necessary:

16
JΩi=0, i=0,1,2,3...N.

After simplifying Eq. (13), we get a homogeneous system of equations in Ωi whose non zero solution gives equation of frequency as:

17
P-λ2Q=0,

where, P=pijN+1 and Q=qijN+1 are square matrices of order (n+1) with i=0,1,2...N and j=0,1,2...N. The time period is calculated using the expression:

18
K=2πλ,

where λ is a frequency obtained from Eq. (14).

3. Numerical results and discussion

This study investigated the impact of parameters such as tapering, thermal gradient, and non-homogeneity on the vibration time period of an orthotropic parallelogram plate. The plate had a fixed aspect ratio of a/b=1.5, a skew angle of θ= 30°, circular thickness, and Poisson’s ratio. The analysis considered specific edge conditions and linear temperature effects, with material parameters sourced from [5]: E2*/E1*=0.01, E*/E1*=0.3, G0/E1*=0.0333, E1*/ρ0= 3.0×105 and νξ0=0.345. The results are presented in Tables 1-3.

Table 1Time period of orthotropic parallelogram plate at SCCC edge condition corresponding to non-homogeneity parameter m

m
α = 0.2
β = 0.0
β = 0.2
β = 0.4
β = 0.6
β = 0.8
K2
K1
K2
K1
K2
K1
K2
K1
K2
K1
SCCC
0.0
0.44915
0.07578
0.44001
0.07564
0.43034
0.07543
0.42022
0.07515
0.40973
0.07480
0.2
0.44909
0.07577
0.43995
0.07563
0.43027
0.07542
0.42013
0.07515
0.40963
0.07479
0.4
0.44903
0.07577
0.43989
0.07563
0.43021
0.07542
0.42006
0.07514
0.40954
0.07479
0.6
0.44897
0.07577
0.43982
0.07563
0.43015
0.07542
0.42000
0.07514
0.40948
0.07479
0.8
0.44893
0.07577
0.43976
0.07563
0.43005
0.07542
0.41991
0.07514
0.40938
0.07478

Table 2Time period of orthotropic parallelogram plate at SCCC edge condition corresponding to thermal gradient α

α
m = 0.2
β = 0.0
β = 0.2
β = 0.4
β = 0.6
β = 0.8
K2
K1
K2
K1
K2
K1
K2
K1
K2
K1
SCCC
0.0
0.43238
0.07212
0.42396
0.07202
0.41501
0.07186
0.40564
0.07163
0.39594
0.07134
0.2
0.44909
0.07577
0.43995
0.07563
0.43027
0.07542
0.42013
0.07515
0.40963
0.07479
0.4
0.46807
0.08004
0.45808
0.07985
0.44752
0.07958
0.43643
0.07923
0.42500
0.07880
0.6
0.48996
0.08512
0.47894
0.08485
0.46725
0.08450
0.45503
0.08406
0.44240
0.08352
0.8
0.51573
0.09130
0.50335
0.09093
0.49022
0.09046
0.47655
0.08989
0.46247
0.08921

Table 3Time period of orthotropic parallelogram plate at SCCC edge condition corresponding to tapering parameter β

β
m=α=0.0
m=α=0.2
m=α=0.4
m=α=0.6
m=α=0.8
K2
K1
K2
K1
K2
K1
K2
K1
K2
K1
SCCC
0.0
0.43244
0.07213
0.44909
0.07577
0.46800
0.08004
0.48984
0.08511
0.51551
0.09129
0.2
0.42402
0.07202
0.43995
0.07563
0.45801
0.07985
0.47878
0.08485
0.50310
0.09092
0.4
0.41510
0.07186
0.43027
0.07542
0.44743
0.07958
0.46706
0.08450
0.48996
0.09045
0.6
0.40574
0.07163
0.42013
0.07515
0.43637
0.07923
0.45484
0.08405
0.47627
0.08988
0.8
0.39600
0.07134
0.40963
0.07479
0.42490
0.07880
0.44221
0.08352
0.46216
0.08920

Table 1 presents time period data for an orthotropic parallelogram plate under the SCCC edge condition and a constant thermal gradient (α=0.2). The table covers a range of non-homogeneity parameter (m) and tapering parameter (β) values from 0.0 to 0.8. The results show that higher β values lead to a decrease in time period K, similar to the effect of increasing m on reducing K.

Table 2 exhibits time period (K) data for an orthotropic parallelogram plate under the SCCC edge condition, with a constant non-homogeneity parameter (m=0.2). The table encompasses β and α values from 0.0 to 0.8. Remarkably, higher thermal gradient (α) values result in elevated time period (K), while increasing the tapering parameter (β) leads to a significant decrease in K.

Table 3 presents time periods for an orthotropic parallelogram plate under the SCCC edge condition, with variable non-homogeneity (m), thermal gradient (α), and tapering (β) parameters ranging from 0.0 to 0.8. Higher values of α and m lead to an increase in the time period (K), while raising β results in a reduction of K.

4. Conclusions

In terms of time periods, Table 1 indicates that β holds greater influence than m under the SCCC edge condition. Additionally, in Tables 2 and 3, β demonstrates a more substantial effect on the time period’s rate of change compared to α and m, respectively.

References

  • A.W. Leissa, “Vibration of plates scientific and technical information division,” National Aeronautics and Space Administration, 1969.
  • M. Di Sciuva, “Bending, vibration and buckling of simply supported thick multilayered orthotropic plates: An evaluation of a new displacement model,” Journal of Sound and Vibration, Vol. 105, No. 3, pp. 425–442, Mar. 1986, https://doi.org/10.1016/0022-460x(86)90169-0
  • A. K. Gupta and Mamta, “Effect of exponential temperature variation on frequencies of a rectangular plate of non-linear varying thickness,” American Journal of Applied Sciences, Vol. 10, No. 2, 2019.
  • A. K. Gupta, A. Kumar, and H. Kaur, “Vibration of visco-elastic orthotropic parallelogram plate with linear thickness variation in both directions,” International Journal of Acoustics and Vibration, Vol. 2, pp. 72–80, 2011.
  • A. K. Gupta, A. Kumar, and D. V. Gupta, “Vibration of visco-elastic orthotropic parallelogram plate with parabolically thickness variation,” Annals of the Faculty of Engineering Hunedoara, Vol. 10, No. 2, pp. 61–70, 2012.
  • A. Khanna and A. Singhal, “Effect of plate’s parameters on vibration of isotropic tapered rectangular plate with different boundary conditions,” Journal of Low Frequency Noise, Vibration and Active Control, Vol. 35, No. 2, pp. 139–151, 2016.
  • N. Lather, R. Bhardwaj, A. Sharma, K. Kumar, and A. L. Schiavo, “Time period analysis of orthotropic skew plate with 2-D circular thickness and 1-D circular density,” Mathematical Problems in Engineering, Vol. 2022, pp. 1–15, 2022.
  • A. Sharma, “Natural vibration of parallelogram plate with circular variation in density,” Acta Technica, Vol. 63, No. 6, pp. 763–774, 2018.
  • A. Sharma, A. Raghav, and V. Kumar, “Mathematical study of vibration on non-homogeneous parallelogram plate with thermal gradient,” International Journal of Mathematical Sciences, Vol. 9, No. 2, 2016.
  • A. K. Sharma and M. K. Dhiman, “Vibration analysis of simply supported parallelogram plate with bi-dimensional thickness and temperature deviation,” International Journal of Research in Advent Technology, Vol. 6, No. 12, 2018.
  • A. K. Sharma and V. Verma, “Inspection of vibration using Rayleigh-Ritz method for orthotropic non-homogeneous plate with varying thickness and thermal effect,” Romanian Journal of Acoustics and Vibration, Vol. 14, No. 1, 2017.
  • A. K. Sharma, M. K. Dhiman, H. Singh, P. Prashar, and P. Sharma, “Effect of non-homogeneity with thickness and temperature variation on vibration of orthotropic parallelogram plate with simply supported edges.,” ARPN Journal of Engineering and Applied Sciences, Vol. 14, No. 15, pp. 2757–2762, 2006.

About this article

Received
31 October 2023
Accepted
11 November 2023
Published
27 November 2023
SUBJECTS
Modal analysis and applications
Keywords
orthotropic parallelogram plate
thermal gradient
circular tapering
nonhomogeneity
Acknowledgements

The authors have not disclosed any funding.

Data Availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflict of interest

The authors declare that they have no conflict of interest.