[pt] Em análises efetuadas através de Métodos de Elementos de Contorno, o procedimento de integração exerce um papel fundamental, já que faz-se necessária a avaliação de integrais singulares e quase-singulares que introduzem erros numéricos nos resultados quando não são devidamente avaliadas. Nos últimos anos, vários pesquisadores sugeriram diferentes técnicas de integração para tratar os problemas de integração de uma forma adequada. Este trabalho inicia apresentando alguns conceitos básicos e uma revisão bibliográfica das principais técnicas sugeridas. Em seguida, apresenta-se uma técnica de integração unificada, que possui uma forma simples e oferece resultados com excelente precisão. A técnica proposta foi aplicada para integrais singulares ou quase- singulares possuindo pólos simples ou múltiplos, sendo que tanto integrais unidimensionais quanto integrais bidimensionais foram consideradas.Paralelamente ao estudo das integrais, foi desenvolvido um programa computacional em linguagem orientada a objetos (C++), destinado a implementar simultaneamente a formulação convencional e as formulações híbridas dos Métodos de Elementos de Contorno. Da forma em que foi planejado, a implementação suporta diferentes aplicações de análises de engenharia. Este programa resultou de um trabalho conjunto realizado com pesquisadores da Universidade de Stuttgart. Por fim, apresentam-se diversos exemplos numéricos e resultados de análises, ressaltando o bom desempenho da técnica proposta e a influência do procedimento de integração em análises através de Métodos de Elementos de Contorno. / [en] In Boundary Element Method analysis, the integration procedure is one of the most important tasks, since one has to deal with singular and quasi-singular integrals which introduce numerical errors in the results, if they are not evaluated adequately. In the last years, several researchers have suggested different techniques with the aim of handling the problem adequately. This work begins presenting some basic concepts and a review of the most important work published before. Following, it introduces a unified integration technique which has a simple form and provides highly accurate results. The proposed scheme also deals with one- or two-dimensional singular or quasi- singular integrals having single or multiple poles. Besides the study of the integrals, a computational code was developed using an objectoriented computer language (C++). This code takes into account the conventional and hybrid formulations of the Boundary Element Method and supports different types of engineering analysis. This computer program was developed in a frame of a joint project with some researchers from the University of Stuttgart.Finally, several numerical examples and analysis results are displayed, showing the good performance of the proposed technique and the influence of the integration task in analysis using Boundary Element Methods.

[pt] O presente trabalho tem como principais objetivos estudar a formulação do método dos elementos de contorno e implementar um programa computacional para análise de tensões de problemas bi-dimensionais (estado plano de deformação e axissimetria) considerando também a possibilidade de acoplamento dos métodos dos elementos de contorno (MEC) e dos elementos finitos (MEF). Dentro deste estudo são abordados as soluções fundamentais para materiais linearmente elásticos, a formulação das integrais no contorno, as técnicas para tratamento de singularidade, a utilização de nós duplos para estudos de canto, etc. Os resultados numéricos obtidos em alguns problemas bi-dimensionais pelos métodos dos elementos de contorno e dos elementos finitos, considerando isoladamente ou acoplados, são comparados mutuamente e com soluções analíticas da teoria da elasticidade linear. Vantagens e desvantagens destes métodos, bem como dificuldades de implementação numérica, são discutidas ao longo deste trabalho. / [en] The main objectives of this work are the study of the boundary element formulation as well as the implementation of computer program for stress analyses of bidimensional problems under axisymmetric or plane strain conditions. The computer program also combines the boundary element method (BEM) with the finite element method (FEM) thus permitting its application to wide range of geotechnical problems. In this study research several mathematical aspects of the boundary element method are reviewed, such as the fundamental solutions for elastic materials, contour integration, singularities, corner problems, etc. The numeric results obtained in some 2D problems the BEM and FEM, in a coupled or isolated form, are compared with analytical solutions provided by classical theory of elasticity. Advantages and shortcomings of both methods, as well as the difficulties in the computer implementation, are also discussed in this work.

