O trabalho tem como objetivo trazer contribuições à simulação numérica pelo método dos elementos de contorno (MEC) de formação e propagação de descontinuidades no campo de deslocamentos (descontinuidades fortes) em sólidos bidimensionais. A formação de descontinuidades fortes caracteriza o processo de falha material, que pode estar associado ao fraturamento em materiais quase frágeis ou a superfícies de deslizamentos de materiais dúcteis. Apresenta-se uma formulação do MEC baseada na incorporação de interfaces de descontinuidade no interior de células internas, que possibilita propagação arbitrária de descontinuidades usando uma malha de células internas fixa, definida antes da análise. Comparam-se diferentes alternativas provenientes do relaxamento dos requisitos de consistência estática e analisa-se a influência do alinhamento da malha. Apresenta-se também um possível esquema de construção adaptativa de células internas com interface incorporada para capturar a trajetória arbitrária da descontinuidade que se propaga durante o processo de carregamento. Este esquema visa aumentar a robustez e reduzir o esforço computacional. As características geométricas das células internas geradas são estabelecidas em função da orientação da descontinuidade fornecida pelo critério de falha, de maneira a proporcionar melhor eficiência numérica. Os estudos são levados a cabo através da simulação numérica de testes experimentais colhidos da literatura. / This work has the objective of bringing contributions to the numeric simulation using the boundary elements method (BEM) to model the initiation and propagation of strong discontinuities in the displacement field in bidimensional solids. The initiation process of strong discontinuities characterizes the failure process of material, which can be associated with the fracture of quasi-brittle materials and slip lines in ductile materials such as metals. The effect of the displacement jump of a discontinuity interface embedded in an internal cell is provided by an equivalent strain field over the cell. This model allows the study of arbitrary crack growth using a fixed mesh defined before the analysis. The dissipative process in the cell interface is described by an isotropic damage model in the continuum approach of strong discontinues. Alternatives that come from relaxing the static consistencies and the influence of the mesh alignment are analyzed. An adaptative algorithm for internal cells creation is also presented to capture the path of the crack growth during the loading process. This algorithm intends to overcome some convergence problems found in models with predefined meshes and also to reduce the computational efforts. The geometric characteristics of the generated internal cells are defined using the crack orientation, given by the failure criterion, to provide a better numerical efficiency. The results obtained with the proposed formulation are compared with the ones obtained with other numerical methods and also from experiments.

Tese (Doutorado) - Universidade Federal de Santa Catarina, Centro Tecnológico. / Made available in DSpace on 2012-10-17T13:30:56Z (GMT). No. of bitstreams: 0Bitstream added on 2014-09-25T19:15:18Z : No. of bitstreams: 1 161846.pdf: 91791177 bytes, checksum: f51a3f5dd4b7584432ec2bf16b4b082c (MD5) / Apresenta-se uma metodologia de modelagem global para caixas de engrenagem, tendo por base os métodos de elementos finitos e elementos de contorno, cujo enfoque principal é para as vibrações e ruídos transportados pela estrutura e que têm como origem a excitação provocada pelo engrenamento, de modo a permitir a predição de parâmetros vibro-acústicos. Assim, são apresentadas as formulações matemáticas que permitem a determinação de matrizes de rigidez generalizadas, para a representação de mancais de elementos rolantes e do engrenamento; o desenvolvimento de um modelo numérico de sistema de eixos engrenados, que permite estimar as reações nos mancais devido ao engrenamento e a influência da modificação do perfil dos dentes das engrenagens; a concepção e criação de uma malha de elementos finitos para representação do alojamento da caixa de engrenagem, bem como a malha acústica correspondente usada para a estimativa de parâmetros acústicos pelo método de elementos de contorno. A metodologia é aplicada a uma caixa de engrenagem de uso veicular real e o comportamento vibro-acústico é predito para condições operacionais distintas.

