Neste trabalho desenvolveu-se uma formulação adaptativa do método de elementos de contorno (MEC) para a análise de problemas de fratura elástica linear. Foi utilizado o método da colocação para a formulação das equações integrais de deslocamento e de tensão. Para a discretização das equações integrais foram utilizados elementos lineares que possibilitaram a obtenção das expressões exatas das integrais (integração analítica) sobre elementos de contorno e fratura. Para a montagem do sistema de equações algébricas foram utilizadas apenas equações de deslocamento, apenas equações de forças de superfície, ou as duas escritas para nós opostos da fratura levando, portanto ao método dos elementos de contorno dual usualmente empregado na análise de fratura. Para o processo de crescimento da trinca foi desenvolvido um procedimento especial objetivando a correta determinação da direção de crescimento da trinca. Os fatores de intensidade de tensão são calculados por meio da conhecida técnica de correlação de deslocamentos a qual relaciona os deslocamentos atuantes nas faces da fissura. Após a determinação dos fatores de intensidade de tensão é utilizada a teoria da máxima tensão circunferencial para a determinação do ângulo de propagação. O modelo adaptativo empregado é do tipo h onde apenas a sub-divisão dos elementos é feita com base em erros estimados. O erro a ser considerado foi estimado a partir de normas onde se consideraram: a variação aproximada dos deslocamentos, a variação das forças de superfície e a variação da energia de deformação do sistema, calculada com a sua integração sobre o contorno. São apresentados exemplos numéricos para demonstrar a eficiência dos procedimentos propostos. / In this work, an adaptative formulation of the boundary element method is developed to analyze linear elastic fracture problems. The collocation point method was used to formulate the integral equations for the displacements and stresses (or tractions). To discretize the integral equations, linear elements were used to obtain the exact expressions of the integrals over boundary elements and fracture. To construct the linear system of equations were used only displacement equations, traction equations or both of them written for opposite nodes of the fracture, leading to the dual boundary element formulation usually employed in the fracture analyses. For the process of growth of the crack a special procedure was developed aiming at the correct determination of the direction of growth of the crack. The stress intensity factors, to calculate he crack growth angle, are calculated through of correlation displacements technique which relates the displacements actuants in the faces of the crack. The employed adaptative model is the h-type where only the sub-division of the elements is done based on error estimate. The error estimates considered in this work are based on the following norms: displacement, traction and strain energy variations, this last considered from the integration over the boundary. Numerical examples are presented to demonstrate the efficiency of the proposed procedures.

O principal objetivo deste trabalho é estudar as singularidades e hipersingularidades existentes nas formulações: singular - clássica - e hipersingular no Método dos Elementos de Contorno (MEC). Também é proposto um esquema residual h-adaptativo para a solução numérica do problema físico governado pela equação de Laplace. Usa-se malha poligonal, juntamente, com funções de interpolação - distribuição - de forma, dos tipos: constantes e lineares. Para controlar o erro a posteriori, é considerado o valor do resíduo, fora dos pontos de colocação. Também é testada uma técnica de quadratura numérica chamada adaptativa, específica para subelementos, no sentido de verificar se a precisão no cálculo das integrais com singularidades é melhorada. O uso de funções hierárquicas é discutido na forma de um algoritmo para atualização da matriz principal do sistema linear. / The main purpose of this work is to study the existing singularities and hypersingularities in the Boundary Element Method (BEM) with singular - classical - and hypersingular formulations. Also, an h-adaptive residual scheme for the numerical solution of the physical problem, driven by Laplace equation, is proposed. Boundary polygonal mesh, with constant and linear interpolation - distribution - shape functions together are used. To control the a posteriori error, is considered the residue value outside the collocation points. Also, a sub-element specific adaptive numerical quadrature technique, in an effort to verify if the precision when dealing with integrals possessing singularities is increased, is tested. The use of hierarchical functions is discussed, as an algorithm to update the linear system main matrix.