[pt] Mais de três décadas atrás, Przemieniecki introduziu uma formulação para análise de elementos de barra e treliça baseada em uma expansão em série de freqüências. Recentemente esta formulação foi generalizada para análise de sistemas elásticos submetidos a carregamento qualquer e deslocamentos iniciais. Baseado no método da superposição modal, um sistema acoplado, com equações diferenciais de movimento de alta ordem, é transformado em um sistema desacoplado com equações diferenciais de segunda ordem, que pode ser resolvido por qualquer método conhecido na literatura. A motivação para este desenvolvimento é o Método Híbrido dos Elementos de Contorno, que tem sido desenvolvido para problemas dependentes do tempo e problemas dependentes da freqüência. Esta formulação, assim como a introduzida por Pian para o Método dos Elementos Finitos, obtém uma matriz de rigidez utilizando apenas integrais de contorno, para um domínio de forma qualquer contendo vários graus de liberdade. O uso de termos com freqüências de alta ordem melhora muito a precisão numérica. A análise modal de um problema dinâmico, conforme se apresenta, é aplicável a qualquer formulação de elementos finitos, em geral, desde que a matriz de rigidez generalizada possa ser obtida. Este trabalho é uma tentativa de consolidação da formulação teórica proposta, em que se faz uso de integrais exclusivamente no contorno, com a discussão de diversos casos particulares e a conseqüente avaliação numérica: estruturas restringidas ou não; consideração de deslocamentos e velocidades iniciais, tanto em termos de valores nodais quanto de campos prescritos no domínio (incluindo deslocamentos de corpo rígido); deslocamentos forçados dependentes do tempo; forças de massa dependentes do tempo; cálculo de resultados em pontos internos. Vários exemplos acadêmicos para problemas de potencial bidimensionais ilustram este trabalho. / [en] More than three decades ago, Przemieniecki introduced a formulation for the free vibration analysis of bar and beam elements based on a power series of frequencies. Recently, this formulation was generalized for the analysis of the dynamic response of elastic systems submitted to arbitrary nodal loads as well as initial displacements. Based on the mode-superposition method, a set of coupled, higher-order differential equations of motion is transformed into a set of uncoupled second order differential equations, which may be integrated by means of standard procedures. Motivation for this theoretical achievement is the hybrid boundary element method, which has been developed for time-dependent as well as frequency-dependent problems. This formulation, as a generalization of Pian`s previous achievements for finite elements, yields a stiffness matrix for which only boundary integrals are required, for arbitrary domain shapes and any number of degrees of freedom. The use of higher-order frequency terms drastically improves numerical accuracy. The introduced modal assessment of the dynamic problem is applicable to any kind of finite element for which a generalized stiffness matrix is available. The present work is an attempt of consolidating this boundary- only theoretical formulation, in which a series of particular cases are conceptually outlined and numerically assessed: Constrained and unconstrained structures; initial displacements and velocities as nodal values as well as prescribed domain fields (including rigid body movement); forced time-dependent displacements; time-dependent body forces; evaluation of results at internal points. Several academic examples for 2D problems of potential illustrate the formulation.

[pt] Este trabalho trata o Método Híbrido dos Elementos de Contorno com vista à análise de problemas que envolvam simetria ou antissimetria. Nestes casos, apenas uma parte da estrutura, que pode ser a metade, um quarto ou um oitavo, deve ser discretizada e capaz de representar o todo. Os métodos de contorno apresentam a vantagem, quando comparados com os de domínio, de não ser necessário nenhum tipo de discretização ao longo dos eixos ou planos de simetria, sem a introdução de mais aproximações, visto que apenas o contorno é discretizado. Embora estas simplificações venham a restringir alguns deslocamentos de corpo rígido (para problemas de elasticidade), no Método dos Elementos de Contorno convencional (colocação ou Galerkin) a ausência de tais deslocamentos não acarreta alterações na sistemática do método. Nos Métodos Híbridos de Elementos de Contorno, por outro lado, os deslocamentos de corpo rígido são necessários direta ou indiretamente para a aplicação de condições de ortogonalidade e avaliação das propriedades espectrais que são essenciais na obtenção da diagonal principal de certas matrizes inerentes ao método, tais como de flexibilidade, de deslocamentos e de tensões. Esta necessidade de avaliação é uma característica de suma importância do método e, quando não houver possibilidade de fazê-la, deve-se procurar uma forma substituta conceitualmente equivalente. Verifica-se que, apesar de este método ser baseado em funções singulares de Green, é capaz de representar estados simples de tensões, tanto por trabalhos virtuais quanto por interpolações no domínio. Como objetivo principal deste trabalho, será demonstrado que para cada deslocamento de corpo rígido perdido, devido às restrições impostas pela simetria ou antissimetria, poderá ser utilizado um estado simples de tensão (constantes na maioria dos casos), que permitirá o estabelecimento de propriedades espectrais apropriadas. De forma a se garantir uma sistemática estruturada para o trabalho, faz-se uma abordagem de conceitos fundamentais aplicados a problemas da elastostática e potencial estacionário, na formulação variacional do Método Híbrido dos Elementos de Contorno com posteriores considerações especiais de estados simples de tensão (representados polinomialmente), para elasticidade tridimensional em geral, visto que para problemas bidimensionais o caso se torna uma particularização. Todas as combinações de simetria e antissimetria são avaliadas com a implementação numérica. Diversos exemplos de problemas bidimensionais ilustram a formulação teórica. / [en] The boundary element methods are suited for the analysis of symmetric and antisymmetric problems - in which only a part (half, quadrant or octant) of the structure needs to be explicitly considered - since, as an additional advantage when compared with a domain discretization method, no interpolation is required along the symmetry axes (for 2D problems) or planes (for 3D problems) and, consequently, no approximations are introduced thereon. Although such computational simplification may prevent some of the structures allowable rigid body movements (elasticity problems considered), this fact may be completely ignored as concerning the implementation of the traditional (collocation or Galerkin) boundary element methods. In the hybrid boundary element methods, on the other hand, special orthogonality conditions, directly or indirectly related to rigid body displacements, are required for the evaluation of elements about the main diagonal of some matrices (flexibility, displacement and stress matrices). Then, a central issue in such methods is the assessment of these matrices spectral properties for any combination of symmetry and antisymmetry and, most important, the investigation of conceptually equivalent, substitutive properties. As presented in this work, the hybrid boundary element methods, although based on singular Green s functions, are able to simulate, in terms of both virtual work and field interpolation, the simplest stress states. Then, one demonstrates that for every missing rigid body displacement - brought about by some symmetry or antisymmetry consideration - one may lay hold of a simple (in most cases constant) stress state, which enables establishing appropriate spectral properties. This work introduces the underlying variational concepts of the hybrid boundary element method and outlines the special consideration of simple (polynomial) stress states, as generally formulated for 3D elasticity, since 2D elasticity and problems of potential may be dealt with as particular cases. All combinations of symmetry and antisymmetry are outlined with the aim of numerical implementation. A series of 2D examples for problems of potential illustrate the theoretical

[pt] Este trabalho faz parte de um projeto para a implementação de um programa que possa simular problemas com milhões de graus de liberdade em um computador pessoal. Para isto, combina-se o Método Fast Multipole (FMM) com o Método Expedito dos Elementos de Contorno (EBEM), além de serem utilizados resolvedores iterativos de sistemas de equações. O EBEM é especialmente vantajoso em problemas de complicada topologia, ou que usem funções fundamentais muito complexas. Neste trabalho apresenta-se uma formulação para o Método Fast Multipole (FMM) que pode ser usada para, virtualmente, qualquer função e também para contornos curvos, o que parece ser uma contribuição original. Esta formulação apresenta um formato mais compacto do que as já existentes na literatura, e também pode ser diretamente aplicada a diversos tipos de problemas praticamente sem modificação de sua estrutura básica. É apresentada a validação numérica da formulação proposta. Sua utilização em um contexto do EBEM permite que um programa prescinda de integrações sobre segmentos – mesmo curvos – do contorno quando estes estão distantes do ponto fonte. / [en] This is part of a larger project that aims to develop a program able to simulate problems with millions of degrees of freedom on a personal computer. The Fast Multipole Method (FMM) is combined with the Expedite Boundary Element Method (EBEM) for integration, in the project s final version, with iterative equations solvers. The EBEM is especially advantageous when applied to problems with complicated topology as well as in the case of highly complex fundamental solutions. In this work, a FMM formulation is proposed for the use with virtually any type of fundamental solution and considering curved boundaries, which seems to be an original contribution. This formulation presents a more compact format than the ones shown in the technical literature, and can be directly applied to different kinds of problems without the need of manipulation of its basic structure, being numerically validated for a few applications. Its application in the context of the EBEM leads to the straightforward implementation of higher-order elements for generally curved boundaries that dispenses integration when the boundary segment is relatively far from the source point.

[pt] Apresenta-se uma formulação particular do método híbrido dos elementos de contorno para a análise de problemas planos de potencial e de elasticidade que, apesar de completamente geral, é apropriada a aplicações de mecânica da fratura. Funções do tipo de Westergaard são usadas como soluções fundamentais, Em uma generalização de uma proposta inicialmente feita por Tada et al. A formulação leva a conceitos de elementos de contorno de deslocamentos semelhante à apresentada por Crouch e Starfield, mas em um contexto variacional que permite interpretações mecânicas bem simples das equações matriciais resultantes. Problemas de topologia geral podem ser modelados, como no caso de domínios infinitos ou multiplamente conexos. A formulação, que é diretamente aplicável a placas com entalhes ou trincas curvas externas ou internas, permite a descrição adequada de altos gradientes de tensão e é uma ferramenta simples de avaliação de fatores de intensidade de tensão, com o que se podem verificar numericamente conceitos estabelecidos por Rice em 1968. Esta dissertação tem foco na fundamentação matemática da formulação para problemas de potencial e de elasticidade. Apresenta-se a implementação da formulação e são discutidos vários exemplos numéricos de validação. / [en] A particular implementation of the hybrid boundary element method is presented for the two dimensional analysis of potential and elasticity problems, which although general in concept, is suited for fracture mechanics applications. Generalized Westergaard stress functions, as proposed by Tada et al in 1993, are used as the problem‘s fundamental solutions. The proposed formulation leads to displacement-based concepts that resemble those presented by Crouch and Starfield, although in a variational framework that leads to matrix equations with sound mechanical meanings. Problems of general topology, such as in the case of unbounded and multiply-connected domains, may be modeled. The formulation, which is directly applicable to notches and generally curved, internal or external cracks, is specially suited for the description of the stress field in the vicinity of crack tips and is an easy means of evaluating stress intensity factors and of checking some basic concepts laid down by Rice in 1968. This dissertation focuses on the mathematical fundamentals of the formulation. Several validating numerical examples are presented.

[pt] A modelagem matemática de microdispositivos, em que estrutura e microestrutura têm aproximadamente a mesma escala de magnitude, assim como de macroestruturas de natureza predominantemente granular ou cristalina, requer uma abordagem não-local de deformações e tensões. Há mais de cem anos os irmãos Cosserat já tinham desenvolvido uma teoria de grãos rígidos. No entanto, e sem detrimento de desenvolvimentos devidos a Toupin e outros pesquisadores, os trabalhos de Mindlin na década de 1960 podem ser considerados a base da chamada teoria gradiente de deformações, que se tornou recentemente objeto de um grande número de investigações analíticas e experimentais, motivadas pelo desenvolvimento de novos materiais estruturais e do crescente uso de dispositivos micro- e nanomecânicos na indústria. Mais recentemente, Aifantis e colaboradores conseguiram desenvolver uma teoria gradiente de deformações mais simplificada, com base somente em duas constantes elásticas adicionais, representativas de comprimentos característicos relacionados às energias de deformação superficial e volumétrica. Uma série de trabalhos recentes desenvolvidos por Beskos e colaboradores estendeu o campo de aplicações da proposta inicial de Aifantis e introduziu uma solução fundamental que de fato remonta aos trabalhos de Mindlin. A equipe de pesquisa de Beskos propôs as primeiras implementações 2D e 3D de elementos de contorno para análises de elasticidade gradiente tanto estáticas quanto no domínio da freqüência, inclusive para problemas da mecânica da fratura. Desde o tempo de Toupin e Mindlin procura-se estabelecer uma base variacional da teoria e uma formulação consistente das condições de contorno cinemáticas e de equilíbrio, o que parece ter tido êxito com os recentes trabalhos de Amanatidou e Aravas. Esta dissertação faz uma revisão da teoria gradiente da deformações e apresenta um estudo didático do problema mais simples que se possa conceber, que é o de uma barra sob diferentes tipos de ações axiais (Aifantis, Beskos). A solução fundamental para problemas 2D e 3D também é apresentada e estudada, tanto em termos de forças pontuais aplicadas, para uma implementação em termos de elementos de contorno, quanto de desenvolvimentos polinomiais (no caso estático), para implementação em termos de elementos finitos. Mostra-se que a teoria gradiente de deformação de Aifantis é adequada a uma formulação no contexto do potencial de Hellinger-Reissner, o que possibilita implementações híbridas de elementos finitos e de contorno. O presente trabalho de pesquisa objetiva o estudo do estado da arte no tema, com uma abordagem dos principais problemas de implementação computacional, inclusive em termos das integrais singulares que surgem. O desenvolvimento completo de programas de análise de elementos híbridos finitos e de contorno, para problemas estáticos e dinâmicos, está planejado para uma tese de doutorado em futuro próximo. / [en] The mathematical modeling of micro-devices in which structure and the microstructure are about the same scale of magnitude, as well as of macrostructure of markedly granular or crystal nature (microcomposites), demands a nonlocal approach for strains and stresses. More than one hundred years ago the Cosserat brothers had already developed a theory for rigid grains. However, and in no detriment due to Toupin and other researchers, Mindlin s work in the 1960s may be accounted the basis of the so-called strain gradient theory, which has recently become the subject of a large number of analytical and experimental investigations motivated by the development of news structural materials together with the increasing use of micro and nano-mechanical devices in the industry. More recently, Aifantis and coworkers managed to develop a simplified strain gradient theory based only on two additional elasticity constants that are representative of material lengths related to surface and volumetric strain energy. A series of very recent works done by Beskos and collaborators extended the field of applications of Aifantis propositions and introduced a fundamental solution that actually remounts to developments already laid down by Mindlin. Beskos workgroup may be regarded as the proponent of the first of the first boundary element 2D and 3D implementations on the subject for both statics and frequency-domain analyses, also including crack problems. Since Toupin and Mindlin`s time, investigations have been under development to establish the variational basis of the theory and to consistently formulate equilibrium and kinematic boundary conditions established by Amanatidou and Aravas. This dissertation makes a revision of the gradient strain elasticity theory and presents a didactic study of the simplest problem that can be conceived, i.e., a bar under different axial actions (Aifantis, Beskos). The fundamental solution for 2D and 3D problems is also presented and studied for an elastic medium submitted to a point force, for boundary methods developments, as well as submitted to polynomial stress fields (for static problems), as in the hybrid finite element method. It is shown that Aifantis strain gradient theory may be developed in the context of the Hellinger-Reissner potential, for the sake of hybrid finite and boundary element implementations. Goal of the present research work is as a detailed study of state art of the theme, which comprises an investigation of the singular integrals one must deal with in a computational implementation. The complete computational development for static and dynamic hybrid boundary/finite analyses is planned for a future doctoral thesis.

