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.

O estudo do comportamento mecânico do complexo sistema advindo da interação entre solo/subestrutura/superestrutura é o tema do trabalho. Neste contexto, a representação do maciço é feita usando-se o método dos elementos de contorno (MEC) em abordagem 3D, de maneira que se possa simular o maciço com características mecânicas não-homogêneas, além de se considerar uma camada de apoio indeslocável a distâncias prescritas a priori e condição de aderência perfeita. A subestrutura também é representada via MEC tridimensional, a qual está imersa dentro deste meio heterogêneo. A infra e a superestrutura são modeladas empregando o método dos elementos finitos (MEF), com o uso de elementos estruturais reticulares e elementos laminares. São apresentados alguns exemplos em que se valida a formulação e outros que demonstram a potencialidade e a necessidade de se empregar a formulação para a melhor análise do complexo fenômeno em estudo. Por fim, demonstra-se a obrigatoriedade de se otimizar a formulação, empregando-se duas grandes ferramentas numéricas: o paralelismo e o emprego de um adequado método de resolução de sistemas esparsos. / The analysis of the soil-structure system interaction is a vast field of interest in the area of civil engineering. A realistic representation of its behaviour. Thus, in the present research, the soil is considered a non-homogeneous continuum supported by a rigid and adhesive interface and modelled by boundary element method via Kelvin solution in 3D space. The foundation is also modelled by this above-mentioned modelling technique. The raft foundation and the superstructure are represented by finite shell and 3D frame elements. In order to estimate the accuracy and the potentiality of the proposed numerical formulation, some examples are validated when compared to similar approaches, and others simulations are presented to stress the necessity of coupling the non-homogeneous soil-foundation-radier-superstructure system as a whole. Finally, to acquire numerical time efficiency, it is shown that it is imperative to apply parallel processing and sparse techniques for the solution of the final system.

O presente trabalho trata do desenvolvimento de formulações numéricas para avaliar o comportamento mecânico de estruturas compostas planas, no contexto de elasticidade linear e mecânica da fratura não linear. As formulações propostas são baseadas no Método dos Elementos de Contorno (MEC), por meio das representações integrais singular e hiper singular dos problemas elastostáticos. A técnica de multi-regiões é considerada para acoplar a interface de sólidos multifásicos. O MEC é uma técnica numérica robusta e precisa para analisar o fenômeno da fratura em sólidos. Esse método numérico apresenta uma natural redução na dimensionalidade do problema, tornando mais simples a modelagem das superfícies de fratura. Além disso, essa redução de dimensionalidade faz também com que o tratamento de interfaces materiais em estruturas compostas seja uma tarefa menos árdua. Com o uso da solução fundamental de Kelvin nas representações integrais, materiais isotrópicos podem ser considerados para constituir as estruturas compostas. Por outro lado, utilizando a solução fundamental de Cruse & Swedlow, também é possível lidar, de maneira geral, com materiais anisotrópicos em estruturas compostas. Nessas estruturas, as fraturas são assumidas como ocorrendo ao longo das interfaces e o comportamento não linear é introduzido pelo modelo coesivo de fratura, o qual é aplicável a materiais quase frágeis. Nessas análises, o sistema não linear de equações pode ser solucionado utilizando dois distintos algoritmos de resolução iterativa. O primeiro sempre leva em consideração a rigidez elástica da estrutura e é, portanto denominado Operador Constante (OC). Já o segundo é denominado Operador Tangente (OT), pois considera uma rigidez tangente à resposta estrutural não linear, o que resulta em melhores taxas de convergência em comparação ao OC. Como aplicações das formulações, estruturas compostas teóricas foram analisadas em regime elástico. Além disso, testes experimentais de fratura em espécimes de concreto e madeira também foram simulados. A comparação dos resultados com as referências demonstrou que, as formulações foram efetivas e precisas para avaliar respostas mecânicas de estruturas, seja em regime elástico linear ou nos testes de fratura quase frágil. / The present work deals the development of numerical formulations to evaluate the mechanical behaviour of plane composed structures, in the context of linear elasticity and nonlinear fracture mechanics. The proposed formulations are based on the Boundary Element Method (BEM), through its classical singular and hyper singular integral equations. The multi-region technique is adopted to couple the interfaces of non-homogeneous multiphase bodies. The BEM is a robust and accurate numerical technique to analyse fracture phenomena in solids. This numerical method presents a mesh dimensionality reduction, which makes easier the modelling of cracks surfaces. Besides, this dimensionality reduction also makes the treatment of interfaces in composed structures a less complex task. Considering the use of Kelvin fundamental solutions at the integrals equations, isotropic materials can be represent as parts of the composed structures. On the other hand, using Cruse & Swedlow fundamental solution it is also possible to deal with general anisotropic materials. At the composed structures, cracks can propagate along the materials interfaces and the cohesive crack model is responsible for the nonlinear structural behaviour of the quasi-brittle failures. The nonlinear system of equations at the fracture analyses is solved using two different algorithms for iterative resolution. The first always takes into account the structure elastic strength and, hence it is called Constant Operator (CO). On the other hand, the second is denominated Tangent Operator (TO) due to the fact that it considers strengths at the tangent directions of the nonlinear structural response. Therefore, convergence rates are faster when compared with the CO. As applications, composed structures were analysed with the developed formulations in linear elastic range. In addition, experimental fracture testes performed in concrete and wood specimens were also analysed. The confront of obtained results with the reference ones show that, the formulation was effective and accurate to evaluate the mechanical responses of composed structures in linear elastic range, and also to perform nonlinear quasi-brittle fracture tests.