[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.

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.

Submitted by Elizabete Silva (elizabete.silva@ufes.br) on 2015-12-03T19:01:45Z No. of bitstreams: 2 license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) Avaliação do desempenho do procedimento recursivo do método dos elementos de contorno aplicado em problemas da elasticidade.pdf: 5162837 bytes, checksum: 7136087fd573dd0565f46bb109ead161 (MD5) / Approved for entry into archive by Morgana Andrade (morgana.andrade@ufes.br) on 2015-12-30T11:42:59Z (GMT) No. of bitstreams: 2 license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) Avaliação do desempenho do procedimento recursivo do método dos elementos de contorno aplicado em problemas da elasticidade.pdf: 5162837 bytes, checksum: 7136087fd573dd0565f46bb109ead161 (MD5) / Made available in DSpace on 2015-12-30T11:42:59Z (GMT). No. of bitstreams: 2 license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) Avaliação do desempenho do procedimento recursivo do método dos elementos de contorno aplicado em problemas da elasticidade.pdf: 5162837 bytes, checksum: 7136087fd573dd0565f46bb109ead161 (MD5) Previous issue date: 2015 / CNPQ / Com o intuito de se obter resultados mais precisos com menor custo computacional, novos procedimentos e técnicas vem sendo avaliadas e testadas exaustivamente no contexto dos mais importantes métodos numéricos. O procedimento recursivo se insere como uma dessas iniciativas que têm sido estudado para melhorar o desem-penho do Método dos Elementos de Contorno (MEC). Basicamente, o processo consiste na introdução, a posteriori, de novos pontos fonte sobre o contorno, com localização distinta dos pontos nodais originais, para neles se calcular valores da variável básica como se faz para determinar valores em pontos internos. O procedimento se baseia nos princípios do Método dos Resíduos Ponderados, de modo que a reaplicação da equação integral para pontos de contorno equivaleria a uma nova minimização de resíduos. A técnica foi testada anteriormente com êxito em problemas de campo escalar, e no presente trabalho é implementada junto aos problemas elásticos governados pela Equação de Navier, destacando particularmente sua aplicação ao recálculo das for-ças de superfície, pois nesse caso uma análise da formulação hipersingular do MEC precisou ser realizada, observando suas peculiaridades e cuidados. Vários problemas de elasticidade foram utilizados para que o procedimento pudesse ser comparado com o método direto e também com as soluções analíticas disponí-veis, para assim melhor avaliar a precisão da técnica proposta e concluir que o pro-cedimento é eficaz em certas situações. / In order to obtain more accurate results with lower computational cost, various pro-cedures and techniques have been evaluated and tested thoroughly. Among them, the recursive procedure is being studied for the boundary element method (BEM). The process consists of introducing a posteriori new source points on the contour, with distinct location of the original nodal points, in order to calculate values of the basic variable as is done to determine values at internal points. The procedure is based on the principles of Weighted Residual Method, so that re-application of the integral equation in the contour points is equivalent to a new resid-ual minimization. The technique was previously tested successfully in scalar field problems, and in this work is implemented with the elastic problems governed by the Navier equation. To use the recursive procedure, an analysis of the BEM hypersingular formulation was performed, noting their peculiarities and cautions. Several problems of elasticity were used so that the procedure could be compared with the direct method. We sought to evaluate the accuracy of the recursive procedure for the displacements and surface forces, showing that the procedure is effective in certain situations.

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.

Submitted by maria angelica Varella (angelica@sibi.ufrj.br) on 2018-02-20T18:07:28Z No. of bitstreams: 1 158139.pdf: 1149685 bytes, checksum: 165bbdaaaec313c6a4bbd28b6c09d515 (MD5) / Made available in DSpace on 2018-02-20T18:07:29Z (GMT). No. of bitstreams: 1 158139.pdf: 1149685 bytes, checksum: 165bbdaaaec313c6a4bbd28b6c09d515 (MD5) Previous issue date: 1983-03 / Apresenta uma análise da aplicação do Método dos Elementos de Contorno ao problema definido pela vibração forçada de um par de cilindros esbeltos imersos em uma região fluida infinita. Na formulação do problema foram utilizadas hipóteses simplificadoras, reduzindo-se o mesmo à determinação de um potencial complexo. Para a realização de tal análise, foi desenvolvida a solução analítica para o problema e, com ela, foram comparados os valores obtidos através da solução numérica acima mencionada. Para a comparação de tais soluções, foram utilizados os valores dos coeficientes de massa adicionais característicos do problema. / The purpose of this work is to analyse the use of the Boundary Element Method (BEM) for solving the problem defined by the forced vibration of a pair of slender circular cylinders immersed in aninfinite fluid domain. The formulation of the problem is done through the assumption of simplifying hipoteses, which reduce its solution to a complex potential determination. In arder to perform the above mentioned analysis, an analytic solution for the problem was developed and compared with the values obtained through the BEM numerical solution. In the comparing of those solutions, the values of the additional mass coeficients for the problem were used.