[pt] Este trabalho visa à implementação de um programa de elementos de contorno para problemas com milhões de graus de liberdade. Isto é obtido com a implementação do Método Fast Multipole (FMM), que pode reduzir o número de operações, para a solução de um problema com N graus de liberdade, de O(N(2)) para O(NlogN) ou O(N). O uso de memória também é reduzido, por não haver o armazenamento de matrizes de grandes dimensões como no caso de outros métodos numéricos. A implementação proposta é baseada em um desenvolvimento consistente do convencional, Método de colocação dos elementos de contorno (BEM) – com conceitos provenientes do Hibrido BEM – para problemas de potencial e elasticidade de larga escala em 2D e 3D. A formulação é especialmente vantajosa para problemas de topologia complicada ou que requerem soluções fundamentais complicadas. A implementação apresentada, usa um esquema para expansões de soluções fundamentais genéricas em torno de níveis hierárquicos de polos campo e fonte, tornando o FMM diretamente aplicável para diferentes soluções fundamentais. A árvore hierárquica dos polos é construída a partir de um conceito topológico de superelementos dentro de superelementos. A formulação é inicialmente acessada e validada em termos de um problema de potencial 2D. Como resolvedores iterativos não são necessários neste estágio inicial de simulação numérica, podese acessar a eficiência relativa à implementação do FMM. / [en] This work aims to present an implementation of a boundary element solver for problems with millions of degrees of freedom. This is achieved through a Fast Multipole Method (FMM) implementation, which can lower the number of operations for solving a problem, with N degrees of freedom, from O(N(2)) to O(NlogN) or O(N). The memory usage is also very small, as there is no need to store large matrixes such as required by other numerical methods. The proposed implementations are based on a consistent development of the conventional, collocation boundary element method (BEM) - with concepts taken from the variationally-based hybrid BEM - for large-scale 2D and 3D problems of potential and elasticity. The formulation is especially advantageous for problems of complicated topology or requiring complicated fundamental solutions. The FMM implementation presented in this work uses a scheme for expansions of a generic fundamental solution about hierarchical levels of source and field poles. This makes the FMM directly applicable to different kinds of fundamental solutions. The hierarchical tree of poles is built upon a topological concept of superelements inside superelements. The formulation is initially assessed and validated in terms of a simple 2D potential problem. Since iterative solvers are not required in this first step of numerical simulations, an isolated efficiency assessment of the implemented fast multipole technique is possible.