No presente trabalho, uma combinação do método dos elementos de contorno (MEC) com o método dos elementos finitos (MEF) é apresentada para a análise da interação entre estacas e o solo, considerado como um meio infinito tridimensional e homogêneo. O meio contínuo tridimensional de domínio infinito é modelado pelo MEC, enquanto as estacas consideradas como elementos reticulares são tratadas pelo MEF. As equações das estacas oriundas do método dos elementos finitos são combinadas com as do meio contínuo obtidas a partir do método dos elementos de contorno, resultando em um sistema completo de equações, que convenientemente tratadas, proporcionam a formulação de coeficientes de rigidez do conjunto solo-estacas. Finalmente, uma formulação para a análise do comportamento não-linear do solo na interface com a estaca é desenvolvida, tornando o modelo mais abrangente. / In the present work a combination of the Boundary Element Method (BEM) and the Finite Element Method (FEM) is used for pile-soil interaction analyses, considering the soil as a homogeneous, three-dimensional and infinite medium. The three-dimensional infinite continuous medium is modeled by the BEM, and the piles are, considered as beam elements, modeled by the FEM. This combination also is used for studying the interaction of plates sitting on a continuous medium. The pile equations generated from the FEM are combined with the medium equations generated from the BEM, resulting a complete equation system. Manipulating properly this equation system, a set of stiffness coefficients for the system soil-pile is obtained. Finally, to make the model more comprehensive, it presented a formulation to take into account the soil nonlinear behavior at the pile interface.

