Numerical experiments with stable versions of the Generalized Finite Element Method

• Autor(a) principal: Sato, Fernando Massami
• Data de Publicação: 2017
• Tipo de documento: Dissertação
• Idioma: eng
• Título da fonte: Biblioteca Digital de Teses e Dissertações da USP
• Texto Completo: http://www.teses.usp.br/teses/disponiveis/18/18134/tde-16102017-101710/
• topic: Generalized Finite Element Method; Higher Order Stable Generalized Finite Element Method; Método dos Elementos Finitos Generalizados; Método dos Elementos Finitos Generalizados Estável; Método dos Elementos Finitos Generalizados Estável de ordem superior; Número de condição escalonado; Rate of convergence; Razão de convergência; Scaled condition number; Stable Generalized Finite Element Method

The Generalized Finite Element Method (GFEM) is essentially a partition of unity based method (PUM) that explores the Partition of Unity (PoU) concept to match a set of functions chosen to efficiently approximate the solution locally. Despite its well-known advantages, the method may present some drawbacks. For instance, increasing the approximation space through enrichment functions may introduce linear dependences in the solving system of equations, as well as the appearance of blending elements. To address the drawbacks pointed out above, some improved versions of the GFEM were developed. The Stable GFEM (SGFEM) is a first version hereby considered in which the GFEM enrichment functions are modified. The Higher Order SGFEM proposes an additional modification for generating the shape functions attached to the enriched patch. This research aims to present and numerically test these new versions recently proposed for the GFEM. In addition to highlighting its main features, some aspects about the numerical integration when using the higher order SGFEM, in particular are also addressed. Hence, a splitting rule of the quadrilateral element area, guided by the PoU definition itself is described in detail. The examples chosen for the numerical experiments consist of 2-D panels that present favorable geometries to explore the advantages of each method. Essentially, singular functions with good properties to approximate the solution near corner points and polynomial functions for approximating smooth solutions are examined. Moreover, a comparison among the conventional FEM and the methods herein described is made taking into consideration the scaled condition number and rates of convergence of the relative errors on displacements. Finally, the numerical experiments show that the Higher Order SGFEM is the more robust and reliable among the versions of the GFEM tested.