[pt] No método dos elementos de contorno tradicional, a modelagem numérica de trincas é usualmente realizada por meio de uma solução fundamental hipersingular. Um procedimento mais natural seria utilizar uma solução fundamental capaz de representar a singularidade 1/raiz quadrada r que surge quando se analisa o campo de tensões próximo à ponta da trinca. Esta representação já foi realizada por Dumont e Lopes em 2003, com alguns refinamentos conseguidos por Dumont e Mamani em 2011, numa formulação do Método Híbrido de Elementos de Contorno, onde as soluções fundamentais são desenvolvidas a partir de funções de tensão generalizadas do tipo Westergaard para problemas de trincas com deslocamento prescrito, conforme proposto por Tada et al, em 1993. O presente trabalho, que é uma continuação das pesquisas de Dumont e Mamani, realiza um estudo sobre o uso destas funções generalizadas para a representação de grandezas na ponta da trinca em problemas de elasticidade e potencial. Os resultados obtidos são comparados conceitualmente com os desenvolvimentos clássicos de Westergaard e Williams. Também foram analisados alguns resultados com funções de tensão generalizadas de trinca com abertura semielíptica e polinomiais, além do uso de funções que representam a rotação relativa das faces da trinca. Além disso, é apresentada a aplicação da função de tensão de Westergaard generalizada como solução fundamental do método dos Elementos de Contorno Convencional, mais especificamente para a obtenção da matriz G do sistema, uma vez que a matriz H já foi desenvolvida, em trabalhos anteriores, com bons resultados. São apresentados alguns exemplos numéricos de aplicação para contornos externos, furos e trincas. / [en] In the traditional boundary element methods, the numerical modelling of cracks is usually carried out by means of a hypersingular fundamental solution. A more natural procedure should make use of fundamental solutions capable of representing the 1/square root of r singularity that arises when one analyses the stress field near the crack tip. This representation has already been made by Dumont and Lopes in 2003, with some refinements achieved by Dumont and Mamani in 2011, in a formulation of the Hybrid Boundary Element Method, where the fundamental solutions are developed from Westergaard-type generalized stress functions for displacement-prescribed crack problems, as proposed by Tada et al in 1993. The present work, which is a continuation of research work done by Dumont and Mamani, carries out a detailed study on the use of these generalized functions for the representation of quantities at the crack tip for problems of elasticity and potential. The results obtained are compared conceptually with the classic developments by Westergaard and Williams. Some results were also analyzed with generalized stress functions for a crack with semi-elliptical opening shapes, besides the use of functions to represent the relative rotation of the crack faces. In addition, the application of the generalized Westergaard stress function as a fundamental solution in the Conventional Boundary Element Method is presented, more specifically for obtaining the matrix G of the system, since the matrix H has already been developed in previous works, with good results. Some numerical examples of application are presented for external boundary, holes and cracks.