O objetivo do presente trabalho é o desenvolvimento de formulações e de códigos computacionais que possibilitem a análise bidimensional estática de meios contínuos anisotrópicos viscoplásticos reforçados ou não por fibras. Especificamente, as análises numéricas envolvem aplicações dos métodos dos elementos de contorno (MEC) e dos elementos finitos (MEF), comparando-se os resultados obtidos com respostas analíticas e experimentais, disponíveis na literatura, buscando-se assim, subsídios teóricos que permitam o entendimento de problemas mais gerais envolvendo meios anisotrópicos. Para tanto são empregados elementos finitos triangulares com aproximações cúbica e quadrática para os deslocamentos na modelagem dos domínios. Na consideração do reforço com fibras, elementos finitos de barras simples são empregados. A formulação desenvolvida proporciona também a consideração de distribuição randômica das fibras imersas no meio sem qualquer aumento dos graus de liberdade do problema analisado, diferindo-se assim, das formulações conhecidas até o momento. Com o MEC, a análise de plasticidade e viscoplasticidade em meios com anisotropia geral é feita de maneira original no trabalho, destacando-se a consideração de lei de fluxo plástico não-associativa e o tratamento de viscosidade apenas com integrais de contorno, sem a utilização de aproximações de domínio. Uma quantidade significativa de exemplos é apresentada, possibilitando a verificação da eficiência das formulações e dos códigos desenvolvidos / The objective of the present work is the development of formulations and computational codes that enable the static bidimensional analysis of the viscoplastic anisotropic medias reinforced, or not, by fibers. Specifically, the numerical analysis involve applications of the boundary elements method (BEM) and finite elements (FEM), comparing the results obtained with analytical and experimental solutions available in the literature, allowing the understanding of general problems in anisotropic media. Two-dimensional finite elements with cubic and quadrate approximations for the displacements are used to model domains. Reinforcements are modeled by truss finite elements. The developed formulation provides the consideration of random distribution of the fibers, without any additional degree of freedom of the problem. With the BEM, the plasticity and viscoplasticity analysis in general anisotropic medias is originally developed in the present work, emphasizing the consideration of non-associative plastic flow and the treatment of viscosity just with boundary integrals, without domain approximation. Various examples are shown in order to verify the efficiency of the proposed formulation and developed computational codes

Neste trabalho, propõe-se uma formulação não linear baseada no método dos elementos de contorno, para representação de domínios poro-elasto-plásticos reforçados. Esta formulação é apresentada para os casos saturado e não saturado. Para o problema poroso enrijecido um acoplamento com o método dos elementos finitos é empregado, e a técnica de mínimos quadrados permite a regularização dos deslocamentos e do vetor de forças de superfície ao longo das interfaces de acoplamento. São empregadas expressões analíticas para o tratamento das integrais de contorno e de domínio presentes na formulação do método dos elementos de contorno. A formulação de Biot é empregada para a descrição de meios porosos saturados e uma formulação energética baseada nos trabalhos de Coussy é adaptada para a extensão ao caso não saturado. Neste caso, a pressão capilar e energia das interfaces são levadas em consideração. O nível de saturação é descrito pelo modelo de Van Genuchten e o comportamento do esqueleto é descrito ou pelo modelo de Drucker-Prager ou pelo modelo de Cam-Clay modificado. O problema não linear obtido por uma descrição temporal associada a discretização espacial é resolvido pelo método de Newton-Raphson. No caso saturado, o operador tangente consistente é definido e utilizado para obtenção da solução do sistema. Exemplos numéricos são apresentados para validar a formulação proposta. / In this work a nonlinear formulation of the boundary element method (BEM) is proposed to deal with saturated and unsaturated poro-elasto-plastic 2D reinforced domains. To model reinforced porous domains a BEM/FEM (Finite Element Method) modified coupling technique is employed. The coupling is made by using the least square method to regularize the displacement and traction distributions along the interfaces. Analytical expressions have been derived for all boundary and domain integrals required for the formulation. The Biot formulation is used for the description of the saturated porous environments and an energetic consistent formulation based on work of Coussy is adopted for its extension to the framework of unsaturated porous media. In this case, the capillar pressure and the interface energy are taken into account. The Van Genuchten model is used for the determination of saturation level in non-saturated poro-elasto-plastic problems. The Drucker-Prager modified model if used for the saturated poro-elasto-plastic problems and the modified Cam-Clay model for the representation of non-saturated poro-elasto-plastic problems. For the saturated case, the consistent tangent operator is derived and employed inside a Newton procedure to solve non-linear problems. Numerical solutions are presented to validate the proposed models.