Neste trabalho desenvolveu-se uma formulação adaptativa do método de elementos de contorno (MEC) para a análise de problemas de fratura elástica linear. Foi utilizado o método da colocação para a formulação das equações integrais de deslocamento e de tensão. Para a discretização das equações integrais foram utilizados elementos lineares que possibilitaram a obtenção das expressões exatas das integrais (integração analítica) sobre elementos de contorno e fratura. Para a montagem do sistema de equações algébricas foram utilizadas apenas equações de deslocamento, apenas equações de forças de superfície, ou as duas escritas para nós opostos da fratura levando, portanto ao método dos elementos de contorno dual usualmente empregado na análise de fratura. Para o processo de crescimento da trinca foi desenvolvido um procedimento especial objetivando a correta determinação da direção de crescimento da trinca. Os fatores de intensidade de tensão são calculados por meio da conhecida técnica de correlação de deslocamentos a qual relaciona os deslocamentos atuantes nas faces da fissura. Após a determinação dos fatores de intensidade de tensão é utilizada a teoria da máxima tensão circunferencial para a determinação do ângulo de propagação. O modelo adaptativo empregado é do tipo h onde apenas a sub-divisão dos elementos é feita com base em erros estimados. O erro a ser considerado foi estimado a partir de normas onde se consideraram: a variação aproximada dos deslocamentos, a variação das forças de superfície e a variação da energia de deformação do sistema, calculada com a sua integração sobre o contorno. São apresentados exemplos numéricos para demonstrar a eficiência dos procedimentos propostos. / In this work, an adaptative formulation of the boundary element method is developed to analyze linear elastic fracture problems. The collocation point method was used to formulate the integral equations for the displacements and stresses (or tractions). To discretize the integral equations, linear elements were used to obtain the exact expressions of the integrals over boundary elements and fracture. To construct the linear system of equations were used only displacement equations, traction equations or both of them written for opposite nodes of the fracture, leading to the dual boundary element formulation usually employed in the fracture analyses. For the process of growth of the crack a special procedure was developed aiming at the correct determination of the direction of growth of the crack. The stress intensity factors, to calculate he crack growth angle, are calculated through of correlation displacements technique which relates the displacements actuants in the faces of the crack. The employed adaptative model is the h-type where only the sub-division of the elements is done based on error estimate. The error estimates considered in this work are based on the following norms: displacement, traction and strain energy variations, this last considered from the integration over the boundary. Numerical examples are presented to demonstrate the efficiency of the proposed procedures.

O principal objetivo deste trabalho é estudar as singularidades e hipersingularidades existentes nas formulações: singular - clássica - e hipersingular no Método dos Elementos de Contorno (MEC). Também é proposto um esquema residual h-adaptativo para a solução numérica do problema físico governado pela equação de Laplace. Usa-se malha poligonal, juntamente, com funções de interpolação - distribuição - de forma, dos tipos: constantes e lineares. Para controlar o erro a posteriori, é considerado o valor do resíduo, fora dos pontos de colocação. Também é testada uma técnica de quadratura numérica chamada adaptativa, específica para subelementos, no sentido de verificar se a precisão no cálculo das integrais com singularidades é melhorada. O uso de funções hierárquicas é discutido na forma de um algoritmo para atualização da matriz principal do sistema linear. / The main purpose of this work is to study the existing singularities and hypersingularities in the Boundary Element Method (BEM) with singular - classical - and hypersingular formulations. Also, an h-adaptive residual scheme for the numerical solution of the physical problem, driven by Laplace equation, is proposed. Boundary polygonal mesh, with constant and linear interpolation - distribution - shape functions together are used. To control the a posteriori error, is considered the residue value outside the collocation points. Also, a sub-element specific adaptive numerical quadrature technique, in an effort to verify if the precision when dealing with integrals possessing singularities is increased, is tested. The use of hierarchical functions is discussed, as an algorithm to update the linear system main matrix.