Funções de interpolação e regras de integração tensorizaveis para o metodo de elementos finitos de alta ordem

Orientador: Marco Lucio Bittencourt Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica Resumo: Este trabalho tem por objetivo principal o desenvolvimento de funções de interpolaçao e regras de integraçao tensorizaveis para o Metodo dos Elementos Finitos (MEF) de alta ordem hp, considerando os sistemas de referencias locais dos elementos. Para isso, primeiramente, determinam-se ponderaçoes especficas para as bases de funçoes de triangulos e tetraedros, formada pelo produto tensorial de polinomios de Jacobi, de forma a se obter melhor esparsidade e condicionamento das matrizes de massa e rigidez dos elementos. Alem disso, procuram-se novas funçoes de base para tornar as matrizes de massa e rigidez mais esparsas possiveis. Em seguida, escolhe-se os pontos de integraçao que otimizam o custo do calculo dos coeficientes das matrizes de massa e rigidez usando as regras de quadratura de Gauss-Jacobi, Gauss-Radau-Jacobi e Gauss-Lobatto-Jacobi. Por fim, mostra-se a construçao de uma base unidimensional nodal que permite obter uma matriz de rigidez praticamente diagonal para problemas de Poisson unidimensionais. Discute-se ainda extensoes para elementos bi e tridimensionais Abstract: The main purpose of this work is the development of tensor-based interpolation functions and integration rules for the hp High-order Finite Element Method (FEM), considering the local reference systems of the elements. We first determine specific weights for the shape functions of triangles and tetrahedra, constructed by the tensorial product of Jacobi polynomials, aiming to obtain better sparsity and numerical conditioning for the mass and stiffness matrices of the elements. Moreover, new shape functions are proposed to obtain more sparse mass and stiffness matrices. After that, integration points are chosen that optimize the cost for the calculation of the coefficients of the mass and stiffness matrices using the rules of quadrature of Gauss-Jacobi, Gauss-Radau-Jacobi and Gauss-Lobatto-Jacobi. Finally, we construct an one-dimensional nodal shape function that obtains an almost diagonal stiffness matrix for the 1D Poisson problem. Extensions to two and three-dimensional elements are discussed. Doutorado Mecânica dos Sólidos e Projeto Mecânico Doutor em Engenharia Mecânica