O objetivo do presente trabalho é o desenvolvimento de um código computacional que possibilite a análise dinâmica de estruturas tridimensionais em regime elástico-linear acopladas ao solo, tratado como meio infinito elastoplástico. As superestruturas são tratadas por elementos finitos simples de casca e de barra geral, as estruturas de fundações são tratadas por elementos de casca que simulam o contato com o solo, modelando radiers, túneis e reservatórios enterrados. Blocos são modelados por elementos de contorno tridimensionais. O solo é modelado de duas maneiras distintas: na região plastificada emprega-se a solução fundamental de Kelvin (estática) e na região não plastificada (elástica) adota-se a solução fundamental do problema de Stokes. O acoplamento entre os meios é feito aplicando-se a técnica de subregiões. Deve ficar claro que todo procedimento estático equivalente foi implementado. Vários exemplos numéricos são apresentados, onde se percebe a eficiência do código computacional desenvolvido / The objective of the present work is the development of a computational code that makes possible dynamic analyses of three-dimensional structures in elastic-linear behavior coupled to the soil, modeled as elastoplastic infinite medium. Simple finite elements, shell and general bars, are used to model elastic structures. The structures of foundations are modeled by shell’s elements which simulate the contact with the soil, modeling radiers, tunnels and buried reservoirs. Blocks are modeled by three-dimensional boundary elements. The soil is modeled in two different ways: in the plastic region Kelvin’s fundamental solution (static) is used and in the elastic region the fundamental solution of the Stoke’s problem is adopted. The coupling among the media is done applying the sub-region technique. It is important to note that the equivalent static procedure has been implemented. Several numerical examples are presented, demonstrating the efficiency of the developed computational code

