A aplicação da análise inversa é objeto de estudo nos mais diversos campos da ciência e da engenharia. A motivação para o tratamento de tais problemas se deve ao fato de que em muitas aplicações dessas áreas do conhecimento, há a necessidade da identificação de parâmetros físicos e geométricos a partir de dados do domínio medidos experimentalmente, já que tais parâmetros de entrada são desconhecidos para uma análise direta do problema. Neste tipo de análise o problema principal está na quantidade e qualidade dos dados experimentais obtidos, que são na maioria das vezes insuficientes para garantir que o sistema gerado apresente solução única, gerando com isto um problema essencialmente mal-posto. Assim, de forma geral o emprego confiável da análise inversa implica na utilização de ferramentas eficientes de aquisição de dados experimentais aliada a técnicas numéricas de regularização que buscam a minimização da função objetiva gerada por algum método numérico, como por exemplo, o método dos elementos de contorno (MEC). Sendo assim, o presente trabalho tem por objetivo apresentar uma formulação para resolução de problemas inversos de valor de contorno e estimativa dos parâmetros do modelo coesivo, através de medidas de campos de deslocamentos, em sólidos bidimensionais com domínio formado por multi-regiões via (MEC), utilizando-se de técnicas tais como: mínimos quadrados, regularização de Tikhonov, decomposição em valor singular (SVD) e filtro de Tikhonov, para regularização do problema. Além disto, são apresentados alguns exemplos de aplicação da formulação desenvolvida. / The application of inverse analysis is nowadays subject of research of many fields in engineering and science. The motivation to consider this problem is due to the fact that in many applications of these knowledge areas, physical and geometric parameters, that are not directly known, can be identified using domain data measured experimentally. In this kind of analysis the main problem is the quantity and the quality of the obtained experimental data, which, many times, are not sufficient to guarantee that the generated system of equations has only one solution, leading therefore to an ill-posed problem. Thus, in general the reliable use of the inverse analysis requires using efficient tools for experimental data acquisition together with the numerical techniques of regularization needed to impose the minimization of the objective function written by using any numerical method, as the boundary element method (BEM) for instance. In this context, the objective of the present work is to derive a formulation for resolution of boundary-value inverse problems and to estimate the material parameters of the cohesive model, by using measured displacements fields, in multi-region two-dimensional solid by BEM, using techniques such as: least squares, Tikhonov regularization, singular value decomposition (SVD) and Tikhonov filtering, for the problem regularization. Some application examples are presented using the developed formulation to illustrate its performance.

Esse trabalho trata da análise de corpos multi-fraturados utilizando o método dos elementos de contorno. Este método numérico é conhecido por ser robusto e preciso neste tipo de problema e também por requerer pequeno esforço computacional na criação da malha de elementos para o crescimento das fissuras. Duas metodologias para a análise do comportamento das fissuras são consideradas. A primeira é a já consagrada metodologia dual. Por meio desta técnica equações integrais distintas são aplicadas às faces da fissura. Estas equações integrais são escritas em termos de deslocamentos e forças de superfície. A segunda metodologia é a que emprega a formulação singular onde a fissura é considerada como um vazio no domínio sendo as faces da fissura separadas por uma pequena distância. No tocante ao crescimento das fissuras foi desenvolvido um procedimento especial para a determinação da direção de crescimento das fissuras o qual mostrou-se muito eficiente levando a resultados precisos. O crescimento das fissuras é efetuado considerando o fator de intensidade de tensão atuante na extremidade de cada fissura. Dessa forma, as fissuras mais solicitadas apresentam maior comprimento de propagação tornando a análise mais realista. Os fatores de intensidade de tensão são calculados por meio de duas técnicas. A primeira é a já conhecida técnica de correlação de deslocamentos a qual relaciona os deslocamentos atuantes nas faces da fissura. Uma técnica alternativa é também utilizada a qual emprega o campo de tensões presente na extremidade da fissura. Após a determinação dos fatores de intensidade de tensão quatro diferentes teorias de interação de modos podem ser utilizadas para a determinação do ângulo de propagação. Foram analisadas estruturas sendo os resultados comparados aos previstos analiticamente e também numericamente. As respostas obtidas foram satisfatórias validando assim a metodologia proposta neste trabalho / This work deals with analysis of multi-fractured bodies using boundary element method. This numerical method is known to be robust and accurate in this kind of problem and by small computational effort to create elements mesh of crack growth. Two methodologies to analyze of crack behavior are considerate. The first is consecrated dual methodology. Through this technique different kind of integral equations are applied to crack boundaries. These integrals equations are written in displacements and traction variables. Second methodology is singular formulation. Through this technique crack is represented like a hole in body and the crack boundaries is separated by a small gap. For crack growth was created a special proceeding to determination crack growth direction. This method is very efficient and your results are accurate. Crack growth is made through the stress intensity factor performed in crack tip. Then the cracks more requested going to propagate with a larger length’s growth turning this model very realistic. The stress intensity factors are calculated through two techniques in this work. First is the known correlation displacement technique which related displacement in crack boundaries. An alternative technique is also used which consider stress field in crack tip. After determination of stress intensity factors four different theories are used to calculate the crack growth angle. In this work were analyzed structures with results are compared with analytical and numerical answers. The results obtained went very satisfactory validating the methodology proposed

