Analysis of free vibration of nonhomogenous trapezoidal plate with 2D varying thickness and thermal effect
Ashish Kumar Sharma^{1} , Amit Sharma^{2} , A. K. Raghav^{3} , Vijay Kumar^{4}
^{1}Department of Mathematics, IEC University, Baddi, H.P., India
^{2, 3, 4}Department of Mathematics, Amity University Haryana, Gurgaon, India
^{1}Corresponding author
Journal of Measurements in Engineering, Vol. 4, Issue 4, 2016, p. 201208.
https://doi.org/10.21595/jme.2016.16378
Received 6 September 2015; received in revised form 28 September 2016; accepted 25 October 2016; published 31 December 2016
JVE Conferences
The present of goal investigation is to study the effect of 2D linearly varying temperature on the vibrations of nonhomogeneous trapezoidal plate whose thickness also varies linearly in both directions and density varies linearly in one direction. Boundary condition which is CSCS with two term deflection is taken into consideration. RayleighRitz method is used to find the solution of the problem. The effect of other plate parameters such as nonhomogeneity constant, taper constant and aspect ratios have also been studied. Results are calculated with accuracy and presented in tabular graphical form.
Keywords: thermal gradient, nonhomogeneity, vibration, trapezoidal plate, thickness.
1. Introduction
Plate theory has been applied to reduce vibration and noise in structures since the end of the 19th century where it began with the work of German physicist Chladni, who discovered various modes of free vibrations experimentally.
Since then, it has developed into an escalating and expansive field with a wide variety of theoretical and empirical techniques, dealing with increasingly complicated problems. In marine and aerospace engineering fields, where lightweight structural elements with orthotropic materials are of primary importance, orthotropic trapezoidal plate has extensive application. They are essentially thin plate structures which provide enhanced stiffness and stability characteristics. Hey also have widespread application in many other branches of modern civil, mechanical and structural engineering and, in these applications, they are regularly subjected to static and time varying loads. Hence, analysis of orthotropic trapezoidal plate under different conditions has always been an area of immense interest to researchers. Research work on the dynamic characteristics of orthotropic trapezoidal plate has a long history as different researchers have employed different methodologies to perform dynamic analysis of such structures.
Plates of variable thickness are used in nuclear reactor structures, naval structures, and aeronautical fields, electro mechanical transducers for the electronic telephones and mirrors and lenses in optical system. he practical importance of these plates has made vibration analysis essential for avoiding resonance excited by internal or external forces.
Also, in recent years, the development of solid propellant rocket motors and the increased use of sot filaments in aerospace structure and high speed runways, many practical solid structure interaction problems such as floor slabs of multistory buildings and buildings activities in cold regions have intensified the need for solution of various problems of plates and other structures supported on elastic media.
Vibrational effect is not only restricted in the field of engineering and science but, our daily life is also affected by it. There are many authors who deal in the study of vibration phenomena of plates by considering different kind of shapes, varying thickness, homogeneity or nonhomogeneity, with or without consideration of temperature effect.
Many authors have reported about the vibrational analysis of different types of plates. Some of them are mentioned here.
Zenkour [1] studied an exact solution for bending thin rectangular plates with uniform, linear, and quadratic thickness variations. Tomar et al. [2] attempted a problem of axisymmetric vibrations of an isotropic elastic nonhomogeneous circular plate of linearly varying thickness. Tomar et al. [3] studied features of the free vibration of an isotropic nonhomogeneous infinite plate of parabolically varying thickness. Srinivas et al. [4] did an exact analysis for vibration of simply supported homogeneous and laminated thick rectangular plates. Saliba [5] considered fully clamped symmetrical trapezoidal plates with free transverse vibration.
The finite element method of the vibration of trapezoidal plates is discussed by Orris and Petyt [6]. Qatu [7] studied the natural frequencies for cantilevered laminated composite right triangular and trapezoidal plates. Kitipornchai et al. [8] discussed the problem of vibrations of thick trapezoidal plates.
Recently, Gupta and Sharma [9] studied the thermal effect on free vibration of orthotropic thin trapezoidal plate of linearly varying thickness. The authors have so far not come across any paper dealing with bi linearly varying temperature, thickness and linearly varying density. In the present work, the effect of nonhomogeneity, taper constant, thermal gradient and aspect ratios has been studied. The frequencies for the first two modes of vibration are obtained for CSCS nonhomogeneous trapezoidal plate by RayleighRitz method and display in the form of graphs.
2. Geometry of the plate
The plate which is used for this paper is shown in Fig. 1. This is thin, symmetric and nonhomogeneous trapezoidal plate. At the left edge, maximum occurring thickness is ${h}_{0}$ and at the right, minimum occurring thickness is $\alpha {h}_{0}$.
Fig. 1. Variable varying thickness of trapezoidal plate
3. Methodology and equation of motion
Consider atrapezoidal plate with bilinearly varying thickness in both $x$ and $y$ directions taken as:
where ${h}_{0}$ is thickness along $\xi =x/a$ and $\eta =y/b$.
Nonhomogeneity of plate arises due to linear variation in density along $x$axis as:
Twodimensional temperature distribution 𝜏 along the length for trapezoidal plate under consideration is:
where $\tau $ denotes the temperature excess above the reference and ${\tau}_{0}$ denotes the temperature excess above the reference temperature.
Modulus of elasticity for elastic materials, is described as:
where ${E}_{0}$ is the value of Young’s modulus along the reference temperature i.e. at $\tau =0$ and $\gamma $ is the slope of the variation of $E$ and $\tau $.
On substituting value $\tau $ of from Eq. (3) into (4):
where $\alpha =\gamma {\tau}_{0}\left(0\le \alpha \le 1\right)$.
4. Equation of motion
The governing differential equation for kinetic energy $T$ and strain energy $V$ is given by:
Flexural rigidity $D\left(\xi \right)$ of the plate is:
Also:
is the plate’s flexural rigidity at the side $\xi =\frac{1}{2}$, A is the area of the plate $\rho $ and is the mass density per unit area of the plate.
Using Eqs. (9) and (5) in (8) flexural rigidity is given by:
$\times \left(1\alpha \left(\frac{1}{2}\xi \right)\left(\frac{1}{2}\eta \right)\right).$
Using Eq. (10) in (7) and (1), (2) in (6):
$\times \left(1\alpha \left(\frac{1}{2}\xi \right)\left(\frac{1}{2}\eta \right)\right)\left\{{\left(\frac{1}{{a}^{2}}\frac{{\partial}^{2}w}{{\partial \xi}^{2}}+\frac{1}{{b}^{2}}\frac{{\partial}^{2}w}{{\partial \eta}^{2}}\right)}^{2}\right.$
$\left.2\left(1\nu \right)\left[\frac{1}{{a}^{2}}\frac{{\partial}^{2}w}{{\partial \xi}^{2}}\frac{{\partial}^{2}w}{{\partial \eta}^{2}}{\left(\frac{1}{ab}\frac{{\partial}^{2}w}{\partial \xi \partial \eta}\right)}^{2}\right]\right\}d\alpha .$
$\left.\times \left[1\left(1{\alpha}_{1}\right)\left(\xi +\frac{1}{2}\right)\right]\right]{\omega}^{2}d\alpha .$
5. Deflection function and boundary condition
A two term deflection function is taken as:
$+{A}_{2}{\left\{\left(\xi +\frac{1}{2}\right)\left(\xi \frac{1}{2}\right)\right\}}^{32+\left(12\right)},$
where ${A}_{1}$ and ${A}_{2}$ are constants to be evaluated. Eq. (13) satisfies boundary conditions and provides good frequency estimation. ClampedSimply supportedClampedSimply supported plate is taken into consideration as shown in Fig. 2.
Fig. 2. Boundary conditions for the plate
Also, considered plate having boundaries which is defined a:
6. Method of solution
For finding the solution, RayleighRitz technique is used. According to its maximum strain energy must be equal to maximum kinetic energy, so for the problem under consideration we have:
Using Eq. (14) in (11) and (12):
$\times \left(1\alpha \left(\frac{1}{2}\xi \right)\left(\frac{1}{2}\eta \right)\right)\left\{{\left(\frac{1}{{a}^{2}}\frac{{\partial}^{2}w}{{\partial \xi}^{2}}+\frac{1}{{b}^{2}}\frac{{\partial}^{2}w}{{\partial \eta}^{2}}\right)}^{2}\right.$
$\left.2\left(1\nu \right)\left[\frac{1}{{a}^{2}}\frac{{\partial}^{2}w}{{\partial \xi}^{2}}\frac{{\partial}^{2}w}{{\partial \eta}^{2}}{\left(\frac{1}{ab}\frac{{\partial}^{2}w}{\partial \xi \partial \eta}\right)}^{2}\right]\right\}d\eta d\xi ,$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\times \left.\left[1\left(1{\alpha}_{1}\right)\left(\xi +\frac{1}{2}\right)\right]\right]d\eta d\xi .$
Now Eq. (15) becomes:
where:
$\times \left(1\alpha \left(\frac{1}{2}\xi \right)\left(\frac{1}{2}\eta \right)\right){\left(\frac{1}{{a}^{2}}\frac{{\partial}^{2}w}{{\partial \xi}^{2}}+\frac{1}{{b}^{2}}\frac{{\partial}^{2}w}{{\partial \eta}^{2}}\right)}^{2}$
$\left.2\left(1\nu \right)\left[\frac{1}{{a}^{2}}\frac{{\partial}^{2}w}{{\partial \xi}^{2}}\frac{{\partial}^{2}w}{{\partial \eta}^{2}}{\left(\frac{1}{ab}\frac{{\partial}^{2}w}{\partial \xi \partial \eta}\right)}^{2}\right]\right\}d\eta d\xi ,$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\times \left.\left[1\left(1{\alpha}_{1}\right)\left(\xi +\frac{1}{2}\right)\right]\right]d\eta d\xi .$
and ${\lambda}^{2}=12{\rho}_{0}{a}^{4}{\omega}^{2}\left({1\nu}^{2}\right)/{{E}_{0}{h}_{0}}^{2}$ is a frequency parameter.
Two unknown ${A}_{1}$ and ${A}_{2}$ are involve in Eq. (18) because of substitution of $w$ from Eq (13). These two unknowns are to be evaluated from Eq. (18), as follows:
On simplifying Eq. (21), we get:
where ${b}_{m1}$, ${b}_{m2}\left(m=1,2\right)$ involve parametric constant and the frequency parameter.
For nonzero solution, the determinant of the coefficient of Eq. (22) must vanish. Therefore, the equation of frequency becomes:
After solving Eq. (23), we get quadratic equation in ${\lambda}^{2}$.
7. Result and discussion
Here numerical results for a nonhomogenous trapezoidal plate of bilinearly varying thickness are calculated. Frequencies are calculated for the first two modes of vibration for different values of taper constant and thermal gradient, for CSCS boundary condition.
Frequency Eq. (23) is quadratic in ${\lambda}^{2}$, so we have tworoots. First two vibrational modes of frequency are calculated for a CSCS, trapezoidal plate, for various values of $a/b$, $c/b$, ${\beta}_{1}$, ${\beta}_{2}$, $\alpha $ and $v=$ 0.33 along with linearly varying thickness. Graphical representations are used for the results.
In Fig. 3 it is observed that the frequency parameter $\lambda $ increases with increasing values of thermal gradient $\alpha $ for different values of taper constant (${\beta}_{1}={\beta}_{2}=$4.0, 0.0), nonhomogeneity constant (${\alpha}_{1}=$0.4, 1.0) and aspect ratios ($a/b=$ 1.0), ($c/b=$0.5) for both the modes of vibration. The rate of increase as well as the value of frequency parameter $\lambda $ is higher in second mode in comparison to the first mode.
In Fig. 4 it is seen that frequency parameter decreases with increasing values of taper constants (${\beta}_{1}$ and ${\beta}_{2}$), for different values of thermal gradient ($\alpha =$0.0, 0.4), taper constants (${\beta}_{1}$ and ${\beta}_{2}$), nonhomogeneity constant (${\alpha}_{1}=$0.4, 1.0) and aspect ratios ($a/b=$ 1.0), ($c/b=$0.5) for both the modes of vibration. The rate of decrease as well as the value of frequency parameter $\lambda $ is higher in second mode in comparison to the first mode.
Fig. 3. Frequency parameter ($\lambda $) for a trapezoidal plate for different values of taper constant (${\beta}_{1}={\beta}_{2}=$4.0, 0.0), thermal gradient ($\alpha $), nonhomogeneity constant (${\alpha}_{1}=$0.4, 1.0) and aspect ratios ($a/b=$ 1.0), ($c/b=$0.5)
Fig. 4. Frequency parameter ($\lambda $) for a trapezoidal plate for different values of thermal gradient ($\alpha =$0.0, 0.4), taper constant (${\beta}_{1}$ and ${\beta}_{2}$), nonhomogeneity constant (${\alpha}_{1}=$0.4, 1.0) and aspect ratios ($a/b=$ 1.0), ($c/b=$0.5)
From Fig. 5 and 6 it is noticed that the frequency parameter decreases with increasing values of aspect ratio $a/b$ and $c/b$ for different combinations of thermal gradient (${\beta}_{1}$ and ${\beta}_{2}$), taper constant ($\alpha $) and fixed value of nonhomogeneity constant (${\alpha}_{1}=$ 0.0) and aspect ratio ($c/b=$ 0.75) for both the modes of vibration. The rate of decrease of frequency parameter $\lambda $ in Fig. 5 is higher in first and second mode in comparison
Fig. 7 clearly shows that frequency parameter $\lambda $ decreases with increasing value of nonhomogeneity constant ${\alpha}_{1}$ for different value of nonhomogeneity constant (${\alpha}_{1}$) with different combinations of thermal gradient (${\beta}_{1}$ and ${\beta}_{2}$), and taper constant ($\alpha $) and aspect ratios ($a/b=$ 1.0), ($c/b=$ 1.0) for both the mode of vibration.
Fig. 5. Frequency parameter ($\lambda $) for a trapezoidal plate for different combinations of thermal gradient (${\beta}_{1}$ and ${\beta}_{2}$), taper constant ($\alpha $) and fixed value of nonhomogeneity constant (${\alpha}_{1}=$ 0.0) and aspect ratio ($c/b=$0.75)
Fig. 6. Frequency parameter ($\lambda $) for a trapezoidal plate for different combinations of thermal gradient (${\beta}_{1}$ and ${\beta}_{2}$), taper constant ($\alpha $) and fixed value of nonhomogeneity constant (${\alpha}_{1}=$ 0.0) and aspect ratio ($a/b=$ 0.75)
Fig. 7. Frequency parameter ($\lambda $) for a trapezoidal plate for different value of nonhomogeneity constant (${\alpha}_{1}$) with different combinations of thermal gradient (${\beta}_{1}$ and ${\beta}_{2}$) and taper constant ($\alpha $) and aspect ratios ($a/b=$ 1.0), ($c/b=$1.0)
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. 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
 Zenkour A. M. An exact solution for the bending of thin rectangular plates with uniform, linear, and quadratic thickness variations. International Journal of Mechanical Sciences, Vol. 45, Issue 2, 2003, p. 295315. [Publisher]
 Tomar J. S., Gupta D. C., Jain N. C. Axisymmetric vibrations of an isotropic elastic nonhomogeneous circular plate of linearly varying thickness. Journal of Sound and Vibration, Vol. 85, Issue 3, 1982, p. 365370. [Publisher]
 Tomar J. S., Gupta D. C., Jain N. C. Free vibrations of anisotropic nonhomogeneous infinite plate of parabolically varying thickness. International Journal of Pure and Applied Mathematics, Vol. 15, Issue 2, 1984, p. 220221. [Search CrossRef]
 Srinivas S., Joga Rao C. V., Rao A. K. An exact analysis for vibration of simplysupported homogeneous and laminated thick rectangular plates. Journal of Sound and Vibration, Vol. 12, Issue 2, 1970, p. 187199. [Publisher]
 Saliba H. T. Transverse free vibration of fully clamped symmetrical trapezoidal plates. Journal of Sound and Vibration, Vol. 126, Issue 2, 1988, p. 237247. [Publisher]
 Orrisand R. M., Petyt M. A finite element study of the vibration of trapezoidal plates. Journal of Sound and Vibration, Vol. 27, Issue 3, 1973, p. 325344. [Publisher]
 Qatu M. S. Natural frequencies for cantilevered laminated composite right triangular and trapezoidal plates. Composites Science and Technology, Vol. 51, Issue 3, 1994, p. 441449. [Publisher]
 Kitipornchai S., Xiang Y., Liew K. M., Lim M. K. A global approach for vibration of thick trapezoidal plates. Computers and Structures, Vol. 53, Issue 1, 1994, p. 8392. [Publisher]
 Gupta A. K., Sharma S. Thermally induced vibration of orthotropic trapezoidal plate of linearly varying thickness. Journal of Vibration and Control, Vol. 17, Issue 10, 2011, p. 15911598. [Publisher]
 Khanna A., Sharma A. K. A computational prediction on vibration of square plate by varying thickness with bidimensional thermal effect. International Journal of Emerging in Engineering and Development, Vol. 2, Issue 3, 2012, p. 191196. [Search CrossRef]
 Khanna A., Sharma A. K. Effect of thermal gradient on vibration of viscoelastic plate with thickness variation. American Journal of Computational and Applied Mathematics, Vol. 2, Issue 1, 2012, p. 3436. [Publisher]
 Khanna A., Sharma A. K. Mechanical vibration of viscoelastic plate with thickness variation. International Journal of Applied Mathematical Research, Vol. 1, Issue 2, 2012, p. 150158. [Search CrossRef]
 Khanna A., Kaur N., Sharma A. K. Effect of varying Poisson ratio on thermally induced vibrations of nonhomogeneous rectangular plate. Indian Journal of Science and Technology, Vol. 5, Issue 9, 2012, p. 32633267. [Search CrossRef]
 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. 312317. [Search CrossRef]
 Sharma Subodh Kumar, Sharma Ashish Kumar Rayleigh‑Ritz method for analyzing free vibration of orthotropic rectangular plate with 2D thickness and temperature variation. Journal of Vibroengineering, Vol. 17, Issue 4, 2015, p. 19892000. [Search CrossRef]
 Sharma Subodh Kumar, Sharma Ashish Kumar Mathematically study on vibration of visco‑elastic parallelogram plate. Mathematical Models in Engineering, Vol. 1, Issue 1, 2015, p. 1219. [Search CrossRef]
Cited By
Journal of Vibroengineering
Amit Sharma, Ashish Kumar Sharma, Vijay Kumar

2018
