Research Article  Open Access
Mohammad M. Elahi, Seyed M. Hashemi, "A Framework for Extension of Dynamic Finite Element Formulation to Flexural Vibration Analysis of Thin Plates", Shock and Vibration, vol. 2017, Article ID 5905417, 10 pages, 2017. https://doi.org/10.1155/2017/5905417
A Framework for Extension of Dynamic Finite Element Formulation to Flexural Vibration Analysis of Thin Plates
Abstract
Dynamic Finite Element formulation is a powerful technique that combines the accuracy of the exact analysis with wide applicability of the finite element method. The infinite dimensionality of the exact solution space of plate equation has been a major challenge for development of such elements for the dynamic analysis of flexible twodimensional structures. In this research, a framework for such extension based on subset solutions is proposed. An example element is then developed and implemented in MATLAB® software for numerical testing, verification, and validation purposes. Although the presented formulation is not exact, the element exhibits good convergence characteristics and can be further enriched using the proposed framework.
1. Introduction
Dynamic stiffness modeling is a wellestablished technique in vibrational analysis of structural elements. These methods seek to propose formulations that have the accuracy of the exact solutions and wider applicability of the finite element methods by incorporating some form of the closed form solutions of governing equations instead of polynomials used by classic finite element method (FEM).
Two of the most famous dynamic stiffness formulations, mainly applied to various beamstructures, are dynamic stiffness matrix (DSM) and Dynamic Finite Element (DFE) methods, which produce accurate results with much coarser mesh compared to FEM formulations.
Extensively developed in 1970s, and pioneered again by Banerjee and his coworkers since 1990s [1–8], the beam DSM formulation dates back to much earlier times [9]. The 1941 work by Kolousek [9] is probably the first to derive the dynamic stiffness matrix for the EulerBernoulli beam. The DSM formulation is performed by developing a dynamic relation between force and displacement matrix for a problem domain, where the displacement functions are combination of trigonometric and hyperbolic frequencydependent functions, which satisfy the EulerBernoulli or Timoshenko beam equations. In this method, the boundary and loading conditions of the problem in hand are taken into consideration in developing the forcedisplacement relationship and thus this process yields a onestep formulation to the problem but is very case specific [1–6]. Using truncated Taylor series expansion and assigning penalty for different boundary conditions, the DSM has been more recently generalized by Pagani and his coworkers [7, 8] to vibrational analysis of composite beams and solid and thinwalled structures subjected to various boundary conditions. This approach was also further extended to use higherorder 1D dynamic stiffness elements for the vibration analysis of composite plates [10].
On the other hand, the DFE technique follows the general procedure of FEM element development by using the general displacement solution of beam equation as basis functions. Next, the element dynamic shape functions are developed and supplied to a residual minimization scheme, such as Galerkin’s technique, to develop the frequencydependent element matrix. The matrix thus developed is not case specific and can be assembled. The boundary and loading conditions are applied in a manner similar to FEM [11–14].
Both DSM and DFE methods generate one matrix per problem which has the effects of mass and stiffness combined. In fact, the mass terms are distributed over the element displacement as a dynamic force using D’Alembert’s principle. The assembly of the element dynamic stiffness matrices and the application of the system boundary conditions then lead to a nonlinear eigenvalue problem, which is then solved using a root finding technique, such as WittrickWilliams algorithm [15].
Both the DSM and DFE formulations have been very successful for onedimensional elements such as rods, beams, and beamstructures under various loading and boundary conditions. The expansion of such techniques to twodimensional (2D) elements such as plates, however, has led to cumbersome mathematics with limited applicability which served as a motive for current research. The particular problems with such extension arise from the fact that the plate governing equations are twodimensional partial differential equations with infinite number of solutions. This is unlike beam and rod situation where the frequencydependent solution of the ordinary differential equation can be used directly as the basis functions of approximation space to develop set of dynamic shape functions [16].
The other problem in the case of DFE formulation comes from evaluating the resulting integral equations, written in terms of the dynamic shape functions, to develop element matrices. Since these functions are transcendentally dependent on frequency, the integrations are difficult to handle given the current computational power of computers especially for infinite dimensional domains. This difficulty is particularly important when extension of such elements to irregular shapes is required. In such situations, the area integrals used in mapping of the element to its natural coordinate system will involve the absolute value of Jacobian mapping matrix that is difficult to handle for such complicated displacement functions.
Because the Dynamic Finite Element matrices are a result of these integrations, use of numerical integration schemes such as Gauss quadrature formula will require introduction of numerous evaluation points for an acceptable result. This is unlike FEM extension to plates of arbitrary shapes where the polynomial shape functions provide descent result with small number of evaluation points. Besides, FEM shape functions are much easier to evaluate numerically as they do not have any dependency on vibrational frequency and produce constant matrices of numerical values that can be manipulated with ease. Therefore, much of the available literature for exact solution of plates was focused on solving these equations for simple geometries and loading cases, as presented by Leissa [17–22], or involved variation error minimization techniques, such as finite element method [23–34].
Recently, there have been more researches to extend the dynamic stiffness matrix scheme to plate vibrations. Since this method involves forcedisplacement relationship development to arrive at the element DSM, it requires developing differential relationships, avoiding the integration problem. Casimir et al. [16] developed the DSM for rectangular thin plates and demonstrated that accurate results can be achieved with limited number of elements. However, their method involved solving infinite dimensional matrices, limited to simple cases, and provided problem specific formulation that must be reformulated for each new configuration. This reformulation for plate elements would therefore involve solving large matrices and is less favourable to be used as a general analysis tool. More recently, Liu and Banerjee [35–37] presented a novel SpectralDSM (SDSM) method to generalize the DSM using truncation of Fourier series and the application of spectral analysis. The SDSM, when applied to the free vibration analysis of complex problems involving isotropic laminated plates under various boundary conditions, demonstrated significantly improved efficiency over the conventional FEM formulations.
In this paper, the focus is to present an alternate dynamic stiffness formulation in the form of DFE for the vibration analysis of Kirchhoff plates. Particularly, the focus is to generate a model that can be easily extended to plates of different shapes without loss of accuracy of the original formulation. The mathematical procedure of beam DFE formulation served as a guideline for development of the current plate DFE design, as presented in the following section. An example 4node 12DOF element is then developed and validated against FEM, where its accuracy and efficiency are demonstrated through the free flexural vibration analysis.
2. Mathematical Modeling
To develop Dynamic Finite Element of thin plates, we begin with the partial differential equation governing the free flexural vibration [38]:where is the flexural displacement, is the mass density, is the thickness, and represents bending coefficient. The terms and denote Young’s modulus and Poisson’s ratio of the plate material, respectively.
Applying Galerkin’s weighted residual minimization scheme [39], the integral form of the governing equation (1) is written aswhere is the residual weighting function and is the differential of the area within closed element boundary.
Green’s theorem is a conversion mechanism between area and line integrals performed over a closed boundary. For functions and Green’s theorem stateshere, and are components of the unit normal vector to the element boundary, and is differential segment along the boundary.
Since Galerkin’s method uses the same approximation for both field variable and weighting function , by applying Green’s theorem (4) to the above area integrals twice and taking the conformity of virtual displacements and loads across element boundaries into consideration, one can change Galerkin’s integral equation (3) to its equivalent weak form:where can also be interpreted as the virtual displacement. Equation (5) is the formulation generally used for finite element development of plates. The DFE derivation will deviate from FEM at this point.
Let be the dynamic flexural displacement, where is a row vector containing the dynamic shape function, derived from the frequencydependent (dynamic) basis functions of approximation space, and is the vector of nodal displacements. The dynamic basis functions, in this case, are defined to be the solutions of the Kirchhoff’s plate governing equation.
Applying Green’s theorem (4) to the above area integrals twice, (5) can be written as
Note that the last 4 area integrals can be combined to generate Kirchhoff’s plate equation in terms of virtual displacement and thus can be written as
Assuming, as mentioned earlier, the dynamic shape functions are derived from the dynamic basis function, which are in turn solutions of Kirchhoff’s equation, the terms in the bracket will yield zero. Therefore, the DFE stiffness matrix can be obtained using only evaluation of line integrals performed over element boundaries, as follows:
Although not explicitly visible, the above matrix is symmetric as it was obtained using mathematical manipulation of the symmetric Galerkin’s weak integral form.
2.1. Extension to Arbitrary Shapes
The absence of area integrals in DFE stiffness matrix evaluation enables exact integration of the boundary line integrals when transformed between different coordinate systems using variable substitution. Consider, for example, the quadrilateral case given in Figure 1(b).
(a)
(b)
The variables and from element coordinate system can be related to the variables and from natural coordinate system using bilinear relations:where the coefficients and can be obtained by evaluating these relations for the nodes of the element. For ,
The DFE stiffness matrix (8) defined in last section was developed for any arbitrary closed area in element domain. Using relationships between element and natural coordinate systems (9) developed above, differentials of DFE matrix can be replaced with appropriate natural coordinate corresponding terms. The first integral, for example, can be evaluated in the natural coordinate system as
In the above integrations, the shape function vector, , and its element coordinate system (Figure 1(b)) derivatives must be converted to the natural coordinate system (Figure 1(a)) as well before the integration is performed. This can be done by evaluating the derivatives in the element coordinate system and then replacing the variables and with the respective definitions in natural coordinate system. Alternatively, one can transfer the shape function to natural coordinate system and then define the relation between derivatives in natural coordinate system and element coordinate system using application of chain rule. It is important to note that the relationship between the natural and element coordinate system does not need to be bilinear and the geometry is not restricted to quadrilaterals. The same procedure is applicable to any geometry and any relationship.
3. Development of Example Element
The governing differential equation of Kirchhoff plates is a twodimensional homogenous partial differential equation. The solution to such equations forms an infinite dimensional Hilbert space, which is usually prescribed in terms of Fourier series for dynamic formulations of thin plate vibrations. Although the DFE mathematical formulation presented does not have any restriction on the dimension of the shape function vector used, to generate a widely applicable scheme, a subset of the infinite solution space is used to develop example element. In this approach, instead of using the general Fourier trigonometric functions, a solution space is generated considering the physics of the element equation rather than a pure mathematical approach. Since in the development of the DFE the entire solution space is not used, the solution will not be exact for one element but would rather converge to an analytical answer. Special focus will be given to solutions that have both the following characteristics:(1)Being symmetric with respect to and (2)Being continuous.
Referring to the differential equation of motion (1), governing a harmonically oscillating Kirchhoff (thin) plate, while there are multiple solutions possible to this equation, one simple subset is obtained by considering the following reduction:
Note that this simplification has made the plate equation similar to that of a thin beam in direction. For this simplification, it is necessary that the solution considered must make the summation equal to zero. One such solution set can be obtained when each of the terms of the summation is zero individually. These solutions can be found using separation of variables and solving the equation above as an ordinary differential equation. Assigning
two solution sets are thus possible:
The first set is identical to beam differential equation solutions in direction. Similar 8 other solutions can be developed for a reduction in direction, written as
Now, the FEM development can be followed to produce the shape functions using these solutions as basis function. It is worth noting that although the above 16 basis functions developed here are linearly independent in a continuous Hilbert space, when evaluated over some finite nodal points, they are found to produce linearly dependent values. This problem is resolved by mixing these solutions to create 12 new basis shape functions which in low frequency values (i.e., when the frequency tends to zero; static behavior) will converge to 12 Hermitian polynomial basis functions commonly used in the conventional FEM formulation of thin plates [38], that is, , , , , , , , .
The new dynamic basis functions are written as
The nonnodal displacement function can be defined as a linear combination of these shape functions, written as
By evaluating these displacements and corresponding derivatives at nodes of an element (Figure 1(a)), the following relations are obtained:
The constant matrix can be obtained aswhere for a 2 by 2 square element, matrix is obtained as
When premultiplied by the basis function matrix , the inverse of matrix generates the shape function vector , containing the 12 new dynamic (frequencydependent) shape/interpolation functions:
Note that although the shape functions obtained through this method are variables of the rotational vibration frequency on the element domain, they still maintain the characteristics of Hermitian shape functions at nodes in that they produce a value of 1 at their corresponding nodes and zero at other nodes. Also, these shape functions converge to their FEM counterparts for small values of vibrational frequency, when tends to zero. The expressions for dynamic shape functions are very lengthy and are omitted here for brevity. For demonstration purposes, the first resulting dynamic interpolation function is presented graphically in Figure 2, for two values of (approaching zero, Figure 2(a)) and (Figure 2(b)), alongside the corresponding first Hermitian shape function (Figure 2(c)). As can be seen from Figures 2(a) and 2(c), as the frequency, , approaches zero, the dynamic shape function, , becomes the same as the Hermitian shape function, [40], that is, .
(a)
(b)
(c)
Using the displacement and weighting functions (16) expressed in terms of these dynamic shape functions and introduced in the DFE boundary integral equation (8), the DFE frequencydependent element dynamic stiffness matrix, , is obtained.
Having this matrix, next an illustrative example of a square plate with three sides simply supported and one side free boundary condition is analyzed. The natural frequency is normalized over plate material and mechanical properties using the relation . The element length and width are equal to 2 units and Poisson’s ratio is taken as 0.3. A conventional FEM analysis using Hermitian shape functions [40] is also performed to highlight the improvements obtained from DFE formulation.
The element matrices are assembled, where applicable, to find the system’s global dynamic stiffness matrix . The system boundary conditions are then applied by modifying the DFE matrix, where the constrained degrees of freedom (DOF) are enforced by eliminating the corresponding rows and columns, in a similar way to the conventional FEM. This produces a nonlinear eigenvalue problem, written aswhere is the vector of system’s DOFs. The system natural frequencies are then evaluated by setting the determinant of DFE matrix equal to zero, . The determinant of is a transcendental function of natural frequency of vibration that must then be swept across the frequency domain to find the values for which determinant function is zero. For a plate with 3 sides simply supported and one side free (SSSSSSF) boundary condition, modeled using a single DFE, the determinant matrix is plotted in Figure 3.
The determinant function also produces undefined regions, known as poles [41], where the determinant value approaches infinity. These regions represent the denominator of the determinant function [12].
The plate’s first five nondimensionalized natural frequencies () were evaluated using FEM and DFE techniques and compared with the analytical data from [17] (see Table 1).