Orientador: Paulo Sollero / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-10T03:16:07Z (GMT). No. of bitstreams: 1 Paiva_SeilaVastiFariade_M.pdf: 1387945 bytes, checksum: 04ad3c76f01698f24ef5fd7ebfb49264 (MD5) Previous issue date: 2006 / Resumo: Nesta dissertação é apresentada uma análise de problemas da mecânica da fratura elástica linear em estruturas bimateriais anisotrópicas. Utilizando o método dos elementos de contorno é possível calcular os fatores de intensidade de tensão em problemas planos (2D) devido à presença de trincas interfaciais entre as lâminas que compõem o material. A estrutura pode estar submetida à carregamento em modo I ou modo misto. O problema é modelado usando-se a técnica de sub-regiões para descrever cada um dos diferentes subdomínios, representado por cada material. Na interface das sub-regiões, em que o domínio é dividido, são impostas condições de equilíbrio de forças e continuidade de deslocamentos, exceto na região que corresponde à trinca. O comportamento singular apresentado pelo campo de tensões próximo à ponta da trinca é modelado com elementos de ponto a um quarto com singularidade de forças de superfície. São apresentados exemplos numéricos de problemas com carregamentos no plano. Foi também apresentada a análise de convergência de malhas, mostrando uma pequena dependência da discretização mesmo quando malhas pouco refinadas foram usadas. Alguns dos exemplos têm correspondentes na literatura, os quais foram utilizados para comparação com os resultados obtidos. Observou-se uma boa concordância na comparação dos resultados / Abstract: This thesis presents an analysis of problems of linear elastic fracture mechanics in anisotropic bimaterial structures. Using the boundary element method, it is possible to evaluate stress intensity factors in plane problems (2D) due to the presence of interfacial cracks between the laminae that constitute the material, when the structure is submitted to a mode I or in mixed mode loading. The problem is modeled using the subregion technique to describe each one of the different subdomains, represented by each material. On the interface of subregions, which the domain is divided, conditions of tractions equilibrium and displacements continuity are imposed, except in the corresponding crack region. The singular behavior presented by the stress field near the crack tip is modeled by traction singular quarter point element. Numerical examples of problems with in-plane loading are presented. Mesh convergence analyses are also presented, showing little dependence on the discretization even when coarse meshes were used. Some of these examples have correspondents in literature, that were used for comparisons with the obtained results. A good agreement in the comparisons of results was observed. / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica

