Spelling suggestions: "subject:"scaled condition number"" "subject:"localed condition number""
1 |
Partições da Unidade flat-top e trigonométricas no Método dos Elementos Finitos Generalizados / Flat-top and trigonometric Partitions of Unity in the Generalized Finite Element MethodRamos, Caio Silva 11 April 2019 (has links)
Atualmente, no que concerne as problemáticas pertinentes à engenharia estrutural, o Método dos Elementos Finitos (MEF) é a principal ferramenta utilizada para obter soluções aproximadas de Problemas de Valor de Contorno (PVC). No entanto, tal metodologia exige um elevado custo computacional ao demandar malhas muito refinadas para solucionar problemas que apresentam singularidades, ou seja, que apresentam regiões onde ocorrem gradientes de deformação fortemente localizados. Para superar esse inconveniente, o Método dos Elementos Finitos Generalizados (MEFG) propõe a expansão do espaço de aproximação do MEF mediante a inserção de funções (conhecidas como funções de enriquecimento) que melhor representem localmente o comportamento da solução procurada. Tais funções podem apresentar características específicas ou mesmo serem geradas numericamente. Neste caso, dispensam-se malhas muito refinadas. Entretanto, o aumento do espaço de aproximação de modo irrestrito pode introduzir dependências lineares no sistema de equações do MEFG, tornando a solução obtida imprecisa ou mesmo impedindo a solução do sistema por métodos diretos. A chamada versão estável do MEFG explora uma modificação imposta às funções de enriquecimento a fim de melhorar o condicionamento da matriz de rigidez. Contudo, tal modificação não se configura como condição suficiente para garantir uma redução efetiva do número de condição. Neste trabalho, considera-se uma proposição recente para a modificação do espaço das funções de forma do MEFG associadas ao enriquecimento: trata-se do emprego de funções do tipo flat-top e trigonométricas como Partição da Unidade (PU), as quais são empregadas exclusivamente na construção das funções de forma enriquecidas (essas partições são definidas para elementos finitos quadrilaterais e triangulares). Exemplos numéricos são selecionados para evidenciar as vantagens dessas novas versões do MEFG em relação às anteriores e ao MEF convencional. Demonstra-se que tanto a PU flat-top quanto a PU trigonométrica, preservam as excelentes propriedades de convergência do MEFG. Além disso, mostra-se que o condicionamento da matriz de rigidez associada é próximo ao apresentado pelo MEF (uma vez que o enriquecimento, mesmo polinomial, não gera dependências) e que a formulação apresenta-se robusta na consideração de descontinuidades fortes. / Currently, regarding structural engineering issues, the Finite Element Method (FEM) is the main tool used to obtain approximate solutions of Boundary Value Problems (BVP). However, such methodology requires very refined meshes to solve problems that have singularities, i.e., that have regions where strongly localized deformation gradients occur, which leads to a high computational cost. To overcome this drawback, the Generalized Finite Element Method (GFEM) proposes the expansion of the FEM approach space by inserting functions (known as enrichment functions) that best represent locally the behavior of the searched solution. Such functions may have specific characteristics or even be generated numerically. In this case, very refined meshes are dispensed. However, the increase of the unrestricted approach space can introduce linear dependencies in the system of equations of the GFEM, making the solution imprecise or even preventing the solution of the system by direct methods. The so-called stable version of the GFEM exploits a modification imposed on the enrichment functions in order to improve the conditioning of the stiffness matrix. However, such a modification is not a sufficient condition to ensure an effective reduction in the condition number. In this work, it is considered a recent proposition to modify the space of the shape functions of GFEM associated with enrichment: the use of flat-top and trigonometric functions such as Partition of Unity (PU), which are used exclusively in the construction of the enriched shape functions (these partitions are defined for finite elements quadrilateral and triangular). Numerical examples are selected to highlight the advantages of these new versions of the GFEM over the previous ones and the conventional FEM. It is demonstrated that both flat-top PU and trigonometric PU preserve the excellent convergence properties of GFEM. In addition, it is shown that the conditioning of the associated stiffness matrix is close to that presented by FEM (since enrichment, even polynomial, does not generate dependencies) and that the formulation is robust in the consideration of strong discontinuities.
|
2 |
Numerical experiments with stable versions of the Generalized Finite Element Method / Experimentos numéricos com versões estáveis do Método dos Elementos Finitos GeneralizadosSato, Fernando Massami 21 August 2017 (has links)
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. / O Método dos Elementos Finitos Generalizados (MEFG) é essencialmente baseado no método da partição da unidade, que explora o conceito de partição da unidade para compatibilizar um conjunto de funções escolhidas para localmente aproximar de forma eficiente a solução. Apesar de suas vantagens bem conhecidas, o método pode apresentar algumas desvantagens. Por exemplo, o aumento do espaço de aproximação por meio das funções de enriquecimento pode introduzir dependências lineares no sistema de equações resolvente, assim como o aparecimento de elementos de mistura. Para contornar as desvantagens apontadas acima, algumas versões aprimoradas do MEFG foram desenvolvidas. O MEFG Estável é uma primeira versão aqui considerada na qual as funções de enriquecimento do MEFG são modificadas. O MEFG Estável de ordem superior propõe uma modificação adicional para a geração das funções de forma atreladas ao espaço enriquecido. Esta pesquisa visa apresentar e testar numericamente essas novas versões do MEFG recentemente propostas. Além de destacar suas principais características, alguns aspectos sobre a integração numérica quando usado o MEFG Estável de ordem superior, em particular, são também abordados. Por exemplo, detalha-se uma regra de divisão da área do elemento quadrilateral, guiada pela própria definição de sua partição da unidade. Os exemplos escolhidos para os experimentos numéricos consistem em chapas com geometrias favoráveis para explorar as vantagens de cada método. Essencialmente, examinam-se funções singulares com boas propriedades de aproximar a solução nas vizinhanças de vértices de cantos, bem como funções polinomiais para aproximar soluções suaves. Ademais, uma comparação entre o MEF convencional e os métodos aqui descritos é feita levando-se em consideração o número de condição do sistema escalonado e as razões de convergência do erro relativo em deslocamento. Finalmente, os experimentos numéricos mostram que o MEFG Estável de ordem superior é a mais robusta e confiável entre as versões do MEFG testadas.
|
3 |
Numerical experiments with stable versions of the Generalized Finite Element Method / Experimentos numéricos com versões estáveis do Método dos Elementos Finitos GeneralizadosFernando Massami Sato 21 August 2017 (has links)
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. / O Método dos Elementos Finitos Generalizados (MEFG) é essencialmente baseado no método da partição da unidade, que explora o conceito de partição da unidade para compatibilizar um conjunto de funções escolhidas para localmente aproximar de forma eficiente a solução. Apesar de suas vantagens bem conhecidas, o método pode apresentar algumas desvantagens. Por exemplo, o aumento do espaço de aproximação por meio das funções de enriquecimento pode introduzir dependências lineares no sistema de equações resolvente, assim como o aparecimento de elementos de mistura. Para contornar as desvantagens apontadas acima, algumas versões aprimoradas do MEFG foram desenvolvidas. O MEFG Estável é uma primeira versão aqui considerada na qual as funções de enriquecimento do MEFG são modificadas. O MEFG Estável de ordem superior propõe uma modificação adicional para a geração das funções de forma atreladas ao espaço enriquecido. Esta pesquisa visa apresentar e testar numericamente essas novas versões do MEFG recentemente propostas. Além de destacar suas principais características, alguns aspectos sobre a integração numérica quando usado o MEFG Estável de ordem superior, em particular, são também abordados. Por exemplo, detalha-se uma regra de divisão da área do elemento quadrilateral, guiada pela própria definição de sua partição da unidade. Os exemplos escolhidos para os experimentos numéricos consistem em chapas com geometrias favoráveis para explorar as vantagens de cada método. Essencialmente, examinam-se funções singulares com boas propriedades de aproximar a solução nas vizinhanças de vértices de cantos, bem como funções polinomiais para aproximar soluções suaves. Ademais, uma comparação entre o MEF convencional e os métodos aqui descritos é feita levando-se em consideração o número de condição do sistema escalonado e as razões de convergência do erro relativo em deslocamento. Finalmente, os experimentos numéricos mostram que o MEFG Estável de ordem superior é a mais robusta e confiável entre as versões do MEFG testadas.
|
Page generated in 0.0961 seconds