[pt] Desde a década de 1990, o Método Fast Multipole (FMM) tem sido usado em conjunto com o Métodos dos Elementos de Contorno (BEM) para a simulação de problemas de grande escala. Este método utiliza expansões em série de Taylor para aglomerar pontos da discretização do contorno, de forma a reduzir o tempo computacional necessário para completar a simulação. Ele se tornou uma ferramenta bastante importante para os BEMs, pois eles apresentam matrizes cheias e assimétricas, o que impossibilita a utilização de técnicas de otimização de solução de sistemas de equação. A aplicação do FMM ao BEM é bastante complexa e requer muita manipulação matemática. Este trabalho apresenta uma formulação do FMM que é independente da solução fundamental utilizada pelo BEM, o Método Fast Multipole Generalizado (GFMM), que se aplica a elementos de contorno curvos e de qualquer ordem. Esta característica é importante, já que os desenvolvimentos de fast multipole encontrados na literatura se restringem apenas a elementos constantes. Todos os aspectos são abordados neste trabalho, partindo da sua base matemática, passando por validação numérica, até a solução de problemas de potencial com muitos milhões de graus de liberdade. A aplicação do GFMM a problemas de potencial e elasticidade é discutida e validada, assim como os desenvolvimentos necessários para a utilização do GFMM com o Método Híbrido Simplificado de Elementos de Contorno (SHBEM). Vários resultados numéricos comprovam a eficiência e precisão do método apresentado. A literatura propõe que o FMM pode reduzir o tempo de execução do algoritmo do BEM de O(N2) para O(N), em que N é o número de graus de liberdade do problema. É demonstrado que esta redução é de fato possível no contexto do GFMM, sem a necessidade da utilização de qualquer técnica de otimização computacional. / [en] The Fast Multipole Method (FMM) has been used since the 1990s with the Boundary Elements Method (BEM) for the simulation of large-scale problems. This method relies on Taylor series expansions of the underlying fundamental solutions to cluster the nodes on the discretised boundary of a domain, aiming to reduce the computational time required to carry out the simulation. It has become an important tool for the BEMs, as they present matrices that are full and nonsymmetric, so that the improvement of storage allocation and execution time is not a simple task. The application of the FMM to the BEM ends up with a very intricate code, and usually changing from one problem s fundamental solution to another is not a simple matter. This work presents a kernel-independent formulation of the FMM, here called the General Fast Multipole Method (GFMM), which is also able to deal with high order, curved boundary elements in a straightforward manner. This is an important feature, as the fast multipole implementations reported in the literature only apply to constant elements. All necessary aspects of this method are presented, starting with the mathematical basics of both FMM and BEM, carrying out some numerical assessments, and ending up with the solution of large potential problems. The application of the GFMM to both potential and elasticity problems is discussed and validated in the context of BEM. Furthermore, the formulation of the GFMM with the Simplified Hybrid Boundary Elements Method (SHBEM) is presented. Several numerical assessments show that the GFMM is highly efficient and may be as accurate as arbitrarily required, for problems with up to many millions of degrees of freedom. The literature proposes that the FMM is capable of reducing the time complexity of the BEM algorithms from O(N2) to O(N), where N is the number of degrees of freedom. In fact, it is shown that the GFMM is able to arrive at such time reduction without resorting to techniques of computational optimisation.