As demonstrated in Table 1, for a single element mesh, the performance of the presented DFE plate surpasses that of its FEM counterpart, especially for the fundamental frequency where DFE produces very small error. It is also worth noting that Hermitian FEM vibrational analysis produces a linear eigenvalue problem in form of a polynomial in terms of frequency. This polynomial, in turn, has finite number of roots (in the case of above boundary condition, two roots only) and therefore, detection of higher modes requires mesh refinement. The DFE element, however, leads to a nonlinear eigenproblem through which one can extract multiple frequencies, even when a singleelement mesh is used.
The 2 × 2 square plate is further analyzed using both FEM and DFE formulations and under three different boundary conditions, for which analytical results were reported by Leissa [17]. These boundary conditions include ClampedClampedFreeFree (CCFF) boundary condition, ClampedFreeSimply SupportedFree (CFSSF) boundary condition, and Simply SupportedFreeSimply SupportedFree (SSFSSF) boundary condition.
Since FEM can capture different number of frequencies for different boundary conditions, only the fundamental natural frequencies are reported here for consistency in comparison. The result of this analysis is reported in Table 2.

Besides the performance of a single DFE element, it is important to analyze the convergence characteristics of the developed DFE plate. The convergence of the DFE element was compared against the 12DOF plate finite element, which utilizes Hermitian polynomials. For this study, the mesh density was increased from a single element to 4 elements (2 by 2), and to 25 elements (5 by 5). For illustration purposes only, the result for the second and fifth mode of vibration of a plate with (SSSSSSF) boundary condition is compared in Figure 4 for DFE and FEM formulations as examples of monotonic and nonmonotonic convergence to the analytical solution reported in [17]. Similar convergence was observed for other modes as well.
(a)
(b)
As can be seen from Figure 4, the DFE plate has comparable convergence rate with FEM. However, as mentioned earlier, since DFE produces a transcendental determinant equation, it is capable of providing potentially infinite number of natural frequencies, while FEM plates create constant mass and stiffness matrices which can only generate a finite number of solutions through linear eigenvalue problem. Also, the plate DFE presented here was produced based on one subset of solution of the plate equation, which reduced the plate behavior to two beamlike equations. Many other solution sets can be found, for example, based on distribution of dynamic load over the element strain, and combined with these solutions to enrich element physics further. The technique can also be extended to generate plate elements with higher number of nodes and DOFs.
4. Conclusion
The wellestablished framework for development of beam Dynamic Finite Element (DFE) formulation, based on exact solutions of governing differential equation, is extended to twodimensional (2D) plate elements. Unlike other dynamic plate formulations proposed in past, which exploit an infinite dimensional solution, an element with limited number of solutions was developed. The solutions obtained were founded on physicsbased simplification of Kirchhoff’s equation to beamlike equations and demonstrated comparable performance with classic FEM formulations. The research is underway to further enrich the presented 2D element by using the proposed framework to potentially achieve a quasiexact plate Dynamic Finite Element (DFE) formulation and further extend DFE plates to arbitrary geometries of various distortions.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors would like to acknowledge the support provided by Ryerson University in carrying out this research.
References
 J. R. Banerjee, “Dynamic stiffness formulation for structural elements: a general approach,” Computers & Structures, vol. 63, no. 1, pp. 101–103, 1997. View at: Publisher Site  Google Scholar
 J. R. Banerjee, “Free vibration of sandwich beams using the dynamic stiffness method,” Computers & Structures, vol. 81, pp. 1915–1922, 2003. View at: Google Scholar
 J. R. Banerjee and A. J. Sobey, “Dynamic stiffness formulation and free vibration analysis of a threelayered sandwich beam,” International Journal of Solids and Structures, vol. 42, no. 8, pp. 2181–2197, 2005. View at: Publisher Site  Google Scholar
 J. R. Banerjee and F. W. Williams, “Coupled bendingtorsional dynamic stiffness matrix of an axially loaded timoshenko beam element,” International Journal of Solids and Structures, vol. 31, no. 6, pp. 749–762, 1994. View at: Publisher Site  Google Scholar
 J. R. Banerjee, C. W. Cheung, R. Morishima, M. Perera, and J. Njuguna, “Free vibration of a threelayered sandwich beam using the dynamic stiffness method and experiment,” International Journal of Solids and Structures, vol. 44, no. 2223, pp. 7543–7563, 2007. View at: Publisher Site  Google Scholar
 J. R. Banerjee, H. Su, and C. Jayatunga, “A dynamic stiffness element for free vibration analysis of composite beams and its application to aircraft wings,” Computers & Structures, vol. 86, no. 6, pp. 573–579, 2008. View at: Publisher Site  Google Scholar
 A. Pagani, E. Carrera, M. Boscolo, and J. R. Banerjee, “Refined dynamic stiffness elements applied to free vibration analysis of generally laminated composite beams with arbitrary boundary conditions,” Composite Structures, vol. 110, no. 1, pp. 305–316, 2014. View at: Publisher Site  Google Scholar
 A. Pagani, M. Boscolo, J. R. Banerjee, and E. Carrera, “Exact dynamic stiffness elements based on onedimensional higherorder theories for free vibration analysis of solid and thinwalled structures,” Journal of Sound and Vibration, vol. 332, no. 23, pp. 6104–6127, 2013. View at: Publisher Site  Google Scholar
 V. Kolousek, “Anwendung des gesetzes der virtuellen verschiebungen und des reziprozitatssatzes in der Stabwerksdynamic,” Ingenieur Archiv, vol. 12, pp. 363–370, 1941. View at: Google Scholar
 A. Pagani, E. Carrera, J. R. Banerjee, P. H. Cabral, G. Caprio, and A. Prado, “Free vibration analysis of composite plates by higherorder 1D dynamic stiffness elements and experiments,” Composite Structures, vol. 118, pp. 654–663, 2014. View at: Publisher Site  Google Scholar
 S. M. Hashemi and E. J. Adique, “A quasiexact dynamic finite element for free vibration analysis of sandwich beams,” Applied Composite Materials, vol. 17, no. 2, pp. 259–269, 2010. View at: Publisher Site  Google Scholar
 S. M. Hashemi and M. J. Richard, “A Dynamic Finite Element (DFE) method for free vibrations of bendingtorsion coupled beams,” Aerospace Science and Technology, vol. 4, no. 1, pp. 41–55, 2000. View at: Publisher Site  Google Scholar
 S. M. Hashemi and M. J. Richard, “Free vibrational analysis of axially loaded bendingtorsion coupled beams: a dynamic finite element,” Computers & Structures, vol. 77, no. 6, pp. 711–724, 2000. View at: Publisher Site  Google Scholar
 S. M. Hashemi and A. Roach, “A Dynamic Finite Element for the Free Vibration Analysis of ExtensionTorsion Coupled Composite Beams,” vol. 1, pp. 221–23, 2010. View at: Google Scholar
 W. H. Wittrick and F. W. Williams, “A general algorithm for computing natural frequencies of elastic structures,” The Quarterly Journal of Mechanics and Applied Mathematics, vol. 24, no. 3, pp. 263–284, 1971. View at: Publisher Site  Google Scholar  MathSciNet
 J. B. Casimir, S. Kevorkian, and T. Vinh, “The dynamic stiffness matrix of twodimensional elements: application to Kirchhoff's plate continuous elements,” Journal of Sound and Vibration, vol. 287, no. 3, pp. 571–589, 2005. View at: Publisher Site  Google Scholar
 W. Leissa, “Vibration of Plates," National Aeronautics and Space Administration (NASASP160), Washington, D.C., 1969”. View at: Google Scholar
 A. Leissa, “Literature review : plate vibration research, 1976  1980: complicating effects,” Shock and Vibration, vol. 13, no. 10, pp. 17–36, 1981. View at: Publisher Site  Google Scholar
 A. W. Leissa, “Literature review: plate vibration research, 1976–1980: classical theory,” Shock and Vibration, vol. 13, no. 9, pp. 11–22, 1981. View at: Publisher Site  Google Scholar
 A. W. Leissa, “Recent studies in plate vibration 1981–1985 part II: complicating effects,” The Shock and Vibration Digest, vol. 19, no. 3, pp. 10–24, 1987. View at: Google Scholar
 A. Leissa, “Recent research in plate vibrations: classical theory,” Shock and Vibration, vol. 9, no. 10, pp. 13–24, 1977. View at: Publisher Site  Google Scholar
 W. Leissa, “Recent research in plate vibrations: complicating effects,” The Shock and Vibration Digest, vol. 9, no. 11, pp. 21–35, 1977. View at: Google Scholar
 R. D. Henshell, D. Walters, and G. B. Warburton, “A new family of curvilinear plate bending elements for vibration and stability,” Journal of Sound and Vibration, vol. 20, no. 3, pp. 381–397, 1972. View at: Publisher Site  Google Scholar
 J. N. Reddy, “Free vibration of antisymmetric, angleply laminated plates including transverse shear deformation by the finite element method,” Journal of Sound and Vibration, vol. 66, no. 4, pp. 565–576, 1979. View at: Publisher Site  Google Scholar
 J. N. Reddy, “A Review of the Literature on Finite Element Modelling of Laminated Composite Plates and Shells,” The Shock and Vibration Digest, vol. 17, no. 4, pp. 3–8, 1985. View at: Google Scholar
 C. W. Bert, J. N. Reddy, W. C. Chao, and V. S. Reddy, “Vibration of thick rectangular plates of bimodulus composite material,” Journal of Applied Mechanics, vol. 48, no. 2, pp. 371–376, 1981. View at: Publisher Site  Google Scholar
 J. N. Reddy, “Large amplitude flexural vibration of layered composite plates with cutouts,” Journal of Sound and Vibration, vol. 83, no. 1, pp. 1–10, 1982. View at: Publisher Site  Google Scholar
 J. N. Reddy and T. Kuppusamy, “Natural vibrations of laminated anisotropic plates,” Journal of Sound and Vibration, vol. 94, no. 1, pp. 63–69, 1984. View at: Publisher Site  Google Scholar
 Y. V. K. Sadasiva Rao and G. Singh, “Vibration of corner supported thick composite plates,” Journal of Sound and Vibration, vol. 111, no. 3, pp. 510–514, 1986. View at: Publisher Site  Google Scholar
 C.N. Chang and F.K. Chiang, “Vibration analysis of a thick plate with an interior cutout by a finite element method,” Journal of Sound and Vibration, vol. 125, no. 3, pp. 477–486, 1988. View at: Publisher Site  Google Scholar
 L.C. Shiau and J.T. Chang, “Transverse shear effect on vibration of laminated plate using higherorder plate element,” Computers & Structures, vol. 39, no. 6, pp. 735–740, 1991. View at: Publisher Site  Google Scholar
 M. Ganapathi, T. K. Varadan, and B. S. Sarma, “Nonlinear flexural vibrations of laminated orthotropic plates,” Computers & Structures, vol. 39, no. 6, pp. 685–688, 1991. View at: Publisher Site  Google Scholar
 Y. K. Cheung and W. L. Kwok, “Dynamic analysis of circular and sector thick, layered plates,” Journal of Sound and Vibration, vol. 42, no. 2, pp. 147–158, 1975. View at: Publisher Site  Google Scholar
 C. LienWen and C. C. Chen, “Asymmetric vibration and dynamic stability of bimodulus thick annular plates,” Computers & Structures, vol. 31, no. 6, pp. 1013–1022, 1989. View at: Publisher Site  Google Scholar
 X. Liu and J. Banerjee, “The SpectralDynamic Stiffness Method: A Novel Approach for Exact Free Vibration Analysis of PlateLike Structures,” in Proceedings of the The Fifteenth International Conference on Civil, Structural and Environmental Engineering, Prague, Czech Republic, September 2015. View at: Publisher Site  Google Scholar
 X. Liu and J. R. Banerjee, “An exact spectraldynamic stiffness method for free flexural vibration analysis of orthotropic composite plate assemblies  Part I: Theory,” Composite Structures, vol. 132, pp. 1274–1287, 2015. View at: Publisher Site  Google Scholar
 X. Liu and J. R. Banerjee, “An exact spectraldynamic stiffness method for free flexural vibration analysis of orthotropic composite plate assemblies  Part II: Applications,” Composite Structures, vol. 132, pp. 1288–1302, 2015. View at: Publisher Site  Google Scholar
 P. Ragesh, V. Mustafa, and T. Somasundaran, “An integrated kirchhoff element by galerkin method for free vibration analysis of plates on elastic foundation,” Procedia Technology, vol. 24, pp. 148–154, 2016. View at: Publisher Site  Google Scholar
 L. Logan, A First Course in the Finite Element Method, Luton, UK, 2007.
 R. J. Melosh, “Structural Analysis of Solids,” ASCE Structural Journal, vol. 89, no. 4, pp. 205–248, 1963. View at: Google Scholar
 N. H. Erdelyi and S. M. Hashemi, “On the finite element free vibration analysis of delaminated layered beams: a new assembly technique,” Shock and Vibration, vol. 2016, Article ID 3707658, 2016. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Mohammad M. Elahi and Seyed M. Hashemi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.