Este trabalho apresenta o desenvolvimento e implementação computacional de técnicas de otimização de topologia para problemas governados pela equação de Poisson. O método numérico utilizado para solução numérica das equações foi o método dos elementos de contorno (MEC). Para tanto, três metodologias foram desenvolvidas. A primeira é direcionada à aplicação de algoritmos genéticos (AG) para investigar como um domínio inicialmente preenchido com cavidades aleatórias evolui durante um processo de otimização e verificar a possibilidade de se extrair topologias ótimas a partir da interpretação da solução encontrada. Os contornos externos permanecem fixos enquanto as posições e as dimensões das cavidades são otimizadas com o objetivo extremizar uma função custo especificada. O desempenho do algoritmo proposto é ilustrada com uma série de exemplos e os resultados são discutidos. A segunda metodologia apresenta um algoritmo numérico para otimização topológica baseado na avaliação da derivada topológica (DT), adotando a energia potencial total como função custo. Este procedimento é uma alternativa às tradicionais técnicas de otimização, evitando assim soluções de projeto com densidade de material intermediária. Sólidos com comportamento anisotrópico são estudados sob condições de contorno de Robin, Neumann e Dirichlet. Uma transformação linear de coordenadas é utilizada para mapear o problema original e suas condições de contorno para um novo domínio equivalente isotrópico, onde o procedimento de otimização é aplicado. A solução otimizada é então transformada de volta ao domínio original. A metodologia proposta mostrou-se particularmente atrativa para resolver esta classe de problemas já que o MEC dispensa o uso de malha no domínio, reduzindo significantemente o custo computacional. Na última parte deste trabalho foi implementada uma formulação de sensibilidade topológica para problemas de otimização de transferência de calor e massa simultâneos. Como as sensibilidades para cada equação diferencial são diferentes, utiliza-se um coeficiente de ponderação para compor a sensibilidade do problema acoplado. Isto permite a imposição de distintos fatores para cada problema, de acordo com uma prioridade especificada. Diversos exemplos são apresentados e seus resultados comparados com os da literatura, quando disponíveis, a fim de validar as formulações propostas. / This work presents the computational development and implementation of topology optimization techniques for problems governed by the Poisson equation. The boundary element method was the numerical technique chosen to solve the equations. Three different methodologies were developed aiming this objective. The first methodology is directed to the application of genetic algorithms to investigate how a domain previously populated with randomly placed cavities evolves during the optimization process, and to verify the resemblance of the final solution with a optimal design. The external boundaries remain fixed during the process, while the location and dimension of the cavities are optimized in order to extremize a given cost function. The performance of the proposed algorithm is verified with a number of examples and the results are discussed. The second methodology presents a numerical algorithm for topology optimization based on the evaluation of topological derivatives, using the total potential energy as the cost function. This procedure is an alternative to the traditional optimization techniques, avoiding design solutions containing intermediary material densities. Solids with anisotropic constitutive behavior are studied under Robin, Neumann and Dirichlet boundary conditions. A linear coordinate transformation approach is used to map the original problem into an isotropic one, where the optimization is carried out. The final solution is then mapped back to the original coordinate system. The proposed method was found to be an attractive way to solve this class of problems, since no interior mesh is necessary, which reduces significantly the computational cost of the analysis. In the last part of the present work the topological derivative approach was further developed to deal with the optimization of problems under simultaneous heat and mass transfer. Since the sensitivities for each differential equation are different, a weighting factor was used to evaluate the final sensitivities of the coupled problem. This allows the imposition of different priorities for each problem Several examples are presented and their results are compared with the literature, when available, in order to validate the proposed formulations.

Orientador: Paulo Sollero / Tese (doutorado) - Universidade Estadual de Campinas. Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-06T09:16:06Z (GMT). No. of bitstreams: 1 Noritomi_PedroYoshito_D.pdf: 4114284 bytes, checksum: 806cfb9e9589c5cd848086ed60178957 (MD5) Previous issue date: 2005 / Resumo: Este trabalho mostra o desenvolvimento de uma metodologia para análise de problemas de bioengenharia, aplicando modelagem numérica elastostática de tensões e deformações, baseada no método dos elementos de contorno com formulação 3D para meios transversalmente isotrópicos lineares, incluindo a capacidade de simulação do comportamento de remodelagem óssea superficial. A implementação do núcleo transversalmente isotrópico sobre a estrutura básica de análise por elementos de contorno 3D utilizou a solução fundamental proposta por Pan & Chou e revisada por Loloi, tendo exigido o cálculo adicional das soluções fundamentais de força de superfície a partir da derivação das soluções fundamentais de deslocamento. O modelo de remodelagem óssea superficial baseou-se na hipótese de estímulo biológico por campo de deformação, partindo de um modelo 2D, adaptado para o espaço 3D com o uso de deformações principais como grandezas de referência. As implementações foram testadas através de análises numéricas de problemas com solução analítica e validações com resultados de aplicações comerciais baseadas em elementos finitos, para problemas padrão de engenharia, bem como comparações com resultados da literatura para problemas de bioengenharia. A análise dos resultados mostrará que, tanto a metodologia quanto as implementações são funcionais, oferecendo uma base sólida para desenvolvimento e teste de novas soluções de bioengenharia / Abstract: This work shows the development of a methodology to analyse bioengineering problems using elastostatic stress-strain numerical modeling based on a 3D transversely isotropic linear boundary element formulation including surface bone remodeling simulation capabilities. The transversely isotropic kernel implementation on the basic 3D boundary element analysis program used the fundamental solution purposed by Pan & Chou and revised by Loloi, with additional fundamental solutions for traction calculation made with the displacement fundamental solution derivatives. The surface bone remodeling model was based on a 2D strain field biological stimulus, extended to the 3D space by using the principal strain as reference values. The implementations were tested through numerical analysis of problems with analytical solution and validation with commercial finite elements applications for standard engineering problems, as well as comparison with literature data for bioengineering problems. The analysis of results will show that both, the methodology and the implementations are fully functional, offering a solid start for development and test of new bioengineering solutions / Doutorado / Mecanica dos Sólidos e Projeto Mecanico / Doutor em Engenharia Mecânica