Submitted by Renata Lopes (renatasil82@gmail.com) on 2016-05-31T15:07:36Z No. of bitstreams: 1 lucasvitorfonsecadeoliveira.pdf: 1420167 bytes, checksum: b5075b7ccc33c587012a46901a9b21af (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2016-07-02T12:44:37Z (GMT) No. of bitstreams: 1 lucasvitorfonsecadeoliveira.pdf: 1420167 bytes, checksum: b5075b7ccc33c587012a46901a9b21af (MD5) / Made available in DSpace on 2016-07-02T12:44:37Z (GMT). No. of bitstreams: 1 lucasvitorfonsecadeoliveira.pdf: 1420167 bytes, checksum: b5075b7ccc33c587012a46901a9b21af (MD5) Previous issue date: 2012-03-30 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Com a abertura legal do setor elétrico brasileiro na década de 90 para investimentos privados, em especial com a implantação sistemática de leilões para definição dos concessionários de transmissão, diversos estudos foram efetuados com o objetivo de viabilizar projetos de linha de transmissão de modo a concorrem nesses leilões. Grandes transmissoras, como Furnas, Chesf e Eletronorte promoveram pesquisas em consórcio com empresas construtoras de linhas e fabricante de ferragens e isoladores, com o objetivo de desenvolverem projetos com baixa relação Reais/MW transmitido. Neste cenário, muitos projetos de linhas de potência natural elevada foram idealizados e implantados, e como fruto desses projetos, publicações foram realizadas descrevendo os resultados obtidos. Todavia, apesar da quantidade, qualidade e riqueza das informações publicadas, a descrição do processo de otimização dos parâmetros elétricos da linha, em especial da impedância característica, não receberam destaque, sendo muita das vezes suprimida nos artigos. Neste sentido, este trabalho propõe a investigação de um método inovador de otimização da capacidade de transmissão de linhas de potência natural elevada, através da análise indireta dos efeitos da variação do posicionamento dos cabos condutores no feixe sobre a impedância característica da linha, por meio de uma abordagem original que utiliza basicamente o vetor de Poynting como função objetivo do Método de Otimização por Enxame de Partículas, sendo os campos elétricos e magnéticos necessários para definição do vetor de Poynting calculados por meio do Método dos Elementos de Contorno. O método foi validado utilizando como exemplos a otimização das configurações de feixes da linha de 500 kV Interligação Norte/SUL III – Trecho 2 e Linha de 500 kV Presidente Dutra / Teresina / Sobral / Fortaleza. Foram encontrados ganhos na capacidade de transmissão de 7% e 22% respectivamente devido à redução da impedância característica calculada após a otimização através da rotina Line Constants do programa ATP/EMTP, e comparando-a com os valores originais. / With the legal opening for private investments in the Brazilian energy sector in the 90's, in particular the systematic implementation of auctions for defining power transmission concessions, several studies were made with the aim of developing transmission line projects in order to compete in these auctions. Major transmission companies such as Furnas, Eletronorte and CHESF promoted consortium research with line builders and hardware and insulators manufacturers, in order to develop projects with low cost/MW transmitted. In this scenario, many projects of high surge-impedance loading lines were developed and implemented, and as a result of these projects, publications were made describing the results. However, despite the quantity, quality and resourcefulness of published information, the descriptions of the optimization process of electrical line parameters, especially regarding characteristic impedance, were not given prominence, being often suppressed from the articles. Thus, this study proposes the investigation of a method for optimizing the transmission capacity of high surge-impedance loading lines, varying the power cables in the bundle, indirectly reducing its characteristic impedance through an original approach that uses basically the Poynting's vector as objective function of the Particle Swarm Optimization method.The electric and magnetic fields needed for defining the Poynting vector were calculated using the Boundary Element Method. The method was validated through the optimization of bundle configuration, using as a model the characteristics of the 500 kV line North / South Interconnection III - Segment 2 and the 500 kV line Presidente Dutra / Teresina / Sobral / Fortaleza. It was found transmission capacity gains of 7% and 22% respectively, by reducing the characteristic impedance, which was calculated after the optimization using the EMTP/ATP Line Constants Program, comparing it with the original values.

Submitted by Renata Lopes (renatasil82@gmail.com) on 2017-02-23T19:26:50Z No. of bitstreams: 1 francianeconceicaopeters.pdf: 4730497 bytes, checksum: 201c60342a8bf9edc9b308fa50fafa54 (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2017-02-24T12:08:43Z (GMT) No. of bitstreams: 1 francianeconceicaopeters.pdf: 4730497 bytes, checksum: 201c60342a8bf9edc9b308fa50fafa54 (MD5) / Made available in DSpace on 2017-02-24T12:08:43Z (GMT). No. of bitstreams: 1 francianeconceicaopeters.pdf: 4730497 bytes, checksum: 201c60342a8bf9edc9b308fa50fafa54 (MD5) Previous issue date: 2010-01-27 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / A construção de imagens associadas à distribuição de condutividades no interior de um meio condutor a partir de injeção de corrente elétrica e medidas de potencial no contorno externo do corpo é uma técnica conhecida como tomografia por impedância elétrica. É um problema inverso que tem sido estudado visando aplicações biomédicas, monitoramento de processos industriais e investigação geofísica. Em alguns casos, é possível levar em consideração informações conhecidas sobre o domínio do corpo no processo de construção da imagem, recaindo no problema da detecção de inclusões que é o problema efetivamente tratado neste trabalho. Este problema pode ser resolvido por meio da minimização de uma função da diferença entre potenciais medidos no contorno e calculados para uma dada distribuição de condutividades. O presente trabalho desenvolve uma estratégia para a solução deste problema baseada na parametrização da geometria do contorno das inclusões cujas formas e dimensões se pretende determinar. O problema de minimização é resolvido por meio do Método de Levenberg-Marquardt e o problema direto via Método dos Elementos de Contorno. Para avaliar o desempenho da estratégia proposta são apresentados resultados numéricos envolvendo contornos definidos por splines, problemas com a presença de ruído nas medidas, avaliação de protocolos de injeção de corrente e medição de potencial elétrico e ainda uma aplicação voltada ao monitoramento cardíaco. / The images reconstruction of the conductivity distribution inside a conductive body based on electrical current injection and potential measurements on the outer boundary of this body is a technique known as electrical impedance tomography. This is an inverse problem that has been studied in biomedical applications, industrial process monitoring and geophysics investigation. In some cases, it is possible to take into account in the reconstruction process, informations about the body, leading to the problem of identifying inclusions, that is the problem actually treated in this work. This inverse problem can be solved by the minimization of a function, defined as the difference between the measured potentials and the computed ones for a given conductivity distribution. The present work describes a strategy to solve this problem based on the parametrization of the inclusions boundary, whose shape and size is intended to be determined. The minimization problem is solved via Levenberg-Marquardt Method and the forward one is solved via Boundary Elements Method. In order to evaluate the performance of the proposed strategy, numerical experiments with inclusions of boundaries defined by splines, problems with noisy data, current injection and potential measurement protocols and an application of the strategy to the cardiac function monitoring are presented.

[pt] Apresenta-se neste trabalho um procedimento de análise de propagação de trincas a partir de um programa de computador baseado na formulação do método consistente dos elementos de contorno para problemas bidimensionais. Este método tem como uma das suas principais características a solução exata dos problemas de singularidade presentes na formulação. Além disso, com esta metodologia é possível representar a geometria da trinca com aberturas micrométricas, de forma semelhante ao observado em ensaios laboratoriais. Neste estudo, são analisados os resultados de propagação em três estruturas com geometrias distintas, cada estrutura submetida a diferentes combinações de carga, com o objetivo de reproduzir modos puros de carregamento I e II, assim como modo misto de carregamento. É realizado um estudo sobre o tamanho dos incrementos utilizados nos modelos e do ângulo de propagação, possibilitando determinar que o tamanho ideal dos elementos de novos trechos deve se limitar à mesma dimensão dos elementos vizinhos, evitando possíveis erros numéricos, enquanto o ângulo de propagação pode ser determinado utilizando os fatores de intensidade de tensão (FIT) KI e KII, empregando o conceito de tensão principal máxima. O FIT é obtido por meio de deslocamentos recíprocos próximos à ponta da trinca, sendo realizado um estudo com um exemplo de referência para medir a confiabilidade da técnica, com diferenças de no máximo 7 por cento. O desempenho observado utilizando a metodologia adotada neste estudo é comparado com outros resultados encontrados na literatura, mostrando caminhos de propagação de trinca semelhantes em todas as simulações. No decorrer do trabalho são explicados os conceitos de mecânica da fratura linearmente elástica e da geometria da trinca adotada, assim como o desenvolvimento do código computacional. / [en] This work presents a crack propagation analysis procedure on a computer program based on the consistent boundary element formulation for two-dimensional problems. This method has as one of its main features the exact solution of the singularity problems present in the formulation. In addition, with this methodology it is possible to represent the crack geometry with micrometric openings, similar to the cracks presented in laboratory tests. In this study, the propagation results in three structures with different geometries are analyzed, each structure subjected to different load combinations, in order to reproduce pure loading modes I and II, as well as mixed loading modes. A study is carried out on the size of the increments used in the models and on the propagation angle, making it possible to determine that the ideal size of the elements of new sections should be limited to the same dimension of the neighboring elements, avoiding possible numerical errors, while the propagation angle can be determined using the stress intensity factors (FIT) KI e KII, employing the concept of maximum principal stress. The FIT is obtained through reciprocal displacements close to the crack tip, and a study is carried out with a reference example to measure the reliability of the technique, with differences of at most 7 per cent. The performance observed using the methodology adopted in this study is compared with other results found in literature, showing similar crack propagation paths in all simulations. In the course of the chapters, the concepts of linearly elastic fracture mechanics and the adopted crack geometry are explained, as well as the development of the computational code.