Made available in DSpace on 2015-05-08T14:59:51Z (GMT). No. of bitstreams: 1 arquivototal.pdf: 2509612 bytes, checksum: f63968d7b069e4d6a16ca4bafd49aa7b (MD5) Previous issue date: 2013-07-26 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES / The laminated composite beams are structural components which have aroused great interest in mechanical industry due to providing components of high strength-to-weight ratio, which make them particularly suitable to the automotive and aerospace industry. In recent decades the solutions of mathematical models for laminated beam have been established in a analytical or numerical forms. The latter have been built based on finite element method (FEM) philosophy. In this work a new solution based on Boundary Element Method (BEM) for both classical and for first order laminated beam theory for static loading is established. In addition, the BEM solution is extended to deal with harmonic loading under classic beams theory hypothesis. Discussions on mathematical steps to write down both integral equations and fundamental solutions for laminated beam problem are properly made. From the formulations here proposed, codes are implemented in C++, providing displacements and efforts at the boundary and domain of the beam. In addition, numerical results for typical cases are presented as well. / As vigas de compósitos laminados são componentes estruturais que têm despertado grande interesse na indústria mecânica, principalmente por prover componentes de alta razão resistência-peso, o que as tornam particularmente aplicável à indústria automobilística e aeronáutica. Nas últimas décadas as soluções dos modelos matemáticos de vigas laminadas têm sido estabelecidas na forma analítica ou numérica. Para o último caso, as soluções têm sido construídas na filosofia do método dos elementos finitos (MEF). Nesta dissertação uma nova solução baseada no Método dos Elementos de Contorno (MEC) para as teorias de vigas laminada clássica e de primeira ordem são estabelecidas para carregamentos estáticos. Além disso, a solução para o MEC é estendida para lidar com carregamento dinâmico harmônico sob a hipótese da Teoria Clássica de viga. Nas discussões sobre as etapas matemáticas são descritas as equações integrais e as soluções fundamentais para os problemas de vigas laminadas. A partir das formulações aqui propostas, códigos são implementados na linguagem C++, capaz de fornecer deslocamentos e esforços no contorno e no domínio da viga. Assim como, as frequências naturais. Além disso, os resultados numéricos, comparados com as soluções analíticas disponíveis na literatura, se mostraram de boa qualidade.

Neste trabalho apresenta-se um modelo numérico para a análise de problemas tridimensionais envolvendo a interação mecânica estaca-solo, acoplando-se o Método dos Elementos de Contorno (MEC) ao Método dos Elementos Finitos (MEF). O solo é modelado com o MEC utilizando-se as soluções fundamentais de Mindlin, assumindo um meio semi-infinito, homogêneo, isotrópico e elástico-linear. As estacas, modeladas com o MEF, consistem em um elemento único, com quatro nós e 14 parâmetros nodais (três deslocamentos em cada nó e mais duas rotações no topo da estaca). Cada uma das estacas é levada em consideração no MEC como uma linha de carga. Considera-se o escorregamento das estacas em relação ao maciço, empregando modelos de aderência para a definição da evolução das tensões tangenciais ao longo do comprimento das estacas. São empregados, como funções de forma, polinômios do quarto grau para os deslocamentos horizontais, cúbicos para os deslocamentos verticais e tensões horizontais ao longo do fuste e quadráticos para as tensões verticais do fuste e escorregamento. A reação da ponta da estaca é calculada assumindo tensão constante na base. / This work presents a method for tri-dimensional pile-soil interaction problems, by coupling the Boundary Element Method (BEM) to the Finite Element Method (FEM). The soil is modeled with BEM, using the Mindlins fundamental solutions, supposing a semi-infinite, homogeneous, isotropic, elastic and linear space. Piles are modeled with FEM and are represented by one element with four nodes and 14 nodal parameters (three displacements in each node and two rotations at the top node). Each pile is represented in BEM as a line load. The pile slip is considered using adherence models to evaluate the evolution of shaft tractions. There are employed fourth grade polynomial shape functions for horizontal displacements, cubic polynomial functions for vertical displacements and horizontal tractions along shaft and quadratic polynomial functions for vertical tractions and slip. Tip reaction is calculated supposing constant traction at the base.

