Analysis of the Behavior of a Square Plate in Free Vibration by FEM in Ansys

In the realization of mechanical structures, achieving stability and balance is a problem commonly encountered by engineers in the field of civil engineering, mechanics, aeronautics, biomechanics and many others. The study of plate behavior is a very sensitive subject because it is part of the structural elements. The study of the dynamic behavior of free vibration structures is done by modal analysis in order to calculate natural frequencies and modal deformations. In this paper, we present the modal analysis of a thin rectangular plate simply supported. The analytical solution of the differential equation is obtained by applying the method of separating the variables. We are talking about the exact solution of the problem to the limit values. However, numerical methods such as the finite element method allow us to approximate these functions with greater accuracy. It is one of the most powerful computational methods for predicting dynamic response in a complex structure subject to arbitrary boundary conditions. The results obtained by MEF through Ansys 15.0 are then compared with those obtained by the analytical method.


Introduction
Structures and buildings are generally subject to increasingly complex excitations. It appears essential to characterize them and then to control their vibratory behavior in order to preserve them against fatigue and rupture [1]. Structural failure is generally caused by structural fatigue due to acoustic radiation and amplitude of vibration amplitudes in response to a deterministic or random forced excitation. The increasing use of the plate element in the various industrial sectors, shows the imperative need of the study of their vibrational behavior which thus becomes of great importance in order to help the engineers to design better structures. Vibration analysis is an important current topic, both from an academic and an industrial point of view. The latter affects many areas, such as space technology, naval and civil engineering, automotive, aeronautics, and bridges, buildings, or nuclear engineering. To ensure the desired and desirable mechanical performance in favor of the supporting structure, the determination of the dynamic characteristics of the plates is essential. Rectangular plates are widely used in various engineering disciplines and, from a technical point of view, it becomes necessary to know the natural frequencies of such structures [2].
As the excellent review articles by Leissa [3] [4], Liew, Xiang and Kitipornchai [5] show, vibration analysis is not a recent science. And yet during these decades, it will know a renewed interest for the needs of design and dimensioning of structures responding to the best operating conditions, safety, economy, aesthetics and sound insulation. In addition, it is to optimize, lighten the structures commonly used and subject to significant levels of excitement. For over a hundred years, plaque behavior has been the subject of exhaustive research. It is clear that this research based on the classical Kirchhoff hypothesis neglects the effect of shear deformation and rotational inertia, unlike the Risner-Mindlin plates which take into account these effects.
Born out of the need to solve problems of elasticity and structural analysis involving complex fields of civil engineering and aeronautics, the finite element method was born between the 1960s and the 1970s by the work of the researchers Argyris, Clough, and Zienkiewicz [6] [7]. As a result, several MEFs will be developed for dynamic analysis of structures such as thin and thick plates. Iso parametric elements are the most widely used in plate dynamics [6] [8]. The finite element method is certainly the most favorable because it is one of the most powerful calculation methods for predicting the dynamic response in a complex structure with arbitrary boundary conditions [9].
In this modest work, we are therefore involved in this vast and important field, whose objective is to determine the eigenfrequencies as well as the modes of the vibrations of homogeneous thin isotropic plates in free dynamics by different methods (analytic, finite element method). To treat these vibratory problems, the general idea is to express the deflection of the plate w by a linear combination of the eigen modes. The whole problem then amounts to expressing the mn A coefficients of the proposed form functions. Either we express them experimentally by replacing in the theoretical expression w by the experimental deflection, or by trying to calculate them theoretically by proposing mathematical models, finally, calculate the matrices of mass and rigidity of the plate and solve the problem with the eigenvalues to obtain the natural frequencies of the plate.
To achieve this, the work done is structured in different sections. Section 2 gives a brief description of the kinematics of rectangular plates and governs the equation of motion. In section 3, the study of frequencies and modes of vibration is done by analytical resolution. Section 3 deals with the formulation of dynamic plate behavior by MEF on Ansys. In Section 4, the numerical results are compared to the analytical results for a rectangular plate simply supported on two opposite sides.

Description of the Problem
We consider a plate having a length a = 1 m, a width b = 1 m and a thickness h = 0.020 m. The properties of the structure illustrated in Figure 1 are given: Material: steel; Young's modulus E in a: 2 × 10 11 ; Density ρ in kg/m 3 : 7850;

Geometry of the Plate
The plate being simply supported (no displacement and no moment) on two edges, the conditions of fixity make it possible to write:

Fields of Displacement
In the classical Kirchhoff model, the normal remains straight and is perpendicular to the average surface after deformation. Kirchhoff's field of displacement is then written: the displacement according to oz or displacement of bending. On the other hand, the rotations are given by: Which allows to write:

Deformation
The state of deformation of a plate can be considered as the state of superposition of membrane deformations and flexural deformations. In small deformation, we know that: What causes these constraints:

Expression of Bending Moments
The bending moments per unit length in the plate Mx, My and Mxy are obtained by integrating the constraints on the thickness of the plate.

Definition of Efforts
Let Nx, Ny and Nxy be these efforts, Hooke allows to write: N S σ = . Proceeding by integration, we have:

Expression of Deformation and Kinetic Energy
The deformation energy of the plate is calculated by integration on the volume of the deformation energy density.
The kinetic energy T is calculated by integration of the volume element:

Movement Equation
The vertical equilibrium of the plate element (dx, dy) gives rise to the relation: We obtain the differential equation of Lagrange which is a partial derivative equation verified by the vertical displacement: Such vibrations are called free or natural transverse vibrations. As previously stated, the natural vibrations depend only on the properties of the material and the geometry of the plate, and are inherent properties of the elastic plate regardless of any charge. Thus, in the case of free or natural vibration (harmonic movement external load ( ) , p x y is equal to zero), the only transverse forces acting on the plate are the inertial forces due to mass ρ per unit area, the equation above becomes: Classical equation of plate theory, which for most technical applications is sufficient for the study of bending problems.

Analytical Study of the Rectangular Plate
The solution to this equation is obtained by looking for the transverse displacement ( ) , , w x y t such that: , , e j t w x y t X x Y y T t w x y ω = = Which is a solution of the form function.  By introducing these elements, Equation (14) becomes: We put: , and we consider unknown β Which results in: Of these two equations, it follows from the equation of complete displacement of the plate simply supported on two opposite edges:  ρ ω The resolution of this equation leads to a natural pulse mn ω of the mode m and n such that: Eigen frequency: The analysis of the first two modes gives us: , 11 603.02 rad s ω = The analysis of natural pulsations by the Rayleigh technique confirms the veracity of these results especially for the first vibratory mode.

Modeling of the Dynamic Behavior of Plates by FEM in Ansys
The finite element method is a widely used and powerful tool for the analysis of complex structures. It consists of a discretization of the element into a finite number of generally triangular or rectangular elements [10]. This resolution technique involves the calculation of the matrix of mass and rigidity of the whole structure: Equation (24) will be solved to give the eigenfrequencies and the eigen mode or modal deformations of the structure using the ANSYS software for a simply supported plate.
Let us always consider our rectangular plate whose characteristics are defined in section 1 and represented under Ansys in Figure 2. The plate is modeled on a volume of 2, e + 007 mm 3 and has a mass of 157 kg.
The dimensions of the plate illustrated in Figure 2 and Figure 3 are given in section 2. The simulations are made on the plate having a referential temperature of 22˚C and discretized in 2401 elements for 17,300 nodes as seen in Figure  3. The bar graph indicates the specific frequency of each calculated mode and whose values are given in Table 1 after modal analysis of the structure.
It is noted that the frequencies grow with vibratory modes as illustrated in Figure 4 and Figure 5. On the other hand the displacements of the plate as shown in Figure 6 do not follow this law, they oscillate around 6 mm.      As shown in Figure 6, the maximum displacement is observed when m = 3, n = 1 and is equal to 7.548 mm. These frequencies are geometrically translated by the following modes of Figure 7 (from (a) to (l)): Table 1 and Table 2 show the values of the eigenfrequencies obtained by two methods considered, analytical and finite element structure calculation, corresponding to the rectangular configuration. From these two tables and in accordance with Figure 8 and Figure 9, we use Equation (25)     The error between the two methods use ın thıs case study is minimized and converges to zero according to the error graph illustrated in Figure 9.

Conclusion
In this paper, we highlight the dynamic analysis of a plate simply supported on two opposite edges in free vibration. Two techniques have been deployed to approach the fundamental eigenfrequencies of the plate to be studied. It is the method of separable variables based on the modeling of the transversal displacement from the characteristic functions of vibrations of the thin plates and the analysis by the method of the finite elements. The responses obtained analytically and numerically from the nonlinear equations developed for the calculation of the eigenvalues of the plate show us that the frequencies increase with the modal deformations. The conditions of fixity and their location impose a great influence on the behavior of the plate structure in vibration. A comparison between the frequencies calculated and those obtained by Ansys for the rectangular plates reveals a convergence of the two calculation techniques.