Este trabalho trata da análise da propagação de fissuras, independente do tempo, em domínios bidimensionais utilizando uma formulação alternativa do método dos elementos de contorno (MEC). O MEC vem sendo utilizado com sucesso na análise de diversos problemas de engenharia. Considerando problemas de mecânica da fratura, o MEC é especialmente eficiente devido à redução da dimensionalidade de sua malha, o que permite a simulação do crescimento das fissuras sem as dificuldades do processo de remalhamento. Nesta pesquisa, desenvolvem-se formulações não lineares do MEC para a análise da propagação de fissuras em materiais quase frágeis. Nesses materiais, a zona de processo à frente da ponta da fissura introduz efeitos fisicamente não lineares no comportamento estrutural. Assim, para a simulação da presença da zona de processo, modelos não lineares são necessários. Classicamente a formulação dual do MEC é utilizada para modelar propagação de fissuras na quais equações singulares e hipersingulares são escritas para elementos definidos ao longo das faces das fissuras. O presente trabalho propõe uma segunda formulação utilizando um campo de tensões iniciais para a representação da zona coesiva. Nesta formulação, o termo de domínio da equação integral clássica do MEC é degenerado, de forma a atuar somente ao longo do caminho de crescimento das fissuras, sendo que esse procedimento dá origem a uma nova variável denominada dipolo, responsável por garantir o atendimento das condições de contorno. Em conjunto com essa nova formulação, se propõe o uso do operador tangente (OT), que é deduzido no trabalho, a fim de acelerar o processo de convergência da solução. Os resultados obtidos, por meio da formulação alternativa, são comparados tanto com dados experimentais quanto com o MEC dual, ambos disponíveis na literatura. As respostas encontradas foram satisfatórias no sentido de conseguir reproduzir o comportamento real da estrutura explorando as vantagens computacionais proporcionadas pelo OT.

This work presents a time-independent crack propagation analysis, in two-dimensional domains, using an alternative boundary element method (BEM) formulation. BEM has been used successfully to analyze several engineering problems. Considering fracture mechanics problems, BEM is especially efficient due to its mesh reduction aspects, which allows the simulation of crack growth without remeshing difficulties. In this research, nonlinear BEM formulations are develop in order to analyze crack propagation in quasi-brittle materials. Considering these materials, the process zone ahead of the crack tip leads to nonlinear effects related to structural behavior. Thus, nonlinear models are required for simulating the presence of the process zone. Classically, the dual BEM is used for modeling the crack propagation, in which singular and hyper-singular equations are written for elements defined along the crack faces. This work proposes an alternative formulation using the initial stress field to represent the cohesive zone. In this formulation, the classic domain integral term is degenerated in order to be non-null only at the crack growth path. This procedure leads the creation of new variable called dipole, which is responsible for ensuring the compliance of the boundary conditions. In addition to this new formulation, it is proposed the use of the tangent operator (TO), which is derived in this work, in order to accelerate the convergence. The results obtained using the new formulation, are compared with experimental data and dual BEM results available in the literature. The responses were found satisfactory in reproducing the behavior of real structures exploiting the computational advantages provided by the TO.