Esta tese apresenta o desenvolvimento de um método de elementos de contorno (BEM) no domínio do tempo baseado em fontes de Rankine para analise linear de comportamento no mar de sistemas oceânicos. O método e formulado por dois problemas de valor inicial de contorno definidos para os potenciais de velocidade e aceleração, sendo este ultimo utilizado para calcular de maneira acurada a derivada temporal do potencial de velocidades. Testes de verificação são realizados para a solução dos problemas de difração, radiação e de corpo livre para flutuar. Uma vez verificada, a ferramenta e aplicada em dois problemas multicorpos considerados no estado-da-arte em termos de modelagem hidrodinâmica utilizando BEM. O primeiro trata do problema envolvendo duas embarcações atracadas a contrabordo. Este é um caso no qual os códigos baseados na teoria de escoamento potencial são conhecidos por apresentarem dificuldades na determinação das soluções, tendendo a superestimar as elevações de onda no vão entre as embarcações e a apresentar problemas de convergência numérica associados a efeitos ressonantes de onda. O problema e tratado por meio do método de damping lid e a convergência das series temporais e investigada avaliando diferentes níveis de amortecimento. Os resultados são comparados com dados experimentais. O segundo problema se refere a analise de sistemas multicorpos com grandes deslocamentos relativos. Neste problema, ferramentas no domínio da frequência nao podem ser utilizadas, por considerarem apenas malhas fixas. Deste modo, o presente método e estendido para considerar um gerador de malhas de paineis e um algoritmo de interpolação de ordem alta no laco de tempo do código, possibilitando a mudança de posições relativas entre os corpos durante a simulação. Os resultados são comparados com dados de experimentos executados especificamente para fins de verificação do código, apresentando uma boa concordância. De acordo com o conhecimento do autor, esta e a primeira vez que certas questões relativas a modelagem numérica destes dois problemas multicorpos são relatadas na literatura especializada em hidrodinâmica computacional. / The development of a time domain boundary elements method (BEM) based on Rankine\'s sources for linear seakeeping analysis of offshore systems is here addressed. The method is formulated by means of two Initial Boundary Value Problems defined for the velocity and acceleration potentials, the latter being used to ensure an accurate calculation of the time derivatives of the velocity potential. Verification tests for solving the difraction, radiation and free floating problems are presented. Once verified, the code is applied for two complex multi-body problems considered to be in the state-of-the-art for hydrodynamic modelling using BEM. The first is the seakeeping problem of two ships arranged in side-by-side, a problem in which all potential flow codes are known to have a poor performance, tending to provide unrealistic high wave elevations in the gap between the vessels and to present numerical convergence problems associated to resonant effects. The problem is here addressed by means of a damping lid method and the convergence of the time series with different damping levels is investigated. Results are compared to data measured in an experimental campaign. The second problem refers to the analysis of multi-body systems composed of bodies undergoing large relative displacements. This is a case that cannot be properly analyzed by frequency domain codes, since they only consider fixed meshes. For this application, the present numerical method is extended to consider a panel mesh generator in the time loop of the code, enabling the change of body relative positions during the computations. Furthermore, a higher order interpolation algorithm designed to recover the solutions of a previous time-step was also implemented, enabling the calculations to progress with reasonable accuracy in time. The numerical results are compared to data of experimental tests designed and executed for verification of the code, and presented a very good agreement. To the author\'s knowledge, this is the first time that certain issues concerning the numerical modelling of these two complex multi-body problems are reported in the literature specialized in hydrodynamic computations.