Formulação do MEC considerando efeitos microestruturais e continuidade geométrica G1: tratamento de singularidade e análise de convergência
Data
2017-11-15
Autores
Título da Revista
ISSN da Revista
Título de Volume
Editor
Biblioteca Digital de Teses e Dissertações da USP
Universidade de São Paulo
Escola de Engenharia de São Carlos
Universidade de São Paulo
Escola de Engenharia de São Carlos
Resumo
Descrição
Neste trabalho, uma abordagem micromecânica com aproximação da geometria dada por funções de Bézier triangulares com continuidade geométrica G1 é inserida ao Método dos Elementos de Contorno, o qual é aplicado em problemas da elastostática tridimensional. Para consideração do efeito microestrutural, foi utilizado a teoria gradiente elástica simplificada de Aifantis, a qual é uma particularização da teoria geral de Mindlin. Nesta teoria, um argumento variacional é estabelecido para determinar todas as possíveis condições de contorno, clássica e não-clássica, para o problema de valor de contorno geral. A partir deste argumento, a solução fundamental da elasticidade gradiente é explicitada e com o auxílio da identidade integral recíproca é construído a representação integral de contorno. Para tornar o problema de valor de contorno bem-posto, em adição à representação integral de contorno para deslocamento, uma segunda representação integral para derivada normal do deslocamento foi utilizada. Expressões integrais para deslocamento e tensão em pontos internos são apresentadas. Todos os núcleos das equações integrais são explicitamente desenvolvidos. Para a discretização do MEC foram utilizados elementos triangulares curvos, aproximados tanto para a geometria quanto para os parâmetros físicos por funções de Proriol (com características espectrais) e por funções aqui chamadas de Polinomiais, onde esta última é construída a partir de uma base nodal equidistante e pela imposição da partição da unidade. Entretanto estas funções aproximadoras garantem apenas continuidade C0 entre os elementos triangulares, ou seja, a garantia da continuidade do plano tangente não necessariamente é satisfeita. Com o objetivo de anular o termo de integral de linha presente na formulação microestrutural, a hipótese de superfície suave se faz necessária e assim funções de Bézier com continuidade geométrica G1, a qual depende apenas da posição e das normais dos nós nos vértices da malha triangular é utilizada. Para auxiliar na obtenção das coordenadas e das normais nodais para geometrias complexas foi utilizado o software de computação gráfica BlenderTM 2.7, o qual foi acoplado ao programa do MEC elastostático gradiente. Na sequência foi verificada, por meio de exemplos, a suavidade na intersecção entre os elementos triangulares G1 e estes foram comparados com as aproximações de Proriol e Polinomial. Em seguida, as singularidades presentes nas soluções fundamentais foram tratadas através da expansão em série de Laurent aplicada à técnica de subtração de singularidade. Condições necessárias e suficientes para a convergência das expansões em série das soluções fundamentais, estimador do erro para estas expansões, assim como, a correlação matemática entre o tamanho da malha e o parâmetro micromecânico g foram estabelecidos. Expressões explicitas da série de Laurent dos núcleos das integrais singulares e hipersingulares do MEC clássico e não clássico foram apresentadas. A verificação do tratamento da singularidade aplicado a elementos triangulares curvos foi realizada, tanto na direção radial quanto na direção angular. E pôde ser observado que ocorre uma perda de eficiência no tratamento da singularidade na direção angular, devida a presença do efeito de camada limite para elementos curvos distorcidos. Entretanto, este efeito de quase singularidade pode ser amenizado por meio da abordagem micromecânica, uma vez que foi observado menor presença do efeito da camada limite à medida que o parâmetro g é diminuído. Por último, foi desenvolvido um programa na linguagem FORTRAN 11.0, o qual contempla as abordagens clássica e micromecânica com continuidade geométrica G1. Sua validação foi feita por meio de exemplos considerados Benchmarks.
In this work, a micromechanical approach with approximation of geometry solved by Bézier triangular functions that guaranty continuity G1 is inserted to the Boundary element Method (BEM). This formulation is applied in three-dimensional elastostatic problems. The simplified elastic gradient theory proposed by Aifantis, which is a particularization of the general theory of Mindlin is used to consider the microstructural effect. In this theory a variational argument is established to determine all possible boundary conditions, classical and non-classical, for the general boundary value problem. From this argument, the fundamental solution of the gradient elasticity is explicited and by the reciprocal integral identity the boundary integral representation is achieved. In addition to the boundary integral representation for dispacement, a second integral representation regarding its normal derivative is used to make the well-posed boundary value problem. Integral expressions for displacement and stress on internal points are also presented. All kernels in the integral equations are explicitly developed. Curved triangular elements are used for the discretization of the BEM. The approximation of both the geometry and physical parameters is performed by Proriol functions (with spectral characteristics) and by Polynomial functions. The last is built from an equidistant nodal basis enforcing the partition of unity. However these approximating functions ensure only C0 continuity between the triangular elements, that is, the tangent plane continuity assurance is not necessarily satisfied. In order to cancel line integral terms in the microstructural approach, the hypothesis of smooth surface is required and thus Bézier function with geometric continuity G1, which depends only on the position and the normal of the nodes at the vertices of the triangular mesh is used. In this study the computer graphics software called BlenderTM 2.7 is used to assist in obtaining coordinates and normal vectors at nodes when complex geometries are analyzed. BlenderTM 2.7 is coupled to the gradient elastic BEM program. The smoothness of the resulting mesh using G1 elements is compared to Proriol and Polynomial approximations by means of simple examples. The singularities present in the fundamental solutions are treated by employing the expansion in Laurent series and the singularity subtraction technique. Necessary and sufficient conditions for the convergence of expansions in series of fundamental solutions, error estimator for these expansions, as well as the mathematical correlation between the size of the mesh and the micromechanical parameter, g, are established. Explicit expressions of Laurent series of the classical and micromechanical kernels forthe singular and hipersingular BEM integrals are presented. Treatment of singularity, both in the radial direction and in the angular direction, applied to curved triangular elements is verified. It can be observed that there is a loss of efficiency in the treatment of singularity in the angular direction, due to the presence of the boundary layer effect for distorted curved boundary elements. However, this nearly singularity effect could be alleviated by micromechanics approach, since minor boundary layer effect was observed as the parameter g is decreased. Finally, using FORTRAN 11.0 language, a computational code is developed, which includes the classic and micromechanics approach with geometric continuity G1, and its results are validated by means of Benchmark examples.
In this work, a micromechanical approach with approximation of geometry solved by Bézier triangular functions that guaranty continuity G1 is inserted to the Boundary element Method (BEM). This formulation is applied in three-dimensional elastostatic problems. The simplified elastic gradient theory proposed by Aifantis, which is a particularization of the general theory of Mindlin is used to consider the microstructural effect. In this theory a variational argument is established to determine all possible boundary conditions, classical and non-classical, for the general boundary value problem. From this argument, the fundamental solution of the gradient elasticity is explicited and by the reciprocal integral identity the boundary integral representation is achieved. In addition to the boundary integral representation for dispacement, a second integral representation regarding its normal derivative is used to make the well-posed boundary value problem. Integral expressions for displacement and stress on internal points are also presented. All kernels in the integral equations are explicitly developed. Curved triangular elements are used for the discretization of the BEM. The approximation of both the geometry and physical parameters is performed by Proriol functions (with spectral characteristics) and by Polynomial functions. The last is built from an equidistant nodal basis enforcing the partition of unity. However these approximating functions ensure only C0 continuity between the triangular elements, that is, the tangent plane continuity assurance is not necessarily satisfied. In order to cancel line integral terms in the microstructural approach, the hypothesis of smooth surface is required and thus Bézier function with geometric continuity G1, which depends only on the position and the normal of the nodes at the vertices of the triangular mesh is used. In this study the computer graphics software called BlenderTM 2.7 is used to assist in obtaining coordinates and normal vectors at nodes when complex geometries are analyzed. BlenderTM 2.7 is coupled to the gradient elastic BEM program. The smoothness of the resulting mesh using G1 elements is compared to Proriol and Polynomial approximations by means of simple examples. The singularities present in the fundamental solutions are treated by employing the expansion in Laurent series and the singularity subtraction technique. Necessary and sufficient conditions for the convergence of expansions in series of fundamental solutions, error estimator for these expansions, as well as the mathematical correlation between the size of the mesh and the micromechanical parameter, g, are established. Explicit expressions of Laurent series of the classical and micromechanical kernels forthe singular and hipersingular BEM integrals are presented. Treatment of singularity, both in the radial direction and in the angular direction, applied to curved triangular elements is verified. It can be observed that there is a loss of efficiency in the treatment of singularity in the angular direction, due to the presence of the boundary layer effect for distorted curved boundary elements. However, this nearly singularity effect could be alleviated by micromechanics approach, since minor boundary layer effect was observed as the parameter g is decreased. Finally, using FORTRAN 11.0 language, a computational code is developed, which includes the classic and micromechanics approach with geometric continuity G1, and its results are validated by means of Benchmark examples.
Palavras-chave
Análise de convergência, Tratamento de singularidade, Série de Laurent, Método dos Elementos de Contorno (MEC) micromecânico, Continuidade geométrica G1, Laurent series, Geometric continuity G1, Micromechanical BEM, Convergence analysis, Treatment of singularity