Orientadores: Eder Lima de Albuquerque, Alberto Luiz Serpa / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-06T09:57:35Z (GMT). No. of bitstreams: 1 Gouvea_AdrianadosReis_M.pdf: 936019 bytes, checksum: 2509e5d6c5977816a1830afe11bcd691 (MD5) Previous issue date: 2006 / Resumo: Neste trabalho é apresentada uma análise, através do uso do método dos elementos de contorno, de problemas estáticos em um modelo de estado plano para materiais anisotrópicos. A formulação de elementos de contorno é desenvolvida usando soluções fundamentais para elasto-estática e são utilizados apenas elementos quadráticos contínuos (3 nós por elementos). A formulação desenvolvida é aplicada no cálculo de tensões e deslocamentos de estruturas planas de materiais compósitos laminados submetidos a cargas no plano da estrutura. São estudados também os critérios de falha para materiais compósitos, abordando dois critérios quadráticos: o critério de Tsai-Hill e o critério de Tsai- Wu. São utilizados métodos de otimização para encontrar os valores ótimos dos ângulos de orientação das fibras que maximizem a resistência da estrutura. Para a otimização de problemas de uma variável é usado o método da Seção Áurea e para problemas de várias variáveis é usado o método Simplex Modificado. A função objetivo para os dois casos é o critério de falha. As variáveis de projeto são os ângulos de orientação das fibras nas diversas lâminas do laminado. Os resultados obtidos são comparados com resultados disponíveis na literatura e mostram uma boa concordância / Abstract: This work presents an analysis, using the boundary element method, of static problems in a plain state model in anisotropic materiaIs. The boundary element formulation is developed using basic solutions for elasto-statics and only continuous quadratic elements (three nodes per element) are used. The developed formulation is applied to compute the stress and displacements of laminated composite structures plain submitted to loads in the structure plane. It is studied also the failure criteria for composite materials, considering two quadratic criteria: the Tsai-Hill criterion and the Tsai-Wu criterion. Optimization methods are used to find optimum values for the fibres orientation angles that maximize the strength of the structure. For optimization problems of one variable, the Golden Section method is used and for problems of several variables the Modified Simplex method is used. The objective function for both cases is the failure criteria. Design variables are fibres orientation angles of all lamina of the laminate. Results are compared with other results from literature and show a good agreement / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica

Orientador: Leandro Palermo Jr / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Civil, Arquitetura e Urbanismo / Made available in DSpace on 2018-08-09T15:27:34Z (GMT). No. of bitstreams: 1 Sakanaka_SandraHiromi_M.pdf: 2300706 bytes, checksum: bf530823f185b6dbee7f56dd83d01a67 (MD5) Previous issue date: 2006 / Resumo: A análise de vibração livre e de instabilidade de placas finas e placas moderadamente espessas é apresentada através do método dos elementos de contorno (MEC) considerando o efeito da deformação pela força cortante e, particularmente para o cálculo de freqüências naturais, o efeito da inércia rotatória é também considerado. A formulação da solução fundamental é baseada na teoria de Mindlin (1951) mas resultados para a teoria de Kirchhoff (1850) também podem ser obtidos [Palermo Jr. (2000)]. O presente trabalho usa a técnica da iteração inversa através do coeficiente de Rayleigh para a determinação das menores freqüências naturais e cargas críticas de instabilidade das placas. A implementação numérica emprega elementos de contorno isoparamétricos lineares contínuos e descontínuos. Elementos constantes de domínio são usados. Os parâmetros nodais são posicionados nos extremos dos elementos e os pontos de carregamento dos elementos descontínuos são deslocados para o interior a uma distância igual a um quarto do comprimento do elemento. Expressões analíticas das integrais de contorno são desenvolvidas para os casos em que o elemento contém o ponto de carregamento e integração numérica de Gauss-Legendre é feita nos outros casos. As integrais de domínio foram transformadas em integrais de contorno para cada célula e foram tratadas como cargas de superfície atualizadas através de um processo iterativo. Os resultados obtidos foram comparados com valores encontrados na literatura para demonstrar a precisão do presente trabalho / Abstract: Free-vibration analysis and static buckling loads analysis of thin and thick plates considering the shear deformation effects using the Boundary Element Method (BEM) is presented. For the calculation of natural frequencies, the rotatory inertia is also counted. The formulation of the fundamental solution considers Mindlin¿s plates but results according to the classic theory can also be obtained [Palermo Jr. (2000)]. The present article makes use of the inverse iteration with Rayleigh coefficient to determine the smallest natural frequencies and the smallest static buckling loads of the plates. The numerical implementation employed continuous or discontinuous isoparametric linear boundary elements according to the characteristics of the problem to be solved. Constant domain elements are used. Nodal parameters have been placed at the ends of the elements and the source point of the discontinuous elements were positioned at a distance equal to one quarter of the element length. Analytical expressions have been employed in the integration on elements containing the source point and Gauss-Legendre numerical integration scheme otherwise. The domain integrals containing the inertia effects or nonlinear effect have been transformed into boundary integrals for each cell and were treated as surface loads updated in an iterative process. The obtained results were compared to those in literature to demonstrate the precision of this proposal / Mestrado / Estruturas / Mestre em Engenharia Civil

