Investigation of Bending Vibrations in Circular Plates at Elasticity
Nimafar M1*, Ganji DD2, Azadi R2 and Akbari MR3
1Department of Mechanical Engineering, Central Tehran Branch, Islamic Azad University, Iran
2Department of Mechanical Engineering, Babol Noshirvani University of Technology, Iran
3Department of civil engineering, University of Tehran, Iran
Submission: September 12, 2017; Published: October 17, 2017
*Corresponding author: Nimafar M, Department of Mechanical Engineering, Central Tehran Branch, Islamic Azad University, Tehran, Iran, Tel: +981133334205, Email: m.nimafar@gmail.com
How to cite this article: Nimafar M, Ganji DD, Azadi R, Akbari MR. Investigation of Bending Vibrations in Circular Plates at Elasticity. Biostat Biometrics Open Acc J. 2017;3(4): 555616. DOI: 10.19080/BBOAJ.2017.03.555616.
Abstract
In this literature, the vibration of circular plates in the presence of dynamic loads has been assessed and analyzed by Elasticity theory in dynamic state and also bending moments and shearing forces instantaneously have been discussed completely. It is necessary to mention that the afore-mentioned procedure has been done by utilizing the method of separation of variables. Then, attempts have been made to obtain modal frequencies by presenting a general formula for circular plates which can be considered as a novelty in the field of vibrations. A couple differential equations for vibrational plates by the separation equation has been solved and consequently modal mass and shearing force have been studied completely. Moreover, modal frequencies and modal shapes has generally been achieved and depicted graphically. In accordance with the obtained displacement charts, it is clear to say that displacement in the first mode of vibration is more than the next ones and also the same result is obtained for bending moment, shearing and stress forces. Furthermore, it is noteworthy that inasmuch as the mentioned circular plate is axisymmetric, the charts of modal shapes have been acquired symmetrically.
Keywords: Bending in Circular Page; Circular Plates Vibration; Modal Shapes; Polar Coordinates
Introduction
Generally speaking, all unknowns of the problems of theory of elasticity in dynamic state are consisted of three components of displacement, six stress components and six strain components along with time and then the mentioned problems will be solved by creating three equilibrium equations, six compatibility equations of strains and finally six equations which are related to the structural behavior. But a wide variety of problems due to their geometric forms and types of loading on them can easily be changed into simpler forms and solved. In the discussion of 2-D problems, it is considered that if the plate thickness is small in comparison with the plate dimensions and also the loading of the problem is limited to the internal forces of the plate, the problem can be solved in the state of 2-D stress with acceptable approximation [1-3]. In this paper, problems of the theory of elasticity in dynamic state have been investigated that they, in fact, are diffusion of stress waves in elastic media. By using Dynamic knowledge in building engineering and earthquake engineering, elastic theory in dynamic state is applied.
Investigating dynamic distributed loading on structural plates and reaction of plate structures as vibration is of great importance for users in civil engineering [4,5]. Plate vibrations in the presence of dynamic forces and loads exerted on them is, in fact, the investigation of the maximum instantaneous strength of the structural elements which the elements are dependent on some parameters such as stresses, bending, shearing forces, and instantaneous displacement due to the vibration of the structures [6-8]. In order to investigate the vibration of the circular plates in regard to dynamic loading [9-12], it is necessary to have the equation governing on the vibrational system which can be non-linear, so in this case, a specific method (AGM) to solve these equations solved [13-15].
Nomenclature
Application
According to Figure 1, in order to gain the motion equations in terms of displacement function, it is essential to write Newton's first law for an optional volume element in a coordinate. Asregards Newton's second law for equilibrium of forces in accordance with Figure 1 which is a circular plate element, the vibrational differential equation governing on a circular plate with dynamic normal load in bending state is achieved as follows [16]:
In partial differential equation (1), parameter is defined as the unit mass per area and is expressed as and and is bending rigidity of the structure that is written as:
Therefore, Eq.(1) is rewritten as follows[13]:
The operator ∇2 in polar coordinates is defined as follows:
In this paper, as regards Figure 2 in the following, a circular plate is considered which its surrounding is fixed that the distributed load of vibration is exerted on it. With regard to Figure 2a & 2b, since loading in different points of the plate is applied homogeneously on the circular plate, the forces exerted on the plate are axisymmetric. In axisymmetric state, all the terms of differential Eq.(3) are independent from the angular . So, all the derivatives in relevance to are equaled to zero and as a result, in axisymmetric state, the governing differential Eq.(3) will be changed as follows:
Therefore in regard to Figure 2a & b q(r,t) is written as follows:
5.1. Solving the differential equation
Differential equation(5) is solved by the separation equation [14-16] as follows:
𝒰z = V(r).z(t) (7)
After substituting Eq.(7) into Eq.(5) and after arranging and substituting square of modal frequency, two differential equations have been yielded in the following forms:
And then, Eq.(9) can be rewritten as:
In Eq. (10), parameters ƒn(𝒯) and Mn are force function and modal mass of the plate vibration, respectively which are defined as equations bellow:
In Eq.(12),vn(r) is the function of modal shapes which by applying boundary conditions will be computed in the forgoing part of the paper.
Computing modal frequencies and modal shapes
By solving Eq. (8) with regard to the boundary conditions of clamping ends, modal frequencies and shapes can be calculated. The answer of differential equation(8) is achieved as the Bessel function of the first and second kinds (J and Y) from zero order as follows:
It is notable that in the above equation, parameter(i) is imaginary partof complex number and (λ) is defined as follows:
Applying boundary conditions
In regard to the circular plate with fixed surrounding and also with axisymmetric loading, boundary conditions can be defined as follows:
The plane slope at the circle center (r = 0) is zero because of the axisymmetric state in the circular plane:
Due to the characteristics of Bessel function from the second kind, (Y) , in (r = 0) we have:
Y0 = ∞ (16)
According to Eq.(16) and the answer of differential equation(13), (c2and c4) coefficients will be zero in order to keep the stability of the system:
C2 = 0 , C4 = 0 (17)
As a result, the answer of differential equation.(13) can be rewritten in accordance with Eq.(17) as follows:
Due to the fixed end of the circular plate, displacement and plane slope in points which joined to the body will be obtained zero as follows:
In order to compute non-zero answer from Eq.(19) and Eq.(20), the matrix determinant of the afore-mentioned set of equations should be zero as follows:
The answer of Eq.(21) is yielded by introducing new variable (βn = λ.R) as:
If Eq.(22) is depicted as a function in Cartesian coordinates, the roots of the equation can be observed in the following chart. (Figure 3)
Finally, the values of (βn) from Eq.(22) is acquired as follows:
On the basis of the obtained values from Table 1 and changing variable (βn = λR) , parameter (2) can be gained as follows:
By substituting Eq.(23) and variable changing(βn = λ.R) in to Eq.(14), modal frequencies can be achieved as follows (Figure 4-7):
Computing the equation of modal shapes
By computing the coefficient 𝒞3 in terms of the coefficient 𝒞3 from Eq.(19) or Eq.(20) and substituting into Eq.(18), the equation of modal shapes can be calculated as follows:
In regards to the values of βn from Table 1 or Eq.(23), the modal shapes can be depicted in different modes for the circular plate as follows:
Solving the differential equation of time index
By determining the function of modal shape 𝒰n(r) from Eq.(25), the functions of modal mass Mn and modal force ƒn (𝒯) from integrated equations(11,12) can easily be achieved. And also by determining modal mass and force, the differential equation of time index which is Eq.(10) can simply be computed. Table 2: The obtained modal frequencies in different modes of vibrations. It is notable that in order to compute the function of time index, Eq.(10), initial conditions (initial displacement and initial velocity) are needed which by supposing the lack of initial conditions, the answer of the time index differential equation can be yielded in the form of:
By determining the displacement of time index zn (t) from Eq.(26) and the modal shape function Vn(r) from Eq.(25), the vibrational displacement of structure element of the circular plate for every mode can be acquired as follows:
Afterwards the displacement of the structure can be gained for all of the modes according to the following equation:
Calculating bending moments, stresses and shearing forces in every mode of vibration
Bending moments for each mode of vibration (n) are achieved according to the following equations:
And shearing forces for each mode of vibration are computed in regard to equations bellow:
In this step, the stresses exerted on the plate for each mode of vibration are acquired in the forms of:
By choosing the following physical values for the structure, we will have:
On the basis of the afore-mentioned formulae and given physical values, the displacement of the circular plate, bending moment, shearing force and stress can be observed in the following charts for each mode of vibration. It is noteworthy that each time period, is divided into 4 time domains as follows (Table 2):
In regard to the afore-mentioned table in the first mode of vibration, the following charts are depicted on the basis of (ω, = 1.067 :
I. The displacement chart of the circular element Figure 8)
II. The charts of bending moment (Figure 9 & 10)
III. Shearing force charts (Figure 11)
IV. Stress charts (Figure 12 & 13)
Eventually, in the second mode of vibration which ω2 = 4.152 , the following charts are illustrated:
V. The displacement chart of the circular element (Figure 14)
It is notable that the other charts of the second mode and upper modes of vibration such as stress, shearing force and bending moment will easily be computed like the above procedure.
Conclusion
In the present paper, the vibration of circular plates in regard to the dynamic loading has been investigated by the method of separation of variables. In order to perform the afore-mentioned procedure, the Newton's first law for an arbitrary circular plate element and the Newton's second law for equilibrium of forces have been applied to achieve the vibrational differential equation governing on a circular plate with dynamic normal load in bending state. Since loading in different points of the mentioned plate is applied homogeneously, the forces exerted on the plate are axisymmetric and as everyone knows in axisymmetric state, all the terms of the obtained differential Equations are independent from the angular . So, all the derivatives in relevance to are equaled to zero and as a result, the equation governing on the present system simplified significantly.
After that by utilizing Separation law, modal frequencies and modal shapes have been acquired in regard to the boundary conditions. Consequently, the functions of modal mass and modal frequencies have been achieved. In this step, the time index function has been yielded in order to compute the vibrational displacement of the circular plate for each mode of vibration on the basis of the most general solution which is the sum of all solutions (Separation Law). Eventually, the values of bending moment, shearing force and shearing stress for the first and second mode of vibrations have been acquired and depicted completely. Results show that displacement in the first mode of vibration is more than the next ones and also the same result is obtained for bending moment, shearing and stress forces. Furthermore, it is noteworthy that in as much as the mentioned circular plate is axisymmetric, the charts of modal shapes have been acquired symmetrically.
References
- Akbari MR, Ganji DD, Nimafar M, Ahmadi AR (2014) Significant progress in solution of nonlinear equations at displacement of structure and heat transfer extended surface by new AGM approach. Frontiers of Mechanical Engineering 9(4): 390-401.
- Rostami AK, Akbari MR, Ganji DD, Heydari S (2014) Investigating Jeffery-Hamel flow with high magnetic field and nanoparticle by HPM and AGM. Central European Journal of Engineering 4(4): 357-370.
- Spencer PMaDE (1961) Field Theory for Engineers. Princeton Nj: Van Nostrand.
- YC Fung (1979) Foundations of Solid Mechanics: Englewood Cliffs NJ : Prentice - Hall, New York.
- M (1953) Methods of Theoretical Physics. New York: McGraw-Hill.
- Pekeris cl (1995) The seismic surface pulse. Proc Natl Acad Sci U S A 41(7): 469-480.
- Dukkipati RV (2009) Mechanical vibrations. New dehli, India: Narosa publishing house, India.
- Graff KF (1991) Wave Motion in Elastic Solids: Dover Publications.
- Dukkipati RV (2006) Advanced mechanical vibrations. Narosa publishing house, New Dehli, India.
- Dukkipati RV (2010) Problem solver in mechanical vibrations. New age international publishers(p)Ltd, New Dehli, India.
- Rand RH (2005) Lecture Notes on Nonlinear vibrations. Cornell University, New York, USA.
- Dettman JW (1962) Mathematical Methods in Physics and Engineering. New York: Mcgraw-Hill, USA.
- Akbari MR, Ganji DD, Rostami AK, Nimafar M (2015) Solving nonlinear differential equation governing on the rigid beams on viscoelastic foundation by AGM. Journal of Marine Science and Application 14(1): 30-38.
- Ganji D D AMR, Goltabar A R (2014) Dynamic Vibration Analysis for Non-linear Partial Differential Equation of the Beam - columns with Shear Deformation and Rotary Inertia by AGM. DAOE 3(1): 22-31.
- Akbari MR, Ganji DD, Ahmadi AR, Kachapi SHH (2014) Analyzing the nonlinear vibrational wave differential equation for the simplified model of Tower Cranes by Algebraic Method. Frontiers of Mechanical Engineering 9(1): 58-70.
- Seyed H, Hashemi Kachapi, Ganji DD (2011) Analysis of Nonlinear Equations in Fluids: Progress in Nonlinear Science.