Geometrically nonlinear polygonal finite element analysis of functionally graded porous plates. (December 2018)

• Record Type: Journal Article
• Authors: Nguyen, Nam V.; Nguyen, Hoang X.; Lee, Seunghye; Nguyen-Xuan, H.
• Is Part Of: Advances in engineering software. Volume 126(2018)
• Page Start: 110
• Page End: 126
• Publication Date: 2018-12
• Subjects: Polygonal finite element method -- Quadratic serendipity shape functions -- Functionally graded materials -- Porosity -- Timoshenko's beam -- Nonlinear dynamics; Computer-aided engineering -- Periodicals; Engineering -- Computer programs -- Periodicals; Engineering -- Software -- Periodicals; Periodicals; 620.0028553
• DOI: 10.1016/j.advengsoft.2018.11.005
• Languages: English
• ISSNs: 0965-9978

Abstract

Highlights: An improved polygonal finite element formulation for functionally graded porous plate analysis is presented. Quadratic serendipity shape functions over polygonal elements are exploited. Shear locking issue is overcome by the Timoshenko's beam element. Both linear and nonlinear anlyses of FG porous plates are shown. Excellent performance of the present approach is found through numerical results. Abstract: In this study, an efficient polygonal finite element method (PFEM) in combination with quadratic serendipity shape functions is proposed to study nonlinear static and dynamic responses of functionally graded (FG) plates with porosities. Two different porosity types including even and uneven distributions through the plate thickness are considered. The quadratic serendipity shape functions over arbitrary polygonal elements including triangular and quadrilateral ones, which are constructed based on a pairwise product of linear shape functions, are employed to interpolate the bending strains. Meanwhile, the shear strains are defined according to the Wachspress coordinates. By using the Timoshenko's beam to interpolate the assumption of the strain field along the edges of polygonal element, the shear locking phenomenon can be naturally eliminated. Furthermore, the C 0 –type higher-order shear deformation theory ( C 0 –HSDT), in which two additional variables are included in the displacement field, significantly improves the accuracy of numerical results. The nonlinear equations of static and dynamic problems are solved by Newton–Raphson iterative procedure and by Newmark's integration scheme in association with the Picard methods, respectively. Through various numerical examples in which complex geometries and different boundary conditions are involved, the proposed approach yields more stable and accurate results than those generated using other existing approaches.