Orientadores: Eder Lima de Albuquerque, Leandro Palermo Junior / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-15T14:43:24Z (GMT). No. of bitstreams: 1 Reis_Adrianados_D.pdf: 2025685 bytes, checksum: f44ef47ff8671c2244a0b67e8e7c3fe4 (MD5) Previous issue date: 2010 / Resumo: Neste trabalho é apresentada uma análise, através do método dos elementos de contorno, de problemas estáticos em placas anisotrópicas sob flexão. As formulações dos elementos de contorno é desenvolvida usando soluções fundamentais para elasto-estática e são utilizados apenas elementos quadráticos descontínuos (3 nós por elementos). As integrais de superfície provenientes das cargas de domínio, são transformadas em integrais de contorno usando o método da integração radial. São estudadas três formulações: a formulação clássica ou de placas finas, a formulação de três parâmetros e a formulação de placas espessas ou placas com efeito da deformação por cortante ou teoria de primeira ordem. São apresentados detalhes das deduções das soluções fundamentais para a teoria clássica e para a teoria de primeira ordem. A solução fundamental da teoria clássica é obtida usando variáveis complexas. Para o caso da teoria de primeira ordem, o operador de Hörmander é usado para transformar os sistema de equações diferenciais parciais que representa as equações de equilíbrio em apenas uma equação diferencial parcial. Usando a transformada de Radon, essa equação diferencial parcial é reduzida a uma equação diferencial ordinária. A formulação desenvolvida é aplicada no cálculo de deslocamentos, tensões e momentos de estruturas planas de materiais compósitos laminados submetidos a cargas distribuídas no domínio da estrutura. Os resultados obtidos são comparados com resultados analíticos disponíveis na literatura, com o método de elementos finitos e com o método sem malhas. Em geral há uma boa concordância com os resultados da literatura / Abstract: This work presents an analysis, using the boundary element method, of static problems of anisotropic plate bending. Boundary element formulations are developed using fundamental solutions for elasto-statics and only discontinuous quadratic elements (three nodes per element) are used. Surface integrals that come from domain loads are transformed into boundary integrals using the radial integration method. Three formulations are studied: the classic formulation or of thin plate theory, the formulation of three parameters, and the formulation of thick plates, or plate with the effect of transversal shear deformation or first order theory. Details of the derivation of the fundamental solutions are presented to the classical and first order theory. The fundamental solution of the classical theory are obtained using complex variables. In the case of thick plates, the Hörmander operator is used to trans- form the partial differencial equation system that represents the equilibrium equation in only one partial differential equation. Using Radon transform, this partial differential equation is reduced to an ordinary differential equation. The developed formulations are applied to compute displacements, stresses and moments of laminated composite structures submitted to transversal loads in the structure domain. Results are compared with analytical results available in literature, with finite element method, and meshless method. In general there is a good agreement between the results obtained in this work and those available in literature / Doutorado / Mecanica dos Sólidos e Projeto Mecanico / Doutor em Engenharia Mecânica

Apresentam-se neste trabalho formulações do método dos elementos de contorno (MEC), visando sua utilização em problemas de fratura e também de modelos de confiabilidade e otimização aplicados na análise de problemas de fadiga. Com relação aos progressos e avanços nas formulações do MEC, apresentam-se modelos que representam o processo de crescimento de fissuras em domínios planos constituídos por materiais frágeis, quase-frágeis e dúcteis. Considerando esses diferentes tipos de materiais, a formulação numérica adotada na análise descreve o comportamento estrutural não linear decorrente do processo de propagação das fissuras e conseqüente degradação estrutural. Nos modelos de fratura é empregada a formulação MEC dual, a qual é mais adequada para a análise da propagação aleatória de fissuras. São também apresentadas as expressões dos operadores tangente para as formulações não lineares que tratam os problemas de fratura elástico linear e coesiva, problemas de contato e os problemas de domínios enrijecidos. Com relação às análises de confiabilidade estrutural, o modelo mecânico de fadiga é acoplado a algoritmos de confiabilidade para a determinação do índice de confiabilidade e do conjunto de valores aleatórios com maior probabilidade de ocorrência. São testados alguns algoritmos de confiabilidade, podendo-se claramente definir um deles como mais eficiente para a análise de problemas de fadiga. A esse modelo é acoplado um algoritmo de otimização para a determinação das dimensões do elemento estrutural e dos intervalos para os procedimentos de manutenção e inspeção, que levam ao mínimo custo estrutural com base nas incertezas determinadas pelo modelo de confiabilidade. São apresentados vários exemplos validando e mostrando a eficiência das formulações desenvolvidas. / This work deals with the development of boundary element method (BEM) formulations to be used in engineering problems. Particular attention is given to using these formulations in development of reliability and optimization models applied to fatigue problems. Contributions to BEM formulations are developed, particularly, models that deal with crack growth in plane domains composed by brittle, quasi-brittle and ductile materials. Taking into account these different types of materials, the proposed formulation properly represents the nonlinear structural behaviour induced by crack growth and the resulting structural damage. The dual BEM formulation is adopted here for the proposed crack model and to analyse random crack propagation. In this thesis tangent operators are used in the non-linear BEM formulations, in order to deal with cohesive crack, contact problems and debonding problems in reinforced domains. Regarding structural reliability analysis, the fatigue mechanical model was coupled with appropriate reliability algorithms to compute the reliability index and other important random values. Several reliability algorithms were tested for this coupled model, in order to find the most efficient in the analysis of fatigue problems. An optimization model was also coupled with the fatigue reliability model, in order to evaluate the optimal structural element dimensions and also to schedule the intervals for maintenance and inspection procedures, taking into account the minimum cost and problem uncertainties. Many examples are presented in order to show the efficiency and accuracy of the proposed formulations in dealing with crack propagation, fatigue reliability analysis and optimization problems.

O presente trabalho desenvolve uma formulação do Método dos Elementos de Contorno para análise de problemas tridimensionais de fraturamento no regime transiente. Utilizam-se as soluções fundamentais da elastostática para obter a matriz de massa, empregando-se o Método da Reciprocidade Dual e a discretização do domínio por células tridimensionais. Para a integração no tempo são utilizados os algoritmos de Newmark e Houbolt. O fenômeno do fraturamento é abordado através da consideração de um campo de tensões iniciais, introduzindo-se o conceito de dipolos de tensão. Os tensores desenvolvidos que se relacionam aos dipolos, derivados das soluções fundamentais, são também apresentados. É utilizado o modelo de fratura coesiva. O contorno é discretizado utilizando-se elementos triangulares planos com aproximação linear, e elementos constantes para a superfície fictícia de fraturamento. São feitas várias aplicações cujos resultados obtidos confirmam a importância e a adequação da formulação apresentada para os problemas propostos. / This work presents a Boundary Element Method (BEM) formulation for analysis of three-dimensional fracture mechanics transient problems. Elastostatics fundamental solutions are considered in order to obtain the mass matrix, using both Dual Reciprocity Method and three-dimensional cell discretization. Newmark and Houbolt algorithms are employed to evaluate the time integrals. The fracture effects are captured by using dipoles of stresses, derived from an initial stress field. The tensors related to those dipoles, developed in the present work, are presented. The cohesive crack is the adopted model. Body boundary is discretized though linear flat triangular elements and the fracture surfaces are approximated by constant flat triangular elements. Some applications are processed to show the efficiency of presented BEM formulations.